Deduccion de Ecuaciones de Flujo de Gas

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 11

UNIVERSIDAD MAYOR DE SAN ANDRES

FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque

DEDUCCION DE LA ECUACION DE FLUJO DE GAS

ECUACION GENERAL DE FLUJO

La ecuación general de flujo, también llamada la ecuación de flujo Fundamental, para flujo
isotérmico en estado estacionario en un gasoducto es la ecuación básica para relacionar la
caída de presión con el caudal. La forma más común de esta ecuación en los EE.UU es el
sistema consuetudinario (USCS) de las unidades, se da en términos del diámetro de la
tubería, propiedades del gas, presiones, temperaturas, y la tasa de flujo.

𝑇𝑏 𝑃12 − 𝑃22 0.5 2.5


𝑄 = 77.54 ( ) ( ) 𝐷 (𝑈𝑆𝐶𝑆 𝑢𝑛𝑖𝑡𝑠)
𝑃𝑏 𝐺𝑇𝑓 𝐿𝑍𝑓

Donde:

Q = caudal de gas, medido en condiciones normales, ft3 / día (SCFD)

Factor f = fricción, sin dimensiones

Pb = presión base, psia

Tb = temperatura base, ºR (460 + ºF)

P = presión aguas arriba, psia

P2 = presión aguas abajo, psia

G = gravedad de gas (aire = 1,00)

T1 = temperatura promedio que fluye el gas, ºR (460 + ºF)

L = longitud del segmento de tubería, in

Z = factor de compresibilidad del gas a la temperatura que fluye, sin dimensiones

D = diámetro de la tubería en el interior, in.


UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque

Debe de tenerse en cuenta que para el segmento de tubería de la sección 1 a la sección 2,


la temperatura del gas T1 se supone que es constante (flujo isotérmico).

En unidades del SI, la ecuación de flujo general se indica como sigue:


0.5
−3
𝑇𝑏 𝑃12 − 𝑃22
𝑄 = 1.1494 × 10 ( )( ) ∗ 𝐷 2.5 (𝑆𝐼 𝑢𝑛𝑖𝑡𝑠)
𝑃𝑏 𝐺𝑇𝑓 𝐿𝑍𝑓

Donde:

Q = caudal de gas, medido en condiciones normales, m3 / día

Factor f = fricción, adimensional

Pb, = presión de base, kPa

Tb = temperatura base, K (273 + ºC)

P1= presión aguas arriba, en kPa

P2 = presión aguas abajo, en kPa

G = gravedad de gas (aire = 1,00)

Tf = temperatura media del gas que fluye, K (273 + ºC)

L = longitud del segmento de tubería, km

Z = factor de compresibilidad del gas a la temperatura que fluye, sin dimensiones

D = diámetro de la tubería, mm en el interior

Debido a la naturaleza de la anterior ecuación, las presiones pueden estar también en MPa
o Bar, siempre y cuando la misma unidad coherente se utilice.

La ecuación se refiere a la capacidad (velocidad de flujo o caudal) de un segmento de tubo


de longitud L, en base a una presión aguas arriba de P1 y una presión aguas abajo de P2
UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque
como se muestra en la Figura. Se asume que no hay diferencia de elevación entre los
puntos aguas arriba y aguas abajo; Por lo tanto, el segmento de tubería es horizontal.

DEDUCIÓN DE LA ECUACIÓN GENERAL DE FLUJO DE GAS:


Por Balance de energia (Del anterior grafico) :
E1 = E2
H1 + ECIN1 + EPOT 2 = H2 + ECIN2 + EPOT 2
ΔH + ΔECIN + ΔEPOT = 0...(1)
Por termodinamica :
ΔH = ΔU + ΔPV = ΔU + PΔV + VΔP....(2)
Primera ley de la termodinamica : ΔU = Q - W...(3)
Reemplazando (3),(2) en (1) :
Q - W + PΔV + VΔP + ΔECIN + ΔEPOT = 0
Donde : W = PΔV - WF ......(4)

