Captulo 6

Flujo en estranguladores

CAPTULO 6: FLUJO EN ESTRANGULADORES

6.1 Conceptos fundamentales
6.1.1 Estrangulador
La mayora de los pozos productores o fluyentes cuenta con algn tipo de
restriccin en la superficie, como lo es un estrangulador superficial instalado cerca de la
cabeza del pozo para controlar los gastos de produccin (en las cantidades deseadas) y/o las
presiones corriente abajo.
Un estrangulador es una herramienta cuya funcin es la de restringir el paso de un
fluido bajo presin, con el objeto de controlar el gasto del pozo en las cantidades deseadas
y con esto:
Ejercer la contrapresin suficiente para evitar la entrada de arena en el pozo
o la conificacin de agua.
Conservar la energa del yacimiento, asegurando una declinacin ms lenta
de su presin.
Proteger el equipo superficial y subsuperficial.
Escoger un estrangulador para un pozo fluyente es una prctica comn en la
industria petrolera de tal manera que pequeas variaciones en la presin corriente abajo
(esto es, en la presin de la lnea de flujo, provocado, por ejemplo, por el uso de un
separador de descarga) no afecten la presin en la cabeza de la TP y por ende el desempeo
del pozo.
El propsito de un estrangulador es proveer un control preciso del gasto de aceite y
gas en la cabeza del pozo. Un estrangulador es una restriccin en una lnea de flujo que
provoca una cada de presin o reduce el gasto a travs de un orificio. Los estranguladores
son capaces de provocar cadas de presin grandes. Por ejemplo, el gas puede entrar en un
estrangulador a 5,000 lb/pg2 abs y salir a 2,000 lb/pg2 abs.
6.1.2 Flujo crtico y subcrtico
El flujo multifsico a travs de restricciones, puede ser evaluado bajo condiciones
de flujo crtico y subcrtico. Como una prctica de campo estndar o comn, el
comportamiento de flujo en la cabeza del pozo es evaluado bajo condiciones de flujo crtico
mientras que el comportamiento de flujo a travs de estranguladores y vlvulas de
seguridad es registrado en flujo subcrtico.
El flujo crtico o snico de un fluido es definido como el flujo del fluido a la
velocidad equivalente a la velocidad (sin friccin) de propagacin de una onda de presin
(sonido) en el medio (fluido).
El flujo crtico o snico ocurre cuando la velocidad relativa de un fluido en
cualquier punto (vf, pie/seg) es equivalente a la velocidad de onda elstica (vp, pie/seg) en
ese punto, o sea:

225

Captulo 6

Flujo en estranguladores

vf
= M = 1 10 5
vp

(6.1)

Donde M = Nmero de Mach

En funcin de este nmero se definen tres diferentes regmenes de flujo:

Para M < 1 el flujo es subsnico.

Para M > 1 el flujo es supersnico.
Para M = 1 el flujo es snico o crtico.

Cuando M = 1, el rea de flujo alcanza su valor mnimo y se dice que se ha logrado

una condicin de garganta a cuyas propiedades se les denomina crticas y a las anteriores
a ella como relaciones crticas.
El flujo crtico para los gases ocurre aproximadamente cuando:

p2
= 0.528
p1

(6.2)

Mach 1 o flujo snico para gases (1100 pie/seg para aire) es diferente que para los
lquidos.
Para flujo multifsico algunos investigadores han observado que la velocidad snica
para una mezcla de gas-lquido es menor que para una sola fase. Olson (Essentials of
Engineering Fluid Mechanics, Intl. Textbook, 1961) estableci que para mezclas gaslquido la velocidad acstica es menor que para cualquiera de las fases por s solas.
Durante la presencia de flujo crtico o supercrtico, el gasto a travs del
estrangulador alcanza un valor mximo con respecto a las condiciones prevalecientes
corriente arriba. La velocidad del fluido bifsico fluyendo a travs de la restriccin alcanza
la presin snica o presin de velocidad de propagacin de la onda para los fluidos en dos
fases. Esto implica que el flujo es estrangulado porque los disturbios corriente abajo no
pueden propagarse corriente arriba. Por lo tanto disminuciones no tan significantes en la
presin corriente abajo no hacen incrementar el gasto, es decir, el flujo crtico o snico es el
flujo en el cual perturbaciones de presin y temperatura corriente abajo no son transmitidas
corriente arriba tal que puedan afectar el gasto, no as en flujo subcrtico. Si la presin
corriente abajo es gradualmente incrementada, esta presin no podra cambiar el gasto o la
presin corriente arriba, esto hasta alcanzar el lmite de flujo crtico-subcrtico. A partir de
este momento si la presin corriente abajo se incrementa ligeramente cerca de las
condiciones lmite, el gasto y la presin corriente arriba sern afectados, entonces la
velocidad de los fluidos pasando a travs del estrangulador caen debajo de la velocidad
snica. Aqu, el gasto depende de la diferencial de presin, o bien, los cambios en la
presin corriente abajo afectan la presin corriente arriba, este comportamiento es
caracterizado como flujo subcrtico.
Los estranguladores superficiales se disean de manera que el flujo sea crtico,
mientras que en las vlvulas de seguridad subsuperficiales el flujo es subcrtico. El flujo a
travs de restricciones en el pozo tambin es subcrtico.

226

Captulo 6

Flujo en estranguladores

Figura 6.1. Gasto de masa vs relacin de presiones

6.1.3 Coeficiente de descarga
El coeficiente de descarga (CD) es la constante requerida para corregir el gasto
terico al gasto real que se presenta en los estranguladores.
Distintos autores han presentado grficas, tablas y ecuaciones de CD, con las cuales
el gasto de gas lquido o multifsico calculado puede ser corregido para obtener el gasto
real para ambos tipos de flujo: crtico y subcrtico. Esto muestra que CD para flujo crtico
no es el mismo que para flujo subcrtico.
El coeficiente de descarga est definido por:
CD =

