PNG 301
Introduction to Petroleum and Natural Gas Engineering

(link is external) (link is external)

5.4.2.1: Derivation of the Diffusivity Equation in Radial-Cylindrical Coordinates for Compressible Gas Flow

PrintPrint

As with the flow of oil, we begin the derivation of diffusivity equation for compressible gas flow with a mass balance on a thin ring or Representative Elemental Volume, REV, in the reservoir as shown in Figure 5.05. This mass balance results in:

Mass InMass Out=Mass Accumulation
Equation 5.26

We can elaborate on the definitions of terms in Equation 5.26 as:

Mass In=ρq]r+dr dt=(ρgqg]r+(ρgqg)rdr)dt
Equation 5.27
Mass Out=(ρgqg]r)dt
Equation 5.28

and,

Mass Accumulation=ϕVb REVρg=2πr dr ϕh ρg
Equation 5.29

Where:

  • ρg is the gas density, lb/ft3
  • qg is the gas rate, ft3/day
  • r is the radial coordinate in a radial-cylindrical coordinate system, ft
  • r is the radius of the representative elemental volume, REV, ft
  • t is time, days
  • ϕ is the porosity of the reservoir, fraction
  • Vb REV is the bulk volume of the representative elemental volume, REV, ft3
  • h is the reservoir thickness, ft
Diagram showing REV in a radial-cylindrical coordinates with arrows showing mass in and mass out pointing toward the center of the circle
Figure 5.05: Representative Elemental Volume, REV, in Radial-Cylindrical Coordinates
Source: Greg King

Substituting Equation 5.27 through Equation 5.29 into Equation 5.26 results in:

[(ρgqg]r+(ρgqg)rdr)ρgqg]r]dt=2πrdr ϕρg
Equation 5.30a

or,

(ρgqg)rdr dt=2πrdr ϕρg
Equation 5.30b

Dividing by the term r dr dt results in:

1r(ρgqg)r=2π ϕhρgdt
Equation 5.31

Now, substituting Darcy’s Law, Equation 5.06 with l=r and without the formation volume factor, Bg, (we want the flow rate in reservoir ft3/day not SCF/day). The unit conversion factor of 5.615 ft3/bbl converts Darcy’s Law from bbl/day to ft3/day:

1rr[(5.615)(0.001127)ρgkg(2πrh)μgpr]=2π ϕhρgdt
Equation 5.32

If we assume that the permeability, kg, and the thickness, h, are uniform, then we have:

1rr(rρgμgpr)=ϕ(5.615)(0.001127)kgρgdt
Equation 5.33a

or,

1rr(rρgμgpr)=ϕ0.006328 kgρgdt
Equation 5.33b

To this point, the derivation for the compressible gas equation is identical to the derivation for a slightly compressible liquid equation. This derivation will begin to deviate now. Applying the definition of density for a real gas (Equation 3.71):

ρg=p MWgZRTr
Equation 5.34

Equation 5.34 is a direct result of the Real Gas Law. Substituting into Equation 5.33b results in:

1rr(rpMWgμgZRTrpr)=ϕ0.006328 kgdt(pMWgZRTr)
Equation 5.35a

or,

1rr(rpμgZpr)=ϕ0.006328 kgdt(pZ)
Equation 5.35b

or,

1rr(rpμgZpr)=ϕ0.006328 kgdp(pZ)pdt
Equation 5.35c

Now, we saw in Lesson 3 that for a real gas, the definition of compressibility is (Equation 3.70b):

cg=pZddp(Zp)]T=Tr=pZddp[(pZ)1]]T=Tr=1p1ZdZdp
Equation 5.36a

or,

ddp(pZ)]T=Tr=cgpZ
Equation 5.36b

Substituting into Equation 5.35c results in:

1rr(rpμgZpr)=ϕcg0.006328 kgpZpdt
Equation 5.37

Equation 5.37 is the nonlinear diffusivity equation describing the transient behavior of compressible (real) gases. We say that it is Nonlinear because of the functions of pressure appearing in the equation. In this nonlinear form, we cannot solve the equation analytically (exactly). In order to obtain analytical solutions to this equation, we must first Linearize it. We do this in the same manner as we linearized the stabilized flow equations: by use of the Pressure, the Pressure-Squared, and the Pseudo-Pressure Formulations.