From fe6ece0df7369d1fa50ff202dff43b22710e5d19 Mon Sep 17 00:00:00 2001 From: Kristian Bendiksen Date: Mon, 24 Aug 2020 13:51:03 +0200 Subject: [PATCH] #6301 Minor tweak to match SFG calculation in WBS plot. --- ...igFemPartResultCalculatorMudWeightWindow.cpp | 17 +++++++++++------ 1 file changed, 11 insertions(+), 6 deletions(-) diff --git a/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorMudWeightWindow.cpp b/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorMudWeightWindow.cpp index d02318927d..6362dcaa6e 100644 --- a/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorMudWeightWindow.cpp +++ b/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorMudWeightWindow.cpp @@ -302,16 +302,15 @@ RigFemScalarResultFrames* RigFemPartResultCalculatorMudWeightWindow::calculate( segmentStress ); caf::Ten3d wellPathStressDouble( wellPathStressFloat ); - RigGeoMechBoreHoleStressCalculator sigmaCalculator( wellPathStressDouble, - porePressureBar, - poissonsRatio, - ucsBar, - 32 ); - // Calculate upper limit float upperLimit = inf; if ( upperLimitParameter == RimMudWeightWindowParameters::UpperLimitType::FG && isSand ) { + RigGeoMechBoreHoleStressCalculator sigmaCalculator( wellPathStressDouble, + porePressureBar, + poissonsRatio, + ucsBar, + 32 ); upperLimit = sigmaCalculator.solveFractureGradient() / hydroStaticPressureForNormalization; } else if ( upperLimitParameter == RimMudWeightWindowParameters::UpperLimitType::SH_MIN ) @@ -351,6 +350,12 @@ RigFemScalarResultFrames* RigFemPartResultCalculatorMudWeightWindow::calculate( } else { + RigGeoMechBoreHoleStressCalculator sigmaCalculator( wellPathStressDouble, + hydroStaticPressureForNormalization, + poissonsRatio, + ucsBar, + 32 ); + double SFG = sigmaCalculator.solveStassiDalia(); lowerLimit = std::max( porePressureBar, static_cast( SFG ) ); }