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

Check if guiderate is violated

Browse Source 
If guiderate is violated change to group controll.
Note that a factor 1.01 is added to minimize oscilations.
Fix missing  multiplication with group efficiency when accumulating guiderates
This commit is contained in:
Tor Harald Sandve
2022-02-15 08:53:59 +01:00
parent 161e72f6de
commit b931c5ef97
1 changed files with 40 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -1040,7 +1040,8 @@ namespace WellGroupHelpers
// We are a group, with default guide rate.
// Compute guide rate by accumulating our children's guide rates.
const Group& group = schedule_.getGroup(name, report_step_);
return guideRateSum(group, always_included_child);
const double eff = group.getGroupEfficiencyFactor();
return eff * guideRateSum(group, always_included_child);
}
} else {
// No group-controlled subordinate wells.
@@ -1164,6 +1165,10 @@ namespace WellGroupHelpers
return tcalc.calcModeRateFromRates(groupTargetReductions);
};
auto localCurrentRate = [&](const std::string& group_name) {
const std::vector<double>& groupSurfaceRates = group_state.production_rates(group_name);
return tcalc.calcModeRateFromRates(groupSurfaceRates);
};
const double orig_target = tcalc.groupTarget(group.productionControls(summaryState));
// Assume we have a chain of groups as follows: BOTTOM -> MIDDLE -> TOP.
// Then ...
@@ -1181,6 +1186,20 @@ namespace WellGroupHelpers
}
}
// check whether guide rate is violated
if (local_reduction_level > 0) {
const auto& guided_group = chain[local_reduction_level];
const double grefficiency
= schedule.getGroup(guided_group, reportStepIdx).getGroupEfficiencyFactor();
const double currentRateFraction = grefficiency * localCurrentRate(guided_group) / (localCurrentRate(chain[local_reduction_level-1]));
const double guiderateFraction = localFraction(guided_group);
// we add a factor here to avoid switching due to numerical instability
const double factor = 1.01;
if (currentRateFraction > (guiderateFraction * factor)) {
return std::make_pair(true, guiderateFraction/currentRateFraction);
}
}
double target = orig_target;
for (size_t ii = 0; ii < num_ancestors; ++ii) {
if ((ii == 0) || guideRate->has(chain[ii])) {
@@ -1323,6 +1342,11 @@ namespace WellGroupHelpers
return tcalc.calcModeRateFromRates(groupTargetReductions);
};
auto localCurrentRate = [&](const std::string& group_name) {
const std::vector<double>& groupSurfaceRates = group_state.injection_surface_rates(group_name);
return tcalc.calcModeRateFromRates(groupSurfaceRates);
};
const double orig_target = tcalc.groupTarget(group.injectionControls(injectionPhase, summaryState), deferred_logger);
// Assume we have a chain of groups as follows: BOTTOM -> MIDDLE -> TOP.
// Then ...
@@ -1339,6 +1363,21 @@ namespace WellGroupHelpers
local_reduction_level = ii;
}
}
// check whether guide rate is violated
if (local_reduction_level > 0) {
const auto& guided_group = chain[local_reduction_level];
const double grefficiency
= schedule.getGroup(guided_group, reportStepIdx).getGroupEfficiencyFactor();
const double currentRateFraction = grefficiency * localCurrentRate(guided_group) / (localCurrentRate(chain[local_reduction_level-1]));
const double guiderateFraction = localFraction(guided_group);
// we add a factor here to avoid switching due to numerical instability
const double factor = 1.01;
if (currentRateFraction > (guiderateFraction * factor)) {
return std::make_pair(true, guiderateFraction/currentRateFraction);
}
}
double target = orig_target;
for (size_t ii = 0; ii < num_ancestors; ++ii) {
if ((ii == 0) || guideRate->has(chain[ii], injectionPhase)) {