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

#6139 MSW IDC : Use length fraction for calculation of parameters

Browse Source 
Fix typo
This commit is contained in:
Magne Sjaastad 2020-06-24 15:53:51 +02:00
parent e935343fab
commit ba6d70eb46
1 changed files with 29 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

42
ApplicationCode/Commands/CompletionExportCommands/RicMswValveAccumulators.cpp
View File

@ -45,11 +45,25 @@ bool RicMswICDAccumulator::accumulateValveParameters( const RimWellPathValve* we
if ( wellPathValve->componentType() == RiaDefines::WellPathComponentType::ICV || if ( wellPathValve->componentType() == RiaDefines::WellPathComponentType::ICV ||
wellPathValve->componentType() == RiaDefines::WellPathComponentType::ICD ) wellPathValve->componentType() == RiaDefines::WellPathComponentType::ICD )
{ {
double lengthFraction = 1.0;
if ( wellPathValve->componentType() == RiaDefines::WellPathComponentType::ICD )
{
std::pair<double, double> valveSegment = wellPathValve->valveSegments()[subValve];
double valveSegmentLength = std::fabs( valveSegment.second - valveSegment.first );
if ( totalValveLengthOpenForFlow > 1.0e-8 )
{
lengthFraction = valveSegmentLength / totalValveLengthOpenForFlow;
}
}
double combinedFraction = contributionFraction * lengthFraction;
double icdOrificeRadius = wellPathValve->orificeDiameter( m_unitSystem ) / 2; double icdOrificeRadius = wellPathValve->orificeDiameter( m_unitSystem ) / 2;
double icdArea = icdOrificeRadius * icdOrificeRadius * cvf::PI_D; double icdArea = icdOrificeRadius * icdOrificeRadius * cvf::PI_D;
m_areaSum += icdArea * contributionFraction; m_areaSum += icdArea * combinedFraction;
m_coefficientCalculator.addValueAndWeight( wellPathValve->flowCoefficient(), icdArea * contributionFraction ); m_coefficientCalculator.addValueAndWeight( wellPathValve->flowCoefficient(), icdArea * combinedFraction );
return true; return true;
} }
return false; return false;
@ -99,15 +113,6 @@ bool RicMswAICDAccumulator::accumulateValveParameters( const RimWellPathValve* w
m_deviceOpen = m_deviceOpen || params->isOpen(); m_deviceOpen = m_deviceOpen || params->isOpen();
if ( params->isOpen() ) if ( params->isOpen() )
{ {
std::array<double, AICD_NUM_PARAMS> values = params->doubleValues();
for ( size_t i = 0; i < (size_t)AICD_NUM_PARAMS; ++i )
{
if ( RiaStatisticsTools::isValidNumber( values[i] ) )
{
m_meanCalculators[i].addValueAndWeight( values[i], contributionFraction );
}
}
std::pair<double, double> valveSegment = wellPathValve->valveSegments()[subValve]; std::pair<double, double> valveSegment = wellPathValve->valveSegments()[subValve];
double valveSegmentLength = std::fabs( valveSegment.second - valveSegment.first ); double valveSegmentLength = std::fabs( valveSegment.second - valveSegment.first );
double lengthFraction = 1.0; double lengthFraction = 1.0;
@ -116,6 +121,17 @@ bool RicMswAICDAccumulator::accumulateValveParameters( const RimWellPathValve* w
lengthFraction = valveSegmentLength / totalValveLengthOpenForFlow; lengthFraction = valveSegmentLength / totalValveLengthOpenForFlow;
} }
double combinedFraction = contributionFraction * lengthFraction;
std::array<double, AICD_NUM_PARAMS> values = params->doubleValues();
for ( size_t i = 0; i < (size_t)AICD_NUM_PARAMS; ++i )
{
if ( RiaStatisticsTools::isValidNumber( values[i] ) )
{
m_meanCalculators[i].addValueAndWeight( values[i], combinedFraction );
}
}
// https://github.com/OPM/ResInsight/issues/6126 // https://github.com/OPM/ResInsight/issues/6126
// //
// flowScalingFactor = 1 / (length_fraction * N_AICDs) // flowScalingFactor = 1 / (length_fraction * N_AICDs)
@ -123,11 +139,11 @@ bool RicMswAICDAccumulator::accumulateValveParameters( const RimWellPathValve* w
// length_fraction = length_COMPSEGS / Sum_lenght_COMPSEGS_for_valve // length_fraction = length_COMPSEGS / Sum_lenght_COMPSEGS_for_valve
// N_AICDs = number of AICDs in perforation interval // N_AICDs = number of AICDs in perforation interval
double divisor = wellPathValve->valveLocations().size() * lengthFraction * contributionFraction; double divisor = wellPathValve->valveLocations().size() * combinedFraction;
m_accumulatedFlowScalingFactorDivisor += divisor; m_accumulatedFlowScalingFactorDivisor += divisor;
m_accumulatedLength += lengthFraction * contributionFraction; m_accumulatedLength += combinedFraction;
} }
} }
return true; return true;