TRANPORTE DE GAS NATURAL – IGPUBA 2020

Ing. Gustavo Cavallo

FLUJO BIFASICO/MULTIFASICO

Este tipo de flujo se presenta con frecuencia en las líneas de yacimientos de petróleo y
se trata del flujo de líquidos hidrocarburos “revueltos” con vapores o gases. Este tipo
de flujo suele establecerse en ductos y gasoductos de captación.

En el transporte de las mezclas, el gas tiende a ubicarse en el centro y el líquido a

mojar las paredes, pero existen varios grados de dispersión de cada fase en la otra, lo
que origina distintos tipos de flujos.

Cabe destacar que Influye en el tipo de régimen de flujo:

• La cantidad relativa de las fases y sus velocidades (relativas y absolutas)

son de gran importancia.
• El tamaño y la inclinación de los tramos de la línea
• La presión y la temperatura del fluido en cada tramo.
• En flujos distribuidos la precisión de los cálculos puede llegar del 5 al
10% y en otros regimenes entre 10 y 15%.

Pueden hacerse tres grandes divisiones de los regimenes/tipos básicos del flujo
bifásico:

– a) Distribuidos: una sola fase “continua” y la otra distribuida (burbujas,

niebla y espuma).

– b) Segregados: dos fases son continuas con la fase gaseosa moviéndose

mas rápido (estratos, ondas, anular).

– c) Intermitentes: la continuidad de ambas fases están interrumpidas y

una serie de fuerzas inerciales y cinéticas se modifican (tapón, olas).

Regimenes de Flujo

1
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo

La caída de presión total se expresa en general por:

(Δp/ΔL)tot=C[(Δp/ΔL)elev +(Δp/ΔL)fric +(Δp/ΔL)acel ]

• En flujos “distribuidos” el ΔP es entre el fluido compuesto y el caño (pero no se

puede tratar como un flujo de una sola fase).
• En flujos “segregados” o “intermitentes” el ΔP es el resultado de la fricción con
el caño y la interacción entre las fases.

Todas son correlaciones matemáticas de resolución del flujo (calculo de la perdida de

carga) son semi-empíricas. El problema es determinar el tipo de flujo para un
determinado “set” de condiciones y aplicar la correlación correcta.

2
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo
Flujo Horizontal* (No hay líneas realmente horizontales, pero se entiende que el largo
es mucho mayor que cualquier cambio de elevación)

El tipo de régimen de flujo depende de:

– Caudales de gas y de líquido

– Relación gas-líquido
– Propiedades físicas del gas y del líquido
– Diámetro de la línea
– Energías de interfases (σ: tensión superficial) y fuerzas de corte (μ:
viscosidades) entre las fases presentes.

Método de LOCKHART y MARTINELLI (1949)

Es un método simplificado que propone cuatro combinaciones de flujos “viscoso” y

“turbulento” y se basa en dos hipótesis:

– Caída de presión igual para el gas y el liquido:

ΔPG = ΔPL
– Volumen total del caño:
Ѵc = ѴG + ѴL

El procedimiento propuesto por L&M fue:

– Calcular el ΔP de cada fase como si fluyeran separadas (dentro de la

cañería como flujo único).
– Determinar la correlación X = (ΔPG / ΔPL)0,5
– Calcular el Re del líquido y del Gas como si fueran flujos individuales y
utilizar la tabla y el gráfico.
– Del gráfico de L&M, conocido X calcular ØL y ØG.
– Estimar ΔPTP = زL ΔPL = زG ΔPG

Gráfico de Lockhart y Martinelli

3
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo

Tabla de Lockhart y Martinelli

Numero de Reynolds
Tipo de Flujo
Liquido φL Gas φG

Liq.Turb.-Gas
> 2000 Ltt > 2000 Gtt
Turb.

Liq.Viscoso-
< 1000 Lvt > 2000 Gvt
Gas Turb.

Liq.Turb.-Gas
> 2000 Ltv < 1000 Gtv
Viscoso

Liq.Viscoso-
< 1000 Lvv < 1000 Gvv
Gas Viscoso

4
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo
El paso inicial (calculo de ∆P involucra considerar en el diámetro D como si cada fluido
fuera por una cañería independiente). Este método se aplica con éxito a líneas
absolutamente horizontales conteniendo dos fases continuas.

