Wells in Niger Delta
Wells in Niger Delta
Wells in Niger Delta
net/publication/310447785
CITATIONS READS
0 1,384
1 author:
Ibianga Sukubo
Shell Petroleum Development Company, Port Harcourt, Nigeria
9 PUBLICATIONS 2 CITATIONS
SEE PROFILE
All content following this page was uploaded by Ibianga Sukubo on 17 November 2016.
By
(DE 97/0971)
August, 2003.
CERTIFICATION
This is to certify that every bit of this project work was planned, executed and
successfully completed in the Department of Petroleum Engineering, Rivers State
University of Science and Technology, Port Harcourt, during the 2001/2002
academic session.
This project work is unique in itself. It has not been presented in any form nor
for the award of any other degree in our Nigerian Universities and it fulfils the
requirement for the award of a Bachelor of Technology (B. Tech) degree in the
Department of Petroleum Engineering, Rivers State University of Science and
Technology, Port Harcourt.
............................................ .......................................
Sukubo Ibianga A. Date
(Researcher)
............................................ .......................................
Prof. D. T. Numbere Date
(Project Supervisor)
........................................... .......................................
Dr. S. A. Amadi Date
(Head of Department)
.............................................. .......................................
External Examiner Date
i
DEDICATION
This work is dedicated to everybody who has similar and/or better zeal for work
– all those persons who are creative, constructive and innovative by
demonstration and at heart. Together we can better the world.
ii
ACKNOWLEDGEMENT
My special thanks to God Almighty for being my only source of hope and
inspiration and for his loving kindness and favour unleashed upon me: like I have
never seen or heard.
My appreciation goes to my project supervisor, Professor D.T Numbere,
University of Missouri-Rolla, U.S.A for being a very strong drive to the success of
this work. Also, my thanks to Mrs. Areje Atike-Odibe, SPDC-East, Port Harcourt,
for providing me with the necessary data I needed.
To my parents, Chief and Mrs. Alu Mark Sukubo Ekwe; my brothers Tonye
Minaibi, Japu, Tekena, Owanaba and Abioton-a and my only sister, Mina; who
have been of immense assistance in innumerable ways. Also, thanks to
Tamunoibelemam Nangi Obu for prayers and moral support.
I must appreciate the effort of every other person and acquaintance towards the
successful completion of this project work. Particularly, my thanks go to Preye
Nelson Allagoa for setting his magic typing-fingers on this work. Thanks to all.
iii
LIST OF FIGURES
Figure 1.1 Typical Performance curves for Solution Gas Drive Reservoirs
Figure 1.2 Typical Performance curves for Water Drive Reservoirs
Figure 1.3 Typical Performance curves for Gas Cap Drive Reservoirs
Figure 1.4 A plot of q Vs Pwf
Figure 3.1 Real gas Pseudo-pressure versus Pressure
Figure 4.1 No Turbulence Drawdown Plot
Figure 4.2 Turbulence Drawdown Plot
Figure 4.3 No Turbulence AOFP Plot
Figure 4.4 Turbulence AOFP Plot
Figure 4.5 Vertical Well Turbulence Effect
Figure 4.6 Horizontal Well Turbulence Effect
iv
LIST OF TABLES
Table 3.1 Development of PVT Data
Table 3.2 PVT Data at 190°F
Table 3.3 Calculation procedure for Gas Pseudo-pressure, m(p)
Table 3.4 Pseudo-steady State IPR Calculations for Vertical Gas Well
Table 3.5 Pseudo-steady State IPR Calculations for Horizontal Gas Well
v
ABSTRACT
Horizontal wells in recent years have become the common feature of the
development plans in several hydrocarbon fields worldwide. But here in the Niger
Delta of Nigeria, due to the high permeability, horizontal wells are only employed
in the oil wells with no single horizontal gas well. Reason being that the high
permeability will cause excessive near wellbore turbulence in the horizontal wells.
However, no study has been conducted to validate this premise. This project
work, therefore, analysed and compared the horizontal and vertical gas well
inflow performance relationships (IPR) for Niger Delta gas fields where
turbulence was accounted for.
Due to the fact that, there are no available horizontal gas well data to work with
– the effective wellbore radius concept, assuming unstimulated vertical well, is
applied and an analogous horizontal gas well deliverability equation was
developed for comparative purposes. This equation is solvable by reducing it to
an equivalent form of the Forchheimer equation. The results, based on a fixed
1000ft of horizontal well length, showed the following:
1. A horizontal gas well will produce considerably more gas than a vertical gas
well with or without turbulence.
2. Turbulence substantially reduces the flow rate of horizontal gas wells. This
reduction increases with decreasing bottom-hole flowing wellbore pressure.
Thus, the decrease in AOFP due to turbulence is much higher in the
horizontal well as compared with the vertical well. This confirms the reason
for the lack of horizontal gas wells in the Niger Delta. The data used was
from a well in one of the gas fields in Niger Delta.
vi
TABLE OF CONTENTS
CERTIFICATION .................................................... i
DEDICATION ..................................................... ii
ACKNOWLEDGEMENT ............................................. iii
LIST OF FIGURES ................................................. iv
LIST OF TABLES ................................................... v
ABSTRACT ....................................................... vi
TABLE OF CONTENTS ............................................. vii
CHAPTER ONE ..................................................... 1
INTRODUCTION ................................................... 1
1.1 TYPES OF RESERVOIRS ......................................... 2
1.1.1 Solution Gas Drive .......................................... 2
1.1.2 Water Drive ............................................... 3
1.1.3 Gas Cap Expansion Drive...................................... 5
1.2 INFLOW PERFORMANCE RELATIONSHIPS ............................ 6
1.2.1 Productivity Index ........................................... 7
1.3 NIGER DELTA FORMATION ....................................... 8
1.4 STATEMENT OF THE PROBLEM .................................... 9
1.5 SCOPE ..................................................... 10
CHAPTER TWO ................................................... 11
LITERATURE REVIEW............................................... 11
2.1 Rate-Pressure Relation for Gas Wells ............................... 14
2.2 Horizontal Wells in Gas Reservoirs ................................. 19
2.3 Horizontal Well Application ...................................... 20
2.3.1 Tight Gas Reservoirs ......................................... 20
2.3.2 High Permeability Reservoirs ................................... 21
2.4 Horizontal Gas Well Equations .................................... 22
2.5 Effective Wellbore Radius, r’w Concept .............................. 23
CHAPTER THREE .................................................. 26
METHODOLOGY ................................................... 26
3.1 Calculation of the Real Gas Pseudo-pressure Function, m(p): ............. 29
3.2 Calculation Of The Gas IPR ...................................... 31
CHAPTER FOUR ................................................... 41
vii
ANALYSIS AND COMPARISON OF RESULTS ............................... 41
4.1 No Turbulence Drawdown ....................................... 41
4.2 Turbulence drawdown .......................................... 42
4.3 No Turbulence AOFP ........................................... 43
4.4 Turbulence AOFP ............................................. 44
4.5 Vertical Well Turbulence Effect ................................... 45
4.6 Horizontal Well Turbulence Effect ................................. 46
CHAPTER FIVE .................................................... 48
CONCLUSION AND RECOMMENDATION ................................. 48
REFERENCES ..................................................... 50
NOMENCLATURE .................................................. 52
APPENDIX "A" .................................................... 55
APPENDIX ”B" .................................................... 56
APPENDIX "C" .................................................... 58
APPENDIX "D" .................................................... 59
APPENDIX "E" .................................................... 61
viii
CHAPTER ONE
INTRODUCTION
The inflow performance of a well represents the ability of that well to give up
fluids. It could be analysed using typical plots of bottom-hole flowing pressure
against the flow rate of the well. The resulting curve gives the inflow
performance relationship (IPR) for the well. IPR curves may differ in shape,
depending upon the drive type. For instance, the flowing pressure versus rate
may be essentially a straight line (for water drive and /or pressure above
saturation pressure) or it may curve (for solution gas drive and flow below the
bubble point). The ability of a well to give up fluids depends to a great extent
upon the type of reservoir and drive mechanism and such variables as reservoir
pressure, permeability, etc.1
It is common practice to assume that inflow into a particular well with constant
conditions is directionally proportional to average reservoir pressure (PR).
Normally, this is true only for flowing pressures above bubble point (single phase
flow). Hence the productivity index of the well, which is the negative inverse of
the first derivative of the inflow performance, is constant. When the value of the
slope of the inflow performance is constant the well is said to have a single
productivity index (PI).
However, it is known that curvature exists in this plot for many wells. Hence, a
well cannot be said to have a linear PI (i.e. straight line), because the slope
varies with the variation in drawdown.
The IPR (Inflow Performance) and PI are not equivalents. The IPR is the
relationship between flowing bottom-hole pressure and flow rate. The PI is the
first differential of the IPR in the special case where the IPR is a straight line, or
is close enough to being straight that curvature can be neglected.
1
and again depends upon the type of reservoir. For a very active drive in which
the pressure remains above the bubble point, the PI remains constant, and for a
solution gas drive in which the flowing pressures are below the bubble point, the
PI changes rapidly.1 If an inflow performance curve is constructed at any time in
the life of a practical reservoir then the q against Pwf curve will probably be a
straight line for the water drive reservoir and curved for the solution gas drive
reservoir. The PI will probably be high for the water drive with very little drop in
average reservoir pressure at the high rates.
The type of reservoir will materially influence the production rate, hence the
shape of inflow performance curve.1 There are generally conceded to be three
basic reservoir types with possible two - and three - way combination of these
three. A brief discussion of each follows.
1. A constant volume. This means that there is no change in the initial size of
the reservoir. There is no water-encroachment for this particular type of
drive mechanism
2. There is two-phase flow at pressures below the bubble point. In other words,
gas coming out of solution flows along with the oil
3. The gas comes out of solution but does not move to form a gas cap. Gas
bubbles formed in the oil phase remain in the oil phase, resulting in
2
simultaneous flow of both oil and gas
4. Oil production is the result of the volumetric expansion of the solution gas
and the volumetric expulsion of oil
In the early production life of the reservoir, oil is replaced by gas on an equal
volume basis, but as pressure declines, a larger gas phase develops. More gas
expansion is then required per unit volume of oil produced because of the free-
flow ability of the gas phase.1 Now, creating excessive drawdown in a depletion
drive reservoir results in an increased gas phase in the reservoir - an increased
permeability to gas and decreased permeability to oil.
Figure (1.1) shows typical performance curves for the type of reservoir. Of
particular importance is the rapid decline in pressure, the rapid decline in PI and
the increase in the gas-oil ratio. Generally, wells in this type of reservoir can
expect to be low rate producers of oil (but high rate producers of gas) in their
latter life unless they have long pay intervals.1
Figure 1.1: Typical Performance curves for Solution Gas Drive Reservoir
3
hydraulic control. Some of the associated characteristics are 1:
1. The reservoir volume for oil does not remain constant. Water encroaches,
changing the initial volume of the reservoir
3. This reservoir type could also have a gas phase, resulting in a combination
water depletion drive
Figure (1.2) shows typical performance curves for this type of reservoir. In a
very active water drive, the pressure decline may be very small and in fact,
pressure may remain constant. Of great importance is the trend of the PI to
remain constant over the life of the well. In return the GOR remains constant.
In a water drive reservoir, the PIs of individual wells are normally more reliable
than those of a depletion drive reservoir. Water encroachment may be such that
there is very little loss in bottom-hole pressure. It is generally conceded that PI
information may be extrapolated linearly for drawdowns necessary to give the
desired production. Probably, in most cases, the pressure remains above the
bubble point. However, there are instances where an increase in water-cut
causes a decrease in PI.
4
1.1.3 Gas Cap Expansion Drive
This type of reservoir drive mechanism may also be referred to as segregation
drive. The reservoir is in a state of segregation – an oil zone overlain by a gas
cap. The drive may be further classified as to whether or not gas coming out of
solution in the reservoir flows to the gas cap. A segregation drive with counter
flow will have gas coming out of solution and moving to the gas cap. As
production proceeds, the gas cap expands and moves down, resulting in gas cap
expansion drive.
Figure 1.3: Typical Performance curve for Gas Cap Drive Reservoir
Summarily, there are many reservoirs having combination drive mechanisms and
5
their performance may differ considerably from the typical curves given for water
drive, solution gas drive, or gas-cap expansion drive mechanisms. An effort
should be made to identify the reservoir drive mechanism in order to permit a
better determination of the ability of the well to give up fluids. There is no
substitute for good data taken during the early life of the reservoir in predicting
future performance.1
All well deliverability equations relate the well production rate and the driving
force in the reservoir, that is, the pressure difference between the initial, outer
boundary or average reservoir pressure and the flowing bottom-hole pressure.
If the bottom-hole pressure is given, the production rate can be obtained readily.
However the bottom-hole pressure is a function of the wellhead pressure which,
in turn, depends upon production engineering decisions, separator or pipelines
etc. Therefore, what a well will actually produce must be the combination of
what the reservoir can deliver and what the imposed wellbore hydraulics can
allow.2
It is then useful to present the well production rate as a function of the bottom-
hole pressure. This type of presentation is known as an "Inflow Performance
Relationship" (IPR) curve. Usually the bottom-hole pressure, Pwf is graphed on
the ordinate and the p reduction rate, q is on the abscissa.
6
1.2.1 Productivity Index
The commonly used term, PI (Productivity Index), represents one point on the
inflow performance curve. The PI is defined as q/∆p in bpd per psi pressure drop
(for oil wells) from static reservoir pressure to flowing bottom-hole pressure.1 It
is common practice to measure one or two PIs in the early life of a well, then use
that same PI to estimate drawdowns necessary for greater production rates as
well as assuming that the same PI exists later in the life of the well. This is
probably fairly safe for wells in a water drive field where the flowing pressure is
above the bubble point. But this may be in error for wells in a solution gas drive
reservoir or at flow below the bubble point.1
The Productivity Index not only changes with time or cumulative production but
is also subject to change with increased drawdown at any one specific time in
the life of the well. If we measure several PIs in a well during a specific time
interval, a relationship will be obtained between rate and flowing pressure which
normally is not linear for a solution gas drive field. This phenomenon may be due
to one or more of the following factors: 1
2. Change from laminar to turbulent flow in some flow capillaries near the
wellbore at increased producing rates.
3. Exceeding critical flow rates through pores at formation face in the wellbore.
These pores act as orifices and when the critical rate is exceeded, increased
drawdowns have a diminished effect on increasing rates.
A plot of q vs Pwf is called Inflow Performance and was first used by Gilbert3 in
describing well performance. A typical plot is noted in figure (1.4) and differs
depending upon the type of reservoir.
7
Figure 1.4: Plots of q vs Pwf
The Niger Delta Basin is a prolific oil (and gas) province with an area covering
some 80,000 sq. km. The bulk of the reserves are onshore although sizeable oil
reserves have also been discovered offshore.
The basin had its origin closely associated with the rift separation of Africa and
South America in early Cretaceous times. The formation of the Niger Delta,
however, can be traced back to the start of Oligocene while the present cone-
shaped front developed during the Miocene. The thick wedge of clastic
sediments beneath the Niger Delta is therefore of early Tertiary to recent age.
Gravity and Magnetic data indicate the sediment fill to be over 40,000 ft thick,
with the thickest part, west of the Niger River. The floor of the basin is
characterised by basement block faulting.4
The tertiary sequence of the Niger Delta is subdivided vertically into three broad
lithofacies units (in descending order) as follows4:
8
2. Agbada Fo rmation (or paralic sequence), consisting of inter-bedded
sandstones, siltstones and shales. This sequence is associated with
synsedimentary growth faulting and contains the bulk of the known oil
accumulations in the Niger Delta. Estimated thickness: about 4,000 - 10,000
ft
In the Paralic and Continental intervals, there are a few anomalous thick clay-
filled bodies interpreted as submarine channel features (e.g. the Afam, the Qua
Iboe and the Opuama).4These features are located on the NW and the SE fringes
of the Niger Delta. One of the most important and conspicuous geological
features of the Niger Delta is its growth fault pattern. Growth faulting is the main
structure and trap forming element of the Delta and it is as a result of density
contrasts and gravity sliding caused by rapid deltaic sedimentation.
It should be of note however that the boundaries between each unit or formation
are not always sharp and transitional situations are common. Typically, a well in
the Niger Delta will penetrate a section of Continental sands on top, grading
downwards through paralic sands and shales into marine shale sediments which,
at depth, are in a state of undercompaction. Only a small number of wells drilled
on the Delta have penetrated into the marine Akata shales.
The statement of problem as well as the purpose of this study is to compare the
IPRs of both vertical and horizontal gas wells using similar data from the Niger
Delta field. Hence, this project, on conclusion, is expected to serve as a simple
reference to vertical and horizontal gas wells IPRs. This would help in enhancing
our understanding of horizontal gas well technology and reasons for its
9
application or disapproval in the Niger Delta gas formations. Also, this would give
a positive contribution into the effects of production optimization of our fields.
1.5 SCOPE
3. Determination of the need for pressure maintenance and plan for future
artificial lift requirements.1
Therefore the obvious reason for the IPR is to determine what the production
rate will be if a certain backpressure is exerted on the wellhead.
10
CHAPTER TWO
LITERATURE REVIEW
Perhaps the simplest and most widely used IPR equation is the straight line IPR
which states that rate is directly proportional to pressure drawdown in the
reservoir. The constant of proportionality, J is called the productivity index;
defined as the ratio of rate to pressure drawdown in the reservoir. This can be
directly obtained from Darcy's radial steady state flow equation as:
11
𝐾𝑜 ℎ(𝑃𝑅 − 𝑃𝑤𝑓 )
𝑞𝑜 = 7.08 𝑥 10−3 (2.1)
𝑟 1
𝜇𝑜 𝐵𝑜 [𝐼𝑛 ( 𝑒⁄𝑟𝑤 ) − ]
2
𝑞𝑜 7.08 𝑥 10−3 𝐾𝑜 ℎ
𝐽= = (2.2)
(𝑃𝑅 − 𝑃𝑤𝑓 ) 𝑟 1
𝜇𝑜 𝐵𝑜 [𝐼𝑛 ( 𝑒⁄𝑟𝑤 ) − ]
2
A limitation on the straight line IPR is the assumption that oil is undersaturated
and so is only slightly compressible. Obviously, this condition does not apply to
gases or saturated oil wells; which evolve considerable amount of gas, both of
which are highly compressible. The effect of compressible gas at two-phase flow
on IPR was observed in the 1920's and 1930's during field testing. Instead of
linear rate increase with pressure drawdown, it was observed that larger-than-
linear pressure drops were required to increase the rate. The rate-pressure
relationship shows curvature pronounced at higher rates. In terms of
Productivity Index, J decreases with increase in drawdown. Note that, J is not
represented by the tangent to the rate-pressure curve but defined as in equation
(2.2).
Gilbert3 called the flow from reservoir into the well "Inflow Performance" and a
plot of producing rate versus bottom-hole flowing pressure is called "Inflow
Performance relationship" or IPR.
Methods for predicting IPR's for both the present, real time and for future times
shall be mentioned. Oil and gas reservoir performance shall be highlighted.
12
𝑞𝑜 𝑃𝑤𝑓 𝑃𝑤𝑓 2
⁄𝑞𝑜,𝑚𝑎𝑥 = 1 − 0.2 𝑃̅ − 0.8 ( 𝑃̅ ) (2.4)
𝑅 𝑅
Although, the method was proposed for saturated dissolved-gas drive reservoir
mainly, it has been observed to apply for any reservoir in which gas saturation
increases as pressure decreases. The reason why equation (2.4) is only valid for
two-phase flow in the reservoir is because, Vogel neglected gravity segregation
in his assumptions and simulation. In fact, Vogel's method did not account for
the effects of a non-zero skin factor.
Fetkovich10 proposed a method for predicting Inflow Performance for oil wells
using the Multipoint testing procedure as for gas wells. The procedure involves
analysing flow-after-flow and Isochronal tests conducted in reservoirs with
permeabilities ranging from 6md to greater than 1000md. Pressure conditions in
the reservoir range from highly undersaturated to saturated at initial pressure
and to partially depleted fields with gas saturation above the critical. He noted
that oil wells producing below the bubble point behave much like gas wells.
13
the surface equipment.
Rawlins and Schellhardt11, engineers from the US bureau of Mines, proposed the
classic backpressure equation relating gas rate, to flowing pressure.
𝑛
𝑞𝑔 = 𝐶(𝑃𝑅 2 − 𝑃𝑤𝑓 2 ) (2.5)
Equation (2.5), though generally accepted for gas wells, has been shown to be
an accurate relation even for saturated oil wells.
Since, Darcy's law breaks down at high flow velocity, there had been several
proposed models to either replace or modify the Darcy's law for high-velocity
flow. The most accepted model was that proposed by Forchheimer12:
𝜕𝑃⁄ = 𝑎𝑣 + 𝑏𝑣 2 (2.6)
𝜕𝑟
where 𝑎 and 𝑏 are constants and 𝑣 = 𝑞/𝐴 is the apparent fluid velocity. Later
14
equation (2.6) was expressed in terms of fluid and rock properties by Green and
Duwez13 and Cornell and Katz14.
where;
𝜇 = gas viscosity
𝜌 = gas density
𝐾 = formation permeability
𝛽 = high veloicty coeffiient; which is a property of the formation rock that accounts
for the deviation from Darcy ′ s law
Note that the deviation from Darcy's flow is more pronounced in gas wells than
in oil wells.
In relation to the radial flow equation, there are two major notable differences
between gas and undersaturated oil flow:
The integral represents the area under the curve of 𝑃/𝜇8 𝑍 versus pressure,
popularly known as the real gas pseudo-pressure, 𝑃/𝜇8 𝑍 plot; which is linear at
15
low pressures of less than 2000psia, and intercepts at the origin5. At high
pressures of greater than 3000psia, the pressure function is nearly constant,
showing some decrease in increasing pressures. And at intermediate pressures –
between 2000psia and 3000psia – the pressure function shows distinct
curvature.
This implies that, 1/𝜇𝑔 𝑍 is constant. Hence, 𝜇g and 𝑍 can be evaluated at any
pressures. But conventionally, they are evaluated at PR. Thus equation (2.8) now
becomes;
0.000703𝐾ℎ 2
Where the skin, 𝑆 = ′
(𝑃𝑤𝑓 − 𝑃𝑤𝑓 2 )
𝑞𝑔 𝑇𝜇𝑔 𝑍
′
𝑃𝑤𝑓 = Wellbore flowing pressure for an ideal well
Using the Forchheimer modification of Darcy's law for high velocity flow;
16
𝑃𝑅
𝑃 𝑃
2∫ 𝑑𝑝 = 2 (𝑃 − 𝑃𝑤𝑓 ) (2.12)
𝑃𝑤𝑓 𝜇𝑔 𝑍 𝜇𝑔 𝑍 𝑅
Note; 𝑃/𝜇8 𝑍 is evaluated at any pressure between 𝑃𝑤𝑓 and 𝑃𝑅 . The resulting
IPR equation for gas wells producing at high flowing and static pressure is:
such that;
𝑃
𝐴(𝑃) = ∫0 (𝑃⁄𝜇 𝑍) 𝑑𝑝 (2.14)
𝑔
𝑃
A special name is given to 2 ∫0 (P/𝜇𝑔 𝑍)𝑑𝑝 = 2𝐴(𝑃).This special name is the
Al-Hussainy and Ramey15 derived the real gas flow equation by solving the
partial differential equation of real gas flow, which comes from the combination
of the continuity equation and Darcy's law for radial flow. They used the
transformation;
P
P
𝑚 (𝑃 ) = 2 ∫ ( ) dp (2.15)
0 𝜇𝑔 𝑍
17
which is analogous to Kirchhoff’s transformation, to solve the equation and called
it the gas pseudo-pressure6; where 𝑚(𝑃) is the analogue to pressure or
pressure-squared in equation (2.9). The differential pseudo-pressure, ∆𝑚(𝑃) =
𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 ), represents the driving force or potential, moving gas toward
the well.
P𝑅
∆𝑚(𝑃) = 2 ∫ (𝑃⁄𝜇 𝑍) dp
P𝑤𝑓 𝑔
P𝑅 P𝑤𝑓
= 2∫ (𝑃⁄𝜇 𝑍) dp − 2 ∫ (𝑃⁄𝜇 𝑍) dp
0 𝑔 0 𝑔
Equation (2.8), and including the effect of high-velocity flow, now becomes
0.000703𝐾ℎ[𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 )]
𝑞𝑔 = (2.17)
𝑟
𝑇 [𝐼𝑛 ( 𝑒⁄𝑟𝑤 ) − 0.75 + 𝑆 + 𝐷𝑞𝑔 ]
Equation (2.17) may be considered the most general gas IPR for stabilized flow.
The relation for skin is then given as:
0.000703𝐾ℎ ′
𝑆= [𝑚(𝑃𝑤𝑓 ) − 𝑚(𝑃𝑤𝑓 )] (2.18)
𝑞𝑔 𝑇
′ ′
where 𝑚(𝑃𝑤𝑓 ) corresponds to the ideal wellbore flowing pressure 𝑃𝑤𝑓 and
18
2. The magnitude of 𝑚(𝑃) is much larger than pressure (usually on the order of
100 times pressure-squared, or 𝑚(𝑃) ≈ 100𝑃2 .
Note that, the constant 𝐷 in the rate-dependent skin factor is related to 𝑏 in the
Forcehheimer model and it is related to the rock properties;
𝛾𝑔 𝐾ℎ
𝐷 = 2.222 𝑥 10−18 𝛽 (2.19)
𝜇𝑔 𝑟𝑤 ℎ𝑝2
𝛾𝑔 ℎ
𝐷 = 3.75 𝑥 10−8 2 , (1⁄𝑆𝐶𝐹𝐷 ) (2.21)
𝜇𝑔 𝑟𝑤 ℎ𝑝
Horizontal wells can be suitable for gas reservoirs. They are applicable in low
permeability reservoirs as well as in high-permeability gas reservoirs (Joshi,
19
S.D).16
20
start off pseudo-steady state flow can be very large), vertical wells can be drilled
at close spacing to efficiently drain the reservoir. But this would require a large
number of vertical wells.
Note that, a fracture job that establishes a link with a natural open fault may
result in a high productivity well, since an open fault itself represents a high-
conductivity fracture.16
In gas wells where near-wellbore turbulence is very high (i.e. high permeability
gas wells) horizontal wells are also applicable. The near wellbore turbulence is
inversely proportional to the wells perforated interval. By drilling a horizontal
21
well, production length can be increased and so, decrease the near-wellbore
turbulence and enhance well productivity.
𝐾𝐻 ℎ(𝑃𝑅 2 − 𝑃𝑤𝑓 2 )
𝑞𝑔 =
𝐿 2
𝑎 + √𝑎2 − ( ) 𝐼𝑎𝑛𝑖 ℎ 𝐼𝑎𝑛𝑖 ℎ
1424𝜇𝑔 𝑍𝑇 𝐼𝑛 2 + [𝐼𝑛 { } + 𝐷𝑞𝑔 ]
𝐿 𝐿 𝑟𝑤 (𝐼𝑎𝑛𝑖 + 1)
2
[ { } ]
(2.22)
𝐾𝐻
𝐼𝑎𝑛𝑖 = √ (2.23)
𝐾𝑉
Equation (2.22) is the non-Darcy steady state equation for horizontal gas wells.
The non-Darcy pseudo-steady state equation for horizontal gas well is2:
𝐾𝐻 ℎ(𝑃𝑅 2 − 𝑃𝑤𝑓 2 )
𝑞𝑔 =
2
𝐿
𝑎 + √𝑎2 − ( ) 𝐼 ℎ 𝐼 ℎ
1424𝜇𝑔 𝑍𝑇 𝐼𝑛 2 + 𝑎𝑛𝑖 [𝐼𝑛 { ( 𝑎𝑛𝑖 } − 0.75 + 𝐷𝑞𝑔 ]
𝐿 𝐿 𝑟𝑤 𝐼𝑎𝑛𝑖 + 1)
2
[ { } ]
(2.25)
Using real gas pseudo-pressures, 𝑚(𝑃) equations (2.22) and (2.25) become:
22
𝐾𝐻 ℎ[(𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 )]
𝑞𝑔 =
𝐿 2
𝑎 + √𝑎 2 − ( ) 𝐼𝑎𝑛𝑖 ℎ 𝐼𝑎𝑛𝑖 ℎ
1424𝑇 𝐼𝑛 2 + [𝐼𝑛 { } + 𝐷𝑞𝑔 ]
𝐿 𝐿 𝑟𝑤 (𝐼𝑎𝑛𝑖 + 1)
2
[ { } ]
(2.26)
𝐾𝐻 ℎ[(𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 )]
𝑞𝑔 =
𝐿 2
𝑎 + √𝑎 2 − ( ) 𝐼 ℎ 𝐼 ℎ
1424𝑇 𝐼𝑛 2 + 𝑎𝑛𝑖 [𝐼𝑛 { ( 𝑎𝑛𝑖 } − 0.75 + 𝐷𝑞𝑔 ]
𝐿 𝐿 𝑟𝑤 𝐼𝑎𝑛𝑖 + 1)
2
[ { } ]
(2.27)
The effective wellbore radius concept is used to represent the well which is
producing at a rate different than that expected from calculations based upon a
drilled wellbore radius. Effective wellbore radius is the theoretical well radius
To calculate the required vertical wellbore diameter to produce gas at the same
rate as that of a horizontal well, equal drainage volumes, 𝑟𝑒ℎ = 𝑟𝑒𝑣 and equal
23
2𝜋𝐾ℎ ℎ 2𝜋𝐾ℎ ℎ
⁄𝜇 𝐵 ⁄𝜇 𝐵
𝑔 𝑔 𝑔 𝑔
=
2
𝑟 𝐿
𝑇 [𝐼𝑛 ( 𝑒⁄ )] 𝑎 + √𝑎2 − ( ) 𝐼 ℎ 𝐼𝑎𝑛𝑖 ℎ
𝑟′𝑤 𝑇 𝐼𝑛 2 + 𝑎𝑛𝑖 {𝐼𝑛 }
[ ]𝑣 𝐿 𝐿 𝑟𝑤 (𝐼𝑎𝑛𝑖 + 1)
[ [ { 2 } ]]ℎ
(2.28)
𝑟𝑒ℎ (𝐿⁄2)
𝑟′𝑤 = 𝐼𝑎𝑛𝑖 ℎ (2.29)
2
𝐿 𝐼𝑎𝑛𝑖 ℎ 𝐿
𝑎 [1 + √1 − ( ) ][ ]
2𝑎 𝑟𝑤 (𝐼𝑎𝑛𝑖 + 1)
𝑟𝑒ℎ (𝐿⁄2)
𝑟′𝑤 = ℎ (2.30)
2
𝐿 ℎ 𝐿
𝑎 [1 + √1 − ( ) ][ ]
2𝑎 2𝑟𝑤
24
general comparison. Hence, the effective wellbore radius concept as employed in
general comparison gives a fair estimation of productivity improvements with
horizontal wells.16
25
CHAPTER THREE
METHODOLOGY
From the single point data and gas composition given by the industry, the PVT
data for the wet-gas reservoir is developed using fluid properties correlations.
Hence, the molecular weight and specific gravity of the gas condensate is
obtained. Then, the pseudo-reduced pressure and pseudo-reduced temperature
relations, together with the Beggs and Brill's correlation (1973) for Z-factor
which, according to Golan M. et al5 serves as one of the best-fit correlation for
the Standing-Katz Z-factor graphical correlation (1942) and the Lee et al
correlation for gas viscosity calculation are used to write a simple programme in
Excel to calculate the PVT properties of the gas condensate (specifically the gas
compressibility factor, Z and the gas viscosity, 𝜇𝑔 ) at various working pressures
but constant temperature.
The trapezoidal rule of integration and its tabular approach is used to calculate
the real gas pseudo-pressures at the various wellbore flowing pressures through
Excel programme; during which the Z-factor and gas viscosities at the various
pressures are used as input variables, and the result displayed in tabular form.
Then, a plot of real gas pseudo-pressure versus pressure is made for the given
gas field in the Niger Delta. This may apply to all wells in the gas field throughout
the production life.
Similarly, the IPR is calculated for Darcy and non-Darcy flows, using the vertical
and horizontal gas well deliverability equations, at the various pseudo-pressures.
26
For ease in computation, a programme in Excel is written and result shown in
both tabular and graphical forms, using the given data by the industry in
Appendix ‘A’, and the gas deliverability equations at the various pseudo-
pressures.
(1−𝐴)
𝑍=𝐴+ + 𝐶𝑃𝑝𝑟 𝐷 (3.1)
𝑒𝑥𝑝 𝐵
27
Where:
0.066 0.32
𝐵 = (0.62 − 0.23𝑇𝑝𝑟 )𝑃𝑝𝑟 + [ − 0.037] 𝑃𝑝𝑟 2 + 𝑃𝑝𝑟 6
(𝑇𝑝𝑟 − 0.86) 109 (𝑇𝑝𝑟 − 1)
𝐶 = (0.132 − 0.32𝑙𝑜𝑔𝑇𝑝𝑟 )
2)
𝐷 = 10(0.3106 – 0.49𝑇𝑝𝑟 + 0.1824𝑇𝑝𝑟
𝑇 𝑃
𝑇𝑝𝑟 = ; 𝑃𝑝𝑟 =
𝑇𝑝𝑐 𝑃𝑝𝑐
Where:
𝑇𝑝𝑟 = pseudo-reduced temperature of the gas mixture,
𝑃𝑝𝑟 = pseudo-reduced pressure of the gas mixture,
T = working temperature of the gas mixture, and
P = working pressure of the gas mixture.
Where:
(9.4 + 0.02𝑀)𝑇1.5
𝐾=
209 + 19𝑀 + 𝑇
986
𝑋 = 3.5 + + 0.01𝑀
𝑇
𝑌 = 2.4 − 0.2𝑋
𝛾𝑔 𝑃
𝜌𝑔 = 0.0433 , 𝑔/𝑐𝑐 (T is in °R)
𝑍𝑇
Table (3.2) shows the generated PVT data at a temperature of 190ºF across the
possible pressure ranges.
28
Table 3.2: PVT Data at 190°F
𝒁−
𝑷𝒊 , 𝒑𝒔𝒊𝒂 𝝁𝒈 , 𝒄𝒑
𝒇𝒂𝒄𝒕𝒐𝒓
0 1.0000 0.0122
14.7 0.9984 0.0122
400 0.9498 0.0126
800 0.8972 0.0135
1200 0.8510 0.0147
1600 0.8166 0.0163
2000 0.7967 0.0183
2250 0.7917 0.0197
2500 0.7921 0.0212
2750 0.7975 0.0228
3000 0.8072 0.0243
3250 0.8206 0.0259
3500 0.8370 0.0274
3750 0.8558 0.0289
4000 0.8765 0.0303
4250 0.8986 0.0317
4500 0.9218 0.0331
4750 0.9458 0.0344
5000 0.9703 0.0357
5150 0.9851 0.0364
5167 0.9868 0.0365
29
7. 2(P/𝜇𝑔 𝑍)𝑎𝑣 ∆P: incremental pseudo-pressure (psi2/cp)
Note that, each step corresponds to a column in table (3.3) below which
tabulates the eight quantities.
a
Column [5] is the arithmetic average of the pressure function, P/μg Z in the
pressure interval, ∆P.
b
1.78E+04 is the notation for 1.78 * 104
30
Below is a plot of the calculated pseudo-pressures against the working pressures
in the above table.
For vertical wells, the IPR for Darcy pseudo-steady state flow is2:
31
𝛽 = 2.73 𝑥 1010 𝐾𝑎−1.1045 (2.20)
Where; 𝛾𝑔 is the gas specific gravity; 𝜇𝑔 is the gas viscosity; 𝑟𝑤 is the wellbore
radius; ℎ𝑝 is the perforation thickness and 𝛽 is the high velocity coefficient.
𝛾𝑔 ℎ
𝐷 = 3.75 𝑥 10−8 (2.21)
𝜇𝑔 𝑟𝑤 ℎ𝑝 2
𝑟
𝑇 [𝐼𝑛 (𝑟𝑒 ) − 0.75 + 𝑆] 𝑇𝐷
𝑤
[𝑚(𝑃𝑅 ) − 𝑚 (𝑃𝑤𝑓 )] = 𝑞𝑔 + 𝑞𝑔 2
7.03 𝑥 10−4 𝐾ℎ −4
7.03 𝑥 10 𝐾ℎ
(3.5)
𝑇𝐷
𝐵=
7.03 𝑥 10−4 𝐾ℎ
32
Example 3.1
Using the data in Appendix ‘A’, calculate the vertical gas well IPR when the gas
viscosity, 𝜇𝑔 = 0.0122cp; the gas specific gravity, 𝛾𝑔 = 0.8925, 𝑚(𝑃𝑅 ) =
1.25E+09psi2/cp and 𝑚(𝑃𝑤𝑓 ) = 1.78E+04psi2/cp [Hint: ℎ=ℎ𝑝 ]
Solution
First, the rate-dependent non-Darcy skin is determined where there is no
damage or stimulation near the wellbore , 𝐾𝑎 = 𝐾𝑔 . Then,
𝛾𝑔 ℎ
𝐷 = 3.75 𝑥 10−8 2 = 1.512 𝑥 10−8 , (1⁄𝑆𝐶𝐹𝐷 )
𝜇𝑔 𝑟𝑤 ℎ𝑝
𝐾ℎ = (420mD)(395ft) = 165,900mD-ft
𝑟
[𝐼𝑛 ( 𝑒⁄𝑟𝑤 ) − 0.75 + 𝑆]= 7.2076
Hence, the Darcy pseudo-steady state gas deliverability equation (3.3) may be
solved by substituting values from the deductions as:
7.03 𝑥 10−4 (165,900)[1.25 𝑥 109 − 1.78 𝑥 104 ]
𝑞𝑔 = , 𝑀𝑠𝑐𝑓/𝑑
(650)(7.2076)
= 31,117.27Mscf/d ≈ 31.117MMscf/d
(650)(7.2076)
𝐴= ≈ 40.17
7.03 𝑥 10−4 (165,900)
(650)(1.512 𝑥 10−8 )
𝐵= −4
≈ 8.43 𝑥 10−8
7.03 𝑥 10 (165,900)
33
Hence,
−8 9 0.5
[(40.17)2 + 4(8.43 𝑥 10 )(1.25 𝑥 10 )] − 40.17
𝑞𝑔 = −8
, 𝑀𝑠𝑐𝑓/𝑑
2(8.43 𝑥 10 )
= 29,314.37Mscf/d ≈ 29.314MMscf/d
For Horizontal wells, the IPR for Darcy pseudo-steady state flow is2:
𝐾𝐻 ℎ[(𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 )]
𝑞𝑔 =
𝐿 2
𝑎 + √𝑎 2 − ( ) 𝐼𝑎𝑛𝑖 ℎ 𝐼𝑎𝑛𝑖 ℎ
1424𝑇 𝐼𝑛 2 + [𝐼𝑛 { } − 0.75]
𝐿 𝐿 𝑟𝑤 (𝐼𝑎𝑛𝑖 + 1)
2
[ { } ]
(3.8)
𝐾𝐻 ℎ[(𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 )]
𝑞𝑔 =
𝐿 2
𝑎 + √𝑎 2 − ( ) 𝐼𝑎𝑛𝑖 ℎ 𝐼𝑎𝑛𝑖 ℎ
1424𝑇 𝐼𝑛 2 + [𝐼𝑛 { } − 0.75 + 𝐷𝑞𝑔 ]
𝐿 𝐿 𝑟𝑤 (𝐼𝑎𝑛𝑖 + 1)
2
[ { } ]
(2.27)
2, 16, 17
Where :
0.5 0.5
2𝑟𝑒𝐻 4 𝐿
𝑎 = 0.5𝐿 {0.5 + [0.25 + ( ) ] } 𝑓𝑜𝑟 < 0.9𝑟𝑒𝐻 (2.24)
𝐿 2
According to Economides2, 𝐼𝑎𝑛𝑖 = √𝐾𝐻 /𝐾𝑉 and 𝐼𝑎𝑛𝑖 ℎ/𝐿 is the reservoir anisotropic
scaled aspect ratio.
34
the reservoir, 𝐼𝑎𝑛𝑖 becomes equal to unity. Hence, the ratio ℎ/𝐿 may be referred
to as the isotropic aspect ratio, in the course of this work.
Using the effective wellbore radius, 𝑟′𝑤 concept according to Joshi, S.D13; for
steady state Darcy flow:
𝑟𝑒ℎ (𝐿⁄2)
𝑟′𝑤 = ℎ (2.30)
2
𝐿 ℎ 𝐿
𝑎 [1 + √1 − ( ) ][ ]
2𝑎 2𝑟𝑤
It can be proved (as shown in Appendix ‘E’) that the effective wellbore radius for
Darcy and non-Darcy pseudo-steady state flows are:
0.472𝑟𝑒ℎ (𝐿⁄2)
𝑟′𝑤𝑝 = ℎ
(3.9)
2
𝐿 0.472ℎ 𝐿
𝑎 [1 + √1 − ( ) ][ ]
2𝑎 2𝑟𝑤
and
0.472𝑟𝑒ℎ (𝐿⁄2)
𝑟′𝑤𝑛𝑝 = ℎ
(3.10)
2 ℎ
𝐿 0.472ℎ 𝐿
𝑎 [1 + √1 − ( ) ][ ] 𝑒 |(𝐷𝐻 𝐿 −𝐷𝑉)|𝑞𝑔
2𝑎 2𝑟𝑤
This implies:
𝑟 ′ 𝑤𝑝
𝑟′𝑤𝑛𝑝 = ℎ (3.11)
𝑒 |(𝐷𝐻 𝐿 −𝐷𝑉)|𝑞𝑔
Where; 𝑟 ′ 𝑤𝑝 is the Darcy pseudo-steady state effective wellbore radius; 𝑟′𝑤𝑛𝑝 is
the non-Darcy pseudo-steady state effective wellbore radius and (𝐷𝐻 ℎ/𝐿 − 𝐷𝑉 ) is
a measure of the horizontal isotropic rate-dependent skin deviation: 𝐷𝐻 ℎ/𝐿 is the
non-Darcy horizontal rate-dependent skin; dependent upon the reservoir
isotropic aspect ratio and 𝐷𝑉 is the vertical rate-dependent skin.
35
It should be of note that 𝐷𝐻 < 𝐷𝑉 . Hence, 𝐷𝐻 ℎ/𝐿 is always less than 𝐷𝑉 . 𝐷𝐻 is
given mathematically as:
𝛾𝑔 ℎ
𝐷𝐻 = 3.75 𝑥 10−8 (3.12)
𝜇𝑔 𝑟𝑤 𝐿2
and 𝐷𝑉 = 𝐷
Now, the difference between 𝐷𝐻 ℎ/𝐿 and 𝐷𝑉 i.e., (𝐷𝐻 ℎ/𝐿 − 𝐷𝑉 ) gives a negative
number. However, since both 𝐷𝐻 ℎ/𝐿 and 𝐷𝑉 are scaler (that is to say, skin is
scalar and so has no direction), their difference will be scaler. Hence, the
negative sign is only an indication of an improvement in the rate-dependent skin
along the horizontal well length – that is, the rate-dependent skin decreases (or
improves) with increasing horizontal well length as can be seen in equation (3.5).
Therefore, absolute values of the difference in skin alone may be considered in
the course of using equations (3.9) and (3.10) since the magnitude of this
difference is a prime concern.
𝐾𝐻 ℎ[𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 )]
𝑞𝑔 = (3.13)
𝑟𝑒
1424𝑇 [𝐼𝑛 ( ) + 𝑆]
𝑟 ′ 𝑤𝑝
and
𝐾𝐻 ℎ[𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 )]
𝑞𝑔 = (3.14)
𝑟𝑒 ℎ
1424𝑇 [𝐼𝑛 ( ′ ) + 𝑆 + |(𝐷𝐻 − 𝐷𝑉 )| 𝑞𝑔 ]
𝑟 𝑤𝑝 𝐿
𝑟𝑒 ℎ
1424𝑇 {𝐼𝑛 ( ) + 𝑆} 1424𝑇 |(𝐷𝐻 − 𝐷𝑉 )|
𝑟 ′ 𝑤𝑝 𝐿
[𝑚(𝑃𝑅 ) − 𝑚 (𝑃𝑤𝑓 )] = 𝑞𝑔 + 𝑞𝑔 2
𝐾𝐻 ℎ 𝐾𝐻 ℎ
(3.15)
36
Equation (3.15) is equivalent to: 𝐶 = 𝐴𝑞𝑔 + 𝐵𝑞𝑔 2
𝑟𝑒
1424𝑇 {𝐼𝑛 ( ) + 𝑆}
𝑟′ 𝑤𝑝
𝐴=
𝐾𝐻 ℎ
ℎ
1424𝑇 |(𝐷𝐻 − 𝐷𝑉 )|
𝐿
𝐵=
𝐾𝐻 ℎ
Example 3.2
Calculate the horizontal gas well IPR for both Darcy and non-Darcy pseudo-
steady state, using the data in Example 3.1. Assume; 𝑟𝑒ℎ = 𝑟𝑒𝑣 and L = 1000ft.
Solution
First, the rate-dependent non-Darcy skins, 𝐷𝐻 and 𝐷𝑉 , are determined. Then, the
half-axis of the drainage ellipsoid formed by the horizontal well of length, 1000ft
is determined so that the pseudo-steady state effective wellbore radius may be
deduced. Also, the skin improvement due to the horizontal drain-hole is
calculated.
𝛾𝑔 ℎ
𝐷𝐻 = 3.75 𝑥 10−8 = 2.359 𝑥 10−9 , (1⁄𝑆𝐶𝐹𝐷 )
𝜇𝑔 𝑟𝑤 𝐿2
0.5 0.5
2𝑟𝑒𝐻 4
𝑎 = 0.5𝐿 {0.5 + [0.25 + ( ) ] }
𝐿
37
0.5
4 0.5
2(1312.336)
= 0.5(1000) {0.5 + [0.25 + ( ) ] }
1000
= 1360.79ft
1000
0.472(1312.336) ( )
𝑟 ′ 𝑤𝑝 = 2
2 0.395
1000 0.472(395)
1360.79 [1 + √1 − ( ) ][ ]
2 𝑥 1360.39 2 𝑥 0.4593
= 14.460ft
ℎ 395
𝐷𝐻 = 2.359 𝑥 10−9 ( ) = 9.32 𝑥 10−10
𝐿 1000
Hence, the IPR for the horizontal pseudo-steady state Darcy flow can be solved
as:
Also, the IPR for the horizontal pseudo-steady state non-Darcy flow is calculated
thus:
1312.336
1424(650) {𝐼𝑛 ( ) − 3.4}
𝐴= 14.460 ≈ 5.9070
165900
1424(650)|(9.32 𝑥 10−10 − 1.513 𝑥 10−8 )|
𝐵= ≈ 7.917 𝑥 10−8
165900
𝐶 = [1.25 𝑥 109 − 1.78 𝑥 104 ] ≈ 1.25 𝑥 109
Hence,
38
[(5.9070)2 + 4(7.917 𝑥 10−8 )(1.25 × 109 )]1/2 − 5.9070
𝑞𝑔 = , 𝑀𝑠𝑐𝑓/𝑑
2(7.917 𝑥 10−8 )
= 93,769.757Mscf/d ≈ 93.770MMscf/d
The pseudo-steady state IPRs are calculated and displayed in tables (3.4) and
(3.5) below for the vertical and horizontal well scenarios respectively.
Table 3.4: Pseudo-steady State IPR Calculations for Vertical Gas Well
𝑷𝒘𝒇 𝒎(𝑷𝒘𝒇 ) ∆𝒎(𝒑) 𝝁𝒑𝒘𝒇 𝑫𝑽 q (scf/d)
(psi) (psi2/cp) (psi2/cp) (cp) 1/(scf/d) No Turb. Turb.
0 0 1.25E+09 0.0122 1.516E-08 31,108,459 29,306,584
14.7 1.78E+04 1.25E+09 0.0122 1.515E-08 31,108,017 29,306,189
400 1.33E+07 1.24E+09 0.0126 1.461E-08 30,776,634 29,010,930
800 5.31E+07 1.20E+09 0.0135 1.369E-08 29,785,667 28,126,036
1200 1.18E+08 1.13E+09 0.0147 1.256E-08 28,170,920 26,677,801
1600 2.04E+08 1.05E+09 0.0163 1.131E-08 26,019,090 24,735,474
2000 3.07E+08 9.42E+08 0.0183 1.007E-08 23,458,133 22,404,998
2250 3.78E+08 8.72E+08 0.0197 9.350E-09 21,708,407 20,800,691
2500 4.51E+08 7.99E+08 0.0212 8.692E-09 19,886,172 19,119,274
2750 5.26E+08 7.24E+08 0.0228 8.101E-09 18,017,996 17,383,990
3000 6.02E+08 6.48E+08 0.0243 7.581E-09 16,124,689 15,613,264
3250 6.78E+08 5.71E+08 0.0259 7.127E-09 14,221,623 13,820,880
3500 7.55E+08 4.95E+08 0.0274 6.731E-09 12,319,600 12,016,656
3750 8.31E+08 4.19E+08 0.0289 6.385E-09 10,425,897 10,207,314
4000 9.06E+08 3.43E+08 0.0303 6.083E-09 8,545,240 8,397,304
4250 9.81E+08 2.68E+08 0.0317 5.816E-09 6,680,597 6,589,501
4500 1.06E+09 1.94E+08 0.0331 5.579E-09 4,833,765 4,785,715
4750 1.13E+09 1.21E+08 0.0344 5.366E-09 3,005,782 2,987,063
5000 1.20E+09 4.81E+07 0.0357 5.173E-09 1,197,200 1,194,208
5150 1.24E+09 4.88E+06 0.0364 5.066E-09 121,469 121,438
5167 1.25E+09 0.00E+00 0.0365 5.055E-09 0 0
39
Table 3.5: Pseudo-steady State IPR Calculations for Horizontal Gas Well
𝑷𝒘𝒇 𝒎(𝑷𝒘𝒇 ) ∆𝒎(𝒑) 𝝁𝒑𝒘𝒇 𝑫𝑯 q (scf/d)
2 2
(psi) (psi /cp) (psi cp) (cp) 1/(scf/d) No Turb. Turb
0 0 1.25E+09 0.0122 2.36E-09 211,551,494 93,752,831
14.7 1.78E+04 1.25E+09 0.0122 2.36E-09 211,548,485 93,751,975
400 1.33E+07 1.24E+09 0.0126 2.28E-09 209,294,929 93,108,895
800 5.31E+07 1.20E+09 0.0135 2.14E-09 202,555,908 91,166,631
1200 1.18E+08 1.13E+09 0.0147 1.96E-09 191,574,907 87,937,268
1600 2.04E+08 1.05E+09 0.0163 1.76E-09 176,941,495 83,499,685
2000 3.07E+08 9.42E+08 0.0183 1.57E-09 159,525,841 77,996,022
2250 3.78E+08 8.72E+08 0.0197 1.46E-09 147,626,913 74,079,485
2500 4.51E+08 7.99E+08 0.0212 1.36E-09 135,234,897 69,848,587
2750 5.26E+08 7.24E+08 0.0228 1.26E-09 122,530,461 65,330,126
3000 6.02E+08 6.48E+08 0.0243 1.18E-09 109,655,125 60,538,225
3250 6.78E+08 5.71E+08 0.0259 1.11E-09 96,713,422 55,472,550
3500 7.55E+08 4.95E+08 0.0274 1.05E-09 83,778,813 50,116,861
3750 8.31E+08 4.19E+08 0.0289 9.96E-10 70,900,783 44,436,746
4000 9.06E+08 3.43E+08 0.0303 9.49E-10 58,111,470 38,375,025
4250 9.81E+08 2.68E+08 0.0317 9.07E-10 45,431,057 31,842,288
4500 1.06E+09 1.94E+08 0.0331 8.70E-10 32,871,770 24,697,176
4750 1.13E+09 1.21E+08 0.0344 8.37E-10 20,440,668 16,702,050
5000 1.20E+09 4.81E+07 0.0357 8.07E-10 8,141,496 7,406,341
5150 1.24E+09 4.88E+06 0.0364 7.90E-10 826,043 817,095
5167 1.25E+09 0.00E+00 0.0365 7.89E-10 0 0
40
CHAPTER FOUR
ANALYSIS AND COMPARISON OF RESULTS
From the results of the previous chapter, particularly from tables (3.4) and (3.5),
several plots are made to enhance effectively, the discussion, analysis and
comparison of the results so obtained. Thus, to compare the IPRs for vertical and
horizontal gas wells, the analysis and comparison are done graphically under the
following headings:
1. No turbulence drawdown (figure 4.1)
2. Turbulence drawdown (figure 4.2)
3. No turbulence Absolute Open Flow Potential (figure 4.3)
4. Turbulence Absolute Open Flow Potential (figure 4.4)
5. Vertical Well Turbulence effect (figure 4.5)
6. Horizontal Well Turbulence effect (figure 4.6)
In similar sense, the higher the wellbore real gas pseudo-pressure, the lower the
drawdown and hence lower the required energy. However, from figure (4.1) it
can be seen that a specified deliverability rate of 20MMscf/d requires a
corresponding change in the real gas pseudo-pressures, ∆𝑚(𝑃𝑤𝑓 ) of 8.0 x
108psi2/cp for the vertical well and 1.21 x 108psi2/cp for the horizontal well.
Hence, it can be deduced that, it requires a lower energy to produce at a
41
particular (required) rate for horizontal wells than it does for vertical wells in non-
turbulent flow.
42
Figure 4.2: Turbulence Drawdown Plot
From figure (4.3), the AOFP for the horizontal well (which is 211.55MMscf/d) is
higher than the AOFP of the vertical well (31.11MMscf/d) in the no turbulence
situation of flow. This means that the horizontal well has better productive
capabilities than does the vertical well.
43
Figure 4.3: No Turbulence AOFP Plot
44
Figure 4.4: Turbulence AOFP Plot
From figure (4.5) below, it can be deduced that, at real gas pseudo-pressure of
45
about 10.0 x 108psi2/cp (corresponding to 6.0MMscf/d) and above, the vertical
gas well exhibits laminar flow regime; below which the vertical gas well
experiences turbulence. The effect of turbulence (causing a reduction in the
flowrate) on the vertical well increases with decreasing real gas pseudo-pressure
around the wellbore and hence increasing drawdown (energy required). At a
given gas rate of 25 MMscf/D, the turbulence effect on the real gas pseudo
pressure is 0.4 x 108psi2/cp.
46
Figure 4.6: Horizontal Well Turbulence Effect
47
CHAPTER FIVE
CONCLUSION AND RECOMMENDATION
Considering the amount of energy, based upon the drawdown required for
producing at a given rate, it is pertinent to conclude that horizontal gas wells
require lower energy and stress to produce. Also, in terms of the gas wells
potential for flow, the horizontal gas well has higher potential to flow than the
vertical gas well. Howbeit, it is important to recap here that, the influence of
turbulence is to increase the pressure drop or pressure drawdown (around the
wellbore) required to produce at the given gas rate. Thus, the presence of
turbulence reduces net production from the well; both for the horizontal and
vertical gas wells. But the effect at a given flowing bottom-hole pressure is more
on the horizontal gas well than on the vertical gas wells. However, the turbulence
effect can be minimized by reducing fluid velocity near the wellbore – since the
highest fluid velocity is expected near the wellbore where flow converges. Hence
if there is no turbulence near the wellbore, there will be no turbulence in the
reservoir. The fluid velocity is reduced by increasing perforation/producing height
(ℎ𝑝 ) in the vertical well or the well length (L) for the horizontal well.
The horizontal-vertical well cost ratio in the Niger Delta, based upon how
sophisticated the available technology used is, is about 1.4 to 3 times; depending
upon drilling method and the completion technique employed and drilling
experience in the Niger Delta area. Obviously, the incremental cost of drilling
horizontal wells over vertical wells has reduced significantly over the last decade.
Presently in the Niger Delta, a typical horizontal wells costs 1.4 times the vertical
well cost. With time and extensive drilling experience, the horizontal well cost
would be reduced even further.
48
compared to the reserves; this is probably the reason why most of the Niger
Delta produced gas are flared rather than processed. Lack of the proper
processing size and handling equipment may also have been one of the reasons
for the very small gas need in this area. Otherwise, the following
recommendations are made for present (if the need arises) and future purposes
of improved technology, cost ratio and infrastructure.
Recommendations
Based on the whole work done, the following recommendations are made. These
recommendations should be considered if the need for gas production increases
with increasing market demand.
Horizontal gas wells should be drilled in the Niger Delta - with high
permeability reservoirs.
49
REFERENCES
1. Brown, K. E., “The Technology of Artificial Lift Methods”, Vol. 1, PennWell
Books, Tulsa, Oklahoma, pp 1-5
3. Gilbert, W.E., "Flowing and Gas-Lift Well Performance", Drilling and Production
Practice, 1954, API p. 143
5. Golan, M. and Whitson, C.H., Well Performance, IHRDC, Boston, MA, 1986, pp
21, 127-147
10. Fetkovich, M.J., "The Isochronal Testing of Oil Wells", SPE Paper No. 4529 -
48th Annual Fall Meeting of SPE of AIME, Las Vegas, Nevada, September 30 to
Octobers 3, 1973
11. Rawlins, E.L., and Schellhardt, M.A., "Back-Pressure Data on Natural gas Wells
and Their Application to Production Practices", U.S Bureau of Mines Monograph
7, 1936
50
12. Forchheimer, P., "Wasserbewegung durch Boston", Zeitz. Ver Deutsch Ing,
Berlin, 45,1901, pp 1781 – 1788
13. Green, L. and Duwez, P., "Fluid Flow through porous materials", J. Appl. Mech.
(March, 1951) 18, 39
14. Cornell, D. and Katz, D. L., "Flow of Gas through Consolidated Porous Media",
Ind. and Eng. Chem. (Oct., 1953) 45,2145
15. Al-Hussainy, R., and Ramey, H.J., Jr. "Application of Real Gas Theory to Well
Testing and Deliverability Forecasting", JPT, May 1966, p.637
16. Joshi, S.D., “Horizontal Well Technology”, PennWell Books, Tulsa, Oklahoma,
1991, pp 90 - 91, 329 - 220 and 345-350
17. Joshi, S.D., "Augmentation of Well Productivity with Slant and Horizontal Wells",
JPT, June 1988, pp.729 – 739
18. Economides, M.J., Deimbacher, F.X. Brand, C.W., and Heinemann, Z.E.,
"Comprehensive Simulation of Horizontal Well Performance", SPE 20717, 1990
and SPEFE, December 1991, pp418 -426
19. Lee, J and Wattenberg, R.A., “Gas Reservoir Engineering”, SPE inc, Richardson
TX, 1996, p.4
51
NOMENCLATURE
a = Half the major axis of an elliptical drainage Area.
A = Area of flow
C = Backpressure Constant
h = Reservoir thickness, ft
J = Productivity Index
= Fluid Viscosity, cp
P = Pressure, psia
52
PR = Initial Reservoir pressure, psia
T = Temperature, °R
= pi = 3.1416
Subscripts
H, h = horizontal
V = vertical
g = gas
53
o = oil
max. = maximum
pc = pseudo-critical
pr = pseudo-reduced
s = surface
Abbreviations
GOR = Gas-Oil Ratio
PI = Productivity Index
Hor. = Horizontal
Vert. = Vertical
Turb. = Turbulence
No-Turb. = No Turbulence
54
APPENDIX "A"
GAS CONDENSATE STUDY DATA
55
APPENDIX ”B"
Stewart et al Method for Calculating Pseudo-Critical Properties of the
C7+ fraction of Hydrocarbon Composition
0.15178 0.15427 3
𝑇𝑏𝑐7+ = (4.5579𝑀𝑐7+ 𝛾𝐶7+ ) , °𝑅 (𝐵. 1)
(𝐵. 2)
Example B.1:
Given the data in Appendix "A"; calculate the pseudo-critical pressure and
temperature of the C7+ fraction.
Solution:
From the data in Appendix "A":
Where:
56
PpcC7+ = Pseudo-critical pressure of the C7+ fraction, psia
57
APPENDIX "C"
Beggs and Brill (1973) Correlation for Calculating Z-factor
(1 − 𝐴) 𝐷
𝑍=𝐴+ + 𝐶𝑃𝑃𝑟 (𝐶. 1)
𝑒𝐵
Where:
0.066 2
0.32 6
𝐵 = (0.63 − 0.23𝑇𝑃𝑟 )𝑃𝑃𝑟 + [ − 0.037] 𝑃𝑝𝑟 + 9 𝑃𝑝𝑟
(𝑇𝑃𝑟 − 0.86) 10 (𝑇𝑝𝑟 − 1)
(𝐶. 3)
Example C.1:
Given, Tpc = 420.62°R; Ppc = 645.36psia; T = 650°R and P = 14.7 psia. Calculate
the gas compressibility factor (z-factor) using Beggs and Brill correlation.
Solution:
Tpr =1.546 and Ppr = 0.023
C ≈ 0.07145
D ≈ 0.97502
A ≈ 0.4422
B ≈ 0.006053768 ≈ 0.00605
Therefore:
Z = 0.9984
58
APPENDIX "D"
Lee et al Correlation for calculating Gas Viscosity
Where:
(9.4 + 0.02𝑀)𝑇1.5
𝐾= (𝐷. 2)
209 + 19𝑀 + 𝑇
986
𝑋 = 3.5 + + 0.01𝑀 (𝐷. 3)
𝑇
Example D.1:
Given T = 650°R; P = 14.7psia; Z = 0.998; M = 25.8514 lbm/lb-mole and 𝛾𝑔 =
0.8925.
Solution:
From equation (D.5);
59
𝜇𝑔 = 121.7197 × 10−4 𝑒𝑥𝑝[5.2754 𝑥 (8.7539 𝑥 10−4 )1.3449 ]
= 0.012177cp ≈ 0.0122cp
60
APPENDIX "E"
The Horizontal well Pseudo-steady State Rate Equation for Horizontal-
Vertical Well Productivity Comparison
(Developed by Sukubo, Ibianga Alu)
Introduction:
The basic concept used in this work of developing the horizontal well pseudo-
steady rate equation for effective horizontal and vertical well productivity
comparison, is the effective wellbore radius concept as discussed in Chapter 2 of
the project work. It should therefore be of note that the effective wellbore radius
concept is used to determine what the productivity would be for horizontal wells
(when horizontal well data are not available) from given vertical well data. This
concept, though theoretical, can be used for practical comparison in fields and
reservoirs where no horizontal well has been drilled.
Hence, the equation may be used to ascertain or justify one's reasons for or
against drilling horizontal wells in oil, gas or both reservoirs or fields. Even
though one's choice as a petroleum production engineer may be based on cost
effectiveness.
However the assumption made in course of developing my formula are the same
as the basic assumption made by Joshi S.D in developing his effective wellbore
radius concept. These assumptions are as follows:
Development:-
Joshi has already developed the effective wellbore radius concept for steady
state radial flow. However, the steady state radial flow system is only an ideal
case of flow. Hence, the equation developed herein is for pseudo-steady
state radial flow for both non-turbulent and turbulent regimes.
61
This implies;
𝑎+√𝑎 −( 2 𝐿⁄ )2
𝑟
] + 𝐼𝑎𝑛𝑖 ℎ⁄𝐿 [𝐼𝑛 (
2 𝐼𝑎𝑛𝑖 ℎ
[𝐼𝑛 ( 𝑒⁄𝑟 ′ ) − 0.75] = 𝐼𝑛 [ 𝐿⁄ ) − 0.75]
𝑤 2 𝑟𝑤 (𝐼𝑎𝑛𝑖 +1)
(𝐸. 2)
2 2
Let 𝐴 = 𝑎 + √𝑎2 − (𝐿⁄2) = 𝑎 (1 + √1 − (𝐿⁄2𝑎) )
Then the effective wellbore radius for the pseudo-steady state Darcy flow, 𝑟′𝑤𝑝 , now
becomes:
62
0.4724 𝑟𝑒ℎ (𝐿⁄2)
𝑟𝑤′ = 𝑟𝑤𝑝
′
= (𝐸. 3)
𝐼𝑎𝑛𝑖 ℎ⁄𝐿
𝐼𝑎𝑛𝑖 ℎ ℎ⁄
𝐴( ) 𝑒 −0.75 𝐼𝑎𝑛𝑖 𝐿
𝑟𝑤 (𝐼𝑎𝑛𝑖 + 1)
Assuming an isotropic reservoir then, 𝐼𝑎𝑛𝑖 = 1
The horizontal well skin or skin due to horizontal drain hole, according to Joshi
now becomes:
′
rwp
S = −In ( ⁄ rw ) (E. 5)
Therefore, the horizontal well rate equation now becomes, in a simplified vertical
well rate form;
2
𝐾𝐻 ℎ(𝑃𝑒2 − 𝑃𝑤𝑓 )
𝑞𝑔 = (𝐸. 7)
𝑟
1424𝜇𝑔 𝑍𝑇 [𝐼𝑛 ( 𝑒⁄𝑟𝑤𝑝 ) + 𝑆]
𝐾𝐻 ℎ[𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 )]
𝑞𝑔 = (𝐸. 8)
𝑟
1424 𝑇 [𝐼𝑛 ( 𝑒⁄𝑟𝑤𝑝 ) + 𝑆]
63
𝑟 𝐴 ℎ
[𝐼𝑛 ( 𝑒ℎ⁄𝑟 ′ ) − 0.75 + 𝐷𝑣 𝑞] = 𝐼𝑛 ( ) + ℎ⁄𝐿 [𝐼𝑛 ( ) − 0.75 + 𝐷ℎ 𝑞]
𝑤 𝐿⁄ 2𝑟𝑤
2
(𝐸. 9)
ℎ
𝑟𝑒ℎ 𝐴 ℎ 𝐿 ℎ ℎ
𝐼𝑛 ( ′ ) + 𝐼𝑛 𝑒 −0.75 − 𝐼𝑛 𝑒 −𝐷𝑣𝑞 = 𝐼𝑛 ( ) + 𝐼𝑛 ( ) + 𝐼𝑛 𝑒 −0.75 𝐿 − 𝐼𝑛 𝑒 −𝐷ℎ 𝑞 𝐿
𝑟𝑤 𝐿⁄ 2𝑟𝑤
2
ℎ
𝑟𝑒ℎ 𝐴 ℎ 𝐿 −0.75 ℎ 𝑒 −𝐷𝑣𝑞
[𝐼𝑛 ( ′ ) 𝑒 −0.75 ] = 𝐼𝑛 [( )( ) 𝑒 𝐿 ] + 𝐼𝑛(
ℎ)
𝑟𝑤 𝐿⁄ 2𝑟𝑤 −𝐷ℎ 𝑞
2 𝑒 𝐿
ℎ
𝑟𝑒ℎ 𝐴 ℎ 𝐿 −0.75 ℎ (−𝐷𝑣𝑞+𝐷ℎ 𝑞 ℎ)
( ′ ) 𝑒 −0.75 = ( )( ) 𝑒 𝐿𝑒 𝐿
𝑟𝑤 𝐿⁄ 2𝑟𝑤
2
Then the effective wellbore radius for the pseudo-steady state non-Darcy
′
flow, 𝑟𝑤𝑛𝑝 , now becomes:
0.4724𝑟𝑒ℎ (𝐿⁄2)
′
𝑟𝑤𝑛𝑝 = ℎ (𝐸. 10)
ℎ 𝐿 ℎ ℎ
𝐴 (2𝑟 ) 𝑒 −0.75 𝐿 𝑒 (𝐷ℎ 𝐿 −𝐷𝑣)𝑞
𝑤
This implies:
′
𝑟𝑤𝑝 ℎ
′
𝑟𝑤𝑛𝑝 = ℎ
′
= 𝑟𝑤𝑝 ∗ 𝑒 −(𝐷ℎ 𝐿 −𝐷𝑣)𝑞
𝑒 (𝐷ℎ 𝐿 −𝐷𝑣)𝑞
ℎ ℎ
But, 𝐷ℎ << 𝐷𝑣 hence, 𝐷ℎ − 𝐷𝑣 shall be negative. However, the rate-dependent
𝐿 𝐿
′
skin cannot be negative and the implication is that 𝑟𝑤𝑛𝑝 cannot be greater
′
than 𝑟𝑤𝑝 . Hence, absolute value is considered. Thus:
ℎ
′
𝑟𝑤𝑛𝑝 ′
= 𝑟𝑤𝑝 ∗ 𝑒 −(|𝐷ℎ 𝐿 −𝐷𝑣|)𝑞
64
7.08 ∗ 10−3 𝐾𝐻 ℎ∆𝑃
𝑞𝑜 = (𝐸. 11)
𝑟 ℎ
𝜇𝑜 𝐵𝑜 [𝐼𝑛 ( 𝑒⁄𝑟𝑤𝑝 ) + 𝑆 + (|𝐷ℎ 𝐿 − 𝐷𝑣 |) 𝑞𝑜 ]
2
𝐾𝐻 ℎ(𝑃𝑒2 − 𝑃𝑤𝑓 )
𝑞𝑔 = (𝐸. 12)
𝑟 ℎ
1424𝜇𝑔 𝑍𝑇 [𝐼𝑛 ( 𝑒⁄𝑟𝑤𝑝 ) + 𝑆 + (|𝐷ℎ 𝐿 − 𝐷𝑣 |) 𝑞𝑔 ]
𝐾𝐻 ℎ[𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 )]
𝑞𝑔 = (𝐸. 13)
𝑟 ℎ
1424 𝑇 [𝐼𝑛 ( 𝑒⁄𝑟𝑤𝑝 ) + 𝑆 + (|𝐷ℎ 𝐿 − 𝐷𝑣 |) 𝑞𝑔 ]
These equations are solvable by reducing them to the equivalent forms of the
Forchheimer equation. For instance, equation (E.13) is thus solvable:
𝑟 ℎ
1424 𝑇 [𝐼𝑛 ( 𝑒⁄𝑟𝑤𝑝 ) + 𝑆] 1424 𝑇 (|𝐷ℎ 𝐿 − 𝐷𝑣 |)
𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 ) = 𝑞𝑔 + 𝑞𝑔2
𝐾𝐻 ℎ 𝐾𝐻 ℎ
(𝐸. 14)
𝑟 ℎ
1424 𝑇[𝐼𝑛( 𝑒⁄𝑟𝑤𝑝)+𝑆] 1424 𝑇(|𝐷ℎ −𝐷𝑣 |)
Let C = 𝑚(𝑃𝑅 ) − 𝑚(𝑃𝑤𝑓 ), A = and B = 𝐿
𝐾𝐻 ℎ 𝐾𝐻 ℎ
𝐶 = 𝐴𝑞𝑔 + 𝐵𝑞𝑔2
𝐵𝑞𝑔2 + 𝐴𝑞𝑔 + 𝐶 = 0 .
√𝐴2 + 4𝐵𝐶 − 𝐴
𝑞𝑔 = (𝐸. 15)
2𝐵
65