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

StimPlanModel export: use linear interpolation as last fallback for EQLNUM

Browse Source
This commit is contained in:
Kristian Bendiksen 2021-05-28 13:32:39 +02:00 committed by Magne Sjaastad
parent 326a1765ef
commit 429d1e4166
4 changed files with 38 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModel.cpp
View File

@ -1111,7 +1111,8 @@ std::deque<RimStimPlanModel::MissingValueStrategy>
else if ( curveProperty == RiaDefines::CurveProperty::EQLNUM )
return { RimStimPlanModel::MissingValueStrategy::DEFAULT_VALUE,
RimStimPlanModel::MissingValueStrategy::CELLS_ABOVE,
RimStimPlanModel::MissingValueStrategy::CELLS_BELOW };
RimStimPlanModel::MissingValueStrategy::CELLS_BELOW,
RimStimPlanModel::MissingValueStrategy::LINEAR_INTERPOLATION };
else
return { RimStimPlanModel::MissingValueStrategy::DEFAULT_VALUE };
}

3
ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelPressureCalculator.cpp
View File

@ -209,6 +209,9 @@ bool RimStimPlanModelPressureCalculator::extractValuesForProperty( RiaDefines::C
return true;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::tuple<std::vector<double>, std::vector<double>, std::vector<double>>
RimStimPlanModelPressureCalculator::interpolateMissingValues( const std::vector<double>& targetTvds,
const std::vector<double>& targetMds,

39
ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelWellLogCalculator.cpp
View File

@ -127,7 +127,13 @@ bool RimStimPlanModelWellLogCalculator::calculate( RiaDefines::CurveProperty cur
{
// K-1 is up
int kDirection = -1;
if ( !replaceMissingValuesWithOtherKLayer( curveProperty, stimPlanModel, timeStep, measuredDepthValues, values, kDirection ) )
if ( !replaceMissingValuesWithOtherKLayer( curveProperty,
stimPlanModel,
timeStep,
measuredDepthValues,
tvDepthValues,
values,
kDirection ) )
{
return false;
}
@ -136,7 +142,13 @@ bool RimStimPlanModelWellLogCalculator::calculate( RiaDefines::CurveProperty cur
{
// K+1 is down
int kDirection = 1;
if ( !replaceMissingValuesWithOtherKLayer( curveProperty, stimPlanModel, timeStep, measuredDepthValues, values, kDirection ) )
if ( !replaceMissingValuesWithOtherKLayer( curveProperty,
stimPlanModel,
timeStep,
measuredDepthValues,
tvDepthValues,
values,
kDirection ) )
{
return false;
}
@ -561,6 +573,7 @@ bool RimStimPlanModelWellLogCalculator::replaceMissingValuesWithOtherKLayer( Ria
const RimStimPlanModel* stimPlanModel,
int timeStep,
const std::vector<double>& measuredDepths,
const std::vector<double>& tvDepthValues,
std::vector<double>& values,
int moveDirection ) const
{
@ -594,6 +607,7 @@ bool RimStimPlanModelWellLogCalculator::replaceMissingValuesWithOtherKLayer( Ria
if ( std::isinf( values[idx] ) )
{
double measuredDepth = measuredDepths[idx];
double tvd = tvDepthValues[idx];
cvf::Vec3d position = wellPathGeometry->interpolatedPointAlongWellPath( measuredDepth );
size_t cellIdx = mainGrid->findReservoirCellIndexFromPoint( position );
@ -606,11 +620,13 @@ bool RimStimPlanModelWellLogCalculator::replaceMissingValuesWithOtherKLayer( Ria
bool isValid = mainGrid->ijkFromCellIndex( cellIdx, &i, &j, &k );
if ( isValid )
{
RiaLogging::info( QString( "Missing value at MD: %1. Cell [%2, %3, %4]" )
.arg( measuredDepth )
.arg( i + 1 )
.arg( j + 1 )
.arg( k + 1 ) );
RiaLogging::debug(
QString( "K-Layer replacement: Replace missing value at MD: %1 TVD: %2. Cell [%3, %4, %5]" )
.arg( measuredDepth )
.arg( tvd )
.arg( i + 1 )
.arg( j + 1 )
.arg( k + 1 ) );
int neighborK = static_cast<int>( k ) + moveDirection;
const int minK = static_cast<int>( 1 );
@ -643,6 +659,15 @@ bool RimStimPlanModelWellLogCalculator::replaceMissingValuesWithOtherKLayer( Ria
}
}
}
else
{
RiaLogging::debug( QString( "K-Layer Replacement: Invalid cell idx: MD=%1 TVD=%2 Pos: [%3 %4 %5]" )
.arg( measuredDepth )
.arg( tvd )
.arg( position.x() )
.arg( position.y() )
.arg( position.z() ) );
}
}
}

1
ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelWellLogCalculator.h
View File

@ -100,6 +100,7 @@ protected:
const RimStimPlanModel* stimPlanModel,
int timeStep,
const std::vector<double>& measuredDepths,
const std::vector<double>& tvDepthValues,
std::vector<double>& values,
int kDirection ) const;