Section 2 Multiphase Flow, Flowing Well Performance

Similar documents
PRODUCTION I (PP 414) ANALYSIS OF MAXIMUM STABLE RATE AND CHOKES RESOLVE OF OIL WELLS JOSE RODRIGUEZ CRUZADO JOHAN CHAVEZ BERNAL

Gas Lift Workshop Doha Qatar 4-88 February Gas Lift Optimisation of Long Horizontal Wells. by Juan Carlos Mantecon

UNIT 2 PATTERNS OF NATURAL

EFFECT OF GAS-OIL-RATIO ON OIL PRODUCTION MASTER OF SCIENCE

Acoustic Liquid-Level Determination of Liquid Loading in Gas Wells

Using the Artificial Gas Lift to Increase the Productivity of Noor Oil Field / Mishrif Formation

New power in production logging

COPYRIGHT. Production Logging Flowmeter Survey Solution Guide

The water supply for a hydroelectric plant is a reservoir with a large surface area. An outlet pipe takes the water to a turbine.

SURFACE CASING SELECTION FOR COLLAPSE, BURST AND AXIAL DESIGN FACTOR LOADS EXERCISE

A Successful Experience in Optimization of a Production Well in a Southern Iranian Oil Field

KICK CIRCULATION ANALYSIS FOR EXTENDED-REACH AND HORIZONTAL WELLS. A Thesis MAXIMILIAN M. LONG

Numerical Analysis of Two Phase Flow Patterns in Vertical and Horizontal Pipes

Simulation of Gas Holdup in a Bubble Column with a Draft Tube for Gas Dispersion into an Annulus

Optimization of Separator Train in Oil Industry

Lecture Outline Chapter 15. Physics, 4 th Edition James S. Walker. Copyright 2010 Pearson Education, Inc.

Environmental Science: An Indian Journal

Air Bubble Departure on a Superhydrophobic Surface

Bubble Dynamics in a Vibrating Liquid. By: James Wymer, Jaggar Henzerling, Aaron Kilgallon, Michael McIntire, Mohammed Ghallab

DETAILED ANALYSIS OF PRESSURE DROP IN A LARGE DIAMETER VERTICAL PIPE THESISPRESENTED TO THE DEPARTMENT OF PETROLEUM ENGINEERING

W I L D W E L L C O N T R O L PRESSURE BASICS AND CONCEPTS

Critical Velocity/ Nodal Stability Gas Well Predictions

OLGA. The Dynamic Three Phase Flow Simulator. Input. Output. Mass transfer Momentum transfer Energy transfer. 9 Conservation equations

and its weight (in newtons) when located on a planet with an acceleration of gravity equal to 4.0 ft/s 2.

Using Gas Lift to Unload Horizontal Gas Wells

Bioreactor System ERT 314. Sidang /2011

In this third article of a three-part series, the results of

A COMPARATIVE ANALYSIS OF CONTINUOUS GAS LIFT VALVE POSITIONING METHODS AND ITS APPLICATION WITH MS EXCEL & VBA

Flow Scanner. Production logging in multiphase horizontal wells

Downhole Gas Separator Selection

Chapter 15 Fluid. Density

4-3 Rate of Change and Slope. Warm Up. 1. Find the x- and y-intercepts of 2x 5y = 20. Describe the correlation shown by the scatter plot. 2.

Permeability. Darcy's Law

Study on the Influencing Factors of Gas Mixing Length in Nitrogen Displacement of Gas Pipeline Kun Huang 1,a Yan Xian 2,b Kunrong Shen 3,c

Flow and Mixing in the Liquid between Bubbles

Design Tapered Electric Submersible Pumps For Gassy Wells Desheng Zhou, SPE, Rajesh Sachdeva, SPE, IHS INC.

Homework of chapter (3)

Numerical Simulations of a Train of Air Bubbles Rising Through Stagnant Water

ANNUAL REPORT UIUC, August 16, Bubble Formation, Breakup and Coalescence in Stopper-rod Nozzle Flow and Effect on Multiphase Mold Flow

EDUCTOR. principle of operation

Beam Pumping to Dewater Horizontal Gas Wells

Two Phase Fluid Flow (ENGINEERING DESIGN GUIDELINE)

Tool Weight Analysis An Engineering Approach to Optimize Pressure Temperature and Spinner Data Acquisition

Hydro-Mech Bridge Plug

Flow transients in multiphase pipelines

GAS KICK MECHANISTIC MODEL. A Thesis RAHEEL ZUBAIRY

