Diseo de Gasoducto

El proceso que se requiere para transportar el gas natural del subsuelo hasta el
destino final es complicado, el movimiento efectivo y eficiente del gas natural desde
las regiones de produccin a las de consumo requiere un sistema de transporte
extensivo y elaborado.
El gas natural generalmente tiene que ser transportado por grandes distancias desde
las reas de produccin hasta las ciudades para su consumo o plantas de
procesamiento.

1. OBJETIVO GENERAL.
Realizar un modelo matemtico para un diseo de gasoducto.

2. OBJETIVOS ESPECFICOS.

Determinar el flujo de entrada.

Determinar el dimetro del ducto.
Determinar la longitud del ducto.
Determinar la presin de entrada y salida.
Anlisis de resultados.

3. MARCO TERICO.
3.1 Gasoducto.
Las tuberas que se utilizan en el transporte de gas natural son denominadas
gasoductos, estas pueden medir desde 6 a 48 pulgadas de dimetro, pero algunos
componentes de las tuberas de transporte.

3.2 Construccin de un gasoducto.

La construccin de un gasoducto requiere planeamiento y preparacin, estudios
extensivos de la ruta son completados, estudios areos y terrestres son realizados
para asegurar que no se presenten contratiempos durante el ensamblado del
gasoducto.
3.3 Ecuacin General de Flujo.
La Ecuacin General de Flujo tambin denominado Ecuacin Fundamental de Flujo
relaciona el caudal (Qg) a las condiciones estndar (60F, 14.7 psia), las propiedades
del gas, la longitud de la tubera (L) y la temperatura del flujo (T) existente entre las
presiones de entrada y de salida (P1, P2) de un segmento de tubera de dimetro
constante (di) como se muestra en la siguiente figura.
FIGURA 3.1 REPRESENTACION DEL FLUJO DE GAS EN UNA TUBERIA

FUENTE: GUA DE AUXILIATURA TRANSPORTE Y ALMACENAJE DE

HIDROCARBUROS II, 2014, Chvez Poma Pablo

La Ecuacin General de Flujo permite determinar la cada de presin debido a la

friccin en el interior de la tubera y explicar adecuadamente los cambios en las
propiedades de los gases debido a la presin y la temperatura para cualquier seccin
de la tubera.
La Ecuacin General de Flujo est representada por la siguiente ecuacin:

2.5

Q g =38.744 E d F

Tb
Pb

( )

( )

Pm
P P 0.0375
( H 1H 2 )
zm T
2
1

2
2

L zm T

0.5

( 3.1)

Dnde:
Qg = Flujo volumtrico del gas a condiciones estndar, PCD
E = Factor de eficiencia del ducto, fraccin
Tb = Temperatura estndar, R (Tb = 60F = 520 R)
Pb = Presin estndar, psia (Pb = 14.7 psia)
H1 = Elevacin respecto al punto de referencia a la entrada, pies
H2 = Elevacin respecto al punto de referencia a la salida, pies
P1 = Presin de entrada del ducto, psia
P2 = Presin de salida del ducto, psia
Pm = Presin promedio en la lnea, psia.
GE = Gravedad especfica del gas, adimensional
T = Temperatura promedio del flujo de gas, R
Le = Longitud del ducto, mi
Zm = Factor de compresibilidad del gas
di = Dimetro interior de la tubera, plg.

Tambin existen otras ecuaciones empricas para predecir la cada de presin o el

flujo de un gas en una tubera resultado de la experimentacin de varios autores.
La existencia de estas ecuaciones se derivan principalmente de la ecuacin general
de flujo, estos modelos permiten determinar el comportamiento en las distintas
condiciones del rgimen de flujo, entre las que se destacan son las siguientes:

Ecuacin de Weymouth
Ecuacin de Panhandle A
Ecuacin de Panhandle B
Ecuacin de Spitzglass

A continuacin solo se explican la ecuacin de Panhandle A, dado a que el presente

proyecto tiene como objeto realizar un diseo de Gasoducto con esta ecuacin.

3.4 Ecuacin de Panhandle A

La ecuacin de Panhandle A es adecuada para el clculo de sistemas tpicos de
distribucin de gas natural y gas licuado de petrleo, para tuberas de dimetros de
6 pulgadas a 24 pulgadas bajo nmeros de Reynolds de 5 a 14 millones.
Tambin es adecuado para sistemas de tuberas de gas de combustin de baja
presin ( 2 psig) para hornos y calderas Para una aproximacin del clculo en
regiones de flujo parcialmente turbulento se ajusta un factor de eficiencia de 0.90,
aunque la mayora de los calculistas tambin consideran un factor de eficiencia de
0.92. El porcentaje de exactitud de la ecuacin es mejor que la proporcionada por la
ecuacin de Weymouth (10%).
En esta ecuacin, la rugosidad de la tubera no es usada.

