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

Merge pull request #1819 from joakim-hove/aicd

Browse Source 
Aicd
This commit is contained in:
Joakim Hove
2020-06-12 09:12:18 +02:00
committed by GitHub
parent e08a3666e3 224911bd10
commit b07bd377cd
9 changed files with 382 additions and 71 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
CMakeLists_files.cmake
View File

@@ -112,6 +112,7 @@ if(ENABLE_ECL_INPUT)
src/opm/parser/eclipse/EclipseState/Schedule/MSW/Compsegs.cpp
src/opm/parser/eclipse/EclipseState/Schedule/MSW/Segment.cpp
src/opm/parser/eclipse/EclipseState/Schedule/MSW/WellSegments.cpp
src/opm/parser/eclipse/EclipseState/Schedule/MSW/AICD.cpp
src/opm/parser/eclipse/EclipseState/Schedule/MSW/SICD.cpp
src/opm/parser/eclipse/EclipseState/Schedule/MSW/Valve.cpp
src/opm/parser/eclipse/EclipseState/Schedule/Network/Branch.cpp
@@ -715,6 +716,7 @@ if(ENABLE_ECL_INPUT)
opm/parser/eclipse/EclipseState/Schedule/MSW/Segment.hpp
opm/parser/eclipse/EclipseState/Schedule/MSW/Segment.hpp
opm/parser/eclipse/EclipseState/Schedule/MSW/WellSegments.hpp
opm/parser/eclipse/EclipseState/Schedule/MSW/AICD.hpp
opm/parser/eclipse/EclipseState/Schedule/MSW/SICD.hpp
opm/parser/eclipse/EclipseState/SimulationConfig/ThresholdPressure.hpp
opm/parser/eclipse/EclipseState/SimulationConfig/BCConfig.hpp

69
opm/parser/eclipse/EclipseState/Schedule/MSW/AICD.hpp Normal file
View File

@@ -0,0 +1,69 @@
/*
Copyright 2020 Equinor 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/>.
*/
#ifndef AICD_HPP_HEADER_INCLUDED
#define AICD_HPP_HEADER_INCLUDED
#include <map>
#include <utility>
#include <vector>
#include <string>
#include <opm/parser/eclipse/EclipseState/Schedule/MSW/SICD.hpp>
namespace Opm {
class DeckRecord;
class DeckKeyword;
class AutoICD : public SICD {
public:
AutoICD() = default;
AutoICD(const DeckRecord& record);
static AutoICD serializeObject();
// the function will return a map
// [
// "WELL1" : [<seg1, aicd1>, <seg2, aicd2> ...]
// ....
static std::map<std::string, std::vector<std::pair<int, AutoICD> > >
fromWSEGAICD(const DeckKeyword& wsegaicd);
bool operator==(const AutoICD& data) const;
template<class Serializer>
void serializeOp(Serializer& serializer)
{
AutoICD::serializeOp(serializer);
}
private:
double m_flow_rate_exponent;
double m_visc_exponent;
double m_oil_density_exponent;
double m_water_density_exponent;
double m_gas_density_exponent;
double m_oil_viscosity_exponent;
double m_water_viscosity_exponent;
double m_gas_viscosity_exponent;
};
}
#endif

32
opm/parser/eclipse/EclipseState/Schedule/MSW/SICD.hpp
View File

