Flujo Multifasico en Tuberias Verticales

Capítulo 3 Flujo multifásico en tuberías verticales
CAPÍTULO 3:
FLUJO MULTIFÁSICO EN TUBERÍAS VERTICALES
3.1 Introducción
El problema de predecir con precisión las caídas de presión en la tubería de
producción (T.P.) se ha ido incrementando a muchas soluciones especializadas para
condiciones limitadas. La razón para estas muchas soluciones es que el flujo multifásico es
complejo y se dificulta su análisis incluso para estas correlaciones estudiadas de
condiciones limitadas.
Al pasar los fluidos provenientes del yacimiento a través de la tubería de
producción, se consume la mayor parte de presión disponible para llevarlos del yacimiento
a las baterías de separación, por lo que es de suma importancia realizar una evaluación
precisa de la distribución de la presión a lo largo de dicha tubería. Al hacerlo
conjuntamente con un análisis integral del sistema de producción, es posible:
a. Diseñar las tuberías de producción y líneas de descarga.

b. Proyectar aparejos de producción artificial (neumático, mecánico, etc).
c. Obtener presión de fondo fluyendo, sin intervenir los pozos.
d. Calcular el efecto de los estranguladores sobre el gasto.
e. Determinar la vida fluyente de los pozos.
f. Corroborar los datos obtenidos con las correlaciones para su ajuste.
Cuando el flujo es vertical las caídas de presión por aceleración son muy pequeñas,
por lo que el gradiente de presión debido a la misma generalmente se desprecia, quedando
la ecuación 2.11 de la siguiente forma:
 Δp   Δp   Δp  
      (3.1)
 Δh T  Δh  e  Δh f
Las diversas correlaciones existentes para el cálculo de distribuciones de presión

con flujo multifásico en tuberías, pueden clasificarse en tres grupos en base al criterio
utilizado en su desarrollo. En la tabla 3.1 se muestran los tres grupos y los criterios
desarrollados.
Hay muchas correlaciones que dan excelentes resultados dependiendo de las
diferentes condiciones de flujo. Basados en las comparaciones hechas por Lawson y Brill,
los métodos vistos en este capítulo están entre las mejores correlaciones para flujo vertical,
ya que pusieron especial cuidado en varios tamaños de tuberías, en las propiedades de los
fluidos y los gastos de flujo. Estas correlaciones son las siguientes:
 Poettmann y Carpenter (1952)

 Duns y Ros (1963)
 Hagedorn y Brown (1965)
83
 Orkiszewski (1967)
 Beggs y Brill (1973)
Tabla 3.1. Criterios en el desarrollo de las correlaciones.
Grupo I Grupo II Grupo III
Resbalamiento
No se considera Si se considera Si se considera
entre las fases
Se calcula en Se calcula Se calcula
Densidad de la función de las considerando el considerando el
mezcla ( m) propiedades de los efecto del efecto del
fluidos colgamiento colgamiento
Se correlaciona con Se correlaciona con
Factor de fricción Se determina las propiedades las propiedades del
(f) empíricamente combinadas del gas fluido en la fase
y el líquido continua
Patrones de flujo No se distinguen No se distinguen Sí se distinguen
Duns y Ros,
Poettmann y
Orkiszewski, Aziz,
Métodos o Carpenter, Fancher y
Hagedorn y Brown Beggs y Brill,
Correlaciones Brown y Baxendell
Chierici, Gould y
y Thomas
Tek entre otras
3.2 Correlaciones
3.2.1 Poettman y Carpenter
Poettmann y Carpenter desarrollaron un método semi empírico utilizando diámetros
de tubería de producción de 2, 2 ½, 3 pulgadas en el que incorporaron la ecuación general
de energía y datos de 34 pozos fluyentes de aceite y 15 con bombeo neumático. El aceite,
gas y agua fueron considerados para estar en una sola fase, y no se trato de determinar la
correlación de colgamiento. Toda la pérdida de energía incluyendo el resbalamiento y
líquidos de emergencia se introdujeron en un factor de perdida de energía la cual se tomó
como constante en toda la longitud de la tubería. También se asumió el flujo vertical de
aceite, gas y agua como turbulento, además de no tomar en cuenta los patrones de flujo.
La correlación de Poettmann y Carpenter sigue siendo usada extensamente en los
campos bajo las siguientes condiciones:
 Tuberías de 2, 2 ½ y 3 pg.  RGL menores a 1,500 pies3/bl.

 Viscosidades menores a 5 cp.  Gastos mayores a 400 bpd.
84
Fancher y Brown obtuvieron una correlación para tuberías de 2 pulgadas de

diámetro interno, y se ajusta con mayor precisión que la correlación de Poettmann y
Carpenter para gastos bajos y altas RGL.
1. Se debe tener u obtener p1, p2, T1, T2, R, γg, γo, γw, qo, qw, h, d así como Bo y Rs a p1
y p2.
2. Graficar h vs p.
3. Determinar la masa de aceite, gas y agua por barril de aceite.
M  350  γo  0.0764  γg  R  350  γw  WOR (3.2)
4. Calcular el gasto másico de la mezcla wm en lbm/día.
5. Suponer una presión inicial en la TP.
6. Calcular el volumen de la mezcla VT a la presión y temperatura establecida:
V  5.61 B  5.61 WOR  R  R 14.7  T  Z  (3.3)

s 
p  520  1 
T o
   

7. Determinar la densidad a la presión supuesta:
M
m  (3.4)
VT
8. De la figura 3.1, determinar el gradiente de presión fluyendo para el gasto másico y

la densidad, calculadas en los pasos 4 y 7 respectivamente.
9. Suponer una segunda presión en la TP para una mayor profundidad y repetir los
pasos 6, 7 y 8.
10. De los dos gradientes obtenidos, determinar un gradiente de presión promedio:
11. De la diferencia de presiones entre p1 y p2, determinar la distancia entre estas dos
presiones con:
p2  p1  dp
 h (3.5)
 dp   dp 
   
 dh  prom dh  prom
12. Repetir el procedimiento hasta completar la profundidad del pozo.
13. Para obtener el gradiente de presión incluyendo las pérdidas por fricción tenemos:
85
1.4737 105  w
d  v    m
(3.6)
d/12
Donde wm es el calculado en el paso 4, en lbm/día
14. De la figura 3.2 determinar el valor del factor de fricción f, con el valor obtenido en
el paso 13.
15. Calcular el gradiente de presión total con:
 dp  1   f w
2


  

 
 m
   m 
  (3.7)


 dh T 144  
7.41310   m  d/12  
10 5
Donde:
m1  m 2
m  (3.8)
2
16. Por último calculamos la profundidad total con:
p1  p2 dp

  Δh (3.9)
   dp 
dp
   
 dh T  dh  T
Figura 3.1. Correlación de la densidad del fluido (Poettmann).
86
Figura 3.2. Correlación de datos de campo en pozos fluyentes y de Bombeo Neumático.
87
3.2.2 Duns y Ros

El método de Duns y Ros es el resultado de una investigación de laboratorio a gran
escala con modificaciones y ajustes usando datos de campo.
Duns y Ros escogieron una aproximación un poco diferente que la mayoría de los
investigadores. El gradiente de presión es expresado como una fracción del gradiente de
líquido hidrostático (  L g ). Ellos definieron arbitrariamente el gradiente de presión estática
como el peso del volumen por la densidad in-situ y desarrollaron correlaciones para la
fricción en la pared de la tubería de sus extensos datos de laboratorio para cada una de las
tres amplias regiones de flujo. Aunque ellos usaron en un punto específico un balance de
presión en lugar de un balance de energía, sus ecuaciones son un balance de energía
termodinámico. El gradiente total incluye un gradiente estático, un gradiente de fricción y
un gradiente por aceleración. Los efectos de resbalamiento entre el gas y el líquido son
incorporados en el gradiente estático y se mantienen separados de los efectos debido a la
fricción. Ellos separaron el flujo dentro de tres tipos de regiones y prepararon correlaciones
separadas para el resbalamiento y fricción en los tres. Las tres regiones son:
 Región I: La fase líquida es continua y el flujo burbuja, flujo tapón y parte

del flujo burbuja existe en este régimen.
 Región II: En esta región las fases de líquido y gas se alternan. La región
por lo tanto cubre del patrón de flujo bache y el resto del flujo burbuja.
 Región III : En esta región el gas es la fase continua por lo que en esta
región se encuentra el flujo neblina.
Procedimiento de cálculo:
1. Calcular la densidad relativa del aceite.
2. Obtener la masa asociada con un barril de líquido a condiciones estándar:

 1   WOR 
M  γ 350  γ 350  0.0764 R γ (3.10)
o   w   g
1 WOR 1 WOR
   

3. Determinar la densidad de la fase líquida:

 
   1   WOR  

 γ w  
1 WOR 
  
(3.11)
L 62.428 γ  
 o 1 WOR     

4. Obtener la presión promedio en psia:
p1  p2
p p atm (3.12)
2
88
5. Obtener la temperatura promedio:
T1  T2
T (3.13)
2
6. Obtener Z
7. Calcular la densidad promedio de la fase de gas:


  γ  0.0764  520   p

  1   (3.14)
 T  460  14.7   Z 

g g
    

