UNIVERSIDAD MAYOR DE SAN ANDRES

FACULTAD DE INGENIERIA
INGENIERIA PETROLERA

TRANSPORTE Y ALMACENAJE DE
HIDROCARBUROS II
INTEGRANTES:
1. CAMPOS ORTIZ YESHENIA VERONICA
2. CHINO ARELLANO RUDY
3. FLORES POMA LILIANA LUPE
4. HUANCA QUISPE NILDA
5. LOPEZ SANCHEZ ALEJANDRA MADELEN
6. MORALES LOPEZ OMAR LUIS
7. POMA PAREDES JESSICA MERY
8. QUENTA MARQUEZ GARY BRANDON
9. TORREZ CALLE JUAN VICTOR

DOCENTE:
ING. SANDALIO CHOQUE POMA
FECHA:

22 de marzo del 2019

1
 DEDUCCION DE LA ECUACION
GENERAL DE FLUJO
Para la deducción de la ecuación general de flujo, consideramos que un flujo de gas
en una tubería circular, como se muestra en la figura.

Donde:
𝑃1 𝑦 𝑃2 = 𝑃𝑟𝑒𝑠𝑖𝑜𝑛𝑒𝑠
𝑑 = 𝐷𝑖𝑎𝑚𝑒𝑡𝑟𝑜
𝐿 = 𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑑𝑒 𝑙𝑎 𝑡𝑢𝑏𝑒𝑟𝑖𝑎
𝑄 = 𝑓𝑙𝑢𝑗𝑜 ( 𝐶𝑎𝑢𝑑𝑎𝑙)

Para obtener la ecuación general de flujo de gas se puede obtener a partir de la

ecuación de Bernoulli que esta se expresa de la siguiente manera:
1
𝑃+ 𝜌 𝑣 2 + 𝜌 𝑔 ℎ = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
2
Remplazando 𝛾 se obtiene:
𝑃 1
+ 𝑣 2 + ℎ = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 … (1)
𝜌 2𝑔
La ecuación (1) se aplica para un flujo ideal, pero en este caso se utilizará la
variación siguiente de la ecuación de Bernoulli para un flujo real:
1 𝑣
𝑑𝑃 + 𝑑𝑣 + 𝑑ℎ − ℎ𝐿 = 0 … (2)
𝜌 2𝑔

2
En este entendido consideramos que a lo largo de la tubería ocurrirá lo siguiente:

1. El Flujo es isotérmico, por lo tanto, la energía interna no varía ∆U = 0.

2. Existe una diferencia de nivel por lo tanto la energía potencial varía con la
altura, debido al peso de la columna del fluido (columna hidrostática).
3. Es despreciable la variación de energía cinética, siempre que la longitud de
la tubería sea suficientemente grande respecto al diámetro (L>>>>>>>D).
4. Solo se considera la variación de la energía de presión, la cual es única que
balancea el trabajo de fricción realizado por el gas.

Tomando en cuenta los puntos descritos anterior mente podemos decir que:

El flujo por ser isotérmico y adiabático, la energía cinética producida por el fluido
parcialmente turbulento es mínimo (3%-8%) razón por la cual no se considera tal
como mencionamos anteriormente

𝑣
( ) 𝑑𝑣 ≅ 0
2𝑔

Por otra parte, en transporte de fluidos en la tubería muestra una sola perdida de
energía que se produce, que es la fricción.

ℎ𝐿 = 𝐹

Por lo que la ecuación (2) se transforma a:

1
𝑑𝑃 + 𝑑ℎ − 𝐹 = 0 … (3)
𝜌
Para describir el factor de fricción, la caída de presión y el factor de fricción de
Faning, consideremos que el factor de fricción (frotamiento entre las partículas del
fluido, así como la pared interna de la tubería) es directamente proporcional a:

∗ 𝐴𝑟𝑒𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒 𝑚𝑜𝑗𝑎𝑑𝑎 = 𝜋 𝐿 𝐷

∗ 𝐷𝑒𝑛𝑠𝑖𝑑𝑎𝑑 𝑑𝑒𝑙 𝑓𝑙𝑢𝑗𝑜 = 𝜌

