129 lines
4.4 KiB
C++
129 lines
4.4 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/Completion.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/Schedule/ScheduleEnums.hpp>
|
|
|
|
namespace Opm {
|
|
|
|
Completion::Completion(int i, int j , int k , CompletionStateEnum state , double CF) {
|
|
m_i = i;
|
|
m_j = j;
|
|
m_k = k;
|
|
|
|
m_state = state;
|
|
m_CF = CF;
|
|
}
|
|
|
|
|
|
bool Completion::sameCoordinate(const Completion& other) const {
|
|
if ((m_i == other.m_i) &&
|
|
(m_j == other.m_j) &&
|
|
(m_k == other.m_k))
|
|
return true;
|
|
else
|
|
return false;
|
|
}
|
|
|
|
/**
|
|
This will break up one record and return a pair: <name ,
|
|
[Completion1, Completion2, ... , CompletionN]>. The reason it
|
|
will return a list is that the 'K1 K2' structure is
|
|
disentangled, and each completion is returned separately.
|
|
*/
|
|
|
|
std::pair<std::string , std::vector<CompletionConstPtr> > Completion::completionsFromCOMPDATRecord( DeckRecordConstPtr compdatRecord ) {
|
|
std::vector<CompletionConstPtr> completions;
|
|
std::string well = compdatRecord->getItem("WELL")->getString(0);
|
|
int I = compdatRecord->getItem("I")->getInt(0);
|
|
int J = compdatRecord->getItem("J")->getInt(0);
|
|
int K1 = compdatRecord->getItem("K1")->getInt(0);
|
|
int K2 = compdatRecord->getItem("K2")->getInt(0);
|
|
CompletionStateEnum state = CompletionStateEnumFromString( compdatRecord->getItem("STATE")->getString(0) );
|
|
{
|
|
DeckItemConstPtr CFItem = compdatRecord->getItem("CF");
|
|
if (CFItem->defaultApplied())
|
|
throw std::invalid_argument("The connection factor item can not be defaulted");
|
|
}
|
|
double CF = compdatRecord->getItem("CF")->getRawDouble(0);
|
|
|
|
for (int k = K1; k <= K2; k++) {
|
|
CompletionConstPtr completion(new Completion(I , J , k , state , CF ));
|
|
completions.push_back( completion );
|
|
}
|
|
|
|
return std::pair<std::string , std::vector<CompletionConstPtr> >( well , completions );
|
|
}
|
|
|
|
/*
|
|
Will return a map:
|
|
|
|
{
|
|
"WELL1" : [ Completion1 , Completion2 , ... , CompletionN ],
|
|
"WELL2" : [ Completion1 , Completion2 , ... , CompletionN ],
|
|
...
|
|
}
|
|
*/
|
|
|
|
|
|
std::map<std::string , std::vector< CompletionConstPtr> > Completion::completionsFromCOMPDATKeyword( DeckKeywordConstPtr compdatKeyword ) {
|
|
std::map<std::string , std::vector< CompletionConstPtr> > completionMapList;
|
|
for (size_t recordIndex = 0; recordIndex < compdatKeyword->size(); recordIndex++) {
|
|
std::pair<std::string , std::vector< CompletionConstPtr> > wellCompletionsPair = completionsFromCOMPDATRecord( compdatKeyword->getRecord( recordIndex ));
|
|
std::string well = wellCompletionsPair.first;
|
|
std::vector<CompletionConstPtr>& newCompletions = wellCompletionsPair.second;
|
|
|
|
if (completionMapList.find(well) == completionMapList.end())
|
|
completionMapList[well] = std::vector<CompletionConstPtr>();
|
|
|
|
{
|
|
std::vector<CompletionConstPtr>& currentCompletions = completionMapList.find(well)->second;
|
|
|
|
for (size_t ic = 0; ic < newCompletions.size(); ic++)
|
|
currentCompletions.push_back( newCompletions[ic] );
|
|
}
|
|
}
|
|
return completionMapList;
|
|
}
|
|
|
|
|
|
int Completion::getI() const {
|
|
return m_i;
|
|
}
|
|
|
|
int Completion::getJ() const {
|
|
return m_j;
|
|
}
|
|
|
|
int Completion::getK() const {
|
|
return m_k;
|
|
}
|
|
|
|
CompletionStateEnum Completion::getState() const {
|
|
return m_state;
|
|
};
|
|
|
|
double Completion::getCF() const {
|
|
return m_CF;
|
|
}
|
|
|
|
}
|
|
|
|
|