Carter-Tracy aquifer object to collect items from keywords AQUCT, AQUTAB and AQUANCON

2018-01-05 10:56:28 +01:00
parent 2b4fa833ea
commit 1edd7cbd92
3 changed files with 241 additions and 0 deletions
--- a/opm/parser/eclipse/EclipseState/AquiferCT.hpp
+++ b/opm/parser/eclipse/EclipseState/AquiferCT.hpp
@@ -0,0 +1,113 @@
+/*
+  Copyright (C) 2017 TNO
+
+  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/>.
+ */
+
+#ifndef OPM_AQUIFERCT_HPP
+#define OPM_AQUIFERCT_HPP
+
+/*
+  The AquiferCT which stands for AquiferCarterTracy is a data container object meant to hold the data for the aquifer carter tracy model.
+  This includes the logic for parsing as well as the associated tables. It is meant to be used by opm-grid and opm-simulators in order to
+  implement the Carter Tracy analytical aquifer model in OPM Flow.
+*/
+
+#include <opm/parser/eclipse/EclipseState/EclipseState.hpp>
+#include <opm/parser/eclipse/Parser/ParserKeywords/A.hpp>
+#include <opm/parser/eclipse/Deck/Deck.hpp>
+#include <opm/parser/eclipse/Deck/DeckItem.hpp>
+#include <opm/parser/eclipse/Deck/DeckRecord.hpp>
+#include <opm/parser/eclipse/Deck/DeckKeyword.hpp>
+#include <opm/parser/eclipse/EclipseState/Tables/Aqudims.hpp>
+#include <opm/parser/eclipse/EclipseState/Tables/TableContainer.hpp>
+#include <opm/parser/eclipse/EclipseState/Tables/AqutabTable.hpp>
+#include <boost/concept_check.hpp>
+
+namespace Opm {
+
+    class AquiferCT {
+        public:
+
+            struct AQUCT_data{
+
+                    // Aquifer ID
+                    int aquiferID;
+                    // Table IDs
+                    int inftableID, pvttableID;
+                    std::vector<int> cell_id;
+                    // Variables constants
+                    double  phi_aq , //aquifer porosity
+                            d0,   //aquifer datum depth
+                            C_t , //total compressibility
+                            r_o , //aquifer inner radius
+                            k_a , //aquifer permeability
+                            c1, // 0.008527 (METRIC, PVT-M); 0.006328 (FIELD); 3.6 (LAB)
+                            h , //aquifer thickness
+                            theta , //angle subtended by the aquifer boundary
+                            c2 ; //6.283 (METRIC, PVT-M); 1.1191 (FIELD); 6.283 (LAB).
+
+                    std::vector<double> td, pi;
+            };
+
+            struct AQUANCON_data{
+
+                    // Aquifer ID
+                    int aquiferID;
+                    // Grid cell box definition to connect aquifer
+                    int i1, i2, j1, j2, k1, k2;
+                    // Variables constants
+                    double  influx_coeff,  //Aquifer influx coefficient
+                            influx_mult;   //Aquifer influx coefficient Multiplier       
+                    // Cell face to connect aquifer to        
+                    std::string face;           
+
+            };
+
+            AquiferCT(const EclipseState& eclState, const Deck& deck);
+
+            const std::vector<AquiferCT::AQUCT_data>& getAquifers() const;
+            const std::vector<AquiferCT::AQUANCON_data>& getAquancon() const;
+            const int getAqInflTabID(size_t aquiferIndex);
+            const int getAqPvtTabID(size_t aquiferIndex);
+            const double getAqInfluxCoeff(size_t aquanconRecord);
+            const double getAqInfluxMult(size_t aquanconRecord);
+
+            std::vector<AquiferCT::AQUCT_data> m_aquct;
+            std::vector<AquiferCT::AQUANCON_data> m_aquancon;
+    
+        private:
+
+            inline void set_default_tables(std::vector<double>& td, std::vector<double>& pi)
+            {
+                std::vector<double> default_pressure_ = { 0.112, 0.229, 0.315, 0.376, 0.424, 0.469, 0.503, 0.564, 0.616, 0.659, 0.702, 0.735, 
+                                                          0.772, 0.802, 0.927, 1.02, 1.101, 1.169, 1.275, 1.362, 1.436, 1.5, 1.556, 1.604, 
+                                                          1.651, 1.829, 1.96, 2.067, 2.147, 2.282, 2.388, 2.476, 2.55, 2.615, 2.672, 2.723, 
+                                                          2.921, 3.064, 3.173, 3.263, 3.406, 3.516, 3.608, 3.684, 3.75, 3.809, 3.86 };
+
+                std::vector<double> default_time_ = { 0.01, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 
+                                                        1.5, 2, 2.5, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 25, 30, 40, 50, 60, 70, 
+                                                        80, 90, 100, 150, 200, 250, 300, 400, 500, 600, 700, 800, 900, 1000 };
+
+                td = default_time_;
+                pi = default_pressure_;
+            }
+
+    };
+}
+
+
+#endif