CONCEPTOS BASICOS

Elevación Energía cinética

Fricción
(Acceleración)

m tp m m
m
m

Factor Fricción: f( m m )
m

Al hablar de flujo multifásico se entiende la producción de aceite, agua y gas , por lo que las propiedades
que intervienen en la ecuación que forman la mezcla de flujo deberán ser la la mejor ponderación que refleje
de la mejor manera los aspectos físicos del fenómeno de flujo por tuberías.
 Fracción volumen del
liquido o colgamiento del
líquido
AL
YL = =HL
AT
Ag
Yg = =HG
AT AL

YL + Yg = 1

Colgamiento sin resbalamiento

qL qL
λL = =
q L + q g qT
Velocidad real

Con colgamiento sin resbalamiento

Con colgamiento
Densidad de la mezcla
Gasto másico
Viscosidad de la Mezcla
Burbuja Bache Bache Anular
Anular/ Niebla
Caotico
Mapas de transporte de Flujo
Diagrama plano que muestra regiones de existencia de diferentes patrones de flujo y la
transferencia de un patrón de flujo a otro.
Calculo de caídas de presión en
tuberías horizontales con flujo
multifasico
• Correlaciones:

1. Bertuzzi, Tek y Poettmann

2. Eaton, Andrews y Knowles
3. Beggs y Brill
4. Dukler
JAGP 1805-LA FUNCION DE PRODUCCION 12

9.1.1.4. Método de Dukler et al

El trabajo para flujo horizontal fue publicado en 1964. Ellos acumularon toda la data
publicada a la fecha sobre flujo bifásico horizontal, esta data consideraba datos de
laboratorio para tuberías de corta longitud y datos de campo correspondiente a
tuberías largas, más de 20000 datos experimentales fueron tomados desde 1959,
solo 2600 fueron consistentes. El trabajo consiste en el desarrollo de una correlación
basado en el concepto de similitud. Primero discute un análisis, donde no existe
deslizamiento entre las fases y flujo homogéneo fue considerado a existir. Luego su
análisis para flujo bifásico. El identificó que en el flujo del fluido actuaban solo cuatro
fuerzas principales, las cuales fueron identificadas como la fuerza debida a la
presión, la fuerza debida al esfuerzo de corte, la fuerza debida a la inercia o
aceleración del fluido y la fuerza debida a la gravedad. Dukler consideró la existencia
de similitud dinámica y cinemática en el flujo del fluido.

Sus métodos no consideran la existencia de patrones de flujo, primero se determina

las propiedades de las fases liquido y gas a condiciones de flujo. Luego usando las
correlaciones de Dukler para flujo bifásico se calcula el número de Reynolds y a
partir de este el factor de fricción es determinado para luego calcular las pérdidas de
fricción.

En 1969 fue desarrollado uno de los métodos más utilizado para amplias
condiciones de flujo en tubería horizontal, sirvió de base para el desarrollo del
Manual de Diseño de la AGA-API.

El procedimiento de cálculo para la caída de presión en flujo horizontal a partir del

conocimiento de la presión en el extremo agua arriba de una tubería.

1.- Asuma un valor para la presión en el extremo agua abajo de la tubería (Pf)sup.

2.- Calcule la presión promedio

12
JAGP 1805-LA FUNCION DE PRODUCCION 13

( PF ) sup + PI
Pprom = (9.1.18)
2

Donde : PI , PF representan valores de presión, psf

3.- Determine las propiedades de cada fase (ρL, ρG,µG, µL,σ), solubilidad Rs, factor
volumétrico Bo, viscosidad de las fases, factor de compresibilidad del gas Z, etc…

4.- Calcule las tasas de flujo volumétrico de líquido y gas

5.615QLSC Bo
QL = (9.1.19)
86400

QLSC (GOR − Rs )PSCTZ

QG = (9.1.20)
86400 PTSC

QL representa el caudal de liquido, cfs

QG representa el caudal de gas, cfs

QLSC representan el caudal de liquido a condiciones de tanque, bbl

Bo representa el factor volumétrico del crudo

GOR representan la relación gas petróleo, scf/sbl

Rs representa la solubilidad del gas , cf/sbl

P PSC representan valores de presión, psf

T TSC representan valores de temperatura, ºR

Z representa el factor de compresibilidad del gas, adimensional

5.- Calcular la fracción volumétrica de la fase líquida

