OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

#11034 Use water density from model.

Browse Source
This commit is contained in:
Kristian Bendiksen 2024-01-18 08:49:16 +01:00
parent 926a9ac850
commit a48faeb6d0
3 changed files with 6 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccess.cpp
View File

@ -105,6 +105,7 @@ RimFaultReactivationDataAccess::RimFaultReactivationDataAccess( const RimFaultRe
property,
porePressureGradient,
seabedDepth,
model.waterDensity(),
model.lateralStressCoefficientX(),
model.lateralStressCoefficientY(),
densities ) );

7
ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccessorStressEclipse.cpp
View File

@ -50,6 +50,7 @@ RimFaultReactivationDataAccessorStressEclipse::RimFaultReactivationDataAccessorS
RimFaultReactivation::Property property,
double gradient,
double seabedDepth,
double waterDensity,
double lateralStressComponentX,
double lateralStressComponentY,
const std::map<RimFaultReactivation::ElementSets, double>& densities )
@ -57,6 +58,7 @@ RimFaultReactivationDataAccessorStressEclipse::RimFaultReactivationDataAccessorS
, m_eclipseCase( eclipseCase )
, m_caseData( nullptr )
, m_mainGrid( nullptr )
, m_waterDensity( waterDensity )
, m_lateralStressComponentX( lateralStressComponentX )
, m_lateralStressComponentY( lateralStressComponentY )
, m_densities( densities )
@ -93,9 +95,6 @@ void RimFaultReactivationDataAccessorStressEclipse::updateResultAccessor()
m_wellPaths = wellPaths;
m_extractors = extractors;
// TODO: get from model
double waterDensity = 1030.0;
for ( auto [gridPart, wellPath] : m_wellPaths )
{
auto extractor = m_extractors[gridPart];
@ -104,7 +103,7 @@ void RimFaultReactivationDataAccessorStressEclipse::updateResultAccessor()
addOverburdenAndUnderburdenPoints( intersections, wellPath->wellPathPoints() );
m_stressValues[gridPart] =
integrateVerticalStress( *wellPath.p(), intersections, *m_model, gridPart, m_seabedDepth, waterDensity, m_densities );
integrateVerticalStress( *wellPath.p(), intersections, *m_model, gridPart, m_seabedDepth, m_waterDensity, m_densities );
}
}

2
ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccessorStressEclipse.h
View File

@ -47,6 +47,7 @@ public:
RimFaultReactivation::Property property,
double gradient,
double seabedDepth,
double waterDensity,
double lateralStressComponentX,
double lateralStressComponentY,
const std::map<RimFaultReactivation::ElementSets, double>& densities );
@ -89,6 +90,7 @@ private:
RigEclipseCaseData* m_caseData;
const RigMainGrid* m_mainGrid;
cvf::ref<RigResultAccessor> m_resultAccessor;
double m_waterDensity;
double m_lateralStressComponentX;
double m_lateralStressComponentY;