real:

gasto medido
gasto calculado

(6.3)

El coeficiente de descarga CD se multiplica por el gasto terico para obtener el gasto

q TPreal = q TPterico C D

(6.4)

Prcticamente, el valor de CD depende del tamao del estrangulador, diseo, tipo de

flujo y rgimen de flujo. En la derivacin de las ecuaciones analticas de flujo a travs de
restricciones, es de uso comn utilizar coeficientes de descarga (CD) como un factor de
modificacin final en la ecuacin de gasto. Normalmente se espera que el empleo de ste
absorba errores debido a suposiciones realizadas mientras se est desarrollando un modelo.
As, los valores de CD dependen de las suposiciones realizadas durante el desarrollo del
modelo. Un modelo perfecto tendr CD = 1.0, esto es, el uso de CD ser redundante si
227

Captulo 6

Flujo en estranguladores

todos los procesos de flujo son tomados en cuenta de manera precisa. Las leyes de la
termodinmica implican que pueden resultar valores de CD menores que la unidad.
Frecuentemente, CD tambin es correlacionado con factores como el nmero de Reynolds,
diferencial de presin, factor de expansin del gas, etc. Obviamente, si un modelo toma en
cuenta variables como stas, tal dependencia no debe de existir.
En el lenguaje comn se ha adoptado decir flujo crtico a lo que estrictamente es
flujo supercrtico, cabe aclarar que el flujo crtico es la condicin de flujo que se alcanza
cuando la velocidad del fluido es igual a la velocidad de propagacin de una onda en la
mezcla del fluido en cuestin. En la prctica difcilmente se puede mantener estable el valor
de la velocidad y se llega a sobrepasar por lo que en realidad es un flujo supercrtico el cual
no tiene mayor diferencia en trminos significativos. En lo posterior se mencionar flujo
crtico para hacer referencia al flujo supercrtico, dejando en claro que el flujo crtico es el
lmite entre flujo supercrtico y flujo subcrtico.

6.2 Modelos para gas

El objetivo de cualquier operacin de produccin de gas es mover el gas de algn
punto en el yacimiento hacia la lnea de venta. Para poder lograr esto, el gas debe de pasar
por muchas reas de cadas de presin, o si es que se utiliza un compresor, la presin se
gana o se pierde.
Aunque todos estos componentes del sistema integral de produccin se pueden
analizar de manera independiente, para poder determinar el desempeo de un pozo, stos
deben de manejarse de manera combinada en un sistema total o en anlisis nodal. sto se
logra ms fcilmente dividiendo el sistema total en dos subsistemas distintos y
determinando los efectos de los cambios realizados en uno o ambos subsistemas en el
desempeo o comportamiento del pozo.
Existen muchas localizaciones en el sistema de produccin de gas donde el gas debe
de pasar a travs de restricciones relativamente pequeas. Algunos ejemplos de estas
restricciones son vlvulas de seguridad subsuperficiales y estranguladores superficiales. El
flujo puede ser crtico o subcrtico.
Una ecuacin general para el flujo a travs de restricciones se puede obtener
combinando la ecuacin de Bernoulli con una ecuacin de estado y asumiendo que no hay
prdidas irreversibles o por friccin. Un coeficiente de descarga emprico se incluye para
tomarse en cuenta por las simplificaciones utilizadas al derivar la ecuacin. La siguiente
ecuacin puede ser utilizada para flujo crtico (snico) y subcrtico (subsnico).

q c.s.

2
k +1

C n p1 d 2 k p 2 k p 2 k
=

g T1 Z1 k 1 p1 p1

(6.5)

C s C d Tcs
p cs

(6.6)

Donde:

Cn =

228

Captulo 6

Flujo en estranguladores

qc.s Gasto de gas.

Cn Coeficiente basado en el sistema de unidades.
d Dimetro interno de agujero abierto para el flujo de gas.
g Densidad especfica del gas (aire = 1.0), adimensional.
k Relacin de calores especficos (Cp/Cv), adimensional.
p1 Presin corriente arriba, unidades absolutas.
p2 Presin corriente abajo, unidades absolutas.
T1 Temperatura corriente arriba, unidades absolutas.
Z1 Factor de compresibilidad a p1 y T1, adimensional.
Cs Coeficiente basado en el sistema de unidades.
Cd Coeficiente de descarga (emprico), adimensional.
Tc.s. Temperatura a condiciones estndar, unidades absolutas.
pc.s. Presin a condiciones estndar, unidades absolutas.
Rpc Relacin de la presin crtica, adimensional.

Las tablas siguientes muestran valores para las constantes en la ecuacin para varios
sistemas de unidades.

Tabla 6.1. Coeficientes y Unidades para la Ecuacin 6.5.

Sistema
Sistema
Sistema
Smbolo
Mtrico
Ingls
Mtrico
Internacional
2
Mpies@cs/d
m /d
m2/d
qsc
pg
mm
mm
d
2
2
lb/pg abs
kg/cm
kPa
pabs
R
K
K
Tabs
27.611
1.6259
1.6259
Cs
Tabla 6.2. Coeficiente para la Ecuacin 6.5
Sistema de
Unidades
Ingls
Mtrico
Mtrico
Internacional

Cd

p @c.s.

T @c.s.

Cn

0.865
0.865
0.865
0.865
0.865
0.865

14.696 lb/pg2
14.696 lb/pg2
1.0332 kg/cm2
1.0332 kg/cm2
101.325 kPa
101.325 kPa

491.68 R
519.68 R
273.16 K
288.72 K
273.16 K
288.72 K

