PRODUCCION I I PET - 208

METODO DE BEGGS AND BRILL

La correlación de Begss and Brill fue desarrollada de datos

experimentales obtenidos en una pequeña escala para la facilidad de las
pruebas. La facilidad consistió de secciones de tubería de 1” y 1.5” de
90 ft de longitud y el material fue de acrílico. La tubería podía ser
inclinada a cualquier ángulo. Los parámetros estudiados y los rangos de
variación fueron:

1. Caudal de flujo de gas de 0 – 30 MPCD

2. Caudal de flujo de liquido de 0 – 30 gpm
3. Presión promedio del sistema de 35 – 95 pisa
4. Diámetro de tubería 1 – 1.5 “
5. Liquido de hold up de 0 – 0.870
6. Gradiente de presión de 0 – 0.8 psi/ft
7. Angulo de inclinación de -90 - 90 º
8. Modelo de flujo horizontal

Los fluidos utilizados fueron aire y agua. Para cada diámetro de tubería.
Los caudales de líquido y gas fueron variando para todos los modelos de
flujo observado con la tubería horizontal.

Después de un rango de caudales de flujo fue seleccionado, el ángulo de

la tubería fue variando así que el efecto del ángulo de Hol Up y el
gradiente de presión fueron observados. El liquido Hola Up y el
gradiente de presión fueron medidos con ángulos respecto a la
horizontal de 5, 10, 15, 20, 35, 55, 75 y 90 º.

Las correlaciones fueron desarrolladas de 584 pruebas medidas.

Diferentes correlaciones para líquidos Hold Up son presentados para
cada uno de los tres regimenes de flujo “horizontal” Hold Up el cual
podrá existir “si la tubería fuese horizontal” es primero calculado y luego
para cualquier ángulo de inclinación de la tubería. Existiendo gráficos de
modelo de flujo horizontal.

Método de Beggs and Brill

1
PRODUCCION I I PET - 208

El máximo Hold Up es encontrado aproximadamente a 50 º respecto a la

horizontal y un mínimo aproximadamente a -50º.

DETERMINACION DE REGIMENES DE FLUJO

Las siguientes variables son utilizadas para determinar cual régimen de

flujo podrá existir si la tubería esta en una posición horizontal.
Estos regimenes de flujo es una correlación de parámetros q no
proporciona información acerca de los regimenes de flujo, a menos q la
tubería este horizontal.

Método de Beggs and Brill

2
PRODUCCION I I PET - 208

FLUJO SEGREGADO (ZONA UNO)

Caracterizado por la estratificación de la fase liquida y gaseosa q

ocasiona un alto resbalamiento, o retraso del líquido frente del gas.

Los límites del flujo son:

Segregado
Limites λ L < 0.01 y N FR < L1
O
λ L > = 0.01 y N FR < L2

FLUJO INTERMITENTE (ZONA DOS)

Caracterizado por la formación de bolsones de gas, de tal forma q se

alternan en la tubería estos bolsones de gas en la masa de líquido.

Los límites del flujo son:

Intermitente
Limites 0.01 < = L < 0.4 y L3 < N FR < = L1
Método de Beggs and Brill

3
PRODUCCION I I PET - 208

o
λ L > = 0.4 y L3 < N FR <= L4

FLUJO DISTRIBUIDO (ZONA TRES)

Este flujo esta dividido en dos zonas:

ZONA 3 A
Caracterizada por el flujo de niebla donde gotas de liquido están
distribuidas en la masa de gas.
ZONA 3 B
Caracterizada por el flujo de burbuja, donde burbujas de gas están
dispersas en la masa de liquido.
Los límites del flujo son:
Distribuido
Limites λ L < 0.4 y N FR >= L1
o
λ L > = 0.4 N FR y > L4

FLUJO DE TRANSICION (ZONA CUATRO)

Método de Beggs and Brill

4
PRODUCCION I I PET - 208

Que es justamente la transición del flujo segregado al intermitente, y

esta caracterizada por la inestabilidad del flujo.

Los límites del flujo son:

Trasmisión
Limites λ L > = 0.01 y L2 < N FR + < L3

Cuando el flujo cae en la zona de régimen de transición el líquido Hold

Up podrá ser calculado usando dos ecuaciones, del segregado e
intermitente e interpolando con los siguientes factores.