Qg =435.87 E d

2.6182

Tb
Pb

1.0778

( )

P2m
2
2
P1 P20.0375
( H 1H 2 )
zm T

( )

0.8538 L z m T

0.5394

(3.2)

Dnde:

Qg = Flujo volumtrico del gas a condiciones estndar, PCD

E = Factor de eficiencia del ducto, fraccin
Tb = Temperatura estndar, R (Tb = 60F = 520 R)
Pb = Presin estndar, psia (Pb = 14.7 psia)
H1 = Elevacin respecto al punto de referencia a la entrada, pies
H2 = Elevacin respecto al punto de referencia a la salida, pies
P1 = Presin de entrada del ducto, psia
P2 = Presin de salida del ducto, psia
Pm = Presin promedio en la lnea, psia.
GE = Gravedad especfica del gas, adimensional
T = Temperatura promedio del flujo de gas, R
Le = Longitud del ducto, mi
Zm = Factor de compresibilidad del gas
4

di = Dimetro interior de la tubera, plg.

De la ecuacin general de Panhandle A, se obtiene:
Dimetro del Gasoducto.

d=

Qg

2.6182

435.87 E

Tb
Pb

1.0778

( )

P21 P220.0375

2
m

( )

P
( H 1H 2 )
zm T

0.8538 L z m T

0.5394

(3.3)

Longitud del Gasoducto.

P21P22 0.0375
L=

P 2m
( H 1H 2 )
zm T

( )

Qg

0.5394

435.87 E d 2.6182

Tb
Pb

1.0788

( )

(3.4)

( 0.8538 z m T )

Presin de entrada.
2
1

Qg

P =0.5394

435.87 E d

2.6182

0.8538

Tb
Pb

1.0788

( )