799.06
344.57
371.83
393.01
3.7915
4.0075

Los valores de k se pueden obtener de:

229

Captulo 6

Flujo en estranguladores

k=

Cp
Cv

M Cp
M C p 1.987

(6.7)

Donde:
M = Peso molecular, lbm/mol.
Cp = Calor especfico, BTU/lbm-R
El valor de k tambin se puede obtener de la figura 6.1:

Figura 6.2. Relacin de calores especficos en funcin de la temperatura y la

densidad relativa
La relacin de presin en donde el flujo se vuelve crtico depende del valor de k
para el gas y est dada por:
k

R pc

2 k 1
=

k +1

(6.8)

En el clculo de valores para Cn dados en la tabla 6.2, se utiliz un coeficiente de

descarga de 0.865. El coeficiente de descarga depende actualmente del Nmero de
Reynolds, de la relacin del dimetro de la tubera al dimetro de la restriccin y de la
geometra de la restriccin.

6.3 Modelos multifsicos

Para predecir el comportamiento de flujo para flujo multifsico a travs de
estranguladores se requiere que primero se haga la prediccin de la frontera entre el flujo
crtico y subcrtico, lo cual es ms difcil que para el flujo monofsico de gas. La eleccin
adecuada depende de que se haga el clculo de relacin de presin crtica, por debajo de la
cual, el gasto msico total es constante, o si se estima la velocidad snica de una mezcla
multifsica.
Como el gasto en flujo crtico no depende de la presin corriente abajo del
estrangulador, adems de presentar un comportamiento lineal con respecto a la presin

230

Captulo 6

Flujo en estranguladores

antes del estrangulador, as como el dimetro, esto hace que las soluciones para flujo crtico
sean ms sencillas hasta cierto punto y por lo tanto ms abundantes en comparacin a
aquellas para flujo subcrtico.
Las soluciones para flujo subcrtico requieren para su aplicacin gran cantidad de
informacin referente a los fluidos en trnsito y por otra parte involucran procesos
iterativos, que traen consigo problemas de convergencia.

6.3.1 Modelos para flujo crtico

6.3.1.1 Correlaciones de Gilbert, Ros y Achong

A partir de datos de produccin, Gilbert desarroll una expresin aplicable al flujo

simultneo gas-lquido a travs de estranguladores. En su trabajo describe en forma
detallada el papel del estrangulador en un pozo y analiza cul es el efecto sobre la
produccin de cambios bruscos en el dimetro del orificio.
Tomando como base la relacin entre las presiones antes y despus de un orificio
para flujo snico de una fase, Gilbert recomend para tener flujo snico, una relacin de
0.588 o menor, entre la presin promedio en el sistema de recoleccin (despus del
estrangulador) y la presin en la boca del pozo (antes del estrangulador).
Utilizando datos adicionales Baxendell actualiz la ecuacin de Gilbert,
modificando los coeficientes.
Ros orient su trabajo al flujo de mezclas con alta relacin gas-aceite, en las que el
gas fue la fase continua. En su desarrollo lleg a una expresin similar a Gilbert; pero con
coeficientes diferentes. Aparentemente su expresin la comprob con datos de campo.
Achong tambin revis la ecuacin de Gilbert y estableci una expresin que valid
con ms de 100 pruebas de campo.
La forma general de las ecuaciones desarrolladas por los investigadores citados es:
A qL R B
p1 =
d C

(6.9)

p1 = Presin corriente arriba (lb/pg2)

qL = Produccin de lquido (bl/da)
R = Relacin gas-lquido (pies3/bl)
d = Dimetro del estrangulador (64avos de pg).
A, B, C = Constantes que dependen de la correlacin y se toman los valores de la tabla
6.3.

Analizando la ecuacin 6.9, se ve claramente cmo el gasto es independiente de la

presin corriente abajo del estrangulador.
La ecuacin 6.9 es tambin muy sensible al tamao del estrangulador. Gilbert
mencion que en un error de 1/128 pulgadas en el tamao del estrangulador, puede causar
errores de 5 a 20% en la estimacin de la presin.

231

Captulo 6

Flujo en estranguladores

Tabla 6.3. Coeficientes empricos para correlaciones de flujo crtico en dos fases.
Correlacin
A
B
C
Gilbert
10
0.546
1.89
Ros
17.4
0.5
2.0
Baxendell
9.56
0.546
1.93
Achong
3.82
0.650
1.88

6.3.1.2 Frmula de Ros (adaptacin de Poettman y Beck)

Este modelo fue establecido a partir del trabajo presentado por Ros en 1960, quien
se bas en el anlisis de la ecuacin de balance de energa, desarroll una frmula de
medidor de flujo a partir de un anlisis terico del flujo simultneo gas-lquido a velocidad
crtica (snica) a travs de restricciones (orificios) y una correlacin para el
comportamiento PVT de los fluidos. No se consider produccin de agua. La precisin de
los resultados obtenidos se comprob comparndolos con 108 datos medidos.
Para que exista flujo crtico se supuso que la presin corriente abajo, debe ser al
menos de 0.55 de la presin en la boca del pozo. Bajo estas condiciones el gasto en el
estrangulador es slo funcin de la presin corriente arriba y de la relacin gas-aceite a
condiciones de flujo.
Despus en 1963 Poettman y Beck convirtieron la ecuacin a unidades de campo y
la redujeron a una forma grfica, lo que result en la siguiente expresin:

9273.6 p1
86,400 C D A o

qL =
5.6142 LS + 0.0764 g R VL (1 + 0.5m L )

0.5

0.4513 r + 0.766

r
+
0.5663

(6.10)

Donde:

r=

0.00504 T1 Z1 (R R s )
p1 B o 1
mL =