8. Calcular la viscosidad promedio del aceite con las ecuaciones 1.51 – 1.57.
9. Calcular la viscosidad promedio del agua con las ecuación 1.110.
10. Calcular la viscosidad de la mezcla del líquido:


μμ  1  
 WOR   (3.15)
o   μw  
  

1  WOR 1  WOR 
L  
   

11. Calcular la tensión superficial de la mezcla de líquido:


   1 
   WOR  

   w
  (3.16)
1  WOR 1 WOR 
L o  
   

12. Calcular Rs a p y T .
13. Calcular Bo a p y T .
14. Calcular el área transversal de la tubería.
15. Calcular el número de viscosidad del líquido:
 1 0.25

NLμ  0.15726 μ L  (3.17)
L σ 
3
16. Calcular la velocidad superficial del líquido vsL (suponiendo Bw = 1.0):
v sL 
5.6142144q L B 1   B  WOR  
86,400  A  o  1  WOR     (3.18)
 1  WOR  
w
p   
89


17. Calcular la número de velocidad del líquido: 


  0.25
NLv  1.938  vsL  L  (3.19)
 σ L 
18. Calcular la velocidad superficial del gas:


144    1  
q L R  Rs 


1  WOR  14.7  T  460  Z 
 
 
v     (3.20)
 
  
sg
86,400  A p  p  520
   1 

19. Calcular el número de velocidad del gas:
 
0.25
Ngv  1.938  vsg  L  (3.21)

 σ L 
20. Calcular el número de diámetro de la tubería: 


120.872  d   0.5
Nd   L  (3.22)
12  σ L 
21. Seleccionar el régimen de flujo adecuado de la figura 3.3.
22. Determinar el factor de resbalamiento dependiendo de la región de flujo obtenida en

el paso 21.
a) Para la Región I:
El resbalamiento se obtiene de la siguiente ecuación:
2

SFFN  F'  Ngv  


(3.23)
1 2 Lv 3
 1 N Lv 

Donde F1, F2, F3, y F4 se obtienen de la figura 3.4 y:
F'  F  F4 (3.24)
3 3
Nd
Para flujo anular Nd se basa en el perímetro mojado, por lo tanto:
90

d  d c  d t  (3.25)
Los límites de la Región I son: 


0  N gv  L1  L 2  N Lv  (3.26)
Donde L1 y L2 se obtienen de la figura 3.5.
b) Para la Región II:

N
S  1 F  
gv
 F '6
0.982
(3.27)
1 F7  NLv 2
5
F5, F6 y F7 se obtienen de la figura 3.6 y donde:
F6'  0.029  N d  F6 (3.28)
Los límites de la Región II son: 


L1  L 2  N Lv   N gv  50  36  N Lv  (3.29)
c) Para la Región III:
S=0 (3.30)
Por lo tanto: 

1
HL  (3.31)
1  vsg /vsL
El límite de la Región III es: 



N gv  75  84  N
0.75
Lv  (3.32)
23. Determinar la velocidad de resbalamiento si el flujo se encuentra dentro de la región

I ó II:
S
vs  (3.33)
1.938 L /σL 0.25
24. Determinar el colgamiento de líquido:
91
H
vs

vsg

vsL  v s  
vsg

vsL 2  4 v s  vsL 
0.5
(3.34)
L
2 vs
25. Determinar el Número de Reynolds para el líquido:
1488 L  vsL  d
N Re  L 12 μ (3.35)
L
26. Determinar el gradiente de fricción de acuerdo a la región de flujo.
a) Para la Región I y II:  



  2 f NLv N Lv  Ngv  (3.36)
fr w
Nd
Donde:
f w  f 1  f 2
f
(3.37)
3
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 f1  R  Nd   2/3
Donde:
vsg
R  (3.38)
vsL
f3 = 1 + f1 (R/50)0.5 (3.39)
El factor de fricción fw es válido para la región I y II. Es bueno para NLv = 0

y arriba del límite dado por Ngv = (50 + 36 NLv).
b) Para la Región III:
En el flujo niebla donde N gv

 75  84  N0.75
Lv
  2 f N N  gv
2
(3.40)
fr w 
Nd
Donde:
92
g 
N 
 
(3.41)
 L
En la Región III fw se toma como f1 y podría ser obtenida de la figura 3.7.

Para ε > 0.05d el valor de f1 es calculado con la ecuación:
f1  1  0.067 ε /d 
1.73
(3.42)
 4  log 0.027  ε /d  2
Para ε > 0,05(d), el valor de d  ε podría ser sustituido por d durante el

cálculo gradiente de fricción, y también esta sustitución puede se hecha por:
vsg  d 2
v  (3.43)
sg
d  ε  2
27. Determinar el gradiente estático: 


 st  H L  1 H L  N  (3.44)
Donde: 
g 
N 
 
(3.45)
 L
28. Determinar el gradiente de presión adimensional total:
a) Para las Regiones I y II:

 T   st  fr (3.46)
b) Para la Región III (tomando en cuenta la aceleración)


 st   fr   (3.47)
1 L  vsl  g  vsg vsg /p
29. Convertir el gradiente a psi/pie:
dp   st   L ; (3.48)
dhst 144
dp   fr   L ; (3.49)
dhfr 144
o:
93
dp  T   L
 (3.50)
dh T 144
30. Determinar la longitud para esa caída de presión:
Δp
Δh  (3.51)
dp
dh total
Figura 3.3. Regiones para la correlación de Duns y Ros.
Figura 3.4. F1, F2, F3, y F4 vs Número de la viscosidad NLμ (por Ros).
94
Figura 3.5. Factores L1 y L2 vs el Número de diámetro de la tubería (por Ros).
Figura 3.6. F5, F6, y F7 vs Número de viscosidad del líquido NLμ (por Ros).
95
Figura 3.7. f1 vs NRe (por Ros).
96
Figura 3.8. Corrección de la Fricción de Burbuja (por Ros).
3.2.3 Hagedorn y Brown (general)

En un esfuerzo hecho por Hagedorn y Bown para determinar una correlación
general la cual incluyera prácticamente todos los rangos de flujo, un amplio rango de gas
liquido, todos los tamaños de tuberías usados ordinariamente y los efectos de las
características de los líquidos. Los datos fueron tomados para rangos de tubería a partir de 1
pulgada a 2 1/2 pulgadas de tubería. El estudio incluyó todo el trabajo anterior hecho por
este equipo de investigadores en el efecto de la viscosidad líquida como “correlaciones
limitadas”. El término de energía cinética fue incorporado en la ecuación de energía por
que fue considerada para ser significante en tuberías de diámetro pequeño en la región
cerca de la superficie donde el fluido tiene baja densidad. Dos ajustes fueron encontrados
para mejorar esta correlación. La correlación Griffith fue usada para que el flujo burbuja
existiera y el colgamiento fue revisado para asegurarse de que excediera el colgamiento sin
resbalamiento y si no, fue usado.
Comenzar con una presión conocida p1, y asumir el valor para p2. Calcular el
incremento de profundidad.
1. Calcular la presión media entre los dos puntos de la presión, psia.
p1  p2
p 14.7 (3.52)
2
Dependiendo de los requisitos del problema, es decir, si o no la presión del fondo del
agujero que no fluye debe ser determinada de la información superficial, o si los
97
cálculos comienzan de la profundidad total y sube por la tubería, la presión inicial debe
ser conocida. Los incrementos o decrementos de la presión deben ser asumidos de entre
los puntos de presión (1) y (2) los cuales serán calculados.
2. Calcular la densidad relativa del aceite  o .
141.5
 o  (3.53)
131.5  API
3. Encuentra la masa total asociada con un barril de líquido a condiciones de tanque.
 0.0764 RGLγ  
 1   WOR 
m  γ  350  γ  350 (3.54)

  w  
1  WOR 1  WOR
o  
g
   

4. Calcular el flujo másico con:
wqm (3.55)
5. Con las ecuaciones del capítulo 1 obtener R s a la p y T .
6. Calcular la densidad de la fase líquida.
 Rs  γg  0.0764
 γo  62.4 
 WOR 
 1    62.4  (3.56)
ρ  5.614  γ
  w  
L  1  WOR 1  WOR 
 BO      
 
 
7. Obtener Z .
8. Calcular la densidad media de la fase gaseosa.



    0.0764  p    520    1 

(3.57)


14.7   T   

 
g g
Z
    

9. Calcular la viscosidad media del aceite con la correlación apropiada (capítulo 1).
10. Calcular la viscosidad media del agua.
11. Calcular la viscosidad líquida de la mezcla.

 1  
μμ   μ   WOR   (3.58)
o   w  
 
1  WOR 1 WOR 
L  
98
   
99



Esto puede ser solamente una aproximación, puesto que la viscosidad de dos
líquidos inmiscibles es absolutamente compleja.
12. Asumiendo la tensión superficial constante en cada punto de presión, calcular la

tensión superficial de la mezcla líquida.
 1  
       WOR   (3.59)
o  w 
 
1  WOR 
L  
1 WOR
   

Esto representa sólo una aproximación de la tensión superficial de la fase líquida.
13. Calcular el número de viscosidad del líquido NL:
 1 14
NL  0.15726  L    (3.60)
  
 L L 