Método de BAKER (1958)

Este investigador introdujo al método de L&M diferentes regímenes (mapa de flujos) y

una corrección por flujo inclinado.

En el gráfico, siendo:

– G: velocidad de masa del gas (área total del caño) en lb/hr/ft.²

– L: velocidad de masa del líquido (área total del caño) en lb/hr/ft.²

λ = [(ρG/0,075)(ρL/62,3)]0,5

Ψ=(73/σ)[μL(62,3/ρL)]1/3

Siendo:

σ : tensión superficial del líquido (dynas/cm)

μL :viscosidad abs. del liquido (ctp)

Para cada tipo de régimen fueron propuestas ecuaciones cuando líquido y gas están en
flujo turbulento.

5
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo
Método de BEGGS & BRILL (1973)

El procedimiento básico incluye:

a) Definir propiedades físicas y términos en tramo Δx.

b) Determinar las ecuaciones del mapa de régimen y calcular el HL.
c) Corregir el HL por inclinación.
d) Determinar el factor de fricción de 2 fases.
e) Resolver la ecuación ΔP/ΔL del Tramo.

• a) cálculo de propiedades a p y T promedios del tramo:

ρG y ρL : densidades del gas y del liquido

VSL = qL/A velocidad superficial del liquido
qL : caudal de volumen del liquido
VSG = qG/A velocidad superficial del gas
qG : caudal de volumen del gas
Vm = (qL + qG )/A velocidad de la mezcla
μG y μL : viscosidades del gas y del liquido
VSL = qL/A velocidad superficial del liquido

6
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo
Gm = GL + GG caudal de masa de la mezcla
σ = tensión superficial del liquido
qG : caudal de volumen del gas
λ = qL / (qL + qG ) velocidad de la mezcla

• b) Determine el tipo de flujo con el gráfico de B&B siendo los limites del
régimen horizontal:

λ < 0,1 y NFR < L1 Segregado

λ ≥ 0,01 y NFR < L2 Segregado
λ ≥ 0,01 y L2 ≤ NFR ≤ L3 Transición
0,01 ≤ λ < 0,4 y L3 < NFR ≤ L1 Intermitente
λ ≥ 0,4 y L3 < NFR ≤ L4 Intermitente
λ < 0,4 y NFR ≥ L1 Distribuido
λ ≥ 0,4 y NFR > L4 Distribuido

• b) Determine el régimen de flujo con el grafico (las curvas llenas indican el

trabajo original, las punteadas una modificación posterior)

Las ecuaciones de las líneas punteadas son:

L1= 316 λ0,302 L2= 0,0009252 λ-2,4684

L3= 0,10 λ-1,4516 L4= 0,52 λ-6,738

NFR= Vm2 /gD λ = VSL / Vm

Gráfico de B&B (revisado)

7
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo

• c) Según sea el régimen que resulte, calcule HL para caño horizontal con:

SEGREGADO

HL= 0,98 λ0,4546 /NFR0,0868

Corrección por Inclinación:

• Hacia arriba: C+= (1-λ) ln[0,011NVL3,539/λ3,768 NFR1,614]
• Hacia abajo: C-= (1-λ) ln[4,7NVL0,1244/λ0,3692 NFR0,5056]

INTERMITENTE

HL= 0,845 λ0,5651 /NFR0,0173

Corrección por Inclinación:

• Hacia arriba: C+= (1-λ)ln[2,46 λ0,305 NFR0,0978 /NVL0,4473]
• Hacia abajo: C-= (1-λ) ln[4,7NVL0,1244/λ0,3692 NFR0,5056]

DISTRIBUIDO

8
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo
HL= 1,065 λ0,5824 /NFR0,0609

Corrección por Inclinación:

• Hacia arriba: C+= sin correción

• Hacia abajo: C-= (1-λ) ln[4,7NVL0,1244/λ0,3692 NFR0,5056]

Nota: una vez calculado el HL y los pasos d) y e); se calculan los cambios causados por
la inclinación usando el C+ o C- combinado con el Ψ.

• d) para aplicar las formulas de c) se tiene:

NVL= 0,1778 VSL [ρL /g.σ]0,25

