BP Exploration Completion Design Manual (Vol 1)
BP Exploration Completion Design Manual (Vol 1)
Contents
INTRODUCTION 1
INTRODUCTION
Page 1
How to use the Manual – Chapter 1 BP Exploration
Completion Design Manual
Page 2
BP Exploration How to use the Manual – Chapter 1
Completion Design Manual
Statement of requirements
Page 3
How to use the Manual – Chapter 1 BP Exploration
Completion Design Manual
Page 4
BP Exploration How to use the Manual – Chapter 1
Completion Design Manual
Detailed completion
design
Prepare technical
specifications
Issue enquiry
Commercial evaluation
of bids
Presenting a recommendation
to a bid committee
Pre-production meeting
with vendors
QA and inspection
Witness onshore
make-up
Run completion
Page 5
How to use the Manual – Chapter 1 BP Exploration
Completion Design Manual
2. Developing a Statement
of Requirements
3. Conceptual
Completion Design a. Reservoir performance
d. Artificial lift
5. Architecture of Completions
a. Optimizing the completion
in the near wellbore region
8. Selection of
Completion Equipment
9. Procurement of Equipment
Page 6
BP Exploration How to use the Manual – Chapter 1
Completion Design Manual
To cope with the large scope and overall length of the manual
and to make the manual manageable, each key topic has
been split into five colour coded subsections. Figure 1.5
illustrates the breakdown. Subdividing the key topics in this
manner allows the reader to quickly locate the subsection
appropriate to his knowledge and interest. It is felt that all
users, irrespective of discipline, should read the summary
and practical guidelines for all topics in order to achieve a
greater awareness of the whole process.
The first section in any chapter will deal with the theory or
background information associated with each topic. This
section will primarily be used by readers new to the topic or
by users with a particular problem which requires further
theoretical investigation. The techniques section details the
actual procedures which should be employed, and includes
details of how to use the appropriate software or company
procedures. The worked examples illustrate how the
techniques are employed in a real situation. The summary
and practical guidelines essentially provide a quick reference
containing rules of thumb associated with any topic. The
final subsection contains a comprehensive bibliography and
list of references together with a list of company specialists
who can be contacted for further assistance. The aim of this
section is to provide a means of obtaining further information
and specialist advice on problems not covered in sufficient
detail in the manual.
Page 7
How to use the Manual – Chapter 1 BP Exploration
Completion Design Manual
Page 8
BP Exploration How to use the Manual – Chapter 1
Completion Design Manual
CONTACTS,
REFERENCES,
ETC
WORKED
EXAMPLE
TECHNIQUES
THEORY
CONTACTS,
REFERENCES,
ETC
SUMMARY AND
Example: Architecture of PRACTICAL
Completions GUIDELINES
WORKED
EXAMPLE
PROCEDURE
BACKGROUND
Page 9
How to use the Manual – Chapter 1 BP Exploration
Completion Design Manual
Each chapter contains the subdivided Key Topic supported as necessary by the
associated Subsidiary Topic(s)
b SUBSIDIARY TOPIC
a SUBSIDIARY TOPIC
CONTACTS,
REFERENCES,
ETC
WORKED EXAMPLE
TECHNIQUES
THEORY
Page 10
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual CONTENTS
BACKGROUND
INTRODUCTION 1
TECHNIQUES
INTRODUCTION 13
WORKED EXAMPLE
SUMMARY 31
CONTACTS 33
REFERENCES 33
BIBLIOGRAPHY 33
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual BACKGROUND
INTRODUCTION
Page 1
Developing a Statement of Requirements – Chapter 2 BP Exploration
BACKGROUND Completion Design Manual
Page 2
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual BACKGROUND
Page 3
Developing a Statement of Requirements – Chapter 2 BP Exploration
BACKGROUND Completion Design Manual
Page 4
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual BACKGROUND
DIFFICULTIES IN
ESTABLISHING
PRECISE AND
AGREED
OBJECTIVES
Optimizing completion design requires the best estimate
available of the individual well behaviour throughout its life.
Traditionally, individual well performance, production
problems and characteristics have not been addressed at an
early enough stage in the field development plan.
Page 5
Developing a Statement of Requirements – Chapter 2 BP Exploration
BACKGROUND Completion Design Manual
Page 6
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual BACKGROUND
Definitions
of System
Terminology
Objective
An objective defines the purpose of the system. It is also a
standard with which to compare performance. The overriding
objective of a well completion is to safely maximize the
profitability of the well throughout its life.
Environment
The environment of a system consists of the natural setting,
physical properties and administrative factors that affect the
Page 7
Developing a Statement of Requirements – Chapter 2 BP Exploration
BACKGROUND Completion Design Manual
Constraints
System constraints are inherent parts of the system which
can limit its performance. Some control can be exerted over
system constraints. Some examples are scale, sand
production and corrosion. The factors that create these
constraints include conditions for carbonate precipitation,
acid gases, and poor or no cementation of sand grains, all of
which are part of the environment since they cannot be
controlled. However, the results of these factors, scale, sand
or corrosion problems, can be controlled depending on the
well system and completion design.
Page 8
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual BACKGROUND
Table 2.1
Page 9
Developing a Statement of Requirements – Chapter 2 BP Exploration
BACKGROUND Completion Design Manual
Table 2.2
Deposits:
Produced oil: Oil, water, and treating fluid
Asphaltenes (acids/frac fluids) compatibility
Paraffins analyses and bottom hole samples
Produced water:
Scale
Salt
Table 2.3
Page 10
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual BACKGROUND
Table 2.4
Resources
System resources include money, materials, manpower,
equipment and certain system characteristics that can be
utilized and controlled to meet the system objectives. Some
resources such as reservoir pressure and temperature and
reservoir drive may seem to be part of the environment.
Resources can often be controlled, changed or improved,
whereas the environment is a given and non-controllable
resource. The more controllable and flexible resources come
from outside the system, ie materials, manpower and money.
Page 11
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual TECHNIQUES
INTRODUCTION
2) Interdisciplinary discussions.
Page 13
Developing a Statement of Requirements – Chapter 2 BP Exploration
TECHNIQUES Completion Design Manual
PHASE 1:
REVIEWING
DATA BASE AND
ESTABLISHING
INITIAL
OBJECTIVES
The first step in developing the SOR is to understand the
initial objectives of the client (geoscience and reservoir
departments).
After the client objectives are known, the next step is to have
each discipline involved make a thorough review of the
available data base and set its initial objectives independent
of the other disciplines. This phase includes:
Page 14
Completion Design Manual
BP Exploration
PHASES FOR DEVELOPING AN SOR
PHASE 1 2 3 4 5
COMPREHENSIVE DISCUSSIONS ANALYSIS FEEDBACK WRITE SOR
REVIEW OF PREVIOUS BETWEEN PLANNING
EXPERIENCE GROUPS
ACTIVITY • Exploration, Reservoir, • Each group presents: • Validate ERCs. • Meet with all groups. Write developed SOR
Drilling, Completion and Objectives, ERC. • Analyse objectives • Present initial SOR and Conceptual
Production Groups • Reconcile obvious for viability. with alternatives. Completion Plan.
establish their initial conflicts. • Identify and quantify • Discuss affects of risk
objectives. Review uncertainty. and uncertainty.
data base. • Determine if • Answer questions
• Identify good ideas, constraints are real. pertaining to SOR.
procedures. • Establish further data • Reconcile differences,
• Avoid previous requirements. make compromises.
mistakes. • Develop IPRs and TPCs
for overall well system
behaviour, considering
Table 2.5
Page 15
future conditions.
• Establish alternatives
and limitations.
GUIDELINES
1. The goal of this process is to arrive at an SOR, not for each group to get its way. Groups must work together.
2. Present everything that is known and germane at the Phase 2 meeting. Each group must do this. If anything is held back it will complicate things later.
3. Avoid wishful thinking: use logical reasoning.
TECHNIQUES
4. Keep uncertainty and limitations reasonable.
5. Bias and experience must be tempered to project at hand.
6. Each group should realize that SOR is a compromise of everybody's ‘must haves’, ‘wants’ and ‘would like to haves’.
Developing a Statement of Requirements – Chapter 2 BP Exploration
TECHNIQUES Completion Design Manual
This review will include data from the discovery well and any
appraisal/development wells in order to establish objectives
for further drilling and recommendations for any future work.
From
Geoscience␣ and
Reservoir Groups
Statement of Objectives
Producing rate
Bottom hole location
Depth to zone(s) of interest
Plateau rate
Future recompletion of other zones
Potential change of well duty in future
Statement of Environment
(Give known values or best estimates possible, along with
assumptions)
Depth to zone(s) of interest
Trap classification(s)
Strike and dip of formation(s)
Rock properties (porosity, permeability, thickness, wettability
and saturations)
Page 16
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual TECHNIQUES
Statement of Resources
(Give best estimates, along with assumptions)
Reservoir drive
Reservoir pressure and temperature
General reservoir description (as expected):
thickness/permeability/porosity/etc.
Reserves: OOIP, OGIP
IPR and/or deliverability of subject well
Alternative completions
Statement of Constraints
Maximum offtake rates
Minimum bottom hole flowing pressure
From
Drilling␣ Group
Statement of Objectives
Well cost and duration
Well type, eg land, platform, subsea
Surface location
Deviation, kick-off point, measured depth,
horizontal displacement
Casing programme
Statement of Environment
Temperature and pressure profiles
Fluid properties
Thief zones
Aquifers
Salt or anhydrite sections
Page 17
Developing a Statement of Requirements – Chapter 2 BP Exploration
TECHNIQUES Completion Design Manual
Statement of Resources
Drilling budget
Type of rig
Statement of Constraints
Cost
Time
Drilling problems, eg hole stability
Environmental factors
Governmental regulations
From
Completion
Group
Statement of Objectives
Well production rate
Field production profile
Completion life
Type of well, eg platform, subsea etc
Statement of Environment
Fluid properties (as above)
Geothermal gradient
Reservoir pressure and temperature
Other zones, aquifers etc
Statement of Resources
Budget
Alternative completions
Artificial lift
Workover and well servicing
Well inflow performance
Reservoir drive
Stimulation
Statement of Constraints
Fluid properties (as above)
Formation damage
Rate sensitive sand production
Governmental Regulations
Page 18
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual TECHNIQUES
From Field
Operations/
Development
Statement of Objectives
Field production profile
Number of wells
Type of wells, eg platform, land etc
Total project cost
Statement of Environment
Surface pressure and temperatures
Fluid properties
Statement of Resources
Project budget
Statement of Constraints
Minimum wellhead pressure
Production facilities, eg topside weight offshore
Simultaneous operations, eg drilling and production
Governmental regulations
PHASE 2:
INTER-
DISCIPLINARY
DISCUSSIONS
Following Phase 1, the necessary discussions between the
exploration, reservoir, drilling, completion, and production
groups can begin. At these meetings each group presents
their objectives and ERC as they have specified them.
Page 19
Developing a Statement of Requirements – Chapter 2 BP Exploration
TECHNIQUES Completion Design Manual
PHASE 3:
ANALYSIS OF
THE COLLECTED
ERC AND
EXAMINATION
OF ALL
OBJECTIVES
In Phase 3 the completion design group takes the information
resulting from Phase 2, and compares the data and
assumptions with what is known of the ERC from the most
recent well(s). The comparison is made to ensure that the
most relevant data is used and estimates provided from
previous wells are␣ applicable.
Page 20
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual TECHNIQUES
PHASE 4:
RECONCILIATION
AND
MODIFICATION
OF INITIAL SOR
WITH THE
OTHER GROUPS
During this phase the completion designer presents the
initial SOR along with alternatives, risks and uncertainties.
The presentation will cover what effects these risks and
uncertainties may have on the completion and what
requirements and limitations are acceptable while still
satisfying the SOR. These discussions will also consider
probable well life and well duty, safety procedures, workover
philosophy, production goals and limitations, well operation
and regulatory requirements.
PHASE 5: WRITE
THE FINAL SOR
The completion group agrees upon the content and writes the
final SOR. This final SOR forms the basis for the conceptual
completion design.
Page 21
Developing a Statement of Requirements – Chapter 2 BP Exploration
TECHNIQUES Completion Design Manual
Page 22
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual WORKED EXAMPLE
OBJECTIVES
Producing rate: 17 000 stb/d
Bottom hole location: Northwest section of reservoir and at least the same distance from
the platform as the discovery well, approximately 2.5 km.
Future recompletion of
other zones: No specification
ENVIRONMENT
Trap classification: Stratigraphy (sandstone fan)
RESOURCES
Refer to the BP Drilling Proposal
Table 2.6
Page 23
Developing a Statement of Requirements – Chapter 2 BP Exploration
WORKED EXAMPLE Completion Design Manual
Bottom hole location: One (1) km NNE of the platform. This location should penetrate zone
up dip of OWC
Future recompletion of
other zones: None
ENVIRONMENT
Reservoir pressure
and temperature: 6100 psi and 250°F at 4085 m subsea
Rock properties: Logs and preliminary side wall core (SWC) analysis show porosity of
18 to 21% and permeability of 30 to 600 md, average of 200 md
Fluid properties: Oil, water and gas properties and analyses (recombined sample from
well 16/7b-24 as reported in the PVT report for this well). H2S
approximately 200 ppm, CO2 negligible There were problems
experienced in collecting the recombined sample, and the PVT report
is suspect.
Wax problems: No wax deposition report has been made, but an analogous field has
a wax appearance temperature of 170°F and a wax content of 13.7%.
Scale potential: No valid water analysis available. First available water sample should
be analysed.
RESOURCES
OOIP: 31 mmstb
CONSTRAINTS
Formation damage, poorer quality reservoir rock in certain areas, limited number of wells for gas injection,
insufficient data
Table 2.7
Page 24
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual WORKED EXAMPLE
Multirate tests and flowing gradient runs for IPR and TPC validation
ENVIRONMENT
Other formations of interest: 180 m salt section immediately above pay, aquifer below pay. No
thief zones. Potential disposal zone at approximately 1980 m.
CONSTRAINTS
Formation damage: Damage experienced with fresh water muds, recommend using oil-
base mud.
Estimated flowline
pressure and temperature: 800 to 1100 psig and 200 to 240°F
Table 2.8
PHASE 2
During Phase 2 of the SOR development, the following
conflicts were resolved:
Page 25
Developing a Statement of Requirements – Chapter 2 BP Exploration
WORKED EXAMPLE Completion Design Manual
Page 26
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual WORKED EXAMPLE
PHASE 5
Phase 5 resulted in a written SOR that included
the␣ following:
Objectives:
• Producing rate is designed for 15 000 stb/d.
Page 27
Developing a Statement of Requirements – Chapter 2 BP Exploration
WORKED EXAMPLE Completion Design Manual
Page 28
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual WORKED EXAMPLE
Environment
• Trap classification: Stratigraphy (sandstone fan).
Page 29
Developing a Statement of Requirements – Chapter 2 BP Exploration
WORKED EXAMPLE Completion Design Manual
Resources
• Refer to the BP Drilling Proposal.
• OOIP: 31 mmstb.
Constraints
• Formation damage: Damage experienced with fresh water
muds, recommend using oil-base mud.
Uncertainties
• Corrosion.
• Scaling potential.
• Oil-water contact.
• Coning potential.
• Gas toxicity.
• Hydrate potential.
• Formation damage and ability to correct.
• Reservoir quality and thickness.
Page 30
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual SUMMARY AND PRACTICAL GUIDELINES
SUMMARY
• An SOR is developed from a list of detailed objectives,
additional data requirements, resources, environment
factors and constraints for a given well completion.
Page 31
Developing a Statement of Requirements – Chapter 2 BP Exploration
SUMMARY AND PRACTICAL GUIDELINES Completion Design Manual
Page 32
BP Exploration Developing a Statement of Requirements – Chapter 2
Completion Design Manual CONTACTS, REFERENCES, ETC
CONTACTS
Phil Murray, SPE Well Technology Studies, Aberdeen x2898.
REFERENCES
2.1 Patton, L D and Abbot, W A: Well Completions and
Workovers – The Systems Approach, 2nd Edition,
Energy Publications, Dallas, TX (1985)
BIBLIOGRAPHY
Murray, P: ‘Systematic Completion Design – the Key to
Increased Well Performance’, presented at the 1989 Joint
Engineering Technology Exchange Conference, Montgomery,
TX (April 2 to 7, 1989)
Page 33
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual CONTENTS
BACKGROUND
BACKGROUND 1
TECHNIQUES
INTRODUCTION 3
DETERMINING THE RANGE OF WELL PERFORMANCE (CHAPTER 4) 3
NEAR WELLBORE PERFORMANCE (CHAPTER 4, SECTIONS 4b AND 5a) 5
ARCHITECTURE OF COMPLETIONS (CHAPTER 5) 6
TUBING MOVEMENT AND STRESS CALCULATIONS (CHAPTER 6) 7
SELECTION OF TUBULARS AND MATERIALS (CHAPTER 7 AND
SECTION 7a) 8
SELECTION OF COMPLETION EQUIPMENT (CHAPTER 8) 9
WELL SERVICING AND WORKOVER PHILOSOPHIES (CHAPTER 11) 9
WORKED EXAMPLE
SUMMARY 23
CONTACTS 25
BIBLIOGRAPHY 25
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual BACKGROUND
BACKGROUND
Page 1
Conceptual Completion Design – Chapter 3 BP Exploration
BACKGROUND Completion Design Manual
Development
DEVELOPING STATEMENT Well
prospect OF REQUIREMENTS prospect
CONCEPTUAL
COMPLETION DESIGN
Complete conceptual
completion design
Detailed well
Field plan including
development plan casing design
Review any changes to field
development plans and
assess impact on the basis
for detailed design
Page 2
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual TECHNIQUES
INTRODUCTION
Inflow
Performance
(Chapter 4 and
Section 4a)
Conceptual Design Consideration
– Predict the range of initial ‘near ideal’ (minimum near
wellbore effects) inflow performance and determine how
this performance is likely to change with time.
Data Required
– All well test and formation evaluation data from appraisal
or analogous wells, including RFTs, cores and logs.
Questions To Be Addressed
– What is the appropriate inflow performance relationship?
– What is the likely range of inflow performance and the
impact on the project?
– How will the reservoir drive affect well performance?
– How will the initial reservoir pressure vary?
– What is the likely change in GLR and water cut?
Page 3
Conceptual Completion Design – Chapter 3 BP Exploration
TECHNIQUES Completion Design Manual
Outflow
Performance
(Chapter 4 and
Section 4c)
Conceptual Design Considerations
– In order to determine production rates, a method of
accurately modelling tubing performance is required. This
model should be used to determine the initial range of
tubing performance and predict how this performance will
change with time.
Data Required
– PVT samples and flash data.
– Flowing gradient surveys from appraisal wells or analogous
wells to validate pressure drop prediction method.
– Initial view on likely wellhead flowing pressure (from SOR).
– Desired completion life (from SOR).
– Answers to questions on inflow performance raised earlier.
Questions To Be Addressed
– Is there an appropriate method for predicting tubing
performance and what are the potential errors and their
impact on the predictions?
– What is the likely range of tubing performance over the
area of the field and with time?
Combining the
Inflow and
Outflow
Performance
(Chapter 4)
Conceptual Design Considerations
– Confirm to an acceptable range of accuracy that the rates
in the SOR can be met. Determine the likely range of well
performance over the field life.
Data Required
– Results from above inflow and outflow analyses.
Page 4
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual TECHNIQUES
Questions To Be Addressed
– Is the range of initial rates in the SOR economically and
technically feasible?
– Is the selected WHFP the economic optimum. If not, what
are the repercussions of changing it?
– What is the appropriate tubing size for the range of inflow
performance and desired production rates?
– Will this tubing size sustain flow at the required rate over
the completion life?
– Is re-completing with a smaller tubing size economically
attractive?
– Is artificial lift likely to be a requirement? If so, what is the
conceptual technique?
NEAR
WELLBORE
PERFORMANCE
(CHAPTER 4,
SECTIONS
4b␣ AND␣ 5a)
Conceptual Design Considerations
– The likely well performance determined in the previous
section was based on achieving ‘near ideal’ inflow
performance (ie near zero skin). In the conceptual phase,
the aim is to identify wells where the near wellbore
performance is likely to be less than ‘ideal’ (positive skin),
or where there is a potential for better than ‘ideal’, ie
stimulation (negative skin).
Data Required
– Appraisal well tests should be designed to acquire all the
data necessary to avoid near wellbore problems in the
development wells.
Questions To Be Addressed
– Identify the potential types of near wellbore damage,
possible stimulation methods and potential sand
production problems?
– Evaluate the sensitivity to perforating techniques?
– Identify if near wellbore performance is a critical factor in
the field development and attempt to quantify the risks?
Page 5
Conceptual Completion Design – Chapter 3 BP Exploration
TECHNIQUES Completion Design Manual
ARCHITECTURE
OF
COMPLETIONS
(CHAPTER 5)
Conceptual Design Considerations
– Based on the above work, alternative designs for the
reservoir/wellbore interface need to be developed.
– Determine the number of zones to be completed and likely
method of production, ie segregated or commingled.
– Based on the tubing stress analysis and an evaluation of
the potential well servicing and workover techniques (see
below) alternatives for the casing/tubing interface should
be evaluated.
– Provide input to casing design.
Data Required
– Formation evaluation data from previous or analogous
wells.
– Results from well performance, near wellbore performance
and tubing stress analyses.
– Required completion life.
– Brief evaluation of potential well servicing and workover
philosophies.
Questions To Be Addressed
– Do all the alternative structures satisfy all the
requirements?
– Are the completions simple, safe, reliable and flexible?
– Do the completions have the ability to cope with
unforeseen problems?
– Does the design life of the completion tie in with the
production profile forecasts? (See Table 3.1.)
Page 6
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual TECHNIQUES
Table 3.1
TUBING
MOVEMENT AND
STRESS
CALCULATIONS
(CHAPTER 6)
Conceptual Design Considerations
– Identify whether standard stock tubulars are suitable for
the application.
– Can tubing movement be eliminated.
Data Required
– Pressures and temperatures associated with all likely
well␣ conditions.
Page 7
Conceptual Completion Design – Chapter 3 BP Exploration
TECHNIQUES Completion Design Manual
Questions To Be Addressed
– Will the selected strength of tubular satisfy all the
operating conditions with acceptable design factors?
SELECTION OF
TUBULARS AND
MATERIALS
(CHAPTER 7
AND␣
SECTION␣ 7a)
Conceptual Design Considerations
– Identify the alternative materials that can provide the
required completion life given the well conditions.
– Is the desired completion life and workover frequency
realistic?
– Evaluate the economic trade-off between corrosion
resistant materials and workover frequency.
Data Required
– Accurate samples and analyses of well fluids or reliable
data from analogous wells.
Questions To Be Addressed
– Is the data representative of conditions throughout
the␣ field?
– Is there any experience of using the same materials in
similar conditions that can be applied to the situation?
– Can material of suitable strength be made in the desired
corrosion resistant material?
– Will conditions change with time?
– Have all the fluids likely to be used been considered?
– Does the selected material and completion life provide the
lowest risk and best economic return?
– Is the lead time on the selected material compatible with
project timing?
Page 8
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual TECHNIQUES
SELECTION OF
COMPLETION
EQUIPMENT
(CHAPTER 8)
Conceptual Design Considerations
– This activity should be performed as part of the detailed
design process. However, if the completion requires special
equipment, care should be taken to ensure that the
equipment can be designed, tested and built within the
time frame of the project.
WELL SERVICING
AND WORKOVER
PHILOSOPHIES
(CHAPTER 11)
Conceptual Design Considerations
– Identify the techniques required to maintain the well
throughout its life.
Data Required
– Likely well duty.
– Likely production chemistry problems.
Questions To Be Addressed
– Are the conceptual completion designs compatible with the
well maintenance philosophy?
Page 9
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual WORKED EXAMPLE
DETERMINING
THE RANGE
OF WELL
PERFORMANCE
Inflow
Performance
The reservoir performance was based on a reservoir
simulation model, and the productivity indices, PIs, were
found to range from less than 10 stb/d/psi to almost
200 stb/d/psi.
Page 11
Conceptual Completion Design – Chapter 3 BP Exploration
WORKED EXAMPLE Completion Design Manual
Page 12
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual WORKED EXAMPLE
Combining
the␣ Inflow
and␣ Outflow
Performance
The decision on tubing size depends on the initial
productivity index, the timing of water breakthrough and the
decline of productivity with increasing water cut.
Formation
Damage
Mechanism
The main near wellbore performance problem in Miller is
believed to be calcium carbonate and barium sulphate scale
deposition in and around the perforation tunnels, although
the exact mechanisms are not yet fully understood.
Page 13
Conceptual Completion Design – Chapter 3 BP Exploration
WORKED EXAMPLE Completion Design Manual
Both salt water, fresh water and oil based mud have been
used to drill the reservoir section in the appraisal wells. The
experience so far shows a trend towards smaller skin factors
where oil based mud was used. This indicates that part of the
formation damage problem in Miller may be barium sulphate
scale deposition in the near wellbore area.
Perforating
Method
A number of perforating options were investigated. The need
for underbalanced perforating to minimize the perforating
skin, and hence the pressure drop across the perforations,
was evident to reduce the amount of calcium carbonate scale
being deposited in the perforations.
Page 14
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual WORKED EXAMPLE
Page 15
Conceptual Completion Design – Chapter 3 BP Exploration
WORKED EXAMPLE Completion Design Manual
Page 16
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual WORKED EXAMPLE
TUBING
MOVEMENT
AND STRESS
CALCULATIONS
The load cases expected to give maximum stresses or tubing
movements are:
Page 17
Conceptual Completion Design – Chapter 3 BP Exploration
WORKED EXAMPLE Completion Design Manual
Page 18
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual WORKED EXAMPLE
Page 19
Conceptual Completion Design – Chapter 3 BP Exploration
WORKED EXAMPLE Completion Design Manual
• Materials.
• Workovers and well servicing.
• Completion equipment.
Materials
Whichever metallurgy is selected for the tubulars, there will
be a trade off between material strength, weight and expense.
Workovers and
Well Servicing
Scale prevention measures, including zone isolation and
inhibitor squeezes, need a more detailed review. Methods to
be investigated include bridge plugs and retrievable straddle
packers, run on either wireline or coiled tubing. Limitations
due to depth, deviation, injectivity, pressure and corrosion
must be considered.
Page 20
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual WORKED EXAMPLE
Page 21
Conceptual Completion Design – Chapter 3 BP Exploration
WORKED EXAMPLE Completion Design Manual
9 5/8" SHOE
THE REQUIRED NIPPLES ARE INDICATED IN THE COMPLETION, ALTHOUGH THE SELECTION
OF NIPPLES AND NIPPLE SIZES NORMALLY BELONGS TO THE DETAILED COMPLETION DESIGN PHASE.
Figure 3.2 – Miller Conceptual Completion Design (7in and 5 1/2in, 5000 psi)
Page 22
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual SUMMARY AND PRACTICAL GUIDELINES
SUMMARY
The early analysis of development well performance as part of
the conceptual completion design is an important part of the
field development planning. This analysis can have a major
impact on the production profile, number of wells, process
design and hence project economics.
Page 23
BP Exploration Conceptual Completion Design – Chapter 3
Completion Design Manual CONTACTS, REFERENCES, ETC
CONTACTS
Phil Murray, SPE Well Technology Studies, Aberdeen x2898.
BIBLIOGRAPHY
Patton, L D and Abbot, W A: Well Completions and Workovers
– The Systems Approach, 2nd Edition, Energy Publications,
Dallas, TX (1985)
Page 25
Chapter 4 - Well Performance
BP Exploration
Jonathan Bellarby
20 October 1998
BP Completion Design Manual - Section 4
Contents
4.1 Introduction .......................................................................................................................3
4.1.1 The Definition of Well Performance ...........................................................................3
4.1.2 Inflow Performance Relationship (IPR) ......................................................................4
4.1.3 Vertical Lift Performance (VLP)..................................................................................5
4.1.4 Overall Pressure Drops ..............................................................................................5
4.1.5 Well Design for the Entire Life of the Well..................................................................5
4.1.6 Check List for Well Performance Modelling Using Prosper .......................................7
4.2 PVT ...................................................................................................................................9
4.2.1 Direct use of Laboratory Data...................................................................................11
4.2.2 Equation of State ......................................................................................................12
4.2.3 Empirical Correlations (Untuned) .............................................................................13
4.2.4 Empirical Correlations (Tuned).................................................................................14
4.3 Inflow Performance .........................................................................................................16
4.3.1 Radial Inflow Equation and Skin...............................................................................16
4.3.2 Vogel Inflow Performance ........................................................................................17
4.3.3 Fetkovich ..................................................................................................................17
4.3.4 Jones ........................................................................................................................17
4.3.5 Hydraulically Fractured Well.....................................................................................17
4.3.6 Horizontal Wells........................................................................................................18
4.3.7 Perforation and Deviation Skin .................................................................................19
4.3.8 Multiple Zone Completions .......................................................................................21
4.3.9 Gravel Pack Completions .........................................................................................21
4.4 Tubing Performance .......................................................................................................23
4.4.1 Which Correlation Should I Use? .............................................................................24
4.4.2 Using Correlations to Predict When a Well Dies......................................................27
4.4.3 Tuning Pressure Drop Correlations ..........................................................................27
4.4.4 Natural Flow .............................................................................................................28
4.4.5 Shut-in Conditions ....................................................................................................28
4.4.6 Artificial Lift ...............................................................................................................29
4.4.7 Heat Transfer ...........................................................................................................36
4.4.8 Erosion .....................................................................................................................37
4.4.9 Tubing Size Selection ...............................................................................................39
4.5 Well and Field Optimisation ............................................................................................44
4.5.1 Practical vs. Theoretical Optimisation ......................................................................46
4.6 Glossary ..........................................................................................................................48
4.1 Introduction
The purpose of this section of the manual is to define well performance and how it is
assessed and modelled. It has been written from a users viewpoint rather than any
theoretical perspective. The more theoretical aspects are covered in sections 4a, 4b, and
4c.
This section of the manual should be used in conjunction with the well performance
software “Prosper” manual. This section does not attempt to cover any of the mechanics
of using Prosper, rather what issues need to be addressed when attempting accurate and
representative well performance modelling. There are some areas of overlap, between the
Prosper manual and this section, however these have been included to emphasise
particular points.
6300
5600
4900
4200
Pressure (psig)
3500
2800
2100
1400
700
0
0 4000 8000 12000 16000 20000
Liquid Rate (STB/day)
Superimposed on the graph are two well properties. The blue line represents what the
reservoir will deliver to the well and is termed the inflow performance relationship (IPR).
The concept was introduced in the 1920’s with the development of the bottomhole
pressure gauge. From the gradient of this line the productivity index, or PI, can be
calculated and is quoted as stock tank barrels per day per psi drawdown. This gives a
means of comparing productivity between wells. The red line defines what the bottomhole
flowing pressure has to be in order for the well to produce at the rate on the bottom axis.
This relationship accounts for the pressure required to get the produced fluids up through
the tubing used to construct the well and is termed the vertical lift performance (VLP) or
tubing performance. If the datum point for each system is the same, the point at which the
two lines intersect defines the rate at which the well will flow and the bottom hole pressure
at this point. This is an example of nodal analysis, so called because the whole system
performance is described by two independent subsystems. The two subsystems are
upstream and downstream of a common node. It is possible for this node to anywhere in
the system, for example when looking at facilities performance as well, it is often useful to
define the node as the wellhead or where the flowline meets the manifold. Such concepts
are explored more in the well optimisation section (4.5).
Having defined a means of describing well performance, the design of the well can be
considered in order to optimise this. In this context, optimisation isn’t necessarily the same
as maximise. Certainly high rate wells can be designed, however, this may have
implications on the performance of the reservoir throughout field life that may reduce
overall hydrocarbon recovery. The key issue is that a means of reviewing well
performance is available to contribute to the overall review of the best production strategy
for the reservoir.
As shown in the graph above, there are two parameters to consider when assessing well
performance:
The Inflow into the well - described by the IPR
The outflow through the well to surface - described by the VLP
An analogous relationship exists for injection wells.
Oil
Pressure
Gas
Separator
Pressure
• The fluids handling and how this may change e.g. there is no point in designing a large
well to cope with high GORs if the wells will have to be constantly choked back
because the gas can’t be managed at surface.
• Production problems. Scales, asphaltenes, wax etc. all impact the well design and
productivity. Asphaltene deposition occurs at a certain pressure and therefore the
bottom hole pressure is maintained above this point in order to prevent asphaltene
depositing in the reservoir or at the bottom of the well (e.g. perforations).
• The likelihood of workovers. This will allow for a tubing size change or implementation
of artificial lift.
PVT
R Bubble point or dew point
R Oil density at standard conditions
R Formation Volume Factor at Reservoir Pressure and bubble point
R GOR (if available for below the bubble point, this must be constant composition
data) or CGR (again constant composition data is a must), alternatively CGR for
gas wells (with CGR values below the dew point, these should again be for
constant composition expansion).
R Oil viscosity at reservoir, bubble point and separator conditions (constant
composition)
R Gas gravity and viscosity
R CO2,, N2 and H2S concentrations
R Water salinity or density
R Emulsion data (viscosity vs. water cut)
R Confidence level in the PVT data i.e. sample method and information on
consistency checks
Reservoir
R Geometry - no flow boundaries and well position
R Well inclination w.r.t. reservoir
R Net reservoir penetration and how this is arranged i.e. partial penetration
R Horizontal permeability, large (reservoir) and small (inches) scale kv/kh
R Relative permeability (oil / water / gas) if available
R Reservoir net thickness
R Porosity and UCS (for perforation performance determination)
R Pressure and datum depth and any expected variation over time
R Expected water cuts and GORs / CGRs - with any variations over time
R Any formation damage tendencies or responses to well kills
Well Test Data
R Completion schematic with dimensions, deviation data and depths
R Fluid rates - water, gas, oil - how they are measured and the confidence in them -
ideally also a variety of flowrates, covering likely rates over field life.
R Surface
and downhole pressures and temperatures - ideally not just spot
measurements, but covering both ranges in time and depth
Completion
R Liner size and weight, tubing size(s) and weights
R Casing sizes and depths (especially if accurate temperature data is required)
R Completion fluid (density and whether oil or water)
R Potential restrictions (e.g. wireline retrievable sub-surface safety valves)
R Deviation Survey
R Water depth (if offshore), RT elevation
R Tubing metallurgy and condition (scale and / or corrosion) or direct roughness data
Artificial Lift
R Gas lift - compressor discharge and wellhead pressure(s) and gas gravity
R ESP - power availability, setting depth, casing size
R Pump performance data (if available)
Surface
R Separator pressure
R Flowline dimensions and lengths (including manifolds) from tree to separator
R Flowline /manifold geometry (90° bends and the like)
R Choke sizes
4.2 PVT
The Pressure Volume Temperature (PVT) relationship describes how a fluid behaves
under changing conditions. With an accurate PVT relationship, the density, viscosity and
gas-oil ratio for the fluid under expected pressures and temperatures can be reliably
extracted. This is then used to determine the inflow performance and more importantly the
tubing performance.
As a hydrocarbon fluid is produced, the temperature and pressure changes. These
changes will initially only cause the oil or gas to change viscosity and density. At a certain
point however the fluid will change from a single phase to two phase. For a black-oil fluid,
the gas will start to come out of solution at the bubble point. For a condensate,
condensate will start to come out of solution. For a given fluid composition, these points
(the saturation points) define the phase envelope. Outside the phase envelope, the fluid is
single phase, inside the fluid is two phase. The critical point is the point on the phase
diagram where to the left the fluid that first comes out of solution is a gas. To the right, the
fluid that first comes out of solution is a liquid. Figure 3 shows an example of a phase
envelope.
Figure 3 - Example Phase Envelope
Pressure (psig)
2400 Cricondentherm
882.2 (degrees F)
2000
858.3 (psig)
1600
1200
Cricondenbar
800 431.5 (degrees F)
400
3913.1 (psig)
0
-200 0 200 400 600 800 1000 1200 1400 1600
Temperature (degrees F)
• A black-oil or volatile oil is one where with a drop in pressure, gas will come out of
solution. The difference between a black oil and a volatile oil is purely arbitrary and
relates to the higher GOR and formation volume factor of a volatile oil compared to a
black oil. The FVF for a volatile oil will be above approximately 1.5., below this it will be
a black-oil
• A retrograde condensate is one where with a drop in pressure, liquid will first come out
of solution. Note most retrograde condensates will exhibit the behavior where as the
pressure is reduced further, the liquid may vaporise again.
• A gas is where a drop in pressure will not result in the phase envelope being crossed.
These definitions apply to the conditions in the reservoir (isothermal). In the tubing,
cooling will occur and therefore the fluid will almost certainly cross the phase envelope at
some point. Therefore a gas will produce some liquid (condensate) before it reaches the
surface.
It is possible to have all of these fluids with the same composition, it is just the initial
pressure and temperature that may change. For example in the phase envelope (Figure
3), a reservoir temperature of 650°F or below would result in a black or volatile oil. A
reservoir temperature of between 650°F and 858°F would be a retrograde condensate
and a reservoir temperature of above 858°F would be a dry gas. Note in this example a
black-oil fluid is likely as most reservoir temperatures are below 650°F!
Below the saturation pressure, the proportions of the phases will change. This can be
examined on the phase envelope plot (as in Figure 4).
Figure 4 - Fluid Behavior From Reservoir to Surface
Black
4000 Oil
3600 Gas
3200
2800
2400
Pressure (psig)
2000
Gas
Gas
40%
1600
60%
Gas
1200
20%
800
Gas
80%
400
0
-200 0 200 400 600 800 1000 1200
Temperature (degrees F)
As the fluid proportions and properties change, this will clearly effect the well
performance. The reservoir or near wellbore performance will be effected by:
1. The viscosity
20 October 1998 Page 10
BP Completion Design Manual - Section 4
2. The expansion or contraction of the fluid (included in the formation volume factor)
3. Any relative permeability effects due to liquid or gas break-out (note Prosper can not
handle gas relative permeabilities and these will have to be estimated separately)
The tubing performance on the other hand will be effected by:
1. The density of the fluid(s) - accounted for by the FVF and oil / gas gravity
2. The proportion of gas to liquid (the GOR or CGR)
3. The viscosity (to a lesser extent)
The PVT data is therefore critical to the well performance predictions. Historically most
errors in well performance prediction have been attributed to poor or inaccurate PVT data.
Regardless of which PVT system is used to describe the fluids within Prosper, it is vital
that the original laboratory PVT data is both representative of the reservoir fluids and
covers well performance issues as well as reservoir performance issues:
1. Ensure that the well has cleaned-up adequately (as evidenced by stable flow with a
steady GOR and water cut).
2. There are two options, either bottom hole single phase samples or separator samples,
recombined according to the GOR. Sampling at the separator introduces more errors
than downhole sampling as the proportions of oil and gas have to be accurately
measured. Which ever method is chosen, ensure that the fluid is single phase at the
sandface otherwise condensate or gas breakout may make the produced fluids
unrepresentative. Sampling of fluids where the reservoir pressure is close to the
bubble or dew point is always going to be error prone and this uncertainty must be
acknowledged.
3. The value of accurate samples is huge. For example on Pompano, approximately US
$20 million could have been saved if the paraffin content had been accurately known
and the expensive TFL completions avoided. The only oil sample was unfortunately
“lost” in the laboratory.
4. Multiple samples should be checked against each other. The bubble point is the most
useful consistency check.
5. The PVT analysis should include a constant composition expansion experiment for
conditions between the bubble point and separator pressure. Note, the reservoir
engineer is more likely to be interested in differential expansion where at each stage of
the expansion, the produced gas is removed and therefore the composition changes.
The two methods will produce different GORs, and Formation Volume Factors and
hence different predictions about well performance. Constant composition is more
valid in the tubing as the oil and gas will be in constant contact with each other.
6. Any PVT model may be ideal for the reservoir or the facilities, but not the tubing. For
example the correct fluid density may not be critical for either the reservoir or the
facilities model. However for the tubing the correct fluid density over a large range of
pressures and temperatures is absolutely vital. Likewise, the reservoir fluids do not
need to include a large sensitivity to temperature and facilities correlations do not
necessarily need to cover a large range in pressure.
Phase Envelope
Pressure (psig)
50000 Critical Point
45000 656.7(degrees F)
35000
30000 Cricondentherm
886.3 (degrees F)
25000
904.4 (psig)
20000
15000
Cricondenbar
10000 317.9 (degrees F)
0
-100 0 100 200 300 400 500 600 700 800 900
Temperature (degrees F)
4. It is known that many EoS models poorly represent the liquid density of fluids and
errors of around 5% are commonplace. Errors in liquid densities would not significantly
effect either reservoir models or most facilities models, but will seriously effect the
tubing performance.
5. An EoS model is not particularly accurate at determining the viscosity of fluids. This
can be improved by the entering of “critical volumes” for each component.
Gas Viscosity (note if a gas is chosen, all the condensate drop out is assumed to occur
in the separator - NOT in the tubing). Therefore any condensate hold-up problems can
NOT be analysed.
• Lee at al
• Carr et al
Any of these correlations can be used directly without tuning so long as the GOR and
densities are known. It would be possible to examine each of these correlations in turn to
see which correlation was developed to model a similar fluid to one in question. However
this implies a complete lack of understanding of the reservoir fluids. At the very least, the
bubble point should be matched and the correlation chosen accordingly.
Retrograde Condensate
Prosper has its own Retrograde condensate empirical model that is capable of predicting
fluid properties and condensate gas ratios (CGRs) below the dew point. Even more than a
black oil model, this model is only realistic if it is matched to real conditions.
0.8
0.7
Untuned
Viscosity (cP)
0.6
Data
0.5
0.4
0.3
0.2
Tuned to
reservoir Tuned to reservoir
0.1 viscosity and wellhead
viscosity
0
800 1200 1600 2000 2400 2800 3200 3600 4000 4400
Pressure (psi)
Note it is not possible to treat any of these effects in isolation. For example a well
perforated in a small proportion of the reservoir will make any poor perforation effects
more pronounced.
More details on the inflow performance is found in the BP Near Wellbore Performance
1
Manual
Pwf Pwf
2
Where Qb and Pb are a flowrate and a corresponding bottom hole bubble point pressure.
In order to use this correlation therefore test data must be available.
4.3.3 Fetkovich
Q = J o ( Pr2 − Pwf2 )n
Jo is the P.I. above the bubble point adjusted for any relative permeability effects
This empirical relationship has similar application to the Vogel IPR and should be used
below the bubble point.
4.3.4 Jones
( Pr − Pwf ) = aQ 2 + bQ
This empirical equation is an expansion of the Darcy inflow equation to include rate
dependent or non Darcy effects and is therefore applicable to oil and gas wells. The b
term is the same for the Darcy equation. The a term accounts for the flow velocity:
0.0005359 Bo2 ρ
a=
k 1.201hp2rw
Where hp is the completed interval.
Note the Jones equation does not account for any flow convergence around perforations
or for any case that is not a simple open hole completion.
1
Near Wellbore Performance Manual WEO-W08 September 1990
Where Kf and K are the permeabilities of the proppant and reservoir respectively, W f is the
fracture width and Xf is the fracture half length (distance from the well to the tip of the
fracture). Fracture widths are usually in the order of 0.15 to 0.75 inch. A fully effective
fracture should have an Fcd > 10. Note the permeability of proppant is much less under
real conditions than measured in the laboratory due to compression, gel residues, fines
entrapment, long term crushing, non Darcy effects or contamination. A useful rule of
thumb is to reduce the lab derived permeability value by a factor of 10. The fracture width
will also be reduced by embedment of the proppant into the rock. In severe cases, this
can lead to a severe reduction in the effectiveness of the fracture and techniques for
extending the width of the fracture to compensate (i.e. tip-screen out techniques) should
be used.
The model used within Prosper is a transient model and accounts for the early time being
dominated by fracture conductivity and the later time being dominated more by overall
reservoir performance. An example of the effect this has is shown in Figure 7.
Figure 7 - Propped Fracture Transient Performance
4000
1:Time (days)
3600
0=0.18
1=0.32
2=0.58
3200 3=1.05
4=1.90
5=3.42
2800 6=6.16
7=11.10
8=20.00
8 9=40.00
2400
Pressure (psig)
9 0
7
2000 5
6
1600 1
4 2
1200 3
800
400
0
0 2000 4000 6000 8000 10000 12000 14000 16000 18000
In this example, only after 20 days, has the well settled down to close to steady state
performance.
The third method allows for pressure drops along the completion interval itself. A number
of horizontal well models are available including the preferred ones of Goode and
2
Wilkinson and Kuckuk and Goode. Note both of these formulations are valid even when
the horizontal wellbore length approaches the length or width of the reservoir. The method
of Goode and Wilkinson is less recommended when frictional pressure drops are
included, but is preferred when friction is ignored (infinite conductivity). Note in order to
Use Goode and Wilkinson with Infinite conductivity, set the roughness to 0 (this turns off
friction, rather than assuming friction for a smooth pipe. The other horizontal well models
are included for completeness and should not be used in practice.
Note the Wong and Clifford model for a deviated well will be available in Prosper v.6 by
late 1998.
2
Inflow Performance of Partially Open Horizontal Wells P.A. Goode and D.J. Wilkinson SPE 19341 August
1991
Perforation Diameter This is the diameter in the rock (not the casing). It can be
extracted from the gun vendor with a correction for the rock
strength. Alternatively the equation
D = EH N 80 [327
. − 0.61 ln(UCS )] 3, for deep penetrating
charges where D is the perforation hole diameter, UCS is the
unconfined compressive rock strength and EHN80 is the
entrance hole diameter in N80 steel (available from perforation
catalogues).
Shots per Foot For underbalanced perforating this will be the gun’s shot per
foot. If overbalance perforating with dirty fluids, the open
perforation may only be 0.1 - 0.25 the nominal shot per foot.
Perforation Length This is a property of the gun and charge, with corrections for
confining stress and rock strength. These corrections are vital
and the equations will soon be incorporated into Prosper but
are available in THoR and from the gun companies. Typical
numbers are between 12 and 24 inches
Damaged Zone This will depend on the mud and rock properties and how the
Thickness wells are drilled. Typical values are 3 to 12 in. The actual
values. The value is not particularly critical - so long as the
perforations extend behind this. The actual thickness can be
determined from filter cake tests on core or estimated from
filtrate depth of invasion calculated from mud losses - difficult
to achieve in practice.
Damaged Zone This is the damage caused by the mud filtrate invasion into a
Permeability reservoir. The ratio between this and the formation permeability
can be anything from 0 to 99% damage with an average of
around 0.8. Core flood tests (return permeability) are vital for
this. Note that the damage profile is also important. For
example with a surfactant adsorption damage mechanism, the
damage is likely to be most severe near the wellbore - where it
has least effect in a perforated completion, by the time the
filtrate reaches the total invasion depth, the surfactant may be
completely adsorbed and therefore no further damage occurs.
Crushed Zone From CAT scans and thin sections of perforated core, this is
Thickness typically 0.5" thick.
Crushed Zone Assuming an optimum perforation underbalance, this should be
Permeability the same as the formation permeability, otherwise a figure of
up to half the formation permeability can be used.
Shot Phasing A choice of the engineer! Low values (60 - 90 degrees) will
normally give optimum perforation efficiency.
Deviation A choice of the engineers!
Penetration This is the fraction of the wellbore that is open to flow. Prosper
assumes that it the top interval only that is completed (i.e. from
the top boundary of the reservoir downwards).
Vertical Permeability There is some uncertainty around what figure should be used.
In THoR and KT Perf, there was a natural division between
short and large scale k v/kh. Short scale Kv/kh applied to the
scale of the perforations i.e. on the scale of a few inches. Large
scale kv/kh applied to the effects of anisotropy on the deviation
or partial penetration skin. Prosper assumes that the large and
small scale kv/kh are the same. In a vertical fully completed
well, the small scale kv/kh should be used, but in a high angle or
partially completed well, the large scale effects may be more
3
Underbalance criteria for Minimum Perforation Damage L.A. Behrmann SPE 30081, May 1995
Note, the criticality of getting a large number of large diameter perforations filled with high
permeability gravel. In the event that the perforations are not effectively packed, they will
fill with reservoir sand. The permeability of the loose reservoir sand is likely to be much
lower than properly sorted and clean gravel.
The gravity term is the weight of the fluid. For a single phase system, it is simply the
density of the fluid. For multiphase systems, the density is a function of the relative
proportions of the different phases. If the phases are travelling at the same speed (no
slip), then the proportion of each phase will depend on the PVT properties and the
pressure. In reality the gas will travel faster than the liquid and this will increase the
proportion of liquid (termed the hold-up). Correlations are used to determine the amount
of slip and therefore the proportion of each phase. The correlations are dependent on
accurate PVT data in order to make valid predictions of the slippage.
The relative contribution of the friction and the density will vary with production rate (Figure
8) and also tubing size. Note that the average density reduces as the flowrate increases.
This is because the slippage will reduce until the hold-up approaches the no slip hold-up.
Below a certain rate the density increases markedly and will approach the liquid density at
zero flowrate. This is because the slippage is excessive and in fact the liquid will not be
produced at a fast enough rate to be produced out of the well. When this happens the
liquid builds up until it is produced out in a large slug. This is called heading and is outside
of the scope of the pressure drop correlations available for steady state flow - this is
because it is a transient or time dependent phenomena.
Total
Pressure Drop (psi / ft)
Density
Friction
Flowrate (bpd)
The importance of this correlation is that it should under predict pressure drop. Therefore
when plotted against real data it should lie to the left of the real (measured) data. If it lies
to the right it indicates either errors in the measured data, errors in friction or more likely
errors in density predictions i.e. PVT problems.
4.4.1.4 Orkiszewski
4
This correlation is an amalgam of various correlations including Duns and Ros . It works
well over a wide range of conditions and often appears as the closest correlation to match
well test results. Unfortunately there is a discontinuity in the predictions at a mixture
velocity of 10 ft/s. This results in a jump in predicted pressures (as seen in Figure 9). This
means that the correlation is of little practical use.
Figure 9 - Orkiszewski VLP Predictions
4400
4200
Pressure (psia)
4000
4.5" Tubing
3800
3600
5.5" Tubing
3400
3200
3000
0 4000 8000 12000
4
Predicting Two Phase Pressure Drops in Vertical Pipe J. Orkiszewski JPT June 1967 SPE 1546
4.4.1.6 Ansari
6
This correlation is recently available in Prosper. It has been available in THoR and
Multiflo for some time. It is a “mechanistic” model in that it predicts pressure drops based
on the mechanism i.e. the flow regime. This is similar to the more sophisticated traditional
methods such as Duns and Ros. However newer mechanistic models take this further by
eliminating more empirical steps and consider the transitional regimes further.
The Ansari correlation should therefore cover all flow regimes, all sizes of tubing, all
deviations and everything from heavy oil to gas. In a comparison of this model across a
7
number of oil and gas fields , at a variety of conditions, this correlation was the most
widely applicable - although by no means perfect, with average absolute errors of around
7%. There is still concern that this correlation (and many others) is not fully applicable for
deviated wells.
This correlation finds greatest application with low GOR, heavy or medium oil fluids such
as Harding, Foinaven and Andrew. It is acceptable for use on higher GOR fluids and
gasses, but correlations such as Gray are often better for gas.
5
A Study of Two-Phase Flow in Inclined Pipes H. Beggs and J. Brill May 1973 JPT
6
A Comprehensive mechanistic Model for Upward Two-phase Flow in Wellbores. A. Ansari, et al SPE 20630
Sept 1990
7
An Evaluation of Recent “Mechanistic” Models of Multiphase Flow for Predicting Pressure Drops in Oil and
Gas Wells J. Pucknell at al SPE 26682 Sept 1993
4.4.1.8 Gray
This correlation was developed for vertical gas condensate wells and was originally used
in the API 14B subsurface safety valve sizing program. The published limits of the
1
correlation are a flow velocity of 50 ft/sec, 3 /2" tubing or less, and a condensate ratio of
under 50 bbl/MMscf. In reality however, Gray has proven accuracy in low to moderate
CGR gas wells and is useful for predicting liquid loading. It is empirical and uses its own
PVT module for condensate and liquid prediction. It is therefore not compatible with
compositional PVT models.
Mukerjee Brill
2500 Ansari
Petroleum Experts
2000 Petroleum Experts 2
1500
1000
500
0
10000 20000
Liquid Rate (STB/day)
In general the following caveats apply:
1. Don’t use the correlation outside of its intended regime i.e. don’t use Duns and Ros for
oil.
2. Avoid using Orkiszewski, Fancher Brown, Beggs and Brill or Hagedorn Brown for
predicting tubing lift minima.
3. In the absence of a verified correlation for the field, Ansari has proven good at
predicting minima over a wide range of fluids.
3000
True Vertical Depth (feet)
4000
5000
Fluid G a s
Flow Tubing
6000 T y p e Producer
7000
10000
0 100 200 300 400 500 600 700 800 900 1000 1100
Clearly this will frequently tend to be an over estimation of the wellhead pressure,
particularly if the fluids are not gassy or a lot of water is being produced. What
frequently happens is that the high gas oil ratios are associated with depletion and the
therefore high gas oil ratios (and therefore a simplification of gas only) do not combine
with high reservoir pressures to give high wellhead pressures. However as a first pass
it is quick, simple and conservative.
2. The next approach is a bit more subtle. If the BHP vs Rate could be extrapolated down
to a zero rate without any of the problems associated with varying hold-up at low rates,
we would be close to extracting a realistic pressure gradient. This is easily achieved
using the Fancher - Brown no slip correlation. This correlation will give the highest
proportion of in-situ gas and therefore will be conservative. This gradient will be
realistic for conditions immediately after a shut-in. Unfortunately phase segregation will
occur - this can result in higher pressures than immediately after . Fortunately the oil
will always drop in comparison to the gas. This means that there will be no further
phase transition from oil to gas. Therefore as the standard volume of gas is known at
the time of shut-in, it is conservative to assume that the amount of space occupied by
the gas remains the same. This then means that the standard gas volume = (1-total
hold-up) x completion volume. The actual gas volume is this standard gas volume
corrected to the pressure in the well.
3. Use a transient model - if you can find one?
9500
1
9000 1
0
8500
0
Liquid Rate (STB/day)
1
8000
0
Gaslift Gas Injection Rate
7500 Curve 0 = 4 (MMscf/day)
Curve 1 = 8
7000 0
1
6500
6000
5500
5000 01
0 2000 4000 6000 8000 10000
Figure 13 - Gas Lift Injection Depth vs. Rate for 5 Untuned Correlations
15000
13500
12000
Orkiszewski
10500 0
0 Hagedorn Brown
Liquid Rate (STB/day)
0
0 0 0
9000 0
0 0 0 0
0
0 0
0 00
7500 0
0 0 0 0
0
0
6000 Ansari 0 0
0
4500
Petroleum 0 0
3000 Experts 2
0
1500 Duns and Ros Modified
0
0 4000 8000 12000 16000
Injection Depth (feet)
11400
1 1
10800 1 1
1
10200 Maximum liquid
1 Optimum gas injection production for deep gas
9600 rate for deep gas lift lift
1
Liquid Rate (STB/day)
9000
Injection Depth
1
Curve 0 = 3000 (feet)
8400 Curve 1 = 12000
7800
1
7200 0 0 0 0
0 0
0 0
6600 0
Optimum gas Maximum liquid
6000 injection rate for production level for
01 shallow gas lift shallow gas lift
5400
0 3 6 9 12 15 18 21 24 27 30
0
1000 GASLIFT DESIGN DATA
Design Situation : New Well
2000 Valve Type : Casing Sensitive
Dome Pressure Correction
(above 1200psig) : Yes
3000
Vertical Lift Correlation : Petroleum Experts 3
True Vertical Depth (feet)
4000
Maximum Gas Available : (MMscf/day) 6.000
5000 Flowing Top Node Pressure : (psig) 500.0
Operating Injection Pressure : (psig) 2700.0
Valve Water Cut : (percent) 90.000
6000 ACTUAL Liquid Rate : (STB/day) 7209.4
ACTUAL Oil Rate : (STB/day) 720.9
ACTUAL Gas Injection Rate : (MMscf/day) 5.975
7000 ACTUAL Injection Pressure : (psig) 2200.0
8000
Valve
9000
10000
Orifice
11000
12000
0 400 800 1200 1600 2000 2400 2800 3200 3600
Pressure (psig)
The design capabilities of Prosper are quite extensive and will not be gone into in detail
here. Some points that may be useful are:
1. There is a choice of valve type. All gas lift valves are sensitive to both tubing and
casing pressure. However casing pressure operated valves are more sensitive to
casing pressure and vice versa. The choice should not be dictated by performance and
uncertainty. A casing sensitive valve requires a drop in casing pressure to close the
valve, therefore potentially more valves may be required or a higher unloading
pressure. Figure 15 shows a plot of a gas lift design using a casing pressure drop to
close the unloading valves. However it is less sensitive to changes in reservoir
properties as the casing pressure is much easier to determine than the tubing
pressure. Tubing pressure operated valve systems are conversely harder to design as
they are more dependent on changes in reservoir properties.
2. The number of valves required is dependent on the separator or wellhead pressure
during the unloading and the assumption about the static gradient of the load fluid. The
gradient may be the completion fluid gradient, however this can produce a design that
is overkill for all other conditions. Consider the possibility of designing the gas lift
system so that the unloading (i.e. removing the completion fluid from the tubing) can
be performed by a nitrogen lift and routine kick-off performed as normal. This means
that the static gradient would become the worst case static gradient of the reservoir
fluid. This can be assumed to be 100% formation water or whatever.
3. Consider the point at which gas lift is required. During the early field life, there may be
little or substantial benefit in gas lift depending on the conditions and especially the
GOR of the produced fluid.
4. Do not design the system just for high rates and high reservoir pressures. This will
push the use of large orifice sizes which may introduce heading problems later on in
the well life.
5. Depending on the access difficulties, it can be useful to consider changing the orifice
valve position with time. As the reservoir depletes, the use of a deeper injection
position may be possible and beneficial. This implies that either two orifice valves are
run from day 1 or the lower mandrel is run with a dummy valve.
6. The pressure drops with gas between the casing and the tubing are usually not
significant. However Prosper has the capability to include these as required.
8 th
Energy Management in ESP Wells A. Simpson 5 European ESP Roundtable January 1998
Diffuser
At
Throat
An
Power Fluid
Reservoir
fluids
Nozzle
9
Fundamentals of Oilwell Jet Pumping A.W. Grupping, J.L.R. Coppes, J.G. Root SPE 15670 1986
2. Assume a constant heat transfer coefficient. Prosper has the capability to model heat
transfer with one overall heat transfer coefficient. Typical values are 4-6 BTU/h/ft2/F.
The heat capacity of the fluids is also required. The heat capacity of the oil will vary
considerably (factor of 2) between a light and a heavy oil. This is probably the best
method to use where there is some temperature data from well tests. The well test
configuration can be modelled and a trial and error approach used to match the heat
transfer coefficient with the measured temperature data. Care should be taken with
offshore wells however as an overall heat transfer coefficient will likely over simplify the
model.
3. Model the well in detail. Either Prosper or Welltemp (part of the Wellcat) package can
be used for detailed temperature modelling. Considerable detail is required including
all casing strings, annular contents, and rock types and properties. Unless all of this
information is rigorously included, the accuracy will not be any higher than option 2
above. This option is the only method that will allow an estimate of the transient heat
transfer effects rather than steady state.
Further details on heat transfer is contained in section 4c.
4.4.8 Erosion
Corrosion and erosion are fundamentally related to the flowrate and flow type. Therefore
Prosper is the ideal tool in which to assess the chances of erosion or erosion-corrosion.
There are no exact or theoretical relationships between the flowrate and the resulting
erosion. However, BP have further developed the empirical and conservative API erosion
limit with corrections for flow type, metallurgy, corrosion and sand.
For multiphase flow, a C factor is used, where:
C
VE =
ρm
Where ρm is the mixture density in lbs/ft , VE is the erosional velocity limit (ft/s) and C
3
10
Erosion Guidelines Revision 2.0 (1996) ESR.97.ER.002 John Martin BP Sunbury
* Sand free is defined as less than 1 parts per thousand barrels (pptb)
** The present velocity limits for "sand containing" conditions should only be applied up
to 275 pptb. Above this level the appropriate specialist should be consulted.
C/S - carbon or low-alloy steels, 13Cr - 13%Cr steels, DSS - duplex stainless steels
An example of extracting the C factor from Prosper is shown in Figure 17.
Figure 17 - C factor vs. Depth Extracted from Prosper for a Variety of Conditions
900
800
10,000 bpd Gas Oil Ratio 2000.0 (scf/STB) Tubing/Pipe Diameter 3.80 (inches)
10,000 bpd Gas Oil Ratio 6000.0 (scf/STB) Tubing/Pipe Diameter 4.80 (inches)
10,000 bpd Gas Oil Ratio 6000.0 (scf/STB) Tubing/Pipe Diameter 3.80 (inches)
700 10,000 bpd Gas Oil Ratio 2000.0 (scf/STB) Tubing/Pipe Diameter 4.80 (inches)
5,000 bpd Gas Oil Ratio 2000.0 (scf/STB) Tubing/Pipe Diameter 4.80 (inches)
5,000 bpd Gas Oil Ratio 6000.0 (scf/STB) Tubing/Pipe Diameter 3.80 (inches)
C=236
600 C=100
500
C factor
400
300
200
100
0
0 500 1000 1500 2000 2500 3000 3500 4000 4500
Depth (m)
Hig
h
Pr
es
su
re,
Mo rve
Lo
de
rat e Cu
w
e nc
PI ma
PI
an erfor
dp P
res ing
su Tub
re
Low Press
ure, High
PI
7500
7000 Initi
al R
6500 Pre eservo
ssu
re ir
6000 /2"
41
/8"
/2"
5500 /2"
27
31
51
Pressure
(psig)
5000
7"
4500
4000
Res
ervo
3500 afte ir Pres
r5y s
Res ears ure
erv
3000 afte oir Pre
r 10 ss
yea ure
2500 rs
2000
0 6000 12000 18000 24000 30000
Liquid Rate (STB/day)
In this example, the reservoir depletion is indicated by the changing inflow performance
curves. The minima of the tubing curves is also very clear, with 7" having a minima at
1
9,000 bpd, whereas 3 /2" tubing has a minima at around 2,000 bpd. The maximum rates
1
for initial reservoir pressures are 25,000 bpd for 7" and 11,000 bpd for 3 /2". However the
constraints of the erosional velocity must also be superimposed on this model.
With initial reservoir pressure, all the tubing sizes are stable, likewise after 5 years.
However after between 5 and 10 years, the wells will cease to flow on their own. In this
instance there is not much difference in when any of the tubing sizes shown become
unstable. This is because the productivities are high. With lower productivities, or for
example gas break out in the reservoir, the benefit of a smaller tubing size may be more
apparent.
For gas wells, the same approach can be used, however great care must be taken with
selecting the PVT representation of the gas. The PVT section (4.2) discusses the pitfalls
in the various Black-oil correlations and their assumption that the gas density can be
increased to represent the condensate.
The most important point about using this method is that care must be taken with the flow
correlation. The correlations section (section 4.4.1) indicates which correlations can be
used to accurately predict the instability region.
The next improvement in this approach that could be made is to examine the flow
behavior to the left of the minima to see if it is operationally acceptable to facilities. This is
covered in section.
One aspect of well performance in producers that is often overlooked is the location of the
bubble point and the impact that this has on tubing size. There will be no minima in the lift
curve if the fluid is single phase (i.e. all liquid or all gas). In these circumstances, the
bigger the tubing the better (from a well performance stand point). For example, in one
field, the bubble point was only 90 psi, therefore the fluid was single phase all the way to
surface. The development plan proposed had a 7" liner, but 4.5" tubing as this was
appropriate to the production rates expected. In this case, 7" tubing would not only have
reduced pressure drops slightly, but would have led to a monobore completion, with the
only downside, the increased cost of the 7" tubing and tree.
The same approach can be used for liner sizing. If the fluids will be single phase in the
liner and a long liner is to be used, then from a well performance stand point, a large liner
(e.g. 7") can be accepted. However this is not the only consideration for liner sizing and a
full discussion in contained in section 5.
Pressure (psig)
11 th
Uses of Reliable Downhole Monitoring In the Captain Field D.J. Cohen, 5 European ESP Round Table 1998
If severe slugging occurs above the pump, there is the potential for the wellbore to fill with
essentially pure liquid. Centrifugal pump stages are designed to be able to handle such a
condition, however this will increase the power requirement of the motor and must be
addressed in the ESP or hydraulic submersible pump design. The slugs are then a
problem for the facilities rather than a potential cause of the well dying. If the slugging
results in the loads on the pump varying, then in a severe case, this could lead to the
pump being periodically operated outside of the upthrust / downthrust envelope and
increased pump wear occurring. In the case of positive displacement pumps, the slugging
may place additional loads on the pump and rods and again the pump unit should be
designed to be able to handle a 100% liquid (water) head.
This phenomena means that much larger tubing can be effectively used. This is of
significant value in the case of electrical or hydraulic submersible pumps deployed inside
5 5
the tubing by means of coiled tubing or cable. The use of 7", 7 /8" or even 9 /8" tubing then
becomes possible.
Well A Oil
14000
Well B Oil
12000 Well C Oil
Well D Oil
10000
Oil Production
8000
6000
4000
2000
0
0 1 2 3 4 5 6 7 8
Gas Lift Injection Rate
20 70000
Well D Gas Injection Rate
Allocated Gas Injection Rate Well C Gas Injection Rate
Well B Gas Injection Rate
Well A Gas Injection Rate
15
Well D Production Rate 60000
Well C Production Rate
Well B Production Rate
Well A Production Rate
10
50000
40000
0 10
12
14
16
18
20
0
30000
-5
10000
-15
-20 0
Gas Injection Rate
In order to do this accurately, the wells and facilities must be accurately modelled and the
performance predictions must be accurate. Particular care should be exercised in tuning
the lift curve correlations on individual wells. Personally I recommend that a generic well
model is constructed to best match all the wells and that if tuning is required it is then
based on a large dataset and individual measurement errors on each well should not
grossly effect the result. See section 4.4.3 for further details on tuning correlations.
Combined Method
A theoretical model such as GAP is flexible in that it can offer predictions and answers to
what-if scenarios and is easily changed. It is vital that the model is however compared
with reality. I particularly recommend the use of multi-rate well tests on gas lifted wells as
frequently there can be discrepancies between theory and reality. If there are
discrepancies then these can be investigated in terms of the model, facilities measuring
errors and the likes.
4.6 Glossary
BHP .................................... Bottom hole pressure
BHT .................................... Bottom hole temperature
CGR ................................... Condensate gas ratio (the volume of dissolved liquid held in
the gas)
FVF..................................... Formation volume factor. This is the ratio between the volume
of a fluid at reservoir conditions and at standard conditions.
The is also called Bo
GOR ................................... Gas Oil Ratio (the volume of dissolved gas held in solution in
the liquid)
head ................................... The height a column of fluid can be supported. It is related to
the pressure by the density of the fluid
Liquid hold-up..................... The fraction of the pipe occupied by the liquid
IPR ..................................... Inflow Performance Relationship
Solution GOR ..................... The ratio of gas remaining in solution to the oil
Slippage ............................. The phenomena whereby two phases travel at different
velocities
Superficial Velocity ............. The velocity of a phase assuming that it takes up the entire
pipe area
THP .................................... Tubing head pressure
THT .................................... Tubing head temperature
UCS.................................... Unconfined compressive strength
VLP..................................... Vertical Lift Performance (or tubing performance)
w/c .............................. Water cut or base sediment and water (BS&W)
CONTENTS
INTRODUCTION 1
REFERENCES 41
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
INTRODUCTION
Page 1
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
RADIAL FLOW
CHARACTER-
ISTICS
Flow from reservoirs is characterized as transient, pseudo-
steady state or steady state flow depending on whether the
pressure response initiated by opening the well has reached
the drainage area boundary and on the type of boundary.
Transient flow occurs when the well is initially opened (or has
a significant rate change), and is a result of the pressure
disturbance moving towards the outer boundary of the
drainage area. During this period, production conditions at
the wellbore change rapidly, and the bottom hole flowing
pressure (BHFP), pwf , decreases exponentially with time.
Most drill stem tests (DSTs) and many well tests are
conducted under transient flow conditions and consequently
the observed productivity will often be greater than that
which will be seen in long-term production. Therefore, to
compute meaningful IPR values from short duration well
tests, appropriate corrections must be applied to compensate
for the transient flow behaviour, as well as the possible
variation in the skin effects.
When the disturbance meets the outer boundary, flow
becomes steady state or pseudo-steady state. If the boundary
is a constant pressure boundary, the reservoir pressure pR
will not change with time, and the flow is termed steady
state. If the boundary is a no-flow boundary, the reservoir
pressure pR will decline as a result of depletion, and the
resulting flow is termed pseudo-steady state. Steady state
equations should be used for oil wells where a waterflood is
in operation, or there is a strong natural water drive. Pseudo-
steady equations are used in the other cases. For a more
detailed discussion on this topic see Sections 2.2, 2.6 and 2.7
(*Reference 4a.1) of Golan & Whitson*.
Stabilization
Time
When the BHFP, pwf , appears to be constant, or declining
only slowly proportionately with time, the reservoir is said to
have achieved stabilized flow conditions. Under these
circumstances, the pseudo-steady state flow equations
presented in this section may be used to predict the long
term deliverability of a well.
Page 2
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
TRANSIENT
TRANSIENT
PSEUDO-STEADY PSEUDO-STEADY
pwf STATE p
wf STATE
log t t
θ µ ct r e2
(Equation 4a.1) t s = 948 (hours)
k
For wells producing from reservoirs of different shapes, the
stabilization time is given by:
θ µ ct A
(Equation 4a.2) t s = 3790 tDApss (hours)
k
where: θ = porosity, fraction
µ = viscosity, cp
ct = total fluid compressibility, psi-1
re = drainage radius, ft
A = drainage area, ft2
k = effective permeability, md
t DApss = the dimensionless time to stabilization for
different geometry and well placements
(*Reference 4a.1)
(see Table 2.4 of Golan and Whitson* or
(*Reference 4a.2) Table C.1 of Earlougher*)
Page 3
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
The Radial
Flow Equation
Consider an idealized well drilled in a fully oil bearing
reservoir rock of constant thickness and uniform
permeability. Assuming the producing interval is completely
penetrated by a circular wellbore, oil production will be
characterized by radial flow geometry. Under pseudo-steady
state conditions, the reservoir pressure declines very slowly
with cumulative production and Darcy’s equation for a
circular drainage area when expressed in practical oilfield
units, becomes:
p R – p wf
(Equation 4a.3) q o = 0.00708 k oh (stb/d)
µo B o
ln .472 rr e + S'
w
Page 4
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
k oh p R – p wf
(Equation 4a.4) qo = (stb/d)
141.2 µ oBo ln rr e – 0.75 + S
w
0.0397 k gh p R – p wf
(Equation 4a.5) qg = (scf/d)
µg Bg ln .472 rr e + S'
w
p
p
(Equation 4.18) m(p) = 2 dp
pb
µZ
Page 5
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
Hence:
k gh m p R – m p wf
(Equation 4a.6) qg = (mscf/d)
1422 T ln rr e – 0.75 + S'
w
Pseudo-
PVT Data Numerical Integration Pressures
p
µ 2p 2p 2p 2p
p Z ∆p × ∆p m = Σo ∆p
µZ µZ µZ µZ
(psia) (cp) (psia)2/cp
400 .01286 .937 66391 33196 400 13.278 x 106 13.278 x 106
800 .01390 .882 130508 98449 400 39.380 x 106 52.658 x 106
1200 .01530 .832 188537 159522 400 63.809 x 106 116.467 x 106
1600 .01680 .794 239894 214216 400 85.686 x 106 202.153 x 106
2000 .01840 .770 282326 261110 400 104.444 x 106 306.597 x 106
2400 .02010 .763 312983 297655 400 119.062 x 106 425.659 x 106
2800 .02170 .775 332986 322985 400 129.194 x 106 554.853 x 106
3200 .02340 .797 343167 338079 400 135.231 x 106 690.084 x 106
3600 .02500 .827 348247 345707 400 138.283 x 106 828.367 x 106
4000 .02660 .860 349711 348979 400 139.592 x 106 967.958 x 106
4400 .02831 .896 346924 348318 400 139.327 x 106 1107.285 x 106
2p is the average of 2p for the current row and the previous row.
µZ µZ
Table 4a.1
Page 6
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
(Equation 4a.7) Bg = ZT
pR + pwf
35.37
2
pR + pwf p 2 – p wf 2
p R – p wf = R
2 2
yields:
k gh p R2 – p wf 2
(Equation 4a.8) qg = (mscf/d)
1422 µ g ZT ln .472 rr e + S'
w
Note that the gas well IPR will always be a curved function
of pwf .
Page 7
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
PARAMETERS
IN RADIAL
FLOW
EQUATION
Permeability
The absolute permeability, k, of a reservoir rock is dependent
primarily on the rock structure and the effective confining
stress. It is a constant so long as one single-phase fluid (oil,
water, gas or air) saturates the entire porous media, flowrates
are low and the rock stress remains constant.
A D Where:
A = k o at S w = 0. Absolute permeability of the rock (k)
EFFECTIVE PERMEABILITY
Page 8
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
In-situ kh estimates derived from well test and DST data are
therefore preferred when making theoretical IPR calculations.
Even then, well test kh values may have to be adjusted to
take into account:
Page 9
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
Net Pay
Thickness
The net pay thickness, h, of a producing formation is the
average over the drainage area. It is usually based on a
combination of the analysis of open hole logs and the
geologist’s isopach maps (lines of equal pay thickness).
Portions of the formation which will not contribute
substantial quantities of hydrocarbons should be excluded.
These can usually be identified based on high water
saturations (from a resistivity log), poorly developed
porosities (from the density, neutron and sonic logs) and
poor sand development (from gamma ray (GR) or
spontaneous potential (SP) logs). Varous ‘net-to-gross’
criteria are used to define these ‘non-pay’ portions. The ratio
of net pay to total or ‘gross’ pay in a well is known as the
‘net-to-gross’ ratio. In a deviated well, or a well intersecting a
steeply dipping formation, the net pay value must be
adjusted to remove the apparent thickening caused by not
penetrating the pay zone perpendicular to the bedding plane.
Oil Formation
Volume Factor
Oil formation volume factor, Bo, is defined as the ratio of the
volume of oil (plus gas in solution) at reservoir pressure and
temperature to the volume of oil at standard conditions in
the stock tank. Hence:
Bo ≥ 1.0
Shrinkage = 1/Bo
Page 10
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
In light and medium oils, the viscosity of live oil (ie containing
dissolved gas) is directly dependent on the volume of
(*Reference 4a.7) dissolved gas (refer to Beal* and MULTIFLO for dead and live
oil viscosity correlations, respectively).
The viscosity of free gas is far below that of liquids. Since gas
is extremely compressible, its viscosity is strongly dependent
on both reservoir temperature and pressure as well as on the
gas composition and gravity. In the absence of laboratory
measured PVT data, µg can be estimated from gas
pseudo-reduced pressure and pseudo-reduced temperature
using the correlations developed by Carr, Kobayashi and
(*Reference 4a.9) Burrows*.
Page 11
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
Reservoir
Temperature
The reservoir temperature, TR, may be extrapolated from the
data recorded on well logs or recorded during testing. The
temperature of a reservoir can also be estimated from the
expression:
°R = °F + 460
Drainage and
Wellbore Radii
The drainage radius, re, is primarily a function of well
spacing. It will also depend, to some degree, on the relative
production rates of the well and its offsets.
Page 12
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
0.00708 k oh p R – p wf
(Equation 4a.9) qo ≈ (stb/d)
µoBo 7 + S'
Hence:
ln .472 rr e ≈ 8
w
k gh p R2 – p wf 2
(Equation 4a.10) qg ≈ (mscf/d)
1422 µ g ZT 8 + S'
Skin Effects
A deviation from the ideal flow model, which predicts a
logarithmic pressure change with distance from the well,
frequently occurs in the near wellbore region and is referred
to as the total skin effect S'.
rs
rw
ACTUAL PRESSURE
PRESSURE WHEN NO
SKIN PRESENT
p
SKIN REGION
Page 13
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
.00708 kh ∆ p skin
(Equation 4a.11) S' =
q µB
Page 14
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
h o = 30ft
6000
k o = 100md
µ o = 0.5cp
5000 B o = 1.5
s = -5
2000 s=
s 0
=
+5
s
s=
s = +10
=
+1
1000 +2
0
0
0
0
0 2500 5000 7500 10 000 12 500 15 000
OIL RATE , stb/d
Page 15
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
Reservoir
Pressure and
Bottom Hole
Flowing Pressure
The reservoir pressure, pR, is either found from an RFT log
(for uncased wells) or from pressure build-up test data, with
test data being the preferred method. For production
engineering calculations, it may be necessary to correct
reservoir pressure data from datum depth to the mid-point of
perforations (or top of perforations) using the fluid gradient
estimated for the reservoir.
OIL WELL
PRODUCTIVITY
EQUATION
Productivity
Index
q
(Equation 4.1) J =
p R – p wf
(stb/d/psi)
Page 16
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
The skin is obtained from well tests along with koh. However,
if the productivity index has changed with time and there is
no change in water cut and the well’s BHFP is above the
bubble point, the new skin can be calculated from the
following equation where Jo is the new oil productivity index.
This equation is derived from 4a.13.
qo
(Equation 4a.15) Specific PIs = (stb/d/psi ⋅ ft)
p R – p wf h
Page 17
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
so that:
0.5
–A + A 2 + 4B m p R – m p wf
(Equation 4a.18) qg = (mscf/d)
2B
Page 18
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
AOF
Page 19
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
LINE THROUGH
STABILIZED POINT
TO OBTAIN A
m (pR ) - m (p wf )
qg LINE THROUGH POINTS
TAKEN AT END OF EACH
FLOW PERIOD
USED TO OBTAIN B
A SLOPE = B
qg
Page 20
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
1422 T ln 0.472 rr e + S
w
(Equation 4a.20) A = psi 2/ c p / m s c f / d
k gh
or
3.16 × 10 –12 β γg T
(Equation 4a.22) B = p s i 2/ c p / m s c f / d2
µwf h p rw
2
Page 21
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
4.11 × 10 1 0
(Equation 4a.23) β =
k p1.33
2.73 × 10 1 0
(Equation 4a.24) β =
k p1.1045
B k gh
(Equation 4a.25) D = (mscf/d)–1
1422 T
or
2.22 × 10 –15 β γg k g h
(Equation 4a.26) D = (mscf/d)–1
h p rw µw f
2
Page 22
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
EFFECTS
OF WATER
PRODUCTION
Water
Production
in Oil Wells
As discussed in the Theory section of Chapter 4, the
production engineer is often more interested in the gross PI
than in the net oil productivity, when water encroachment
can be expected as a result of natural water drive or in
response to a waterflood scheme.
Page 23
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
OWC
SEGREGATED FLOW
TRANSITION
OWC
FRACTIONAL FLOW
+ 0.00708 k w h o (stb/d/psi)
µw Bw ln .472 rr e + S' o
w
The first two terms describe the segregated flow of oil and
water through separate reservoir intervals, while the first and
third term characterize the fractional flow of formation water
due to saturation changes in an oil zone or a transition zone.
Page 24
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
1.0
WATER CUT
GROSS
p wf 0.8
0.6
WATER
0.4
OIL 0.2
0
q GROSS q
Figure 4a.8 – Segregated Flow IPR Curve and Water Cut Plot
Page 25
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
ko + kw
µo B o µw B w
Figure 4a.9 shows how to determine the gross PI, Jgr, under
fractional flow conditions. This figure demonstrates how the
PI initially decreases with increasing water cuts, but at
higher water saturations the gross PI may start to rise.
Sw k ro k rw k ro / µ o Bo k ro / µ w B w PI Sw / PI Sw = o Water
Cut
0.30 1.00 0.000 1.282 — 1.000 0.000
0.40 0.92 0.010 1.179 0.033 0.945 0.027
0.45 0.75 0.025 0.962 0.083 0.815 0.079
0.50 0.40 0.050 0.641 0.167 0.630 0.207
0.55 0.15 0.070 0.192 0.233 0.332 0.548
0.60 0.06 0.100 0.077 0.333 0.315 0.812
0.65 0.02 0.135 0.026 0.450 0.371 0.945
0.70 0.00 0.180 — 0.600 0.468 1.000
1.0 1.0
PI Sw / PI ( S w = 0)
0.8 0.8
WATER CUT
0.6 0.6
qg
0.4 0.4
0.2 0.2
0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0
WATER CUT SW
Page 26
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
(d)
(a) DURING LATER PRODUCTION
BEFORE PRODUCTION WATER PERCOLATES PAST SHALES AGAIN AND
AND DURING VERY EARLY PRODUCTION FINGERS ALONG HIGH PERMEABILITY ZONES —
VERTICAL EQUILIBRIUM — IDEAL SITUATION VERTICAL EQUILIBRIUM LOST
K K
J J
H H
F F
E E
ORIGINAL OWC 2217 m ORIGINAL OWC 2217 m
(b) (e)
DURING LATER PRODUCTION DURING LATER PRODUCTION
WATER PERCOLATES PAST SHALES BRIDGE PLUGS SET AND VERTICAL
AND VERTICAL EQUILIBRIUM LOST EQUILIBRIUM RE-ESTABLISHED
K K
J J
H H
F F
(f)
(c) CURRENT POSITION
DURING LATER PRODUCTION WATER PERCOLATES PAST SHALES AND FINGERS
VERTICAL EQUILIBRIUM ALONG HIGH PERMEABILITY ZONES AGAIN.
RE-ESTABLISHED BY BRIDGE PLUGS FINAL LOSS OF VERTICAL EQUILIBRIUM
K K
J J
H H
F F
Page 27
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
Page 28
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
p wf
GAS
WATER
0
10 20 30 40 50
q g (mmscf/d)
0
100 200 300 400 500
q w (bbl/d)
85
FLOWING WGR = = 4.7 bbl/mmscf
18
CONDENSED WGR = 1.3 bbl/mmscf
EFFECT OF
PRESSURE
DEPLETION
Oil Wells
To develop meaningful production forecasts and appropriate
artificial lift designs, estimates of expected IPR curves
throughout the reservoir life are frequently required. These
IPR predictions must account for declining reservoir pressure
as well as the effects of changing fluid saturations.
Page 29
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
• Standing’s Method
To make this generally applicable, and to avoid confusion
(*Reference 4a.14) with measured PI, J', Standing (1970)* defined a term J*,
as the PI at low drawdowns derived from the pseudo-
steady state radial flow equation, Equation 4a.3. This
means the theoretical PI with no two-phase flow effects,
which is the case at low drawdowns.
• Eickmeier/Fetkovich
(*Reference 4a.15) Fetkovich* further simplified this by suggesting that the
change in relative permeability was directly proportional
to the change in the static reservoir pressures. Since the
viscosity and Bo effects are often relatively small and tend
to cancel out, then:
(Equation 4a.28)
J *f u t u r e = p R future
J *initial p R initial
(*Reference 4a.16)
Eickmeier* reportedly modified this concept and
suggested that a more realistic expression was:
J *future pR future
Page 30
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
• Empirical Methods
Empirical correlations have been used successfully to
predict future IPR curves in solution gas drive reservoirs.
The most common one is now often attributed to Kermit
(*Reference 4a.18) Brown*.
(*Reference 4a.19) Uhri and Blount* have also proposed a graphical, pivot
point method by which two well tests at different reservoir
pressures can be used to develop a forecast of the ongoing
IPR decline trend.
(Equation 4a.30) A = A
µgZ future µgZ initial
(Equation 4a.31)
B = B
Zf u t u r e Zi n i t i a l
Page 31
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
0.07
0.06
µg Z
0.05
µ g or µ g Z, cp
0.04 1.6
1.5
0.03 1.4
Z FACTOR
µg 1.3
0.02 Z 1.2
1.1
0.01 1.0
0.9
0 0.8
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10 000
PRESSURE, psi
EFFECT OF GAS
PRODUCTION
IN OIL WELLS
Vogel’s Two-
phase Flow IPR
Once the BHFP of a producing oil well drops below the
bubble point pressure, pb, the liberated gas causes a
build-up␣ of gas saturation in the vicinity of the wellbore.
Moreover, the volume occupied by this gas will be a function
of the flowing pressure. The resulting reduction of oil relative
Page 32
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
J *p R
(Equation 4a.32) q max = (stb/d)
1.8
0.8
0.6
p wf / p R
LINEAR PI
0.4
VOGEL IPR
100%
0.2 DRAWDOWN
0
0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
q / q max 9 q max
q max 5
q pwf pwf 2
(Equation 4.3)
q max = 1 – 0.2 pR – 0.8 pR
Page 33
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
1.00
p wf
pR
0.80
=
BOTTOM HOLE FLOWING PRESSURE
AVERAGE RESERVOIR PRESSURE
0.60
0.40
0.20
0
0 0.20 0.40 0.60 0.80 1.00
PRODUCING RATE q o
=
MAXIMUM PRODUCING RATE q o max
Figure 4a.14 –␣ Vogel’s Reference Curve for Solution Gas Drive Reservoirs
Page 34
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
1.0
p wf / p R , BOTTOM HOLE WELL PRESSURE, AS
A FRACTION OF RESERVOIR PRESSURE
0.8
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1.0
q/ q m , PRODUCING RATE AS A
FRACTION OF MAXIMUM RATE
Page 35
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
qi = qA – qb
qtmax = qb + qmax
ln rr e – 3
p – p' w f w 4
(Equation 4a.33) FE = R = ≈ 7
pR – pwf 7+S
ln rr e – 3 + S
w 4
Page 36
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
1.0
FRACTION OF RESERVOIR PRESSURE p wf / p R
0.8 FLOW
BOTTOM FLOWING PRESSURE,
FLOW 0.6
0.6 EFFICIENCY
0.5
1.1
1.5
1.3
1.7
2.1
2.3
2.5
1.9
0.
0.4
9
0.
1.2
1.0
1.4
1.6
2.0
2.2
1.8
2.4
0.4
0.5
1.0
0.
0.3
0.6
0.2 0.2
0.1
0 0
0 0.2 0.4 0.6 0.8 1.0 0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8
Generalized
IPR Curves
Because of the problems of using the Vogel/Standing model
in practical situations where the flow efficiency is uncertain,
it may be better to use the generalized IPR curves presented
in the Theory section of Chapter 4 for oil wells producing
below the bubble point, if multirate test data is available.
or
Page 37
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
Wells Producing
Free Gas
If the producing gas-oil ratio (GOR) is greater than that
expected on the basis of the solution GOR at that specific
reservoir pressure, then free gas is being produced from
the␣ reservoir.
p wf
Page 38
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
p R1
p R2
p R3
p wf
COMPOSITE IPR
Page 39
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
When the well is shut-in, cross flow will occur between the
higher pressured, lower permeability layers and the lower
pressured, higher permeability sections, in response to the
pressure differentials.
Page 40
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4a – RESERVOIR PERFORMANCE
REFERENCES
Page 41
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4a – RESERVOIR PERFORMANCE Completion Design Manual
Page 42
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
CONTENTS
INTRODUCTION 1
NON-DARCY SKIN 22
ONSET OF NON-DARCY FLOW 23
NON-DARCY SKIN IN GAS WELLS 23
NON-DARCY SKIN IN OIL WELLS 25
PRACTICAL CONSIDERATIONS WITH RESPECT TO NON-DARCY SKIN 25
REFERENCES 34
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
INTRODUCTION
In most wells, the total skin, S', is the only IPR parameter
that can be changed by drilling, completion, stimulation,
workover and production practices, and is therefore a primary
concern of production engineers.
Page 1
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
Lp DRILLING DAMAGE
CEMENT (Part of Formation Damage)
CASING
CONVERGING
FLOW FROM
h
h p
PARTIAL
COMPLETION
rd
DEEP FORMATION DAMAGE
rs
rw = WELLBORE RADIUS
rd = RADIUS OF DRILLING DAMAGE
rs = RADIUS OF FORMATION DAMAGE
Lp = EFFECTIVE LENGTH OF PERFORATION
The total skin, S', measured in a well test may include both
low rate, Darcy skin effects, S, and high rate, non-Darcy flow
effects, Dq , as discussed in Section 4a – Reservoir
Performance.
Page 2
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
(Equation 4b.1) S = h Sp + Sc
hp
Page 3
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
pb DARCY
IDEAL IPR
SKIN EFFECT(S)
ACTUAL IPR 2
p R - p wf = aq + bq
p wf
q
PRODUCTION LOST
DUE TO EFFECTIVE SKIN (S')
IDEAL IPR
RADIAL FLOW EQUATION
PR m(p r ) — m(p wf ) = Aq
P wf TURBULENCE
EFFECTS (S TB )
DARCY
ACTUAL IPR SKIN EFFECT(S)
2
m(p R ) - m(p wf ) = A q+ B q
Q
PRODUCTION LOST
DUE TO EFFECTIVE SKIN (S')
Page 4
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
Page 5
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
PARTICLES IN WELLBORE
FLUIDS BLOCK PORE THROATS
CEMENT
ku
CASING
ks SOLIDS IN
WELLBORE
FLUIDS
SHALE
SWELLING
SAND GRAINS
SAND GRAINS
ADJACENT TO
WELLBORE SHALE SWELLING
SOLIDS IN
WELLBORE TRAPPED
FLUIDS HIGH WATER
CAN BLOCK
PERFORATIONS SATURATION
PRODUCED FINES
TRAPPED
HIGH WATER
SATURATION
rw
rs
PRODUCED
FINES
k u – ks
(Equation 4b.2) Sd = ln rr s
ks w
Page 6
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
PERFORATION
SKIN
Perforation skin, Sp, is caused by additional pressure drops
as the radial flow of reservoir fluids deviates to spherical/
cylindrical flow into the perforations and crosses the crushed
zone around the perforation tunnel to enter the wellbore. The
perforation skin factor, Sp, will increase significantly where the
shot density is low, or effectively low because of perforations
Page 7
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
Page 8
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
0.33
15 0.5 7.9
5.3
12 0.6
SKIN FACTOR
3.4
o
0 PHASING
o
PRODUCTIVITY RATIO
6 90 PHASING 2.0 DAMAGED ZONE
DAMAGED ZONE THICKNESS (in.)
0.8
DIAMETER
3 0 0.9 0.9
d
b' - c'
0 3 o
1.0 0.0 OPEN HOLE
120
DIAMETER
0 o 180 o 90 o
6 1.1 -0.7
Ø = PHASE ANGLE
b-c 6" BOREHOLE CRUSHED
9 ANGULAR
1.2 -1.3
ZONE DIAMETER
0.4 0.1 (k d/ ku ) PERFORATION
PHASING
SPACING
o
8 spf, 0 PHASING (DEPENDENT ON
(k c / k u )
SHOT DENSITY)
CRUSHED ZONE PERFORATION
(0.1)
RELATIVE DIAMETER
PERMEABILITY
(k c / ku ) (0.2)
PERFORATION
LENGTH
(0.4)
(1.0)
4
2 1
Page 9
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
Table 4b.1
Page 10
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
PARTIAL
COMPLETION
SKIN
In situations where less than about 85% of the total net pay
thickness is open, flow will converge to enter the perforated
interval and cause a significant additional pressure drop. The
resultant skin, Sc, is variously called the partial completion
skin, partial penetration skin, geometric skin or pseudo-skin.
Page 11
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
28
24
20
10 0
1 000 0
Sc
00
16
100
30
50
12
h' h
rw rw
8
20
10
4
5
2
0 1 hp
0 0.2 0.4 0.6 0.8 1.0 b=
h
30ft
15ft
0.25ft 6ft
75ft
Examples of partial well completion showing; (a) Well only partially penetrating the formation;
(b) Well producing from only the central portion of the formation; (c) Well with 5 intervals open to production.
h p cos θ
=
h
Page 12
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
where: θ =
well angle measured from the normal to the
bedding plane (usually near vertical). This
is illustrated by the three examples of
Figure 4b.5.
b = 30/150 = 0.2 (for all 3 cases)
h/rw = 150/0.25 = 600 and Sc = 18 for (a)
h/rw = 75/0.25 = 300 and Sc = 16 for (b)
h/rw = 15/0.25 = 60 and Sc = 9.5 for (c)
(*Reference 4b.9) More rigorous analyses have been presented by Odeh* (1976)
(*Reference 4b.10) and Streltsova-Adams* which confirmed that the partial
completion skin was strongly dependent on the open interval
length, the ratio of horizontal to vertical permeability, other
skin effects and more weakly on the location within the pay.
0.825
S'c = 1.35 h –1 ln h k / k v + 7 – 1.95
h p cos θ
rwc = rw (ft)
(Equation 4b.5)
Page 13
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
Page 14
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
2.06 1.865
(Equation 4b.6) S"c = – θ ' – θ' log h k
41 56 100 rw kv
0
15°
-1
DEVIA 30°
TION
ANGL
E θ
-2
45°
-3
S c"
-4
60°
-5
-6
75°
-7 4
2 3
10 10 10 10
h
hD = rw
Page 15
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
Page 16
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
STANDARD FRACTURE
x f < 15% r e
S = In 2r w
xf
LARGE FRACTURE
xf
x f > 15% r e
FLOW IS LINEAR
Page 17
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
(Equation 4b.9) Ss ≈ 7 –7
J f/ J o
SKIN DUE TO
RESERVOIR
HETEROGENEITY
Many formations, especially carbonates, are naturally
fractured and will exhibit a negative skin once the drilling
damage is bypassed by acidization or by effective
perforations. A similar effect is seen when high permeability
streaks act as a feeder for the rest of the pay. An unexplained
negative skin is therefore generally attributed to reservoir
Page 18
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
Page 19
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
0.8 A
PP 25% DRAWDOWN S m = + 0.875
A
R
EN VO
T G 50% DRAWDOWN S m = + 2.0
P / P avg
0.6 LI EL
N
EA IP
R R IDE
IP (S AL
R o= LIN 75% DRAWDOWN S m = + 3.5
0.4 ,J EA
'(S 0. RI
0) PR
m ,J
= * (S
5.
6) o= 100% DRAWDOWN S m = + 5.6
0.2 0.0
)
0
0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
q / q max
Figure 4b.8 – Multiphase Flow Skin in Oil Wells Below the Bubble Point
Page 20
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
Page 21
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
NON-DARCY SKIN
+ 75 -3
D = 2.1 x 10 / (stb/d)
∆ SKIN
SLOPE = =D
∆ RATE
APPARENT SKIN, S'
+ 50
FLOWRATE, stb/d
Page 22
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
ONSET OF NON-
DARCY FLOW
(*Reference 4b.19) As early as 1937, Muskat* showed that deviations from Darcy
flow could be expected if the Reynold’s number exceeded 1.0.
(*Reference 4b.18) Golan* has presented a more practical limit as the point
where the high rate flow effects exceed 5% of the total
pressure drop, given by the following equations:
19 × 10 3 µo rw hp2 ln rr e + S
w
(Equation 4b.11) q o, HVF > stb/d
ρ o Bo h k / k a
1.4 × 10 3 µg rw hp2 ln rr e + S
w
(Equation 4b.12) q g, HVF > mscf/d
γg h k / k a
where: qo, HVF = critical oil rate for high velocity flow
effects, stb/d
qg ,HVF = critical gas rate for high velocity flow
effects, mscf/d
rw = wellbore radius, ft
hp = perforated interval open to flow, ft
h = total net pay thickness, ft
µ = gas or oil viscosity, cp
γg = reservoir gas gravity, (air = 1)
ρo = live oil density at flowing pressure, lb/ft3
Bo = oil FVF at flowing pressure, bbl/stb
k = formation permeability, md
ka = near wellbore permeability, md
S = Darcy skin factor
From this, it can be seen that gas wells, and oil wells
producing large volumes of free gas at the formation face,
will␣ quickly reach non-Darcy flow conditions. However,
non-Darcy flow effects do not become significant in oil wells
with normal GLRs, unless rates exceeds some 100 stb/d/ft,
or an internal gravel pack is installed (as discussed later).
NON-DARCY SKIN
IN GAS WELLS
The gas well productivity discussion in Section 4a – Reservoir
Performance presents the derivation and general form of the
non-Darcy skin component, D, in terms of the non-Darcy flow
coefficient, B.
Page 23
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
B kg h
(Equation 4a.24) D = m s c f / d- 1
1422 T
2.22 × 10 –15 γg k g h β
(Equation 4a.25) D = ×r m s c f / d- 1
h p2 µ w f w
2.22 × 10 –15 γg k g h β d r d – rw
(Equation 4b.13) Dd = × (mscf/d)–1
h p2 µ w f r d rw
2.22 × 10 –15 γg k g h r –r
(Equation 4b.14) Dr = × β r re r d (mscf/d)–1
h p µw f
2 e w
Page 24
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
NON-DARCY SKIN
IN OIL WELLS
A similar set of equations can be developed for the non-Darcy
skin components (D in [stb/d]–1) in a high rate oil well:
1.635 × 10 –16 ρo k oh Bo β d r d – r w
(Equation 4b.15) Dd = × r d rw (stb/d)–1
µo h p2
Page 25
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
Page 26
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
PERFORATIONAL TUNNEL
FORMATION SAND
FILTER SIZE
SAND
GRAVEL
Page 27
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
µo L Bo qo Bo q o 2
(Equation 4b.17) ∆p G = + 9.1 × 10 –13 βG L ρ o psi
1.127 × 10 –3 k G Ap Ap
(Equation 4b.19) SG = 2 π kh L
k G Ap
µo Ap
Four main parameters influence the magnitude of gravel pack
skin:
Page 28
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
π EH 2
(Equation 4b.22) Ap = h p Νp Ep ft 2
576
3.4 × 10 1 2
(Equation 4b.23) β G= (ft–1 )
k G1.54
Page 29
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
SKIN EFFECTS
FROM AN
INTERNAL
GRAVEL PACK
IN GAS WELLS
The pressure losses and skin effect from an internal gravel
pack in a gas well are strongly dependent on flowrate, since
high velocity flow effects will occur in the reservoir rock as
well as in the gravel filled perforations. Moreover the
computations are complicated by the dependence of gas
(*Reference 4b.22) properties on pressure. McLeod* shows that the pressure
drop across the gravel pack is:
8936 µ g T Z L q g
p 2wfs – p 2 w f =
k G Ap
2
qg
(Equation 4b.24) + 1.25 × 10 –10 γg Z T L β G psi 2
Ap
2 π kh 8.79 × 10 –14 γg L βG kh q g
(Equation 4b.25) S'G = +
k G Ap Ap2 µ g
so
γ g L βG kh
(Equation 4b.27) DG = 8.79 × 10 –14 m s c f / d –1
A µg
2
p
Page 30
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
Page 31
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
Page 32
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
Page 33
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4b – NEAR WELLBORE PERFORMANCE Completion Design Manual
REFERENCES
Page 34
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4b – NEAR WELLBORE PERFORMANCE
Page 35
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
CONTENTS
INTRODUCTION 1
COMPUTING ALGORITHM 37
TEMPERATURE PREDICTION 40
WELLBORE HEAT TRANSFER 41
TEMPERATURE PREDICTION 42
REFERENCES 50
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
INTRODUCTION
The accurate design of oil and gas well tubing strings requires
the ability to predict flow behaviour in the wells. All wells
produce a mixture of gas and liquids, regardless of whether
they are classified as oil wells or gas wells. This multiphase
flow is significantly more complex than single-phase flow.
However, the technology to predict multiphase flow behaviour
has improved dramatically in the past decade. It is now
possible to select tubing sizes, predict pressure drops and
calculate flowrates in wells with acceptable engineering
accuracy.
Page 1
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
FLOWING SEPARATOR
WELLHEAD
PRESSURE TO SALES
STOCK
TANK
INCLINED FLOW
HORIZONTAL FLOW
IN UBIN
VERTICAL
CL
T
IN G
INTAKE ED FLOW THROUGH
POROUS MEDIA
pR , k, IPR
Page 2
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Page 3
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Page 4
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Page 5
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
4000
OIL GAS-CONDENSATE GAS
RESERVOIRS RESERVOIRS RESERVOIRS
CRICONDENTHERM = 250 °F
A
Tc = 127 °F
3500 B
CRITICAL
POINT
TION
RESERVOIR PRESSURE, psia
3000
DUC
B1
DE
PRO
W
NT PO
POI
2500 IN
H OF
B LE T
UB
PAT
B
% B2
80
2000
%
40
ID %
QU ME 20
1500 I
L LU
VO
%
10
5%
1000
B3 0%
A2
500
0 50 100 150 200 250 300 350
RESERVOIR TEMPERATURE, °F
Page 6
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Once the black oil model parameters are known, oil density
and other physical properties of the two phases can be
calculated. Methods for predicting these properties are given
in the BP MULTIFLO Manual.
Page 7
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Volumetric
Flowrates
After mass transfer calculations are completed, it is possible
to calculate the in situ volumetric flowrates of each phase.
For the black oil model, volumetric flowrates are determined
from:
Page 8
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
p SC Z T
(Equation 4c.4) Bg =
p ZSC T SC
w t 1–x
(Equation 4c.5) qL =
ρL
qg = wt x
(Equation 4c.6)
ρg
DEFINITION
OF VARIABLES
When performing multiphase flow calculations for wells,
single-phase flow equations are often modified to account for
the presence of a second phase. This involves defining
mixture expressions for velocities and fluid properties that
use mixture weighting factors based on either volume or
Page 9
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Mixture
Weighting
Factors
When gas and liquid flow together up a well, the gas normally
travels faster than the liquid. Under steady state conditions
this results in a reduced area fraction for the gas phase and
an expanded area fraction for the liquid phase. Thus,
slippage of gas past liquid results in an accumulation of
liquid in the pipe, a phenomenon known as liquid hold-up.
Liquid hold-up can be defined as the fraction of a pipe cross
section or volume increment that is occupied by the liquid
phase. The difference in velocities of the gas and liquid
phases can be estimated from a variety of equations
developed from experimental data for a limited range of flow
conditions. A more common procedure is to develop empirical
correlations for predicting liquid hold-up for a broad range of
flow conditions.
q
(Equation 4c.7) f Lns = q +L q
L g
and where qL is the sum of the oil and water flowrates for the
black oil model, or is given by Equation 4c.5 for the
compositional model. If free water exists when using the
compositional model, the water flowrate must be added to the
oil or condensate flowrate to account for all the liquid. Since
the no-slip liquid hold-up can be determined rigorously, it is
commonly used as a correlating parameter for other
multiphase flow parameters, such as fL.
Page 10
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
q
(Equation 4c.8) f o = q +o q
o w
The water cut, fw, based on in situ rather than stock tank
flowrates, is simply 1–fo .
VMg
(Equation 4c.9) x =
VMg+ LML
Thus,
qL
(Equation 4c.10) vSL =
Ap
Page 11
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
qg
(Equation 4c.11) vSg =
Ap
qL + qg
(Equation 4c.12) vm = = vSL + vSg
Ap
vSg
(Equation 4c.14) vg =
1 – fL
(Equation 4c.15) vs = vg – vL
Page 12
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Mixture Fluid
Properties
Regardless of whether the black oil or the compositional
model is used, it is possible to calculate the density and
viscosity of each phase (oil, water and gas) and the gas-oil
and gas-water surface tensions. A variety of methods have
then been used to define mixture fluid properties.
Oil-Water Mixture
For the black oil model case, and for the compositional model
case when free water exists, oil and water properties are
combined as follows:
(Equation 4c.16) ρ L = ρ o fo + ρ w fw
(Equation 4c.17) σL = σo fo + σw fw
(Equation 4c.18) µL = µo fo + µw fw
Page 13
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
0.1
0.09
LAMINAR CRITICAL TRANSITION
FLOW ZONE ZONE FULLY ROUGH ZONE
0.08
0.05
0.07
0.04
0.06
0.03
LAM
0.05 0.02
RELATIVE ROUGHNESS, D
e
INAR 64
f = Re
0.015
FRICTION FACTOR, f
FLO
0.04
0.01
W
0.08
0.006
0.03 Re
cr 0.004
0.025
0.002
0.02 0.0001
0.0008
0.0006
0.0004
0.015
0.0002
SMOOTH
PIPES 0.0001
0.000,05
0.01
0.009
0.008 0.000,01
3 2 3 4 5 6 8 4 2 3 4 5 6 8 5 2 3 4 5 6 8 6 2 3 4 5 6 8 7 2 3 4 5 6 8 8
10 10 10 10 10 10
e e
= 0.000,001 = 0.000,005
pVD D D
REYNOLDS NUMBER, Re µ
Page 14
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
4
10
LAMINAR FLOW
REGIME LINE
3
10
OIL VISCOSITY, cp
2
10
10
1
0 0.2 0.4 0.6
INPUT WATER FRACTION REQUIRED
TO INVERT THE MIXTURE
Page 15
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Gas-Liquid Mixture
Numerous equations have been proposed for describing PVT
properties of gas-liquid mixtures. In general, these equations
are referred to as ‘slip’ or ‘no-slip’ properties, depending upon
whether fL or fLns is used as the volumetric weighting factor.
Thus, for the case of two-phase viscosity:
(Equation 4c.20) µs = µL fL + µg 1 – f L
(Equation 4c.21) µs = µL fL × µ g 1 – fL
Page 16
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
(Equation 4c.23) ρ s = ρ L fL + ρ g 1 – f L
(Equation 4c.25) f 2
1 – f Lns 2
ρk = ρL L ns
+ ρg
fL 1 – fL
(Equation 4c.26) h = h L 1 – x + h gx
Page 17
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
FLOW
PATTERNS
As a multiphase mixture flows up a well, the gas and liquid
can be distributed in a variety of configurations called flow
patterns. Predicting the flow pattern that occurs at a given
location in a well is extremely important since the empirical
correlation or mechanistic model used to predict pressure
drop varies with flow pattern. Numerous investigators have
described flow patterns in wells and made attempts to predict
when they occur. Consequently, a concensus of how to
classify flow patterns does not exist.
Flow Pattern
Classification
For upward multiphase flow of gas and liquid, most
investigators now recognize the existence of four flow
patterns:
• Bubble flow.
• Slug flow.
• Annular-slug transition (churn) flow.
• Mist (annular) flow.
• Bubble Flow
Bubble flow is characterized by a uniformly distributed
gas phase as discrete bubbles in a continuous liquid
phase. Based on the presence or absence of slippage
between the two phases, bubble flow is further classified
into bubbly and dispersed bubble flows. In bubbly flow,
relatively fewer and larger bubbles move faster than the
liquid phase due to slippage. On the other hand, in
dispersed bubble flow numerous tiny bubbles are
transported by the liquid phase, causing no relative
motion between the two phases. Dispersed bubble flow is
often referred to as froth flow.
Page 18
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
• Slug Flow
Slug flow is characterized by a series of slug units, each of
which is composed of a gas pocket called a Taylor bubble,
a plug of liquid called a slug and a film of liquid around
the Taylor bubble flowing downwards. For vertical flow,
the Taylor bubble is an axially symmetrical bullet-shaped
gas pocket that occupies almost the entire cross-sectional
area of the pipe. The liquid slug, carrying distributed gas
bubbles, bridges the pipe and thus separates two
consecutive Taylor bubbles.
Page 19
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Flow Pattern
Occurrence
The following is a description of a typical sequence of flow
patterns in an oil wellbore. Starting with a flowing bottom
hole pressure above the bubble point pressure, only a liquid
phase exists at the bottom. As the liquid moves upward, it
experiences pressure decline resulting in the liberation of
some of the gas dissolved in the liquid phase. The liberated
gas appears as small bubbles in the continuous liquid phase,
which characterizes the bubble flow pattern. As the flow
continues upward, further decrease in pressure and
temperature occurs, resulting in gas expansion and the
liberation of more solution gas from the oil phase. This
creates more and larger bubbles that start coalescing with
each other. The coalescence creates larger Taylor bubbles
separated by the continuous liquid phase. The slug flow
pattern thus occurs. Further continuation of the upward
motion of the flow into the region of lower pressure causes
the expansion of Taylor bubbles along with the liberation of
more gas from liquid slugs. This creates a chaotic two-phase
flow, defined earlier as churn flow. The churn flow continues
to exist until the gas flowrate is sufficiently high to push the
liquid against the pipe wall. This characterizes the existence
of annular flow.
Page 20
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
ANNULAR FLOW
(MIST FLOW)
CHURN FLOW
(ANNULAR-SLUG FLOW)
SLUG FLOW
BUBBLE FLOW
LIQUID FLOW
Page 21
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
2
10
10
FROTH FLOW
5
N
TIO
Ls
NSI
1
TRA
BUBBLE
5 FLOW
Lm
2
PLUG HEADING SLUG
-1 FLOW FLOW MIST 1
10
-1 2 5 2 5 2 5 2 2 5 3
10 1 10 10 10
N Gv
ρL 0.25 ρL 0.25
N Lv = 1.938 VSL N Gv = 1.938 VSG
σL σL
ρG
1 1
ρ σ /4 /3
L WA
Y =
ρ σ
X = ρA
W L
100
10
ft / sec
BUBBLE
1
YVSL
SLUG ANNULAR
0.1 MIST
FROTH
0.01
0.1 1 10 100 1000
XVSG ft / sec
Page 22
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
• Bubble-Dispersed Bubble
Turbulence due to high liquid velocities causes a
homogeneous, no-slip behaviour.
• Bubble-Slug
Increased gas superficial velocity results in the liquid
hold-up falling to a critical value of 0.75, below which gas
bubbles coalesce to form Taylor bubbles.
• Dispersed Bubble-Slug
Superficial gas velocity reaches a sufficiently high value
for no-slip liquid hold-up to have minimum possible value
for dispersed bubble flow of 0.48. Below this, gas bubbles
must coalesce regardless of turbulence.
• Annular Mist
Gas velocity is just high enough to drag largest liquid
droplets and the liquid film on the pipe wall upward.
20
10
DISPERSED BUBBLE
SUPERFICIAL LIQUID VELOCITY, m/s
C
B
1 BARNEA
BUBBLY TRANSITION
ANNULAR
0.1 D
A
SLUG OR CHURN
0.01
0.002
0.02 0.1 1 10 100
Page 23
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
LIQUID HOLD-UP
Liquid hold-up is the most critical variable in the calculation
of pressure drop in wells. This is because of its importance in
calculating mixture density in Equation 4c.23 for use in the
hydrostatic head component in the pressure gradient
equation.
DIMENSIONAL
ANALYSIS
The use of dimensional analysis to solve engineering
problems that cannot be solved analytically is a powerful tool
that has been used successfully in many engineering fields,
including single-phase fluid mechanics and heat transfer.
The increased complexity of multiphase flow has logically
resulted in even broader use of the same principles. When
performed properly, and with the availability of appropriate
experimental data, use of dimensionless groups to scale flow
behavior beyond the range of data used to develop empirical
correlations can be very successful. Failure to include all
pertinent variables will always result in empirical correlations
with limited accuracy.
Page 24
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
ρL 0.25
(Equation 4c.27) N Lv = 1.938 vSL
σL
ρL 0.25
(Equation 4c.28) N g v = 1.938 vSg
σL
ρL 0.5
(Equation 4c.29) N D = 120.872 D
σL
Page 25
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
CONSERVATION
PRINCIPLES
The basis for calculating pressure drop during fluid flow in
pipes is conservation of mass and linear momentum. This is
true for both single-phase and multiphase flow. Following is
the development of the pressure gradient equation for steady
state single-phase flow.
Conservation
of Mass
Conservation of mass simply means that for a given control
volume such as the segment of pipe shown in Figure 4c.10,
the mass in minus mass out must equal the mass
(*Reference 4c.2 Knudsen) accumulation*.
dL
dZ
dX
∂ρ ∂
(Equation 4c.31) + ρυ = 0
∂t ∂L
Page 26
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Conservation
of Linear
Momentum
Page 27
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Page 28
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
dp ρ n g sin θ f TP ρ n vm 2
(Equation 4c.38) = gc +
dL 2 gc D
Slip Flow
Methods
Slip flow methods recognize that the gas and liquid do not
travel at the same velocity, but are independent of flow
pattern. The most commonly used methods in this category
(*References 4.21, 4.20) were presented by Hagedorn and Brown* and Gray*. For this
method, the pressure gradient is obtained from the following
equation:
ρ s g sin θ f TP ρ TP vm 2
gc +
dp 2 gc D
(Equation 4c.39) =
dL 1 – EK
Page 29
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
ρ TP d vTP
g
(Equation 4c.40) EK = c dL
dp
dL
(*Reference 4.14) • If Duns and Ros* predict mist flow or the transition
region␣ to occur, then the Duns and Ros correlation should
be used.
Page 30
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Gray Correlation
(*Reference 4.20) Gray* developed a slip model based on Equation 4c.39 for use
in the API 14B Subsurface Controlled Subsurface Safety
Valve sizing computer program. Empirical correlations were
presented for predicting liquid hold-up and pseudo wall
roughness when mist flow occurs. A revised kinetic energy
term may also have been used.
Flow Pattern
Dependent
Methods
All of the more modern correlations for predicting pressure
gradients in wells recognize the importance of both flow
pattern and slippage. Empirical maps or correlations are first
used to predict the flow pattern that exists at a given location
in the tubing. Equation 4c.39 is then used to calculate
pressure gradient, with different liquid hold-up and friction
factor correlations being used for each flow pattern.
vm vSG ρn
(Equation 4c.41) EK = gc p
Page 31
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Page 32
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Orkiszewski Correlation
After concluding that none of the existing pressure gradient
correlations was sufficiently accurate, Orkiszewski combined
the Griffith bubble flow correlation, the Duns and Ros mist
flow correlation and a new development for slug flow into a
comprehensive pressure gradient prediction method. The slug
flow development combined the Griffith and Wallis Taylor
bubble rise velocity correlation with a new parameter called
the liquid distribution coefficient to obtain the slip mixture
density and the friction component. Empirical equations for
the liquid distribution coefficient, depending on the
Page 33
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
(*Reference 4.32) • Chierici et al* recorrelated the Griffith and Wallis Taylor
bubble rise velocity data in an attempt to improve the
Orkiszewski correlation. The statistical performance of the
(*Reference 4.24) result was worse than the original correlation*. (Not
available in MULTIFLO.)
Page 34
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Caetano (TUFFP)
(*Reference 4c.17) Caetano* conducted an experimental investigation of
multiphase flow in a 52.5 ft long vertical annulus consisting
of 3.0in ID casing and 1.66in OD tubing. Tests were
conducted with both air-water and air-kerosene mixtures and
with both concentric and fully eccentric annuli. Modifications
to the Taitel et al flow pattern prediction model were made for
flow in an annulus. Mechanistic models were formulated
based on experimental observations and on existing models
for specific flow patterns in pipe flow. The resulting
mechanistic model should be valid for vertical wells in which
production is through the casing-tubing annulus. Future
modifications may be necessary for deviated wells. This
model will shortly become available in MULTIFLO.
Ozon et al
In a joint study in Boussens, France by ELF Aquitaine, IFP
and Total-CFP, a comprehensive mechanistic model was
developed for predicting pressure drop in wells at any
inclination angle. Tests were conducted in 98 ft long pipes
with diameters of 3.0 and 6.0in using condensate, crude oil
Page 35
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Ansari (TUFFP)
(*Reference 4c.20)
Following an approach similar to Hasan and Kabir, Ansari*
developed a comprehensive mechanistic model that was
based primarily on combining the studies of Barnea* for flow
(References 4.18, 4c.14-
4c.16) pattern prediction and Ansari and Sylvester*, Alves et al *
and Sylvester* for pressure gradient prediction. The TUFFP
well databank was expanded to include high capacity wells
from the Prudhoe Bay Field. The resulting databank of 1775
flowing well surveys was used to evaluate the comprehensive
model in comparison to various empirical correlations. The
Ansari model is available in MULTIFLO.
Page 36
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
COMPUTING ALGORITHM
n
(Equation 4c.42) ∆p =
dp
dL
dL ≈ ∑ dp
dL i
∆ Li
i =1
0
where: L = length, ft
Page 37
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
START
GUESS p i +1
YES TEMP
PROFILE
KNOWN
NO
CALCULATE
Ti + 1 NO
ENTHALPY
BALANCE
YES
CALCULATE T i + 1
GUESS T i +1 BY COULTER &
BARDON EQUATION
CALCULATE
dq / dL & dh / dL
h i + 1 = h i ± dh
dL
. ∆L
FIND T i + 1 THAT
MATCHES h + FROM
h i + 1 = hL (1 - x) + hg x
NO CONVERGE
ON T i + 1
YES
CALCULATE dp / dL
p i + 1 = p i ± dp
dL
. ∆L
NO CONVERGE
ON p i + 1
YES
NEXT INCREMENT
Page 38
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Page 39
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
TEMPERATURE PREDICTION
HEAT EXCHANGER
U'1 dE Q U'2
dt
p1V1 p2V2
mv12 mv22
2 gc Z2 2 gc
PUMP OR TURBINE
mgZ 1 mgZ 2
Z1
gc gc
REFERENCE PLANE
W'
Page 40
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
p
dU = dh – 1 d
(Equation 4c.45)
J ρ
WELLBORE
HEAT TRANSFER
When hot reservoir fluids enter a wellbore, they immediately
begin losing heat to the cooler surrounding rock as they flow
to the surface. The surrounding rock gradually heats up,
reducing the temperature difference and thus the heat
transfer between the fluids and the rock. Eventually, for a
constant mass flowrate, the entire earth surrounding the well
reaches a steady state temperature distribution. Prediction
of the fluid temperature in the wellbore as a function of depth
(and time) is necessary to determine fluid physical properties
required to calculate pressure gradients.
A simple expression for the total heat loss from the fluids in
(*Reference 4c.21)
the tubing can be estimated from Newton’s Law of cooling*:
(Equation 4c.47) Q = 2π r t o ∆ L U ∆T
Page 41
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
dQ' U π dt o T f – Te
(Equation 4c.48) = w
dL
Page 42
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
T L , t = T R – gG L sin θ + gG A sin θ 1 – e –L / A +
(Equation 4c.49) T bh – T R e –L /A
wc p 1
(Equation 4c.50) A = .
2 π r t oU
g
T L , t = T R – gG L sin θ + A sin θ – gG e –L /A – 1 +
gcJ cp
dp –L /A
(Equation 4c.51) T bh – T R e –L /A – A η e – 1
dL
where h is the Joule-Thompson coefficient that describes
cooling (or heating) by expansion.
Page 43
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Page 44
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
–k
2R 1 R1 1 – X cb – X c + 1 X –e
c –1
–1/k b
k 1 + R1 X c
(Equation 4c.53) Xc =
R1
where: b = k–1
k
e = k+1
k
R1 = vS g 1
vS L 1
Figure 4c.13 shows a plot of the solution to Equation 4c.53 for
different values of specific heat ratio, k. It is clear that, as
single phase gas flow is approached at high values of vSg, the
accepted value of 0.50 to 0.55 for Xc is predicted. However, as
the in situ gas/liquid ratio declines, much lower pressure
ratios (higher pressure drops) are required to reach critical
flow. The region above the curves represents subcritical flow,
and region below the curve is critical flow.
Page 45
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
0.7
k si 04 k=1.04
0.6 k si 20 k=1.20
SUB-CRITICAL
0.4 ASHFORD
CRITICAL FLOW
0.3
0.2
0.1
0.0
0.001 0.01 0.1 1 10 100 1000 10 000 100 000
C p1
(Equation 4c.55) q gsc =
γg T 1
Page 46
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Table 4c.1
ρ n1 vmB1 2
(Equation 4c.56) ∆p =
2 gc CD2
Page 47
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Page 48
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
6
DEPTH x 1000 ft
12 000 ft
10
12
14 50
100
GAS/LIQUID RATIO, SCF/STB
16
200
18 300
400
20
30
20 0
15
12
10
80
60
50
0
00
00
00
00
3. From 13 840 ft on the vertical scale, proceed horizontally to the 400 scf/stb line
FIND THE FLOWING BOTTOM-HOLE PRESSURE. and read a flowing pressure of 3500 psi
Page 49
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
REFERENCES
Page 50
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4c – MULTIPHASE FLOW
Page 51
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4c – MULTIPHASE FLOW Completion Design Manual
Page 52
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
CONTENTS
INTRODUCTION 1
GAS LIFT 8
GAS LIFT PROCESS AND APPLICATIONS 8
IMPLICATIONS FOR COMPLETION DESIGN 12
COMMON PROBLEMS 13
PLUNGER LIFT 40
REFERENCES 42
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
INTRODUCTION
Page 1
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
Page 2
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
Page 3
SECTION 4d – ARTIFICIAL LIFT
Determining the Optimum Well Performance – Chapter 4
MAIN ARTIFICIAL LIFT SYSTEMS
Lift System Common Applications Normal Depth Rates Min BHFP Problems
Range (ft) (stb/d)
Continuous gas lift Offshore wells 12 000 10–50 000+ (0.05)–0.15 Gas supply
Deviated wells psi/ft Hydrates
Deep wells Single wells
High GLR wells Casing design
Moderately productive wells Depletion
Electrical submersible High rate, low GLR wells 12 000 500–20 000 100 psi Casing size
pumps (ESPs) Water supply wells (5000) (500–10 000) Handling
High WOR wells BHT >300°F
Alternative to high rate Sand
Gas lift High GLR
Table 4d.1
Hydraulic jet pump Mod-low GLR wells 12 000 100–8000+ 20% net lift Casing size
Page 4
Screw pumps Heavy oil wells (<18° API) 5000 1–500 Wellhead BHT >250°F
Shallow/light wells (1–1200) pressure High GLR
BP Exploration
Sandy wells (1–400)
Plunger lift Dewatering gas wells 12 000 1–500 100-500 psi Low GLR
High GLR/low rate oil wells Solids
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
(*Reference 4d.1) According to Brown*, the percentage use of these lift methods
in the USA on wells producing more than 10 stb/d is:
Page 5
SECTION 4d – ARTIFICIAL LIFT
Determining the Optimum Well Performance – Chapter 4
TYPICAL ARTIFICIAL LIFT PERFORMANCE SUMMARY
Service
Frequency Service
(yrs/job) Cost Rating of Ability to Handle Problems
High High
>250°F GLR Wax Sand Corrosion Scale Visc Doglegs
ESP
– <150°F 2
– 150-250°F 1 H 2 3 4 2 2 3 2 3
– >250°F 0.5
Hydraulic
– Jet 1.5 L 4 4 5 4 5 5 3 4
– Piston 0.5 M 4 2 4 2 4 2 5 4
Plunger 1 VL 4 5 4 1 2 1 1 5
BP Exploration
1 – Poor M – Medium VL – Very Low
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
Page 7
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
GAS LIFT
GAS LIFT
PROCESS AND
APPLICATIONS
Continuous gas lift simply supplements the natural flow
process by adding additional gas into the produced fluid to
reduce the hydrostatic head component and hence, the back
pressure on the formation.
Page 8
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
The lift gas is normally injected down the annulus and into
the tubing through a gas lift valve (GLV), which contains a
non-return check. In some wells, the gas is supplied via the
tubing with production being taken either up the annulus, or
through a second tubing string (concentric or parallel to the
supply string).
Page 9
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
As shown in the figure, gas passes through the first GLV and
lightens the fluid column to surface, according to the total
gas liquid ratio (TGLR) curve shown. As the fluid gradient in
the tubing changes, the gas in the casing moves down to the
second GLV, unloading kill fluid from casing. As the gas
reaches the second GLV and begins to lighten the fluid
column, the first GLV closes, so that all gas is passing
through the second GLV.
Page 10
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
Since the GLR requirements are very much subject to the law
of diminishing returns, most gas lift systems are designed
either on the basis of a near optimum GLR (NOGLR) that
achieves a BHFP within 20 to 50 psi of the minimum, or on
the basis of available injection gas supply volumes, Qi, (see
Figure 4d.6). It should be apparent that:
3) IGLR = Qi/q
Page 11
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
IMPLICATIONS
FOR
COMPLETION
DESIGN
Higher lift gas supply pressures are now commonly used to
reduce or eliminate the requirement for unloading valves to
kick-off wells. However, this greatly increases the risk of gas
leakage through the tubing and casing couplings. Couplings
with improved gas sealing capacity are therefore generally
preferred. Because of the large spiral void area, buttress
threaded couplings (BTC) commonly give problems in high
pressure gas lift operations and should be avoided.
Page 12
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
COMMON
PROBLEMS
The most common problem with gas lift design is spacing the
valves too far apart because of underestimating the pressure
losses in the gas supply system and the back-pressure
caused by slugging. In operation, the most common problem
is gas leakage through the upper valves, and/or through
tubing or casing couplings.
Page 13
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
ESP EQUIPMENT
AND
APPLICATIONS
ESPs have their greatest application in moving large volumes
(>500 stb/d) of low GLR fluid (<100 scf/stb). They are
particularly attractive for water supply wells, high water cut
producers and high deliverability undersaturated oil wells.
However, units are available at pumping rates as low as
90␣ stb/d. Where pump intake pressure is below the bubble
point, modified design procedures and improved downhole
gas separators/compressors allow effective operation at up to
1000 scf/stb.
Page 14
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
(Equation 4d.3)
∆pp = 0.433 (SG) Ns Hs
Page 15
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
It can be seen from the above that the ESP delivery capacity
will vary with:
Page 16
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
Page 17
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
PUMP AND
SYSTEM
PERFORMANCE
It is usual to select the largest pump that will fit into the
production casing or to size the production casing to
accommodate the pump. Small casing or liners will limit the
maximum feasible horsepower and pump capacity (see
Table␣ 4d.3 ).
Page 18
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
Table 4d.3
Page 19
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
Page 20
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
Page 21
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
The power supply cable has low inherent strength and must
be properly supported by the tubing. Two cable clamps per
joint are recommended. The material used for the clamps
must be compatible with the operating environment and
armour material. Centralizers or protective cable clamps
should be used to provide crush protection through doglegs
and build-up sections.
Page 22
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
COMMON
PROBLEMS
The main problem with ESP operations is short run times
and consequently high service costs and production losses.
These are commonly caused by:
Page 23
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
Page 24
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
BASIC
CONCEPTS AND
APPLICATIONS
Since it is not particularly sensitive to temperature, depth,
deviation or severe operating environments, hydraulic
transmission has obvious attractions as an alternative to
electrical or mechanical means of transferring energy
downhole. The major problems are, however, the need for at
least two reasonably large conduits to minimize fluid friction
losses, the difficulty of maintaining a clean, solids-free power
fluid, and the high capital and maintenance costs for high
pressure pump units.
The two main systems, the jet pump and piston pump, are
interchangeable in most installations, providing a high degree
of flexibility to cope with changing well conditions. Moreover,
in many cases, the bottom hole pump assembly can be
circulated in or out of the well using a set of swab cups as
a␣ pump-down motor. Alternatively, the units can be run and
pulled with wireline. This permits low downhole service costs
and minimal downtime for pump servicing or replacement.
Page 25
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
Page 26
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
The annulus may also be used for gas venting, especially for
the piston pumps, which are sensitive to gas locking.
Occasionally, in low rate installations, a third conduit is
installed to permit the use of a closed power fluid system.
Page 27
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
0.45 0.5
0.25 1.0
0.17 1.5
0.10 2.0
Table 4d.4
Page 28
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
Page 29
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
IMPLICATIONS
FOR
COMPLETION
DESIGN If multiple strings are to be used in a hydraulic pump
installation, the casing size is obviously a critical limitation.
Concentric completions are generally favoured because they
are more space efficient and provide a large return conduit
size. If the annulus is used for the produced fluid, special
attention must be paid to corrosion protection, materials
selection, and de-oxygenization of the power fluid.
Page 30
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
BASIC
EQUIPMENT AND
APPLICATIONS
On land, the most common artificial lift system for low rate
wells is rod or beam pumping. Typically this is used for
stripper wells (<10 stb/d) and for shallow wells (<8000 ft)
producing less than 500 stb/d, although rod pumps are
sometimes used to produce up to 1500 stb/d.
The annulus is normally left open and used to vent any free
gas that can be separated downhole. Tubing is used as the
production conduit to protect the casing from wear and
corrosion. The tubing is generally anchored to the casing and
pulled into tension to minimize tubing movement, buckling
and rod wear.
(Equation 4d.4)
PD = Ct x Sp x N x D2 x Ep
Page 31
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
Page 32
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
(Equation 4d.5) Fo = Ct x SG x D2 x H
Page 33
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
The surface loads on the polished rod will include the fluid
loads, the rod weight, and the dynamic forces. Rod loads are
significant in deeper wells. A tapered rod string is therefore
often employed which involves using thicker rods on the
upper sections, where loads are higher. Buoyancy will vary
throughout the stroke, but is normally calculated on the
downstroke when the travelling valve is open.
Page 34
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
The API design method is not reliable for heavy oil wells
(<20°␣ API), unless proper corrections are applied for the fluid
friction load caused by the high viscosities and the lack of rod
weight on the downstroke. Rod fall problems generally
restrict pump rates to 1.5 to 2.5 SPM and favour the use of
long stroke, variable speed pumping units, such as the
hydraulic pump jacks. The high viscosities also lead to sand
production problems, poor valve performance, poor gas
separation and poor lubrication at standard clearances
(0.003in).
Page 35
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
Figure 4d.17 – API Bulletin 11L4, Composite Design Curves for
Rod Pump Installations
Page 36
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
IMPLICATIONS
FOR
COMPLETION
DESIGN
If free gas will be present downhole at the intended BHFP, it
is highly desirable to have an open annulus through which
gas can be vented. Therefore, packers should not be installed
in wells that can be expected to be converted to rod pump
within a short period of time (<2 years), unless they are
essential for zonal isolation or safety reasons.
Page 37
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
Page 38
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
Page 39
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
PLUNGER LIFT
Plunger lift is a lift system for high GLR wells that produce
liquids at relatively low rates (< 500 stb/d). It requires an
area where gas energy can be stored and then supplied to the
tubing at very high rates. Usually the tubing-casing annulus
is used for this purpose, but a natural or induced fracture
system may be adequate. The gas energy is used to drive a
piston, or plunger, carrying a small slug of liquid up to the
surface. After production of the tail gas, the well is closed-in
and the plunger falls back to bottom (see Figure 4d.19).
Page 40
BP Exploration Determining the Optimum Well Performance – Chapter 4
Completion Design Manual SECTION 4d – ARTIFICIAL LIFT
Page 41
Determining the Optimum Well Performance – Chapter 4 BP Exploration
SECTION 4d – ARTIFICIAL LIFT Completion Design Manual
REFERENCES
Page 42