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

Merge pull request #3452 from joakim-hove/constraint-bug

Browse Source 
BUG? Make sure WellState injection control is updated
This commit is contained in:
Joakim Hove 2021-08-03 14:38:20 +02:00 committed by GitHub
commit 6542a62cc0
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/simulators/wells/WellInterfaceFluidSystem.cpp
View File

@ -92,7 +92,7 @@ checkIndividualConstraints(WellState& well_state,
if (well.isInjector()) { if (well.isInjector()) {
const auto controls = well.injectionControls(summaryState); const auto controls = well.injectionControls(summaryState);
auto currentControl = well_state.currentInjectionControl(well_index); const auto currentControl = well_state.currentInjectionControl(well_index);
if (controls.hasControl(Well::InjectorCMode::BHP) && currentControl != Well::InjectorCMode::BHP) if (controls.hasControl(Well::InjectorCMode::BHP) && currentControl != Well::InjectorCMode::BHP)
{ {
@ -149,7 +149,7 @@ checkIndividualConstraints(WellState& well_state,
current_rate += well_state.wellReservoirRates(well_index)[ pu.phase_pos[BlackoilPhases::Vapour] ]; current_rate += well_state.wellReservoirRates(well_index)[ pu.phase_pos[BlackoilPhases::Vapour] ];
if (controls.reservoir_rate < current_rate) { if (controls.reservoir_rate < current_rate) {
currentControl = Well::InjectorCMode::RESV; well_state.currentInjectionControl(well_index, Well::InjectorCMode::RESV);
return true; return true;
} }
} }
@ -159,7 +159,7 @@ checkIndividualConstraints(WellState& well_state,
const auto& thp = getTHPConstraint(summaryState); const auto& thp = getTHPConstraint(summaryState);
double current_thp = well_state.thp(well_index); double current_thp = well_state.thp(well_index);
if (thp < current_thp) { if (thp < current_thp) {
currentControl = Well::InjectorCMode::THP; well_state.currentInjectionControl(well_index, Well::InjectorCMode::THP);
return true; return true;
} }
} }