OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-11 07:56:08 -06:00
Code Issues Packages Projects Releases Wiki Activity

#2035 Remove the unused Pseudo Length calculator

Browse Source
This commit is contained in:
Jacob Støren 2017-11-06 10:28:13 +01:00
parent 94c6bc763d
commit 1089060539
2 changed files with 0 additions and 143 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

93
ApplicationCode/ReservoirDataModel/RigSimWellBranchPseudoLengthCalculator.cpp
View File

@ -1,93 +0,0 @@
/////////////////////////////////////////////////////////////////////////////////
//
// Copyright (C) 2017 Statoil ASA
//
// ResInsight is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// ResInsight is distributed in the hope that it will be useful, but WITHOUT ANY
// WARRANTY; without even the implied warranty of MERCHANTABILITY or
// FITNESS FOR A PARTICULAR PURPOSE.
//
// See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html>
// for more details.
//
/////////////////////////////////////////////////////////////////////////////////
#include "RigSimWellBranchPseudoLengthCalculator.h"
#include "RigSimulationWellCoordsAndMD.h"
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
RigSimWellBranchPseudoLengthCalculator::RigSimWellBranchPseudoLengthCalculator(const std::vector< std::vector <cvf::Vec3d> >& pipeBranchesCLCoords, const std::vector< std::vector <RigWellResultPoint> >& pipeBranchesCellIds): m_pipeBranchesCLCoords(pipeBranchesCLCoords),
m_pipeBranchesCellIds(pipeBranchesCellIds)
{
m_pipeBranchesCellMDs.resize(pipeBranchesCLCoords.size());
calculatePseudoLength(0, 0.0);
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RigSimWellBranchPseudoLengthCalculator::calculatePseudoLength(size_t branchIdx, double startPseudoLengthFromTop)
{
const std::vector<RigWellResultPoint>& branchCells = m_pipeBranchesCellIds[branchIdx];
const std::vector <cvf::Vec3d>& branchClPoints = m_pipeBranchesCLCoords[branchIdx];
RigSimulationWellCoordsAndMD mdCalculator(branchClPoints);
// Set the pseudo lengths for this branch
m_pipeBranchesCellMDs[branchIdx] = mdCalculator.measuredDepths();
for ( double& depth : m_pipeBranchesCellMDs[branchIdx] )
{
depth += startPseudoLengthFromTop;
}
// Find the downstream branches, and calculate for them as well, using the correct start offset
int clSegIdx = static_cast<int>(branchCells.size()) - 1;
while ( clSegIdx >= 0 )
{
RigWellResultPoint previousResultPoint = branchCells[clSegIdx];
--clSegIdx;
while ( clSegIdx >= 0 && previousResultPoint.isEqual(branchCells[clSegIdx]) ) { --clSegIdx; }
int cellUpperPointIndex = clSegIdx + 1;
double pseudoLengthFromTop_upper = mdCalculator.measuredDepths()[cellUpperPointIndex] + startPseudoLengthFromTop;
std::vector<size_t> downStreamBranchIndices = findDownStreamBranchIdxs(branchCells[cellUpperPointIndex]);
for ( size_t dsBidx : downStreamBranchIndices )
{
if ( dsBidx != branchIdx && !m_pipeBranchesCellMDs[dsBidx].size() ) // Not this branch or done before
{
calculatePseudoLength(dsBidx, pseudoLengthFromTop_upper);
}
}
}
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::vector<size_t> RigSimWellBranchPseudoLengthCalculator::findDownStreamBranchIdxs(const RigWellResultPoint& connectionPoint)
{
std::vector<size_t> downStreamBranchIdxs;
for ( size_t bIdx = 0; bIdx < m_pipeBranchesCellIds.size(); ++bIdx )
{
if ( m_pipeBranchesCellIds[bIdx][0].m_gridIndex == connectionPoint.m_gridIndex
&& m_pipeBranchesCellIds[bIdx][0].m_gridCellIndex == connectionPoint.m_gridCellIndex
&& m_pipeBranchesCellIds[bIdx][0].m_ertBranchId == connectionPoint.m_ertBranchId
&& m_pipeBranchesCellIds[bIdx][0].m_ertSegmentId == connectionPoint.m_ertSegmentId )
{
downStreamBranchIdxs.push_back(bIdx);
}
}
return downStreamBranchIdxs;
}

50
ApplicationCode/ReservoirDataModel/RigSimWellBranchPseudoLengthCalculator.h
View File

@ -1,50 +0,0 @@
/////////////////////////////////////////////////////////////////////////////////
//
// Copyright (C) 2017 Statoil ASA
//
// ResInsight is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// ResInsight is distributed in the hope that it will be useful, but WITHOUT ANY
// WARRANTY; without even the implied warranty of MERCHANTABILITY or
// FITNESS FOR A PARTICULAR PURPOSE.
//
// See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html>
// for more details.
//
/////////////////////////////////////////////////////////////////////////////////
#pragma once
#include <vector>
#include <vector>
#include "cvfBase.h"
#include "cvfVector3.h"
#include "RigSimWellData.h"
class RigSimWellBranchPseudoLengthCalculator
{
public:
RigSimWellBranchPseudoLengthCalculator (const std::vector< std::vector <cvf::Vec3d> >& pipeBranchesCLCoords,
const std::vector< std::vector <RigWellResultPoint> >& pipeBranchesCellIds);
const std::vector< std::vector <double> >& pseudoLengthsPrBranch() { return m_pipeBranchesCellMDs; }
private:
void calculatePseudoLength(size_t branchIdx, double startPseudoLengthFromTop);
std::vector<size_t> findDownStreamBranchIdxs(const RigWellResultPoint& connectionPoint);
private:
const std::vector< std::vector <RigWellResultPoint> >& m_pipeBranchesCellIds;
const std::vector< std::vector <cvf::Vec3d> >& m_pipeBranchesCLCoords;
std::vector< std::vector <double> > m_pipeBranchesCellMDs;
};