OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Merge pull request #2131 from jalvestad/swel-cor

Browse Source 
Correction of SWEL to handle zero rates in WCONHIST and WELOPEN SHUT
This commit is contained in:
Joakim Hove 2020-11-26 14:07:32 +01:00 committed by GitHub
commit f395629dd5
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with 69 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

29
src/opm/output/eclipse/AggregateWellData.cpp
View File

@ -517,17 +517,27 @@ namespace {
const auto& pc = well.productionControls(smry);
const auto& predMode = well.predictionMode();
if (pc.oil_rate != 0.0) {
if (predMode) {
if ((pc.oil_rate != 0.0)) {
sWell[Ix::OilRateTarget] =
swprop(M::liquid_surface_rate, pc.oil_rate);
}
if ((pc.water_rate != 0.0)) {
sWell[Ix::WatRateTarget] =
swprop(M::liquid_surface_rate, pc.water_rate);
}
if ((pc.gas_rate != 0.0)) {
sWell[Ix::GasRateTarget] =
swprop(M::gas_surface_rate, pc.gas_rate);
sWell[Ix::HistGasRateTarget] = sWell[Ix::GasRateTarget];
}
} else {
sWell[Ix::OilRateTarget] =
swprop(M::liquid_surface_rate, pc.oil_rate);
}
if (pc.water_rate != 0.0) {
sWell[Ix::WatRateTarget] =
swprop(M::liquid_surface_rate, pc.water_rate);
}
if (pc.gas_rate != 0.0) {
sWell[Ix::GasRateTarget] =
swprop(M::gas_surface_rate, pc.gas_rate);
sWell[Ix::HistGasRateTarget] = sWell[Ix::GasRateTarget];
@ -579,6 +589,11 @@ namespace {
sWell[Ix::ResVRateTarget] = getRateLimit(units, M::rate, pc.resv_rate);
//}
}
if ((well.getStatus() == Opm::Well::Status::SHUT)) {
sWell[Ix::OilRateTarget] = 0.;
sWell[Ix::WatRateTarget] = 0.;
sWell[Ix::GasRateTarget] = 0.;
}
}
else if (well.isInjector()) {
const auto& ic = well.injectionControls(smry);

48
tests/test_AggregateWellData.cpp
View File

@ -364,7 +364,22 @@ COMPDAT
WCONPROD
'OP_6' 'OPEN' 'ORAT' 20000 4* 1000 /
/
TSTEP -- 8
DATES -- 8
18 OCT 2014 /
/
WELSPECS
'OP_6' 'OP' 9 9 1* 'OIL' 1* 1* 1* 1* 1* 1* 1* /
/
COMPDAT
'OP_6' 9 9 3 9 'OPEN' 1* 32.948 0.311 3047.839 1* 1* 'X' 22.100 /
/
WCONHIST
'OP_6' 'OPEN' 'RESV' 275. 0. 34576. /
/
WELOPEN
'OP_5' 'SHUT' /
/
TSTEP -- 9
10 /
)~" };
@ -724,6 +739,37 @@ BOOST_AUTO_TEST_CASE (Declared_Well_Data)
BOOST_CHECK_EQUAL(zwell[i1 + Ix::WellName].c_str(), "OP_2 ");
}
// SWEL (OP_6)
// Report Step 8: 2014-10-18 --> 2014-10-28
const auto rptStep_8 = std::size_t{8};
const auto ih_8 = MockIH {
static_cast<int>(simCase.sched.getWells(rptStep_8).size())
};
BOOST_CHECK_EQUAL(ih_8.nwells, MockIH::Sz{6});
//smry = sim_state();
awd = Opm::RestartIO::Helpers::AggregateWellData{ih_8.value};
awd.captureDeclaredWellData(simCase.sched,
simCase.es.getUnits(), rptStep_8, action_state, smry, ih_8.value);
{
using Ix = ::Opm::RestartIO::Helpers::VectorItems::SWell::index;
const auto i1 = 4*ih_8.nswelz;
const auto& swell = awd.getSWell();
BOOST_CHECK_CLOSE(swell[i1 + Ix::OilRateTarget], 0.0f, 1.0e-7f);
BOOST_CHECK_CLOSE(swell[i1 + Ix::WatRateTarget], 0.0f, 1.0e-7f);
BOOST_CHECK_CLOSE(swell[i1 + Ix::GasRateTarget], 0.0f, 1.0e-7f);
const auto i2 = 5*ih_8.nswelz;
BOOST_CHECK_CLOSE(swell[i2 + Ix::OilRateTarget], 275.f, 1.0e-7f);
BOOST_CHECK_CLOSE(swell[i2 + Ix::WatRateTarget], 0.0f, 1.0e-7f);
BOOST_CHECK_CLOSE(swell[i2 + Ix::GasRateTarget], 34576.0f, 1.0e-7f);
}
}