Q - PΔV - WF + PΔV + VΔP + ΔECIN + ΔEPOT = 0....(5)


Donde el trabajo por factor de friccion es :
v2 * ΔL
WF = 4 * f *
2g * D
 ΔECIN = v * Δv
 ΔEPOT = gΔz
Reemplazando en (5)
v2 * ΔL
Q +4 * f * + VΔP + v * Δv + gΔz = 0...(6)
2g * D
Q = 0  Debido a que no hay transferencia ni perdida de calor
gΔz = 0  Por estar al mismo nivel

Expresando en su forma diferencial :


v2 * L
4* f * + VP + v * v = 0
2g * D
q γ m m
consideremso : v = * = * V  v = * V..en(6)
A γ A g*A
2
m 
 * V  * L
A
4* f *   m m 1
+ VP + * V * * V = 0 / / * 2
2g * D A g*A V

 m
2
* L 1 m
2
1 R *T*z
4* f * + P + 2 * * V = 0  Donde : V =
2g * A2 * D V A *g V Mg * P
Integrando :
UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque

 m
2
L * L P2
Mg * P
m 2
1
2 V


0
4* f * 2
+ 
2g * A * D P1 R * T * z
P + 2
* 
A * g V1 V
* V = 0

 m Mg *  P22 - P12 
2
2
*L m V 
4* f * 2
+ + 2 * ln  2  = 0
2g * A * D R *T *z *2 A *g  V1 
Ordenando :

 m  V2  Mg *  P1 - P2 
2
2 2 2
*L
m P1 * V1 P2 * V2 V P
4* f * + * ln   =  Donde : = T1 = T2 2 = 1
2g * A2 * D A2 * g  V1  R *T *z *2 T1 T2 V1 P2
Reemplazando y despejando flujo masico :
Mg *  P12 - P22 

 m
2
= R *T *z *2
L 1 P 
4* f * 2
+ 2 * ln  1 
2g * A * D A * g  P2 
Donde flujo masico :
m P1 * Mg q * P1 * Mg
A*v=q= ρ= m=
ρ R *T R *T
Reemplazando :
Mg *  P12 - P22 
q * P1 * Mg R *T *z *2 P 
=  Donde : ln  1  = 0 (considerando despreciable)
R *T L 1 P   P2 
4* f * 2
+ 2 * ln  1 
2g * A * D A * g  P2 
Mg *  P12 - P22  * g * A2 * D R *T π
q= *  Mg = SGg * MAIR ; A = * D2
4* f *R *T*z*L P1 * Mg 4
Reemplazando :

q=
P 2
1 - P22  * g *  π / 4 * D5
2

*
R *T
4 * f * R * T * z * L * SGg * MAIR P1
Finalmente factorizando las constantes :

q=K*
T P
2
1 - P22  * D5
l.q.q.d .
P1 f * T * z * L * SGg
DONDE :
ΔECIN = Variacion de energia cinetica.
ΔEPOT =Variacion de energia portencial.
ΔH = Variacion de entalpia.
UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque

ΔU = Variacion de energia interna.


Q = Calor
W = Trabajo
WF = Trabajo debido al factor de fricción.
f = Factor de fricción.
Δv = Varicaion de velocidad.
v = Velocidad.
g = Constante gravitacional.
Δz = Variacion de elevación.
Z = Elevación.
ΔL = Variación de longitud.
L = Longitud.
D = DIametro.
q = Caudal.
A = Area del ducto.
m= Flujo masico.
γ = Peso especifico.
ρ = Densidad.
V = Volumen.
ΔV = Variación de volumen.
R = Constante universal de los gases.
T = Temperaturta.
z = Factor de compresibilidad.
Mg = Peso molecular del gas.
P = Presión.
MAIR = Peso molecular.
K = Constante.
UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque

Al examinar la ecuación general de flujo, vemos que para un segmento de tubo de longitud
L y diámetro D, Q el caudal de gas (en condiciones normales) depende de varios factores.
Q depende de las propiedades del gas representados por la gravedad G y el factor de
compresibilidad Z. Si la gravedad de gas se incrementa (gas más pesado), el caudal se
reducirá. Del mismo modo, como los Z factor de compresibilidad aumenta, la velocidad de
flujo disminuirá. Además, como el gas que fluye aumenta la temperatura T1, el rendimiento
disminuirá. Por lo tanto, cuanto más caliente esté el gas, menor será la velocidad de flujo.
Cuanto mayor sea el diámetro, mayor es la velocidad de flujo. El término i2 P - P i2
representa la fuerza motriz que hace que la tasa de flujo desde el extremo corriente arriba
hasta el extremo aguas abajo. A medida que la presión P2 aguas abajo se reduce,
manteniendo la presión constante aguas arriba P1, la velocidad de flujo aumentará. Es
obvio que cuando no hay caudal, P1 es igual a P2. Es debido a la fricción entre las paredes
de tubos de gas y que la caída de presión (P1 - P2) se produce desde el punto 1 aguas
arriba hasta el punto 2 aguas abajo. El factor de fricción f depende de la condición interna
de la tubería, así como el tipo de flujo (laminar o turbulento).

A veces la ecuación de flujo general se representa en términos del factor de transmisión F


en lugar del factor de fricción f Esta forma de la ecuación es como sigue.

0.5
𝑇𝑏 𝑃12 − 𝑃22
𝑄 = 38.77𝐹 ( ) ∗ ( ) ∗ 𝐷 2.5 (𝑈𝑆𝐶𝑆 𝑢𝑛𝑖𝑡𝑠)
𝑃𝑏 𝐺𝑇𝑓 𝐿𝑍𝑓

Donde el factor de transmisión F y factor de fricción f están relacionadas por:

2
𝐹=
√𝑓

Y en unidades del sistema internacional:

0.5
−4
𝑇𝑏 𝑃12 − 𝑃22
𝑄 = 5.474 × 10 𝐹 ( ) ∗ ( ) ∗ 𝐷 2.5 (𝑆𝐼 𝑢𝑛𝑖𝑡𝑠)
𝑃𝑏 𝐺𝑇𝑓 𝐿𝑍
UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque
Vamos a discutir varios aspectos de la ecuación de flujo general antes de pasar a las otras
fórmulas para el cálculo de caída de presión.

EFECTO DE AREAS DE LA TUBERIA

Cuando se incluye la diferencia de elevación entre los extremos de un segmento de tubería,


laEcuación de flujo general se modifica de la siguiente manera:
0.5
𝑇𝑏 𝑃12 − 𝑒 5 𝑃22
𝑄 = 38.77𝐹 ( ) ∗ ( ) ∗ 𝐷 2.5 (𝑈𝑆𝐶𝑆 𝑢𝑛𝑖𝑡𝑠)
𝑃𝑏 𝐺𝑇𝑓 𝐿𝑒 𝑍

En unidades del SI:


0.5
−4
𝑇𝑏 𝑃12 − 𝑒 5 𝑃22
𝑄 = 5.747 × 10 ( )∗( ) ∗ 𝐷 2.5 (𝑆𝐼 𝑢𝑛𝑖𝑡𝑠)
𝑃𝑏 𝐺𝑇𝑓 𝐿𝑒 𝑍

Donde:

𝐿(𝑒 𝑠 − 1)
𝐿𝑒 =
𝑠
La longitud equivalente, Le , y el término es toma en cuenta la diferencia de elevación
entre los extremos aguas arriba y aguas abajo del segmento de tubería. El parámetro s
depende de la gravedad gas, factor de compresibilidad del gas, la temperatura que fluye, y
la diferencia de altitud. Se define como sigue en USCS unidades:

𝐻2 − 𝐻1
𝑠 = 0.0375 ( ) ( 𝑈𝑆𝐶𝑆 𝑢𝑛𝑖𝑡𝑠)
𝑇𝑓

Donde

s = parámetro de ajuste de elevación, adimensional

