Three-phase relative permeability model¶

Overview¶

For the simulation of three-phase flow in porous media, it is common to use a specific treatment (i.e., different from the typical two-phase procedure) to evaluate the oil relative permeability. Specifically, the three-phase oil relative permeability is obtained by interpolation of oil-water and oil-gas experimental data measured independently in two-phase displacements.

Let $k_{rw,wo}$ and $k_{ro,wo}$ be the water-oil two-phase relative permeabilities for the water phase and the oil phase, respectively. Let $k_{rg,go}$ and $k_{ro,go}$ be the oil-gas two-phase relative permeabilities for the gas phase and the oil phase, respectively. In the current implementation, the two-phase relative permeability data is computed analytically using the Brooks-Corey relative permeability model.

The water and gas three-phase relative permeabilities are simply given by two-phase data and only depend on $S_w$ and $S_g$ , respectively. That is,

$k_{rw,wog}(S_w) = k_{rw,wo}(S_w),$

$k_{rg,wog}(S_g) = k_{rg,go}(S_g).$

The oil three-phase relative permeability is obtained using a variant of the saturation-weighted interpolation procedure initially proposed by Baker. Specifically, we compute:

$k_{ro,wog}(S_w,S_g) = \frac{ (S_w - S_{w,\textit{min}}) k_{ro,wo}(S_w) + S_g k_{rg,go}(S_g) }{ (S_w - S_{w,\textit{min}}) + S_g }.$

This procedure provides a simple but effective formula avoiding the problems associated with the other interpolation methods (negative values).

Parameters¶

The relative permeability constitutive model is listed in <Constitutive> block of the input XML file. The relative permeability model must be assigned a unique name via name attribute. This name is used to assign the model to regions of the physical domain via a materialList attribute of the <ElementRegion> node.

The following attributes are supported:

Name	Type	Default	Description
gasOilRelPermExponent	real64_array	{1}	Rel perm power law exponent for the pair (gas phase, oil phase) at residual water saturation
gasOilRelPermMaxValue	real64_array	{0}	Maximum rel perm value for the pair (gas phase, oil phase) at residual water saturation
name	string	required	A name is required for any non-unique nodes
phaseMinVolumeFraction	real64_array	{0}	Minimum volume fraction value for each phase
phaseNames	string_array	required	List of fluid phases
waterOilRelPermExponent	real64_array	{1}	Rel perm power law exponent for the pair (water phase, oil phase) at residual gas saturation
waterOilRelPermMaxValue	real64_array	{0}	Maximum rel perm value for the pair (water phase, oil phase) at residual gas saturation

Below are some comments on the model parameters.

phaseNames - The number of phases should be 3. Supported phase names are:

Value	Phase
oil	Oil phase
gas	Gas phase
water	Water phase

phaseMinVolFraction - The list of minimum volume fractions $S_{\ell,min}$ for each phase is specified in the same order as in phaseNames. Below this volume fraction, the phase is assumed to be immobile.
waterOilRelPermExponent - The list of exponents $\lambda_{\ell,wo}$ for the two-phase water-oil relative permeability data, with the water exponent first and the oil exponent next. These exponents are then used to compute $k_{r \ell,wo}$ in the Brooks-Corey relative permeability model.
waterOilRelPermMaxValue - The list of maximum values $k_{\textit{r} \ell,wo,\textit{max}}$ for the two-phase water-oil relative permeability data, with the water max value first and the oil max value next. These exponents are then used to compute $k_{r \ell,wo}$ in the Brooks-Corey relative permeability model.
gasOilRelPermExponent - The list of exponents $\lambda_{\ell,go}$ for the two-phase gas-oil relative permeability data, with the gas exponent first and the oil exponent next. These exponents are then used to compute $k_{r \ell,go}$ in the Brooks-Corey relative permeability model.
gasOilRelPermMaxValue - The list of maximum values $k_{\textit{r} \ell,go,\textit{max}}$ for the two-phase gas-oil relative permeability data, with the gas max value first and the oil max value next. These exponents are then used to compute $k_{r \ell,go}$ in the Brooks-Corey relative permeability model.

Example¶

<Constitutive>
   ...
 <BrooksCoreyBakerRelativePermeability name="relperm"
                                       phaseNames="oil gas water"
                                       phaseMinVolumeFraction="0.05 0.05 0.05"
                                       waterOilRelPermExponent="2.5 1.5"
                                       waterOilRelPermMaxValue="0.8 0.9"
                                       gasOilRelPermExponent="3 3"
                                       gasOilRelPermMaxValue="0.4 0.9"/>
   ...
</Constitutive>