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

Some refactor before #1203

Browse Source
This commit is contained in:
Jacob Støren 2017-03-09 17:04:43 +01:00
parent 5d0df1c86f
commit 8b3029c3a3
2 changed files with 18 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
ApplicationCode/ReservoirDataModel/RigAccWellFlowCalculator.cpp
View File

@ -59,7 +59,8 @@ RigAccWellFlowCalculator::RigAccWellFlowCalculator(const std::vector< std::vecto
m_pipeBranchesCellIds(pipeBranchesCellIds), m_pipeBranchesCellIds(pipeBranchesCellIds),
m_tracerCellFractionValues(&tracerCellFractionValues), m_tracerCellFractionValues(&tracerCellFractionValues),
m_cellIndexCalculator(cellIndexCalculator), m_cellIndexCalculator(cellIndexCalculator),
m_smallContributionsThreshold(smallContribThreshold) m_smallContributionsThreshold(smallContribThreshold),
m_isProducer(isProducer)
{ {
m_connectionFlowPrBranch.resize(m_pipeBranchesCellIds.size()); m_connectionFlowPrBranch.resize(m_pipeBranchesCellIds.size());
m_pseudoLengthFlowPrBranch.resize(m_pipeBranchesCellIds.size()); m_pseudoLengthFlowPrBranch.resize(m_pipeBranchesCellIds.size());
@ -98,7 +99,8 @@ RigAccWellFlowCalculator::RigAccWellFlowCalculator(const std::vector< std::vecto
m_pipeBranchesCellIds(pipeBranchesCellIds), m_pipeBranchesCellIds(pipeBranchesCellIds),
m_tracerCellFractionValues(nullptr), m_tracerCellFractionValues(nullptr),
m_cellIndexCalculator(RigEclCellIndexCalculator(nullptr, nullptr)), m_cellIndexCalculator(RigEclCellIndexCalculator(nullptr, nullptr)),
m_smallContributionsThreshold(smallContribThreshold) m_smallContributionsThreshold(smallContribThreshold),
m_isProducer(true)
{ {
m_connectionFlowPrBranch.resize(m_pipeBranchesCellIds.size()); m_connectionFlowPrBranch.resize(m_pipeBranchesCellIds.size());
m_pseudoLengthFlowPrBranch.resize(m_pipeBranchesCellIds.size()); m_pseudoLengthFlowPrBranch.resize(m_pipeBranchesCellIds.size());
@ -331,7 +333,9 @@ void RigAccWellFlowCalculator::calculateAccumulatedFlowPrConnection(size_t branc
// Accumulate the connection-cell's fraction flows // Accumulate the connection-cell's fraction flows
std::vector<double> flowPrTracer = calculateFlowPrTracer(branchCells, clSegIdx); const RigWellResultPoint& wellCell = branchCells[clSegIdx];
std::vector<double> flowPrTracer = calculateFlowPrTracer(wellCell);
for (size_t tIdx = 0; tIdx < flowPrTracer.size(); ++tIdx) for (size_t tIdx = 0; tIdx < flowPrTracer.size(); ++tIdx)
{ {
@ -400,7 +404,7 @@ void RigAccWellFlowCalculator::calculateFlowPrPseudoLength(size_t branchIdx, dou
currentSegmentIndex = cellUpperPointIndex; currentSegmentIndex = cellUpperPointIndex;
} }
std::vector<double> flowPrTracer = calculateFlowPrTracer(branchCells, currentSegmentIndex); std::vector<double> flowPrTracer = calculateFlowPrTracer(branchCells[currentSegmentIndex]);
double pseudoLengthFromTop_lower = mdCalculator.measuredDepths()[cellBottomPointIndex] + startPseudoLengthFromTop; double pseudoLengthFromTop_lower = mdCalculator.measuredDepths()[cellBottomPointIndex] + startPseudoLengthFromTop;
double tvd_lower = -mdCalculator.wellPathPoints()[cellBottomPointIndex][2]; double tvd_lower = -mdCalculator.wellPathPoints()[cellBottomPointIndex][2];
@ -496,16 +500,16 @@ void RigAccWellFlowCalculator::storeFlowOnDepthWTvd(BranchFlow *branchFlow, doub
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
/// ///
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
std::vector<double> RigAccWellFlowCalculator::calculateFlowPrTracer(const std::vector<RigWellResultPoint>& branchCells, std::vector<double> RigAccWellFlowCalculator::calculateFlowPrTracer(const RigWellResultPoint& wellCell ) const
int clSegIdx) const
{ {
std::vector<double> flowPrTracer(m_tracerNames.size(), 0.0); std::vector<double> flowPrTracer(m_tracerNames.size(), 0.0);
if ( m_tracerCellFractionValues ) if ( m_tracerCellFractionValues )
{ {
if ( branchCells[clSegIdx].isCell() && branchCells[clSegIdx].m_isOpen ) if ( wellCell.isCell() && wellCell.m_isOpen )
{ {
size_t resCellIndex = m_cellIndexCalculator.resultCellIndex(branchCells[clSegIdx].m_gridIndex, size_t resCellIndex = m_cellIndexCalculator.resultCellIndex(wellCell.m_gridIndex,
branchCells[clSegIdx].m_gridCellIndex); wellCell.m_gridCellIndex);
size_t tracerIdx = 0; size_t tracerIdx = 0;
double totalTracerFractionInCell = 0.0; double totalTracerFractionInCell = 0.0;
for ( const auto & tracerFractionIt: (*m_tracerCellFractionValues) ) for ( const auto & tracerFractionIt: (*m_tracerCellFractionValues) )
@ -513,7 +517,7 @@ std::vector<double> RigAccWellFlowCalculator::calculateFlowPrTracer(const std::v
double cellTracerFraction = (*tracerFractionIt.second)[resCellIndex]; double cellTracerFraction = (*tracerFractionIt.second)[resCellIndex];
if ( cellTracerFraction != HUGE_VAL && cellTracerFraction == cellTracerFraction ) if ( cellTracerFraction != HUGE_VAL && cellTracerFraction == cellTracerFraction )
{ {
double tracerFlow = cellTracerFraction * branchCells[clSegIdx].flowRate(); double tracerFlow = cellTracerFraction * wellCell.flowRate();
flowPrTracer[tracerIdx] = tracerFlow; flowPrTracer[tracerIdx] = tracerFlow;
totalTracerFractionInCell += cellTracerFraction; totalTracerFractionInCell += cellTracerFraction;
@ -522,13 +526,13 @@ std::vector<double> RigAccWellFlowCalculator::calculateFlowPrTracer(const std::v
} }
double reservoirFraction = 1.0 - totalTracerFractionInCell; double reservoirFraction = 1.0 - totalTracerFractionInCell;
double reservoirTracerFlow = reservoirFraction * branchCells[clSegIdx].flowRate(); double reservoirTracerFlow = reservoirFraction * wellCell.flowRate();
flowPrTracer[tracerIdx] = reservoirTracerFlow; flowPrTracer[tracerIdx] = reservoirTracerFlow;
} }
} }
else else
{ {
flowPrTracer[0] = branchCells[clSegIdx].flowRate(); flowPrTracer[0] = wellCell.flowRate();
} }
return flowPrTracer; return flowPrTracer;

3
ApplicationCode/ReservoirDataModel/RigAccWellFlowCalculator.h
View File

@ -94,7 +94,7 @@ private:
std::vector<double> calculateFlowPrTracer(const std::vector<RigWellResultPoint> &branchCells, int clSegIdx) const; std::vector<double> calculateFlowPrTracer(const RigWellResultPoint& wellCell) const;
void sortTracers(); void sortTracers();
void groupSmallContributions(); void groupSmallContributions();
@ -116,6 +116,7 @@ private:
RigEclCellIndexCalculator m_cellIndexCalculator; RigEclCellIndexCalculator m_cellIndexCalculator;
std::vector<QString> m_tracerNames; std::vector<QString> m_tracerNames;
double m_smallContributionsThreshold; double m_smallContributionsThreshold;
bool m_isProducer;
struct BranchFlow struct BranchFlow
{ {