HL (Transición) = A * HL (segregado) + B * HL (intermitente)

Donde: A = (L3 - N FR ) / (L3 – L2)

B = 1- A
Las mismas ecuaciones son utilizadas para calcular el líquido Hold Up
para todos los regimenes de flujo. Los coeficientes y exponentes
utilizados en la ecuación son diferentes para cada régimen de flujo.

DENSIDAD BIFASICO

HL (0) = HL (0) * ψ
Donde:

HL (0) es el Hold Up a las condiciones de una tubería horizontal. Esto es

calculado con la siguiente formula:
HL (0) = (a * (λ L)b )/ NFR c

Donde:

Método de Beggs and Brill

5
PRODUCCION I I PET - 208

a , b , c son determinados para cada modelo de flujo en la tabla con la

restricción que HL (0) > = λ L .

MODELO DE a b c
FLUJO
Segregado 0.98 0.4846 0.0868
Intermitente 0.845 0.5351 0.0173
distribuido 1.065 0.5824 0.0609

El factor para corregir el Hold Up por efectos de la inclinación de la

tubería es la siguiente.

Ψ = 1+ C (sen(1.8*Ф)-0.333*sen3(1.8* Ф))

Donde: Ф es el Angulo de la tubería respecto a la horizontal. Para flujo

vertical ascendente, el ángulo de 90º y Ψ llegara a ser:

Ψ = 1+ 0.3 * C

