Diseo de Gasoductos

Generalidades
Configuracin (Layout) Hidrulica Materiales Otros Evaluacin Econmica

Configuracin
Ducto Vlvulas Estaciones de Compresin Medidores y reguladores Bridas aislantes Cruces Proteccin catdica Recubrimiento Telecomunicaciones
Trampa de lanzamiento Trampa de recepcin

Cruces

Hidrulica (Flujo de Fluidos)

Teorema de Bernoulli H= Ze + 144 P + V2 2g Ze1 + 144 P1 + V12 = Ze2 + 144 P2 + V22 + HL 1 2g 2 2g
A lo largo del flujo de un fluido se producen cadas de presin, que slo pueden ser explicadas adecuadamente conociendo las propiedades fsicas del mismo. Las propiedades que ms frecuentemente intervienen en la solucin de problemas de flujo de fluidos son la viscosidad y la densidad La viscosidad es dependiente de la temperatura. La viscosidad en la mayora de los lquidos decrece con el incremento de la temperatura mientras que la de los gases se incrementa. La presin casi no tiene efecto en la viscosidad de lquidos o gases perfectos La gravedad especfica es el inverso del volumen especfico. La gravedad especfica es la relacin entre la densidad de un lquido con la densidad del agua a 60F. En un gas es la relacin entre el peso molecular del gas y el del aire. Densidad relativa: MW/28.9625

Prdidas de friccin

Las Prdidas de Presin

La friccin resulta en una prdida de energa disponible para trabajo Ecuacin de Darcy-Weisbach: hL= (fm)(L)(V2) (pies) hL= ()(fm)(L)(V2) (lb/pulg2) (144) D (2g) 2gD
fm factor de friccin de Moody ff factor de Fanning=0.25 fm

