OilfieldToolsNet/opm-simulators
4
0 Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2026-10-25 15:10:09 -05:00
Code Issues Packages Projects Releases Wiki Activity

make sure phase is active before getting phase_pos

Browse Source
This commit is contained in:
Kai Bao committed 2025-10-24 13:32:19 +02:00
1 parent 235a36da1f
commit e718ce8359
1 file changed
+13 -6
Download Patch File Download Diff File Expand all files Collapse all files
opm/simulators/wells/WellConstraints.cpp
+13 -6
View File
@@ -220,21 +220,28 @@ activeProductionConstraint(const SingleWellState<Scalar, IndexTraits>& ws,
}
if (controls.hasControl(Well::ProducerCMode::LRAT) && currentControl != Well::ProducerCMode::LRAT) {
const int oil_pos = pu.canonicalToActivePhaseIdx(IndexTraits::oilPhaseIdx);
const int water_pos = pu.canonicalToActivePhaseIdx(IndexTraits::waterPhaseIdx);
Scalar current_rate = -ws.surface_rates[oil_pos];
current_rate -= ws.surface_rates[water_pos];
Scalar current_rate = 0.;
if (pu.phaseIsActive(IndexTraits::oilPhaseIdx)) {
const int oil_pos = pu.canonicalToActivePhaseIdx(IndexTraits::oilPhaseIdx);
current_rate -= ws.surface_rates[oil_pos];
}
if (pu.phaseIsActive(IndexTraits::waterPhaseIdx)) {
const int water_pos = pu.canonicalToActivePhaseIdx(IndexTraits::waterPhaseIdx);
current_rate -= ws.surface_rates[water_pos];
}
bool skip = false;
if (controls.liquid_rate == controls.oil_rate) {
if (controls.liquid_rate == controls.oil_rate && pu.phaseIsActive(IndexTraits::waterPhaseIdx)) {
const int water_pos = pu.canonicalToActivePhaseIdx(IndexTraits::waterPhaseIdx);
const Scalar current_water_rate = ws.surface_rates[water_pos];
if (std::abs(current_water_rate) < 1e-12) {
skip = true;
deferred_logger.debug("LRAT_ORAT_WELL", "Well " + well_.name() + " The LRAT target is equal the ORAT target and the water rate is zero, skip checking LRAT");
}
}
if (!skip && controls.liquid_rate < current_rate)
if (!skip && controls.liquid_rate < current_rate) {
return Well::ProducerCMode::LRAT;
}
}
if (controls.hasControl(Well::ProducerCMode::RESV) && currentControl != Well::ProducerCMode::RESV) {