PNG 301
Introduction to Petroleum and Natural Gas Engineering

5.4.1.3: Stabilized Flow of Gas to a Vertical Production Well in Terms of the Real Gas Pseudo-Pressure, m(p)

PrintPrint

In the two previous developments of the stabilized inflow performance relationship for gas (in terms of pressure and the pressure-squared), we made simplifying assumptions on the pressure dependency of the group, p μ g  Z . I stated that within certain pressure ranges these assumptions were valid. However, there may be cases where these assumptions are not valid over the entire pressure range of interest. In these cases, we need to develop a general equation that is valid for all pressures. We can achieve this with the Real Gas Pseudo-Pressure Formulation. We define the real gas pseudo-pressure, m( p ) , as:

m( p )=2 14.7 p p μ g Z dp
Equation 5.19

The terminology, m( p ) , implies that the real gas pseudo-pressure is a function of pressure. The units of m( p ) are ps i 2 / cp . We can see from Equation 5.19 that:

m( p e )m( p wf )=2( 14.7 pe p μ g   Z dp 14.7 p wf p μ g   Z dp )=2 p wf p e p μ g   Z dp
Equation 5.20a

or,

p wf p e p μ g   Z dp = [ m( p e )m( p wf ) ] 2
Equation 5.20b

The relationships expressed in Equation 5.20 are illustrated in Figure 5.01. In this figure, p μ g  Z , is plotted on the y-axis and p is plotted on the x-axis. The blue area under the curve represents m( p wf ) 2 = 14.7 p wf p μ g  Z dp while the total area (blue plus green areas) under the curve represents m( p e ) 2 = 14.7 p e p μ g  Z dp . These are two of the integrals in Equation 5.20. The green area under the curve then represents [ m( p e )m( p wf ) ] 2 = p wf p e p μ g  Z dp .

Graph showing the pseudo pressure concept described above.
Figure 5.01: Illustration of the Pseudo-Pressure, m( p ) , Concept
Source: Greg King

Substituting Equation 5.20b into Equation 5.10 results in:

r w r e 1 r dr = 0.001127( 2π )( 5.615 ) T SC   k g   h 2   q g   T r   p SC [ m( p e )m( p wf ) ]
Equation 5.21

Performing the integration on the left-hand side of Equation 5.21 results in:

[ log e ( r ) ] r w r e = 0.001127( 2π )( 5.615 ) T SC   k g   h 2   q g   T r   p SC [ m( p e )m( p wf ) ]
Equation 5.22

Rearranging Equation 5.22 results in:

q g = 0.001127( 2π )( 5.615 ) T SC   k g   h 2   T r   p SC [ m( p e )m( p wf ) ] log e ( r e r w )
Equation 5.23a

or,

q g = 0.01988   T SC   k g   h T r   p SC [ m( p e )m( p wf ) ] log e ( r e r w )
Equation 5.23b

or, after substituting the normal U.S. definitions of p SC and T SC :

q g = 0.70325   k g   h T r [ m( p e )m( p wf ) ] log e ( r e r w )
Equation 5.23c

The only thing remaining for the pseudo-pressure formulation is the construction of the m( p ) function. To do this, we perform numerical integration on the p μ g  Z function. There are many ways to perform numerical integration; however, one of the most common methods is to approximate the area under the curve with a series of rectangles and to sum the areas of these rectangles. This is illustrated in Figure 5.02.

graph showing the numerical integration of the p/(ugZ) functiion described above.
Figure 5.02: Numerical Integration of the p μ g  Z Function
Source: Greg King

Using this approach, we can approximate the integral (area under the curve) with:

m(p)=2 14.7 p p μ g   Z dp2 n=2 N [ ( p μ g   Z ) n1 + ( p μ g   Z ) n ] 2 ( p n p n1 )
Equation 5.24

The Real Gas Pseudo-Pressure Formulation is the most rigorous and general formulation for the inflow performance relationship for gas and is valid over all pressure ranges. Even in pressure ranges where another formulation may be valid, the real gas pseudo-pressure formulation is more rigorous.

As I just mentioned, there are pressure ranges where the Pressure Formulation and the Pressure-Squared Formulation are valid. If we were to plot the p μ g  Z group (on the y-axis) versus p on the x-axis, then we would develop the plot shown in Figure 5.03. This plot illustrates that the p μ g  Z group is approximately constant for p > 3,000 psi and that the Pressure Formulation is valid when all pressures in the pressure range of interest are greater than 3,000 psi.

graph showing a plot of p/(ugz) where all lines rise until 3000 psi and then they all remain approximately constant above 3000 psi.
Figure 5.03: Plot of p μ g  Z and the Pressure Range Where it is Approximately Constant
Source: Greg King

We can also plot μ g  Z (on the y-axis) versus p on the x-axis to develop the plot shown in Figure 5.04. This plot illustrates that the μ g  Z product is approximately constant for p < 2,000 and that the Pressure-Squared Formulation is valid when all pressures are less than 2,000 psi.

graph showing a plot of ugZ highlighting where the approximately constant pressure range below 2000 psi. described above.
Figure 5.04: Plot of μ g  Z and the Pressure Range Where it is Approximately Constant
Source: Greg King