Classes | Namespaces | Functions
blackoilmodel.hh File Reference

A fully-implicit black-oil flow model. More...

#include <opm/material/localad/Math.hpp>
#include "blackoilproblem.hh"
#include "blackoilindices.hh"
#include "blackoilextensivequantities.hh"
#include "blackoilprimaryvariables.hh"
#include "blackoilintensivequantities.hh"
#include "blackoilratevector.hh"
#include "blackoilboundaryratevector.hh"
#include "blackoillocalresidual.hh"
#include "blackoilnewtonmethod.hh"
#include "blackoilproperties.hh"
#include <ewoms/models/common/multiphasebasemodel.hh>
#include <ewoms/io/vtkcompositionmodule.hh>
#include <ewoms/io/vtkblackoilmodule.hh>
#include <opm/material/fluidsystems/BlackOilFluidSystem.hpp>
#include <sstream>
#include <string>
Include dependency graph for blackoilmodel.hh:

Go to the source code of this file.

Classes

class  Ewoms::BlackOilModel< TypeTag >
 A fully-implicit black-oil flow model. More...
 
singleton  Ewoms::EclGridManager< TypeTag >
 
class  Ewoms::BlackOilModel< TypeTag >
 A fully-implicit black-oil flow model. More...
 

Namespaces

 Ewoms
 
 Ewoms::Properties
 

Functions

 Ewoms::Properties::NEW_TYPE_TAG (BlackOilModel, INHERITS_FROM(MultiPhaseBaseModel, VtkBlackOil, VtkComposition))
 The type tag for the black-oil problems. More...
 
 Ewoms::Properties::SET_TYPE_PROP (BlackOilModel, LocalResidual, Ewoms::BlackOilLocalResidual< TypeTag >)
 Set the local residual function. More...
 
 Ewoms::Properties::SET_TYPE_PROP (BlackOilModel, NewtonMethod, Ewoms::BlackOilNewtonMethod< TypeTag >)
 Use the black-oil specific newton method. More...
 
 Ewoms::Properties::SET_TYPE_PROP (BlackOilModel, Model, Ewoms::BlackOilModel< TypeTag >)
 The Model property. More...
 
 Ewoms::Properties::SET_TYPE_PROP (BlackOilModel, BaseProblem, Ewoms::BlackOilProblem< TypeTag >)
 The Problem property. More...
 
 Ewoms::Properties::SET_PROP (BlackOilModel, BlackOilFluidState)
 The BlackOilFluidState property. More...
 
 Ewoms::Properties::SET_TYPE_PROP (BlackOilModel, RateVector, Ewoms::BlackOilRateVector< TypeTag >)
 the RateVector property More...
 
 Ewoms::Properties::SET_TYPE_PROP (BlackOilModel, BoundaryRateVector, Ewoms::BlackOilBoundaryRateVector< TypeTag >)
 the BoundaryRateVector property More...
 
 Ewoms::Properties::SET_TYPE_PROP (BlackOilModel, PrimaryVariables, Ewoms::BlackOilPrimaryVariables< TypeTag >)
 the PrimaryVariables property More...
 
 Ewoms::Properties::SET_TYPE_PROP (BlackOilModel, IntensiveQuantities, Ewoms::BlackOilIntensiveQuantities< TypeTag >)
 the IntensiveQuantities property More...
 
 Ewoms::Properties::SET_TYPE_PROP (BlackOilModel, ExtensiveQuantities, Ewoms::BlackOilExtensiveQuantities< TypeTag >)
 the ExtensiveQuantities property More...
 
 Ewoms::Properties::SET_TYPE_PROP (BlackOilModel, Indices, Ewoms::BlackOilIndices< 0 >)
 The indices required by the model. More...
 
 Ewoms::Properties::SET_PROP (BlackOilModel, FluidSystem)
 Set the fluid system to the black-oil fluid system by default. More...
 
 Ewoms::Properties::SET_INT_PROP (BlackOilModel, BlackoilNumChoppedIterations, 4)
 

Detailed Description

A fully-implicit black-oil flow model.

The black-oil model is a three-phase, three-component model widely used for oil reservoir simulation. The phases are denoted by lower index $\alpha \in \{ w, g, o \}$ ("water", "gas" and "oil") and the components by upper index $\kappa \in \{ W, G, O \}$ ("Water", "Gas" and "Oil"). The model assumes partial miscibility:

  • Water and the gas phases are immisicible and are assumed to be only composed of the water and gas components respectively-
  • The oil phase is assumed to be a mixture of the gas and the oil components.

The densities of the phases are determined by so-called formation volume factors:

\[ B_\alpha := \frac{\varrho_\alpha(1\,\text{bar})}{\varrho_\alpha(p_\alpha)} \]

Since the gas and water phases are assumed to be immiscible, this is sufficint to calculate their density. For the formation volume factor of the the oil phase $B_o$ determines the density of saturated oil, i.e. the density of the oil phase if some gas phase is present.

The composition of the oil phase is given by the gas dissolution factor $R_s$, which defined as the volume of gas at atmospheric pressure that is dissolved in a given amount of oil at reservoir pressure:

\[ R_s := \frac{\varrho_{o}^G}{\varrho_o^O}\;. \]

This allows to calculate all quantities required for the mass-conservation equations for each component, i.e.

\[ \sum_\alpha \frac{\partial\;\phi c_\alpha^\kappa S_\alpha }{\partial t} - \sum_\alpha \mathrm{div} \left\{ c_\alpha^\kappa \mathbf{v}_\alpha \right\} - q^\kappa = 0 \;, \]

where $\mathrm{v}_\alpha$ is the filter velocity of the phase $\alpha$.

By default $\mathrm{v}_\alpha$ is determined by using the standard multi-phase Darcy approach, i.e.

\[ \mathbf{v}_\alpha = - \frac{k_{r\alpha}}{\mu_\alpha} \mathbf{K} *\left(\mathbf{grad}\, p_\alpha - \varrho_{\alpha} \mathbf{g} \right) \;, \]

although the actual approach which is used can be specified via the FluxModule property. For example, the velocity model can by changed to the Forchheimer approach by

SET_TYPE_PROP(MyProblemTypeTag, FluxModule, Ewoms::ForchheimerFluxModule<TypeTag>);

The primary variables used by this model are:

  • The pressure of the phase with the lowest index
  • The two saturations of the phases with the lowest indices