3
14. Determinar CNL de la figura. 3.9.
15. Calcular el área de la Tp.
Π d 2
Ap  (3.61)
4
16. Obtener Bo a p y T .
pie
17. Asumiendo BW  1.0 , calcular la velocidad superficial del líquido vsL ,
seg
  1 
5.61 q L  B   B  WOR   (3.62)
 
v  O  1 WOR  W 

sL
86400  A 

1 WOR
p     

18. Calcular el número de velocidad del líquido NLV.
 L  4
1
NLV  1.938  VsL    (3.63)

 L 

19. Calcular la velocidad superficial del gas, vsg.
10
0
  1 
q  RGL  R  

L   
1 WOR   14.7   T

S
  Z 
 
 
v      
   
 
(3.64)
sg
86400  A P  p  520  1
     

20. Calcular el número de velocidad del gas, NGv.
 L  4
1
NGv  1.938  vsg    (3.65)

 L 

21. Revisar el régimen de flujo para determinar si se continúa con la correlación

Hagedorn y Brown o se procede con la correlación de Griffith para el flujo burbuja.
El numero “A” se calcula con la siguiente fórmula:


A  1.071- 0.2218  v sL
 vsg  2
 (3.66)
d
Si A  0.13 entones utiliza ese valor, y, si es menor a 0.13 usa A = 0.13

El número “B” se calcula con la siguiente fórmula:
vsg
B (3.67)
vsL  vsg
Si (B-A) es positivo o cero, se continua con la correlación Hagendorn y Brown.

Si (B-A) es negativo, sigue con la correlación de Griffith, el procedimiento se
encuentra en el método de Orkiszewski.
Nota: este procedimiento sigue de cerca los incrementos de presión. Es posible que
el procedimiento del cálculo cambie de un régimen a otro.
22. Encuentra el número de diámetro de la tubería, Nd.

L
Nd  120.872  d  (3.68)
L
23. Calcular el colgamiento en función de  :
 N   p 0.10  CN 
   Lv       L  (3.69)
 gv   14.7   Nd 


 N 


HL
24. Obtener de la figura 3.10, el valor de .
ψ
100
25. Determina la segunda corrección en función de  :
 NLv N0.380 

  L (3.70)
 N 2.14 
 d 

26. Obtener ψ de la figura 3.11.
27. Calcular el valor de HL:

 HL 

HL     ψ (3.71)
 ψ 
Para bajas viscosidades no habrá corrección y ψ = 1.00
28. Determinar el valor para dos fases del número de Reynolds, (NRe)TP. 

(N Re ) TP  2.2 102 w
 
d   μ H  μ 1H 
L
L g
L
 (3.72)

29. Determinar el valor para 
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.
30. Obtener el factor de fricción con la figura 3.12.
31. Calcular la densidad media de la mezcla  m :
(a) Utilizando el valor de HL del paso 27, calcular  m con la ecuación 2.84
32. Repetir los pasos 5, 7, 16, 17, y 19 para p1 y p2.
33. Calcular la velocidad de la mezcla para p1 y p2.
vm1  vsL1  vsg1 (3.73)
vm2  vsL2  vsg2 (3.74)
34. Determinar el valor para  vm2  

  
 v2  v2  v2
m m1 m2
 (3.75)
101
35. Calcular  h correspondiente a p  p1 - p2.
 v2 
144 Δp  ρm Δ   m 
 2 
h   (3.76)
gc
fw2
 m 
2.9652 1011 d5  m

36. Comenzando con p2, y la profundidad conocida a la p2, asumir otro punto de presión
y repetir el procedimiento hasta alcanzar la profundidad total o hasta alcanzar la
superficie.
Figura 3.9. Correlación para el coeficiente C del número de viscosidad (por Hagedorn).
Figura 3.10. Correlación para el factor de colgamiento (por Hagedorn).
102
Figura 3.11. Correlación para la segunda corrección (por Hagedorn)..
Figura 3.12. Correlación del factor de fricción (por Hagedorn).
103
3.2.4 Orkiszewski
La correlación de Orkiszewski fue resultado de un análisis de muchos de los
métodos publicados para determinar si alguno de ellos son lo bastante extensos para
predecir con precisión la caída de presión para una amplia gama de condiciones del pozo.
Orkiszewski enfatizó que el colgamiento de líquido fuera derivado del fenómeno
físico observado y que el gradiente de presión fuera relacionado a la distribución
geométrica de la fase de líquido y gas, por lo que determinó que la densidad de la mezcla se
determinara mediante el colgamiento, considerando en ella el resbalamiento entre las fases.
El factor de fricción se correlacionó con las propiedades del fluido en la fase continua.
Él reconoció cuatro tipos de patrones de flujo e hizo correlaciones separadas para
establecer la velocidad de resbalamiento y fricción para cada uno de los patrones. Los
patrones son burbuja, bache, transición (bache-niebla) y niebla.
Considerando similitudes en conceptos teóricos y diferentes categorías, Orkiszewski
comparó cinco métodos seleccionados y determinó la cantidad de desviación entre las
caídas de presión medidas y predichas. Orkiszewski observó que los mejores resultados,
bajo ciertas condiciones de flujo, se obtenían con los métodos de Griffith y Wallis, y Duns
y Ros, por lo que tomó estas correlaciones como base para desarrollar su método,
combinándolas para los diferentes patrones de flujo considerados.
1. Selecciona el punto donde se va a comenzar, ya sea en la cabeza del pozo o en el

fondo del pozo.
2. Determinar el gradiente de temperatura del pozo.
3. Fijar una Δp que puede ser la de la superficie o la del fondo del incremento.
Encontrar la p de ese incremento.
4. Suponer un incremento de profundidad Δh, y determinar la profundidad promedio

h del incremento.
5. Del gradiente de temperatura, determinar la T .
6. Determinar las propiedades de los fluidos a p y T , así como

L , g , vsL , vsg , vm , μ L , μg , Ngv , NLv con las ecuaciones de los capítulos anteriores.
7. Determine LB, LS y LM para poder obtener el tipo de régimen de flujo de la siguiente

tabla:
104
Tabla 3.2. Límites de los patrones de flujo

Limites Régimen de Flujo
vsg
L B Burbuja
vm
vsg
 L B, N gv  LS Bache
vm
LM  Ngv  LS Transición bache-niebla
Ngv  LM Niebla
Donde LB, LS, y LM son los límites de burbuja-bache, bache-transición, transición-

niebla respectivamente, adimensionales, y Ngv se calcula con la ecuación 2.72.
Donde:
 vm2 
LB  1.071  2.6616  (3.77)
 d h 
Con el límite LB  0.13
Y:
dh = diámetro hidráulico de la tubería = área transversal/perímetro mojado, pg.
LS  50  36  NLv (3.78)
L M  75  84  NLv
0.75
(3.79)
8. Basándose en el régimen de flujo obtenido en el paso 7 determinar el gradiente de

elevación y de fricción:

a) Régimen de burbuja
 
    H L   g 1 H L 
 Δp  1
(3.80)
ΔL 144 L
 e
Donde HL se calcula con:
C1  C 2
HL  1  (3.81)
2
105
Y: 

 2 4 
0.5
C  C  v  (3.82)
2 1 sg 
 0.8 

vm 1
C1 (3.83)
0.8
El gradiente de fricción se calcula con:
 Δp  1  12  f    v 
2
 
  L 
L  (3.84)
 ΔL f 144  2 g  d h 

El factor de fricción se calcula utilizando un proceso iterativo con las

ecuaciones 2.20 ó 2.21 según sea el caso y con:
124  dh  vL  L
N Re  L μ (3.85)
L
b) Régimen de bache
El gradiente de elevación se obtiene de acuerdo al procedimiento delineado

por Griffith y Wallis:
 Δp   C3 
   1    L  δ (3.86)
ΔL  
 e 144  v m  v b 

Donde:
C 3   L v sL  v b    g  v sg (3.87)
El término δ se conoce como el coeficiente de distribución del líquido, el

cual considera los siguientes fenómenos físicos:
1) El líquido está distribuido en tres espacios: el bache, la película alrededor

de la burbuja de gas y dentro de la misma como gotas atrapadas. Un
cambio en su distribución cambiará las pérdidas netas por fricción.
2) Las pérdidas por fricción están constituidas esencialmente por dos

componentes, una corresponde al bache de líquido y la otra a la película
del mismo.
106
3) La velocidad de elevación de la burbuja se aproxima a cero conforme el

flujo tiende al tipo burbuja.
El coeficiente de distribución de líquido (δ) se calcula como se indica en la

siguiente tabla:
Tabla 3.3. Relación de la ecuación a aplicar, con la velocidad de la mezcla y

la fase continua.
ECUACIÓN A
FASE CONTINUA vm
APLICAR
Agua < 10 (3.88)
fw > 0.75 >10 (3.89)
Aceite <10 (3.90)
fo >0.25 >10 (3.91)
Las ecuaciones para obtener δ son:
1.38  d 
δ  0.681  0.013  h 
d
 log μ L  0.232  log v m  0.428  log  h  (3.88)
 12   12 
 dh 
0.799  d 
δ  0.709  0.0451   log μ L  0.162  log v m  0.888  log  h  (3.89)
 12   12 
1.415   dh 
 dh  
 log
  
μ L  1  0.167  log v m  0.113  log 


δ  0.284  0.0127  12  

