My Project
Loading...
Searching...
No Matches
EclipseState.hpp
1/*
2 Copyright 2013 Statoil ASA.
3
4 This file is part of the Open Porous Media project (OPM).
5
6 OPM is free software: you can redistribute it and/or modify
7 it under the terms of the GNU General Public License as published by
8 the Free Software Foundation, either version 3 of the License, or
9 (at your option) any later version.
10
11 OPM is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 GNU General Public License for more details.
15
16 You should have received a copy of the GNU General Public License
17 along with OPM. If not, see <http://www.gnu.org/licenses/>.
18*/
19
20#ifndef OPM_ECLIPSE_STATE_HPP
21#define OPM_ECLIPSE_STATE_HPP
22
23#include <opm/input/eclipse/EclipseState/Aquifer/AquiferConfig.hpp>
24#include <opm/input/eclipse/EclipseState/Compositional/CompositionalConfig.hpp>
25#include <opm/input/eclipse/EclipseState/EclipseConfig.hpp>
26#include <opm/input/eclipse/EclipseState/Grid/EclipseGrid.hpp>
27#include <opm/input/eclipse/EclipseState/Grid/FIPRegionStatistics.hpp>
28#include <opm/input/eclipse/EclipseState/Grid/FaultCollection.hpp>
29#include <opm/input/eclipse/EclipseState/Grid/FieldPropsManager.hpp>
30#include <opm/input/eclipse/EclipseState/Grid/LgrCollection.hpp>
31#include <opm/input/eclipse/EclipseState/Grid/NNC.hpp>
32#include <opm/input/eclipse/EclipseState/Grid/TransMult.hpp>
33#include <opm/input/eclipse/EclipseState/MICPpara.hpp>
34#include <opm/input/eclipse/EclipseState/Runspec.hpp>
35#include <opm/input/eclipse/EclipseState/SimulationConfig/SimulationConfig.hpp>
36#include <opm/input/eclipse/EclipseState/Tables/TableManager.hpp>
37#include <opm/input/eclipse/EclipseState/TracerConfig.hpp>
38#include <opm/input/eclipse/EclipseState/Co2StoreConfig.hpp>
39#include <opm/input/eclipse/EclipseState/WagHysteresisConfig.hpp>
40
41#include <opm/input/eclipse/Units/UnitSystem.hpp>
42
43#include <cstddef>
44#include <memory>
45#include <optional>
46#include <vector>
47
48namespace Opm {
49 class Deck;
50 class DeckKeyword;
51 class InitConfig;
52 class IOConfig;
53 class DeckSection;
54} // namespace Opm
55
56namespace Opm { namespace RestartIO {
57 class RstAquifer;
58 class RstNetwork;
59}} // namespace Opm::RestartIO
60
61namespace Opm {
62
64 public:
65 enum EnabledTypes {
66 IntProperties = 0x01,
67 DoubleProperties = 0x02,
68
69 AllProperties = IntProperties | DoubleProperties
70 };
71
72 EclipseState() = default;
73 explicit EclipseState(const Deck& deck);
74 virtual ~EclipseState() = default;
75
76 const IOConfig& getIOConfig() const;
78
79 const InitConfig& getInitConfig() const;
81
82 const SimulationConfig& getSimulationConfig() const;
83 virtual const EclipseGrid& getInputGrid() const;
84
85 const FaultCollection& getFaults() const;
86 const TransMult& getTransMult() const;
87 TransMult& getTransMult();
88
91 const NNC& getInputNNC() const;
92 void appendInputNNC(const std::vector<NNCdata>& nnc);
93 void setInputNNC(const NNC& nnc);
94 bool hasInputNNC() const;
95
96 // The potentially parallelized field properties
97 virtual const FieldPropsManager& fieldProps() const;
98 // Always the non-parallel field properties
99 virtual const FieldPropsManager& globalFieldProps() const;
100
101 virtual void computeFipRegionStatistics();
102 const FIPRegionStatistics& fipRegionStatistics() const;
103
104 const TableManager& getTableManager() const;
105 const EclipseConfig& getEclipseConfig() const;
106 const EclipseConfig& cfg() const;
107 const GridDims& gridDims() const;
108
109 const LgrCollection& getLgrs() const;
110 bool hasInputLGR() const;
111
112 // the unit system used by the deck. note that it is rarely needed
113 // to convert units because internally to opm-parser everything is
114 // represented by SI units.
115 const UnitSystem& getDeckUnitSystem() const;
116 const UnitSystem& getUnits() const;
117
118 std::string getTitle() const;
119
120 void apply_schedule_keywords(const std::vector<DeckKeyword>& keywords);
121
122 const Runspec& runspec() const;
123 const AquiferConfig& aquifer() const;
124 const TracerConfig& tracer() const;
125 const MICPpara& getMICPpara() const;
126 const WagHysteresisConfig& getWagHysteresis() const;
127 const Co2StoreConfig& getCo2StoreConfig() const;
128
129 void reset_actnum(const std::vector<int>& new_actnum);
130 void pruneDeactivatedAquiferConnections(const std::vector<std::size_t>& deactivated_cells);
131 void loadRestartAquifers(const RestartIO::RstAquifer& aquifers);
132 // TODO: it is possible that the aquifer are opened through SCHEDULE and not specified in the SOLUTION section
133 // For the ease of the implementation, we create inactive aquifer in the AquiferConfig.
134 // At the moment, it only works for Constant Flux Aquifers(AQUFLUX) as indicated by the function name
135 // When we know and decide to handle the same for AQUFETP and AQUCT, this part will be refactored
136 void appendAqufluxSchedule(const std::unordered_set<int>& ids);
137
138 void loadRestartNetworkPressures(const RestartIO::RstNetwork& network);
139 const std::optional<std::map<std::string, double> >& getRestartNetworkPressures() const { return this->m_restart_network_pressures; }
140
141 template<class Serializer>
142 void serializeOp(Serializer& serializer)
143 {
144 // FieldPropsManager is handled through a different mechanism.
145 // Do not add the member (i.e., field_props) to this list.
146 serializer(m_tables);
147 serializer(m_runspec);
148 serializer(m_eclipseConfig);
149 serializer(m_deckUnitSystem);
150 serializer(m_inputNnc);
151 serializer(m_gridDims);
152 serializer(m_lgrs);
153 serializer(m_simulationConfig);
154 serializer(aquifer_config);
155 serializer(compositional_config),
156 serializer(m_transMult);
157 serializer(m_faults);
158 serializer(m_title);
159 serializer(tracer_config);
160 serializer(m_micppara);
161 serializer(wag_hyst_config);
162 serializer(co2_store_config);
163 serializer(this->fipRegionStatistics_);
164 }
165
166 static bool rst_cmp(const EclipseState& full_state, const EclipseState& rst_state);
167
168 private:
169 void initIOConfigPostSchedule(const Deck& deck);
170 void assignRunTitle(const Deck& deck);
171 void reportNumberOfActivePhases() const;
172 void initLgrs(const Deck& deck);
173 void conveyNumericalAquiferEffects();
174 void applyMULTXYZ();
175 void initFaults(const Deck& deck);
176 void initPara(const Deck& deck);
177
178 void setMULTFLT(const Opm::DeckSection& section, bool edit = false);
179
180 void complainAboutAmbiguousKeyword(const Deck& deck,
181 const std::string& keywordName);
182
183 protected:
184 TableManager m_tables;
185 Runspec m_runspec;
186 EclipseConfig m_eclipseConfig;
187 UnitSystem m_deckUnitSystem;
188 EclipseGrid m_inputGrid;
189 NNC m_inputNnc;
190 GridDims m_gridDims;
191 FieldPropsManager field_props;
192 LgrCollection m_lgrs;
193 SimulationConfig m_simulationConfig;
194 AquiferConfig aquifer_config;
195 CompositionalConfig compositional_config;
196 TransMult m_transMult;
197 TracerConfig tracer_config;
198 MICPpara m_micppara;
199 WagHysteresisConfig wag_hyst_config;
200 Co2StoreConfig co2_store_config;
201
202 std::string m_title{};
203 FaultCollection m_faults{};
204
205 std::optional<std::map<std::string, double> > m_restart_network_pressures{std::nullopt};
206
207 std::optional<FIPRegionStatistics> fipRegionStatistics_{std::nullopt};
208 };
209} // namespace Opm
210
211#endif // OPM_ECLIPSE_STATE_HPP
Definition AquiferConfig.hpp:46
Definition Co2StoreConfig.hpp:30
Definition CompositionalConfig.hpp:34
Definition DeckSection.hpp:47
Definition Deck.hpp:49
Definition EclipseConfig.hpp:32
About cell information and dimension: The actual grid information is held in a pointer to an ERT ecl_...
Definition EclipseGrid.hpp:55
Definition EclipseState.hpp:63
const EclipseConfig & getEclipseConfig() const
[[deprecated]] — use cfg()
Definition EclipseState.cpp:257
const IOConfig & getIOConfig() const
[[deprecated]] — use cfg().io()
Definition EclipseState.cpp:238
const NNC & getInputNNC() const
non-neighboring connections the non-standard adjacencies as specified in input deck
Definition EclipseState.cpp:301
const InitConfig & getInitConfig() const
[[deprecated]] — use cfg().init()
Definition EclipseState.cpp:248
Basic descriptive statistics about a model's fluid-in-place regions.
Definition FIPRegionStatistics.hpp:40
Definition FaultCollection.hpp:35
Definition FieldPropsManager.hpp:42
Definition GridDims.hpp:31
Definition IOConfig.hpp:147
Definition InitConfig.hpp:32
Definition LgrCollection.hpp:33
Definition MICPpara.hpp:26
Definition NNC.hpp:102
Definition aquifer.hpp:45
Definition network.hpp:39
Definition Runspec.hpp:481
Class for (de-)serializing.
Definition Serializer.hpp:84
Definition SimulationConfig.hpp:38
Definition TableManager.hpp:66
Definition TracerConfig.hpp:33
Definition TransMult.hpp:52
Definition UnitSystem.hpp:34
Definition WagHysteresisConfig.hpp:30
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30