OCW Drilling Hydraulics Lecture
OCW Drilling Hydraulics Lecture
OCW Drilling Hydraulics Lecture
Contents
Dye
This one?
or
this one?
or
Transition flow
this one?
Turbulent flow
Laminar flow
Transition flow
Turbulent flow
vd
2000
laminar flow
vd
2000
4000 transition flow
vd
4000
turbulent flow
vd
Re
dimensionless
vd
where:
v
d
fluid density
fluid average velocity
pipe inside diameter
fluid absolute viscosity
vd
Re
Laminar flow
Turbulent flow
vmaks
v vavg 12 vmaks
v vavg
v vavg
Qin
element ( mg )1 @ Station 1
HL
v1
Turbine
HE
z1
element ( mg )2 @ Station 2
Pump
v2
HA
Control volume
Qout
z2
2
Datum
The energy possessed by a flowing fluid consists of internal energy and energies due to
pressure, velocity, and position
energy at energy energy energy energy at
section 1 added
lost
extracted section 2
This equation, for steady flow of incompressible fluids in which the change in internal
energy is negligible, simplifies to
p1 v12
p2 v22
z
1
A
L
E
2
2g
2g
10
p2 v22
2 g z1 H A H L H E 2 g z2
z1 H P
z2 H T H L major H L minor
2g
2g
p1
11
2
64
v
L
HL
Re d 2 g
12
b. Turbulent flow
From Darcy-Weisbach equation for laminar flow
2
64
v
L
HL
Re d 2 g
2
v
L
HL f
d 2g
64
a simple mathematical relationship.
Re
For turbulent flow, f has to be solved empirically experiment need to be done.
13
Friction Factor
a.
Laminar flow
Darcy-Weisbach equation
2
HL f L v
d 2g
b.
where f 64
Re
Turbulent flow
In the literature (from 1900s current date), there are many studies that have been
conceded by various researchers.
14
Laminar
Transition
f 64
Re
Turbulent
Complete Turbulent
vd
Re vd
Moody Chart
15
Absolute roughness,
ft
mm
0.003 - 0.03
0.001- 0.01
0.0006 - 0.003
0.00085
0.0005
0.0004
0.00015
0.000005
0.0 (smooth)
0.9 - 9.0
0.3 - 3.0
0.18 - 0.9
0.26
0.15
0.12
0.045
0.0015
0.0 (smooth)
16
Attention
1. In this subject, SKM1043, the f that we are using, is the American friction factor, fAmerican.
2. The value of fAmerican is different to the one that used by the British
f American 4 f British
64
Re
16
Re
f American 4 f British
17
Since the mud enters the drill string and leaves the annulus
at essentially the same elevation, the only pressure
required is to overcome the frictional losses in the system.
Hence, the discharge pressure at the pump is defined by:
pt p s p p p c pb p ac p ap
....... (3.1)
where:
pt =
ps =
pp =
pc =
pb =
pac =
pap =
18
19
N Re
where:
928 vd
N Re = Reynoldss number
v = average velocity of flow, ft/sec
.......... (3.2)
20
21
22
23
24
Hose
Swivel, etc.
Kelly
Eq. length,
3.826 ID
ID
Length
ID
Length
ID
Length
ID
Length
40 ft.
2.5
45 ft.
20 ft.
2.25
40 ft.
2,600 ft.
2.5 x 10-4
3.5
40 ft.
2.5
55 ft.
2.5
25 ft.
3.25
40 ft.
946 ft.
9.6 x 10-5
45 ft.
55 ft.
2.5
25 ft.
3.25
40 ft.
610 ft.
5.3 x 10-5
45 ft.
55 ft.
30 ft.
40 ft.
424 ft.
4.2 x 10-5
25
.......... (3.5)
where:
ps =
q =
m =
E =
p
4
144p Yt mv
3
where:
144p = pressure drop, lb/ft2
4
3
Yt = Yb, lb/ft2
m = L/(1500d2), slope of linear portion (from Eq. (3.3))
26
27
p vL
LYb
p
300d 1500d 2
L
p
300d
where:
pv
Yb
5
d
laminar flow
.......... (3.6)
28
Avg. velocity
inside the pipe
Avg. velocity
in the annulus
1 ft 3
1 min
q gal/min
7.48 gal
60 sec
( /4)(d /12)2
q
....... (3.7a)
v
2
2.45 d
OD
ID
Hole
q
....... (3.7b)
v
2
2
2.45 (d h - d p )
where :
v average velocity, ft/sec.
q = flow rate, gpm
d = diameter, in.
Drill pipe
Annulus
Annulus Area Ah Ap
(d h2 d p2(OD ) )
29
p
v
Substituting the above Eq. for in the Reynoldss number of Eq.
(3.2), equating the resulting equation to 2000, and solving for v
gives:
1.08 p 1.08 p2 9.3 d 2Yb
.......... (3.8)
vc
d
where:
vc = critical velocity, ft/sec, above which turbulent flow
exists and below which the flow is laminar.
Eq. (3.8) assumes that turbulence occurs at NRe = 2000. Therefore,
if:
v vc , flow is laminar
v vc , flow is turbulent
30
NRe
NRe
928 vd
t
2,970vd
.......... (3.10)
31
32
Example 3.1
Mud is flowing through 4 1/2 inch OD, internal flush drill pipe.
Calculate the frictional pressure drop per 1000 ft of pipe.
Mud properties
Mud density, m
Pipe ID
Bingham yield, Yb
Circulating rate, q
Plastic viscosity, p
=
=
=
=
=
10 lb/gal
3.640 in.
10 lb/100 ft2
400 gal/min
30 cp
33
Solution 3.1
Eq. (3.7a) : v
q
2.45d 2
Eq. (3.8) : v
c
400
12.3 ft/sec
2.45(3.64)2
(1)
(2)
(2,970)(10)(12.3)(3.64)
44,300
30
(0.0062)(10)(1000)(12.3)2
(c) p p
100 psi/1000 ft
(25.8)(3.64)
34
hydraulic radius, rh =
(r12 r22 ) r1 r2
for an annulus
rh =
2 (r1 r2 )
2
r2 r
r2
r
re = r1 r2
or
de = d1 d2
.......... (3.11)
r1
35
v12 p2 v22
2g 2g
p1
where:
v1 , v2
= density, lb/ft3
= velocities at points 1 and 2, ft/sec
.......... ( a)
36
or
p1 v12 p2 v22
w 2g w 2g
.......... (a )
v22 v12
2g
Practically, v22 v12 v22 , therefore:
p
.......... (b)
v22 2 g
37
By substituting Eq. (c) into Eq. (b), and rearranging it, the equation
becomes:
q2
.......... (3.12)
p
2 2
2 gC A2
Altering Eq. (3.12) to practical units for mud flow, we:
q2
pb
7, 430 C 2 d e 4
.......... (3.13)
38
Multiple Nozzles
39
q
For parallel flow, q1 , where n = number of nozzles.
n
therefore:
pm q12 A2
q12 A2
2 2 2 2 2
p q A1 n q1 A1
A2
1
2 2
n A1
A2 n2 A12
or
A nA1
.......... (3.14)
40
.......... (3.15a)
41
Example 3.2
A 10 lb/gal mud is being circulated at the rate of 500 gal/min.
through a tri-cone bit having three 3/8 in. diameter jets. What
is the pressure drop across the bit?
Solution 3.2
Drill string
p1
Hole
Nozzle
vn1
vn2 vn3
p2
42
43
Solution 3.3
5,500 ft.
2.45 7 8 4 2 (3)
307 gpm
500 ft.
Table 3.1
E 9.6 10-5
q
307
8.56 ft/sec
2
2
2.45d
2.45(3.826)
m d
(10)(3.826)
4.25 ft/sec
45
NRe
2,970 vd
N Re 32, 400
Curve II
(2,970)(10)(8.58)(3.826)
f Lv 2 (0.0066)(10)(5,500)(8.56)2
pp
269 psi
25.8 d
(25.8)(3.826)
46
15.84 ft/sec
2
2
2.45d
2.45(2.813)
The critical velocity:
vc
m d
(10)(2.813)
4.64 ft/sec
47
NRe
2,970 vd
N Re 44,100
Curve II
(2,970)(10)(15.84)(2.813)
44,112 44,100
30
Fig. 7.1 f 0.0062
f Lv 2 (0.0062)(10)(500)(15.84)2
pc
107 psi
25.8 d
(25.8)(2.813)
48
1
3
q
307 / 3
0.41 in.
2.45v
(2.45)(250)
nozzle diameter of
13
32
in. is chosen
102
v
252 ft/sec
13 2
2.45( 32 )
1
3
1
3
1
3
49
Eq. (3.15a) d e nd 2
Eq. (3.15b) de ad12 bd 22 etc.
Using Eq. (3.15) or (3.15a), the actual nozzle diameter:
2
d 3 ( 13
)
0.704 in.
32
q 2 m
Eq. (3.13) pb
7, 430C 2 d 4
Pressure drop across the bit, pb :
(307)2 (10)
pb
573 psi
2
4
7, 430(0.95) (0.704)
50
(10)(1 )
7.26 ft/sec
51
NRe
2,970 vd
(2,970)(7.62)(1 18 )
8, 487 8,500
30
N Re 8, 400
Curve IV (for annuli
in uncased hole)
f Lv 2 (0.0098)(10)(500)(7.62)2
pac
98 psi
1
25.8 d
(25.8)(1 8 )
52
v 3 ft/sec
d a d1 d 2
d 7 78 4 12 3 83 in.
The critical velocity:
vc
v vc
3
8
(10)(3 )
4.39 ft/sec
5,500
pap
300 (3 83 )
30 (3)
83 psi
10
3
5(3 8 )
pv
Yb
5
d
53
q p
HP
1,714 v m
.......... (3.17)
where :
q = flow rate, gpm
v = volumetric efficiency
m = mechanical efficiency
Assuming volumetric and mechanical efficiencies of the pump are 90% and 85%
respectively:
HP
307 (1,166)
273 horsepower
1,714(0.90)(0.85)
54
Summary
Bingham Plastic Model: Calculation Steps
pt p s p p p c pb p ac p ap
p
vc
Eq. (3.8)
E q . (3 .1 6 ) o r F ig . 3 .3
v
Eq. (3.7a ) or (3.7b)
de
Eq. (3.15a) or (3.15b)
pb
No
(lam in ar)
p p , pc , pac , pap
Eq. (3.6)
if
v vc
(Eqn. 3.13)
Yes
(tu rb u len t)
N Re
Eq. (3.2)
f
(Fig. 7.1)
p p , pc , pac , pap
Eq. (3.4)
55
p
p pu
9.5 3.2(3)
Example 3.4
Solution 3.4
From Example 3.3: q = 307 gpm, bit = 3 13/32 in. nozzles
(a) Surface equipment losses ( ps)
q = 307 gpm
Curve type 2
Fig. 7.3
pu 27 psi
0.14
10 30
ps 27
33 psi
9.5 3.2(3)
(b) Pressure losses inside drill pipe ( pp)
q = 307 gpm
Fig. 7.5 (for 4.5 d/p)
Curve 7
(assume ID = 3 )
10 30
p p 176
9.5 3.2(3)
32
pu
5,500 176 psi
1,000
0.14
217 psi
0.14
56
p
p pu
9.5 3.2(3)
Fig. 7.7
(assume ID = 2 )
10 30
pc 75
9.5 3.2(3)
13 "
32
93 psi
Fig. 7.9
nozzle
pb 550
15
500 75 psi
100
0.14
pu
10
579 psi
9.5
p pu
pu 550 psi
m
9.5
0.14
57
p
p pu
9.5 3.2(3)
Fig. 7.10
(bit size = 7 7/8)
10 30
pac 125
9.5 3.2(3)
25
500 125 psi
100
pu
0.14
154 psi
Fig. 7.10
(bit size = 7 7/8)
10 30
pap 77
9.5 3.2(3)
pu
1.4
5,500 77 psi
100
0.14
95 psi
0.14
58
Surface connections, ps
36
33
269
217
107
107
Bit nozzles, pb
573
579
98
154
83
95
1,166
1,185
Total circulating
pressure, pt
59
Additional Information
Besides Newtonian and Bingham Plastic Models, there are
several other model used to predict pressure losses in mud
circulating systems.
Generally, each model is based on a set of assumptions which
cannot be completely fulfilled in any drilling situation.
Power law, Herschel-Bulkley (Yield Power Law @ API Power
Law) models are the most widely used in the oil industry.
Table 3.3 shows a summary of pressure loss equations