QL
λL = (9.1.21)
QL + QG

13
JAGP 1805-LA FUNCION DE PRODUCCION 14

6.- Cálculo de la densidad de la fase líquida en la tubería

62.4γ L + 0.0764γ G RS / 5.615

ρL = (9.1.22)
Bo

donde: γ L γ g representa la gravedad especifica del gas, adimensional

9.- Cálculo de la densidad de la fase gaseosa en la tubería

0.0764γ G PTSC P
ρG = = 2.701γ G (9.1.23)
PSCTZ TZ

Donde :
P representa el valor promedio de la presión, psia
T representa el valor promedio de la temperatura, ºR

8.- Cálculo de la velocidad superficial de la mezcla en la tubería

QL + QG Q +Q
Vm = = 576 L 2 G (9.1.24)
A πd

Donde:
d representa el diámetro de la tubería, in

9.- Cálculo de la viscosidad de la mezcla en la tubería

µ m = λµ L + (1 − λ )µ G (9.1.25)

10.- Estime el holdup HL de la fase líquida en la tubería

11.- Cálculo de la densidad de la mezcla en la tubería

14
JAGP 1805-LA FUNCION DE PRODUCCION 15

ρm = ρL
λ2
+ ρG
(1− λ)
2
(9.1.26)
HL (1 − H L )

12.- Cálculo del número de Reynolds

ρ mVm D
Re = (9.1.27)
µm

13.- Con los valores de λ y Re determine el holdup de líquido haciendo uso de la

figura 9.3.

14.- Compare los valores de los holdup asumidos y calculados. Si la diferencia

absoluta entre ellos es mayor que una tolerancia permitida, repita los pasos 10 a 14,
en caso contrario continúe con el proceso

15.- Determine el factor de fricción de Darcy para un flujo monofásico. Existen varias
correlaciones para determinar el factor de fricción para tubería lisa o rugosa. Dukler
en su desarrollo recomendó la correlación siguiente

0.50
f n = 0.00560 + (9.1.28)
(Re )0.32

16.- Determine la relación entre el factor de frición bifásico y monofásico de la figura

fTP
9.4, .
fn

17.- Calcule el factor de fricción para la mezcla bifasica

fTP
fTP = fn (9.1.29)
fn

18.- Calcule la caída de presión debida a la fricción

15
JAGP 1805-LA FUNCION DE PRODUCCION 16

f LV ρ
2
∆Pf = TP m m (9.1.30)
2 gc D

Figura 9.3 Correlación de holdup de líquido (Dukler)

Dukler correlacionó el factor de fricción bifasico mediante la relación

fTP = Fβf D (9.1.31)

γ
F = 1+ (9.1.31a)
1.218 − 0.478γ + 0.444γ 2 − 0.094γ 3 + 0.00843γ 4

γ = − ln λL (9.11b)

ρ L λ2L ρ G (1 − λL )2
β= + (9.11c)
ρ ns H L ρ ns (1 − H L )

16
JAGP 1805-LA FUNCION DE PRODUCCION 17

Fig. 8.4 Curva normalizada para el factor de fricción

19.- Calcule la caída de presión debida a la energía cinética

1 ⎧⎪⎡ ρ Q 2 ρ Q 2 ⎤ ⎡ρ Q 2 ρ Q 2 ⎤ ⎫⎪
∆Pacc = ⎨⎢
G G
+ L L ⎥ −⎢ G G + L L ⎥ ⎬ (9.1.32)
⎪⎩⎢⎣ 1 − H L H L ⎥⎦ agua − abajo ⎢⎣ 1 − H L
2
144 g c H L ⎥⎦ agua − arriba ⎪
⎭

20.- Cálculo de la caída de presión debida al cambio de energía potencial

20a.- Cálculo de la velocidad superficial de la fase gaseosa en la tubería

576QG
VGSP = (9.1.33)
πd 2

20b.- Determine el factor de corrección HLF usando la figura 9.14.

20c.- Calcule las caída de presión correspondiente a la diferencia de energía

potencial

17
JAGP 1805-LA FUNCION DE PRODUCCION 18

ρ L H LF ∑ ∆h
∆Pelev = (9.1.34)
144

21- Calcule la caída de presión total

∆P = ∆Pf + ∆Pacc + ∆Pelev (9.1.35)

22- Calcule la presión en el extremo final (Pf)calc de la tubería o segmento de la