H1 = elevación de aguas arriba, ft

H2 = elevación aguas abajo, pies

e = base de los logaritmos naturales (e = 2.718 ...)

En las unidades del SI, el parámetro de ajuste de elevación s se define como sigue:
UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque

𝐻2 − 𝐻1
𝑠 = 0.0375 ( ) ( 𝑆𝐼 𝑢𝑛𝑖𝑡𝑠)
𝑇𝑓

Donde

H1 = elevación de aguas arriba, m

H2 = elevación aguas abajo, m

En la ecuación de flujo general, se utiliza el factor de compresibilidad Z. Esto debe ser


calculado a la temperatura y presión de gas media que fluye en el segmento de tubería.

Por lo tanto, es importante para calcular primero la presión media en un segmento de


tubería, considerar un segmento de tubería con presión P1, P2 aguas arriba y la presión
aguas abajo. Una presión media para este segmento debe ser utilizado para calcular el
factor de compresibilidad del gas en el medio de la temperatura Tf gas. Como una primera
aproximación podemos utilizar una media aritmética de (P1 + P2) / 2. Sin embargo, se ha
encontrado que un valor más exacto de la presión de gas promedio en un segmento de
tubería es;

2 𝑃1 𝑃2
𝑃𝑎𝑣𝑔 = (𝑃1 + 𝑃2 − )
3 𝑃1 + 𝑃2

Otra forma de la presión media en un segmento de tubería es:

2 𝑃1 3 −𝑃2 3
𝑃𝑎𝑣𝑔 = ( 3 )
3 𝑃1 −𝑃2 3

Debe tenerse en cuenta que las presiones utilizadas en la ecuación de flujo general están
todos en unidades absolutas. Por lo tanto, las unidades de presión de vía deben ser
convertidos a presión absoluta mediante la adición de la presión base.

VELOCIDAD DE GAS EN UN OLEODUCTO

La velocidad del flujo de gas en una tubería representa la velocidad a la que las moléculas
de gas pasan de un punto a otro. A diferencia de una tubería de líquido, debido a la
UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque
compresión, la velocidad del gas depende de la presión y, por tanto, variará a lo largo de la
tubería incluso si el diámetro de la tubería es constante. La velocidad más alta será en el
de aguas abajo extremo, donde la presión es lo menos. Correspondientemente, la menor
velocidad será a aguas arriba, donde la presión es más alta.

Consideremos un gas de transporte de tubería desde el punto A al punto B, como se


muestra en la Figura. Bajo flujo de estado estacionario, en A, el caudal másico de gas se
designa como M y la voluntad sera la misma que la tasa de flujo de masa en el punto B, si
entre A y B no hay inyección o entrega de gas. La masa es el producto del volumen y la
densidad, podemos escribir la siguiente relación para el punto A:

𝑀 = 𝑄𝑟

La tasa de flujo Q se puede expresar en términos de la velocidad de flujo transversal de la


tubería y área de la sección A de la siguiente manera:

𝑄 =𝑢∗𝐴

Por lo tanto, la combinación de las dos últimas ecuaciones y la aplicación de la conservación


de la masa de los puntos A y B, obtenemos:

𝑀1 = 𝑢1𝐴1𝑟1 = 𝑀2 = 𝑢2𝐴2𝑟2

Puesto que el flujo de gas en una tubería puede dar lugar a variación de la temperatura del
punto A al punto B, la densidad del gas también variará con la temperatura y la presión. Si
la densidad y la velocidad en un momento dado son conocidos, la velocidad
correspondiente en el otro punto se puede calcular utilizando la ecuación

𝑢1 𝜌1 = 𝑢2 𝜌2

Si las condiciones de entrada están representados por el punto A y el caudal Q de volumen


en condiciones estándar de 60 ° F y 14.7 psia son conocidos, podemos calcular la velocidad
en cualquier punto a lo largo de la tubería a la que son la presión y la temperatura del gas
P y T, respectivamente.

