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

#6301 Minor tweak to match SFG calculation in WBS plot.

Browse Source
This commit is contained in:
Kristian Bendiksen 2020-08-24 13:51:03 +02:00
parent 5cd94eed47
commit fe6ece0df7
1 changed files with 11 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorMudWeightWindow.cpp
View File

@ -302,16 +302,15 @@ RigFemScalarResultFrames* RigFemPartResultCalculatorMudWeightWindow::calculate(
segmentStress ); segmentStress );
caf::Ten3d wellPathStressDouble( wellPathStressFloat ); caf::Ten3d wellPathStressDouble( wellPathStressFloat );
RigGeoMechBoreHoleStressCalculator sigmaCalculator( wellPathStressDouble,
porePressureBar,
poissonsRatio,
ucsBar,
32 );
// Calculate upper limit // Calculate upper limit
float upperLimit = inf; float upperLimit = inf;
if ( upperLimitParameter == RimMudWeightWindowParameters::UpperLimitType::FG && isSand ) if ( upperLimitParameter == RimMudWeightWindowParameters::UpperLimitType::FG && isSand )
{ {
RigGeoMechBoreHoleStressCalculator sigmaCalculator( wellPathStressDouble,
porePressureBar,
poissonsRatio,
ucsBar,
32 );
upperLimit = sigmaCalculator.solveFractureGradient() / hydroStaticPressureForNormalization; upperLimit = sigmaCalculator.solveFractureGradient() / hydroStaticPressureForNormalization;
} }
else if ( upperLimitParameter == RimMudWeightWindowParameters::UpperLimitType::SH_MIN ) else if ( upperLimitParameter == RimMudWeightWindowParameters::UpperLimitType::SH_MIN )
@ -351,6 +350,12 @@ RigFemScalarResultFrames* RigFemPartResultCalculatorMudWeightWindow::calculate(
} }
else else
{ {
RigGeoMechBoreHoleStressCalculator sigmaCalculator( wellPathStressDouble,
hydroStaticPressureForNormalization,
poissonsRatio,
ucsBar,
32 );
double SFG = sigmaCalculator.solveStassiDalia(); double SFG = sigmaCalculator.solveStassiDalia();
lowerLimit = std::max( porePressureBar, static_cast<float>( SFG ) ); lowerLimit = std::max( porePressureBar, static_cast<float>( SFG ) );
} }