Deduccion de Ecuaciones de Flujo de Gas
Deduccion de Ecuaciones de Flujo de Gas
Deduccion de Ecuaciones de Flujo de Gas
FACULTAD DE INGENIERIA
TRANSPORTES DE HIDROCARBUROS II PET - 214 Doc: Ing.Sandalio Choque
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.
Donde:
Donde:
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.
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
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).
0.5
𝑇𝑏 𝑃12 − 𝑃22
𝑄 = 38.77𝐹 ( ) ∗ ( ) ∗ 𝐷 2.5 (𝑈𝑆𝐶𝑆 𝑢𝑛𝑖𝑡𝑠)
𝑃𝑏 𝐺𝑇𝑓 𝐿𝑍𝑓
2
𝐹=
√𝑓
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.
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
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
2 𝑃1 𝑃2
𝑃𝑎𝑣𝑔 = (𝑃1 + 𝑃2 − )
3 𝑃1 + 𝑃2
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.
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.
𝑀 = 𝑄𝑟
𝑄 =𝑢∗𝐴
𝑀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
𝑄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
𝑄𝑏 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:
𝑍𝑅𝑇
𝑢𝑚𝑎𝑥 = 100√ (𝑈𝑆𝐶𝑆 𝑢𝑛𝑖𝑡𝑠)
29𝐺𝑃
Donde