OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Avoid negative targets and currentWellRates

Browse Source
This commit is contained in:
Tor Harald Sandve 2020-03-26 11:46:17 +01:00 committed by Atgeirr Flø Rasmussen
parent d2dac6c19c
commit 21b98b5573
2 changed files with 11 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
opm/simulators/wells/WellGroupHelpers.hpp
View File

@ -408,8 +408,14 @@ namespace Opm {
const auto& it = wellState.wellMap().find(wellName); const auto& it = wellState.wellMap().find(wellName);
if (it != end) { // the well is found on this node if (it != end) { // the well is found on this node
int well_index = it->second[0]; int well_index = it->second[0];
const auto& wellTmp = schedule.getWell(wellName, reportStepIdx);
int sign = 1;
// production wellRates are negative. The users of currentWellRates uses the convention in
// opm-common that production and injection rates are positive.
if (!wellTmp.isInjector())
sign = -1;
for (int phase = 0; phase < np; ++phase) { for (int phase = 0; phase < np; ++phase) {
rates[phase] = wellStateNupcol.wellRates()[well_index*np + phase]; rates[phase] = sign * wellStateNupcol.wellRates()[well_index*np + phase];
} }
} }
wellState.setCurrentWellRates(wellName, rates); wellState.setCurrentWellRates(wellName, rates);
@ -1204,7 +1210,8 @@ namespace Opm {
} }
target *= localFraction(chain[ii+1]); target *= localFraction(chain[ii+1]);
} }
const double target_rate = std::max(0.0, target / efficiencyFactor); // Avoid negative target rates comming from too large local reductions.
const double target_rate = std::max(1e-12, target / efficiencyFactor);
return std::make_pair(current_rate > target_rate, target_rate / current_rate); return std::make_pair(current_rate > target_rate, target_rate / current_rate);
} }

5
opm/simulators/wells/WellInterface_impl.hpp
View File

@ -2263,10 +2263,9 @@ namespace Opm
target -= localReduction(chain[ii]); target -= localReduction(chain[ii]);
target *= localFraction(chain[ii+1]); target *= localFraction(chain[ii+1]);
} }
const double target_rate = target / efficiencyFactor; // Avoid negative target rates comming from too large local reductions.
const double target_rate = std::max(0.0, target / efficiencyFactor);
const auto current_rate = -tcalc.calcModeRateFromRates(rates); // Switch sign since 'rates' are negative for producers. const auto current_rate = -tcalc.calcModeRateFromRates(rates); // Switch sign since 'rates' are negative for producers.
control_eq = current_rate - target_rate; control_eq = current_rate - target_rate;
} }