misma.

23- Compare los valores la presión en el extremo final (Pf)calc y (Pf)sup.

a.- Si la diferencia es mayor que el valor de la tolerancia permitida, repita los
pasos 1 al 23.
b.- En caso contrario finalice los cálculos o continúe con otro segmento de la
tubería.

9.1.1.5. Método de Beggs y Brill

Esta correlación fue desarrollada considerando la inclinación de la tubería

comprendida entre 0 y 90° a – 90°. La data experimental fue generada en un circuito
de prueba de 1 y 1.5 pulgadas utilizando aire y agua como fluido de trabajo. Tres
regímenes de flujo fueron considerados, correlaciones para el holdup y el factor de
fricción correspondiente a cada patrón de flujo fueron desarrolladas. Los patrones de
flujo utilizados corresponden a los patrones de flujo observado en tuberías
horizontales. El holdup se determina primero para flujo horizontal y luego es
corregido para el ángulo de inclinación. Los patrones de flujo usados por esta
correlación se muestran en la siguiente figura 9.6.

Los parámetros estudiados y su intervalo fueron:

1- Tasa de gas entre 0 a 300 MMscfd
2- Tasa de liquido entre 0 a 30 gal/min
3- Presión promedio entre 35 a 95 psia

18
JAGP 1805-LA FUNCION DE PRODUCCION 19

4- Diámetro de la tubería 1 y 1.5 in

5- Holdup de liquido entre 0 y 0.870
6- Gradiente de presión entre 0 y 0.8 psi/ft
7- Angulo de inclinación entre -90 y 90º
8- Fluidos utilizados agua y aire

Fig. 9.6 Patrones de flujo tubería horizontal (Beggs-Brill)

Para cada diámetro de la tubería, las tasas de líquido y gas fueron variadas, tal que
todos los patrones de flujo fueron observados. Luego un conjunto particular de tasas
de flujo fueron seleccionadas y el ángulo de la tubería fue variado dentro del
intervalo correspondiente a la inclinación de la tubería tal que el efecto de la
inclinación de la tubería sobre el holdup y gradiente de presión fuese observado. El
holdup y gradiente de presión fueron medidos para la posición horizontal ángulos
más y menos de 5, 10, 15, 20, 35, 55, 75 y 90º. Las correlaciones fueron
desarrolladas para resultados correspondientes a 584 pruebas.

Diferentes correlaciones fueron presentadas para los tres patrones de flujo. El

holdup de líquido, el cual existiría si la tubería fuese horizontal es calculado y luego
es corregido para la inclinación de la tubería. El holdup fue encontrado a tener un

19
JAGP 1805-LA FUNCION DE PRODUCCION 20

máximo en + 50º de la horizontal y un mínimo a – 50º. El mapa original de los

patrones de flujo fue modificado para incluir una zona de transición entre los
regímenes de flujo segregado e intermitente. La superposición de los mapas de los
patrones de flujo se presenta en la figura 9.7. Un factor de fricción bifásico es
calculado usando ecuaciones, las cuales son independiente del régimen de flujo
pero dependiente del holdup. Un grafico de un factor de fricción normalizado como
una función del holdup y del contenido de liquido fue elaborado, como se muestra en
la figura 9.7 a.

Las variables siguientes N FR número de Fraude y λ fracción volumétrica de la fase

líquida son usadas para determinar el régimen de flujo si la tubería tuviese en una
posición horizontal. Este régimen de flujo es solo un parámetro correlaciónate y no
da ninguna información sobre el régimen de flujo actual excepto que la tubería fuese
horizontal.

El procedimiento de cálculo para la caída de presión en flujo horizontal a partir del

conocimiento de la presión en el extremo agua arriba de una tubería.

1.- Suponga un valor para la presión en el extremo agua abajo de la tubería (Pf)sup.

2.- Calcule la presión promedio

(Pf )sup + PI
Pprom = (9.1.36)
2

3.- Determine las propiedades de cada fase (ρL, ρG,µG, µL,σ), solubilidad Rs, factor
volumétrico Bo, viscosidad de las fases, factor de compresibilidad del gas Z, etc…

En la identificación del régimen de flujo requiere del cálculo de varios números

adimensionales, incluyendo el número de Froude correspondiente al flujo bifásico.
Esta correlación es utilizada para cualquier ángulo de inclinación.

