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

Add INVALID_WELL context for WCONINJ(E/H)

Browse Source
This commit is contained in:
Lars Petter Øren Hauge
2018-04-26 10:39:54 +02:00
parent 183fd4eabb
commit ba107498a1
3 changed files with 64 additions and 32 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
opm/parser/eclipse/EclipseState/Schedule/Schedule.hpp
View File

@@ -154,12 +154,12 @@ namespace Opm
void handleCOMPLUMP( const DeckKeyword& keyword, size_t currentStep );
void handleWELSEGS( const DeckKeyword& keyword, size_t currentStep);
void handleCOMPSEGS( const DeckKeyword& keyword, size_t currentStep);
void handleWCONINJE( const SCHEDULESection&, const DeckKeyword& keyword, size_t currentStep);
void handleWCONINJE( const SCHEDULESection&, const DeckKeyword& keyword, size_t currentStep, const ParseContext& parseContext);
void handleWPOLYMER( const DeckKeyword& keyword, size_t currentStep);
void handleWSOLVENT( const DeckKeyword& keyword, size_t currentStep);
void handleWTEMP( const DeckKeyword& keyword, size_t currentStep);
void handleWINJTEMP( const DeckKeyword& keyword, size_t currentStep);
void handleWCONINJH( const SCHEDULESection&, const DeckKeyword& keyword, size_t currentStep);
void handleWCONINJH( const SCHEDULESection&, const DeckKeyword& keyword, size_t currentStep, const ParseContext& parseContext);
void handleWELOPEN( const DeckKeyword& keyword, size_t currentStep );
void handleWELTARG( const SCHEDULESection&, const DeckKeyword& keyword, size_t currentStep);
void handleGCONINJE( const SCHEDULESection&, const DeckKeyword& keyword, size_t currentStep);

70
src/opm/parser/eclipse/EclipseState/Schedule/Schedule.cpp
View File