(3.90)
  12
 

1.371   dh 
 dh  

 log μ
 
 1  0.569  log  

δ  0.161  0.0274  12 

L
   12
 
 1.571   dh 
 dh  

  
log μ  1  0.63  log   log v

 0.397  0.01  12 

L m (3.91)
12
    

El valor de δ debe estar dentro de los límites siguientes:
Para vm < 10
δ  0.065  vm (3.92)
Para vm > 10
107
vb   m 

δ   1   (3.93)
vm  v b  L 
El valor de vb (velocidad de elevación de la burbuja o velocidad del bache,

pies/seg) se determina por ensaye y error con las ecuaciones siguientes:
124 dh  vm  L
N Re 
L μ (3.94)
L
124  dh  vb  L
N Re 
b μ (3.95)
L
El procedimiento para calcular la velocidad del bache se inicia suponiendo

un valor de vb = 1.75 pies/seg, y se compara con el calculado de la ecuación
3.96, 3.97, 3.99 ó 3.100 (según sea el caso) hasta que el supuesto sea igual al
calculado:
Cuando NRe  6000 tenemos los siguientes tres casos para

L
NRe :
b
1. Si NRe  3000 :
b


  g  dh 
0.5
6
vb  0.546  8.74 10  NRe   (3.96)


 12 
L
2. Si 3000  NRe  8000 : 


b
 13.59 μ 0.5
L 
v b 0.5 α   α 
2
(3.97)
  L d h /12 
0.5

Donde:

 6   g  dh 
0.5
α  0.251  8.74 10  NRe   (3.98)


 12 
L
3. Si NRe  8000 :
b


 g  d 
0.5
h
6
vb  0.350  8.74 10  NRe   (3.99)


 12 
L
108
Cuando NRe  6000 y NRe  32.5 se utiliza la siguiente ecuación para

L b
determinar vb:
 g  d 0.5
v b  C1 C 2 h   (3.100)
 12 

Donde:
C 2 1.36  C 5 C C
6  C
3 C
2
7 3 (3.101)
C 7 1.16110-3  4.6 10-5 C  2.954 10-3 C  5.5 10-4 C  6.67 10-4 C

2 3 4
4 4 4 4
(3.102)
C 6 0.0413  0.01122 C 4 0.012 C 1.110 C 1.118 10-3 C

-32 3 4
4 4 4 (3.103)
C 5 0.220623  0.03408 C 4 9.549999 103 C 4 8.283001103 C 4 2.645 103 C

2 3 4
4
(3.104)
N Re  5,500
C4  b
(3.105)
1,000
N Re  3,000
C3  L
(3.106)
1,000
C1 0.013805  0.4246 C 8 0.1753 C 8 0.02363 C

2 3
8 (3.107)
N Re
C8  b
(3.108)
10
Y si NRe  6000 y NRe  32.5 ; C1 = 0.351 y calculamos C2 con las

L b
ecuaciones 3.101 a la 3.106 y vb con la ecuación 3.100.
El gradiente por fricción se obtiene con la ecuación:
 Δp  f  v  v  v 
2
    
 m L
 sL b
 δ (3.109)
 ΔL  f 12   2  g  d h  v m  v b 

El factor de fricción se puede calcular con la ecuación 2.21 y el número de

Reynolds con la ecuación:
109
124 dh  vm  L
N Re 
L μ (3.110)
L
c) Régimen de transición bache-niebla.
Para este caso Orkiszewski adoptó el método de interpolación propuesto por

Duns y Ros que consiste en calcular (Δp/ΔL)e y ( Δp/LΔ)f en las fronteras
para flujo bache y flujo niebla, para luego ponderar linealmente cada término
respecto al valor de Ngv.
La zona de transición está definida por:
Lm > Ngv >Ls (3.111)
Donde:
L m  84  NLv  75
0.75
(3.112)
El valor del término por elevación, está dado por:
 Δp   Lm  Ngv  Δp  
Ngv - Ls  Δp   (3.113)
  
  L  L   LL  
  
ΔL ΔL  ΔL 
 e m s  e m s  e
bache niebla
Y el término por fricción, por:
 Δp   Lm  Ngv  Δp  
Ngv - Ls  Δp   (3.114)
 
  L  L   LL  
 
ΔL ΔL f ΔL 
 f m s  m s  f
bache niebla
Se obtiene un valor más preciso del factor de fricción en la región de niebla,

si el gasto de gas se obtiene con la siguiente ecuación:
0.25
  
q  A p  L mL g (3.115)
g σ L 
d) Régimen de niebla.
Para calcular el gradiente de presión correspondiente a esta región se aplica

el método de Duns y Ros.
La región de niebla está definida para:
Ngv > Lm (3.116)
110
El gradiente por elevación, dado que el líquido va en suspensión dentro de la

corriente de gas y no existe diferencia de velocidad entre las fases, se
calcular:
 Δp  1  L  vsL  g  vsg 
   144  v  (3.117)
ΔL
 e  m 

En el gradiente por fricción, se considera que la mayor parte de las caídas de

presión por fricción se deben al flujo de gas por la tubería.
 Δp  12  f  v  
2
  

sg g (3.118)
 ΔL  f 2  g  dh
El factor de fricción se puede calcular con la ecuación 2.21 y con el número

de Reynolds de la ecuación:
124  dh  vsg  g
N Re 
L μ (3.119)
g
En este caso la rugosidad relativa ε/d se determina a través de una función del
número de Weber según los lineamientos establecidos por Duns y Ros, quienes señalan que
sólo será significativo cuando su valor esté comprendido entre 1 x 10 -3 y 0.5. Entre estos
límites, ε/d se calcula de la siguiente manera:
Si:
Nw  0.005 (3.120)
 0.8994 σ
d  g  vsg2  d L (3.121)
h
Y si:
Nw  0.005 (3.122)
 4.6242  σ
 N 0.302
L w (3.123)
d g  vsg 2  d h
Donde:
111
 μ  
2
Nw  0.093 g 
 L 
 (3.124)
vsg
L  σ L 

9. Calcular ΔL de la siguiente ecuación:
 Δp 
  Δp 
  
 Δp     

 Δh   Δh 
   e f (3.125)
 Δh  T 1 E k
Donde: 
Wm  vsg
Ek  (3.126)
4,637  p  Ap
Ek es el término donde se incluyen las caídas de presión por aceleración.
10. Si la Δh calculada en el paso 9 no es igual a la supuesta en el paso 4, entonces

repetir el procedimiento a partir del paso 4 tomando como supuesta la Δh calculada
en el paso 9. Este procedimiento se repite hasta que Δhsup = Δhcal.
11. Determinar los valores de p y Z para ese incremento Δh.
12. Repetir el procedimiento del paso 3 hasta completar la profundidad total del pozo.
3.2.5 Beggs y Brill

La correlación de Beggs y Brill (1973) fue desarrollada en 584 pruebas tomadas de
datos obtenidos experimentalmente de una prueba de arreglo a pequeña escala. La prueba
consistió en una sección de tubería de acrílico de 1 pg y 1.5 pg de diámetro y de 90 pies de
longitud, la cual tenia un mecanismo que podía inclinar la tubería de horizontal a vertical y
los fluidos utilizados eran aire y agua. Los parámetros estudiados y sus rangos de variación
son:
 Gasto de gas, 0 a 300 Mpies3/día;
 Gasto de líquido, 0 a 30 gal/min (0 a 1.635 x 106 litros/día);
 Presión promedio del sistema, 35 a 95 psia;
 Diámetro de la tubería, 1 y 1.5 pg;
 Colgamiento de líquido, 0 a 0.870;
 Gradiente de presión, 0 a 0.8 psi/pie;
 Ángulo de inclinación, -90o a +90o;
 Patrón de flujo horizontal.
Para cada diámetro de tubería, los gastos de líquido y gas variaban por lo que se
pudieron observar todos los patrones de flujo cuando la tubería estaba en posición
112
horizontal. Una vez establecido cada patrón de flujo se procedió a variar el ángulo de
inclinación, así que se pudo observar como el ángulo de inclinación afectaba el colgamiento
y el gradiente de presión. El colgamiento y el gradiente de presión 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 máximo en +50 grados y a su valor mínimo en -50 grados.
El mapa de patrones de flujo original que obtuvieron Beggs y Brill fue ligeramente
modificado para poder incluir la zona de transición entre el patrón 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 ecuación para determinar el gradiente de presión es:
g fT  M m  v m 12 
   m  sen θ  
Δp g 2 g d 
  c c 
 (3.127)
Δh   m  v m  vsg 

144 1- 
 g c  p  144 


Figura 3.13. Mapa de patrón de flujo horizontal modificado.
Resolviendo la ecuación 3.127 para el incremento en la profundidad, Δh:
113

Δp    m  vm  vsg
144 1 g  p 14.7144

Δh  c  (3.128)
g  m f  Mm  vm 12
senθ  T
gc 2 g c d
Comenzando con una presión conocida p1, seleccionar el valor para p2 y suponer un
valor Δh entre los dos puntos de modo que la T , que está en función de la profundidad,
pueda ser evaluada. Calcular el incremento en la profundidad, que debe coincidir con el
valor supuesto.
1. Calcular la presión (psia) y profundidad promedio entre los dos puntos de presiones
dadas:
p1  p2
p p atm (3.129)
2
Δh
hh (3.130)
2
para p1
2. Determinar la temperatura promedio T , a la profundidad promedio. Este valor se

