From 98450b5aa843eaa344dff629ce222d5b1388abc4 Mon Sep 17 00:00:00 2001 From: Kai Bao Date: Fri, 4 Jan 2019 15:10:18 +0100 Subject: [PATCH] fixing the way in handling the RESV rates for WCONHIST --- opm/autodiff/BlackoilWellModel_impl.hpp | 32 ++++++++++++++----------- 1 file changed, 18 insertions(+), 14 deletions(-) diff --git a/opm/autodiff/BlackoilWellModel_impl.hpp b/opm/autodiff/BlackoilWellModel_impl.hpp index de999dfd8..63f65b027 100644 --- a/opm/autodiff/BlackoilWellModel_impl.hpp +++ b/opm/autodiff/BlackoilWellModel_impl.hpp @@ -1536,22 +1536,26 @@ namespace Opm { well_controls_iset_distr(ctrl, rctrl, & distr[0]); - // for the WCONHIST wells, we need to calculate the rates - const WellMap::const_iterator i = wmap.find(wells()->name[*rp]); + // for the WCONHIST wells, we need to calculate the RESV rates since it can not be specified directly + if (is_producer) { + const WellMap::const_iterator i = wmap.find(wells()->name[*rp]); - if (i == wmap.end()) { - OPM_THROW(std::runtime_error, "Failed to find the well " << wells()->name[*rp] << " in wmap."); + if (i == wmap.end()) { + OPM_THROW(std::runtime_error, "Failed to find the well " << wells()->name[*rp] << " in wmap."); + } + const auto* wp = i->second; + const WellProductionProperties& production_properties = wp->getProductionProperties(step); + // historical phase rates + std::vector hrates(np); + SimFIBODetails::historyRates(phase_usage_, production_properties, hrates); + + std::vector hrates_resv(np); + rateConverter_->calcReservoirVoidageRates(fipreg, pvtreg, hrates, hrates_resv); + + const double target = -std::accumulate(hrates_resv.begin(), hrates_resv.end(), 0.0); + + well_controls_iset_target(ctrl, rctrl, target); } - // not handling injector for now - const auto* wp = i->second; - const WellProductionProperties& production_properties = wp->getProductionProperties(step); - // historical phase rates - std::vector hrates(np); - SimFIBODetails::historyRates(phase_usage_, production_properties, hrates); - - const double target = - std::inner_product(distr.begin(), distr.end(), - hrates.begin(), 0.0); - well_controls_iset_target(ctrl, rctrl, target); } } }