Donde: C = (1- λ L)* Ln (d* λ Lc * (NLV )F (NFR )g

Donde: d, e, f y g son determinados para cada condición de flujo de

la sigt. Tabla con la restricción q C>= 0

MODELO DE d e f g
FLUJO
Segregado 0.011 -3.768 3.539 -1.614
ascendente
Intermitente 2.96 0.305 -0.4473 0.0978
ascendente
Distribuido - - - -
ascendente

Distribuido ascendente : C= 0. Ψ = 1, HL distinto de f(Ф)

FACTOR DE FRICCION

 dp  f tp   o  Vm 2
  
 dz  f 2 * gc * d

Método de Beggs and Brill

6
PRODUCCION I I PET - 208

donde:  o   L*  L   g *  g

f tp
f tp  f o *
fo

El factor de fricción de no resbalamiento es determinado de la curva de

la tubería lisa del diagrama de Mordí o de:

1
fo  2
  Nre  
 2 * log  
 
  4.2523 * log( Nre  3.8215)  

utilizando el siguiente numero de Reynolds:

 M * Vm* d
N Re 
o
donde:  o   L * LL   g * L g

la razón de los factores de fricción bifásico y de no desplazamiento es

calculada de la siguiente manera:

f tp
 eS
fo

donde:

Lny
S
 0.0523  3.182 * Lny  0.8725 * ( Lny ) 2  0.01853 * ( Lny ) 4

L
y
Y H    L
2

El valor de S llegaría a ser un punto en el intervalo 1 < y <1.2; y para

este intervalo, la función S es calculada por la ecuación:
Método de Beggs and Brill

7
PRODUCCION I I PET - 208

S  Ln 2.2 y  1.2 

TERMINO DE LA ACELERACION

Aunque el gradiente de presión de aceleración es muy pequeño ecepto

para alta velocidad de flujo, el término podrá ser incluido si es posible.

 dp   S * Vm * Vsg  dp 
   * 
 dz  acc gc * P  dz 

si definimos en términos de la aceleración como:

s * Vm * Vsg
Ek 
gc * P

El gradiente de presión total puede ser calculado con la siguiente

ecuación.

 dp   dp 
   
 dp   dz  cl  dL  f
  
 dz  total 1  E K 

donde:

S  L H L  g H g

 dp   g
 dz    * S
cl  gc 

Método de Beggs and Brill

8
PRODUCCION I I PET - 208

Método de Beggs and Brill

Ejercicio:

Datos: -
T cabeza = T2 = 120 ºF
-
Ps = 500 PSI -
σo = 30 dina/cm
-
P fondo = 1000 PSI -
σw = 70 dina/cm
Tubería: -
Qg = 0.5 MPCD
-
ID = 1.995” -
Qo = 400 BPD
-
Ф = 30º -
Qw = 600 BPD
Fluidos: -
Bw = 1.01 Bbl/Bf
-
SGg = 0.65 -
ε/d = 0.0012
-
SGw = 1.07 -
e = 0.0001995
-
ºAPI = 22 -
Atub = 0.0217 pies2
-
SGo = 0.92 -
RGP = 1250 Pcs/Bf
-
T fondo = T1 = 154 ºF

Qw
WOR   1.5
Qo

Qg
GLR   500 Pcs
Qo  Qw Bf

Método de Beggs and Brill

9
PRODUCCION I I PET - 208

- Masa asociada con un barril de tanque:

 1   WOR 
m  Yo  350     Yw  350     0.0764  GLR  Yg
 1  WOR   1  WOR 

 1   1.5 
m  0.92  350     1.07  350     0.0764  500  0.65  378.3 Lb Bfliquido
 1  1.5   1  1.5 

_ _ _
P , T , Z , para un ΔZ = 2000

_
P1  P2
P  14.7
2
_
500  1000
P  14.7  764.7 PSI
2

_
T1  T2
T
2

_
154  120
T  137 º F
2

_
ΔZ
Z  Z
2

_
2000
Z  0  1000
2

- Calculo de Rs , Bo y µo:

-
Rs = 92 Pc/Bf
-
Bo = 1.065 Bbl/Bf
-
µod = 6.53 cp
-
µob = 4.8 cp

- Calculo del factor Z , Bg , ρg , µg:

-
Z = 0.91
-
Bg = 0.0195 Pc/Pcs
-
ρg = 0.0395 gr/cc

Método de Beggs and Brill

10
PRODUCCION I I PET - 208

-
µg = 0.0131 cp

- Calculo del factor ρw , µw:

350  SGw
ρw 
5.615  Bw

350  1.07
ρw   66.7 Lb
5.615  1 Pc

μw  e 1.0031.47910 
2
T 1.98210 5 T 2

μw  e 1.0031.47910 137 1.98210 137   0.63 cp

2 5 2

- Calculo del Q gas:

Qo  _ 
Qg  3.27  10 7  Zg  __
  R  Rs    T  460 
P  

400
Qg  3.27  10 7  0.91   1250  92   137  460   0.1076 Pc
764.7 seg

- Calculo de QL:

QL  6.49  10 5   Qo  Bo  Qw  Bw 

QL  6.49  10 5   400  1.063  600  1  0.06653 Pc

seg

- Calculo de VSL, VSG , VM:

QL
VSL 
At

0.06653
Vsl   3.066 pies
0.0217 seg

Qg
VSG 
At

Método de Beggs and Brill

11
PRODUCCION I I PET - 208

0.1076
VSG   4.96 pies
0.0217 seg

VM  3.066  4.96  8.026 pies

seg

- Calculo de Fo y Fw:

1
Fo 
1  WOR

1
Fo   0.4
1  1.5

WOR
Fw 
1  WOR

1.5
Fw   0 .6
1  1.5

- Calculo de la ρL:

ρL  ρo  Fo  ρw  Fw

350  SGg  0.0764  Rs  SGg

ρo 
5.615  Bo

350  0.92  0.0764  92  0.65

ρo   54.71 Lb
5.615  1.063 Pc

ρL  54.71  0.4  66.7  0.6  61.9 Lb

Pc

- Calculo de las proporciones de masa de flujo:

GL  ρL  VSL

GL  61.9  3.066  189.8 Lb

Pc 2  seg
Gg  ρg  VSG

Método de Beggs and Brill

12
PRODUCCION I I PET - 208

Gg  2.473  4.96  12.3 Lb

Pc 2  seg

Gh  189.8  12.3  202.1 Lb

Pc 2  seg

Utilizando el método:

- Cálculo del No Slip Hold Up:

QL
λ
QL  Qg

0.06653
λ  0.382
0.06653  0.1076

- Calculo del numero de Fraude, µL , µmezcla y σo :

2
VM
NFR 
gd

8.026 2  12
NFR   12.03
32.2  1.995 
μL  μo  Fo  μw  Fw

μL  17  0.4  0.63  0.6  7.18 cp

μmez  6.27  10 4   μL  λ  μg  1  λ    0.0017 cp

σL  σo  Fo  σw  Fw

σL  30  0.4  70  0.6  54 dina

cm

- Cálculo del Nre, No Slip y Nº vel.L:

Gm  d
Nre 
μm

Método de Beggs and Brill

13
PRODUCCION I I PET - 208

202.1  1.995
Nre   19764.2
12  0.0017

0.25
 ρL 
NLv  1.938  VSL   
 σL 

0.25
 61.90 
NLv  1.938  3.066     6.15
 54 

- Calculo de L1, L2 , L3 y L4:

L1  316  λ 0.302  236.2

L2  0.0009252  λ -2.4684  0.010

L3  0.10  λ -1.4516  0.404

L4  0.5  λ -6.738  327.4

- Caudal del modelo de flujo:

Los límites del flujo horizontal son:

Segregado
Limites λ L < 0.01 y N FR < L1 O λ L > = 0.01 y N FR < L2

Trasmisión
Limites λ L > = 0.01 y L2 < N FR + < L3

Intermitente
Limites 0.01 < = L < 0.4 y L3 < N FR < = L1 o λ L > = 0.4 y
L3 < N FR <= L4

Distribuido
Limites λ L < 0.4 y N FR >= L1 o λ L > = 0.4 N FR y > L4

Si λ = 0.382
NFR = 12.03
Como 0.01 < = L < 0.4 y L3 < N FR < = L1 o λ L > = 0.4 y L3 < N FR
<= L4
El modelo de flujo es intermitente.

- Calculo del estancamiento horizontal HL (Horiz.) = HL(0)

Método de Beggs and Brill

14
PRODUCCION I I PET - 208

a  b
HL 
NFR C

Donde a, b, c son determinados del cuadro al modelo de flujo (cuadro

antes descrito)

Si el modelo es intermitente a = 0.845 b = 0.5351 c = 0.0173

0.845  0.382 0.5351

H L(O)   0.483
12.03 0.0173

- Calculo del Factor de corrección del líquido estancado:


C  1     Ln d  e  N LV  NFR g
F

Donde d , e . F y g son obtenidos del cuadro de acuerdo al modo de
flujo.

Como el flujo es intermitente d = 2.96 e = 0.305 F = -0.4473 g =

0.0978

 
C  1  0.382   Ln 2.96  0.382 0.305  6.15 0.4473  12.03 0.0978  0.137

- Calculo del factor de corrección de inclinación del líquido estancado:

ψ1 C Sen 18  Φ   0333  Sen 3 18 Φ 

Si el pozo es vertical ψ = 1+0.3*C

En el caso del ejercicio Ф = 30º

ψ1 0.137 Sen 18  30   0333  Sen 3 18 30   1.09

- Calculo del estancamiento de líquido y la densidad de las dos fases:

HL(30)  HL(0)  

HL(30)  0.483  1.09  0.526

Método de Beggs and Brill

15
PRODUCCION I I PET - 208

tp  L  HL(30)  g  1  HL(30) 

tp  61.9  0.526  2.473  1  0.526   33.73 Lb

Pc

- Calculo de la razón del factor de fricción:

  
Y  2 
 HL(30) 

 0.382 
Y   1.38
 0.526 
2

Ftp
si  eS
Fns

Ln Y 
S
 
 0.0523  3.182  Ln(Y )  0.8725  Ln Y   0.01853  Ln Y 
2

4

Ln(Y) = Ln (1.38) =0.32

0.32
S
 
 0.0523  3.182  0.32  0.8725  0.32 2  0.01853  0.32 4    0.365
Ftp
 e S  1.441
Fns

- Calculo del factor de fricción No Slip:

 0.5 
Fns  0.0056   0.32 
 Nrens 

 0.5 
Fns  0.0056   0.32   0.0267
 19764.5 

- Calculo del factor de fricción de las 2 fases:

Fns  FTp
Ftp 
Fns

Ftp  0.0267  1.441  0.03845

Método de Beggs and Brill

16
PRODUCCION I I PET - 208

- Calculo del Gradiente de Presión:

0.03845  202.1  8.026

33.73  Sen30 
1.995
2  32.2 
P 12
  0.1576 PSI
Z 33.73  8.026  4.96 pie
1
32.2  764.7 144 g Ftp  Gm
 tp  Sen 
P gc 2  gc  d

Z tp  Vm  Vsg
1 _
gc  P 144

P
 0.1576 PSI
Z pie

Método de Beggs and Brill

17