1
g1

1 + r
L1

VL =

mL

o1

(6.11)

(6.12)

(6.13)

Siendo:
qL = Gasto en barriles @ c.s./da.
CD = Coeficiente de descarga (1.03).
Ao = rea transversal de la garganta en pulgadas cuadradas (la garganta es el rea
transversal mnima de un estrangulador).
232

Captulo 6

Flujo en estranguladores

Figura 6.3. Nomograma de flujo en estranguladores (20o API) de Poettmann y Beck

233

Captulo 6

Flujo en estranguladores

Figura 6.4. Nomograma de flujo en estranguladores (30o API) de Poettmann y Beck

234

Captulo 6

Flujo en estranguladores

Figura 6.5. Nomograma de flujo en estranguladores (40o API) de Poettmann y Beck

235

Captulo 6

Flujo en estranguladores

LS = densidad del lquido en lb/pie3 @ 60F y 14.7 lb/pg2.

g = Densidad especfica del gas referida al aire @ 60F y 14.7 lb/pg2.
R = Relacin gas-aceite en pie3 de gas @cs/barril de aceite @ cs.
pwh = Presin en la cabeza de la TP en lb/pg2.
p1 = Presin corriente arriba, lbf/pie2.
r = Relacin gas libre-aceite a condiciones de flujo, pies3 de gas/pies3 de lquido.
VL = Volumen especfico del lquido (pies3 de lq./lbm de mezcla).
mL = Masa de lquido por unidad de masa de mezcla (adimensional).
T = Temperatura absoluta de la TP, se asume que es 85F (545R).
Z = factor de compresibilidad del gas a la presin de la TP y 85F.
Rs = Relacin de solubilidad a la presin de la TP y 85F.
Bo = Factor de volumen de formacin del crudo a la presin de la TP y 85F.
3
L = densidad del crudo a la presin p y 85F en lbm/pie .
3
g = densidad del gas a la presin p y 85F en lbm/pie .
Utilizando las ecuaciones 6.10 a la 6.13 y correlaciones empricas para determinar
Rs1 y Bo1, Poettman y Beck construyeron grficas para densidades del crudo de 20, 30 y
40API (Figuras.6.3, 6.4 y 6.5), asumiendo una densidad del gas de 0.6, la temperatura de
la TP de 85F y una presin corriente abajo menor a 0.55 veces la presin corriente arriba.
Las tablas son para 20, 30 y 40API respectivamente. Mediante las grficas de las figuras
6.3, 6.4 y 6.5 se pueden resolver los siguientes problemas:
1. Diseo de estranguladores para pozos nuevos.
2. Estimados de las relaciones gas-aceite y gastos de produccin de gas a partir de
pozos existentes, conociendo las presiones en la TP y los gastos de aceite.
3. Prediccin del desempeo de un estrangulador dado conociendo la relacin gasaceite producida.
4. Revisar por obstruccin por parafinas o corte del estrangulador por gas o arena.
Se obtienen buenos resultados a partir de las tablas si es que no se presenta
produccin de agua y si el flujo es bifsico a condiciones de flujo crticas.
6.3.1.3 Correlacin de Omaa R.

En 1968 Omaa desarroll una correlacin (para flujo crtico) entre el gasto, la
presin corriente arriba del estrangulador, la relacin gas-lquido, la densidad de los
lquidos y el tamao del orificio. Dicha correlacin se obtuvo a partir de datos
experimentales de campo tomados en las instalaciones de la Union Oil Company del campo
Tigre Lagon de California en Louisiana para revisar las correlaciones existentes y
desarrollar la suya. Los experimentos de campo se realizaron con agua y gas natural. En
vista de que estos datos estuvieron dentro de rangos muy limitados, esta correlacin no es
aceptada ampliamente debido a:
1. Limitaciones en el tamao del estrangulador (4, 6, 8, 10 12 y 14/64 de pg).
2. Limitaciones en el gasto (0 a 800 bl/da de agua).
236

Captulo 6

Flujo en estranguladores

3. Limitaciones en la presin corriente arriba p1 (de 400 a 1,000 lb/pg2 manomtricas).

4. Uso de agua en vez de aceite o mezcla agua-aceite en los experimentos de campo.
5. Gastos de gas de 0 a 7 MMpies3/da ( g = 0.611)
Sin embargo, para estranguladores con tamao de hasta 14/64 pg, se considera que
esta correlacin es muy precisa.
Las condiciones de flujo crtico se fijaron para una relacin de presiones igual o
menor de 0.546 y una relacin gas-lquido mayor de 1.0.
La ecuacin que se estableci, mediante un anlisis de regresin mltiple, es:
qL =

N qL

1.84 L
L

(6.14)

1.25

Donde:
N q = 0.263 N 3.49 (N p )

3.19

g
L

N=

1
N p = p1
L g
Q=
1+

g
Nd = d L

0.5

(Q )0.657 (N d )1.8

(6.15)
(6.16)

0.5

0.0174 p1

(6.17)

( L L )0.5

1
(R Rs )Bg

(6.18)

5.615 Bo

= 120.872 d L

0.5

(6.19)

La secuencia de clculo para aplicar la correlacin de Omaa puede sintetizarse en

los pasos siguientes:
1. Calcular g, L, y a la presin y temperatura existentes antes del estrangulador.
2. Evaluar N, Np, Q y Nd a las condiciones prevalecientes corriente arriba del
estrangulador.
3. Obtener Nq con la ecuacin (6.15) y qL con la ecuacin (6.14).

237

Captulo 6

Flujo en estranguladores

Antes de usar la ecuacin (6.14) es conveniente comprobar su validez y ajustarla

