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

#6526 Fracture Model: Make TVD/MD values consistent in Deviation.FRK

Browse Source 
The change in measured depth (MD) must be greater or equal to change in
true vertical depth (TVD), and this was not always the case due to
interpolation imprecision. Fixed by adjusting TVD.
This commit is contained in:
Kristian Bendiksen 2020-09-18 14:17:02 +02:00
parent 8bb460c2db
commit fad6256d9e
1 changed files with 20 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
ApplicationCode/ReservoirDataModel/RigWellPathGeometryExporter.cpp
View File

@ -72,6 +72,9 @@ void RigWellPathGeometryExporter::exportWellPathGeometry( const RigWellPath* w
double currMd = wellPathGeom->measureDepths().front() - mdStepSize; double currMd = wellPathGeom->measureDepths().front() - mdStepSize;
double endMd = wellPathGeom->measureDepths().back(); double endMd = wellPathGeom->measureDepths().back();
bool isFirst = true;
double prevMd = 0.0;
double prevTvd = 0.0;
while ( currMd < endMd ) while ( currMd < endMd )
{ {
currMd += mdStepSize; currMd += mdStepSize;
@ -80,9 +83,26 @@ void RigWellPathGeometryExporter::exportWellPathGeometry( const RigWellPath* w
auto pt = wellPathGeom->interpolatedPointAlongWellPath( currMd ); auto pt = wellPathGeom->interpolatedPointAlongWellPath( currMd );
double tvd = -pt.z(); double tvd = -pt.z();
// The change in measured depth (MD) must be greater or equal to change in
// true vertical depth (TVD), and this was not always the case due to
// interpolation imprecision. Fix by adjusting TVD.
if ( !isFirst )
{
double deltaMd = currMd - prevMd;
double deltaTvd = tvd - prevTvd;
if ( deltaMd < deltaTvd )
{
tvd = prevTvd + deltaMd;
}
}
xValues.push_back( pt.x() ); xValues.push_back( pt.x() );
yValues.push_back( pt.y() ); yValues.push_back( pt.y() );
tvdValues.push_back( tvd ); tvdValues.push_back( tvd );
mdValues.push_back( currMd + rkbOffset ); mdValues.push_back( currMd + rkbOffset );
prevMd = currMd;
prevTvd = tvd;
isFirst = false;
} }
} }