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

#1858 Fix missing connections from fishbone laterals

Browse Source
This commit is contained in:
Jacob Støren
2017-12-12 17:40:32 +01:00
parent b8b5707ebc
commit dfd25f7b63
1 changed files with 18 additions and 37 deletions
Download Patch File Download Diff File Expand all files Collapse all files

55
ApplicationCode/Commands/CompletionCommands/RicWellPathExportCompletionDataFeature.cpp
View File

@@ -876,47 +876,28 @@ void RicWellPathExportCompletionDataFeature::assignLateralIntersections(const Ri
lateralMDs.push_back(coordMD.second);
}
std::vector<WellPathCellIntersectionInfo> intersections = RigWellPathIntersectionTools::findCellIntersectionInfosAlongPath(caseToApply->eclipseCaseData(),
lateralCoords,
lateralMDs);
std::vector<WellPathCellIntersectionInfo> intersections = RigWellPathIntersectionTools::findCellIntersectionInfosAlongPath(caseToApply->eclipseCaseData(),
lateralCoords,
lateralMDs);
double previousExitMD = lateralMDs.front();
double previousExitTVD = lateralCoords.front().z();
auto intersection = intersections.cbegin();
double mdFromPreviousIntersection = 0;
double tvdChangeFromPreviousIntersection = 0;
cvf::Vec3d startPoint = lateralCoords[0];
int attachedSegmentNumber = location->icdSegmentNumber;
for (size_t i = 1; i < lateralCoords.size() && intersection != intersections.cend(); ++i)
for (const auto& cellIntInfo: intersections)
{
if (isPointBetween(startPoint, lateralCoords[i], intersection->endPoint))
{
mdFromPreviousIntersection += (intersection->endPoint - startPoint).length();
tvdChangeFromPreviousIntersection += intersection->endPoint.z() - startPoint.z();
++(*segmentNum);
++(*segmentNum);
WellSegmentLateralIntersection lateralIntersection((*segmentNum),
attachedSegmentNumber,
cellIntInfo.globCellIndex,
cellIntInfo.endMD - previousExitMD,
cellIntInfo.endPoint.z() - previousExitTVD,
cellIntInfo.intersectionLengthsInCellCS);
lateral.intersections.push_back(lateralIntersection);
WellSegmentLateralIntersection lateralIntersection((*segmentNum),
attachedSegmentNumber,
intersection->globCellIndex,
mdFromPreviousIntersection,
tvdChangeFromPreviousIntersection,
intersection->intersectionLengthsInCellCS);
lateral.intersections.push_back(lateralIntersection);
mdFromPreviousIntersection = 0;
tvdChangeFromPreviousIntersection = 0;
startPoint = intersection->endPoint;
attachedSegmentNumber = (*segmentNum);
++intersection;
}
else
{
mdFromPreviousIntersection += (lateralCoords[i] - startPoint).length();
tvdChangeFromPreviousIntersection += lateralCoords[i].z() - startPoint.z();
startPoint = lateralCoords[i];
}
attachedSegmentNumber = (*segmentNum);
previousExitMD = cellIntInfo.endMD;
previousExitTVD = cellIntInfo.endPoint.z();
}
}
}