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

adding a fix for the zero rate control for StandardWell

Browse Source 
to avoid under-determined problem.
This commit is contained in:
Kai Bao 2018-05-16 15:35:58 +02:00
parent efbac8a628
commit ce849fe18e
1 changed files with 15 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
opm/autodiff/StandardWell_impl.hpp
View File

@ -668,14 +668,22 @@ namespace Opm
assert(number_phases_under_control == 1); // only handles single phase injection now
control_eq = getGTotal() - target_rate;
} else if (well_type_ == PRODUCER) {
EvalWell rate_for_control(0.);
const EvalWell& g_total = getGTotal();
for (int phase = 0; phase < number_of_phases_; ++phase) {
if (distr[phase] > 0.) {
rate_for_control += g_total * wellVolumeFractionScaled(flowPhaseToEbosCompIdx(phase));
}
if (target_rate != 0.) {
EvalWell rate_for_control(0.);
const EvalWell& g_total = getGTotal();
for (int phase = 0; phase < number_of_phases_; ++phase) {
if (distr[phase] > 0.) {
rate_for_control += g_total * wellVolumeFractionScaled(flowPhaseToEbosCompIdx(phase));
}
}
control_eq = rate_for_control - target_rate;
} else {
// there is some special treatment for the zero rate control well
// 1. if the well can produce the specified phase, it means the well should not produce any fluid
// 2. if the well can not produce the specified phase, it cause a under-determined problem, we
// basically assume the well not producing any fluid as a solution
control_eq = getGTotal() - target_rate;
}
control_eq = rate_for_control - target_rate;
}
break;
}