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

#2046 Add line from well head to new branches. Handle branch lines with only one cell

Browse Source
This commit is contained in:
Jacob Støren 2017-11-02 11:09:13 +01:00 committed by Jacob Støren
parent 956b73c0c9
commit daeb7c3bdc
1 changed files with 9 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
ApplicationCode/ReservoirDataModel/RigSimulationWellCenterLineCalculator.cpp
View File

@ -718,6 +718,12 @@ private:
{
m_branchedWell.m_wellResultBranches.push_back(RigWellResultBranch());
branchIdx = static_cast<int>( m_branchedWell.m_wellResultBranches.size()) - 1;
RigWellResultPoint wellHeadAsPoint;
const RigCell& whCell = m_eclipseCaseData->cellFromWellResultCell(m_orgWellResultFrame.m_wellHead);
cvf::Vec3d whStartPos = whCell.faceCenter(cvf::StructGridInterface::NEG_K);
wellHeadAsPoint.m_bottomPosition = whStartPos;
m_branchedWell.m_wellResultBranches[branchIdx].m_branchResultPoints.push_back(wellHeadAsPoint);
}
RigWellResultBranch& currentBranch = m_branchedWell.m_wellResultBranches[branchIdx];
@ -782,6 +788,8 @@ private:
const RigCell& c1 = m_eclipseCaseData->cellFromWellResultCell(orgWellResultPoints[cIdx]);
m_cellsWithNeighbors[cIdx]; // Add an empty set for this cell, in case we have no neighbors
for ( size_t idxToCloseCell : closeCells )
{
if ( idxToCloseCell != cIdx && m_cellsWithNeighbors[cIdx].count(idxToCloseCell) == 0 )
@ -824,7 +832,7 @@ private:
{
m_unusedWellCellIndices.erase(it);
// Create a new branchline containing with the cell itself.
// Create a new branchline containing the cell itself.
m_branchLines.push_back(std::make_pair(true, std::deque<size_t>()));
auto currentBranchLineIt = std::prev(m_branchLines.end());
auto& branchList = currentBranchLineIt->second;