𝑣2
∗ 𝐶𝑢𝑎𝑑𝑟𝑎𝑑𝑜 𝑑𝑒 𝑙𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 =
2𝑔

Bajo la consideración anterior está dado por lo siguiente:

𝑣2
𝐹 ∝ 𝑑𝐿𝜋𝑑𝜌 ( )
2𝑔

3
Pero la longitud se considera como variación 𝐿 = 𝑑𝐿

𝑣2
𝐹 ∝ 𝑑𝐿 𝜋 𝑑 𝜌 ( ) 𝑑𝐿 … (4)
2𝑔

Donde el trabajo realizado debido a la fricción resulta:

𝑣2
𝑑𝑤𝐹 = 𝑓 𝑑𝐿 𝜋 𝑑 𝜌 ( ) 𝑑𝐿 … (5)
2𝑔

Donde:

𝑓 = Factor de fricción de Fanning

𝑊𝑓 = trabajo realizado debido al factor de fricción

Dividiendo la ecuación 5 entre:

𝜋 2
∗ 𝐴𝑟𝑒𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑒𝑐𝑐𝑖𝑜𝑛 𝑐𝑖𝑟𝑐𝑢𝑙𝑎𝑟 𝑑𝑒 𝑙𝑎 𝑡𝑢𝑏𝑒𝑟𝑖𝑎 = 𝑑
4

∗ 𝐷𝑒𝑛𝑠𝑖𝑑𝑎𝑑 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 = 𝜌

∗ 𝐷𝑖𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎 𝑑𝑒 𝑙𝑜𝑛𝑔𝑖𝑡𝑢𝑑 = 𝑑𝐿

se obtiene la siguiente expresión:

𝑣2
𝑓 𝑑𝐿 𝜋𝑑𝜌 (2𝑔) 𝑑𝐿
𝑑𝑤𝐹 = 𝜋 2 … (6)
4 𝑑 𝜌 𝑑𝐿
Simplificando:
2𝑓𝑣 2 𝑑𝐿
𝑑𝑤𝐹 =
𝑔𝑑
Integrando de 0 a L:
𝐿
2𝑓𝑣 2
𝑑𝑤𝐹 = ∫ 𝑑𝐿
0 𝑔𝑑

2𝑓𝑣 2 𝐿
𝑤𝐹 = … (7)
𝑔𝑑

Donde:
∆𝑃
𝑤𝐹 = … (8)
𝜌

4
Igualando las ecuaciones 7 y 8:

∆𝑃 2𝑓𝑣 2 𝐿
=
𝜌 𝑔𝑑

Despejando la variación de la Presión se tiene:

2𝑓𝑣 2 𝜌𝐿
∆𝑃 = … (9)
𝑔𝑑

Donde:
∗ ∆𝑃 = 𝐶𝑎𝑖𝑑𝑎 𝑑𝑒 𝑝𝑟𝑒𝑠𝑖𝑜𝑛 𝑒𝑛 𝑒𝑙 𝑖𝑛𝑡𝑒𝑟𝑖𝑜𝑟 𝑑𝑒 𝑙𝑎 𝑡𝑢𝑏𝑒𝑟𝑖𝑎 (𝑃𝑠𝑖)
∗ 𝑣 = 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑓𝑙𝑢𝑗𝑜 (𝑓𝑡/𝑠𝑒𝑔)
𝑙𝑏
∗ 𝜌 = 𝐷𝑒𝑛𝑠𝑖𝑑𝑎𝑑 𝑑𝑒𝑙 𝑓𝑙𝑢𝑗𝑜 ( 3 )
𝑓𝑡
∗ 𝐿 = 𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑑𝑒 𝑙𝑎 𝑡𝑢𝑏𝑒𝑟𝑖𝑎 (𝑓𝑡)
∗ 𝑑 = 𝐷𝑖𝑎𝑚𝑒𝑡𝑟𝑜 𝑖𝑛𝑡𝑒𝑟𝑛𝑜 𝑑𝑒 𝑙𝑎 𝑡𝑢𝑏𝑒𝑟𝑖𝑎 (𝑓𝑡)
∗ 𝑓 = 𝐹𝑎𝑐𝑡𝑜𝑟 𝑑𝑒 𝑓𝑟𝑖𝑐𝑐𝑖𝑜𝑛 𝑓𝑎𝑛𝑖𝑛𝑔 (𝐴𝑑𝑖𝑚𝑒𝑛𝑠𝑖𝑜𝑛𝑎𝑙)
∗ 𝑔 = 𝑐𝑜𝑛𝑡𝑎𝑛𝑡𝑒 𝑔𝑟𝑎𝑣𝑖𝑡𝑎𝑐𝑖𝑜𝑛𝑎𝑙

