OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-01-24 15:26:48 -06:00
Code Issues Packages Projects Releases Wiki Activity

Wbs: Make sure last LAS-value is used if within a threshold

Browse Source
This commit is contained in:
Gaute Lindkvist 2020-05-26 11:37:52 +02:00
parent 08f3789c01
commit 920ebd081d
1 changed files with 16 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
ApplicationCode/ReservoirDataModel/RigGeoMechWellLogExtractor.cpp
View File

@ -103,7 +103,7 @@ void RigGeoMechWellLogExtractor::performCurveDataSmoothing( int
if ( !mds->empty() && !values->empty() ) if ( !mds->empty() && !values->empty() )
{ {
std::vector<std::vector<double>*> dependentValues = { tvds, &interfaceShValuesDbl, &interfacePorePressuresDbl }; std::vector<std::vector<double>*> dependentValues = {tvds, &interfaceShValuesDbl, &interfacePorePressuresDbl};
std::vector<unsigned char> smoothOrFilterSegments = determineFilteringOrSmoothing( interfacePorePressuresDbl ); std::vector<unsigned char> smoothOrFilterSegments = determineFilteringOrSmoothing( interfacePorePressuresDbl );
@ -296,7 +296,6 @@ std::vector<RigGeoMechWellLogExtractor::WbsParameterSource>
double waterDensityGCM3 = m_userDefinedValues[RigWbsParameter::waterDensity()]; double waterDensityGCM3 = m_userDefinedValues[RigWbsParameter::waterDensity()];
#pragma omp parallel for
for ( int64_t intersectionIdx = 0; intersectionIdx < (int64_t)m_intersections.size(); ++intersectionIdx ) for ( int64_t intersectionIdx = 0; intersectionIdx < (int64_t)m_intersections.size(); ++intersectionIdx )
{ {
// Loop from primary source and out for each value // Loop from primary source and out for each value
@ -1171,6 +1170,21 @@ double RigGeoMechWellLogExtractor::getWellLogIntersectionValue( size_t intersect
} }
} }
} }
// If we found no match, check first and last value within a threshold.
if ( !wellLogValues.empty() )
{
const double relativeEps = 1.0e-3 * std::max( 1.0, intersection_md );
if ( std::abs( wellLogValues.front().first - intersection_md ) < relativeEps )
{
return wellLogValues.front().second;
}
else if ( std::abs( wellLogValues.back().first - intersection_md ) < relativeEps )
{
return wellLogValues.back().second;
}
}
return std::numeric_limits<double>::infinity(); return std::numeric_limits<double>::infinity();
} }