20
JAGP 1805-LA FUNCION DE PRODUCCION 21

Fig. 9.7 Mapa de patrones de flujo en tubería horizontal (Beggs y Brill)

Fig. 9.7a Factor de fricción bifásico (Beggs y Brill)

21
JAGP 1805-LA FUNCION DE PRODUCCION 22

4.- Cálculo de la densidad de la fase líquida en la tubería.

62.4 * 5.615γ L + 0.0764γ G RS

ρL = (9.1.37)
5.615 Bo

5.- Cálculo de la densidad de la fase gaseosa en la tubería.

0.0764γ G PTSC P
ρG = = 2.701γ G (9.1.38)
PSCTZ TZ

6.- Calcule las tasas de flujo volumétrico de líquido y gas en los gasoductos.

5.615QLSC Bo
QL = (9.1.39)
86400

QLSC (GOR − Rs )PSCTZ

QG = (9.1.40)
86400 PTSC

7.- Calcular la fracción volumétrica de la fase líquida

QL
λL = (9.1.41)
QL + QG

8.- Calcule las velocidades superficiales para cada fase y la mezcla

QL
VSL = (9.1.42)
AT

QG
VSG = (9.1.43)
AT

QM QL + QG
VSM = = (9.1.44)
AT AT

22
JAGP 1805-LA FUNCION DE PRODUCCION 23

9.- Determine el patrón de flujo

9a.- Calcule los números adimensionales

2
VSM
N FR = (9.1.45)
gD

0.25
⎛ρ ⎞
N LV = 1.938VSL ⎜⎜ L ⎟⎟ (9.1.46)
⎝σ L ⎠
Donde:
VSL representa la velocidad superficial de la fase liquida, ft/seg

ρ L representa la densidad de la fase liquida, lbm/ft3

σ L representa la tensión superficial de la fase liquida, dina/cm

9b.- Calcule los contornos para los patrones de flujo

L1 = 316λL
0.302
(9.1.47)

−2.4684
L2 = 0.0009252λL (9.1.48)

−1.4516
L3 = 0.10λL (9.1.49)

−6.738
L4 = 0.5λL (9.1.50)

Contornos de los patrones de flujo

Segregado

λL < 0.01 N FR < L! (9.1.51)

λL > 0.01 N FR < L2 (9.1.52)

23
JAGP 1805-LA FUNCION DE PRODUCCION 24

Transición

λL ≥ 0.01 L2 < N FR < L3 (9.1.53)

Intermitente

0.01 ≤ λL < 0.4 L3 < N FR ≤ L! (9.1.54)

λ L ≥ 0 .4 L3 < N FR ≤ L4 (9.1.55)

Distribuido

λL < 0.4 N FR ≥ L! (9.1.56)

λL ≥ 0.4 N FR > L4 (9.1.57)

Cuando el patrón de flujo se ubica en la región de transición, el holdup de líquido

debe ser calculado usando las ecuaciones de flujo segregado e intermitente,
interpole utilizando la siguiente relación:

H L (transición ) = AH L ( segregado ) + BH L (int ermitente) (9.1.58)

L3 − N FR
A= (9.1.59)
L3 − L2

B =1− A (9.1.60)

9c.- Determine el patrón de flujo basándose en los valores del número de Froude y la
fracción volumétrica de la fase líquida

24
JAGP 1805-LA FUNCION DE PRODUCCION 25

10.- Calcule el holdup de líquido una vez conocido el patrón de flujo haciendo uso de
las relaciones

H L (θ ) = H L (0 )ϕ (9.1.61)

aλL
b
H L (0) = c
(9.1.62)
N FR

donde:

Patrón de flujo a b c
Segregado 0.98 0.4846 0.0868
Intermitente 0.845 0.5351 0.0173
Distribuidos 1.065 0.5824 0.0609

El holdup debe satisfacer la siguiente condición

H L (0) ≥ λL (9.1.63)

El factor de corrección ϕ para el holdup, que permite considerar el efecto de la

inclinación de la tubería está dado por

ϕ = 1 + C [sen(1.8θ ) − 0.333sen3 (1.8θ )] (9.1.64)

con
αλL e N LV f
C = (1 − λL ) ln( g
) (9.1.65)
N FR

donde:

Patrón de flujo α e f g
Segregado ascendente 0.011 -3.768 3.539 -1.614

