Flujo Multifásico en Tuberías
Flujo Multifásico en Tuberías
Flujo Multifásico en Tuberías
FACULTAD DE INGENIERA
DIVISIN EN CIENCIAS DE LA TIERRA
Tesis:
Flujo multifsico en tuberas
Ingeniero Petrolero
Presentan:
Torres Coria Juan Manuel
Trauwvitz Rosas Edgar Fernando
NDICE
INTRODUCCIN.5
CAPTULO 1: PROPIEDADES DE LOS FLUIDOS
1.1 Introduccin..7
1.2 Propiedades del aceite saturado..15
1.2.1 Correlacin de Standing...15
1.2.2 Correlacin de Vzquez...16
1.2.3 Correlacin de Oistein..18
1.2.4 Correlacin de J. A. Lasater.22
1.2.5 Densidad del aceite saturado23
1.2.6 Viscosidad del aceite saturado.24
1.2.7 Tensin superficial del aceite saturado24
1.3 Propiedades del aceite bajosaturado24
1.3.1 Compresibilidad...24
1.3.2 Densidad del aceite bajosaturado.24
1.3.3 Viscosidad del aceite bajosaturado..25
1.3.4 Factor de volumen del aceite bajosaturado..25
1.3.5 Presin de burbuja25
1.4. Propiedades del gas natural25
1.4.1 Densidad relativa del gas..25
1.4.2 Factor de volumen del gas27
1.4.3 Densidad del gas...28
1.4.4 Factor de compresibilidad28
1.4.5 Propiedades de gases que contienen N2, CO2 y H2S33
1.4.6 Viscosidad del gas34
1.4.7 Correccin de la viscosidad del gas por presencia de gases contaminantes34
1.5 Propiedades del agua saturada34
1.5.1 Factor de volumen35
1.5.2 Densidad...37
1.5.3 Viscosidad37
1.5.4 Tensin superficial agua-gas39
1.5.5 Solubilidad del gas en agua..39
1.6 Propiedades del agua bajosaturada..39
1.6.1 Compresibilidad...40
1.6.2 Factor de volumen40
1.7 Ejemplos..41
CAPTULO 2: FUNDAMENTOS DE FLUJO MULTIFSICO
2.1Introduccin.53
2.2 Ecuaciones fundamentales..53
2.2.1 Ecuacin General de Energa...54
2.2.2 Prdidas de presin por friccin...56
2.2.3 Flujo de lquidos por tuberas...59
2.2.3.1 Ecuacin general59
2.2.3.2 Nmero de Reynolds para el lquido.61
2.2.3.3 Eficiencia de flujo..61
2.2.4 Flujo de gas por tuberas..62
2.2.4.1 Ecuacin general62
2.2.4.2 Nmero de Reynolds para el gas...65
2.2.4.3 Eficiencia de flujo..66
2.3 Colgamiento de lquido...66
2.3.1 Colgamiento sin resbalamiento70
2.3.2 Velocidades superficiales.70
2.3.3 Velocidad real...71
2.3.4 Densidad de la mezcla de los fluidos...71
2.3.5 Gasto de masa...73
2.3.6 Viscosidad de la mezcla.......74
2.3.7 Tensin superficial de la mezcla de lquidos74
2.3.8 Densidad de la mezcla de lquidos...74
2.4 Patrones de flujo..75
2.4.1 Patrones de flujo en tuberas horizontales75
2.4.1.1 Correlacin de Taitel y Dukler para predecir la transicin de los
patrones de flujo77
2.4.2 Patrones de flujo en tuberas verticales79
2.5 Ejemplos..80
CAPTULO 3: FLUJO MULTIFSICO EN TUBERAS VERTICALES
3.1Introduccin.83
3.2 Correlaciones..84
3.2.1 Poettman y Carpenter...84
3.2.2 Duns y Ros...88
3.2.3 Hagedorn y Brown...97
3.2.4 Orkiszewski104
3.2.5 Beggs y Brill...112
3.3 Ejemplos119
3.3.1 Mtodo Poettman y Carpenter119
3.3.2 Mtodo Duns y Ros122
3.3.3 Mtodo Hagedorn y Brown126
3.3.4 Mtodo Orkiszewski...132
3.3.5 Mtodo Beggs y Brill.137
CAPTULO 4: FLUJO MULTIFSICO EN TUBERAS HORIZONTALES
4.1 Introduccin..143
4.2 Correlaciones.144
4.2.1 Lockhart y Martinelli..144
4.2.2 Baker..147
4.2.3 Dukler.152
4.2.4 Eaton...158
4.2.5 Beggs y Brill...161
4.3 Ejemplos167
4.3.1 Mtodo Lockhart y Martinelli167
4.2.2 Mtodo Baker.169
4.3.3 Mtodo Dukler...172
4.3.4 Mtodo Eaton.177
4.3.5 Mtodo Beggs y Brill.183
CAPTULO 5: FLUJO MULTIFSICO EN TUBERAS INCLINADAS
5.1 Introduccin..187
5.2 Correlaciones187
5.2.1 Correlacinde Flanigan..187
5.2.2 Correlacin de Beggs y Brill..190
5.3 Ejemplos195
5.3.1 Mtodo Flanigan195
5.3.2 Mtodo Beggs y Brill.197
5.4 Modelos Mecansticos..201
5.4.1 Modelo para predecir la transicin de los patrones de flujo segn Xiao...202
5.4.2 Modelo para predecir la transicin de los patrones de flujo segn Barnea204
5.4.3 Caracterizacin hidrodinmica de los patrones de flujo208
5.4.3.1 Modelo para flujo estratificado208
5.4.3.2 Modelo para flujo intermitente213
5.4.3.3 Modelo para flujo anular.219
5.4.3.4 Modelo para flujo burbuja dispersa.223
CAPTULO 6: FLUJO EN ESTRANGULADORES
6.1 Conceptos fundamentales..225
6.1.1 Estrangulador..225
6.1.2 Flujo crtico y subcrtico.225
6.1.3 Coeficiente de descarga..227
6.2 Modelos para gas...228
6.3 Modelos multifsicos....230
6.3.1 Modelos para flujo crtico..231
6.3.1.1 Correlacin de Gilbert, Ros y Achong....231
6.3.1.2 Frmula de Ros (adaptacin de Poettman y Beck).232
6.3.1.3 Correlacin de Omaa R.236
6.3.1.4 Ecuacin de Ashford...238
6.3.1.5 Velocidades snicas de Nguyen..238
6.3.1.6 Velocidad snica de Wallis.239
6.3.1.7 Ashford y Pierce..239
6.3.1.8 Sachdeva..240
Captulo 1
CAPTULO 1:
PROPIEDADES DE LOS FLUIDOS
1.1 Introduccin
La evaluacin de propiedades fsicas de mezclas de hidrocarburos es un importante
paso para el diseo de diferentes etapas de operacin de un campo petrolero.
Las propiedades de los fluidos varan conforme cambia la presin y la temperatura, la
mayora de estas propiedades deben ser evaluadas por la ingeniera de yacimientos e
ingeniera de produccin.
El clculo del gradiente de presin en un sistema de tuberas con dos o ms fases
fluyendo, requiere de la prediccin de propiedades de los fluidos como el gas disuelto, los
factores de formacin, compresibilidad del aceite, la viscosidad en distintos puntos de la
tubera. Incluso las mediciones en laboratorio para estas propiedades pueden estar
disponibles como funcin de la presin, usualmente dichas medidas son realizadas a
temperatura del yacimiento. El comportamiento de estas propiedades cambia con la
temperatura y son predichas empricamente.
A continuacin daremos algunas definiciones de las propiedades de los fluidos, as
como algunos conceptos bsicos empleados en relacin con dichas propiedades.
Aceite Estabilizado.- Es el aceite que ha sido sometido a un proceso de separacin con el
objeto de ajustar su presin de vapor y reducir su vaporizacin al quedar expuesto,
posteriormente, a las condiciones atmosfricas.
Aceite residual.- Es el lquido que permanece en la celda pVT al terminar un proceso de
separacin en el laboratorio. Generalmente, el aceite residual se determina a condiciones
estndar.
Aceite en el tanque de almacenamiento.- Es el lquido que resulta de la produccin de los
hidrocarburos de un yacimiento a travs del equipo superficial empleado para separar los
componentes gaseosos. Las propiedades y la composicin del aceite dependen de las
condiciones de separacin empleadas, como son: nmero de etapas de separacin,
presiones y temperaturas. El aceite en el tanque de almacenamiento se acostumbra a
reportarlo a condiciones estndar.
Anlisis PVT.- Son requisito indispensable para conocer las propiedades de los fluidos,
consiste en simular en el laboratorio el comportamiento de los fluidos en el yacimiento a
temperatura constante. Estos estudios son absolutamente necesarios para llevar a cabo
actividades de ingeniera de yacimientos, anlisis nodales y diseo de instalaciones de
produccin.
Para tener la certeza de que el muestreo es representativo, se hace una validacin
exhaustiva tomando en cuenta todos los parmetros del yacimiento medidos durante la
toma de muestras como son:
Captulo 1
Del anlisis PVT se obtiene el comportamiento de los fluidos mediante grficas que
nos ayudan a determinar las propiedades de los fluidos.
Condiciones Estndar.- Son definidas por los reglamentos de pases o estados; en Mxico
las condiciones estndar son de P = 14.69 lb/pg2 (101.33 kPa, 1.033kg/cm2, 1 atm.) y T =
60 oF (15.56 o C, 288.71 K)
Correlacin.- Es un conjunto de datos experimentales obtenidos mediante parmetros de
campo o de laboratorio al que se le ajusta una ecuacin o un conjunto de ecuaciones
empricas o experimentales, (tambin se pueden presentar en forma grfica) la cual
representar un determinado fenmeno.
Se aplica para obtener una aproximacin de las propiedades de los fluidos sin tener
que realizar un anlisis PVT o mediciones directas
Densidad del aceite.- En la prctica comn, el aceite crudo se clasifica de acuerdo a su peso
especfico, y es expresado en una escala normalizada por el Instituto Estadounidense del
Petrleo (American Petroleum Institute). Esta escala es llamada densidad API, o
comnmente conocida como grados API (oAPI).
La densidad en oAPI esta definida como:
o
API =
141.5
131.5
(1.1)
donde:
o = densidad relativa del aceite a condiciones estndar o aceite muerto.
En la tabla 1.1 se muestra una clasificacin del aceite crudo en trminos de su
densidad, aunque hay que enfatizar que la produccin de diversos yacimientos productores
est asociada a un solo tipo de aceite crudo.
Captulo 1
p V = Z n R (T + 460)
(1.2)
Donde:
Captulo 1
Z = Factor de compresibilidad
R = Constante universal de los gases, 10.71 (pies/mole-oR)
n = Nmero de moles de gas, (lb-mole)
p = Presin, psia
V = Volumen (pies3)
T = Temperatura, oF
10
Captulo 1
Bo =
(Vo )p ,T
(Vo )c.s.
(1.3)
Bg =
(V )
(V )
g p ,T
(1.4)
g c.s .
Figura 1.4. Comportamiento tpico del factor de volumen del gas en funcin de la
presin del yacimiento a temperatura constante.
Factor de volumen del agua (Bw).- Se define como el volumen que ocupa un barril de agua
medido a condiciones de yacimiento por cada barril de agua medido a condiciones estndar.
La expresin matemtica del Bw es:
Bw =
(Vw )p ,T
(Vw )c.s.
(1.5)
Factor de volumen total (Bt).- Es tambin llamado factor de volumen de las dos fases, es
decir, cuando la presin est por debajo de la presin de burbujeo. El Bt, a diferencia del
Bo, considera tambin el gas que se liber en el yacimiento al bajar la presin por debajo de
la presin de burbujeo.
La ecuacin que representa al Bt es:
11
Captulo 1
B t = Bo + Bg (R sb R s )
(1.6)
Donde:
Bo = factor de volumen del aceite a una presin menor a la pb, en [email protected]. / [email protected]
Bg = factor de volumen del gas, en [email protected]. / [email protected].
Rsb = relacin de solubilidad a la presin de burbujeo, en [email protected]. / [email protected].
Rs = relacin de solubilidad a una presin menor a la pb, en [email protected]. / [email protected].
Al analizar la ecuacin podemos observar que al restar la Rsb de la Rs nos da la
cantidad de gas liberado en el yacimiento debido a la disminucin de la presin, y se
multiplica por el factor de volumen del gas para que ese gas liberado este a condiciones de
yacimiento.
La figura 1.5 nos muestra la comparacin del factor de volumen total con el factor de
volumen del aceite. Los dos factores de volumen son idnticos cuando se encuentran arriba
de la presin de burbujeo, ya que no se ha liberado gas todava; la diferencia ocurre cuando
la presin disminuye por debajo de la pb, y es debido a que el Bt si considera el gas que ha
sido liberado del aceite en el yacimiento, por lo cual el Bt aumenta.
Figura 1.5. Comportamiento del factor de volumen total y del aceite en funcin de la
presin del yacimiento.
Mole.- Es el peso molecular de cualquier sustancia. Por ejemplo el peso molecular del
metano es de 16.04 lb = 1 lb-mole. En igual forma 16.04 gr = 1 gr-mole del metano. Una
mole-lb de un gas ocupa 379 pies3 a condiciones estndar.
Presin de vapor.- Es la presin que ejerce el vapor de una sustancia cuando sta y su
vapor estn en equilibrio. El equilibrio se establece cuando el ritmo de evaporacin de la
sustancia es igual al ritmo de condensacin de su vapor.
Presin de vapor Reid.- Es la medida de la presin de vapor de una muestra a 100 oF, en
una celda especial, y siguiendo la norma de evaluacin as denominada.
12
Captulo 1
13
Captulo 1
Viscosidad del aceite.- La viscosidad es una medida de la resistencia al flujo ejercida por
un fluido, y sus unidades usualmente son los centipoise ( grm / 100 seg cm ) .
La viscosidad, como otras propiedades fsicas de los lquidos, es afectada por la
presin y la temperatura. Un incremento en la temperatura causa una disminucin en la
viscosidad. Una disminucin en la presin ocasiona una disminucin en la viscosidad, con
la condicin de que el nico efecto de la presin es comprimir el lquido. En los
yacimientos de aceite la viscosidad tambin se ve afectada, debido a que una disminucin
en la cantidad de gas en solucin en el aceite causa un aumento de la viscosidad de aceite, y
por supuesto, la cantidad de gas en solucin esta en funcin directa de la presin.
Figura 1.7. Forma tpica del comportamiento de la viscosidad del aceite en funcin de la
presin, a temperatura de yacimiento constante.
En la grfica de la figura 1.7 se puede observar como la viscosidad del gas disminuye
si la presin del yacimiento disminuye, esto es debido a que como la presin es menor las
molculas se encuentran ms separadas entre s por lo que se pueden mover ms fcilmente
unas con otras. Tambin se puede observar el efecto que tiene la temperatura sobre la
14
Captulo 1
viscosidad (figura 1.8), a presiones bajas la viscosidad del gas aumenta conforme aumenta
la temperatura aumenta, sin embargo, a altas presiones la viscosidad disminuye conforme la
temperatura se incrementa.
gd
0.83
(1.7)
R s = gd
18
(1.8)
El factor de volumen del aceite fue relacionado con la relacin gas disuelto-aceite,
temperatura, densidad relativa del gas y densidad del aceite.
El factor de volumen por debajo de la pb (Bo) obtenido fue:
B o = 0.972 + 0.000147(F)1.175
(1.9)
Donde:
gd
F = R s
o
0.5
+ 1.25 T
(1.10)
15
Captulo 1
(1.11)
Haciendo nfasis de que hay veces que se requieren propiedades fsicas de los fluidos
del yacimiento que slo se obtendrn en un laboratorio de anlisis pVT. Muchas veces estas
propiedades son requeridas cuando solamente se tiene disponible la densidad del aceite y
gas, adems de presin y temperatura del yacimiento. El propsito de este estudio fue la de
utilizar una extensa base de datos pVT para el desarrollo emprico de correlaciones.
Muchas de las correlaciones usadas fueron desarrolladas haces muchos aos con un nmero
limitado de datos y fueron utilizados mas all del rango para el que fueron intentados. Estos
estudios fueron realizados por la reunin de ms de 600 anlisis PVT de campos de todo el
mundo.
En el desarrollo de la correlacin se encontr que todas las propiedades de los fluidos
podan ser correlacionadas como funcin de la presin, temperatura, densidad del gas y
densidad del aceite. La densidad del gas es un fuerte parmetro a correlacionar y,
desafortunadamente, es una de las variables mas cuestionadas por su veracidad, porque
depende de las condiciones a las que se separa el gas del aceite. Por dicha razn, se decidi
usar la densidad del gas obtenida de condiciones particulares de presin de separacin en
todas las correlaciones que requieren de la densidad del gas. El valor de 100 psig se utiliz
como referencia de presin por ser la presin mnima a la que el aceite se encoge por una
prueba de separacin disponible. Tambin se consider la presin de 0 psig, pero el valor
de 100 psig representa ms cercanamente la condicin del yacimiento y se tendr menos
error s la presin real del separador se obtiene de una muestra de gas.
Por lo tanto, es el objetivo de este estudio utilizar las dems correlaciones, la
densidad del gas se debe corregir primero de un resultado de separacin a 100 psig. El
anlisis de regresin dio como resultado la siguiente ecuacin para corregir la densidad del
gas para las condiciones del separador.
p
gs = gp 1 + 5.912 10 5 o API Ts log s
114.7
(1.12)
16
C3 API
(
T
460
)
+
(1.13)
Captulo 1
- C3 API C 2
T + 460
R e
pb =
C1 gs
(1.14)
Los valores de los coeficientes en las ecuaciones 1.13 y 1.14 se obtienen de la tabla
1.3.
API 300
0.0362
1.0937
25.724
0.0178
1.187
23.931
La expresin que se obtuvo para determinar el factor de volumen para aceite saturado
fue una correlacin como funcin del gas disuelto, temperatura, densidad del aceite y
densidad del gas:
( p pb )
o API
o
+ C 3 R s (T 60) API
B o = 1 + C1 R s + C 2 (T 60)
gs
gs
(1.15)
Los valores de los coeficientes para Bo (ecuacin 1.15) se obtienen de la tabla 1.4.
Tabla 1.4. Coeficientes de Vzquez para Bo.
0
0
Coeficiente
API 300
API > 300
4.67710-4
4.6710-4
C1
-5
1.75110
1.110-5
C2
-1.81110-8
1.33710-9
C3
(1.16)
a 1 + a 2 R s + a 3 T + a 4 gs + a 5 o API
p 10 5
(1.17)
17
Captulo 1
Donde:
a1 = -1443
a2 = 5
a3 = 17.2
a4 = -1180
a5 = 12.61
(1.18)
Donde:
m = C1 p C2 e (C3 +C4 P )
(1.19)
y:
C1 = 2.6
C2 = 1.187
C3 = -11.513
C 4 = 8.98 10 5
Temperatura
18
Captulo 1
p b = f R, g , o API, T
(1.20)
0.816
T 0.172
o
API 0.989
(1.21)
Usando esta ecuacin, se graficaron (escala log-log) los datos obtenidos de cada una
de las muestras, obteniendo una curva parablica a la cual se le ajust una ecuacin por
anlisis de regresin, y con la cual podemos calcular la presin de burbujeo. La ecuacin es
la siguiente:
p b = 10
b
b
(1.22)
Donde:
R
pb =
g
*
0.816
Tn
0.989
API
(1.23)
19
Captulo 1
(p b )c = p b C N
C CO 2 C H 2S
(1.24)
[(
(1.25)
Y:
a5 = 1.954 10 11
a1 = 2.65 104
a2 = 5.5 103
a3 = 0.0931
a4 = 0.8295
a6 = 4.699
a7 = 0.027
a8 = 2.366
C CO 2 = 1.0 693.8 y CO 2 T 1.553
(1.26)
)(
(1.27)
(1.28)
Donde:
a = 3.184 11 T 3.444
(1.29)
b = [10.213(log T ) 36.447]
(1.30)
Donde:
20
API corr =
API* corr o
API
o
API*
(1.31)
Captulo 1
o API0.989
R s = g
p b
n
(1.32)
p b = 10[2.8869(14.18113.3093log p ) ]
(1.33)
a
Bob = f 2 R g + b T
(1.34)
Bob = 1 + 10 A
(1.35)
21
Captulo 1
Donde:
(1.36)
y:
B ob
g
= R
o
0.526
+ 0.968
(1.37)
T 0.5
2.910 0.00027 R
p 1.1089
B t = k 4 R 0.3 o
(1.38)
Bt = 10 A
(1.39)
Donde:
(1.40)
y:
T 0.5 2.9100.00027R 1.1089
p
B t = R 0.3 o
g
(1.41)
22
p f (T + 460 )
gd
(1.42)
Captulo 1
(1.43)
Rs
379.3
yg =
Rs
350 o
+
o
379.3
(1.44)
El peso molecular del aceite en el tanque (Mo) se correlacion con los API del aceite
en el tanque de almacenamiento, a cuya curva se le ajustaron las siguientes expresiones.
Si 15 0API <40,
(1.45)
Mo = (1048.33 / 0API)1.6736
(1.46)
Si 40 0API <55,
yg o
R s = 132755
(1 y )M
g
o
(1.47)
o =
62.4 o + 0.01362 Rs gd
Bo
(1.49)
Donde:
o =
141.5
131.5+ oAPI
(1.50)
23
Captulo 1
o = a( om ) b
(1.51)
Donde:
a = 10.715(R s + 100) -0.515
(1.52)
(1.53)
om = 10 x 1
(1.54)
X = Y T -1.163
(1.55)
Y = 10 z
(1.56)
Donde:
(1.57)
(1.58)
La compresibilidad del aceite bajosaturado puede ser calculada con la ecuacin 1.17
de la correlacin de Vzquez.
1.3.2 Densidad del aceite bajosaturado
o = ob e( c (p p
o
Donde:
24
b ))
(1.59)
Captulo 1
ob =
62.4 o + 0.01362 R gd
Bob
(1.60)
(1.61)
Donde:
(1.62)
(1.63)
En los clculos de las propiedades de los fluidos se utilizan tres tipos de densidades
relativas del gas, por lo que es conveniente distinguirlas.
25
Captulo 1
g =
i =1
gi
q gi
(1.64)
qg
Donde:
n = nmero de etapas de separacin
gi = densidad relativa del gas en la salida del separador i
qgi = gasto en la salida del separador i (pies3g a c.s. / da)
La densidad relativa del gas disuelto puede obtenerse con la correlacin de Katz.
(1.65)
gd gf 0.55
El valor de la densidad relativa del gas libre gf puede obtenerse de un balance
msico:
w g = w gd + w gf
(1.66)
Donde:
wg =
0.0764 R q o g
86400
(1.67)
lb g
3
pie c.s.
lb aire pie g c.s.
bl
c.s.
1
da
o
(q )
( g )
(0.0764 ) 3
(1.68)
w g = R
bl c.s. o dia 86400seg
lb aire
pie
c.s.
o
aire
3
g
w gd =
26
0.0764 R s q o gd
86400
(1.69)
Captulo 1
w gd
lb gd
3
pie c.s.
pie
c.s.
bl
c.s.
lb
1
da
gd
(q ) o
(0.0764 ) 3 aire ( gd )
(1.70)
= Rs
bl c.s. o dia 86400seg
lb aire
pie aire c.s.
o
3
gd
w gf =
w gf
0.0764 (R R s ) q o gf
86400
(1.71)
lb gl
3
pie c.s.
pie
c.s.
bl
c.s.
lb
1
da
gl
(q ) o
(0.0764 ) 3 aire ( gf )
= (R R s )
bl c.s. o dia 86400seg
lb aire
pie aire c.s.
(1.72)
3
gl
R g R s gd
R Rs
(1.73)
nR =
p V
p V
= cte.
=
Z (T + 460) z (T + 460) c.s.
(1.74)
y:
Bg =
p Z (T + 460)
V
= c.s .
Vc.s. p Z c.s. (Tc.s. + 460)
Bg =
(14.7 ) Z (T + 460)
(520) p (1.0)
B g = 0.02827
Z (T + 460)
p
(1.75)
(1.76)
(1.77)
27
Captulo 1
(1.78)
m
M a gf
(1.79)
g =
28.96 p gf
p gf
M
=
= 2.7
V
Z 10.72 (T + 460)
Z (T + 460)
(1.80)
g =
Donde:
(1.81)
c.s.
Bg
= gf (0.0764)
c.s.
g =
0.0764 gf
Bg
(1.82)
(1.83)
g =
2.7044 p gf
Z (T + 460)
(1.84)
Existen diferentes correlaciones para calcular las propiedades pseudocrticas del gas.
28
Captulo 1
Las curvas correspondientes a los gases han sido establecidas utilizando gases de los
separadores y vapores obtenidos en los tanques de almacenamiento.
La ecuacin para obtener la temperatura y presin pseudocrtica para gases
superficiales es:
Tpc = 167 + 316.67 gf
(1.85)
p pc = 702.5 50 gf
(1.86)
(1.87)
p pc = 740 100 gf
(1.88)
T + 460
Tpc
(1.89)
Ppr =
p
p pc
(1.90)
r =
0.27 p pr
(1.91)
Z sup Tpr
A A 5 A 2
2
A
A
A
Z = 1 + A1 + 2 + 33 r + A4 + 5 r2 + 5 6 r + 7 3 r 1 + A8 r2 e ( A8 r )
T
T pr T pr
T pr
T pr
pr
(1.92)
Donde:
A1 = 0.31506
A2 = 1.0467
A3 = 0.5783
A4 = 0.5353
A5 = 0.6123
A6 = 0.6123
A7 = 0.68157
A8 = 0.68446
29
Captulo 1
Figura 1.9. Factores de compresibilidad para gases naturales. (por Standing y Katz)
30
Captulo 1
31
Captulo 1
32
Captulo 1
Si:
(T
pc
p pc Tpc
(1.93)
(1.94)
) )
+ y H 2S 1 y H 2S 3
(1.95)
Donde:
A = y H 2S + y CO2
33
Captulo 1
Los valores de T' pc y p' pc as obtenidos, se usan en vez de Tpc y Ppc para el clculo
del factor de compresibilidad Z.
1.4.6 Viscosidad del gas
g = K 10 4 e
g
X
62.428
(1.96)
Donde:
gf
(1.97)
gf
Y = 2.4 0.2 X
X = 3.5 +
(1.98)
986
+ 0.2896 gf
(T + 460)
(1.99)
La viscosidad del gas natural corregida, se obtiene con las siguientes expresiones
(correlacin de Carr):
g = gc + C N + C CO + C H S
(1.100)
(1.101)
(1.102)
(1.103)
34
Captulo 1
inters particular. La densidad relativa o densidad del agua es una importante propiedad y
esta normalmente disponible de los anlisis de agua de los campos petroleros.
1.5.1 Factor de volumen
Figura 1.12. Tpica grfica del factor de volumen del agua contra presin
B w = Bw + y B w
Donde:
Bw = Factor de volumen del agua pura.
y
(1.104)
(Eichelberg)
(1.105)
35
Captulo 1
Donde:
TX = T 60
(1.106)
y:
T=o F
Bw .
Mtodo 2
Mc Cain sugiri el mtodo de Katz para calcular Bo se puede utilizar para calcular
Bw = (1 + Vwp )(1 + wT )
(1.107)
Mc Cain indica que el mtodo 2 es vlido para aguas del campo petrolero con
concentraciones de salmuera que varan, puesto que los efectos de la concentracin de la
salmuera en Vwp y VwT se estn compensando.
Los factores del volumen del agua de formacin sobre la presin del punto de
burbujeo se calculan la misma manera que para los sistemas de gas-aceite.
36
Captulo 1
(1.108)
1.5.2 Densidad
Se puede utilizar la misma ecuacin que para la densidad del aceite, excepto que los
efectos de gas en solucin normalmente no se toman en cuenta.
w =
w c.s
Bw
w c.s. 62.4
Bw
(1.109)
La densidad del agua a condiciones estndar esta en funcin de los slidos disueltos y
puede ser obtenida con la figura 1.15.
37
Captulo 1
T +1.98210 5 T 2
(1.110)
Figura 1.16. Viscosidad del agua a presin y temperatura de yacimiento (por Van Wingen)
38
Captulo 1
w1 = 52.5 0.006 (P )
(1.111)
w 2 = 76 e (0.00025P )
(1.112)
280 T
( w 2 w1 ) + w1
206
w =
(1.113)
R sw = A + B p + C p 2
Donde:
(1.114)
) (
(1.115)
) (
) (
(1.116)
(1.117)
(Rsw )b = Rsw CS
C s = 1.0 (0.0753 0.000173) S
(1.118)
(1.119)
Donde (Rsw)b = solubilidad de salmuera en agua scf/STBW, Rsw = solubilidad del gas
en agua pura, CS = correccin por salinidad = factor de correccin por salinidad, donde S es
la salinidad del agua porcentaje en peso de NaCl.
39
Captulo 1
considera 1.0. Actualmente en los campos hay agua que muestra solubilidad al gas, y as, se
contrae bajo una cada de presin.
1.6.1 Compresibilidad
Meehan present una correlacin emprica para determinar el coeficiente isotrmico
de compresibilidad del agua, libre de gas y con gas en solucin.
Para agua libre de gas, la relacin que obtuvo fue la siguiente:
(c w )f = (A + BT + CT 2 ) 10-6
Donde:
(1.120)
Y:
A = 3.8546 (0.000134) p
(1.121)
(1.122)
(1.123)
(c w )g = (c w )f (1.0 + 8.9 10 3 R sw )
Donde:
(1.124)
40
(1.125)
Captulo 1
1.7 Ejemplos
Ejemplo 1
De un fluido de hidrocarburos se tiene la siguiente informacin:
API = 35
g = 0.75 @ c.s.
Bo = 1.23 @ 1500 psia y 250 F
Rs = 277 pie3/bl @ 1500 psia y 250 F
Wgd=1.13 lbm/seg
Wgf=3.85 lbm/seg
Obtener la densidad del aceite o (lbm/pie3) y del gas libre gf (lbm/pie3) @1500 psia y
250 F
Solucin:
Para calcular la densidad del aceite saturado utilizamos la ecuacin 1.49:
o =
62.4 o + 0.01362 R s gd
Bo
o =
141.5
141.5
=
= 0.85
o
131.5+ API 131.5 + 35
Con la correlacin de Katz (ecuacin 1.65) obtenemos la densidad relativa del gas
disuelto.
gd = 0.25 + 0.02(35) + 277 10 6 (0.6874 3.5864(35))
gd = 0.915
Sustituyendo valores en o
o =
lb
62.4(0.85) + 0.01362(277 )(0.92)
= 45.94 m3
1.23
pie
o = 45.94
lb m
@ 1500 psia y 250 oF
3
pie
Ahora para obtener la densidad del gas libre con la ecuacin 1.84, calculamos:
41
Captulo 1
qo =
w gd
0.0764 R s gd
1.13
bl
= 0.058
0.0764(277 )(0.915)
seg
wg
0.0764 q o g
4.98
pie 3
= 1498.5
0.0764 (0.058)(0.75)
bl
gf =
Ppc =
p
1500
=
= 2.249
p pc 666.93
42
gf =
(2.7044)(0.7114)(1500) = 4.458
(0.91)(250 + 460)
lb m
pie 3
Captulo 1
gf = 4.458
lb m
@ 1500 psia y 250 oF
pie 3
Ejemplo 2
Con los siguientes datos:
Gas hmedo
g =0.85
g = gp (1 etapa de separacin)
% NaCl en el agua = 0.18
y CO 2 = 0.0287
43
Captulo 1
30 1.0937
- 25.724
220 + 460
750 e
pb =
= 3,626 psia
(0.0362 )(0.8522 )
pb
750
=
0.85
0.816
220 0.172
0.989 = 22.165
30
2
p b = 10 (1.7669+1.7447log (22.165 )0.30218(log 22.165 ) ) = 3,694 psia
1.553
= 0.9954
(p b )c = (3694)(0.9954)(0.7947 ) = 2,922
psia
o =
141.5
= 0.88
30 + 131.5
750
379.3
= 0.68
yg =
750 + (350 )(0.88)
336.41
379.3
pb =
44
psia
Captulo 1
10
bl
18
1.0937
30
25.724
( 220+ 460 )
= 675
pies 3
bl
1.2255
30 0.989
19.68
R s = 0.85
0.172
220
= 648
pies 3
bl
(0.6849 )(0.88)
pies 3
= 752
R s = 132755
bl
(1 0.6849 )(336.41)
Al comparar cada valor obtenido de Rs y el de laboratorio (Rs,lab), se observa que la
correlacin que se ajusta de manera ms aproximada es la de Lasater.
Para obtener el factor de volumen a la pb mediante las correlaciones vistas utilizamos:
0.88
0.5
+ (1.25)(220 ) = 1013.718
45
Captulo 1
30
B o = 1 + 4.677 10 4 (750 ) + 1.75110 5 (220 60 )
0.85
30
1.81110 8 750(220 60)
= 1.3730
0.85
0.85
= 750
0.88
0.526
+ (0.968)(220 ) = 951.096
R s,lab
R s
(1.127)
C Bo =
Bo,lab
Bo
(1.128)
Donde pb , R s y Bo son los valores obtenidos con las correlaciones vistas para
determinar el factor de correccin.
Utilizando las ecuaciones 1.126, 1.127 y 1.128 para los valores de p b , R s y Bo de
cada correlacin, se obtuvieron los resultados mostrados en la tabla 1.6.
Con base en los resultados obtenidos, se aplicarn las correlaciones de Lasater para
calcular pb y Rs, y la de Standing para calcular Bo.
46
Captulo 1
Tabla 1.6. Relacin de resultados obtenidos de las correlaciones con sus factores de
correccin y con los datos del laboratorio
Correlacin
Rs
CRs
Standing
Vzquez
Oistein
Lasater
Laboratorio
736
675
648
752
750
1.0192
1.1114
1.1568
0.9979
-
C Bo
Bo
1.4723 0.9882
1.3730 1.0598
1.4324 1.0158
1.4550
-
Pb
C pb
3,352
3,634
2,922
3,299
3,300
0.9844
0.9080
1.1294
1.0004
-
La expresin para la densidad relativa del gas disuelto, tambin se debe ajustar para el
valor proporcionado como dato para lo cual se calcula primero gd mediante la ecuacin
1.65:
gd = 0.25 + (0.02 )(30 ) + 750 10 6 (0.6874 (3.5864 )(30 )) = 0.7698
g
0.85
=
= 1.1042
gd 0.7698
(1.129)
(1.130)
gd =0.7698
2. Mediante la correlacin de Lasater se resuelve la ecuacin 1.43.
Pf =
y g = 0.5288
47
Captulo 1
Rs = 457.4468
gf =
48
Captulo 1
A = y H 2S + y CO2
Tpc = 419.05 29.04 = 390.01 o R
Ppc =
662.20 390.01
= 690.24 psia
419.05 + 0.2327 (1 0.2327 ) 29.04
200 + 460
= 1.6923
390.01
Ppr =
2000
= 3.2853
609.01
Y por ltimo obtenemos Z de forma iterativa con las ecuaciones 1.91 y 1.92:
Z = 0.87527
Con las ecuaciones 1.77 y 1.83 obtenemos Bg, y g :
Bg =
g =
0.0764 0.796
lb
= 7.4529 m3
0.00816
pie
La viscosidad del gas se obtiene con las ecuaciones 1.99, 1.98, 1.97 y 1.96:
X = 3.5 +
986
+ 0.2896 0.7960 = 5.2245
(200 + 460)
g = 127.9289 10 e
4
1.3551
5.2245 7.4529
62
.
428
= 0.0172 cp
Para corregir la g por presencia de H2S y CO2 utilizamos las ecuaciones 1.103,
1.102:
49
Captulo 1
0.8762
0.5
Para obtener la viscosidad del aceite muerto utilizamos las ecuaciones 1.57, 1.56,
1.55 y 1.54:
Z = 3.0324 0.02023 30 = 2.4255
Y = 10 2.4255 = 266.379
om = 10 0.5616 1 = 2.6439 cp
Calculamos los factores a y b con las ecuaciones 1.53 y 1.52:
b = 5.44(457 + 150) -0.338 = 0.6234
a = 10.715(457 + 100) -0.515 = 0.4128
50
Captulo 1
El factor de volumen del agua saturada, se obtiene con las ecuaciones 1.106 y 1.105:
TX = 200 60 = 140
Bw = 1.0 + 1.2 10 4 140 + 1.0 10 6 140 2 3.33 10 6 2000 = 1.0297
w =
62.43
lb
= 60.627 m3
1.0297
pie
) = 0.3128
w 2 = 76 e (0.000252000 ) = 46.0963
dina
280 200
(46.0963 40.5) + 40.5 = 42.6733
cm
206
w =
La solubilidad del gas en el agua, se determina con las ecuaciones 1.115, 1.116, 1.117
y 1.114.
51
Captulo 1
) (
La solubilidad del gas es corregida por efecto de salinidad se obtiene con las
ecuaciones. 1.119 y 1.118:
Cs = 1.0-[0.0753-0.000173]S = 1.0-[0.0753-0.000173](0.18) = 0.9865
52
pie 3 /bl
Captulo 2
CAPTULO 2:
FUNDAMENTOS DE FLUJO MULTIFSICO
2.1 Introduccin
En las tuberas, el flujo de gas y lquido ocurre frecuentemente y la precisin del
clculo de la cada de presin es muy importante en la industria del petrleo. Las mezclas
de gas y lquido son transportadas a grandes distancias lo que ocasiona cadas de presin
que influyen en el diseo del sistema.
Las cadas de presin en el flujo multifsico son diferentes al de una sola fase, ya que
en la mayora de los casos existe una interfase, el gas se desliza dejando atrs el lquido lo
que ocasiona superficies de diferentes tipos de rigidez, dependiendo del patrn de flujo.
Cada fase fluye a travs de un rea ms pequea, provocando grandes cadas de presin
comparado con el flujo en una sola fase.
Prdidas
=
Prdidas
Prdidas
+
+
por
por
por
Elevacin
Friccin
Aceleracin
(2.1)
dp
dp
dp
dp
=
+
+
dL total dL elev. dL friccin dL acel.
(2.2)
Captulo 2
E1 + Wf + Ws = E 2
(2.3)
Donde:
Wf = Prdidas de energa por friccin.
Ws = Prdidas de energa por trabajo externo.
E1 = Energa por unidad de masa, en la posicin 1, (lbf pie / lb m ).
E 2 = Energa por unidad de masa, en la posicin 2, (lbf pie / lb m )
54
(2.4)
Captulo 2
Donde:
V = Volumen especfico, pie3 / lbm
Energa potencial ( Ep )
lb pie
pie 1
= g
E p = f
2
lb
seg
gc
m
lb f seg 2
lb m pie
g
h (pie ) =
gc
(2.5)
Donde:
gc = factor de conversin en la segunda ley de Newton = 32.174 (lbm-pie/lbf-seg2)
Energa cintica ( Ec )
E c = f
=
2
lb m 2 seg g c lb m pie 2 g c
(2.6)
v 2
g
+ Wf + WS = 0
h +
2 gc
gc
(2.7)
(2.8)
Donde:
1
Ve = Volumen especfico medio del fluido Ve =
(2.9)
(2.10)
55
Captulo 2
p p p
p
=
+
+
L T L e L ac L f
(2.11)
Ecuacin de Darcy
=
L f 2 g c d
(2.12)
Ecuacin de Fanning
=
L f 2 g c rh
(2.13)
Donde:
d2
4
rh =
d
=d
4
(2.14)
Y:
rh = radio hidrulico =
=
gc d
L f
56
Factor de friccin
(2.15)
Captulo 2
(2.16)
dv
(2.17)
pi
i =1
(2.18)
pi
Li
Ai
i =1
n
Donde:
n
pi = P
i =1
Ps
(2.19)
(pg)
Estriada
0.00006
Produccin o
perforacin
0.0006
Escurrimiento
0.0007
Galvanizada
0.006
57
Captulo 2
Flujo laminar
N Re < 2300
Flujo turbulento
N Re > 3100
Para flujo laminar de una sola fase, el factor de friccin depende exclusivamente del
nmero de Reynolds y est dado por:
64
(2.20)
f=
N Re
Para flujo turbulento (NRe >3100), el factor de friccin est dado por la ecuacin de
Colebrook y White
2.514
+
f = 2 log
f N Re
3.71 d
(2.21)
Para NRe < 2300 (flujo laminar) f depende exclusivamente del nmero de
Reynolds.
A partir de NRe = 3100 se inicia la zona de transicin. Dentro de sta, f
depende tanto de NRe como de / d (rugosidad relativa).
La zona francamente turbulenta se inicia a diferentes valores de NRe,
dependiendo del valor de / d. En esta zona f es independiente de NRe y
varia nicamente con la rugosidad relativa, El valor de f puede obtenerse,
para flujo turbulento, con:
f = 2 log
3.71 d
d)
N Re 2300
1.3521
f=
2300
2.514
2.3026 log 3.71 d + 3100 f
58
(2.22)
+ 0.032
(2.23)
Captulo 2
(2.24)
(2.25)
v2 L
2 g c d
(2.26)
Como:
59
Captulo 2
lb m
3
pie
= 62.428 L
(2.27)
Y:
v=
4 q pies
d2 seg
(2.28)
4 2 q2
2 d4
(2.29)
Entonces:
v2 =
f L q2 L
d5
(2.30)
1 pies 5
125 pg 5
pies 3
pies 3 1 da
bl
= q L
q L
5.6142
da
seg
bl 86400 seg
q L2 =
q 2L (5.6142) 2
86400 2
pg 2
lb
lb
pf f 2 = p f f2 12 2 2
pie
pg
pie
(2.31)
(2.32)
(2.33)
(2.34)
f L q 2L L
d5
(2.35)
Donde:
1.1476 10 -5 = factor de conversin para poder introducir los valores de q, L, d y
pf en unidades prcticas.
60
Captulo 2
Y por ltimo al sustituir las ecuaciones 2.35 y 2.25 en la ecuacin 2.24 obtenemos:
p L = 0.433 L h + 1.1476 10
f L q 2L L
5
d
(2.36)
Que es la ecuacin que determina el gradiente o cada de presin total para el flujo de
lquidos por tuberas en unidades de campo.
2.2.3.2 Nmero de Reynolds para el lquido
N ReL = 92.2
qL L
d L
(2.37)
Donde:
92.2 = Factor de conversin para introducir q, y d en unidades prcticas.
2.2.3.3 Eficiencia de flujo
p L = 0.433 L h + 1.1476 10
2
5
E
d
(2.38)
Para el diseo de tuberas se puede utilizar la ecuacin 2.38, para obtener el dimetro
para un gasto y cada de presin dados. El procedimiento de solucin es iterativo, ya que el
nmero de Reynolds, y por lo tanto, el factor de friccin es funcin del dimetro.
61
Captulo 2
1.1476 10 5 f L q 2 L
d= 2
E (p L 0.433 L h )
0.2
(2.39)
p L 0.433 L h
q = E d 5
5
1
.
1476
10
f
0.5
(2.40)
Pf = f
L v2
(2.41)
2 g c d
g = 0.0764
(2.42)
Bg
p o Z(T + 460)
B g =
p
To + 460
(2.43)
q 4 q Bg
=
A
d 2
(2.44)
v=
v = 1.6211
q 2 Bg
d4
(2.45)
(0.0764 ) (1.6211q
=f
g
2 B g d
B g L
2
(2.46)
Simplificando:
g q 2 B g L
f
Pf =
519.98
d 5
Sustituyendo la ecuacin 2.43 en 2.47:
62
(2.47)
Captulo 2
2
f p o Z g q (T + 460) L
Pf =
d5 p
519.98 T o +460
Los trminos de la ecuacin 2.48 excepto p , estn en unidades consistentes:
(
)
(lb/pg abs)
(F )
p f lb/pie 2
p
T
q pies 3 /seg
L (pies )
d (pies )
(2.48)
p=
(p1 p 2 )
2
(2.49)
y:
p = p 1 - p 2
(2.50)
(p
2
1
p2
2
f
p o Z q g (T + 460) L
=
d5
33513.5316 To + 460
(2.51)
Ecuacin que permite evaluar las prdidas de presin por friccin en gasoductos
horizontales con unidades consistentes de campo: q = pies3/da, L = pies, d = pg y pf = psi.
De la ecuacin 2.5, las prdidas de presin debido al cambio de elevacin se expresan
de la siguiente forma:
p e = h
(2.52)
2.7044
p e =
144 Z (T + 460 )
(2.53)
Y como
p e = p1 p2
(2.54)
63
Captulo 2
p1 p 2 = 0.01878
p g h
(p
(2.55)
Z (T + 460 )
(p1 + p 2 )
2
se obtiene finalmente:
2
1
p2
p g h
) = 0.03756 Z (T + 460)
(2.56)
Sumando las ecuaciones 2.51 y 2.56 se obtiene la siguiente ecuacin para calcular p1
o p2 en unidades prcticas y bajo las siguientes condiciones supuestas:
(p
2
1
p2
2
2
p g h
Z q g g (T + 460 ) L
f
p o
+ 0.03756
2
Z (T + 460 )
d5
33513.5316 To + 460
(2.57)
Para determinar el gasto:
64
Captulo 2
2
0.03756 p g h
K4 2
2
q g = 0.5 p1 p 2
Z (T + 460 )
L
0.5
d 2.5
(2.58)
L
d =
K
4
1
2
qg
2
0.03756
h
2
g
2
p1 p 2
Z (T + 460 )
0.2
(2.59)
Donde:
T + 460
K 4 = 33513.5316 o
po
0.5
f Z (T + 460 )
g
0.5
(2.60)
N Re g =
g =
v=
Bg =
d v g
(2.61)
0.0764 g
(2.62)
Bg
4 q g Bg
(2.63)
d2
0.02825 Z (T + 460 )
p
(2.64)
N Re g
g
d
(
2.07044 ) (4) (0.02825) q g
=
(2.65)
65
Captulo 2
(2.66)
De igual forma que en las tuberas con flujo de lquidos, la eficiencia es un factor de
ajuste para compensar los efectos de corrosin, erosin, rugosidad e incrustaciones, que no
se consideran en la deduccin de las ecuaciones de flujo, por lo tanto los resultados
obtenidos se tienen que corregir y as obtener un gasto ms real.
Los valores ms comunes de eficiencia E para flujo de gas son:
Tabla 2.2. Valores comunes de Eficiencia de Flujo.
Lnea
Contenido de
E
lquidos
Gas seco
0.1
0.92
Gas hmedo
7.2
0.77
Gas y condensado
800
0.60
66
Captulo 2
VL
VP
(2.67)
Donde:
HL = Colgamiento de lquido.
VL = Volumen de lquido en la seccin de tubera.
Vp = Volumen de la seccin de tubera.
Resbalamiento.- Se usa para describir el fenmeno natural del flujo a mayor
velocidad de una de las fases. El resbalamiento entre fases en el flujo multifsico en
tuberas es inevitable a cualquier ngulo de inclinacin. Las causas de este fenmeno son
diversas, a continuacin mencionaremos algunas de las causas:
67
Captulo 2
Para cada ngulo ascendente o ascendente, el espacio vaco se traz como una
funcin de velocidad superficial de gas para la velocidad superficial del lquido.
Para gas a muy altos gastos, las curvas casi se vuelven asintticas en 0% de
colgamiento.
Las formas generales de estas graficas incitaron la seleccin de una ecuacin de
regresin no lineal.
La correlacin es la siguiente:
H L = eA
(2.68)
Donde:
N Cgv5
2
2
A = C1 + C 2sen + C 3 sen + C 4 N L C6
N Lv
2
L
Y:
C = Constante Emprica.
68
= 0.15726 L
3
L
0.25
(2.70)
0.25
N Lv
= 1.938 v sL L
L
0.25
N gv
= 1.938 v sg L
L
(2.69)
(2.71)
(2.72)
Captulo 2
En la tabla 2.3 aparecen los coeficientes de la ecuacin 2.69 para diferentes patrones
y direcciones de flujo, usando programas de regresin no lineal.
Las ecuaciones de 2.68 y 2.69 muestran que el colgamiento debe incrementar al
mismo tiempo que el ngulo ascendente se incrementa. Esto se demuestra grficamente
comparando los valores de colgamiento del lquido, obtenidos de la grafica en la figuras 2.5
y 2.6, donde la fraccin de vaco es graficada para algn aceite en 3 velocidades
superficiales de lquido similares para tubera horizontal y ascendente con una inclinacin
de 30.
La fraccin de vaco o gas, se define como:
g = (1 H L )
Tabla 2.3. Coeficientes de la ecuacin 2.69 para diferentes patrones y direcciones de flujo.
Direccin
Patrn de
C2
C3
C4
C5
C6
C1
del flujo
flujo
Horizontal
o
Todos
-0.38011 0.12988 -0.11979 2.34323 0.47569 0.28866
Ascendente
Estratificado -1.33028
4.80814
4.17156
56.2623
0.07995
0.50489
0.78981
0.55163
15.5192
0.37177
0.39395
Descendente
Otros
-0.51664
Figura 2.5. Fraccin de vaco vs velocidad superficial del gas a una velocidad superficial
del lquido fija, para flujo horizontal.
69
Captulo 2
Figura 2.6. Fraccin de vaco vs velocidad superficial del gas a una velocidad superficial
del lquido fija, para flujo ascendente a 30o.
2.3.1 Colgamiento sin resbalamiento ()
qL
=
qL + qg
1+
1
q o (R R s )Bg
(2.73)
5.615(q o Bo + q w B w )
qL 6.4979 10 5 (q o B o + q w B w )
=
Ap
A p / 144
v sg =
q g
Ap
1.1574 10 5 q o (R R s )B g
A p /144
(2.74)
(2.75)
o:
v sg =
Entonces:
70
qg
Ap
(2.76)
Captulo 2
vm =
qL + qg
Ap
= v sL + v sg
(2.77)
v sL
vm
(2.78)
2
d ci d 2te
4
(2.79)
v sg =
0.01191(q o Bo + q w B w )
d ci2 d 2te
0.002122 q o (R R s ) Bg
(d
2
ci
d 2te
(2.80)
(2.81)
vg =
qg
Ag
v
q L
q L
=
= sL
AL Ap HL HL
=
qg
A p (1 H L )
v sg
(1 H L )
(2.82)
(2.83)
m = L H L + g (1 H L )
(2.84)
71
Captulo 2
ns = L + g (1 )
(2.85)
ns =
M
Vm
(2.86)
Donde:
Los valores de M y Vm se obtienen con las ecuaciones siguientes:
M = Mo + Mg + Mw
lb
(lb o /pie 3o )
62.48 w3
M o = o
3
pie w
(lb w /pie w )
(2.87)
pie 3
5.615 o
bl o
(2.88)
M o = 350.5 o
(2.89)
3
(lb g /pie 3g )
lb a pie gp a c.s.
0.0764 3 R
Mg = g
(lb /pie 3 )
pie
bl
o
a
a
a
(2.90)
M g = 0.0764 g R
(lb /pie 3w )
lb
62.428 w3
M w = w w
3
(lb w /pie w )
pie w
(2.91)
pie 3
5.615 w
bl w
M w = 350.5 w WOR
(2.92)
(2.93)
(2.94)
pies3o + g + w @ c.s.
Vm
bl
@
c.s.
o
Y:
72
(2.95)
Captulo 2
Vmo
Vmg
pies 3o
= 5.615
bl o
pies 3o a c.esc.
B o
pies 3 a c.s.
o
bl
Vmw = WOR w
bl o
pies 3w
5.615
bl w
pies 3g @ c.esc.
pie 3 @ c.s.
g
Bw
pies 3w @ c.esc.
pie 3 @ c.s.
w
(2.96)
(2.97)
(2.98)
Por lo tanto:
(2.99)
ns =
(2.100)
lb m de lquido + gas
segundo
(2.101)
qo M
86400
wm = wo + ww + wg
wo =
ww =
o q o B o
15388
w q w Bw
15388
(2.102)
(2.103)
(2.104
(2.105)
73
Captulo 2
wg =
g q o (R R s ) Bg
86400
(2.106)
Dependiendo del mtodo que se aplique, se usan las siguientes ecuaciones para
obtener la viscosidad de la mezcla de dos fases:
ns = L + g (1 )
(2.107)
m = HL L (1gH L )
(2.108)
(2.109)
fo =
q o Bo
q w B w + q o Bo
(2.110)
fo =
Bo
Bo + WOR B w
(2.111)
Donde:
f w = 1 fo
(2.112)
(2.113)
L = o fo + w f w
74
(2.114)
Captulo 2
Tericos
Ventaja
Desventaja
Los patrones de flujo en tuberas horizontales descritos por Beggs son los siguientes:
Flujo Segregado Estratificado: Este patrn de flujo se presenta relativamente a bajos
gastos de gas y lquido, para el cual las dos fases son separadas por efecto de la gravedad,
donde el lquido fluye en el fondo de la tubera y el gas en la parte superior.
75
Captulo 2
Flujo Segregado Ondulado: ste rgimen de flujo se presenta a gastos mas altos que
en el estratificado, con presencia de ondas estables en la interfase.
Flujo Segregado Anular: El flujo anular se presenta a muy altos gastos de flujo de
gas. La fase gaseosa fluye como un ncleo a alta velocidad, el cual puede llevar gotas de
lquido atrapadas. La fase lquida fluye como una pelcula delgada pegada a la pared interna
de la tubera, generalmente, esta pelcula es ms gruesa en el fondo que en la pared superior
de la tubera, dependiendo de la magnitud relativa de los gastos de flujo de gas y lquido.
Flujo Intermitente : Este flujo es caracterizado por el flujo alternado de lquido y
gas, fluyendo sucesivamente tapones o baches de lquido ocupando completamente el rea
transversal de la tubera, separados por bolsas o burbujas de gas, el cual contiene una capa
estratificada de lquido que a su vez se desplaza en el fondo de la tubera.
Este tipo de flujo es consecuencia de la inestabilidad hidrodinmica de una interfase
estratificada lquido-gas bajo ciertas condiciones, donde el mecanismo de flujo consiste de
una pelcula de lquido que fluye a lo largo del fondo de la tubera a una velocidad menor
que la del bache, y que crece al arrastrar lquido de la pelcula en su parte delantera, pero a
su vez pierde lquido en la parte trasera en una misma proporcin
El lquido en el cuerpo del bache puede ser aireado por pequeas burbujas que son
concentradas hacia el frente del bache y en la parte superior de la tubera. El flujo
intermitente es dividido en tapn y bache, y el comportamiento del flujo bache y tapn son
los mismos con respecto al mecanismo de flujo, y por lo tanto, no existe distincin entre
ellos generalmente.
Patrones de Flujo en Tuberas Horizontales/Inclinadas
Estratificado
Flujo Segregado
Ondulado
Anular
Tapn
Flujo Intermitente
Bache
Burbuja
Flujo distribuido
Niebla
76
Captulo 2
Flujo Intermitente Tapn: Es considerado como el caso lmite del flujo bache,
cuando el bache de lquido est libre de burbujas, lo cual ocurre a gastos de gas
relativamente bajos cuando el flujo es menos turbulento.
Flujo Intermitente Bache: A altos gastos de gas, donde el flujo en el frente del bache
est en forma de un remolino (debido al recogimiento del movimiento lento de la pelcula)
se denomina flujo bache.
Flujo Burbuja o Burbujas Dispersas: En este tipo de patrn de flujo la tubera se
encuentra casi llena de lquido y la fase de gas libre es pequea. El gas est presente en
pequeas burbujas distribuidas aleatoriamente, al igual que sus dimetros. Las burbujas se
mueven a diferentes velocidades dependiendo de sus respectivos dimetros, el lquido se
mueve a una velocidad bastante uniforme y, a excepcin de la densidad, la fase de gas tiene
un efecto mnimo en el gradiente de presin. Existe sin embargo, condiciones donde hay
pequeas burbujas discretas a bajos gastos, que son a veces designadas como flujo burbuja.
La diferencia entre flujo burbuja y burbujas dispersas no siempre es claramente visible. El
flujo de burbujas dispersas se observa sobre un rango completo de inclinacin de tubera,
mientras que el patrn de flujo burbujeante es observado solamente en vertical y tuberas de
dimetro relativamente grandes.
Niebla o Neblina: En este patrn la fase continua es el gas, el cual arrastra y
transporta al lquido. El lquido deja una pelcula en la pared de la tubera, pero sus efectos
son secundarios, el gas es el factor predominante.
2.4.1.1 Correlacin de Taitel y Dukler para predecir la transicin de los patrones de flujo
Taitel y Dukler en 1976 obtuvieron un modelo terico para predecir con exactitud la
transicin entre los patrones de flujo basndose en modelos de mecanismos fsicos. Los
patrones considerados son el intermitente (bache y tapn), estratificado liso, estratificado
ondulado, distribuido burbuja y anular niebla. La correlacin predice los lmites de
transicin del flujo y el efecto que stos tienen en el dimetro de la tubera, las propiedades
de los fluidos y el ngulo de inclinacin.
Los clculos realizados por Taitel y Dukler arrojaron los siguientes grupos
adimensionales:
dp
dx L
X=
dp
dx g
(2.115)
77
Captulo 2
T=
( L
2
dp
dx L
g ) g cos
(2.116)
2
g
F = vg
(d g cos )( L g )
(2.117)
2
g vg 2 vL
K=
( L g ) g L cos
Y=
g ) g sen
dp
dx g
(2.118)
(2.119)
78
Captulo 2
Figura 2.8. Mapa de patrones de flujo generalizado para flujo multifsico horizontal
(por Taitel y Dukler)
Donde:
Curva
AyB
Coordenadas
F vs X
K vs X
T vs X
79
Captulo 2
Tapn: Aqu la fase de gas es ms pronunciada, aunque la fase lquida sigue siendo
continua las burbujas de gas se unen y forman burbujas estables de aproximadamente el
mismo tamao y forma que la tubera (que estn rodeadas por una pelcula de lquido), y
son separadas por tramos de lquido. La velocidad de la burbuja es mayor que la del lquido
y puede ser predicho en relacin a la velocidad del bache de lquido. La velocidad del
lquido no es constante mientras el tramo o bache de lquido se mueva siempre hacia arriba
(en la direccin del flujo), el lquido de la pelcula que rodea a la burbuja podra moverse
hacia arriba, aunque posiblemente a una menor velocidad, o incluso podra moverse hacia
abajo. Esta variacin de las velocidades del lquido, podra resultar no slo en una variacin
de las prdidas por friccin en la pared, sino tambin en el colgamiento de lquido, que est
influenciado por la densidad de la mezcla que fluye por la tubera. A mayores velocidades
de flujo, el lquido puede incluso ser arrastrado dentro de las burbujas de gas.
Transicin: El cambio de una fase continua de lquido a una continua de gas ocurre
en este patrn de flujo. El bache de lquido entre las burbujas virtualmente desaparece, y la
fase gaseosa arrastra una cantidad significativa de lquido, aunque los efectos del lquido
son significativos, el gas es el que predomina.
Neblina: En este patrn la fase continua es el gas, el cual arrastra y transporta al
lquido. El lquido deja una pelcula en la pared de la tubera, pero sus efectos son
secundarios, el gas es el factor predominante.
MONOFSICO
BURBUJA
TAPN
TR ANSICIN
2.5 Ejemplos
Dado en un pozo vertical:
qo = 1000 bl/da
R = 1000 pies3/bl
Rs =450 pie3/bl
Bo = 1.27
80
NEBLINA
Captulo 2
d = 1.995 pg
Z = 0.85
gf = 0.7
o
API = 30
p = 2,500 psia
T = 200 oF
a)
b)
c)
12 v m
gd
(2.120)
(0.845)L 0.5351
N Fr
0.0173
(2.121)
Calcular:
dp
dz e
Si:
g
dp
m sen
=
dz e g c
d)
(2.122)
Solucin:
a)
v sg =
6.4979 10 5 (1000)(1.27 )
pies
= 3.8016
2
1.995
seg
(4)(144)
81
Captulo 2
=
b)
v sL 3.8026
=
= 0.6717
v m 5.6608
o =
lb
62.4(0.8762) + 0.01362(450 )(0.8019)
= 46.9210 m3
1.27
pie
g =
lb
2.7044(2500)(0.7 )
= 8.4352 m3
0.85(200 + 460 )
pie
N Fr =
HL =
= 0.6621
d)
Las velocidades reales del lquido y del gas las obtenemos con las
ecuaciones 2.82 y 2.83
vL =
vg =
82
v sL 3.8016
pies
=
= 5.7417
H L 0.6621
seg
v sg
1.8592
pies
=
= 5.5022
(1 H L ) (1 0.6621)
seg
Captulo 3
CAPTULO 3:
FLUJO MULTIFSICO EN TUBERAS VERTICALES
3.1 Introduccin
El problema de predecir con precisin las cadas de presin en la tubera de
produccin (T.P.) se ha ido incrementando a muchas soluciones especializadas para
condiciones limitadas. La razn para estas muchas soluciones es que el flujo multifsico es
complejo y se dificulta su anlisis incluso para estas correlaciones estudiadas de
condiciones limitadas.
Al pasar los fluidos provenientes del yacimiento a travs de la tubera de
produccin, se consume la mayor parte de presin disponible para llevarlos del yacimiento
a las bateras de separacin, por lo que es de suma importancia realizar una evaluacin
precisa de la distribucin de la presin a lo largo de dicha tubera. Al hacerlo
conjuntamente con un anlisis integral del sistema de produccin, es posible:
a.
b.
c.
d.
e.
f.
Cuando el flujo es vertical las cadas de presin por aceleracin son muy pequeas,
por lo que el gradiente de presin debido a la misma generalmente se desprecia, quedando
la ecuacin 2.11 de la siguiente forma:
p
p p
=
+
h T h e h f
(3.1)
83
Captulo 3
Orkiszewski (1967)
Beggs y Brill (1973)
Tabla 3.1. Criterios en el desarrollo de las correlaciones.
Grupo I
Grupo II
Grupo III
Resbalamiento
entre las fases
No se considera
Si se considera
Si se considera
Densidad de la
mezcla ( m)
Se calcula en
funcin de las
propiedades de los
fluidos
Factor de friccin
(f)
Se determina
empricamente
Se calcula
considerando el
efecto del
colgamiento
Se correlaciona con
las propiedades
combinadas del gas
y el lquido
Se calcula
considerando el
efecto del
colgamiento
Se correlaciona con
las propiedades del
fluido en la fase
continua
Patrones de flujo
No se distinguen
No se distinguen
S se distinguen
Mtodos o
Correlaciones
Poettmann y
Carpenter, Fancher y
Brown y Baxendell
y Thomas
Hagedorn y Brown
Duns y Ros,
Orkiszewski, Aziz,
Beggs y Brill,
Chierici, Gould y
Tek entre otras
3.2 Correlaciones
3.2.1 Poettman y Carpenter
Poettmann y Carpenter desarrollaron un mtodo semi emprico utilizando dimetros
de tubera de produccin de 2, 2 , 3 pulgadas en el que incorporaron la ecuacin general
de energa y datos de 34 pozos fluyentes de aceite y 15 con bombeo neumtico. El aceite,
gas y agua fueron considerados para estar en una sola fase, y no se trato de determinar la
correlacin de colgamiento. Toda la prdida de energa incluyendo el resbalamiento y
lquidos de emergencia se introdujeron en un factor de perdida de energa la cual se tom
como constante en toda la longitud de la tubera. Tambin se asumi el flujo vertical de
aceite, gas y agua como turbulento, adems de no tomar en cuenta los patrones de flujo.
La correlacin de Poettmann y Carpenter sigue siendo usada extensamente en los
campos bajo las siguientes condiciones:
84
Tuberas de 2, 2 y 3 pg.
Viscosidades menores a 5 cp.
Captulo 3
(3.2)
(3.3)
m =
M
VT
(3.4)
Captulo 3
dv =
1.4737 10 5 w m
d/12
(3.6)
f wm
1
dp
=
+
m
5
10
(
)
7
.
413
10
d/12
dh T 144
m
Donde:
m =
m1 + m 2
2
(3.7)
(3.8)
86
(3.9)
Captulo 3
87
Captulo 3
WOR
M = o (350 )
+ w (350 )
+ (0.0764 ) R g
1 + WOR
1 + WOR
(3.10)
WOR
+ w
1 + WOR
1 + WOR
L = 62.428 o
(3.11)
p=
88
p1 + p 2
+ p atm
2
(3.12)
Captulo 3
T1 + T2
2
(3.13)
6. Obtener Z
7. Calcular la densidad promedio de la fase de gas:
520 p 1
T + 460 14.7 Z
g = g (0.0764)
(3.14)
8. Calcular la viscosidad promedio del aceite con las ecuaciones 1.51 1.57.
9. Calcular la viscosidad promedio del agua con las ecuacin 1.110.
10. Calcular la viscosidad de la mezcla del lquido:
WOR
L = o
+ w
1 + WOR
1 + WOR
(3.15)
WOR
+ w
1 + WOR
1 + WOR
L =o
(3.16)
12. Calcular Rs a p y T .
13. Calcular Bo a p y T .
14. Calcular el rea transversal de la tubera.
15. Calcular el nmero de viscosidad del lquido:
N L
= 0.15726 L
3
L
0.25
(3.17)
v sL =
(5.6142)(144) q L
86,400 A p
1
WOR
B o 1 + WOR + B w 1 + WOR
(3.18)
89
Captulo 3
= 1.938 v sL L
L
0.25
(3.19)
144 q L R R s
v sg =
86,400 A p
p 520 1
(3.20)
= 1.938 v sg L
L
0.25
(3.21)
Nd =
12
L
0.5
(3.22)
Para la Regin I:
El resbalamiento se obtiene de la siguiente ecuacin:
N gv
S = F1 + F2 N Lv + F
1 + N Lv
'
3
(3.23)
F3' = F3
F4
Nd
90
(3.24)
Captulo 3
(3.25)
(3.26)
(N )
)
0.982
+ F6'
gv
(1 + F7 N Lv )2
(3.27)
F6' = 0.029 N d + F6
(3.28)
(3.29)
(3.30)
Por lo tanto:
HL =
1
1 + v sg /v sL
(3.31)
N gv > 75 + 84 N Lv
0.75
(3.32)
1.938 ( L / L )
0.25
(3.33)
91
Captulo 3
HL =
v s v sg v sL + (v s v sg v sL ) + 4 v s v sL
2
0.5
2 vs
(3.34)
1488 L v sL d
12 L
(3.35)
fr = 2 f w
N Lv (N Lv + N gv )
Nd
(3.36)
Donde:
f w = (f 1 )
f2
f3
(3.37)
Y:
f1 se obtiene de la figura 3.7.
f2 se obtiene de la figura 3.8 y para obtenerla se debe determinar la
( )
abscisa f 1 R N d
2/3
Donde:
R=
v sg
(3.38)
v sL
f3 = 1 + f1 (R/50)0.5
(3.39)
(N )
fr = 2 f w N
Donde:
92
gv
Nd
(3.40)
Captulo 3
N =
g
L
(3.41)
1
1.73
+ 0.067 ( /d )
2
[4 log (0.027 /d )]
(3.42)
v sg d 2
(d )2
(3.43)
st = H L + (1 H L ) N
(3.44)
Donde:
N =
g
L
(3.45)
T = st + fr
(3.46)
st + fr
1 ( L v sl + g v sg )(v sg /p )
(3.47)
(3.48)
dp fr L
;
=
dh fr
144
(3.49)
o:
93
Captulo 3
dp T L
=
dh T
144
30. Determinar la longitud para esa cada de presin:
h =
p
dp
dh total
94
(3.50)
(3.51)
Captulo 3
Figura 3.6. F5, F6, y F7 vs Nmero de viscosidad del lquido NL (por Ros).
95
Captulo 3
96
Captulo 3
Comenzar con una presin conocida p1, y asumir el valor para p2. Calcular el
incremento de profundidad.
1. Calcular la presin media entre los dos puntos de la presin, psia.
p=
p1 + p 2
+ 14.7
2
(3.52)
Dependiendo de los requisitos del problema, es decir, si o no la presin del fondo del
agujero que no fluye debe ser determinada de la informacin superficial, o si los
97
Captulo 3
clculos comienzan de la profundidad total y sube por la tubera, la presin inicial debe
ser conocida. Los incrementos o decrementos de la presin deben ser asumidos de entre
los puntos de presin (1) y (2) los cuales sern calculados.
2. Calcular la densidad relativa del aceite o .
o =
141.5
131.5 + API
(3.53)
WOR
m = o (350)
+ w (350)
+ (0.0764) (RGL ) ( g )
1 + WOR
1 + WOR
(3.54)
w = qm
(3.55)
o (62.4 ) +
1
WOR
5.614
L =
+ w (62.4 )
BO
1 + WOR
1 + WOR
(3.56)
7. Obtener Z .
8. Calcular la densidad media de la fase gaseosa.
p 520 1
14.7 T Z
g = g (0.0764 )
(3.57)
9. Calcular la viscosidad media del aceite con la correlacin apropiada (captulo 1).
10. Calcular la viscosidad media del agua.
11. Calcular la viscosidad lquida de la mezcla.
1
WOR
L = o
+ w
1 + WOR
1 + WOR
98
(3.58)
Captulo 3
Esto puede ser solamente una aproximacin, puesto que la viscosidad de dos
lquidos inmiscibles es absolutamente compleja.
12. Asumiendo la tensin superficial constante en cada punto de presin, calcular la
tensin superficial de la mezcla lquida.
1
WOR
+ w
1 + WOR
1 + WOR
L = o
(3.59)
1 4
N L = 0.15726 L
3
L L
(3.60)
(3.61)
16. Obtener Bo a p y T .
17. Asumiendo B W = 1.0 , calcular la velocidad superficial del lquido v sL ,
v sL =
5.61 q L
86400 A p
WOR
B O
+ BW
1 + WOR
1 + WOR
pie
seg
(3.62)
N LV
= 1.938 VsL L
L
(3.63)
99
Captulo 3
q L RGL R S
1 + WOR 14.7 T Z
v sg =
86400 A P
p 520 1
(3.64)
N Gv = 1.938 v sg L
L
1
4
(3.65)
[0.2218 (v
A = 1.071 -
+ v sg )
sL
(3.66)
v sg
(3.67)
v sL + v sg
L
L
(3.68)
=
14.7
N
N
d
gv
24. Obtener de la figura 3.10, el valor de
100
HL
.
(3.69)
Captulo 3
=
N d2.14
(3.70)
H
H L = L ( )
(3.71)
2.2 10 2 w
=
(d ) HL L (g1HL )
( )(
(3.72)
d
Si el valor para no es conocido, un buen valor para usar es 0.00015 pies, el cual
es un valor promedio dado para el acero comercial.
v m1 = v sL1 + v sg1
(3.73)
v m2 = v sL2 + v sg2
(3.74)
( )
( ) [
v 2m = v 2m1 v 2m2
(3.75)
101
Captulo 3
v2
144 p m m
2g
c
h =
f w2
m +
2.9652 1011 d 5 m
(3.76)
36. Comenzando con p2, y la profundidad conocida a la p2, asumir otro punto de presin
y repetir el procedimiento hasta alcanzar la profundidad total o hasta alcanzar la
superficie.
Figura 3.9. Correlacin para el coeficiente C del nmero de viscosidad (por Hagedorn).
Captulo 3
103
Captulo 3
3.2.4 Orkiszewski
Procedimiento de clculo:
104
Captulo 3
v sg
< LB
Burbuja
> L B , N gv < L S
Bache
vm
v sg
vm
Rgimen de Flujo
L M > N gv > L S
Transicin bache-niebla
N gv > L M
Niebla
Donde LB, LS, y LM son los lmites de burbuja-bache, bache-transicin, transicinniebla respectivamente, adimensionales, y Ngv se calcula con la ecuacin 2.72.
Donde:
2
v
L B = 1.071 2.6616 m
dh
(3.77)
0.75
(3.78)
(3.79)
1
p
( L H L + g (1 H L ))
=
L e 144
(3.80)
C1 C 2
2
(3.81)
105
Captulo 3
Y:
4
2
C 2 = C1
v sg
0.8
C1 = 1 +
0.5
vm
0.8
(3.82)
(3.83)
L f 144 2 g d h
(3.84)
124 d h v L L
L
(3.85)
b) Rgimen de bache
1 C3
p
+ L
=
L e 144 v m + v b
(3.86)
C 3 = L (v sL + v b ) + g v sg
(3.87)
Donde:
106
Captulo 3
<10
(3.90)
>10
(3.91)
1.38
d
= 0.709 + 0.0451 h
12
d
= 0.284 + 0.0127 h
12
d
= 0.161 + 0.0274 h
12
1.371
d
log L + 0.232 log v m 0.428 log h
12
0.799
1.415
(3.88)
d
log L 0.162 log v m 0.888 log h
12
(3.89)
d
log ( L + 1) + 0.167 log v m + 0.113 log h
12
(3.90)
d
log ( L + 1) + 0.569 log h
12
1.571
d
dh
log ( L + 1) + 0.63 log h log v m
0.397 + 0.01
12
12
(3.91)
(3.92)
Para vm > 10
107
Captulo 3
vb
vm + vb
m
1
(3.93)
124 d h v m L
L
(3.94)
124 d h v b L
L
(3.95)
N Re b =
12
0.5
(3.96)
2. Si 3000 N Re b 8000 :
13.59 L
v b = 0.5 + 2 +
0.5
(
)
d
/12
L
h
0.5
(3.97)
Donde:
= 0.251 + 8.74 10 N Re L
g dh
12
0.5
(3.98)
3. Si N Re b 8000 :
g dh
v b = 0.350 + 8.74 10 6 N Re L
12
108
0.5
(3.99)
Captulo 3
12
0.5
(3.100)
Donde:
C 2 = 1.36 + C 5 + C 6 C 3 + C 7 C 3
(3.101)
3
(3.103)
C4 =
C3 =
N Re b 5,500
(3.105)
1,000
N Re L 3,000
(3.106)
1,000
2
C8 =
N Re b
(3.107)
(3.108)
10
=
L f 12 (2 ) g d h
2
v sL + v b
+
vm + vb
(3.109)
109
Captulo 3
N Re L =
124 d h v m L
L
(3.110)
(3.111)
Donde:
L m = 84 N Lv
0.75
+ 75
(3.112)
L m N gv
p
=
L e L m L s
L e
+
bache
N gv - L s p
L m L s L e
(3.113)
niebla
=
Lm Ls
L f
L f
+
bache
N gv - L s p
L m L s L f
(3.114)
niebla
q g = A p L m L
g L
0.25
(3.115)
d) Rgimen de niebla.
Para calcular el gradiente de presin correspondiente a esta regin se aplica
el mtodo de Duns y Ros.
La regin de niebla est definida para:
Ngv > Lm
110
(3.116)
Captulo 3
=
vm
L e 144
(3.117)
=
2 g dh
L f
2
(3.118)
124 d h v sg g
g
(3.119)
(3.120)
0.8994 L
(3.121)
g v sg 2 d h
Y si:
N w 0.005
4.6242 L N w
g v sg 2 d h
(3.122)
0.302
(3.123)
Donde:
111
Captulo 3
g v sg L
N w = 0.093
L L
(3.124)
p
h e h f
=
1 Ek
h T
(3.125)
Donde:
Ek =
Wm v sg
4,637 p A p
(3.126)
112
Captulo 3
horizontal. Una vez establecido cada patrn de flujo se procedi a variar el ngulo de
inclinacin, as que se pudo observar como el ngulo de inclinacin afectaba el colgamiento
y el gradiente de presin. El colgamiento y el gradiente de presin fueron medidos en
ngulos que variaban de 5, 10, 15, 20, 35, 55, 75 y 90 grados, y se encontr que el
colgamiento llegaba a su valor mximo en +50 grados y a su valor mnimo en -50 grados.
El mapa de patrones de flujo original que obtuvieron Beggs y Brill fue ligeramente
modificado para poder incluir la zona de transicin entre el patrn de flujo segregado y el
intermitente. El mapa de patrones de flujo modificado fue sobrepuesto al original y se
muestra en la figura 3.13.
La ecuacin para determinar el gradiente de presin es:
g
f T M m v m 12
m sen +
2 gc d
p g c
=
h
m v m v sg
144 1
g c p (144)
(3.127)
113
Captulo 3
m v m v sg
p (144 )1
g c (p + 14.7 )(144 )
h =
g m
f M m v m 12
sen + T
gc
2 gc d
(3.128)
Procedimiento de clculo:
Comenzando con una presin conocida p1, seleccionar el valor para p2 y suponer un
valor h entre los dos puntos de modo que la T , que est en funcin de la profundidad,
pueda ser evaluada. Calcular el incremento en la profundidad, que debe coincidir con el
valor supuesto.
1. Calcular la presin (psia) y profundidad promedio entre los dos puntos de presiones
dadas:
p=
p1 + p 2
+ p atm
2
h=h+
h
2
(3.129)
(3.130)
para p1
2. Determinar la temperatura promedio T , a la profundidad promedio. Este valor se
debe conocer de la relacin temperatura vs profundidad (gradiente de temperatura).
3. Del anlisis PVT o las correlaciones apropiadas, calcular a la p y T :
R s , Bo , B w , o , w , g , o , w , Z
(3.131)
141.5
131.5 + o API
(3.132)
+ w
1 + WOR
1 + WOR
L = o
o =
114
(350
+ 0.0764 R s g )
5.6146 B o
(3.133)
(3.134)
Captulo 3
w =
g =
350 w
5.6146 B w
(3.135)
0.0764 g p 520
(3.136)
3.27 10 7 Z q o (R - R s ) (T + 460 )
p
q L = 6.49 10 5 (q o B o + q w B w )
(3.137)
(3.138)
Donde:
qL y qg = pies3/seg
7. Calcular las velocidades superficiales del gas, lquido y la mezcla in situ.
v sL =
v sg =
144 q L
Ap
144 q g
Ap
v m = v sL + v sg
(3.139)
(3.140)
(3.141)
(3.142)
G g = g v sg
(3.143)
GT = GL + Gg
(3.144)
N FR =
vm
g d/12
(3.145)
115
Captulo 3
WOR
L = o
+ w
1 + WOR
1 + WOR
m = L + g (1 )
WOR
+ w
1 + WOR
1 + WOR
L =o
(3.146)
(3.147)
(3.148)
N Lv
M m d/12
m 6.72 10 4
= 1.938 v sL L
L
(3.149)
0.25
(3.150)
12. Para determinar el patrn de flujo que existe en el flujo horizontal, calcular los
parmetros correlacionados, L1, L2, L3, y L4.
L1 = 316 0.302
(3.151)
L 2 = 0.0009252 2.4684
(3.152)
L 3 = 0.10 -1.4516
(3.153)
L 4 = 0.5 -6.738
(3.154)
Transicin
Intermitente
Distribuido
116
Captulo 3
a b
(N FR )c
(3.155)
Segregado
0.98
0.4846
0.0868
Intermitente
0.845
0.5351
0.0173
Distribuido
1.065
0.5824
0.0609
C = (1 )ln d ( ) (N Lv ) (N FR )
e
(3.156)
Segregado ascendente
0.011
- 3.768
3.539
- 1.614
Intermitente ascendente
2.96
0.305
- 0.4473
0.0978
Distribuido ascendente
Todos los patrones de
flujo descendente
Sin Correccin (C = 0)
4.70
- 0.3692
0.1244
- 0.5056
= 1 + 0.3 C
(3.157)
(3.158)
(3.159)
117
Captulo 3
m = L H L + g (1 H L )
(3.160)
18. Calcular la relacin del Factor de Friccin de las dos fases (fT) con respecto al
Factor de Friccin sin resbalamiento (fns).
fT
= eS ,
f ns
(3.161)
Donde:
S=
ln (y )
2
4
0.0523 + 3.182 ln (y ) 0.8725[ln (y )] + 0.01853 [ln (y )]
(3.162)
Y:
y=
[H L ( )]2
(3.163)
S se indetermina en un punto del intervalo 1 < y < 1.2; para y en este intervalo, la
funcin S se calcula de:
S = ln (2.2 y 1.2 )
(3.164)
N Re
2 log
4.5223
log
N
3
.
8215
Re
(3.165)
o:
f ns = 0.0056 +
0.5
(N Re )0.32
(3.166)
118
fT
f ns
(3.167)
Captulo 3
m v m v sg
p (144 )1
g c (p + 14.7 )(144 )
h =
g m
f M m v m 12
sen + T
gc
2 gc d
(3.168)
3.3 Ejemplos
3.3.1 Mtodo Poettman y Carpenter.
1. Dados los siguientes datos:
d = 1.995 pg.
T1 = 120 oF
T2 = 150 oF
Presin en la cabeza = p1 =
Presin de fondo = p2 =
Presin atmosfrica =
qo =
qw =
R=
o =
g =
w =
w =
o =
g =
500 psig
1,000 psig
14.7 psia.
400 bpd
600 bpd
500 pies3/bl
22 oAPI
0.65
1.07
70 dinas/cm. (cte)
30 dinas/cm. (cte)
0.018 cp. (cte)
119
Captulo 3
lb m
bl O
lb m
dia
(0.935)
514.7 520
= 5.88 + 8.43 + 35.69
= 50 pie 3
7. Densidad a 500 psig.
m1 =
947
lbm
= 18.96 3
50
pie
pie 3
bl
Z = 0.90
1014
.
7
520
1
120
Captulo 3
m2 =
947
lbm
= 29.8 3
31.8
pie
1.4737 10 5 (378,800 )
= 33.6
1.995/12
m =
m1 + m 2
2
lb
18.96 + 29.8
= 24.38 m3
2
pie
(0.008)(378,800)2
psi
1
dp
= 0.205
+
24
.
38
=
5
10
7.413 10 (24.38)(1.995 / 12)
pie
dh T 144
16.
p 2 p1 1000 500
=
= 2,439 pies
0.205
dp
dh T
121
Captulo 3
1.
141.5
141.5
=
= 0.9218
o
131.5 + API 131.5 + 22
2.
WOR =
q w 600
=
= 1.5
q o 400
lb
1.5
1
M = (0.9218)(350 )
+ (0.0764 )(500)(0.65) = 378.582 m
+ (1.07 )(350 )
bl o c.s.
1 + 1.5
1 + 1.5
3.
L = 62.4280.9218
4.
p=
lb
1
1.5
+ 1.07
= 63.0972 m3
2.5
2.5
pie
p1 + p 2
500 + 1000
+ p atm =
+ 14.7 = 764.7 psia
2
2
5.
T=
T1 + T2 120 + 150
=
= 135 0 F
2
2
6.
Z = 0.912
122
Captulo 3
7.
lb
520 764.7 1
= 2.4755 m3
pie
135 + 460 14.7 0.912
g = 0.65(0.0764 )
WOR
L = o
= 3.4313 cp
+ 0.5313
+ w
= 7.7813
2.5
2.5
1 + WOR
1 + WOR
11.
dinas
1
1.5
1
WOR
= 54
+ 70
+ w
= 30
cm
2.5
2.5
1 + WOR
1 + WOR
L =o
12.
Rs = 107.22 pies3/bl
13.
Bo = 1.0736
14.
d 2 (1.38)
=
= 1.4957 pg 2
4
4
2
Ap =
15.
N L
1
= 0.15726 L
3L
L
0.25
= 0.15726(3.4313)
3
(63.0972 )(54 )
0.25
= 9.6112 10 3
(86,400)(1.4957 )
seg
2.5 2.5
123
Captulo 3
17.
N Lv
= 1.938 v sL L
L
0.25
63.0972
= 1.938(6.44 )
54
0.25
= 12.9761
18.
v sg =
(144)(1000)500 107.22
pies
1 + 1.5 14.7 135 + 460 0.912
= 10.2179
86,400 (1.4957 )
seg
764.7 520 1
19.
N gv = 1.938 v sg L
L
20.
0.25
63.0972
= 1.938(10.2179 )
54
120.872(1.38) 63.0972
Nd =
12
54
0.25
= 20.5882
0.5
= 15.0256
21. De la figura 3.3 y con Ngv y NLv se tiene que el flujo se encuentra en la Regin II.
22. Para la Regin II el factor de resbalamiento se obtiene de la siguiente manera:
De la figura 3.6 se tiene que:
F5 = 0.19
F6 = - 0.1
F7 = 0.092
Y:
(N )
)
0.982
gv
+ F6'
(1 + F7 N Lv )2
= (1 + 0.19)
23.
vs =
24.
1.938 ( L / L )
0.25
4.7379
pie
= 2.3514
0.25
seg
1.938(63.0972 / 54)
124
0.5
Captulo 3
H L = 0.4210
Este valor de HL puede revisarse haciendo:
vs =
v sg
1 HL
v sL
10.2179
6.44
pie
=
= 2.35
H L 1 0.4210 0.4210
seg
( )
f1 R N d
2/3
v sg
v sL
10.2179
= 1.5866
6.44
= 0.0064(1.5866)(15.0256)
2/3
= 0.0618
f w = (f1 )
f2
1.05
= 0.0064
= 6.7126 10 3
f3
1.0011
Por lo tanto:
fr = 2 f w
N Lv (N Lv + N gv )
Nd
= 2 6.7126 10 3
+ 20.5882 )
)12.9761(1215.9761
= 0.3891
.0256
fr = 0.3891 Adimensional
27. El gradiente esttico es:
125
Captulo 3
g
2.4755
=
= 0.03923
L 63.0972
h =
p 1000 500
=
= 1,370.24 pies
dp
0.3649
dh T
126
120 oF
150 oF
500 psig
1,000 psig
14.7 psia.
400 bpd
600 bpd
500 pies3/bl
22 oAPI
0.65
1.07
70 dinas/cm. (cte)
30 dinas/cm. (cte)
0.018 cp. (cte)
Captulo 3
Se requiere la distancia entre los puntos de presin 500 psi y 1000 psi.
Solucin:
1. Peso especfico del aceite.
o =
141.5
= 0.9218
131.5 + 22
lbm
da
1
0.9218 (62.4) +
lb
1.5
5.614
L =
+ 1.07 (62.4)
= 62.01 m3
1.0736
pie
1 + 1.5
1 + 1.5
5. Calcular la presin.
Presin media = p =
500 + 1000
+ 14.7 = 764.7 psia
2
6. Calcular Z.
Temperatura media = T =
120 + 150
= 135 o F
2
Z = 0.91
g = (0.65) (0.0764)
127
Captulo 3
8. Rs media es:
Rs = 107.22
o = 7.7813 cp
10. La viscosidad media del agua es:
w = 0.5313 cp
11. Determinar la viscosidad lquida de la mezcla:
1.5
1
L = 7.7813
= 3.4313 cp
+ 0.5313
1 + 1.5
1 + 1. 5
4
1
= 0.0097
N L = 0.15726 3.4313
62.0121 543
Ap =
16. Obtener Bo a p y T . Captulo 1
B O = 1.0736
128
Captulo 3
v sL =
pie
seg
N LV
62.0121 4
= 1.938 3.08
= 6.1520
54
v sg =
pie
1 + 1.5 14.7 595 0.91
= 4.8784
86400 0.0217
seg
764.7 520 1
N Gv
62.0121 4
= 1.938 4.8784
= 9.7870
54
A = 1.071 -
0.166
129
Captulo 3
62.021
= 21.53
54
9.787 14.7
0.10
0.0020
= 2.29X10^ 4
21.53
2
1
.
53
(N Re ) TP =
130
2 . 2 10 2 (378585 )
1.995
0.2
0.018
3.4313
12
)(
(1 0 .2 )
= 9 . 74 x10 5
Captulo 3
0.00015
= 0.000902
1.995
12
Si el valor para no es conocido, un buen valor para usar es 0.00015 pies, que es un
valor promedio dado para el acero comercial.
30. De la figura 3.12:
f = 0.022
lb m
pie 3
T1 = 120F
Z1 = 0.9335
T2 = 150F
Z1 = 0.8976
B O1 = 1.0513
B O2 = 1.0974
R S1 = 73.0731 R S2 = 75.2569
VSL1 = 3.05
VSL2 = 3.11
VSG1 = 7.6572
VSG1 = 3.7278
( )
( ) [
131
Captulo 3
67.9754
144 (1000 - 500 ) 14.3872
64.4
h =
= 3555 pies
0.022 378585 2
14.3872 +
5
11 1.995
2.9652 10
14.3872
12
120 + 150
= 135 o F
2
132
Captulo 3
a)
p=
p1 + p 2
500 + 1,000
+ p atm =
+ 14.7 = 764.7 psia
2
2
Bo = 1.0736
Rs = 107.223 pies3/bl@ c.s.
Z = 0.912
b)
qg =
3.27 10 7 Z q o (R - R s ) (T + 460 )
p
q L = 6.49 10 5 (q o B o + q w B w )
q L = 6.49 10 5 [(400)(1.0736 ) + (600 )(1.0)] = 0.0668
pies 3
seg
c)
w L = 4.05 10 3 ((400 )(0.9218) + (600 )(1.07 )) + 8.85 10 7 (1,000 )(0.65)(107.223) = 4.1551 lb m /seg
g =
wg
qg
lb
0.2259
= 6.1958 m3
0.03646
pie
133
Captulo 3
L =
e)
lb
w L 4.1551
=
= 62.2021 m3
qL
0.0668
pie
Velocidades superficiales:
dh
(1.38)
Af =
=
= 1.4957 pg 2
4
4
2
v sg =
144 q L 144(0.0668)
pies
=
= 6.4312
Ap
1.4957
seg
144 q g
Ap
144(0.03646 )
pies
= 3.5102
1.4957
seg
pies
seg
La viscosidad del aceite (Beggs y Robinson), del agua y del lquido es:
o = 7.7813 cp
w = 0.5313 cp
1
1.5
WOR
1
L = o
+ 0.5313
= 3.4313 cp
= 7.7813
+ w
1 + WOR
2.5
2.5
1 + WOR
g)
= 30
+ 70
= 54
+ w
cm
1 + WOR
2.5
2.5
1 + WOR
L =o
h)
134
N Lv
= 1.938 v sL L
L
N gv
= 1.938 v sg L
L
0.25
62.2021
= 1.938(6.4312 )
54
0.25
62.2021
= 1.938(3.5102 )
54
= 12.9121
0.25
= 7.0475
Captulo 3
7.
2
2
vm
(
9.9414 )
L B = 1.071 2.6616
= 1.071 2.6616 1.38 = 189.5448
d
h
3.5102
= 0.3530
9.9414
L S = 50 + 36 N Lv = 50 + 36(12.9121) = 514.8356
Como
a)
v sg
vm
1
1
=
= 0.4
1 + WOR 1 + 1.5
fw = 1 fo =1-0.4 = 0.6
Por lo que el principal lquido que fluye en la tubera es agua (fase continua),
y como vm < 10 utilizamos la ecuacin 3.88:
d
= 0.681 + 0.013 h
12
1.38
= 0.681 + 0.013
12
1.38
1.38
d
log L + 0.232 log v m 0.428 log h
12
1.38
log (3.4313) + 0.232 log (9.9414 ) 0.428 log
= 1.4521
12
0.065(9.9414 )
135
Captulo 3
N Re L =
124 d h v b L 124(1.38)(1.75)(62.2021)
=
= 5,429
L
3.4313
13.59 L
v b = 0.5 + 2 +
0.5
(
)
d
/12
L
h
0.5
32.174(1.38)
= 0.251 + 8.74 10 (30,839 )
12
0.5
= 1.0012
Entonces:
13.59(3.4313)
v b = 0.5(1.0 ) + 1.0 2 +
0.5
62.2021(1.38/12 )
0.5
= 1.9868
124 d h v b L 124(1.38)(1.9868)(62.2021)
=
= 6,163.13
L
3.4313
1 C3
1 545.3658
psi
p
+ L =
+ 62.2021(1.4521) = 0.9447
pie
L e 144 v m + v b
136
Captulo 3
Con NRe =30,839 y con la rugosidad = 6 x 10-5 se obtiene el factor de friccin del
diagrama de Moody:
f = 0.02346
Y por lo tanto el gradiente de friccin es:
2
p 12 f v m L v sL + v b
+
2 (12) g d h v m + v b
h f
9.
Ek =
Wm v sg
4,637 p A p
4.381(3.5102 )
= 2.8996 10 6
4,637(764.7 )(1.4957 )
p p
psi
p
h e h f 0.9447 + 0.2920
=
= 1.2367
=
6
1 Ek
pie
1 2.8996 10
h T
Entonces L es:
h =
p
500
=
= 404.3018 pies
1.2367 1.2367
Captulo 3
WOR = 1.5
R = 500 pies3/bl
o = 22 oAPI
g = 0.65
w = 1.07
w = 70 dinas/cm.
o = 30 dinas/cm.
g = 0.018 cp.
1. Suponiendo una h = 2000 pies, entre p1 y p2 tenemos que:
p=
p1 + p 2
500 + 1000
+ p atm =
+ 14.7 = 764.7 psia
2
2
h=h+
2.
h
2,000
=0+
= 1,000 pies
2
2
w = 0.5214 cp
Z = 0.9131
Bw = 1.0
4.
o =
5.
o =
(350
141.5
141.5
=
= 0.9218
o
131.5 + API 131.5 + 22
+ 0.0764 R s g )
5.6146 B o
w =
lb
350(0.9218) + 0.0764(106.8419 )(0.65)
= 54.3529 m3
5.6146(1.0746 )
pie
350 w
lb
350(1.07 )
=
= 66.7011 m3
5.6146 B w 5.6146 (1.0)
pie
lb
1
WOR
1
1.5
+ w
= 54.3529
+ 66.7011
= 61.7618 m3
pie
1 + WOR
1 + WOR
2.5
2.5
L = o
138
Captulo 3
g =
6.
qg =
0.0764 g p 520
lb
0.0764(0.65)(764.7 )(520 )
= 2.4643 m3
14.7(137 + 460 )(0.9131)
pie
7.
pies 3
seg
d 2 (1.995)
=
= 3.1259 pg 2
4
4
2
Ap =
v sL =
v sg =
144 q L 144(0.0668)
pies
=
= 3.0773
seg
Ap
3.1259
144 q g
Ap
pies
144(0.03666 )
= 1.6888
seg
3.1259
pies
seg
9.
=
10.
qL
0.0668
=
= 0.6457
q L + q g 0.0668 + 0.03666
vm
(12)(4.7661) = 12(4.7661) = 4.2433
=
g d/12
32.2(1.995)
32.2(1.995)
2
N FR =
139
Captulo 3
1
1.5
1
WOR
L = o
= 3.2906 cp
+ 0.5214
= 7.4445
+ w
2.5
2.5
1 + WOR
1 + WOR
WOR
1
1.5
+ w
= 30
+ 70
= 54
cm
1 + WOR
1 + WOR
2.5
2.5
L =o
11.
N Re =
N Lv
12.
M m d/12
194.2213(1.995)
=
= 22,546.826
4
(2.1311) 6.72 10 4 (12)
m 6.72 10
= 1.938 v sL L
L
0.25
61.7618
= 1.938(3.0773)
54
0.302
0.25
= 6.1674
= 276.8958
2.4684
1.4516
6.738
= 2.7236 10 3
= 0.1887
= 9.5276
13. Dado que 0.4 y L3 < NFR L4, el patrn de flujo es intermitente.
14. El colgamiento horizontal para flujo intermitente es:
(0.845)(0.6457 )
a b
H L (0 ) =
=
c
(4.2433)0.0173
(N FR )
0.5351
= 0.6521
0.305
140
Captulo 3
lb m
pie 3
18.
y=
S=
0.6457
=
= 1.432
2
[H L (90)] (0.6715)2
ln(1.432)
= 0.3672
2
4
0.0523 + 3.182 ln(1.432) 0.8725[ln(1.432)] + 0.01853 [ln(1.432)]
fT
= e S = e 0.3672 = 1.4437,
f ns
19.
f ns =
22,546.826
2 log
(
)
4.5223
log
22,546.826
3
.
8215
= 0.0264
fT
= (0.0264 )(1.4437 ) = 0.0381
f ns
21.
m v m v sg
p (144 )1
g c (p + 14.7 )(144 )
h =
g m
f M m v m 12
sen + T
gc
2 gc d
h =
= 1,579.6 pies
(
0.0381) (194.2213) (4.7661)(12 )
(42.2826) (sen 90) +
(2) (32.2) (1.995)
141
Captulo 3
142
Captulo 4
CAPTULO 4:
FLUJO MULTIFSICO EN TUBERAS HORIZONTALES
4.1 Introduccin
Hay mucha literatura disponible sobre flujo multifsico en tuberas horizontales, lo
que hace muy complicado determinar cul de estas publicaciones ha contribuido ms al
desarrollo de informacin relacionada con este tipo de flujo.
Numerosos autores han presentado mtodos experimentales de clculo, conocidos
tambin como correlaciones para evaluar el gradiente de presin en tuberas horizontales. El
primer trabajo publicado sobre este tema fue en 1830, posteriormente ha habido
innumerables trabajos publicados dentro de los cuales hay 5 correlaciones generales que se
consideran las mejores:
De estas 5 correlaciones las mejores para todos los rangos de gastos y dimetros de
tubera son las de Dukler, Eaton y la de Beggs y Brill con la limitante de que para la de
Eaton se requieren viscosidades menores a 12 centipoise. Adicionalmente mencionaremos
que en la correlacin de Beggs y Brill puede ser usada para cualquier ngulo de flujo.
Debido a que para el flujo horizontal no se tiene el gradiente de elevacin es posible
que se piense que el colgamiento no sea necesario determinarlo, pero eso no es cierto, ya
que ste es necesario para calcular las velocidades verdaderas para el trmino de la
aceleracin, adems de que el colgamiento tambin est involucrado en la determinacin del
factor de volumen para algunas correlaciones
La mayora de las condiciones de flujo multifsico horizontal son en la regin de
flujo turbulento.
Para flujo horizontal, el gradiente de presin debido al cambio de elevacin es igual
a cero por lo que la ecuacin general de energa vista en el captulo 2 queda:
p
p
p
=
+
L T L ac L f
(4.1)
f v
v
p
+
=
L T 2 g c L 2 g c d
2
(4.2)
143
Captulo 4
La mayora de los investigadores han adoptado la ecuacin anterior para evaluar las
caractersticas del flujo de dos fases y posteriormente determinar el gradiente de presin
total. El problema de la variacin de las caractersticas de flujo se elimina al suponer que la
mezcla gas-lquido es homognea en un intervalo pequeo de la tubera. As la ecuacin 4.2
se puede escribir como:
v m 2 f tp m v m
p
+
= m
2 gc d
L T 2 g c L
(4.3)
4.2 Correlaciones
4.2.1 Lockhart y Martinelli
= L
L T
L L
(4.4)
p
2 p
= g
L T
L g
(4.5)
y:
X=
(p/L)L
(p/L)g
(p )L
(p )g
(4.6)
Donde:
p
144
Captulo 4
Procedimiento de clculo:
1. Calcular la cada de presin suponiendo que solamente fluye lquido en la tubera.
Con la ecuacin 2.35 y despreciando las prdidas por elevacin se puede calcular
dicha cada de presin en unidades prcticas.
f L q2 L
d5
p L = 1.1476 10 5
(4.7)
0.32
d1/3
(4.8)
El factor de friccin puede ser tambin obtenido del diagrama de Moody, o de las
ecuaciones vistas en el captulo 2 dependiendo del tipo de flujo que se tenga y de la
rugosidad de la tubera.
2. Calcular la cada de presin suponiendo que fluye solamente la fase gaseosa en la
tubera, que se puede calcular mediante la ecuacin de Weymouth que fue
modificada para incluir el factor de compresibilidad.
(T + 460)
q = 5,615.44 c.s.
p c.s.
(p
2
p2 d5
g (T + 460 ) L Z f
2
0.5
(4.9)
2
q 2 g (T + 460 ) L f Z
p 2 = p1 2.5343 10 11
d5
0.5
(4.10)
145
Captulo 4
(p/L)L
(p/L)g
(p )L
(p )g
(4.11)
4. Determine el nmero de Reynolds para ambas fases, suponiendo que cada una fluye
sola en la tubera, mediante la ecuacin 2.37 y la 2.66 vistas en el captulo 2.
5. Determinar el tipo de flujo de la siguiente tabla4.1:
Tabla 4.1. Tipo de flujo de acuerdo al nmero de Reynolds.
Nmero de Reynolds
Lquido
Gas
> 2000
> 2000
Lquido turbulento gas turbulento (tt)
< 1000
> 2000
Lquido laminar gas turbulento (vt)
> 2000
< 1000
Lquido turbulento gas laminar (tv)
< 1000
< 1000
Lquido laminar gas laminar (vv)
146
Captulo 4
= L
L T
L L
(4.12)
p
2 p
= g
L g
L T
(4.13)
4.2.2 Baker
Baker public una serie de artculos relacionados con el tema de flujo multifsico en
tuberas horizontales e inclinadas. En su trabajo inicial Baker describi siete diferentes
patrones de flujo y present un mtodo para predecir estos patrones. Usando su mtodo las
regiones de flujo ms precisas son la de bache y anular. Su mtodo en general es mejor para
dimetros de tubera mayores a 6 pulgadas, y la mayora de sus datos los tom de una
tubera de 8 y 10 pulgadas de dimetro. De un patrn de flujo a otro ocurre una
discontinuidad algo abrupta.
Bsicamente Baker present un acercamiento similar al de Lockhart y Martinelli, la
principal diferencia entre los dos es que Baker us el concepto de patrones de flujo y
tambin present diferentes ecuaciones para cada patrn.
Debido a que el cambio del patrn de flujo puede suceder en cualquier lugar de la
lnea, es mejor trabajar el problema escogiendo pequeos decrementos de presin, o sea,
escoger pequeos tramos de tubera en lugar de tomar la tubera completa, lo que nos
arrojara una presin media diferente y no tan precisa.
Procedimiento de clculo:
147
Captulo 4
f L q L
p
=
5
L L 1,359,947 d
5. Calcular el nmero de Reynolds para el gas con la ecuacin 2.66.
2
(4.14)
=
2 1010 d 5 p
L g
2
(4.15)
(p/L )L
(p/L )g
(4.16)
9. Calcular con:
g L
=
0.075 62.4
0.5
(4.17)
10. Calcular el flujo msico del gas en lbm/hr pie2 con la ecuacin:
Gg =
11. Determine
Gg
(0.0754) g q g (144)
(24) A p
(4.18)
q L (5.6146) L (144)
24 A p
(4.19)
13. Calcular .
73 62.428
L
=
L L
148
2 1/ 3
(4.20)
Captulo 4
14. Calcular:
15. Usando
GL
Gg
Gg
(4.21)
GL
obtener el patrn de flujo de la figura 4.3
Gg
gtt =
1402 X 0.75
L0.1
(4.22)
b) Flujo Tapn:
gtt =
27.315 X 0.855
L0.17
(4.23)
c) Flujo Estratificado:
gtt =
15,400 X
0.8
GL
(4.24)
p T =
fT Gg L
193.2 d g
(4.25)
y:
f TP = f g g
(4.26)
Donde:
g = factor de correlacin de Lockhart y Martinelli .
e) Flujo Bache:
gtt =
f) Flujo Anular:
1,190 X 0.815
0.5
GL
(4.27)
(4.28)
Esta ecuacin podr ser utilizada para dimetros de tubera menores a 12 pg.
Cuando d es mayor a 10 pg, siempre use 10 pg en la ecuacin.
149
Captulo 4
150
Captulo 4
(4.29)
(4.30)
3. Para 1.0 X 10
(4.31)
4. Para 10 X 100
(4.32)
p T = gtt p g
2
(4.33)
18. Calcular L.
L =
p1 p 2
p
L T
(4.34)
151
Captulo 4
Para todas las cadas de presin (p ) , los valores de L son calculados siguiendo
los pasos del 1 al 18. Todos los L son sumados hasta obtener la longitud total de la lnea.
4.2.3 Dukler
En 1964 Dukler public su trabajo sobre flujo multifsico horizontal y
posteriormente en 1969 un manual. Acumul todos los datos publicados sobre este tema y
form lo que ellos llaman un banco de datos, los cuales consistan en datos de laboratorio de
tubera corta y datos de campo de largos tramos de tubera con aceite.
La correlacin presentada por Dukler consiste esencialmente en dos partes: caso I y
caso II.
Caso I
1. Suponer la cada de presin corriente abajo que puede ser para toda la longitud de la
lnea o solo para una distancia corta, y con sta, calcular la presin promedio entre p1
y p2.
2. Obtener R s , B o , Z.
3. Calcular el gasto de lquido y gas en pies3/da.
q L = q L B o 5.6142
(4.35)
p (T + 460 )
q g = q L (R R s ) c.s.
Z
p (Tc.s + 460 )
(4.36)
L =
(4.37)
T + 460 p c.s. Z
g = g aire
c .s .
152
(4.38)
Captulo 4
wm =
L qL + g qg
86,400
(4.39)
w m (144 )
Ap
(4.40)
m = L + g (1 )
(4.41)
(4.42)
wT
(N Re )T = 1,488
d
(4)(12)
(4.43)
0.125
(N Re )T 0.32
(4.44)
2 fT GT
p
=
L f g c m d (12)
(4.45)
(20736)(16) w T w G p
2 g c d 4 p1 p 2 g
(4.46)
153
Captulo 4
L f
p
=
L T (1 a )
(4.47)
(4.48)
Caso II
1. Suponer la cada de presin p y calcular la presin promedio, p .
p=
p1 + p 2
2
(4.49)
2. Obtener R s , B o , Z.
3. Calcular el gasto de lquido y gas en pies3/da.
q L = q L B o 5.61
(4.50)
p (T + 460 )
q g = q L (R R s ) c.s.
Z
p (Tc.s + 460 )
(4.51)
154
Captulo 4
wL
qL
L
=
qL + qg w L w g
+
(4.52)
L =
(4.53)
T + 460 p c.s. Z
g = g aire
c .s .
(4.54)
(q
+ q g ) (144 )
(4.55)
d2
86,400
4
m = L + g (1 )
(4.56)
m = L
+ g
(1 )2
1 HL
(4.57)
(N Re )T =
d vm m
(12) 6.72 10 -4 T
(4.58)
12. Con el colgamiento sin resbalamiento calculado en el paso 4 y (NRe)T del paso 11,
ir a la figura 4.4 y leer el valor del colgamiento H L .
155
Captulo 4
fT
.
fo
0.125
(N Re )T 0.32
(4.59)
fT
fo
fo
(4.60)
156
(4.61)
Captulo 4
Esto puede cambiarse para calcular p/L o resolverse para L , sobre una cada de
presin. Si la p1 es conocida, el valor de p2 puede ser supuesto, y L solucionarse
directamente, despus sumar todas las L hasta completar el largo de la lnea
L = londitud de la lnea .
18. La cada de presin debido a la aceleracin puede ser despreciable dentro de la
tubera, pero puede ser considerada en procesos de instalacin de tuberas. Baker dio
la siguiente ecuacin:
1
p a =
144 g c A 2
q 2
g q g 2 q 2
L qL2
g g
L
L
+
cos
(
H L
1 HL )
H L
(1 H L )
para p 2
para p1
(4.62)
(4.63)
20. Si sucede algn cambio en la elevacin, agregue la componente por prdida por
elevacin en el paso 19.
157
Captulo 4
4.2.4 Eaton
Tres lquidos fueron probados en cada lnea. El gasto de lquido vari de 50 a 2,500
barriles por da en la lnea de 2pg y de 50-5,500 bpd en la de 4 pg, y para cada gasto de
lquido la relacin gas-aceite se vari desde cero a el mximo permitido por el sistema.
La precisin del mtodo para determinar la presin en algn punto de la tubera
depender de las magnitudes de los decrementos de presin tomadas, entre ms pequeos
sean los decrementos de presin aumentar la precisin del clculo.
Procedimiento de clculo:
(N Lv )0.575
N gv N d
0.0277
pb
0.05
NL
N LB
0.1
(4.63)
Donde:
N Lv = 1.938 v sL L
158
0.25
(4.64)
Captulo 4
N gv = 1.938 v sg L
0.25
120.872 d L
Nd =
12
(4.65)
0.5
N L = 0.15726 L
3
L
(4.66)
0.25
(4.67)
p
p
=
p c.s. 14.65
(4.68)
(4.69)
Es importante notar que el gas en solucin y el gas libre deben determinarse para
poder evaluar correctamente vsL y vsg.
5. Obtener HL1 y HL2 de la figura 4.6.
6. Evaluar v L1 , v L2 , v L , v g1 , v g2 , v g .
7. Obtener el valor del Factor de Friccin de la figura 4.7, calculando primero:
(GR )
1.25
0.5
dB
d
GT d
g (12) 6.72 10 4
(4.70)
Donde:
GR =
wm
(4.71)
wL
wm
(4.72)
wm
A p / 144
(4.73)
LR =
GT =
wg
159
Captulo 4
Con este valor obtenga el valor de f(LR)0.1 de la figura 4.7, y conociendo LR calcular
el Factor de Friccin f.
w L w g w L v L 2 + w g v g 2
+
144 p
2 gc
L g
(4.74)
160
Captulo 4
161
Captulo 4
horizontal. Una vez establecido cada patrn de flujo se procedi a variar el ngulo de
inclinacin, as que se pudo observar como el ngulo de inclinacin afectaba el colgamiento
y el gradiente de presin. El colgamiento y el gradiente de presin fueron medidos en
ngulos que variaban de 5, 10,15, 20, 35, 55, 75 y 90 grados, y se encontr que el
colgamiento llagaba a su valor mximo en +50 grados y a su valor mnimo en -50 grados. El
mapa de patrones de flujo original que obtuvieron Beggs y Brill fue ligeramente modificado
para poder incluir la zona de transicin entre el patrn de flujo segregado y el intermitente.
El mapa de patrones de flujo modificado fue sobrepuesto al original y se muestra en la
figura 4.8.
g
f G v
m sen + T m m
2 g c d 12
p g c
=
m v m v sg
L
1g c (p + 14.7 )(144)
Donde p/L esta en: psi/pie, y para flujo horizontal sen = 0.
Procedimiento de clculo:
162
(4.75)
Captulo 4
1. Comenzando con una presin conocida p1, estimar el valor para la cada de presin
p.
2. Calcular la presin promedio en el intervalo:
p = p1 +
p
si p1 es la presin corriente abajo.
2
(4.76)
p = p1
p
si p1 es la presin corriente arriba.
2
(4.77)
R s , Bo , B w , o , w , g , o , w , Z
4. Calcular la densidad relativa del aceite:
o =
141.5
131.5 + o API
(4.78)
WOR
+ w
1 + WOR
1 + WOR
L = o
(4.79)
o =
(350
+ 0.0764 R s g )
5.6146 B o
w =
g =
350 w
5.615 B w
(4.80)
(4.81)
(4.82)
3.27 10 7 Z q o (R - R s ) (T + 460)
(p + 14.7 )
(4.83)
163
Captulo 4
(4.84)
Donde:
qL y qg = pies3/seg
7. Calcular las velocidades superficiales del gas, lquido y la mezcla in situ.
v sL =
v sg =
144 q L
Ap
(4.85)
144 q g
(4.86)
Ap
v m = v sL + v sg
(4.87)
(4.88)
G g = g v sg
(4.89)
GT = GL + Gg
(4.90)
N FR
Vm
=
g d/12
(4.91)
WOR
L = o
+ w
1 + WOR
1 + WOR
m = L + g (1 )
WOR
+ w
1 + WOR
1 + WOR
L =o
(4.92)
(4.93)
(4.94)
Captulo 4
N Re =
N Lv
G T d/12
m 6.72 10 4
= 1.938 v sL L
L
(4.95)
0.25
(4.96)
12. Para determinar el patrn de flujo que existe en el flujo horizontal, calcular los
parmetros correlacionados, L1, L2, L3, y L4 con:
L1 = 316 0.302
(4.97)
L 2 = 0.0009252 2.4684
(4.98)
L 3 = 0.10 -1.4516
(4.99)
L 4 = 0.5 -6.738
(4.100)
Transicin
Intermitente
Distribuido
a b
(N FR )c
(4.101)
165
Captulo 4
0.98
0.4846
0.0868
Intermitente
0.845
0.5351
0.0173
Distribuido
1.065
0.5824
0.0609
(4.102)
Donde:
S=
ln(y )
2
4
0.0523 + 3.182 ln(y ) 0.8725[ln(y )] + 0.01853 [ln(y )]
(4.103)
Y:
y=
[H L (0)]2
(4.104)
(4.105)
N Re
2 log
4.5223
log
N
3
.
8215
Re
(4.106)
o:
f ns = 0.0056 +
18. Calcular el factor de friccin de las dos fases.
166
0.5
(N Re )0.32
(4.107)
Captulo 4
f T = f ns
fT
f ns
(4.108)
19. Calcular:
f G v
L T m m
2 g c d 12
p =
m v m v sg
1g c (p + 14.7 ) (144 )
(4.109)
4.3 Ejemplos
4.3.1 Mtodo Lockhart y Martinelli.
Determinar la presin corriente abajo para una presin corriente arriba conocida.
Datos disponibles:
Dimetro interno de la tubera = 2 pg
Longitud de la lnea = 1,500 pies
qL = 2,000 bpd de agua.
w = 1.07
RGL = 1,000 pies3/bl
g = 0.65
Presin corriente arriba = p1 = 850 psia.
w = 66.7 dinas/cm.
g = 0.015 cp.
T = 120 oF.
L = 1.0 cp.
1. Calcular la cada de presin suponiendo que solamente fluye lquido:
f=
p L = 1.1476 10 5
0.032
= 0.0254
21/3
167
Captulo 4
2 10 6
p 2 = 850 2 2.5343 10 11
) (0.65)(120 + 460)(1500)(0.0254)(0.8937)
2
25
0.5
= 825.7 psia
Entonces:
p g = p 2 p1 = 850 825.7 = 24.3 psia
3. Calcular el parmetro X:
X=
58.48
= 1.55
24.3
4. Los nmeros de Reynolds para ambas fases, suponiendo que cada una fluye sola en
la tubera son:
(2000)(1.07 ) = 98,654
(2)(1.0)
a)
N ReL = 92.2
b)
N Reg = 0.0201056
(2 10 )(0.65) = 871,242
6
(2)(0.015)
L = 3.5
g = 5.4
7. A cada de presin de las dos fases es:
168
Captulo 4
p
2
p
2
T
Por lo tanto:
p 2 = 850 708.59 = 141.41 psia.
4.3.2 Mtodo Baker
500 + 460
= 480 psi = 494.7 psia
2
Adems se obtuvieron:
g = 1.6
lb m
,
pie 3
R s = 124.84
pies 3
,
bl
B o = 1.0762,
Z = 0.9359
Captulo 4
L =
Donde:
gd =
o + gd
Bo
o =
Por lo que:
L =
141.5
141.5
=
= 0.8155
API + 131.5 42 + 131.5
pie 3
1.0764
y:
N ReL = 1.4781
(2000)(48.9739) = 36,194
(4)(1)
5
pie
1,359,947 (4)
L L
2
170
qg g
d g
Captulo 4
Por lo tanto:
N Reg = 0.0201056
(4)(0.02)
=
5
pie
2 1010 (4) (494.7 )
L g
2
8.
X=
(p/L )L
(p/L )g
9.5163 10 4
= 1.3909
4.9185 10 4
9.
1.6 48.9739
=
0.075 62.4
0.5
= 4.0918
Mg =
lb m
0.0764 (0.65) 1.750320 10 6 (144)
= 41,517.60
(24)(12.5637 )
hr pie 2
11.
Gg
41510.49
= 10,144.80
4.0918
lb m
2000 (5.6146 )(48.9739 )(144)
= 262,631.733
24(12.5637 )
hr pie 2
13.
171
Captulo 4
2
73 62.428
= 1
30 48.9739
14.
1/ 3
= 2.8607
G L (262631.733)(4.0918(2.8607 ))
=
= 74.046
Gg
41517.6
gtt = (4.8 0.3125 d ) X 0.343-0.021d = (4.8 0.3125 (4 ))1.3909 0.343 - 0.021(4 ) = 3.8667
17.
p
p
4
3 psi
=
gtt = 4.9185 10 (3.8667 ) = 1.9018 10
pie
L T L g
Entonces la cada de presin total para toda la lnea es:
p1 p 2
500 460
=
= 21,032 pies
1.9018 10 3
p
L T
Captulo 4
g = 0.65
o = 42 oAPI
Presin corriente arriba = 500 psi.
o = 30 dinas/cm.
o = 1 cp.
g = 0.02 cp.
/d = 1.0 x 10-4
Relacin agua aceite WOR = 0
1. Suponiendo una p = 40 psi. tenemos que:
p 2 = p1 p = 500 - 40 = 460 psi
p=
500 + 460
= 480 psi = 494.7 psia
2
2.
g = 1.6
lb m
,
pie 3
R s = 124.84
pies 3
,
bl
B o = 1.0762,
Z = 0.9359
3.
q L = (2000 )(1.0762)(5.6142 ) = 12,084
pies 3
da
pies 3
14.7 120 + 460
q g = 2000(1000 124.84 )
0.9359 = 54,293.39
da
494.7 60 + 460
4.
=
qL
12,084
=
= 0.1820
q L + q g 12,084 + 54,293.39
5. Calcular la relacin del flujo msico del lquido y del gas total, WT.
a) Densidad del lquido:
o =
o =
141.5
141.5
=
= 0.8155
API + 131.5 42 + 131.5
173
Captulo 4
= 1.601 m3
pie
(120 + 460) 14.7 0.9359
g = (0.65)(0.0764 )
L qL + g qg
86,400
seg
6.
Ap =
Gm =
d2
42
=
= 12.5664 pg 2
4
4
7.
lb m
pie 3
9.
wT
(N Re )T = 1,488
d
(4)(12)
(1,488)(7.7657 ) = 222,583
(4)
(0.1983)
(4)(12)
10.
f T = 0.00140 +
11.
0.125
0.125
= 0.00140 +
= 2.4306 10 3
0.32
0.32
(N Re )T
(222,583)
2 fT GT
psi
2 2.4306 10 3 (88.9881)
p
=
=
= 2.4645 10 3
(32.2)(10.1059)(4)(12)
pie
L f g c m d
2
174
Captulo 4
13.
p
3
psi
p
L f 2.4645 10
=
= 2.568 10 3
=
(1 a ) (1 0.0403)
pie
L T
14. La cada de presin en toda la lnea es:
p
3
p = L
= (3000 ) 2.568 10 = 7.7039 psi
L T
Como la p calculada no es igual a la p supuesta se debe suponer la nueva p
calculada y comenzar desde el paso nmero 1.
o = 52.2 lb m /pie 3
g = 1.45 lb m /pie 3
o = 4.2 cp
g = 0.0105 cp
o = 22.3 dinas/cm
d = 12 pg
L = 134,370 pies
p1 = 424.7 psia
1. Suponiendo p2 = 319.7 psia:
p=
p1 + p 2 424.7 + 319.7
=
= 372.2 psia
2
2
175
Captulo 4
qg = 105,600 pies3/da
4.
=
qL
40000
=
= 0.2747
q L + q g 40000 + 105600
5.
(q
+ q g ) (144 )
86,400
8.
d
4
m = L
+ g
11.
(N Re )T =
(1 0.2747 )2
(0.2747 )2
(1 )2
lb
= 14.219 m3
+ 1.45
= 52.2
1 0.30
0.30
1 HL
pie
d vm m
(12)(2.1456)(14.219) = 39,093
=
-4
(12) 6.72 10 m (12) 6.72 10 4 (1.1613)
2
m = L
HL
176
+ g
(1 )2
(0.2747 )2
(1 0.2747 )2
lb m
= 52.2
+ 1.45
1 HL
0.43
1 0.43 = 10.4987 pie 3
Captulo 4
(N Re )T =
d vm m
(12)(2.1456)(10.4987 ) = 28,865
=
-4
(12) 6.72 10 m (12) 6.72 10 4 (1.1613)
Usando la figura 4.4 con estos nuevos valores obtenidos tenemos que HL = 0.45, que
es el valor que se utilizar.
14. De la figura 4.5 se obtuvo:
fT
= 2.00
fo
15. Calcular:
f o = 0.00140 +
0.125
0.125
= 0.00140 +
= 6.0729 10 3
0.32
0.32
(N Re )T
(28,865)
16.
fT =
17.
fT
f o = (2.00) 6.0729 10 3 = 0.0121458
fo
Captulo 4
p=
850 + 650
= 750 psia
2
T = 120 oF. (dada)
3.
= 2.5239 m3
pie
(120 + 460) 14.7 0.90
b) L = (1.07 )(62.428) = 66.7979
c) Rs = 0 (se supone cero debido a que es agua)
d) Bw = 1.0 (se supone que no hay compresibilidad)
e) w = 66.7 dinas/cm
f) g = 0.015 cp
g) w = 1.0 cp
h) w L =
i) w g =
Captulo 4
p (T + 460 )
pies 3
14.7 580
q g = q L (R R s ) c.s.
Z = 2 10 6
(0.92 ) = 35,492.85
da
850 520
p (Tc.s + 460 )
v sg =
Donde:
q g 144
86,400 A p
(35,492.85)(144) = 18.83
86,400 (3.1416 )
pies
seg
d 2 (2 )
Ap =
=
= 3.1416 pg 2
4
4
2
v sL =
0.25
= 1.938 v sg L
Nd =
0.25
120.872 d L
12
N L = 0.15726 L
3
L
0.25
66.7979
= 1.938(5.9573)
66.7
66.7979
= 1.938(18.83)
66.7
0.5
120.872(2) 66.7979
12
66.7
0.25
= 11.5494
0.25
= 36.5059
0.25
= 20.1527
= 0.15726(1.0)
3
(66.7979 )(66.7 )
0.25
= 0.002357
p
p
850
=
=
= 57.8231
p c.s. 14.7 14.7
NLB = constante = 0.00226
Por lo tanto:
(N Lv )0.575
N gv N d
0.0277
pb
0.05
NL
N LB
0.1
Captulo 4
da
650 520
p (Tc.s + 460)
v sg =
q g 144
86,400 A p
Donde:
(46,010.13)(144) = 24.409
86,400 (3.1416 )
pies
seg
d 2 (2 )
=
= 3.1416 pg 2
4
4
2
Ap =
v sL =
N Lv
= 1.938 v sL L
N gv
= 1.938 v sg L
Nd =
0.25
120.872 d L
12
N L = 0.15726 L
3
L
0.25
66.7979
= 1.938(5.9573)
66.7
66.7979
= 1.938(24.409 )
66.7
0.5
120.872(2) 66.7979
12
66.7
180
= 11.5494
0.25
= 47.322
0.25
= 20.1527
= 0.15726(1.0)
3
(66.7979 )(66.7 )
p
p
650
=
=
= 44.2176
p c.s. 14.7 14.7
Por lo tanto:
0.25
0.25
= 0.002357
Captulo 4
(N Lv )0.575
0.0277
N gv N d
pb
0.05
NL
N LB
0.1
1
1
pies
= 5.9573
= 19.8576
H L1
0.30
seg
v L2 = v sL2
1
1
pies
= 5.9573
= 22.0641
H L2
0.27
seg
v g2 = v sg2
1
1
pies
= 18.83
= 26.9
(1 - H L1 )
0.70
seg
1
1
pies
= 24.409
= 33.4369
(1 - H L2 )
0.73
seg
wg
wT
1.1495
= 0.1168
9.837
dB 1
= = 0.5
d 2
GT =
wm
lb m
9.837(144)
=
= 450.894
A p / 144
3.1416
seg - pie 2
Por lo tanto:
(GR )
0.5
1.25
dB
d
(450.894)(2)
Gm d
0.5
1.25
= (0.1168) (0.5)
4
(0.015)(12)(6.72 104 )
g (12) (6.72 10 )
181
Captulo 4
(GR )
0.5
dB
d
1.25
Gm d
= 1.071 10 6
g (12 ) 6.72 10 4
w L 8.6815
=
= 0.8825
w T 9.837
L =
0.022
= 0.0223
0.9875
v m1 + v m2 24.7873 + 30.3663
pies
=
= 27.5768
2
2
seg
2 gc d
2
12 w T v m f
2
2
w
w g w L v L + w g v g
144 p L +
2
g
L
g
c
L =
2(32.2)(2 )
2
12(9.837 )(27.5768) (0.0223)
2
2
2(32.2 )
66.7979 2.5239
L = 1,084.66 pies
Como se puede observar se considera despreciable el trmino de la energa
cinemtica.
Este valor obtenido al compararlo con la longitud total de la tubera (L = 1500 pies)
se nota que es diferente, por lo que se debe suponer otra cada de presin. El valor
correcto final es de 580 psia.
Se puede obtener un resultado ms preciso si se suponen decrementos de presin ms
pequeos.
182
Captulo 4
p = p1
p
= 500 50 = 450 psi
2
3.
Rs = 100, Bo = 1.065, Z = 0.94
4.
o =
5.
o =
g =
6.
qg =
141.5
= 0.82
131.5 + 42
183
Captulo 4
2 2
d = 4 = 12.5663 pg 2
4
4
Ap =
v sL =
pies 3
seg
v sg =
144(0.6905)
pies
= 7.9126
12.5663
seg
pies
seg
8.
G L = L v sL = (48.8274 )(1.5814 ) = 77.2156
G g = g v sg = (1.4973)(7.9126 ) = 11.8475
lb m
seg pie 2
lb m
seg - pie 2
lb m
seg - pie 2
9.
=
qL
0.1382
=
= 0.1667
q L + q g 0.1382 + 0.6905
10.
Vm
(9.494) = 8.3977
=
g d (32.2)(4 / 12)
2
N FR =
184
M m d/12
(89)(4 / 12)
=
= 2.4084 10 5
4
m 6.72 10
0.1833 6.72 10 4
Captulo 4
N Lv
12.
= 1.938 v sL L
L
0.25
48.8274
= (1.938)(1.5814 )
30
0.302
0.25
= 3.4616
= 183.955
2.4684
1.4516
6.738
= 0.077
= 1.3472
= 87,411.37
13. Como 0.01 < 0.4 y L3 < NFR L1, el flujo es intermitente.
14.
H L (0) =
15.
a b
0.845(0.1667 )
=
c
(N FR )
(8.3977 )0.0173
0.5351
= 0.3123
16.
y=
S=
0.1667
=
= 1.7092
2
[H L (0)] (0.3123)2
ln (1.7092 )
= 0.3817
2
4
0.0523 + 3.182 ln(1.7092 ) 0.8725[ln (1.7092 )] + 0.01853 [ln(1.7092 )]
fT
= e S = e 0.3817 = 1.4648
f ns
17.
f ns = 0.0056 +
0.5
0.5
= 0.0056 +
0.32
N Re
2.4084 10 5
0.32
= 0.015
18.
f T = f ns
fT
= 0.015(1.4648) = 0.022
f ns
185
Captulo 4
19.
f Mm vm
(0.022 )(89.0631)(9.494 )
Z T
3000
2 g c d 12
2(32.2)(4)(12 )
= 18.064 lb f
p =
=
(16.2784)(9.494)(7.9126)
m v m v sg
pg 2
1
186
Captulo 5
CAPTULO 5:
FLUJO MULTIFSICO EN TUBERAS INCLINADAS
5.1 Introduccin
El flujo inclinado se define como el flujo a travs de tuberas que se desvan a partir
de la horizontal o como el flujo a travs del terreno accidentado. El flujo direccional se
define como el flujo a travs de tuberas que se desvan con respecto a la vertical, y es
referida, como el flujo a travs de tubera en pozos de perforacin direccional. Ambos
ofrecen problemas similares pero se discuten por separado.
La prdida de presin total en la tubera de descarga debido a terreno montaoso o
accidentado, es la suma de las prdidas por friccin, aceleracin y por elevacin necesarias
para transportar los fluidos a lugares con mayor elevacin a cualquier distancia.
Para este captulo se puede aplicar la misma ecuacin general (ecuacin 2.11)
p p p
p
=
+
+
L T L e L ac L f
(5.1)
o:
p
g v v m m f v m
= m sen + m m
+
L
gc
gc
z
2 gc d
(5.2)
5.2 Correlaciones
5.2.1 Correlacin de Flanigan
187
Captulo 5
fueron de 1 a 12 pies por segundo y el flujo de lquido a gas fue de 20 a 1200 bbl/MMcf.
Los fluidos usados fueron gas natural y condensado. Una vez obtenida la eficiencia en la
tubera, podemos usar la ecuacin de cada de presin similar a la de Panhandle para
determinar el componente de la cada de presin.
Examinando los datos de prueba para tubera de 16 pulgadas Flanigan not lo
siguiente:
a)
b)
c)
d)
e)
El procedimiento sugerido por Flanigan para calcular la cada de presin total es como
sigue.
1. Calcula la velocidad superficial del gas, v sg aplicando la frmula siguiente:
qg Z T
v sg = 0.031194 2
d p TSC
(5.3)
4. Con la eficiencia del paso anterior calcular la cada de presin por friccin,
utilizando la ecuacin de Panhandle o una similar.
1.07881
2
2
Tatm + 460
P1 P2
q g = 435.87
L
P
atm
Z (T + 460 )
5277
188
0.5394
d2
0.4606 E (5.4)
Captulo 5
P =
L HF H
144
(5.5)
189
Captulo 5
Beggs y Brill condujeron una investigacin experimental de dos fases, flujo de gaslquido en tuberas inclinadas para determinar que efecto tiene el ngulo de inclinacin en
el colgamiento del lquido y las perdidas de presin.
El estudio experimental consisti en mediciones de colgamiento del lquido y cadas
de presin en tuberas de 1 y 1.5 pulgadas de dimetro interno. Los fluidos usados fueron
aire y agua con flujos de variaron de 0 a 4 pies 3 /min . De lquido y de 0 a 300 Mpies 3 /dia
gas. Los datos tomados fueron para las inclinaciones de ngulos de:
90, 85, 75, 55, 35, 20, 15, 10, 5 y 0 a partir de la horizontal.
Las correlaciones empricas fueron desarrolladas para el colgamiento del lquido y
factor de friccin como funcin de las propiedades de flujo y ngulo de inclinacin. El
colgamiento del lquido es fuertemente influenciado por el ngulo de inclinacin y el factor
de friccin es influenciado por el colgamiento del lquido. Se desarrollaron diferentes
correlaciones para el colgamiento cerca de los tres regmenes de flujo. La correlacin del
factor de friccin se hizo para ser independiente del rgimen de flujo, pero requiere de un
valor para del colgamiento del lquido.
Beggs y Brill llegaron a las siguientes conclusiones.
El ngulo de inclinacin de una tubera con flujo de dos fases, ocurren
efectos de colgamiento de lquido y cadas de presin.
En el flujo inclinado de dos fases, el mximo colgamiento del lquido es con
b)
un ngulo aproximado de +50 y el mnimo aproximadamente -50 respecto a la
horizontal. El hecho de que el colgamiento sea aproximadamente igual en los
ngulos de +90 y +20 explica por que las correlaciones de colgamiento
vertical pueden ser usadas para algunos grados de flujo horizontal.
La recuperacin de presin cuesta abajo en la seccin de la tubera con dos
c)
fases que esta sobre el terreno existe y puede ser considerada en el diseo de la
tubera.
a)
g
f T G m v m 12
g m sen +
2 gc d
p c
=
L
v v
144 1 - m m sg
g c p (144 )
m = L H L + g (1 - H L )
(5.6)
(5.7)
190
Captulo 5
p = p1 +
(5.8)
141.5
131.5 + API
(5.9)
WOR
+ w
= o fo + w f w
1 + WOR
1 + WOR
L = o
o =
350 o + 0.0764 R s g
5.615 B o
w =
g =
350 w
5.615 B w
0.0764 g p (520 )
(14.7 ) (T + 460) Z g
(5.10)
(5.11)
(5.12)
(5.13)
3.27 10 7 Z g q o (R Rs ) (T + 460 )
p
q L = 6.49 10 5 (q o B o + q w B w )
(5.14)
(5.15)
Donde: qL y qg = pies3/seg
7. Calcular las velocidades superficiales del gas, lquido y la mezcla:
191
Captulo 5
v sL =
v sg =
144 q L
Ap
144 q g
Ap
v m = v sL + v sg
(5.16)
(5.17)
(5.18)
(5.19)
G g = g v sg
(5.20)
Gm = GL + Gg
(5.21)
qL
qL + qg
(5.22)
N FR
v
= m
d
g
12
1
WOR
L = o
= O fO + W fW
+ w
1 + WOR
1 + WOR
m = ( L + g (1 ))
1
WOR
L = o
= O fO + W fW
+ w
1 + WOR
1 + WOR
(5.23)
(5.24)
(5.25)
(5.26)
N Rens
192
d
12
=
m 6.72 10 4
Gm
(5.27)
Captulo 5
N LV
= 1.938 v sL L
L
0.25
(5.28)
12. Para determinar el patrn de flujo que existe en el flujo horizontal, calcular los
parmetros correlacionados L1, L2, L3, y L4.
L1 = 316 0.302
(5.29)
L 2 = 0.0009252 2.4684
(5.30)
L 3 = 0.10 -1.4516
(5.31)
L 4 = 0.5 -6.738
(5.32)
13. Determinar el patrn de flujo usando los siguientes lmites de la tabla 5.1.
Segregado
0.01 y NFR < L2
0.01 y L2 < NFR L3
Transicin
Intermitente
Distribuido
H L (0) =
a b
(5.33)
(N FR )c
0.98
0.4846
0.0868
Intermitente
0.845
0.5351
0.0173
Distribuido
1.065
0.5824
0.0609
193
Captulo 5
C = (1 ) ln d e (N Lv
)f (N FR )g )
(5.34)
Segregado ascendente
0.011
- 3.768
3.539
- 1.614
Intermitente ascendente
2.96
0.305
- 0.4473
0.0978
Distribuido ascendente
Todos los patrones de
flujo descendente
Sin Correccin (C = 0)
4.70
- 0.3692
0.1244
- 0.5056
(5.35)
(5.36)
18. Calcular la relacin del Factor de Friccin de las dos fases (fT) con respecto al
Factor de Friccin sin resbalamiento (fns).
fT
= eS ,
f ns
(5.37)
Donde:
S=
ln (y )
+ 0.01853 [ln (y )]
(5.38)
Y:
y=
[H L ( )]2
(5.39)
S se indetermina en un punto del intervalo 1 < y < 1.2; para y en este intervalo, la
funcin S se calcula de:
194
Captulo 5
S = ln (2.2 y 1.2)
(5.40)
N Rens
2 log
4.5223 log N Rens 3.8215
(5.41)
fT
f ns
(5.42)
21. Calcular:
g
f G m v m 12
L m sen + T
gc
2 gc d
p =
m v m v sg
(144) 1
g c p (144 )
(5.43)
5.3 Ejemplos.
5.3.1 Mtodo Flanigan
Dados los siguientes datos:
d=
E=
L=
P1 =
RGL =
qL =
w =
w =
g =
2 pg.
Eficiencia de Panhandle, %
1500 pies
850 psia.
1000 pie3/bl.
2000 bpd
1.07
66.7dinas/cm
0.015 cp
T = 120 F
195
Captulo 5
L = 1.0 cp
66.7 lbm/pie^3
L =
g = 0.65
# Colina
1
2
3
4
5
6
altura
15
10
160
100
55
En este ejemplos solo se mostrar los efectos de las colinas y la prdida de presin
total no ser calculada.
1. Calcular Vsg.
qg =
Suponiendo
Tenemos
seg
p2 = 650 psia
p = 750 psia
Z = 0.9
pie
14.7 120 + 460 0.9
q g = (23.1481)
= 0.4554
seg
750 520 1
v sg =
qg
Ap
0.4554
pie 3
=
2
0
.
874
seg
(2)2
4(144 )
2. De la figura 5.2
HF = 0.13
Por efecto de las colinas:
P =
Se calcular la prdida de presin por friccin para 1500 pies por el mtodo de
Eaton (del mismo ejemplo al del captulo 4).
196
Captulo 5
p
= 850 580 = 270 psia
L
Entonces la prdida total es = 270 + 21.6775 = 291.6775 psia
Finalmente la presin corriente abajo p2 = 850 291.6775 = 558.3225 psia
5.3.2 Mtodo Beggs y Brill
Datos:
O = 20 API
g = 0.70
d = 12 pg
L = 1 milla = 5280 pies
p1 = 425 psia
p = 100 psig
2.
p = 425 50 = 375 psia.
3. Las propiedades de los fluidos se calcularon con las ecuaciones vistas en el captulo
1.
R s = 51.3871 B o = 1.0279 o = 4.2 cp
g = 0.0113
o = 25.2504
Z g = 0.9413.
141.5
= 0.9340
131.5 + 20
197
Captulo 5
o =
g =
pie 3
(0.0764)(0.70)(375)(520) = 1.3703 lb m
(14.7 )(90 + 460)(0.9413)
pie 3
375
seg
pie 3
seg
pie
seg
0.4763
= 0.3238
0.4763 + 0.9948
198
Captulo 5
N FR =
1.8730 2
= 0.1090
12
32.2
12
L = o = 25.2504
11.
N Rens =
N LV
36.3734 (1)
= 39,578.1539
(1.3676) 6.72 10 4
57.1146
= (1.938)( 0.6065)
25.2504
0.25
= 1.4414
L 1 = (316
= 224 . 73
L 2 = (0.0009252 )(0.3238)
= 0.015
2.4684
L 3 = (0.10)(0.3238)
-1.4516
L 4 = (0.5)(0.3238)
-6.738
= 0.5139
= 997.0676
(0.98)(0.3238)0.4846
(0.1090)0.0868
= 0.6878
(0.845)(0.3238)0.5351
(0.1090)0.0173
= 0.4802
199
Captulo 5
Como el colgamiento se encuentra entre estos dos valores, tomaremos el valor de:
H L (0 ) = 0.58
15. El coeficiente del factor de correccin por inclinacin es:
0.305
-0.4473
(0.1090)0.0978 = 0.2441
C = (1 0.3238)ln (2.96 )(0.3238) (1.4414 )
0.3238
= 0.9202
[0.5932]2
ln (y ) = ln(0.9202 ) = 0.08316
S=
- 0.08316
= 0.2575
2
4
0.0523 + 3.182 (- 0.08316 ) 0.8725 [- 0.08316 ] + 0.01853 [- 0.08316 ]
fT
= eS = e 0.2575 = 1.2937
f ns
19. Calcular el Factor de Friccin sin considerar el resbalamiento.
f ns =
39578.1539
2 log
(4.5223) log (39578.1539) 3.8215
21. Calcular:
200
= 0.0220
Captulo 5
(0.0285)(36.3734)(1.8730)(12)
52801(57.1146)(0.0523) +
(2)(32.174)(12)
= 110.6418 psia
p =
(57.1146)(1.8730)(1.2666 )
(144)1
(32.174)(375)(144)
201
Captulo 5
La gran ventaja sobre los modelos tradicionales es que, cuando esto se consigue, es
posible realizar una experimentacin intensiva, sistemtica y automatizable sobre el modelo
en vez de sobre el sistema natural.
NOTA: Para ste tema de Modelos Mecansticos las variables de las ecuaciones tendrn
que estar en unidades consistentes y su nomenclatura se indicar en su momento.
5.4.1 Modelo para predecir la transicin de los patrones de flujo segn Xiao
h L ( L g ) g A g cos
v g > 1
d
A
g L
h L
0.5
(5.44)
Cuando las ondas son inestables, el flujo puede cambiar ya sea a flujo intermitente o
a anular, dependiendo si existe o no el suficiente suministro de lquido.
202
Captulo 5
(5.45)
Considera que los mecanismos que gobiernan esta transicin son por un lado, el
proceso de turbulencia el cual separa las burbujas de gas evitando su coalescencia, y por
otro, los efectos de la fuerza de tensin superficial actuando para preservar la configuracin
de flujo burbuja. El criterio de transicin es expresado como:
4 A g g cos g
1
vL >
S
f
i
L
L
0.5
(5.46)
Donde:
Si = permetro mojado de la interfase
fL = factor de friccin de Fanning
Transicin Estratificado LisoEstratificado Ondulado
Ce L g v L
0.5
(5.47)
vL
(g h L )0.5
> 1.5
(5.48)
203
Captulo 5
5.4.2 Modelo para predecir la transicin de los patrones de flujo segn Barnea
Barnea (1987), utiliz los siguientes criterios y mecanismos para identificar y
modelar las principales fronteras de transicin entre los patrones de flujo.
Transicin para Burbuja o Burbuja Dispersa
v gS 2 f m 3
D C 0.725 + 4.15
vm
v m L d
(5.49)
= tensin superficial
d = dimetro de la tubera o dimetro hidrulico
La frontera de transicin es vlida para 0 < g < 0.52, ya que en el lmite superior
se alcanza la mxima densidad de empacamiento volumtrico de las burbujas y la
coalescencia ocurre an a altos niveles de turbulencia. Esta condicin puede expresarse
como:
v sL = v gS
Donde:
1 g
(5.50)
D CD
0.4
= 2
g
L
D CB =
204
0.5
2
3 L fm vm
8 ( L g ) g cos
(5.51)
(5.52)
Captulo 5
Donde:
DCD = tamao de burbuja crtica arriba del cual la burbuja es deformada
DCB = tamao de burbuja crtica abajo del cual se previene la migracin de
las burbujas a la parte superior
Transicin Estratificado No Estratificado
Taitel y Dukler consideran flujo estratificado con ondas finitas en la superficie sobre
el cual fluye gas, y como el gas se acelera sobre la cresta de la onda, la presin en la fase
gas disminuye debido al efecto de Bernoulli y la onda tiende a crecer.
Para inclinaciones ascendentes despreciables, el flujo estratificado se encoge
substancialmente y prcticamente desaparece para ngulos de inclinacin de 30o. Para
cambios de inclinacin descendente se tiene un efecto profundo en el rgimen de flujo
estratificado, ya que se expande considerablemente a medida que el ngulo de inclinacin
incrementa hasta ngulos grandes aproximadamente 80o, para desaparecer totalmente a
ngulos cercanos a -90o. Este criterio puede expresarse en trminos adimensionales como:
2 A L
vg
h L
2
1
Fr
2
(1 h L ) A g
(5.53)
0.5
v gS
0.5
(d g cos )
(5.54)
vL >
g d (1 h L )cos
fL
(5.55)
205
Captulo 5
p
2
f
AL
x LS
2
Z=
(1 h L ) LS
L g cos
fL
A
(5.56)
Donde:
fL = factor de friccin de la fase lquida empleando el dimetro del lquido y
la velocidad real del lquido
fLS = factor de friccin considerando que solo fluye lquido
(p / x )LS = gradiente de presin si fluyera solo lquido
Transicin Anular Intermitente
b)
1 + 75 H L
X
3
2.5
(1 H L ) H L H L
(5.57)
3
2 HL
2
Y
X2
3
3
H L 1 H L
(5.58)
Donde:
2
4
f LS L v LS
2
X = d
v 2
4 f gS g gS
2
d
206
0.5
p
x LS
=
p
x
gS
0.5
(5.59)
Captulo 5
Y=
Siendo:
(p / x )LS
(p / x )gS
g ) g sen
p
x gS
(5.60)
(5.61)
Los criterios para flujo estratificado liso y estratificado ondulado son los mismos a
los utilizados en el modelo de Xiao (1990), pero para el rango completo del ngulo de
inclinacin.
Subregiones en Flujo Intermitente
207
Captulo 5
gb
2 f v 3 0.4 0.6
m
m
L
= 1 H b = 0.058d C
0
.
725
d
(5.62)
208
Captulo 5
interfase plana de Taitel y Dukler (1976), para describir el flujo ondulado estratificado en la
horizontal y cercana a la horizontal.
Chen, X. Cai y James P. Brill (1996), desarrollaron para el factor de friccin
interfacial un modelo mecanstico denominado "doble crculo" para predecir el colgamiento
del lquido y la cada de presin durante el flujo estratificado ondulado en tuberas
horizontales.
El modelo de Taitel y Dukler (1976), mostrado en la figura 5.4, es ampliamente
usado para predecir el flujo gas - lquido estratificado ondulado en tuberas horizontales y
cercanas a la horizontal, de acuerdo con Andritsos y Hanratty (1987), Cheremisinoff y
Davis (1979), Kowalski (1987), Oliemans (1987) y Shoham y Taitel (1984). Este modelo
de interfase plana puede desviar ampliamente la configuracin interfacial real y causar una
baja estimacin del rea interfacial, la cual puede ser de importancia.
Hamersma - Hart (1987) y Hart - colaboradores (1989) propusieron el modelo de
"superficie rugosa aparente", mostrado en la figura 5.4, el cual asume que el permetro
mojado de la pared es cubierto por una capa lquida de espesor constante. Esto no es
realista, puesto que la fuerza gravitacional siempre causar que el lquido se acumule
mayormente en el fondo que en cualquier otra parte de la tubera.
= 0.52 H L
0.374
+ 0.26 Fr 0.58
(5.63)
209
Captulo 5
Donde:
vL 2
( L g ) g d
Fr =
(5.64)
Fr =
(5.65)
g ) Di g cos -1 (R )
SL
=
d
(5.66)
H L
i =
2
sen ( ) tan (i )
(5.67)
210
d sen ( )
sen (i )
(5.68)
Captulo 5
(5.69)
v gt =
v
e
g
L
0.5
(5.70)
0.20
sg 1
v
gt
0.08
(5.71)
211
Captulo 5
L = fL
i = fi
g vg 2
2
L vL2
2
g (vg v L )2
2
(5.72)
(5.73)
(5.74)
Los factores de friccin de Fanning para lquido - pared y gas - pared; fL y fg;
respectivamente, se obtienen con:
Para Re < 2000
f=
16
Re
(5.75)
1
9.35
2
= 3.48 4 log
+
f
Re f
d
(5.76)
dL =
4 Ag
sg + si
4 AL
sL
(5.77)
(5.78)
212
dp
- A L + i s i wLs L AL L g sen ( ) = 0
dx
(5.79)
dp
- A g i s i wg s g Ag g g sen ( ) = 0
dx
(5.80)
Captulo 5
s g
sL
wg
AL
Ag
i
+
wg
si
s
+ i + ( L + g ) g sen ( ) = 0
A
L Ag
(5.81)
dp wL s L + wg s g A L
=
+
L +
g g sen ( )
A
A
dx
A
(5.82)
Captulo 5
La gran mayora de los modelos para flujo bache consideran flujo estable, el cual es
ordenado y con baches relativamente cortos (menores a 100d), y un gasto promedio de
lquido y gas constante, sobre un perodo de tiempo de un ciclo, llamado unidad bache,
adems, no se considerar transferencia de calor en el flujo.
En la zona de pelcula, el gas y el lquido se separan por efecto de gravedad como
en el caso de flujo estratificado, variando la altura del lquido, para este modelo se
considera una altura constante a lo largo de toda la zona.
(5.83)
1.0
(5.84)
1.39
1.0 + v s
8.66
Si H S p 0.48 H S = 0.48
(5.85)
g ( L g ) 4
n
v b = 1.2 v S + 1.53
H S sen ( )
2
L
(5.86)
214
(5.87)
Captulo 5
v g = v f H f + v g (1 H f )
(5.88)
vL =
v S v b (1 H S )
HS
(5.89)
Velocidad de traslacin.
(5.90)
(5.91)
(5.92)
215
Captulo 5
Donde:
Nmero de Reynolds:
(5.93)
Re =
vS d
vL
(5.94)
Fr =
vS
gd
(5.95)
Froude :
L g d2
(5.96)
Dumitrescu (1943)
C = 0.35
Tuberas Verticales
Benjamn (1968)
C = 0.54
Tuberas Horizontales
C = 2.27
Flujo Laminar
C = 1.2
Flujo Turbulento
C=
2.27
Re
1+
Rec
1.2
Re
1+
Rec
(5.97)
216
S vS d
S
(5.98)
Captulo 5
(5.99)
S = H S L + (1 H S ) g
(5.100)
(5.101)
Donde
HS, y Hf son los colgamientos de lquido en, el cuerpo del bache y en la zona de la
pelcula, respectivamente.
De un balance de masa aplicado a dos secciones de rea transversal con relacin a
un sistema coordenado movindose a la velocidad de traslacin V t , para la fase liquida
resulta:
(v t v L ) H S = (v t v f ) H f
(5.102)
Despejando Vf:
vf = vt
(v t v L ) H S
Ef
(5.103)
v t H S + v b (1 H S ) v SG
vt
(5.104)
sg s
sf
s
g
+ i i + i + ( L g ) sen ( ) = 0
Af
A g g A f A g
(5.105)
Donde:
217
Captulo 5
f = ff
g = fg
i = fi
L vf vf
(5.106)
g vg vg
(5.107)
g v g v f (v g - v f )
(5.108)
2
Re f =
Re G =
L vf d L
(5.109)
G vG dG
(5.110)
s vs 2
(5.111)
S = H S L + (1 + H S ) g
(5.112)
S = H S L + (1 + H S ) g
(5.113)
0.1
(5.114)
Para tuberas de dimetro menor a 0.0381 mm. (1.5 pg.) Esta longitud puede
calcularse de manera aproximada por:
218
Captulo 5
LS = (32)d
(5.115)
L u =LS
vL vf Hf
v SL v f H f
(5.116)
L u = LS + L f
(5.117)
L f = L u LS
(5.118)
Longitud de la pelcula.
Del balance de fuerza para una unidad de bache se llega a la expresin para calcular
el gradiente de presin promedio para flujo intermitente.
1
dp
= u g sen ( ) +
Lu
dx
i sf + g sg
L f
S
L S +
A
(5.119)
u = H m L + (1 H m ) g
(5.120)
En flujo anular la fase lquida fluye como una pelcula sobre la pared de la tubera y
como gotas de lquido atrapadas en el ncleo gaseoso. El flujo anular se presenta a altos
gastos de gas y de bajas a medios gastos de lquido, como se muestra en la Figura 5.7.
219
Captulo 5
(5.122)
c = H c L + (1 H c ) g
(5.123)
v sL FE
v sg + v sL FE
(5.124)
Para el clculo del factor de arrastre de lquido (FE), se encuentra las correlaciones
de Wallis (1969). Whalley y Hewit (1978). Y Olierrians y Col. (1986). Xiao (1990), expone
220
Captulo 5
que el factor de arrastre propuesto por Oliemaris y Col. (1986) da mejores resultados,
donde la expresin para el factor de arrastre de lquido es:
FE
2
4
7
8
= 100 1 g 3 g 5 d 6 v SL v SG g 9
1 FE
(5.125)
Donde los parmetros beta () con los coeficientes de regresin obtenidos por
Oliemans y Col. (1986), empleando la base de datos de Harwell. Con los siguientes valores:
0 = 2.52 1 = 1.08 2 = 0.18 3 = 0.27 4 = 0.28
5 = 1.80 6 = 1.72 7 = 0.70 8 = 1.44 9 = 0.46
Suponiendo condiciones de equilibrio entre la fase lquida y el ncleo de gas y
combinando las dos ecuaciones anteriores se llega a:
wL
1
sL
1
+ ( L c ) g sen ( ) = 0
i s i
+
Af
Af Ac
(5.126)
H L = 1 1 2
d v sg + v sL FE
(5.127)
AC =
(5.128)
A F = (d 2 )
(5.129)
SL = d Si = (d 2 )
(5.130)
4 (d )
d hC = 2 (d )
d
(5.131)
221
Captulo 5
4 (d )
(5.132)
(v
sg
+ v sL FE ) d 2
(d 2 )2
(5.133)
L vf 2
(5.134)
Donde:
d v
f f = C f hF F
L
(5.135)
Esta ecuacin es una forma de la ecuacin de Blasius, los coeficientes usados son:
CF
Flujo
16
Laminar
0.046
0.2
Turbulento
(5.136)
222
(5.137)
Captulo 5
(5.138)
f i = f c 1 + 2250
c (v c v f )2
Donde el factor de friccin interfacial del ncleo fc, puede calcularse con las
ecuaciones (5.75) (5.76) empleando la siguiente definicin del nmero de Reynolds.
Rec =
Donde
c vc dc
(5.139)
c = E c L + (1 E c ) G
(5.140)
d c = d 2
(5.141)
dp s
= wL L + vL f + c c g sen ( )
A
A
A
dx
(5.142)
Velocidad de la mezcla.
La velocidad de la mezcla se define mediante la ecuacin (5.18); utilizando el
modelo de flujo deriva para la velocidad de las burbujas dispersas:
vG =
v sg
(5.143)
1 HL
g ( L g ) 4
0.1
v G 1.2 v S 1.53
E L sen ( ) = 0
2
(5.144)
223
Captulo 5
(5.145)
m = H L L + (1 H L ) g
(5.146)
m vs d
m
(5.147)
Gradiente de presin. El gradiente de presin puede ser llevado a cabo como flujo
de una sola fase con propiedades promedio de la mezcla:
2
dp 2 f m m v s
=
+ m g sen ( )
d
dx
224
(5.148)
Captulo 6
Flujo en estranguladores
225
Captulo 6
Flujo en estranguladores
vf
= M = 1 10 5
vp
(6.1)
p2
= 0.528
p1
(6.2)
Mach 1 o flujo snico para gases (1100 pie/seg para aire) es diferente que para los
lquidos.
Para flujo multifsico algunos investigadores han observado que la velocidad snica
para una mezcla de gas-lquido es menor que para una sola fase. Olson (Essentials of
Engineering Fluid Mechanics, Intl. Textbook, 1961) estableci que para mezclas gaslquido la velocidad acstica es menor que para cualquiera de las fases por s solas.
Durante la presencia de flujo crtico o supercrtico, el gasto a travs del
estrangulador alcanza un valor mximo con respecto a las condiciones prevalecientes
corriente arriba. La velocidad del fluido bifsico fluyendo a travs de la restriccin alcanza
la presin snica o presin de velocidad de propagacin de la onda para los fluidos en dos
fases. Esto implica que el flujo es estrangulado porque los disturbios corriente abajo no
pueden propagarse corriente arriba. Por lo tanto disminuciones no tan significantes en la
presin corriente abajo no hacen incrementar el gasto, es decir, el flujo crtico o snico es el
flujo en el cual perturbaciones de presin y temperatura corriente abajo no son transmitidas
corriente arriba tal que puedan afectar el gasto, no as en flujo subcrtico. Si la presin
corriente abajo es gradualmente incrementada, esta presin no podra cambiar el gasto o la
presin corriente arriba, esto hasta alcanzar el lmite de flujo crtico-subcrtico. A partir de
este momento si la presin corriente abajo se incrementa ligeramente cerca de las
condiciones lmite, el gasto y la presin corriente arriba sern afectados, entonces la
velocidad de los fluidos pasando a travs del estrangulador caen debajo de la velocidad
snica. Aqu, el gasto depende de la diferencial de presin, o bien, los cambios en la
presin corriente abajo afectan la presin corriente arriba, este comportamiento es
caracterizado como flujo subcrtico.
Los estranguladores superficiales se disean de manera que el flujo sea crtico,
mientras que en las vlvulas de seguridad subsuperficiales el flujo es subcrtico. El flujo a
travs de restricciones en el pozo tambin es subcrtico.
226
Captulo 6
Flujo en estranguladores
CD =
real:
gasto medido
gasto calculado
(6.3)
(6.4)
Captulo 6
Flujo en estranguladores
todos los procesos de flujo son tomados en cuenta de manera precisa. Las leyes de la
termodinmica implican que pueden resultar valores de CD menores que la unidad.
Frecuentemente, CD tambin es correlacionado con factores como el nmero de Reynolds,
diferencial de presin, factor de expansin del gas, etc. Obviamente, si un modelo toma en
cuenta variables como stas, tal dependencia no debe de existir.
En el lenguaje comn se ha adoptado decir flujo crtico a lo que estrictamente es
flujo supercrtico, cabe aclarar que el flujo crtico es la condicin de flujo que se alcanza
cuando la velocidad del fluido es igual a la velocidad de propagacin de una onda en la
mezcla del fluido en cuestin. En la prctica difcilmente se puede mantener estable el valor
de la velocidad y se llega a sobrepasar por lo que en realidad es un flujo supercrtico el cual
no tiene mayor diferencia en trminos significativos. En lo posterior se mencionar flujo
crtico para hacer referencia al flujo supercrtico, dejando en claro que el flujo crtico es el
lmite entre flujo supercrtico y flujo subcrtico.
q c.s.
C p d2
= n 1
g T1 Z1
2
k +1
k p 2 k p 2 k
k 1 p1 p1
(6.5)
C s C d Tcs
p cs
(6.6)
Donde:
Cn =
228
Captulo 6
Flujo en estranguladores
Las tablas siguientes muestran valores para las constantes en la ecuacin para varios
sistemas de unidades.
Tabla 6.1. Coeficientes y Unidades para la Ecuacin 6.5.
Sistema
Sistema
Sistema
Smbolo
Mtrico
Ingls
Mtrico
Internacional
2
Mpies@cs/d
m /d
m2/d
qsc
pg
mm
mm
d
2
2
lb/pg abs
kg/cm
kPa
pabs
R
K
K
Tabs
27.611
1.6259
1.6259
Cs
Tabla 6.2. Coeficiente para la Ecuacin 6.5
Sistema de
Unidades
Ingls
Mtrico
Mtrico
Internacional
Cd
p @c.s.
T @c.s.
Cn
0.865
0.865
0.865
0.865
0.865
0.865
14.696 lb/pg2
14.696 lb/pg2
1.0332 kg/cm2
1.0332 kg/cm2
101.325 kPa
101.325 kPa
491.68 R
519.68 R
273.16 K
288.72 K
273.16 K
288.72 K
799.06
344.57
371.83
393.01
3.7915
4.0075
229
Captulo 6
Flujo en estranguladores
k=
Cp
Cv
M Cp
M C p 1.987
(6.7)
Donde:
M = Peso molecular, lbm/mol.
Cp = Calor especfico, BTU/lbm-R
El valor de k tambin se puede obtener de la figura 6.1:
R pc
2 k 1
=
k +1
(6.8)
230
Captulo 6
Flujo en estranguladores
antes del estrangulador, as como el dimetro, esto hace que las soluciones para flujo crtico
sean ms sencillas hasta cierto punto y por lo tanto ms abundantes en comparacin a
aquellas para flujo subcrtico.
Las soluciones para flujo subcrtico requieren para su aplicacin gran cantidad de
informacin referente a los fluidos en trnsito y por otra parte involucran procesos
iterativos, que traen consigo problemas de convergencia.
p1 =
A qL R B
d C
(6.9)
231
Captulo 6
Flujo en estranguladores
Tabla 6.3. Coeficientes empricos para correlaciones de flujo crtico en dos fases.
Correlacin
A
B
C
Gilbert
10
0.546
1.89
Ros
17.4
0.5
2.0
Baxendell
9.56
0.546
1.93
Achong
3.82
0.650
1.88
Este modelo fue establecido a partir del trabajo presentado por Ros en 1960, quien
se bas en el anlisis de la ecuacin de balance de energa, desarroll una frmula de
medidor de flujo a partir de un anlisis terico del flujo simultneo gas-lquido a velocidad
crtica (snica) a travs de restricciones (orificios) y una correlacin para el
comportamiento PVT de los fluidos. No se consider produccin de agua. La precisin de
los resultados obtenidos se comprob comparndolos con 108 datos medidos.
Para que exista flujo crtico se supuso que la presin corriente abajo, debe ser al
menos de 0.55 de la presin en la boca del pozo. Bajo estas condiciones el gasto en el
estrangulador es slo funcin de la presin corriente arriba y de la relacin gas-aceite a
condiciones de flujo.
Despus en 1963 Poettman y Beck convirtieron la ecuacin a unidades de campo y
la redujeron a una forma grfica, lo que result en la siguiente expresin:
9273.6 p1
86,400 C D A o
qL =
5.6142 LS + 0.0764 g R VL (1 + 0.5m L )
0.5
0.4513 r + 0.766
r + 0.5663
(6.10)
Donde:
r=
0.00504 T1 Z1 (R R s )
p1 Bo1
mL =
1
g1
1 + r
L1
VL =
mL
o1
(6.11)
(6.12)
(6.13)
Siendo:
qL = Gasto en barriles @ c.s./da.
CD = Coeficiente de descarga (1.03).
Ao = rea transversal de la garganta en pulgadas cuadradas (la garganta es el rea
transversal mnima de un estrangulador).
232
Captulo 6
Flujo en estranguladores
233
Captulo 6
Flujo en estranguladores
234
Captulo 6
Flujo en estranguladores
235
Captulo 6
Flujo en estranguladores
En 1968 Omaa desarroll una correlacin (para flujo crtico) entre el gasto, la
presin 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 de campo tomados en las instalaciones de la Union Oil Company del campo
Tigre Lagon de California en Louisiana para revisar las correlaciones existentes y
desarrollar la suya. Los experimentos de campo se realizaron con agua y gas natural. En
vista de que estos datos estuvieron dentro de rangos muy limitados, esta correlacin no es
aceptada ampliamente debido a:
1. Limitaciones en el tamao del estrangulador (4, 6, 8, 10 12 y 14/64 de pg).
2. Limitaciones en el gasto (0 a 800 bl/da de agua).
236
Captulo 6
Flujo en estranguladores
qL =
N qL
1.84 L
L
(6.14)
1.25
Donde:
N q = 0.263N 3.49 (N p )
3.19
N=
g
L
1
N p = p1
L g
Q=
1+
g
Nd = d L
0.5
(Q )0.657 (N d )1.8
(6.15)
(6.16)
0. 5
0.0174 p1
( L L )0.5
(6.17)
1
(R Rs )Bg
(6.18)
5.615 Bo
= 120.872 d L
0.5
(6.19)
237
Captulo 6
Flujo en estranguladores
0.5
qo =
(6.20)
donde:
d = Dimetro del estrangulador (64avos de pg).
P1 = presin corriente arriba, psia
6.3.1.5 Velocidades Snicas de Nguyen
Nguyen estudi la velocidad snica en sistemas de dos fases como funcin del
patrn de flujo.
Para flujo estratificado, una velocidad snica combinada no existe porque cada fase
es continua en la direccin axial. Una velocidad snica efectiva existe en cada fase que est
influenciada por la otra fase. Si se somete el lquido y el gas a un cambio de presin de
manera simultnea, la perturbacin se propaga con velocidades diferentes en ambas fases
en la direccin axial.
La ecuacin 6.21 es la velocidad snica efectiva para la fase gaseosa y muestra que
la velocidad snica efectiva est gobernada primariamente por la velocidad snica del gas
porque el segundo trmino en el denominador es pequeo.
vE*2, g =
1
1 H L g 1
+
v*g2 1 H L L vL*2
238
(6.21)
Captulo 6
Flujo en estranguladores
vE*2, L =
1
1 1 H L L 1
+
vL*2 H L g v*g2
(6.22)
vL* v*g
H L v*g + (1 H L )vL*
(6.23)
Para flujo homogneo, Nguyen combin expresiones para las velocidades snicas
de cada fase fluyendo dentro de una frontera elstica con el concepto de que el frente de la
onda pasa de manera secuencial a travs de zonas de lquido y gas dentro de la mezcla
homognea. La ecuacin 6.24 da la expresin resultante.
v* =
(1 )
1 L
(1 ) g
*2 +
+ *2 +
*
2
v
g vg
vg
LvL*2
L
(6.24)
En la ecuacin 6.25 Wallis present una expresin para calcular la velocidad snica
o de onda de compresibilidad de una mezcla homognea.
vm* = ( g g + L L ) g*2 + L*2
v
LvL
g g
0.5
(6.25)
Ashford y Pierce desarrollaron una expresin para el gasto msico total de una
mezcla multifsica. Ellos asumieron flujo isoentrpico a travs de la restriccin, lquido
239
Captulo 6
Flujo en estranguladores
k 1 + R1 yc
yc =
R1
e
b
1 yc yc + 1 yc 1
R1
(6.26)
Donde:
b=
k 1
k
(6.27)
e=
k +1
k
(6.28)
N k 1
yc =
D
(6.29)
Donde:
N=
Y:
240
(1 x g1 ) g1 (1 yc )
k
+
k 1
x g 1 L
(6.30)
Captulo 6
Flujo en estranguladores
D=
k
n n(1 x g1 ) g 2 n (1 x g1 ) g 2
+
+ +
2 x g1 L
k 1 2
x g1 L
(6.31)
x g1 (C pg Cvg )
x g1 C vg + (1 x g1 )C L
(6.32)
Y:
xg 1 =
wg1
wg1 + wL1
(6.33)
Donde:
wg1 = 0.0764 g q Lsc (R p f o Rs1 )
(6.34)
(6.35)
6.3.2.1 Fortunati
Fortunati present un mtodo emprico que puede ser utilizado para calcular el flujo
multifsico crtico y subcrtico a travs de estranguladores. El asumi una mezcla
homognea y aclar que esa suposicin era vlida siempre y cuando vm sea mayor que 32.8
pie/seg y que el nmero Froude de la mezcla sea mayor a 600, adems asumi que no haba
resbalamiento entre las fases, aunque reconoci que el resbalamiento existe incluso para
lquidos inmiscibles.
Utilizando datos experimentales, Fortunati desarroll la figura 6.6, la cual puede ser
utilizada para definir la frontera entre el flujo crtico y el subcrtico. Las curvas de la figura
241
Captulo 6
Flujo en estranguladores
6.6 se basaron en una presin de 0.137 MN/m2 (19.8 psia) corriente abajo. Se utiliza la
ecuacin 6.36 para calcular la velocidad Fortunati de la mezcla a partir de la velocidad
actual de la mezcla y de la presin corriente abajo, en la cual sus variables deben estar en
unidades congruentes.
vm 2
Donde:
p2
= v
p2 F
= (1 3 )
0.38
Y:
242
(6.36)
(6.37)
Captulo 6
Flujo en estranguladores
o w
v2 =
2
A (x g + L1 )
1 2
A1 x y n +
g
c
L1
0.5
(6.38)
243
Captulo 6
288 g c p1
Wi = A2
v1
Flujo en estranguladores
n 1
1 y c n + L1 (1 y c )
1
(x g + L1 )
A2
1
x
y
g c n
1
A1 x g y c n + L1
2
+ L1
(6.39)
n 1
n
1
y
c
+ 2 L1 (1 yc )
2
(1+ n )
2
2 x +
A2
x g1 (1+ n ) A2 x g1 (x g1 + L1 )yc n
g1
L1
n
y
+
c
1
2
A1 n x y 1n +
A1 x g1 yc n + L1 n
g1 c
L1
1
A2 x g1 + L1
= 1
x g1 yc n + L1
1 x g1 yc n + L1
(6.40)
n 1 1n
yc + L1
Donde:
xo1
L1 = n1
o1
xw1
w1
(6.41)
= x g + (x g C vg + xo C vo + x w Cvw )
Z R
(6.42)
244
x g k C vg + xo C vo + x w C vw
x g C vg + xo C vo + x w C vw
(6.43)
Captulo 6
Flujo en estranguladores
( p1 p4 )
d
1 c
d d
1.85
(6.44)
El mtodo para distinguir entre flujo crtico y subcrtico y como manejar flujo
crtico y subcrtico es como sigue:
1. Calule p2 (la presin en la garganta del estrangulador cuando el flujo es
crtico), resolviendo la ecuacin 6.40.
2. Calcule p3 (la presin de descarga justo corriente abajo del estrangulador),
resolviendo la ecuacin 6.44.
3. Si p2 > p3, entonces el flujo es crtico y yc = p2/p1, se deber usar las
ecuaciones 6.38 y 6.39 para calcular la velocidad isoentrpica y gasto
msico.
4. Si p2 = p3, entonces el flujo est en el lmite entre flujo crtico y subcrtico.
(Use yc = p2/p1 = p3/p1 para calcular la velocidad isoentrpica y gasto
msico).
5. Si p2 < p3, entonces el flujo es subcrtico. (Use yc = p3/p1 en las ecuaciones
6.38 y 6.39 para calcular velocidad isoentrpica y gasto msico).
En todos los casos, el gasto real es determinado por multiplicar el gasto isoentrpico
por el coeficiente de descarga CD.
6.3.2.3 Modelo de Ashford y Pierce
Ashford y Pierce establecieron una ecuacin que describe la dinmica de las cadas
de presin y capacidades de flujo en condiciones de flujo multifsico. Este modelo
relaciona el comportamiento del estrangulador en ambos regmenes de flujo, crtico y no
crtico.
La capacidad y cadas de presin que se presentan en la restriccin se han
relacionado con sus dimensiones y las propiedades de los fluidos manejados.
Los datos usados, reflejan el comportamiento de una vlvula de seguridad OTIS tipo
J-22J037. Sin embargo, el modelo puede usarse para estimar las cadas de presin a travs
de cualquier dispositivo que restrinja el flujo.
Para la validacin del modelo, se dise una prueba de campo en un pozo fluyente.
Tanto las cadas de presin como el gasto se midieron directamente y luego se compararon
con datos anlogos obtenidos del modelo. Esta informacin se us para determinar el
245
Captulo 6
Flujo en estranguladores
coeficiente de descarga del orificio, definido por la relacin de gasto medido entre el gasto
calculado.
La ecuacin obtenida por los autores es:
qo = 1.9706 C d 2
(6.45)
donde:
1
= (Bo + WOR ) 2
(6.46)
k
0.5
k
T1Z1
k
(
)
198
.
6
+
R
R
E
o + 0.000217 g Rs + WOR w
s
p
1
k=
E=
Cp
]
(6.47)
(6.48)
Cv
p1
p2
(6.49)
m = k 1
(6.50)
559
615
0.9089
16
484
402
1.2039
16
334
C=
0.9089 + 1.2039
= 1.0564
2
246
Captulo 6
Flujo en estranguladores
1.1511
16
1.0564
20
0.9760
C = 2.398 0.477 ln (d )
(6.51)
p TP = p L L + p g g
qL
p L =
2 g c 144 C L A c
p g =
2 g c 144 Z C g A c
(6.52)
2
(6.53)
d 2 1 p 2 p1
d1 k p1
C=
CD
d1
1
d2
(6.54)
(6.55)
(6.56)
247
Captulo 6
Flujo en estranguladores
C
DL
= p L 1 + g
Z DDg
(6.57)
Donde:
d 4
qm
p L = L 1 + 1
2
d 2 8083 d1 C DL
Y:
(6.58)
6.5 Ejemplos
Ejemplo 1
Dados los siguientes datos mostrados abajo, calcula el tamao de estrangulador en
64avos de pulgada, usando las correlaciones de Gilbert, Ros, Baxandell, Achong, Poettman
y Beck, Omana y Ashford para flujo critico.
qo=
qg=
RGL =
dt =
Ap =
CD =
API =
g =
T=
P1 =
L1 =
Rs1 =
Bo1 =
Z1 =
g1 =
g1 =
248
1800
1.116
620
2.994
0.0487
1
30
0.7
128
1663
9
335
1.10
0.78
6.856
52.59
blp/d
Mmpc/d
Pg
Pie2
F
Psia
Dinas/cm
Lbm/pie3
Lbm/pie3
Captulo 6
Flujo en estranguladores
A qL R BP C
d =
P1
Correlacin
d
1
1663
Gilbert
1663
Ros
1663
Baxendell
1663
Achong
v sgl =
q o (R p R s1 ) B gl
Ap 86400
v sL1 =
lb m
bl
128
= 0.0077
1663
R=
v sgl
v sLl
0.9415
= 0.3564
2.6418
249
Captulo 6
Flujo en estranguladores
g s.c.
gl =
B gl
mL =
1+ R
gl
q L m
86400 C D
(1800)(341.23)
(86400)(1)
lb
0.7 0.765
= 6.9346 m3
0.0077
pie
L1
vL =
Ao =
1
6.856
1 + 0.3564
52.59
= 0.9556
m L 0.7865
=
= 0.0182
L1
52.59
1
=
0.4513 R + 0.766
9273.6 p1
v L (1 + 0.5 m L )
R + 0.5663
= 5.7247 10 4 ft 2
5.7247 10 4 =
d2
4
Omaa
R=
N gv
N Lv
QD =
250
20.7344
pg
64
v sg1
v sL1
0.9415
= 0.3564
2.6418
1
1
=
= 0.7373
1 + R 1 + 0.3564
Captulo 6
Flujo en estranguladores
Np = 1.74 10 p1
L1 L1
0.5
= 1.74 10 1663
52.59 9
0.5
N qL
= 1.84 q L L1
L1
N d = 120.872 d O
N qL = (0.263) N
3.49
)(
= 1.3301
1.25
= 30089.5820
L1
= 290.0759
L1
) (N ) (Q
3.19
0.657
D
) (N )
1.80
)(
)(
d 1.80
=
O
1.80
1.80
30089.5820
= 0.0017
1.773 10 7
d O = 0.0289 ft
d O = 0.0289 12 = 0.3468 pg
dO =
22.20
pg
64
Ashford
A = 5.04 10 3 T1 z 1 (R P R s1 ) = 658.7905
B = 62.4
145.7
= 62.4 0.8793 = 54.8683
135.7 + 30
1
1
=
= 0.9535
B O1
1.1
251
Captulo 6
Flujo en estranguladores
dc =
qO + 0.56
P1
A + 0.76 p1
(B + 0.01353 g R p )2 + (B + 0.01353 g R s1 )1
658.7905
+ 0.56
1800
1663
= 18.5899
=
252
18.59
pg
64
CONCLUSIONES Y RECOMENDACIONES
A pesar de que ha aumentado el uso de nuevas tcnicas para la prediccin y
comportamiento de los fluidos dentro de las tuberas, no deben dejar de utilizarse los
mtodos analticos existentes, como las correlaciones. Esto se debe a que son ms prcticos
para ciertas aplicaciones especficas.
Se debe tomar mucho en cuenta el flujo maultifasco en tuberas pues directamente
su aplicacin es para el diseo de las instalaciones y tuberas que se utilizan en la industria
petrolera, ya que en el flujo multifsico siempre se tiene presente en la explotacin de
hidrocarburos, por lo que es importante determinar adecuadamente, mediante las
propiedades de los fluidos, las cadas de presin en las tuberas.
Se debe tener un continuo seguimiento a los textos utilizados para la materia as
como su actualizacin continua y bsqueda de nuevos mtodos para el estudio de flujo
multfasico en tuberas.
253
NOMENCLATURA
NOTA: sta nomenclatura es vlida para toda la tesis a excepcin de las ecuaciones en que
se indique su propia nomenclatura.
Smbolo
Descripcin
Unidades
Ag
(pg2)
AL
(pg2)
Ap
(pg2)
Bg
Bo
Bob
Bw
Bwb
CD
Adimensional
Co
(pg2/lb)
Cw
(pg2/lb)
C N2
C CO2
C H2S
Dimetro
(pulgadas)
dci
(pulgadas)
dci
(pulgadas)
dB
dh
Dimetro hidrulico
(pulgadas)
(64avos de pg)
(psi/pie)
(psi/pie)
(psi/pie)
Fraccin
255
Adimensional
fo
Adimensional
Adimensional
fw
Adimensional
Gg
(lbm/seg-pie2)
GL
(lbm/seg-pie2)
(lbm/seg-pie2)
(pie/seg2)
gc
(lbm-pie/lbf-seg2)
Profundidad
(pies)
Diferencia de altura
(pies)
Suma de elevaciones
(pies)
HL
Colgamiento de lquido
Adimensional
HF
Adimensional
Longitud
(pies)
Ma
Mg
Mo
Mw
Masa
(lbm)
Nd
Nmero de dimetro
Adimensional
NFR
Nmero de Fraude
Adimensional
Ngv
Adimensional
NL
Adimensional
NLv
Adimensional
NRe
Nmero de Reynolds
Adimensional
NReb
Adimensional
NReg
Adimensional
NReL
Adimensional
NRens
Adimensional
256
(lb-mole)
Presin
(lbf/pg2 o psia)
Presin media
(lbf/pg2 o psia)
patm
Presin atmosfrica
(lbf/pg2 o psia)
pb
(lbf/pg2 o psia)
ppc
Presin Pseudocrtica
(lbf/pg2 o psia)
p' pc
(lbf/pg2 o psia)
ppr
Presin Pseudoreducida
Adimesional
ps
(lbf/pg2 o psia)
pth
(lbf/pg2 o psia)
pwf
(lbf/pg2 o psia)
pws
Presin esttica
(lbf/pg2 o psia)
p1
(lbf/pg2 o psia)
p2
(lbf/pg2 o psia)
Cada de Presin
(lbf/pg2 o psia)
pa
(lbf/pg2 o psia)
pe
(lbf/pg2 o psia)
pf
(lbf/pg2 o psia)
pL
(lbf/pg2 o psia)
pT
(lbf/pg2 o psia)
(psi/pie)
(psi/pie)
(psi/pie)
(psi/pie)
qg
(pies3g a c.s./da)
qg
qL
qL
(blL a c.s./da)
qo
(blo a c.s./da)
qT
(bl a c.s./da)
257
qw
(blw a c.s./da)
RGL
Relacin gas-lquido
Rs
Relacin de solubilidad
Rsw
rh
Radio hidrulico
(pulgadas)
Temperatura
(oF)
Temperatura media
(oF)
T.P.
Tubera de Produccin
T.R.
Tubera de Revestimiento
Tatm
Temperatura atmosfrica
(oF)
Tpc
Temperatura pseudocrtica
(oR)
T' pc
(oR)
Tpr
Temperatura pseudoreducida
Adimensional
Ts
(oF)
Volumen
(pies3)
Vmg
Vmo
Vmw
v, v
Velocidad de flujo
(pies/seg)
v, v
(pies/seg)
vb
(pies/seg)
vf
(pies/seg)
vg
(pies/seg)
258
vL
(pies/seg)
vm
Velocidad de la mezcla
(pies/seg)
vp
vsg
(pies/seg)
vsL
(pies/seg)
wg
(lbmg/seg)
wgd
(lbmgd/seg)
wgf
(lbmgf/seg)
(lbm/seg)
wo
(lbmo/seg)
ww
(lbmw/seg)
yg
y CO 2
y H 2S
y N2
Adimensional
Fraccin de vaco
Adimensional
Adimensional
gd
Adimensional
gf
Adimensional
gp
Adimensional
gs
Adimensional
Adimensional
Adimensional
Adimensional
wg
Adimensional
Adimensional
(pg)
(oR)
259
Adimensional
(cp)
(cp)
(cp)
ns
(cp)
(cp)
oD
(cp)
ob
(cp)
om
(cp)
(cp)
(lbm/pie3)
(lbm/pie3)
(lbm/pie3)
(lbm/pie3)
ns
(lbm/pie3)
(lbm/pie3)
Densidad pseudoreducida
(lbm/pie3)
(lbm/pie3)
Tensin superficial
(dinas/cm)
(dinas/cm)
(dinas/cm)
(dinas/cm)
fr
Adimensional
st
Adimensional
Adimensional
260
BIBLIOGRAFA
1. Brown Kermit, E., et.al., The Technology of Artificial Lift Methods, Vol. 4. Tulsa
Oklahoma. Pennwell Books.
2. Brown Kermit, E., et.al., The Technology of Artificial Lift Methods, Vol. 1. Tulsa
Oklahoma. Pennwell Books, 1977.
3. McCain Jr. William D, The Properties of Petroleum Fluids, Second Edition
PennWell Publishing Company Tulsa, Oklahoma.
4. Garaicochea Petrirena, Francisco; Bernal
Hidrocarburos por Ductos, Mxico: CIPM.
Huicochea;
Csar,
Transporte
de
261
15. Tesis de Maestra de Luis Carlos Leal Lemus, Modelo Mecanstico del Perfil de
Presiones para Flujo Multifsico en Tuberas Horizontales e Inclinadas. Mxico:
UNAM, Divisin de Estudios de Posgrado, Facultad de Ingeniera.
16. James P. Brill, Hemanta Mukherjee, Multiphase in Flow Wells, Monograph Volume
17 SPE, Richardson, Texas 1999.
17. Rey Meneses Hegel Anlisis de modelos de flujo multifsico a travs de
estranguladores. Mxico:. UNAM. 2004
18. T. E. W. Nind, Principles of Oil Well Production. McGraw Hill, USA, 1981.
19. H. Dale Beggs, Gas Production Operations. Oil & Gas, Consultants International
Publications, Inc. Tulsa, Oklahoma, 1991.
20. Ken Arnold, Maurice Stewart, Surface Production Operations. Vol 2. Design of GasHandling Systems and Facilities. Gulf Publishing Company, Houston, Texas, 1989.
21. George V. Chilingarian, John O. Robertson, Sanjay Kumar, Surface Operations in
Petroleum Production I. Elsevier Science Publishers B. V., Netherlands, 1987.
22. George V. Chilingarian, John O. Robertson, Sanjay Kumar, Surface Operations in
Petroleum Production II. Elsevier Science Publishers B. V., Netherlands, 1989.
23. R. Sachdeva, Z. Schmidt, J. P. Brill, R. M. Blais, Two-Phase Flow Through Chokes
Paper SPE 15657.
24. H. Elshahawi, M. Hashem, P. Hegeman, O. C. Mullins, G. Fujisawa, S. Betancourt
Insitu Characterization of Formation Fluid Samples Case Studies. Paper SPE
90932.
25. Katz D. L., Handbook of Natural Gas Engineering. McGraw Hill Book Co., 1959
262