Reemplazando la ecuación 9 en la ecuación 8:

2𝑓𝑣 2 𝐿
𝐹= … (10)
𝑔𝑑
Reemplazando la ecuación 10 en la ecuación 4:

1 2𝑓𝑣 2 𝑑𝐿
𝑑ℎ + ( ) 𝑑𝑃 + = 0 … (11)
𝛾 𝑔𝑑

Multiplicando por 𝛾 2 a la ecuación 11:

1 2𝑓𝑣 2 𝑑𝐿
𝜌2 𝑑ℎ + 𝜌2 ( ) 𝑑𝑃 + 𝜌2 = 0 … (12)
𝛾 𝑔𝑑

Tomando en cuenta:
𝑣
𝑉=
𝑚
𝑣
𝑉=𝑚 𝑣 =1
𝜌
𝑣
1
𝜌𝑏 =
𝑉𝑏
1
𝜌 = … (13)
𝑉

5
Donde: 𝐴 = 𝐴𝑏

𝑄 ∙ 𝜌 = 𝑄𝑏 𝜌𝑏 … (14)

Reemplazando 13 en 14:
1 1
𝑄 ( ) = 𝑄𝑏 ( )
𝑉 𝑉𝑏
𝑉
𝑄 = 𝑄𝑏 ( ) … (15)
𝑉𝑏

Pero, según la ecuación de la continuidad:

𝑄 𝑄
𝑄 =𝑣∙𝐴→𝑣 = → 𝑣 = 𝜋 … (16)
𝐴 2
4𝑑

Reemplazando la ecuación 16 en la ecuación 12:

2 2
1 2
2𝑓 4𝑄 2
𝜌 𝑑ℎ + 𝜌 ( ) 𝑑𝑃 + 𝜌 ( ) ( 2 ) 𝑑𝐿 = 0
𝜌 𝑔𝑑 𝜋𝑑

32𝑓𝑄 2
𝜌2 𝑑ℎ + 𝜌𝑑𝑃 + 𝜌2 𝑑𝐿 = 0
𝜋 2 𝑔𝑑 5
Donde:
𝑉
𝑄 = 𝑄𝑏 ( )
𝑉𝑏
Reemplazando:
2
2
32𝑓 2
𝑉
𝜌 𝑑ℎ + 𝜌𝑑𝑃 + 𝜌 5
(𝑄𝑏 ( )) 𝑑𝐿 = 0. . . (17)
𝜋𝑔𝑑 𝑉𝑏

Sabemos que la ecuación de estado de los gases reales es:

𝑚
𝑃𝑣=𝑍 𝑅𝑇
𝑀
𝑚
𝑃𝑀=𝑍 𝑅𝑇
𝑣
𝑃𝑀
𝜌=
𝑍𝑅𝑇

6
 𝑣 1 1 1
𝑉= =𝑚= = … (18)
𝑚 𝜌 𝑃𝑀
𝑣 𝑍𝑅𝑇

Reemplazando la ecuación 18 en la ecuación 17 obtenemos:

2
2 2
𝑃𝑀 𝑃𝑀 𝑃𝑀 32𝑓 1
( ) 𝑑ℎ + 𝑑𝑃 + ( ) 2 5 (𝑄𝑏 ( )) 𝑑𝐿 = 0
𝑍𝑅𝑇 𝑍𝑅𝑇 𝑍𝑅𝑇 𝜋 𝑔𝑑 𝑃𝑀
𝑉𝑏 ∗ (𝑍𝑅𝑇)