COPYRIGHT. Reservoir Fluid Fundamentals. Reservoir Brine Basic Workflow Exercise Review. Brine Density at Standard Conditions

Downhole Gas Separator Selection

The Mechanism Study of Vortex Tools Drainage Gas Recovery of Gas Well

CVEN 311 Fluid Dynamics Fall Semester 2011 Dr. Kelly Brumbelow, Texas A&M University. Final Exam

LOW PRESSURE EFFUSION OF GASES revised by Igor Bolotin 03/05/12

Production Optimization Using Bean Size Selection for Niger Delta Oil Wells

Scott Denning CSU CMMAP 1

GAS CONDENSATE RESERVOIRS. Dr. Helmy Sayyouh Petroleum Engineering Cairo University

Plunger Fall Velocity Model

Plunger Fall Velocity Model

Air Pollution Dispersion

Continuous Flow Gas Lift Design By Jack R. Blann 2/9/2008 1

Old-Exam.Questions-Ch-14 T072 T071

Example: Calculate the density of methane at 50 psig and 32 ⁰F. H.W. In previous example calculate the density of methane in gm/m 3.

Energy and mass transfer in gas-liquid reactors.

Petr 311 Wellbore pressures and Kick Control

Well Control Modeling Software Comparisons with Single Bubble Techniques in a Vertical Well

Flow in a shock tube

B04 Guided Interpretation #4

TUTORIAL. NPSHA for those who hate that stuffy word. by Jacques Chaurette p. eng. copyright 2006

G.C.E (A/L) Examination March In collaboration with

Fully coupled modelling of complex sand controlled completions. Michael Byrne

Static Fluids. **All simulations and videos required for this package can be found on my website, here:

PMI Pulse Decay Permeameter for Shale Rock Characterization Yang Yu, Scientist Porous Materials Inc., 20 Dutch Mill Road, Ithaca NY 14850

Advanced Intermittent Gas Lift Utilizing a Pilot Valve

STUDY OF SLUG CONTROL TECHNIQUES IN PIPELINE SYSTEMS

Irrigation &Hydraulics Department lb / ft to kg/lit.

Appalachian Basin Gas Well. Marietta College, Marietta, Ohio June 7-8, Development. Robert McKee, P.E. Design Engineer Multi Products Company

Robust Autonomous Flow Control Developments in Norway

L 15 Fluids [4] The Venturi Meter. Bernoulli s principle WIND. Basic principles of fluid motion. Why does a roof blow off in high winds?

COPYRIGHT. Reservoir Fluid Core. Single Phase, Single Component Systems. By the end of this lesson, you will be able to:

Fluid Mechanics - Hydrostatics. Sections 11 5 and 6

Assumptions 1 At specified conditions, air behaves as an ideal gas. 2 The volume of the tire remains constant.

Casing Design. Casing Design. By Dr. Khaled El-shreef

32 nd Gas-Lift Workshop. Optimizing Subsea Well Kick-off Operations Case Study Using FlowLift2

From the Reservoir Limit to Pipeline Flow: How Hydrocarbon Reserves are Produced. 11/10/

Well Performance and Nodal TM Analysis Fundamentals COPYRIGHT. Skill Module Overview

Understanding pressure and pressure

Workshop 1: Bubbly Flow in a Rectangular Bubble Column. Multiphase Flow Modeling In ANSYS CFX Release ANSYS, Inc. WS1-1 Release 14.

LOW PRESSURE EFFUSION OF GASES adapted by Luke Hanley and Mike Trenary

SPE Copyright 2012, Society of Petroleum Engineers

Write important assumptions used in derivation of Bernoulli s equation. Apart from an airplane wing, give an example based on Bernoulli s principle

FLOWING GRADIENTS IN LIQUID LOADED GAS WELLS

Numerical simulation of an intermediate sized bubble rising in a vertical pipe

Gas Lift Optimization Using Smart Gas Lift Valve

Chapter 13 Fluids. Copyright 2009 Pearson Education, Inc.

Intermittent Gas Lift Utilizing a Pilot Valve

Multiphase Flow Prof. Gargi Das Department of Chemical Engineering Indian Institute of Technology, Kharagpur

4 RESERVOIR ENGINEERING

MODELING AND SIMULATION OF VALVE COEFFICIENTS AND CAVITATION CHARACTERISTICS IN A BALL VALVE

3.3 USING A SIMPLE PARCEL MODEL TO INVESTIGATE THE HAINES INDEX

Two interconnected rubber balloons as a demonstration showing the effect of surface tension

Dynamic Surface Tensiometer Bubble Pressure Method BP-D5, BP-D5L. Kyowa Interface Science Co., Ltd.

Distribution of the Sizes of Rock Cuttings in Gas Drilling At Various Depths

Transcription:

Section 2 Multiphase Flow, Flowing Well Performance Multiphase Vertical Flow When producing an oil or gas well, the flow of the fluids up the tubing will be in most cases be 2 phase, liquid and gas. The pressure at any point in the tubing string is influenced by many factors. Such as type of flow at that point, the viscosity s of the oil and gas, the PVT characteristics of the fluid which are pressure and temperature dependent, the bottom hole and tubing pressures. In the tubing there is a pressure drop and a temperature variation which causes several of the fluid properties to vary along with the pressure drop that allows the solution gas to come out of the oil at different rates in various parts of the tubing, causing holdup behavior. The drag of the flowing fluids (friction) also plays a part along with gravity in resisting the flow. Flow Regimes As the oil travels up the tubing different types of flow can occur as the gas is liberated from the liquid phase. If the well is producing above the bubble point pressure the flow is single phase but as it rises up the tubing the pressure will drop below this pressure and the multiphase regimes will start to form. Bubble Flow The gas starts to come out of solution and form in the liquid medium in an even dispersion.

Slug or Plug Flow As the fluid moves upward the gas bubbles move faster than the liquid, because of the buoyancy effect, and collect in large bubbles. These bubbles grow to a size where extend across the diameter of the tubing separating slugs of liquid containing the smaller bubbles Churn Flow With increasing gas velocity, the larger bubbles become unstable and collapse, resulting in a highly turbulent flow pattern. Churn flow is characterized by oscillatory, up and down motions of the liquid. Annular Flow As the pressure is lowered these large pockets of case break through the slugs and form a continuos phase in the center of the tubing carrying droplets of the liquid with a film of oil flowing along the walls of the tubing.

Mist Flow Eventually the gas volume is so large that the film of liquid on the wall of the tubing disappears and the only liquid moving is the droplets or mist in the gas phase. Not all of these regimes are not necessarily present in a well but usually more than one is present. As can be seen the flow of oil and gas up a string of tubing has many complex components that must be considered in designing a flowing well. For a given bottom hole pressure there is a certain amount of reservoir fluids that will enter the well bore, what is needed to be known is whether the pressure drop in the tubing will allow flow to the surface. Several methods have been developed to describe this flow so pressure curves of the tubing can be drawn to use in the planning of the well and its equipment.

Poettmann & Carpenter Method Based on the energy equation with the assumptions that the difference in the kinetic energy at the top and bottom is negligible and that external work done by the flowing fluid is also negligible. The energy loss W f is caused by the slippage and frictional effects. 2 Δh W f = 4 fu (2-1) 2gD where u is the average mixture velocity over Δh and D is the inside diameter of the tubing, g is gravity constant and f is found empirically. With the given assumptions and field units equation 2-1 can be reduced to 144 Δ p Δh _ = ρ+ K _ (2-2) ρ where Δp = pressure drop over Δh, psi ρ = average density of the fluid, lb/cu ft _ K = 2 2 fq M 10 5 ( 7. 413x10 D ) (2-3) q = liquid production rate, STB/day M = Total mass of gas and liquid associated with 1 STB, lb D = ID of tubing, ft f = energy loss factor defined by equation 2-1 The way that this is used is to start at the bottom of the tubing with the known p wf which can produce a certain production of reservoir fluids that is known from the productivity index of the well. All the PVT and reservoir parameters such as FVF, water-cut, gas solubility are known or can be estimated. Divide the tubing into equal parts, H 1, H 2, H 3, H 4, each a length of Δh and let the pressures, p 1, p 2, p 3, p 4, at the H n points. At H 1 the pressure is p wf, so this pressure can be used for H 1 H 2 for the first approximation for the factors needed for equation 2-2, such as the density and Mass. If a better

approximation is needed use the average between of p wf and the first approximation for the next calculations. This process is repeated for all the intervals in the tubing, to the surface or till a negative pressure is calculated. Doing this for several starting p wf a plot can be generated for any size tubing. Figure 2-1 This graph can be used to find the bottom hole pressure and production rates for any producing tubing head pressure. The maximum rate that the well can flow is indicated by point E on the graph, this point must always be lower than q since there is a positive pressure loss in the tubing.

Figure 2-2 Plotting the tubing head pressures for several sizes of tubing is used to find the optimal tubing size for different flow rates. In this case at 200 psi THP the 2 tubing is better but at 100 psi the 3 ½ tubing produces at a higher rate. The only factor in equation 2-3 that is not easily found is the energy loss factor f. In an effort to establish values for f as a function of the variables involved in the flow of oil wells a reverse calculation was made. This involved taking field data in the form of production, PVT data, and pressure and temperature transverses. It was found that f correlated best as a function of tubing ID and mass velocity of the fluid in the tubing. In terms of field units 5 Mq D v ρ = (1.4737x10 ) (2-4) D From this the following plot is used to find f for different flow rates and tubing sizes.

Figure 2-3 It is noted here that this plotted on log-log paper and that the scatter of the points is greater that its appears on this graph. Also Baxendell extended this work to in include flow in an annulus, and with Thomas plotted f over higher flow rates. The weakness of this approach is the inaccuracy of finding f. Work by others have been to overcome this problem.

Using the Poettman & Carpenter Method Calculating M, mass of the fluid in STB M = wt of Oil + wt of Water + wt of Gas in lb m /STB wt of Oil = Oil cut x γ o x 350 lb/bbl wt of Water = Water cut x γ w x 350 lb/bbl wt of gas = GLR cf/bbl x γ g x.0764 lb/cf w = Mass flow rate, = Mq, lb m / day 141.5 Note: γ o = γ o w and γ g and GOR should be given & 131.5 + API the oil and water cut calculated from the WOR. Calculating ave. Density, ρ, lb m /cf ρ = M/volume Volume = vol of oil + vol of water + vol of free gas in cf Using the pressure of interval being calculated, lab or correlations of PVT data needed. β o, β g, R s Volume bbls = Oil cut x β o + Water cut + β g { GLR Oil cut x R s } Volume cf = Volume bbls x 5.614 cf/bbl ρ = M/(5.614[oil cut x β o + water cut] + β g {GLR Oil cut x R s }) With values of M and ρ the Δp/Δh can be calculated using equations 2-2 and 2-3 by finding f from figure 2-3. By dividing the pressure increment by

this pressure gradient gives the length of tubing over which the increment occurs. This way a plot of the pressure can be obtained. Example: Given Tubing 2 q o. = 400 bpd 22 o API oil q w.= 600 bpd γ w = 1.07 GLR = 500 cf/bbl γ g =.65 γ o = 141.5/131.5 + 22 =.922 M =.4*.922*350 +.6*1.07*350 + 500*.65*.0764 = 378.6 lb/bbl w = M*q = 378.6*1000 = 378610 lb/day D ρv = 1.4737E-5 x w/d = 1.4737E-5 * 378610/1.995/12 = 33.56 Form fig. 2-3 f =.008 K = f w 2 /(7.413E 10 D 5 ) =.008*378610^2 / (7.413E+10*.166^5) = 122.7 For pressure point 1, 500 psi βo = 1.043 Rs = 59 scf/bbl z =.93 T = 105F ρ = M/5.614[oil cut x β o + water cut] + β g {GLR Oil cut x R s } 378/(5.614[.4*1.043 +.6]+.0298{500 -.4*59}) = 19 lb/cf For pressure point 2, 1000 psi βo = 1.083 Rs = 120 scf/bbl z =.9 T = 145F ρ = 378.6/(5.614[.4*1.083 +.6] +.0153{500 -.4*120}) = 29.8 lb/cf ave ρ = ρ 1 + ρ 2 / 2 = 24.4 lb/cf Use this average density to find Δp/Δh in equation 2-2

Δp/Δh = [ρ + K/ρ]/144 = [24.4 + 122.7/24.4]/144 =.204 psi/ft Now solve for /Δh for the Δp Δh = 1000 500 /.204 = 2447 ft. Plot this point, say THP is 500 psi then it 1000 psi at a depth of 2447. Then repeat the process for the next interval. A large interval was used here in the example, the interval used when doing a pressure profile of a well the pressure intervals will be much smaller, like about 25 psi. Homework Given Tubing 2.5 q o. = 600 bpd 28 o API oil q w.= 400 bpd γ w = 1.07 GLR = 750 cf/bbl γ g =.65 For pressure point 1, 500 psi βo = 1.040 Rs = 60 scf/bbl z =.93 T = 105F For pressure point 2, 1000 psi βo = 1.080 Rs = 125 scf/bbl z =.9 T = 145F