La ecuacin de Darcy-Weisbach es vlida para flujo laminar y turbulento de cualquier lquido. La ecuacin debe aplicarse en segmentos de lnea suficientemente cortos de manera que la densidad del fluido sea constante. Flujo laminar: (Poiseuille's law) fm= 64/Re o ff= 16/Re fm=64 e = 64 12 Pf = 0.000668 LV DV V 1488 d d2 Flujo turbulento: (Colebrook) 1 = - 2 log10 ( + 2.51 ) /D rugosidad relativa fm 3.7D Refm viscosidad , cp

Factores de Friccin
crtica laminar turbulenta rugosas

Lisas

Re Reynold Number = DV

La Prdida de Presin en gases

Una diferencia considerable entre los gases y los lquidos es el elevado cambio en el volumen que puede lograrse en los gases, usualmente definido mediante el factor de compresibilidad (Z)

Z= 1 / (Fpv)2

Fpv = 1 + (avg)(3.444)(105)(10 (1.785) (S))

Tavg 3.825
S gravedad especfica (aire=1)

FLUJO EN UNA SOLA FASE

Flujo Isotrmico Q= 38.77 Tb E Pb Factor de eficiencia 1 ff P12 - P22 S Lm Tavg Zavg
0.5

d 2.5

pies cbicos/da

Tb = 520R Pb = 14.73 psia

En la ecuacin general, todas las irreversibilidades y no idealidades, excepto las cubiertas por leyes expresas, se incluyen en el factor de friccin

Las Ecuaciones AGA (American Gas Asociation)

Flujo Totalmente turbulento
1 = 4 log10 3.7 D ff

Flujo parcialmente turbulento

1 = 4 log10 Re ff 1 ff - 0.6

La Ecuacin Weymouth (1912)

ff = 0.008 d1/3 1 = 11.18 d 1/6 ff Gran aproximacin en ductos cortos y sistemas de distribucin El grado de error se incrementa con la presin Una forma de corregir el caudal es : Q final = Q calculado . 1/ Z

Panhandle A (1940) Panhandle Eastern Pipeline Co

Re = 1.934 QS Q= pc/d d d= pulg 1 = 7.211 (QS) 0.07305 ff d 5'000,000 < Re < 11'000,000 Basado en experiencias medidas Tuberas lisas Recomendacin: E= 0.90 - 0.92 Es menos precisa con flujos altos

Panhandle B (1956) Panhandle Eastern Pipeline Co

1 = 16.70 (QS) 0.01961 ff d El valor E en la prctica es un factor de ajuste Este factor es un factor prctico resultado de experiencias El mantenimiento ayuda a E de diseo

Panhandle A parcialmente turbulento Panhandle B totalmente turbulento

Modelo de Clculo
Condiciones Generales Tb = 520R = 60F = 16C 20C Pb = 14.73 psia =1 atm = 1.033 kgf/cm2 Caractersticas del Ducto Dimetro: Espesor: Rugosidad: Clculo Caudal: Presin Inicial P1: Presin Final P2: Longitud: Dimetro externo: Eficiencia (E): Caractersticas del Fluido Densidad:

t= PD 2Sadm

(ANSI acero)

Opciones Colebrook AGA Panhandle B Weymouth

Ejemplo:
CALCULO DE GASODUCTOS MONOFASICOS PRESSAO BASE: . 1.033 kgf/cm2 abs TEMPERATURA BASE: . 20 graus centigrados TEMPERATURA DO FLUXO: ... 25 graus centigrados DENSIDADE DO GAS: (densidad relativa) ESPESSURA DA TUBULACO:.... pol RUGOSIDADE DA TUBULACO:.. 0.045 mm VAZZAO: PRESSAO INICIAL: . PRESSAO FINAL: COMPRIMENTO: . DIAMETRO EXTERNO: .. EFICIENCIA DO GASODUTO: .. x 1000 m3/dia kgf / cm2 man. kgf / cm2 man. km pol 1

qu deseas calcular?
CALCULO DE GASODUCTOS MONOFASICOS PRESSAO BASE: . 1.033 kgf/cm2 abs TEMPERATURA BASE: . 20 graus centigrados TEMPERATURA DO FLUXO: ... 25 graus centigrados DENSIDADE DO GAS: 0.60 (densidad relativa) ESPESSURA DA TUBULACO:.... 0.50 pol RUGOSIDADE DA TUBULACO:.. 0.045 mm SELECO DO CALCULO 0. FIM DOS CALCULOS VAZZAO: 1. VAZAO PRESSAO INICIAL: . 2. PRESSAO INICIAL PRESSAO FINAL: 3. PRESSAO FINAL COMPRIMENTO: . 4. COMPRIMENTO DIAMETRO EXTERNO: .. 5. DIAMETRO EXTERNO EFICIENCIA DO GASODUTO: .. 6. MUDAR AS PROPIEDADES 7. IMPRIMIR A TELA x 1000 m3/dia kgf / cm2 man. kgf / cm2 man. km pol 1

con qu mtodo?
CALCULO DE GASODUCTOS MONOFASICOS PRESSAO BASE: . 1.033 kgf/cm2 abs TEMPERATURA BASE: . 20 graus centigrados TEMPERATURA DO FLUXO: ... 25 graus centigrados DENSIDADE DO GAS: 0.60 (densidad relativa) ESPESSURA DA TUBULACO:.... 0.50 pol RUGOSIDADE DA TUBULACO:.. 0.045 mm VAZZAO: PRESSAO INICIAL: . 150 FORMULA DE CALCULO PRESSAO FINAL: 30 1. COLEBROOK COMPRIMENTO: . 150 2. AGA DIAMETRO EXTERNO: .. 24 3. PANHANDLE B DO GASODUTO: .. 0.92 EFICIENCIA 4. WEYMOUTH x 1000 m3/dia kgf / cm2 man. kgf / cm2 man. km pol

Resultados
VOLMENES MXIMOS
ECUACIN 1000 M3/D COLEBROOK 23059 AGA 23059 PANHANDLE 28174 WEYMOUTH 23216 MMPCD 814 814 995 820

Tubo de 24"; esp: 0.5 " P1: 150 bar; P2: 30 bar Longitud: 150 km
30000 25000 20000

km3/d

15000 10000

Cada 100 Mw de generacin hidrulica equivalen a 20 MMpcd 30,000 familias consumen 1 MMpcd

5000 0

ol eb ro o

ey m ou th

A AG Pa

nh

Escenarios de Demanda
YEAR 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 GENERATION MMPCD 113 148 148 149 149 148 162 176 201 227 254 282 312 INDUSTRY MMPCD 47 63 75 83 87 91 95 99 104 108 113 118 124 COM/DOM MMPCD 2 3 5 9 12 15 19 22 26 28 29 30 30 TOTAL MMPCD 162 215 228 241 248 254 276 297 330 363 396 430 466 MAX FLOW MMPCD 227 303 324 348 365 344 378 413 461 508 551 596 644

MERCADO GENERATION INDUSTRY COMMERCIAL/RESID

FACTOR DE CARGA
5 AOS SGTES

ELCTRICO

70.00% 80.00% 26.00%

77.00% 90.00% 30.00%

Clculo de L

Bases de Diseo Tb = 20C P Pb = 1.033 kgf/cm2 Para un mismo caudal Temperatura media del gas= 25C Densidad= 0.60 Rugosidad de la tubera = 0.045 mm Espesor de la tubera = 0.625 " Caudal = 325 MMpcd Mximo = 9207 x 103 m3 P inicial = 150 kgf/cm2 L P final = 30 kgf/cm2 Optimizacin (Colebrook) Para 600 km P inicial necesaria= 126 kgf/cm2 = 1840 psi Resultados Para P inicial = 100 kgf/cm2 Longitud = 889 km slo cubre 355 kms

Tubo de 24"; esp: 0.625 " P1: 150 bar; P2: 30 bar Longitud: ? km

Diferentes distancias para distintos caudales

Tubo de 24"; esp: 0.625 " P1: 150 bar; P2: 30 bar Longitud: ? Km Distintos caudales
AOS OPER 3 5 10 15 CAUDAL MAX MMPCD 1000 M3/D 324 9,178 365 10,340 508 14,391 644 18,244 LONGITUD KMS 895 705 364 226

El Diseo a futuro
Horizonte de diseo Desarrollo por etapas Caudal mximo econmico Mnimos costos de operacin Mnima afectacin de medio ambiente y poblacin

150

300

600

Kms

Estaciones de Compresin
P
El diseo del gasoducto Camisea-Lima en 1988 segn Petroper consider tres estaciones de compresin para un flujo de diseo mximo de 450 MMPCD y una presin inicial de 100 bar
1 bar = 1.01972 kgf/cm2 = 0.98692 atm.

150

300

400

600
Kms

Estaciones intermedias de compresin

Estaciones de Compresin
P
p2
K = Cp / Cv

65

real --> politrpico isotrmico

MW

p1
15

V
Simplificado HP= 22 . (R) . N . (MMscfd) . F R : relacin de compresin por etapas N : nmero de etapas R recomendable menor a 4.5 Mayores valores llevan a temperaturas elevadas Para gases con gravedad = 0.8 - 1.0, usar 20

1.04 1.08 1.12 1.16 1.20 1.24 1.28 1.32 k

En clculos ms sofisticados es comn utilizar 1.25 - 1.26 como coeficiente politrpico

Segn fabricantes (14.7 PSI, 60F)

BHP / MMPCD PARA COMPRIMIR GAS NATURAL
PRESION DE DESCARGA (psig)
25 0 10 25 50 75 100 150 200 250 300 350 400 450 500 550 600 750 1000 1250 1500 1750 2000 3000 4000 35 50 65 40 75 85 55 25 100 145 105 70 40 20 150 170 130 95 70 45 30 200 190 145 120 80 55 45 20 350 235 185 150 120 95 75 55 40 25 20 500 250 215 175 140 120 100 75 60 45 35 25 20 15 600 260 230 185 155 130 115 85 70 55 45 35 30 20 15 15 750 275 235 205 170 145 130 105 80 65 55 50 40 30 30 25 20 800 280 240 210 175 150 135 110 85 70 60 55 45 35 35 30 25 15 900 290 250 220 180 155 140 115 95 80 70 60 50 40 40 35 30 20 1000 300 255 230 190 165 145 120 100 85 75 65 55 50 45 40 35 25 1250 315 275 240 205 180 160 135 115 105 80 75 70 60 55 50 45 35 20 1500 330 290 255 215 195 170 150 125 115 85 90 80 70 60 55 55 40 30 25 1750 345 300 265 230 205 180 160 140 125 100 100 90 80 75 70 65 55 40 30 20 2000 350 310 275 245 215 195 165 150 135 120 110 105 100 95 90 85 75 50 40 30 20 320 285 255 225 205 175 165 145 130 120 110 105 100 95 90 80 60 50 35 30 20 330 295 265 235 215 185 175 160 140 130 120 115 110 110 100 90 70 60 50 40 35 350 320 280 260 235 210 195 180 160 150 140 130 125 125 115 105 90 75 65 55 45 20 360 330 285 270 245 220 205 190 170 160 150 140 135 130 125 115 105 90 75 65 55 35 20 390 325 300 275 260 240 225 205 195 175 165 160 155 140 125 110 90 80 70 50 30 185 2500 3000 4000 4500 6000 65 100 130

PRESION DE SUCCIN (psig)

Knight Industries

Una etapa

Dos etapas

Tres etapas

370

Compresin: Resultados
CAUDAL MMPCD 324 365 508 644 METODO SIMPLIFICADO KNIGHT 34,427 34,020 38,783 38,325 53,977 53,340 68,428 67,620 PROGRAMA 28,822 32,471 45,192 57,292

Un factor de correccin : (TS/520) * (FGESP) Para: TS = 25C=537 FGESP = 0.965 (GPA) Factor= 0.9965

CAUDAL MMPCD 324 365 80000 70000 508 60000 644

HP

METODO COMPARACION RESULTADOS SIMPLIFICADO KNIGHTDEPROGRAMA 34427 34020 28822 38783 38325 32471 53977 53340 45192 68428 67620 57292

50000 40000 30000 20000 10000 0 300

SIMPLIFICADO KNIGHT PROGRAMA

400

500 MMPCD

600

700

ASME B 31.8 - 1995

Diseo Fabricacin Instalacin Inspeccin Pruebas Seguridad Operacin y Mantenimiento Materiales Tuberas de acero Gasoductos Redes de distribucin Ductos menores Estaciones de compresin Estaciones de Medicin y regulacin Lneas de Servicio

: API 5L : ASTM A 53 : ASTM A 106 Hierro dctil : ANSI A21.52 Tuberas plsticas y accesorios : ASTM D 2513-2517 Componentes y detalles de fabricacin : ANSI B16, API...

Criterios de diseo
Expansin y flexibilidad Esfuerzos combinados Soportes y anclajes Esfuerzos circunferenciales Factor de localizacin (F) Factor de tubera (E) Factor de temperatura (T)

t=

PD . 2 (SMYS) . FET

Material y Dimensionamiento: Dimetro (D) Espesor (t)

P presin de diseo SMYS specified minimun yield stenght

S ADM = (SMYS) . FET esfuerzo admisible

t=

PD

2 S ADM

Factores
CLASE DE LOCALIZACIN CLASE 1, DIV 1 CLASE 1, DIV 2 CLASE 2 CLASE 3 CLASE 4 F 0.80 0.72 0.60 0.50 0.40 N EDIFIC. PH <= 10 1.25 <= 10 1.1 10 < # < 46 > 46 > 46 ALTO TRFICO

CLASE DE TUBERA A 53 SIN COSTURA A 53 RESISTENCIA ELCTRICA A 53 SOLDADURA A TOPE - HORNO A 106 SIN COSTURA API 5L SIN COSTURA API 5L RESISTENCIA ELCTRICA API 5L ELECTRIC FLASH WELDED API 5L ARCO SUMERGIDO A 106 SOLDADURA A TOPE - HORNO
TEMPERATURA F 250 O MENOS 300 350 400 450

E 1.00 1.00 0.60 1.00 1.00 1.00 1.00 1.00 0.60

T 1.000 0.967 0.933 0.900 0.867

Usual: 0.72 Ejemplo: API 5LX 60 SMYS=60000 psi E=1

SADM = 43,200 psi

Si: D = 24" P = 150 bar = 2250 psi Luego:

t = 0.625 pulg
(no adic. x corrosin)

Adems...
Tuberas de plstico
General: menor o igual a 100 psi

Diseo termoplstico operacin mayor a -20F menor a "long term temp" y 140F

P= (2 S) t (0.32) (D-t)
P presin de diseo SMYS long term hydrostatic stenght

Diseo reforzado -20F <operacin<150F t>0.060"

Estaciones de Compresin Localizacin . Espacios abiertos Edificio albergue. ANSI/NFPA 220 Salidas >= 2 reas cercadas Facilidades elctricas. ANSI/NFPA 70 Equipos. Tratamiento de gas y otros Tuberas en la estacin

Vlvulas (transporte y distribucin) Espaciamiento Localizacin

LINEAS DE TRANSPORTE
LOCATION CLASE 1 CLASE 2 CLASE 3 CLASE 4 DISTANCIA MX 20 mi 15 mi 10 mi 5 mi

Evaluacin econmica
OPCION 1 Ducto 24" Presin mxima : 150 bar t= 0.688" Peso = 255 kg/m Total en 600 km : 153,000 ton Tubera : US$ 450 millones OPCION 2 Ducto 24" Presin mxima : 100 bar t= 0.438" Peso = 164 kg/m Total en 600 km : 98,400 ton Tubera : US$ 340 millones

Ducto de 600 kms Demanda indicada No std by compressor Mdulos de 15000 HP US$ 1500 / HP US$ 2000 / ton+10000/km-"

Costos de compresin
YEAR TOTAL MMPCD 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 162 215 228 241 248 254 276 297 330 363 396 430 466 MAX FLOW MMPCD 227 303 324 348 365 370 378 413 461 508 551 596 644 HP 150 BAR 23,804 31,802 34,026 36,524 38,315 38,850 39,723 43,320 48,441 53,336 57,902 62,595 67,622 HP COMER 30,000 45,000 45,000 45,000 45,000 45,000 45,000 45,000 60,000 60,000 60,000 67,000 67,000 NO NO NO 1 1 1 1 1 1 3 3 3 3 REQ INV mm$ 0 0 68 0 0 0 0 23 180 0 32 0 0 HP 100 BAR 15,869 21,202 22,684 24,350 25,543 25,900 26,482 28,880 32,294 35,557 38,602 41,730 45,081 HP COMER 15,000 30,000 30,000 30,000 30,000 30,000 30,000 30,000 45,000 45,000 45,000 45,000 45,000 1 2 2 2 2 5 5 5 5 5 5 5 5 REQ INV mm$ 23 68 0 0 0 135 0 0 113 0 0 0 0 0

Comparacin de NPV
YEAR TOTAL MMPCD 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 TASA 0% 10% 12% 162 215 228 241 248 254 276 297 330 363 396 430 466 NO NO NO 1 1 1 1 1 1 3 3 3 3 REQ INV mm$ 450 0 0 68 0 0 0 0 23 180 0 32 0 0 NPV 752 544 519 9 9 9 11 11 11 11 12 12 14 14 15 15 NPV 154 79 71 227 303 324 348 365 370 378 413 461 508 551 596 644 NPV 5,488 2,721 2,417 1,502 745 662 62 83 89 95 100 101 104 113 126 139 151 163 176 1 2 2 2 2 5 5 5 5 5 5 5 5 MANTTO mm$ MAX MMPCD INGRESOS REQ mm$ INV mm$ 363 68 0 0 0 135 0 0 113 0 0 0 0 0 NPV 678 510 487 11 13 13 13 13 19 19 19 20 20 20 20 20 NPV 221 112 100 MANTTO mm$

US$0.75 / mpc
MAX MMPCD 227 303 324 348 365 370 378 413 461 508 517 517 517 NPV 5,248 2,647 2,358 1,437 725 645 INGRESOS mm$ 62 83 89 95 100 101 104 113 126 139 142 142 142

GASODUCTO A 150 BAR NPV@10% ING-EG = $122 mm NPV@12% ING-EG = $72 mm

GASODUCTO A 100 BAR NPV@10% ING-EG = $103 NPV@12% ING-EG = $58