@@ -27,12 +27,11 @@
#include <string>
#include <opm/parser/eclipse/EclipseState/Schedule/MSW/icd.hpp>
#include <opm/parser/eclipse/Deck/DeckRecord.hpp>
#include <opm/parser/eclipse/Deck/DeckKeyword.hpp>
namespace Opm {
class DeckRecord;
class DeckKeyword;
class SICD {
public:
@@ -91,6 +90,33 @@ namespace Opm {
serializer(m_scaling_factor);
}
template<class ICD>
static std::map<std::string, std::vector<std::pair<int, ICD> > >
fromWSEG(const DeckKeyword& wseg) {
std::map<std::string, std::vector<std::pair<int, ICD> > > res;
for (const DeckRecord &record : wseg) {
const std::string well_name = record.getItem("WELL").getTrimmedString(0);
const int start_segment = record.getItem("SEG1").get<int>(0);
const int end_segment = record.getItem("SEG2").get<int>(0);
if (start_segment < 2 || end_segment < 2 || end_segment < start_segment) {
const std::string message = "Segment numbers " + std::to_string(start_segment) + " and "
+ std::to_string(end_segment) + " specified in WSEGSICD for well " +
well_name
+ " are illegal ";
throw std::invalid_argument(message);
}
const ICD spiral_icd(record);
for (int seg = start_segment; seg <= end_segment; seg++) {
res[well_name].push_back(std::make_pair(seg, spiral_icd));
}
}
return res;
}
private:
double m_strength;
double m_length;

10
opm/parser/eclipse/EclipseState/Schedule/MSW/Segment.hpp
View File

@@ -26,6 +26,7 @@
#include <opm/parser/eclipse/EclipseState/Schedule/MSW/Valve.hpp>
#include <opm/parser/eclipse/EclipseState/Schedule/MSW/SICD.hpp>
#include <opm/parser/eclipse/EclipseState/Schedule/MSW/AICD.hpp>
namespace Opm {
@@ -105,10 +106,12 @@ namespace Opm {
bool operator!=( const Segment& ) const;
const SICD& spiralICD() const;
const AutoICD& autoICD() const;
const Valve& valve() const;
void updatePerfLength(double perf_length);
void updateSpiralICD(const SICD& spiral_icd);
void updateAutoICD(const AutoICD& aicd);
void updateValve(const Valve& valve, const double segment_length);
void updateValve(const Valve& valve);
void addInletSegment(const int segment_number);
@@ -123,6 +126,11 @@ namespace Opm {
return std::holds_alternative<SICD>(this->m_icd);
}
inline bool isAICD() const
{
return std::holds_alternative<AutoICD>(this->m_icd);
}
inline bool isValve() const
{
return std::holds_alternative<Valve>(this->m_icd);
@@ -194,7 +202,7 @@ namespace Opm {
bool m_data_ready;
std::optional<double> m_perf_length;
std::variant<RegularSegment, SICD, Valve> m_icd;
std::variant<RegularSegment, SICD, AutoICD, Valve> m_icd;
// We are not handling the length of segment projected onto the X-axis and Y-axis.
// They are not used in the simulations and we are not supporting the plotting.

69
src/opm/parser/eclipse/EclipseState/Schedule/MSW/AICD.cpp Normal file
View File

@@ -0,0 +1,69 @@
/*
Copyright 2020 Equinor 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 <opm/parser/eclipse/EclipseState/Schedule/MSW/AICD.hpp>
#include <opm/parser/eclipse/Parser/ParserKeywords/W.hpp>
namespace Opm {
AutoICD AutoICD::serializeObject() {
AutoICD aicd;
return aicd;
}
using AICD = ParserKeywords::WSEGAICD;
AutoICD::AutoICD(const DeckRecord& record) :
SICD(record),
m_flow_rate_exponent(record.getItem<AICD::FLOW_RATE_EXPONENT>().get<double>(0)),
m_visc_exponent(record.getItem<AICD::VISC_EXPONENT>().get<double>(0)),
m_oil_density_exponent(record.getItem<AICD::OIL_FLOW_FRACTION>().get<double>(0)),
m_water_density_exponent(record.getItem<AICD::WATER_FLOW_FRACTION>().get<double>(0)),
m_gas_density_exponent(record.getItem<AICD::GAS_FLOW_FRACTION>().get<double>(0)),
m_oil_viscosity_exponent(record.getItem<AICD::OIL_VISC_FRACTION>().get<double>(0)),
m_water_viscosity_exponent(record.getItem<AICD::WATER_VISC_FRACTION>().get<double>(0)),
m_gas_viscosity_exponent(record.getItem<AICD::GAS_VISC_FRACTION>().get<double>(0))
{
}
// the function will return a map
// [
// "WELL1" : [<seg1, aicd1>, <seg2, aicd2> ...]
// ....
std::map<std::string, std::vector<std::pair<int, AutoICD> > >
AutoICD::fromWSEGAICD(const DeckKeyword& wsegaicd) {
return SICD::fromWSEG<AutoICD>(wsegaicd);
}
bool AutoICD::operator==(const AutoICD& other) const {
return SICD::operator==(other) &&
this->m_flow_rate_exponent == other.m_flow_rate_exponent &&
this->m_visc_exponent == other.m_visc_exponent &&
this->m_oil_density_exponent == other.m_oil_density_exponent &&
this->m_water_density_exponent == other.m_water_density_exponent &&
this->m_gas_density_exponent == other.m_gas_density_exponent &&
this->m_oil_viscosity_exponent == other.m_oil_viscosity_exponent &&
this->m_water_viscosity_exponent == other.m_water_viscosity_exponent &&
this->m_gas_viscosity_exponent == other.m_gas_viscosity_exponent;
}
}

27
src/opm/parser/eclipse/EclipseState/Schedule/MSW/SICD.cpp
View File

@@ -101,32 +101,13 @@ namespace Opm {
return result;
}
std::map<std::string, std::vector<std::pair<int, SICD> > >
SICD::fromWSEGSICD(const DeckKeyword& wsegsicd)
{
std::map<std::string, std::vector<std::pair<int, SICD> > > res;
for (const DeckRecord &record : wsegsicd) {
const std::string well_name = record.getItem("WELL").getTrimmedString(0);
const int start_segment = record.getItem("SEG1").get<int>(0);
const int end_segment = record.getItem("SEG2").get<int>(0);
if (start_segment < 2 || end_segment < 2 || end_segment < start_segment) {
const std::string message = "Segment numbers " + std::to_string(start_segment) + " and "
+ std::to_string(end_segment) + " specified in WSEGSICD for well " +
well_name
+ " are illegal ";
throw std::invalid_argument(message);
}
const SICD spiral_icd(record);
for (int seg = start_segment; seg <= end_segment; seg++) {
res[well_name].push_back(std::make_pair(seg, spiral_icd));
}
}
return res;
return SICD::fromWSEG<SICD>(wsegsicd);
}
double SICD::maxAbsoluteRate() const {

12
src/opm/parser/eclipse/EclipseState/Schedule/MSW/Segment.cpp
View File

@@ -229,6 +229,9 @@ namespace {
if (this->isSpiralICD())
return SegmentType::SICD;
if (this->isAICD())
return SegmentType::AICD;
if (this->isValve())
return SegmentType::VALVE;
@@ -274,6 +277,15 @@ namespace {
return std::get<SICD>(this->m_icd);
}
void Segment::updateAutoICD(const AutoICD& aicd) {
this->m_icd = aicd;
}
const AutoICD& Segment::autoICD() const {
return std::get<AutoICD>(this->m_icd);
}
void Segment::updateValve(const Valve& input_valve) {
// we need to update some values for the vale

41
src/opm/parser/eclipse/share/keywords/000_Eclipse100/W/WSEGAICD
View File

@@ -9,62 +9,63 @@
"value_type": "STRING"
},
{
"name": "SEGMENT1",
"value_type": "INT",
"default": 0
"name": "SEG1",
"value_type": "INT"
},
{
"name": "SEGMENT2",
"value_type": "INT",
"default": 0
"name": "SEG2",
"value_type": "INT"
},
{
"name": "AICD_STRENGTH",
"name": "STRENGTH",
"value_type": "DOUBLE",
"default": 0
"dimension": "Pressure*Time*Time/GeometricVolume*GeometricVolume"
},
{
"name": "ICD_LENGTH",
"name": "LENGTH",
"value_type": "DOUBLE",
"dimension": "Length",
"default": 12
},
{
"name": "RHO",
"name": "DENSITY_CALI",
"value_type": "DOUBLE",
"dimension": "Mass/Length*Length*Length",
"dimension": "Density",
"default": 1000.25
},
{
"name": "VISCOSITY",
"name": "VISCOSITY_CALI",
"value_type": "DOUBLE",
"dimension": "Viscosity",
"default": 1.45
},
{
"name": "WATER_LIMIT",
"name": "CRITICAL_VALUE",
"value_type": "DOUBLE",
"dimension": "1",
"default": 0.5
},
{
"name": "TRANSITION_WIDTH",
"name": "WIDTH_TRANS",
"value_type": "DOUBLE",
"dimension": "1",
"default": 0.05
},
{
"name": "MAX_SOMETHING",
"name": "MAX_VISC_RATIO",
"value_type": "DOUBLE",
"dimension": "1",
"default": 5
},
{
"name": "SCALING_METHOD",
"name": "METHOD_SCALING_FACTOR",
"value_type": "INT",
"default": -1
},
{
"name": "MAX_QICD",
"name": "MAX_ABS_RATE",
"value_type": "DOUBLE",
"dimension": "ReservoirVolume/Time"
"dimension": "GeometricVolume/Time"
},
{
"name": "FLOW_RATE_EXPONENT",
@@ -75,7 +76,7 @@
"value_type": "DOUBLE"
},
{
"name": "ICD_FLAG",
"name": "STATUS",
"value_type": "STRING",
"default": "OPEN"
},
@@ -95,7 +96,7 @@
"default": 1
},
{
"name": "OIL_VSIC_FRACTION",
"name": "OIL_VISC_FRACTION",
"value_type": "DOUBLE",
"default": 1
},

191
tests/parser/MultisegmentWellTests.cpp
View File

@@ -46,6 +46,148 @@
#include <opm/parser/eclipse/Parser/ParseContext.hpp>
BOOST_AUTO_TEST_CASE(AICDWellTest) {
auto dir = Opm::Connection::Direction::Z;
const auto kind = Opm::Connection::CTFKind::DeckValue;
Opm::WellConnections connection_set(Opm::Connection::Order::TRACK, 10,10);
Opm::EclipseGrid grid(20,20,20);
connection_set.add(Opm::Connection( 19, 0, 0,grid.getGlobalIndex(19,0,0), 1, 0.0, Opm::Connection::State::OPEN , 200, 17.29, 0.25, 0.0, 0.0, 0, dir, kind, 0, true) );
connection_set.add(Opm::Connection( 19, 0, 1,grid.getGlobalIndex(19,0,1), 1, 0.0, Opm::Connection::State::OPEN , 200, 17.29, 0.25, 0.0, 0.0, 0, dir, kind, 0, true) );
connection_set.add(Opm::Connection( 19, 0, 2,grid.getGlobalIndex(19,0,2), 1, 0.0, Opm::Connection::State::OPEN , 200, 17.29, 0.25, 0.0, 0.0, 0, dir, kind, 0, true) );
connection_set.add(Opm::Connection( 18, 0, 1,grid.getGlobalIndex(18,0,1), 1, 0.0, Opm::Connection::State::OPEN , 200, 17.29, 0.25, 0.0, 0.0, 0, Opm::Connection::Direction::X, kind, 0, true) );
connection_set.add(Opm::Connection( 17, 0, 1,grid.getGlobalIndex(17,0,1), 1, 0.0, Opm::Connection::State::OPEN , 200, 17.29, 0.25, 0.0, 0.0, 0, Opm::Connection::Direction::X, kind, 0, true) );
connection_set.add(Opm::Connection( 16, 0, 1,grid.getGlobalIndex(16,0,1), 1, 0.0, Opm::Connection::State::OPEN , 200, 17.29, 0.25, 0.0, 0.0, 0, Opm::Connection::Direction::X, kind, 0, true) );
connection_set.add(Opm::Connection( 15, 0, 1,grid.getGlobalIndex(15,0,1), 1, 0.0, Opm::Connection::State::OPEN , 200, 17.29, 0.25, 0.0, 0.0, 0, Opm::Connection::Direction::X, kind, 0, true) );
BOOST_CHECK_EQUAL( 7U , connection_set.size() );
const std::string compsegs_string = R"(
WELSEGS
'PROD01' 2512.5 2512.5 1.0e-5 'ABS' 'HF-' 'HO' /
2 2 1 1 2537.5 2537.5 0.3 0.00010 /
3 3 1 2 2562.5 2562.5 0.2 0.00010 /
4 4 2 2 2737.5 2537.5 0.2 0.00010 /
6 6 2 4 3037.5 2539.5 0.2 0.00010 /
7 7 2 6 3337.5 2534.5 0.2 0.00010 /
8 8 3 7 3337.6 2534.5 0.2 0.00015 /
/
COMPSEGS
PROD01 /
20 1 1 1 2512.5 2525.0 /
20 1 2 1 2525.0 2550.0 /
20 1 3 1 2550.0 2575.0 /
19 1 2 2 2637.5 2837.5 /
18 1 2 2 2837.5 3037.5 /
17 1 2 2 3037.5 3237.5 /
16 1 2 3 3237.5 3437.5 /
/
WSEGAICD
'PROD01' 8 8 0.002 -0.7 1* 1* 0.6 1* 1* 2* 1.0 1.0 'SHUT' /
/
)";
Opm::Parser parser;
Opm::Deck deck = parser.parseString(compsegs_string);
const Opm::DeckKeyword compsegs = deck.getKeyword("COMPSEGS");
BOOST_CHECK_EQUAL( 8U, compsegs.size() );
const Opm::DeckKeyword welsegs = deck.getKeyword("WELSEGS");
Opm::WellSegments segment_set(welsegs);
BOOST_CHECK_EQUAL(7U, segment_set.size());
Opm::ErrorGuard errorGuard;
Opm::ParseContext parseContext;
parseContext.update(Opm::ParseContext::SCHEDULE_COMPSEGS_INVALID, Opm::InputError::THROW_EXCEPTION);
parseContext.update(Opm::ParseContext::SCHEDULE_COMPSEGS_NOT_SUPPORTED, Opm::InputError::THROW_EXCEPTION);
const auto& [new_connection_set, new_segment_set] = Opm::Compsegs::processCOMPSEGS(compsegs, connection_set, segment_set, grid, parseContext, errorGuard);
// checking the ICD segment
const Opm::DeckKeyword wsegaicd = deck.getKeyword("WSEGAICD");
const auto aicd_map = Opm::AutoICD::fromWSEGAICD(wsegaicd);
BOOST_CHECK_EQUAL(1U, aicd_map.size());
const auto it = aicd_map.begin();
const std::string& well_name = it->first;
BOOST_CHECK_EQUAL(well_name, "PROD01");
const auto& aicd_vector = it->second;
BOOST_CHECK_EQUAL(1U, aicd_vector.size());
const int segment_number = aicd_vector[0].first;
const Opm::AutoICD& aicd0 = aicd_vector[0].second;
BOOST_CHECK_EQUAL(8, segment_number);
Opm::Segment segment = segment_set.getFromSegmentNumber(segment_number);
segment.updateAutoICD(aicd0);
BOOST_CHECK(Opm::Segment::SegmentType::AICD == segment.segmentType());
auto aicd = segment.autoICD();
BOOST_CHECK_GT(aicd.maxAbsoluteRate(), 1.e99);
BOOST_CHECK(aicd.status()==Opm::ICDStatus::SHUT);
// 0.002 bars*day*day/Volume^2
BOOST_CHECK_EQUAL(aicd.strength(), 0.002*1.e5*86400.*86400.);
BOOST_CHECK_EQUAL(aicd.length(), -0.7);
BOOST_CHECK_EQUAL(aicd.densityCalibration(), 1000.25);
// 1.45 cp
BOOST_CHECK_EQUAL(aicd.viscosityCalibration(), 1.45 * 0.001);
BOOST_CHECK_EQUAL(aicd.criticalValue(), 0.6);
BOOST_CHECK_EQUAL(aicd.widthTransitionRegion(), 0.05);
BOOST_CHECK_EQUAL(aicd.maxViscosityRatio(), 5.0);
BOOST_CHECK_EQUAL(aicd.methodFlowScaling(), -1);
// the scaling factor has not been updated properly, so it will throw
BOOST_CHECK_THROW(aicd.scalingFactor(), std::runtime_error);
const int outlet_segment_number = segment.outletSegment();
const double outlet_segment_length = segment_set.segmentLength(outlet_segment_number);
// only one connection attached to the outlet segment in this case
const Opm::Connection& connection = new_connection_set.getFromIJK(15, 0, 1);
const auto& perf_range = connection.perf_range();
const auto connection_length = perf_range->second - perf_range->first;
aicd.updateScalingFactor(outlet_segment_length, connection_length);
// updated, so it should not throw
BOOST_CHECK_NO_THROW(aicd.scalingFactor());
BOOST_CHECK_EQUAL(0.7, aicd.scalingFactor());
BOOST_CHECK_EQUAL(7U, new_connection_set.size());
const Opm::Connection& connection1 = new_connection_set.get(0);
const int segment_number_connection1 = connection1.segment();
const double center_depth_connection1 = connection1.depth();
BOOST_CHECK_EQUAL(segment_number_connection1, 1);
BOOST_CHECK_EQUAL(center_depth_connection1, 2512.5);
const Opm::Connection& connection3 = new_connection_set.get(2);
const int segment_number_connection3 = connection3.segment();
const double center_depth_connection3 = connection3.depth();
BOOST_CHECK_EQUAL(segment_number_connection3, 3);
BOOST_CHECK_EQUAL(center_depth_connection3, 2562.5);
const Opm::Connection& connection5 = new_connection_set.get(4);
const int segment_number_connection5 = connection5.segment();
const double center_depth_connection5 = connection5.depth();
BOOST_CHECK_EQUAL(segment_number_connection5, 6);
BOOST_CHECK_CLOSE(center_depth_connection5, 2538.83, 0.001);
const Opm::Connection& connection6 = new_connection_set.get(5);
const int segment_number_connection6 = connection6.segment();
const double center_depth_connection6 = connection6.depth();
BOOST_CHECK_EQUAL(segment_number_connection6, 6);
BOOST_CHECK_CLOSE(center_depth_connection6, 2537.83, 0.001);
const Opm::Connection& connection7 = new_connection_set.get(6);
const int segment_number_connection7 = connection7.segment();
const double center_depth_connection7 = connection7.depth();
BOOST_CHECK_EQUAL(segment_number_connection7, 8);
BOOST_CHECK_EQUAL(center_depth_connection7, 2534.5);
}
BOOST_AUTO_TEST_CASE(MultisegmentWellTest) {
@@ -64,30 +206,31 @@ BOOST_AUTO_TEST_CASE(MultisegmentWellTest) {
BOOST_CHECK_EQUAL( 7U , connection_set.size() );
const std::string compsegs_string =
"WELSEGS \n"
"'PROD01' 2512.5 2512.5 1.0e-5 'ABS' 'HF-' 'HO' /\n"
"2 2 1 1 2537.5 2537.5 0.3 0.00010 /\n"
"3 3 1 2 2562.5 2562.5 0.2 0.00010 /\n"
"4 4 2 2 2737.5 2537.5 0.2 0.00010 /\n"
"6 6 2 4 3037.5 2539.5 0.2 0.00010 /\n"
"7 7 2 6 3337.5 2534.5 0.2 0.00010 /\n"
"8 8 3 7 3337.6 2534.5 0.2 0.00015 /\n"
"/\n"
"\n"
"COMPSEGS\n"
"PROD01 / \n"
"20 1 1 1 2512.5 2525.0 /\n"
"20 1 2 1 2525.0 2550.0 /\n"
"20 1 3 1 2550.0 2575.0 /\n"
"19 1 2 2 2637.5 2837.5 /\n"
"18 1 2 2 2837.5 3037.5 /\n"
"17 1 2 2 3037.5 3237.5 /\n"
"16 1 2 3 3237.5 3437.5 /\n"
"/\n"
"WSEGSICD\n"
"'PROD01' 8 8 0.002 -0.7 1* 1* 0.6 1* 1* 2* 'SHUT' /\n"
"/\n";
const std::string compsegs_string = R"(
WELSEGS
'PROD01' 2512.5 2512.5 1.0e-5 'ABS' 'HF-' 'HO' /
2 2 1 1 2537.5 2537.5 0.3 0.00010 /
3 3 1 2 2562.5 2562.5 0.2 0.00010 /
4 4 2 2 2737.5 2537.5 0.2 0.00010 /
6 6 2 4 3037.5 2539.5 0.2 0.00010 /
7 7 2 6 3337.5 2534.5 0.2 0.00010 /
8 8 3 7 3337.6 2534.5 0.2 0.00015 /
/
COMPSEGS
PROD01 /
20 1 1 1 2512.5 2525.0 /
20 1 2 1 2525.0 2550.0 /
20 1 3 1 2550.0 2575.0 /
19 1 2 2 2637.5 2837.5 /
18 1 2 2 2837.5 3037.5 /
17 1 2 2 3037.5 3237.5 /
16 1 2 3 3237.5 3437.5 /
/
WSEGSICD
'PROD01' 8 8 0.002 -0.7 1* 1* 0.6 1* 1* 2* 'SHUT' /
/
)";
Opm::Parser parser;
Opm::Deck deck = parser.parseString(compsegs_string);