2
𝑃𝑀 2 𝑃𝑀 32𝑓 1
( ) 𝑑ℎ + 𝑑𝑃 + 2 5 (𝑄𝑏 ( )) 𝑑𝐿 = 0
𝑍𝑅𝑇 𝑍𝑅𝑇 𝜋 𝑔𝑑 𝑉𝑏

𝑅
Multiplicamos por ( ):
𝑀

2
𝑃2 𝑀 𝑃 32𝑓𝑅 1
2 2
𝑑ℎ + 𝑑𝑃 + 2 5
(𝑄𝑏 ( )) 𝑑𝐿 = 0
𝑍 𝑅𝑇 𝑍𝑇 𝜋 𝑔𝑀𝑑 𝑉𝑏
Considerando la ecuación 18:

1
𝑉𝑏 =
𝑃𝑀
𝑍𝑅𝑇
Reemplazando y desarrollando:

𝑃𝑚2 𝑀 ℎ2 1 𝑃2
32𝑓𝑀 𝑃𝑏 𝑄𝑏2 𝐿
∫ 𝑑ℎ + ∫ 𝑑𝑃 + 2 5 ∫ 𝑑𝐿 = 0
𝑍𝑚2 𝑅𝑇 2 ℎ1 𝑍𝑚2 𝑇 𝑃1 𝜋 𝑔𝑑 𝑇𝑏 0

Resolviendo la ecuación diferencial:

𝑀 ℎ2
𝑃 𝑃2
32𝑓𝑀 𝑃𝑏 𝑄𝑏2 𝐿
∫ 𝑑ℎ + ∫ 𝑑𝑃 + ∫ 𝑑𝐿 = 0
𝑧𝑅𝑇 2 ℎ1 𝑍𝑚2 𝑇 𝑃1 𝜋 2 𝑔𝑑 5 𝑇𝑏 0

𝑀 ℎ2
𝑃 𝑃2
32𝑓𝑀 𝑃𝑏 𝑄𝑏2 𝐿
∫ 𝑑ℎ + ∫ 𝑑𝑃 + 2 5 ∫ 𝑑𝐿 = 0
𝑧𝑅𝑇 2 ℎ1 𝑍𝑚2 𝑇 𝑃1 𝜋 𝑔𝑑 𝑇𝑏 0

𝑀 𝑃 𝑃22 𝑃12 32𝑓𝑀 𝑃𝑏 𝑄𝑏2

(ℎ − ℎ1 ) + ( − )+ 2 5 (𝐿 − 0) = 0
𝑧𝑅𝑇 2 2 𝑍𝑚2 𝑇 2 2 𝜋 𝑔𝑑 𝑇𝑏

7
 32𝑓𝑀 𝑃𝑏 𝑄𝑏2 𝑀 𝑃𝑚2 𝑃12 𝑃22
𝐿 = (ℎ − ℎ ) + ( − )=0
𝜋 2 𝑔𝑑5 𝑇𝑏 𝑍𝑚𝑅𝑇 2 1 2
𝑍𝑚2 𝑇 2 2

𝑀 𝑃𝑚2 𝑃12 𝑃22 𝜋 2 𝑔𝑑 5 𝑇𝑏

𝑄𝑏2=( (ℎ )
− ℎ2 + ( − )) ∗
𝑍𝑚𝑅𝑇 2 1 𝑍𝑚2 𝑇 2 2 𝐿 ∗ 32𝑓𝑀 𝑃𝑏
Donde el factor de transmisión es:
2
𝐹=
√𝑓
4
𝑓=
𝐹
Reemplazando 𝑓 obtenemos:

𝑀 𝑃𝑚2 𝑃12 𝑃22 𝜋 2 𝑔𝑑 5 𝑇𝑏

𝑄𝑏2 = ( (ℎ − ℎ ) + ( − )) ∗
𝑍𝑚𝑅𝑇 2 1 2
𝑍𝑚2 𝑇 2 2 4
𝐿 ∗ 32 ∗ 𝐹 𝑀 𝑃𝑏
Desarrollando:

𝜋2𝑔 𝑑5 𝑇𝑏 𝑀 𝑃𝑚2 𝑃12 − 𝑃22

