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 
34 
35 #if HAVE_ECL_INPUT
36 #include <opm/input/eclipse/EclipseState/EclipseState.hpp>
37 #include <opm/input/eclipse/Schedule/Schedule.hpp>
38 #endif
39 
40 namespace Opm {
45 template <class Scalar>
47 {
49  typedef std::vector<std::pair<Scalar, Scalar> > SamplingPoints;
50 
51 public:
53  ConstantCompressibilityOilPvt(const std::vector<Scalar>& oilReferenceDensity,
54  const std::vector<Scalar>& oilReferencePressure,
55  const std::vector<Scalar>& oilReferenceFormationVolumeFactor,
56  const std::vector<Scalar>& oilCompressibility,
57  const std::vector<Scalar>& oilViscosity,
58  const std::vector<Scalar>& oilViscosibility)
59  : oilReferenceDensity_(oilReferenceDensity)
60  , oilReferencePressure_(oilReferencePressure)
61  , oilReferenceFormationVolumeFactor_(oilReferenceFormationVolumeFactor)
62  , oilCompressibility_(oilCompressibility)
63  , oilViscosity_(oilViscosity)
64  , oilViscosibility_(oilViscosibility)
65  { }
66 
67 #if HAVE_ECL_INPUT
68 
76  void initFromState(const EclipseState& eclState, const Schedule&)
77  {
78  const auto& pvcdoTable = eclState.getTableManager().getPvcdoTable();
79  const auto& densityTable = eclState.getTableManager().getDensityTable();
80 
81  assert(pvcdoTable.size() == densityTable.size());
82 
83  size_t numRegions = pvcdoTable.size();
84  setNumRegions(numRegions);
85 
86  for (unsigned regionIdx = 0; regionIdx < numRegions; ++ regionIdx) {
87  Scalar rhoRefO = densityTable[regionIdx].oil;
88  Scalar rhoRefG = densityTable[regionIdx].gas;
89  Scalar rhoRefW = densityTable[regionIdx].water;
90 
91  setReferenceDensities(regionIdx, rhoRefO, rhoRefG, rhoRefW);
92 
93  oilReferencePressure_[regionIdx] = pvcdoTable[regionIdx].reference_pressure;
94  oilReferenceFormationVolumeFactor_[regionIdx] = pvcdoTable[regionIdx].volume_factor;
95  oilCompressibility_[regionIdx] = pvcdoTable[regionIdx].compressibility;
96  oilViscosity_[regionIdx] = pvcdoTable[regionIdx].viscosity;
97  oilViscosibility_[regionIdx] = pvcdoTable[regionIdx].viscosibility;
98  }
99 
100  initEnd();
101  }
102 #endif
103 
104  void setNumRegions(size_t numRegions)
105  {
106  oilReferenceDensity_.resize(numRegions);
107  oilReferencePressure_.resize(numRegions);
108  oilReferenceFormationVolumeFactor_.resize(numRegions);
109  oilCompressibility_.resize(numRegions);
110  oilViscosity_.resize(numRegions);
111  oilViscosibility_.resize(numRegions);
112 
113  for (unsigned regionIdx = 0; regionIdx < numRegions; ++regionIdx) {
114  setReferenceFormationVolumeFactor(regionIdx, 1.0);
115  setReferencePressure(regionIdx, 1.03125);
116  }
117  }
118 
122  void setReferenceDensities(unsigned regionIdx,
123  Scalar rhoRefOil,
124  Scalar /*rhoRefGas*/,
125  Scalar /*rhoRefWater*/)
126  { oilReferenceDensity_[regionIdx] = rhoRefOil; }
127 
131  void setViscosity(unsigned regionIdx, Scalar muo, Scalar oilViscosibility = 0.0)
132  {
133  oilViscosity_[regionIdx] = muo;
134  oilViscosibility_[regionIdx] = oilViscosibility;
135  }
136 
140  void setCompressibility(unsigned regionIdx, Scalar oilCompressibility)
141  { oilCompressibility_[regionIdx] = oilCompressibility; }
142 
146  void setReferencePressure(unsigned regionIdx, Scalar p)
147  { oilReferencePressure_[regionIdx] = p; }
148 
152  void setReferenceFormationVolumeFactor(unsigned regionIdx, Scalar BoRef)
153  { oilReferenceFormationVolumeFactor_[regionIdx] = BoRef; }
154 
158  void setViscosibility(unsigned regionIdx, Scalar muComp)
159  { oilViscosibility_[regionIdx] = muComp; }
160 
164  void initEnd()
165  { }
166 
170  unsigned numRegions() const
171  { return oilViscosity_.size(); }
172 
176  template <class Evaluation>
177  Evaluation internalEnergy(unsigned,
178  const Evaluation&,
179  const Evaluation&,
180  const Evaluation&) const
181  {
182  throw std::runtime_error("Requested the enthalpy of oil but the thermal option is not enabled");
183  }
184 
189  template <class Evaluation>
190  Evaluation viscosity(unsigned regionIdx,
191  const Evaluation& temperature,
192  const Evaluation& pressure,
193  const Evaluation& /*Rs*/) const
194  { return saturatedViscosity(regionIdx, temperature, pressure); }
195 
199  template <class Evaluation>
200  Evaluation saturatedViscosity(unsigned regionIdx,
201  const Evaluation& temperature,
202  const Evaluation& pressure) const
203  {
204  Scalar BoMuoRef = oilViscosity_[regionIdx]*oilReferenceFormationVolumeFactor_[regionIdx];
205  const Evaluation& bo = saturatedInverseFormationVolumeFactor(regionIdx, temperature, pressure);
206 
207  Scalar pRef = oilReferencePressure_[regionIdx];
208  const Evaluation& Y =
209  (oilCompressibility_[regionIdx] - oilViscosibility_[regionIdx])
210  * (pressure - pRef);
211  return BoMuoRef*bo/(1.0 + Y*(1.0 + Y/2.0));
212  }
213 
217  template <class Evaluation>
218  Evaluation inverseFormationVolumeFactor(unsigned regionIdx,
219  const Evaluation& temperature,
220  const Evaluation& pressure,
221  const Evaluation& /*Rs*/) const
222  { return saturatedInverseFormationVolumeFactor(regionIdx, temperature, pressure); }
223 
230  template <class Evaluation>
231  Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx,
232  const Evaluation& /*temperature*/,
233  const Evaluation& pressure) const
234  {
235  // cf. ECLiPSE 2011 technical description, p. 116
236  Scalar pRef = oilReferencePressure_[regionIdx];
237  const Evaluation& X = oilCompressibility_[regionIdx]*(pressure - pRef);
238 
239  Scalar BoRef = oilReferenceFormationVolumeFactor_[regionIdx];
240  return (1 + X*(1 + X/2))/BoRef;
241  }
242 
246  template <class Evaluation>
247  Evaluation saturatedGasDissolutionFactor(unsigned /*regionIdx*/,
248  const Evaluation& /*temperature*/,
249  const Evaluation& /*pressure*/) const
250  { return 0.0; /* this is dead oil! */ }
251 
255  template <class Evaluation>
256  Evaluation saturatedGasDissolutionFactor(unsigned /*regionIdx*/,
257  const Evaluation& /*temperature*/,
258  const Evaluation& /*pressure*/,
259  const Evaluation& /*oilSaturation*/,
260  const Evaluation& /*maxOilSaturation*/) const
261  { return 0.0; /* this is dead oil! */ }
262 
269  template <class Evaluation>
270  Evaluation saturationPressure(unsigned /*regionIdx*/,
271  const Evaluation& /*temperature*/,
272  const Evaluation& /*Rs*/) const
273  { return 0.0; /* this is dead oil, so there isn't any meaningful saturation pressure! */ }
274 
275  template <class Evaluation>
276  Evaluation diffusionCoefficient(const Evaluation& /*temperature*/,
277  const Evaluation& /*pressure*/,
278  unsigned /*compIdx*/) const
279  {
280  throw std::runtime_error("Not implemented: The PVT model does not provide a diffusionCoefficient()");
281  }
282 
283  const Scalar oilReferenceDensity(unsigned regionIdx) const
284  { return oilReferenceDensity_[regionIdx]; }
285 
286  const std::vector<Scalar>& oilReferenceFormationVolumeFactor() const
287  { return oilReferenceFormationVolumeFactor_; }
288 
289  const std::vector<Scalar>& oilCompressibility() const
290  { return oilCompressibility_; }
291 
292  const std::vector<Scalar>& oilViscosity() const
293  { return oilViscosity_; }
294 
295  const std::vector<Scalar>& oilViscosibility() const
296  { return oilViscosibility_; }
297 
298  bool operator==(const ConstantCompressibilityOilPvt<Scalar>& data) const
299  {
300  return this->oilReferenceDensity_ == data.oilReferenceDensity_ &&
301  this->oilReferencePressure_ == data.oilReferencePressure_ &&
302  this->oilReferenceFormationVolumeFactor() == data.oilReferenceFormationVolumeFactor() &&
303  this->oilCompressibility() == data.oilCompressibility() &&
304  this->oilViscosity() == data.oilViscosity() &&
305  this->oilViscosibility() == data.oilViscosibility();
306  }
307 
308 private:
309  std::vector<Scalar> oilReferenceDensity_;
310  std::vector<Scalar> oilReferencePressure_;
311  std::vector<Scalar> oilReferenceFormationVolumeFactor_;
312  std::vector<Scalar> oilCompressibility_;
313  std::vector<Scalar> oilViscosity_;
314  std::vector<Scalar> oilViscosibility_;
315 };
316 
317 } // namespace Opm
318 
319 #endif
This file provides a wrapper around the "final" C++-2011 statement.
Class implementing cubic splines.
Implements a linearly interpolated scalar function that depends on one variable.
Implements a scalar function that depends on two variables and which is sampled uniformly in the X di...
This class represents the Pressure-Volume-Temperature relations of the oil phase without dissolved ga...
Definition: ConstantCompressibilityOilPvt.hpp:47
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:177
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:270
void setViscosity(unsigned regionIdx, Scalar muo, Scalar oilViscosibility=0.0)
Set the viscosity and "viscosibility" of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:131
void initEnd()
Finish initializing the oil phase PVT properties.
Definition: ConstantCompressibilityOilPvt.hpp:164
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:256
void setReferencePressure(unsigned regionIdx, Scalar p)
Set the oil reference pressure [Pa].
Definition: ConstantCompressibilityOilPvt.hpp:146
void setCompressibility(unsigned regionIdx, Scalar oilCompressibility)
Set the compressibility of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:140
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:190
void setReferenceDensities(unsigned regionIdx, Scalar rhoRefOil, Scalar, Scalar)
Initialize the reference densities of all fluids for a given PVT region.
Definition: ConstantCompressibilityOilPvt.hpp:122
void setViscosibility(unsigned regionIdx, Scalar muComp)
Set the oil "viscosibility" [1/ (Pa s)].
Definition: ConstantCompressibilityOilPvt.hpp:158
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:218
Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the formation volume factor [-] of gas saturated oil.
Definition: ConstantCompressibilityOilPvt.hpp:231
void setReferenceFormationVolumeFactor(unsigned regionIdx, Scalar BoRef)
Set the oil reference formation volume factor [-].
Definition: ConstantCompressibilityOilPvt.hpp:152
Evaluation saturatedGasDissolutionFactor(unsigned, const Evaluation &, const Evaluation &) const
Returns the gas dissolution factor [m^3/m^3] of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:247
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:200
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition: ConstantCompressibilityOilPvt.hpp:170
Implements a linearly interpolated scalar function that depends on one variable.
Definition: Tabulated1DFunction.hpp:47