55 static constexpr int Water = BlackoilPhases::Aqua;
56 static constexpr int Oil = BlackoilPhases::Liquid;
57 static constexpr int Gas = BlackoilPhases::Vapour;
58 static constexpr int NUM_PHASES = 3;
59 static constexpr double ALQ_EPSILON = 1
e-8;
62 using GLiftSyncGroups = std::set<int>;
63 using Rate = GasLiftGroupInfo::Rate;
71 double alq_,
bool alq_is_limited_) :
74 oil_is_limited{oil_is_limited_},
76 gas_is_limited{gas_is_limited_},
78 water_is_limited{water_is_limited_},
80 alq_is_limited{alq_is_limited_} {}
86 double new_water_rate;
87 bool water_is_limited;
94 const std::string& name()
const {
return well_name_; }
96 std::optional<GradInfo> calcIncOrDecGradient(
double oil_rate,
double gas_rate,
104 std::unique_ptr<GasLiftWellState> runOptimize(
const int iteration_idx);
120 const Parallel::Communication& comm,
131 bhp_is_limited{rates.bhp_is_limited}
143 bhp_is_limited = rates.bhp_is_limited;
150 double operator[](Rate
rate_type)
const {
159 return this->oil + this->water;
161 throw std::runtime_error(
"This should not happen");
165 double oil, gas, water;
171 enum class LimitType {well, group, none};
174 bool oil_is_limited_,
bool gas_is_limited_,
180 oil_is_limited{oil_is_limited_},
181 gas_is_limited{gas_is_limited_},
182 water_is_limited{water_is_limited_},
186 set_initial_limit_type_();
191 bool oil_is_limited_,
bool gas_is_limited_,
192 bool water_is_limited_
195 oil_is_limited{oil_is_limited_},
196 gas_is_limited{gas_is_limited_},
197 water_is_limited{water_is_limited_}
199 set_initial_limit_type_();
202 bool limited()
const {
203 return oil_is_limited || gas_is_limited || water_is_limited;
207 LimitType limit_type;
210 bool water_is_limited;
211 std::optional<Rate> oil_limiting_target;
212 std::optional<Rate> water_limiting_target;
214 void set_initial_limit_type_() {
215 limit_type = limited() ? LimitType::well : LimitType::none;
225 stop_iteration{
false},
235 std::pair<std::optional<double>,
bool> addOrSubtractAlqIncrement(
double alq);
236 double calcEcoGradient(
double oil_rate,
double new_oil_rate,
237 double gas_rate,
double new_gas_rate);
238 bool checkAlqOutsideLimits(
double alq,
double oil_rate);
239 bool checkEcoGradient(
double gradient);
240 bool checkOilRateExceedsTarget(
double oil_rate);
241 bool checkRatesViolated(
const LimitedRates& rates)
const;
242 void debugShowIterationInfo(
double alq);
243 double getBhpWithLimit();
244 void warn_(std::string
msg) {parent.displayWarning_(
msg);}
249 std::pair<std::optional<double>,
bool> addOrSubtractAlqIncrement_(
250 double alq,
bool increase)
const;
251 double calcEcoGradient_(
double oil_rate,
double new_oil_rate,
252 double gas_rate,
double new_gas_rate,
bool increase)
const;
253 bool checkALQequal_(
double alq1,
double alq2)
const;
254 bool checkGroupTargetsViolated(
256 bool checkInitialALQmodified_(
double alq,
double initial_alq)
const;
257 virtual bool checkThpControl_()
const = 0;
258 virtual std::optional<double> computeBhpAtThpLimit_(
double alq,
bool debug_output =
true)
const = 0;
259 std::pair<std::optional<double>,
double> computeConvergedBhpAtThpLimitByMaybeIncreasingALQ_()
const;
260 std::pair<std::optional<BasicRates>,
double> computeInitialWellRates_()
const;
261 std::optional<LimitedRates> computeLimitedWellRatesWithALQ_(
double alq)
const;
262 virtual BasicRates computeWellRates_(
double bhp,
bool bhp_is_limited,
bool debug_output =
true)
const = 0;
263 std::optional<BasicRates> computeWellRatesWithALQ_(
double alq)
const;
264 void debugCheckNegativeGradient_(
double grad,
double alq,
double new_alq,
265 double oil_rate,
double new_oil_rate,
266 double gas_rate,
double new_gas_rate,
267 bool increase)
const;
268 void debugPrintWellStateRates()
const;
269 void debugShowAlqIncreaseDecreaseCounts_();
270 void debugShowBhpAlqTable_();
271 void debugShowLimitingTargets_(
const LimitedRates& rates)
const;
272 void debugShowProducerControlMode()
const;
273 void debugShowStartIteration_(
double alq,
bool increase,
double oil_rate);
274 void debugShowTargets_();
275 void displayDebugMessage_(
const std::string&
msg)
const override;
276 void displayWarning_(
const std::string& warning);
277 std::pair<double, bool> getBhpWithLimit_(
double bhp)
const;
278 std::pair<double, bool> getGasRateWithLimit_(
280 std::pair<double, bool> getGasRateWithGroupLimit_(
282 std::pair<std::optional<LimitedRates>,
double> getInitialRatesWithLimit_()
const;
284 std::tuple<double,double,bool,bool> getLiquidRateWithGroupLimit_(
285 const double new_oil_rate,
const double oil_rate,
286 const double new_water_rate,
const double water_rate,
const std::string&
gr_name_dont_limit)
const;
287 std::pair<double, bool> getOilRateWithGroupLimit_(
289 std::pair<double, bool> getOilRateWithLimit_(
const BasicRates& rates)
const;
290 std::pair<double, std::optional<Rate>> getOilRateWithLimit2_(
292 double getProductionTarget_(Rate
rate)
const;
294 std::pair<double, std::optional<Rate>> getRateWithLimit_(
296 std::tuple<double, const std::string*, double> getRateWithGroupLimit_(
298 std::pair<double, bool> getWaterRateWithGroupLimit_(
299 double new_water_rate,
double water_rate,
const std::string&
gr_name_dont_limit)
const;
300 std::pair<double, bool> getWaterRateWithLimit_(
const BasicRates& rates)
const;
301 std::pair<double, std::optional<Rate>> getWaterRateWithLimit2_(
304 bool hasProductionControl_(Rate
rate)
const;
305 std::pair<LimitedRates, double> increaseALQtoPositiveOilRate_(
307 std::pair<LimitedRates, double> increaseALQtoMinALQ_(
309 void logSuccess_(
double alq,
311 std::pair<LimitedRates, double> maybeAdjustALQbeforeOptimizeLoop_(
312 const LimitedRates& rates,
double alq,
bool increase)
const;
313 std::pair<LimitedRates, double> reduceALQtoGroupAlqLimits_(
315 std::pair<LimitedRates, double> reduceALQtoGroupTarget(
317 std::pair<LimitedRates, double> reduceALQtoWellTarget_(
319 std::unique_ptr<GasLiftWellState> runOptimize1_();
320 std::unique_ptr<GasLiftWellState> runOptimize2_();
321 std::unique_ptr<GasLiftWellState> runOptimizeLoop_(
bool increase);
323 std::unique_ptr<GasLiftWellState> tryIncreaseLiftGas_();
324 std::unique_ptr<GasLiftWellState> tryDecreaseLiftGas_();
325 void updateGroupRates_(
331 void updateWellStateAlqFixedValue_(
const GasLiftWell& well);
333 void debugInfoGroupRatesExceedTarget(
335 void warnMaxIterationsExceeded_();
337 const Well& ecl_well_;
341 GLiftSyncGroups& sync_groups_;
359 std::string well_name_;
364 bool debug_limit_increase_decrease_;
365 bool debug_abort_if_decrease_and_oil_is_limited_ =
false;
366 bool debug_abort_if_increase_and_gas_is_limited_ =
false;