𝑄𝑏2 =( ) ( (ℎ − ℎ2 ) + ( ))
32 ∗ 4 𝐿 𝑀 𝑃𝑏 𝑍𝑚𝑅𝑇 2 1 𝑍𝑚2 𝑇 2

𝜋2𝑔 𝑑 5 𝑇𝑏 𝑀 𝑃𝑚2 𝑃12 − 𝑃22

𝑄 = √( ) ( (ℎ )
− ℎ2 + ( ))
32 ∗ 4 𝐿 𝑀 𝑃𝑏 𝑍𝑚𝑅𝑇 2 1 𝑍𝑚2 𝑇 2

Operando:
0.5
2 2 𝐺𝐸 ∗ 𝑃𝑚2 (𝐻2 − 𝐻1 )
𝑇𝑏 𝑃1 − 𝑃2 − 0.0375 ( )
𝑍𝑚 𝑇
𝑄 = 38.744 𝐹[ ] 𝑑 2.5
𝑃𝑏 𝐺𝐸𝐿𝑇𝑍𝑚

Multiplicamos por E:

0.5
𝐺𝐸 ∗ 𝑃𝑚2 (𝐻2 − 𝐻1 )
𝑇𝑏 𝑃12 − 𝑃22 − 0.0375 ( )
𝑍𝑚 𝑇
𝑄 = 38.744 𝐹[ ] 𝑑 2.5 𝐸
𝑃𝑏 𝐺𝐸 ∗ 𝐿𝑇𝑍𝑚

8
Donde:
𝑄𝑏 =(𝑓𝑡 3 /día) flujo de gas a condiciones de base
𝑇𝑏 = (R) temperatura base
𝑃𝑏 = (psia) presión base
𝐹 = factor de transmisión
𝑃1 = (psia) presión de entrada
𝑃2 = (psia) presión de salida
𝑃𝑚 = (psia) presión promedio en la línea
𝐺𝐸 = gravedad especifica del gas
ℎ1 = (ft) elevación con referencia al punto de entrada
ℎ2 = (ft) elevación con referencia al punto de salida
𝑍𝑚 = factor de compresibilidad del gas
𝑇 = (R) temperatura promedio de la línea
𝐿 = (millas) longitud de la línea
𝑑 = (pulg) diámetro interno de la tubería
𝐸 = factor de eficiencia del ducto

VALOR DE ( E ) CONDICIONES DE LA TUBERIA

1.0 Completamente nueva
0.95 En buenas condiciones
0.92 En condición promedio
0.85 En condiciones no favorables

 ECUACIÓN GENERAL DEL FLUJO DE GAS

0.5
2 2 𝐺𝐸𝑃𝑚2 (𝐻2 − 𝐻1 )
𝑇𝑏 𝑃1 − 𝑃2 − 0.0375 ( )
𝑍𝑚 𝑇
𝑄 = 38.744 𝐹[ ] 𝑑 2.5 𝐸
𝑃𝑏 𝐺𝐸𝐿𝑇𝑍𝑚

De la ecuación anterior, se puede despejar diferentes variables, entre ellas las más
importantes:

 DIÁMETRO INTERNO DE LA TUBERÍA

1⁄
2.5

𝑄
𝑑= 0.5
2 2 𝐺𝐸𝑃𝑚2 (𝐻2 − 𝐻1 )
𝑃1 − 𝑃2 − 0.0375 ( )
𝑇𝑏 𝑍𝑚 𝑇
38.744 𝐹∙𝐸[ ]
𝑃𝑏 𝐺𝐸𝐿𝑇𝑍𝑚
{ }

9
  LONGITUD DEL DUCTO

2 2 𝐺𝐸𝑃𝑚2 (𝐻2 − 𝐻1 )
𝑇𝑏 𝑃1 − 𝑃
22 − 0.0375 ( )
𝑍𝑚 𝑇
𝐿 = 1501.09 ( 𝐹) ( ) 𝑑 4.5 𝐸 2
𝑃𝑏 𝐺𝐸𝑇𝑍𝑚 𝑄 2

 PRESIÓN DE ENTRADA

