1.

1 Primary Reservoir Characteristics general, reservoir fluids are classified into

three groups:

Flow in porous media is a very complex (1) Incompressible fluids;

phenomen and cannot be described as (2) Slightly compressible fluids;
explicitly as flow through pipes or conduits. It (3) Compressible fluids.
is rather easy to measure the length and
The isothermal compressibility coefficient is
diameter of a pipe and compute its flow
described mathematically by the following
capacity as a function of pressure; however,
two equivalent expressions:
in porous media flow is different in that
there are no clear-cut flow paths which lend In terms of fluid volume:
themselves to measurement.
−1 𝜕𝑉
𝑐= 𝑉 𝜕𝑝
[1.1.1]
The analysis of fluid flow in porous media has
evolved throughout the years along two
fronts: the experimental and the analytical.
Physicists, engineers, hydrologists, and the In terms of fluid density:
like have examined experimentally the −1 𝜕𝜌
𝑐= [1.1.2]
behavior of various fluids as they flow 𝜌 𝜕𝑝

through porous media ranging from sand

packs to fused Pyrex glass. On the basis of
their analyses, they have attempted to Where
formulate laws and correlations that can be 𝑉 = 𝑓𝑙𝑢𝑖𝑑 𝑣𝑜𝑙𝑢𝑚𝑒
utilized to make analytical predictions for
𝜌 = 𝑓𝑙𝑢𝑖𝑑 𝑑𝑒𝑛𝑠𝑖𝑡𝑦
similar systems.
𝑝 = 𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒, 𝑝𝑠𝑖 −1
The main objective of this chapter is to
𝑐 = 𝑖𝑠𝑜𝑡ℎ𝑒𝑟𝑚𝑎𝑙 𝑐𝑜𝑚𝑝𝑟𝑒𝑠𝑠𝑖𝑏𝑖𝑙𝑖𝑡𝑦 𝑐𝑜𝑒𝑓𝑓𝑖𝑐𝑖𝑒𝑛𝑡, 𝜓 −1
present the mathematical relationships that
are designed to describe the flow behavior of In compressible fluids
the reservoir fluids. The mathematical forms
An incompressible fluid is defined as the fluid
of these relationships will vary depending
whose volume or density does not change
upon the characteristics of the reservoir.
with pressure. That is
These primary reservoir characteristics that
must be considered include: 𝜕𝑉
= 0 and
𝜕𝜌
=0
𝜕𝑝 𝜕𝑝
 Types of fluids in the resevoir;
Incompressible fluids do not exist; however,
 Flow regimes;
this behavior may be asumed in some cases
 Reservoir geometry;
to simplify the derivation and the final form
 Number of flowing fluids in the
of many flow equations.
reservoir.
Slightly compressible fluids
1.1.1 Types of fluids
These ´´ Slightly`` compressible fluids exhibit
The isothermal compressibility coefficient is small changes in volume, or density, with
essentially the controlling factor in changes in pressure. Knowing the volume of
identifying the type of the reservoir fluid. In a slightly compressible liquid at a reference
(initial) pressure, the changes in the
volumetric behavior of this fluid as a function 𝑝𝑟𝑒𝑓 = 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 𝑖𝑛𝑖𝑡𝑖𝑎𝑙 ( 𝑟𝑒𝑓𝑒𝑟𝑒𝑛𝑐𝑒)𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒 𝑝𝑟𝑒𝑓
of pressure can be mathematically describe It should be pointed out that crude oil and water
by integrating Equation 1.1.1, to give: systems fi into this category.
𝑝 𝑉
𝑑𝑉 Compressible fluids
−𝑐 ∫ 𝑑𝑝 = ∫
𝑝𝑟𝑒𝑓 𝑉𝑟𝑒𝑓 𝑉
These are fluids that experience large changes in
𝑉 volume as a function of pressure. All gases are
𝑒𝑥𝑝[𝑐(𝑝𝑟𝑒𝑓 − 𝑝)] = considered compressible fluids. The truncation of
𝑉𝑟𝑒𝑓
the series expansion as given by Equation 1.1.5 is
𝑉 = 𝑉𝑟𝑒𝑓 𝑒𝑥𝑝[𝑐(𝑝𝑟𝑒𝑓 − 𝑝)] [1.1.3] not valid in this category and the complete
expansion as given by Equation 1.1.4 is used.

The isothermal compressibility of any

Where: compressible fluid is described by the following
expression:
𝑝 = 𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒, 𝑝𝑠𝑖𝑎
1 1 𝜕𝑍
𝑐𝑔 = − ( ) [1.1.8]
𝑉 = 𝑣𝑜𝑙𝑢𝑚𝑒 𝑎𝑡 𝑝𝑟𝑒𝑠𝑠𝑢𝑒 𝑝, 𝑓𝑡 3 𝑝 𝑍 𝜕𝑝 𝑇

𝑝𝑟𝑒𝑓 = 𝑖𝑛𝑖𝑡𝑖𝑎𝑙 ( 𝑟𝑒𝑓𝑒𝑟𝑒𝑛𝑐𝑒)𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒, 𝑝𝑠𝑖𝑎 Figures 1.1 and 1.2 show schematic illustrations
of the volume and density changes as a function
𝑉𝑟𝑒𝑓 of pressure for the three types of fluids.
= 𝑓𝑙𝑢𝑖𝑑 𝑣𝑜𝑙𝑢𝑚𝑒 𝑎𝑡 𝑖𝑛𝑖𝑡𝑖𝑎𝑙 ( 𝑟𝑒𝑓𝑒𝑟𝑒𝑛𝑐𝑒) 𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒, 𝑝𝑠𝑖𝑎
1.1.2 Flow regimes
The exponential ℯ 𝑥 may be represented by a
series expansion as: There are basically three types of flow regimes
that must be recognized in order to describe the
𝑥2 𝑥2 𝑥𝑛 fluid flow behavior and reservoir pressure
𝑒𝑥 = 1 + 𝑥 + + + ⋯+ [1.1.4]
2! 3! 𝑛!
distribution as a function of time. These three
Because the exponent 𝑥 (which represents flow regimes are:
the term 𝑐(𝑝𝑟𝑒𝑓 − 𝑝) is very small, the 𝑒 𝑥 (1) Steady-state flow;
term can be approximated by truncating (2) Unsteady-state flow;
Equation 1.1.4 to: (3) Pseudosteady-state flow.

𝑒𝑥 = 1 + 𝑥 [1.1.5]

Combining Equation 1.1.5 with 1.1.3 gives:

𝑉 = 𝑉𝑟𝑒𝑓 [1 + 𝑐(𝑝𝑟𝑒𝑓 − 𝑝)] [1.1.6]

A similar derivation is applied to Equation

1.1.2, to give:

𝑝 = 𝑝𝑟𝑒𝑓 [1 − 𝑐(𝑝𝑟𝑒𝑓 − 𝑝)] [1.1.7]

Where:

𝑉 = 𝑣𝑜𝑙𝑢𝑚𝑒 𝑎𝑡 𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒 𝑝
𝜌 = 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 𝑎𝑡 𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒 𝑝
𝑉𝑟𝑒𝑓
= 𝑣𝑜𝑙𝑢𝑚𝑒 𝑎𝑡 𝑖𝑛𝑖𝑡𝑖𝑎𝑙 ( 𝑟𝑒𝑓𝑒𝑟𝑒𝑛𝑐𝑒)𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒 𝑝𝑟𝑒𝑓
 with respect to time at every
position is constant, or:
Steady-state flow
The flow regime is identified as a 𝜕𝑝
steady-state flow if the pressure at ( 𝜕𝑡 ) = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡 [1.1.11]
𝑖
every location in the reservoir It should be pointed out that
remains constant, i.e., does not pseudosteady-state flow and
change with time. Mathematically, quasisteady-state flow.
this condition is expressed as: Figure 1.3 shows a schematic
comparison of the pressure declines
as a function of time of the three
𝜕𝑝 flow regimes.
( 𝜕𝑡 ) = 0 [1.1.9]
𝑖
This equation states that are rate of
change of pressure with respect to
time at any location is zero. In
reservoirs, the steady-state flow
condition can only occur when the
reservoir is completely recharged
and supported by strong aquifer or
pressure maintenance operations.

Unsteady-state flow

Unsteady-state flow (frequently

called transient flow) is defined as
the fluid flowing condition at which
the rate of change of pressure with
respect to time at any position in the
reservoir is not zero or constant. This
definition suggests that the pressure
derivative with respect to time is
essentially a function of both
position and time, thus:
𝜕𝑝
( 𝜕𝑡 ) = 𝑓(𝑖, 𝑡) [1.1.10]
1.1.3 Reservoir geometry
Pseudosteady-state flow
When the pressure at different The shape of a reservoir has a
locations in the reservoir is declining significant effect on its flow
linearly as a function of time, i.e., at behavior. Most reservoirs have
a constant declining rate, the flowing irregular boundaries and a rigorous
condition is characterized as mathematical description of their
pseudosteady-state flow. geometry is often possible only with
Mathematically, this definition states the use of numerical simulators.
that the rate of change of pressure However, for many engineering
 purposes, the actual flow geometry
may be represented by one of the
following flow geometrics:

 Radial flow;
 Linear flow;
 Spherical and hemispherical
flow.

Radial flow

In the absence of severe reservoir

heterogeneities, flow into or away from a
wellbore will follow radial flow lines a
substantial distance from the wellbore.
Because fluids move toward the well from all
directions and coverage at the wellbore, the
term radial flow is used to characterize the
flow of fluid into the wellbore. Figure 1.4
shows idealized flow lines and isopotential
lines for a radial flow system.

Linear flow

Linear flow occurs when flow paths are

parallel and the fluid flows in a single Area to flow must be constant. Figure 1.5
direction. In addition, the cross-sectional shows and idealize linear flow system. A
common application of linear flow equations
is the fluid flow into vertical hydraulic
fractures as illustrated in figure 1.6.

Spherical and hemispherical flow

Depending upon the type of wellbore

completion configuration, it is possible to
have spherical or hemispherical flow near
the wellbore. A well with a limited
perforated interval could result in spherical
flow in the vicinity of the perforations as
illustrated in figure1.7. A well which only
partially penetrates the pay zone, as shown
in figure 1.8, could result in hemispherical
flow. The condition could arise where coning
of bottom water is important.

1.1.4 Number of flowing fluids in the

reservoir
The mathematical expressions that are used 𝑞 𝑘 𝑑𝑝
𝑣= = −
to the volumetric performance and pressure 𝐴 𝜇 𝑑𝑥
behavior of a reservoir vary in form and
complexity depending upon the number of
mobile fluids in the reservoir. There are 𝑣 is the apparent velocity in centimeters per
𝑞
generally three cases of flowing system: second and is equal to 𝐴 , where 𝑞 is the
(1) Single-phase flow ( oil, water, or volumetric flow rate in cubic centimeters per
gas); second and 𝐴 is the total cross-sectional area
(2) Two-phase flow ( oil-water, oil-gas, of the rock in square centimeters. In other
or gas-water); words, 𝐴 includes the area of the rock
(3) Three-phase flow (oil, water; and material as well as the area of the pore
gas). channels. The fluid viscosity,𝜇, is expressed
in centipoise units, and the pressure
The description of fluid flow and subsequent 𝑑𝑝
gradient, 𝑑𝑥
, is in atmospheres per
analysis of pressure data becomes more
difficult as the number of mobile fluids centimeter, taken in the same direction as 𝑣
and 𝑞. The proportionality constant,𝑘, is the
increases.
permeability of the rock expressed in Darcy
1.2 Fluid Flow Equations units.
The fluid flow equations that are used to The negative sign in equation 1.2.1a is added
describe the flow behavior in a reservoir can 𝑑𝑝
because the pressure gradient 𝑑𝑥
is negative
take many forms depending upon the
in the direction of flow as showm in figure
combination of variables presented
1.9.
previously (i.i., types of flow, types of fluid,
etc.). by combining the conservation of mass
equation with the transport equation(
Darcy`s equation) and various equations of
state, the necessary flow equations can be
developed. Since all flow equations to be
considered depend on Darcy`s law, it is
important to consider this transport
relationship first.

1.2.1 Darcy`s law

The fundamental law of fluid motion in

porous media is Darcy`s law. The
mathematical expression developed by Darcy
in 1956 states that the velocity of a
homogeneous fluid in a porous medium is
For a horizontal- radial system, the pressure
proportional to the pressure gradient, and
gradient is positive ( see figure 1.10) and
inversely proportional to the fluid viscosity.
Darcy`s equation can be expressed in the
For a horizontal linear system, this
following generalized radial form:
relationship is:
 𝑞𝑟 𝑘 𝜕𝑝
𝑣= = ( )
𝐴𝑟 𝜇 𝜕𝑟 𝑟

Where:

𝑞𝑟 = 𝑣𝑜𝑙𝑢𝑚𝑒𝑡𝑟𝑖𝑐 𝑓𝑙𝑜𝑤 𝑟𝑎𝑡𝑒 𝑎𝑡 𝑟𝑎𝑑𝑖𝑢𝑠 𝑟

𝐴𝑟 = 𝑐𝑟𝑜𝑠𝑠 − 𝑠𝑒𝑐𝑡𝑖𝑜𝑛𝑎𝑙 𝑎𝑟𝑒𝑎 𝑡𝑜 𝑓𝑙𝑜𝑤 𝑎𝑡 𝑟𝑎𝑑𝑖𝑢𝑠 𝑟

𝜕𝑝
( ) = 𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒 𝑔𝑟𝑎𝑑𝑖𝑒𝑛𝑡 𝑎𝑡 𝑟𝑎𝑑𝑖𝑢𝑠 𝑟
𝜕𝑟 𝑟  Radial flow of compressible fluids;
 Multiphase flow.
𝑣 = 𝑎𝑝𝑝𝑎𝑟𝑒𝑛𝑡 𝑣𝑒𝑙𝑜𝑐𝑖𝑡𝑦 𝑎𝑡 𝑟𝑎𝑑𝑖𝑢𝑠 𝑟
Linear flow of incompressible fluids
The cross-sectional are at radius r is
essentially the surface area of a cylinder. For In a linear system, it is assumed that the flow occurs
through a constant cross-sectional area 𝐴, where both
a fully penetrated well with a net thickness
end are entirely open to flow. It is also assumed that
of ℎ, the cross-sectional area 𝐴, is given by: no flow crosses the sides, top, or bottom as shown in
figure 1.11. I f an incompressible fluid is flowing across
𝐴𝑟 = 2𝜋𝑟ℎ
the element dx, then the fluid velocity 𝑣 and the flow
Darcy`s law applies only when the following conditions rate 𝑞 are constants at all points. The flow behavior in
exist: this system can be expressed by the differential form
of Darcy`s equation, i.e., Equation 1.2.1a. Separating
 Laminar ( viscous) flow; the variables of Equation 1.12.1a and integrating over
 Steady-state flow; the length of the linear system:
 Incompressible fluids;
 Homogeneous formation. 𝑞 𝐿 𝑘 𝑝2
∫ 𝑑𝑥 = − ∫ 𝑑𝑝
𝐴 0 𝑢 𝑝1
For turbulent flow, which occurs at higher velocities,
the pressure gradient increases at a greater rate than Which results in:
does the flow rate and a special modification of Darcy`s
equation is needed. When turbulent flow exists, the
application of Darcy`s equation can result in serious
errors. Modifications for turbulent flow will be
discussed later in this chapter. 𝑘𝐴(𝑃1 − 𝑃2 )
𝑞=
𝜇𝐿
1.2.2 Steady-state flow
It is desirable to express the above relationship in
As defined previously, steady-state flow represents the customary field units, or:
condition that exists when the pressure throughout the
reservoir does not change with time. The applications 0.001127𝑘𝐴(𝑝1 − 𝑝2 )
𝑞=
of steady-state flow to describe the flow behavior of 𝜇𝐿
several types of fluid in different reservoir geometries
Where:
are presented below. These include:
𝑏𝑏𝑙
 Linear flow of incompressible fluids; 𝑞 = 𝑓𝑙𝑜𝑤 𝑟𝑎𝑡𝑒,
𝑑𝑎𝑦
 Linear flow of slightly compressible fluids;
𝑘 = 𝑎𝑏𝑠𝑜𝑙𝑢𝑡𝑒 𝑝𝑒𝑟𝑚𝑒𝑎𝑏𝑖𝑙𝑖𝑡𝑦, 𝑚𝑑
 Linear flow of compressible fluids;
 Radial flow of incompressible fluids; 𝑝 = 𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒, 𝑝𝑠𝑖𝑎
 Radial flow of slightly compressible fluids;
𝜇 = 𝑣𝑖𝑠𝑐𝑜𝑠𝑖𝑡𝑦, 𝑐𝑝

𝐿 = 𝑑𝑖𝑠𝑡𝑎𝑛𝑐𝑒, 𝑓𝑡
 𝐴 = 𝑐𝑟𝑜𝑠𝑠 − 𝑠𝑒𝑐𝑡𝑖𝑜𝑛 𝑎𝑟𝑒𝑎 𝑓𝑡 2
𝑞 (1.6905)(5.615) 𝑓𝑡
Example 1.1 An incompressible fluid flows in a linear 𝑣= = = 0.0105
𝐴 (0.15)(6000) 𝑑𝑎𝑦
porous media with the following properties:

𝐿 = 2000 𝑓𝑡 ℎ = 20 𝑓𝑡 Width = 300 ft

𝑘 = 100 𝑚𝑑 𝜙 = 15% 𝜇 = 2 𝑐𝑝 (c) Calculate the actual fluid velocity :

𝑝1 = 2000 𝑝𝑠𝑖 𝑝2 = 1990 𝑝𝑠𝑖 𝑞 (1.6905)(5.615)

𝑣= = = 0.0105 𝑓𝑡/𝑑𝑎𝑦
𝜙𝐴 (0.15)(6000)
Calculate:
The difference in the pressure (𝑝1 − 𝑝2 ) in
(a) Flow rate in bbl /day
Equation 1.2.2 is not only driving force in a tilted
(b) Apparent fluid velocity in ft / day;
(c) Actual fluid velocity in ft / day. reservoir. The gravitational force is the other
important driving force that must be accounted
Solution Calculate the cross-sectional area A: for to determine the direction and rate of flow.
𝐴 = (ℎ)(𝑤𝑖𝑑𝑡ℎ) = (20)(100) = 6000 𝑓𝑡 2
The fluid gradient force (gravitational force) is
always directed vertically downward while the
(a) Calculate the flow rate from Equation1.2.2: force that results from applied pressure drop may
be in any direction. The force causing flow would
0.001127𝑘𝐴(𝑝1 − 𝑝2 ) then be the vector sum of these two. In practice
𝑞= we obtain this result by introducing a new
𝜇𝐿
parameter, called ´´fluid potential`` which has the
(0.001127)(100)(6000)(2000 − 1990)
𝑞= same dimensions as pressure, e.g., psi. Its symbol
(2)(2000)
is
= 1.6905 𝑏𝑏𝑙/𝑑𝑎𝑦
Eq
(b) Calculate the apparent velocity:
uation 1.2.15 can be arranged to solve for the Solution
pressur p at any radius r, to give:

Step 1. Rearrange Equation 1.2.15arid solve

Example 1.5 An oil well in the Nameless Field forthe pressure , p at radius r:
is pro ducing at a stabilized rate of 600
STB/day at a stabilized bottom-hole flowing
pressure of 1800 psi. Analysis of the pressure Step 2. Calculate the pressure at the
buildup test data indicates that the pay zone designated radii
is characterized by a permeability of 120 md
and a uniform thickness of 25 ft. The well
drains an area of approximately 40 acres. The
following additional data is available:
Figure 1.14 shows the pressure profile as a
function of radius for the calculated data.
Calculate the pressure profile (distribution)
and list the pres- sure drop across 1 ft
intervals fromr to 1.25 ft, 4 to 5 ft, 19 to 20 ft, Results of the above example reveal that the
99 to 100 ft, and 744 to 745 ft. pressure drop just around the wellbore (i.e.,
142 psi) is 7.5 times greater than at the 4 to 5
interval, 36 times greater than at 19-20 ft,
and 142 times than that at the 99-100 ft fluids. If this equation is substituted into the
interval. The reason for this large pressure radial form of Darcy's law, the following is
drop around the wellbore is that the fluid obtained:
flows in from a large drainage area of 40
acres.
where gref is the flow rate at some reference
pressure pref. Separating the variables and
The external pressure pe used in Equation assuming a constant com pressibility over the
1.2.15cannot be measured readily, but pe entire pressure drop, and integrating over the
does not deviate length of the porous medium:

substantially from the initial reservoir

pressure if a strong and active aquifer is
where qref is the oil flow rate at a reference
present
pressure Pref Choosing the bottom-hole flow
pressure wt as the reference pressure and
expressing the flow rate in STB/day gives:
Several authors have suggested that the
average reservoir pressure pr, which often is
reported in well test results, uld be used in
Example 1.6 The following data is available on
performing material balance calcula- ons and
a well in the Red River Field
flow rate prediction. Craft arid Hawkins
(1959) showed that the average pressure is
located at about 61 % of the drainage radius
re for a steady-state flow condition. Assuming a slightly compressible fluid,
calculate the oil low rate. Compare the result
with that of an incompressible fluid.
Substituting 0.61re in Equation 1.2.17 gives:

Solution For a slightly compressible fluid, the

oil flow rate can be calculated by applying
Golan and Whitson (1986) suggested a
Equation 1.2.21:
method for approximating the drainage area
of wells producing from a common reservoir.
These authors assume that the volume
drained by a single well is proportional to its
rate of flow. Assuming constant reservoir Assuming an incompressible fluid, the flow
properties and a uniform thickness, the rate can be estimated by applying Darcy's
approximate drainage area of a single well A equation, i.e. Equation 1.2.15:
is:

Radial flow of compressible gases

Radial flow of slightly compressible fluids
The basic differentiak form of Darcy's law for
Terry and co-authors (1991) used Equation a horizontal laminar flow is valid for
1.2.6 to express the dependency of the flow describing the flow of both gas and liquid
rate on pressure for slightly compressible
systems. For a radial gas flow, Darcy's under the curve from p 0 to any pressure p
equation takes the form: represents the value of corresponding to p.
The following example will illustrate the
procedure.
can be expanded to give:

Example 1.7 The PVT data from a gas well in

Replacing the integral in Equation 1.2.24 with the Anaconda Gas Field is given below
the abone expanded form yields:

The well is producing at a stabilized bottom-

The integral l 2p/(Z)dp is called the "realgas hole flowing pressure of 3600 psi. The
pseudo- potential" or "real-gas wellbore radius is 0.3 ft. The following
pseudopressure" and it is usually represented additional data is available:
by m(p) or v. Thus:

Calculate the gas flow rate in Mscf/day.

Equation 1.2.27 can be written in terms of
the real-gas pseudopressure as:
Solution

Equation 1.2.28 indicates that a graph of v

vs. In(r/r.) yields a straight line with a slope of Step 1.
Q T/0.703kh and an intercept value of vw as
shown in Figure 1.15. The exact flow rate is
then given by Calculate the term 2p/ugZ for each pressure
as shown below

Because the gas flow rate is commonly

expressed in Mscf/day, Equation 1.2.30 can Step 2. Plot the term 2p/usZ versus pressure
be expressed as: as shown in Figure 1.16.

Equation 1.2.31 can be expressed in terms of Step 3. Calculate numerically the area under
the average reservoir pressure p, instead of the curve for each value of p. These areas
the initial reservoir pressure Pe as: correspond to the real gas pseudopressure at
each pressure. These v values are tabulated
below; notice that 2p/gZ vs p is also plotted in
To calculate the integral in. Equation 1.2.31, the figure.
the values of 2p/u Z are calculated for several
values of pressure p. Them 2p/p Z vs. p is
plotted on a Cartesian scale and the area Step 4. Calculate the flow rate by applying
under the curve is calculated either Equation 1.2.30
numerically or graph ically, where the area
In the approximation of the gas flow rate, the Solution
exact ga flow rate as expressed by the
Step 1. Calculate the arithmetic average
different forms of Darcy's las i.e., Equations
pressure:
1.2.25 through 1.2.32, can be approximated
moving the term 2/H Z outside the integral as
a constant should be pointed out that the
product of Zpg is considere constant only Step 2. Determine the gas viscosity and gas
under a pressure range of less than 2000 ps compressibility factor at 4020 psi:
Equation 1.2.31 can be rewritten as

Step 3. Apply Equation 1.2.33

where:

Step 4. Results show that the pressure-

The term (uZg is evaluated at an average squared method approximates the exact
pressure p that is defined by the following solution of 37 614 with an absolute error of
expression: 1.86 % . This error is due to the lim- ited
applicability of the pressure-squared method
to a pressure range of less than 2000 psi.
The above approximation method is called
the pressure squared method and is limited
to flow calculations when the reservoir Horizontal multiple-phase flow
pressure is less that 2000 psi. Other
approximation methods are discussed in
Chapter 2. When several fluid phases are flowing
simultaneously in a horizontal porous system,
the concept of the effective perme- ability of
Example 1.8 Using the data given in Example each phase and the associated physical
1.7, re- solve the gas flow rate by using the properties must be used in Darcy's equation.
pressure-squared method. Compare with the For a radial system, the generalized form of
exact method (i.e., real-gas pseudopressure Darcy's equation can be applied to each
solution). reservoir as follows