308 lines
10 KiB
C++
308 lines
10 KiB
C++
/*
|
|
Copyright 2013 Statoil ASA.
|
|
|
|
This file is part of the Open Porous Media project (OPM).
|
|
|
|
OPM is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OPM is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
#include <stdexcept>
|
|
#include <iostream>
|
|
#include <boost/filesystem.hpp>
|
|
|
|
#define BOOST_TEST_MODULE GroupTests
|
|
#include <boost/test/unit_test.hpp>
|
|
#include <boost/date_time/posix_time/posix_time.hpp>
|
|
|
|
#include <opm/parser/eclipse/Deck/Deck.hpp>
|
|
#include <opm/parser/eclipse/Parser/Parser.hpp>
|
|
#include <opm/parser/eclipse/Parser/ParseContext.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/Util/Value.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/Grid/EclipseGrid.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/Schedule/Schedule.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/Schedule/Group/Group2.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/Schedule/Group/GuideRateModel.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/Schedule/SummaryState.hpp>
|
|
|
|
#include <opm/parser/eclipse/EclipseState/Schedule/Well/WellProductionProperties.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/Schedule/Well/WellInjectionProperties.hpp>
|
|
|
|
using namespace Opm;
|
|
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(CreateGroup_CorrectNameAndDefaultValues) {
|
|
Opm::Group2 group("G1" , 1, 0, 0, UnitSystem::newMETRIC());
|
|
BOOST_CHECK_EQUAL( "G1" , group.name() );
|
|
}
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(CreateGroupCreateTimeOK) {
|
|
Opm::Group2 group("G1" , 1, 5, 0, UnitSystem::newMETRIC());
|
|
BOOST_CHECK_EQUAL( false, group.defined( 4 ));
|
|
BOOST_CHECK_EQUAL( true, group.defined( 5 ));
|
|
BOOST_CHECK_EQUAL( true, group.defined( 6 ));
|
|
}
|
|
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(CreateGroup_SetInjectorProducer_CorrectStatusSet) {
|
|
Opm::Group2 group1("IGROUP" , 1, 0, 0, UnitSystem::newMETRIC());
|
|
Opm::Group2 group2("PGROUP" , 2, 0, 0, UnitSystem::newMETRIC());
|
|
|
|
group1.setProductionGroup();
|
|
BOOST_CHECK(group1.isProductionGroup());
|
|
BOOST_CHECK(!group1.isInjectionGroup());
|
|
|
|
group2.setInjectionGroup();
|
|
BOOST_CHECK(!group2.isProductionGroup());
|
|
BOOST_CHECK(group2.isInjectionGroup());
|
|
}
|
|
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(ControlModeOK) {
|
|
Opm::Group2 group("G1" , 1, 0, 0, UnitSystem::newMETRIC());
|
|
Opm::SummaryState st;
|
|
const auto& prod = group.productionControls(st);
|
|
BOOST_CHECK_EQUAL( Opm::GroupInjection::NONE , prod.cmode);
|
|
}
|
|
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(GroupChangePhaseSameTimeThrows) {
|
|
Opm::Group2 group("G1" , 1, 0, 0, UnitSystem::newMETRIC());
|
|
Opm::SummaryState st;
|
|
const auto& inj = group.injectionControls(st);
|
|
BOOST_CHECK_EQUAL( Opm::Phase::WATER , inj.phase); // Default phase - assumed WATER
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(GroupDoesNotHaveWell) {
|
|
Opm::Group2 group("G1" , 1, 0, 0, UnitSystem::newMETRIC());
|
|
|
|
BOOST_CHECK_EQUAL(false , group.hasWell("NO"));
|
|
BOOST_CHECK_EQUAL(0U , group.numWells());
|
|
}
|
|
|
|
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(createDeckWithGEFAC) {
|
|
Opm::Parser parser;
|
|
std::string input =
|
|
"START -- 0 \n"
|
|
"19 JUN 2007 / \n"
|
|
"SCHEDULE\n"
|
|
|
|
"WELSPECS\n"
|
|
" 'B-37T2' 'PRODUC' 9 9 1* 'OIL' 1* 1* 1* 1* 1* 1* 1* / \n"
|
|
" 'B-43A' 'PRODUC' 8 8 1* 'OIL' 1* 1* 1* 1* 1* 1* 1* / \n"
|
|
"/\n"
|
|
|
|
"COMPDAT\n"
|
|
" 'B-37T2' 9 9 1 1 'OPEN' 1* 32.948 0.311 3047.839 1* 1* 'X' 22.100 / \n"
|
|
" 'B-43A' 8 8 2 2 'OPEN' 1* 46.825 0.311 4332.346 1* 1* 'X' 22.123 / \n"
|
|
"/\n"
|
|
|
|
"GEFAC\n"
|
|
" 'PRODUC' 0.85 / \n"
|
|
"/\n";
|
|
|
|
auto deck = parser.parseString(input);
|
|
EclipseGrid grid(10,10,10);
|
|
TableManager table ( deck );
|
|
Eclipse3DProperties eclipseProperties ( deck , table, grid);
|
|
Runspec runspec (deck );
|
|
Opm::Schedule schedule(deck, grid, eclipseProperties, runspec);
|
|
|
|
auto group_names = schedule.groupNames("PRODUC");
|
|
BOOST_CHECK_EQUAL(group_names.size(), 1);
|
|
BOOST_CHECK_EQUAL(group_names[0], "PRODUC");
|
|
|
|
const auto& group1 = schedule.getGroup2("PRODUC", 0);
|
|
BOOST_CHECK_EQUAL(group1.getGroupEfficiencyFactor(), 0.85);
|
|
BOOST_CHECK(group1.getTransferGroupEfficiencyFactor());
|
|
}
|
|
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(createDeckWithWGRUPCONandWCONPROD) {
|
|
|
|
/* Test deck with well guide rates for group control:
|
|
GRUPCON (well guide rates for group control)
|
|
WCONPROD (conrol data for production wells) with GRUP control mode */
|
|
|
|
Opm::Parser parser;
|
|
std::string input =
|
|
"START -- 0 \n"
|
|
"19 JUN 2007 / \n"
|
|
"SCHEDULE\n"
|
|
|
|
"WELSPECS\n"
|
|
" 'B-37T2' 'PRODUC' 9 9 1* 'OIL' 1* 1* 1* 1* 1* 1* 1* / \n"
|
|
" 'B-43A' 'PRODUC' 8 8 1* 'OIL' 1* 1* 1* 1* 1* 1* 1* / \n"
|
|
"/\n"
|
|
|
|
"COMPDAT\n"
|
|
" 'B-37T2' 9 9 1 1 'OPEN' 1* 32.948 0.311 3047.839 1* 1* 'X' 22.100 / \n"
|
|
" 'B-43A' 8 8 2 2 'OPEN' 1* 46.825 0.311 4332.346 1* 1* 'X' 22.123 / \n"
|
|
"/\n"
|
|
|
|
|
|
"WGRUPCON\n"
|
|
" 'B-37T2' YES 30 OIL / \n"
|
|
" 'B-43A' YES 30 OIL / \n"
|
|
"/\n"
|
|
|
|
"WCONPROD\n"
|
|
" 'B-37T2' 'OPEN' 'GRUP' 1000 2* 2000.000 2* 1* 10 200000.000 5* / / \n"
|
|
" 'B-43A' 'OPEN' 'GRUP' 1200 2* 3000.000 2* 1* 11 0.000 5* / / \n"
|
|
"/\n";
|
|
|
|
|
|
|
|
auto deck = parser.parseString(input);
|
|
EclipseGrid grid(10,10,10);
|
|
TableManager table ( deck );
|
|
Eclipse3DProperties eclipseProperties ( deck , table, grid);
|
|
Runspec runspec (deck );
|
|
Opm::Schedule schedule(deck, grid, eclipseProperties, runspec);
|
|
const auto& currentWell = schedule.getWell2("B-37T2", 0);
|
|
const Opm::Well2::WellProductionProperties& wellProductionProperties = currentWell.getProductionProperties();
|
|
BOOST_CHECK(wellProductionProperties.controlMode == Opm::Well2::ProducerCMode::GRUP);
|
|
|
|
BOOST_CHECK_EQUAL(currentWell.isAvailableForGroupControl(), true);
|
|
BOOST_CHECK_EQUAL(currentWell.getGuideRate(), 30);
|
|
BOOST_CHECK_EQUAL(currentWell.getGuideRatePhase(), Opm::GuideRate::OIL);
|
|
BOOST_CHECK_EQUAL(currentWell.getGuideRateScalingFactor(), 1.0);
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(createDeckWithGRUPNET) {
|
|
Opm::Parser parser;
|
|
std::string input =
|
|
"START -- 0 \n"
|
|
"31 AUG 1993 / \n"
|
|
"SCHEDULE\n"
|
|
|
|
"GRUPNET \n"
|
|
" 'FIELD' 20.000 5* / \n"
|
|
" 'PROD' 20.000 5* / \n"
|
|
" 'MANI-B2' 1* 8 1* 'NO' 2* / \n"
|
|
" 'MANI-B1' 1* 8 1* 'NO' 2* / \n"
|
|
" 'MANI-K1' 1* 9999 4* / \n"
|
|
" 'B1-DUMMY' 1* 9999 4* / \n"
|
|
" 'MANI-D1' 1* 8 1* 'NO' 2* / \n"
|
|
" 'MANI-D2' 1* 8 1* 'NO' 2* / \n"
|
|
" 'MANI-K2' 1* 9999 4* / \n"
|
|
" 'D2-DUMMY' 1* 9999 4* / \n"
|
|
" 'MANI-E1' 1* 9 1* 'NO' 2* / \n"
|
|
" 'MANI-E2' 1* 9 4* / \n"
|
|
"/\n";
|
|
|
|
auto deck = parser.parseString(input);
|
|
EclipseGrid grid(10,10,10);
|
|
TableManager table ( deck );
|
|
Eclipse3DProperties eclipseProperties ( deck , table, grid);
|
|
Runspec runspec (deck );
|
|
Opm::Schedule schedule(deck, grid, eclipseProperties, runspec);
|
|
|
|
const auto& group1 = schedule.getGroup2("PROD", 0);
|
|
const auto& group2 = schedule.getGroup2("MANI-E2", 0);
|
|
const auto& group3 = schedule.getGroup2("MANI-K1", 0);
|
|
BOOST_CHECK_EQUAL(group1.getGroupNetVFPTable(), 0);
|
|
BOOST_CHECK_EQUAL(group2.getGroupNetVFPTable(), 9);
|
|
BOOST_CHECK_EQUAL(group3.getGroupNetVFPTable(), 9999);
|
|
}
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(Group2Create) {
|
|
Opm::Group2 g1("NAME", 1, 1, 0, UnitSystem::newMETRIC());
|
|
Opm::Group2 g2("NAME", 1, 1, 0, UnitSystem::newMETRIC());
|
|
|
|
BOOST_CHECK( g1.addWell("W1") );
|
|
BOOST_CHECK( !g1.addWell("W1") );
|
|
BOOST_CHECK( g1.addWell("W2") );
|
|
BOOST_CHECK( g1.hasWell("W1"));
|
|
BOOST_CHECK( g1.hasWell("W2"));
|
|
BOOST_CHECK( !g1.hasWell("W3"));
|
|
BOOST_CHECK_EQUAL( g1.numWells(), 2);
|
|
BOOST_CHECK_THROW(g1.delWell("W3"), std::invalid_argument);
|
|
BOOST_CHECK_NO_THROW(g1.delWell("W1"));
|
|
BOOST_CHECK_EQUAL( g1.numWells(), 1);
|
|
|
|
|
|
BOOST_CHECK( g2.addGroup("G1") );
|
|
BOOST_CHECK( !g2.addGroup("G1") );
|
|
BOOST_CHECK( g2.addGroup("G2") );
|
|
|
|
// The children must be either all wells - or all groups.
|
|
BOOST_CHECK_THROW(g1.addGroup("G1"), std::logic_error);
|
|
BOOST_CHECK_THROW(g2.addWell("W1"), std::logic_error);
|
|
}
|
|
|
|
BOOST_AUTO_TEST_CASE(createDeckWithGCONPROD) {
|
|
Opm::Parser parser;
|
|
std::string input = R"(
|
|
START -- 0
|
|
31 AUG 1993 /
|
|
SCHEDULE
|
|
|
|
GRUPTREE
|
|
'G1' 'FIELD' /
|
|
'G2' 'FIELD' /
|
|
/
|
|
|
|
GCONPROD
|
|
'G1' 'ORAT' 10000 3* 'CON' /
|
|
'G2' 'RESV' 10000 3* 'CON' /
|
|
/)";
|
|
|
|
auto deck = parser.parseString(input);
|
|
EclipseGrid grid(10,10,10);
|
|
TableManager table ( deck );
|
|
Eclipse3DProperties eclipseProperties ( deck , table, grid);
|
|
Runspec runspec (deck );
|
|
Opm::Schedule schedule(deck, grid, eclipseProperties, runspec);
|
|
SummaryState st;
|
|
|
|
const auto& group1 = schedule.getGroup2("G1", 0);
|
|
const auto& group2 = schedule.getGroup2("G2", 0);
|
|
|
|
auto ctrl1 = group1.productionControls(st);
|
|
auto ctrl2 = group2.productionControls(st);
|
|
|
|
BOOST_CHECK(ctrl1.exceed_action == Group2::ExceedAction::RATE);
|
|
BOOST_CHECK(ctrl2.exceed_action == Group2::ExceedAction::CON);
|
|
}
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE(TESTGuideRateModel) {
|
|
Opm::GuideRateModel grc_default;
|
|
BOOST_CHECK_THROW(Opm::GuideRateModel(0.0,GuideRateTarget::NONE, -5,0,0,0,0,0,true,1,true), std::invalid_argument);
|
|
BOOST_CHECK_THROW(grc_default.eval(1,0.50,0.50), std::invalid_argument);
|
|
|
|
Opm::GuideRateModel grc_delay(10, GuideRateTarget::NONE, 1,1,0,0,0,0,true,1,true);
|
|
BOOST_CHECK_NO_THROW(grc_delay.eval(1.0, 0.5, 0.5));
|
|
}
|