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

GasLiftStage2: template Scalar type

Browse Source
This commit is contained in:
Arne Morten Kvarving 2024-02-20 08:24:50 +01:00
parent 09065dd130
commit ababeb79c4
2 changed files with 397 additions and 316 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

445
opm/simulators/wells/GasLiftStage2.cpp
View File

@ -39,21 +39,20 @@
namespace Opm { namespace Opm {
GasLiftStage2::GasLiftStage2( template<class Scalar>
const int report_step_idx, GasLiftStage2<Scalar>::GasLiftStage2(const int report_step_idx,
const Parallel::Communication& comm, const Parallel::Communication& comm,
const Schedule& schedule, const Schedule& schedule,
const SummaryState& summary_state, const SummaryState& summary_state,
DeferredLogger &deferred_logger, DeferredLogger &deferred_logger,
WellState<double>& well_state, WellState<Scalar>& well_state,
const GroupState<double>& group_state, const GroupState<Scalar>& group_state,
GLiftProdWells &prod_wells, GLiftProdWells &prod_wells,
GLiftOptWells &glift_wells, GLiftOptWells &glift_wells,
GasLiftGroupInfo<double>& group_info, GasLiftGroupInfo<Scalar>& group_info,
GLiftWellStateMap &state_map, GLiftWellStateMap &state_map,
bool glift_debug bool glift_debug)
) : : GasLiftCommon<Scalar>(well_state, group_state, deferred_logger, comm, glift_debug)
GasLiftCommon(well_state, group_state, deferred_logger, comm, glift_debug)
, prod_wells_{prod_wells} , prod_wells_{prod_wells}
, stage1_wells_{glift_wells} , stage1_wells_{glift_wells}
, group_info_{group_info} , group_info_{group_info}
@ -80,14 +79,11 @@ GasLiftStage2::GasLiftStage2(
// currently has the largest weighted incremental gradient. The // currently has the largest weighted incremental gradient. The
// procedure takes account of any limits on the group production rate // procedure takes account of any limits on the group production rate
// or lift gas supply applied to any level of group, including the FIELD level group. // or lift gas supply applied to any level of group, including the FIELD level group.
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::runOptimize()
runOptimize()
{ {
const auto& group = this->schedule_.getGroup("FIELD", this->report_step_idx_); const auto& group = this->schedule_.getGroup("FIELD", this->report_step_idx_);
optimizeGroupsRecursive_(group); optimizeGroupsRecursive_(group);
} }
@ -100,17 +96,19 @@ runOptimize()
// saved in "grad_map") // saved in "grad_map")
// INPUT: grad_map : map of incremental (if "add" is true) or decremental // INPUT: grad_map : map of incremental (if "add" is true) or decremental
// (if "add" is false) GradInfo structs for each well name. // (if "add" is false) GradInfo structs for each well name.
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
addOrRemoveALQincrement_(GradMap &grad_map, const std::string& well_name, bool add) addOrRemoveALQincrement_(GradMap &grad_map,
const std::string& well_name,
bool add)
{ {
// only applies to wells in the well_state_map (i.e. wells on this rank) // only applies to wells in the well_state_map (i.e. wells on this rank)
auto it = this->well_state_map_.find(well_name); auto it = this->well_state_map_.find(well_name);
if (it == this->well_state_map_.end()) if (it == this->well_state_map_.end())
return; return;
GasLiftWellState<double>& state = *(it->second.get()); GasLiftWellState<Scalar>& state = *(it->second.get());
const GradInfo &gi = grad_map.at(well_name); const GradInfo& gi = grad_map.at(well_name);
if (this->debug) { if (this->debug) {
auto new_alq = gi.alq; auto new_alq = gi.alq;
auto old_alq = state.alq(); auto old_alq = state.alq();
@ -120,32 +118,35 @@ addOrRemoveALQincrement_(GradMap &grad_map, const std::string& well_name, bool a
this->displayDebugMessage_(msg); this->displayDebugMessage_(msg);
} }
state.update(gi.new_oil_rate, gi.oil_is_limited, state.update(gi.new_oil_rate, gi.oil_is_limited,
gi.new_gas_rate, gi.gas_is_limited, gi.new_gas_rate, gi.gas_is_limited,
gi.alq, gi.alq_is_limited, gi.new_water_rate, gi.water_is_limited, add); gi.alq, gi.alq_is_limited,
gi.new_water_rate, gi.water_is_limited, add);
this->well_state_.setALQ(well_name, gi.alq); this->well_state_.setALQ(well_name, gi.alq);
const auto& pu = this->well_state_.phaseUsage(); const auto& pu = this->well_state_.phaseUsage();
std::vector<double> well_pot(pu.num_phases, 0.0); std::vector<Scalar> well_pot(pu.num_phases, 0.0);
if(pu.phase_used[BlackoilPhases::PhaseIndex::Liquid]) if (pu.phase_used[BlackoilPhases::PhaseIndex::Liquid])
well_pot[pu.phase_pos[BlackoilPhases::PhaseIndex::Liquid]] = gi.new_oil_rate; well_pot[pu.phase_pos[BlackoilPhases::PhaseIndex::Liquid]] = gi.new_oil_rate;
if(pu.phase_used[BlackoilPhases::PhaseIndex::Aqua]) if (pu.phase_used[BlackoilPhases::PhaseIndex::Aqua])
well_pot[pu.phase_pos[BlackoilPhases::PhaseIndex::Aqua]] = gi.new_water_rate; well_pot[pu.phase_pos[BlackoilPhases::PhaseIndex::Aqua]] = gi.new_water_rate;
if(pu.phase_used[BlackoilPhases::PhaseIndex::Vapour]) if (pu.phase_used[BlackoilPhases::PhaseIndex::Vapour])
well_pot[pu.phase_pos[BlackoilPhases::PhaseIndex::Vapour]] = gi.new_gas_rate; well_pot[pu.phase_pos[BlackoilPhases::PhaseIndex::Vapour]] = gi.new_gas_rate;
this->well_state_[well_name].well_potentials = well_pot; this->well_state_[well_name].well_potentials = well_pot;
} }
std::optional<GasLiftStage2::GradInfo> template<class Scalar>
GasLiftStage2:: std::optional<typename GasLiftStage2<Scalar>::GradInfo>
calcIncOrDecGrad_( GasLiftStage2<Scalar>::
const std::string well_name, const GasLiftSingleWell &gs_well, const std::string& gr_name_dont_limit, bool increase) calcIncOrDecGrad_(const std::string well_name,
const GasLiftSingleWell& gs_well,
const std::string& gr_name_dont_limit,
bool increase)
{ {
// only applies to wells in the well_state_map (i.e. wells on this rank) // only applies to wells in the well_state_map (i.e. wells on this rank)
if(this->well_state_map_.count(well_name) == 0) if (this->well_state_map_.count(well_name) == 0)
return std::nullopt; return std::nullopt;
GasLiftWellState<double>& state = *(this->well_state_map_.at(well_name).get()); GasLiftWellState<Scalar>& state = *(this->well_state_map_.at(well_name).get());
if (checkRateAlreadyLimited_(well_name, state, increase)) { if (checkRateAlreadyLimited_(well_name, state, increase)) {
return std::nullopt; return std::nullopt;
} }
@ -179,16 +180,18 @@ calcIncOrDecGrad_(
} }
} }
bool template<class Scalar>
GasLiftStage2:: bool GasLiftStage2<Scalar>::
checkRateAlreadyLimited_(const std::string& well_name, checkRateAlreadyLimited_(const std::string& well_name,
GasLiftWellState<double>& state, GasLiftWellState<Scalar>& state,
bool increase) bool increase)
{ {
auto current_increase = state.increase(); auto current_increase = state.increase();
bool do_check = false; bool do_check = false;
if (current_increase) { if (current_increase) {
if (*current_increase == increase) do_check = true; if (*current_increase == increase) {
do_check = true;
}
} }
else { else {
// If current_increase is not defined, it means that stage1 // If current_increase is not defined, it means that stage1
@ -198,10 +201,13 @@ checkRateAlreadyLimited_(const std::string& well_name,
// cause the new rates not to be limited. However, if // cause the new rates not to be limited. However, if
// "increase" is false, subtracting ALQ can make the new rates // "increase" is false, subtracting ALQ can make the new rates
// not limited. // not limited.
if (increase) do_check = true; if (increase) {
do_check = true;
}
} }
if (do_check) { if (do_check) {
if (state.gasIsLimited() || state.oilIsLimited() || state.alqIsLimited() || state.waterIsLimited()) { if (state.gasIsLimited() || state.oilIsLimited() ||
state.alqIsLimited() || state.waterIsLimited()) {
const std::string msg = fmt::format( const std::string msg = fmt::format(
"Well {} : alq = {} : skipping {} gradient since {} was limited in previous step", "Well {} : alq = {} : skipping {} gradient since {} was limited in previous step",
well_name, well_name,
@ -216,9 +222,10 @@ checkRateAlreadyLimited_(const std::string& well_name,
return false; return false;
} }
GasLiftStage2::GradInfo template<class Scalar>
GasLiftStage2:: typename GasLiftStage2<Scalar>::GradInfo
deleteGrad_(const std::string &name, bool increase) GasLiftStage2<Scalar>::
deleteGrad_(const std::string& name, bool increase)
{ {
GradMap &map = increase ? this->inc_grads_ : this->dec_grads_; GradMap &map = increase ? this->inc_grads_ : this->dec_grads_;
auto value = map.at(name); auto value = map.at(name);
@ -226,56 +233,58 @@ deleteGrad_(const std::string &name, bool increase)
return value; return value;
} }
GasLiftStage2::GradInfo template<class Scalar>
GasLiftStage2:: typename GasLiftStage2<Scalar>::GradInfo
deleteDecGradItem_(const std::string &name) GasLiftStage2<Scalar>::
deleteDecGradItem_(const std::string& name)
{ {
return deleteGrad_(name, /*increase=*/false); return deleteGrad_(name, /*increase=*/false);
} }
GasLiftStage2::GradInfo template<class Scalar>
GasLiftStage2:: typename GasLiftStage2<Scalar>::GradInfo
deleteIncGradItem_(const std::string &name) GasLiftStage2<Scalar>::
deleteIncGradItem_(const std::string& name)
{ {
return deleteGrad_(name, /*increase=*/true); return deleteGrad_(name, /*increase=*/true);
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
displayWarning_(const std::string &msg, const std::string &group_name) displayWarning_(const std::string& msg, const std::string& group_name)
{ {
const std::string message = fmt::format("GROUP: {} : {}", group_name, msg); const std::string message = fmt::format("GROUP: {} : {}", group_name, msg);
displayWarning_(message); displayWarning_(message);
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
displayWarning_(const std::string &msg) displayWarning_(const std::string& msg)
{ {
logMessage_(/*prefix=*/"GLIFT2", msg, MessageType::WARNING); this->logMessage_(/*prefix=*/"GLIFT2", msg, MessageType::WARNING);
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
displayDebugMessage_(const std::string &msg) const displayDebugMessage_(const std::string& msg) const
{ {
if (this->debug) { if (this->debug) {
logMessage_(/*prefix=*/"GLIFT2", msg); this->logMessage_(/*prefix=*/"GLIFT2", msg);
} }
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
displayDebugMessage2B_(const std::string &msg) displayDebugMessage2B_(const std::string& msg)
{ {
if (this->debug) { if (this->debug) {
logMessage_(/*prefix=*/"GLIFT2B", msg); this->logMessage_(/*prefix=*/"GLIFT2B", msg);
} }
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
displayDebugMessage_(const std::string &msg, const std::string &group_name) displayDebugMessage_(const std::string& msg, const std::string& group_name)
{ {
if (this->debug) { if (this->debug) {
const std::string message = fmt::format( const std::string message = fmt::format(
@ -284,9 +293,10 @@ displayDebugMessage_(const std::string &msg, const std::string &group_name)
} }
} }
std::tuple<double, double, double, double> template<class Scalar>
GasLiftStage2:: std::tuple<Scalar, Scalar, Scalar, Scalar>
getCurrentGroupRates_(const Group &group) GasLiftStage2<Scalar>::
getCurrentGroupRates_(const Group& group)
{ {
return {this->group_info_.oilRate(group.name()), return {this->group_info_.oilRate(group.name()),
this->group_info_.gasRate(group.name()), this->group_info_.gasRate(group.name()),
@ -294,21 +304,23 @@ getCurrentGroupRates_(const Group &group)
this->group_info_.alqRate(group.name())}; this->group_info_.alqRate(group.name())};
} }
std::optional<double> template<class Scalar>
GasLiftStage2::getGroupMaxALQ_(const Group &group) std::optional<Scalar>
GasLiftStage2<Scalar>::getGroupMaxALQ_(const Group& group)
{ {
if (this->glo_.has_group(group.name())) { if (this->glo_.has_group(group.name())) {
const auto &gl_group = this->glo_.group(group.name()); const auto& gl_group = this->glo_.group(group.name());
return gl_group.max_lift_gas(); return gl_group.max_lift_gas();
} }
return std::nullopt; // If GLIFTOPT is missing from schedule, assume unlimited alq return std::nullopt; // If GLIFTOPT is missing from schedule, assume unlimited alq
} }
std::optional<double> template<class Scalar>
GasLiftStage2::getGroupMaxTotalGas_(const Group &group) std::optional<Scalar>
GasLiftStage2<Scalar>::getGroupMaxTotalGas_(const Group& group)
{ {
if (this->glo_.has_group(group.name())) { if (this->glo_.has_group(group.name())) {
const auto &gl_group = this->glo_.group(group.name()); const auto& gl_group = this->glo_.group(group.name());
return gl_group.max_total_gas(); return gl_group.max_total_gas();
} }
return std::nullopt; // If GLIFTOPT is missing from schedule, assume unlimited alq return std::nullopt; // If GLIFTOPT is missing from schedule, assume unlimited alq
@ -322,22 +334,23 @@ GasLiftStage2::getGroupMaxTotalGas_(const Group &group)
// NOTE: This means that wells are located at the leaf nodes of the tree, and // NOTE: This means that wells are located at the leaf nodes of the tree, and
// groups are located at the other nodes (not leaf nodes) of the tree // groups are located at the other nodes (not leaf nodes) of the tree
// //
std::vector<GasLiftSingleWellGeneric<double>*> template<class Scalar>
GasLiftStage2:: std::vector<GasLiftSingleWellGeneric<Scalar>*>
getGroupGliftWells_(const Group &group) GasLiftStage2<Scalar>::
getGroupGliftWells_(const Group& group)
{ {
std::vector<GasLiftSingleWell *> wells; std::vector<GasLiftSingleWell*> wells;
getGroupGliftWellsRecursive_(group, wells); getGroupGliftWellsRecursive_(group, wells);
return wells; return wells;
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
getGroupGliftWellsRecursive_(const Group &group, getGroupGliftWellsRecursive_(const Group& group,
std::vector<GasLiftSingleWell *> &wells) std::vector<GasLiftSingleWell*> &wells)
{ {
for (const std::string& group_name : group.groups()) { for (const std::string& group_name : group.groups()) {
if(this->schedule_.back().groups.has(group_name)) { if (this->schedule_.back().groups.has(group_name)) {
const Group& sub_group = const Group& sub_group =
this->schedule_.getGroup(group_name, this->report_step_idx_); this->schedule_.getGroup(group_name, this->report_step_idx_);
getGroupGliftWellsRecursive_(sub_group, wells); getGroupGliftWellsRecursive_(sub_group, wells);
@ -345,15 +358,15 @@ getGroupGliftWellsRecursive_(const Group &group,
} }
for (const std::string& well_name : group.wells()) { for (const std::string& well_name : group.wells()) {
if (this->stage1_wells_.count(well_name) == 1) { if (this->stage1_wells_.count(well_name) == 1) {
GasLiftSingleWell *well_ptr = this->stage1_wells_.at(well_name).get(); GasLiftSingleWell* well_ptr = this->stage1_wells_.at(well_name).get();
wells.push_back(well_ptr); wells.push_back(well_ptr);
} }
} }
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
mpiSyncGlobalGradVector_(std::vector<GradPair> &grads_global) const mpiSyncGlobalGradVector_(std::vector<GradPair>& grads_global) const
{ {
if (this->comm_.size() == 1) if (this->comm_.size() == 1)
return; return;
@ -367,10 +380,10 @@ mpiSyncGlobalGradVector_(std::vector<GradPair> &grads_global) const
mpiSyncLocalToGlobalGradVector_(grads_local, grads_global); mpiSyncLocalToGlobalGradVector_(grads_local, grads_global);
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
mpiSyncLocalToGlobalGradVector_( mpiSyncLocalToGlobalGradVector_(const std::vector<GradPair>& grads_local,
const std::vector<GradPair> &grads_local, std::vector<GradPair> &grads_global) const std::vector<GradPair>& grads_global) const
{ {
assert(this->comm_.size() > 1); // The parent should check if comm. size is > 1 assert(this->comm_.size() > 1); // The parent should check if comm. size is > 1
using Pair = std::pair<int, double>; using Pair = std::pair<int, double>;
@ -402,14 +415,14 @@ mpiSyncLocalToGlobalGradVector_(
for (std::size_t i = 0; i < grads_global_tmp.size(); ++i) { for (std::size_t i = 0; i < grads_global_tmp.size(); ++i) {
grads_global.emplace_back( grads_global.emplace_back(
std::make_pair( std::make_pair(
well_state_.globalIdxToWellName(grads_global_tmp[i].first), this->well_state_.globalIdxToWellName(grads_global_tmp[i].first),
grads_global_tmp[i].second)); grads_global_tmp[i].second));
} }
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
optimizeGroup_(const Group &group) optimizeGroup_(const Group& group)
{ {
const auto& group_name = group.name(); const auto& group_name = group.name();
const auto prod_control = this->group_state_.production_control(group_name); const auto prod_control = this->group_state_.production_control(group_name);
@ -433,9 +446,9 @@ optimizeGroup_(const Group &group)
} }
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
optimizeGroupsRecursive_(const Group &group) optimizeGroupsRecursive_(const Group& group)
{ {
for (const std::string& group_name : group.groups()) { for (const std::string& group_name : group.groups()) {
if(!this->schedule_.back().groups.has(group_name)) if(!this->schedule_.back().groups.has(group_name))
@ -447,15 +460,16 @@ optimizeGroupsRecursive_(const Group &group)
optimizeGroup_(group); optimizeGroup_(group);
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
recalculateGradientAndUpdateData_( recalculateGradientAndUpdateData_(GradPairItr& grad_itr,
GradPairItr &grad_itr, const std::string& gr_name_dont_limit, bool increase, const std::string& gr_name_dont_limit,
bool increase,
//incremental and decremental gradients, if 'grads' are incremental, then //incremental and decremental gradients, if 'grads' are incremental, then
// 'other_grads' are decremental, or conversely, if 'grads' are decremental, then // 'other_grads' are decremental, or conversely, if 'grads' are decremental, then
// 'other_grads' are incremental // 'other_grads' are incremental
std::vector<GradPair> &grads, std::vector<GradPair> &other_grads) std::vector<GradPair>& grads,
std::vector<GradPair>& other_grads)
{ {
// NOTE: We make a copy of the name string instead of taking a reference // NOTE: We make a copy of the name string instead of taking a reference
// since we may have to erase grad_itr (in the "else" condition below) // since we may have to erase grad_itr (in the "else" condition below)
@ -549,10 +563,12 @@ recalculateGradientAndUpdateData_(
// just do it once for the topmost group "PLAT-A" and then skip redistribution for // just do it once for the topmost group "PLAT-A" and then skip redistribution for
// all sub groups of "PLAT-A" // all sub groups of "PLAT-A"
// //
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
redistributeALQ_(std::vector<GasLiftSingleWell *> &wells, const Group &group, redistributeALQ_(std::vector<GasLiftSingleWell*>& wells,
std::vector<GradPair> &inc_grads, std::vector<GradPair> &dec_grads) const Group& group,
std::vector<GradPair>& inc_grads,
std::vector<GradPair>& dec_grads)
{ {
OptimizeState state {*this, group}; OptimizeState state {*this, group};
if (this->comm_.size() == 1) { if (this->comm_.size() == 1) {
@ -615,10 +631,11 @@ redistributeALQ_(std::vector<GasLiftSingleWell *> &wells, const Group &group,
// Lift gas increments are removed in turn from the well that currently has // Lift gas increments are removed in turn from the well that currently has
// the smallest weighted decremental gradient, until there is no surplus // the smallest weighted decremental gradient, until there is no surplus
// lift gas in the group. // lift gas in the group.
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
removeSurplusALQ_(const Group &group, removeSurplusALQ_(const Group& group,
std::vector<GradPair> &inc_grads, std::vector<GradPair> &dec_grads) std::vector<GradPair>& inc_grads,
std::vector<GradPair>& dec_grads)
{ {
if (dec_grads.empty()) { if (dec_grads.empty()) {
displayDebugMessage_("no wells to remove ALQ from. Skipping"); displayDebugMessage_("no wells to remove ALQ from. Skipping");
@ -644,8 +661,12 @@ removeSurplusALQ_(const Group &group,
displayDebugMessage_(msg); displayDebugMessage_(msg);
} }
SurplusState state {*this, group, oil_rate, gas_rate, water_rate, alq, SurplusState state {*this, group, oil_rate, gas_rate, water_rate, alq,
min_eco_grad, controls.oil_target, controls.gas_target, controls.water_target, static_cast<Scalar>(min_eco_grad),
controls.liquid_target, max_glift, max_totalgas }; static_cast<Scalar>(controls.oil_target),
static_cast<Scalar>(controls.gas_target),
static_cast<Scalar>(controls.water_target),
static_cast<Scalar>(controls.liquid_target),
max_glift, max_totalgas };
while (!stop_iteration) { while (!stop_iteration) {
if (dec_grads.size() >= 2) { if (dec_grads.size() >= 2) {
@ -709,9 +730,9 @@ removeSurplusALQ_(const Group &group,
} }
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
saveGrad_(GradMap &map, const std::string &name, GradInfo &grad) saveGrad_(GradMap& map, const std::string& name, GradInfo& grad)
{ {
if (auto it = map.find(name); it == map.end()) { if (auto it = map.find(name); it == map.end()) {
[[maybe_unused]] auto result = map.emplace(name, grad); [[maybe_unused]] auto result = map.emplace(name, grad);
@ -722,23 +743,23 @@ saveGrad_(GradMap &map, const std::string &name, GradInfo &grad)
} }
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
saveDecGrad_(const std::string &name, GradInfo &grad) saveDecGrad_(const std::string& name, GradInfo& grad)
{ {
saveGrad_(this->dec_grads_, name, grad); saveGrad_(this->dec_grads_, name, grad);
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
saveIncGrad_(const std::string &name, GradInfo &grad) saveIncGrad_(const std::string& name, GradInfo& grad)
{ {
saveGrad_(this->inc_grads_, name, grad); saveGrad_(this->inc_grads_, name, grad);
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
sortGradients_(std::vector<GradPair> &grads) sortGradients_(std::vector<GradPair>& grads)
{ {
auto cmp = [](GradPair a, GradPair b) { auto cmp = [](GradPair a, GradPair b) {
return a.second < b.second; return a.second < b.second;
@ -746,11 +767,12 @@ sortGradients_(std::vector<GradPair> &grads)
std::sort(grads.begin(), grads.end(), cmp); std::sort(grads.begin(), grads.end(), cmp);
} }
std::optional<GasLiftStage2::GradInfo> template<class Scalar>
GasLiftStage2:: std::optional<typename GasLiftStage2<Scalar>::GradInfo>
updateGrad_(const std::string &name, GradInfo &grad, bool increase) GasLiftStage2<Scalar>::
updateGrad_(const std::string& name, GradInfo& grad, bool increase)
{ {
GradMap &map = increase ? this->inc_grads_ : this->dec_grads_; GradMap& map = increase ? this->inc_grads_ : this->dec_grads_;
std::optional<GradInfo> old_value = std::nullopt; std::optional<GradInfo> old_value = std::nullopt;
if (map.count(name) == 1) { if (map.count(name) == 1) {
old_value = map.at(name); old_value = map.at(name);
@ -759,9 +781,11 @@ updateGrad_(const std::string &name, GradInfo &grad, bool increase)
return old_value; return old_value;
} }
void template<class Scalar>
GasLiftStage2:: void GasLiftStage2<Scalar>::
updateGradVector_(const std::string &name, std::vector<GradPair> &grads, double grad) updateGradVector_(const std::string& name,
std::vector<GradPair>& grads,
Scalar grad)
{ {
for (auto itr = grads.begin(); itr != grads.end(); itr++) { for (auto itr = grads.begin(); itr != grads.end(); itr++) {
if (itr->first == name) { if (itr->first == name) {
@ -778,14 +802,15 @@ updateGradVector_(const std::string &name, std::vector<GradPair> &grads, double
* Public methods declared in OptimizeState * Public methods declared in OptimizeState
***********************************************/ ***********************************************/
void template<class Scalar>
GasLiftStage2::OptimizeState:: void GasLiftStage2<Scalar>::OptimizeState::
calculateEcoGradients(std::vector<GasLiftSingleWell *> &wells, calculateEcoGradients(std::vector<GasLiftSingleWell*>& wells,
std::vector<GradPair> &inc_grads, std::vector<GradPair> &dec_grads) std::vector<GradPair>& inc_grads,
std::vector<GradPair>& dec_grads)
{ {
for (auto well_ptr : wells) { for (auto well_ptr : wells) {
const auto &gs_well = *well_ptr; // gs = GasLiftSingleWell const auto& gs_well = *well_ptr; // gs = GasLiftSingleWell
const auto &name = gs_well.name(); const auto& name = gs_well.name();
auto inc_grad = this->parent.calcIncOrDecGrad_(name, gs_well, group.name(), /*increase=*/true); auto inc_grad = this->parent.calcIncOrDecGrad_(name, gs_well, group.name(), /*increase=*/true);
if (inc_grad) { if (inc_grad) {
inc_grads.emplace_back(std::make_pair(name, inc_grad->grad)); inc_grads.emplace_back(std::make_pair(name, inc_grad->grad));
@ -799,10 +824,9 @@ calculateEcoGradients(std::vector<GasLiftSingleWell *> &wells,
} }
} }
template<class Scalar>
bool bool GasLiftStage2<Scalar>::OptimizeState::
GasLiftStage2::OptimizeState:: checkAtLeastTwoWells(std::vector<GasLiftSingleWell*>& wells)
checkAtLeastTwoWells(std::vector<GasLiftSingleWell *> &wells)
{ {
int numberOfwells = 0; int numberOfwells = 0;
for (auto well : wells){ for (auto well : wells){
@ -821,18 +845,20 @@ checkAtLeastTwoWells(std::vector<GasLiftSingleWell *> &wells)
return true; return true;
} }
void template<class Scalar>
GasLiftStage2::OptimizeState:: void GasLiftStage2<Scalar>::OptimizeState::
debugShowIterationInfo() debugShowIterationInfo()
{ {
const std::string msg = fmt::format("redistribute ALQ iteration {}", this->it); const std::string msg = fmt::format("redistribute ALQ iteration {}", this->it);
displayDebugMessage_(msg); displayDebugMessage_(msg);
} }
std::pair<std::optional<GasLiftStage2::GradPairItr>, template<class Scalar>
std::optional<GasLiftStage2::GradPairItr>> std::pair<std::optional<typename GasLiftStage2<Scalar>::GradPairItr>,
GasLiftStage2::OptimizeState:: std::optional<typename GasLiftStage2<Scalar>::GradPairItr>>
getEcoGradients(std::vector<GradPair> &inc_grads, std::vector<GradPair> &dec_grads) GasLiftStage2<Scalar>::OptimizeState::
getEcoGradients(std::vector<GradPair>& inc_grads,
std::vector<GradPair>& dec_grads)
{ {
if (!inc_grads.empty() && !dec_grads.empty()) { if (!inc_grads.empty() && !dec_grads.empty()) {
this->parent.sortGradients_(inc_grads); this->parent.sortGradients_(inc_grads);
@ -871,11 +897,12 @@ getEcoGradients(std::vector<GradPair> &inc_grads, std::vector<GradPair> &dec_gra
// a new decremental gradient given the new ALQ. The new incremental gradient // a new decremental gradient given the new ALQ. The new incremental gradient
// for this well is set equal to the current decremental gradient // for this well is set equal to the current decremental gradient
// (before the ALQ is subtracted) // (before the ALQ is subtracted)
void template<class Scalar>
GasLiftStage2::OptimizeState:: void GasLiftStage2<Scalar>::OptimizeState::
recalculateGradients( recalculateGradients(std::vector<GradPair>& inc_grads,
std::vector<GradPair> &inc_grads, std::vector<GradPair> &dec_grads, std::vector<GradPair>& dec_grads,
GradPairItr &min_dec_grad_itr, GradPairItr &max_inc_grad_itr) GradPairItr& min_dec_grad_itr,
GradPairItr& max_inc_grad_itr)
{ {
this->parent.recalculateGradientAndUpdateData_( this->parent.recalculateGradientAndUpdateData_(
max_inc_grad_itr, this->group.name(), /*increase=*/true, inc_grads, dec_grads); max_inc_grad_itr, this->group.name(), /*increase=*/true, inc_grads, dec_grads);
@ -888,9 +915,10 @@ recalculateGradients(
} }
// Take one ALQ increment from well1, and give it to well2 // Take one ALQ increment from well1, and give it to well2
void template<class Scalar>
GasLiftStage2::OptimizeState:: void GasLiftStage2<Scalar>::OptimizeState::
redistributeALQ( GradPairItr &min_dec_grad, GradPairItr &max_inc_grad) redistributeALQ(GradPairItr& min_dec_grad,
GradPairItr& max_inc_grad)
{ {
const std::string msg = fmt::format( const std::string msg = fmt::format(
"redistributing ALQ from well {} (dec gradient: {}) " "redistributing ALQ from well {} (dec gradient: {}) "
@ -908,16 +936,16 @@ GasLiftStage2::OptimizeState::
* Private methods declared in OptimizeState * Private methods declared in OptimizeState
**********************************************/ **********************************************/
void template<class Scalar>
GasLiftStage2::OptimizeState:: void GasLiftStage2<Scalar>::OptimizeState::
displayDebugMessage_(const std::string &msg) displayDebugMessage_(const std::string& msg)
{ {
this->parent.displayDebugMessage_(msg, this->group.name()); this->parent.displayDebugMessage_(msg, this->group.name());
} }
void template<class Scalar>
GasLiftStage2::OptimizeState:: void GasLiftStage2<Scalar>::OptimizeState::
displayWarning_(const std::string &msg) displayWarning_(const std::string& msg)
{ {
this->parent.displayWarning_(msg, this->group.name()); this->parent.displayWarning_(msg, this->group.name());
} }
@ -926,9 +954,11 @@ displayWarning_(const std::string &msg)
* Public methods declared in SurplusState * Public methods declared in SurplusState
**********************************************/ **********************************************/
void template<class Scalar>
GasLiftStage2::SurplusState:: void GasLiftStage2<Scalar>::SurplusState::
addOrRemoveALQincrement(GradMap &grad_map, const std::string& well_name, bool add) addOrRemoveALQincrement(GradMap& grad_map,
const std::string& well_name,
bool add)
{ {
if (this->parent.debug) { if (this->parent.debug) {
const std::string msg = fmt::format("group: {} : well {} : {} ALQ increment", const std::string msg = fmt::format("group: {} : well {} : {} ALQ increment",
@ -938,12 +968,12 @@ addOrRemoveALQincrement(GradMap &grad_map, const std::string& well_name, bool ad
this->parent.addOrRemoveALQincrement_(grad_map, well_name, add); this->parent.addOrRemoveALQincrement_(grad_map, well_name, add);
} }
bool template<class Scalar>
GasLiftStage2::SurplusState:: bool GasLiftStage2<Scalar>::SurplusState::
checkALQlimit() checkALQlimit()
{ {
if (this->max_glift) { if (this->max_glift) {
double max_alq = *(this->max_glift); Scalar max_alq = *(this->max_glift);
if ((max_alq) < (this->alq) ) { if ((max_alq) < (this->alq) ) {
if (this->parent.debug) { if (this->parent.debug) {
const std::string msg = fmt::format("group: {} : " const std::string msg = fmt::format("group: {} : "
@ -955,8 +985,8 @@ checkALQlimit()
} }
} }
if (this->max_total_gas) { if (this->max_total_gas) {
double max_total = *(this->max_total_gas); Scalar max_total = *(this->max_total_gas);
double total_gas_rate = (this->alq + this->gas_rate); Scalar total_gas_rate = (this->alq + this->gas_rate);
if ((max_total) < total_gas_rate ) { if ((max_total) < total_gas_rate ) {
if (this->parent.debug) { if (this->parent.debug) {
const std::string msg = fmt::format("group: {} : " const std::string msg = fmt::format("group: {} : "
@ -970,9 +1000,9 @@ checkALQlimit()
return false; return false;
} }
bool template<class Scalar>
GasLiftStage2::SurplusState:: bool GasLiftStage2<Scalar>::SurplusState::
checkEcoGradient(const std::string &well_name, double eco_grad) checkEcoGradient(const std::string& well_name, Scalar eco_grad)
{ {
if (eco_grad < this->min_eco_grad) { if (eco_grad < this->min_eco_grad) {
if (this->parent.debug) { if (this->parent.debug) {
@ -988,9 +1018,9 @@ checkEcoGradient(const std::string &well_name, double eco_grad)
} }
} }
bool template<class Scalar>
GasLiftStage2::SurplusState:: bool GasLiftStage2<Scalar>::SurplusState::
checkGasTarget(double delta_gas) checkGasTarget(Scalar delta_gas)
{ {
if (this->group.has_control(Group::ProductionCMode::GRAT)) { if (this->group.has_control(Group::ProductionCMode::GRAT)) {
// the change in gas rate is added to the gas rate to make sure the // the change in gas rate is added to the gas rate to make sure the
@ -1009,9 +1039,10 @@ checkGasTarget(double delta_gas)
} }
return false; return false;
} }
bool
GasLiftStage2::SurplusState:: template<class Scalar>
checkLiquidTarget(double delta_liquid) bool GasLiftStage2<Scalar>::SurplusState::
checkLiquidTarget(Scalar delta_liquid)
{ {
if (this->group.has_control(Group::ProductionCMode::LRAT)) { if (this->group.has_control(Group::ProductionCMode::LRAT)) {
// the change in liquid rate is added to the liquid rate to make sure the // the change in liquid rate is added to the liquid rate to make sure the
@ -1032,9 +1063,9 @@ checkLiquidTarget(double delta_liquid)
return false; return false;
} }
bool template<class Scalar>
GasLiftStage2::SurplusState:: bool GasLiftStage2<Scalar>::SurplusState::
checkOilTarget(double delta_oil) checkOilTarget(Scalar delta_oil)
{ {
if (this->group.has_control(Group::ProductionCMode::ORAT)) { if (this->group.has_control(Group::ProductionCMode::ORAT)) {
// the change in oil rate is added to the oil rate to make sure the // the change in oil rate is added to the oil rate to make sure the
@ -1054,10 +1085,10 @@ checkOilTarget(double delta_oil)
return false; return false;
} }
bool template<class Scalar>
GasLiftStage2::SurplusState:: bool GasLiftStage2<Scalar>::SurplusState::
checkWaterTarget(double delta_water) checkWaterTarget(Scalar delta_water)
{ {
if (this->group.has_control(Group::ProductionCMode::WRAT)) { if (this->group.has_control(Group::ProductionCMode::WRAT)) {
// the change in water rate is added to the water rate to make sure the // the change in water rate is added to the water rate to make sure the
// group still can produce its target // group still can produce its target
@ -1076,21 +1107,22 @@ checkWaterTarget(double delta_water)
return false; return false;
} }
std::array<double, 4> template<class Scalar>
GasLiftStage2::SurplusState:: std::array<Scalar, 4>
computeDelta(const std::string &well_name) GasLiftStage2<Scalar>::SurplusState::
computeDelta(const std::string& well_name)
{ {
std::array<double, 4> delta = {0.0, 0.0, 0.0, 0.0}; std::array<Scalar, 4> delta = {0.0, 0.0, 0.0, 0.0};
// compute the delta on wells on own rank // compute the delta on wells on own rank
if (this->parent.well_state_map_.count(well_name) > 0) { if (this->parent.well_state_map_.count(well_name) > 0) {
const GradInfo &gi = this->parent.dec_grads_.at(well_name); const GradInfo& gi = this->parent.dec_grads_.at(well_name);
GasLiftWellState<double>& state = *(this->parent.well_state_map_.at(well_name).get()); GasLiftWellState<Scalar>& state = *(this->parent.well_state_map_.at(well_name).get());
GasLiftSingleWell &gs_well = *(this->parent.stage1_wells_.at(well_name).get()); GasLiftSingleWell& gs_well = *(this->parent.stage1_wells_.at(well_name).get());
const WellInterfaceGeneric &well = gs_well.getWell(); const WellInterfaceGeneric& well = gs_well.getWell();
// only get deltas for wells owned by this rank // only get deltas for wells owned by this rank
if (this->parent.well_state_.wellIsOwned(well.indexOfWell(), well_name)) { if (this->parent.well_state_.wellIsOwned(well.indexOfWell(), well_name)) {
const auto &well_ecl = well.wellEcl(); const auto& well_ecl = well.wellEcl();
double factor = well_ecl.getEfficiencyFactor(); Scalar factor = well_ecl.getEfficiencyFactor();
auto& [delta_oil, delta_gas, delta_water, delta_alq] = delta; auto& [delta_oil, delta_gas, delta_water, delta_alq] = delta;
delta_oil = factor * (gi.new_oil_rate - state.oilRate()); delta_oil = factor * (gi.new_oil_rate - state.oilRate());
delta_gas = factor * (gi.new_gas_rate - state.gasRate()); delta_gas = factor * (gi.new_gas_rate - state.gasRate());
@ -1105,9 +1137,9 @@ computeDelta(const std::string &well_name)
return delta; return delta;
} }
void template<class Scalar>
GasLiftStage2::SurplusState:: void GasLiftStage2<Scalar>::SurplusState::
updateRates(const std::array<double, 4>& delta) updateRates(const std::array<Scalar, 4>& delta)
{ {
const auto& [delta_oil, delta_gas, delta_water, delta_alq] = delta; const auto& [delta_oil, delta_gas, delta_water, delta_alq] = delta;
this->oil_rate += delta_oil; this->oil_rate += delta_oil;
@ -1116,5 +1148,6 @@ updateRates(const std::array<double, 4>& delta)
this->alq += delta_alq; this->alq += delta_alq;
} }
template class GasLiftStage2<double>;
} // namespace Opm } // namespace Opm

268
opm/simulators/wells/GasLiftStage2.hpp
View File

@ -30,8 +30,7 @@
#include <tuple> #include <tuple>
#include <vector> #include <vector>
namespace Opm namespace Opm {
{
class DeferredLogger; class DeferredLogger;
class GasLiftOpt; class GasLiftOpt;
@ -42,86 +41,113 @@ class Schedule;
class WellInterfaceGeneric; class WellInterfaceGeneric;
template<class Scalar> class WellState; template<class Scalar> class WellState;
class GasLiftStage2 : public GasLiftCommon<double> { template<class Scalar>
using GasLiftSingleWell = GasLiftSingleWellGeneric<double>; class GasLiftStage2 : public GasLiftCommon<Scalar>
{
using GasLiftSingleWell = GasLiftSingleWellGeneric<Scalar>;
using GLiftOptWells = std::map<std::string,std::unique_ptr<GasLiftSingleWell>>; using GLiftOptWells = std::map<std::string,std::unique_ptr<GasLiftSingleWell>>;
using GLiftProdWells = std::map<std::string,const WellInterfaceGeneric*>; using GLiftProdWells = std::map<std::string,const WellInterfaceGeneric*>;
using GLiftWellStateMap = std::map<std::string,std::unique_ptr<GasLiftWellState<double>>>; using GLiftWellStateMap = std::map<std::string,std::unique_ptr<GasLiftWellState<Scalar>>>;
using GradPair = std::pair<std::string, double>; using GradPair = std::pair<std::string, Scalar>;
using GradPairItr = std::vector<GradPair>::iterator; using GradPairItr = typename std::vector<GradPair>::iterator;
using GradInfo = typename GasLiftSingleWellGeneric<double>::GradInfo; using GradInfo = typename GasLiftSingleWellGeneric<Scalar>::GradInfo;
using GradMap = std::map<std::string, GradInfo>; using GradMap = std::map<std::string, GradInfo>;
using MPIComm = typename Dune::MPIHelper::MPICommunicator; using MessageType = typename GasLiftCommon<Scalar>::MessageType;
static const int Water = BlackoilPhases::Aqua; static const int Water = BlackoilPhases::Aqua;
static const int Oil = BlackoilPhases::Liquid; static const int Oil = BlackoilPhases::Liquid;
static const int Gas = BlackoilPhases::Vapour; static const int Gas = BlackoilPhases::Vapour;
public: public:
GasLiftStage2( GasLiftStage2(const int report_step_idx,
const int report_step_idx, const Parallel::Communication& comm,
const Parallel::Communication& comm, const Schedule& schedule,
const Schedule& schedule, const SummaryState& summary_state,
const SummaryState& summary_state, DeferredLogger& deferred_logger,
DeferredLogger& deferred_logger, WellState<Scalar>& well_state,
WellState<double>& well_state, const GroupState<Scalar>& group_state,
const GroupState<double>& group_state, GLiftProdWells& prod_wells,
GLiftProdWells& prod_wells, GLiftOptWells& glift_wells,
GLiftOptWells& glift_wells, GasLiftGroupInfo<Scalar>& group_info,
GasLiftGroupInfo<double>& group_info, GLiftWellStateMap& state_map,
GLiftWellStateMap& state_map, bool glift_debug);
bool glift_debug
);
void runOptimize(); void runOptimize();
protected: protected:
void addOrRemoveALQincrement_( void addOrRemoveALQincrement_(GradMap& grad_map,
GradMap& grad_map, const std::string& well_name, bool add); const std::string& well_name,
std::optional<GradInfo> calcIncOrDecGrad_( bool add);
const std::string name, const GasLiftSingleWell& gs_well, const std::string& gr_name_dont_limit, bool increase);
std::optional<GradInfo> calcIncOrDecGrad_(const std::string name,
const GasLiftSingleWell& gs_well,
const std::string& gr_name_dont_limit,
bool increase);
bool checkRateAlreadyLimited_(const std::string& well_name, bool checkRateAlreadyLimited_(const std::string& well_name,
GasLiftWellState<double>& state, GasLiftWellState<Scalar>& state,
bool increase); bool increase);
GradInfo deleteDecGradItem_(const std::string& name); GradInfo deleteDecGradItem_(const std::string& name);
GradInfo deleteIncGradItem_(const std::string& name); GradInfo deleteIncGradItem_(const std::string& name);
GradInfo deleteGrad_(const std::string& name, bool increase); GradInfo deleteGrad_(const std::string& name, bool increase);
void displayDebugMessage_(const std::string& msg) const override; void displayDebugMessage_(const std::string& msg) const override;
void displayDebugMessage2B_(const std::string& msg); void displayDebugMessage2B_(const std::string& msg);
void displayDebugMessage_(const std::string& msg, const std::string& group_name); void displayDebugMessage_(const std::string& msg,
void displayWarning_(const std::string& msg, const std::string& group_name); const std::string& group_name);
void displayWarning_(const std::string& msg,
const std::string& group_name);
void displayWarning_(const std::string& msg); void displayWarning_(const std::string& msg);
std::tuple<double, double, double, double> getCurrentGroupRates_(const Group& group);
std::optional<double> getGroupMaxALQ_(const Group &group); std::tuple<Scalar, Scalar, Scalar, Scalar>
std::optional<double> getGroupMaxTotalGas_(const Group &group); getCurrentGroupRates_(const Group& group);
std::vector<GasLiftSingleWell *> getGroupGliftWells_(
const Group& group); std::optional<Scalar> getGroupMaxALQ_(const Group& group);
void getGroupGliftWellsRecursive_( std::optional<Scalar> getGroupMaxTotalGas_(const Group& group);
const Group& group, std::vector<GasLiftSingleWell *>& wells);
std::vector<GasLiftSingleWell*> getGroupGliftWells_(const Group& group);
void getGroupGliftWellsRecursive_(const Group& group,
std::vector<GasLiftSingleWell*>& wells);
void optimizeGroup_(const Group& group); void optimizeGroup_(const Group& group);
void optimizeGroupsRecursive_(const Group& group); void optimizeGroupsRecursive_(const Group& group);
void recalculateGradientAndUpdateData_(
GradPairItr& grad_itr, const std::string& gr_name_dont_limit, bool increase, void recalculateGradientAndUpdateData_(GradPairItr& grad_itr,
std::vector<GradPair>& grads, std::vector<GradPair>& other_grads); const std::string& gr_name_dont_limit,
void redistributeALQ_( bool increase,
std::vector<GasLiftSingleWell *>& wells, const Group& group, std::vector<GradPair>& grads,
std::vector<GradPair>& inc_grads, std::vector<GradPair>& dec_grads); std::vector<GradPair>& other_grads);
void removeSurplusALQ_(
const Group& group, void redistributeALQ_(std::vector<GasLiftSingleWell*>& wells,
std::vector<GradPair>& inc_grads, std::vector<GradPair>& dec_grads); const Group& group,
std::vector<GradPair>& inc_grads,
std::vector<GradPair>& dec_grads);
void removeSurplusALQ_(const Group& group,
std::vector<GradPair>& inc_grads,
std::vector<GradPair>& dec_grads);
void saveGrad_(GradMap& map, const std::string& name, GradInfo& grad); void saveGrad_(GradMap& map, const std::string& name, GradInfo& grad);
void saveDecGrad_(const std::string& name, GradInfo& grad); void saveDecGrad_(const std::string& name, GradInfo& grad);
void saveIncGrad_(const std::string& name, GradInfo& grad); void saveIncGrad_(const std::string& name, GradInfo& grad);
void sortGradients_(std::vector<GradPair>& grads); void sortGradients_(std::vector<GradPair>& grads);
std::optional<GradInfo> updateGrad_(
const std::string& name, GradInfo& grad, bool increase);
void updateGradVector_(
const std::string& name, std::vector<GradPair>& grads, double grad);
void mpiSyncGlobalGradVector_(std::vector<GradPair>& grads_global) const;
void mpiSyncLocalToGlobalGradVector_(
const std::vector<GradPair>& grads_local,
std::vector<GradPair>& grads_global) const;
std::optional<GradInfo> updateGrad_(const std::string& name,
GradInfo& grad, bool increase);
void updateGradVector_(const std::string& name,
std::vector<GradPair>& grads,
Scalar grad);
void mpiSyncGlobalGradVector_(std::vector<GradPair>& grads_global) const;
void mpiSyncLocalToGlobalGradVector_(const std::vector<GradPair>& grads_local,
std::vector<GradPair>& grads_global) const;
GLiftProdWells& prod_wells_; GLiftProdWells& prod_wells_;
GLiftOptWells& stage1_wells_; GLiftOptWells& stage1_wells_;
GasLiftGroupInfo<double>& group_info_; GasLiftGroupInfo<Scalar>& group_info_;
GLiftWellStateMap& well_state_map_; GLiftWellStateMap& well_state_map_;
int report_step_idx_; int report_step_idx_;
@ -133,12 +159,14 @@ protected:
int max_iterations_ = 1000; int max_iterations_ = 1000;
//int time_step_idx_; //int time_step_idx_;
struct OptimizeState { struct OptimizeState
OptimizeState(GasLiftStage2& parent_, const Group& group_ ) : {
parent{parent_}, OptimizeState(GasLiftStage2& parent_, const Group& group_)
group{group_}, : parent{parent_}
it{0} , group{group_}
, it{0}
{} {}
GasLiftStage2& parent; GasLiftStage2& parent;
const Group& group; const Group& group;
int it; int it;
@ -147,70 +175,90 @@ protected:
using GradPair = typename GasLiftStage2::GradPair; using GradPair = typename GasLiftStage2::GradPair;
using GradPairItr = typename GasLiftStage2::GradPairItr; using GradPairItr = typename GasLiftStage2::GradPairItr;
using GradMap = typename GasLiftStage2::GradMap; using GradMap = typename GasLiftStage2::GradMap;
void calculateEcoGradients(std::vector<GasLiftSingleWell *>& wells,
std::vector<GradPair>& inc_grads, std::vector<GradPair>& dec_grads); void calculateEcoGradients(std::vector<GasLiftSingleWell*>& wells,
bool checkAtLeastTwoWells(std::vector<GasLiftSingleWell *>& wells); std::vector<GradPair>& inc_grads,
std::vector<GradPair>& dec_grads);
bool checkAtLeastTwoWells(std::vector<GasLiftSingleWell*>& wells);
void debugShowIterationInfo(); void debugShowIterationInfo();
std::pair<std::optional<GradPairItr>,std::optional<GradPairItr>> std::pair<std::optional<GradPairItr>,std::optional<GradPairItr>>
getEcoGradients( getEcoGradients(std::vector<GradPair>& inc_grads,
std::vector<GradPair>& inc_grads, std::vector<GradPair>& dec_grads); std::vector<GradPair>& dec_grads);
void recalculateGradients(
std::vector<GradPair>& inc_grads, std::vector<GradPair>& dec_grads, void recalculateGradients(std::vector<GradPair>& inc_grads,
GradPairItr& min_dec_grad_itr, GradPairItr &max_inc_grad_itr); std::vector<GradPair>& dec_grads,
void redistributeALQ( GradPairItr& min_dec_grad, GradPairItr& max_inc_grad); GradPairItr& min_dec_grad_itr,
GradPairItr &max_inc_grad_itr);
void redistributeALQ( GradPairItr& min_dec_grad,
GradPairItr& max_inc_grad);
private: private:
void displayDebugMessage_(const std::string& msg); void displayDebugMessage_(const std::string& msg);
void displayWarning_(const std::string& msg); void displayWarning_(const std::string& msg);
}; };
struct SurplusState { struct SurplusState
SurplusState( GasLiftStage2& parent_, const Group& group_, {
double oil_rate_, double gas_rate_, double water_rate_, double alq_, SurplusState(GasLiftStage2& parent_,
double min_eco_grad_, const Group& group_,
double oil_target_, double gas_target_, double water_target_, double liquid_target_, Scalar oil_rate_,
std::optional<double> max_glift_, std::optional<double> max_total_gas_) : Scalar gas_rate_,
parent{parent_}, Scalar water_rate_,
group{group_}, Scalar alq_,
oil_rate{oil_rate_}, Scalar min_eco_grad_,
gas_rate{gas_rate_}, Scalar oil_target_,
water_rate{water_rate_}, Scalar gas_target_,
alq{alq_}, Scalar water_target_,
min_eco_grad{min_eco_grad_}, Scalar liquid_target_,
oil_target{oil_target_}, std::optional<Scalar> max_glift_,
gas_target{gas_target_}, std::optional<Scalar> max_total_gas_)
water_target(water_target_), : parent{parent_}
liquid_target{liquid_target_}, , group{group_}
max_glift{max_glift_}, , oil_rate{oil_rate_}
max_total_gas{max_total_gas_}, , gas_rate{gas_rate_}
it{0} , water_rate{water_rate_}
, alq{alq_}
, min_eco_grad{min_eco_grad_}
, oil_target{oil_target_}
, gas_target{gas_target_}
, water_target(water_target_)
, liquid_target{liquid_target_}
, max_glift{max_glift_}
, max_total_gas{max_total_gas_}
, it{0}
{} {}
GasLiftStage2 &parent; GasLiftStage2 &parent;
const Group &group; const Group &group;
double oil_rate; Scalar oil_rate;
double gas_rate; Scalar gas_rate;
double water_rate; Scalar water_rate;
double alq; Scalar alq;
const double min_eco_grad; const Scalar min_eco_grad;
const double oil_target; const Scalar oil_target;
const double gas_target; const Scalar gas_target;
const double water_target; const Scalar water_target;
const double liquid_target; const Scalar liquid_target;
std::optional<double> max_glift; std::optional<Scalar> max_glift;
std::optional<double> max_total_gas; std::optional<Scalar> max_total_gas;
int it; int it;
void addOrRemoveALQincrement( void addOrRemoveALQincrement(GradMap &grad_map,
GradMap &grad_map, const std::string& well_name, bool add); const std::string& well_name,
bool add);
bool checkALQlimit(); bool checkALQlimit();
bool checkEcoGradient(const std::string& well_name, double eco_grad); bool checkEcoGradient(const std::string& well_name, Scalar eco_grad);
bool checkGasTarget(double delta_gas); bool checkGasTarget(Scalar delta_gas);
bool checkLiquidTarget(double delta_liquid); bool checkLiquidTarget(Scalar delta_liquid);
bool checkOilTarget(double delta_oil); bool checkOilTarget(Scalar delta_oil);
bool checkWaterTarget(double delta_water); bool checkWaterTarget(Scalar delta_water);
std::array<double, 4> computeDelta(const std::string& name); std::array<Scalar, 4> computeDelta(const std::string& name);
void updateRates(const std::array<double, 4>& delta); void updateRates(const std::array<Scalar, 4>& delta);
}; };
}; };