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

Merge pull request #3526 from plgbrts/well-control

Browse Source 
Allow NONE for item 7 in GCONPROD
This commit is contained in:
Tor Harald Sandve
2023-06-09 10:50:57 +02:00
committed by GitHub
parent 78c527fb53 db4c2d7f38
commit ba0845c5f3
4 changed files with 90 additions and 17 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
src/opm/input/eclipse/Schedule/Group/Group.cpp
View File

@@ -130,6 +130,7 @@ namespace {
production.exceed_action = Opm::Group::ExceedActionFromInt(rst_group.exceed_action);
production.guide_rate_def = Opm::Group::GuideRateProdTargetFromInt(rst_group.prod_guide_rate_def);
//TODO allow also for ExceedAction::NONE (item 7 of GCONPROD)
if ((production.cmode == Opm::Group::ProductionCMode::ORAT) ||
(production.cmode == Opm::Group::ProductionCMode::WRAT) ||
(production.cmode == Opm::Group::ProductionCMode::GRAT) ||

74
src/opm/input/eclipse/Schedule/KeywordHandlers.cpp
View File

@@ -572,27 +572,69 @@ File {} line {}.)", wname, location.keyword, location.filename, location.lineno)
production.resv_target = resv_target;
production.available_group_control = availableForGroupControl;
if ((production.cmode == Group::ProductionCMode::ORAT) ||
(production.cmode == Group::ProductionCMode::WRAT) ||
(production.cmode == Group::ProductionCMode::GRAT) ||
(production.cmode == Group::ProductionCMode::LRAT))
production.exceed_action = Group::ExceedAction::RATE;
else
production.exceed_action = exceedAction;
production.exceed_action = exceedAction;
production.production_controls = 0;
if (!apply_default_oil_target)
// GCONPROD
// 'G1' 'ORAT' 1000 100 200 300 NONE => constraints 100,200,300 should be ignored
//
// GCONPROD
// 'G1' 'ORAT' 1000 100 200 300 RATE => constraints 100,200,300 should be honored
if (production.cmode == Group::ProductionCMode::ORAT){
production.production_controls += static_cast<int>(Group::ProductionCMode::ORAT);
if (!apply_default_gas_target)
if (exceedAction == Group::ExceedAction::RATE) {
if (!apply_default_gas_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::GRAT);
if (!apply_default_water_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::WRAT);
if (!apply_default_liquid_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::LRAT);
}
}
else if (production.cmode == Group::ProductionCMode::GRAT){
production.production_controls += static_cast<int>(Group::ProductionCMode::GRAT);
if (!apply_default_water_target)
if (exceedAction == Group::ExceedAction::RATE) {
if (!apply_default_oil_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::ORAT);
if (!apply_default_water_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::WRAT);
if (!apply_default_liquid_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::LRAT);
}
}
else if (production.cmode == Group::ProductionCMode::WRAT){
production.production_controls += static_cast<int>(Group::ProductionCMode::WRAT);
if (!apply_default_liquid_target)
if (exceedAction == Group::ExceedAction::RATE) {
if (!apply_default_oil_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::ORAT);
if (!apply_default_gas_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::GRAT);
if (!apply_default_liquid_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::LRAT);
}
}
else if (production.cmode == Group::ProductionCMode::LRAT){
production.production_controls += static_cast<int>(Group::ProductionCMode::LRAT);
if (exceedAction == Group::ExceedAction::RATE) {
if (!apply_default_oil_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::ORAT);
if (!apply_default_gas_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::GRAT);
if (!apply_default_water_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::WRAT);
}
}
else {
if (exceedAction == Group::ExceedAction::RATE) {
if (!apply_default_oil_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::ORAT);
if (!apply_default_water_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::WRAT);
if (!apply_default_gas_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::GRAT);
if (!apply_default_liquid_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::LRAT);
}
}
if (!apply_default_resv_target)
production.production_controls += static_cast<int>(Group::ProductionCMode::RESV);

7
tests/parser/GroupTests.cpp
View File

@@ -259,11 +259,13 @@ BOOST_AUTO_TEST_CASE(createDeckWithGCONPROD) {
GRUPTREE
'G1' 'FIELD' /
'G2' 'FIELD' /
'G3' 'FIELD' /
/
GCONPROD
'G1' 'ORAT' 10000 3* 'CON' /
'G1' 'ORAT' 10000 3* 'RATE' /
'G2' 'RESV' 10000 3* 'CON' /
'G3' 'ORAT' 10000 3* 1* /
/)";
auto schedule = create_schedule(input);
@@ -271,12 +273,15 @@ BOOST_AUTO_TEST_CASE(createDeckWithGCONPROD) {
const auto& group1 = schedule.getGroup("G1", 0);
const auto& group2 = schedule.getGroup("G2", 0);
const auto& group3 = schedule.getGroup("G3", 0);
auto ctrl1 = group1.productionControls(st);
auto ctrl2 = group2.productionControls(st);
auto ctrl3 = group3.productionControls(st);
BOOST_CHECK(ctrl1.exceed_action == Group::ExceedAction::RATE);
BOOST_CHECK(ctrl2.exceed_action == Group::ExceedAction::CON);
BOOST_CHECK(ctrl3.exceed_action == Group::ExceedAction::NONE);
}

25
tests/parser/ScheduleTests.cpp
View File

@@ -1509,6 +1509,18 @@ BOOST_AUTO_TEST_CASE(createDeckModifyMultipleGCONPROD) {
GCONPROD
'G*' 'ORAT' 2000 0 0 0 'NONE' 'YES' 148 'OIL'/
/
DATES -- 3
10 DEC 2008 /
/
GCONPROD
'G*' 'ORAT' 2000 1000 0 0 'NONE' 'YES' 148 'OIL'/
/
DATES -- 4
10 JAN 2009 /
/
GCONPROD
'G*' 'ORAT' 2000 1000 0 0 'RATE' 'YES' 148 'OIL'/
/
)";
const auto& schedule = make_schedule(input);
@@ -1530,6 +1542,19 @@ BOOST_AUTO_TEST_CASE(createDeckModifyMultipleGCONPROD) {
BOOST_CHECK_EQUAL(g.productionControls(st).guide_rate, 148);
BOOST_CHECK_EQUAL(true, g.productionControls(st).guide_rate_def == Group::GuideRateProdTarget::OIL);
}
{
auto g = schedule.getGroup("G1", 3);
BOOST_CHECK_CLOSE(g.productionControls(st).oil_target, 2000 * siFactorL, 1e-13);
BOOST_CHECK(g.has_control(Group::ProductionCMode::ORAT));
BOOST_CHECK(!g.has_control(Group::ProductionCMode::WRAT));
}
{
auto g = schedule.getGroup("G1", 4);
BOOST_CHECK_CLOSE(g.productionControls(st).oil_target, 2000 * siFactorL, 1e-13);
BOOST_CHECK(g.has_control(Group::ProductionCMode::ORAT));
BOOST_CHECK_CLOSE(g.productionControls(st).water_target, 1000 * siFactorL, 1e-13);
BOOST_CHECK(g.has_control(Group::ProductionCMode::WRAT));
}
auto g2 = schedule.getGroup("G2", 2);
BOOST_CHECK_CLOSE(g2.productionControls(st).oil_target, 2000 * siFactorL, 1e-13);