OilfieldToolsNet/opm-simulators
4
0 Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2026-10-26 23:50:08 -05:00
Code Issues Packages Projects Releases Wiki Activity

Set thp from network in well testing

Browse Source
This commit is contained in:
Tor Harald Sandve committed 2025-03-31 14:12:29 +02:00
1 parent 563481ce03
commit c0431ce4f9
2 files changed
+22 -13
Download Patch File Download Diff File Expand all files Collapse all files

No files matched your search

opm/simulators/wells/BlackoilWellModel_impl.hpp
+12
View File
@@ -565,6 +565,18 @@ namespace Opm {
well->setPrevSurfaceRates(this->wellState(), this->prevWellState());
}
const auto& network = this->schedule()[timeStepIdx].network();
if (network.active() && !this->node_pressures_.empty()) {
if (well->isProducer()) {
const auto it = this->node_pressures_.find(well->wellEcl().groupName());
if (it != this->node_pressures_.end()) {
// The well belongs to a group which has a network nodal pressure,
// set the dynamic THP constraint based on the network nodal pressure
const Scalar nodal_pressure = it->second;
well->setDynamicThpLimit(nodal_pressure);
}
}
}
try {
using GLiftEclWells = typename GasLiftGroupInfo<Scalar>::GLiftEclWells;
GLiftEclWells ecl_well_map;
opm/simulators/wells/WellInterface_impl.hpp
+10 -13
View File
@@ -397,6 +397,15 @@ namespace Opm
WellState<Scalar> well_state_copy = well_state;
auto& ws = well_state_copy.well(this->indexOfWell());
const auto& summary_state = simulator.vanguard().summaryState();
const bool has_thp_limit = this->wellHasTHPConstraints(summary_state);
if (has_thp_limit) {
ws.production_cmode = Well::ProducerCMode::THP;
}
else {
ws.production_cmode = Well::ProducerCMode::BHP;
}
updateWellStateWithTarget(simulator, group_state, well_state_copy, deferred_logger);
calculateExplicitQuantities(simulator, well_state_copy, deferred_logger);
updatePrimaryVariables(simulator, well_state_copy, deferred_logger);
@@ -715,19 +724,8 @@ namespace Opm
DeferredLogger& deferred_logger)
{
OPM_TIMEFUNCTION();
// keep a copy of the original well state
const WellState<Scalar> well_state0 = well_state;
const double dt = simulator.timeStepSize();
const auto& summary_state = simulator.vanguard().summaryState();
const bool has_thp_limit = this->wellHasTHPConstraints(summary_state);
bool converged;
if (has_thp_limit) {
well_state.well(this->indexOfWell()).production_cmode = Well::ProducerCMode::THP;
}
else {
well_state.well(this->indexOfWell()).production_cmode = Well::ProducerCMode::BHP;
}
converged = iterateWellEquations(simulator, dt, well_state, group_state, deferred_logger);
const bool converged = iterateWellEquations(simulator, dt, well_state, group_state, deferred_logger);
if (converged) {
deferred_logger.debug("WellTest: Well equation for well " + this->name() + " converged");
return true;
@@ -735,7 +733,6 @@ namespace Opm
const int max_iter = this->param_.max_welleq_iter_;
deferred_logger.debug("WellTest: Well equation for well " + this->name() + " failed converging in "
+ std::to_string(max_iter) + " iterations");
well_state = well_state0;
return false;
}