debe conocer de la relación temperatura vs profundidad (gradiente de temperatura).
3. Del análisis PVT o las correlaciones apropiadas, calcular a la p y T :
Rs , Bo , Bw , μo , μw , μg , σo , σw , Z (3.131)
4. Calcular la densidad relativa del aceite.
141.5
γo  (3.132)
131.5  o API
5. Calcular las densidades del líquido y del gas en lbm/pie3 a p y T .


  1   WOR  
   (3.133)
   
1  WOR 1 WOR 
L o w  
   
350  γ  0.0764  R s  γ g 

  o
(3.134)
5.6146  B o
o
114
350 γw
 w  (3.135)
5.6146 Bw
0.0764  γg  p  520
 g 
14.7T  460Z
(3.136)
6. Calcular los gastos de gas y líquido in situ (a condiciones de escurrimiento).
3.27 107  Zq o R - R s T  460 

q  (3.137)

g
p

q L 6.49105 q  oB oq  Bw w (3.138)
Donde:
qL y qg = pies3/seg
7. Calcular las velocidades superficiales del gas, líquido y la mezcla in situ.
144  q L
v sL  (3.139)
A p
144  q g
v sg  (3.140)
Ap
vm  vsL  vsg (3.141)
8. Calcular el flujo másico del gas, líquido y total.
GL  L  vsL (3.142)
Gg  g  vsg (3.143)
GT  GL  Gg (3.144)
9. Calcular el colgamiento de líquido sin resbalamiento con la ecuación 2.73.
10. Calcular el Número de Froude, NFR, la viscosidad del líquido y de la mezcla, y la

tensión superficial del líquido.
v m2

N FR  (3.145)
g  d/12
115
 1
μ μ  
 WOR   (3.146)
o  μ w  
1  WOR 1 WOR 
L 
   


μ m  μ L  λ  μ g1 λ   (3.147)

  1   WOR  

   w 
(3.148)
o 

1 WOR 
L  
1 WOR
   

11. Calcular el Número de Reynolds sin resbalamiento y el número de velocidad del

líquido.
M m  d/12
NRe  (3.149)
μ m 6.72 104
NLv  1.938 vsL  L

0.25
 
(3.150)

 L 

12. Para determinar el patrón de flujo que existe en el flujo horizontal, calcular los
parámetros correlacionados, L1, L2, L3, y L4.
L1  316  λ0.302 (3.151)
L2  0.0009252 λ2.4684 (3.152)
L
3
 0.10 λ-1.4516 (3.153)
L4  0.5 λ-6.738 (3.154)
13. Determine el patrón de flujo usando los siguientes límites:
λ < 0.01 y NFR < L1

Segregado ó
λ  0.01 y NFR < L2
Transición λ  0.01 y L2 < NFR  L3
0.01  λ < 0.4 y L3 < NFR  L1
Intermitente ó
λ  0.4 y L3 < NFR  L4
λ < 0.4 y NFR  L1
Distribuido ó
λ  0.4 y NFR > L4
116
14. Calcular el colgamiento horizontal, HL (0). Si el patrón de flujo es transición, es

necesario interpolar entre los valores de flujo segregado y el intermitente.
a  λb

H L 0  N c (3.155)
FR
Donde a, b y c son determinados para cada patrón de flujo de la tabla siguiente:
Tabla 3.4. Coeficientes para determinar el colgamiento según el patrón de flujo.

Patrón de Flujo a b c
Segregado 0.98 0.4846 0.0868
Intermitente 0.845 0.5351 0.0173
Distribuido 1.065 0.5824 0.0609

15. Calcular el coeficiente del factor de corrección por inclinación:


C  1  λ ln d  λ   N Lv   N FR 
e f g
 (3.156)
Donde d, e, f y g se determinan para cada condición de flujo de la tabla:
Tabla 3.5. Coeficientes para la corrección por inclinación.

Patrón de Flujo d e f g
Segregado ascendente 0.011 - 3.768 3.539 - 1.614
Intermitente ascendente 2.96 0.305 - 0.4473 0.0978
Distribuido ascendente Sin Corrección (C = 0)
Todos los patrones de
4.70 - 0.3692 0.1244 - 0.5056
flujo descendente
16. Calcular el factor de corrección del colgamiento de líquido debido a la inclinación:



ψ  1  C sen1.8 θ   0.333  sen3 1.8 θ  (3.157)
Para pozos verticales se usa:

ψ  1 0.3 C (3.158)
17. Calcular el colgamiento de líquido corregido y la densidad de la mezcla con:
H L θ   H L 0   ψ (3.159)
117

m  L  HL  g 1 HL  (3.160)
18. Calcular la relación del Factor de Fricción de las dos fases (fT) con respecto al
Factor de Fricción sin resbalamiento (fns).
fT S
e, (3.161)
fns
 
Donde:

lny 

S 
 0.0523  3.182  ln y  0.8725lny  2
 0.01853 lny 
4
 (3.162)
Y:
y λ (3.163)
 2
HL θ

S se indetermina en un punto del intervalo 1 < y < 1.2; para “y” en este intervalo, la
función S se calcula de:

S  ln 2.2 y 1.2 (3.164)
19. Calcular el Factor de Fricción sin considerar el resbalamiento.
1
f ns  2 (3.165)
 Re N 
2 log   
  4.5223 log N Re  3.8215  

 o: 
 
0.5
fns  0.0056  (3.166)
N Re 0.32
20. Calcular el factor de fricción de las dos fases.
f
fT  f ns f T (3.167)
ns
21. Calcular Δh:
118

Δp      m  v m  vsg 
144 1 g  p 14.7144

Δh  c  (3.168)
g  m f  Mm  vm 12
sen θ  T
gc 2 g c d
Si el valor supuesto en el paso 1 y el calculado en el paso 21 no son suficientemente

cercanos, el valor calculado es tomado como el nuevo valor supuesto de Δh y el
procedimiento se repite hasta que los valores sean iguales. Un nuevo incremento de presión
es escogido y el procedimiento continúa hasta que la suma de todas las Δh’s sea igual a la
profundidad del pozo.
3.3 Ejemplos
3.3.1 Método Poettman y Carpenter.
1. Dados los siguientes datos:
d = 1.995 pg.
T1 = 120 oF
T2 = 150 oF
Presión en la cabeza = p1 = 500 psig
Presión de fondo = p2 = 1,000 psig
Presión atmosférica = 14.7 psia.
qo = 400 bpd
qw = 600 bpd
R= 500 pies3/bl
o = 22 oAPI
γg = 0.65
γw = 1.07
σw = 70 dinas/cm. (cte)
σo = 30 dinas/cm. (cte)
μg = 0.018 cp. (cte)
2. Para este ejemplo no se realiza grafica.
3. Determina la masa en barriles de aceite a condiciones estándar.
119
M  3500.922 0.076412500.65 3501.071.5

M  323  62  562
lb
M  947 m
blO
4. El gasto total de masa por día (wm) es:
lbm
w m  400947  378,800
dia
5. Comienza los cálculos con p1 = 500 psig.
6. Calcula el volumen del aceite, gas y agua a 500 psig. (Bo y Rs calculados con
ecuaciones capitulo1)
pies3
Bo  1.043 Rs  59
bl  14.7  580 
V  5.611.043 5.611.5 1250  59 0.935

T   
 514.7  520 
 5.88  8.43  35.69
 50 pie3
7. Densidad a 500 psig.

947 lbm
   18.96
m1
50 pie3
8. El gradiente de presión, de la figura 3.1

dp psi
 0.176 a 500 psig
dh pie
9. Determinar el gradiente de presión a 1000 psig.
Bo  1.083
pie3
Rs  120
bl
Z  0.90
V  1.0835.61 1.55.61 1130  14.7  610  0.90 
 31.8pie3
T    
1014.7  520  1 
120


947 lbm
   29.8
m2
31.8 pie3
De la figura 3.1 se obtuvo:

dp psi
 0.236 a 1000 psig
dh pie
10. Gradiente de presión promedio.
 dp  0.176  0.236 psi

  0.206
  pie
 dh  prom
2
11. La distancia entre 500 y 1000 psi es:
p2  p1 1000  500
  2,427 pies
 dp  0.206
 
 dh  prom
12. La profundidad del pozo sin incluir las pérdidas por fricción es de 2,425 pies.
13.

1.4737 105 378,800
d  v     33.6 
1.995/12
14. De la figura 3.2 se obtuvo el valor de f = 0.0080.
15.
m1  m2 18.96  29.8 lb

m    24.38 m3
2 2 pie
 dp   1 24.38   0.008 378,800 2  

5 
psi
  0.205
144 7.41310 24.381.995 /12
10
pie
dh
 T   
16.
p2  p1 1000  500
  2,439 pies
 dp  0.205
 
 dh  T
121
3.3.2 Método Duns y Ros

Calcular la longitud para una Δp dada con los siguientes datos:
Diámetro nominal = 1.25 pg.

Diámetro de flujo = 1.38 pg
Temperatura superficial = T1 = 120 oF
Temperatura en el fondo = T2 = 150 oF
qo = 400 bpd
qw = 600 bpd
R = 500 pies3/bl
 o= 22 oAPI