𝑇𝑏 2 2 𝐺𝐸𝑃𝑚2 (𝐻2 − 𝐻1 )
1501.09 ( 𝐹) [𝑃2 + 0.0375 ( )] + 𝐺𝐸𝑇𝑍𝑚 𝑄 2
𝑃𝑏 𝑍𝑚 𝑇
𝑃1 = √ 2 ∙ 𝑑 2.25 𝐸
𝑇𝑏
1501.09 ( ) 𝐺𝐸 𝑇 𝑍𝑚 𝑄 2
𝑃𝑏𝑓

 PRESIÓN DE SALIDA

𝑇𝑏 2 2 𝐺𝐸𝑃𝑚2 (𝐻2 − 𝐻1 )
1501.09 ( 𝐹) [𝑃1 − 0.0375 ( )] + 𝐺𝐸𝑇𝑍𝑚 𝑄 2
𝑃𝑏 𝑍𝑚 𝑇
𝑃2 = √ 2 ∙ 𝑑 2.25 𝐸
𝑇𝑏
1501.09 ( ) 𝐺𝐸 𝑇 𝑍𝑚 𝑄 2
𝑃𝑏𝑓

La ecuación general de flujo de gas en tuberías, es más conocida cuando se

considera que no existen cambios de nivel en el trayecto de la tubería, y en ese
caso, la ecuación general de flujo sería:
0.5
2 2 𝐺𝐸𝑃𝑚2 (𝐻2 − 𝐻1 )
𝑇𝑏 𝑃1 − 𝑃2 − 0.0375 ( 𝑍𝑚 𝑇 )
𝑄 = 38.744 𝐹[ ] 𝑑 2.5 𝐸
𝑃𝑏 𝐺𝐸𝐿𝑇𝑍𝑚

Donde:
(𝐻2 − 𝐻1 ) = 0

Tomando en cuenta el factor de fricción:

𝑇𝑏 𝑃12 − 𝑃22
𝑄 = 77.54 ( ) ( ) 𝑑 2.5 𝐸 (𝑢𝑛𝑖𝑑𝑎𝑑𝑒𝑠 𝑈𝑆𝐶𝑆)
𝑃𝑏 𝐺𝐸 𝑇𝑓 𝐿 𝑍 𝑓

Tomando en cuenta el factor de transición:

𝑇𝑏 𝑃12 − 𝑃22
𝑄 = 38.744 ( ) 𝐹 ( ) 𝑑 2.5 𝐸 (𝑢𝑛𝑖𝑑𝑎𝑑𝑒𝑠 𝑈𝑆𝐶𝑆)
𝑃𝑏 𝐺𝐸 𝑇𝑓 𝐿 𝑍

10
FACTORES RELEVANTES DE LA ECUACION GENERAL DE FLUJO
 EFECTO DE LA ELEVACION DE LA TUBERIA.
Cuando se incluye la diferencia de elevación entre los extremos de un segmento de
tubería, la Ecuación de flujo general se modifica de la siguiente manera:
0.5
𝑇𝑏 𝑃12 − 𝑒 5 𝑃22
𝑄 = 37.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 𝑒 𝑠 toma en cuenta la diferencia de
elevación entre los extremos aguas arriba y aguas abajo del segmento de tubería.
Los parámetros depende de la gravedad gas, factor de compresibilidad del gas, la
temperatura del 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, ft
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:

𝐻2 − 𝐻1
𝑠 = 0.0684 ∗ 𝐺𝐸 ( ) (𝑢𝑛𝑖𝑑𝑎𝑑𝑒𝑠 𝑆𝐼)
𝑇𝑓 ∗ 𝑍

11
Donde:
H1 = elevación de aguas arriba, m
H2 = elevación aguas abajo, m
 PRESION PROMEDIO
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 presión media en un segmento de tubería es:
2 𝑃13 − 𝑃23
𝑃𝑚 = ( )
3 𝑃12 − 𝑃22

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 convertidas a presión absoluta mediante la adición de la presión base.
 FACTOR DE FRICCION