para las condiciones de flujo observadas en un campo, introduciendo el coeficiente de
descarga.
6.3.1.4 Ecuacin de Ashford

A partir de un balance de energa y considerando que el fluido se expande

politrpicamente al pasar por el estrangulador, Ashford deriv una ecuacin que describe el
flujo multifsico, bajo condiciones snicas, a travs de un orificio.
Para compensar la ecuacin por las suposiciones incluidas en su desarrollo, se
introdujo en ella un coeficiente de descarga. Sin embargo, al evaluarla, comparando sus
resultados con datos medidos en 14 pozos, se encontr que el coeficiente de descarga
resultaba muy cercano a la unidad.
En su derivacin Ashford supuso una relacin de calores especficos k = 1.04 y una
relacin de presiones, para obtener flujo snico en el orificio de 0.544.
La ecuacin propuesta por Ashford en unidades de campo es:

1.53 d2 C D p1 [(T1 + 460 ) Z1 (R Rs ) + 151 p1 ]( o + 0.000217 g Rs + WOR w )

0.5

qo =

(Bo + WOR )0.5 [(T1 + 460) Z1 (R Rs ) + 111 p1 ]( o + 0.000217 g

R + WOR w )
(6.20)

donde:
d = Dimetro del estrangulador (64avos de pg).
P1 = presin corriente arriba, psia
6.3.1.5 Velocidades Snicas de Nguyen

Nguyen estudi la velocidad snica en sistemas de dos fases como funcin del
patrn de flujo.
Para flujo estratificado, una velocidad snica combinada no existe porque cada fase
es continua en la direccin axial. Una velocidad snica efectiva existe en cada fase que est
influenciada por la otra fase. Si se somete el lquido y el gas a un cambio de presin de
manera simultnea, la perturbacin se propaga con velocidades diferentes en ambas fases
en la direccin axial.
La ecuacin 6.21 es la velocidad snica efectiva para la fase gaseosa y muestra que
la velocidad snica efectiva est gobernada primariamente por la velocidad snica del gas
porque el segundo trmino en el denominador es pequeo.
vE*2, g =

1
1 H L g 1

+
v*g2 1 H L L vL*2

La expresin paralela para la velocidad snica efectiva en la fase lquida es:

238

(6.21)

Captulo 6

Flujo en estranguladores

vE*2, L =

1
1 1 H L L 1

+
vL*2 H L g v*g2

(6.22)

g y v*g2 en el segundo trmino del denominador de la ecuacin 6.22 son

relativamente pequeos, dando una mayor influencia en la velocidad snica efectiva en la
fase lquida.
En contraste con el flujo estratificado, una expresin combinada de la velocidad
snica se desarroll para una unidad slug idealizada. La ecuacin 6.23 da el resultado.

v* =

v L* v *g

(6.23)

H L v *g + (1 H L )v L*

Para flujo homogneo, Nguyen combin expresiones para las velocidades snicas
de cada fase fluyendo dentro de una frontera elstica con el concepto de que el frente de la
onda pasa de manera secuencial a travs de zonas de lquido y gas dentro de la mezcla
homognea. La ecuacin 6.24 da la expresin resultante.

v* =

(1 )

1 L
(1 ) g

+ +
+
*
2
*
2
*
2
v
g v g
vg
L v L*2
L

(6.24)

6.3.1.6 Velocidad Snica de Wallis

En la ecuacin 6.25 Wallis present una expresin para calcular la velocidad snica
o de onda de compresibilidad de una mezcla homognea.

g
L
vm* = ( g g + L L )
+
v *2 v *2

L L
g g

0.5

(6.25)

La velocidad snica de la mezcla homognea no cae necesariamente entre las

velocidades snicas de cada fase y en algunas circunstancias puede ser mucho menor que
ambas. Por ejemplo, una mezcla agua/aire a presin atmosfrica va a tener una velocidad
snica de 1,100 pie/seg, una relacin de densidad de 0.0012 y una velocidad snica mnima
de la mezcla de slo 75 pie/seg. Wallis dijo que la velocidad snica de una mezcla
homognea pasa a travs de un mnimo a una fraccin invlida sin resbalamiento de 0.5.
6.3.1.7 Ashford y Pierce

Ashford y Pierce desarrollaron una expresin para el gasto msico total de una
mezcla multifsica. Ellos asumieron flujo isoentrpico a travs de la restriccin, lquido
239

Captulo 6

Flujo en estranguladores

incompresible, el lquido no se flashea" (separacin flash) en la restriccin y una mezcla

homognea. La ecuacin 6.26 asume que la derivada del gasto con respecto a la relacin de
presin es cero en la frontera crtica.
2 R1

k 1 + R1 yc
yc =

R1
e
b
1 yc yc + 1 yc 1

R1

(6.26)

Donde:

b=

k 1
k

(6.27)

e=

k +1
k

(6.28)

yc es definido como la relacin de la presin en la garganta del estrangulador a la

presin corriente arriba, p2/p1.
La ecuacin anterior requiere de un proceso iterativo para determinar los valores de
yc como funcin de la relacin gas/lquido in situ para diferentes valores de k. La relacin
gas/lquido in situ a condiciones corriente arriba, R1, se puede calcular fcilmente como la
relacin de las velocidades superficiales del gas y del lquido que se determinan a
condiciones inmediatamente corriente arriba del estrangulador.
6.3.1.8 Sachdeva

Sachdeva llev a cabo un estudio combinado experimental y terico que result en

las siguientes ecuaciones para determinar yc, el cual es definido como la relacin de la
presin en la garganta del estrangulador a la presin corriente arriba, p2/p1.
k

N k 1
yc =
D

(6.29)

Donde:

N=
Y:

240

(1 xg1 ) g1 (1 yc )
k
+
k 1
x g1 L