@@ -172,7 +172,7 @@ namespace Opm {
handleWCONPROD(keyword, currentStep, parseContext);
else if (keyword.name() == "WCONINJE")
handleWCONINJE(section, keyword, currentStep);
handleWCONINJE(section, keyword, currentStep, parseContext);
else if (keyword.name() == "WPOLYMER")
handleWPOLYMER(keyword, currentStep);
@@ -187,7 +187,7 @@ namespace Opm {
handleWINJTEMP(keyword, currentStep);
else if (keyword.name() == "WCONINJH")
handleWCONINJH(section, keyword, currentStep);
handleWCONINJH(section, keyword, currentStep, parseContext);
else if (keyword.name() == "WGRUPCON")
handleWGRUPCON(keyword, currentStep);
@@ -612,11 +612,15 @@ namespace Opm {
}
void Schedule::handleWCONINJE( const SCHEDULESection& section, const DeckKeyword& keyword, size_t currentStep) {
void Schedule::handleWCONINJE( const SCHEDULESection& section, const DeckKeyword& keyword, size_t currentStep, const ParseContext& parseContext) {
for( const auto& record : keyword ) {
const std::string& wellNamePattern = record.getItem("WELL").getTrimmedString(0);
for( auto* well : getWells( wellNamePattern ) ) {
auto wells = getWells(wellNamePattern);
if (wells.empty())
InvalidWellPattern(wellNamePattern, parseContext, keyword);
for( auto* well : wells) {
WellInjector::TypeEnum injectorType = WellInjector::TypeFromString( record.getItem("TYPE").getTrimmedString(0) );
WellCommon::StatusEnum status = WellCommon::StatusFromString( record.getItem("STATUS").getTrimmedString(0));
@@ -800,9 +804,9 @@ namespace Opm {
}
}
void Schedule::handleWCONINJH( const SCHEDULESection& section, const DeckKeyword& keyword, size_t currentStep) {
void Schedule::handleWCONINJH( const SCHEDULESection& section, const DeckKeyword& keyword, size_t currentStep, const ParseContext& parseContext) {
for( const auto& record : keyword ) {
const std::string& wellName = record.getItem("WELL").getTrimmedString(0);
const std::string& wellNamePattern = record.getItem("WELL").getTrimmedString(0);
// convert injection rates to SI
WellInjector::TypeEnum injectorType = WellInjector::TypeFromString( record.getItem("TYPE").getTrimmedString(0));
@@ -811,35 +815,41 @@ namespace Opm {
WellCommon::StatusEnum status = WellCommon::StatusFromString( record.getItem("STATUS").getTrimmedString(0));
auto& well = this->m_wells.get( wellName );
updateWellStatus( well, currentStep, status );
WellInjectionProperties properties(well.getInjectionPropertiesCopy(currentStep));
auto wells = getWells( wellNamePattern );
properties.injectorType = injectorType;
if (wells.empty())
InvalidWellPattern( wellNamePattern, parseContext, keyword);
const std::string& cmodeString = record.getItem("CMODE").getTrimmedString(0);
WellInjector::ControlModeEnum controlMode = WellInjector::ControlModeFromString( cmodeString );
if (!record.getItem("RATE").defaultApplied(0)) {
properties.surfaceInjectionRate = injectionRate;
properties.addInjectionControl(controlMode);
properties.controlMode = controlMode;
}
properties.predictionMode = false;
for (auto* well : wells) {
updateWellStatus( *well, currentStep, status );
WellInjectionProperties properties(well->getInjectionPropertiesCopy(currentStep));
if ( record.getItem( "BHP" ).hasValue(0) )
properties.BHPH = record.getItem("BHP").getSIDouble(0);
if ( record.getItem( "THP" ).hasValue(0) )
properties.THPH = record.getItem("THP").getSIDouble(0);
properties.injectorType = injectorType;
if (well.setInjectionProperties(currentStep, properties))
m_events.addEvent( ScheduleEvents::INJECTION_UPDATE , currentStep );
const std::string& cmodeString = record.getItem("CMODE").getTrimmedString(0);
WellInjector::ControlModeEnum controlMode = WellInjector::ControlModeFromString( cmodeString );
if (!record.getItem("RATE").defaultApplied(0)) {
properties.surfaceInjectionRate = injectionRate;
properties.addInjectionControl(controlMode);
properties.controlMode = controlMode;
}
properties.predictionMode = false;
if ( ! well.getAllowCrossFlow() && (injectionRate == 0) ) {
std::string msg =
"Well " + well.name() + " is an injector with zero rate where crossflow is banned. " +
"This well will be closed at " + std::to_string ( m_timeMap.getTimePassedUntil(currentStep) / (60*60*24) ) + " days";
OpmLog::note(msg);
updateWellStatus( well, currentStep, WellCommon::StatusEnum::SHUT );
if ( record.getItem( "BHP" ).hasValue(0) )
properties.BHPH = record.getItem("BHP").getSIDouble(0);
if ( record.getItem( "THP" ).hasValue(0) )
properties.THPH = record.getItem("THP").getSIDouble(0);
if (well->setInjectionProperties(currentStep, properties))
m_events.addEvent( ScheduleEvents::INJECTION_UPDATE , currentStep );
if ( ! well->getAllowCrossFlow() && (injectionRate == 0) ) {
std::string msg =
"Well " + well->name() + " is an injector with zero rate where crossflow is banned. " +
"This well will be closed at " + std::to_string ( m_timeMap.getTimePassedUntil(currentStep) / (60*60*24) ) + " days";
OpmLog::note(msg);
updateWellStatus( *well, currentStep, WellCommon::StatusEnum::SHUT );
}
}
}
}

22
tests/parser/ParseContextTests.cpp
View File

@@ -526,6 +526,28 @@ BOOST_AUTO_TEST_CASE( test_invalid_wtemplate_config ) {
)";
testSamples.push_back(testSample);
// Invalid well name in WCONINJE
testSample = R"(
COMPDAT
'INJ' 10 10 3 3 'OPEN' 1* 1* 0.5 /
/
WCONINJE
'SOMETHINGELSE' 'WATER' 'OPEN' 'RATE' 20000 4* 1000 /
/
)";
testSamples.push_back(testSample);
// Invalid well name in WCONINJH
testSample = R"(
COMPDAT
'INJ' 10 10 3 3 'OPEN' 1* 1* 0.5 /
/
WCONINJH
'SOMETHINGELSE' 'WAT' 'OPEN' 20000 /
/
)";
testSamples.push_back(testSample);
std::string deckinput;
for (std::string sample : testSamples) {