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

MswRollUp: Removed obsolete code

Browse Source 
This is an intermediate commit and does not compile
p4#: 22220
This commit is contained in:
Jacob Støren 2013-08-26 14:27:39 +02:00
parent 9f2bc4f1d0
commit b807a0e0ee
1 changed files with 14 additions and 74 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

88
ApplicationCode/ModelVisualization/RivWellPipesPartMgr.cpp
View File

@ -47,7 +47,6 @@
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
/// ///
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
@ -219,28 +218,8 @@ void RivWellPipesPartMgr::calculateWellPipeCenterline( std::vector< std::vector
cvf::Vec3d whIntermediate = whStartPos; cvf::Vec3d whIntermediate = whStartPos;
whIntermediate.z() = (whStartPos.z() + whCell.center().z()) / 2.0; whIntermediate.z() = (whStartPos.z() + whCell.center().z()) / 2.0;
const RigWellResultPoint* prevWellResPoint = NULL; const RigWellResultPoint* prevWellResPoint = NULL;
#if 0
// Use well head if branch head is not specified
if (!wellResults->isMultiSegmentWell())
{
// Create a new branch from wellhead
pipeBranchesCLCoords.push_back(std::vector<cvf::Vec3d>());
pipeBranchesCellIds.push_back(std::vector <RigWellResultPoint>());
prevResCell = whResCell;
pipeBranchesCLCoords.back().push_back(whStartPos);
pipeBranchesCellIds.back().push_back(*prevResCell);
pipeBranchesCLCoords.back().push_back(whIntermediate);
pipeBranchesCellIds.back().push_back(*prevResCell);
}
#endif
CVF_ASSERT(wellResults->isMultiSegmentWell() || resBranches.size() <= 1); CVF_ASSERT(wellResults->isMultiSegmentWell() || resBranches.size() <= 1);
// The centerline is calculated by adding a point when the pipe enters a cell, // The centerline is calculated by adding a point when the pipe enters a cell,
@ -254,6 +233,7 @@ void RivWellPipesPartMgr::calculateWellPipeCenterline( std::vector< std::vector
{ {
// Skip empty branches. Do not know why they exist, but they make problems. // Skip empty branches. Do not know why they exist, but they make problems.
bool hasValidData = false; bool hasValidData = false;
for (size_t cIdx = 0; cIdx < resBranches[brIdx].m_branchResultPoints.size(); ++cIdx) for (size_t cIdx = 0; cIdx < resBranches[brIdx].m_branchResultPoints.size(); ++cIdx)
{ {
@ -268,66 +248,26 @@ void RivWellPipesPartMgr::calculateWellPipeCenterline( std::vector< std::vector
prevWellResPoint = NULL; prevWellResPoint = NULL;
// Find the start the MSW well-branch centerline. Normal wells are started "once" at wellhead in the code above // Find the start the MSW well-branch centerline. Normal wells are started "once" at wellhead in the code above
// if (wellResults->isMultiSegmentWell()) pipeBranchesCLCoords.push_back(std::vector<cvf::Vec3d>());
// { pipeBranchesCellIds.push_back(std::vector <RigWellResultPoint>());
pipeBranchesCLCoords.push_back(std::vector<cvf::Vec3d>());
pipeBranchesCellIds.push_back(std::vector <RigWellResultPoint>());
if (brIdx == 0) if (brIdx == 0)
{ {
// The first branch contains segment number 1, and this is the only segment connected to well head // The first branch contains segment number 1, and this is the only segment connected to well head
// See Eclipse documentation for the keyword WELSEGS // See Eclipse documentation for the keyword WELSEGS
prevWellResPoint = whResCell; prevWellResPoint = whResCell;
pipeBranchesCLCoords.back().push_back(whStartPos); pipeBranchesCLCoords.back().push_back(whStartPos);
pipeBranchesCellIds.back().push_back(*prevWellResPoint); pipeBranchesCellIds.back().push_back(*prevWellResPoint);
pipeBranchesCLCoords.back().push_back(whIntermediate); pipeBranchesCLCoords.back().push_back(whIntermediate);
pipeBranchesCellIds.back().push_back(*prevWellResPoint); pipeBranchesCellIds.back().push_back(*prevWellResPoint);
} }
#if 0 // Branch is supposed to contain its start point except the
else
{
const RigWellResultPoint* leafBranchHead = staticWellFrame.findResultCellFromOutletSpecification(resBranches[brIdx].m_outletBranchIndex_OBSOLETE, resBranches[brIdx].m_outletBranchHeadCellIndex_OBSOLETE);
if (leafBranchHead && leafBranchHead->hasGridConnections())
{
// Create a new branch and use branch head as previous result cell
prevResCell = leafBranchHead;
const RigCell& cell = rigReservoir->cellFromWellResultCell(*leafBranchHead);
cvf::Vec3d branchHeadStartPos = cell.faceCenter(cvf::StructGridInterface::NEG_K);
pipeBranchesCLCoords.back().push_back(branchHeadStartPos);
pipeBranchesCellIds.back().push_back(*prevResCell);
}
else if (leafBranchHead)
{
cvf::Vec3d interpolatedCoord = leafBranchHead->m_bottomPosition;
CVF_ASSERT(interpolatedCoord != cvf::Vec3d::UNDEFINED);
if (interpolatedCoord != cvf::Vec3d::UNDEFINED)
{
pipeBranchesCLCoords.back().push_back(interpolatedCoord);
pipeBranchesCellIds.back().push_back(RigWellResultPoint());
}
}
else
{
// No branch head found: Possibly main branch
CVF_ASSERT(false);
}
}
#endif
// }
// Loop over all the resultCells in the branch // Loop over all the resultCells in the branch
const std::vector<RigWellResultPoint>& resBranchCells = resBranches[brIdx].m_branchResultPoints; const std::vector<RigWellResultPoint>& resBranchCells = resBranches[brIdx].m_branchResultPoints;
for (int cIdx = 0; cIdx < static_cast<int>(resBranchCells.size()); cIdx++) // Need int because cIdx can temporarily end on -1 for (int cIdx = 0; cIdx < static_cast<int>(resBranchCells.size()); cIdx++) // Need int because cIdx can temporarily end on -1