γg = 0.65
γw = 1.07
σw = 70 dinas/cm.
σo = 30 dinas/cm.
μg = 0.02 cp.
1.
141.5 141.5
γo    0.9218
131.5  API 131.5  22
o
2.
qw 600
WOR    1.5
qo 400
 1   
M  0.9218350  1.07350  1.5  0.07645000.65  378.582 lbm
   
11.5 11.5 bl
    
o c.s.
3.
 1 1.5  lbm
  62.428 0.9218  1.07  63.0972
L  2.5 2.5  pie3
4.
p1  p2 500 1000
p p  14.7  764.7 psia
atm
2 2
5.
T1  T2 120 150
T   135 0 F
2 2
6.
Z  0.912
122
7.    
  0.650.0764 
1 lbm
 764.7 
520  2.4755
  

135  460
g


14.7 pie3
0.912 
   

8. Con las ecuaciones 1.51 a la 1.57 se obtuvo μ o :
μo  7.7813 cp
9. Con la ecuación 1.110 se obtuvo:
μw  0.5313 cp
10.     1   1.5 
μμ   1 
 WOR  7.7813  0.5313  3.4313 cp
 μ
o  w     
  

1  WOR 1  WOR 2.5 
L 
2.5
       

11.  1    
   
 WOR  30  1   1.5  dinas
o   w     70   54
1  WOR 1  WOR
L
2.5 2.5 cm
       

12.
Rs = 107.22 pies3/bl
13.
Bo = 1.0736
14. 
π  d 2 π 1.382
Ap    1.4957 pg 2
4 4
15.  
 0.25
 0.157263.4313
NLμ  0.15726 μ L  1 
0.25
  1   9.6112 103
  L σ 3 
 63.0972543 
 
L
16. Suponiendo Bw = 1.0
v sL 
5.61421441000   1   1.5  pies
86,4001.4957 1.0736 2.5   1 2.5   6.44 seg
    
123


17. 0.25
 63.0972 

0.25
NLv  1.938  vsL  L   1.9386.44   12.9761

 σ L   54 
18. 
 1 
  
107.22 
 

144 1000 500 1 1.5  14.7  135  460  0.912 pies
        10.2179
v  

86,400 1.4957 764.7 
sg
520  1  seg
19. 0.25
 63.0972 
 L 
0.25
Ngv  1.938  vsg    1.93810.2179   20.5882

σ
 L   54 
20.
120.8721.38 63.0972 
0.5
Nd     15.0256
12  54 
21. De la figura 3.3 y con Ngv y NLv se tiene que el flujo se encuentra en la Región II.
22. Para la Región 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:
F6'  0.029  N d  F6  0.029 15.0256 0.1  0.3357
Por lo que : N
    20.5882   0.3357
0.982
 F' 
0.982
S  1 F5  1 0.19
gv 6
 4.7379
1 F7  NLv  2
1  0.09212.97612
23. 
S 4.7379 pie
vs    2.3514
1.938  L /σL  1.93863.0972 / 54 
0.25 0.25
seg
24.

2.3514 10.2179  6.44  2.3514 10.2179  6.44   42.3514 6.44 
2

0.5
HL  
2 2.3514 
124



HL  0.4210
Este valor de HL puede revisarse haciendo:
vsgvsL 10.2179  6.44  2.35 pie

vs  H  
1 HL L 1  0.4210 0.4210 seg
Como es igual al vs calculado en el paso 23, el HL calculado es el correcto.
25.
1488  L  vsL  d 148863.09726.441.38
N    20,264
12 μ L 123.4313
ReL
26. El gradiente de fricción para la Región II es:

De la figura 3.7 obtenemos:
f1 = 0.0064
Para poder obtener f2 primero calculamos:
vsg 10.2179
R   1.5866
vsL 6.44
f1  R  N d  2/3
 0.00641.5866 15.0256 
2/3
 0.0618
Con este valor, de la figura 3.8 obtenemos f2 = 1.05
f3 = 1 + f1 (R/50)0.5 = 1 + (0.0064)(1.5866/50)0.5 = 1.0011
Y:
f  f  f2  0.0064 1.05  6.7126 103
w 1
f3 1.0011
Por lo tanto:
  2 f N Lv N Lv  N gv   26.7126 10 3 12.976112.9761 20.5882  0.3891

fr w
Nd 15.0256
 fr  0.3891 Adimensional
27. El gradiente estático es:
125
 st  H L  1 H L  N   0.4210  1 0.42100.03923  0.4437 Adimensional
Donde:
g
N  
2.4755
 0.03923

L 63.0972
28. El gradiente de presión total es:
 T   st  fr  0.4437  0.3891  0.8328
29. El gradiente de presión total en psi/pie es:
dp

 T  L

0.832863.0972  0.3649 psi
dh T 144 144 pie
30. Determinar la longitud: 


Δp  1000  500  1,370.24 pies
Δh 
dp 0.3649
dh T
3.3.3 Método Hagedorn y Brown (general).

El siguiente problema se trabajó por el método Hagedorn y Brown para tubería de 2
pulgadas con un diámetro interno de 1.995 pulgadas. Como se muestra en el
procedimiento se determinará la distancia entre 500 psi y 1000 psi.
d = 1.995 pg.
Presión atmosférica = 14.7 psia.
qo = 400 bpd
qw = 600 bpd
R= 500 pies3/bl
o = 22 oAPI
γg = 0.65
γw = 1.07
σw = 70 dinas/cm. (cte)
σo = 30 dinas/cm. (cte)
μg = 0.018 cp. (cte)
126
Se requiere la distancia entre los puntos de presión 500 psi y 1000 psi.
Solución:
1. Peso específico del aceite.
141.5
o   0.9218
131.5  22
2. Encuentra la masa asociada con un barril de líquido a condiciones de tanque.

 1   1.5 
M  0.09218 350  1.07 350  0.0764 500 0.65  378.5854
   
11.5 11.5 
   

3. Calcular el flujo másico
lbm
w m  (400  600)  (378.5854)  378545
día
4. Calcular la densidad de la fase líquida.


 107.220.6510.0764  1.5 
L  5.614
0.9218  62.4   
     62.4    62.01 lbm
   
1.0736 11.5  1.07 11.5  pie3
     
 
5. Calcular la presión.
500  1000
Presión media = p   14.7  764.7 psia
2
6. Calcular Z.
120  150
Temperatura media = T   135o F
2
Z = 0.91
7. Calcular la densidad media de la fase gaseosa.

 764.7   520   1 
  0.65 0.0764    2.4810 lbm/cu pie
   
0.91 
g
14.7 595
     
127



8. Rs media es:
Rs  107.22
9. Calcular la viscosidad media del aceite:
 o  7.7813 cp
10. La viscosidad media del agua es:
 w  0.5313 cp
11. Determinar la viscosidad líquida de la mezcla:

 1   1.5 
μ  7.7813   0.5313    3.4313 cp

1  1.5  1  1.5

   

12. Asumiendo la tensión superficial constante en cada punto de presión, calcular la

tensión superficial de la mezcla líquida.
 1   1.5 
σ  30   70   54 dinas/cm
   
11.5 
L
11.5
   

Esto representa sólo una aproximación de la tensión superficial de la fase líquida.
13. Calcular el número de viscosidad del líquido NL.

1
 4  0.0097

NL  0.15726  3.4313  

1
3

 
62.0121 54 

14. Determinar CNL de la figura 3.9.
CNL  0.002

15. Calcular el área de la Tp.
 
Π  1.995
2
Ap   0.0217 pie2
144  4
16. Obtener Bo a p y T . Capítulo 1
128
BO  1.0736
129
pie
17. Asumiendo BW  1.0 , calcular la velocidad superficial del líquido vsL ,
seg
5.61 600  400   1.5 

v sL  86400  0.0217  1.0736 
1  pie
  1   3.08
  11.5   11.5 seg

18. Calcular el número de velocidad del líquido NLV.
 62.012114
N  1.938 3.08   6.1520
LV  
 54 

19. Calcular la velocidad superficial del gas, vsg.

   1
600  400   107.22  

 
500 1  1.5   14.7    595   0.91 pie
  
v sg   
  
    4.8784
 
86400  0.0217  764.7 

  520   1  seg
20. Calcular el número de velocidad del gas, NGv.
 62.0121 14
N  1.938  4.8784   9.7870
Gv  
 54 

21. Revisar el régimen de flujo para determinar si se continúa con la correlación

Hagedorn y Brown o se procede con la correlación de Griffith para el flujo burbuja.
El numero “A” se calcular con la siguiente fórmula:
A  1.071-
0.2218 3.08  4.8784   83.41
2
0.166
Si A  0.13 entones utiliza ese valor, y, si es menor a 0.13 usa A = 0.13
El número “B” se calcula con la siguiente fórmula:
4.8784
B  0.61
3.08  4.8784
Como A  0.13 se toma el valor indicado de 0.13
B  A  0.61 0.13  0.48
 Continuamos con la correlación de Hagendorn y Brown.
130
Si (B-A) es positivo o cero, se continúa con la correlación Hagedorn y Brown.

