Apuntes de Transporte de Hidrocarburos
Apuntes de Transporte de Hidrocarburos
Apuntes de Transporte de Hidrocarburos
...... b.o (/
._~
./
/
APUNTES DE
TRANSPORTE DE
HIDROCARBUROS
FRANCISCO GARAICOCHEA P.
Noviembre de 1983.
TRANSPORl'E DE HIDROCARBUROS
CONTENIDO
,1.1.- Introduccin.
r'Vg
R
P = 18" ( _'s_,_ )
.'
0.83
, (1.2.1)
Por lo que ~espejando la reacin gas disuelto aceite {Rs ' de 1,11. aCU!-
cian anterior' se tiene:
100.0125 ( Yo)'
R (1.2.2)
s 100 .00091 (T)
Donde:,
F (1.2.4)
Donde:
R
s
e x p (C 3 ( >;,/tr + 460) (1.2.6)
1-2
Los valores de loa coeficientes son:
Coeficientes . 30 API
Yo
Cl 0.0362 0.0178
C2 1.0937 1.1870
C3 25.724 23.931
(1.2.7)
(lag, Pl3
Rs Yg
3.- Calcule B *
o
t' 0"1
P
T 0.130
con:
Yo
0.;816
0.526
B *
o
R
s 'i Yro ) + 0.968 T (1.2.10)
I-3
4. - Determine B con:
o
log (BO - 1) *
- 6.58511 + 2.91329 log Bo * - 0.27683(10g Bo )
2
(1.2.11)
b 5.44 (R
s + 150)-0.338 (1.2.15)
10 x - 1 (1.2.16)
.l1om
(1.2.181
Donde:
1-4
;L. 3. - Propidades de;L aceite bajosaturado.
Donde: al .= '-
1433, a 2 = 5, a 3 = 17.2, a 4 = 1180, - aS 12,61, a 6 = 105
1.3"2.- Densidad
., del aceite bajosaturado
(1.3.2)
Donde:
PObdensidlid del aceite a Pb
~. . 2
1.3.3.- Correlacin para calcular la viscosidad del aceite bajosaturado.
c
m = el p2 e~ (e 3 + e 4 p) (1. 3.4)
Donde:
el = 2..:~
-5
Ci = 1.167 e4 = - '8.98 x 10
}lob - viscosid<i.d
I-S
-_._-~-~--~~---'
.. ,....
Es obvio que el primer paso para obtener ~.as propiedades del aceite bajo-
saturado es la determinaci5n de la presin de saturacin del aceite.
Esta presi5n puede obtenerse con la Ec. (1.2.1) cuendo el aceite es de
bajo encogimiento. Para aceites vOltiles se puede ''',<ir
la correlacin
de 0i stiIi :
* R 0.816 T 0.130
Donde: Pb =(-v:) (1.3.7)
gd 1'0.989
o
El valor de PJ:> .as obtenido puede correqirse usando una ecuaein similr
.a la (1.19). Es conveniente indicar adems que generalmente la presin
de saturacin del<gua se considera igual a la del aceite.
dad relativa del gas libre ( P9f) Y la densidad relativa del gas disuel-
1-6
PMp.i.e.da.de de. ta.6 Fttddo.6
2.7044 P
460) (1.4.5)
I-7
(1.4.11)
(1.4.12)
2
Z = til + Aa / Tpr + A3/ Tp ; 1 Pr + (A4 + AS/Tprl Pr
, 5 ' 2 3 ' 2
+ (As A6 Pr 1 / T pr +CA7' Pr /Tpr ) (1 + A8 Pr )
Donde:
Al "0.31506 AS =- 0.6123
- 1.0467 A6 =- 0.10489
A3 .. - 0.5783 A7 = 0.68157
(1.4.14)
P
pcm = {1 - yN2- yC0 2 - y H2S1 P
pc
+ 493 yN 2
+ 1071 YC02 + 1306 yH 2S (1.4.151
1-a
Estos valores debende ajustarse, pr su contenido de CO 2 ')' H2 S, median-
te las ecuaciones de wichert y ~iz:
(l;~.17)
T' = 'r - C
pcm pcm 'lila
Bw'" 1.0 + 1.2 : 10-4 (T-60) + 1.9 x 10~, (,+,,-pO! - 3.33 x 10-6 P
(1.5.1)
I-9
1.5.3.- Correlaci~,paradetezminar .la viscoslda.d del a~ saturada. 9
. '. .. -li i
f(p, T) = 1+ 3.5xlO (P) (T-40) (1.5.6)
60F ~ T ~400F
(l.s.S)
=,[. .
,
(28Q -'1\ )
206 !7wl ) + (1.5.9),
-lO
1.S.S.-Correlacin para determinar la solubilidad del gas en el agua. 8
donde:
.6~21'
(1.6.5)
Donde:
1-11
El ranqo de aplicaciGn de e!lta correlaci6n es:
(io6.6)
(1.6.7)
Ejemplo 1.1.
Soluc;ioo:
. p -18
_C.i
el de Rs' Aplicando la Ec. (1.2;1)
( Pq ) '
0.83 10.0009l(T>l
100.0125J Vol j
[
358 0.83 100~00091(200) J
p= 18 . ( 0.75) 100.0125(3O>J
1-12
Ejemplo 1.2.-
100,0125 (38)
R
s
o
.L
8[4500
18 100.00091 (240)
R = 12113 . 3 gd
pJ.es ~ ca C.s.
s 'bl a.c.s.
o
F, = Rs ,( Jl / Yo ) 1/2
C( + 1.25T
B 1.1.815
o
I-13
Se observa que el factor de volumen del aceite obtenido, al aplicar la
correlacin de Standing, resulta muy bajo y no corresponde al de un -
aceite supuestamente voltil. Esto se debe a que la correlacin de Stan-
ding fu desarrollada utilizando muestras de aceite de bajo encogimien-
to. Por lo tanto, estos resultados (los de R y B ) no son confiables.
Ser necesario utilizar otras coorrelacionessparaocalcular los valores
de Rs y Bo'
o,a028
R
s
. 3
R 1136.54 pl.es gd a c.S.
s
bl a c.s.
o
B 1.6245
o
De la Ec. (1.2.8)
r-14
p. = .27.901
r' = 141.5
ro ],31. 5 +' API
141.5 0.835
131.5 + 38
0.526
B
+ 0.968 T
o
0.526
BO = 1622,.5412 (. g::35 + 0.968 (240)
Be * 1819
iog (B o - 1) =- 0.0300085
BO '" 1.9332
I-15
A continuacin se muestra una tabla comparativ,a de los valores obtenidos
con cada una de las ,correlaciones utilizadas y el reportado de labor ato-
rio.
R B
o
s
R (1.19)
sr
Donde:
R valor de R corregido
sr s
R valor de Rs obtenido con la correlaci6n
sc
Rslab valor de Rs obtenido ,del laboratorio.
Ejemplo 1.3.-
Calcular la densidad y la viscosidad del aceite del ejemplo anterior ~
utilizando los valores ,Obtenidos en el laboratorio.
36.54: lb/pie)
I-16
'.';'
10 2 2636, = 183.4847
x = y T -1.163
Pam
103129 1 = LO~565 cp
Pam
b = 5.44(R + 150)-0338
s
.. = 0.2241 U.05565)0.4259
"'0 '
r.:. ,;
Po' ='0.2293 cp
EjelllPlo 1. 4.-
C
o
= (-14JJ + (5) (1725) + (17.2) (240)' + (-1180) (0.8028; + (12.61) (38) )
/ 105 (4500)
C
o
= 2.4115 X 10-5 pg2/ lb a p = 4500 lb/pg2 Y T = 240F
C
o
= (- l4J3 + (51 (17251 + (17.2) (240) + (-1180) (0.8028) + (12.61) (38))
/ 105 (6000)
C
o
= 1.8086 x 10-5 pg2 /lb a p = 6000 lb/pg2 Y T = 240F
P a J7.888 lb/pie 3
O
1-18
que es la densidad del aceite a p '" 6000 lb/pl Y T 240F
1 1~7 . '-5 .
ro ;= (2.6) (6000)' exp (-;1.1.513 + (,..8,98 x 10.(6000))
m = 0.4630
!lo
Po = (0.22933) (6000/4500)0.4630
1-19
REFERENCIAS CAP. 1
5.- Baker, O.: "eesigning Pipelines for Simu1taneous F10w of Oil and
Gas" Pipeline Engineer, Feb. 1960.
10.- Wichert, E., y Aziz, K.: Ca1cu1ate Zs for Sour Gases". Hydro-
carbon Processing, Mayo 1972.
11.- Benedict, M., etal, "An Empirical Equation for Thermodinamic
l'roperties of Light Hydrocarbons and Their nixtures". J . Chem.
Fhys. Vol. 8. 1940.
1-20
------------- -------------------------'--------
2.-FUNDAMENTOS DE FLUJO A TRAVES DE TUBERIAS.
CAt!BIADOR DE CAIDR
llb
m
2-1
Fundamen.tO-6 de F.t..u.jo tt Tltav6de TubflJa.6
. . . . . . . . . . . . . 11 (2.1)
Donde:
Donde:
3
V.- volmen especfico t~)
lb
m
Energa potencial(Ep
(2.3)
2-2
Fundamen-t:o.6 de Flujo a. tlta.vA de TubeJQ.6
Donde:
V. - velocidad (pie l
seg
.lh + + dW =o (2.6)
s
Donde:
+ p .+ P dWf
+ --zrI;" =O (2.7)
+ (2.8)
Donde:
2-3
Fundame.n.tM de Flu.jo a TlUlV de Tu.bviU>
(~ (2.10)
.1L I f
(~) (2.11)
.1L f
Donde:
7r d d/4
(2.12)
2-4
2.2.3 Factor de friccin.
El valor del factor d/;!, fri~cin (f l' e,s funcin de la' rugosidad .de, la tu-
beda ( E ) Y del nUmero de~ynolds (NRe) , esto 'es:
N (2.13)
Re
Para flujo laminar de una sola' fase" el factor de fricc.6n dep~n~e exdu
sivamente del nmero de Reyrlolds y est 'dado por:'
64 ,
(2.14)
lii
Re
Para flujo turbulento (NRe :=::' 3100), el. factor de friccin est dado pOr
la ecuacin de Col.ebrook y White:
2-!?
Fundamento~ de Flujo a Tnavl6 de Tub~
-2 ,
21.25 )] "
f = [1.14 - 2 109<+ + (2.18)
~
Re
2-6
o._tO
o.o~ ~~__.___ ~_L____ :l_~_ -
Flujo Zona lona d.
O.oe lominar I crtica transicin Turbulencia completa, tubos- rUl)osos -
0.07 \ I //'/1 - - -
~ ~ r--- .2:..'\.
\ 0.04
0.06
U
e
'O
:E
..,
Q)
~
e
Q)
:~
0.04
0.03
I
I(
W ~' ~RCT
Flujo lamlnar\
I:~\
/;/f;,
A%~ ~ ~ ~
1---
~
'" K
0.002
+~
]
o
o
"ti
.;
g,
r--
N
Qj
o
U 0.02
Re \ ~==-- " "- 0.001
::l
n::
~
~T"b'"'~
\ 0;0008
.0906
Ex10 5 (pg)
r-::
~"-
Material 0.0004
VIdrio. cobreplstico, hule 5.905 .
Fierro fundido n8VO '
Fierro fundido semioxldado
19.7 a 39.4
39.4 059.0
"- o,po02
Fierro fundido oxidado 59.0 a 118 "-...
Cementa liso 11.8 o 31.5 0.0001
0.00008
Acero 157 o 394 0.00006
-.:::
~~04~. 7 870 mas
Albllto comento
0.01 r-- Concreto
0.00004'
+ (3.1)
(3.2)
como:
p = 62.428 Yr (p~s3 t (3:4)
y
4 q' ( pies)
v' (3.5)
seg
" d,2
v,2 42 q'2
(3.6)
~
" d
Sustituyendo (3.4) Y (3.6) en (3.3), se tiene:
1 pies 5 ) (3.8)
5
12 5 pq
3.1
Flujo de L(quido~ po~ Tub~
. 3 . J
q' ( E!.L
seg
i z q bl
(dIa)
5.6146 ( p;.,es ) _1_- (d.a 1
bl 96400 seg
q2 {S.614612
(J.9)
86400:1
dp'
f d (. lb )
Pf - 2 - ' (3.11)
pg
(3.12)
0.06056 f)~ q2 L
(3.13)
dS
(3.17)
f d',.. ,
3.2
Sea:
. 3 bl . 3
q' (Pl.::cj ) .. q (d!a ) 5.6142 tPl.:~, ) "=,,,l~__ , (di~ ) (3.18)
, 36400seg
_1_ _
d' (pies) = d(pg) (~) (3.21)
12 pg ,
N ' 92.2' ~
Re (3.22)
dft
3.3
Fi..uJo de. LZqtUdolJ pOJr. rub~
Soluc,in:
64 ~
f 0.05456
NRe 1173
0.06056 f ,\:"r q2 L
d5
P f = 11.7 lb/pg2
3.4
Ftujo de Uqu..i.dO.6 pOlt TubeJa.6
0'
f
c'
,+ 2~14
.J f N
)J-2
'
(3.24)
s Re
.f
c [-2 log (0.000041 + 0.0004286
"fs
Iteracin f
c
1 0.03700
2 0.03576
3 0.03596
4 0.03592
5 0.03593
Por lo que:
0.06056 f Yr l L
192.8
La: ecuacin (3.23) puede aplicarse tambin para obtener' el dimetro para
un gasto y cada de presin dados. El procedimiento de solucin es itera-
tivo, ya'que el nmerq de Reynolds (y por lo, tanto el factor defriqcin)
es funcin del dimetro.
3.5
Ftujo de Uqu.i.do.6 poI!. TubeJ1:t6
d
0.06056 f Y r q2 L
[ --=-=-==2~=---
E (dPT - 0.433 Y r dh)
.J 2
(3.25)
q
E[ .=d5-!..(_d..;:;P...T__-...;0:..:.:...:4:.::3.::.3_Y;.,lr;...:dh~),-- JO. 5 (3.26)
0.06056 f Yr L
a) La estacin de bombeo
b) La l!nea superficial de distribucin
c) La tuber!a vertical de inyeccin
d) El yacimiento
3.6
Estacin
de bombeo
;'
Tuberlo vertical
de inyeccin
3.7
(J~27)
Para simplificar los clcu,l.os es conveniente graficar los ,valOres dia '
los gastos supuestos contra los de, la presi6ri esttica calculada. Co-
mo se muestra' en :la Pig,. 3,2.
3.8
P REAL
ws
3.9
Para la lnea superficial:
r rw =1.00
.
L = 2 millas.. p.=.1 cp
f= 0.0006 pcj
Solucin:
GASTO SUPUESTO
(pl/dia) Pth Pwf P;"s
(lb/pg2)
3.:1.0
Le
Flujo de. Llqui.d04 poI!. TubeJa.6
L.D. q T.P. q Pp 2 q
(pg) (bl/dal (pg) (bl/da) (lb/pg 1 (bl/cta
Los codos, las vlvUiasy las conexiones, as como otros tipos de. adita-
mentos, incrementan las prdidas de presin por friccin en el sistema.
En el mejor de los casos, stas prdidas adicionales de presin slo pue-
~en ser estimadas en forma aproximada. Ellas deben de incluirse en el an!
iisis de un sistema, considerando que cada conexin es sustituida por una
longitud equivalente de tubera recta (Le), la cual producir la misma
prdida de presin por friccin que la conexin real. La longitud equiva-
lente de cada conexin se agrega a la longitud axial de la tubera (L),
antes de calcular la prdida de presin total en el sistema. Los valores
aproximados de Le, para las vlvulas comunes y otras conexiones, para flu-
jo turbulento, pueden obtenerse mediante el nomograma mostrado en la Fig.
(3.7). Por jemplo un codo de 2 pg., en ngulo recto,es equivalente a 12
pies de ima tubera recta.
3.11
q (blIda)
8000
7000
6000
5000
3 6 L.O. (pg)
Fig :3.4 Relacin entre el gasto de inyeccin y el dimetro de la
lnea de descarga
q (bl/da)
13000
11000
9000
7000
,
11000
9000
7000
5000L--~---~----~-~----...I.--.....
1000 2000 3000 4000 5000 6000
Pp ( Ib/pg 2)
Fig :3.6 Relacin entre el gasto de inyeccin y la presin de
bombeo
3.12
-
-----------------------------------------------------,
8
.(
...
!,!
o:
lo.
'"'"
~!
:>0:
>'"
'~~ 3000
ZOOO
1000
50
VLVULA ANGULAR,
AmE1tTA
CONEXiN REGULAR
EN T
SOO 42
48
36
@ / U? CODO EN NGULO
RECTO
300
200
30
22
24
30
16 o
14
ff] 50
10
IZ
~~~~l~d~ t~ ciR~- 30
9
lO
FlelO LATERAL
DE DESCARGA 20 8
'"'"a:: -_o .
.~
~tB.t
Z -6--
10 '"
~ ".-'
"Z 5 S
...
Uf
Q
~
CODO ORDINARIO, U
... 4t
O CONEXiN EN T
REDUCIDA A LA '"
o:
".::>-'
1~
MITAD S ,~
4
rJ
~::O ~EGULAR
a: >- Jt z
3 '"
ID
Z
'"
;;: J 3 '"
~
.
~tB.t
a: ~
r~
o
2 ..J U
2~ ...1-o:
~~oo~O~~~~~N~
~TRACCIN SUelTA
'"o
'"
1-
.
,~
en
2 2 "
oa:
T REDUCIDA A 1/4
...z :>
...
1-
d/D-I/4
d/D-I/2 ..J
.
1- I~
,~
~
.
J
.-;" d/O-3.14
:;
..J
1\4 Q
0.5 a o
\$
= gEc'b~;la~AR- ~
1-
0.3 ;
z
0.2 o
'":>o .
..J
z
~
o
z ~
REGULAR EN T cOOO DE 45
..J oa:
...
1-
t
0.1
~
Q
Fig 3.7 0.5
-Cuadro de alineaciul utilizable en la determinaciu de la re;;;islel1cia al fl~o~ ofrecida por las vlvulas y
conexiones. en trmDos de la lonqitud equivaleute de una tubera reda. (Cortesa de la Creme Co.)
3.13
4. FLUJO DE GAS POR TUBERIAS Y ESTRANGULADOPES.
De la Ec. (1.36):
p'
2.7044 p ro (4.2)
Z (T + 460)
" i
Bg. segn se vi (Ec. (1. 34 es:
p (lb/pg21 , T ( F)
4.1
<-4.81
,', 3 3
qi2 (....e!!!L)2 .q2 ( Pdiieas,' "1 2 '~2 ( d!a2 1 . , (4,.~)
seg 86400 seg
P l + P2
Ademis: p .. 2
, y p = Pl - P2
Donde:
Pl = presin a la entrada del dueto
f "g q 2 Z L . T + 460)
(4.11)
lb 1, .i 2
( --'- ) """""i44 (~) h (pies) (4.12)
pie 3 pg
(4.13)
4.2
_2
0.03756 P I'a h
Z (T + 460)
(4.16)
h (pies), T( F)
q = 2744 (..--;.)
1 05 2 5
d'
['. . .p, ~"'" P~
- - (0.0375
...:::.".....--'~_ _ _ _ _-'2.._-,-_--'_ _-I
V h p,2/ z (T + 460) .0.5 , ( 4.1n
,f V.. L Z(T + 460)
Asimismo se tiene:
f Y t Z T +460.) (4.1B)
d = .:[. ':l2 9
:2 2 .
( Pl- P . - -2 - ' "
274i .. 2 (0.0375 V.. h P Z ,(T + 460,
n d,2 . P
(4.21)
nd
NRe .. 0.09727
q' Vg '
(4.22)
d' ",'
4 ..::
procediendQ ahora al cambio deuni~des:
3
pies 3 ......J.-(~)
q'(~) q ( dia) (4.23)
seg 8640.0. seg
lb
m
fl (ep) 0..0.0.0.67197 ( p1e-seq ep) (4.25) .
p. '( pie - seq
4.3.,.. OBSERVACIONES:
4.4
F.f.LLjo de Ga.6 polt TubeJLu Y Etltangula.dolte&
NODO
SEPAI'.AOOR
ESTRANGULADOR
INTERCONEXION
FONDO DEL POZO
RADIO DE l)P.ENE
4.5
Flujo de. Ga6 polL TubeJa.6 !J E&:tJr.angut.adolLeA
al Suponga un gasto
Pwf =J P;s
_ ( 51. ) lln
C
(4.30)
4.6
Flujo de GIU pOlt TubVlU> IJ E.tltanguladolte.6
1 , ,,
/'
/'
/'
\.
,\
\
\
\
\
\
\
\
'\nax
4.7
Flujo de Ga6 '{XJ1t TubelLl.6 1jEt.t:Jumguladoltu
2 rl-p-w-S-=-5-0.... ~....:...lb-/~
Pws - 1000 lb/pg
q Pth
l--lq-X-1-0-6-'1 :: I
2 x 106 834
q Pe
2 x 10 6 153
4 x 10 6 292
10 x 10 6 698
12 x 10 6 828
16 x 10 6 1081
4.8
p{lb/pgl)
Pth
1200
\
800
\
\ \
\ \
\ \. \
Pws = 2000 (lb/pg2)
\ Pws = 1500
\p.ws =1000
\Pws =500
4 8 12 16 20
Qg xl(f6 (pies!dlo)
2000
1500
1000
500
%~------~--------~------~15-------'
qgx 10--(pies'/da)
4.9
500r---------------======~===
Pwf
400
," "
300
200 \
\
100
~--------------- Ps =50(lb/pg2 )
0~-----L------L------L------4L-----~
qg X lO-6(pies/da)
0.90..----,--........,.........,............--,----,..........--...,
Yg
0.80 I----+----'-\-l-\--\-~ri"'--'-:..:::..::..-~--_I
0.70 I----+---=.::-=--:....--+-~*~--I-----I
0.601----+-----l----4<-~~~l----,---1
(4.31)
(4.32)
Donde:
e calor especifico a presin constante
K~
ev -calor especfico a volumen constante
[~ P2 2/K P2 . (K+l)/KJ]0.S
<SSSO,,".
64.34 K
I'g(T + 460) (K-l)
--)
Pl f
(-)
Pl
Donde:
2
A . are a del estrangulador, pg
C coeficiente de descarga del estrangulador.
------- --------------------~----------------
T = temperatura en el estrangulador F
= presin corriente a.;-r~a .del estrangulador (Pth)
Donde:
0.6 t<:=1.31
l'g
- 2
p= 700 lh/pl P2 = 600 lh/W
1
Solucil. 1
2 ,K/(K-l) .. 0.545
K+ 1
P2 / Pl .. 0.857
4.12
Ftujo de Ga poJt TubeJM tJ E~.tJtangu..f.a.doJtu
CA = __________-.~------~8~0~0~0~0~~----------------------
(155500) 700[ (0.6)
(64.34) (1.31)
(610) (0.31)
[(O . 857)1.527- (0.857)1.763]]5
CA = 0.005068 pg2
4,13
REFERENCIAS. CAP. 4-
4.14
5. - FLUJO MULTIFASleO EN 'rUBElUAS.
(5.1)
:; 1
5.1
Estratificado
Flujo Ondulado
segregado
Anular
Tapn
Flojo
intermitente
Bache
'>!,'>jujo
distribuido
N
LV
< 10 (0.321 - 0.017 N<Jv - 4.267 sen e - 2.972 N L
5.2
nmero de la viscosidad del lquido O 15726 .. '1/ PLtl3) 0.25
." r-L '
N
gv nmero de la velocidad del gas
Direcci6n Tipo de
del Flujo Flujo Cl C2 C3 C4 C5 > C6
(5.2)
VSL (pies/seg) = ~~ =
0.01191 ('& Bo + qw Bw> (5.3)
d2
q' 0.002122 q:o (R ~ Rs) Bg
vsg (pies/seg) '" -?-- ,=
p
d2 (5.4)
5.3
-------------------------------------------------------------------------------------,
Donde:
v
,1 = sL
(5.6)
v
m
q' q' v
L L ~
vL (5.7)
~ ApH L HL
q' q' v
v ~ 9: =~ (5.8)
9 A A U-H L ) (l-HLI
9 p
Pm eL HL +
P9 (1 - HL) (5.9)
e
9
(1 - ,1 ) (5.10)
Donde:
5.4
v = volumen de la mezcla a c. esc. por barril de aceite producido a c.s.
m Pies~ a c. esc./blo a c.s.
M=M +M +M (5.12)
o 9 w
. 3
(lbo/pie~ (lb", l pl.eo
M
o Pro x 62.428 . 3 x 5.615 (-;;;- l
pl.e",
M 350.5 F (5.13)
o ro
3
(lb Ipie 3 lb pies a e.s
M
"g
!Z SI x 0.0764 e_a_ l x R ( 212
9
(lb Ipie 3
, 3
pl.e bl a c.s
a a o
a
M 0.0764 F R {5.14l
9 9
3
lb lb pie
I piew x 62.428 ( ~) x 5.615
M
w "w
lb
'" I pie 3 l
' 3
pl.e
bf-:l
w
w
w
'"
bl
x WOR ' w
'~
o
M = 350.5
"ro + 0.0764 R "9 + 350.5 I'w WOR (5.16)
pieso3 + 9 + W a c. ese. )
Clilcul0 de V
m
bl0 a c.s.
V
'mo 5.615 B
o
pies
3
!Zl a c. ese. a e.s )
pies 3 c. ese.
V (R - Rsl ( B
m9 bl0 a c.s. 9 3
pie a c.s.
5.5
. 3
bl PJ.e"w pies 3
( w a c.esc.)
~ l 5.615 ( b/w . 3
bl0 pJ.ew a c.s.
(5.17)
w (5.20)
m 86400
w = W + w + w (5.21)
m o w 9
w = Po qo B I 15391 (5.22)
o o
W = Pw B I 15391 (5.23)
w '\. w
+ (1 - A ) (5.25)
"ns= "L 4
5.6
.
,-,., =
(S.261
.
Donde:
P'
,L = (5~27) "'
BO
f
~
= B
o +WORB' (S.28)
w
(5.29)
~
es 'la relacin a~ite-lquido y fw es la' relacin ag1UL-lquido.
"
Donde fe)
+
(S.30)
5.7
Flujo MuUi.646i.c.o en Tube.Wt6
REFERENCIAS. CAP. 5.
10- Beggs, Ho Do y Brill, JoPo: "A Study of Two Phase Flow in Inclined
Pipes". Transo AIME . 1973.
508
6. FLUJO MULTlFASlCO ,EN TUBERIAS HORIZONTALES.
6.1INTRODUCClON.
(
AL
.1L )T ("""*' ) f + ( 1t) ac (6.1)
O sea:
( -1f-) f e v
2
+
e (v2
(6.2)
T 2 gc d 2 gc L
6.1
l
Flujo Mu.Ui.64.6..i.c.o en TubeJLu Holzonta.l.e.6
<lP
6.2
6. 3~ - CORRii:LaCION DE B~TIIZZI, TE!{ y. POET'l1WiN.( 1)
6.3.1.. Introducci6n.
(6.4)
(6.5)
En donde:
y ,,= (6.9.)
.l?.P
174 . 158ft .
m
,i
(6.10)
A L e ' n s d5
6.3
0.100
S---' KA
0.080
0.060
...... ......
" Curvas Ay A'
'"
Relacin de masa oos-acelte
0.6-1.0
o.
K Curvas S y S' 0.4-0.6
..:: 0.040
' ........... '"'" ........ Curvas C yC' 0.2-0.4
~
-15
"'" --
Curvas O yO' O -0.2
'C) r--
-
.~ 0.Q20 """ 1'-...
........
f'-
""~
1.....
G) r-- t'-..
"t:I
~I'... r-.. :: 1- ...... .... -
'"
l5 0.010
.....
- - - -
"'
o 0.008
.
I.L 0.006 -- - --- - - B'
le' .-
O'F
0.004~
0.002 .'
..
0.00
102 103 104 105 106
4> = Funcin del nmer() de Reynolds. (NRe9)a(NR.~)b
Fig 6.1 Grfica del factor de disipacin de energa ( tomada de Sertuzzi y otros)(1)
ElnGlilerode Reynolds del. l!quido y ~. gas lile obtienen con las expre- .
. siones siguientes. .
~
(6.U)
w
N =.' 2.2737 ---L- (6.12
Reg d I'g
;:l. tactor de friccin puede obtenerse con la 1'1g. 6.1 o empleando las
ecuaciones si9w.entes: ( 2)
Para:
Para:
11 .:: 10 000
Para:
500 .~. 11 * 10. 000
En" donde:
<.6~l6)
6.5
6.3.3. Procedimiento de clculo.
4.- Cbtenga wm con las Ecs. (5.16) Y (5.20) o con las Ecs.(5.2l),
(5.22) Y (5.23). Obtenga .A con la Ec. (5.2)
d 6 pg. T 80F
).10m = 15.18 cp
Yg 0.'7
Solucin:
Bo 1.1867
R 460.42
s
6.6
Fiujo MuUi.6tU.ieo en Tubva.6 /falzonta.lu
Z = 0.7786
De la correlacin de Lee:
"'g = 0.01396 cp
w 4.7648 lbm/seg
g
WL 17.6468 lbm/seg
wm 22.4116 lbm/seg
..:1.=0.3172
NReL 29891
N 1 293 548
Reg
5.b.- Los valores de ~,a y b, de acuerdo con las Ecs. (6.7), {6.Sl
y (6.9) son
~ 0.27001
a = 0.2126
b 0.9734
5.c.-
fI = (1293548)0.2126 (29891)9734
fI = 452 598
6.7
De la Fig. 6.1 con ~ y W, se obtiene f tp = 0.005209
El mismo resultado se obtiene con la Ec. (6.14)
6. - De la Ec. (6.4)
JO
\ns 18.4359 lbm/pie 3
0.003179 lb/pg2/pie
(6.19)
Donde:
(6.20)
el de la ordenada:
0.1
y = ] (6.22)
ftp
6.8
- o
2" agua
4" agua
17" destilado
"
(3) .. ,
Fig 6.2 Correlacin de Eaton para el factor de perdidas de
energ(a
6.9
Para obtener las velocidades reales del lquido (v ) y del gas (v ), es
preciso conocer el colgamiento del lquido (H ) enLcualquier partil de
la tubera. Esto es necesario slo cuando lasLprdidas de presin por
aceleracin (Ec. (6.20no son despreciables.
Para: 60000
Para:
Y = tllOO (6.25)
S log d (6.28)
6.10
Para X ~ 819194 - .39981.7 d + 2838.8 d 2 - 73~26 d 3
y .. (21.525 - 1.5~34 d + 0.02278 d2 + 0.00131 d 3) X-o 49 . (6.29)
w
m
8.- Aplicando la Ec. (6.19) calcule el valor de (p/ ,4,L) Y con ste el
valor de AL correspondiente a la .4p sup~esta"
P R B
3 s o
(lb/pq 2 (pie Ibl) (bl/bl) ~
1)
-2
.. 1300 - 50 . 1250
2
lb/pq ..
p Z
"L B
q
a
p e, L.'.
VSL v
sg
6:i1
._ . . __ ._---------~--'--~-------'
Flujo Mu..Ui6cf.6.ic.o en TubvM HM.;On.tale-b
w = 4.7648
9
wL " 17.647
w
m
= 22.412
.\ = 0.3172
(5).- Clculo de a L1 ya L2
NL = 0.15726
N
L1
= (0.15726) (2.199) (1/(45.664) (10.68)3)0.25
NL1 0.022518
NLv =.
1 938 VSL ( \~L/cr)0.25
,
NLv1 = 5.4813
N
gv
11.429
2 (l1.429)0~47569
a Ll = exp ( (-.38011 + (2.34323) (0.022518) ) (5.4813)0.28866
0.4777
6.12
De la misma ~era se ObtieneHL2 '
. NL2 .. 0.022678
N~2 .. 12.0529
HL2 : 0.4673
(6).-'c!lcu10 de
v SL1
1.9669
v =- 4.1174 pies/seg
Ll HL1 0.~777
1-.9539
.. 4.1813 pies/seg
0.4673
"~ 4.1011
v .. '.' = 7.852,0 pie!\l/seg
gl (1 - HL1 0.5223
0.00006131
6.13.
w
De la Fig. 6.2 (_L_ )0.1
w
0.016
m
y 0.016
(7.b) f
tp WL 0.1 = 0.01639.
(--)
W
m
(8) .-
(:) 18.436
\ ns
-1f-
A = (43.539) (22.412)2 (0.01639)
,!,,:,=':=,!,-l.,:.:c.:..;:,:=::;..(_.>.:::.:.::.:=::.!..._ 0.0025 lb/pg2/pie
(6)5 (18.436) (1 - 0.00006131)
6.5.1.- Introduccin.
f w2
= 43.539 Q
tu
,
m'
dS (6.31)
'ns
Definiendo:
wg.
7.2557 fm wm (6.32)
La Ec .. (6.24) se transforma a:
(6.33)
6.14
El factor de friccin se calcula de:
f
ftp = (T) fn (6.34)
n
En donde:
124 d Vm fns
NRi= = /lns (6.36)'
/l ns + /'g(l-.\) (6.37)
-~ s
e
f
(6.38)
n
'ln y
s = -,.0,.0523 + 3.182 lny - 0.8725(lny)2 + 0.01853(lny)4 (6.39)
Donde,
y= -;r-
..\
(6.40)
L
Las tuberas transparentes permitieron observar y clasificar tres pa-
trones de flujo, De acuerdo a esta clasificacin se desarroll el ma-
pa de patrones de flujo, mostrado en la Fig. (6.3), en funcin de
y del nmero de Froude. El patrn de flujo para cualquier condicin
de flujo se puede determinar de este mapa o de la tabla siguiente:
Transicin
6.15
100 XII
...
a::
Z
ti
'Q
-
::J
o
....
IL.
lO
ti
Mapa revisado
'Q
o
.... S{mbolo Rgimen de flujo
ti
E 1 Sec;rec;ado
~::J
Z II Intermitente 1
m Distribuido
:m: Transicin
O.l .....-------'-------'--------'---~----lIII
0.0001 0.001 0.01 0.1
Colgamiento sin resbalamiento, X
6.16
Intermitente 6
Distribuido
En donde:
Ll . 316
0.302
L2 0.0009552 ..\-2.4684
L3 . 0.1
-1.4516
L4 = 0.5 ,\
-6.738
6.17
FluJo Mu.tU6tfUc.c en. Tubelll4 HolLizontai.e.t.
R 460.42
s
B 1.1867 b1/b1
o
z = 0.7786
B .0.009316 pi 3 /pie 3
g
f'L 2.237 cp
f'g 0.01396 cp
C1 10.87 dinas/cm.
(4) 45.763 lb 3
fL ,m/pie
fg
w
g
W
m'
..\. = '0.3172
f ns =18.436 U>m/pie3
(5)
, ,,3
0.15726 "f'L (1/ eL (5 ) 0.25
0.25
(0.15726) (2.237)"'(1 / (45.763) (10.87i 3 )
0.022593
Ct:J /'C1)0.25
NLv '1.938 VsL , \ L'
" C " , _ 0.25
NLV (1.938) .(1.9614) (45.163/10.8,7)
N = 5.4441
Lv
Ngv 1.938 V Sg ( fL/O)0.25
6.18
N (1~938) (4.'2271) ,(45.763/ 10~81}0.25,
'1v
Nqv p.7346
HL = 0.4725
fm .. ~\ HL + (1 - HL )
,'1
e
e m .. (45.763) (0.4725) '+ (0.5275) (5.74b7)
~ = 0.0001092
"L "
A
+
fn 0.017393
.. 0.3172 .. 1.4202
y =
0.47262
$= 0.36662
6.19
~
f = 1 .442.8 .
n
ftp= f n . eS = (0.0173931 (1.4428)= 0.025095
(8)
~
.4L
= 0.003829 lb/pq2/pie
(9)
.4L - (.4p! .4P ! .4 L) =50/0.00382. = 13 058
6.20
Ftu.jo l.fulti.6t!6.i.c.o en rubIM Holzim.taleA
REFERENCIAS CAP. 6
6.21
---------------------------.,---------c---------------,
7.1. - Introduqci6n.
7.1
las prdidas .de presin, Fig. 7.1 (al. Este fenmeno se explica de la -
manera siguiente: :
7.2
Relacin gas-Ilqu ido .Didmetro delconducfo
(a) 4p
Gostode I{quido
Je)
Fig 7.1 Relacin entre la cada de presin y qL ,R Y P, en
tuberas verticales
.So"che ..
.. A:
. '
5 10 50100
7.3
desarrolladas experimentalmente en base a la distribucin de las fases
en la corriente. Por 10 general, de las correlaciones citadas se deter-
mina el valor del colgamiento "HOLD-UP" y de ste, fcilmente pueden ob
tenerse los volmenes de gas y lquido en el flujo.
7.4
100r----------;r---------~-------->--_r----------,
Espuma
Burbuja
---'
Niebla
Bache
0.11...-----....I.--'-_..L-_ _ _......._ _ _-'-....I.'-_ _ _ _ _...
0.1 10 100 1000
7.5
de este tipo estn incluidos los mt01gf de Poettmann y Carpenter (3) , Fa!!.
cher y Brown(4) y Baxendell y Thomas.
7.4.1. Introduccin.
(7.2)
(13)
La correlacin as obtenida fue'extendida por Baxende;l y Thomas ,p~
ra ser aplicab'le a pozos con altos c:astos y flujo por~' el espacio anular.
En la Fig. 7.4 se muestran los resultados de ambas correlaciones.
7.6
0.5
0.1
0.05
0.01
0.005
00020~--:20:-~4~0-~60:-~BO~~I~OO:-~'20
dvPns
_ . - Calculado de datos
de campo
Factor de
friccin
f'oelfmony
Carpenfer
1500<G/L<3000
0.005
0.001 I......I.....L..J...j,..o..r..u.I......JI-I....L.J..i..l.U.l.-...L.......I..I...I.U..u
1 5 10 50 100 500 1000
7.7
Donde:
a .= (7.4)
2.- Calcular para las cOndiciones. inedias del. intervalo i?,.:r (sta temp~
ratura generalmente estimadal, los valores de Z,. Bo, Rs y Bg'
3.- Calcular
..
Pns a P cn la Ec. (5.18}
qo 943.5bl/da d = 1.995 pg
R = 1122 pies3~1 P = 960 lb/pl
T = 153F Vg = 0 . 65
"ro = '0.85
R
s
= 212.19
De las otras cOrrelaciones, se obtiene:
Z = 0.903 B
g
0.01628
qo ' R ,d,WOR
NO
7.9
3.- De la Ee. (5.18):
a = 5.98
Con la Ee. (7;3)
ftp 0.00938
.~~ - ": r
I
, 6
- 5.723 x 10- 4 d JI: to ) + d x 10 6
( 1.848 x 10- 4 ( )2
'\, M qo M
+ 3.5843 x 19-'6 d x 10 6
'\,M
)3 J ] (7.5)
d
d
dpl dh)
qo ~1 ) 144 X 2.9~9 le: 10 5 aS- P
ns
[2
/
x 5.415 x 10- 3 q
op
M
7.10
- (5.723 x 10 -4 ) (d x 10 6 ) _
,(q
-6
(3.5843 x 10) (d xlQ,-) '= O
UZ "
63] " (7.6)
op ",
Simplificando ,
"
d= 3 pg. qo M 275913
Como El1 gasto real es mucho menor que 'el ptimo, se infiere que el flujo
es ineficiente, dEibido a un colgamiento excesivo que incrementa el gra-
dientEl de presin. '""
7.
d = q'oMI, 91970 = 1.073" pg."
Aunque el mtodo es aproximado, permite en el caso analizado, con-
cluir la con.wriieni d usr tuberas de menordiametro, af.in de
prolongar la vida fluyente de los pozos .
7.S.1.-Introduccin.
2
f t2 y,
m
'ns
[t m +
5.362 d
~ 1
<lb = - ' y, y
(7.11)
m Si
1 - '111
~cp
En esta ecuacin se identifican los gradientes, por dlnsidad y' por fri~
cin, as! como el trmino de aceleracin, o sea.
(~) + (~)
<lb e <lb , f
(7.12)
1 -
4;- ,calcular Y S9 ' Y SL , Ym y J , con las Ecs. (5.3), (5.4), (5.5) Y (S.6)
7.13
EjemPlo.- calClil~ elqriadiente de presi6n para las mismas condiciones _
.def1ujo del ejemplo anterior. Considere despreciable el efecto por _
aCttleracian. El valor de p._ a f es de 10cp.
BO . 1.1002 z
Rs ... 2l2.1~
Bq
.. 0.903
= 0.01628
#L
. 3~07S u 13.21
#9 = 0.068.8
Po 49.916
Pg 3.0495
NL .. 0.026256
NLil .. 8.3933
Nqli 20.l3~
Slistittiyendo:
+ fA (1.,.A.)
q
7:14
N
124 d vm P ns (124) (1.995) (10.557) (16.825)
Re /L ns (0.9530)
N
Re
= 46 Hi7
f = 0.021293
n
s = 0.31298
De la Ec. (6.38):
0.31298
f 0.021293 e 0.030919
tp
0.20052 lb/pi/pie
7.15
Para un conducto anular.,
.. d . - d
C1 te!'
E .. +
te
Donde:
(
ci' - rugosidad de
la superficie 'interior de la T.R., i'g.
= 0.00181
( ~ci
Es evidente que para calcular las cadas de presin por las tuberas, -
se seleccionar~ el mtodo .que permita reproducir, con mayor aproxima~
cin, los valores medidos en el campo, de dichas cardas de presin.
Sin embargo generalmente es necesario realizar algunos ajustes para 10-
grr mayor precisin en los resultados. Los valores de las propiedades
de los fluidos deben corregirse, segn se indic en el captulo 1. Si
se observa una desviacin sistemtica en cierto sentido, entre los va-
lores calculados y los medidos, entonces'el ajuste es sencillo. Basta-
r en este caso incluir un factor de eficiencia o modificar el valor
de la rugosidad de la tubera, hasta lograr la mejor aproximacin.
7 .16
Respecto a los patrones de flujo,determin que la correlacin de Beggs
y Brill es aceptable para las tuberas horizontales. Los mapas propues-
tos por Duns y Ros(l) y Aziz(8), fueron los mejores para predecir el
tipo de flujo en los, "risers".
7.17
R ~ F E R'E Ne I AS CAP. ' 7
1.- Duns H. y Ros N.C.J.: "Vertical Flow'of Gas and Liquid Mixtures in
Wells", .l'rqc. 6thWor;Ld petroleum.congress. 1963.
8.- AZiz j{,;> 'Gi>V'iEii G;W, y F,ogarasi M.: "P:-essure Orop in Well P:-odu-
ing Oil andG.as:'JerUl:!'al of caladian' Pt. Tech., julio-septiembre
1972'" '
.a.
11.- GouJ.;el '1' ):.., Tek' M ''"1 Katz D. i..: "Two-phase Flow Through V"'rtiGl t
Im::lind 61: curvedPipll", .JeUJ:Da1 of, P",t. Tech., agosto 1974.
7. 18
a.- FLUJO MULTIFASICO A TRAVES DE ESTRANGULADORES.
a.l.- Introduccin.
a.l
Presi6ri corriente arriba, Ub/pl)
d .-
c
Dimetro del estrangulador, {64 avos de pgl
lores siguientes:
Correlacin A B e
Gilbert 10.0 0.546 1.89
Ros 17.40 0.500 2.00
Baxendell 9.56 0.546 1.9.,
Achong- 3.82 0.650 J .88,
Para que exista flujo crtico se supuso que la presin corriente abajo
debe ser ,al menos de 0.55 de, la presi6ri en la boca del pozO. Bajo es-
tas cOO.diC:iones ei g~o en el ertrangulador es slo funcioo de la pre-'
si6ri COrriente arriPa y _de la ,relacin gas-aceite a condiciones de 'flu-
jo.
La ecu,aci6ri de P Y B es:
qo
1.549 di
73.856 Po +~
1'9 R V/Vl
92 73.6 Pl
(l +'0.5 m}
I 0,.4513
r
Ir + 0.766
+ 0.5663
]ls.2)
Donde:
(8.5)
l_
8.2
Siendo:
[DT l
0.5
+ 460) Zl (R-Rs ) F
ro + 0.000217 Vg R s +WOR P ])
w
1.53 d~ PI
(8.6)
Donde:
14/64 1.1510
16/64 1.0564
20/64 0.9760
8.3
La ecuacin obtenida por los autores es,
(8.71
a = (B o + WOR) - 1/2
(8.8)
k -k--- 0.5
k - 1
0.00504 (~ )(T l + 460lZ 1 (R - Rsl (1 - X + Pl (1 - Xl
~=~~----------------------~~--------~ (8.9)
(T, + 460)Z -l/k
1 + 0.00504
Pl
Donde:
k Cp / Cv
X P2 / P l
Omaa desarroll una correlacin (para flujo crticol entre el gasto, la pr~
Slon corriente arriba del estrangulador, la relacin gas lquido, la
densidad de los lquidos y el tamao del orificio. Dicha correlacin se
obtuvo a partir de datos experimentales.-En vista de que estos datos estu
vieron dentro de rangos muy limitados, su aplicacin slo se recomienda -
para orificios hasta de 14/64 de pg. y gastos mximos de 800 bl/da.
(8.10)
Donde:
N (8.11)
q
N (8.12)
8.4
0.0174 Pl
N = (8.13)
P
(.1' u10. 5
L
1 (8.14)
Q
(R- Rsl Bg
1 +
5.615 Bo
(8.15)
8.5
REFERENCIAS CAP. 8
1.- Gilbert, W.E.: "Flowing ond Gas Lift Well Performance". Dril1.
and Prado Pract. 1954.
7. - Omaa R.: "Mu1tiphase F10w Through Chokes". SPE 2682. 44th Fall
Meeting. Sep-Oct. 1969.
8.6
9.- COMPORTAt4IENTO DE POZOS FLUYENTES
9.1.- Introduccin.
El mtodo que se ver mas adelante pennite obtener .las curvas de compor
tamiento del sisteIl'A constituido por el yacimiento, la tubera verti~al.
la tubera horizontal y el estr~~gulador. Estas curvas (Fig. 9.3) repre-
sentan la variacin de la presin en los nodos y la cada de presin en
los elementos del sistema, con el gasto.
Se observa en esta figura que elgasto mximo que puede obtenerse corres
ponde al del flujo sin estrangulador (punto B). se aprecia tambien que -
al estrangular el pozo el gasto disminuye; sin embargo la presin en
la boca decrece (en vez de continuar amnentando) a partir del. punto A,
en el que se obtiene el gasto ptimo. Tambin se nota el incremento c~
siderable en la cada de presin por la T.P., provOcado, como se indi-
c, p,?r el colgamiento del lquido.
(9.1)
9.1
Aceite
Nodo
<D Separador
o Estrangulador
o Boco del pozo
@ FORdadel pozo
Radio de drene
/fr.:,.,..,..~~.,.,
9.2
Compo~ento de Pozo~ Ftuyent~
(9.2)
qL
Pwf p (9.4)
'ws K
.toi B
ro oi
J. (1 + WOR)
1 K B
roi l' o o
ES evidente que para resolver la Ec. (9.4) se requiere qe la prediccin
del comportamiento del yacimiento, ya que los valores d Kro y de WOR -
son funcin de la saturacin de aceite yde la sat~racin de agua.
En dicha figura se observa que cuando Pth = Pe ,se tiene el gasto mxi-
mo, correspondiente al flujo sin estrangulador. Los gastos inferiores se
pueden obtener usando estranguladores en el cabezal del pozo. El tamao
del estrangulador se puede calcular mediante las ecuaciones presentadas
en el captulo anterior. '
9.3
Compo~ento de Paza~ Ftuyenteh.
lb
p(=r)
pg
1 - = : : : : : - - - - - - - - - - - - - - - - - __
"
"
__ - Pe
q (bl/dia)
'\nax o
9. 4
CompoJr.t:am.nto de POZ04 FR.uye.ntu
Fig. 9.4.- GASTO MAXIMO OBTENIBLE EN FUNCION DEL/\. PRESION ESTATICA DEL
POZO.
9.5
..
! "1.... '~~,
. '".;.l
.J
9.6 ..
Compo4tam~ento de Pozo~ Ftuyent~
La Fig. 9.6 muestra la variaci6n del gasto mximo al usar lmeas de des-
carga de diferentes dimetros. Se observa que para una tuber!a de produc
ci6n dada, existe un dimetro de lmea de descarga para el cual se ob- -
tiene el mximo gasto. Incrementos adicionales en el dimetro de la l-
nea de descarga ya no proporcionan mayor producci6n.
T.
Pws
P.}. CONSTANTES
I
GASTO MAXIMO
9. 7
CompolLtamiento de POZO.6 Flu.yentu.
L.O.} CONSTANTES
Pws
9.8
- gasto ptimo de aceite, a partir del cual la presin de fondo aumen-
t?- al disminuir el gasto.
2. - Todos los. mtodos indican que los gastos ptiIOOs decrecen al dis-
minuir el dimetro de la tubera de
.
produccin. ,
3. - Los valores d.e los gastos ptimos son dife~nte.s para cada corre-
lacin.
9.5.-0ptimizacinde un sistemadeproducC!n.
I .
9.9
d) Cambiar la tubera de produccin.
e) Cambiar la lnea de descarga o instalar ~a lnea adicional.
f) Instalar un sistema artificial de produccin.
9.10
REFSRENCI.AS CAP. 9
2.:", Fetkovich, 'M. ,J. "Isoeromi.l., Testing of Oil ,Wells,n. SPE ~ 4!29.
48 th FaIl, Meeting of IM!;;, 1973.
9.11
10.- CALCULO DE I:,A DISTRIBUCION DE LA TEMPERATURA EN TUBERIAS.
10.1. Introduccin.
I I
~mumvH/ffrH#n/nrlffH/UW//AiH/Il7Ullj'll?WJ7//lWl/d-:
I I '
1
) T dQ T - dT
,
I I
Ta~ dx ~Ta
I I
10-1
'O.(J,tJt.ibuc6n de .ea. Tempvr.a.twr.a. en TubVll.6
En donde:
(10.2)
Que es' la prdiJa de calor del flu~do al pasar de (x) a (x + dx); donde:
Agrupando:
dT
T - Ta
=- a ale (10.3)
Donde:
a= lI' d U / 12 Wf C f (10.4)
a = (10.5)
Integrando:
In := - ax (10.6)
10.2
Despejando Tlt, queda:'
Donde:
Eje1l1Plo.- Se de'sea calc~lar li.. ,tE!1l1Peratura del gas, en :un:.' gasdduc:o ente-
rrado, del que se tiene la sigu~entE;l informacin:
d = 12pq
Cfg = 0.8 BtU/lb, - F
10.3
0.40
......
N'
. p =100lb/pg2
01
Q.
......
.........
.Q
lI...
..!..."
c:
O
f
O
.s::.
0.28
....
CII
:; 0~24
O
'")
-CII
c:
CII
"O
0.20
;:
'!
(,) 0.16
OJ2
2000;,...
10.4
Solucioo.-
W= '0.0764 ~ Pg
W = 2.4 x 10 6 lb l<tia
'. m .
, ;a.= +,.5708. x 10
-5
.1'.
- 5.4288
'1'
'x
= 7l.68F
La ~olucin rigurosa de. la Ec. (10.8) cOmprende la determinacin previa
de", de ac:uerd~,con; (8) .s H
.5'1' . ' 'OS ,)'1'
" = (-P)a /J.H (1,0.9),
,( /J'1' lp
.' .
DOnde H es la 'entalpia esp~cfic~, cuyas valores se. calculan mediante la
aplicacin de )lila .ecul/-cin, de estado, a diferentes presiones y t~mpElra-'
. turas; Las deri vades de la er;talpia 'pon l; prsin y ,cenla: t.emperatura
se abtienendel calculo anterior. ColIio H yen' consecuencia" dependen de
la cOIlIposicin del fluido, laprsinyla temperatUra, para res.olver
la J:c. (10.8) esnecesariotilizarunproceso iterativo, ya, que
dp/dXdepende 'taiibil de'. la temperatura: El proceso iterativo consiste en:
10. 5
al - A pa:tir de, un punto d~, pre,si6n ~ temperatura conocidos ,(Pl' '1') r,
Cl9Jlsl.derar un incremento de lonql.tud A x y suponer una calda ,~ ..;
, presi6n y una oalda, de temperatura.
a=0.2618Ud/(W
, ' q Cf q + W,o Cf o + Ww C~)
...,. (10.10)
En esta ecuaciSn Cfq' Cfo y Cfw' son ios calores especficos del qas, -
del aeite y del agua. Los qastos mSsic:qs pueden obtenerse' ~n:
W (R-Rs)Bq (10.11)
q P9 io
"'o 5.615 Po io B
o
(10.12)
Ww = 5.615 Pw ~ Bw (10.13)
10.6
21.10
o
---------~-
0.50
0.40 .
30 50 100 150 250 300 350
Temperatura (OF)
Fig 10.3 Calor epeefico del gas nafural con densidad relativa de
0~60 0.7 '
10.7
V.i.WbucUn de. la. TempeJUt.twuz en TubeILf.aA
Donde:
(10.151
y
,':.
~=l +
U h
La 'funcin del tiempO, f(tl puede calcularse, para tiempos menores, de, 4PO
dias, con:
lag f(t) = 0.31333 lag Y - 0.06 (lag y)2 + 0.006666 (lag n,3
Donde:
y 552 t / d ce 2 ~(lO.::&)
Donde: ,
F 1 + B fet)
tlO.,~O)
2 B
B
zti U/~e , (lO.j
10.8
i gasto de inyeccin, bl/dia
t tiempo, dias
Tbh temperatura en el fondo del pozo, F
10.9
Para pozos inyecto'res a travs de la T.R. (sin T.P.), ~e d ci ' por lo
que:
(10.27)
Solucin:
f(t) = 3.5732
De la Ec. (10.26)
T 163F
Los pozos productores con gastos altos pueden sfrir elongaciones consi
derab1es en sus T. P. Tambin es importante estimar l~ temperatura en la
boca del pozo para definir el procesamiento de los fluidos producidos.
10.1.0
V.utJbuc.i.n de .ea. TempelU1.tww. en TubeJL.C,a.6
SoluCin:
U= 45.50/(dc i - d te )
33.6
~e
q M
o
1.896
Se tiene:
qoM C f (dci - d te ) , (33.6 + 1.896 d f(t) / ( dci - d te ) ) (10.32)
A =
400.27 dt i
'El valor de f(t) de acuerdo con las &cs. (10.17) y (10.18) ,.para t = 30
das y d ce = 9.626 pq, resulta:
f(t) = 3.0025
10.11
Sustituyendo valores en la Ec. (10.32):
A = 118834.7
Tth 205.6 F
da:
Tth = 142.8 F
k' (Z - ~)
Donde:
El valor de k', para po~oS productores en Arabia con una relacin gas-
aceite de 540 pies 3 /bl, result de 0.0015 F/pie; sin embargo para al-
tas relaciones gas aceite puede ser mayor dicho valor.
10.12
.f),l6Wbuc,i;n de. .ea. Te.mpe.J!.a;tWr.a. e.n Tubef.it4
1 1
-u=:- +.2!... ~+ h1 (10.34)
l)C+ . l)s o
Donde~
I
Generalmente h. puede despreciarse. CUando el recubrimiento de concreto
es mayor de. 4 pg., tambin son despreciables Kns y he. La conductividad
t~rmicadel concreto puede variar substanciall\lente con el contenido de
hUlnedad y su integridad estructuraL La conductividad trmica: d la pe-
lc.ula exterior es funcin de la. velocidad dl fluido ambiental, normal
a la tuberla. . . .
En esta ecuacin:
lO.
------ -------c~_:_-----------------------,
124d VfPf
R
e = (10.;17)
"'f
Pr '" 58.06 (10.38)
. SoluCin,
Pw = 0.65 C.p.
K=
w S.4 Btu/di;t
' , ' -pie - V , - - -lb
Pf , '"' 61.5 - ...
, , 3
pie
10,14
V~tJbutU611 de .ea. Tempvr.a.twta en TubelCt
1 _1_ x~
h
+ -k-
U
o hs
Sustituyendo:
_1_ 1 0.02604
U 3838 + 600
Donde los valores de Nu ' Rey Pr, corresponden al del flujo en el int~
rior de la tubera.
1 ID (48Z/d)
(10.41)
U 24 ~e/d
Donde:
10.15
Xr espesordelrecubriti\iento, pies
(Btu/dia pie F)
Arenoso Seco 10 ~ 20
Arenoso HIedo 22 ~ 29
~
Arenoso Saturado 53 - 63
Arcilloso Seco 9 - 14
Arcilloso HU1l\edo 19 24
Arcilloso Saturado 29 - 44
Po!iestireno 0 . 54
.Concreto 24.0.0
10.16
V-ibbuc6n de .e.a Tempvr.a:tuJta. en TubVlU.
REFERENCIAS: CAP. 10
1.- Coulter, D. M.: et al.: "Revised Equation Improves Flowing Gas Tem-
perature Prediction". Oil and Gas Journal, Feb. 26, 1979.
2.- Hein .M.: "Here are Methods for Sizing Offshore Pipelines". Oil and
Gas Journa1, Mayo 2, 1983.
4.- Ramey, H.J. Jr.: "Wellbore Heat Transmission", Trans. AIME (1962)
225.
12.- Lawton, L.L. : "Curves give Temperature Drop for Expanding Gases".
Wor1d Oi1, Enero 1984.
10.17
------------------------------.."
NObIENCI;A'rtlRA ,
2
A rea, pg
C compresibilidlid, pg2/lb
diiimetr~, hidrulico, pg
dimetro exterior de la T.P., pg
profundidad, pies
ge grldintegeotrmico, F/pie
L loriqitud, millas
L. o: lnea de descarq~
m masa, lbm
iI
nilmero,de la viscosidad del lquido, adimensional
P presin, lb/pg2
q gast, bl/da
Ili
.-~-----~-~----------------- --------
T temperatura. F
v velocidad, pie!seg
rugosidad,pg
2
~ coeficiente de Joule-Thompson, P/lbf/pg
J
densidad del gas, lbm/pie
/,
'L
viscosida4 del 'l!quido, :.:p
viscosidad de la mezcla, cp
.1
SUBINDlCES
aire
burbujeo
cs condiciones atmosfricas
f friccin
g gas
L lquido
m masa
m mezcla
, ns sin resbalamiento
o aceite
s separacin
w agua
o"
CONSTANTES Y FACTORES DE CONVERSION
Constantes
Conversiones
Longitud Densidad
1 pg 2.54 cm
141. 5
1 pie 30.48 cm
~'a = 131. 5 + API
1 minil 5277 pies
141.5
Presin - 131. 5
~
1 lb 453.59. gr
1 kg =2.2046 lb
ViScosidad
1 cp 0.04031 1bm/pie~in
1 cp = 6 .. 7197 x 10-4 lb m/ pie-seg
,..---~~~~~~~~~~~~~~~~~~--~~~-~ ~~~~