OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-05 22:13:00 -06:00
Code Issues Packages Projects Releases Wiki Activity

scale the gas equation, not the control

Browse Source
This commit is contained in:
Tor Harald Sandve 2019-10-24 10:54:50 +02:00
parent e1f09d89cd
commit 9a3ce14885
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/simulators/wells/MultisegmentWell_impl.hpp
View File

@ -1650,7 +1650,7 @@ namespace Opm
switch(current) {
case Well2::InjectorCMode::RATE:
{
control_eq = getSegmentGTotal(0) * efficiencyFactor - controls.surface_rate*scaling;
control_eq = getSegmentGTotal(0) * efficiencyFactor / scaling - controls.surface_rate;
break;
}
@ -1682,7 +1682,7 @@ namespace Opm
}
control_eq = coeff*getSegmentGTotal(0)*efficiencyFactor - controls.reservoir_rate*scaling;
control_eq = coeff*getSegmentGTotal(0)*efficiencyFactor / scaling - controls.reservoir_rate;
break;
}
@ -1881,7 +1881,7 @@ namespace Opm
case Group2::InjectionCMode::RATE:
{
control_eq = getSegmentGTotal(0) - fraction * scaling * (groupcontrols.surface_max_rate / efficiencyFactor - groupTargetReduction);
control_eq = getSegmentGTotal(0) / scaling - fraction * (groupcontrols.surface_max_rate / efficiencyFactor - groupTargetReduction);
break;
}
case Group2::InjectionCMode::RESV:
@ -1890,7 +1890,7 @@ namespace Opm
Base::rateConverter_.calcCoeff(/*fipreg*/ 0, Base::pvtRegionIdx_, convert_coeff);
double coeff = convert_coeff[phasePos];
double target = std::max(0.0, (groupcontrols.resv_max_rate/coeff/efficiencyFactor - groupTargetReduction));
control_eq = getSegmentGTotal(0) - fraction * scaling * target;
control_eq = getSegmentGTotal(0) / scaling - fraction * target;
break;
}
case Group2::InjectionCMode::REIN:
@ -1899,7 +1899,7 @@ namespace Opm
productionRate += wellGroupHelpers::sumWellRates(group, schedule, well_state, current_step_, phasePos, /*isInjector*/false);
productionRate /= efficiencyFactor;
double target = std::max(0.0, (groupcontrols.target_reinj_fraction*productionRate - groupTargetReduction));
control_eq = getSegmentGTotal(0) - fraction * scaling * target;
control_eq = getSegmentGTotal(0) / scaling - fraction * target;
break;
}
case Group2::InjectionCMode::VREP:
@ -1913,7 +1913,7 @@ namespace Opm
voidageRate += wellGroupHelpers::sumWellResRates(group, schedule, well_state, current_step_, pu.phase_pos[BlackoilPhases::Vapour], /*injector*/false);
voidageRate /= efficiencyFactor;
double target = std::max(0.0, ( groupcontrols.target_void_fraction*voidageRate/coeff - groupTargetReduction));
control_eq = getSegmentGTotal(0) - fraction * scaling * target ;
control_eq = getSegmentGTotal(0) / scaling - fraction * target ;
break;
}
case Group2::InjectionCMode::FLD: