Lecture 1 Choke Performance
Lecture 1 Choke Performance
Lecture 1 Choke Performance
Engineering Faculty
Petroleum and Mining Department
Wellhead chokes are used to limit production rates for regulations, protect
surface equipment from slugging, avoid sand problems due to high
drawdown, and control flow rate to avoid water or gas coning.
Placing a choke at the wellhead means fixing the wellhead pressure and, thus,
the flowing bottom-hole pressure and production rate. For a given wellhead
pressure, by calculating pressure loss in the tubing the flowing bottomhole
pressure can be determined. If the reservoir pressure and productivity index
is known, the flow rate can then be determined on the basis of inflow
performance relationship (IPR).
FRZAN ALI 2
1.2 Sonic and Subsonic Flow
Both sound wave and pressure wave are mechanical waves. When the fluid flow velocity in
a choke reaches the traveling velocity of sound in the fluid under the in situ condition, the
flow is called ‘‘sonic flow.’’
Under sonic flow conditions, the pressure wave downstream of the choke cannot go
upstream through the choke because the medium (fluid) is traveling in the opposite direction
at the same velocity.
Therefore, a pressure discontinuity exists at the choke, that is, the downstream pressure does
not affect the upstream pressure. Because of the pressure discontinuity at the choke, any
change in the downstream pressure cannot be detected from the upstream pressure gauge. Of
course, any change in the upstream pressure cannot be detected from the downstream
pressure gauge either.
This sonic flow provides a unique choke feature that stabilizes well production rate and
separation operation conditions.
FRZAN ALI 3
1.2 Sonic and Subsonic Flow
FRZAN ALI 4
Figure 1-1: A typical choke performance curve
FRZAN ALI 5
1.3 Single-Phase Liquid Flow
2 g c P
q = CD A (1.2)
where
p
q = 8,074C D d 2
(1.3)
2
where
q = flow rate, bbl/d
d2 = choke diameter, in.
p = pressure drop, psi
FRZAN ALI 7
The following correlation has been found to give reasonable
accuracy for Reynolds numbers between 104 and 106 for nozzle-
type chokes (Guo and Ghalambor, 2005):
d 2 0.3167
CD = + 0 .6
+ 0.025[log( N Re ) − 4]
d1 d (1.4)
2
d1
where
d1 = upstream pipe diameter, inch
d2 = choke diameter, inch
NRe= Reynolds number based on d2.
FRZAN ALI 8
1.2
d 2 /d 1
1.15 0.75
0.725
0.7
1.1 0.675
0.65
0.625
CD
1.05
0.6
0.575
0.55
1
0.5
0.45
0.95 0.4
0.9
1,000 10,000 100,000 1,000,000 10,000,000
Reynolds Number
0.45
0.65 0.4
0.3
0.2
0.6
0.55
1,000 10,000 100,000 1,000,000 10,000,000
Reynolds Number
FRZAN ALI 11
The Reynolds number for determining CD is
expressed as:
20q sc g
N Re = (1.6)
d 2
FRZAN ALI 12
Gas velocity under subsonic flow conditions is less than
the sound velocity in the gas at the in-situ conditions:
k −1
z up p down
k
v = vup + 2 g c C p Tup 1 −
2
(1.7)
z dn pup
FRZAN ALI 13
1.4.2 Sonic Flow
Under sonic flow conditions the gas passage rate reaches its
maximum value. Gas passage rate is expressed in the
following equation for ideal gases:
k +1
k 2 k −1
Qsc = 879C D Apup (1.8)
T k + 1
g up
FRZAN ALI 14
Gas velocity under sonic flow conditions is equal to
sound velocity in the gas under the in-situ
conditions:
The choke flow coefficient CD is not sensitive to the Reynolds number for Reynolds
number values greater than 10 6 . Thus, the CD value at the Reynolds number of 10 6 can be
assumed for CD values at higher Reynolds numbers.
z up 2
v = vup + 2 g c C p Tup 1 −
2
(1.9)
z outlet k + 1
or
FRZAN ALI 15
1.4.3 Temperature at Choke
Assuming an isentropic process for an ideal gas flowing through chokes, the
temperature at the choke downstream can be predicted using the following
equation: k −1
z up poutlet k
Tdn = Tup
z outlet pup (1.11)
The outlet pressure is equal to the downstream pressure in subsonic flow
conditions.
FRZAN ALI 16
1.4.4 Applications
• Downstream temperature;
FRZAN ALI 17
To estimate gas passage rate at given upstream and
downstream pressures, the following procedure can be taken:
FRZAN ALI 18
Example Problem 1.1:
FRZAN ALI 19
Solution: (a)
k 1.3
Poutlet 2 k −1 2 1.3−1
= = = 0.5459
P
up c k + 1 1.3 + 1
Pdn 200
= = 0.25 < 0.5459
Pup 800
k +1
k 2 k −1
q sc = 879C D APup
T k + 1
g up
.
1.3+1
1 .3 2 1.3−1
q sc = (879)(0.62)[ (1) / 4](800)
2
(0.6)(75 + 460) 1.3 + 1
FRZAN ALI 21
Check NRe:
= 0.01245
cp by the Carr-Kobayashi-Burrows correlation.
FRZAN ALI 22
(b)
k −1
z up Poutlet k 1.3−1
Tdn = Tup = (75 + 460)(1)(0.5459) = 465
z outlet Pup
1.3
oR= 5 oF < 32 oF
Therefore, heating is needed to prevent icing.
(c)
Poutlet
Poutlet = Pup = (800)(0.5459) = 437psia
P
up
FRZAN ALI 23
Example Problem 1.2:
A 0.65 specific gravity natural gas flows from a 2-in pipe through a
1.1-in nozzle-type choke. The upstream pressure and temperature
are 100 psia and 70 oF, respectively. The down stream pressure is
80 psia (measured 2 ft from the nozzle). The gas specific heat ratio
is 1.21.
FRZAN ALI 24
Solution: (a)
k 1.25
Poutlet
= 2 2
k −1 1.25−1
= = 0.5549
P
up c k + 1 1.25 + 1
Pdn 80
= = 0 .8 > 0.5549
Pup 100
Subsonic flow exists.
d 2 1.5"
= = 0.75
d1 2"
FRZAN ALI 25
Assuming NRe > 106, Figure 1.1 gives C D = 1.2
2 k +1
Pdn Pdn
k k
k
q sc = 1,248C D APup −
. (k − 1) g Tup Pup P
up
2 1.25+1
1 . 25
80 1.25 80 1.25
q sc = (1,248)(1.2)[ (1.5) 2 / 4](100) −
(1.25 − 1)(0.65)(530) 100 100
FRZAN ALI 26
Check NRe:
= 0.0108 cp by the Carr-Kobayashi-Burrows correlation.
d (0.0108)(1.5)
FRZAN ALI 27
(b)
k −1
z up Poutlet k 1.25−1
Tdn = Tup = (70 + 460)(1)(0.8) = 507
z outlet Pup
1.25
oR = 47 oF > 32 oF
Heating may not be needed.
But hydrate curve may need to be checked.
(c)
Poutlet = Pdn = 80 psia for subcritical flow.
FRZAN ALI 28
1.5 Multiphase Flow
When the produced oil reaches the wellhead choke, the wellhead
pressure is usually below the bubble-point pressure of the oil. This means
that free gas exists in the fluid stream flowing through choke. Choke
behaves differently depending on gas content and flow regime (sonic or
subsonic flow).
Tangren et al. (1949) performed the first investigation on gas-liquid two-
phase flowthrough restrictions. They presented an analysis of the
behavior of an expanding gas liquid system. They showed that when gas
bubbles are added to an incompressible fluid, above a critical flow
velocity, the medium becomes incapable of transmitting pressure change
upstream against the flow
FRZAN ALI 29
1.5.1 Critical (Sonic) Flow
Several empirical choke flow models have been developed in the past
half-century. They generally take the following form for sonic flow:
CR m q
p wh = n
(1.12)
S
Pwh = upstream (wellhead) pressure, psia
q = gross liquid rate, bbl/day
R = producing gas-liquid ratio, Scf/bbl
S = choke size, 1⁄64 in.
FRZAN ALI 30
1.5.2 Subcritical (Subsonic) Flow
k
k (1− x1 )V L (1− y c )
+
k −1
k −1
x1VG 1
yc = 2 (1.13)
n n (1− x1 )V L n (1− x1 )V L
k −1 + 2 + x1VG 2 + 2 x1VG 2
k
(1.14)
FRZAN ALI 31
1.5.2 Subcritical (Subsonic) Flow
(1.15)
The critical pressure ratio Yc can be solved from Eq. (1.13) numerically.
p2
ya = (1.16) ya = actual pressure ratio
p1 p1 = upstream pressure, psia
p2 = downstream pressure, psia.
FRZAN ALI 32
1.5.2 Subcritical (Subsonic) Flow
The total mass flux can be calculated using the following equation:
0.5
2 (1 − x1 )(1 − y ) x1k
G2 = CD 2 g c 144 p1 m 2 + (VG1 − yVG 2 ) (1.17)
L k −1
FRZAN ALI 33
1.5.2 Subcritical (Subsonic) Flow
Once the mass flux is determined from Eq. (5.17), mass flow rate can be
calculated using the following equation:
M 2 = G 2 A2 (1.19)
A2 = choke cross-sectional area, ft2
M2 = mass flow rate at down stream, lbm/sec.
M L 2 = (1 − x2 ) M 2 (1.20)
FRZAN ALI 34
1.5.2 Subcritical (Subsonic) Flow
M G 2 = x2 M 2 (1.21)
Gas volumetric flow rate at choke downstream can then be determined using
gas law based on downstream pressure and temperature.
FRZAN ALI 35
1.5.2 Subcritical (Subsonic) Flow
2) For predicting liquid rates of oil wells and gas rates of gas condensate
wells, a discharge coefficient of CD =1.08 should be used.
FRZAN ALI 36
An Example Calculation with Sachdeva’s Choke Model
Input Data:
Choke Diameter (d2): 24 1/64th in.
FRZAN ALI 37
Pre-Calculations:
Gas Specific Heat ratio
(k=Cp/Cv): 1.4
FRZAN ALI 38
Critical Pressure Ratio Computation:
k/(k-1) = 3.5
(1-x1)/x1 = 999
n/2 = 0.500043
VL/VG1 = 0.004811
VL/VG2 = 0.002287
FRZAN ALI 39
Flow Rate Calculations:
Pressure Ratio (yactual): 0.625
1.1 A well is producing 40o API oil at 200 stb/d and no gas. If the beam
size is 1”, estimate pressure drop across the choke.
1.2 A well is producing at 200 stb/d of liquid along with a 900 scf/stb of
gas. If the beam size is 1/2”, assuming sonic flow, calculate the
flowing wellhead pressure using Gilbert’s formula.
1.3 A 0.65 specific gravity gas flows from a 2-in pipe through a 1-in
orifice-type choke. The upstream pressure and temperature are 850
psia and 85 oF, respectively. The down stream pressure is 210 psia
(measured 2 ft from the orifice). The gas-specific heat ratio is 1.3. (a)
What is the expected daily flow rate? (b) Does heating need to be
applied to assure that the frost does not clog the orifice? (c) What is
the expected pressure at the orifice outlet?
FRZAN ALI 41
1.4 A 0.70 specific gravity natural gas flows from a 2-
in pipe through a 1.5-in nozzle-type choke. The
upstream pressure and temperature are 120 psia
and 75 oF, respectively. The down stream
pressure is 90 psia (measured 2 ft from the
nozzle). The gas specific heat ratio is 1.25.
FRZAN ALI 42
1.5 For the following given data, estimate
upstream gas pressure at choke:
Flowline ID: 2 in
Choke discharge
coefficient: 0.95
FRZAN ALI 43
1.6 For the following given data, estimate
downstream gas pressure at choke:
Flowline ID: 2 in
FRZAN ALI 44
1.7 For the following given data, assuming
subsonic flow, estimate liquid and gas
production rate:
Choke Diameter: 32 1/64th in.
Discharge Coefficient: 0.85
Downstream Temperature: 30 oF
FRZAN ALI 45