Fundamentals of Reservoir Engineering & Characterization 1 RESERVOIR ENGINEERING The purpose of Reservoir engineering is economic optimization of the development and production of hydrocarbon reservoirs. This requires most representative solutions to the following aspects: • Quantity of hydrocarbon in place • Recoverable hydrocarbons reserves • Rate of exploitation The determination of these three quantities is the crux of reservoir engineering. RESERVOIR A reservoir is a porous and permeable subsurface formation containing hydrocarbon accumulation. For a reservoir to be commercially exploitable, three basic requirements must be fulfilled: • Sufficient void space generally called porosity to store oil and gas. • Adequate connectivity, i.e. permeability to allow hydrocarbon fluids movement over large distances under pressure gradients. • Accumulation in a trap of impervious cap rock, which should prevent upward migration of the oil and gas. Accumulation of oil and gas in a reservoir (After “Reservoir and Production Fundamentals”, Schlumberger, 1982)
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Fundamentals of Reservoir Engineering & Characterization 1
RESERVOIR ENGINEERING
The purpose of Reservoir engineering is economic optimization of the development and
production of hydrocarbon reservoirs. This requires most representative solutions to the
following aspects:
• Quantity of hydrocarbon in place
• Recoverable hydrocarbons reserves
• Rate of exploitation
The determination of these three quantities is the crux of reservoir engineering.
RESERVOIR
A reservoir is a porous and permeable subsurface formation containing hydrocarbon
accumulation. For a reservoir to be commercially exploitable, three basic requirements
must be fulfilled:
• Sufficient void space generally called porosity to store oil and gas.
• Adequate connectivity, i.e. permeability to allow hydrocarbon fluids movement
over large distances under pressure gradients.
• Accumulation in a trap of impervious cap rock, which should prevent upward
migration of the oil and gas.
Accumulation of oil and gas in a reservoir
(After “Reservoir and Production Fundamentals”, Schlumberger, 1982)
Fundamentals of Reservoir Engineering & Characterization 2
RESERVOIR ROCKS
These rocks are generally sedimentary rocks. Sedimentary rocks are rocks made up of
sediments formed at the earth’s surface by debris or chemical precipitations.
Sedimentary rocks are classified into two groups: clastic (the rocks of detrital origin) and
non-clastic (sediments of biochemical or chemical precipitate origin.)
Clastics rocks
Rock type Particle diameter
Conglomerate Pebbles: 2 to 64 mm
Sandstone Sand: 0.06 to 2 mm
Siltstone Silt: 0.003 – 0.06 mm
Shale Clay: < 0.003 mm
Non-clastic
Rock type Composition
Limestone Calcite –CaCo3
Dolomite Dolomite Ca Mg( Co3)
Sandstone Reservoirs
These reservoir rocks consist of quartz (Silica SiO2). These quartz grains cemented
together form sandstone. Sandstones are very often stratified in a superimposed pattern.
This results from successive deposition at the shore-line or in the form of fluvial or
deltaic alluvia. A vertical cross section generally exhibits alternation deposits of sands,
shaly sands, silts and shales.
Sandstone reservoirs are the widest spread hydrocarbon pools.
Carbonate Reservoirs
The carbonate rocks limestone, dolomite, and chalk comprise about 20% of sedimentary
rocks. Limestone composed mainly of the mineral calcite is concentrated by
accumulation of the shells and skeletons of marine animals or by direct precipitation
from mineral saturated waters. Dolomite is the double carbonate of calcium and
magnesium. When dolomitization (replacement of calcium by magnesium) occurs,
shrinkage of matrix is observed. Matrix porosities and permeabilities of carbonate rocks
are typically low. But formation of vugs, channels, and other cavities add to storage
capacity.
Fundamentals of Reservoir Engineering & Characterization 3
The most prolific hydrocarbon bearing carbonates are highly fractured.
TRAPS
The trap is the place where oil and gas are barred from further movement. The
traps can be classified as
• Structural traps: These traps are formed by uplifting and folding of the strata.
When viewed from above, the dome is circular in shape, whereas the anticline is
in an elongated fold.
• Stratigraphic traps: In these traps, trapping is due to variation in facies, The rock
becomes impermeable laterally. Sandstone lenses, pinch outs and carbonate
reefs are some of examples.
• Combination traps: A combination trap has a two or three elements
- a stratigraphic element causing the edge of permeability of the reservoir rock.
- a structural element causing the deformation that combines with the
stratigraphic element to complete rock portion of the trap
- a down dip flow of formation water increasing the trapping effect.
Examples: eroded anticlines, traps associated with salt dome,
Structural trap: oil and gas accumulation in a dome structure
(After “Reservoir and Production Fundamentals”, Schlumberger, 1982)
Fundamentals of Reservoir Engineering & Characterization 4
Structural trap: oil and gas accumulation in an anticline
(After “Reservoir and Production Fundamentals”, Schlumberger, 1982)
Oil accumulation in a stratigraphic trap formed by a change in permeability
(After “Reservoir and Production Fundamentals”, Schlumberger, 1982)
Combination trap
Oil
Fundamentals of Reservoir Engineering & Characterization 5
RESERVOIR PRESSURE
Reservoir pressure is a dominant variable condition that affects every petroleum
reservoir. It is in the form of stored and available energy. It is one of the most important
parameters of reservoir engineering calculations.
The fluids confined in the pores of the reservoir rock occur under certain degree of
pressure, generally called reservoir pressure, fluid pressure or formation pressure. Since
all the fluids are in contact with one another, they transmit pressures freely, and
pressures measured on fluid are actually the pressures on all fluids. Reservoir pressure
unless otherwise stated is generally thought of as the original or virgin pressure – the
pressure that existed before the natural pressure equilibrium of the formation has been
disturbed by any production. The original pressure can be measured directly only by the
first producing the well drilled into the reservoir, for the pressure begins to decline as
soon as oil and gas are withdrawn. When a producing well is shut in, the reservoir
pressure begins to rise. This rise is rapid at first, and then gradually slows until finally the
maximum pressure is reached. The maximum pressure is called the static bottom hole
pressure, the shut in pressure or static formation pressure. The normal pressure
distribution from surface through a reservoir structure is shown
below:
ABNORMAL PRESSURES
Under certain depositional conditions, or because of earth movements, close to reservoir
structure, fluid pressures may depart substantially from the normal range. Abnormal
Fundamentals of Reservoir Engineering & Characterization 6
pressures can occur, when some part of the overburden load is transmitted to the
formation fluids. Abnormal pressures corresponding to gradients of 0.8 psi/ft to 0.9 psi/ft
and approaching geostatic gradient (1.0 psi/ft) can be considered dangerously high.
RESERVOIR TEMPERATURE The
computation of primary recovery of hydrocarbon reservoirs is based on the assumption
that the reservoir temperature remains constant. Thus, hydrocarbon recovery during
primary phase is an isothermal process.
The average reservoir temperature is needed for laboratory analyses carried at reservoir
conditions. Determining reservoir fluid properties such as viscosity, density, formation
volume factor, and gas in solution, and reservoir rock-fluid interaction properties like
capillary, relative permeability and resistivity measurements require a value for reservoir
temperature. For EOR techniques such as chemical and miscible processes,
temperature affects the phase behavior of injected and produced fluids, and thus the
recovery. The feasibility of these processes must be determined by laboratory tests
carried out at reservoir temperature. In EOR processes that employ heat injection, such
as steam or in-situ combustion, the reservoir temperature is not constant and
hydrocarbon recovery is not an isothermal process.
Reservoir temperature is usually measured at the bottom of the well or wells in a
reservoir using a wireline temperature gauge. If a variation in temperature is detected
across a reservoir after correcting for depth, an average value can be used for the
constant reservoir temperature.
Fundamentals of Reservoir Engineering & Characterization 7
TEXTURE
The rock texture is related to those properties of rocks that concerns with grain to grain
relations. Some of these properties are chemical composition, grain shape, grain
roundness, grain size, sorting and grain orientation. The rock texture influence porosity,
permeability, and the interstitial water saturation. Texture is studied by thin section
analysis and visual inspection of hand specimens.
GOOD SORTING POOR SORTING
(After “Fundamentals of Core analysis”, Core Lab, USA, 1989)
Porosity of a rock is the ratio of the pore volume to the bulk volume. In hydrocarbon
reservoirs, the pore volume is the space available for oil, gas and water storage.
Porosity is generally expressed as a percentage of bulk volume.
100xVbVp
=φ
Fundamentals of Reservoir Engineering & Characterization 8
100xVb
VgVb
−=φ
Where Vp = pore volume
Vg = grain volume
Vb = bulk volume
Total or Absolute Porosity: It is the ratio of the volume of all the pores to the bulk
volume of the material, regardless of whether or not, all the pores are interconnected.
Effective porosity: It is the ratio of the interconnected pore volume to the bulk volume
of the rock. The value of this parameter is used in all reservoir engineering calculations.
(After
“Fundamentals of Core analysis”, Core Lab, USA, 1989)
Porosity types
Basic porosity falls into two classes; one that relates to fabric or texture of the rock and
other independent of it. Porosity in sands and sandstone varies primarily with grain size
distribution and grain shapes and packing. Porosity in carbonate rock is much more
variable in magnitude and depends largely on the post depositional processes of
dolomitization, dissolution or cementation.
The porosity types identified in sandstones and carbonate are as follows:
• Fabric related pores are present at time of sediment accumulation and formed
later by fabric controlled.
Fundamentals of Reservoir Engineering & Characterization 9
Fracture porosity both in sandstones and carbonates
These types are common in most of the reservoirs.
Sandstone Reservoirs: There are four basic types of porosity:
• Intergranular porosity: The interstitial pore spaces between the sand grains are
the intergranular porosity which all sandstones possess initially. It ranges from
5% to 40%.
• Intragranular porosity: It is a product of dissolution of soluble material, principally
carbonate particles, unstable rock fragments, feldspar and sulphate within the
formation.
• Microporosity: Microporosity exists as small pores which are commonly
associated with clay minerals
• Fracture porosity: It is generally artificially created in sandstones to improve the
deliverability of any reservoir.
The factors which control sandstone porosity are:
• Mineralogical composition
• Burial history
• Grain size and sorting
• Paleotemperature
• Pressure history
• Pore water composition
• Carbonate cementation
• Secondary porosity
Carbonate Reservoirs:
• Interparticle: Carbonate with a grain supported framework has a large (30%-
40%) initial porosity.
• Intraparticle: These pores are the body cavities which may become sites of
internal sedimentation and crystal filling.
Fundamentals of Reservoir Engineering & Characterization 10
• Moulding: The cavities are formed by solution of shells or destruction of other
original components of the rock, creating moldic porosity.
• Intercrystalline: Coarse dolomites may show intercrystalline porosity caused by
solution of non-replaced calcite.
• Fracture porosity
Determination of porosity
The porosity is determined by core analysis or by well logging.
Core analysis
In porosity any two of Vp, Vb, Vg are determined. In core analysis, the cylindrical plugs
of either 1.0 inch or 1.5 inch diameter are cut from whole core and then first cleaned and
dried.
Measurement of bulk volume
• Caliper method. The length and diameter of core plug is measured at different points
of the core and averaged values are determined.
Vb = 4
2ldπ
• Measurement of the buoyancy exerted by mercury on the samples immersed
in it.
The mercury based methods are not used for rocks containing fissures or macropores
because of possibility of mercury penetration.
Measurement of pore volume
The pore volume can be measured:
• Helium expansion in the interconnected pores
• Measurement by weighing in a fluid filling the effective pores
• Measurement by mercury injection
The grain volume can also be determined by Helium expansion method.
Effect of pressure on porosity
Porosity decreases with increasing net overburden pressure. Reservoir rocks experience
the lithostatic pressure and fluids pressure in the pores. The production of hydrocarbons
Fundamentals of Reservoir Engineering & Characterization 11
causes a decline in the fluid pressure in the pores resulting in compression of the rock,
until a new equilibrium is attained.
Averaging of porosity
Arithmetic averaging of thickness average porosity: This method is used in cases
when the reservoir rock shows large variation on porosity vertically but does not show
great variations in porosity parallel to the bedding planes.
Arithmetic Average porosity Ø = njφ
Thickness weighted porosity Ø = J
jj
h
hφ
Areal weighted or volumetric weighted average porosity: These averages are used
in cases where the porosity in one portion of the reservoir is greatly different from that in
another area because of sedimentation or depositional changes.
Areal weighted average porosity Ø = j
jj
A
Aφ
Volumetric weighted average porosity Ø = jj
jjj
hA
hAφ
Where n = total number of core samples
hj = thickness of core sample j or reservoir area j
Øj = porosity of core sample j or reservoir area j
Aj = reservoir area j
GRAIN DENSITY
The grain density of a rock is defined as the weight of the rock (exclusive of the weight of
fluids contained in the pore space) divided by the volume of the solid rock material
(exclusive of pore space). The density varies with the mineral composition of the rock
and the state of hydration of the minerals. In complex lithologies containing inter-mixed
limestone, dolomite, sandstones, and heavy minerals, grain density will vary vertically
and horizontally. Even in formations described as homogeneous, measured densities
often vary considerably from published values for pure components as tabulated below.
Fundamentals of Reservoir Engineering & Characterization 12
Minor amounts of secondary cement, such as calcite or siderite, will cause grain
densities to exceed values shown in the table.
Component Approximate grain density (g/cm3)
Sandstone 2.65
Limestone 2.71
Dolomite 2.85-2.87
Anhydrite 2.98
Gypsum 2.3
Pyrite 5.0
Siderite 3.9
Clays 2.2-2.9
Grain density is important in core analysis on the account that it can be used as a quality
control check of the core analysis measurements themselves.
PERMEABILITY
Permeability is a measure of the capacity of formation to transmit fluids. Its unit is Darcy,
named after a French scientist Henry Darcy in 1856. One Darcy equals permeability that
will permit a fluid of one centipoise viscosity to flow at a rate of one cubic centimeter per
second through a cross-sectional area of one square centimeter when the pressure
gradient is one atmosphere per centimeter. Generally permeabilities are given in
millidarcies which is equal to (1/1000) of a Darcy. Its dimension is L2.
K A ∆P ∆∆∆∆P = Press. Differential, atm q = ------------ A = Cross Sectional Area, cm2 µ * L K = Permeability, darcy q = Outlet Flow Rate, cc/sec µ = Fluid Viscosity, cp L = System Length, cm
∆∆∆∆ P
q A
L
Fundamentals of Reservoir Engineering & Characterization 13
Darcy law is used to determine permeability when the following conditions exist:
• Laminar flow
• No reaction between fluid and rock
• One phase present at 100 percent pore space saturation.
The measured permeability at 100% saturation of a single phase is called the absolute
permeability of the rock.
The following terms are generally used to specify the permeability:
<1mD = Very low
1to 10 mD = Low
10 to 50 mD = Medium
50 to 200 mD =Good
200 to 500 mD = Very Good
>500 mD = Excellent
The factors which control magnitude of permeability are:
• Shape and size of sand grains
• Lamination
• Cementation
• Fracturing and solution
Permeability Anisotropy
Permeability is a directional quantity. The long axis of the grains aligns parallel in the
direction of maximum velocity during the process of sediments deposition, thus providing
the maximum cross-sectional area of the grains in a horizontal plane. This results in
highest permeability parallel to long axis of the grains.
In most of reservoir rocks, permeability like porosity is reduced by increase in net
overburden pressure.
Measurement of Permeability The permeability is measured by flowing a fluid of known
viscosity µ through a core plug of measured dimensions (A and L) and then measuring
flow rate q and pressure drop p. Darcy equation becomes
pALq
k∆
= µ
Fundamentals of Reservoir Engineering & Characterization 14
Absolute permeability is usually determined by flowing air through the core plug because
of its convenience and to minimize rock-fluid interaction.
In using dry gas in measuring the permeability, the gas volumetric rate q varies with the
pressure because the gas is a highly compressible fluid. Hence, the equation becomes
bgsc Lp
ppkAQ
µ2)( 2
22
1 −=
Where k = absolute permeability, Darcies
µg = gas viscosity, cp
pb = base pressure ( atmospheric pressure), atm
p1 = inlet pressure (upstream), atm.
p2 = outlet (down stream), atm.
L = length of the core plug, cm
A = cross-sectional area, cm2
Qsc = gas flow rate at standard conditions, cm3/sec.
Klinkenberg effect
Klikenberg (1941) compared the permeability results of measurements made with air as
the flowing fluid as well as with a liquid as the flowing fluid. He observed that the air
permeability is always greater than the liquid permeability. Klinkenberg postulated that
liquids had a zero velocity at the sand grain surface while gases exhibited some finite
velocity at the sand grain surface. And this slippage at the sand grain surface has
resulted in higher flow rate for the gas at a given pressure differential. Further, he also
found that as the mean pressure increased, the calculated permeability of the porous
medium decreased. The magnitude of Klinkenberg effect varies with the core
permeability and the type of gas used in the experiment. The resulting straight
relationship can be expressed as:
Ka = KL + b[1/pm]
Where Ka = measured gas permeability
pm = mean pressure
KL = equivalent liquid permeability
b = slope of line
Further b = c KL where c is a constant which depends on the size of the pore openings
and is inversely proportional to the radius of capillaries.
Fundamentals of Reservoir Engineering & Characterization 15
Klinkenberg effect
A comparison of absolute permeability and Klinkenberg permeability is given below:
Gas Permeability, mD (Ka)
Klinkenberg Permeability, mD (KL)
Ratio of KL/ Ka
0.18 0.12 0.66 1.00 0.68 0.68 10.0 7.80 0.78
100.0 88.0 0.88 1000.0 950.0 0.95
Averaging of absolute permeabilities
An adequate understanding of permeability distribution is critical to the reservoir
performance prediction. Homogeneous reservoirs seldom exist. Because of existence of
small scale heterogeneities, laboratory measured core plug permeabilities needs proper
averaging for flow characteristics representation of the entire reservoir or its individual
reservoir units.
There are three commonly used techniques:
• Weighted average permeability
• Harmonic average permeability
• Geometric average permeability
Weighted Average Permeability: Used to determine the average permeability of
layered – parallel beds with different permeabilities.
Gas
Per
mea
bilit
y
Liquid permeability
Fundamentals of Reservoir Engineering & Characterization 16
=
==n
jj
n
jjj
avg
h
hkk
1
1
Where hj = thickness of layer j
Kj = absolute permeability of layer j
Harmonic Average Permeability: Used to average permeabities where permeability
variations can occur laterally in a reservoir.
=
=
=
n
j j
n
jj
av
kL
Lk
1
1
Where Lj = length of each bed
kj = absolute permeability of each bed
Geometric Average Permeability: Most representative averaging technique for a
heterogeneous formation:
=
=
=n
jj
n
jjj
avg
h
khk
1
1
))ln((exp
Where kj = permeability of core sample j
hi = thickness of core sample j
n = total number of samples
If the thickness of all the core samples is same, then the above equation becomes:
( )nnavg kkkkk1
321 ...=
SATURATION
Fluid saturation is defined as the fraction of pore volume occupied by a particular fluid.
Hence for reservoir fluids, mathematical expressions can be:
Oil saturation, VolumePore
oilofVolumeSo ..
....=
Gas saturation, VolumePore
gasofVolumeS g ..
....=
Fundamentals of Reservoir Engineering & Characterization 17
Water Saturation, VolumePore
waterofVolumeSw ..
....=
Sg + So + Sw = 1.0
Determination of saturation
Fluid saturation in the laboratory is one of the least reliable reservoir property
measurements. Factors that are likely to introduce errors into these measurements
include invasion of the core by mud or mud filtrate during coring process, gas expansion
during core recovery, and handling of the core during preservation and measurement.
Some of the methods generally used for laboratory determination of fluid saturations are:
Soxhlet distillation extraction/Dean-Stark method: In this method, oil is removed from
the sample by extraction i.e. dissolved in suitable solvent; most commonly used toluene
and xylene. A mixture of 80%acetone+20% methanol is frequently employed. Water is
removed from the sample by distillation, then condensed to liquid which is caught in a
trap and measured.
Retort method. It is atmospheric distillation in which rock sample is heated in stages to
1200oF. All the reservoir fluids are vaporized. The most commonly used system employs
electric heating and counter-current cooling with water.
Averaging of saturation data
The representative averaging of saturation data requires that the saturation values be
weighted by both the interval thickness hj and interval porosity øj
=
==n
jjj
n
jjjj
h
SohSo
1
1
φ
φ
=
==n
jjj
n
jjjj
h
SwhSw
1
1
φ
φ
Fundamentals of Reservoir Engineering & Characterization 18
=
==n
jjj
n
jjjj
h
SghSg
1
1
φ
φ
Where the subscript j refers any individual measurement and hj represents the depth
interval to which Øj, Soj, Sgj, Swj apply.
WETTABILITY
Wettability is defined as the tendency of one fluid to adhere or spread on a solid surface
in presence of other immiscible fluids. The varying wetting characteristics of liquids for
the solid can be observed by placing small drop of three liquids namely mercury, oil and
water on clean glass plate.
The spreading tendency is expressed by measuring the angle of contact at the liquid-
solid interface. This angle is called the contact angle . Wettablity can be determined in
the laboratory by measuring the contact angle between a droplet of fluid and a flat
surface of mineral crystal. Wettability has profound influence on distribution of fluids in
the porous media and affects the ultimate recovery. Because of the attractive forces, the
wetting phase tends to occupy the smaller pores of the rock and non-wetting phase
occupies the more open channels.
In reservoirs, generally water is considered to be wetting fluid. However, the oil may be
wetting especially for limestones.
The laboratory studies have indicated that preferentially wettability of the rock is largely
controlled by the compounds adsorbed at the surface of the rock.
CAPILLARY PRESSURE
Surface and interfacial tension result from molecular forces that cause the surface of a
liquid to assume the smallest possible size and to act like a membrane under tension.
Fundamentals of Reservoir Engineering & Characterization 19
Capillarity is the rise or depression of liquids in a fine tube resulting from surface tension
and wetting preferences. Consider a capillary tube of radius ‘r’ placed in large open
vessel containing water. The water will rise in tube, until the total force acting to pull
liquid upward is balanced by the weight of the column of liquid being supported in the
tube.
Fup = 2r.
gw. cos
Fdown = r2.h. (ρw - ρair).g
Since density of air is negligible in comparison to density of water
Fdown = r2.h.ρw.g
At equilibrium Fup = Fdown
2r. gw. cos = r2.h.ρw.g
gw = θρ
cos2... wghr
In porous medium, even when two or more fluids are present at the same subsea
elevation and are in state of pressure equilibrium, they are not at the same pressure.
This is primarily because of differences in the mutual attraction between rock and fluids
(adhesion tension). This difference in pressure between two phases in equilibrium at the
same subsea elevation is the capillary pressure between two phases. The fluid with the
greatest tendency to wet the rock will have the lowest pressure.
Fundamentals of Reservoir Engineering & Characterization 20
Pc = pnw - pw
The pressure excess in the non-wetting fluid is the capillary pressure, and this quantity is
a function of saturation.
Gas - Liquid system
Pc = r
gw θσ cos2
h = )(
cos2
gw
gw
rg ρρθσ
−
Oil - Water System
Pc = r
ow θσ cos2
h = )(
cos2
ow
ow
rg ρρθσ
−
Where ρw = water density, gm/cm3
ρ0 = oil density, gm/cm3
gw = gas-water surface tension, dynes/cm
ow = oil-water surface tension, dynes/cm
r = capillary radius, cm
= contact angle
h = capillary rise, cm
g =- acceleration due to gravity, cm/sec2
Pc = capillary pressure, dynes/cm2
Laboratory determination of capillary pressure data
Three methods are generally used for determination of capillary data on rock samples:
• Purcell’s method/Mercury injection method: The core plug cleaned and dried
with pore volume determined is first subjected to vacuum after pore volume
determination. The mercury is injected into it in increasing pressure stages. At
each stage, the volume of mercury intruded is recorded. The capillary pressure
is the absolute pressure of mercury.
• Restored state method: The core plug saturated with brine in placed on a porous
plate saturated with brine. Air is injected at increasing pressure stages. A
Fundamentals of Reservoir Engineering & Characterization 21
capillary tube is used to measure the volume of water expelled from the core.
The capillary pressure is the relative pressure of the air.
• Centrifuge method: In the centrifuge method, an artificial gravity using the density
difference between the two fluids creates a capillary pressure gradient all along the
plug and thus a saturation variation from the top to the bottom.
Converting laboratory capillary pressure data to reservoir conditions
Since the laboratory measurements are not conducted using reservoir fluids, the lab
results must be corrected to reservoir condition using the relationship:
PcRes = PcLab(Res/ Lab)
Where
PcRes = capillary pressure at reservoir conditions, psi
PcLab = capillary pressure at laboratory conditions, psi
Res = interfacial tension at reservoir conditions, dynes/cm
Lab = interfacial tension at laboratory conditions, dynes/cm
Averaging of capillary pressure data
Leverett (1942) proposed a means of converting all capillary – pressure data to a
universal curve using the dimensionless function of saturation known as J-
function,
φσkPc
J Sw 21645.0)( =
Where J(sw) = Leverett J-function
Pc = capillary pressyre, psi
= interfacial tension dynes/cm
k = permeability, mD
Ø= porosity, fraction
Each capillary pressure curve gives a J-function curve. The average J-curve is
generated. Using this curve, a Pc-Sw can be plotted for a given sample if its k and Ø are
known.
Fundamentals of Reservoir Engineering & Characterization 22
.
The Leverett J-Function for unconsolidated sands (After Leverett,1941)
Initial saturation distribution in a reservoir
An important application of capillary pressure data relates to the fluid distribution in a
virgin reservoir. The capillary pressure - saturation data can be converted into height –
saturation relation ship as given below:
h = ρ∆Pc144
where Pc = capillary pressure,psia
ρ = density difference between wetting phase and non-wetting phase at
reservoir conditions, lb/ft3
h = height above the free water level, ft
Fundamentals of Reservoir Engineering & Characterization 23
Distribution of saturation in the reservoir
(After “Reservoir and Production Fundamentals”, Schlumberger, 1982)
The transition is the vertical thickness over which the water saturation changes 100%
saturation to irreducible water saturation, Swi.
The water oil contact is the uppermost depth in the reservoir where a 100% water
saturation exists. At free water level, there is zero capillary pressure from reservoir
engineering standpoint.
Irreducible water saturation
Irreducible water saturation is the minimum saturation that can be induced by
displacement. At this stage, the wetting phase becomes discontinuous. This minimum
saturation corresponds to smallest mean radius of curvature and maximum capillary
pressure.
Grain size has remarkable influence on irreducuible water saturation:
Cap
illar
y P
ress
ure
H
eigh
t Abo
ve O
il-W
ater
Con
tact
Fundamentals of Reservoir Engineering & Characterization 24
Effect of grain size on Irreducible water saturation
(After “Reservoir and Production Fundamentals”, Schlumberger, 1982)
RELATIVE PERMEABILITY
Production of hydrocarbons involves simultaneous flow of two or three fluids in the
reservoir rock. In this multiphase flow, each fluid tends to interfere with the flow of the
others.
Absolute permeability relates to permeability with one fluid present at 100 percent
saturation. It is also called as specific permeability or base permeability.
Effective permeability is the permeability to a given phase when more than one phase
saturates the porous medium. The effective permeability is a function of saturation.
Relative permeability to a given phase is defined as the ratio of effective
permeability to the absolute or, in some cases, a base permeability. Relative
permeability is also a function of saturation.
Relative permeability = typermeabiliBase
typermeabiliEffective⋅
.
Fundamentals of Reservoir Engineering & Characterization 25
kro =kko
krg =k
k g
krw =k
kw
It is a dimensionless term and generally reported in fraction or percentage.
For an oil-water reservoir, the base permeability, k is taken as effective permeability to
oil at irreducible water saturation. For a gas reservoir, the base permeability will be
effective permeability to gas in the presence of irreducible water.
Imbibition versus Drainage
In relative permeability studies, the terms imbibition and drainage are commonly
referred.
If the wetting phase is decreasing, that phase is draining and the curve is called a
drainage curve.
If the wetting phase is increasing or being imbibed during the test, the curve is referred
to as an imbibition curve.
Water – Oil relative permeability curve
leaving In oil-water system, oil and water relative permeabilities are plotted as functions
of water saturation. At irreducible water saturation, Swi, the relative permeability to
water, krw is zero and oil permeability with respect to oil kro is a value less than unity. This
is due to reduction in oil permeability due to presence of water. At Swi, only oil can flow.
As water saturation increases, relative permeability to water increases and oil
permeabity decreases. The maximum water saturation is reached at the residual oil
saturation (Sor). Residual oil saturation is left in the smaller channels when the
interfacial tension causes the thread of oil to break; behind oil in droplets which tend to
assume spherical form and when gradient pressure is not sufficient to deform the bubble
enough to pass through the smaller pore openings.
Fundamentals of Reservoir Engineering & Characterization 26
0
0 Swi Sor 1
Sw
Gas–Oil relative permeability curve
The gas relative permeability krg remains zero until the critical gas saturation Sgc is
reached. At Sgc, there is enough accumulation of gas for its mobility. As gas saturation
increases, the gas relative permeability increases. The gas relative permeability will
achieve maximum value at residual oil saturation. The oil permeability decreases from
unity to lesser values as the gas saturation increases finally reaching a value of zero at
the residual oil saturation plus irreducible water saturation.
0
0 Sgc Sorg 1 Sg
Laboratory methods for measuring relative permeability
Two major laboratory methods have evolved to measure relative permeability. These are
referred to as the steady-state and nonsteady-state techniques.
krel
1
krel
1
Fundamentals of Reservoir Engineering & Characterization 27
STEADY STATE: The steady-state test, the older of the two methods, is made at low
flow rates. Most research groups prefer data obtained from this test. Two fluids are
injected simultaneously into a core sample and the water saturation is increased slowly.
This simulates the slow increase in water saturation that would occur in the formation
between the injection and producing wells. Saturation increase is monitored by
measuring the gain in weight occurring in the sample or by X-ray technique.
NONSTEADY STATE: The nonsteady-state technique uses viscous oil and is normally
made at a higher flow rate than that present in the reservoir. It is this higher rate that
sometimes yields pessimistic estimates of recovery from rocks of intermediate
wettability.
Normalisation and averaging of relative permeability data
There is a wide variation observed in relative permeability results of experiments
conducted on core plugs of a reservoir rock. To use this data for reservoir engineering
calculations, the proper averaging or normalization of relative permeability data obtained
on individual rock samples is essential so that the effects of different water saturations
and residual oil saturations are removed. The normalized relative permeability data is
then denormalised for different portions of reservoir as per the measured irreducible
water saturation and residual oil saturation.
The following steps are required for averaging of oil and water relative permeability
curves.
1. Starting with Swi, chose several values of Sw and calculate Sw* for each set of
relative permeability curve
Sw* = orwi
wiw
SSSS−−
−1
Calculate the normalized relative permeability for the oil phase at different water
saturation
kro* = Swiro
ro
kk
)(
3. Calculate the normalized relative permeability of the water phase at different
water saturation
Fundamentals of Reservoir Engineering & Characterization 28
krw*= orSrw
rw
kk
)(
4. Make a linear plot of the normalized kro*, krw* versus Sw* for all the core samples.
Obtain a single pair of normalized relative permeability curve by selecting
arbitrary values of Sw* and calculate the average of kro* and krw* using the
following relationships
(kro*)avg =
=
=n
jj
n
jjro
hk
hkk
1
1
)(
*)(
(krw*)avg =
=
=n
jj
n
ijrw
hk
hkk
1
1
)(
*)(
Where n = total number of core samples
hj = thickness of sample j
kj = absolute permeability of sample j
6. Denormalise the average curve to reflect actual reservoir and conditions of Swi
and Sor; using the following equations:
wiorwww SSSSS +−−= )1(*
( k ro)Swi =
[ ]
=
=n
jj
n
jjSwiro
hk
khk
1
1
)(
)(
( k rw)Sor =
[ ]
=
=n
jj
n
jjSorrw
hk
khk
1
1
)(
)(
Where (kro)Swi and (krw)Sor are the average relative permeability of oil and water at
irreducible water saturation and residual oil saturation respectively.
WELL LOGGING
A well log is the continuous recording of the characteristics of the hole drilled formation,
as a function of depth.
Fundamentals of Reservoir Engineering & Characterization 29
Well logs are recorded at the various stages in well under drilling. The drilling is
interrupted during the log recording. The data is recorded and transmitted to the surface
instantaneously. Well logs are essential tools for enhanced reservoir evaluation.
Electric Logs
Spontaneous potential
The SP log is the difference in electric potential between a fixed electrode at the surface
and a moving electrode in the borehole. It is measured in millivolts, and there is no
absolute zero; only changes in potential are recorded.
Two types of potential may contribute to the SP effect. These are the electrochemical
potential (Ec) and the electro kinetic potential (Ek).
The electro kinetic potential (Ek) is produced by the flow of mud filtrate through a porous
and permeable formation. The electrochemical potential (Ec) results from the transfer of
ions from a more concentrated electrolyte (usually the uninvaded zone in the formation)
to a less concentrated electrolyte (usually mud in the bore hole).
The SP log is used in the identification of permeable beds and the location of their
boundaries, and for determination of formation water resistivity in the uninvaded zone
(Rw).
A deflection is observed opposite the reservoir rock compared with a “base line” of the
shale. These deflections are due to different salinities of the reservoir water and the
drilling mud.
Fundamentals of Reservoir Engineering & Characterization 30
Resistivity log
Resistivity logs measure and record the resistance offered by the rocks surrounding the
bore hole to the passage of the electric current. A system of electrodes sends an electric
current into the formation. The apparent resistivity of the reservoir is measured in ohms
per meter. The resistivity logs may be divided into conventional or non focused devices,
focused tools and induction systems.
The Laterolog systems contain an array of electrodes to focus the survey current and
force it to flow laterally into the formations surrounding the borehole. The effective depth
of laterolog investigation is controlled by the extent to which the surveying current is
focused.
The induction log measures the conductivity of the rocks surrounding the borehole by
inducing an electric current through them. The tool consists of a transmitter and a
receiver coil. A constant, high frequency alternating current is sent through the
transmitter coils. This generates an alternating magnetic field which induces secondary
currents (also known as eddy currents) in the rocks surrounding the borehole. These
currents flow in circular paths coaxial with the transmitter coils through the surrounding
rocks. The resulting magnetic field, in turn induces signals in the receiver coils. These
signals are proportional to the conductivity of the formations from which resistivity is
derived and recorded on the log.
The resistivity recorded is a function of the porosity and saturation (water/hydrocarbons).
The rock matrices are insulating and the hydrocarbons have high resistivity, whereas the
resistivity of the water decreases with increasing salinity. The resistivity can differentiate
the water from hydrocarbons.
Empirical equations
m
aRwRo
Fφ
==
n
RtRw
S w φ1=
Fundamentals of Reservoir Engineering & Characterization 31
Where
Ro = resistivity of rocks 100% saturated with water of resistivity Rw.
F = formation factor
a = tortuosity coefficient
m = cementation factor
n = saturation exponent
Rt = calculated resistivity of rock at water saturation Sw.
Radioactivity Logs
Gamma ray log (GR)
This log records the natural radioactivity of formations. The radioactivity arises from the
presence of uranium(U), thorium(Th) and potassium (K40) in the rocks. These elements
continuously emit gamma rays, which are short bursts of high energy radiation similar to
x-rays. Gamma rays are capable of penetrating a few inches of rock, and a fraction of
Fundamentals of Reservoir Engineering & Characterization 32
those that originate close to the borehole traverse the hole and can be detected by a
suitable gamma-ray sensor. The detector gives a discrete electrical pulse for each
gamma ray detected, and the parameter logged is the number of pulses recorded per
unit of time by the detector. The GR log is useful in detecting shale beds. Non
radioactive minerals like coal may be detected by their characteristically low gamma
response. This log is used for correlation of formations in cased holes.
Neutron log
In neutron logging the formations surrounding the borehole are bombarded by high
energy neutrons from an artificial source carried on the device. Neutrons are electrically
neutral particles with a mass almost identical to that of a hydrogen atom. Upon leaving
the source the neutrons enter the formations and collide with nuclei in the rocks forming
the borehole wall. With each collision a neutron loses some of its energy. The amount of
energy lost per collision depends on the relative mass of the nucleus with which the
neutron collides. The greatest energy loss occurs when the neutron strikes a nucleus of
practically equal mass, ie. a hydrogen nucleus. Collisions with heavy nuclei do not slow
the neutron down very much. Thus the slowing down of neutrons depends largely on the
amount of hydrogen in the formation. The sonde emits fast neutrons which bombard the
formation giving rise to slow neutrons, The neutron count rates increase with decreasing
hydrogen content (low porosity in clean formations) and decrease with increasing
hydrogen content (high porosity in clean formations).
Formation Density Compensated (FDC) Log
The Formation Density Compensated (FDC) Log records the bulk density (b) of the
formation surrounding the borehole. Gamma rays are beamed at the formations by the
source. These enter the formations and undergo multiple collisions with the electrons in
the frocks, as the result of which they energy and become scattered in all directions.
This is known as Compton scattering. Some of the scattered gamma rays return to the
borehole and are recorded by the detectors on the device. The intensity of the returned
radiation is proportional to the number of electrons in the formation, and provides a
measure of the electron density of the material. Electron density is approximately is
equal to the bulk density of the rocks and this is recorded in gm/cm3.
Fundamentals of Reservoir Engineering & Characterization 33
mf DDD )1(. φφ −+=
Where D = total density read on log
Df = fluid density (filtrate)
Dm = density of rock matrix
The Borehole Compensated Sonic Log (BHC)
The sonic or acoustic log provides a continuous record of the time taken, in milliseconds
per foot (µsec/ft), by a compressional sound wave to travel through one foot of formation.
Known as the interval transit time, this is the reciprocal of the compressional wave
velocity (Vp).
The velocity of sound through a given formation is a function of its lithology and porosity.
Dense, low porosity rocks are characterized by high matrix velocities (Vm), while porous
and less dense formations are characterized by low Vm, values.
Fundamentals of Reservoir Engineering & Characterization 34
mf VVVφφ −+= 11
or
tmtt f ∆−+∆=∆ ).1(. φφ
Where t = travel time in the transmitter/receiver interval
Some other logs
Caliper log
This log system with arms furnishes the borehole diameter and helps in identifications of
caving, constrictions etc.
Dipmeter log
This is the simultaneous recording of four microlaterolog curves along four 90 degree
generating lines in a plane normal to the bore hole axis. The difference in the four curves
gives the value of dip and its direction.
Fundamentals of Reservoir Engineering & Characterization 35
Cement bond log (CBL)
This log system is a continuous cased hole recording of the amplitude of the acoustic
signal versus depth. The analysis of signals provides information on the presence and
bonding the cement to the casing and to the formation
• In the presence of cement, the signal is weak because cement attenuates the
vibrations of the metal.
• In the absence of cement, the casing vibrates freely generating a strong signal.
Cement Bond Log
Fundamentals of Reservoir Engineering & Characterization 36
Introduction The chemistry of hydrocarbon reservoir fluids is very complex. Some estimates suggest
that perhaps 3,000 organic compounds can exist in a single reservoir fluid. These
compounds contain a variety of substance of diverse chemical nature that includes
hydrocarbons and non hydrocarbons. Hydrocarbons range from methane to substances
that may contain more than 100 carbon atoms. Non-hydrocarbons include substances
such as N2, CO2, H2S, S, H2O, He and even traces of Hg, etc.
The physical properties of these mixtures depend primarily on composition and
temperature & pressure conditions. Reservoir temperature usually can be assumed
constant, however as the oil and gas are produced, reservoir pressure decreases and
the remaining hydrocarbon mixtures change in composition, volumetric properties, and
phase behaviour. Understanding of this behaviour is very important for a petroleum
engineer as it is of prime consideration in the development and management of
reservoirs that would maximize the profits.
With this objective, this particular chapter on reservoir fluid behaviour would familiarize
the reader about reservoir fluid composition, phase behaviour properties, types of
reservoir fluid, various reservoir fluid characteristics and empirical methods for its
determination, various types of laboratory experiment, application of reservoir fluid
characteristics and equation of state. At the end of this chapter reader should be able to
apply these concepts in solving practical engineering problems.
Reservoir Fluid Composition
The empirical formula CnH2n+hSaNbOc can be used to classify nearly all compounds found
in crude oil. The largest portion of crude oil is composed of hydrocarbons with carbon
number n, ranging from 1 to about 60, and h numbers ranging from +2 for low molecular
weight paraffin hydrocarbons to -20 for high molecular-weight organic compounds.
Fundamentals of Reservoir Engineering & Characterization 37
Occasionally, sulphur, nitrogen and oxygen substitutions occur in high molecular weight
organic compounds with a, b and c usually ranging from 1 to 3.
Those hydrocarbons which contain only two elements, hydrogen and carbon are of two
types aliphatic and aromatic. Aliphatic hydrocarbons are further divided into alkanes
(CnH2n+2), alkenes (CnH2n), alkynes (CnH2n-2), and their cyclic analogs.
The series of straight chain alkanes show a smooth gradation of physical properties. As
molecular size increases, each additional CH2 group contributes a fairly constant
increment to boiling point and specific gravity. The boiling and melting points of alkanes
are fairly low because of symmetrical nature of molecules. Chemically, alkanes are
unreactive at ordinary temperature. Hence, naturally occurring petroleum deposits
mainly consist of alkanes.
The physical properties of alkenes and alkynes are very much like the physical
properties of alkanes. However, because of double and triple bonds, alkenes and
alkynes are more reactive than alkanes. Hence, alkenes and alkynes are not usually
found in naturally occurring hydrocarbon deposits.
Cycloalkanes and cycloalkenes are about as reactive chemically as their open chain
analogs. Different members of cyclic group exhibit different chemical reactivities.
Aromatic hydrocarbons show gradation in their physical properties with increase in
molecular weight and they have the same stability as the carbon-carbon single bond
found in alkanes.
There are many families of organic compounds other than alkanes, alkenes, alkynes and
their cyclic analogs which, contains atoms other than carbon and hydrogen e.g. sulphur,
nitrogen and oxygen etc. Mercaptans, alkyl sulphides, aldehydes, ketones, resins and
asphaltenes belong to this category of organic compounds.
Following table lists classification of organic compounds according to functional groups.
Fundamentals of Reservoir Engineering & Characterization 38
Phase Behaviour of a Multicomponent Mixture The phase behaviour of a multi component mixture is not as simple as that of a pure
component. It is more elaborate than that of a pure component. The complexity
compounds as component with widely different structures and molecular sizes comprise
the system. However, reservoir fluids are mainly composed of hydrocarbons with similar
structures. Hence their phase behaviour is not generally complex. Two important
differences between pure and multicomponent systems are (i) the saturated P-T diagram
is represented by a phase envelope rather than by a vapor-pressure curve as the
separation between bubble point and dew point increases with the contrast of system
component, and (ii) the critical temperature and critical pressure no longer define the
extent of the two phase region. Two phase can exist upto cricondentherm and
cricondenbar beyond critical temperature and critical pressure.
Fundamentals of Reservoir Engineering & Characterization 47
Pure Component Multicomponent
Pressure-volume diagram of a multicomponent reservoir fluid is schematically shown
below;
Contrary to a pure system, in a multicomponent system the system pressure decreases
during an isothermal expansion between its bubble and dew points. At the bubble point
(A), the composition of the liquid is essentially equal to the overall composition of the
mixture. However, the infinitesimal amount of gas which is liberated is richer in the more
volatile component. Similarly at the dew point the composition of the vapour is
essentially equal to the over all composition of the mixture with infinitesimal amount of
liquid is richer in the least volatile component.
Critical Point
Temperature
Liquid Solid
Pressure
Vapour
C
Two Phase
Pressure
Temperature
T2
T1
T3
A B
Pressure
Volume
Fundamentals of Reservoir Engineering & Characterization 48
Phase diagram of a mixture is determined by its composition. Figure shown below is that
of ethane-heptane system. The critical temperature of different mixture lies between the
critical temperature of the two pure compounds. However, the critical pressure exceeds
the value of both components as pure, in most cases.
Fundamentals of Reservoir Engineering & Characterization 49
The greater the difference between the critical points of the two components, the higher
the mixture critical pressure would be.
Retrograde Condensation: In a multicomponent phase diagram as shown below,
vapour and liquid phases coexist at any pressure-temperature condition within the phase
envelope. The different liquid/mixture volumetric ratios are conventionally shown as
dashed lines which are called quality lines. The quality lines come very close towards
each other near critical point of the system. Hence small pressure or temperature
changes at a region near the critical point cause major phase changes.
If the reservoir hydrocarbon system is at point A, reduction of pressure for vapor like fluid
at point A, forms the first drop of liquid at point B. Further reduction of pressure will result
in further condensation, as indicated by quality lines. This phenomenon of condensation
with decrease in pressure is called retrograde condensation. The condensation will
cease at some point, point D, and the condensed phase will re-vaporize again. The
shaded region of the phase diagram is called retrograde region. This is an important
phenomenon which is generally observed in gas condensate wells.
Classification of Reservoirs and Reservoir Fluids A typical phase diagram of a reservoir hydrocarbon system can be used to describe
various types of reservoir fluids. Identification of types of reservoir fluids is necessary
and must for production and reservoir engineer, as different types of fluid require
different approaches for exploitation.
How to classify reservoir types? Location of reservoir temperature on the phase diagram
can be used to classify reservoir fluids. There are five types of reservoir fluids; dry gas,
wet gas, gas condensate (retrograde gas), volatile oil and black oil.
Fundamentals of Reservoir Engineering & Characterization 50
Dry Gas Dry gas is primarily composed of methane and some intermediates such as nitrogen and
carbon dioxide. A typical phase diagram of a dry gas is given below
As evident from phase diagram the phase envelop is relatively tight and mostly located
below the ambient temperature. A dry gas does not contain any enough of heavier
molecule to form hydrocarbon liquid at the surface. Gas remains in single phase from
reservoir to separator. Water, however may condense at surface condition due to
cooling effect.
Wet Gas A wet gas exists solely as gas in the reservoir throughout the reduction in reservoir
pressure. However, liquid may form at separator due to its position within the phase
region. A typical phase diagram of wet gas system is given below;
Pressure
Temperature
Reservoir
Separator
Reservoir
Separator
C
Fundamentals of Reservoir Engineering & Characterization 51
The surface liquid is normally called as condensate. As no condensate is formed in the
reservoir, material balance equation for a dry gas is equally suitable for a wet gas.
Producing gas to condensate ratios are typically above 10,000 v/v.
Gas Condensate or Retrograde Gas In a typical gas condensate reservoir the reservoir temperature lies between critical point
and cricondentherm. The gas will drop out liquid by liquid by retrograde condensation in
the reservoir, when the pressure falls below the dew point. A typical gas condensate
phase diagram is shown below;
The phase diagram of a retrograde gas is somewhat smaller than that for oils because of
presence of less heavy hydrocarbons. However, presence of heavy hydrocarbons
(compared to a wet gas system) expands the phase envelope relative to a wet gas
phase envelope. Material balance equation developed for dry gases can be used for a
gas condensate reservoir as long as the reservoir pressure remains above the dew point.
Lowering of reservoir pressure, below dew point, results in gas to form free liquid in the
reservoir. The liquid will normally not flow and can not be produced. Hence,
condensation and loss of valuable compounds in reservoirs could be avoided by
maintaining the reservoir pressure above the fluid dew point by gas recycling. A
compositional material balance method should be used for gas condensate reservoir
system where pressure has fallen below dew point.
Volatile Oil Volatile oil contains relatively higher heavy molecules than a gas condensate system
that makes it to behave liquid-like at reservoir conditions. The phase envelope, as per
phase rule, is relatively wider than that of a gas condensate system, with a higher critical
temperature due to its larger concentration of heavy compounds.
Pressure
Temperature
C
Reservoir
Separator
Fundamentals of Reservoir Engineering & Characterization 52
A typical volatile oil phase diagram is shown below;
The vertical line shows the path taken by the isothermal pressure reduction during
depletion. A small reduction in the pressure below the bubble point causes the release of
a large amount of gas in the reservoir. Saturation pressure of volatile oils is high. Gases
produced below the bubble point, therefore, are quite rich and behave as retrograde
gases. The amount of liquid recovered from the gas constitutes a significant portion of
the total oil recovery. Compositional methods should be applied generally to study
volatile oil reservoirs.
Black Oil Black oil is the most common type of oil reserves. It consists of wide variety of chemical
species including large, heavy, nonvolatile molecules. Therefore, its phase envelope is
the widest of all the types of reservoir fluids, with its critical temperature well above the
reservoir temperature. A typical phase diagram of black oil is shown below
Reserv
2
3
Separator
C 1
Pressure
Temperatu
Separato
Reservoir Pressure
Temperature
Fundamentals of Reservoir Engineering & Characterization 53
In a black oil system the quality lines are broadly spaced at reservoir condition with
separator condition lying on relatively high quality lines. In atypical black oil system GOR
may decrease initially when the reservoir pressure falls below bubble point, as the
evolved gas remains immobile at very low saturation. The saturation pressure of black
oil is relatively low. Contribution of heavy compounds present in evolved gases in
reservoir to the total liquid recovery is not significant.
PVT Properties of Oil and Gas Knowledge of PVT is the first step in the study of any oil field as the information helps in
evaluating reserves, developing optimum recovery plan, and determining the quantity
and quality of produced fluids. In fact PVT parameters are required in every aspect of
reservoir engineering. Hence, accurate and reliable phase behaviour and volumetric
data are essential elements for proper management of petroleum reservoirs.
Most commonly information from black oil PVT tests are the oil formation volume factor,
Solution GOR, gas formation volume factor, as they are used to simplify engineering
calculations. Specifically, they allow for the introduction of surface volumes of gas, oil
and water into material-balance equation. Hence it would help if these terms are defined.
Oil Formation Volume Factor: It is defined as the number of reservoir barrels of oil and
dissolved gas that must be produced to obtain one stock tank barrel of stable oil at the
surface condition. Its unit is reservoir barrel/stock tank barrel.
Solution Gas Oil Ratio: It is defined as the number of standard cubic feet of gas
produced with each stock tank barrel of oil that was dissolved in the oil in the reservoir.
It’s unit is standard cubic feet/stack tank barrel.
Gas Formation Volume Factor:It is defined as volume in barrels that one standard
cubic foot of gas at the surface occupies as free gas in the reservoir. Its unit is reservoir
barrel/standard cubic feet.
Compositional studies are often conducted for gas condensate and volatile oil reservoirs,
where detailed informations on the fluid constituents are used to estimate fluid properties.
Only in special cases such as gas injection or miscible displacement, the compositional
approach is used for black oil reservoirs.
. Methods of Obtaining PVT data : There are primarily three methods by which PVT data are derived
Fundamentals of Reservoir Engineering & Characterization 54
1. Laboratory Measurements
2. Empirical PVT Correlations
3. EOS fluid Characterization
There are several PVT tests that are routinely conducted in the laboratory to study and
quantify the phase behaviour and properties of a reservoir fluid at simulated recovery
conditions.
Empirical correlations and charts, mainly reminiscence of days when hand calculations
were norm to predict PVT data, are still in vogue and much sought after.
A compositional phase behaviour model (EOS), can predict all the PVT data using only
the composition of the original reservoir fluid. However, the model first needs to be
evaluated and tuned against the measured PVT data prior to being used in reservoir
studies with confidence. With the advent of fast computers and robust algorithms
compositional data model are becoming very popular
Routine laboratory tests The majority of laboratory tests are depletion experiments, where the pressure of the
single phase test fluid is lowered in successive steps. The reduction of pressure results
in formation of a second phase, except in dry and wet gas mixtures. Hence type of
laboratory experiment to be conducted also depends upon the type of reservoir fluid.
Determination of fluid compositions is an important test on all reservoir fluid samples.
The gas and liquid phases are commonly analysed by gas chromatography and
distillation respectively.
Laboratory Tests for Dry Gas:In a dry gas reservoir system, no phase change occurs.
Hence, its composition remains the same. The only PVT test required for a dry gas is the
pressure-volume relation at the reservoir temperature and determination of specific
gravity, gas formation volume factor and isothermal compressibility.
Specific gravity = 96.28
Mg where Mg=molecular wt of gas
Bg = 0.02827 (Z T)/ P where Z is the compressibility factor, T and P is the reservoir temperature and pressure A typical gas formation volume factor graph is given below;
Fundamentals of Reservoir Engineering & Characterization 55
Isothermal Compressibility
TTg PZ
ZPPv
VC )(
11)(
1∂∂−=
∂∂−=
Laboratory Tests for Wet Gas:
PVT tests for a wet gas at reservoir conditions are similar to those for a dry gas.
Separate, tests are, however, needed to determine the amount and properties of
produced fluids at the surface conditions. The formation volume factor of a wet gas is
defined as the volume of the gas at reservoir conditions required to produce one unit
volume of the stock-tank liquid. The molecular weight and specific gravity of produced
condensate are also measured in the laboratory.
Laboratory Tests for Black Oil:
The phase transition of undersaturated oil during depletion can be best depicted as
given below. Let the reservoir pressure assumed to be higher than the bubble point
pressure. As the well is opened, pressure drops and as per theory of line solution most
of the pressure drop occurs near to the well bore. Away from the well bore, at Zone A,
which is far away from the well bore, the pressure is still higher than the bubble point
pressure. Hence oil expands as a single phase liquid.
Bg
Pressure
Bg Vs Pressure
Fundamentals of Reservoir Engineering & Characterization 56
The pressure at zone B, is just below the bubble point. Two phase region is formed.
However, the gas saturation is too small to allow its mobilization. The gas is assumed to
be in equilibrium with oil. This reservoir process is simulated in the laboratory at
reservoir temperature as constant composition flash vaporization. In the flash
vaporization the overall phase composition remains constant.
Zone C which is just ahead of the zone B, the gas bubbles coalesce together to form
bigger bubble, thereby, increasing gas saturation. Gas which was immobile in zone B,
now starts moving towards well bore in zone C. In the zone C the overall phase
composition doesn’t remain the same as gas moves out of the mixture. This reservoir
process is simulated in the laboratory at reservoir temperature as differential
vaporization. A series of flash tests at surface temperatures are also carried out to
simulate surface condition phase separation.
Constant Composition Expansion Tests: Constant composition expansion tests are carried out at reservoir temperature on gas
condensate or black oil to simulate the zone A and zone B as shown in above diagram.
The following important PVT properties are determined through these tests;
• Saturation Pressure (bubble point or dew point)
• Isothermal compressibility of the single phase fluid above bubble point.
• Compressibility factor of gas phase
• Total hydrocarbon volume as function of pressure.
Constant Composition expansion test can be schematically shown as following;
Fundamentals of Reservoir Engineering & Characterization 57
A typical PVT test data as reported by a laboratory for constant composition tests is
given below;
Pressure, psig Relative Volume Y function Density Reservoir pressure
Fundamentals of Reservoir Engineering & Characterization 59
• Relative volume is defined as volume of oil at indicated pressure per volume of residual oil at the standard condition
• Relative total volume is defined as volume of total oil plus liberated gas at indicated pressure per volume of residual oil at the standard condition.
• Bg is defined as the volume of gas at indicated pressure per volume at the standard conditions.
The differential liberation test is considered to better describe the separation process
taking place in the reservoir and is also considered to simulate the flowing behaviour of
hydrocarbon system at conditions above critical gas saturation
The test is started from bubble point pressure and the pressure is depleted till the
system pressure becomes atmospheric pressure. The test procedure can be
schematically shown as below;
Separator Test In the separator test, a known volume of the reservoir oil at its bubble point is flashed at
two or more stages, where the last stage represents the stock tank. For oils with high
gas in solution, more than one intermediate separator is often used. A field average
temperature is used for the separator tests. The experiment is carried out at number of
stages to determine the optimum field separation condition which gives lowest formation
volume factor and maximum stock tank oil. Operational limitation may, however, dictate
other pressure conditions in the field.
The behaviour of a reservoir oil during depletion is simulated by a combination of all
there types of tests discussed above. The reservoir oil remains single phase as long as
the pressure is above bubble point. The gas evolved just below the bubble point initially
Expelled Gas Expelled Gas
Oil
Gas
P>Pb P=Pb P<<Pb P<Pb
Fundamentals of Reservoir Engineering & Characterization 60
remains immobile in pores. These processes are simulated by constant composition
expansion test i.e. flash vaporization. The evolved gas begins to move away from the oil
as the gas saturation exceeds a critical value. The process then becomes similar to
differential liberation tests. A part of the gas, however, remains in contact with the oil in
well bore and their subsequent separation in the separator. A flash separation can only
simulate this process.
In material balance calculations and black oil simulation, the properties of fluid produced
at the surface are related to those at reservoir conditions by the results of separator tests,
and not those of differential liberation. The differential liberation test data are based on
the residual oil in the reservoir, whereas the volume factor and solution gas data are
based on the stock tank oil must be used in material balance equation and black oil
simulation. The corresponding values by differential test are almost always higher and
can lead to errors of 10 to 20% in the calculated oil in place and recoverable oil. Hence
following corrections are made in the Bo and Rs values.
Odb
ObFODO B
BBB =
Where BObF = Flash formation volume factor at bubble point BOdb = Differential formation volume factor at bubble point
−−=
Odb
ObFSidSifS B
BRsdRRR )(
Where Rsif= Solution GOR at Bubble point from flash Rsid = Solution GOR at bubble point from differential liberation Constant Volume Depletion Test: Constant volume depletion (CVD) experiments are performed on gas condensate and
volatile oil systems to simulate reservoir depletion performance and compositional
variation. It is commonly assumed that the condensate dropped out in pores remains
immobile. The depletion process is, therefore, simulated by CVD. The test consists of a
series of expansion followed by expelling the excess gas at constant pressure in such a
way that cell volume remains constant at the end of each stage. The expelled gas at
each pressure stage is collected and its composition, volume and compressibility factor
are determined. The condensate volume is also measured.
A schematic diagram of CVD test is given below
Fundamentals of Reservoir Engineering & Characterization 61
Laboratory Test for Volatile Oil
Pressure depletion in volatile oil is associated with high gas liberation. This gas phase
almost immediately becomes mobile. The differential test seems to simulate the process.
However, the mobile gas which is produced with oil behaves as a rich retrograde gas
and contributes significantly to the liquid production. None of the pressure depletion tests
commonly conducted in laboratories can simulate the fluid behaviour as occurs in the
field. Hence, PVT tests for volatile well are not well defined till now. However
compositional model after tuning with pressure volume data and amount of condensate
may mimic the phase behaviour to some extent.
Empirical Correlations
When a reservoir fluid study is unavailable, the engineer must rely on correlations to
estimate values of the physical properties of interest. The main properties which are
determined from empirical correlations are the bubble point, gas solubility, volume
factors, density, compressibility, and viscosity. The correlation typically matches the
employed experimental data with an average deviation of less than a few percent. It is
not unusual, however, to observe a much higher deviation when applied to other fluids.
There are many fluid property correlations. A number of these correlations have used
data of certain localities; hence, their application is limited. Some correlations have
received higher attention and acceptability than others. However no correlation has clear
superiority over other. Some of them have shown their reliability in various comparative
studies. Following table provides information on the range of data used in the
correlations to help selecting correlation for specific purpose.
Gas
Gas Gas
Gas
Condensate
P> Pdew P= Pdew P< Pdew P<< Pdew
Gas Gas
Gas
Fundamentals of Reservoir Engineering & Characterization 62
Correlation Standing Lasater Vasquez-Beggs
Glaso Marhoun
Pb, psia 130-7000 48-5780 15-6055 165-7142 130-3573 Temperature, OF
Before different correlations for the determination of physical properties are given, it is
advised to the reader to use only those empirical relations which satisfy the limitation
given in the above detail and also matches with the experimentally determined results.
Black Oil Physical Properties Correlation for Bubble point Pressure
Bubble point pressure Pb is defined as the highest pressure at which gas is first liberated
from the oil. The correlation to determine Pb are based on the fact that bubble point
pressure is a strong function of solution GOR RS, gas gravity gγ , Oil gravity OAPI, and
temperature T
Standing correlations Standing initially introduced a graphical correlation for determining the bubble point
pressure for Californian crude, and later expressed the graph by the following correlation
( )
−⊗
= 4.1102.1883.0
a
g
sb
RP γ
Where a = 0.00091T – 0.0125 (OAPI) Pb = Bubble Point Pressure, psia Rs = Solution GOR, SCF/STB T = Temperature, OF McCain suggested to replace gas gravity with separator gas i.e. excluding the gas from
stock tank would improve the accuracy of the equation.
Fundamentals of Reservoir Engineering & Characterization 63
Limitations:
• This correlation should be used with caution if non-hydrocarbon component are
also present.
• A deviation of about 15% is expected from this correlation.
Vasquez and Beggs Correlation Vasquez and Beggs pointed out that gas gravity depends upon separator pressure.
Hence, it used the gas gravity ( gnγ ) normalized to a separator pressure of 100 psig.
TS = separator temperature,OF PS = separator pressure, psia This method has an absolute error of 12.7%. Glaso’s Correlation Glaso developed the correlation from studying 45 North Sea crude oil samples
[ ]2** )log(30218.0)log(7447.17669.1)log( bbb PPP −+= Where *
bP is a correlating number and is defined by
cObagsb APITRP )()()/(* γ=
Glaso’s correction for correlating number to account for non-hydrocarbon component
and stock-tank-oil paraffinicity is not widely used. As large variation was observed in the
tuned equation of state model prediction with that of Glaso’s corrected and uncorrected
bubble points as shown in the Figure Below;
Fundamentals of Reservoir Engineering & Characterization 64
Paraffinicity is characterized by Watson’s characterization factor. Marhoun’s Correlation Marhoun used 160 experimentally determined bubble point pressure from the PVT
analysis of 69 Middle Eastern oils to develop a correlation for estimating Pb.
edO
cg
bsb TaRP γγ= psia
Where, T = temperature, OR Oγ = Stock tank oil Sp. Gravity
gγ = Gas specific gravity
a = 5.38088E-3 b = 0.715082 c = -1.87784 d = 3.1437 e = 1.32657
An average absolute relative error of 3.66% is observed. 5. The Petrosky-Farshad Correlation The gas solubility equation of Petrosky-Farshad can be solved for bubble point pressure
051.1391)10(
727.1128439.0
577421.0
−
=
Xg
Sb
RP
γpsia
Where X = (7.916E-4)(OAPI)1.541-(4.561E-5)(T-460)1.3911
T = temperature, OR 6. Lasater Correlation Lasater used mole fraction yg of solution gas in the reservoir oil as main correlating
parameter;
C7+ Watson Characterisation Factor
Pb, psia PR EOS
Glaso uncorrected
Glaso Corrected
Fundamentals of Reservoir Engineering & Characterization 65
g
b
TAP
γ= psia
Where T= temperature in OR 57246.017664.11083918.0 g
In summary significant variation will not be observed for most of the correlations for Pb.
However, Lasater and Standing correlations are recommended for general use and as a
starting point for developing reservoir-specific correlation.
Correlation for Solution GOR 1. Standing Correlation Standing correlation for solution GOR is given as follows;
2048.1
104.12.18
⊗
+= xgS
PR γ scf/STB
Where
X=0.0125 (OAPI)-0.00091(T-460) T=temperature, OR P = System Pressure, psia
gγ =solution gas specific gravity
This equation is valid for application at and below bubble point pressure. 2. Vasquez-Beggs correlation Vasquez-Beggs presented an improved empirical correlation for solution GOR using
5008 measured gas solubility data points. The correlation proposed is as follows
Fundamentals of Reservoir Engineering & Characterization 66
RS = Solution GOR, scf/STB The correlation is based on Californian crude sample. 2. Vasquez-Beggs Correlation Vasquez-Beggs developed a relationship for determining BO as a function of RS, Oγ , gγ
and T. The proposed correlation is based on 6000 measurements of BO at various pressures.
3. Glaso’s Correlation Glaso proposed the following relation based on North Sea crude oil BO = 1.0+10A bbl/STB
Where A= -6.58511+2.91329 log( *
obB ) – 0.27683 ( log( *obB ))2
*obB is a correlating number as is defined as *obB = [Rs( gγ / Oγ )0.526 ]+0.968(T-460)
Fundamentals of Reservoir Engineering & Characterization 68
T = temperature, OR
gγ = solution gas specific gravity
Oγ = Specific gravity of Stock tank oil 4. Marhoun’s Correlation Marhoun developed an equation for BO by the use of non-linear multiple regression analysis on 160 experimental data points and is given as below; B0 = 0.497069+0.862963E-3 T+0.182594E-2 F+0.318099E-5(F)2 bbl/STB Where c
Obg
aSRF γγ=
a = 0.742390 b = 0.323294 c = -1.202040 5. Material Balance Equation BO is defined as
O
gSOO
RB
ργγ 0136.04.62 +
=
where O =density of the oil at the specified pressure and temperature, lb/ft3
gγ = solution gas specific gravity
Oγ = Specific gravity of Stock tank oil RS = Solution GOR, scf/STB Error in calculating BO using material balance equation will depend upon the accuracy of
input variables ( gγ , Oγ and Rs) and the method of calculating O. All the correlations for
BO determination give approximately the same accuracy.
Sutton and Farshad’s comparative study of these correlation suggests that Standing
correlation is slightly better for Bob <1.4 and Glaso’s correlation is best for Bob >1.4.
The Standing and the Vasquez-Beggs correlations suggest that a plot of BO vs. RS
should correlate almost linearly. Hence this plot should be used for checking the
consistency of reported PVT data from a differential liberation plot.
Correlation for isothermal Compressibility Coefficient for Crude Oil Isothermal compressibility coefficients are required to solve many reservoir engineering
problems, including transient fluid flow problems; material balance equation and they are
Fundamentals of Reservoir Engineering & Characterization 69
also required in the determination of physical properties of undersaturated oil.
Compressibility is defined as;
T
O pv
VC
=
δδ1
Strictly speaking, the compressibility of an oil mixture is defined only for pressures
greater than the bubble point pressure. If oil is at bubble point pressure, the
compressibility can be determined and defined only for a positive change in pressure.
Implicit in the definition of compressibility is the fact that mass remains the same.
However, as the pressure is reduced below bubble point pressure gas comes out of oil
and as a result mass of the original system for which compressibility is to be determined
doesn’t remain the same.
Compressibility Factor for Saturated Oil Perrine introduced a definition for the compressibility of a saturated oil that include the
shrinkage effect of saturated-oil, p
BO
δδ
, and the expansion effect of gas coming out of
solution, T
sg p
RB
δδ
T
S
O
g
T
O
OO p
RB
B
PB
BC
+
−=δ
δδ
δ615.511
Compressibility Factor for Under Saturated Oil 1. The Vasquez-Beggs Correlation After studying a total of 4036 experimental data point for compressibility. Vasquez-
Beggs gave following equation for undersaturated compressibility of oil
Co = A/p psi-1
Where A=10-55Rsb + 17.2(T-460)-1180gs+12.61 OAPI
( )( )
+= −
7.114log10912.51 5 S
So
ggs
PTAPIγγ
T = Temperature, OR P = Pressure above bubble point, psia Rsb = Slolution GOR at the bubble point pressure TS = Separator temperature, OR PS = Separator pressure, psia gγ = Solution gas specific gravity
Fundamentals of Reservoir Engineering & Characterization 70
2. Standing Correlation Standing gave a graphical correlation for undersaturated Co that can be represented by
−−−−−+
= −
938.12))(4141.7(1.79)(004347.0
exp10 6
b
bobO PPE
PPC
ρ psi-1
Where ob = oil density, lbm/ft3
Any of the above correlations can be used for C0 determination. However, caution
should be exercised while calculating C0 for volatile oil where Co > 20X10-6 . In case of
volatile oil simple polynomial fit of the relative volume data, Vro = Vo/Vob from a PVT
report should be used for an accurate Co rather than using correlations.
Polynomial fit should be done as follows;
Vro = Ao+A1P+A2P2
2
21
21 )2(pApAA
pAAC
OO ++
+−=
Correlation for Oil Viscosity The live oil viscosity depends upon the solution gas content. Oil Viscosity decreases with
rising pressure as the solution gas increases, upto the bubble point pressure. There are
few empirical correlations to determine the viscosity of saturated and undersaturated
crude oil which accounts for the effect of dissolved gas and pressure on the viscosity of
dead oil.
Correlation for determining dead Oil Viscosity (od) For empirical correlation, the dead oil viscosity is determined first. The dead oil is
defined at atmospheric pressure and at any fixed system temperature without dissolved
gas.
1. Beal Correlation Beal presented a graphical correlation to determine dead oil viscosity, if the oAPI gravity
of the crude oil and the temperature are known. Standing presented this in the form of
mathematical equation for determining dead oil viscosity, od, at 14.7 psia and
temperature T, in oR;
a
APIod T
E
−
++=260
36078.132.0
53.4γµ cp
Fundamentals of Reservoir Engineering & Characterization 71
Where
+
= APIa µ33.8
43.0
10 2. Beggs and Robinson Correlation Beggs and Robinson presented an empirical relation for dead oil viscosity based on 460
dead-oil viscosity measurement
110 −= xodµ
where
163.1
)02023.00324.3(10T
xAPIγ−
=
Temperature T is in oF 3. Glaso Correlation Glaso developed empirical relation for dead oil viscosity base on crude oil samples of
North Sea.
( ) aAPIod TE )(log10141.3 444.3 γµ −+=
Where A = 10.313 log(T)-36.447 Temperature T is in oF 4. Kartoatmodjo and Schmidt Correlation In its empirical form this correlation is a combination of all three previous ones and can
be expressed as; )9718.26)log(7526.5(8177.2 )(log()8160( −−+= T
APIod TE µµ Temperature T is in oF Dead-oil viscosity is one of the most unreliable properties to predict with correlations
primarily because of the large effect that oil type (paraffinicity, aromaticity, and
asphaltene content) has on viscosity.
Correlation for determining live saturated Oil Viscosity The original approach by Chew and Connaly for correlating saturated oil viscosity in
terms of dead oil viscosity and solution GOR is still widely used. Most of other
correlations have in fact used the same concept in development of the relationship;
Chew and Connally gave the empirical relation as follows
( ) 21
Aodo A µµ = cp
Fundamentals of Reservoir Engineering & Characterization 72
This correlation is valid for GOR less than 1000 scf/STB. The functional relations for A1
and A2 reported by various authors differ somewhat, but most are best fit equations of
Chew and Connally’s tabulated results
1. Beggs and Robinson Correlation
A1=10.715(RS+100)-0.515
A2= 5.44(RS+150)-0.338
2. Bergman log(A1)= 4.768-0.8359 log(RS+300)
A2= 0.555 + 3005.133
+sR
3. Standing
2)72.2()44.7(
1 10 ss REREA −+−−=
SSS RERERE
A)374.3()51.1()562.8(2 10
062.010
25.010
68.0−−− ++=
4. Al-Khfaji et. Al. This correlation extended the Chew-Connally correlation to high GOR (upto 2000
scf/STB) 40
30
201 0631.04065.05657.02824.0247.0 AAAAA O +−++=
40
30
202 01008.00736.007667.00546.0894.0 AAAAA O +−++=
Where AO= log10(Rs) RS = Solution GOR, scf/STB Correlation for determining live undersaturated Oil Viscosity Oil viscosity at pressure above the bubble point is estimated by first calculating the oil
viscosity at its bubble point and then adjusting the bubble point viscosity to higher
pressure.
1. Vasquez-Beggs Correlation Vasquez-Beggs proposed follwing correlation for determining live viscosity above bubble
point pressure by analyzing 3593 data points
m
bob P
P
= µµ 0
Where
Fundamentals of Reservoir Engineering & Characterization 73
critical temperature adjustment factor. Following three empirical equations have been widely used in the oil industry for the
determination of Z factor.
• Hall – Yarborough method
• Dranchuk-Abu-Kassem
• Dranchuk-Purvis-Robinson
These empirical equations require iterative processes to get the z factor. Hence, it is
best used through computer programming. In this section only solution of Hall-Yarborugh
method would be explained.
Hall Yarborough method proposed the following mathematical equation to calculate z
factor;
))1(2.1exp(06125.0 2t
Y
tPZ pr −−
=
Where Ppr = Pseudo reduced pressure T= reciprocal of pseudo reduced temperature, i.e. Tpc/T And Y is the reduced density which is obtained as the solution of following equation;
Fundamentals of Reservoir Engineering & Characterization 75
0)3()2()1(
1)( 423
432
=+−−
++++= XYXYXY
YYYYXYF
Where ])1(2.1exp[06125.01 2ttPX pr −−−=
)58.476.976.14(2 32 tttX +−=
)4.422.2427.90(3 32 tttX +−= X4=(2.18+2.82t)
Procedure followed to solve the equation; Step 1. An appropriate initial guess for Yn is made. Where n is an iteration counter.
An appropriate initial guess is given as
[ ]2)1(2.1exp0125.0 ttPY prn −−=
Step 2: Initial value of Y is substituted in the function F(Y). Unless the initial value is the
correct solution the function will have non zero value.
Step 3: A new improved estimate of Y i.e. Yn+1 is calculated from the following
expression
)()((
'1
K
Knn
YFYF
YY −=+
Step 4; The procedure is repeated several times till absolute value of (Yn-Yn+1) becomes
smaller than 10-12
Step 5: The correct value of Y is then used to evaluate z Correlation for Gas viscosity 1. Lee et.al. Lee et. al. presented a semi empirical equation to calculate gas viscosity. The equation
can not be used for sour gas is given below;
= −
Yg
g XK4.62
exp10 4 ρµ cp
Where )19209/()02.04.9( 5.1 TMTMK gg +++=
gMTX 01.0)/986(5.3 ++=
XY 2.04.2 −=
Fundamentals of Reservoir Engineering & Characterization 76
g = Gas density at reservoir pressure and temperature, lbm/ft3
T = Reservoir Temperature, oR Mg = Apparent molecular weight of the gas mixture 2. Carr et. al. This correlation requires the knowledge of the gas composition and the viscosity of each
component at atmospheric pressure and reservoir temperature
=
==n
jjj
n
JJjj
ga
My
My
1
1
µµ
Where n is the number of component in the gas yj=mole fraction of component j
j = viscosity of component j Mj=Molecular weight of component j EOS Fluid Characterization
An equation of state is an algebraic expression that can represent the phase behaviour
of a multi-component mixture both in the two phase envelope and outside the phase
envelop i.e. outside the bimodal curve. The same EOS can be used to calculate the
properties of all the phase. Phase equilibriums are calculated with an EOS by satisfying
the condition of chemical equilibrium. For a two phase system, the chemical potential of
each component in the liquid phase must be equal to the chemical potential each
component in the liquid phase. A component material balance is also required to solve
vapour/liquid equilibrium problems. Solving phase equilibrium with an EOS is a trial and
error procedure, requiring considerable computation. With the advent of powerful
convergence techniques, nevertheless, solution of EOS has become robust and faster.
Before we deal with the subject on EOS, knowledge of few concepts like equilibrium
constant flash calculations etc. are must
Equilibrium Constant For a multicomponent system, such as petroleum fluids, the composition, pressure, and
temperature uniquely define the system phase behaviour. The equilibrium constant Ki, of
Fundamentals of Reservoir Engineering & Characterization 77
a component i is defined as the ratio of the mole fraction of the component in the gas
phase, yi, to the mole fraction of the same component in the liquid phase, xi
Ki = yi/xi
For a real solution, the equilibrium constant are not only function of pressure and
temperature but also the composition of the hydrocarbon phases. In compositional
modeling, the engineering objective is to determine the physical properties of the
individual gas and liquid phases. Consequently, the equilibrium constants which indicate
partitioning of each component between the liquid and gas phases must be known.
Hence it is appropriate to introduce flash calculation, which is the workhorse of most
EOS application
Two-Phase Flash Calculation The two-phase calculation consists of defining the amounts and composition of
equilibrium phases, given the pressure and temperature, and overall composition. An
inherent obstacle to solving the problem, is not knowing, whether mixture may exist as
single phases or split into two or more phases.
The two phase split calculation (Rachford-Rice procedure) can result in either a solution
yielding equilibrium phase composition or a trivial solution. Even when the results appear
physically consistent, a rigorous check of the solution with the phase stability test may
be required. Mathematically the two phase flash calculation is solved by satisfying the
equal fugacity and material balance constraint with a successive substitution or Newton
Raphson algorithm.
The component and phase material balance constraints state that n total moles of feed
with composition zi distribute into nv moles of vapour with composition yi and nL moles of
liquidwith composition xi
The material balance constraint can be written as
n = nv+nL nzi = nvyi +nLxi
Let FV = vapour mole fraction = )( VL
v
nnn+
Hence Zi = FVyi + (1-FV) xi Additionally , the mole fraction of equilibrium phases and the overall mixture sum to unity
Fundamentals of Reservoir Engineering & Characterization 78
===
==n
ii
n
ii
n
ii Zxy
111
This constraint can be expressed as
0)(1
=−=
n
iii xy
Since ki = yi/xi
H(Fv) = =−=
n
iii xy
1
)( 0)1(1
)1(
1
=−+
−
=
n
i iv
ii
kFkz
The above equation is referred to Rachford-Rice equation. With feed composition and k
values known, the only remaining unknown is Fv . H(FV) has asymptotes at Fv =1/(1-ki).
This can be shown in the following graphs;
Muskat and McDowell proposed a solution to the phase split calculation by assuming Ci
= 1/(ki-1), wher Ci = ∞ for ki=1. They proposed following form of the function H(Fv)
H(Fv) = =+
0iv
i
CFZ
Where +−=
2)( iv
i
v CFz
Fh
δδ
Using modified regula falsi method solution converges for FV.
Having solved for Fv, phase compositions in different phases are calculated as
1)1( +−
=iv
ii KF
zx
FvMin FV Max
H(Fv)
Fundamentals of Reservoir Engineering & Characterization 79
1)1( +−
=iv
iii KF
kzy
The composition calculation requires value of ki at the pressure, temperature, and
composition of each phase. There are several methods to determine phase equilibrium
constants, including use of charts. Recommended way, however, is to determine
through equation of state by rigorously checking for stability by minimum fugacity energy
constraint for individual component.
The need for EOS rose when it was established that the equality of fugacity of each
component throughout all phases to be the requirement for chemical equilibrium in multi
component systems. , The fugacity coefficient iφ is defined as;
∞
−
−
=
≠V nVTii ZdVVRT
nP
RTj
ln/1
ln1,,
δδφ
The fugacity coefficient can, therefore be determined from the above with an aid of
equation of state (EOS). Equation of states is basically developed for pure components.
However by employing some mixing rules to determine their parameters for mixtures, it
can be used for multi component mixtures. The mixing rules are considered to describe
the prevailing forces between molecules of different substances forming the mixtures. It
is the capability of EOS and the associated mixing rules determines the success of
phase equilibrium prediction.
Before deliberating on procedure for determining phase behaviour at different pressure
through equation of state, it would be pertinent here to deliberate on types of EOS.
Equation of States Vander Waal’s first proposed the following equation of state by considering the
intermolecular attractive and repulsive forces;
( ) RTbvVa
P =−
+ 2
Where, a/v2 and b represent the attractive and repulsive terms respectively and v is the
molar volume. As the pressure approached infinity, the molar volume becomes equal to
Fundamentals of Reservoir Engineering & Characterization 80
b. Hence, b is also considered as an apparent volume of the molecules and is less than
molar volume v.
The above equation in terms of compressibility factor takes the cubic equation form; 0)1( 23 =−++− ABAZZBZ Where The dimensionless parameter A and B are defined as
2)(RTaP
A ≡
RTbP
B ≡
Hence Vander Waals type of EOS is referred as cubic EOS. A typical response of Van
der Waals EOS is shown below;
Based on the response equation of state can be divided into two main group: cubic and
non cubic. Cubic equation have three roots when T≤ Tc (Critical temperature) and only
one root when T> Tc. At T=Tc, there are three equal roots.
Following figure depicts the deficiency which most of the cubic equation of state exhibit.
Volume
Tc
T1
T3
Pressure
Fundamentals of Reservoir Engineering & Characterization 81
As can be seen the EOS are poor in predicting the under saturated liquid density.
Whereas, they can predict the gas phase volume and density remarkably well. A number
of EOS has been proposed by different authors. Notable among them are, Peng
Robinson (PR) EOS, Zudkevitch-Joffe modification of Redlih-Kwong (ZJRK) EOS,
Soave-Redlich-Kwong (SRK) EOS etc.. However, none of them can be singled out as
the most superior equation to best predict all properties at all conditions. A number of
comparative studies have, however, showed that certain equations exhibit a higher
overall accuracy.
Peng Robinson (PR) EOS and Soave-Redlich-Kwong (SRK) EOS take the general form
of ;
22 )1( cbcvbva
bvRT
p−++
−−
=
When c=1, the above equation becomes the Peng Robinson (PR) EOS and when c=0, it
becomes the Soave-Redlich-Kwong (SRK) EOS.
22 2 bvbva
bvRT
p−+
−−
= - PR EOS
vbv
abv
RTp
+−
−= 2 - SRK EOS
C
Volume
Tc
T1
T3
Pressure EOS Observed
Fundamentals of Reservoir Engineering & Characterization 82
For pure components, the parameter a and b are expressed in terms of the critical
properties and accentric factor ();
αcaa =
( ) ccac pRTa /Ω=
CT
Tk −+= 1(1α
ccb pRTb /Ω= Let A = ap/(RT)2 B= bp/RT By putting compressibility factor Z= pv/RT, THE general EOS becomes [ ] ( )[ ] 0)21()()1( 23223 =+−−+−+−+−− BBcABcBciBAZcBZZ For multi component system parameters a and b are defined using the following mixing
rule;
−== j
jijji
ii adxaxa )1(1
=i
iibxb
where dij is an empirically determined interaction coefficient. Fugacity coefficient is given by
∞
−
−
=
≠V nVTii ZdVVRT
nP
RTj
ln/1
ln1,,
δδφ
Solving and using EOS results in following general expression for fugacity
++
−−
−−−−−=
BZBZ
bb
a
adxa
BA
BZZbb i
jijjii
i1
2
12
ln)1(2
1)ln()1(ln
δδ
δδφ
T The derivation is complicated and beyond the scope of the training programme. Parameters ba ΩΩ , are determined from the critical condition.
At the critical point, the compressibility factor will have three real and equal roots.
(Z-Zc)3 = 0 EOS c δδδδ1 δδδδ2 ΩΩΩΩa ΩΩΩΩb Zc
Peng Robinson 1 1- 2 1+ 2 0.45724 0.07780 0.307
SRK 0 0 1 0.42747 0.08664 0.333
Fundamentals of Reservoir Engineering & Characterization 83
The cubic EOS yields three real roots in the two phase region. The one having lowest z
factor is taken for liquid fraction and highest is taken for vapour phase. The in between z
factor is discarded.
To make the EOS more reliable in predicting liquid phase volume and densities,
Peneloux et.al proposed volume translation technique. A volume translation technique
modifies the molar volume of the system v predicted by the equation of state as follows;
Vcor = v-c
= ii rxc
cii Btr = ti is the dimensionless individual translation value for each component.
c
cbc P
RTB Ω=
Volume translation is found to have no effect on the equilibrium conditions. Therefore, it
doesn’t alter saturation pressure, saturation temperature, equilibrium composition etc.
However, it will modify the molar volumes, compressibility factors and densities of the
fluid.
Solution Algorithm for Phase Split Calculation through EOS To know the mole fraction of a component “i” in the liquid and vapour phase,
compressibility factor for liquid and vapour at that pressure and temperature is must. To
determine the compressibility factor, Z, in the liquid or gas phase, the appropriate EOS
can be solved either by direct or iterative methods. These equations are cubic equations
that yield a single root in the single phase region and three real roots in the two phase
region. The largest root of the cubic equation corresponds to the vapour phase and the
smallest root corresponds to that liquid phase.
The following is the step by step procedure to calculate equilibrium constant and hence
split mole fraction in vapour and liquid phase;
1. The input data required for this calculation are the system pressure, p, temperature,
T and over all system composition.
2. The flash calculation is initialized by estimating a set of k values for each component
in the mixture.
Fundamentals of Reservoir Engineering & Characterization 84
( )( )[ ]
( )PciP
TTwk ciiOld
i
)/(1137.5exp −+= ---- Wilson Equation
3. With estimated ki value Rachford-Rice equation is solved for Fv, with the search for
Fv always lying between Fvmin and Fvmax.
Fvmin = 01
1
max
<− K
and FVmax = 11
1
min
>− K
4. With the determined value of Fv composition of each component in liquid and vapour
phase is determined. Now the next step would be know whether the determined
phase composition is stable or not.
5. Having calculate xi and yi , cubic equation is solved . [ ] ( )[ ] 0)21()()1( 23223 =+−−+−+−+−− BBcABcBciBAZcBZZ 6. Out of the three roots, the middle z value is discarded. Lowest z value is designated
as liquid phase compressibility factor ZL and highest Z value is designated as vapour
phase compressibility factor, ZV . With the help of ZL and ZV liquid phase fugacity and
vapour phase fugacity is determined.
++
−−
−−−−−=
BZBZ
bb
a
adxa
BA
BZZbb
L
Lijijji
LLi
Li1
2
12
ln)1(2
1)ln()1(ln
δδ
δδφ
++
−−
−−−−−=
BZBZ
bb
a
adxa
BA
BZZbb
V
Vijijji
VVi
Vi1
2
12
ln)1(2
1)ln()1(ln
δδ
δδφ
7. Using fugacity value, new equilibrium constant, k, is determined
Vi
Li
i
iNewi x
yk
φφ
==
8. Following convergence criteria is tested
=
−≤−n
iOldi
Newi
kk
1
122 10)1(
9. If the condition at step 7 is not satisfied. With the new k value the procedure from
step 3 is repeated. Till the condition is met
10. Following three types of converged solution we can get
Fundamentals of Reservoir Engineering & Characterization 85
a. A physically acceptable solution is found with o ≤ FV ≤ 1. When Fv=0 , it
correspondence to bubble point condition, when Fv=1 it correspondence to
dew point condition when o ≤ FV ≤ 1, it corresponds two phase condition.
b. A physically unacceptable solution is found, when FV <0 or Fv>1, when the
calculated equilibrium constant satisfy the equal fugacity condition and the
mathematical material balance equation.
c. A so called trivial solution when k value equals one i.e. xi = yi =zi
The solution “a” is usually correct solution. However, two phase condition stability should
be further analysed with the Gibbs tangent plane criterion for minima of Gibbs energy.
Gibbs tangent plane criterion is very complex and computer intensive and beyond the
scope of this training programme. The trivial solution at c should be checked with phase
stability test to find whether the mixture is in single phase.
Fundamentals of Reservoir Engineering & Characterization 86
Introduction The aim of well testing is to get information about a well and about a reservoir. Once the presence of hydrocarbon-bearing-formation is established and obtained its
porosity, thickness and hydrocarbon saturation. Well test analysis helps to get the
answer of three most important questions;
a. What is the volume of hydrocarbon in the reservoir system?
b. At what rate the available hydrocarbon fluid should be produced at the surface?
c. How much of the fluid can be recovered?
Besides it also provides information about following reservoir parameters;
1. Interwell flow capacity of a reservoir
2. static well pressure
3. Extent of well damage
4. Distance to nearest boundary
5. Detecting heterogeneity with in the pay zone.
Answer to the questions and information of reservoir parameters will establish the
commercial viability of the prospect and is the task of reservoir engineer.
Well test analysis is a branch of reservoir engineering. It uses the pressures and rates
under a standard condition for the determination of parameters which influences the fluid
flow through porous media e.g. permeability, fault, fluid contacts etc. Measuring the
variation in pressure versus time and interpreting them give data on the reservoir and
well.
There are few special pressure transient tests, which, can be used to determine the
areal extent of a reservoir and to estimate the volumes of fluid in place. In case of
composite systems like in-situ combustion, steam flooding or polymer flooding, these
Fundamentals of Reservoir Engineering & Characterization 87
well tests can accurately predict the swept zone parameters, enabling the engineers to
determine the efficacy of the processes. Pressure measurements can also be
interpreted to yield quantitative estimation of the well condition, so the efficacy of
stimulation treatments on well productivity can be evaluated.
Well tests do not directly provide estimates of permeability, well condition, pore volume
etc. Measurement must be analyzed and interpreted using a number of laws of fluid
mechanics to arrive at the desired result.
Diffusivity Equation The fundamental basis of transient flow theory is the diffusivity equation, a differential
equation that must be satisfied when fluid flows through porous body under isothermal
conditions. When Darcy’s law is applied to continuity equation, an equation, which, is
developed from conservation of mass principle, it gets transformed to diffusivity equation
that governs the pressure distribution for flow through porous media. The derivation of
the equation is based on two laws and one equation of state which are;
• Darcy law
pgradk
Vµ
−=
It is assumed that Darcy’s law governs the fluid flow. Darcy’s law is not applicable
macroscopically throughout the flow period. It is applicable microscopically during the
time interval when the various parameters and flow rate can be considered constant.
The gravitational forces are neglected.
• Material Balance
It is assumed that mass of fluid contained in the reservoir volume unit is equal to the
difference between the amount of fluid input and the output during the time interval.
0)( =+∇ OSt
V ρφδδρ
• Equation of State
The gravity of the fluid varies with pressure and the variation is equivalent to the
compressibility of the fluid
Fundamentals of Reservoir Engineering & Characterization 88
Te pC )(
1∂∂= ρ
ρ
The following section will derive first the continuity equation based on material balance.
Before that an attempt is made in the following section to derive a mathematical
formulation of what actually happens in a reservoir when a well is flowed, following
simple model and assumptions are needed. It is assumed that;
• A vertical well of radius rw intercepts a horizontal formation of constant thickness
h and of infinite extent.
• The formation is having uniform porosity φ and isotropic permeability K.
• Constant viscosity .
• Constant total compressibility Ct
• The rock properties are not time dependent.
Under these conditions the flow is radial.
Let us assume two dimensional flow in the x-y plane and consider a control volume of
infinitesimal dimensions shown in fig. below. Let us assume that the dimensions of the
control volume are x and y with unit depth perpendicular to the plane of the paper. Let
us also assume that gravitational effect is negligible here.
Let G be the mass flux and Gx and Gy be the component of G in X and Y direction,
respectively.
The various equations reflecting the conservation of mass principle for the volume
element shown in fig above are;
Outflow= (Gy + xy
yGy δδ )∂
∂ + (Gx + yx
xGx δδ )∂
∂
yxt
δδφρ∂
∂ )(
(Gy + xyy
Gy δδ )∂
∂
(Gx + yxx
Gx δδ )∂
∂
xδ
yδ
yδ
xδ
y
x
Fundamentals of Reservoir Engineering & Characterization 89
Inflow = Increase of Storage = Application of the conservation of mass principle will yield
0)( =
∂∂+
∂∂+
∂∂
tyGy
xGx φρ
For three dimension space this equation can be written as
0)( =
∂∂+
∂∂+
∂∂+
∂∂
tzGz
yGy
xGx φρ
For steady flow the above equation can be written in vector notation as
0=∇G For constant density and no pore volume change in time (G = .v) the above equation
can be written as;
0=∇v Where v is the Darcy flow, normal to unit cross section area of the flow. The above equation is called as continuity equation. The appropriate differential equation is obtained by combining the continuity equation,
the flux (Darcy) Law, and an equation of state. Ignoring gravitational effects Darcy law is
given by
v = - pk ∇µ
Substituting this in the continuity equation we have
typk
yxpk
x ∂∂=
∂∂
∂∂+
∂∂
∂∂ )(φρ
µρ
µρ ---- eqn. (A)
Assuming constant compressibility we can write
xp
cxp
∂∂=
∂∂
)(1ρ
Substituting this in eqn. (A) we get
tp
cyp
xpk
cypk
yxpk
x ∂∂=
∂∂+
∂∂+
∂∂
∂∂+
∂∂
∂∂ )(
)(22
φµµµ
yδ xδ
yxt
δδφρ∂
∂ )(
Fundamentals of Reservoir Engineering & Characterization 90
The above equation is a non linear partial differential equation. If we ignore the second
degree term because of its very-very small value then becomes
tp
cypk
yxpk
x ∂∂=
∂∂
∂∂+
∂∂
∂∂ )(
)(φµµ
This equation is a linear equation provided that is a constant. If K and are constant,
then classical methods of solution can be used to obtain the pressure distribution. The
above equation, under these assumptions can be written as
tp
kc
yp
xp
∂∂=
∂∂+
∂∂
)(2
2
2
2 φµ
The above equation is popularly known as diffusivity equation and defines the movement
of fluid into, out of and through the rock pore spaces. The expression suggests that the
pressure disturbance or perturbation diffuses rather than propagates. Had the
perturbation effect propagated in the reservoir, the expression would have been the
second order differential equation versus time.
In a cylindrical coordinate system, the diffusivity equation is represented as;
tp
kcp
rrp
rrr ∂
∂=∂∂+
∂∂
∂∂ φµ
θ 2
2
2
1)(
1
In a radial system, to which most of practical field solutions are arrived at is given as
below;
tp
kc
rp
rrp t
∂∂=
∂∂+
∂∂ φµ
0002637.011
2
2
in field unit
The term tc
kφµ
0002637.0 is called the diffusivity of the medium. It is a measure of how fast
the pressure perturbation will diffuse in the reservoir.
The same set of differential equation arises in many other contexts, and is not unusual to
obtain solution for flow through porous media by mere change in notation. The notable
other contexts are diffusion, diffusion and chemical reactions, and electrical problems
etc..
The solution of the above equation relies on the concepts of dimensionless pressure and
dimensionless time. The basic advantage of these groups is that they permit us to
Fundamentals of Reservoir Engineering & Characterization 91
understand the structure of the solution of interest without consideration of the specific
values of formation properties, fluid properties, or flow rates. The general solution of
diffusivity equation in dimensionless form is given as
)()(),( 0 DoDDD rsBIrsAKsrp += Where; s is the Laplace transform variable with respect to tD and I0(x) and Ko(x) is the modified
Bessel functions of the first and second kind of order 0 respectively. A and B are
constants.
The line source solution of the above equation in the dimensionless form can be written
as
−−=
D
DiD t
rEp
421 2
Using the following dimensionless factors;
)],([2.141
trppqB
khp iD −=
µ
w
D rr
r =
trc
kt
wtD ∆=
2
0002637.0φµ
Selecting a definition of dimensionless pressure is a difficult task. Basically the selection
depends on the wellbore condition-constant rate or constant well bore pressure.
Commonly used dimensionless groups follow from the seminal work of van Everdingen
and Hurst.
The general expression for pressure which is given as;
( ) )4
(4
,2
ktr
Ekh
qBtrpp ii
−−=−π
µ
The above equation suggests that the whole reservoir is affected due to perturbation
created by flow of fluid. However, the practicality of the above expression lies in the fact
that it used to locate the compressible zone w.r.t time within the reservoir. The pressure
drop in the well mainly reflects the reservoir properties in the compressible zone.
Fundamentals of Reservoir Engineering & Characterization 92
That is what a well test enables us to i.e.;
• characterize the average properties far away from the well
• detect facies heterogeneity
• identify permeability barriers.
• define composite system
Before we delve on further on well test analysis it would be appropriate to define various
regime and type of flows. Based on their boundary condition of the flow the diffusivity
equation has been solved and solution has been provided w.r.t time.
Transient State Condition Transient state which is also named as unsteady state as the condition at which rate of
change of pressure with respect to time at position with in the compressible zone due to
perturbation effect is not zero or constant.
Mathematically it can be expressed as
p = f (r,t)
),( trftp =
∂∂
Both are function of time and distance Pseudo Steady State or Semi Steady State Condition It’s the condition which follows transient state. The compressible zone has reached to
the reservoir boundary (no flow boundary), and due to no support outside the boundary
of the reservoir in the form of fluid influx, the pressure declines linearly w.r.t time.
Mathematically this is defined as
0=∂∂
rp
at r= re i.e. no influx
tConstp
tan=∂∂
for all r and t
Steady State Condition In a steady state condition the pressure at every location in the reservoir remains
constant i.e. it does not change with time
Mathematically this can be expressed as
0=∂∂
tp
for all r and t
p= Pe constant at r= re
Fundamentals of Reservoir Engineering & Characterization 93
Radius of Investigation Since the pressure variation in a well test represents the properties of the part of
reservoir involved in the compressible zone. It is important to locate the compressible
zone and this is what is involved in the concept of radius of investigation. However,
caution should be exercised in finding out radius of investigation due to the fact that
radius of investigation is actually a circular system with a pseudo-steady pressure
distribution hence is not meant for locating compressible zone with in transient state
where pressure p and rate of change of pressure w.r.t. time is a function of time and
distance.
The expression for radius of investigation is given as;
t
i ckt
rφµ
032.0= in field units
Well Bore Storage A well test begins with a change in production rates. Because the flow rate is usually
controlled at the well head, the compressible fluids in the wellbore do not allow an
immediate transmission of the pressure disturbance down to the surface. As a result of
compressibility of the fluid column, inequality in the surface and sandface flow rates
occurs resulting in accumulation of mass in the wellbore. The phenomenon is known as
wellbore storage. This can be graphically shown as below;
For a typical drawdown and buildup tests, the well bore storage phenomenon is known
as unloading and afterflow, respectively. In either case, during the initial stages of the
Well Head Flow Rate
Bottom Hole Flow rate qB
qt
Wellbore Storage Effect
Time
Fundamentals of Reservoir Engineering & Characterization 94
test, there is a variable rate at the sandface which invalidates the assumption of constant
production rate for which the solution of the diffusivity equation has been arrived.
Mathematical model
For a radial and constant flow rate in an infinitely large reservoir the flow model in
dimensionless pressure term can be given as;
D
D
D
DD
DD tp
rp
rrr ∂
∂=
∂∂
∂∂1
Near to the well bore
1][ 1 −=∂∂
=DrD
D
rp
However because of well bore storage effect the expression near to sand face becomes
DrD
D qsfrp
D−=
∂∂
=1][
During the initial stages of well testing, the ratio of sand face flow rate and well head flow
rate from material balance and compressibility factor can be given as;
DD
wDD
wh
sf qsfdt
dpC
q
q=−= 1
where
2
894.0
wt
sD rch
CC
φ=
This makes the dimensionless pressure expression near to the well bore;
11
−∂
∂=
∂∂
= D
wDD
rD
D
tp
Crp
D
In case of a drawdown test, the initially produced fluids are being unloaded from a well
bore with very or no flow at the sand face. That makes 0][ 1 =∂∂
=DrD
D
rp
a sense. i.e.
1=∂
∂
D
wDD t
pC
Integrating and taking logs both side of equality, we get DwDD tpC logloglog =+ in a dimensional form
Fundamentals of Reservoir Engineering & Characterization 95
C
qBtp
24=∆
Significance of the above relation is that, should the early data points (plotted in terms of
coordinates Dplog and Dtlog ) exhibit a unit slope line, then most fluid produced
originates from the well bore. As the test progresses, the sand face rate is significantly
increased the term log pwD increases. As the storage effect diminishes log pwD is
described by the flow equation for constant production rate for a drawdown. For a build-
up similar methodology and expression exists.
Hence, it can be deduced that the deviation from unit slope line marks the end of well
bore storage effect. This is a very important observation for a reservoir engineer, as it
would be possible to devote more on quality data representing solution of flow model.
In a linear plot of del(p) vs. time, slope of the straight line is used to compute well bore
storage effect
i.e.
slope
qBc
⊗=
24
The straight line should pass through the origin. There can be several reasons for it not
to pass through the origin;
• a shut in pressure error
• a shut in time error
t
p∆
End of Well Bore Storage
Fundamentals of Reservoir Engineering & Characterization 96
If such error happens, should be corrected by offsetting the data. However, caution
should be exercised in case of following cases;
• Time duration between two consecutive measurements is very high.
• Variable well bore storage due to gas
• Fluid segregation in the well.
Skin Factor Originally the skin factor response was introduced to incorporate the noted difference
between measured pressure response and predicted pressure response by diffusivity
equation model. The measured pressure responses were usually lower than the
predicted pressure. Van-Everdingen and Hurst suggested that the extra pressure drop
reflects a small region of low permeability (damage) around the well bore. In fact they
are credited for introducing the term skin factor to the oil industry.
Skin factor makes the vicinity of well bore characteristics different from those in the
reservoir as a result of drilling and well treatment operation. It reflects the connection
between the reservoir and the well. The difference in pressure drop in the vicinity of the
well bore can be interpreted in several ways;
• By using infinitesimal skin and is defined by s. If sp∆ is the pressure drop due
to skin
+−
−= 23.3log151.1 2
wt
wfi
rck
m
ppS
φµ in field units
µqB
pkhS S
2.141∆
=
• By finite Thickness Skin
KS
K
rs
rw
Ks
Fundamentals of Reservoir Engineering & Characterization 97
W
S
S rr
kk
S ln1
−=
Above equation shows that damage corresponds to a positive skin, and, improvement in
the flow due to well bore treatment corresponds to negative skin. It should be noted that
the negative value of skin may go upto maximum -5. Reporting of skin lower than this
value would have to be doubly sure by reservoir engineers from model verification.
Secondly, +ve skin are reported in the literature as high as +500. However, it is
cautioned that any value greater than +10 should be doubly verified from model
response.
• Effective Radius
The effective radius method consists in replacing the real well with a radius rw and skin
by S by a fictitious well with radius rws
)exp( Srr wws −=
In case of a gravel pack, the effective radius of the well should normally fall between the
screen radius and the under reaming radius. An effective radius that is less than the liner
radius would mean that the gravel pack is particularly ineffective.
As, could be understood from previous section that the skin represented an additional
pressure drop located in the vicinity of the well bore. The additional pressure drop was
explained for understanding due to variation in permeability in an area near to the well
bore. However, the skin concept could be generalized in more practical aspects of well
bore flow phenomena. For example;
P real
P with effective radius
rS rW
Ks < k
Fundamentals of Reservoir Engineering & Characterization 98
• Skin can be used in representing pressure drop due to partial perforation
• Inclination of a well improves the flow in the vicinity of the well bore, which can be
represented as negative skin.
• It can be used as –ve skin to represent the flow characteristic improvement in
hydraulically fractured well.
• Injection of fluid (water or Polymer) into the reservoir creates a composite zone of
different mobility ratio. It causes additional pressure drop that can also
considered as skin.
• In gas well, Darcy law breaks down at high velocity of flow. At high flow velocities,
pressure drop in porous media increases more than predicted linearly with
increasing rate. This extra pressure drop is accounted for by rate dependent skin.
Interpretation Methods Well test can be effectively used to know reservoir flow characteristics, as the pressure
variation near a well bore reflects the reservoir properties in the compressible zone.
Hence well tests are used by reservoir engineers to know a number of reservoir flow
parameters which decides the exploration, development and exploitation of a reservoir.
A number of different methods are used to analyze a well test. However, they can be
classified broadly into two heads;
• Methods using the Type curves
• Conventional methods.
Inside each of the groups the above methods depend on the type of well, reservoir and
reservoir boundaries
Type Curves Type curves are basically a graphical representation of the theoretical response during a
test of an interpretation model that represents the well and the reservoir being tested.
Most of the type curves are for drawdown well test response.
These type curves first appeared in seventies in the form of sets of curves using
dimensionless parameters. From 1983 on, type curve methods were greatly improved as
Fundamentals of Reservoir Engineering & Characterization 99
they were used in conjunction with the pressure derivative. With the advent of powerful
computers and programming use of Type curve has become very easy to use and
interpret. We would limit the introduction of Type curves only for vertical wells completed
in an infinite reservoir.
There are several kinds of Type curves commercially available. Few notable among
them are;
• Agrawal et. al. Type curve
• Mckinley Type curve
• Earlougher and Kersch Type curves
• Gringarten et. al. Type curve.
However, the most widely used Type curves in the oil industry is that of Gringarten et. al.
curves as they are most complete and practical to use. Hence a brief introduction about
Gringarten et. al. curves and methodology for use is given;
Gringarten et. al. Curves Gringarten et. al. Type curve represents the variation in pressure versus time for a
specific reservoir-well configuration. It is calculated using an analytical model and
expressed in dimensionless variables. In the form of
)],([2.141
trppqB
khp iD −=
µ
w
D rr
r =
trc
kt
wtD ∆= 2
0002637.0φµ
2
894.0
wt
sD rch
CC
φ=
Fundamentals of Reservoir Engineering & Characterization 100
In a vertical well in an infinite homogeneous reservoir the dimensionless pressure
variations depend on three factors: time, wellbore storage and skin. i.e.
( )SCtpp DDDD ,,= with Gringarten using the form below;
))2exp(,( SCCt
pp DD
DDD =
Interpretation Method The interpretation method using type curves involves the following steps;
• The pressure drop with respect to the initial pressure should be plotted on a
tracing paper lying on the type curve, using the same scale as that of Type curve.
Keeping the two coordinates axes parallel, the tracing paper is shifted to a
X axis
Fundamentals of Reservoir Engineering & Characterization 101
position on the type curve that represents the best fit of the measurement. Only
translational movement is allowed keeping the two grid parallel.
• To evaluate reservoir parameters, a match point is selected anywhere on the
overlapping position of the curves, and the coordinates of the common point on
both sheets are recorded. Once the match is obtained, the coordinates of the
match point are used to compute formation flow capacity, kh and storavity
constant tChφ . The specification of the type curve where the measured points
match, they correspond to a value of )2exp( SCD
( )( )M
MD
Pp
qBkh∆
= µ2.141
( )
MD
D
M
Ct
tkhC
∆=
µ000295.0
( )
D
D
CSC
S)2exp(
ln21=
What about Build-up Interpretation Using Type Curve? Types curves were established for constant flow rate production i.e. drawdown. Hence
Type curve analysis for build-up should be done with caution.
In following condition only, type curve should be used to interpret build-up data.
Condition 1: Build-up duration should be very-very small compared to production
duration of the well.
Condition2 : Build-up duration should be smaller than the duration of the last production
period before shut-in.
Other than the above condition, it would be incorrect to use for Build-up interpretation
without incorporating certain changes. The effect of short production time can be seen
in a flattening out of the type curve, the build-up curve under the drawdown curve. Force
match between the build-up data and a draw down curve would result in a type curve
located too high on the set of curves and therefore in inaccurate results.
Fundamentals of Reservoir Engineering & Characterization 102
The most useful method of using drawdown type curves for build-up is Agrawal’s method.
It consists of plotting each measurement versus an equivalent time et∆ as defined below
instead of t∆ .
p
e
tt
tt
∆+
∆=∆1
There is a condition to be satisfied before which this equivalent time can be used in
Gringarten Type curve or any Type curve. The condition is; the semi-log straight line
should have reached during the previous drawdown before build-up.
Advantages of Type Curves lie in the fact that it allows the interpreter to make a
diagnosis about the type of reservoir and understand the flow regimes. It also allows the
interpreter to use the flow concept regime in a conventional interpretation method with
ease and confidence. However, assumption of constant well bore storage effect in a
Type curve puts severe limitation to the interpretation.
Conventional Method of Well Test Interpretation This particular section will study the response of flow/pressure behaviour at constant
rate (drawdown) or when rate is zero (Build up).
Drawdown test The solution of diffusivity equation in the transient pressure regime is given as;
Type curve from Type curve set
Type curve calculated for a shut-in well
Fundamentals of Reservoir Engineering & Characterization 103
+−+=− S
rck
tkh
qBpp
wtwfi 87.023.3loglog
6.1622φµ
µ
Hence, if pressure measured at the bottom of well bore is plotted against log of time, it
would result in a straight line with slope, m
kh
qBm
µ6.162= and
+−
−= 23.3log151.1
21
wt
hri
rck
mpp
Sφµ
This slope m can be used for calculating flow capacity, kh of a reservoir and skin. P1hr is
the pressure at 1hr from the start of drawdown test, read from the straight line equation.
Transient state is of short duration. If the test is extended and the compressible zone is
allowed to travel and reach the boundary of the reservoir, the flow regime changes to
pseudo steady state regime in absence of any support from the outer boundary. The
solution of pressure response in a drawdown test in a pseudo-steady state is given as;
+++=− S
CrA
khqB
thAc
qBPp
Awtwfi 87.0
2458.2loglog
6.162234.02
µφ
The above equation suggests that a plot of pressure against time in the pseudo-steady
state region would result in a straight line whose slope is given as;
hAc
qBm
tφ234.0=
hAφ is nothing but pore volume of the reservoir. If this is represented as Vp, then
mCqB
Vt
p
234.0−= ft3
This particular test is called Reservoir Limit Test (RLT). Point should be remembered
that RLT is valid for only for pseudo-steady state only and not for steady state. There are
different ways to calculate the pore volume in a steady state condition.
Build-up Test: Horner’s Method Most of the information from a well test comes from the interpretation of a pressure
build-up. The reason is the fluctuation in the production rate which is inherent to the
production. Fluctuation may cause large variation in bottom hole pressure during draw
down test. This is not the case in a build-up test. In a build up test, the well is allowed to
flow at sufficiently large time to allow the flowing pressure almost constant.
Fundamentals of Reservoir Engineering & Characterization 104
Subsequently, the well is closed and the continuous recording of bottom hole shut in
pressure is done till the surface tubing shut-in pressure stabilizes.
The equation and analysis method was given by Horner.
The following expression for pressure is given ;
)log(6.162
)(t
tT
khqB
tpp pwsi ∆
∆+=∆− µ
The value of pressure measured at the bottom is plotted versus the logarithm of t
tTp
∆∆+
,
on a graph, once the wellbore storage effect has ended a straight line with a slope of m
can be observed
kh
qBm
µ6.162=
This helps to know the flow capacity (kh) of the well. The thickness h is called the
effective thickness and is obtained by subtracting the noncontributing length from gross
thickness of the formation encountered in the well. Skin is determined from the following
expression;
+−
∆∆+
+−
= 23.3log)log(151.12
1
wt
pwfHr
rck
t
tT
m
ppS
φµ
time
Tp
Pws
t∆
Fundamentals of Reservoir Engineering & Characterization 105
log of t
tTp
∆∆+
is considered negligible while determining skin through Horner’s method
and p1Hr, must be calculated from the Horner straight line at hrt 1=∆
Extrapolated Pressure
If the slope of the Horner’s straight line is extrapolated at t
tTp
∆∆+
=1 (i.e. when ∞∆t ),
the value of the pressure read, is called initial reservoir pressure (P*) in most initial tests,
where amount of fluid produced before shut-in is usually negligible compared with the
amount in place. The idea is that if the build-up would have been continued for infinitely
long time, the pressure would have stabilized to initial reservoir pressure. However,
when substantial amount of oil has been produced, the value of P* is not the reservoir
pressure existing at that point of time, rather this value is used to calculate the average
reservoir pressure. There are conditions, when the value of P* is found to be less than
average reservoir pressure −−
P ! So reservoir engineers should use this value with great
caution and understanding.
Miller Dyes and Hutchinson Method of Build-Up
Horner showed that build-up varies linearly with log(t
tTp
∆∆+
). When Tp >>> t∆ , the
term tTp ∆+ can be approximated as PT . Physically it means that the pressure drop due
to previous production is neglected.
Fundamentals of Reservoir Engineering & Characterization 106
Hence the Horner’s equation becomes )ln(ln6.162
pwsi Ttkh
qBpp −∆−=− µ
This equation was proposed by Miller, Dyes and Hutchinson and the particular method
of build up is called MDH method.
Pressure Shapes and Interpretation Methods in Various Characteristic Boundaries When compressible zone created by the perturbation reaches reservoir boundary, it is
perceived as a characteristic response in the pressure at the well. This nature of the
response in the well bore pressure depends upon the characteristics of the boundary.
Few of the characteristic responses observed in different types of boundaries are
explained below;
Linear Sealing Fault The boundary condition corresponding to linear fault is the linear no-flow boundary.
Linear sealing fault and disappearing facies, unconformities are few of the examples of
the characteristic boundaries. In such type of situation two different straight line
segments are seen with slopes having approximate ratio of 2:1 in the semi-log straight
line. The flow capacity and the skin should be calculated on the basis of first line.
However, P* should be calculated based on the second straight line in case of only one
fault. Flow capacity in both the drawdown and the build-up should be calculated based
on the following data;
m
qBkh
µ6.162=
Skin in drawdown
time
Tp
Pws
t∆
p∆ pMDH∆
Fundamentals of Reservoir Engineering & Characterization 107
+−
−= 23.3log151.1 2
wt
ihri
rck
mpp
Sφµ
Skin in build-up
+−
−= 23.3log151.1
2
1
wt
wfHr
rck
m
ppS
φµ
Gray suggested that the distance to the fault or barrier can be approximated using the
following equation.
t
t
ctk
Dφµ
∆= 0328.0
Where; D = Distance to the barrier or fault, ft K = Formation permeability, mD φ = Porosity, fraction = Fluid viscosity, cP ct =Total compressibility, psi-1
=∆ tt End of first straight line segment, hr If the two barrier/faults are approximately the same distance, the characteristic doubling
of slope will not be seen in the plot. In such case after the initial straight line is seen, the
slope of the second line would increase to more than two times. In such case the second
line suggests presence of more than one fault.
In a type curve the derivative of the slopes goes up from 0.5 to 1.
Pressure Build-up Data from a Well Producing from a Long Narrow Reservoir Such as Channel Sand The pressure transient data collected from a well producing from a long narrow reservoir
as shown below have characteristics that show combination of radial flow and linear flow.
During radial flow the pressure varies as logarithm of time. In a linear flow the pressure
varies linearly with square root of time. The channel can be due to number cases such
as;
1. Two parallel sealing faults.
2. a sedimentary deposit channels.
3. two parallel lateral variations in facies. etc.
w
d
Fundamentals of Reservoir Engineering & Characterization 108
The channel is defined by its width w and by the distance, d, from the well to one of its
edges.
During a well test inside a channel, following characteristics in the pressure patterns are
observed;
• A semi-log straight line with stabilization of derivative at 0.5 is observed.
• As the compressible zone reaches the first edge of the channel, fault effect is
seen. The boundary has exactly the same effect as sealing fault in an infinite
reservoir. The slope of the line doubles. This is observed only when the well is
very off centered in the channel.
• When the compressible zone reaches the two edges of the channel, it expands
linearly parallel to the edges of the channel. The pressure varies linearly with
square root time. Plot of pressure vs. t shows a straight line suggesting of a
channel.
A plot of P vs. ttt ∆−∆+ should be made in case of build-up. If the late time data
becomes a straight line on this plot it along with doubling of slope in radial flow indicates
channel reservoir. P* is determined from the linear plot by extrapolating ttt ∆−∆+
to 0.
Linear flow is used to determine the width of the channel and the eccentricity of the well
The width of linear channel can be calculated by
tch
qBmm
wφ
1
2
638.0= ft :for a oil well
m1 = slope of ( pws vs.
∆∆+
t
tt plog ) psi/cycle
m2 = slope of (pws vs. ttt ∆−∆+ ) q = oil flow rate, STB/D B = oil formation volume factor, bbl/STB H = net pay thickness, ft φ = porosity, fraction Ct = total compressibility, psi-1
Fundamentals of Reservoir Engineering & Characterization 109
)1(
02.2 .1
2 w
avg
Sh
PqTZm
mw
−=
φft :for a gas well
Where
Q = Gas rate, mscf/day T = Reservoir Temperature, oR Pavg = Average pressure in the neighborhood of the well Z = Gas deviation factor Sw = connate water saturation
Pressure Build-up Data from a Hydraulically Fractured Well Natural fractures are distributed homogeneously in the reservoir. Artificially fractures are,
however, located in the vicinity of the well bore. They are created by the operations
carried out on the well. They are an effective technique for increasing the productivity of
damaged wells or wells producing from low-flow-capacity formation.
Fractures can be created both in vertical and horizontal direction. At depths of less than
1000m it is possible to create horizontal fractures. However, at great depths, the
overburden weight makes the fractures develop only along vertical planes.
Flow around an Artificially Fractured Well The presence of an artificial fracture modifies the flows near the well bore considerable.
However because of the short distance extension of the fracture, these fractures have
finite conductivities, unlike natural fractures which have infinite conductivities.
In an artificially fractured well, initially, there is a fracture linear flow. This period is quite
short and is normally dominated by wellbore storage. Flow from the reservoir causes the
matrix to contribute to the flow of fluid to the fracture. This period is featured by linear
flows in both fractures and the formation and the fracture tip still has not affected the flow
behaviour of the well. These bilinear flow regimes are experienced only by fractures of
finite conductivity. Bilinear flow is followed by linear flow. During the start of this flow
period, the flow behaviour starts getting affected by the fracture tip. There is a linear flow
from matrix to the fracture. This flow is very often seen during testing of artificially
fractured wells. Finally, at, long times the pseudo-radial flow is reached by all fractured
systems regardless of the fracture conductivity or damage. The system developed for
Fundamentals of Reservoir Engineering & Characterization 110
the radial homogeneous system is equally applicable for interpreting data of this flow
period, albeit with minor modification.
Flow Model for Each Flow Pattern a. Linear flow in the Fracture The flow exists theoretically at the very beginning of the test. During this flow most of the
fluids produced at the well come from expansion in the fracture. The flow is linear. The
pressure varies linearly with t
The variation can be expressed as
Dxfr
D tC
P ηπ2=
or
ftf
wfi Ckt
whqB
pp)(
128.8φµ=−
Fracture
Xf
Fundamentals of Reservoir Engineering & Characterization 111
where
2
0002637.0
ftDxf xc
ktt
φµ=
η is the ratio of diffusivity inside the fracture and diffusivity in the reservoir.
And Cr is the relative conductivity and is expressed as
kx
wkC
f
fr =
The greater the relative conductivity of the fracture more effective and pronounce this
flow regime is seen on the plot. A fracture with relative conductivity of over 100 behaves
as if it had infinite conductivity. At low fracture conductivity, linear flow regime is not seen.
The concept of relative conductivity explains why the smaller the formation permeability,
the more effective the hydraulic fracturing is.
b. Bilinear Flow It is called bilinear because it corresponds to two simultaneous linear flows;
• an incompressible linear flow in the fracture
• a compressible linear flow in the formation
Bilinear flow lasts as long as the ends of the fracture do not affect the flows. This
flow period occurs only in case of finite fracture conductivity cases and where there is no
well bore storage distortion. In this flow regime, the pressure behaviour is featured by
the linear relationship when data are plotted by using the pw and t1/4 coordinates
4/14/12/ )()(
1.44t
kcChqB
pfWt
irf φµ
µ=∆
The equation suggests that slope of the bilinear plot would lead to the estimation of kfw
and fracture half length. However, it should also be noted that determination of fracture
characteristics by this method requires knowledge of reservoir properties.
c. Linear Flow in the Formation This flow is very often visible during testing of artificially fractured well. It is an integral
part of the conventional analysis methods of these tests. The flow regime occurs in the
Fundamentals of Reservoir Engineering & Characterization 112
fracture itself and in the formation proper. This type flow is exhibited by only highly
conductive fractures (Cr > 100). This flow period if exists, should be used for calculation
of fracture properties. It is characterized by a linear variation of the pressure versus t
The flow is characterized by following expression
DxfD tp π=
or kx
tch
qBpp
ftwfi φ
µ064.4=−
d. Pseudo Radial Flow At long time and end of bilinear and linear flows, pseudo-radial flow regime starts. The
reason why it is called pseudo radial flow is that flow period is not fully radial (Russel and
Truitt). Nevertheless, all curves approach a common value of maximum slope which is
dependent on the length of the fracture penetration. Raghvan et. al. constructed a graph
of correction factors fc , which must be used to obtain the correct permeability factor.
i.e. cH fkk = Russel-Truitt Method of Permeability Determination from Pseudo Radial Flow Russel Truitt method for the determination of true permeability of the reservoir is given
below;
The graph is a plot between R=(Measured slope of build-up data/theoretical slope of
build up data) versus ( Lf/D=fracture length/spacing between wells). A prototype of slope
is shown below and is not to the scale;
Procedure
R
Lf/D
Russel-Truitt Plot for slope Correction
Fundamentals of Reservoir Engineering & Characterization 113
For an oil well, the equation relating fracture length with reservoir test parameter is given
as;
Rch
qBmm
Lt
f φ1
2
638.0=
where
Lf = fracture length (tip to tip), ft M1 = slope of Pwf Vs log t∆ plot psi/cycle for drawdown
slope of Pws Vs log
∆∆+
t
tt plog plot psi/cycle for drawdown
M2 = slope of the pwf vs t∆ for drawdown
slope of the pws vs ttt ∆−∆+ for build up φ = porosity, fraction R = Correction factor from Russel-Truitt plot
Steps to solve: Step 1. : Assume an Lf value
Step 2. : Calculate Lf/D value
Step 3. : From the graph of Russel-Truitt calculate R
Step 4. : Calculate Lf from the equation
Step 5 : The assumed value and the calculated value if found equal gives the correct
value of fracture length. Otherwise, repeat the iterative process.
Step 6: Put the value of R in following equation to know correct value of permeability of
the formation
1
6.162m
RqBkh
µ=
Flow Pattern in a Closed Reservoir When the reservoir is limited by no-flow boundaries it is called closed reservoir i.e. when;
0=∂∂
rp
at r= re i.e. no influx
tConstp
tan=∂∂
for all r and t
Fundamentals of Reservoir Engineering & Characterization 114
The beauty of this pseudo-steady state regime is that it helps to define the drainage area
of a well. Drainage area may be due to;
• physical barriers: sealing faults, disappearing facies, etc.
• production from neighboring wells: The boundary between two wells is
proportional to the pore volume drained by each well
Solution of diffusivity equation in the pseudo-steady state regime is given as below;
)2458.2
ln(21
ln21
22
AwDAD Cr
Atp ++= π
or
+++=− S
CrA
khqB
thAc
qBpp
Awtwfi 87.0
2458.2loglog
6.162234.02
µφ
Where
Ac
ktt
tDA φµ
0002637.0=
A = is the drainage area of the well CA = is a shape factor that depends on the shape of reservoir and
the position of the well in it. A table with shape factor corresponding to different well configuration is shown below;
Fundamentals of Reservoir Engineering & Characterization 115
Method to Calculate Shape Factor As evident from the pseudo-steady state equation, a plot of pressure vs time on linear
scale would result in a straight line with slope mL as shown in the figure below;
The slope ML is used to determine the drainage area or the pore volume drained hφA.
Lt hMC
qBA
234.0=
The value of CA can be calculated in the following way;
[ ]mPPmm
C inithrL
A /)(303.2exp456.5 1 −−⊗=
where,
m is from radial transient flow and mL is determined from linear plot
CA value should be compared with the chart to find out the drainage shape. Determination of Average Reservoir Pressure When the compressible zone reaches real no-flow physical boundaries during build-up,
the pressure in the drainage area becomes uniform and constant. The pressure is called
the average pressure of the drainage area.
Matthews, Brons and Hazebroek (MBH) method
Matthews, Brons and Hazebroek calculated AC
ktt
t
ppDA φµ
0002637.0= for various reservoir
well configuration and plotted against m
PPPDMBH
)(303.2 * −= . One of the plot for
rectangular area is shown below;
Pressure
Slope=ML
Elapsed time
Pinit
Fundamentals of Reservoir Engineering & Characterization 116
Where
tp is the production time
P* is the extrapolated pressure from the semi log straight line
m= slope of the semi-log straight line
A= drainage area.
Steps to calculate average reservoir pressure Step 1: From the known drainage area, tpDA is calculated
AC
ktt
t
ppDA φµ
0002637.0=
Step2 : From the curve shown above m
PPPDMBH
)(303.2 * −= is calculated based on the
reservoir well configuration, which can be known from CA.
Fundamentals of Reservoir Engineering & Characterization 117
Step 3: From known m and P* average reservoir pressure is calculated. Deitz Method to Calculate Average Reservoir Pressure Following steps have been suggested for the calculation of average reservoir pressure.
Step 1: P* is calculated from semi-log staright line.
Step 2.: Vi is calculated based on following equation
TOT
TOTii Q
VQV =
Where,
VTOT can be calculated from geological structure map.
QTOT is the total rate from the sand
Step 3: From Vi area A can be calculate by dividing by average thickness.
Step 4: tDA can be calculated as follows
AC
ktt
tDA φµ
0002637.0=
Step 5: Value of St∆ i.e. start of pseudo-steady state time can be calculated from the
following equation;
DAAS
sp tCt
tt=
∆∆+
Step 6: The value of pressure at St∆ read from Horner plot gives the value of average
reservoir pressure.
P* P
Fundamentals of Reservoir Engineering & Characterization 118
Russel Method for BUP Interpretation Some time well bore storage effect affects or distorts the Horner plot so that we don’t get
the Horner straight line which is characteristic of radial flow. This makes the BUP
interpretation useless. Russel suggested a solution to this problem by following method.
Step-1 : Plot
tC
p
∆−
∆1
1 versus log of t∆
Where )()( tPtpp wfws −∆=∆
Step 2: Vale of C should be chosen such that the plot we get a straight line
Step 3 : Slope of the straight Russel line would give the permeability value.
slope
qBkh
µ6.162=
Step 4: Skin should be calculated as follows
+−
∆−
−= 23.3
)(log
)
11
(
)(151.1
2
1
wt
wfhr
rCk
slopetC
PPS
φµ
Horizontal Well Testing Procedures Horizontal well testing is complex and on many occasions it is difficult to interpret.
Detailed discussion would require an exhaustive treatise. Hence, the discussion would
Correct C
C too small
C too large
Fundamentals of Reservoir Engineering & Characterization 119
be limited to only the fundamentals so that one can apply in solving actual field related
problem. There are four transient flow regimes that are theoretically possible with a
build-up or drawdown test in a horizontal well. They are as follows;
Early Time Radial Flow The flow is radial and is equivalent to that of a fully penetrating vertical well in an infinite
reservoir.
Intermediate Time Linear Flow A horizontal well will generally be long compared to the formation thickness; a period of
linear flow may develop once the pressure transient reaches the upper and lower
boundaries.
Late Time Radial Flow If the horizontal well length is sufficiently small as compared to the reservoir size, a
second radial flow known as a pseudo-radial flow will develop at late times.
Late Time Linear Flow This flow period occurs when the pressure transient reaches the lateral extremities of the
reservoir. The intermediate time linear flow and late-time linear flow period develops only
for reservoir of finite width. The identification of these flow regimes is critical to the
proper interpretation of a horizontal well test.
Pressure Response Equations for Different Flow Regime Early time Radial
+−
=− s
rC
tkk
Lkk
qBPP
wt
yv
yv
wfi 868.023.3log6.162
2φµµ
for Drawdown
+
∆∆+=− 1log
6.162 γµt
tT
Lkk
qBPP
WyZ
wsi for BUP
Where
Fundamentals of Reservoir Engineering & Characterization 120
SrC
KKt
LCtk
KK
hL
wt
yv
wt
x
X
v 869.0227.3log)log(023.2log221 ++
−−
−
=
φµφµγ
Intermediate Time Linear FLow
++
=− )(
2.141128.8Ss
KKL
qBCK
tLh
qBPP Z
Vytywfi
µφµ
for Drawdown
Where 838.1180lnln25.0ln −
⊗−
+
=
hZ
SinK
K
rh
S w
V
y
wZ
+∆=− 3
128.8 γφ
µty
wsi CKt
hLqB
PP for Build-up
Where
+
= 023.2log
6.16223 LC
tk
kkh
qB
t
x
yx φµµγ
Late Time Radial Flow
++−
=− )(
2.141]023.2[
6.1622
SsKKL
qBLC
tkhkk
qBPP Z
Vyt
x
yxwfi
µφµ
µ for drawdown
∆∆+=−t
tT
hkk
qBPP
yx
wsi log6.162 µ
for Build Up
Late Time Linear FLow
)(2.141
2128.8
SSSkkL
qBckt
hxqB
PP zx
vytyewfi +++=− µ
φµ
for Draw down
∆−=− )(
128.8tt
CKhhqB
PPtyx
wsi φµ
for Build-up
Gas Well Testing
Fundamentals of Reservoir Engineering & Characterization 121
The gas well testing differs from a well testing fundamentally. The basis on which the
diffusivity equation was derived doesn’t hold good for gas. Unlike oil, gas viscosity and
compressibility vary widely with pressure. Darcy equation for gas flow can be written as;
pZ
p
rr
TP
khTq
R
Wf
P
P g
w
eSC
scg ∂⊗= µ
π2
)ln(
2
The general trend for Z
p
gµ versus pressure is given as follows;
Region I :
Region I which is less than 2000 psi, the pressure function Z
p
gµ shows linear
relationship with pressure. Hence Zgµ
1 can be taken as constant at low pressure, in
that case
Z
PPp
Zp
g
wfRP
P g
R
wfµµ
22
2−
=∂
Hence Darcy equation for gas at low pressure <2000 psi becomes
( )
)75.0(ln(
703.0 22
srr
ZT
PPkhq
w
eg
wfRg
+−
−=
µ
At high velocity flow Forcheimer modified the above Darcy equation to include the rate
dependent skin
Zp
gµ
P, Psi
I
III II
2000 3000
Fundamentals of Reservoir Engineering & Characterization 122
( )
)75.0(ln(
703.0 22
sDqrr
ZT
PPkhq
w
eg
wfRg
++−
−=
µ
Region II In region II, where the pressure is in between 2000 to 3000 psi the pressure function
shows distinct curvature. In this region, the concept of pseudo pressure should be used.
Pseudo pressure is defined as ;
∆= PZ
PP
gµψ 2
)(
Flow equation becomes
( )
)75.0(ln(
703.0
sDqrr
T
khq
w
e
wfRg
++−
−=
ψψ
Region III
Which is a high pressure region, higher than 3000 psi, the pressure function Z
p
gµ is
constant. Hence
)(2
2 wfRg
P
P g
PPZ
Pp
ZpR
wf
−=∂ µµ
Darcy equation becomes
( )
)75.0(ln(
406.1
sDqrr
T
PPZ
Pkh
q
gw
e
wfRg
g
++−
−
=µ
Hence, while interpreting gas well test data, use of correct type of pressure function
must be remembered.
Absolute Open Flow Potential To know the absolute open flow potential (AOFP) of a gas well is one of the most
important parameter for predicting gas production profile. The AOFP is defined as the
theoretical production flow rate of a well reached with the bottom hole pressure equal to
Fundamentals of Reservoir Engineering & Characterization 123
atmospheric pressure. This measurement is done in the pseudo-steady state region.
Well is allowed to flow through three different gas rates. AOFP is calculated in the
following way;
In the pseudo-steady state region, the flow equation can be expressed as 1. Pressure less than 2000 222
ggWfR BqAqPP +=−
Where
+−= s
rr
kh
ZTA
w
eg 75.0ln703.0
µ
Dkh
ZTB g
703.0
µ=
2. Pressure greater than 3000 2
ggWfR BqAqPP +=−
Where
+−= s
rr
khP
ZTA
w
e
Avg
g 75.0ln406.1
µ
DkhP
ZTB
Avg
g
.406.1
µ=
2
wfRavg
PPp
+=
3. Pressure in between 2000 to 3000 psi 2)()( ggwfR BqAqPP +=−ψψ
Where
+−= s
rr
khT
Aw
e 75.0ln703.0
Dkh
TB
703.0=
Fundamentals of Reservoir Engineering & Characterization 124
AOFP is defined as
B
BAAAOFP a
2
)(42 ψψ −++−=
Composite System In-situ combustion process, steam flooding process, polymer flooding process etc., give
rise to composite system, where mobility contrast exist within the reservoir. A steam
flood or in-situ combustion process is modeled as a two-region reservoir, with an inner
swept region surrounding the injection well and an infinitely large unswept region beyond
the front. Figure below shows a typical composite system.
Mobility contrast which exists in between the region I & II is used to model the fluid flow
in such type of composite system. This acts as basis for determination of the swept
volume using pressure transient model. The swept volume adjacent to an injection well
is considered to have both well bore storage and a skin effect. This swept volume has
different permeability, porosity and compressibility of the reservoir fluid then the zone II
ahead of it. A complication of non uniform temperature i.e. adiabatic condition does arise
while modeling a combustion system or a steam flooding, which is just opposite to the
assumption used for derivation of diffusivity equation. In an in-situ combustion process
the temperature in the region adjacent to the well bore will be that of injected air whereas
near to the combustion front it would be as high as 10000 F whereas, in case of steam
flood temperature adjacent to well bore, the temperature would be that of steam while,
that of farthest away near to outer boundary of the swept region, it would be equal to the
reservoir temperature. The solution to modeling such type of problem is found by
assuming swept region to exist at some mean temperature.
I
II Well Bore
Top View
Fundamentals of Reservoir Engineering & Characterization 125
Modeling of the Fluid Flow in a Composite System In any composite system, as explained in the previous paragraph, there will be swept
region from the injection sand face to the displacement front as shown in fig above.
Region-I will be dominated by the injected fluid which can be steam in case steam
flooding, injected air in case of forward in-situ combustion or polymer slug in case of
polymer flooding. Region-II is the zone representing the zone ahead of the displacement
front.
Modeling of the fluid flow in a composite system includes following assumptions;
1. The formation is horizontal, uniform thickness and is homogeneous.
2. The front is of infinitesimal thickness in the radial direction.
3. Flow is radial, and gravity and capillarity effects are negligible.
4. During the well test (Fall-off or Injectivity test) the front is considered to be
stationary.
5. The region behind the front contains only gas in case of in-situ combustion or
steam in case of steam flooding.
6. The fluid is slightly compressible.
The diffusivity equation derivation methodology is same as that of homogeneous system
as explained earlier while, deriving the expression for homogeneous system. In
dimensionless form, the diffusivity equation for two different regions can be written as
shown in the following paragraphs. The reason why the diffusivity equation is written in
dimensionless form is that, it permits to understand the structure of the solutions of
interest without consideration of the specific values of formation properties, fluid
properties or flow rate. The objective here is to obtain a solution that contains no
parameters.
Region-I
D
D
D
DD
DD tP
rP
rrr ∂
∂=
∂∂
∂∂ 111
--- (1)
Region-II
D
D
D
DD
DD tP
rP
rrr ∂
∂=
∂∂
∂∂ 221 η --- (2)
Where :-
Fundamentals of Reservoir Engineering & Characterization 126
PD1 = )(2
11
1 PPqB
hki −
µπ
& PD2 = )(2
22
2 PPqB
hki −
µπ
211
1
1 WtD
rC
tkt
µφ= & rD =
Wrr
21
=tt c
kc
kφµφµ
η
Eqn. 1 and 2 along with initial and boundary conditions, can be solved analytically in
cylindrical coordinates using, Laplace inverter to generate dimensionless bottom hole
pressure data, Pwd, as a function of dimensionless time, tD. The simulation of Pwf function
against time t shows a semi log straight line on a semilog plot for region – I (Swept
region) followed by a break at the front with another straight line having different slope
for region-II. The semilog slope of the first line as shown in fig below gives an idea
about the permeability and skin of the swept zone.
mqB
khµ6.162=
and
+
−
−= 23.3log1513.1
21
wt
hrw
rck
mPP
sφµ
Where m is the slope and P1hr is the extrapolated pressure at one hour shut-in time for
the semilog straight line. The parameter B, , Ø and Ct corresponds to rock and fluid
properties in the swept volume and rw is the radius of the well.
After the semilog straight line, the system starts to react to the radial discontinuity at a
distance where the front lies or where the mobility contrast is highest. The zone near to
front behaves like an impermeable boundary due to the high mobility contrast. As a
result pressure rises above the semilog straight line as seen in fig above. During the
Fundamentals of Reservoir Engineering & Characterization 127
transition period, the system behaves or approximates pseudo steady state flow as
shown by the Cartesian plot in fig.below . This region of pseudo steady state can be
used to calculate the swept zone rock pore volume. The pore volume is related to the
slope of the pseudo steady state Cartesian straight line as follows;
mC
BIV
t
ga
*
**234.0−=
Where Ia is injection rate, Bg is the air or steam formation volume factor, Ct is the total
compressibility of the swept region and m is the Cartesian slope. The above calculation
of swept region for in-situ combustion or steam flooding is easy said than done. Because,
the above calculation is for isothermal condition whereas, in in-situ combustion or steam
flooding the process is non-isothermal. To solve for the swept pore volume, permeability
and skin formation volume factor Bg and the total compressibility Ct which are
temperature and pressure dependent, have to be estimated from the average pressure
and average temperature of the swept zone(region-I). Once the swept volume is known,
calculation of the fuel concentration for a combustion process or the cumulative heat
loss from a steam flood is possible.
Procedure for calculation of K, s and fuel concentration in a swept zone (In-situ combustion process) The basic procedure is to calculate the average reservoir pressure, temperature and the
swept volume simultaneously. Also the permeability thickness, and skin factor, s, can not
be calculated because the air properties, Bg and Cg are not known until the average
reservoir temperature and pressure is found. The procedure for in-situ combustion
process is explained as given below;
Pressure
Slope=ML
Elapsed time
Pinit
Fundamentals of Reservoir Engineering & Characterization 128
• Plot the pressure and time data on a semi-log and Cartesian plot.
• Estimate the average reservoir pressure behind the front from the early time
flattening of the semi-log plot.
• From the Cartesian slope find out the slope of pseudo-steady state straight line,
m1.
• Calculate dimensionless time, t2θ for in-situ combustion process from the
following formula;
( )
tCh
Kt
bOB
OBb
2
2
)1(
2
−=
ρφαθ
where KOB = overburden thermal conductivity h = Thickness of the pay zone, ft α = Thermal diffusivity of the cap rock t = total injection time bC)(ρ = effective specific heat of the swept region.
• Calculate thermal heat efficiency as follows
+−= )(1
21)( 2
2
22 2
terfcet
ttE t
h θπ
θθ
θ θ
• Assume average temperature behind the front
• Find Bg and Cg.
• Calculate swept volume as follows
g
ga
Cm
BIV
1
=
• Calculate total volume behind the front φV
Vb =
• Calculate the average temperature, −T as follows
rb
faaaF
cah
f CV
TTCFH
tIE
TT))(1(
)()(
ρφ
ρ
−
−+
+=−
• Assume different average temperature and repeat the process till the subsequent
swept volume do not change
• With the value of Bg and Cg calculate permeability thickness and skin
Fundamentals of Reservoir Engineering & Characterization 129
m
qBkh
µ6.162=
+
−
−= 23.3log1513.1
21
wt
hrw
rck
mPP
sφµ
• Calculate fuel concentration
1VF
tIC
aF
am
φ= where FaF is air fuel ratio.
Procedure for calculation of K, s and cumulative heat loss in a swept zone (Steam flood process) An important factor for a steam flood is the amount of heat that has been lost to the
overburden. Knowledge of the steam swept volume from a pressure transient well test
enables calculation of the heat loss. The procedure for the steam flood is simpler than
the in-situ combustion because the average reservoir temperature is known.
• Plot the pressure-time data on a semi log graph and Cartesian graph.
• Find the average reservoir pressure behind the front from early time flattening of
the semi log curve. From the steam table estimate the average swept zone
temperature. Find the slope of the early time semi log straight line. Calculate
the permeability thickness in the swept region
m
qBkh
µ6.162=
• From the semi log graph calculate the skin factor
+
−
−= 23.3log1513.1 2
1
wt
hrw
rck
mPP
sφµ
• From the Cartesian plot find the slope m1 of the pseudo steady state straight line
• Calculate swept volume, V1
g
gss
Cm
BIV
11 =
• Heat loss can be calculated as follows
( )
[ ]swws
Th HHtq
VHE
Γ+Γ−=
)1(1
φρρ
where Γ =Steam quality qs = steam feed water H = Enthalpy
Fundamentals of Reservoir Engineering & Characterization 130
( Hρ )T = Total heat content of the swept zone This brings to an end to the well test concept applied to the in-situ and steam flooding
processes. It is evident from the above analysis that the pressure transient fall off well
test of thermal injection wells based on the above mentioned model produces potentially
useful results. However, it is pertinent to mention here that the accuracy of the result will
depend upon how accurately we identify the transient and pseudo steady state period in
the swept zone. The case mentioned above is an ideal one. Accurate determination of
different region requires derivative analysis of the pressure w.r.t. time.
MANAGEMENT OF OIL WELL TEST After having discovered oil/gas pool, it becomes critical to know reliable information
about in Situ reservoir conditions. A proper understanding of the reservoir and fluid
properties is essential for cost effective and efficient development planning. Having
spent enormous amount on exploratory and a. delineation drilling activity to prove up the
reserves, it is negligible to leave the well without establishing data that will be required
for planning the exploitation of the reserves. There ate numerous cases where operators
have had to reenter or redrill a well or, worse still, have installed ill-designed facilities
and preceded with an uneconomic development as a result of inadequately planned,
insufficiently long, poorly supervised or misinterpreted well tests. It is pertinent here to
note the difference between conventional major fields and frontier marginal fields is that
large capital investment has to be made for frontier/marginal field’s development based
almost completely on exploration and delineation well data. While the conventional field
development cases, the data can be refined in a phased manner using the latest drilling
and production results and for frontier marginal field development cases require a
commitment to spend majority of funds long before any production history is available. It
amounts to that the data obtained from exploration and appraisal wells must be
comprehensive and of the best quality possible.
WELL TESTS GENERAL Well testing is a process used by the petroleum industry to solve problems and answer
questions related to the operations and economic evaluation of hydrocarbon reservoirs
and their associated wells. Two general conditions exist within the industry with respect
to the nature of well testing activities.
Fundamentals of Reservoir Engineering & Characterization 131
One most popular connotation in terms of the type and frequency of test occurrence, is
that a well test is an observation of a well's productivity i.e. production or injection rate as
a function of bottom hole or surface flowing pressure.
The second connotation of well testing, as seen mainly through the eyes of engineering
segment of the industry is that a well test is a definition and quantification of the
parameters which control a well's productivity, i.e. static drainage. area, pressure,
permeability, skin, etc. The advantages of second approach to well testing includes:
The ability to determine the accuracy of a well's observed productivity.
The ability to determine the stability of a well's observed productivity.
The ability to determine the impact of changing the parameters which control the
productivity of a well or an entire reservoir.
RESERVOIR ROCK PROPERTIES Well tests can give reliable estimates of reservoir rock properties such as :
Capacity (Kh) : For predicting well productivity, estimating net pay open to flow,
correlating with core data, predicting reservoir stratification and establishing fracture,
stimulation requirements.
Skin(s) : used for estimating well Bore damage and essential for predicting well
productivity and evaluating stimulation potential and results.
Drawdown (Delta P) : used for defining productivity index of the well and evaluating well
bore conditions.
Production Characteristics: These are needed for production forecasting, designing
well completions and sizing top side facilities in particular the following data is needed.
• Inflow Performance Curve or Absolute Flow Potential: For gas wells essential for
production forecasting.
• Tubing Performance Curve: needed to size production tubing and gathering
system.
• Sand Production: Important in designing production and injection well
completions specifically gravel packs.
Fundamentals of Reservoir Engineering & Characterization 132
• Potential Problems: Waxes, sulphur, scaling, corrosion, and hydrates needed for
designing well completions and facilities.
The types of information available from pressure transient tests along with
economicaJ1y significant benefits of obtaining this information are presented in Tables.
Most of these tests are of productivity observation variety, but could be easily and
economically converted to pressure transient test variety with significant potential value
to the industry.
To recapitulate and summarize what has been talked of in the preceding sections, the
data generated from well tests and their utility is summed up below.
DATA REQUIREMENT AND DATA GENERATED FROM WELL TESTS
Most of the data required for evaluation and valuation of a reservoir would be generated
from well tests. The main data requirement expected from a production test programme
is summarized below with their utility and relative importance of such data.
FLUIDS
It is of utmost importance to identify and obtain representative samples of fluid contents
of the reservoir be they oil, gas, condensate or water. These are needed for geological
modeling, predicting fluid contacts, recovery prediction, and formulation of reservoir
depletion plan, production facility design and PVT behavior of the reservoir fluids.
RESERVOIR BOUNDARIES AND HETEROGENITIES
Comprehensive well test data sometimes can provide valuable information about nature
and size of the reservoir being tested. Specific information obtainable from well tests is
fractures, limit of reservoir like pinch outs, nearby gas cap, nearby faults, nearby aquifer,
stratification and inter-block communication. These are the areas of uncertainty can
usually be estimated by an extended production testing by investigating for several days.
When there is doubt about the size of the reserves, extended production testing is the
only answer to gain confidence on the reserves for development decision.
Fundamentals of Reservoir Engineering & Characterization 133
COST EFFECTIVENESS AND PROPOSED MANAGEMENT OF WELL TESTS
It is frequently impractical and not at all times to get all of the data indicated above owing
to various logistic problems. Certain guidelines as per their rank in importance is
indicated in Table. A technical recommendation and management decision has to be
made as to whether to spend the time and money needed to obtain certain items of
information. The recommendations have to be purely based on the need of the situation.
For instance a reservoir boundary is suspected from seismic and other geological
information which is critical to estimate minimum reserves size needed for development,
an extended test should be considered. It would be difficult to calculate cost
effectiveness for petroleum engineer to quantify the cost of not knowing the correct
reservoir fluid compositional analysis, because this missing data will have an impact on
the recovery predictions. The depletion plan the facilities design and ultimately on the
project cash flow. The development plan may turnout to be either too optimistic or
pessimistic. The facilities accordingly will be either under designed or over-designed.
This situation would result into non-optimization
of exploitation strategy. Now-a-days, sophisticated computer modeling tools are
available which would help ' in checking sensitivity of the project cash flow to certain key
assumptions. This can help to quantify cost effectiveness of obtaining certain data but
will not provide the total answer. The bad development would mean recovering less oil
and gas than what would have been expected but how much and at what cost? Under
these conditions, a judicious decision has to be taken depending on the situation, as to
what data is a must and rank remaining information as needs. Thus meeting the
requirement of cost effectiveness.
The problems arise for testing sour oil/gas wells because of concern associated with
high costs and risks in testing. Normally, there will be reluctance to test these wells even
if they are tested, the duration will be for a short time because the completions might not
be designed to overcome the bad effects,. This situation would result into missing of vital
reservoir information which would result into more assumptions.
TESTING GUIDELINES
Having been convinced of the importance of the data generated from well tests, the
following guidelines are given for obtaining the data through various means.
Fundamentals of Reservoir Engineering & Characterization 134
WIRELINE FORMATION TESTING
The repeat formation tester (RFT) is a well tried and proven testing tool which can
provide valuable information quicker at less time than DST or conventional production
tests. The pressures are very useful in identifying different reservoirs, depletion levels of
the reservoirs and geological zones.
DRILLSTEM AND SHORT TERM PRODUCTION TESTING
It is a short term test conducted in a well. These can be run in open hole under
cemented casing under tubing and permanent packers. Successful well testing in frontier
wells consists of finding the correct balance between two opposing needs - obtaining
maximum collection of relevant data with minimum amount of expensive rig and support
costs.
FLOW AND BUILD-UP PERIODS
Adequate pretest planning is required for estimating number and length of flow and
build-up periods. If log, core, wire line formation test or nearby offset well data is
sufficient, flowing and build-up periods can be specifically specified using fluid flow
equation, i.e. by determining stabilization time. While designing the test period the
following should be kept in mind.
• The time required to eliminate well bore storage effects for both drawdown and
build-up testing.
• The time required for semi-log analysis techniques to be applicable.
• The time when flow conditions change from transient to semi-steady-state,
expected flow rates under both flow regimes and radius of investigations at
different times.
In the absence of any specifically designed tests, the following guidelines are suggested.
Initial flow of 15-30 minutes is required to allow equalization of the filtrate invaded
zone back to static reservoir pressure.
It should be followed by 1.0 to 2.0 hours shut-in to" obtain reliable estimates of
initial reservoir pressure and temperature gradients should be seen.
Fundamentals of Reservoir Engineering & Characterization 135
Clean up period should continue until the tubing head pressures and
temperatures, gas-oil ratios, water rate are reasonable stable.
If high drawdown are required to get intended test rates the choke size should be
progressively increased to safeguard against sand production.
Highly productive zones can be produced at high rates immediately to obtain
high tubing head temperatures to minimize hydrate formation and to accelerate
clean-up.
Clean-up rate should be more than the planned test rates to facilitate opening up
of maximum number of perforations.
The response of the well to choke sizes should be well conceived during clean
up, so that a suitable choke size can be chosen before putting the flow through
separators.
Frequent changes of chokes should be avoided which would make analysis
difficult if not useless.
OIL WELL TESTING
• Three flow periods are ideal to maximize reservoir data if there is time
constraint, two rates may be adequate.
• The drawdown to be created should be up to 40% - 50% of reservoir
pressure.
• At least four hours of stabilized flow rate should be adequate to get reliable
data.
• If specific information is needed like sand failure, casing, etc. the drawdown
should be higher to know the sensitivity of drawdown to sand cut.
• If due to operational constraints, the pressure. Build-up study is not amenable
for Horner’s Method, data should be interpreted by log-log curve matching
technique to get the feel of reservoir properties.
GAS WELL TESTING
Gas well testing should be essentially multi-rate flow tests (4 chokes) to obtain
reasonable estimates of flow performance and rate dependent skin effect.
Fundamentals of Reservoir Engineering & Characterization 136
Flow after flow tests or Back-Pressure tests be preferred if the reservoir
permeability is large.
Modified isochronal be chosen if permeability is low.
The build-up time should be approximately twice the cumulative flow time of flow
and clean up time.
If enough details available the time needed for applying semi-log analysis
technique can be applied.
If possible, the well should be dosed-in down hole to minimize well bore storage
effects.
CONCLUDING REMARKS
Well test engineering is the process of successfully deriving useful valuable information
from well tests in the form of problem diagnosis and or reservoir valuation. The tasks
required to perform well test engineering can be grouped into three categories:
l. Planning and Designing
2. Monitoring and Control
3. Interpretation and Diagnosis.
These activities have to be carefully and judiciously planned, executed and Interpreted,
the task of finding a model which adequately represents the physical situation existing in
the wells and reservoirs being tested and quantifying the parameters which are critical
parameters in planning,
development, predicting reservoir depletion and managing the reservoir during the
producing life to get best out of the reservoir.
To conclude, well tests would be able to generate very useful information by which the
"Definition" and "Evaluation" of the reservoir could be accomplished in a very meaningful
manner thereby leading to draw a rational development strategy. As it is amply clear that
reservoir is unbelievably complex and impossible to define completely, to arrive at a
diagnosis of the system, one has to rely on –
A few physical deterministic facts.
Production statistics often of doubtful reliability.
Fundamentals of Reservoir Engineering & Characterization 137
Samples representing approximately one billilanths of the reservoir.
The above equation is called the material balance equation. To condense the material balance into more understandable form Havlena and Odeh
employed following methods
Fundamentals of Reservoir Engineering & Characterization 214
( ) WfwgO WeBEmEENF +++=
In which [ ] WpgSPOP BWBRRBNF +−+= )(
( ) gSSiOiOO BRRBBE )( −+−=
−= 1
gi
gOig B
BBE
( )Wc
fWCWOifw S
pCSCBmE
−∆+
+=1
)(1
Material Balance Above the Bubble Point pCNBBWBN effectiveOiWPOP ∆=+
where
wc
fWCWOOeffective S
CSCSCC
−++
=1
Material Balance for Depletion Below Bubble Point Once the pressure falls below the bubble point , solution gas is liberated from the oil.
Morris Muskat presented the performance prediction of a depletion below bubble point
pressure.
Let us consider an initially gas saturated reservoir from which NP stb of oil has been
produced. Then the oil remaining in the reservoir would be
Nremaining = N-Np= O
O
BVS
Where V is the pore volume (rb0. The change in this volume with pressure is
P
BBS
VP
SB
VP
N O
O
OO
O
R
∂∂
−∂
∂=
∂∂
2
1
The total volume of dissolved and free gas in the reservoir is;
Fundamentals of Reservoir Engineering & Characterization 215
Gr = ( )g
WcOO
SO
BV
SSB
RSV −−+ 1
Its change in volume with pressure is given by
∂∂−−
−∂
∂−
∂∂
−∂
∂+
∂∂
=∂
∂p
B
BSS
pS
BpB
BSR
pS
BR
pR
BS
Vp
G g
g
wcoO
g
O
O
OSO
O
SS
O
Or22
)1(1
Hence producing GOR expression can be given as;
pB
BS
pS
B
p
B
BSS
pS
BpB
BSR
pS
BR
pR
BS
ROoO
O
g
g
WCoO
g
O
O
OSO
O
SS
O
O
∂∂
−∂
∂∂
∂−−−
∂∂
−∂
∂−
∂∂
+∂
∂
=
20
22
1
)1(1
The producing GOR can be approximated with Darcy’s law for GOR
Sgro
OO
g
rg RkB
B
kR +=
µµ
The above two equation can be equated to give
g
O
ro
rg
g
g
WcOO
g
o
ro
rg
O
OS
O
gO
o
k
kp
B
BSS
pB
k
k
BS
pR
B
BS
pS
µµ
µµ
+
∂∂−−
−∂
∂+
∂∂
=∂∂
1
)1(
The above equation can be calculated for change in SO for depletion below bubble point.
At any stage of depletion the oil saturation is related to recovery in following way
)1()(
wcOi
OPO S
NBBNN
S −−
=
Giving
o
oi
WC
OP
BB
SS
NN
−−=
11
Material Balance Method in a Gas Reservoir
Fundamentals of Reservoir Engineering & Characterization 216
The initial gas in place can be determined without even knowing A, h, φ or SW provided
enough pressure production history is available. This can be done due to following
relationship;
As per mole balance on the gas;
Moles of gas produced = Initial moles of gas - moles of gas remaining
As per real gas law following substitution can be done in the above relation;
[ ]
ZRT
WWVP
RTZVP
RTGP pe
i
i
SC
PSC)( −−
−=
Assuming no water production for a volumetric reservoir the above relation reduces to;
VZTP
VTZ
PT
GP
i
i
SC
PSC
−
=
Or
PSC
SC
i
i GVTTP
ZP
ZP
−=
The above equation is an equation of straight line in terms of (P/Z) vs GP Havlena-Odeh Interpretation Havlena – Odeh expressed the material balance in terms of gas production, fluid
expansion, and water influx as;
Underground withdrawal = Gas Expansion + Water expansion/pore compaction + water influx =>
With water influx
Without influx
Gas produced
P/Z
GIIP
Fundamentals of Reservoir Engineering & Characterization 217
( ) ( )We
WC
fWCWgigigwpgp BWp
S
CSCGBBBGBWBG +∆
−+
+−=+1
F = G (Eg + Ef,w) + We BW Where; F= Underground fluid withdrawal= GPBg + WP BW Eg= Gas expansion=Bg - Bgi
Ef,w = Water and rock expansion = ( )
wi
fWiWgi S
CSCB
−+
1
Defining Drive Mechanism Dake (1994) presented an excellent discussion of the strengths and weaknesses of the
material balance equation as straight line. The Havlena-Odeh equation can be
expressed as following equation;
wfo
We
wfO EEBW
NEE
F
,, ++=
+ for oil reservoir
wfg
We
wfg EEBW
GEE
F
,, ++=
+ for gas reservoir
The classical approach to find out the type of drive mechanism, is to plot the right hand
side expression of oil reservoir/gas reservoir as shown in the above equation against oil
production/gas production. Dake suggested that typically two type of trend is observed.
• In a particular case all the points of wfO EE
F
,+ or
wfg EEF
,+ may lie on a
horizontal straight line, as shown in the plot as trend line A. Line A on the plot implies
that the reservoir can be classified as a volumetric reservoir, i.e. We=0. This defines
purely depletion type of reservoir. The energy for production from the reservoir purely
derives from the expansion of the rock, connate water and the oil. Furthermore, the
ordinate value of the plateau determines the initial oil in place N, or initial gas in place
GIIP.
Alternately the values of the ordinate term as shown in the plot below may rise as
illustrated by the trend B and C. Both the curve suggest of aquifer energy. Plot C
represents an active aquifer, whereas plot B represents weak reservoir. However, it
Fundamentals of Reservoir Engineering & Characterization 218
should be remembered that the trend is highly rate dependent. Higher rate than the
water influx into the reservoir may lead to dipping of the trend suggesting otherwise low
strength reservoir, whereas lower withdrawal rate than aquifer influx may suggest active
water drive reservoir. Hence, it may give wrong impression. The conclusion should be
made based on other reservoir parameters viz. aquifer volume, KV/Kh ratio, water
production trend etc.
Abnormally Pressured Gas Reservoir High pressured gas reservoirs, usually do not show the typical straight line relationship
between P/Z and gas production value. It is generally observed that they typically exhibit
two slopes. The second slope is steeper than the first slope. The initial/first slope is due
to gas expansion and significant pressure maintenance brought about by formation
compaction, water expansion. Hence, GIIP calculated from the first slope would be
erroneously very high. At approximately normal pressure gradient, the formation
compaction is essentially complete and the reservoir assumes the characteristics of
normal gas expansion reservoir. This accounts for the second slope.
Roach (1981) proposed a graphical technique for analyzing abnormally pressured gas
reservoirs. He put forward following equation for determining GIIP
αααα =(1/G) * ββββ -ER
Where;
NP or GP
A
B
C wfEEF
,+
Fundamentals of Reservoir Engineering & Characterization 219
( )
( )PP
ZP
ZP
i
i
i
−
−
=1/
α
and
( )
)(
/
Pp
ZPZP
i
i
i
−
=β
Fundamentals of Reservoir Engineering & Characterization 220
# % &
Technology to increase oil recovery from a porous formation beyond that obtained by
conventional means. Conventional oil recovery technologies produce an average of
about one-third of the original oil in place in a formation. Conventional technologies are
primary or secondary. Primary technologies rely on native energy, in the form of fluid
and rock compressibility and natural aquifers, to produce oil from the formation to wells.
Secondary technologies supplement the native energy to drive oil to producing wells by
injecting water or low-pressure gas at injection wells. The target of enhanced recovery
technologies is that large portion of oil that is not recovered by primary and secondary
means.
Many of the challenges encountered by secondary technologies are identical to those
encountered by enhanced recovery technologies. Those challenges include reducing
residual oil saturation, improving sweep efficiency, fitting the technology to the reservoir
heterogeneities, and minimizing up-front and operating costs.
Residual oil remains trapped in a porous rock after the rock has been swept with water,
gas, or any other recovery fluid. The residual oil saturation is the percentage of the pore
space occupied by the residual oil. The residual oil saturation depends on the pore size
distribution and connectivity, the interfacial tension between a recovery agent and the oil,
the relative wettability of the rock surfaces with respect to the recovery agent and the oil,
the viscosity of the fluids, and the rate at which the fluids are moving through the rock.
The sweep efficiency specifies that portion of a reservoir that is contacted by a recovery
fluid. Sweep efficiency increases with volume of injected fluid. It also depends on the
pattern of injection and production wells in a formation, on the mobility of the oil and the
recovery fluid, and on heterogeneities in the formation.
A wide variety of processes have been considered for enhancing oil recovery: thermal
processes, high-pressure gas processes, and chemical processes. Specifically, low
residual oil saturation can be obtained by selecting a recovery fluid that provides a very
low interfacial tension between the oil and the fluid. With very low interfacial tension, the
'
Fundamentals of Reservoir Engineering & Characterization 221
capillary number is large. And high sweep efficiency can be obtained by selecting a
recovery agent with low mobility or by increasing the mobility of the oil.
Different Phases in Field Development
There are broadly three phases in the development of a field. The phases are defined
as;
• Primary recovery phase
• Secondary recovery phase
• Tertiary recovery phase
Primary Recovery Phase
Primary oil recovery phase describes the production of hydrocarbons under the natural
driving mechanism present in the reservoir. The sources of natural reservoir energy are
fluid and rock expansion, solution gas drive, gravity drainage, and the influx of water
from aquifers. Based on the principal source of reservoir energy, the reservoirs are
classified as (1) Water drive, (2) solution gas drive, (3) fluid expansion, (4) gas-cap drive,
and (5) gravity drainage. These natural sources of energy displace oil towards the
producer without supplementary help from injected fluids such as water or gas.
Secondary Recovery Phase
Lack of sufficient natural drive in most reservoirs has led to the practice of
supplementing the natural reservoir energy by introducing some form of artificial drive,
the most basic method being the injection of gas or water.
Primary
Secondary
Tertiary
QO
Time
Fundamentals of Reservoir Engineering & Characterization 222
Water flooding, called secondary recovery because the process yields a second batch of
oil after a field is depleted by primary production
The practice of Water flooding apparently began accidentally as early as 1890, when
operators realized that water entering the productive formation was stimulating
production. The practice of Water flooding expanded rapidly after 1921. The earlier slow
growth of Water flooding was due to several factors. The oil demand was less and
impact of Water flooding on oil production was immense. However, after 1921 demand
of oil picked up and interest for Water flooding grew many folds. Gas injection developed
about the same time as the Water flooding and was a competing process in some
reservoirs.
Water or gas is pumped into the reservoir to produce more pressure on the oil, when
natural pressure is too low to bring the oil to the well.” Typical recoveries are 25-45%
after primary recovery (average 32%) of the total oil in place.
The four basic possibilities in such recovery are:
• Mining
• Squeezing
• Pushing
• Sucking
Mining involves removing the oil bearing rock from its position several hundreds or
thousands of feet below the surface of the earth. It is brought to the surface for
processing as an ore. This method proved to be uneconomical because of the
concentration of ore is low and the depths of most deposits make mining difficult.
Squeezing has to do with the pressing out the oil from the rock by force.
Pushing is the most successful secondary recovery. This is done by displacing the oil
from the rock with some other substance.
Sucking is a type of variation of pushing. The air in the atmosphere is used as a pusher.
The primary techniques are supplemented by the injection of water or gas in the
secondary recovery technique. They do not displace all of the oil. That which is trapped
by capillaries force in the pores is called residual oil.
Fundamentals of Reservoir Engineering & Characterization 223
WATERFLOODING:
In a water flood, water is injected in a well or pattern of wells to displace oil towards
producer. Initially, oil alone is produced as the part of the reservoir at the irreducible
water saturation is swept. When the leading edge of the capillary transition zone reaches
the producer breakthrough occurs (the first appearance of water in the produced
fluids).After breakthrough, both oil and water are produced and the watercut increases
progressively. Eventually the trailing edge of the capillary zone reaches the producer
and only water is produced. Because water is readily available and inexpensive, the
oldest secondary recovery method is water flooding, pumping water through injection
wells in to the reservoir.
The water is forced from injection wells through the rock pores, sweeping the oil ahead
of it towards production wells. This is practical for light to medium crude. Over time, the
percentage of water in produced fluids-the water cut-steadily increases. Some wells
remain economical with water cut as high as 99%.But at some point, the cost of
removing and disposing of water exceeds the income from oil production and secondary
recovery is then halted. While deciding suitability of a candidate reservoir for Water
flooding following reservoir characteristics should be considered;
Flood Pattern
The areal geometry of the reservoir will influence the location of wells and that will
essentially decide the flooding pattern (injection-production well arrangements) to be
Fundamentals of Reservoir Engineering & Characterization 224
used if the reservoir is to be produced through water-injection practices. The commonly
used flood patterns are given in the following figures;
Fundamentals of Reservoir Engineering & Characterization 225
The characteristics of the different flood patterns are given in the following table.
Pattern P/I
Regular
P/I
Inverted
d/a EA, %
Direct Line Drive 1 - 1 56
Staggered Line
Drive
1 - 1 76
4-Spot 2 1/2 0.866 -
5-Spot 1 1 1/2 72
7-Spot 1/2 2 0.866 -
9-Spot 1/3 3 1/2 80
P = number of Production wells
I = number of injection wells
d= distance from an injector to the line containing two producing wells
a = distance between wells in line in regular pattern
EA = Areal sweep efficieny at water break through for W/O = 1
• Mobility Ratio
Mobility ratio, which is the ratio of the displacing phase and the displaced phase, is an
important parameter for the selection of water flooding process. Mobility ratio less than 1
suggests that the water moves slower than the oil. This leads to piston type of
displacement leading to better sweep efficiency than cases where mobility ratio is
greater than 1. Low oil viscosity is preferred for water flooding. The reason is; at
abandonment areal sweep efficiency would be very high.
Mobility ratio = mobility of water in the water contacted portion / mobility of oil in
the oil bank
O
ro
W
rw
K
K
M
µ
µ=
• Recovery Efficiency
A simplistic model for estimating overall recovery involves factoring the recovery
efficiency into individual process efficiencies.
ER = EA * EV * ED * EM
Fundamentals of Reservoir Engineering & Characterization 226
Where;
ER = Overall recovery efficiency
EA = areal sweep efficiency
EV = Vertical sweep efficiency
ED = Displacement efficiency
EM = mobilization efficiency
Areal Sweep Efficiency
It is defined as the fractional area of the field that is invaded by an injected fluid. The
major factors determining areal sweep are fluid mobility, pattern type, areal
heterogeneity, extent of field development, and total volume of fluid injected.
Vertical Sweep Efficiency
It is defined as the fraction of the vertical section that is contacted by injected fluids and
is primarily a function of the vertical heterogeneity and the degree of vertical segregation.
Displacement Efficiency
It is the fraction of the mobile oil in the swept zone that has been displaced and is a
function of the volume injected, the fluid viscosities and the relative permeability curves
of the rock. Displacement efficiency will continually increase with increasing water
saturation in the reservoir. Buckley and Leverett developed a well established theory
called frontal displacement theory to determine the relation ship between the increase in
the average water saturation in the swept area as a function of cumulative water injected.
The theory will be discussed in the subsequent section.
Mobilization Efficiency
It is defined as the fraction of the oil in place at the start of a recovery process that
ultimately could be recovered by that process and is given as
oi
oi
oforpoi
oi
M
BS
BSBS
E/−
=
Soi = oil saturation at start of project
Boi = Oil formation factor at start of project
Sorp = residual oil to process
Bof = Oil formation volume factor at the end of process
Fundamentals of Reservoir Engineering & Characterization 227
Buckley and Leverett Theory of Frontal Displacement The Buckley and Leverett model was developed by application of the law of
conservation of mass to the flow of two fluids (Oil + Water) in one direction. The classic
theory consists of two equations;
Fractional Flow Equation
Frontal Advance Equation
The following three assumptions are made while deriving the frontal displacement
expression.
1. Incompressible flow
2. Fractional flow of water is a function of only water saturation
3. No mass transfer between phases takes place.
For a linear system mass flux rate in x direction both for oil and water can be written as;
( ) ( )φρρ OOoxO St
ux ∂
∂=∂∂− for oil
( ) ( )φρρ WWWxW St
ux ∂
∂=∂∂− for water
The above equation can also be written in volumetric form as;
( ) ( )φρρ OOoO St
Aqx ∂
∂=∂∂− for oil
( ) ( )φρρ WWWW St
Aqx ∂
∂=∂∂− for water
In the Buckley-Leverett model, water and oil are considered incompressible and thus O
and W are constant. Hence the above equation becomes;
t
SA
xq OO
∂∂
=∂
∂− φ
Fundamentals of Reservoir Engineering & Characterization 228
t
SA
xq WW
∂∂
=∂
∂− φ
The sum of above two equations gives;
( )WOWO SS
tA
xqq
+∂∂=
∂+∂
− φ)(
Since SO+SW =1.0
0)(
=∂+∂
−x
qq WO
or qO +qW =qt = constant Saturation qO and qW vary with distance x. However, because of oil and water are
assumed to be incompressible, the total volumetric flow rate at any time t is constant for
every position of x in the linear system.
The fractional flow of a phase is defined as the volume fraction of the phase that is
flowing at x, t.
For oil and water phases;
WO
OO qq
qf
+=
WO
WW qq
qf
+=
=>
t
SqA
xf W
t
W
∂∂
=∂
∂−
φ
The water saturation in a porous rock is a function of distance and time
Hence,
),( txSS WW =
or
Fundamentals of Reservoir Engineering & Characterization 229
dtt
Sdx
xS
Sx
W
t
WW
∂∂
+
∂∂
=∂
To know the value at particular instant of water saturation SW we can put dSw = 0
=>
t
W
x
W
S
xS
tS
dtdx
W
∂∂
∂∂
−=
The term WSdt
dx
is the velocity at which the saturation, SW , moves through the porous
media.
FW happens to be a function of saturation hence;
t
W
tW
W
tW
W
xS
Sf
Sf
∂∂
∂∂
=
∂∂
=>
tW
Wt
S Sf
Aq
dtdx
W
∂∂
=
φ
The above equation is called Buckley-Leverett equation, which states that in a linear
displacement process, each water saturation moves through the porous rock at a
velocity that can be computed from the derivative of the fractional flow with respect to
water saturation. For two-phase flow, the total flow rate qt is essentially equal to the total
injection rate, iW
WW SW
WW
S Sf
Ai
dtdx
∂∂
=
φ615.5
Where;
IW= water injection rate, bbl/day A=cross sectional area, ft2
Fundamentals of Reservoir Engineering & Characterization 230
The total distance specified water saturation will travel during a total time t,
( )W
W
S
WWS dt
dfA
tix
=
φ615.5
Where;
t= time, day
(x)Sw = distance from the injection for any given saturation SW, ft
Tertiary Recovery/EOR Phase
Tertiary recovery involves injecting other gases (such as carbon dioxide), or heat (steam
or hot water) to stimulate oil and gas flow to produce remaining fluids that were not
extracted during primary or secondary recovery phases. Typical recoveries are 5-20% of
OIP after primary and secondary recovery (average13%).The third type of recovery is
tertiary or enhanced. This “can sometimes be achieved if the viscosity of the oil is
lowered so that it flows more easily, either by heating the oil (by injecting steam, for
example) or by injecting chemicals into the reservoir.” The tertiary recovery is also a
supplementation of natural reservoir energy; however it is defined as that additional
recovery over and above what could be recovered from primary and secondary recovery
methods. Various types of tertiary or EOR recovery processes are given as follows;
EOR Processes
Thermal EOR
Processes
Chemical
EOR
Miscible EOR
Processes
Immiscible
EOR
Microbial EOR
Processes
• In-situ
combustion
• Air injection
• Steam
flooding
• Alkali-
Surfactant-
Polymer
• Polymer
• Hydrocarbon
miscible
• CO2 miscible
• N2 miscible
• Flue gas
• Hydrocarbon
immiscible
• CO2
immiscible
• N2
immiscible
• Flue gas
• Consortium
of Bacteria
used for
insitu
generation of
suphonates,
CO2,etc. for
profile
modification
Fundamentals of Reservoir Engineering & Characterization 231
As evident from the above other than Water flooding process, all other natural reservoir
energy supplementation processes have been considered as EOR process. Broadly,
EOR are essentially designed to recover oil, commonly described as residual oil.
Within a broad context, the applicability of the assortment of IOR and EOR technologies
depends by and large on two factors: the API gravity of the oils and the depth of the
reservoirs. In reality, the proper technical selection parameters are the oil viscosity and
the reservoir pressure. These, however, are related empirically to oil gravity and
reservoir depth, respectively.
Enhanced Oil Recovery
The ever increasing demand of hydrocarbon has led to vigorous E&P efforts in finding
new oil reserves. Newer oil reservoirs discovered are in general found to be less efficient
than their predecessors, making future recovery efficiency abysmally low. Conventional
exploitation methods in these reservoirs in general do not perform well. In physical
sense lot of oil remain In-situ. This has led to R&D efforts for improving recovery of oil,
producible beyond primary and secondary methods.
In the last decade or so, many techniques have been investigated in the laboratory to
improve the technology and methods for development and production of oil reservoirs
beyond primary and secondary recovery processes. These techniques for improving the
recovery beyond primary and secondary process have got its appellation Enhanced Oil
Recovery.
An Enhanced Oil Recovery (EOR) process involves supplementation of natural reservoir
energy externally to produce incremental oil that cannot be produced techno-
economically by conventional means.
Application of an EOR process in a particular reservoir involves four important steps- (i)
identification of suitable EOR process, (ii) laboratory studies, (iii) pilot testing, and (iv)
commercialization. Selection of the appropriate EOR process is the single most crucial
factor for success of any EOR project.
There is not a single process that can be considered a “cure-all” for recovering additional
oil from all types of reservoirs. Each process has its specific application, as they not only
depend upon reservoir rock and fluid properties but also on past production history.
Fundamentals of Reservoir Engineering & Characterization 232
Crude oil recovery by EOR processes is rather difficult and high risk operation, and the
likelihood of its success is influenced by a great many factors. However, technical and
economic criteria still dictates the selection of a process. The problem faced by the
reservoir engineer is to identify all the EOR process applicable to a candidate reservoir
or to check the suitability of a particular process in the light of all information available
about the reservoir under study. Prior information on reservoirs similar to the candidate
reservoir may also influence the choice of an EOR method. This leads to crucial need for
experts in this area of EOR process selection.
To understand the type of EOR processes and the range of reservoir rock and fluid
parameters suited for the process, the processes are discussed one by one.
In this process, steam is continuously introduced into an injection well. When steam is
injected into the reservoir, heat is transferred to the oil bearing formation, reservoir fluids
and some of the adjacent cap and base rock. The heat reduces the oil viscosity. This
increases the mobility of oil. As the steam loses heat energy it condenses to yield a
1.1a Steam Flooding Process
1.1 Thermal EOR Process
Fundamentals of Reservoir Engineering & Characterization 233
mixture of steam and hot water. Because of pressure gradient towards producing well,
an oil bank is formed ahead of steam zone. This enables the immobile oil to get
produced from the reservoir. In general steam reduces the oil saturation in the steam
zone to very low value (about 10±%. Some oil is also transported by steam distillation.
Technical Screening Guidelines
Crude Oil
Gravity : <36o API
Viscosity : > 20cP
Composition: Not critical but some light ends for steam distillation will help.
Reservoir
Type of formation : Sand or sand stone with high porosity and
Permeability preferred
Net Thickness : >15 feet
Average permeability : >10mD
Depth : 300-5000 ft
Temperature : Not Critical
Limitations :
• Oil saturation must be high, and the pay zone should be more than 15ft thick to
minimize the heat losses to adjacent formations.
• Lighter, less viscous crude oils can be steam flooded, but normally they are not if
reservoir responds to an ordinary water flood.
Fundamentals of Reservoir Engineering & Characterization 234
• Steam flooding is primarily applicable to viscous oil in massive, high permeability
sandstones or unconsolidated sands.
• Steam flooded reservoirs should be as shallow as possible as long as pressure
for sufficient injection rates can be maintained. This is to avoid excessive heat
losses in the well bore.
• Steam flooding is not normally used in carbonate reservoir.
• Cost per incremental barrel of oil is high.
• Low percentage of water-sensitive clays is desired for good injectivity.
• Adverse mobility ratio and channeling of steam may make this process
unattractive.
This process involves starting a fire in the reservoir and injecting air to sustain the
burning of some of the crude oil, usually in combination with water. A combustion front is
formed at which the injected air burns a small portion of the reservoir oil. The process
combustion can be achieved through low temperature oxidation and high temperature
oxidation. Low temperature oxidation is suited for light oil. Hot flue gas and steam
resulting from combustion and water vaporization displace the oil ahead of the
combustion front. Vaporization of the light ends and thermal cracking also occur. Ahead
of the combustion front, the vaporized light ends condense, providing some assistance
to displacement by solvent dilution of the virgin crude.
1.1 b In-situ Combustion Process
Fundamentals of Reservoir Engineering & Characterization 235
Technical Screening Guidelines
Crude Oil
Gravity : <48o API
Viscosity : < 100cP
Composition: Some asphaltic components to aid coke deposition.
Reservoir
Type of formation : Sand or sand stone with high porosity and
permeability preferred
Net Thickness : >15 feet
Average flow capacity : > 20 mD-ft
Depth : > 500 ft
Temperature : 150oF preferred
Limitations
• The process will not sustain if sufficient coke is not deposited.
• Excessive deposition of coke will lead to low advancement of combustion front
and eventually killing of the process in the presence of sufficient quantity of air.
• Oil saturation and porosity must be high to minimize heat loss to rock.
• Produced flue gases can pose environmental problem.
• Operation problems such as severe corrosion caused by low pH hot water,
increased sand production, pipe failures as a result of high temperature and
adverse mobility ratio makes this process complicated and difficult.
• Heterogeneous formation can result in poor sweep efficiency.
The process also called as micellar or micro emulsion flooding, consists of injecting a
slug that contains surfactant, co-surfactants, oil, water and other chemicals. The function
of the surfactant is to reduce oil/water interfacial tension, but it may also cause
interphase mass transfer of reservoir oil and water. Both the interphase mass transfer
and reduction of IFT increase recovery of oil. Surfactant slug is often followed by
polymer thickened water to improve sweep efficiency.
1.2 Chemical EOR Process 1.2a Alkali Surfactant Polymer Process
Fundamentals of Reservoir Engineering & Characterization 236
In surfactant flooding, surfactant molecules generally are injected along with water to
reduce the oil/ water interfacial tension (IFT), which reduces capillary forces that may
trap oil in rock pores . Normally in chemical flooding processes, inclusion of a viscosifier
(usually a water-soluble polymer) is required to provide an efficient sweep of the
expensive chemicals through the reservoir
Technical Screening Guidelines
Crude Oil
Gravity : > 20o API
Viscosity : < 100 cP
Composition: Light intermediates are desirable
Reservoir
Type of formation : Sand stone preferred
Net Thickness : > 10ft
Average permeability : > 10mD
Depth : 950 to 9000 ft( temperature!)
Temperature : <200o F
Limitations
• Adsorption of chemicals can be detrimental to the process.
• High temperature leads to degradation of chemicals.
• Best results are obtained when alkaline material reacts with crude oil. The oil
should have acid number more than 0.2 mg KOH/g of oil.
Fundamentals of Reservoir Engineering & Characterization 237
• High amounts of anhydrite, gypsum or clays are undesirable.
• Formation water chloride of < 20,000 ppm and divalent ions (Ca++ and Mg++)
<500 ppm are desirable.
• High heterogeneity may lead to the failure of the process.
• Vertical fractures may lead to gravity segregation.
The objective of polymer flooding is to provide better displacement and volumetric
sweep efficiencies during a water flood. They do not lower the residual oil saturation.
They improve recovery by increasing the viscosity of water, decreasing the mobility of
water, contacting a large volume of the reservoir and reducing the injected fluid mobility
to improve aerial and vertical sweep efficiencies. Because, polymer flooding inhibits
fingering, the oil displacement is more efficient in the early stages as compared to a
conventional water-flood.
Technical Screening Guidelines
Crude Oil
Gravity : > 15o API
Viscosity : <200 cP
Composition : Not critical
Reservoir
Type of formation : Sand stone preferred
Net Thickness : > 10ft
1.2b Polymer Process
Fundamentals of Reservoir Engineering & Characterization 238
Average permeability : > 10mD
Depth : < 9000 ft( temperature!)
Temperature : < 229o F
Limitations and Facts
• There are two types polymer synthetically produced polymers (Acrylamide
type) and bio-polymers (Xanthum gum).
• Factors which degrade polymers are salinity, temperature, time, shear rate
and presence of divalent ions.
• Bio-polymers suffer from bacterial degradation and cause well bore plugging.
• Polymers may be ineffective in a mature water flood because of low mobile
oil saturation.
• High adsorption on reservoir rock may kill the process.
• Optimum temperature is a key selection criterion for polymers. Clay increase
polymer adsorption.
• If oil viscosities are high, a higher polymer concentration is needed to achieve
the desired mobility control.
• Some heterogeneity is acceptable but for conventional polymer flooding,
reservoirs with extensive fractures should be avoided.
1.3. EOR by Gas Injection Processes:
Gas Injection is the second largest enhanced oil recovery process, next only to thermal
processes used in heavy oil fields. The residual oil saturations in gas swept zones have
been found to be quite low, however, the volumetric sweep of the flood has always been
a cause of concern. The mobility ratio, which controls the volumetric sweep, between the
injected gas and displaced oil bank in gas processes, is typically highly unfavorable due
to the relatively low viscosity of the injected phase. This difference makes mobility and
consequently flood profile control the biggest concerns for the successful application of
this process.
Hydrocarbon miscible flooding consists of injecting light hydrocarbons through the
reservoir to form a miscible flood. The process recovers crude oil by generating
1.3.a Hydrocarbon Miscible Flooding
Fundamentals of Reservoir Engineering & Characterization 239
miscibility (in the condensing and vaporizing gas drive), increasing oil volume by swelling
and decreasing the viscosity of oil. Three different methods of attaining miscibility are
used. One method uses about 5% PV slug of liquefied petroleum gas such as propane,
followed by natural gas or gas and water. A second method called enriched(condensing)
gas drive , consists of injecting 10 to 20% PV slug of natural gas that is enriched with
ethane through hexane(C2 to C6) , followed by lean gas and possible water. The
enriching components are transferred from gas to the oil. The third method, called High
pressure (Vaporizing) gas drive, consists of injecting lean gas at high pressure to
vaporize C2 to C6 components from the crude oil being displaced.