blackoilmodel.hh File Reference
#include <opm/material/densead/Math.hpp>
#include "blackoilproblem.hh"
#include "blackoilindices.hh"
#include "blackoiltwophaseindices.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 "blackoilsolventmodules.hh"
#include "blackoilpolymermodules.hh"
#include "blackoilfoammodules.hh"
#include "blackoilbrinemodules.hh"
#include "blackoilextbomodules.hh"
#include "blackoildarcyfluxmodule.hh"
#include "blackoilmicpmodules.hh"
#include <opm/models/common/multiphasebasemodel.hh>
#include <opm/models/io/vtkcompositionmodule.hh>
#include <opm/models/io/vtkblackoilmodule.hh>
#include "blackoildiffusionmodule.hh"
#include "blackoildispersionmodule.hh"
#include <opm/models/io/vtkdiffusionmodule.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

struct  Opm::Properties::TTag::BlackOilModel
 The type tag for the black-oil problems. More...
 
struct  Opm::Properties::LocalResidual< TypeTag, TTag::BlackOilModel >
 Set the local residual function. More...
 
struct  Opm::Properties::NewtonMethod< TypeTag, TTag::BlackOilModel >
 Use the black-oil specific newton method. More...
 
struct  Opm::Properties::Model< TypeTag, TTag::BlackOilModel >
 The Model property. More...
 
struct  Opm::Properties::BaseProblem< TypeTag, TTag::BlackOilModel >
 The Problem property. More...
 
struct  Opm::Properties::RateVector< TypeTag, TTag::BlackOilModel >
 the RateVector property More...
 
struct  Opm::Properties::BoundaryRateVector< TypeTag, TTag::BlackOilModel >
 the BoundaryRateVector property More...
 
struct  Opm::Properties::PrimaryVariables< TypeTag, TTag::BlackOilModel >
 the PrimaryVariables property More...
 
struct  Opm::Properties::IntensiveQuantities< TypeTag, TTag::BlackOilModel >
 the IntensiveQuantities property More...
 
struct  Opm::Properties::ExtensiveQuantities< TypeTag, TTag::BlackOilModel >
 the ExtensiveQuantities property More...
 
struct  Opm::Properties::FluxModule< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::Indices< TypeTag, TTag::BlackOilModel >
 The indices required by the model. More...
 
struct  Opm::Properties::FluidSystem< TypeTag, TTag::BlackOilModel >
 Set the fluid system to the black-oil fluid system by default. More...
 
struct  Opm::Properties::EnableSolvent< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnableExtbo< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnablePolymer< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnablePolymerMW< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnableFoam< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnableBrine< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnableVapwat< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnableDisgasInWater< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnableSaltPrecipitation< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnableMICP< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnableTemperature< TypeTag, TTag::BlackOilModel >
 By default, the blackoil model is isothermal and does not conserve energy. More...
 
struct  Opm::Properties::EnableEnergy< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::EnableDiffusion< TypeTag, TTag::BlackOilModel >
 disable diffusion by default More...
 
struct  Opm::Properties::EnableDispersion< TypeTag, TTag::BlackOilModel >
 disable disperison by default More...
 
struct  Opm::Properties::BlackOilEnergyScalingFactor< TypeTag, TTag::BlackOilModel >
 
struct  Opm::Properties::BlackoilConserveSurfaceVolume< TypeTag, TTag::BlackOilModel >
 
class  Opm::BlackOilModel< TypeTag >
 A fully-implicit black-oil flow model. More...
 

Namespaces

namespace  Opm
 
namespace  Opm::Properties
 
namespace  Opm::Properties::TTag
 The generic type tag for problems using the immiscible multi-phase model.
 

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 sufficient 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

template<class TypeTag>
struct FluxModule<TypeTag, TTag::MyProblemTypeTag> { using type = Opm::ForchheimerFluxModule<TypeTag>; };
Specifies a flux module which uses the Forchheimer relation.
Definition: forchheimerfluxmodule.hh:62

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