Si (B-A) es negativo, sigue con la correlación de Griffith, el procedimiento se
encuentra en el método de Orkiszewski.
22. Encontrar el número de diámetro de la tubería, Nd.
62.021
Nd  120.872  0.166   21.53
54
23. Calcular el colgamiento en función de 


 6.1520  0.10
764.7  0.0020 

           2.29X10^4
 9.787   14.7   21.53 

24. De la figura 3.10 se obtuvo:
HL
 0.2
ψ
25. Determinar la segunda corrección en función de  :
 6.1520  0.0097 0.380


 
  0.0024 
 21.532.14 
 

26. De la figura 3.11 se determinó:
ψ  1.00
27. Calcular el valor de HL:
HL  0.2 1.00  0.2
Para bajas viscosidades no habra corrección y ψ = 1.00
28. Determinar el valor para dos fases del numero de Reynolds, (NRe)TP.

(N Re ) TP 2.2  10 2  378585   9.74 x10 5
  1.995 

 12


3.4313
0.2

 0.018 1 0.2  

130


29. Determinar el valor para 
d
 0.00015
  0.000902
d 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
31. Calcular la densidad media de la mezcla  m por dos métodos.
(a) Utilizando el valor de HL del paso 27, calcular  m como sigue:
lb m
 m  62.0121
3
0.2  2.4810  1  0.2  14.3872
pie
32. Repitiendo los pasos 5, 7, 16, 17, y 19 para p1 y p2.
T1  120F T2  150F
Z1  0.9335 Z1  0.8976
BO1  1.0513 BO2  1.0974
R S1  73.0731 R S2  75.2569
VSL1  3.05 VSL2  3.11
VSG1  7.6572 VSG1  3.7278
33. Calcular la velocidad de la mezcla en dos fases para p1 y p2.
vm1  3.05  7.6572  10.7098

vm2  3.11  3.7278  6.8354
34. Determinar el valor para  vm2  

   
Δ vm2  10.7098 2  6.8354 2  67.9754
35. Calcular  h correspondiente a p  p1 - p2.
131
 67.9754 
144  1000 - 50014.3872  

 
 
h  64.4   3555 pies
0.022  37858525
14.3872   1.995 
2.9652 1011    14.3872
 12 


3.3.4 Método Orkiszewski.
Calcular la caída de presión a través de la TP con los siguientes datos:
Diámetro nominal, = 1.25 pg.

Diámetro de flujo = 1.38 pg
qo = 400 bpd
qw = 600 bpd
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. Se selecciona la presión de 500 psi como el punto de inicio (cabeza del pozo).
2. La temperatura en cada punto de presión es obtiene de que a 120 oF se tiene una p =

500 psi y a 150 oF se tiene p = 1,000 psi.
3.
Δp = 1,000 – 500 = 500 psi
4. No es necesario asumir un Δh debido a que la p2 es conocida.
5. La temperatura promedio del incremento es:
120  150
T  135 o F
2
6. Se calculan las propiedades de los fluidos a p y T con las correlaciones vistas en el

primer capítulo:
132
a)
p1  p2 500 1,000
p p  14.7  764.7 psia
atm
2 2
Bo = 1.0736
Rs = 107.223 pies3/bl@ c.s.
Z = 0.912
b) Los gastos de gas y líquido a condiciones de escurrimiento son:
3.27 107  Zq o R - R s T  460

q 
g
p

3.27 107 0.912400500 107.223135  460 pies3
qg   0.03646 
764.7 seg
q L 6.49 105 q  oB  oq  Bw w
q  6.49 105        pies3

L 400 1.0736 600 1.0 0.0668
seg
q T q g q  0.03646
L  0.0668  0.1033 pies3/seg
c) Los gastos de masa son:
wg 8.85 107  q  γ gR  R 

s
wg 8.85 107 1,000 0.65500 107.223  0.2259 lbm/seg
w  4.05 103 oq  γo  q w γ w 8.85 107  q  γ  R

L g s
w L  4.05 103 4000.9218  6001.07  8.85 107 1,0000.65107.223  4.1551 lbm /seg
wm  wg  w L  0.2259  4.1551  4.381lbm /seg
d) Las densidades son:
wg 0.2259  6.1958 lbm

g  
qg 0.03646 pie3
133
w L 4.1551 lb
   62.2021 m
L
qL 0.0668 pie3
e) Velocidades superficiales:
π d π1.38
2 2

Af  h
 1.4957 pg2
4 4
Donde Af = área de flujo de la TP.
144  q L 1440.0668 pies

v sL    6.4312
Ap 1.4957 seg
144  q g 1440.03646 pies

v sg    3.5102
Ap 1.4957 seg
pies
vm  vsL  vsg  6.4312  3.5102  9.9414
seg
f) La viscosidad del aceite (Beggs y Robinson), del agua y del líquido es:
μo = 7.7813 cp
μw = 0.5313 cp
 1   1  1.5 
μ μ  
 WOR  7.7813   0.5313  3.4313 cp
o  μ w     
2.5 
L
1 WOR 1 WOR 2.5
       

g) La tensión superficial del líquido es:

  
   1 
   WOR  30  1 
  1.5  dinas
o
  w     70   54
1  WOR
L 
1 WOR 2.5 2.5 cm

       

h) Los números de velocidad del líquido y gas son:

0.25
 62.2021
 L 
0.25
NLv  1.938  vsL    1.9386.4312   12.9121

 σ L   54 
0.25
 62.2021

0.25
Ngv  1.938  vsg  L

  σ L   1.9383.5102
134
  
54  7.047
5
135
7. 
 v 2   9.9414 2 
LB  1.071  2.6616   1.071   2.6616
m   189.5448
 dh   1.38 

Como el límite de LB debe ser mayor o igual a 0.13, utilizamos:
LB = 0.13
vsg 3.5102
  0.3530
vm 9.9414
LS  50  36  NLv  50  3612.9121  514.8356

vsg
Como  LB, Ngv  LS , el patrón de flujo es bache.
vm
a) Para determinar δ hacemos:
Calculamos el flujo fraccional de aceite o de agua:
1 1
fo    0.4
1  WOR 11.5
fw = 1 – fo =1-0.4 = 0.6
Por lo que el principal líquido que fluye en la tubería es agua (fase continua),
y como vm < 10 utilizamos la ecuación 3.88:
1.38  d 
δ  0.681  0.013  h 
d
 log μ L  0.232  log v m  0.428  log  
h
 12   12 

   
1.38
 1.38   1.38 
δ  0.681  0.013   log 3.4313  0.232 log 9.9414  0.428 log    1.4521

 12   12 
δ  0.0659.9414
Como δ  0.6462 , entonces el valor de δ es correcto.
b) Para determinar vb hacemos lo siguiente:
Calculamos el número de Reynolds:
136
124 dh  vm  L 1241.389.941462.2021

N Re    30,839
μL
L
3.4313
Ahora suponiendo un valor de vb =1.75 pies/seg
124  dh  vb  L 1241.381.7562.2021
N Re    5,429
μL
b
3.4313
Como NRe  6,000 y 3,000  NRe  8,000 tenemos que:

L b
v  0.5 α α2  13.59

0.5
 μ L 

 L d h /120.5 
b


   32.1741.38 
0.5
 6
 1.0012
α  0.251  8.74 10 30,839 
 12 


Entonces:
 
 13.593.4313
v b  0.51.0 1.02  0.5  1.9868
 

62.20211.38/12  
0.5
Como 1.75  1.9868, empezamos ahora suponiendo un valor de vb = 1.9868
124  dh  vb  L 1241.381.986862.2021
N Re    6,163.13
μL
b
3.4313
Como NRe  6,000 y

L
3,000  NRe  8,000 b
cae en el mismo intervalo
tenemos que:
vb = 1.9868 pies/seg
8. Calcule el gradiente de elevación y de fricción con:
C 3   L v sL  v b    g  v sg  62.20216.4312 1.9868  6.19583.5102   545.3658
 
 Δp  1 
 C3   Lδ  1 
 545.3658  62.20211.4521  psi
 e 0.9447 
   144  vm  vb  144 9.9414 1.9868  pie
ΔL
137
Con NRe =30,839 y con la rugosidad ε = 6 x 10-5 se obtiene el factor de fricción del
diagrama de Moody:
f = 0.02346
Y por lo tanto el gradiente de fricción es:
 Δp  12  f  v m    v  v 
2
     δ 
 Δh f
  L sL
2  12 g dh  v m  vb
b


 Δp  0.02346 9.9414  62.2021  6.4312 1.9868 

2
psi
 
 21232.1741.38   1.4521  0.2920 pie
Δh 9.9414 1.9868
 f  

9. 
Wm  vsg 4.3813.5102
Ek    2.8996 106
4,637  p  A p 4,637764.71.4957

  Δp 
 Δp   Δp 
     0.9447  0.2920 psi
 Δh Δh

 e  f   1.2367
   6
 Δh T 1 Ek 1 2.899610 pie
Entonces ΔL es: 

Δp 500
Δh    404.3018 pies
1.2367 1.2367
10. Como la p2 es conocida la Δh calculada es la correcta.
3.3.5 Método Beggs y Brill