(6.30)

Captulo 6

Flujo en estranguladores

D=

k
n n ( 1 x g1 ) g 2 n ( 1 x g1 ) g 2
+ +
+

2 x g1 L
k 1 2
x g1 L

(6.31)

En las ecuaciones 6.30 y 6.31 el parmetro n y la fraccin o calidad de masa de gas

in situ corriente arriba, xg1, se determinan a partir de las ecuaciones 6.32 y 6.33
respectivamente:

n = 1+

x g1 (C pg Cvg )

x g1 Cvg + (1 x g1 ) C L

(6.32)

Y:

xg1 =

wg1
wg1 + wL1

(6.33)

Donde:

wg1 = 0.0764 g q Lsc (R p f o Rs1 )

wL1 = 5.6142 q Lsc ( f o Bo1 o1 + f w Bw1 w1 )

(6.34)
(6.35)

La ecuacin 6.26 es adimensional, as que se puede utilizar cualquier tipo de

unidades consistentes. La determinacin de yc a partir de la ecuacin 6.26 requiere de un
procedimiento iterativo. Se asume primero un valor de yc en la ecuacin 6.29. Esto permite
un clculo de yc. Un mtodo directo de substitucin es adecuado cuando se utiliza el valor
calculado de yc para la siguiente suposicin hasta que los valores calculados y supuestos de
yc que converjan de acuerdo a una tolerancia predeterminada. Se recomienda un valor de
0.5 para la primera suposicin.

6.3.2 Modelos para flujo crtico-subcrtico

6.3.2.1 Fortunati

Fortunati present un mtodo emprico que puede ser utilizado para calcular el flujo
multifsico crtico y subcrtico a travs de estranguladores. El asumi una mezcla
homognea y aclar que esa suposicin era vlida siempre y cuando vm sea mayor que 32.8
pie/seg y que el nmero Froude de la mezcla sea mayor a 600, adems asumi que no haba
resbalamiento entre las fases, aunque reconoci que el resbalamiento existe incluso para
lquidos inmiscibles.
Utilizando datos experimentales, Fortunati desarroll la figura 6.6, la cual puede ser
utilizada para definir la frontera entre el flujo crtico y el subcrtico. Las curvas de la figura
241

Captulo 6

Flujo en estranguladores

6.6 se basaron en una presin de 0.137 MN/m2 (19.8 psia) corriente abajo. Se utiliza la
ecuacin 6.36 para calcular la velocidad Fortunati de la mezcla a partir de la velocidad
actual de la mezcla y de la presin corriente abajo, en la cual sus variables deben estar en
unidades congruentes.

Figura 6.6. Velocidad de mezclas de gas-aceite a travs de estranguladores.

vm 2

= v

p2
p2 F

(6.36)

Donde:

= (1 3 )

0.38

Y:

242

(6.37)

Captulo 6

Flujo en estranguladores

vm2 = velocidad de la mezcla o fase correspondiente a la presin actual p2, m/seg.

v = velocidad de la mezcla o fase correspondiente a la p2 = 0.137 MN/m2, la cual se
lee de la figura 6.6 m/seg.
p2F = 0.137 MN/m2, presin corriente abajo del estrangulador usada para graficar las
curvas experimentales de la figura 6.6.
p2 = presin corriente abajo del estrangulador actual, MN/m2
6.3.2.2 Perkins

El mtodo de anlisis fue probado comparando gastos medidos y calculados de

1432 datos de la literatura comprendiendo flujo crtico y subcrtico de aire-agua, airekeroseno, gas natural, gas natural-aceite, gas natural-agua y nicamente agua.
Perkins no consider sistemas en los cuales se involucra condensacin de la fase
gaseosa, tales como vapor de agua, as como condiciones en las cuales se pueden alcanzar
presiones de garganta igual y menor a la presin de vapor del agua.
Iniciando con la ecuacin general de energa, la relacin entre las variables en
cualquier punto del sistema fluyendo, puede ser determinada considerando las siguientes
suposiciones:
La temperatura vara con la posicin, pero en cualquier punto todas las fases
estn a la misma temperatura.
La velocidad vara con la posicin, pero en cualquier punto todos los
componentes estn movindose con la misma velocidad.
El factor de compresibilidad del gas es constante.
Los lquidos tienen una compresibilidad despreciable comparada a la del gas
Los cambios de elevacin son despreciables
El proceso de flujo es adiabtico y sin friccin.
De acuerdo con la ecuacin general de energa, la velocidad de descarga
isoentrpica de un estrangulador puede ser calculada con la ecuacin 6.38 y el gasto msico
isoentrpico puede ser calculado con la ecuacin 6.39. El trmino yc es definido como la
relacin de la presin en la garganta del estrangulador a la presin corriente arriba, p2/p1. El
valor de yc el cual conduce al flujo crtico, y el gasto mximo posible, es obtenido de
resolver la ecuacin 6.40.

288 g p 1 y nn1 + xo1 + xw1 p (1 y )

c
1
1
c
1
c

o w

v2 =
2

(xg + L1 )
A

1 2

A1 x y n +
g
c
L1

0.5

(6.38)

243

Captulo 6

Flujo en estranguladores

n 1

(
)

1
y
+
1

c
L1
c

288 g c p1

Wi = A2
2

v1

2
1

A2 (x g + L1 )

x g y c n + L1
1

1

A1 x g yc n + L1

(6.39)

n 1

n
1

y
c

+ 2 L1 (1 yc )

2
(1+ n )
2

2 x +

A2
x g1 (1+ n ) A2 x g1 (x g1 + L1 )yc n

g1
L1
n

yc
1
+
1
2

A1 n x y 1n +
A1 x g1 yc n + L1 n
g1 c
L1

(6.40)