La velocidad del gas en la sección 1 está relacionada con la velocidad de flujo Q1 en la


sección 1 y área de la tubería de sección transversal A de la siguiente ecuación:

𝑄1 = 𝑢1𝐴

El flujo de masa M en la sección 1 y 2 es el mismo para el flujo en estado estacionario. Por


Consiguiente:
UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque
𝑀 = 𝑄1 𝜌1 = 𝑄2 𝜌2 = 𝑄𝑏 𝜌𝑏

Donde Qb es el caudal de gas en condiciones estándar y ρb es la densidad del gas


correspondiente.

Por lo tanto, la ecuación se simplifica a:


𝜌𝑏
𝑄1 = 𝑄𝑏 ( )
𝜌1

𝑃1
= 𝑍1 𝑅𝑇1
𝜌1

𝑃1
𝜌1 =
𝑍1 𝑅𝑇1

𝑃𝑏
𝜌𝑏 =
𝑍𝑏 𝑅𝑇𝑏

𝑃𝑏 𝑇1 𝑍1
𝑄1 = 𝑄𝑏 ( ) ( ) ( )
𝑃𝑏 𝑃1 𝑍𝑏

𝑍𝑏 = 100 aproximadamente

𝑃𝑏 𝑇1
𝑄1 = 𝑄𝑏 ( ) ( ) 𝑍1
𝑃𝑏 𝑃1

𝑄𝑏 𝑍1 𝑃𝑏 𝑇1 4 ∗ 144 𝑃𝑏 𝑇1
𝑢1 = ( )( ) = 2
𝑄𝑏 𝑍1 ( ) ( )
𝐴 𝑃𝑏 𝑃1 3.14𝐷 𝑃𝑏 𝑃1

𝑄𝑏 𝑃𝑏 𝑍1 𝑇1
𝑢1 = ( 2 ) ( ) (𝑈𝑆𝐶𝑆 𝑢𝑚𝑖𝑡𝑠)
𝐷 𝑇𝑏 𝑃1

Donde

u, = velocidad del gas aguas arriba, pies / s

Qb = caudal de gas, medido en condiciones normales, ft3 / día (SCFD)

D = diámetro de la tubería en el interior, in.

Pb = presión base, psia

Tb = temperatura base, ºR (460 + ºF) P = presión aguas arriba, psia

T, = temperatura aguas arriba de gas, ºR (460 + ºF)

Z1 = factor de compresibilidad de gas en las condiciones aguas arriba, sin dimensiones

𝑄𝑏 Pb Z2 T2
𝑢2 = 0.002122 ( )( )( ) (USCS units)
𝐷 2 Tb P2
UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque

𝑄𝑏 Pb ZT
𝑢2 = 0.002122 ( )( )( )
𝐷 2 Tb p

𝑄𝑏 Pb ZT
u = 14.7349 ( 2 ) ( ) ( ) (SI units)
𝐷 Tb p

Donde:

u = velocidad del gas, m/s

Qb = caudal de gas, medido en condiciones normales, m3 / día

D = diámetro de la tubería, mm en el interior

Pb = presión de base, kPa

Tb = temperatura base, K (273 + ºC) P = Presión, kPa

T = temperatura media del gas que fluye, K (273 + ºC)

Z = factor de compresibilidad del gas a la temperatura que fluye, adimensional

Puesto que el lado derecho de la ecuación contiene relaciones de presiones, cualquier


unidad coherente se puede utilizar, como kPa, MPa, o bar.

𝑍𝑅𝑇
𝑢𝑚𝑎𝑥 = 100√ (𝑈𝑆𝐶𝑆 𝑢𝑛𝑖𝑡𝑠)
29𝐺𝑃

Donde

Z = factor de compresibilidad del gas, sin dimensión constante

R = gas = 10,73 psia ft3 / lb-Moler

T = temperatura del gas, ºR

G = gravedad de gas (aire = 1,00) P = presión del gas, psia

Por lo general, una velocidad operacional aceptable es 50% de lo anterior

También podría gustarte