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

#7597 MSW : Connect fractures a segment with lower measured depth

Browse Source
This commit is contained in:
Magne Sjaastad 2021-04-20 15:50:31 +02:00
parent 4f89e38334
commit dcbfd4cc27
1 changed files with 35 additions and 53 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

88
ApplicationLibCode/Commands/CompletionExportCommands/RicWellPathExportMswCompletionsImpl.cpp
View File

@ -575,68 +575,43 @@ bool RicWellPathExportMswCompletionsImpl::generateFracturesMswExportInfo(
std::vector<WellPathCellIntersectionInfo> filteredIntersections = std::vector<WellPathCellIntersectionInfo> filteredIntersections =
filterIntersections( cellIntersections, initialMD, wellPath->wellPathGeometry(), eclipseCase ); filterIntersections( cellIntersections, initialMD, wellPath->wellPathGeometry(), eclipseCase );
// Create a dummy perforation interval
RimPerforationInterval perfInterval;
bool foundSubGridIntersections = false; bool foundSubGridIntersections = false;
createWellPathSegments( branch, filteredIntersections, { &perfInterval }, wellPath, -1, eclipseCase, &foundSubGridIntersections );
// Main bore // Check if fractures are to be assigned to current main bore segment
const RigMainGrid* grid = eclipseCase->eclipseCaseData()->mainGrid(); for ( RimWellPathFracture* fracture : fractures )
for ( const auto& cellIntInfo : filteredIntersections )
{ {
if ( cellIntInfo.globCellIndex >= grid->globalCellArray().size() ) double fractureStartMD = fracture->fractureMD();
if ( fracture->fractureTemplate()->orientationType() == RimFractureTemplate::ALONG_WELL_PATH )
{ {
continue; double perforationLength = fracture->fractureTemplate()->perforationLength();
fractureStartMD -= 0.5 * perforationLength;
} }
size_t localGridIdx = 0u; auto segment = branch->findClosestSegmentWithLowerMD( fractureStartMD );
const RigGridBase* localGrid = if ( segment )
grid->gridAndGridLocalIdxFromGlobalCellIdx( cellIntInfo.globCellIndex, &localGridIdx );
QString gridName;
if ( localGrid != grid )
{ {
gridName = QString::fromStdString( localGrid->gridName() ); std::vector<RigCompletionData> completionData =
foundSubGridIntersections = true; RicExportFractureCompletionsImpl::generateCompdatValues( eclipseCase,
wellPath->completionSettings()->wellNameForExport(),
wellPath->wellPathGeometry(),
{ fracture },
nullptr,
nullptr );
assignFractureCompletionsToCellSegment( eclipseCase,
wellPath,
fracture,
completionData,
segment,
&foundSubGridIntersections );
} }
size_t i = 0u, j = 0u, k = 0u;
localGrid->ijkFromCellIndex( localGridIdx, &i, &j, &k );
auto segment = std::make_unique<RicMswSegment>( "Main stem segment",
cellIntInfo.startMD,
cellIntInfo.endMD,
cellIntInfo.startTVD(),
cellIntInfo.endTVD() );
// Check if fractures are to be assigned to current main bore segment
for ( RimWellPathFracture* fracture : fractures )
{
double fractureStartMD = fracture->fractureMD();
if ( fracture->fractureTemplate()->orientationType() == RimFractureTemplate::ALONG_WELL_PATH )
{
double perforationLength = fracture->fractureTemplate()->perforationLength();
fractureStartMD -= 0.5 * perforationLength;
}
if ( cvf::Math::valueInRange( fractureStartMD, cellIntInfo.startMD, cellIntInfo.endMD ) )
{
std::vector<RigCompletionData> completionData =
RicExportFractureCompletionsImpl::generateCompdatValues( eclipseCase,
wellPath->completionSettings()->wellNameForExport(),
wellPath->wellPathGeometry(),
{ fracture },
nullptr,
nullptr );
assignFractureCompletionsToCellSegment( eclipseCase,
wellPath,
fracture,
completionData,
segment.get(),
&foundSubGridIntersections );
}
}
branch->addSegment( std::move( segment ) );
} }
exportInfo->setHasSubGridIntersections( foundSubGridIntersections );
exportInfo->setHasSubGridIntersections( exportInfo->hasSubGridIntersections() || foundSubGridIntersections );
branch->sortSegments(); branch->sortSegments();
std::vector<RimModeledWellPath*> connectedWellPaths = wellPathsWithTieIn( wellPath ); std::vector<RimModeledWellPath*> connectedWellPaths = wellPathsWithTieIn( wellPath );
@ -662,6 +637,7 @@ bool RicWellPathExportMswCompletionsImpl::generateFracturesMswExportInfo(
return true; return true;
} }
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
/// ///
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
@ -1499,6 +1475,12 @@ void RicWellPathExportMswCompletionsImpl::assignFractureCompletionsToCellSegment
for ( const RigCompletionData& compIntersection : completionData ) for ( const RigCompletionData& compIntersection : completionData )
{ {
const RigCompletionDataGridCell& cell = compIntersection.completionDataGridCell(); const RigCompletionDataGridCell& cell = compIntersection.completionDataGridCell();
if ( !cell.isMainGridCell() )
{
*foundSubGridIntersections = true;
}
cvf::Vec3st localIJK( cell.localCellIndexI(), cell.localCellIndexJ(), cell.localCellIndexK() ); cvf::Vec3st localIJK( cell.localCellIndexI(), cell.localCellIndexJ(), cell.localCellIndexK() );
auto intersection = std::make_shared<RicMswSegmentCellIntersection>( cell.lgrName(), auto intersection = std::make_shared<RicMswSegmentCellIntersection>( cell.lgrName(),