A x g1 + L1

n 1 1n
n +
= 1 2
x
y
yc + L1

g1
c
L1
1

A
n

1 x g1 yc n + L1

Donde:
xo1

L1 = n1

o1

xw1

w1

(6.41)

= x g + (x g Cvg + xo Cvo + x w Cvw )

Z R

(6.42)

En la expresin de Perkins, x es la fraccin en peso de una fase dada en la corriente

de flujo y Cv es la capacidad calorfica a volumen constante, (pie-lbf)/(lbm-F). se puede
evaluar presumiblemente ya sea a condiciones corriente arriba o corriente abajo. El
exponente politrpico de expansin para mezclas utilizado en la ecuacin 6.40 se defini
como sigue:

n=

244

x g k Cvg + xo Cvo + xw Cvw

x g Cvg + xo Cvo + x w Cvw

(6.43)

Captulo 6

Flujo en estranguladores

Para el caso de Ashford y Pierce as como de Sachdeva, se requiere de un proceso

iterativo para determinar yc.
La condicin de flujo crtico puede ser determinada por diferenciacin de la
ecuacin general de energa y para distinguir entre gasto crtico y subcrtico, es necesario
para conocer la magnitud de la presin de descarga en la garganta del estrangulador.
Para flujo subcrtico la relacin de Perry puede ser aproximada por:
p3 = p1

( p1 p4 )
d
1 c
d d

1.85

(6.44)

El mtodo para distinguir entre flujo crtico y subcrtico y como manejar flujo
crtico y subcrtico es como sigue:
1. Calule p2 (la presin en la garganta del estrangulador cuando el flujo es
crtico), resolviendo la ecuacin 6.40.
2. Calcule p3 (la presin de descarga justo corriente abajo del estrangulador),
resolviendo la ecuacin 6.44.
3. Si p2 > p3, entonces el flujo es crtico y yc = p2/p1, se deber usar las
ecuaciones 6.38 y 6.39 para calcular la velocidad isoentrpica y gasto
msico.
4. Si p2 = p3, entonces el flujo est en el lmite entre flujo crtico y subcrtico.
(Use yc = p2/p1 = p3/p1 para calcular la velocidad isoentrpica y gasto
msico).
5. Si p2 < p3, entonces el flujo es subcrtico. (Use yc = p3/p1 en las ecuaciones
6.38 y 6.39 para calcular velocidad isoentrpica y gasto msico).
En todos los casos, el gasto real es determinado por multiplicar el gasto isoentrpico
por el coeficiente de descarga CD.
6.3.2.3 Modelo de Ashford y Pierce

Ashford y Pierce establecieron una ecuacin que describe la dinmica de las cadas
de presin y capacidades de flujo en condiciones de flujo multifsico. Este modelo
relaciona el comportamiento del estrangulador en ambos regmenes de flujo, crtico y no
crtico.
La capacidad y cadas de presin que se presentan en la restriccin se han
relacionado con sus dimensiones y las propiedades de los fluidos manejados.
Los datos usados, reflejan el comportamiento de una vlvula de seguridad OTIS tipo
J-22J037. Sin embargo, el modelo puede usarse para estimar las cadas de presin a travs
de cualquier dispositivo que restrinja el flujo.
Para la validacin del modelo, se dise una prueba de campo en un pozo fluyente.
Tanto las cadas de presin como el gasto se midieron directamente y luego se compararon
con datos anlogos obtenidos del modelo. Esta informacin se us para determinar el

245

Captulo 6

Flujo en estranguladores

coeficiente de descarga del orificio, definido por la relacin de gasto medido entre el gasto
calculado.
La ecuacin obtenida por los autores es:
qo = 1.9706 C d2

(6.45)

donde:
1

= (Bo + WOR ) 2

(6.46)

k
0.5
k

T1 Z1 (R Rs ) 1 E + 198.6 (1 E ) o + 0.000217 g Rs + WOR w

T1 Z1
(R Rs )E k o + 0.000217 g Rs + WOR w
198.6 +
p
1

]
(6.47)

k=

E=

Cp

(6.48)

Cv

p1
p2

(6.49)

m = k 1

(6.50)

Para un dimetro dado de estrangulador, su coeficiente de descarga se muestra en la

tabla 6.4 y obtuvo de la siguiente forma:

Tabla 6.4. Coeficientes de descarga para un dimetro dado.

Dimetro del
Gasto
Gasto
Coeficiente de
estrangulador
medido
calculado
descarga
(1/64 pg)
(lb/da @cs)
(lb/da @cs)
16

559

615

0.9089

16

484

402

1.2039

16

334

C=

246

0.9089 + 1.2039
= 1.0564
2

Captulo 6

Flujo en estranguladores

Los resultados obtenidos en las pruebas, para dimetros de estrangulador de 4, 16 y

20/64 de pulgada, se muestran en la tabla 6.5:

Tabla 6.5. Valores de los coeficientes de descarga para

estranguladores de 4, 16 y 20/64 de pg.
Dimetro del
Coeficiente de
estrangulador
descarga
(1/64 pg)
14

1.1511

16

1.0564

20

0.9760

Para dimetros menores de 20.81/64 de pulgada puede aproximarse el coeficiente de

descarga con la siguiente ecuacin, que es el resultado del ajuste de la relacin entre el
dimetro del estrangulador y el coeficiente de descarga:

C = 2.398 0.477 ln (d )

(6.51)

Para valores mayores, el valor de C es constante e igual a 0.95.

6.4 Modelo Mecanstico de la Universidad de Tulsa por Beggs y Brill.

El modelo mecanstico del estrangulador por Beggs y Brill es vlido para flujo
crtico y subcrtico. Es puramente terico y est basado en la combinacin de la ecuacin de
Bernoulli con una ecuacin de continuidad.

