368 lines
15 KiB
C++
368 lines
15 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 <opm/parser/eclipse/EclipseState/Schedule/Schedule.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/Schedule/TimeMap.hpp>
|
|
#include <iostream>
|
|
|
|
namespace Opm {
|
|
|
|
Schedule::Schedule(DeckConstPtr deck) {
|
|
if (deck->hasKeyword("SCHEDULE")) {
|
|
initFromDeck(deck);
|
|
} else
|
|
throw std::invalid_argument("Deck does not contain SCHEDULE section.\n");
|
|
}
|
|
|
|
void Schedule::initFromDeck(DeckConstPtr deck) {
|
|
createTimeMap(deck);
|
|
addGroup( "FIELD", 0 );
|
|
initRootGroupTreeNode(getTimeMap());
|
|
iterateScheduleSection(deck);
|
|
}
|
|
|
|
void Schedule::initRootGroupTreeNode(TimeMapConstPtr timeMap) {
|
|
m_rootGroupTree.reset(new DynamicState<GroupTreePtr>(timeMap, GroupTreePtr(new GroupTree())));
|
|
}
|
|
|
|
void Schedule::createTimeMap(DeckConstPtr deck) {
|
|
boost::gregorian::date startDate(defaultStartDate);
|
|
if (deck->hasKeyword("START")) {
|
|
DeckKeywordConstPtr startKeyword = deck->getKeyword("START");
|
|
startDate = TimeMap::dateFromEclipse(startKeyword->getRecord(0));
|
|
}
|
|
|
|
m_timeMap.reset(new TimeMap(startDate));
|
|
}
|
|
|
|
void Schedule::iterateScheduleSection(DeckConstPtr deck) {
|
|
DeckKeywordConstPtr scheduleKeyword = deck->getKeyword("SCHEDULE");
|
|
size_t deckIndex = scheduleKeyword->getDeckIndex() + 1;
|
|
size_t currentStep = 0;
|
|
while (deckIndex < deck->size()) {
|
|
|
|
DeckKeywordConstPtr keyword = deck->getKeyword(deckIndex);
|
|
|
|
if (keyword->name() == "DATES") {
|
|
handleDATES(keyword);
|
|
currentStep += keyword->size();
|
|
}
|
|
|
|
if (keyword->name() == "TSTEP") {
|
|
handleTSTEP(keyword);
|
|
currentStep += keyword->getRecord(0)->getItem(0)->size(); // This is a bit weird API.
|
|
}
|
|
|
|
if (keyword->name() == "WELSPECS") {
|
|
handleWELSPECS(keyword, currentStep);
|
|
}
|
|
|
|
if (keyword->name() == "WCONHIST")
|
|
handleWCONHIST(keyword, currentStep);
|
|
|
|
if (keyword->name() == "WCONPROD")
|
|
handleWCONPROD(keyword, currentStep);
|
|
|
|
if (keyword->name() == "WCONINJE")
|
|
handleWCONINJE(keyword, currentStep);
|
|
|
|
if (keyword->name() == "WCONINJH")
|
|
handleWCONINJH(keyword, currentStep);
|
|
|
|
if (keyword->name() == "COMPDAT")
|
|
handleCOMPDAT(keyword, currentStep);
|
|
|
|
|
|
if (keyword->name() == "GRUPTREE")
|
|
handleGRUPTREE(keyword, currentStep);
|
|
|
|
if (keyword->name() == "GCONINJE")
|
|
handleGCONINJE( keyword , currentStep );
|
|
|
|
if (keyword->name() == "GCONPROD")
|
|
handleGCONPROD( keyword , currentStep );
|
|
|
|
deckIndex++;
|
|
}
|
|
}
|
|
|
|
void Schedule::handleDATES(DeckKeywordConstPtr keyword) {
|
|
m_timeMap->addFromDATESKeyword(keyword);
|
|
}
|
|
|
|
void Schedule::handleTSTEP(DeckKeywordConstPtr keyword) {
|
|
m_timeMap->addFromTSTEPKeyword(keyword);
|
|
}
|
|
|
|
bool Schedule::handleGroupFromWELSPECS(const std::string& groupName, GroupTreePtr newTree) const {
|
|
bool treeUpdated = false;
|
|
if (!newTree->getNode(groupName)) {
|
|
treeUpdated = true;
|
|
newTree->updateTree(groupName);
|
|
}
|
|
return treeUpdated;
|
|
}
|
|
|
|
void Schedule::handleWELSPECS(DeckKeywordConstPtr keyword, size_t currentStep) {
|
|
bool needNewTree = false;
|
|
GroupTreePtr newTree = m_rootGroupTree->get(currentStep)->deepCopy();
|
|
|
|
for (size_t recordNr = 0; recordNr < keyword->size(); recordNr++) {
|
|
DeckRecordConstPtr record = keyword->getRecord(recordNr);
|
|
const std::string& wellName = record->getItem("WELL")->getString(0);
|
|
const std::string& groupName = record->getItem("GROUP")->getString(0);
|
|
|
|
if (!hasGroup(groupName)) {
|
|
addGroup(groupName , currentStep);
|
|
}
|
|
|
|
if (!hasWell(wellName)) {
|
|
addWell(wellName, record, currentStep);
|
|
}
|
|
|
|
WellConstPtr currentWell = getWell(wellName);
|
|
checkWELSPECSConsistency(currentWell, record);
|
|
|
|
addWellToGroup( getGroup(groupName) , getWell(wellName) , currentStep);
|
|
needNewTree = handleGroupFromWELSPECS(record->getItem(1)->getString(0), newTree);
|
|
}
|
|
if (needNewTree) {
|
|
m_rootGroupTree->add(currentStep, newTree);
|
|
}
|
|
}
|
|
|
|
void Schedule::checkWELSPECSConsistency(WellConstPtr well, DeckRecordConstPtr record) const {
|
|
if (well->getHeadI() != record->getItem("HEAD_I")->getInt(0)) {
|
|
throw std::invalid_argument("Unable process WELSPECS for well " + well->name() + ", HEAD_I deviates from existing value");
|
|
}
|
|
if (well->getHeadJ() != record->getItem("HEAD_J")->getInt(0)) {
|
|
throw std::invalid_argument("Unable process WELSPECS for well " + well->name() + ", HEAD_J deviates from existing value");
|
|
}
|
|
if (well->getRefDepth() != record->getItem("REF_DEPTH")->getSIDouble(0)) {
|
|
throw std::invalid_argument("Unable process WELSPECS for well " + well->name() + ", REF_DEPTH deviates from existing value");
|
|
}
|
|
}
|
|
|
|
void Schedule::handleWCONProducer(DeckKeywordConstPtr keyword, size_t currentStep, bool isPredictionMode) {
|
|
for (size_t recordNr = 0; recordNr < keyword->size(); recordNr++) {
|
|
DeckRecordConstPtr record = keyword->getRecord(recordNr);
|
|
const std::string& wellName = record->getItem("WELL")->getString(0);
|
|
WellPtr well = getWell(wellName);
|
|
double orat = record->getItem("ORAT")->getSIDouble(0);
|
|
double wrat = record->getItem("WRAT")->getSIDouble(0);
|
|
double grat = record->getItem("GRAT")->getSIDouble(0);
|
|
|
|
well->setOilRate(currentStep, orat);
|
|
well->setWaterRate(currentStep, wrat);
|
|
well->setGasRate(currentStep, grat);
|
|
well->setInPredictionMode(currentStep, isPredictionMode);
|
|
}
|
|
}
|
|
|
|
void Schedule::handleWCONHIST(DeckKeywordConstPtr keyword, size_t currentStep) {
|
|
handleWCONProducer(keyword, currentStep, false);
|
|
}
|
|
|
|
void Schedule::handleWCONPROD(DeckKeywordConstPtr keyword, size_t currentStep) {
|
|
handleWCONProducer(keyword, currentStep, true);
|
|
}
|
|
|
|
void Schedule::handleWCONINJE(DeckKeywordConstPtr keyword, size_t currentStep) {
|
|
for (size_t recordNr = 0; recordNr < keyword->size(); recordNr++) {
|
|
DeckRecordConstPtr record = keyword->getRecord(recordNr);
|
|
const std::string& wellName = record->getItem("WELL")->getString(0);
|
|
WellPtr well = getWell(wellName);
|
|
double surfaceInjectionRate = record->getItem("SURFACE_FLOW_TARGET")->getSIDouble(0);
|
|
double reservoirInjectionRate = record->getItem("RESV_FLOW_TARGET")->getSIDouble(0);
|
|
|
|
well->setSurfaceInjectionRate( currentStep , surfaceInjectionRate );
|
|
well->setReservoirInjectionRate( currentStep , reservoirInjectionRate );
|
|
well->setInPredictionMode(currentStep, true);
|
|
}
|
|
}
|
|
|
|
void Schedule::handleWCONINJH(DeckKeywordConstPtr keyword, size_t currentStep) {
|
|
for (size_t recordNr = 0; recordNr < keyword->size(); recordNr++) {
|
|
DeckRecordConstPtr record = keyword->getRecord(recordNr);
|
|
const std::string& wellName = record->getItem("WELL")->getString(0);
|
|
WellPtr well = getWell(wellName);
|
|
double injectionRate = record->getItem("RATE")->getSIDouble(0);
|
|
|
|
well->setSurfaceInjectionRate( currentStep , injectionRate );
|
|
well->setInPredictionMode(currentStep, false );
|
|
}
|
|
}
|
|
|
|
|
|
void Schedule::handleGCONINJE(DeckKeywordConstPtr keyword, size_t currentStep) {
|
|
for (size_t recordNr = 0; recordNr < keyword->size(); recordNr++) {
|
|
DeckRecordConstPtr record = keyword->getRecord(recordNr);
|
|
const std::string& groupName = record->getItem("GROUP")->getString(0);
|
|
GroupPtr group = getGroup(groupName);
|
|
|
|
{
|
|
PhaseEnum phase = PhaseEnumFromString( record->getItem("PHASE")->getString(0) );
|
|
group->setInjectionPhase( currentStep , phase );
|
|
}
|
|
{
|
|
GroupInjection::ControlEnum controlMode = GroupInjection::ControlEnumFromString( record->getItem("CONTROL_MODE")->getString(0) );
|
|
group->setInjectionControlMode( currentStep , controlMode );
|
|
}
|
|
group->setSurfaceMaxRate( currentStep , record->getItem("SURFACE_TARGET")->getSIDouble(0));
|
|
group->setReservoirMaxRate( currentStep , record->getItem("RESV_TARGET")->getSIDouble(0));
|
|
group->setTargetReinjectFraction( currentStep , record->getItem("REINJ_TARGET")->getRawDouble(0));
|
|
group->setTargetVoidReplacementFraction( currentStep , record->getItem("VOIDAGE_TARGET")->getRawDouble(0));
|
|
}
|
|
}
|
|
|
|
|
|
void Schedule::handleGCONPROD(DeckKeywordConstPtr keyword, size_t currentStep) {
|
|
for (size_t recordNr = 0; recordNr < keyword->size(); recordNr++) {
|
|
DeckRecordConstPtr record = keyword->getRecord(recordNr);
|
|
const std::string& groupName = record->getItem("GROUP")->getString(0);
|
|
GroupPtr group = getGroup(groupName);
|
|
{
|
|
GroupProduction::ControlEnum controlMode = GroupProduction::ControlEnumFromString( record->getItem("CONTROL_MODE")->getString(0) );
|
|
group->setProductionControlMode( currentStep , controlMode );
|
|
}
|
|
group->setOilTargetRate( currentStep , record->getItem("OIL_TARGET")->getSIDouble(0));
|
|
group->setGasTargetRate( currentStep , record->getItem("GAS_TARGET")->getSIDouble(0));
|
|
group->setWaterTargetRate( currentStep , record->getItem("WATER_TARGET")->getSIDouble(0));
|
|
group->setLiquidTargetRate( currentStep , record->getItem("LIQUID_TARGET")->getSIDouble(0));
|
|
{
|
|
GroupProductionExceedLimit::ActionEnum exceedAction = GroupProductionExceedLimit::ActionEnumFromString(record->getItem("EXCEED_PROC")->getString(0) );
|
|
group->setProductionExceedLimitAction( currentStep , exceedAction );
|
|
}
|
|
|
|
}
|
|
}
|
|
|
|
void Schedule::handleCOMPDAT(DeckKeywordConstPtr keyword , size_t currentStep) {
|
|
std::map<std::string , std::vector< CompletionConstPtr> > completionMapList = Completion::completionsFromCOMPDATKeyword( keyword );
|
|
std::map<std::string , std::vector< CompletionConstPtr> >::iterator iter;
|
|
|
|
for( iter= completionMapList.begin(); iter != completionMapList.end(); iter++) {
|
|
const std::string wellName = iter->first;
|
|
WellPtr well = getWell(wellName);
|
|
well->addCompletions(currentStep, iter->second);
|
|
}
|
|
}
|
|
|
|
void Schedule::handleGRUPTREE(DeckKeywordConstPtr keyword, size_t currentStep) {
|
|
GroupTreePtr currentTree = m_rootGroupTree->get(currentStep);
|
|
GroupTreePtr newTree = currentTree->deepCopy();
|
|
for (size_t recordNr = 0; recordNr < keyword->size(); recordNr++) {
|
|
DeckRecordConstPtr record = keyword->getRecord(recordNr);
|
|
const std::string& childName = record->getItem("CHILD_GROUP")->getString(0);
|
|
const std::string& parentName = record->getItem("PARENT_GROUP")->getString(0);
|
|
newTree->updateTree(childName, parentName);
|
|
}
|
|
m_rootGroupTree->add(currentStep, newTree);
|
|
}
|
|
|
|
boost::gregorian::date Schedule::getStartDate() const {
|
|
return m_timeMap->getStartDate();
|
|
}
|
|
|
|
TimeMapConstPtr Schedule::getTimeMap() const {
|
|
return m_timeMap;
|
|
}
|
|
|
|
GroupTreePtr Schedule::getGroupTree(size_t timeStep) const {
|
|
return m_rootGroupTree->get(timeStep);
|
|
}
|
|
|
|
void Schedule::addWell(const std::string& wellName, DeckRecordConstPtr record, size_t timeStep) {
|
|
int headI = record->getItem("HEAD_I")->getInt(0);
|
|
int headJ = record->getItem("HEAD_J")->getInt(0);
|
|
double refDepth = record->getItem("REF_DEPTH")->getSIDouble(0);
|
|
WellPtr well(new Well(wellName, headI, headJ, refDepth, m_timeMap , timeStep));
|
|
m_wells[ wellName ] = well;
|
|
}
|
|
|
|
size_t Schedule::numWells() const {
|
|
return m_wells.size();
|
|
}
|
|
|
|
bool Schedule::hasWell(const std::string& wellName) const {
|
|
return m_wells.find(wellName) != m_wells.end();
|
|
}
|
|
|
|
WellPtr Schedule::getWell(const std::string& wellName) const {
|
|
if (hasWell(wellName)) {
|
|
return m_wells.at(wellName);
|
|
} else
|
|
throw std::invalid_argument("Well: " + wellName + " does not exist");
|
|
}
|
|
|
|
std::vector<WellConstPtr> Schedule::getWells() const {
|
|
return getWells(m_timeMap->size()-1);
|
|
}
|
|
|
|
std::vector<WellConstPtr> Schedule::getWells(size_t timeStep) const {
|
|
if (timeStep >= m_timeMap->size()) {
|
|
throw std::invalid_argument("Timestep to large");
|
|
}
|
|
|
|
std::vector<WellConstPtr> wells;
|
|
for (auto iter = m_wells.begin(); iter != m_wells.end(); ++iter) {
|
|
WellConstPtr well = (*iter).second;
|
|
if (well->hasBeenDefined(timeStep)) {
|
|
wells.push_back(well);
|
|
}
|
|
}
|
|
return wells;
|
|
}
|
|
|
|
void Schedule::addGroup(const std::string& groupName, size_t timeStep) {
|
|
if (!m_timeMap) {
|
|
throw std::invalid_argument("TimeMap is null, can't add group named: " + groupName);
|
|
}
|
|
GroupPtr group(new Group(groupName, m_timeMap , timeStep));
|
|
m_groups[ groupName ] = group;
|
|
}
|
|
|
|
size_t Schedule::numGroups() const {
|
|
return m_groups.size();
|
|
}
|
|
|
|
bool Schedule::hasGroup(const std::string& groupName) const {
|
|
return m_groups.find(groupName) != m_groups.end();
|
|
}
|
|
|
|
GroupPtr Schedule::getGroup(const std::string& groupName) const {
|
|
if (hasGroup(groupName)) {
|
|
return m_groups.at(groupName);
|
|
} else
|
|
throw std::invalid_argument("Group: " + groupName + " does not exist");
|
|
}
|
|
|
|
void Schedule::addWellToGroup( GroupPtr newGroup , WellPtr well , size_t timeStep) {
|
|
const std::string currentGroupName = well->getGroupName(timeStep);
|
|
if (currentGroupName != "") {
|
|
GroupPtr currentGroup = getGroup( currentGroupName );
|
|
currentGroup->delWell( timeStep , well->name());
|
|
}
|
|
well->setGroupName(timeStep , newGroup->name());
|
|
newGroup->addWell(timeStep , well);
|
|
}
|
|
|
|
|
|
}
|