.. _formulation-of-black-oil: The Mathematical Formulation of the Black Oil Model =================================================== Introduction ---------------------------------- The black oil model is a generalization of the simple two-phase, two-component :ref:TOIL_IMS model. The black oil model is a three phase, three component model such that gas can dissolve in the reservoir oil phase. This implies that, in the case of the reservoir oil phase, more than one phase can contribute to the flow of a component, which extends the structure of the underlying equations, introduced for the formulation of :ref:thermal oil: .. math:: \frac{\partial}{\partial t} \phi \left(s_\alpha\eta_\alpha\right) + \nabla\cdot\left(\mathbf{q_\alpha}\eta_\alpha\right) = Q_i :label: molar_balance_bo .. math:: \frac{\partial}{\partial t} \left[ \phi \sum_{\alpha}{s_\alpha\eta_\alpha U_\alpha} + (1-\phi)\rho_r c_r T \right] + \nabla \cdot \sum_{\alpha}{\left[ q_\alpha\eta_\alpha H_\alpha + \kappa \nabla T \right]} = Q_e :label: eng_balance_bo .. math:: q_\alpha = \frac{K k_\alpha}{\mu_\alpha} \nabla(P_\alpha - \rho_\alpha g z) :label: darcy_bo The model assumes full instant equilibrium between the components in a cell, so that the oil-gas system can exist in two states possible states: - Saturated – there is free gas in the cell, and the oil is saturated, so that its bubble point (the pressure at which gas will come out of solution) is equal to the oil phase pressure. - Undersaturated – there is no free gas in the cell, and the bubble point pressure is less than the oil phase pressure. In each case, the system of oil, water and gas may be parameterized as a function of three primary variables, these being (:math:P_{o}, :math:S_{o}, :math:S_{g}) for the saturated case and (:math:P_{o}, :math:S_{o}, :math:P_{b}) for the undersaturated case. In the main solution arrays, one or other of these two possibilities is held for each cell, and a state variable is maintained to distinguish between the two cases. The amount of dissolved gas in the oil is commonly parameterized by the solution gas-oil ratio (normally denoted :math:R_{s}) which is the volume of dissolved gas (measured at surface conditions) divided by the volume of oil (also measured at surface conditions) in the reservoir oil phase. Oil properties such as formation volume factor (volume of oil at reservoir conditions/volume of oil at surface conditions), enthalpy and viscosity will generally be functions of pressure and bubble point. The variable substitution method was selected in preference to the alternative of using mass or molar density variables, which avoid the need to switch variables but introduce a flash calculation. This was driven by the wish to make best use of Pflotran’s saturation variable formulation, the reduced equation set in the absence of a extra volume balance condition, and greater ease of solution damping with respect to changes in the relative permeability functions. Phases and Components ---------------------------------- When oil and gas are produced from a reservoir, the process can be pictured as: .. image:: ../manual/images/oil_gas.png Hydrocarbon from the reservoir is produced through the wells, and is then run through separators, to yield stock tank oil (STO) and surface gas (SG). The black oil model assumes that the nature of the STO and SG remain fixed, and that, as the production always splits into these two fluids, the entire reservoir may be regarded as a combination of these two components. In the reservoir, of course, hydrocarbon can exist in liquid (reservoir oil, RO) and vapor (reservoir gas, RG) phases. In the basic dissolved gas picture, we further assume that reservoir gas contains only surface gas, whilst reservoir oil contains both stock tank oil and surface gas. As well as the extra bubble point pressure dependence of the oil properties, the black oil model requires, as data, the saturated solution gas-oil ratio as a function of the saturation pressure, :math:Rssat(P_b). Variable Switching ---------------------------------- In PFLOTRAN the state of a cell is defined by a number of primary variables, from which other secondary variables are derived. In the case of a saturated cell, these can be (:math:P_{o}, :math:S_{o}, :math:S_{g}) and for an undersaturated cell we can use (:math:P_{o}, :math:S_{o}, :math:P_{b}). Other options are, of course, possible, for example (:math:P_{o}, :math:S_{o}, :math:R_{s}). However, using the bubble point has a number of advantages: | It saves Rssat lookups | :math:P_{b} is a useful output variable | :math:P_{b} appears in the data (e.g. the properties of saturated oil). With :math:P_{b}, the switching rule is simple: | if(:math:`S_o>\varepsilon \; \& \; P_b