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

Subtract reduction rate when checking higher group constraints

Browse Source
This commit is contained in:
Stein Krogstad 2024-11-12 13:51:17 +01:00
parent b502206ad7
commit cf877332ab
2 changed files with 29 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/simulators/wells/BlackoilWellModelGeneric.cpp
View File

@ -612,6 +612,10 @@ checkGroupHigherConstraints(const Group& group,
std::vector<Scalar> resv_coeff_inj(phase_usage_.num_phases, 0.0);
calcInjResvCoeff(fipnum, pvtreg, resv_coeff_inj);
// checkGroupConstraintsInj considers 'available' rates (e.g., group rates minus reduction rates).
// So when checking constraints, current groups rate must also be subtracted it's reductioin rate
const std::vector<Scalar> reduction_rates = this->groupState().injection_reduction_rates(group.name());
for (int phasePos = 0; phasePos < phase_usage_.num_phases; ++phasePos) {
const Scalar local_current_rate = WellGroupHelpers<Scalar>::sumWellSurfaceRates(group,
schedule(),
@ -620,7 +624,7 @@ checkGroupHigherConstraints(const Group& group,
phasePos,
/* isInjector */ true);
// Sum over all processes
rates[phasePos] = comm_.sum(local_current_rate);
rates[phasePos] = comm_.sum(local_current_rate) - reduction_rates[phasePos];
}
const Phase all[] = { Phase::WATER, Phase::OIL, Phase::GAS };
for (Phase phase : all) {
@ -665,6 +669,10 @@ checkGroupHigherConstraints(const Group& group,
if (!isField && group.isProductionGroup()) {
// Obtain rates for group.
// checkGroupConstraintsProd considers 'available' rates (e.g., group rates minus reduction rates).
// So when checking constraints, current groups rate must also be subtracted it's reductioin rate
const std::vector<Scalar> reduction_rates = this->groupState().production_reduction_rates(group.name());
for (int phasePos = 0; phasePos < phase_usage_.num_phases; ++phasePos) {
const Scalar local_current_rate = WellGroupHelpers<Scalar>::sumWellSurfaceRates(group,
schedule(),
@ -673,7 +681,7 @@ checkGroupHigherConstraints(const Group& group,
phasePos,
/* isInjector */ false);
// Sum over all processes
rates[phasePos] = -comm_.sum(local_current_rate);
rates[phasePos] = -comm_.sum(local_current_rate) - reduction_rates[phasePos];
}
std::vector<Scalar> resv_coeff(phase_usage_.num_phases, 0.0);
calcResvCoeff(fipnum, pvtreg, this->groupState().production_rates(group.name()), resv_coeff);

31
opm/simulators/wells/WellGroupHelpers.cpp
View File

@ -1181,6 +1181,8 @@ checkGroupConstraintsProd(const std::string& name,
// part of. Later it is the accumulated factor including the group efficiency factor
// of the child of group.
// NOTE: if name is a group, the rates-argument is the portion of the 'full' group rates that
// is potentially available for control by an ancestor, i.e. full rates minus reduction rates
const Group::ProductionCMode& currentGroupControl = group_state.production_control(group.name());
if (currentGroupControl == Group::ProductionCMode::FLD || currentGroupControl == Group::ProductionCMode::NONE) {
@ -1263,7 +1265,7 @@ checkGroupConstraintsProd(const std::string& name,
// Assume we have a chain of groups as follows: BOTTOM -> MIDDLE -> TOP.
// Then ...
// TODO finish explanation.
const Scalar current_rate
const Scalar current_rate_available
= -tcalc.calcModeRateFromRates(rates); // Switch sign since 'rates' are negative for producers.
const auto chain = groupChainTopBot(name, group.name(), schedule, reportStepIdx);
// Because 'name' is the last of the elements, and not an ancestor, we subtract one below.
@ -1298,6 +1300,7 @@ checkGroupConstraintsProd(const std::string& name,
}
}
}
// Compute portion of target corresponding to current_rate_available
Scalar target = orig_target;
for (std::size_t ii = 0; ii < num_ancestors; ++ii) {
if ((ii == 0) || guideRate->has(chain[ii])) {
@ -1309,18 +1312,18 @@ checkGroupConstraintsProd(const std::string& name,
}
// Add my reduction back at the level where it is included in the local reduction
if (local_reduction_level == ii ) {
target += current_rate * efficiencyFactor;
target += current_rate_available * efficiencyFactor;
}
}
target *= localFraction(chain[ii + 1]);
}
// Avoid negative target rates comming from too large local reductions.
const Scalar target_rate = std::max(Scalar(1e-12), target / efficiencyFactor);
const Scalar target_rate_available = std::max(Scalar(1e-12), target / efficiencyFactor);
Scalar scale = 1.0;
if (current_rate > 1e-12)
scale = target_rate / current_rate;
if (current_rate_available > 1e-12)
scale = target_rate_available / current_rate_available;
return std::make_pair(current_rate > target_rate, scale);
return std::make_pair(current_rate_available > target_rate_available, scale);
}
template<class Scalar>
@ -1350,6 +1353,9 @@ checkGroupConstraintsInj(const std::string& name,
// part of. Later it is the accumulated factor including the group efficiency factor
// of the child of group.
// NOTE: if name is a group, the rates-argument is the portion of the 'full' group rates that
// is potentially available for control by an ancestor, i.e. full rates minus reduction rates
auto currentGroupControl = group_state.injection_control(group.name(), injectionPhase);
if (currentGroupControl == Group::InjectionCMode::FLD || currentGroupControl == Group::InjectionCMode::NONE) {
@ -1435,7 +1441,7 @@ checkGroupConstraintsInj(const std::string& name,
// Assume we have a chain of groups as follows: BOTTOM -> MIDDLE -> TOP.
// Then ...
// TODO finish explanation.
const Scalar current_rate
const Scalar current_rate_available
= tcalc.calcModeRateFromRates(rates); // Switch sign since 'rates' are negative for producers.
const auto chain = groupChainTopBot(name, group.name(), schedule, reportStepIdx);
// Because 'name' is the last of the elements, and not an ancestor, we subtract one below.
@ -1473,6 +1479,7 @@ checkGroupConstraintsInj(const std::string& name,
}
}
// Compute portion of target corresponding to current_rate_available
Scalar target = orig_target;
for (std::size_t ii = 0; ii < num_ancestors; ++ii) {
if ((ii == 0) || guideRate->has(chain[ii], injectionPhase)) {
@ -1485,18 +1492,18 @@ checkGroupConstraintsInj(const std::string& name,
// Add my reduction back at the level where it is included in the local reduction
if (local_reduction_level == ii ) {
target += current_rate * efficiencyFactor;
target += current_rate_available * efficiencyFactor;
}
}
target *= localFraction(chain[ii + 1]);
}
// Avoid negative target rates comming from too large local reductions.
const Scalar target_rate = std::max(Scalar(1e-12), target / efficiencyFactor);
const Scalar target_rate_available = std::max(Scalar(1e-12), target / efficiencyFactor);
Scalar scale = 1.0;
if (current_rate > 1e-12)
scale = target_rate / current_rate;
if (current_rate_available > 1e-12)
scale = target_rate_available / current_rate_available;
return std::make_pair(current_rate > target_rate, scale);
return std::make_pair(current_rate_available > target_rate_available, scale);
}
template<class Scalar>