OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-01-03 20:57:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

StimPlanModelPlot: Remove gaps in stress gradient and temperature curves

Browse Source
This commit is contained in:
Kristian Bendiksen 2021-03-23 14:56:02 +01:00 committed by GitHub
parent 5331e84d77
commit 53dbb33e86
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with 46 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

47
ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelStressCalculator.cpp
View File

@ -19,6 +19,7 @@
#include "RiaDefines.h" #include "RiaDefines.h"
#include "RiaEclipseUnitTools.h" #include "RiaEclipseUnitTools.h"
#include "RiaInterpolationTools.h"
#include "RiaLogging.h" #include "RiaLogging.h"
#include "RiaStimPlanModelDefines.h" #include "RiaStimPlanModelDefines.h"
@ -105,10 +106,12 @@ bool RimStimPlanModelStressCalculator::calculate( RiaDefines::CurveProperty curv
else if ( curveProperty == RiaDefines::CurveProperty::STRESS_GRADIENT ) else if ( curveProperty == RiaDefines::CurveProperty::STRESS_GRADIENT )
{ {
values = m_stimPlanModelCalculator->calculateStressGradient(); values = m_stimPlanModelCalculator->calculateStressGradient();
addDatapointsForBottomOfLayers( tvDepthValues, values );
} }
else if ( curveProperty == RiaDefines::CurveProperty::TEMPERATURE ) else if ( curveProperty == RiaDefines::CurveProperty::TEMPERATURE )
{ {
m_stimPlanModelCalculator->calculateTemperature( values ); m_stimPlanModelCalculator->calculateTemperature( values );
addDatapointsForBottomOfLayers( tvDepthValues, values );
} }
if ( eclipseCase ) if ( eclipseCase )
@ -119,15 +122,21 @@ bool RimStimPlanModelStressCalculator::calculate( RiaDefines::CurveProperty curv
// Generate MD data by interpolation // Generate MD data by interpolation
const std::vector<double>& mdValuesOfWellPath = wellPathGeometry->measuredDepths(); const std::vector<double>& mdValuesOfWellPath = wellPathGeometry->measuredDepths();
std::vector<double> tvdValuesOfWellPath = wellPathGeometry->trueVerticalDepths(); const std::vector<double>& tvdValuesOfWellPath = wellPathGeometry->trueVerticalDepths();
if ( mdValuesOfWellPath.empty() ) if ( mdValuesOfWellPath.empty() )
{ {
RiaLogging::error( "Well path geometry had no MD values." ); RiaLogging::error( "Well path geometry had no MD values." );
return false; return false;
} }
measuredDepthValues = // The thickness direction "fake" well path is always a straight line:
RigWellPathGeometryTools::interpolateMdFromTvd( mdValuesOfWellPath, tvdValuesOfWellPath, tvDepthValues ); // measured depth can be interpolated linearly
measuredDepthValues.clear();
for ( double tvd : tvDepthValues )
{
double md = RiaInterpolationTools::linear( tvdValuesOfWellPath, mdValuesOfWellPath, tvd );
measuredDepthValues.push_back( md );
}
CVF_ASSERT( measuredDepthValues.size() == tvDepthValues.size() ); CVF_ASSERT( measuredDepthValues.size() == tvDepthValues.size() );
} }
@ -166,3 +175,35 @@ void RimStimPlanModelStressCalculator::addDatapointsForBottomOfLayers( std::vect
stress = valuesWithBottomLayers; stress = valuesWithBottomLayers;
tvDepthValues = tvdWithBottomLayers; tvDepthValues = tvdWithBottomLayers;
} }
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RimStimPlanModelStressCalculator::addDatapointsForBottomOfLayers( std::vector<double>& tvDepthValues,
std::vector<double>& values )
{
std::vector<double> tvdWithBottomLayers;
std::vector<double> valuesWithBottomLayers;
for ( size_t i = 0; i < values.size(); i++ )
{
// Add the data point at top of the layer
double topLayerDepth = tvDepthValues[i];
double value = values[i];
tvdWithBottomLayers.push_back( topLayerDepth );
valuesWithBottomLayers.push_back( value );
// Add extra data points for bottom part of the layer
if ( i < values.size() - 1 )
{
double bottomLayerDepth = tvDepthValues[i + 1];
double bottomValue = value;
tvdWithBottomLayers.push_back( bottomLayerDepth );
valuesWithBottomLayers.push_back( bottomValue );
}
}
values = valuesWithBottomLayers;
tvDepthValues = tvdWithBottomLayers;
}

2
ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelStressCalculator.h
View File

@ -48,5 +48,7 @@ private:
std::vector<double>& stress, std::vector<double>& stress,
const std::vector<double>& stressGradients ); const std::vector<double>& stressGradients );
static void addDatapointsForBottomOfLayers( std::vector<double>& tvDepthValues, std::vector<double>& values );
RimStimPlanModelCalculator* m_stimPlanModelCalculator; RimStimPlanModelCalculator* m_stimPlanModelCalculator;
}; };