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

update well potential for all well after gaslift

Browse Source
This commit is contained in:
Tor Harald Sandve committed 2025-12-02 12:15:33 +01:00
1 parent abab405b97
commit 0a7501ddac
2 files changed
+22 -16
Download Patch File Download Diff File Expand all files Collapse all files

No files matched your search

opm/simulators/wells/BlackoilWellModelGasLift_impl.hpp
+4
View File
@@ -136,6 +136,10 @@ maybeDoGasLiftOptimize(const Simulator& simulator,
num_wells_changed = glift_wells.size();
}
num_wells_changed = simulator.vanguard().gridView().comm().sum(num_wells_changed);
if (num_wells_changed > 0) {
updateWellPotentials(simulator, well_container, node_pressures, wellState, deferred_logger);
}
return num_wells_changed > 0;
}
opm/simulators/wells/GasLiftSingleWellGeneric.cpp
+18 -16
View File
@@ -160,24 +160,26 @@ runOptimize(const int iteration_idx)
Scalar alq = state->alq();
if (this->debug)
logSuccess_(alq, iteration_idx);
this->well_state_.well(this->well_name_).alq_state.set(alq);
auto& ws = this->well_state_.well(this->well_name_);
ws.alq_state.set(alq);
const auto& pu = this->well_state_.phaseUsageInfo();
std::vector<Scalar> well_pot(pu.numActivePhases(), 0.0);
if (pu.phaseIsActive(IndexTraits::oilPhaseIdx)) {
const int oil_pos = pu.canonicalToActivePhaseIdx(IndexTraits::oilPhaseIdx);
well_pot[oil_pos] = state->oilRate();
const bool isThp = (ws.production_cmode == Well::ProducerCMode::THP);
if (isThp) {
std::vector<Scalar> well_pot(pu.numActivePhases(), 0.0);
if (pu.phaseIsActive(IndexTraits::oilPhaseIdx)) {
const int oil_pos = pu.canonicalToActivePhaseIdx(IndexTraits::oilPhaseIdx);
well_pot[oil_pos] = state->oilRate();
}
if (pu.phaseIsActive(IndexTraits::waterPhaseIdx)) {
const int water_pos = pu.canonicalToActivePhaseIdx(IndexTraits::waterPhaseIdx);
well_pot[water_pos] = state->waterRate();
}
if (pu.phaseIsActive(IndexTraits::gasPhaseIdx)) {
const int gas_pos = pu.canonicalToActivePhaseIdx(IndexTraits::gasPhaseIdx);
well_pot[gas_pos] = state->gasRate();
}
ws.well_potentials = well_pot;
}
if (pu.phaseIsActive(IndexTraits::waterPhaseIdx)) {
const int water_pos = pu.canonicalToActivePhaseIdx(IndexTraits::waterPhaseIdx);
well_pot[water_pos] = state->waterRate();
}
if (pu.phaseIsActive(IndexTraits::gasPhaseIdx)) {
const int gas_pos = pu.canonicalToActivePhaseIdx(IndexTraits::gasPhaseIdx);
well_pot[gas_pos] = state->gasRate();
}
this->well_state_[this->well_name_].well_potentials = well_pot;
}
}
}