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

addressing reviewing comments for PR#3760

Browse Source
This commit is contained in:
Kai Bao 2022-01-10 11:33:11 +01:00
parent dec3100247
commit 2b9b6e5155
1 changed files with 17 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
opm/simulators/wells/MultisegmentWellGeneric.cpp
View File

@ -483,9 +483,12 @@ computeBhpAtThpLimitProd(const std::function<std::vector<double>(const double)>&
auto fflo = [&flo, &frates](double bhp) { return flo(frates(bhp)); };
auto bhp_max = this->bhpMax(fflo, controls.bhp_limit, maxPerfPress, table.getFloAxis().front(), deferred_logger);
// could not solve for the bhp-point, we could not contrinue to find the bhp
// could not solve for the bhp-point, we could not continue to find the bhp
if (!bhp_max.has_value()) {
return {};
deferred_logger.warning("FAILED_ROBUST_BHP_THP_SOLVE_INOPERABLE",
"Robust bhp(thp) solve failed due to not being able to "
"find bhp-point where production becomes non-zero for well " + baseif_.name());
return std::nullopt;
}
// Define the equation we want to solve.
auto eq = [&fbhp, &frates](double bhp) {
@ -506,11 +509,14 @@ computeBhpAtThpLimitProd(const std::function<std::vector<double>(const double)>&
}
if (!finding_bracket) {
deferred_logger.debug(" trying the brute force way for last attempt ");
deferred_logger.debug(" Trying the brute force search to bracket the bhp for last attempt ");
finding_bracket = this->bruteForceBracket(eq, range, low, high, deferred_logger);
}
if (!finding_bracket) {
deferred_logger.warning("FAILED_ROBUST_BHP_THP_SOLVE_INOPERABLE",
"Robust bhp(thp) solve failed due to not being able to "
"bracket the bhp solution with the brute force search for " + baseif_.name());
return std::nullopt;
}
@ -566,7 +572,7 @@ bhpMax(const std::function<double(const double)>& fflo,
// empty optional.
deferred_logger.warning("FAILED_ROBUST_BHP_THP_SOLVE_INOPERABLE",
"Robust bhp(thp) solve failed due to inoperability for well " + baseif_.name());
return {};
return std::nullopt;
} else {
// Still producing, even at high bhp.
assert(f_high < 0.0);
@ -590,7 +596,13 @@ bhpMax(const std::function<double(const double)>& fflo,
}
++it;
}
bhp_max = low;
if (it < maxit) {
bhp_max = low;
} else {
deferred_logger.warning("FAILED_ROBUST_BHP_THP_SOLVE_INOPERABLE",
"Bisect did not find the bhp-point where we produce for well " + baseif_.name());
return std::nullopt;
}
}
deferred_logger.debug("computeBhpAtThpLimitProd(): well = " + baseif_.name() +
" low = " + std::to_string(low) +