El termino factor de fricción es un parámetro adimensional que depende del número
de Reynolds del flujo. En la literatura de ingeniería podemos encontrar 2 diferentes
factores de fricción: factor de fricción de Darcy y factor de fricción Fanning. La
relación que se tiene es la siguiente:
𝑓𝑑
𝑓𝑓 =
4
Donde:
𝑓𝑓 = Factor de fricción de Fanning

𝑓𝑑 = Factor de fricción de Darcy

Se utiliza el factor de fricción de Darcy (f), y se determina dependiendo del tipo de
flujo.

12
Para flujo laminar se utiliza la siguiente ecuación:
64
𝑓=
𝑅𝑒
Para flujo turbulento se utiliza la ecuación de Colebrook White. Que es una relación
entre el factor de fricción y el número de Reynolds, la rugosidad de la tubería y el
diámetro interno de la tubería. La ecuación es:
1 𝜀𝑟 2.51
= −2𝑙𝑜𝑔 ( + )
√𝑓 3.71 𝑅𝑒√𝑓

Esta ecuación se utiliza siempre y cuando el número de Reynolds sea mayor a 4000.
Donde:
f= Factor de fricción
D= Diámetro interno de la tubería (in)
𝜀𝑟 = Rugosidad absoluta (in)
Re = Numero de Reynolds

 FACTOR DE TRANSMISION
Este factor de transmisión F es considerado el opuesto del factor f. Mientras que
factor de fricción indica cuán difícil es mover una cierta cantidad de gas por la
tubería, el factor de transmisión es una medida directa de cuanto gas puede ser
transportado por la tubería. Si aumenta el factor de fricción, el factor de transmisión
disminuye.
La relación que existe para ambos es la siguiente:
2
𝐹=
√𝑓
Tomemos en cuenta que solo hay un factor de transmisión.
Si reemplazamos esta ecuación, en la ecuación de Colebrook White se modifica
por lo cual tenemos lo siguiente:
𝜀𝑟 1.255𝐹
𝐹 = −4𝑙𝑜𝑔 ( + )
3.71 𝑅𝑒

13
  VELOCIDAD DEL 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 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. Bajo
flujo de estado estacionario, en A, el caudal másico de gas se designa como M será
igual a la tasa de flujo por la densidad del fluido, si entre A y B no hay inyección o
entrega de gas por lo tanto en ambos puntos será el mismo caudal másico de gas.
Por lo tanto la masa es el producto será igual 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:
𝑄 = 𝑣 ∗ 𝐴
El flujo de masa M en la sección 1 y 2 es el mismo para el flujo en estado
estacionario. Por Consiguiente:
𝑀 = 𝑄1 𝜌1 = 𝑄2 𝜌2 = 𝑄𝑏 𝜌𝑏
Donde:
M = Flujo másico
𝑄𝑏 =es el caudal de gas en condiciones estándar
𝜌𝑏 = es la densidad del gas correspondiente.
Por lo tanto, la ecuación se simplifica a:
𝜌𝑏
𝑄1 = 𝑄𝑏 ( )
𝜌1
Pero teniendo en cuenta de la ecuación de estado si despejamos la densidad
tenemos lo siguiente:
𝑃1
𝜌1 =
𝑍1 𝑅𝑇1
𝑃𝑏
𝜌𝑏 =
𝑍𝑏 𝑅𝑇𝑏

14
Reemplazando en la ecuación del caudal 1 ( 𝑄1):
𝑃𝑏 𝑇1 𝑍1
𝑄1 = 𝑄𝑏 ( ) ( ) ( )
𝑃1 𝑇𝑏 𝑍𝑏
𝑍𝑏 es 1 aproximadamente
𝑃𝑏 𝑇1
𝑄1 = 𝑄𝑏 ( )( )𝑍
𝑃1 𝑇𝑏 1

Reemplazando entonces el caudal con: 𝑄1 = 𝑣1 𝐴

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

𝜋
Sabemos que el área es: 𝐴 = 4 𝐷2 reemplazando y llevándolo al sistema ingles

4𝑄𝑏 𝑍 𝑃𝑏 𝑇 144𝑖𝑛2 1𝑑𝑖𝑎 1ℎ𝑟 1 𝑚𝑖𝑛

