Constitutive Relations

Capillary Pressure - Saturation Functions

van Genuchten Saturation Function

Capillary pressure is related to saturation by various phenomenological relations, one of which is the van Genuchten (1980) relation

(1)se=[1+(pcpc0)n]m,

where pc represents the capillary pressure [Pa], and the effective saturation se is defined by

(2)se=ssrs0sr,

where sr denotes the residual saturation, and s0 denotes the maximum saturation. The inverse relation is given by

(3)pc=pc0(se1/m1)1/n.

The quantities m, n and pc0 are impirical constants determined by fitting to experimental data.

Brooks-Corey Saturation Function

The Brooks-Corey saturation function is a limiting form of the van Genuchten relation for pc/pc01, with the form

(4)se=(pcpc0)λ,

with λ=mn and inverse relation

(5)pc=pc0se1/λ.

Relative Permeability Functions

Two forms of the relative permeability function are implemented based on the Mualem and Burdine formulations. The quantity n is related to m by the expression

(6)m=11n,     n=11m,

for the Mualem formulation and by

(7)m=12n,     n=21m,

for the Burdine formulation.

Mualem Relative Permeability

For the Mualem relative permeability function based on the van Genuchten saturation function is given by the expression

(8)kr=se{1[1(se)1/m]m}2.

The Mualem relative permeability function based on the Brooks-Corey saturation function is defined by

(9)kr=(se)5/2+2/λ=(pc/pc0)(5λ/2+2).

Burdine Relative Permeability

For the Burdine relative permeability function based on the van Genuchten saturation function is given by the expression

(10)kr=se2{1[1(se)1/m]m}.

The Burdine relative permeability function based on the Brooks-Corey saturation function has the form

(11)kr=(se)3+2/λ=(pcpc0)(3+2λ).

Modified Brooks Corey Relative Permeability

The modified Brooks Corey relative permeability function can be associated with any saturation function. The liquid relative permeability is defined as

(12)sel=(slsrl1srlsrg)krl=krl,maxselnl

The gas phase relative permeability is defined as

(13)seg=(1slsrg1srlsrg)krg=krg,maxsegng

where krα,max, nα and seα are the maximum relative permeability, modified Brooks Corey exponent and effective saturation for phase α.

Smoothing

At the end points of the saturation and relative permeability functions it is sometimes necessary to smooth the functions in order for the Newton-Raphson equations to converge. This is accomplished using a third order polynomial interpolation by matching the values of the function to be fit (capillary pressure or relative permeability), and imposing zero slope at the fully saturated end point and matching the derivative at a chosen variably saturated point that is close to fully saturated. The resulting equations for coefficients ai, i=03, are given by

(14)a0+a1x1+a2x12+a3x13=f1,a0+a1x2+a2x22+a3x23=f2,a1x1+2a2x1+3a3x12=f1,a1x2+2a2x2+3a3x22=f2,

for chosen points x1 and x2. In matrix form these equations become

(15)[1x1x12x131x2x22x23012x13x12012x23x22][a0a1a2a3]=[f1f2f1f2].

The conditions imposed on the smoothing equations for capillary pressure f=se(pc) are x1=2pc0, x2=pc0/2, f1=(se)1, f2=1, f1=(se)1, f2=0. For relative permeability f=kr(se), x1=1, x2=0.99, f1=1, f2=(kr)2, f1=0, f2=(kr)2.