OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-26 03:00:17 -06:00
Code Issues Packages Projects Releases Wiki Activity

Bugfix. Fix setting initial rates in updateWellControls

Browse Source 
- The initial rates are only set to target values for single phase
producers (orat, wrat, grat).
- For injectors compi is used to determine the initial target rates.
This commit is contained in:
Tor Harald Sandve 2016-03-18 13:16:23 +01:00
parent c1b1b9296f
commit 94aa360fe6
1 changed files with 25 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
opm/autodiff/BlackoilModelBase_impl.hpp
View File

@ -1610,11 +1610,33 @@ namespace detail {
break;
case SURFACE_RATE:
for (int phase = 0; phase < np; ++phase) {
if (distr[phase] > 0.0) {
xw.wellRates()[np*w + phase] = target * distr[phase];
// assign target value as initial guess for injectors and
// single phase producers (orat, grat, wrat)
const WellType& well_type = wells().type[w];
if (well_type == INJECTOR) {
for (int phase = 0; phase < np; ++phase) {
const double& compi = wells().comp_frac[np * w + phase];
if (compi > 0.0) {
xw.wellRates()[np*w + phase] = target * compi;
}
}
} else if (well_type == PRODUCER) {
// for single phase producers sum of distr should be 1.0
double sumdistr = distr[0];
for (int phase = 1; phase < np; ++phase) {
sumdistr += distr[phase];
}
std::cout << sumdistr << std::endl;
for (int phase = 0; phase < np; ++phase) {
if (distr[phase] > 0.0 && sumdistr == 1.0 ) {
xw.wellRates()[np*w + phase] = target * distr[phase];
}
}
} else {
OPM_THROW(std::logic_error, "Expected PRODUCER or INJECTOR type of well");
}
break;
}