Determinar la distancia (Δh) entre las presiones p1 y p2.
Datos:
d = 1.995 pg
Temperatura = T = 120 oF + 0.017 z
p1 = 500 psig
p2 = 1,000 psig
Presión atmosférica = patm = 14.7 psia
θ = 0o para vertical = 90o para horizontal
qo = 400 bpd
qw = 600 bpd
138
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:
p1  p2 500 1000
p p  14.7  764.7 psia
atm
2 2
Δh 2,000
hh 0  1,000 pies
2 2
2.
T  120  0.017  z  120  0.0171,000  137 o F
3. Las propiedades de los fluidos a p y T se obtuvieron de las correlaciones vistas en

el primer capítulo:
Rs = 106.8419 pies3/bl μw = 0.5214 cp

Bo = 1.0746 Z = 0.9131
μo =7.4445 cp Bw = 1.0
4. 
141.5 141.5
γo    0.9218
131.5  API 131.5  22
o
350  γ  0.0764  R s  γ g 
5.
o 3500.9218 0.0764106.84190.65 54.3529 lb
  
o
5.6146  B o 
m 5.61461.0746 pie3
350  γw 3501.07  66.7011 lbm

w  
5.6146  Bw 5.6146 1.0 pie3
   1  1.5  lb
   1 
 WOR  54.3529   66.7011  61.7618 m
o
    w     
1 WOR  1 WOR  2.5  2.5 
L 
pie3
       
139

0.0764  γ g  p  520 0.07640.65764.7520 lbm

    2.4643
g
14.7T  460Z 14.7137  4600.9131 pie3
6.
3.27 107 0.9131400500 106.8419137  460 pies3
qg   0.03666 
764.7 seg
  5
      pies3
q L 6.49 10 400 1.0746 600 1.0 0.0668
seg
7. 
π  d 2 π1.995 
2
Ap    3.1259 pg 2
4 4
144  q L 1440.0668 pies

v sL    3.0773
Ap 3.1259 seg
144  q g 1440.03666 pies

v sg    1.6888
Ap 3.1259 seg
pies
vm  vsL  vsg  3.0773  1.6888  4.7661
seg
8.
GL  L  vsL  61.76183.0773  190.0596 lbm /seg  pie2
G g  g  vsg  2.46431.6888  4.1617 lb m/seg  pie2
G T G L G g 190.0596  4.1617  194.2213 lb /seg

m  pie2
9.
qL 0.0668
λ   0.6457
q L  qg 0.0668  0.03666
10.
v
2
124.76612 124.7661
2
NFR m    4.2433

g  d/12 32.21.995 32.21.995
140
   1  1.5 
μμ   1 
 WOR  7.4445   0.5214  3.2906 cp
o
   μ w     
1  WOR 1  WOR 2.5 
L 
2.5
       
 
μ m  μL  λ  μ g1 λ  3.29060.6457 0.0181 0.6457  2.1311cp
  
   1 
   WOR  30  1 
  1.5  dinas
L o
  w     70   54
1 WOR 1 WOR 2.5 2.5 cm
       

194.22131.995
11.
M m  d/12 

  22,546.826
N 
μ
Re
m  6.72 104 
2.1311 6.72 104 12

0.25
 61.7618 
L 
0.25
NLv  1.938 vsL     1.9383.0773   6.1674

  54 
 L 

12.
L1  316 λ 0.302  3160.6457   276.8958
0.302
L2  0.0009252 λ 2.4684  0.00092520.6457 

2.4684
 2.7236 10 3
L3  0.10 λ -1.4516  0.100.6457 

1.4516
 0.1887
L4  0.5 λ -6.738  0.50.6457 

6.738
 9.5276
13. Dado que λ  0.4 y L3 < NFR  L4, el patrón de flujo es intermitente.
14. El colgamiento horizontal para flujo intermitente es:
 a  λb 0.8450.6457 0.5351
H L 0  N c   0.6521
FR 4.24330.0173
15. El coeficiente del factor de corrección por inclinación es:
C  1 0.6457 ln 2.960.6457   0.305

6.16740.4473 4.24330.0978   0.0989
16. Debido a que la tubería es vertical, el factor de corrección del colgamiento de
líquido debido a la inclinación es:
140
ψ  1  0.3 C  1 0.30.0989  1.0297
141
17. El colgamiento de líquido corregido para flujo vertical es:
H L 90   H L 0   ψ  0.65211.0297   0.6715
La densidad de la mezcla es:
m  L  H L  g 1 H L   61.76180.6715 2.46431 0.6715  42.2826

lbm
pie3
18.
λ 0.6457
y   1.432
HL 90 2
0.67152
ln1.432
 0.0523  3.182  ln1.432 0.8725ln1.432  0.01853 ln1.432
S  0.3672 2 4
fT
 eS  e0.3672  1.4437,
fns
19. 
1
f ns  2  0.0264
2 log 
 22,546.826   
4.5223 log 22,546.826 3.8215 
  

20. El factor de fricción de las dos fases es:
f
f T  fns f T  0.02641.4437  0.0381
ns
21.  
Δp    m  vm  vsg

144 1 g  p 14.7144

Δh  c 
g  m f  Mm  vm 12
senθ  T
gc 2 g c d
144   42.2826 4.7661 1.6888

1000 - 500 1  32.2764.7 144  
Δh     1,579.6 pies
   
42.2826 sen 90 0.0381 194.2213 4.766112


2 32.2 1.995
142


Como se puede observar el valor calculado de Δh difiere del supuesto en el primer

paso, por lo que se tiene que repetir el procedimiento iniciando con un nuevo valor
supuesto de Δh = 1,579 pies.
143

Flujo Multifasico en Tuberias Verticales

Cargado por

Copyright:

Formatos disponibles

Flujo Multifasico en Tuberias Verticales

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Flujo Multifasico en Tuberias Verticales

Cargado por

Copyright:

Formatos disponibles

Capítulo 3 Flujo multifásico en tuberías verticales

a. Diseñar las tuberías de producción y líneas de descarga.

Las diversas correlaciones existentes para el cálculo de distribuciones de presión

 Poettmann y Carpenter (1952)

Tabla 3.1. Criterios en el desarrollo de las correlaciones.

Grupo I Grupo II Grupo III

Patrones de flujo No se distinguen No se distinguen Sí se distinguen

 Tuberías de 2, 2 ½ y 3 pg.  RGL menores a 1,500 pies3/bl.

Fancher y Brown obtuvieron una correlación para tuberías de 2 pulgadas de

3. Determinar la masa de aceite, gas y agua por barril de aceite.

M  350  γo  0.0764  γg  R  350  γw  WOR (3.2)

4. Calcular el gasto másico de la mezcla wm en lbm/día.

5. Suponer una presión inicial en la TP.

6. Calcular el volumen de la mezcla VT a la presión y temperatura establecida:

V  5.61 B  5.61 WOR  R  R 14.7  T  Z  (3.3)

7. Determinar la densidad a la presión supuesta:

8. De la figura 3.1, determinar el gradiente de presión fluyendo para el gasto másico y

10. De los dos gradientes obtenidos, determinar un gradiente de presión promedio:

12. Repetir el procedimiento hasta completar la profundidad del pozo.

Donde wm es el calculado en el paso 4, en lbm/día

15. Calcular el gradiente de presión total con:

16. Por último calculamos la profundidad total con:

Figura 3.1. Correlación de la densidad del fluido (Poettmann).

Figura 3.2. Correlación de datos de campo en pozos fluyentes y de Bombeo Neumático.

3.2.2 Duns y Ros

 Región I: La fase líquida es continua y el flujo burbuja, flujo tapón y parte

1. Calcular la densidad relativa del aceite.

2. Obtener la masa asociada con un barril de líquido a condiciones estándar:

3. Determinar la densidad de la fase líquida:

5. Obtener la temperatura promedio:

7. Calcular la densidad promedio de la fase de gas:

9. Calcular la viscosidad promedio del agua con las ecuación 1.110.

10. Calcular la viscosidad de la mezcla del líquido:

11. Calcular la tensión superficial de la mezcla de líquido:

14. Calcular el área transversal de la tubería.

15. Calcular el número de viscosidad del líquido:

16. Calcular la velocidad superficial del líquido vsL (suponiendo Bw = 1.0):

17. Calcular la número de velocidad del líquido: 

18. Calcular la velocidad superficial del gas:

Ngv  1.938  vsg  L  (3.21)

20. Calcular el número de diámetro de la tubería: 

21. Seleccionar el régimen de flujo adecuado de la figura 3.3.

22. Determinar el factor de resbalamiento dependiendo de la región de flujo obtenida en

El resbalamiento se obtiene de la siguiente ecuación:

Donde F1, F2, F3, y F4 se obtienen de la figura 3.4 y:

Para flujo anular Nd se basa en el perímetro mojado, por lo tanto:

Los límites de la Región I son: 

Donde L1 y L2 se obtienen de la figura 3.5.

b) Para la Región II:

F5, F6 y F7 se obtienen de la figura 3.6 y donde:

F6'  0.029  N d  F6 (3.28)

Los límites de la Región II son: 

c) Para la Región III:

El límite de la Región III es: 

23. Determinar la velocidad de resbalamiento si el flujo se encuentra dentro de la región

24. Determinar el colgamiento de líquido:

25. Determinar el Número de Reynolds para el líquido:

26. Determinar el gradiente de fricción de acuerdo a la región de flujo.

a) Para la Región I y II:  