OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

Use opm parser to find all VFP data based on a *.DATA file

Browse Source
This commit is contained in:
Magne Sjaastad
2024-05-03 09:52:43 +02:00
parent 24aa7fca3e
commit 6d478ddb3c
2 changed files with 125 additions and 47 deletions
Download Patch File Download Diff File Expand all files Collapse all files

151
ApplicationLibCode/ProjectDataModel/RiaOpmParserTools.cpp
View File

@@ -30,6 +30,7 @@
#include "opm/input/eclipse/Parser/Parser.hpp"
#include "opm/input/eclipse/Parser/ParserKeywords/I.hpp"
#include "opm/input/eclipse/Parser/ParserKeywords/P.hpp"
#include "opm/input/eclipse/Parser/ParserKeywords/T.hpp"
#include "opm/input/eclipse/Parser/ParserKeywords/V.hpp"
#include "opm/input/eclipse/Parser/ParserKeywords/W.hpp"
@@ -37,10 +38,55 @@
#include <set>
namespace RiaOpmParserTools
{
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::vector<Opm::VFPInjTable> RiaOpmParserTools::extractVfpInjectionTables( const std::string& filename )
Opm::VFPInjTable createInjectionTable( const Opm::DeckKeyword& keyword )
{
Opm::UnitSystem unitSystem;
{
const auto& header = keyword.getRecord( 0 );
if ( header.getItem<Opm::ParserKeywords::VFPINJ::UNITS>().hasValue( 0 ) )
{
std::string units_string;
units_string = header.getItem<Opm::ParserKeywords::VFPINJ::UNITS>().get<std::string>( 0 );
unitSystem = Opm::UnitSystem( units_string );
}
}
return { keyword, unitSystem };
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
Opm::VFPProdTable createProductionTable( const Opm::DeckKeyword& keyword )
{
Opm::UnitSystem unitSystem;
{
const auto& header = keyword.getRecord( 0 );
if ( header.getItem<Opm::ParserKeywords::VFPPROD::UNITS>().hasValue( 0 ) )
{
std::string units_string;
units_string = header.getItem<Opm::ParserKeywords::VFPPROD::UNITS>().get<std::string>( 0 );
unitSystem = Opm::UnitSystem( units_string );
}
}
bool gaslift_opt_active = false;
return { keyword, gaslift_opt_active, unitSystem };
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::vector<Opm::VFPInjTable> extractVfpInjectionTables( const std::string& filename )
{
std::vector<Opm::VFPInjTable> tables;
@@ -60,21 +106,7 @@ std::vector<Opm::VFPInjTable> RiaOpmParserTools::extractVfpInjectionTables( cons
for ( auto kw : keywordList )
{
auto name = kw->name();
Opm::UnitSystem unitSystem;
{
const auto& header = kw->getRecord( 0 );
if ( header.getItem<Opm::ParserKeywords::VFPINJ::UNITS>().hasValue( 0 ) )
{
std::string units_string;
units_string = header.getItem<Opm::ParserKeywords::VFPINJ::UNITS>().get<std::string>( 0 );
unitSystem = Opm::UnitSystem( units_string );
}
}
Opm::VFPInjTable table( *kw, unitSystem );
auto table = createInjectionTable( *kw );
tables.push_back( table );
}
}
@@ -88,7 +120,7 @@ std::vector<Opm::VFPInjTable> RiaOpmParserTools::extractVfpInjectionTables( cons
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::vector<Opm::VFPProdTable> RiaOpmParserTools::extractVfpProductionTables( const std::string& filename )
std::vector<Opm::VFPProdTable> extractVfpProductionTables( const std::string& filename )
{
std::vector<Opm::VFPProdTable> tables;
@@ -106,22 +138,7 @@ std::vector<Opm::VFPProdTable> RiaOpmParserTools::extractVfpProductionTables( co
for ( auto kw : keywordList )
{
auto name = kw->name();
Opm::UnitSystem unitSystem;
{
const auto& header = kw->getRecord( 0 );
if ( header.getItem<Opm::ParserKeywords::VFPPROD::UNITS>().hasValue( 0 ) )
{
std::string units_string;
units_string = header.getItem<Opm::ParserKeywords::VFPPROD::UNITS>().get<std::string>( 0 );
unitSystem = Opm::UnitSystem( units_string );
}
}
bool gaslift_opt_active = false;
Opm::VFPProdTable table( *kw, gaslift_opt_active, unitSystem );
auto table = createProductionTable( *kw );
tables.push_back( table );
}
}
@@ -135,7 +152,65 @@ std::vector<Opm::VFPProdTable> RiaOpmParserTools::extractVfpProductionTables( co
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::map<std::string, std::vector<std::pair<int, int>>> RiaOpmParserTools::extractWseglink( const std::string& filename )
std::pair<std::vector<Opm::VFPProdTable>, std::vector<Opm::VFPInjTable>> extractVfpTablesFromDataFile( const std::string& dataDeckFilename )
{
if ( !std::filesystem::exists( dataDeckFilename ) ) return {};
std::vector<Opm::VFPProdTable> prodTables;
std::vector<Opm::VFPInjTable> injTables;
try
{
Opm::Parser parser( false );
// Required to include the TUNING keyword to avoid parsing error of a Norne DATA file containing the TUNING keyword
// The TUNING keyword is not required nor related to VFP data
std::vector<Opm::ParserKeyword> parserKeywords = { Opm::ParserKeywords::VFPPROD(),
Opm::ParserKeywords::VFPINJ(),
Opm::ParserKeywords::INCLUDE(),
Opm::ParserKeywords::TUNING() };
for ( const auto& kw : parserKeywords )
{
parser.addParserKeyword( kw );
}
Opm::ParseContext parseContext( Opm::InputError::Action::WARN );
auto deck = parser.parseFile( dataDeckFilename, parseContext );
{
std::string prodKeyword = "VFPPROD";
auto keywordList = deck.getKeywordList( prodKeyword );
for ( auto kw : keywordList )
{
auto table = createProductionTable( *kw );
prodTables.push_back( table );
}
}
{
std::string injKeyword = "VFPINJ";
auto keywordList = deck.getKeywordList( injKeyword );
for ( auto kw : keywordList )
{
auto table = createInjectionTable( *kw );
injTables.push_back( table );
}
}
}
catch ( ... )
{
QString text = QString( "Error detected when parsing '%1'. Imported data might be missing or incomplete." )
.arg( QString::fromStdString( dataDeckFilename ) );
RiaLogging::warning( text );
}
return { prodTables, injTables };
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::map<std::string, std::vector<std::pair<int, int>>> extractWseglink( const std::string& filename )
{
if ( !std::filesystem::exists( filename ) ) return {};
@@ -205,7 +280,7 @@ std::map<std::string, std::vector<std::pair<int, int>>> RiaOpmParserTools::extra
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::vector<RiaOpmParserTools::AicdTemplateValues> RiaOpmParserTools::extractWsegAicd( const std::string& filename )
std::vector<RiaOpmParserTools::AicdTemplateValues> extractWsegAicd( const std::string& filename )
{
if ( !std::filesystem::exists( filename ) ) return {};
@@ -291,7 +366,7 @@ std::vector<RiaOpmParserTools::AicdTemplateValues> RiaOpmParserTools::extractWse
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::vector<RiaOpmParserTools::AicdTemplateValues> RiaOpmParserTools::extractWsegAicdCompletor( const std::string& filename )
std::vector<RiaOpmParserTools::AicdTemplateValues> extractWsegAicdCompletor( const std::string& filename )
{
QFile file( QString::fromStdString( filename ) );
if ( !file.open( QFile::ReadOnly ) ) return {};
@@ -346,7 +421,9 @@ std::vector<RiaOpmParserTools::AicdTemplateValues> RiaOpmParserTools::extractWse
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::string RiaOpmParserTools::aicdTemplateId()
std::string aicdTemplateId()
{
return "ID_NUMBER";
}
} // namespace RiaOpmParserTools

21
ApplicationLibCode/ProjectDataModel/RiaOpmParserTools.h
View File

@@ -28,17 +28,18 @@
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
class RiaOpmParserTools
namespace RiaOpmParserTools
{
public:
static std::vector<Opm::VFPInjTable> extractVfpInjectionTables( const std::string& filename );
static std::vector<Opm::VFPProdTable> extractVfpProductionTables( const std::string& filename );
std::vector<Opm::VFPInjTable> extractVfpInjectionTables( const std::string& filename );
std::vector<Opm::VFPProdTable> extractVfpProductionTables( const std::string& filename );
static std::map<std::string, std::vector<std::pair<int, int>>> extractWseglink( const std::string& filename );
std::pair<std::vector<Opm::VFPProdTable>, std::vector<Opm::VFPInjTable>> extractVfpTablesFromDataFile( const std::string& dataDeckFilename );
using AicdTemplateValues = std::map<std::string, double>;
static std::vector<AicdTemplateValues> extractWsegAicd( const std::string& filename );
static std::vector<AicdTemplateValues> extractWsegAicdCompletor( const std::string& filename );
std::map<std::string, std::vector<std::pair<int, int>>> extractWseglink( const std::string& filename );
static std::string aicdTemplateId();
};
using AicdTemplateValues = std::map<std::string, double>;
std::vector<AicdTemplateValues> extractWsegAicd( const std::string& filename );
std::vector<AicdTemplateValues> extractWsegAicdCompletor( const std::string& filename );
std::string aicdTemplateId();
}; // namespace RiaOpmParserTools