25
JAGP 1805-LA FUNCION DE PRODUCCION 26

Intermitente ascendente 2.96 0.305 -0.4473 0.0978

Distribuido ascendente ninguna corrección C=0, ϕ=1 HL=f(φ)
Patrones descendente 4.70 -0.3692 0.1244 -0.5058

11.- Cálculo de la densidad y la viscosidad de la mezcla

ρ s = ρ L H L + ρG H G (9.1.66)

ρ n = ρ L λL + ρ G λG (9.1.67)

µ n = µ L λL + µG λG (9.1.68)

12.- Cálculo del número de Reynolds

ρ nVm D
Re = (9.1.69)
µn

13.- Determine el factor de fricción de Darcy para un flujo monofásico. Existen varias
correlaciones para determinar el factor de fricción para tubería lisa o rugosa.
Originalmente Beggs y Brill recomendaron la correlación siguiente:

−2
⎡ ⎛ Re ⎞⎤
f n = ⎢4 log⎜⎜ ⎟⎟⎥ (9.1.70)
⎣ ⎝ 4.5223 log Re− 3.8215 ⎠⎦

14.- Determine el factor de frición bifásico

fTP
fTP = fn (9.1.71)
fn

El factor de fricción originalmente recomendado por Beggs y Brill es la que

representa el factor de fricción para una tubería lisa, razón por la cual es
independiente de la rugosidad de la tubería. Se recomienda que el factor de fricción

26
JAGP 1805-LA FUNCION DE PRODUCCION 27

incluya este efecto, de allí que pueda ser calculado a partir de la correlacion de
Colebrook o determinado a partir del diagrama de Moody.

1 ⎛ ε 2.51 ⎞⎟
= −2 log⎜ + (9.1.72)
⎜ 3.7 D R f ⎟
f ⎝ e ⎠

fTP
= es (9.1.73)
fn

λL
Y= (9.1.74)
[H L (θ )]2

ln(Y )
S=
{
− 0.0523 + 3.182 ln(Y ) − 0.8725[ln(Y )] + 0.01853[ln(Y )]
2 4
} (9.1.75)

Con la condición:

1 < Y < 1 .2 S = ln(2.2Y − 1.2) (9.1.76)

15.- Calcule el gradiente de presión debida a la fricción

f TPVm ρ n
2
⎛ dP ⎞
⎜ ⎟ = (9.1.77)
⎝ dL ⎠ f 2gc D

16.- Calcule el gradiente de presión debida a la elevación

⎛ dP ⎞ g
⎜ ⎟ = ρ s senφ (9.1.78)
⎝ dL ⎠ elev g c

17.- Calcule el gradiente de presión debida a la energía cinética

⎛ dP ⎞ ρ sVmVsg ⎛ dP ⎞
⎜ ⎟ = ⎜ ⎟ (9.1.79)
⎝ dL ⎠ acc g c P ⎝ dL ⎠ T

27
JAGP 1805-LA FUNCION DE PRODUCCION 28

ρ sVmVsg
EK = (9.1.80)
gc P

20.- Cálculo de la caída de presión total

⎛ dP ⎞ ⎛ dP ⎞
⎜ ⎟ +⎜ ⎟
⎛ dP ⎞ ⎝ dL ⎠ elev ⎝ dL ⎠ f
⎜ ⎟ = (9.1.81)
⎝ dL ⎠ T 1 − Ek

⎛ dP ⎞
∆PT = ⎜ ⎟ L (9.1.82)
⎝ dL ⎠ T

21- Calcule la presión en el extremo final (Pf)calc de la tubería o segmento de la

misma.

22- Compare los valores la presión en el extremo final (Pf)calc y (Pf)sup.

a- Si la diferencia es mayor que el valor de la tolerancia permitida, repita los
pasos 1 al 22.
b- En caso contrario finalice los cálculos o continúe con otro segmento de la
tubería.

9.1.1.6 Método de Oliemans

Oliemans propuso un nuevo esquema para el análisis de flujo de gas en presencia

de líquido. Estableció que cierta cantidad de líquido una vez formado se mantenía
estático en la tubería, esta tesis está dada por la diferencia entre HL y λL. A partir de
esta consideración el área, el diámetro de flujo en la tubería fue redefinido. La
cantidad de líquido reduce el área efectiva de flujo.

A f = [1 − (H L − λL )]Ap
(9.1.83)

28