p TP = p L L + p g g
qL

p L =
2 g c 144 C L A c

p g =
2 g c 144 Z C g A c

(6.52)

(6.53)
2

d 2 1 p 2 p1

Z = 1.0 0.41 + 0.35

d1 k p1

(6.54)

(6.55)

247

Captulo 6

Flujo en estranguladores

C=

CD
d1
1
d2

(6.56)

La cada de presin total de un sistema de dos fases est dado por:

2

(6.57)

d 4

qm
p L = L 1 + 1

2
d 2 8083 d1 C DL

(6.58)

pTP

C
= p L 1 + g DL

Z DDg

Donde:

Y:

AC = rea de la seccin transversal del estrangulador, pie2

C = coeficiente de flujo
Z = factor de compresibilidad
d1 = dimetro de la tubera corriente arriba, pg
d1 = dimetro del orificio
El subndice TP se refiere a dos fases

6.5 Ejemplos
Ejemplo 1
Dados los siguientes datos mostrados abajo, calcula el tamao de estrangulador en
64avos de pulgada, usando las correlaciones de Gilbert, Ros, Baxandell, Achong, Poettman
y Beck, Omana y Ashford para flujo critico.
qo=
qg=
RGL =
dt =
Ap =
CD =
API =
g =
T=
P1 =
L1 =
Rs1 =

248

1800
1.116
620
2.994
0.0487
1
30
0.7
128
1663
9
335

blp/d
Mmpc/d
Pg
Pie2

F
Psia
Dinas/cm

Captulo 6

Flujo en estranguladores

Bo1 =
Z1 =
g1 =

1.10
0.78
6.856

lbm/pie3

L1 =

52.59

lbm/pie3

Gilbert, Ros, Baxendell, Achong

Utilizando A, B, C de la tabla 6.3 tenemos que:
1

A q 'L R BP C

d =
p
1

Correlacin

d
1

(10 )(1800 )(620 )0.546 1.89

= 22.5934
=

1663

Gilbert

(17.4 )(1800 )(620 )0.5 2

= 21.6552
=

1663

Ros

Baxendell

(9.56 )(1800 )(620 )0.546 1.93

= 20.6916
=

1663

(3.82 )(1800 )(620 )0.65 1.88

= 19.6494
=

1663

Achong

Poettman and Beck

m = 5.614 L + 0.0765 g R p
s.c.

m = 5.614 (0.8793)(62.4) + (0.0765)(0.7 )(620) = 341

B gl = 0.028 (0.78)

v sg1 =

q 'o (R p R s1 ) B g1
A p (86400 )

lb m
bl

128
= 0.0077
1663

1800 (620 335) 0.0077

pie
= 0.9415
0.0487 (86400 )
seg

249

Captulo 6

Flujo en estranguladores

v sL1 =

5.614 q o B o1 5.614 (1800 )(1.10 )

=
= 2.6418
86400
86400
v sgl 0.9415
R=
=
= 0.3564
v sLl 2.6418

g c.s.

g1 =

mL =

B g1
1

1 + R g1
L1

qL m
86400 C D

(1800)(341.23)
(86400)(1)

lb
0.7 (0.765)
= 6.9346 m3
0.0077
pie
1

6.856
1 + 0.3564

52.59

mL

= 0.9556

0.7865
= 0.0182
L1 52.59
1
=
0.4513 R + 0.766
9273.6 p1

v L (1 + 0.5 m L )
R + 0.5663

vL =
Ao =

1
=
(0.4513) 0.3564 + 0.766
(9273.6)(1663)

(0.0182)(1 + (0.5)(0.9556))
0.3564 + 0.5663

= 5.7247 10 4 pie 2

5.7247 10 4 =

d2
4

d = 0.0270 pie = 0.0270 (12 ) = 0.3240 pg

d=

Omaa
R=

250

20.7344
pg
64

N gv
N Lv

v sg1
v sL1

0.9415
= 0.3564
2.6418

Captulo 6

Flujo en estranguladores

QD =

1
1
=
= 0.7373
1 + R 1 + 0.3564

N p = 1.74 10

p1
L1 L1

= 1.74 10 1663
(
)
52.59
9

0.5

0.5

= 1.3301

1.25

N qL

= 1.84 q L L1
L1

= 30089.5820
1.25

52.59
N qL = 1.84 (1800)

L1
= 290.0759 d o
L1

N d = 120.872 d o

L1 6.856
=
= 0.1304
g1 52.59

N =

N qL = (0.263) N

= 30089.5820

3.49

)(

)(N )(Q
3.19

0.657
D

)(N )

)(

N qL = (0.263) 0.1304 3.49 1.33013.19 0.73730.657

30089.5820 = 1.77 10 7 d o

d 1.80
=
O

1.80

)((290.0759 d ) )
1.80

1.80|

30089.5820
= 0.0017
1.773 10 7

d o = 0.0289 pie
d o = 0.0289 (12 ) = 0.3468 pg

do =

22.20
pg
64

251

Captulo 6

Flujo en estranguladores

Ashford
A = 5.04 10 3 T1 Z1 (R p R s1 ) = 658.7905

B = 62.4

145.7
= 62.4 (0.8793) = 54.8683
135.7 + 30
=

1
1
=
= 0.9535
B O1
1 .1
1

dc =

q o + 0.56
P1

A + 0.76 p1
(B + 0.01353 g R p )2 + (B + 0.01353 g R s1 )1

658.7905

1800
+ 0.56
1663

658.7905 + 0.76 (1663)

(54.8683 + 0.01353 (0.7 )(620))2 + (54.8683 + 0.01353 (0.7 )(335))1

= 18.5899

252

18.59
pg
64