OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Output WGLIR=0 Unless ALQ=GRAT

Browse Source 
This reduces the risk of users misinterpreting the curves since the
unit of measure when plotting the vector will unconditionally be
surface volume rate for gas.

We will revisit this decision if we gain the ability to capture the
actual lift gas injection rate independent of ALQ type.

Suggested by: [at]tskille
This commit is contained in:
Bård Skaflestad 2021-06-16 13:32:44 +02:00
parent 2ad81c9364
commit b638ba0f0e
1 changed files with 2 additions and 28 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
src/opm/output/eclipse/Summary.cpp
View File

@ -576,15 +576,6 @@ alq_type(const Opm::ScheduleState& sched_state,
return sched_state.vfpprod(pc.vfp_table_number).getALQType();
}
inline bool is_glir_permissible(const Opm::VFPProdTable::ALQ_TYPE alqType)
{
using ALQ_TYPE = Opm::VFPProdTable::ALQ_TYPE;
return (alqType == ALQ_TYPE::ALQ_GRAT)
|| (alqType == ALQ_TYPE::ALQ_IGLR)
|| (alqType == ALQ_TYPE::ALQ_TGLR);
}
inline quantity glir( const fn_args& args ) {
if (args.schedule_wells.empty()) {
return { 0.0, measure::gas_surface_rate };
@ -592,14 +583,6 @@ inline quantity glir( const fn_args& args ) {
const auto& sched_state = args.schedule[args.sim_step];
auto alqType = std::optional<Opm::VFPProdTable::ALQ_TYPE>{};
if (args.schedule_wells.size() > std::vector<Opm::Well*>::size_type{1}) {
// Summing to group and field levels is only possible if the ALQ
// type is gas flow rate.
alqType = Opm::VFPProdTable::ALQ_TYPE::ALQ_GRAT;
}
double alq_rate = 0.0;
for (const auto* well : args.schedule_wells) {
if (well->isInjector()) {
@ -619,12 +602,7 @@ inline quantity glir( const fn_args& args ) {
}
const auto thisAlqType = alq_type(sched_state, production);
if (alqType.has_value() && (thisAlqType != alqType.value())) {
continue;
}
alqType = thisAlqType;
if (! is_glir_permissible(thisAlqType)) {
if (thisAlqType != Opm::VFPProdTable::ALQ_TYPE::ALQ_GRAT) {
continue;
}
@ -632,11 +610,7 @@ inline quantity glir( const fn_args& args ) {
alq_rate += eff_fac * xwPos->second.rates.get(rt::alq, production.alq_value);
}
const auto unit = !alqType.has_value() || (alqType.value() == Opm::VFPProdTable::ALQ_TYPE::ALQ_GRAT)
? measure::gas_surface_rate
: measure::gas_oil_ratio;
return { alq_rate, unit };
return { alq_rate, measure::gas_surface_rate };
}
inline quantity wwirt( const fn_args& args ) {