21 #ifndef OPM_TABLE_MANAGER_HPP
22 #define OPM_TABLE_MANAGER_HPP
28 #include <opm/input/eclipse/EclipseState/Tables/DenT.hpp>
29 #include <opm/input/eclipse/EclipseState/Tables/PvtgTable.hpp>
30 #include <opm/input/eclipse/EclipseState/Tables/PvtgwTable.hpp>
31 #include <opm/input/eclipse/EclipseState/Tables/PvtgwoTable.hpp>
32 #include <opm/input/eclipse/EclipseState/Tables/PvtoTable.hpp>
33 #include <opm/input/eclipse/EclipseState/Tables/PvtsolTable.hpp>
34 #include <opm/input/eclipse/EclipseState/Tables/RocktabTable.hpp>
35 #include <opm/input/eclipse/EclipseState/Tables/Rock2dTable.hpp>
36 #include <opm/input/eclipse/EclipseState/Tables/Rock2dtrTable.hpp>
37 #include <opm/input/eclipse/EclipseState/Tables/PlyshlogTable.hpp>
38 #include <opm/input/eclipse/EclipseState/Tables/PvtwsaltTable.hpp>
39 #include <opm/input/eclipse/EclipseState/Tables/RwgsaltTable.hpp>
40 #include <opm/input/eclipse/EclipseState/Tables/BrineDensityTable.hpp>
41 #include <opm/input/eclipse/EclipseState/Tables/SolventDensityTable.hpp>
42 #include <opm/input/eclipse/EclipseState/Tables/StandardCond.hpp>
43 #include <opm/input/eclipse/EclipseState/Tables/FlatTable.hpp>
44 #include <opm/input/eclipse/EclipseState/Tables/SorwmisTable.hpp>
45 #include <opm/input/eclipse/EclipseState/Tables/SgcwmisTable.hpp>
46 #include <opm/input/eclipse/EclipseState/Tables/MiscTable.hpp>
47 #include <opm/input/eclipse/EclipseState/Tables/PmiscTable.hpp>
48 #include <opm/input/eclipse/EclipseState/Tables/MsfnTable.hpp>
49 #include <opm/input/eclipse/EclipseState/Tables/JFunc.hpp>
50 #include <opm/input/eclipse/EclipseState/Tables/Tabdims.hpp>
51 #include <opm/input/eclipse/EclipseState/Tables/TableContainer.hpp>
52 #include <opm/input/eclipse/EclipseState/Tables/Aqudims.hpp>
53 #include <opm/input/eclipse/EclipseState/Tables/PlymwinjTable.hpp>
54 #include <opm/input/eclipse/EclipseState/Tables/SkprwatTable.hpp>
55 #include <opm/input/eclipse/EclipseState/Tables/SkprpolyTable.hpp>
56 #include <opm/input/eclipse/EclipseState/Tables/Eqldims.hpp>
57 #include <opm/input/eclipse/EclipseState/Tables/Regdims.hpp>
58 #include <opm/input/eclipse/EclipseState/Tables/TLMixpar.hpp>
71 const TableContainer& getTables(
const std::string& tableName )
const;
72 const TableContainer& operator[](
const std::string& tableName)
const;
73 bool hasTables(
const std::string& tableName )
const;
75 const Tabdims& getTabdims()
const;
76 const Eqldims& getEqldims()
const;
77 const Aqudims& getAqudims()
const;
78 const Regdims& getRegdims()
const;
83 size_t numFIPRegions()
const;
133 const JFunc& getJFunc()
const;
135 const std::vector<PvtgTable>& getPvtgTables()
const;
136 const std::vector<PvtgwTable>& getPvtgwTables()
const;
137 const std::vector<PvtgwoTable>& getPvtgwoTables()
const;
138 const std::vector<PvtoTable>& getPvtoTables()
const;
139 const std::vector<PvtsolTable>& getPvtsolTables()
const;
140 const std::vector<Rock2dTable>& getRock2dTables()
const;
141 const std::vector<Rock2dtrTable>& getRock2dtrTables()
const;
145 const DenT& WatDenT()
const;
146 const DenT& GasDenT()
const;
147 const DenT& OilDenT()
const;
149 std::size_t gas_comp_index()
const;
151 const std::vector<PvtwsaltTable>& getPvtwSaltTables()
const;
152 const std::vector<RwgsaltTable>& getRwgSaltTables()
const;
153 const std::vector<BrineDensityTable>& getBrineDensityTables()
const;
154 const std::vector<SolventDensityTable>& getSolventDensityTables()
const;
166 const std::map<int, PlymwinjTable>& getPlymwinjTables()
const;
167 const std::map<int, SkprwatTable>& getSkprwatTables()
const;
168 const std::map<int, SkprpolyTable>& getSkprpolyTables()
const;
169 const std::map<std::string, TableContainer>& getSimpleTables()
const;
186 double rtemp()
const;
188 double salinity()
const;
192 template<
class Serializer>
195 auto simpleTables = m_simpleTables;
196 auto split = splitSimpleTable(simpleTables);
197 serializer.map(simpleTables);
198 serializer(split.plyshMax);
199 serializer.map(split.plyshMap);
200 serializer(split.rockMax);
201 serializer.map(split.rockMap);
202 serializer.vector(m_pvtgTables);
203 serializer.vector(m_pvtgwTables);
204 serializer.vector(m_pvtgwoTables);
205 serializer.vector(m_pvtoTables);
206 serializer.vector(m_pvtsolTables);
207 serializer.vector(m_rock2dTables);
208 serializer.vector(m_rock2dtrTables);
209 m_pvtwTable.serializeOp(serializer);
210 m_pvcdoTable.serializeOp(serializer);
211 m_densityTable.serializeOp(serializer);
212 m_diffCoeffTable.serializeOp(serializer);
213 m_plyvmhTable.serializeOp(serializer);
214 m_rockTable.serializeOp(serializer);
215 m_plmixparTable.serializeOp(serializer);
216 m_shrateTable.serializeOp(serializer);
217 m_stone1exTable.serializeOp(serializer);
218 m_viscrefTable.serializeOp(serializer);
219 m_watdentTable.serializeOp(serializer);
220 serializer.vector(m_pvtwsaltTables);
221 serializer.vector(m_rwgsaltTables);
222 serializer.vector(m_bdensityTables);
223 serializer.vector(m_sdensityTables);
224 serializer.map(m_plymwinjTables);
225 serializer.map(m_skprwatTables);
226 serializer.map(m_skprpolyTables);
227 m_tabdims.serializeOp(serializer);
228 m_regdims.serializeOp(serializer);
229 m_eqldims.serializeOp(serializer);
230 m_aqudims.serializeOp(serializer);
231 serializer(hasImptvd);
232 serializer(hasEnptvd);
233 serializer(hasEqlnum);
234 serializer(hasShrate);
236 oilDenT.serializeOp(serializer);
237 gasDenT.serializeOp(serializer);
238 watDenT.serializeOp(serializer);
239 stcond.serializeOp(serializer);
240 serializer(m_gas_comp_index);
242 serializer(m_salinity);
243 m_tlmixpar.serializeOp(serializer);
244 if (!serializer.isSerializing()) {
245 m_simpleTables = simpleTables;
246 if (split.plyshMax > 0) {
248 for (
const auto& it : split.plyshMap) {
249 container.addTable(it.first, it.second);
251 m_simpleTables.insert(std::make_pair(
"PLYSHLOG", container));
253 if (split.rockMax > 0) {
254 TableContainer container(split.rockMax);
255 for (
const auto& it : split.rockMap) {
256 container.addTable(it.first, it.second);
258 m_simpleTables.insert(std::make_pair(
"ROCKTAB", container));
264 TableContainer& forceGetTables(
const std::string& tableName ,
size_t numTables);
266 void complainAboutAmbiguousKeyword(
const Deck& deck,
const std::string& keywordName);
268 void addTables(
const std::string& tableName ,
size_t numTables);
269 void initSimpleTables(
const Deck& deck);
270 void initRTempTables(
const Deck& deck);
271 void initDims(
const Deck& deck);
272 void initRocktabTables(
const Deck& deck);
273 void initGasvisctTables(
const Deck& deck);
275 void initPlymaxTables(
const Deck& deck);
276 void initPlyrockTables(
const Deck& deck);
277 void initPlyshlogTables(
const Deck& deck);
279 void initPlymwinjTables(
const Deck& deck);
280 void initSkprwatTables(
const Deck& deck);
281 void initSkprpolyTables(
const Deck& deck);
285 template <
class TableType>
286 void initRockTables(
const Deck& deck,
const std::string& keywordName, std::vector<TableType>& rocktable );
288 template <
class TableType>
289 void initPvtwsaltTables(
const Deck& deck, std::vector<TableType>& pvtwtables );
291 template <
class TableType>
292 void initRwgsaltTables(
const Deck& deck, std::vector<TableType>& rwgtables );
294 template <
class TableType>
295 void initBrineTables(
const Deck& deck, std::vector<TableType>& brinetables );
297 void initSolventTables(
const Deck& deck, std::vector<SolventDensityTable>& solventtables);
302 template <
class TableType>
303 void initSimpleTableContainerWithJFunc(
const Deck& deck,
304 const std::string& keywordName,
305 const std::string& tableName,
308 template <
class TableType>
309 void initSimpleTableContainer(
const Deck& deck,
310 const std::string& keywordName,
311 const std::string& tableName,
314 template <
class TableType>
315 void initSimpleTableContainer(
const Deck& deck,
316 const std::string& keywordName,
319 template <
class TableType>
320 void initSimpleTableContainerWithJFunc(
const Deck& deck,
321 const std::string& keywordName,
324 template <
class TableType>
325 void initSimpleTable(
const Deck& deck,
326 const std::string& keywordName,
327 std::vector<TableType>& tableVector);
329 template <
class TableType>
330 void initFullTables(
const Deck& deck,
331 const std::string& keywordName,
332 std::vector<TableType>& tableVector);
334 void checkPVTOMonotonicity(
const Deck& deck)
const;
336 void logPVTOMonotonicityFailure(
const Deck& deck,
337 const std::size_t tableID,
338 const std::vector<PvtoTable::FlippedFVF>& flipped_Bo)
const;
340 std::map<std::string , TableContainer> m_simpleTables;
341 std::vector<PvtgTable> m_pvtgTables;
342 std::vector<PvtgwTable> m_pvtgwTables;
343 std::vector<PvtgwoTable> m_pvtgwoTables;
344 std::vector<PvtoTable> m_pvtoTables;
345 std::vector<PvtsolTable> m_pvtsolTables;
346 std::vector<Rock2dTable> m_rock2dTables;
347 std::vector<Rock2dtrTable> m_rock2dtrTables;
348 PvtwTable m_pvtwTable;
349 PvcdoTable m_pvcdoTable;
350 DensityTable m_densityTable;
351 DiffCoeffTable m_diffCoeffTable;
352 PlyvmhTable m_plyvmhTable;
353 RockTable m_rockTable;
354 PlmixparTable m_plmixparTable;
355 ShrateTable m_shrateTable;
356 Stone1exTable m_stone1exTable;
357 ViscrefTable m_viscrefTable;
358 WatdentTable m_watdentTable;
359 std::vector<PvtwsaltTable> m_pvtwsaltTables;
360 std::vector<RwgsaltTable> m_rwgsaltTables;
361 std::vector<BrineDensityTable> m_bdensityTables;
362 std::vector<SolventDensityTable> m_sdensityTables;
363 std::map<int, PlymwinjTable> m_plymwinjTables;
364 std::map<int, SkprwatTable> m_skprwatTables;
365 std::map<int, SkprpolyTable> m_skprpolyTables;
373 bool hasImptvd =
false;
374 bool hasEnptvd =
false;
375 bool hasEqlnum =
false;
376 bool hasShrate =
false;
377 std::optional<JFunc> jfunc;
383 std::size_t m_gas_comp_index = 77;
384 double m_rtemp {288.7056};
385 double m_salinity {0.0};
387 struct SplitSimpleTables {
390 std::map<size_t, std::shared_ptr<PlyshlogTable>> plyshMap;
391 std::map<size_t, std::shared_ptr<RocktabTable>> rockMap;
394 SplitSimpleTables splitSimpleTable(std::map<std::string,TableContainer>& simpleTables);
Definition: Aqudims.hpp:34
Definition: Eqldims.hpp:32
Definition: Regdims.hpp:36
Definition: Serializer.hpp:38
Definition: TLMixpar.hpp:55
Definition: Tabdims.hpp:36
Definition: TableContainer.hpp:31
Definition: TableManager.hpp:64
bool useEqlnum() const
deck has keyword "EQLNUM" — Equilibriation region numbers
bool useJFunc() const
deck has keyword "JFUNC" — Use Leverett's J Function for capillary pressure
bool useEnptvd() const
deck has keyword "ENPTVD" — Saturation end-point versus depth tables
bool useImptvd() const
deck has keyword "IMPTVD" — Imbition end-point versus depth tables
bool useShrate() const
deck has keyword "SHRATE"
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition: Exceptions.hpp:29
Definition: FlatTable.hpp:45
Definition: FlatTable.hpp:122
Definition: FlatTable.hpp:249
Definition: FlatTable.hpp:283
Definition: FlatTable.hpp:224
Definition: FlatTable.hpp:159
Definition: FlatTable.hpp:187
Definition: FlatTable.hpp:308
Definition: StandardCond.hpp:24
Definition: FlatTable.hpp:333
Definition: FlatTable.hpp:389
Definition: FlatTable.hpp:420