OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-07-07 04:53:03 -05:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #5409 from bska/collect-well-limits-for-smry-output

Browse Source 
Collect Well Level Control Limits/Targets for Summary Output
This commit is contained in:
Bård Skaflestad 2024-06-06 14:51:22 +02:00 committed by GitHub
commit e876fa86e9
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
3 changed files with 83 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/simulators/wells/BlackoilWellModel.hpp
View File

@ -264,11 +264,17 @@ template<class Scalar> class WellContributions;
return this->wasDynamicallyShutThisTimeStep(well_index);
});
const auto& tracerRates = simulator_.problem().tracerModel().getWellTracerRates();
this->assignWellTracerRates(wsrpt, tracerRates);
{
const auto& tracerRates = this->simulator_.problem()
.tracerModel().getWellTracerRates();
this->assignWellTracerRates(wsrpt, tracerRates);
}
BlackoilWellModelGuideRates(*this).assignWellGuideRates(wsrpt, this->reportStepIndex());
BlackoilWellModelGuideRates(*this)
.assignWellGuideRates(wsrpt, this->reportStepIndex());
this->assignWellTargets(wsrpt);
this->assignShutConnections(wsrpt, this->reportStepIndex());
return wsrpt;

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

@ -43,6 +43,7 @@
#include <opm/input/eclipse/Schedule/Network/Balance.hpp>
#include <opm/input/eclipse/Schedule/Network/ExtNetwork.hpp>
#include <opm/input/eclipse/Schedule/Schedule.hpp>
#include <opm/input/eclipse/Schedule/ScheduleTypes.hpp>
#include <opm/input/eclipse/Schedule/Well/WellConnections.hpp>
#include <opm/input/eclipse/Schedule/Well/WellTestConfig.hpp>
@ -890,6 +891,74 @@ setWsolvent(const Group& group,
}
}
template<class Scalar>
void BlackoilWellModelGeneric<Scalar>::
assignWellTargets(data::Wells& wsrpt) const
{
auto pwInfo = this->local_parallel_well_info_.begin();
for (const auto& well : this->wells_ecl_) {
if (! pwInfo++->get().isOwner()) {
continue;
}
// data::Wells is a std::map<>
auto& limits = wsrpt[well.name()].limits;
if (well.isProducer()) {
this->assignProductionWellTargets(well, limits);
}
else {
this->assignInjectionWellTargets(well, limits);
}
}
}
template<class Scalar>
void BlackoilWellModelGeneric<Scalar>::
assignProductionWellTargets(const Well& well, data::WellControlLimits& limits) const
{
using Item = data::WellControlLimits::Item;
const auto ctrl = well.productionControls(this->summaryState());
limits
.set(Item::Bhp, ctrl.bhp_limit)
.set(Item::OilRate, ctrl.oil_rate)
.set(Item::WaterRate, ctrl.water_rate)
.set(Item::GasRate, ctrl.gas_rate)
.set(Item::ResVRate, ctrl.resv_rate)
.set(Item::LiquidRate, ctrl.liquid_rate);
}
template<class Scalar>
void BlackoilWellModelGeneric<Scalar>::
assignInjectionWellTargets(const Well& well, data::WellControlLimits& limits) const
{
using Item = data::WellControlLimits::Item;
const auto ctrl = well.injectionControls(this->summaryState());
limits
.set(Item::Bhp, ctrl.bhp_limit)
.set(Item::ResVRate, ctrl.reservoir_rate);
if (ctrl.injector_type == InjectorType::MULTI) {
// Not supported
return;
}
auto rateItem = Item::WaterRate;
if (ctrl.injector_type == InjectorType::GAS) {
rateItem = Item::GasRate;
}
else if (ctrl.injector_type == InjectorType::OIL) {
rateItem = Item::OilRate;
}
limits.set(rateItem, ctrl.surface_rate);
}
template<class Scalar>
void BlackoilWellModelGeneric<Scalar>::
assignShutConnections(data::Wells& wsrpt,
@ -901,7 +970,7 @@ assignShutConnections(data::Wells& wsrpt,
auto& xwel = wsrpt[well.name()]; // data::Wells is a std::map<>
xwel.dynamicStatus = this->schedule()
.getWell(well.name(), reportStepIndex).getStatus();
.getWell(well.name(), reportStepIndex).getStatus();
const auto wellIsOpen = xwel.dynamicStatus == Well::Status::OPEN;
auto skip = [wellIsOpen](const Connection& conn)
@ -913,7 +982,7 @@ assignShutConnections(data::Wells& wsrpt,
!this->wasDynamicallyShutThisTimeStep(wellID))
{
xwel.dynamicStatus = well.getAutomaticShutIn()
? Well::Status::SHUT : Well::Status::STOP;
? Well::Status::SHUT : Well::Status::STOP;
}
auto& xcon = xwel.connections;

3
opm/simulators/wells/BlackoilWellModelGeneric.hpp
View File

@ -349,6 +349,9 @@ protected:
void assignShutConnections(data::Wells& wsrpt,
const int reportStepIndex) const;
void assignWellTargets(data::Wells& wsrpt) const;
void assignProductionWellTargets(const Well& well, data::WellControlLimits& limits) const;
void assignInjectionWellTargets(const Well& well, data::WellControlLimits& limits) const;
void assignGroupControl(const Group& group,
data::GroupData& gdata) const;
void assignGroupValues(const int reportStepIdx,