donde:
• NVL: grupo adimensional
• VSL [m/s]
• ρL [Kg/m3]
• σ [N/m]
• g = 9,81 m/s2
• d) para aplicar las formulas de c) en el caso de tener régimen de
Transición se interpola

HL= I₁ HL (segregado) + I₂ HL (intermitente)

donde

I₁ = L₃ - NFR /L₃ - L₂ I₂ = 1 - I₁

• e) Influencia de la Inclinación: Calcule HL para la sección inclinada:

HL(θ)= HL(0).Ψ

Siendo
Ψ=[1+ C (sen 1.8 θ - ⅓sen³ 1.8 θ]

θ : ángulo real desde la horizontal del tramo

C : factor calculado con las formulas del punto c

• f) Calcule el factor de fricción fncon la ecuación*:

fn=[2 log(Ren/(4,5223 log Ren-3,8215)]-2

Siendo
Ren = D.Vm.ρn / µn
fn : factor de fricción (sin resbalar)

9
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo
* Nota: el factor fn se puede encontrar asimismo del diagrama de Moody o formula de
Churchill

• g) Calcule el factor de fricción del tramo ftp con la ecuación:

ftp = fn . Es

Siendo:

S= ln y.[-0,0523+3,182 ln y-0,8725(lny)²+0,01853 (lny)⁴]¯¹

y = λ . [HL(θ) ]¯²

• h) Calcule la perdida de presión por unidad de longitud ΔP/ΔL para el tramo,

con:

f tp . n .Vm2
g.sen  L .H L  G 1  H L  
P 2D

L  .H  G 1  H L .Vm .VSG
1 L L
P

• Haga la sumatoria de las caídas de presión para todos los tramos para conocer
el ΔP total.

• i) Conclusiones:

• El procedimiento puede usarse para cada tramo Δx (en subida o en

bajada) y luego sumar los resultados para obtener el ΔP total.
• El trabajo de B&B se realizó sobre caños cortos y de pequeño diámetro,
usando agua y aire.
• B&B usa una carga de liquido constante. En caso de una línea de gas
natural y condensado, que incrementa el líquido a medida que se enfría,
el comportamiento de fase del sistema debe incluirse en el cálculo.

10
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo

EJEMPLO DE FLUJO BIFASICO (Metodologia Beggs & Brill)

Se trata de transportar por una cañería de diámetro D = 20” (e = 8,74 mm) un caudal
de petróleo qL= 3000 m3/dia en conjunto con un caudal de gas natural qG = 9 x 106
Sm3/dia de L = 100 km, partiendo de una presión P 1 = 100 bars (abs) y con presión de
entrega P2 = 40 bars (abs). La cañería recorre en subida 60 km llegando hasta 1000
metros de altura y luego sigue en bajada los siguientes 40 km hasta llegar a 500 metros
en el punto de entrega.

Las características de los fluidos son:

ρL = 885 kg/m3 μL = 250 ctp

ρG = 0,62 (rel. al aire) μG = 0,011 ctp
σ = 24.10-3 N/m (tensión superficial petróleo-gas natural)

Cálculos

= 67,25 Kg/m3

Método de Beggs & Brill

Se resolverá por tramos ΔL ≤ 10% L. En el presente ejemplo se elegirán tramos de 10

km

11
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo

Tipo de Flujo (según grafico B&B)

Entrando al Grafico de B&B (rectificado) para λ y Fr se determina el de tipo de flujo

para el tramo que resulta en zona de Transición.

Comprobación matemática:

Resulta λ > 0,1 y L2 < Fr < L3 lo que comprueba que el régimen del flujo en el primer
tramo es Transición

Calculo de HL (“Liquid Hold up”) para caño horizontal y flujo de transición

12
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo

Calculo de HL para tramo de cañería inclinada

Resulta:

Calculo del factor de fricción

Utilizando el diagrama de Moody fn = 0,026

Siendo el factor corregido para flujo bifásico

El exponente s se calcula con

13
 TRANPORTE DE GAS NATURAL – IGPUBA 2020
Ing. Gustavo Cavallo

(se observa que el factor de fricción aumenta)

Cálculo de la perdida de presión del tramo

Para el tramo de 10 km el ΔP será de aprox. 3,4 bars de modo que el tramo siguiente
se deberá comenzar con 96,6 bars y así continuar.

14