P2m
L z m T ) +0.0375
H 1 H 2) + P 22( 3.5)
(
zm T

( )

Presin de salida.

P2= P10.5394

Qg
435.87 E d 2.6182

Tb
Pb

1.0788 (

( )

0.8538

L z m T ) 0.0375

( )

Pm
( H 1H 2 ) (3.6)
zmT

3.5 Carga Esttica.

Este parmetro de correccin se debe a que en la mayora de los casos la diferencia
de altura es considerada despreciable ya que la densidad de los vapores y gases es
tan baja que incluso cambios de elevacin de 200 ft, de altura representan solo una
pequea variacin de presin.
s=0.0375

P2m
( H 2 H 1 ) (3.7)
zm T

( )

Dnde:
s = Carga esttica.
H1 = Elevacin respecto al punto de referencia a la entrada, pies
H2 = Elevacin respecto al punto de referencia a la salida, pies
Pm = Presin promedio, psia
T = Temperatura promedio del flujo de gas, R
Zm = Factor de compresibilidad

3.6 Temperatura Promedio.

Esta dada por la siguiente ecuacin:
T=

T 1 +T 2
(3.8)
2

Dnde:
T = Temperatura promedio del flujo de gas, R
T1 = Temperatura de ingreso, R
T2 = Temperatura de salida, R

3.7 Presin Promedio.

La presin promedio del ducto est dada por la siguiente ecuacin.

P m=

P P2
2
P 1+ P 2) 1
(3.9)
(
3
P 1+ P 2

Dnde:
Pm = Presin promedio en la lnea, psia
P1 = Presin de entrada del ducto, psia
P2 = Presin de salida del ducto, psia
3.8 Factor de Compresibilidad z.
El clculo del valor del factor de compresibilidad est en funcin de la Presin
psudoreducida y Temperatura pseudoreducida..
z=f ( P pr , T pr )
Dnde:
P pr=

P
(3.10)
P pc

T pr =

T
(3.11)
T pc

Dnde:
T pc =168+325 g12.5 2g (3.12)
2

P pc =677+15.0 g 37.5 g (3.13)

Dnde:

= Gravedad especifica del gas

Para determinar el factor de compresibilidad existen varias correlaciones como ser:

Peng Robinson, Yarborough - Hall, Dranchuk-Abu-Kassem, etc. No obstante para el
desarrollo del presente proyecto se empleara el mtodo de Yarborough Hall.
3.8.1 Mtodo de Yarborough Hall.
Hall y Yarborough (1973) presentaron una ecuacin de estado que representa
exactamente la grfica del factor z de Standing y Katz. La expresin propuesta se
basa en las ecuaciones de estado de Starling-Carnahan. Los coeficientes de la
correlacin fueron determinados ajustndolos a los datos tomados de la grfica del
factor z de Standing y Katz. Hall y Yarborough propusieron la siguiente frmula
matemtica:
0.06125 P pr t exp [ 1.2 ( 1t )2 ]
z=
(3.14)
Y
Dnde:
Ppr = Presin pseudo reducida
t

= Reciproco de la temperatura pseudo reducida, 1/T pr

Y = densidad reducida que puede obtenerse como solucin de la siguiente

ecuacin:
2

Y +Y +Y y
2
D
F ( Y )= A P pr +
B Y +C Y (3.15)
4
( 1Y )
Dnde:

A=0.06125 t exp [1.2 ( 1t )2 ](3.16)

B=14.76 t9.76 t 2+ 4.58 t 3 (3.17)
C=90.7 t242.2 t 2 +42.4 t 3 (3.18)

D=2.18+ 2.82t (3.19)

La ecuacin 3.14 es una ecuacin no lineal y puede ser convenientemente resuelta
para la densidad reducida Y usando la tcnica de iteracin de Newton-Raphson.
Aplicando Newton Raphson.
Y cal=Y as

F (Y )
(3.20)
F ' (Y )

La derivada de la Funcin.
F ' ( Y )= A P pr +

1+4 Y + 4 Y 2 + 4 Y 3 y 4
2
( D1 )
2 B Y + CD Y
(3.21)
4
(1Y )

3.9 Viscosidad del Gas.

La viscosidad del gas es la medida de la friccin interna del fluido o resistencia al
flujo que afecta a la cada de presin por influjo del reservorio al agujero del pozo y a
lo largo de las instalaciones.
La viscosidad de un fluido generalmente se define como la relacin de la fuerza
distribuida por unidad de rea al gradiente de viscosidad local. Las viscosidades se
expresan en trminos de poises, centipoises o micro-poises.
La viscosidad del gas comnmente no se mide en laboratorio porque puede
estimarse con precisin de correlaciones empricas. Como todas las propiedades
intensivas, la viscosidad del gas natural es descrita completamente por la siguiente
funcin:
g =f ( P , T , Y i )
La relacin anterior simplemente establece que la viscosidad es una funcin de la
presin, temperatura, y composicin.

Existen dos mtodos populares que son comnmente usados en la industria

petrolera, los cuales son:

Mtodo de Lee-Gonzalez-Eakin
Mtodo de Carr-Kobayashi-Burrows

Para el desarrollo del presente proyecto se empleara el mtodo de Lee-GonzalezEakin

3.9.1 Mtodo de Lee-Gonzalez-Eakin
La relacin anterior simplemente establece que la viscosidad es una funcin de la
presin, temperatura, y composicin.
Lee, Gonzalez Eakin (1966) presentaron una relacin semi-emprica para calcular la
viscosidad de los gases naturales. Los autores expresaron la viscosidad del gas en
trminos de la temperatura del reservorio, densidad del gas y el peso molecular del
gas.
Su ecuacin propuesta est dada por:
Y

g =104 K e X (3.22)
Dnde:
K=

( 9.4 +0.02 M a ) T 1.5

209+ 19 M a T

X =3.5+

(3.23)

986
+0.01 M a (3.24)
T

Y =2.40.2 X ( 3.25)

g = Densidad del gas a presin y temperatura del reservorio, lb/pc

Ma = Peso molecular aparente, lb/lb-mol
La densidad del gas esta dad por:

10

g =

2.7 g P m
(3.26)
T zm

Dnde:
g = Gravedad especifica del gas
Pm = Presin promedio, psia
T = Temperatura promedio del flujo de gas, R
Zm = Factor de compresibilidad

El peso molecular aparente:

M a= g ( 28.96 ) (3.27)

3.10 Nmero de Reynolds.

Re =0.0004778

Pb
Tb

( )( )

Q
(3.28)
d

Dnde:
Re = Nmero de Reynolds,
Pb = Presin estndar, psia.
Tb = Temperatura estndar, R.
GE = Gravedad especifica del gas.
Q = Flujo del Gas, condiciones estndar, PCD.
d = Dimetro interno, plg.
= Viscosidad del gas, cp.

11

3.11 Factor de Friccin.

Colebrook y While en 1939, propusieron una ecuacin aplicable a tuberas lisas
como tubera de flujo de transicin y totalmente rugosas en las zonas de flujo
turbulento.
F=

0.0768
(3.28)
0.1461
Re

Dnde:
F = Factor de Friccion.
Re = Numero de Reynods

4. ALGORITMO
MEN.
1.
2.
3.
4.
5.

Caudal de gas
Dimetro del Gasoducto
Longitud del Gasoducto.
Presin de Entrada del Gas P1
Presin de salida del Gas P2

4.1 Caudal de Gas

Paso 1. Inicio
Paso 2. Crear una pestaa para el caudal de gas, en la hoja de Excel.
Paso 3. Leer Pb, Tb, GE, L, P1, P2, H1, H2, P, d, E.
Paso 4. Leer T1 y T2 T promedio.
Paso 5. Si se tiene T1 y T2, calcular la temperatura promedio, ecuacin (3.8)
Paso 6. Calcular la presin promedio, ecuacin (3.9)
Paso 7. Calcular Ppc y Tpc, ecuaciones (3.12) y (3.13) respectivamente
12

Paso 8. Calcular Ppr y Tpr, con presin y temperatura promedio.

Paso 9. Calcular el factor z, subrutina para el clculo del factor compresibilidad del
gas z, mtodo Yaoboruogh-Hall.
Paso 10. Calcular el caudal de gas, mtodo de Panhandle A. ecuacin (3.2)
Paso 11. Calcular la viscosidad del gas, subrutina para el clculo de la viscosidad
del gas z, mtodo Lee Gonzales Eiken.
Paso 12. Calcular nmero de Reynolds
Paso 13. Calcular factor de friccin.
Paso 14. Calcular caudal de gas con el factor de friccin incluido.
Paso 15. Fin.
4.2 Dimetro del Gasoducto.
Paso 1. Inicio
Paso 2. Crear una pestaa para el dimetro del gasoducto, en la hoja de Excel.
Paso 3. Leer Pb, Tb, GE, Qg, P1, P2, H1, H2, P, L, E.
Paso 4. Leer T1 y T2 T promedio.
Paso 5. Si se tiene T1 y T2, calcular la temperatura promedio, ecuacin (3.8)
Paso 6. Calcular la presin promedio, ecuacin (3.9)
Paso 7. Calcular Ppc y Tpc, ecuaciones (3.12) y (3.13) respectivamente
Paso 8. Calcular Ppr y Tpr, con presin y temperatura promedio.
Paso 9. Calcular el factor z, subrutina para el clculo del factor compresibilidad del
gas z, mtodo Yaoboruogh-Hall.
Paso 10. Despejas d, de la ecuacin de general de Panhandle A, ecuacin (3.2)
Paso 11. Calcular el dimetro del gasoducto, ecuacin (3.3).

13

Paso 12. Calcular la viscosidad del gas, subrutina para el clculo de la viscosidad
del gas z, mtodo Lee Gonzales Eiken.
Paso 13. Calcular nmero de Reynolds
Paso 14. Calcular factor de friccin.
Paso 15. Calcular el dimetro del gasoducto con el factor de friccin incluido.
Paso 16. Fin.
4.3 Longitud del Gasoducto.
Paso 1. Inicio
Paso 2. Crear una pestaa para la longitud del gasoducto, en la hoja de Excel.
Paso 3. Leer Pb, Tb, GE, Qg, P1, P2, H1, H2, P, d, E.
Paso 4. Leer T1 y T2 T promedio.
Paso 5. Si se tiene T1 y T2, calcular la temperatura promedio, ecuacin (3.8)
Paso 6. Calcular la presin promedio, ecuacin (3.9)
Paso 7. Calcular Ppc y Tpc, ecuaciones (3.12) y (3.13) respectivamente
Paso 8. Calcular Ppr y Tpr, con presin y temperatura promedio.
Paso 9. Calcular el factor z, subrutina para el clculo del factor compresibilidad del
gas z, mtodo Yaoboruogh-Hall.
Paso 10. Despejas L, de la ecuacin de general de Panhandle A, ecuacin (3.2)
Paso 11. Calcular la longitud del gasoducto, ecuacin (3.4).
Paso 12. Calcular la viscosidad del gas, subrutina para el clculo de la viscosidad
del gas z, mtodo Lee Gonzales Eiken.
Paso 13. Calcular nmero de Reynolds
Paso 14. Calcular factor de friccin.

14

Paso 15. Calcular longitud del gasoducto con el factor de friccin incluido.
Paso 16. Fin.
4.4 Presin de Entrada del Gas P1
Paso 1. Inicio
Paso 2. Crear una pestaa para la de salida del gas, en la hoja de Excel.
Paso 3. Leer Pb, Tb, GE, Qg, L, P2, H1, H2, P, d, E.
Paso 4. Leer T1 y T2 T promedio.
Paso 5. Si se tiene T1 y T2, calcular la temperatura promedio, ecuacin (3.8)
Paso 6. Calcular Ppc y Tpc, ecuaciones (3.12) y (3.13) respectivamente
Paso 7. Despejar P1, de la ecuacin de general de Panhandle A, ecuacin (3.2)
Paso 8. Calcular la presin promedio, asumiendo un valor de P 1 ecuacin (3.9)
Paso 9. Calcular Ppr y Tpr, con presin y temperatura promedio.
Paso 10. Calcular el factor z, subrutina para el clculo del factor compresibilidad del
gas z, mtodo Yaoboruogh-Hall.
Paso 11. Calcular presin de entrada del gas, ecuacin (3.5).
Paso 12. Calcular la diferencia entre P1(cal) y P1(as). Si |P1(cal) P1(as)| 1 psia se
procede al siguiente paso. Si no P1(cal) = P1(as), y volver al paso 8
Calcular la viscosidad del gas, subrutina para el clculo de la viscosidad del gas z,
mtodo Lee Gonzales Eiken.
Paso 13. Calcular nmero de Reynolds
Paso 14. Calcular factor de friccin.
Paso 15. Calcular presin de entrada del gas con el factor de friccin incluido.
Paso 16. Fin.

15

4.5 Presin de salida del Gas P2

Paso 1. Inicio
Paso 2. Crear una pestaa la presin de salida del gas, en la hoja de Excel.
Paso 3. Leer Pb, Tb, GE, Qg, L, P1, H1, H2, P, d, E.
Paso 4. Leer T1 y T2 T promedio.
Paso 5. Si se tiene T1 y T2, calcular la temperatura promedio, ecuacin (3.8)
Paso 6. Calcular Ppc y Tpc, ecuaciones (3.12) y (3.13) respectivamente
Paso 7. Despejar P2, de la ecuacin de general de Panhandle A, ecuacin (3.2)
Paso 8. Calcular la presin promedio, asumiendo un valor de P 2 ecuacin (3.9)
Paso 9. Calcular Ppr y Tpr, con presin y temperatura promedio.
Paso 10. Calcular el factor z, subrutina para el clculo del factor compresibilidad del
gas z, mtodo Yaoboruogh-Hall.
Paso 11. Calcular presin de entrada del gas, ecuacin (3.6).
Paso 12. Calcular la diferencia entre P2(cal) y P2(as). Si |P2(cal) P2(as)| 1 psia se
procede al siguiente paso. Si no P2(cal) = P2(as), y volver al paso 8
Calcular la viscosidad del gas, subrutina para el clculo de la viscosidad del gas z,
mtodo Lee Gonzales Eiken.
Paso 13. Calcular nmero de Reynolds
Paso 14. Calcular factor de friccin.
Paso 15. Calcular presin de salida del gas con el factor de friccin incluido.
Paso 16. Fin.
4.6 Subrutina para el clculo de z, mtodo de Yamborough-Hall
Paso 1. Inicio.

16

Paso 2. Leer o calcular Ppr y Tpr.

Paso 3. Invertir la temperatura pseudo reducida, t = 1/Tpr
Paso 4. Asumir primer valor para Y, Y(as) = 0.001
Paso 5. Calcular los coeficientes A, B, C y D; ecuaciones (3.16), (3.17), (3.18), y
(3.19) respectivamente.
Paso 6. Calcular la funcin con Y (asumido), ecuacin (3.15)
Paso 7. Calcular la derivada de la funcin Y (asumido), ecuacin (3.21)
Paso 8. Calcular Y(cal), ecuacin (3.20)
Paso 9. Calcular la diferencia entre Y(cal) y Y(as). Si |Y(cal) Y(as)| 10-6 se procede al
siguiente paso. Si no Y(cal) = Y(as), y volver al paso 5
Paso 10 Calcular z, ecuacin (3.14)
Paso 11. Fin
4.7 Subrutina para el clculo de la viscosidad de gas, mtodo de Lee
Gonzales Eiken
Paso 1. Inicio.
Paso 2. Leer presin, temperatura y gravead especifica del gas (GE)
Paso 3. Calcular, Ppc, Tpc, ecuaciones (3.12) y (3.13), respectivamente.
Paso 4. Calcular Ma, ecuacin (3.27).
Paso 5. Calcular Ppr y Tpr ecuaciones (3.10) y (3.11) respectivamente.
Paso 6. Calcular el factor z, subrutina para el clculo del factor compresibilidad del
gas z, mtodo Yaoboruogh-Hall.
Paso 7. Calcular la densidad del gas, ecuacin (3.26)
Paso 8. Calcular los coeficientes K, X, Y, ecuaciones (3.23), (3.24) y (3.25)
respectivamente.
17

Paso 9. Calcular la viscosidad del gas, ecuacin (3.22)

Paso 10. Fin.

5. APLICACIN PRCTICA.
5.1 Datos de Entrada
Los datos con los que se trabaj para el presente proyecto son los siguientes 1:

GE =0,6
Tb = 60
F
Pb = 14,7 psia
T = 100 F
L = 20
millas
P1 = 2000 psia
P2 = 1500 psia
d = 4,026 plg
E = 1
H = 100 pies
Q = 16691560.12 pcs/d

5.2 Programa Excel.

En primer lugar se cre las pestaas en la hoja de Excel.

1 Anexo I
18

Men para el diseo de Gasoducto

Con los datos de entrada, se procedi a los siguientes clculos generales para
determinar el caudal de gas, dimetro del gasoducto y longitud del gasoducto.
Obteniendo los siguientes resultados:
PRESIN
PROMEDIO
Pm =

1761,905 psia

RESULTADOS
PARCIALES
Ppc =
672,500
Tpc =
358,500
M=
17,376 lb/lb-mol

Ppr =
Tpr =

2,61993
1,56206

FACTOR DE COMPRESIBILIDAD z, calculado por la

subrutina del factor z
z=
0,81885
VISCOSIDAD DEL GAS,
calculado por la subrutina
de viscosidad.
=
0,015293938

Subrutina para el clculo de z, mtodo de Yamborough-Hall

CALCULO DEL FACTOR "z" POR EL METODO DE
YARBOROUGH-HALL
DATOS INGRESADOS
Ppr
2,61993
Tpr
1,56206
Pr
0,001

19

RESULTADOS
PARCIALES
t = 1/Tpr
0,64018
A
0,03357
B
6,65073
C
-30,07208
D
3,98530
ITER
1
2
3
4
5
6
7
8
9
10

Y(as)
0,001
0,08841
0,10730
0,10740
0,10740
0,10740
0,10740
0,10740
0,10740
0,10740

F(Yas)
-0,08695
-0,01398
-0,00007
0,00000
0,00000
0,00000
0,00000
0,00000
0,00000
0,00000
z=

F'(Yas)
0,99473
0,73977
0,73516
0,73519
0,73519
0,73519
0,73519
0,73519
0,73519
0,73519

Ycal
0,08841
0,10730
0,10740
0,10740
0,10740
0,10740
0,10740
0,10740
0,10740
0,10740

(Yas-Ycal)
8,74108E-02
1,88932E-02
9,92387E-05
-1,73624E-09
0,00000E+00
0,00000E+00
0,00000E+00
0,00000E+00
0,00000E+00
0,00000E+00

0,81885

Subrutina para el clculo de la viscosidad de gas, mtodo de Lee Gonzales

Eiken
CALCULO DE LA VISCOCIDAD DEL
GAS POR EL METODO DE LEE GONZALES
g =

0,099751444

K=
X=
Y=

117,5226517
5,434474286
1,313105143

0,015293938
20

5.2.1 Caudal de Gas

Datos de entrada.
GE =
Tb
Pb
T=
L=

0,6
60
14,7
100
20

F
psia
F
millas

P1 =

2000 psia

P2 =
=

1500 psia
cp

H1 =
H2 =
d=
E=
H
=

520 R
560 R

100
200
4,026 plg
1
100

Resultados parciales para determinar el caudal de gas

Qg =435.87 E d

2.6182

Tb
Pb

1.0778

( )

P2m
2
2
P1 P20.0375
( H 1H 2 )
zm T

( )

0.8538 L z m T
B

0.5394

(3.2)

s
P12-P22*s*H GE0,8538*L*T*z
B/C
152,3188709 1734768,113 5929,378828 292,5716442

E
D0,5394

F
435,87*(Tb/Pb)1,0778*d2
,1682
*E

21,39368841

780209,5554

Caudal de gas
21

G
E*F
16691560,1
2

16691560
Q=
,12 pcs/d
Q=
16,69 MMpcs/d
Caudal de gas, incluyendo factor de friccin
NUMERO DE
REYNOLDS
Re =

2196,924
914

F=

0,024952
957

FACTOR DE
FRICCIN

Resultados parciales para determinar el caudal de gas

A

s
P12-P22*s*H GE0,8538*L*T*z
B/C
152,3188709 1734768,113 5929,378828 292,5716442

E
D0,5394

F
435,87*(Tb/Pb)1,0778*d2
,1682
*E

21,39368841

19468,53545

G
E*F
416503,781
1

Caudal de gas.

Q=
Q=

416503,78
11 pcs/d
0,417 MMpcs/d

5.2.2 Dimetro del Gasoducto.

Datos de entrada.
GE =
Tb
Pb

0,6
60 F
14,7 psia
22

520 R

T=
L=

100 F
20 millas

P1 =

2000 psia

P2 =

1500 psia

H1 =

100

H2 =

200
16691560,
12 pcs/d
1
100

Q=
E=
H =

560 R

MMpcs/
16,69 d

Resultados parciales para determinar el dimetro del gasoducto.

d=

Qg

2.6182

435.87 E

Tb
Pb

1.0778

( )

P2m
2
2
P1 P20.0375
( H 1H 2 )
zm T

( )

0.8538 L z m T

0.5394

B
C
2
0,8538
P1 GE
*L*
s
P22*s*H
T*z
152,31887 1734768,1 5929,3788
09
13
28

(3.3)

E
D0,5394
21,393688
41

F
435,87*(Tb/Pb)1,0
778
*E
20348,42308

D
B/C
292,57164
42

E*F

Qb*G

435327,823

4,026

Dimetro del gasoducto.

d=

4,026 plg.

Dimetro del gasoducto, incluyendo factor de friccin

NUMERO DE
REYNOLDS
23

Re =

2638,9028
54

F=

0,0242935
66

FACTOR DE
FRICCION

Dimetro del gasoducto.

d = 16,654 plg
5.2.3 Longitud del Gasoducto.
Datos de entrada.
GE =
Tb
Pb
T=
Q=

0,6
60
14,7
100
16691560
,12

F
psia
F

560 R

pcs/d

P1 =

2000 psia

P2 =

1500 psia

H1 =

520 R

100

H2 =
d=
E=
H =

200
4,026 plg
1
100

Resultados parciales para determinar la longitud del gasoducto.

2

L=

( )

Pm
P P 0.0375
( H 1H 2 )
zm T
2
1

2
2

Qg

0.5394

Tb
Pb

1.0788

435.87 E d 2.6182

( )

A
s

B
P12-

(3.4)

( 0.8538 z m T )

C
GE

0,8538

*T*z

24

D
435,87*(Tb/Pb)1,0778*d2,16

P22*s*H
152,318870 1734768,1
296,4689414
9
13

82

*E

780209,5554

E
F
G
(1/0,5364)
Qb/D
E
C*F
21,393688
86738,40
41
292,5716442
56

H
B/G
20

Longitud del gasoducto.

L=

20,0 millas

Longitud del gasoducto, incluyendo factor de friccin

NUMERO DE
REYNOLDS
Re =

2638,90285
4

F=

0,02429356
6

FACTOR DE
FRICCIN

Longitud del gasoducto.

L=

0,0203 millas

5.2.4 Presin de entrada del Gas.

Datos de entrada.
GE =
Tb
Pb
T=
Q=
L=
P2 =

0,6
60
14,7
100
1669156
0

F
psia
F
pcs/d

20 millas
1500 psia
25

520 R
560 R

H1 =

100

H2 =
d=
E=
H =

200
4,026 plg
1
100

Resultados parciales para determinar la presin de entrada del gas.

Qg

P21=0.5394

435.87 E d 2.6182

Tb
Pb

0.8538 L z m T ) +0.0375
1.0788 (

( )

A
435,87*(Tb/Pb)1,0778*d2
,1682
*E
780209,5554

P2m
2
( H 1 H 2) + P 2( 3.5)
zm T

( )

Qb/A

B(1/0,5364)

21,393688
41

292,571644
2

Determinacin de la presin de entrada del gas por prueba y error

ITER
1
2
3
4
5
6
7
8
9
10

P1 (as)

Pm

Ppr

P1 (cal)

|P1 (cal)-P1
(as)|

1550,537
6
2,305632 0,831762 2005,924
405,92410
1765,130
2005,924
0
2,624729 0,818682 1999,925
5,99948
1761,863
1999,925
7
2,619872 0,818855 2000,001
0,07634
1761,905
2000,001
3
2,619934 0,818852 2000,000
0,00097
1761,904
2000,000
8
2,619933 0,818852 2000,000
0,00001
1761,904
2000,000
8
2,619933 0,818852 2000,000
0,00000
1761,904
2000,000
8
2,619933 0,818852 2000,000
0,00000
1761,904
2000,000
8
2,619933 0,818852 2000,000
0,00000
1761,904
2000,000
8
2,619933 0,818852 2000,000
0,00000
1761,904
2000,000
8
2,619933 0,818852 2000,000
0,00000
Tabla de las iteraciones para determinar la presin de entrada del gas.
1600,000

26

 Presin de entrada del gas

P1 =

2000,0 psia

Presin de entrada, incluyendo factor de friccin

NUMERO DE
REYNOLDS
Re =

2638,90285
4

F=

0,02429356
6

FACTOR DE
FRICCIN

ITER

P1 (as)

Pm

Ppr

P1 (cal)

|P1 (cal)-P1
(as)|

1550,54

2,305632

0,831762

41675,31

40075,30732

1600,00
41675,3
1

27818,28

41,365473

0,818682

2789,67

38885,63459

2789,67

2209,46

3,285441

0,803362

1994,08

795,58973

1994,08

1758,69

2,615146

0,819023

2000,08

5,99270

2000,08

1761,95

2,619994

0,818850

2000,00

0,07666

2000,00

1761,90

2,619932

0,818852

2000,00

0,00098

2000,00

1761,90

2,619933

0,818852

2000,00

0,00001

2000,00

1761,90

2,619933

0,818852

2000,00

0,00000

2000,00

1761,90

2,619933

0,818852

2000,00

0,00000

10

2000,00

1761,90

2,619933

0,818852

2000,00

0,00000

Presin de entrada del gas

P1 =

2000,0 psia

27

5.2.5 Presin de salida del Gas.

Datos de entrada.
GE =
Tb
Pb
T=

0,6
60
14,7
100
1669156
0
2000
20
8,12E-06
100

Q=
P1 =
L=
=
H1 =
H2 =
d=
E=
H =

F
psia
F

520 R
560 R

pcs/d
psia
psia
cp

200
4,026 plg
1
100

Resultados parciales para determinar la presin de salida del gas.

Qg

P2= P10.5394

435.87 E d 2.6182

0.8538

Tb
Pb

1.0788 (

( )

L z m T ) 0.0375

A
435,87*(Tb/Pb)1,0778*d2
,1682
*E
780209,5554

P 2m
( H 1H 2 ) (3.6)
zmT

( )

Qb/A

B(1/0,5364)

21,393688
41

292,57164
42

Determinacin de la presin de entrada del gas por prueba y error

ITER

P1 (as)

Pm

Ppr

P1 (cal)

1
2

1900,000
1504,810

1950,4274
1764,0659

2,900264
2,623146

0,810316
0,818738

1504,810
1500,067

28

|P1 (cal)-P1
(as)|
395,18971
4,74288

3
1500,067 1761,9350 2,619978 0,818851 1500,001
4
1500,001 1761,9052 2,619933 0,818852 1500,000
5
1500,000 1761,9048 2,619933 0,818852 1500,000
6
1500,000 1761,9048 2,619933 0,818852 1500,000
7
1500,000 1761,9048 2,619933 0,818852 1500,000
8
1500,000 1761,9048 2,619933 0,818852 1500,000
9
1500,000 1761,9048 2,619933 0,818852 1500,000
10
1500,000 1761,9048 2,619933 0,818852 1500,000
Tabla de las iteraciones para determinar la presin de entrada del gas.

0,06647
0,00093
0,00001
0,00000
0,00000
0,00000
0,00000
0,00000

Presin de salida del gas

P2 =

1500,0

psia

5.3 Anlisis de resultados.

Se puede observar que este factor afecta en gran proporcin a los valores
Resultados:

Caudal
Diametro

Sin factor Friccion

16 MMpcs/d
4.026 plg

Con factor Friccion

0.417 MMpcs/d
16.65 plg.

Longitud

20 millas

0.020 millas

Presion de entrada

2000 psia

2000 psia

Presion de salida

1500 psia

error

En la tabla anterior se puede notar la diferencia que existe entre los resultados, esta
variacin de resultados dependera de diferentes razones entre las cuales:
El factor de friccin afecta al resultado de una manera drstica.
La ecuacin de factor de Friccin, no sea la apropiada.
La ecuacin para determinar el nmero de Reynolds no es la apropiada.

29

Para el clculo de la presin de salida con el factor de friccin, se observo que los
valores son muy altos y el momento de resolver la raz, se produce un error en el
programa.

6. CONCLUSIONES Y RECOMENDACIONES.
Con base a los resultados obtenidos en el trabajo se concluye que la metodologa
empleada en el diseo de gasoductos es adecuada. No obstante, se recomienda
efectuar el uso de ecuaciones apropiadas para determinar algunos parmetros, tales
como el factor de friccin.
Adems, la informacin contenida en el trabajo como los criterios de diseo y las
opciones de clculo facilita en gran medida a la solucin de problemas relacionados
al transporte de gas. Siendo de este modo, una herramienta alternativa de clculo de
gran utilidad, ya que permite de alguna forma comprender, analizar y determinar el
comportamiento del gas en el ducto..
Para las contribuciones a trabajos a futuros, aprovechando la base de datos en
Excel sera interesante continuar con su desarrollo, con la finalidad de ampliar sus
capacidades e integrar mayor informacin al programa. Por ejemplo; la opcin de
poder predecir propiedades fsicas, disear un gasoducto a partir de la composicin
del gas, etc mediante el uso de ecuaciones de estado y modelos correlativos,
ampliando al mismo tiempo la base de datos del programa.
30

31