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

WBS Plot: use PP (from any source) for shale for SFG.

Browse Source
This commit is contained in:
Kristian Bendiksen
2024-04-15 12:05:56 +02:00
parent 9dc5e99be0
commit b181c55db8
1 changed files with 10 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
ApplicationLibCode/ReservoirDataModel/RigGeoMechWellLogExtractor.cpp
View File

@@ -722,6 +722,10 @@ void RigGeoMechWellLogExtractor::wellBoreWallCurveData( const RigFemResultAddres
CAF_ASSERT( pp.size() == intersections().size() );
std::vector<double> ppShaleValues( intersections().size(), std::numeric_limits<double>::infinity() );
calculateWbsParameterForAllSegments( RigWbsParameter::PP_NonReservoir(), 0, 0, &ppShaleValues, true );
CAF_ASSERT( ppShaleValues.size() == intersections().size() );
#pragma omp parallel for
for ( int64_t intersectionIdx = 0; intersectionIdx < static_cast<int64_t>( intersections().size() ); ++intersectionIdx )
{
@@ -730,6 +734,12 @@ void RigGeoMechWellLogExtractor::wellBoreWallCurveData( const RigFemResultAddres
double hydroStaticPorePressureBar = hydroStaticPorePressureForSegment( intersectionIdx );
double porePressureBar = ppSandAllSegments[intersectionIdx];
if ( resAddr.fieldName == RiaResultNames::wbsSFGResult().toStdString() )
{
// SFG needs PP for shale.
porePressureBar = ppShaleValues[intersectionIdx] * hydroStaticPorePressureBar;
}
if ( porePressureBar == std::numeric_limits<double>::infinity() )
{
porePressureBar = hydroStaticPorePressureBar;