𝑣= ( ) ( ) × × × ×
𝜋𝐷2 𝑃 𝑇𝑏 1𝑓𝑡 2 24ℎ𝑟 60 𝑚𝑖𝑛 60 𝑠𝑒𝑔
𝑄𝑏 𝑍 𝑃𝑏 𝑇
𝑣 = 0.002122 ( )( ) (𝑢𝑛𝑖𝑑𝑎𝑑𝑒𝑠 𝑈𝑆𝐶𝑆)
𝐷2 𝑃 𝑇𝑏
Donde:
v= Velocidad del gas, (ft/s)
Q= Caudal en condiciones estándar, (ft3/día)
D= Diámetro interno de la tubería, (in)
Pb = Presión base, (psia)
Tb = Temperatura base, °R (460 +°F)
Z= Factor de compresibilidad del gas, (adimensional).
T= Temperatura promedio del flujo de gas, °R (460 +°F)
P=Presión, (psia)

15
La velocidad en el sistema internacional es:
𝑄𝑏 𝑍 𝑃𝑏 𝑇
𝑣 = 14.7349 ( )( ) (𝑢𝑛𝑖𝑑𝑎𝑑𝑒𝑠 𝑆𝐼)
𝐷2 𝑃 𝑇𝑏
Donde:
v= Velocidad del gas, (m/s)
Q= Caudal en condiciones estándar, (m3s/día)
D= Diámetro interno de la tubería, (mm)
Pb = Presión base, (kPa)
Tb = Temperatura base, K (273+°C)
Z= Factor de compresibilidad del gas, (adimensional)
T= Temperatura promedio del flujo de gas, K (273+°C)
P=Presión, (kPa)

 VELOCIDAD EROSIONAL.

El límite máximo de la velocidad del gas es usualmente calculado aproximadamente

por la siguiente ecuación:
𝐶
𝑣𝑒 = (𝑢𝑛𝑖𝑑𝑎𝑑𝑒𝑠 𝑈𝑆𝐶𝑆)
√𝜌
Donde:
C = una constante empírica, que toma un valor de 100 para servicio continuo y 125
para servicio intermitente (adimensional).
ve= Velocidad erosionar de flujo, (ft/s)
𝜌 = Densidad delgas a temperatura de flujo, lb/ft3.
Tenemos también en unidades de sistema internacional lo cual vendría a ser:
𝐶
𝑣𝑒 = 1.22 (𝑢𝑛𝑖𝑑𝑎𝑑𝑒𝑠 𝑆𝐼)
√𝜌
Donde:
Ve =Velocidad erosional de flujo, (m/s)
C= Constante empírica, 100 –125, (adimensional)

16
𝜌 = Densidad del gas a temperatura de flujo, (kg/m3).
La velocidad erosional también puede ser expresada en términos de presión y
temperatura:

𝑍𝑅𝑇
𝑣𝑒 = 𝐶 √ (𝑢𝑛𝑖𝑑𝑎𝑑𝑒𝑠 𝑈𝑆𝐶𝑆)
29𝐺𝑃

Donde:
ve= Velocidad erosional, (pie/s)
Z= Factor de compresibilidad del gas a temperatura de flujo, (adimensional)
R= Constante universal de los gases, (10,73 psia pie3/lb mol °R)
T = Temperatura del gas, °R (460 +°F)
G= Gravedad especifica del gas, (adimensional)
P= Presión absoluta del gas, (psia)

En unidades de sistema internacional tenemos:

𝑍𝑅𝑇
𝑣𝑒 = 1.22𝐶 √ (𝑢𝑛𝑖𝑑𝑎𝑑𝑒𝑠 𝑆𝐼)
29𝐺𝑃

Donde:
ve= Velocidad erosional, (m/s)
Z= Factor de compresibilidad del gas a temperatura de flujo, (adimensional)
R= Constante universal de los gases, (8.314,4472 Nm/kmol K)
T= Temperatura del gas, K (273+°C)
G= Gravedad especifica del gas, (adimensional)
P= Presión absoluta del gas, (Pa).

17