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

#7678 Fishbones MSW : Always create an intersection for an ICD segment

Browse Source
This commit is contained in:
Magne Sjaastad 2021-05-13 20:19:49 +02:00
parent 460a387021
commit 9a6073eb7b
1 changed files with 17 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -530,7 +530,8 @@ void RicWellPathExportMswCompletionsImpl::generateFishbonesMswExportInfo(
}
// Find cell intersections closest to each sub location
std::map<size_t, std::vector<size_t>> subAndCellIntersectionIndices;
std::map<size_t, std::vector<size_t>> closestSubForCellIntersections;
std::map<size_t, size_t> cellIntersectionContainingSubIndex;
{
auto fishboneSectionStart = subs->startMD();
auto fishboneSectionEnd = subs->endMD();
@ -549,6 +550,11 @@ void RicWellPathExportMswCompletionsImpl::generateFishbonesMswExportInfo(
{
double subMD = subs->measuredDepth( sub.first );
if ( ( cellIntersection.startMD <= subMD ) && ( subMD <= cellIntersection.endMD ) )
{
cellIntersectionContainingSubIndex[sub.first] = intersectionIndex;
}
auto distanceCandicate = std::abs( subMD - intersectionMidpoint );
if ( distanceCandicate < closestDistance )
{
@ -557,7 +563,7 @@ void RicWellPathExportMswCompletionsImpl::generateFishbonesMswExportInfo(
}
}
subAndCellIntersectionIndices[closestSubIndex].push_back( intersectionIndex );
closestSubForCellIntersections[closestSubIndex].push_back( intersectionIndex );
}
}
}
@ -597,7 +603,15 @@ void RicWellPathExportMswCompletionsImpl::generateFishbonesMswExportInfo(
{
const RigMainGrid* mainGrid = eclipseCase->mainGrid();
for ( auto intersectionIndex : subAndCellIntersectionIndices[sub.first] )
std::set<size_t> indices;
for ( auto intersectionIndex : closestSubForCellIntersections[sub.first] )
{
indices.insert( intersectionIndex );
}
indices.insert( cellIntersectionContainingSubIndex[sub.first] );
for ( auto intersectionIndex : indices )
{
auto intersection = filteredIntersections[intersectionIndex];
if ( intersection.globCellIndex >= mainGrid->globalCellArray().size() ) continue;