My Project
ConstantCompressibilityOilPvt.hpp
Go to the documentation of this file.
1 // -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2 // vi: set et ts=4 sw=4 sts=4:
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 2 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  Consult the COPYING file in the top-level source directory of this
20  module for the precise wording of the license and the list of
21  copyright holders.
22 */
27 #ifndef OPM_CONSTANT_COMPRESSIBILITY_OIL_PVT_HPP
28 #define OPM_CONSTANT_COMPRESSIBILITY_OIL_PVT_HPP
29 
30 #if HAVE_ECL_INPUT
31 #include <opm/input/eclipse/EclipseState/EclipseState.hpp>
32 #include <opm/input/eclipse/Schedule/Schedule.hpp>
33 #endif
34 
35 namespace Opm {
40 template <class Scalar>
42 {
43 public:
45  ConstantCompressibilityOilPvt(const std::vector<Scalar>& oilReferenceDensity,
46  const std::vector<Scalar>& oilReferencePressure,
47  const std::vector<Scalar>& oilReferenceFormationVolumeFactor,
48  const std::vector<Scalar>& oilCompressibility,
49  const std::vector<Scalar>& oilViscosity,
50  const std::vector<Scalar>& oilViscosibility)
51  : oilReferenceDensity_(oilReferenceDensity)
52  , oilReferencePressure_(oilReferencePressure)
53  , oilReferenceFormationVolumeFactor_(oilReferenceFormationVolumeFactor)
54  , oilCompressibility_(oilCompressibility)
55  , oilViscosity_(oilViscosity)
56  , oilViscosibility_(oilViscosibility)
57  { }
58 
59 #if HAVE_ECL_INPUT
60 
68  void initFromState(const EclipseState& eclState, const Schedule&)
69  {
70  const auto& pvcdoTable = eclState.getTableManager().getPvcdoTable();
71  const auto& densityTable = eclState.getTableManager().getDensityTable();
72 
73  assert(pvcdoTable.size() == densityTable.size());
74 
75  size_t numRegions = pvcdoTable.size();
76  setNumRegions(numRegions);
77 
78  for (unsigned regionIdx = 0; regionIdx < numRegions; ++ regionIdx) {
79  Scalar rhoRefO = densityTable[regionIdx].oil;
80  Scalar rhoRefG = densityTable[regionIdx].gas;
81  Scalar rhoRefW = densityTable[regionIdx].water;
82 
83  setReferenceDensities(regionIdx, rhoRefO, rhoRefG, rhoRefW);
84 
85  oilReferencePressure_[regionIdx] = pvcdoTable[regionIdx].reference_pressure;
86  oilReferenceFormationVolumeFactor_[regionIdx] = pvcdoTable[regionIdx].volume_factor;
87  oilCompressibility_[regionIdx] = pvcdoTable[regionIdx].compressibility;
88  oilViscosity_[regionIdx] = pvcdoTable[regionIdx].viscosity;
89  oilViscosibility_[regionIdx] = pvcdoTable[regionIdx].viscosibility;
90  }
91 
92  initEnd();
93  }
94 #endif
95 
96  void setNumRegions(size_t numRegions)
97  {
98  oilReferenceDensity_.resize(numRegions);
99  oilReferencePressure_.resize(numRegions);
100  oilReferenceFormationVolumeFactor_.resize(numRegions);
101  oilCompressibility_.resize(numRegions);
102  oilViscosity_.resize(numRegions);
103  oilViscosibility_.resize(numRegions);
104 
105  for (unsigned regionIdx = 0; regionIdx < numRegions; ++regionIdx) {
106  setReferenceFormationVolumeFactor(regionIdx, 1.0);
107  setReferencePressure(regionIdx, 1.03125);
108  }
109  }
110 
114  void setReferenceDensities(unsigned regionIdx,
115  Scalar rhoRefOil,
116  Scalar /*rhoRefGas*/,
117  Scalar /*rhoRefWater*/)
118  { oilReferenceDensity_[regionIdx] = rhoRefOil; }
119 
123  void setViscosity(unsigned regionIdx, Scalar muo, Scalar oilViscosibility = 0.0)
124  {
125  oilViscosity_[regionIdx] = muo;
126  oilViscosibility_[regionIdx] = oilViscosibility;
127  }
128 
132  void setCompressibility(unsigned regionIdx, Scalar oilCompressibility)
133  { oilCompressibility_[regionIdx] = oilCompressibility; }
134 
138  void setReferencePressure(unsigned regionIdx, Scalar p)
139  { oilReferencePressure_[regionIdx] = p; }
140 
144  void setReferenceFormationVolumeFactor(unsigned regionIdx, Scalar BoRef)
145  { oilReferenceFormationVolumeFactor_[regionIdx] = BoRef; }
146 
150  void setViscosibility(unsigned regionIdx, Scalar muComp)
151  { oilViscosibility_[regionIdx] = muComp; }
152 
156  void initEnd()
157  { }
158 
162  unsigned numRegions() const
163  { return oilViscosity_.size(); }
164 
168  template <class Evaluation>
169  Evaluation internalEnergy(unsigned,
170  const Evaluation&,
171  const Evaluation&,
172  const Evaluation&) const
173  {
174  throw std::runtime_error("Requested the enthalpy of oil but the thermal option is not enabled");
175  }
176 
181  template <class Evaluation>
182  Evaluation viscosity(unsigned regionIdx,
183  const Evaluation& temperature,
184  const Evaluation& pressure,
185  const Evaluation& /*Rs*/) const
186  { return saturatedViscosity(regionIdx, temperature, pressure); }
187 
191  template <class Evaluation>
192  Evaluation saturatedViscosity(unsigned regionIdx,
193  const Evaluation& temperature,
194  const Evaluation& pressure) const
195  {
196  Scalar BoMuoRef = oilViscosity_[regionIdx]*oilReferenceFormationVolumeFactor_[regionIdx];
197  const Evaluation& bo = saturatedInverseFormationVolumeFactor(regionIdx, temperature, pressure);
198 
199  Scalar pRef = oilReferencePressure_[regionIdx];
200  const Evaluation& Y =
201  (oilCompressibility_[regionIdx] - oilViscosibility_[regionIdx])
202  * (pressure - pRef);
203  return BoMuoRef*bo/(1.0 + Y*(1.0 + Y/2.0));
204  }
205 
209  template <class Evaluation>
210  Evaluation inverseFormationVolumeFactor(unsigned regionIdx,
211  const Evaluation& temperature,
212  const Evaluation& pressure,
213  const Evaluation& /*Rs*/) const
214  { return saturatedInverseFormationVolumeFactor(regionIdx, temperature, pressure); }
215 
222  template <class Evaluation>
223  Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx,
224  const Evaluation& /*temperature*/,
225  const Evaluation& pressure) const
226  {
227  // cf. ECLiPSE 2011 technical description, p. 116
228  Scalar pRef = oilReferencePressure_[regionIdx];
229  const Evaluation& X = oilCompressibility_[regionIdx]*(pressure - pRef);
230 
231  Scalar BoRef = oilReferenceFormationVolumeFactor_[regionIdx];
232  return (1 + X*(1 + X/2))/BoRef;
233  }
234 
238  template <class Evaluation>
239  Evaluation saturatedGasDissolutionFactor(unsigned /*regionIdx*/,
240  const Evaluation& /*temperature*/,
241  const Evaluation& /*pressure*/) const
242  { return 0.0; /* this is dead oil! */ }
243 
247  template <class Evaluation>
248  Evaluation saturatedGasDissolutionFactor(unsigned /*regionIdx*/,
249  const Evaluation& /*temperature*/,
250  const Evaluation& /*pressure*/,
251  const Evaluation& /*oilSaturation*/,
252  const Evaluation& /*maxOilSaturation*/) const
253  { return 0.0; /* this is dead oil! */ }
254 
261  template <class Evaluation>
262  Evaluation saturationPressure(unsigned /*regionIdx*/,
263  const Evaluation& /*temperature*/,
264  const Evaluation& /*Rs*/) const
265  { return 0.0; /* this is dead oil, so there isn't any meaningful saturation pressure! */ }
266 
267  template <class Evaluation>
268  Evaluation diffusionCoefficient(const Evaluation& /*temperature*/,
269  const Evaluation& /*pressure*/,
270  unsigned /*compIdx*/) const
271  {
272  throw std::runtime_error("Not implemented: The PVT model does not provide a diffusionCoefficient()");
273  }
274 
275  const Scalar oilReferenceDensity(unsigned regionIdx) const
276  { return oilReferenceDensity_[regionIdx]; }
277 
278  const std::vector<Scalar>& oilReferenceFormationVolumeFactor() const
279  { return oilReferenceFormationVolumeFactor_; }
280 
281  const std::vector<Scalar>& oilCompressibility() const
282  { return oilCompressibility_; }
283 
284  const std::vector<Scalar>& oilViscosity() const
285  { return oilViscosity_; }
286 
287  const std::vector<Scalar>& oilViscosibility() const
288  { return oilViscosibility_; }
289 
290  bool operator==(const ConstantCompressibilityOilPvt<Scalar>& data) const
291  {
292  return this->oilReferenceDensity_ == data.oilReferenceDensity_ &&
293  this->oilReferencePressure_ == data.oilReferencePressure_ &&
294  this->oilReferenceFormationVolumeFactor() == data.oilReferenceFormationVolumeFactor() &&
295  this->oilCompressibility() == data.oilCompressibility() &&
296  this->oilViscosity() == data.oilViscosity() &&
297  this->oilViscosibility() == data.oilViscosibility();
298  }
299 
300 private:
301  std::vector<Scalar> oilReferenceDensity_;
302  std::vector<Scalar> oilReferencePressure_;
303  std::vector<Scalar> oilReferenceFormationVolumeFactor_;
304  std::vector<Scalar> oilCompressibility_;
305  std::vector<Scalar> oilViscosity_;
306  std::vector<Scalar> oilViscosibility_;
307 };
308 
309 } // namespace Opm
310 
311 #endif
This class represents the Pressure-Volume-Temperature relations of the oil phase without dissolved ga...
Definition: ConstantCompressibilityOilPvt.hpp:42
Evaluation internalEnergy(unsigned, const Evaluation &, const Evaluation &, const Evaluation &) const
Returns the specific enthalpy [J/kg] of oil given a set of parameters.
Definition: ConstantCompressibilityOilPvt.hpp:169
Evaluation saturationPressure(unsigned, const Evaluation &, const Evaluation &) const
Returns the saturation pressure of the oil phase [Pa] depending on its mass fraction of the gas compo...
Definition: ConstantCompressibilityOilPvt.hpp:262
void setViscosity(unsigned regionIdx, Scalar muo, Scalar oilViscosibility=0.0)
Set the viscosity and "viscosibility" of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:123
void initEnd()
Finish initializing the oil phase PVT properties.
Definition: ConstantCompressibilityOilPvt.hpp:156
Evaluation saturatedGasDissolutionFactor(unsigned, const Evaluation &, const Evaluation &, const Evaluation &, const Evaluation &) const
Returns the gas dissolution factor [m^3/m^3] of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:248
void setReferencePressure(unsigned regionIdx, Scalar p)
Set the oil reference pressure [Pa].
Definition: ConstantCompressibilityOilPvt.hpp:138
void setCompressibility(unsigned regionIdx, Scalar oilCompressibility)
Set the compressibility of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:132
Evaluation viscosity(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &) const
Returns the dynamic viscosity [Pa s] of gas saturated oil given a pressure and a phase composition.
Definition: ConstantCompressibilityOilPvt.hpp:182
void setReferenceDensities(unsigned regionIdx, Scalar rhoRefOil, Scalar, Scalar)
Initialize the reference densities of all fluids for a given PVT region.
Definition: ConstantCompressibilityOilPvt.hpp:114
void setViscosibility(unsigned regionIdx, Scalar muComp)
Set the oil "viscosibility" [1/ (Pa s)].
Definition: ConstantCompressibilityOilPvt.hpp:150
Evaluation inverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &) const
Returns the formation volume factor [-] of the fluid phase.
Definition: ConstantCompressibilityOilPvt.hpp:210
Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the formation volume factor [-] of gas saturated oil.
Definition: ConstantCompressibilityOilPvt.hpp:223
void setReferenceFormationVolumeFactor(unsigned regionIdx, Scalar BoRef)
Set the oil reference formation volume factor [-].
Definition: ConstantCompressibilityOilPvt.hpp:144
Evaluation saturatedGasDissolutionFactor(unsigned, const Evaluation &, const Evaluation &) const
Returns the gas dissolution factor [m^3/m^3] of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:239
Evaluation saturatedViscosity(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure) const
Returns the dynamic viscosity [Pa s] of gas saturated oil given a pressure.
Definition: ConstantCompressibilityOilPvt.hpp:192
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition: ConstantCompressibilityOilPvt.hpp:162