Análisis de Modelos de Optimización de Redes de Transmisión Del G

Universidad de La Salle
Ciencia Unisalle
Ingeniería Eléctrica Facultad de Ingeniería
2014
Análisis de modelos de optimización de redes de transmisión del

gas natural con énfasis en el sistema colombiano
Andrés Mauricio Macías Valencia
Universidad de La Salle, Bogotá
Camilo Andrés Moya Garzón

Universidad de La Salle, Bogotá
Follow this and additional works at: https://ciencia.lasalle.edu.co/ing_electrica
Part of the Electrical and Electronics Commons
Citación recomendada
Macías Valencia, A. M., & Moya Garzón, C. A. (2014). Análisis de modelos de optimización de redes de
transmisión del gas natural con énfasis en el sistema colombiano. Retrieved from
https://ciencia.lasalle.edu.co/ing_electrica/19
This Trabajo de grado - Pregrado is brought to you for free and open access by the Facultad de Ingeniería at
Ciencia Unisalle. It has been accepted for inclusion in Ingeniería Eléctrica by an authorized administrator of Ciencia
Unisalle. For more information, please contact [email protected].
ANÀLISIS DE MODELOS
DE OPTIMIZACIÓN DE
REDES DE
TRANSMISIÓN DEL GAS
NATURAL CON ÈNFASIS
EN EL SISTEMA
COLOMBIANO
Andrés Mauricio Macías Valencia.
Camilo Andrés Moya Garzón.
Este documento contiene modelos de optimización
en transmisión de gas natural según los distintos
escenarios, haciendo énfasis en la posible aplicación
de estos en el sistema colombiano. Se plantea como
guía para optimización de redes existentes y futuras.
UNIVERSIDAD DE LA SALLE
XX DE 03 DEL 2014
AGRADECIMIENTOS
Dedico este trabajo principalmente a Dios, por haberme dado la vida y permitirme el
haber llegado hasta este momento tan importante de mi formación profesional. A mi
familia (Rafael Moya, Hermelinda Garzón y Laura Moya), por ser el motor que me
impulso para llegar a la meta y demostrarme siempre su cariño y apoyo incondicional
sin importar nuestras diferencias de opiniones. También agradezco a los profesores
pues ellos me enseñaron y me dieron fuerza para lograr cada uno de los objetivos que
me proponía en el semestre. Un reconocimiento especial para el ingeniero Adrián
Correa, director de tesis, por haberse interesado en nuestro trabajo, por sus aportes y
paciencia a lo largo del último año. Y agradezco a mi amigo incondicional Andrés
Macías pues él ha sido un apoyo muy importante, compartimos gran parte de la
carrera y con el compartimos este proyecto y lo sacamos a flote, Gracias.
Finalmente, a mis compañeros, compañeras, amigos y amigas que fueron una ayuda
constante para sortear cada uno de los problemas que se nos presentaron durante la
carrera.
Att: Camilo Andrés Moya Garzón
En primer lugar, doy gracias a Dios, que me ha permitido tener éxito en todos mis
esfuerzos de la vida.
Agradezco todo mi cariño y mi amor a mis padres Luz Dary Valencia y Rafael Dúfran
Macías que hicieron todo en la vida para que yo pudiera lograr mis sueños, por
motivarme y darme la mano cuando sentía que el camino se terminaba, a ustedes por
siempre mi corazón y mi agradecimiento.
A mi compañero y amigo de Camilo Moya por toda la dedicación, esfuerzo y

continuidad no solo en la realización del proyecto de grado sino también en el
crecimiento de esta amistad, a él le deseo los mejores éxitos en su vida personal y
profesional.
A mis profesores que en este andar por la vida, influyeron con sus lecciones y
experiencias en formarme como una persona de bien y preparada para los retos que
pone la vida, a todos y cada uno de ellos les dedico cada una de estas páginas de mi
tesis.
Y por último pero no menos importante a mi director de Tesis el ingeniero Adrián

Correa por todo el apoyo y motivación brindada en el proceso de realización de este
proyecto de grado, gracias por ser consejero y guía.
Att: Andrés Mauricio Macías Valencia.
1
INDICE
INTRODUCCIÓN………………………………………………………………………………..4
PLANTEAMIENTO DEL PROBLEMA……………………………..………………………….4
OBJETIVO GENERAL……………………………………………………..……………………….5
OBJETIVOS ESPECIFICOS…………………………………..……………………………………5
JUSTIFICACION………………………………………………………………………………………5
CAPÍTULO 1…………………………………………………………………………………………..6
GENERALIDADES DEL MERCADO DEL GAS NATURAL
1.1. Historia del gas en Colombia …………………………………………….………...….….6
1.2. Generalidades de la transmisión de Gas Natural………………………..……….6
1.2.1 Gasoductos entre departamentos…………………………………………….7
1.3. Usos del Gas Natural…………………………………………………………….………..…..7
1.4. Estructura de Gasoductos en Sudamérica y Centro América….……….….8
1.5. Empresas transportadoras de gas en Colombia…………………….……….…..9
1.5.1 Transportadora de Gas Internacional – TGI……………….……….…..9

1.5.2 Promigas……………………………………………………………….…….….....10
1.5.3 ProGasur……………………………………………………………….…….…….11
1.5.4 Transgastol…………………………………………………………….….…..….13
1.5.5 Transoccidente ……………………………………………………….….……..14
1.6. Mapa de conexiones de SNT……………………………………………….…….…...15
1.7. Estructura del sector del gas……………………………………………….…….…...17
CAPÍTULO 2…………………………………………………………………………………….…..19
ACERCAMIENTO MATEMÁTICO
2.1. Definiciones……………………………………………………………………….……….…19
2.2. Sistema de Ecuaciones en Estado Estacionario…………………….…………..23
2.2.1. Formulación de ecuaciones de flujo de tubería………....…………23

2.2.2. Análisis cuando la tubería tiene un lazo auxiliar…….……………25
2.2.3. Formulación de las ecuaciones de balance de masa….……….…25
2
CAPÍTULO 3……………………………………………………………………………………….27
MODELOS DE OPTIMIZACIÓN DE TRANSMISIÓN DE GAS NATURAL

3.1. ESTUDIOS ASOCIADOS CON LA OPERACIÓN DE LA RED DE GAS
(MINIMIZACIÓN DE DIFERENTES COSTOS OPERATIVOS)………………………27
3.1.1 Minimización de los costos de combustible al final de la red en las

estaciones compresoras……………………………………………...…28
3.1.2. Optimización en el transporte de gas natural con base en el
consumo………………………………………………………………………32
3.2. ESTUDIOS ASOCIADOS CON LA EXPANSIÓN Y REFUERZO DE REDES

EXISTENTES (MINIMIZACIÓN DE COSTOS DE INVERSIÓN PARA
EXPANSIÓN A MEDIANO Y LARGO PLAZO DE LA RED)……………..…35
3.2.1 Minimizar el consumo de combustible en una red de tuberías…35
3.3. ESTUDIOS ASOCIADOS CON EL DISEÑO ÓPTIMO DE NUEVAS REDES

(MIN. COSTOS DE INVERSIÓN…………………………………………………..38
3.3.1 Localización óptima de las instalaciones de transmisión….38
3.4. ESTUDIOS INTEGRADOS GAS-ELECTRICIDAD………………………………41

3.5. OBSERVACIÓN SOBRE LOS MODELOS……………………………………….42
3.6. MODELOS DE PRUEBA PARA TRANSMISIÓN DE GAS NATURAL….44
3.6.1. Modelo de prueba 1………………………………………………………….44

3.6.2. Modelo de prueba 2………………………………………………………48
3.6.3. Modelo de prueba 3………………………………………………………50
INTERPRETACION Y CONCLUSIONES………………………………………………...52
REFERENCIAS……………………………………………………………………………………54
ANEXOS……..……………………………………………………………………………………56
3
INTRODUCCIÓN
En este proyecto de grado analiza los modelos de optimización de la transmisión del
gas natural buscando conocerlos, puesto que este tema en Colombia no se estudia.
Uno de los objetivos ocultos de este proyecto es dar a conocer estos modelos llegando
a ser utilizados en el sistema colombiano. Para esto se crearon 3 capítulos enfocados al
tema propuesto.
En el primer capítulo se muestra lo que hay en el sistema colombiano con respecto a la

red de transporte de gas natural. Damos a conocer los principales operadores de red
que maniobran el sistema hoy en día.
El segundo capítulo demuestra el origen de las restricciones y función objetivo a partir

de las ecuaciones en estado estacionario de las tuberías, compresores entre otros con
el fin de demostrar principalmente las restricciones.
El tercer capítulo muestra los modelos estudiados y analizados, donde muestra la

descripción de estos, la explicación de las restricciones y la aplicabilidad en el sistema
colombiano. Aparte de lo antes mencionado muestra los modelos de prueba donde se
muestran los datos necesarios para correr una simulación y aparte evidencia la
aplicabilidad en cualquier contexto.
Al final este documento analiza y demuestra que aplicando estos modelos en nuestro
sistema podríamos mejorar el sistema de Colombia. Sin embargo se espera dar un
primer paso dando a conocer este documento.
PLANTEAMIENTO DEL PROBLEMA

Colombia a pesar de ser un país rico en gas natural, la aplicación de modelos de
optimización en la cadena del gas natural y haciendo énfasis en la transmisión es casi
nula. Buscando con este análisis a los modelos de optimización de transporte en el gas
natural ya existentes, conlleve a mas investigaciones y la adaptación y creación de
nuevos modelos de optimización que apliquen completamente a la transmisión de gas
natural colombiana.
Este análisis se hace con la idea de entender los modelos de optimización en la

transmisión del gas natural y darlos a conocer ya que pocas empresas tienen
conocimiento sobre ellos; en algunos modelos de optimización de gas natural se
conocen cuáles son las variables y restricciones pero no en qué consiste cada una de
las variables, como se usan las restricciones, en qué consisten las restricciones, que
función tiene el modelo, etc. Con este análisis se pretende redactar y mostrar una
monografía con cada uno de los conocimientos mencionados, de esta forma incentivar
a la aplicación de los modelos de optimización de gas natural en el contexto
colombiano.
Se analizaran los modelos de optimización en transmisión de gas natural más

importantes y que se consideran factibles al uso en Colombia.
4
OBJETIVO GENERAL
Realizar un compendio de los modelos de optimización en la transmisión de las redes
de gas natural aplicados al contexto de los sistemas de transporte del gas natural en
Colombia.
OBJETIVOS ESPECÍFICOS
 Describir detalladamente las variables que componen cada uno de los modelos de
optimización.
 Puntualizar los modelos completos identificando su posible aplicación.
 Detallar las restricciones asociadas a los modelos que son importantes en cada uno
de estos.
 Realizar un análisis de los modelos aplicables al caso del sistema de transporte gas
natural colombiano.
JUSTIFICACIÓN
Se estudiarán y analizarán los modelos de optimización de trasmisión de gas natural
con la intención de poder identificar cada uno de sus variables, la determinación de
estas en el contexto colombiano y con esto dejar un camino abierto a las diferentes
aplicaciones y soluciones que se puedan crear.
5
CAPITULO 1
GENERALIDADES DEL MERCADO DEL GAS NATURAL
1.1. Historia del gas en Colombia
En Colombia el desarrollo de la industria del gas natural es reciente. Sin embargo en el

año 1950 se empezaron a realizar algunos estudios logrando usos aislados de este
combustible, pero fue a mediados de la década de 1970 cuando comenzó su verdadero
desarrollo gracias al descubrimiento de un pozo de gas en la Guajira donde empezó a
ser explotado en 1977.
En 1986 luego de un largo período de bajo crecimiento, inició el programa “Gas para el
cambio” que permitió ampliar el consumo de gas en las ciudades y con esto realizar
parte de la interconexión nacional que conocemos hoy y de este modo descubrir
nuevos pozos.
En 1993 el Gobierno Nacional definió que ECOPETROL liderara la interconexión

nacional, para lo cual dos años después comenzaron las conexiones entre los
principales yacimientos y centros de consumo, mediante la construcción de más de
2.000 km de gasoductos que pasaron por el departamento de la Guajira, el centro y
suroccidente del país y los Llanos orientales.
Entre 1997 y 1998 se otorgaron concesiones de áreas de distribución exclusiva de gas

para extender la cobertura del servicio a los Departamentos de Quindío, Caldas,
Risaralda, Valle y Tolima.
El Gobierno Nacional interesado en promover el desarrollo de este combustible en

todo el país y de masificar su uso, estableció en el 2003 las “Estrategias para la
dinamización y consolidación del gas natural en Colombia”, donde se formularon
algunas estrategias y recomendaciones para lograr este objetivo.
Un año después se hizo lo mismo para masificar el Gas Natural Vehicular y se ordenó
ofrecer condiciones económicas especiales (especialmente descuentos y bonos) para
beneficiar a quienes utilicen este combustible.
En el 2007 Ecopetrol, PDVSA (petrolera venezolana) y Chevron suscribieron un

contrato mediante el cual determinaron las condiciones para compra y venta de gas
natural entre Colombia y Venezuela durante los próximos 20 años [1].
1.2. Generalidades de la transmisión del Gas
Dentro de las principales generalidades de las redes de transmisión de gas natural

encontramos lo siguiente:
Hay tres tipos principales de tuberías a lo largo de la ruta de transporte: el sistema de

tuberías para recolección, el sistema de tuberías departamental y el sistema de
tuberías para distribución. El sistema de tuberías para recolección consiste en baja
6
presión, donde las tuberías son de gran diámetro puesto que transporta el gas natural
crudo desde el pozo hasta la planta de procesamiento. En casos particulares los pozos
tienen un alto contenido de azufre y de dióxido de carbono (gas ácido), para estos
casos debe ser instalado un tubo para la recolección de gas amargo especializado. El
gas contaminado (gas ácido) es corrosivo, por este motivo para el transporte desde el
pozo a la planta de endulzamiento debe hacerse con cuidado y con materiales
especiales para el transporte de gas [2].
Los otros dos tipos de tuberías se pueden clasificar como departamentales y

municipales o de distribución. Los gasoductos departamentales son similares a las vías
que comunican los departamentos puesto que estos llevan el gas de departamento en
departamento, otro ejemplo son las líneas de transmisión de energía eléctrica, estas
llevan la energía por largos trayectos a través de las fronteras entre departamentos.
Los gasoductos de distribución, por otro lado, son aquellos que transportan el gas
natural dentro de un municipio en particular y tal vez municipios aledaños. Este
proyecto de grado cubre sólo los fundamentos de gasoductos departamentales, sin
embargo, los detalles técnicos y operativos descritos son esencialmente los mismos
para gasoductos municipales.
1.2.1. Gasoductos entre departamentos
La red de gasoductos departamentales transporta gas natural procesado en las plantas

de procesamiento en las regiones productoras de las zonas con altos requerimientos
de gas natural, las zonas urbanas particularmente grandes. Más adelante en este
proyecto veremos cómo se distribuyen las redes de tuberías a través de todo el país.
[2]
El gas natural que se transporta a través de gasoductos departamentales viaja a altas

presiones en las tuberías, con presiones entre 200 y 1.500 libras por pulgada cuadrada
(psi). Esto reduce el volumen del gas natural transportado (hasta en 600 veces), así
como propulsión de gas natural a través de la tubería.
El modelo matemático del sistema de transmisión de gas es representado por un

conjunto de ecuaciones que describen los procesos térmicos e hidráulicos que se
originan en la tubería y en los elementos de la estación de compresores.
1.3. Usos del Gas Natural
El gas natural en el sector industrial y petroquímico puede ser utilizado como

combustible o materia prima.
Como combustible se emplea en varios tipos de equipos, por ejemplo: hornos,

secadores y calderas. En las industrias de cerámicas, cemento, metales, y otras donde
se requieren hornos, el aprovechamiento energético y el ahorro en el consumo son
notorios cuando se utiliza gas natural.
Como materia prima se utiliza en las industrias que requieren metano (principal
componente del gas natural) en sus procesos. Algunos de los subproductos del metano
son: Monóxido de carbono, Hidrógeno, Metanol, Ácido acético, Anhídrido acético,
entre otros [3].
7
1.4. Estructura de Gasoductos en Sudamérica
Como se muestra en la Figura 1, el sistema de conexión en Sudamérica es poco denso

en su estructura pues tiene largos tramos entre estaciones compresoras y las redes
existentes son muy diferentes para cada territorio o país.
Figura 1. Mapa de Sudamérica con la estructura de gasoductos [4].
Por ejemplo Chile en su largo territorio tiene pequeñas redes que están
interconectados con Argentina y este a su vez con Uruguay. Algunas de estas redes
siguen su interconexión con Brasil, Bolivia y Perú. Otro ejemplo es la no existencia de
redes de gas en Paraguay. Uno de los países que no tiene conexión con otro país es
Ecuador, es decir que esté país se abastece con sus propios fosos. Yendo al norte del
8
territorio sudamericano está Colombia interconectado con Venezuela, sin embargo,
Colombia tiene un pequeño porcentaje conectado con Venezuela de su total de redes.
Podemos ver en detalle podemos que las redes de transmisión tienen poca extensión
entre estaciones compresoras. Pero la explicación más acertada es por sus cordilleras y
montañas puesto que hacer fluir el gas con estos cambios de presión y de temperatura
se complica y es por esto que se ve mucho más denso en la imagen antes vista.
1.5. Empresas transportadoras de gas en Colombia
En Colombia hay varios campos de producción de gas natural, las cuencas de la Guajira
y de los Llanos Orientales son las de mayor producción; en el primer trimestre de 2013
contribuyeron con alrededor del 50% y del 33% de la producción nacional
respectivamente.
El gas natural se transporta desde las zonas de producción hasta las zonas de consumo
por medio del Sistema Nacional de Transporte (SNT). Según el decreto 1493 de 2003
del Ministerio de Minas y Energía, el SNT "es el conjunto de gasoductos localizados en
el territorio nacional, excluyendo conexiones y gasoductos dedicados, que vinculan los
centros de producción de gas del país con las Puertas de Ciudad, Sistemas de
Distribución, Usuarios No Regulados, Interconexiones Internacionales o Sistemas de
Almacenamiento". [5]
Las empresas transportadoras de gas natural en Colombia son:
 Transportadora de Gas Internacional – TGI.

 Promigas.
 Promotora de Gases del Sur – Progasur.
 Transportadora de Metano – Transmetano.
 Transportadora Colombiana de Gas – Transcogas.
 Sociedad Transportadora de Gas del Oriente – Transoriente.
 Transportadora Gasoducto del Tolima - Transgastol.
 Transoccidente.
 Coinobras.
A continuación encontraremos una descripción breve de algunas de las distribuidoras

antes mencionadas:
1.5.1 Transportadora de Gas Internacional – TGI.
En la actualidad es la transportadora de gas más grande en Colombia puesto que tiene

una capacidad de transporte promedio de 420 MPCD y cuenta con 3774 kilómetros de
gasoductos.
Su infraestructura se extiende desde el Departamento de la Guajira en el norte sobre

la Costa Caribe hasta los Departamentos de Tolima y Huila en el Centro y Sur del país.
También lo hace desde los Departamentos de Casanare y Meta en el oriente hasta el
Departamento del Valle en el Occidente, pasando por la sabana donde se encuentra
Bogotá en el centro del país. [6]
9
Actualmente desarrolla un proyecto de expansión que implica un aumento de su
capacidad transportadora en más del 60% desde el norte de Colombia (Guajira), por
USD $142 millones, y desde el oriente (Cusiana – Casanare) por USD $450 millones de
dólares. [7]
1.5.2 Promigas
La infraestructura de transporte del país está conformada por aproximadamente 7000

km de gasoductos, de los cuales 2896 corresponden a Promigas en conjunto con las
transportadoras de su portafolio. La capacidad de transporte total del sistema es 737
MPCD.
En la actualidad transportan en conjunto por su sistema de gasoductos alrededor de
50% del gas natural consumido en el país. [8]
Figura 2. Estructura de Promigas del 2014
Tabla 1. Participación de Promigas en algunos tramos de red de transmisión de gas natural.
10
Como se muestra en la Figura 2 y Tabla 1, Promigas está compuesto por varias
empresas que distribuyen el gas para gran parte de Colombia y Perú, sin embargo para
este caso veremos las tres subempresas mas importantes llamadas Transmetano,
Transoccidente y Transoriente o Promioriente puesto que estas son las que más gas
transportan para Colombia. Viendo los ítems que componen la Tabla 1 se encuentran
ítems técnicos y económicos. Dentro de los ítems técnicos encontramos la longitud,
presión máxima y capacidad máxima, entre los técnicos vemos que las redes de
transporte de Transoriente son más extensas que las demás y cuentan con gran
presión para su transporte, sin embargo, la capacidad máxima de transporte es de 42
millones de pies cúbicos por día lo que limita su producción para el consumo. En ese
caso la empresa que cuenta con mayor capacidad es Transmetano.
En la parte económica Promigas controla el 99,64% de la empresa Transmetano, el

69% de la empresa Transocciednete y el 73,3% de la empresa Transoriente o
Promioriente.
Al ver la estructura de Promigas (Figura 2) también se puede ver que controla varias
empresas por completo y que participa de varias, siendo una de las principales
transportadoras y distribuidoras de gas para Colombia.
Tabla 2. Sectores de consumo de gas natural.
En la Tabla 2 se observa el consumo de gas natural por sectores presentado por

Promigas, que podría usarse como ejemplo para notar el consumo general de
Colombia. Vemos que la parte eléctrica para generación de energía es bastante
grande, sin embargo, el consumo ha bajado, en el 2010 el consumo de gas fue mayor
al del 2011 y siguió disminuyendo en el 2012. Viendo esto de forma general el
consumo de gas ha disminuido considerablemente.
1.5.3 ProGasur
A finales de 1994 de obras importantes como el Gasoducto al Sur Neiva - Hobo (Huila),
en el año 2003 el Gasoducto Flandes (Tolima) – Girardot – Ricaurte (Cundinamarca), en
el año 2006 el gasoducto Guando Melgar (Tolima) – Fusagasuga (Cundinamarca) y
recientemente en el 2.009 del Gasoducto ramal Pradera/Jamundi (Valle del Cauca) –
Popayán (Cauca). [9]
11
Figura 3. Zonas donde ProGasur realizar su operación de transporte.
La Figura 3 muestra donde opera la empresa ProGasur. Esta empresa transporta el gas
a través de tuberías tipo departamental y a su vez tipo distribución en los
departamentos antes mencionados. Esta empresa usa un SCADA para monitorear y
hacer maniobras en sus redes de transporte. Los siguientes son dos ejemplos del
esquema SCADA presentado por la ProGasur.
Figura 4. SCADA de Barrancabermeja – Payoa – Bucaramanga
12
El SCADA muestra las estaciones compresoras y el sistema de tuberías que existen en
determinada parte del sistema de transmisión. Vemos que existen varias estaciones
compresoras que impulsan el gas a través de las tuberías. Las estaciones compresoras
tienen un bypass para que su mantenimiento sea más fácil y no tenga que parar la
transmisión de gas, sin embargo, podrían bajar la presión.
Figura 5. SCADA de Gibraltar – Bucaramanga
Esta imagen es otro ejemplo de un sistema de transmisión de gas menos complejo,

donde podemos ver que existen varias estaciones compresoras seguida lo que supone
que estos trayectos son para pasar algo que cambia de altitud y a su vez presión o que
son tramos largos. Acá vemos que pasa lo mismo con el bypass, usado para el
mantenimiento de la estación compresora sin interrupción del transporte.
1.5.4Transgastol
Transgastol es una empresa que opera desde 1994 transportando gas natural,
presentaron pruebas para prestación del servicio de transporte de gas durante 4 años,
donde ganaron varias de los estudios realizados por el Ministerio de Minas y Energía,
dando como resultado la construcción y operación de varios gasoductos, troncales del
subsistema norte Huila - Tolima, con el propósito de obtener la concesión de la zona.
[10] A continuación se ven esquemas de los estudios realizados:
13
Figura 6. Zona donde opera Transgastol transportando gas natural.
La imagen muestra que Transgastol es una pequeña empresa que tiene una pequeña
parte, donde comprende algunos tramos de la zona. Se observa también que tiene
algunas estaciones compresoras en puntos estratégicos que impulsan el gas hasta otro
punto estratégico. Las líneas punteadas son las tuberías que transportan el gas y la
línea gruesa continua es la vía que comunica cada una de las ciudades mencionadas en
la imagen.
1.5.5 Transoccidente
Transoccidente es una empresa vinculada de Promigas, nace como resultado de un

proceso de secesión de la compañía Gases de Occidente S.A. E.S.P. el 20 de Marzo de
1998. Los motivos que originaron dicho proceso obedecen a la evolución y dinámica al
interior del mercado de los servicios públicos. El objeto principal es la realización de la
actividad de transporte de gas combustible mediante la construcción, operación y
mantenimiento de sistemas y subsistemas de transporte.
Para efectos regulatorios y tarifarios, y amparada por las disposiciones de la CREG en

las Resoluciones 202 de 1997, 009 y 077 de 1998, el sistema que representa la
empresa se distingue con el nombre de Subsistema de Transporte del Valle 01, este
subsistema está comprendido entre la puerta de ciudad de Cali y la estación
reguladora y de medición ubicada en la Avenida 3N con calle 70N de la nomenclatura
de Cali, conocido como sector ACOPI, con sus conexiones y ramales asociados. El
subsistema se conecta del sistema de Transporte de Occidente (Mariquita-Cali) en el
CITY GATE de Cali, y presta el servicio de transporte de gas natural a las empresas de
distribución domiciliaria como Gases de Occidente S.A. E.S.P. y a Gases del Norte del
Valle S.A. E.S.P. [11]
14
1.6 Mapa de conexiones de SNT
Estos mapas muestran la distribución de las estaciones compresiones. [6]
Figura 7. Mapa de los oleoductos existentes, indicando las principales estaciones compresoras
15
Figura 8. Mapa de los oleoductos existentes, indicando los principales puntos de consumo.
En las anteriores imágenes se encuentran dos mapas sobre los oleoductos existentes
en Colombia. La primera imagen muestra las principales estaciones de compresión que
impulsan el gas a través de las tuberías. Dentro de las principales estaciones
16
encontramos las de Bucaramanga, Bogotá, Barrancabermeja, Miraflores entre otros.
Estas estaciones compresoras deben estar en puntos estratégicos y en puntos
necesarios para el correcto funcionamiento del sistema. En la Figura 8, vemos los
principales puntos de consumo del gas natural, entre los que encontramos los
departamentos más importantes económicamente para Colombia y las ciudades que
tienen empresas generadoras de energía eléctrica a partir del gas natural.
1.7. Estructura del sector del gas
Con el fin de que el sector funcione y se desarrolle de la mejor forma posible, se

estableció un esquema que involucra a las entidades que producen la energía, las que
la trasportan, las que la venden, las que coordinan a todas las anteriores, las que
establecen las políticas generales, las que hacen las normas para entregar productos
de buena calidad a un precio razonable y las que vigilan que todos cumplan las normas
existentes. [12]
Figura 9. Estructura de sector del gas según la CREG
Al ver la Figura 9 podemos deducir que; la parte política esta para diseñar las reglas
que se deben cumplir para la transmisión y distribución del gas natural siendo el
encargado el ministerio de minas y energía. La parte regulatoria es la encargada de
aplicar los reglamentos correspondientes según la actividad, en este caso el que le
compete al gas natural. Estos reglamentos son aplicados tanto a usuarios como a las
empresas, sin embargo, para el fin de esta tesis, los reglamentos que aplican son los
que tiene que ver con el transporte de gas a nivel de transmisión y no de distribución.
La entidad encargada es la Comisión de Regulación de Energía y Gas. El mercado
compete la parte te abastecer a los usuarios, estos se dividen en dos, los regulados y
los no regulados. Para la parte de transmisión, las empresas manejan directamente el
abastecimiento de las generadoras de energía eléctrica como se ve en la nueva
estructura de la industria (Figura 10). Como último esta la supervisión y el control de
las empresas que distribuyen el gas natural, esta parte está encomendada a la
superintendencia de servicios públicos domiciliarios que determinan quien cumple y
17
quien no cumple las reglas implementadas, y del mismo modo dan las sanciones
pertinentes para que el servicio sea lo más óptimo posible.
Figura 10. Nueva estructura de la industria
18
CAPÍTULO 2
MODELO MATEMÁTICO DE COMPONENTES

Con la problemática existente con la producción de gasolina y ACPM, la solución a
partir del gas natural para el abastecimiento de energía o combustible que supla la
necesidad de los usuarios, es la optimización del transporte en las redes de
transmisión que podría lograr inversiones más bajas pero con el mismo nivel de
seguridad como la disposición de la infraestructura.
El sistema de transporte de gas natural consiste en una compleja red de tuberías,

destinados al transporte de manera rápida y eficiente de gas natural desde su origen, a
las zonas de alta demanda de gas natural. El transporte de gas natural está
estrechamente relacionado con su almacenamiento en instalaciones generalmente
subterráneas.
Con este capítulo buscamos afianzar los conceptos específicos presentados en los
modelos de optimización de la transmisión de gas, esto con el objetivo de poder
facilitar la comprensión y entendimiento de los temas tratados allí.
2.1. DEFINICIONES
Para comprender lo modelos matemáticos y algunas de las partes que comprende este
proyecto de grado, se deben entender las siguientes definiciones:
Sistema en estado estable: Esta definición debe tenerse en cuenta, puesto que es así
como se desarrollaron los modelos de optimización, de esta forma se llega a las
diferentes soluciones para los diferentes problemas de optimización que tiene la red
de transmisión de gas. Una definición básica sobre el sistema en estado estable es que
los valores característicos del flujo de gas en el sistema no dependen del tiempo. Por lo
tanto, el desarrollo de estos modelos no tiene trascendencia a lo largo de su vida útil
para el desarrollo de los modelos de optimización. En este caso, el sistema de
restricciones, particularmente aquellas que describen el flujo del gas a través de los
ductos pueden ser descritas usando ecuaciones algebraicas no lineales. Existe otro tipo
de sistema cuando en el análisis se considera el tiempo, este tipo de sistema es
considerado en estado transitorio. En un sistema como este las variables del sistema,
están en función del tiempo, que afectan directamente las presiones y los flujos a
través de los elementos que comprender el sistema.
Para modelarlo se tiene que hacer uso de ecuaciones diferenciales parciales en función
del tiempo que describen la relación que existe entre estas variables.
19
Balance de masa: El balance de masa se realiza a partir de la ley de la conservación de
la materia basándose en la siguiente ecuación. [13]
𝑀𝑠𝑖𝑠𝑡𝑒𝑚𝑎 = 𝑀𝑒𝑛𝑡𝑟𝑎 − 𝑀𝑠𝑎𝑙𝑒 (1)
La anterior ecuación da a entender que un sistema está en equilibrio si solo si la masa

que entra en el sistema es igual a la masa que sale. Sin embargo podría pasar que la
masa que entra no es igual a la que sale por razones diferentes. Entre algunas, fugas
por falta de mantenimiento, materiales defectuosos, entre otros.
La conservación de la masa dice que la tasa de flujo de masa a través de un tubo es

una constante. Podemos determinar el valor de la tasa de flujo de masa de las
condiciones de flujo.
Si el fluido pasa inicialmente a través de un área A con velocidad V, podemos definir un

volumen de la masa para ser barrido en una cierta cantidad de tiempo t. El volumen v
es:
𝑣 = 𝐴 ∗ 𝑉 ∗ 𝑡 (2)
𝑙𝑜𝑛𝑔𝑖𝑡𝑢𝑑
Una verificación de unidades da á𝑟𝑒𝑎 𝑥 𝑥 𝑡𝑖𝑒𝑚𝑝𝑜 = á𝑟𝑒𝑎 𝑥 𝑙𝑜𝑛𝑔𝑖𝑡𝑢𝑑 =
𝑡𝑖𝑒𝑚𝑝𝑜
𝑣𝑜𝑙𝑢𝑚𝑒𝑛. La masa M que contiene este volumen es simplemente la densidad r veces
el volumen.
𝑀 = 𝑟∗𝐴∗𝑉∗𝑡 (3)
Para determinar el caudal másico 𝑀𝑑𝑜𝑡 , se divide la masa por el tiempo. La definición
resultante de la relación de flujo de masa se muestra en la siguiente ecuación.
𝑀𝑑𝑜𝑡 = 𝑟 ∗ 𝐴 ∗ 𝑉 (4)
Temperatura ambiente constante: Se considera un proceso isotérmico, es decir, que

la temperatura en el sistema se mantiene constante en cada ducto. Esto se hace con la
finalidad de hacer más tratable la restricción de flujo en los ductos.
Número fijo de compresores en cada estación compresora: El número de

compresores operando en cada estación se considera como un parámetro conocido.
Todas las unidades compresoras en una estación compresora son centrifugas, idénticas
e instaladas en paralelo. El hecho de considerar idénticas y conectadas en paralelo, nos
permite suponer que el flujo que entra a la estación compresora se reparte
equitativamente entre el número de compresores que tenga la estación lo cual facilita
el análisis. Este tipo de configuración con compresores centrífugos es también de los
más utilizados en el medio industrial.
Procesos Isentrópicos: Muchos sistemas de ingenierías o aparatos como bombas,

turbinas, toberas y difusores son esencialmente adiabáticos en su operación, y ellos se
desarrollan mejor cuando irreversibilidades, como fricción asociada con el proceso,
son minimizadas. Por ello, un proceso isentrópico puede servir como un modelo para
un proceso real, y también nos ayuda a definir las eficiencias de los procesos
20
comparando el desarrollo real de esos aparatos con el desarrollo bajo condiciones
idealizadas. [14]
Estación Compresora de Gas Natural: En una estación compresora el gas natural es

altamente presurizado a medida que viaja a través de una tubería de una estación a
otra. Para asegurarse que el gas natural fluya a través de cualquier tubería se mantiene
a presión, este requiere la compresión del gas natural periódicamente a lo largo de la
tubería. Por lo general estas estaciones están situadas en intervalos de 70 a 160
kilómetros a lo largo de la tubería. El gas natural entra en la estación de compresión,
donde se comprime mediante turbinas o motores.
Los compresores de turbinas obtienen su energía mediante el uso de una pequeña

proporción del gas natural que comprimen. La turbina en sí sirve para operar un
compresor centrífugo, que contiene un tipo de ventilador que comprime y bombea el
gas natural a través de la tubería. Algunas estaciones de compresores hacen funcionar
mediante el uso de un motor eléctrico para girar el mismo tipo de compresor
centrífugo. Este tipo de compresión no requiere el uso del gas natural de la tubería, sin
embargo, requiere una fuente fiable de electricidad cerca. Motores de gas natural
alternativos también se utilizan para algunas estaciones de compresión. Estos motores
se asemejan a un gran motor de automóvil, y son alimentados por gas natural de la
tubería. La combustión de los pistones de gas en el exterior del motor, sirve para
comprimir el gas natural. [2]
Además de la compresión de gas natural, las estaciones compresoras por lo general

contienen algún tipo de separador de líquido, al igual que los que se utilizan para
deshidratar el gas natural durante su procesamiento. Por lo general, estos separadores
se componen de depuradores y filtros que capturan los líquidos u otras partículas no
deseadas del gas natural en la tubería. Aunque el gas natural en tuberías de gas se
considera "seco", que no es raro que una cierta cantidad de agua y de hidrocarburos
para condensar fuera de la corriente de gas mientras está en tránsito. Los separadores
de líquidos en las estaciones de compresión aseguran que el gas natural en la tubería
es tan puro como sea posible, y por lo general filtrar el gas antes de la compresión.
Gasoducto: Son todas las partes de las instalaciones físicas a través de las cuales se
mueve el gas en su transporte, incluyendo tuberías válvulas, accesorios, bridas
(incluyendo el empernado y las empaquetaduras), reguladores, recipiente a presión,
amortiguadores de pulsación, válvulas de desfogue, y otros accesorios instalados en la
tubería, unidades de compresión, estaciones de medición, estaciones de regulación, y
conjuntos fabricados.
Se incluyen en esta definición las líneas de transporte y recolección de gas, incluyendo

sus complementos o accesorios que se halla instalados costa fuera para el transporte
de gas desde las instalaciones de producción a localidades en tierra y equipo s de
21
almacenamientos de gas del tipo de tubería cerrada, que se fabrican o se forjan de
tubería o se fabrican con tubería y accesorios.
La siguiente ilustración muestra las partes que componen una la red de transmisión o
gasoducto:
Figura 11. Estructura del sistema de transmisión de gas: EC – Estación de Compresión, CG - Compresor
de Gas, TG – Turbina de Gas, EA - Enfriador de Aire (después de más frío), T – Tubería, NE – Nodo de
entrega. [15]
Compresor: Es una máquina que fundamentalmente se utiliza para elevar la presión de

un gas pasando de presión baja a otra más alta. Existen dos tipos principales de
compresores que se explicaran a continuación.
Compresor Centrífugo: El compresor centrífugo es una turbomáquina que consiste en

un rotor que gira dentro de una carcasa provista de aberturas para el ingreso y egreso
del fluido. El rotor es el elemento que convierte la energía mecánica del eje en
cantidad de movimiento y por tanto energía cinética del fluido. En la carcasa se
encuentra incorporado el elemento que convierte la EC en energía potencial de
presión (el difusor) completando así la escala de conversión de energía. [16]
Compresor Axial: Los compresores axiales tienen ciertas ventajas y desventajas con
respecto a los compresores centrífugos. Entre las ventajas se pueden citar menor área
frontal (importante para usos aeronáuticos) y mayores relaciones de compresión y
eficiencias, aunque estas últimas ventajas no son tan grandes si se consideran
compresores centrífugos de varias etapas y de diseño moderno.
Las principales desventajas del compresor axial son su costo y su relativamente menor
robustez, dada la fragilidad de los álabes (comparando con el rotor centrífugo de una
sola pieza).
El compresor axial consiste en un rotor de forma cilíndrica que gira dentro de una
carcasa o estator. El fluido de trabajo circula por el espacio anular entre el rotor y el
estator, pasando por hileras de álabes fijos y móviles. [17]
Turbina de Gas: La turbina de gas es un equipo capaz de transformar la energía

química contenida en un combustible en energía mecánica, ya sea para su
22
aprovechamiento energético o como fuerza de impulso de aviones, automóviles o
barcos. En este caso la turbina esta para generar la energía que necesita el compresor
para impulsar el gas a su destino. Las tuberías pueden clasificarse según el origen de su
desarrollo, por el diseño de su cámara de combustión y por su número de ejes. Como
se ha indicado antes, una parte importante de esta potencia es absorbida
directamente por el compresor. Los gases, que entran a la turbina a una temperatura
de 1200-1400ºC y una presión de 10 a 30 bar., salen a unos 450-600ºC. Esa alta
temperatura hace que la energía que contienen pueda ser aprovechada bien para
mejorar el rendimiento de la turbina.
Tubería: Las tuberías son hechas con materiales especiales según su uso, son usados
para el transporte de fluidos tales como petróleo, agua, gas natural entre otros. Los
sistema de transporte son usualmente interconectados para conformar una red, que
transportan gas de un sistema de recolección, desde la salida de una planta de
procesamiento, o un campo de almacenamiento, hacia un sistema de distribución de
alta o baja presión, un cliente que compra un gran volumen, u otro campo de
almacenamiento.
2.2. Sistema de Ecuaciones en estado estacionario
Para el análisis de los modelos que serán mostrados en el capítulo 3 se usarán las
ecuaciones que se explicaran a continuación, es bueno saber que el modelo
matemático se desarrolló a partir del conocimiento de las características de
rendimiento de los compresores, las ecuaciones que gobiernan el flujo de gas a través
de tuberías, y los principios de conservación de la masa. La formulación matemática y
los tipos de ecuaciones incorporados en las ecuaciones de simulación que rigen
dependen de las configuraciones de la red, la naturaleza del gas y la corrosión interna.
[18]
2.2.1. Formulación de ecuaciones de flujo por tubería
La ecuación de flujo de la tubería es una de las ecuaciones que rigen la mayoría de las
restricciones de los modelos contemplados en esta tesis. Se deriva basado en el
principio del análisis de flujo de gas en las tuberías. El flujo de gas a través de las
tuberías puede verse afectada por varios factores tales como las propiedades del gas,
factor de fricción y la geometría de los tubos. La relación entre la presión de aguas
arriba, la presión aguas abajo y el flujo del gas en las tuberías puede ser descrita por
diversas ecuaciones. Para este análisis, la ecuación de flujo general se utiliza debido a
su frecuente aplicación en la industria de gas para describir la relación entre la
diferencia de presión y flujo de gas en las tuberías. [19]
La ecuación de flujo para la fase de la tubería se obtiene a partir de la siguiente figura

con la relación de presión aguas arriba 𝑃𝑖 (Punto inicial) y la presión aguas abajo 𝑃𝑗
23
(Punto final), gracias a estos cambios de presiones se produce el flujo a través del tubo
𝑄𝑖𝑗 , este se puede expresar como:
Figura 12. Ilustración unifilar de una tubería
𝑃𝑖2 − 𝑃𝑗2 = 𝐾𝑖𝑗 ∗ 𝑄𝑖𝑗

2
(5)
Donde 𝐾𝑖𝑗 toma diferentes formas dependiendo de las dimensiones de los parámetros
utilizados en las ecuaciones de flujo. Por ejemplo, Presión (P) [kPa], Temperatura (T)
[K], Longitud (L) [km], Caudal (Q) [m3/h] y Diámetro (D) [mm], la expresión para 𝐾𝑖𝑗
toma la forma:
𝑓𝐺𝑍𝑇 𝑃𝑛 2
𝐾𝑖𝑗 = 4.3599 𝑥 108 (𝑇 ) 𝐿 (6)
𝐷5 𝑛
Donde,
𝑓= Factor de fricción.
𝐺= Gravedad del gas.
𝑍= Compresibilidad del gas.
𝑇= Temperatura del gas.
𝐷= Diámetro.
𝑃𝑛 = Condiciones de presión estándar.
𝑇𝑛 = Condiciones de temperatura estándar.
𝐿= Longitud.
La ecuación 6 puede ser representada en forma de función, donde esta consta sólo de
parámetros que son desconocidos en esa ecuación. Si las presiones nodales y el caudal
no se conocen, la representación funcional toma la siguiente forma:
𝐹(𝑃𝑖 , 𝑃𝑗 , 𝑄𝑖𝑗 ) = 0 (7)
El enfoque del modelado de flujo monofásico puede no ser adecuadas para predecir la
capacidad de transporte de los gasoductos necesarios para mover los fluidos
mezclados.
Un análisis de flujo de dos fases o análisis de flujo de una sola fase con el factor de
fricción modificado puede ser necesario para predecir adecuadamente las capacidades
24
de transporte de este tipo de sistema. Por otra parte, la corrosión en las industrias de
petróleo y gas es uno de los graves problemas que afectan el desempeño de estas. A
medida que la edad de la tubería aumenta la rugosidad de la tubería tenderá a
aumentar debido a la acumulación de varios elementos alrededor de la superficie
interna de la tubería. Por lo tanto, la ecuación 7 debería ser modificada para tener en
cuenta el efecto de flujo multifásico y corrosión en las tuberías.
2.2.2. Análisis cuando la tubería tiene un lazo auxiliar
Con el fin de reducir la caída de presión en una cierta sección de la tubería debido a la
limitación de la presión o para aumentar la tasa de flujo en las secciones de cuello de
botella, la red de tuberías con un lazo puede ser construido. El sistema de tuberías con
un lazo se compone de dos o más tuberías conectadas de tal manera que el flujo de
gas se divide entre los tubos de bifurcación y, finalmente, se combinan aguas abajo en
un único tubo como se muestra en la siguiente figura.
Figura 13. Ilustración unifilar de una tubería con lazo auxiliar
De la anterior figura se deduce que la tubería se divide, esto se hace principalmente

para generar un punto de mantenimientos, puesto que este puede ser una tubería con
punto crítico donde el flujo o la presiona afecte más este trayecto. También se usa
cuando hay una estación compresora y se genera un lazo auxiliar para no parar de
transmitir y hacer mantenimiento a la máquina.
Estas ecuaciones adicionales se obtienen a partir del lazo. La condición del lazo
establece que para cada uno dentro del sistema de red cerrado la caída de presión es
cero. En la base de la condición de lazo, la caída de presión en el conducto de
derivación 1-2-4 debe ser igual a la caída de presión en el ramal de la tubería 1-3-4.
Esto es debido al hecho de que las dos ramas de tubería tienen un punto de partida
común (nodo 1) y el punto final común (nodo 4). La condición del lazo se puede
expresar con la siguiente ecuación:
𝑄2 𝐿 0.5 𝐷 2.5
= ( 𝐿2 ) 1
(𝐷 ) (8)
𝑄3 1 2
Donde, 𝑄2 , 𝐿1 y 𝐷1 son las partes de la ecuación donde se muestra el flujo a través del
tubo, la longitud y el diámetro para la tubería 1-2-4, respectivamente. Del mismo
modo, 𝑄3 , 𝐿2 y 𝐷2 son las partes de la ecuación donde se muestra el flujo a través del
tubo, la longitud y el diámetro para la tubería 1-3-4, respectivamente.
25
2.2.3. Formulación de las ecuaciones de balance de masa
Además de las ecuaciones de flujo de tuberías y compresores discutidos

anteriormente, el balance de masa proporciona las restantes ecuaciones básicas para
tener una formulación matemática completa para el análisis de los modelos. Las
ecuaciones de balance de masa se obtienen en base al principio de conservación de la
masa. En cualquier nodo, la ecuación de balance de masa generalizada de tuberías
entrantes (𝑡), tuberías salientes (𝑢) y tuberías de carga (𝑤) pueden ser resumidas
como:
𝑖=𝑡 𝑗=𝑢 𝑙=𝑤
∑ 𝑄𝑖 − ∑ 𝑞𝑗 − ∑ 𝐷𝑙 = 0
𝑖=1 𝑗=1 𝑙=1
(9)
Donde 𝑄1, 𝑄2 ,…, 𝑄𝑡 son flujo a través de las tuberías de entrada al nodo 𝑐, 𝑞1 , 𝑞2 ,…, 𝑞𝑢
son el flujo a través de tuberías de salida del nodo 𝑐 y 𝐷1 , 𝐷2 ,..., 𝐷𝑤 son la carga de la
salida del nodo 𝑐. Si todos los caudales a través de las tuberías de entrada y de salida
son desconocidos, entonces, la representación funcional de la ecuación 10 toma la
forma:
𝑓(𝑄1 , 𝑄2 , … , 𝑄𝑡 , 𝑞1 , 𝑞2 , … , 𝑞𝑢 ) = 0 (10)
26
CAPÍTULO 3
MODELOS DE OPTIMIZACIÓN DE TRANSMISIÓN DE GAS NATURAL

En la industria energética como en cualquier otro proceso es de vital importancia
lograr sacar el máximo provecho de los recursos, con el fin de ver reflejado esto en
ganancias económicas siendo más eficientes en algún proceso. La transmisión de gas
natural es común e importante en el proceso de industrial energético, ya sea el
escenario en que se transporte gas para centrales térmicas, centros de distribución
urbanos o inclusive para el uso como combustible en generadores con sistemas Bi-
Fuel.
Cuando se trabaja con redes de gas natural, cada topología se puede clasificar en uno
de los dos tipos siguientes:
 Radial: Sistema con más tiempo de desarrollo y como resultado de esto, es el
más documentado.
En una configuración radial es posible obtener el conjunto óptimo de los
caudales másicos de antemano, por lo que el problema se enfoca en encontrar
el conjunto óptimo de presiones en cada nodo. La mayoría de diseños de ellos
están basados en técnicas de programación dinámica (DP) por sus siglas en
inglés y han tenido un relativo éxito en su aplicación.
 Enmallado: Los sistemas enmallados el conjunto de variables, las tasas de flujo
de masa en los ductos y presiones en los mismos, deben ser tratadas de
manera explícita, así que la DP ya no puede proporcionar soluciones eficientes
como su tiempo de cálculo explota exponencialmente. Por lo tanto, el objetivo
de este trabajo es proporcionar soluciones de buena calidad para sistemas
enmallados.
Este capítulo se enfoca en identificar y caracterizar los diferentes modelos de

optimización encontrados en términos de la función objetivo, restricciones
involucradas y los tipos de variables en cada modelo de optimización de la
transmisión del gas.
3.1. ESTUDIOS ASOCIADOS CON LA OPERACIÓN DE LA RED DE GAS

(MINIMIZACIÓN DE DIFERENTES COSTOS OPERATIVOS).
Generalmente, las estaciones de compresores consumen en combustible

alrededor de 3 a 5% del total del gas que fluye a través de la red. Este valor llega
27
a ser significativo, sobre los miles de millones de pies cúbicos de gas que son
transportados cada día.
El gas natural es transportado por la presión a través de un sistema de tuberías,

esta transmisión está sujeta a pérdidas de energía causadas por:
 Fricción existente entre el gas y la pared interior de la tubería.

 Transferencia de calor entre el gas y el medio ambiente.
De ahí la importancia de encontrar una mejor manera de ubicar estas estaciones

de compresión a través de un sistema de tuberías.
3.1.1 Minimización de los costos de combustible al final de la red en las

estaciones compresoras.[20]
Se considera un problema de minimización de costos del combustible incurridos

en el consumo por estaciones compresoras a través de las redes de transmisión
del gas.
Esta red es representada en un esquema de distribución, donde los ductos
corresponden a tuberías, los trapecios a estaciones compresoras, y nodos
corresponden a sus puntos de interconexión física. Las variables de decisión son
el flujo de masa a través de cada ducto, y el nivel de presión de gas en cada
nodo. En cada estación compresora, hay un costo que depende de la entrada de
presión, la salida (descarga) de presión y del flujo de masa a través del
compresor.
La función que representa este costo es típicamente no convexa y no lineal.
Descripción del modelo:
La función objetivo del problema es la suma de los costos de combustible al final

de las estaciones compresoras en la red.
Este problema implica las siguientes restricciones:
i) Balance de flujo de masa en cada nodo.

ii) Relación de flujo de gas a través de cada tubo.
iii) Presión límite implicada a cada nodo.
iv) Límites de operación para cada estación compresora.
Las primeras dos restricciones mencionadas anteriormente también son

llamadas ecuaciones de flujo de red de estado estacionario. Se hace énfasis que
(i) es lineal y (ii) es no lineal.
28
Para flujos de presiones medianas y altas, cuando toman en cuenta el hecho que
un cambio de la dirección de flujo de corriente del gas puede tomar lugar en la
red, la ecuación de flujo de tubo toma la forma:
𝒑𝟐𝒊 − 𝒑𝟐𝒋 = 𝒄𝒊𝒋 ∗ 𝒖|𝒖|.𝜶 (11)
Donde pi y pj son presiones de los nodos del ducto (i,j), u es la relación de flujo
de masa a través del ducto, α es una constante (α ≈ 𝟏), y la resistencia de la
tubería 𝒄𝒊𝒋 es un cantidad positiva dependiendo de los atributos físicos del ducto
(i,j).
Las ecuaciones de flujo pueden ser estables en una muy concisa forma usando la
matriz de incidencia. Se considera una red con n nodos, L ductos, y m estaciones
compresoras. Cada ducto tiene asignada una dirección, 1 si tiene una dirección o
cero si no hay dirección del flujo de gas a través del ducto. Al dejar Al ser la
matriz nx1 cuyos elementos están dados por:
𝟏, 𝒔𝒊 𝒆𝒍 𝒕𝒖𝒃𝒐 𝒋 𝒔𝒂𝒍𝒆 𝒅𝒆𝒍 𝒏𝒐𝒅𝒐 𝒊;

𝒂𝒍𝒊,𝒋 = {−𝟏, 𝒔𝒊 𝒆𝒍 𝒕𝒖𝒃𝒐 𝒋 𝒆𝒏𝒕𝒓𝒂 𝒆𝒏 𝒆𝒍 𝒏𝒐𝒅𝒐 𝒊 ;
𝟎, 𝒐𝒕𝒓𝒐 𝒄𝒂𝒔𝒐. ;
Al es llamada la matriz incidencia del nodo-tubo. De forma similar, se deja Am

ser la matriz nxm cuyos elementos son dados por:
𝟏, 𝒔𝒊 𝒆𝒍 𝒏𝒐𝒅𝒐 𝒊, 𝒆𝒔 𝒆𝒍 𝒏𝒐𝒅𝒐 𝒅𝒆 𝒅𝒆𝒔𝒄𝒂𝒓𝒈𝒂 𝒅𝒆 𝒍𝒂 𝒆𝒔𝒕𝒂𝒄𝒊𝒐𝒏 𝒌;

𝒂𝒍𝒊,𝒋 = {−𝟏, 𝒔𝒊 𝒆𝒍 𝒏𝒐𝒅𝒐 𝒊 𝒆𝒔 𝒆𝒍 𝒏𝒐𝒅𝒐 𝒅𝒆 𝒔𝒖𝒄𝒄𝒊𝒐𝒏 𝒅𝒆 𝒍𝒂 𝒆𝒔𝒕𝒂𝒄𝒊𝒐𝒏 𝒌;
𝟎, 𝒐𝒕𝒓𝒐 𝒄𝒂𝒔𝒐.
Am es llamada la matriz de incidencia nodo-estación. La matriz formada

añadiendo Am a el lado derecho de Al esta denotado por A, i.e., A=(Al Am), la
cual es una matriz n x (l + m).
Se define que u=(u1……..ul)T y O=(O1……Om)T son la tasa de flujo de masa a través

de tubos y estaciones compresoras, respectivamente.
Define w=(uT, OT)T el componente uj o Ok es positivo si la dirección de flujo

coinciden con el ducto asignado o la dirección de la estación, negativa, o
cualquiera que sea. Decimos pi es la presión de nodo i.
P=(p1,…..,pn)T, y s=(s1,….,sn)T son el vector fuente, donde la fuente si en el nodo i

es positiva(negativa) si el nodo es un nodo suministro (Entrega). Un nodo que no
es ni de suplencia o entrega es llamado un nodo de transición y si es igual a
cero. Se asume, sin perdidas generalmente, la suma de las fuentes es cero:
∑𝒏𝒊=𝟏 𝒔𝒊 = 𝟎 (12)
29
La ecuación de flujo de red ahora puede ser declarada en la sentencia:
𝑨𝒘 = 𝒔, 𝒚
{ (13)
𝑨𝑻𝒍 ∗ 𝒑𝟐 = 𝝓(𝒖),
Donde P2=(p12 ,….,pn2)T, 𝝓(𝒖) = (𝝓(𝒖𝟏 ), … . , 𝝓𝒍 (𝒖𝒍 ))T, con 𝝓𝒋 (𝒖𝒋 ) = 𝒄𝒋 𝒖𝒋 |𝒖𝒋 |.𝒂
siendo la ecuación de flujo en el tubo j.
Ahora suponiendo que el vector fuente s satisfaciendo la condición de sumatoria

cero, y los limites pl, pv de presiones en todos los nodos son dados. El problema
es poder determinar el vector presión p y el vector flujo w, entonces el total del
consumo de combustible es minimizado, es decir:
Minimizar F(w1p)=∑𝒎
𝒌=𝟏 𝒈𝒌 (𝒗𝒌 , 𝒑𝒊𝒏(𝒌) , 𝒑𝒐𝒖𝒕(𝒌) ), (13a)
Sujeto a 𝑨𝒘 = 𝒔, (13b)
𝑨𝑻𝒍 𝒑𝟐 = 𝝓(𝒖𝟏 ), (13c)
𝒑 ∈ [𝒑𝑳 , 𝒑𝑼 ], (13d)
(𝒗𝒌 , 𝒑𝒊𝒏(𝒌) , 𝒑𝒐𝒖𝒕(𝒌) ) ∈ 𝑫𝒌, 𝒌 = 𝟏, 𝟐, … . . , 𝒎 (13e)
Donde 𝒗𝒌 , 𝒑𝒊𝒏(𝒌) y 𝒑𝒐𝒖𝒕(𝒌) son el flujo de rata de masa, presión de succión y

presión de descarga en la estación compresora k, respectivamente. Esta es en (k)
y out(k) denotando los índices con nodos definiendo las estaciones compresoras
k. En la figura 14, por ejemplo, el compresor 2 es formado por el ducto (3,4) los
cuales hacen mención in(2)=3 y out(2)=4. La función gk es la correspondiente a
la función costo, y Dk es el dominio factible en la estación k para el vector
(𝒗𝒌 , 𝒑𝒊𝒏(𝒌) , 𝒑𝒐𝒖𝒕(𝒌) ), para un estudio profundo de la estructura y propiedades de
Dk y gk. Note lo siguiente:
1. El dominio factible Dk son generalmente no convexos.

2. La función de minimización de combustible gk son no lineales, no convexos y
en ocasiones discontinuos, es decir, que la variable independiente no es de
grado 1, del mismo modo es cóncavo y discontinuo para cualquier k.
3. La ecuación de flujo de tubo (13c) define un conjunto no convexo.
30
Figura 14. Ejemplo de una red simple.
La figura muestra un ejemplo de red simple con n=10 nodos, l=6 tubos y m =3
estaciones.
Se asignan dirección a las tuberías, note que todos los nodos, tubos y
estaciones han sido etiquetados separadamente. Las matrices Al y Am para esta
red es dada por:
Figura 15. Matriz de la red simple
Donde la i-ésima fila en cada matriz corresponde al nodo i, la columna Al

corresponde a los ductos (2,3), (4,5), (5,7), (8,9) y (9,10), respectivamente y la
columna K-ésima en Am a estaciones compresoras K(CSk) en la red. Note que
en cada matriz en cualquier columna contiene exactamente dos elementos no
ceros (1 y -1), los cuales corresponden a los dos nodos finales del ducto o
compresor.
31
Aplicabilidad del modelo al contexto colombiano:
Este modelo puede ser aplicado al contexto colombiano puesto que este se
enfoca a la forma ramal la cual tiene Colombia en su tipología. Para poder
aplicarlo tendríamos que tener los datos necesarios para correr esta simulación
en un software especializado. Los datos necesarios son las variables mostradas
en las restricciones mencionadas anteriormente.
3.1.2. Optimización en el transporte de gas natural con base en el consumo. [21]
La formulación del problema se aplica en una situación donde la empresa

encargada de la transmisión es la misma que comercializa, obteniendo la
posibilidad de decidir la cantidad de gas a comprar de varias fuentes
satisfaciendo la demanda en diferentes nodos, garantizando la presión mínima
en los mismos. Estas adaptaciones sencillas permiten reformular el modelo
para una empresa de transporte de gas puro.
La red de una empresa de gas consta de varios puntos de suministro en el que

el gas se inyecta en el sistema, varios puntos de demanda de gas del sistema y
otros nodos intermedios donde el gas es simplemente redirigido.
Las tuberías están representadas por ductos que unen los nodos. Se utiliza la
siguiente notación matemática. La red se define como el par (N’, A’), donde N’
es el conjunto de nodos y A’ C N’ × N’ es el conjunto de ductos que conectan
estos nodos. La red se puede representar como en la siguiente figura.
Figura 16. Red de conjunto de ductos
32
Dos variables están asociadas a cada i nodo de la red: pi representa la presión
del gas en este nodo y si el suministro de gas neto en el nodo i. Si si es positivo
corresponde a un suministro de gas en el nodo i. Si si es negativa implica un gas
Demanda de di = -si en el nodo i.
El flujo de gas fij está asociado con cada arco (i , j) de i a j .

Hay dos tipos de arcos:
 Arcos pasivos (cuyo conjunto se observa Ap) corresponden a las

tuberías.
 Arcos activos corresponden a tuberías con un compresor (cuyo conjunto
es señalado Aa).
Las restricciones del modelo son los siguientes:
 En un nodo de alimentación i, el gas si de entrada debe permanecer

dentro de los límites de asimilación especificadas en los contratos.
 Una empresa de gas contrata una cantidad específica media diaria que
debe ser tomada por la compañía de transmisión de la compañía
productora.
Dependiendo de la flexibilidad del contrato, la empresa de transmisión

tiene la posibilidad de aumentar una cantidad que oscila entre un límite
inferior y superior de la fracción (por ejemplo, entre 0,85 y 1,15), de la
cantidad promedio contratada.
Matemáticamente:
si ≤ 𝒔i ≤ 𝒔̅𝒊 (14)
En un nodo de la demanda, el gas de salida −𝒔𝒊 debe ser mayor o igual a

𝒅𝒊 , lo que corresponde a la demanda en este nodo.
La compañía de transporte de gas no puede llevar el gas a una presión

superior a la asegurada por el proveedor en el punto de entrada. Por el
contrario, en cada punto de salida, la demanda debe ser satisfecha a una
presión mínima garantizada para el usuario industrial o de la compañía de
distribución local.
Matemáticamente:
pi ≤ 𝒑i ≤ 𝒑̅𝒊 (15)
La ecuación de conservación de flujo en el nodo i (ver figura 16) asegura el

equilibrio de gas en el nodo i.
33
Figura 17. Suministro del nodo i
Matemáticamente:
∑ 𝑓𝑖𝑗 = ∑ 𝑓𝑖𝑗 + 𝑠𝑖
𝑗|(𝑖,𝑗)∈𝐴 𝑗|(𝑖,𝑗)∈𝐴
(16)
Ahora, teniendo en cuenta las limitaciones de los ductos, se distingue entre

los arcos pasivos y activos. Para un arco pasivo, la relación entre el flujo 𝑓𝑖𝑗
en el arco (𝑖, 𝑗) y las presiones 𝑝𝑖 y 𝑝𝑗 es de la siguiente forma:
𝑠𝑖𝑔𝑛(𝑓𝑖𝑗 )𝑓𝑖𝑗2 = 𝐶𝑖𝑗2 (𝑝𝑖2 − 𝑝𝑗2 ), ∀(𝑖, 𝑗) ∈ 𝐴𝑝 (17)
Donde 𝑐𝑖𝑗 es una constante que depende de la longitud, el diámetro y la

rugosidad absoluta de la tubería y en la composición del gas. Tenga en
cuenta que 𝑓𝑖𝑗 tienen restricciones en el signo. Si el flujo 𝑐𝑖𝑗 < 0, −𝑓𝑖𝑗 va
desde el nodo 𝑗 al nodo 𝑖.
Para un arco activo correspondiente a una tubería con un compresor,

existe un aparato que aumenta la presión a lo largo del tubo que permite
un flujo mayor que la que corresponde a la disminución de la presión en la
tubería:
𝑠𝑖𝑔𝑛(𝑓𝑖𝑗 )𝑓𝑖𝑗2 ≥ 𝐶𝑖𝑗2 (𝑝𝑖2 − 𝑝𝑗2 ), ∀(𝑖, 𝑗) ∈ 𝐴𝑎 (18)
Hay un límite superior en la presión en la salida del compresor:
𝑝𝑖 ≤ 𝑝̅𝑗 , ∀(𝑖, 𝑗) ∈ 𝐴𝑎 (19)
En la práctica, también hay un costo asociado al aumento de presión. No se

considera este costo en este caso. Para ductos activos, la dirección del flujo
se fija:
𝑓𝑖𝑗 ≥ 0 , ∀(𝑖, 𝑗) ∈ 𝐴𝑎 (20)
La función objetivo características del problema planteado es minimizar el

coste total de los suministros. Esto se puede escribir:
34
min 𝑧 = ∑ 𝑐𝑗 𝑠𝑗
𝑗∈𝑁𝑠
(21)
Donde 𝑐𝑗 es el precio de compra del gas entregado en el nodo 𝑗 y denota a

𝑁𝑠 el conjunto de nodos de suministro.
Definición de funciones
De este modo el problema de transmisión de gas se puede formular de la

siguiente manera:
min 𝑧 (𝑓, 𝑠, 𝑝) = ∑ 𝑐𝑗 𝑠𝑗
𝑗∈𝑁𝑠
(22)
∑ 𝑓𝑖𝑗 = ∑ 𝑓𝑖𝑗 + 𝑠𝑖 ∀𝑖 ∈ 𝑁′ (22.1)

𝑗|(𝑖,𝑗)∈𝐴 𝑗|(𝑖,𝑗)∈𝐴
𝑠𝑖𝑔𝑛(𝑓𝑖𝑗 )𝑓𝑖𝑗2 = 𝐶𝑖𝑗 (𝑝𝑖2 − 𝑝𝑗2 )
2
∀(𝑖, 𝑗) ∈ 𝐴𝑝 (22.2)
𝑠. 𝑡 𝑠𝑖𝑔𝑛(𝑓𝑖𝑗 )𝑓𝑖𝑗2
≥ 𝐶𝑖𝑗2 (𝑝𝑖2 − 𝑝𝑗2 ) ∀(𝑖, 𝑗) ∈ 𝐴𝑎 (22.3)
𝑠⏟𝑖 ≤ 𝑠𝑖 ≤ 𝑠̅𝑖 ∀𝑖 ∈ 𝑁′ (22.4)
𝑝⏟𝑖 ≤ 𝑝𝑖 ≤ 𝑝̅𝑖 ∀𝑖 ∈ 𝑁′ (22.5)
{ 𝑓𝑖𝑗 ≥ 0 ∀ (𝑖, 𝑗) ∈ 𝐴𝑎 (22.6)
Un simple cambio de variables nos permite eliminar las no linealidades

debidas a las variables de la presión.
Definiendo:
𝜋𝑖 = 𝑝𝑖2 , ∀𝑖 ∈ 𝑁 (23)
La restricción (22.5) es remplazada por:
𝝅i ≤ 𝝅i ≤ ̅̅̅
𝝅𝒊 (24)
𝝅𝒊 son, respectivamente el cuadrado de pi y 𝒑̅𝒊 . Las únicas no

Donde 𝝅i y ̅̅̅
linealidades restantes aparecen en las expresiones 𝑠𝑖𝑔𝑛(𝑓𝑖𝑗 )𝑓𝑖𝑗2 en las
restricciones (22.2) y (22.3).
35
3.2. ESTUDIOS ASOCIADOS CON LA EXPANSIÓN Y REFUERZO DE REDES

EXISTENTES (MINIMIZACIÓN DE COSTOS DE INVERSIÓN PARA EXPANSIÓN A
MEDIANO Y LARGO PLAZO DE LA RED).
3.2.1. Minimizar el consumo de combustible en una red de tuberías. [22]
Existen diversas variaciones de este problema en función de los supuestos del

modelo y de las decisiones que se harán. Uno de los supuestos del modelo es
que el número de unidades compresoras que funciona dentro de cada estación
de compresión es fija. En este modelo se considera esto como una variable de
decisión por lo tanto, el modelo se convierte en un problema no lineal entero-
mixto (MINLP).
El problema se modela normalmente como un problema de red de flujo no

lineal donde las variables de decisión son la relación de flujo de masa en cada
ducto y de la presión en cada nodo.
En este trabajo se presenta un modelo MINLP para el problema de minimizar el

consumo de combustible en una red de tuberías. Las variables de decisión son:
(i) La presión en cada nodo de la red.

(ii) El flujo de masa a través de la tubería.
(iii) El número de unidades de compresor que tienen que estar en dentro de
cada estación.
Las restricciones del modelo son los siguientes:
 Las restricciones limitan el flujo, la presión que pasa a través de las tuberías.
Los cual a su vez limita el uso de combustible usado por las estaciones
compresoras para su óptimo uso al impulsar el gas natural.
 Estas también pueden determinar el número de estaciones compresoras
que pueden ser colocadas a lo largo de la red de transmisión de gas natural.
36
Estos son los supuestos del modelo.
• Se supone que el problema es en estado estable. Esto significa, que nuestro

modelo proporcionará una solución para los sistemas que han estado operando
por una cantidad relativamente grande de tiempo. El Análisis transitorio
debería necesariamente aumentar el número de variables y la complejidad de
este problema.
• La red está equilibrada. Esto significa que la suma de todos los flujos netos en
cada nodo de la red es igual a cero. En otras palabras, el flujo de suministro
total es impulsado completamente al flujo de la demanda total, sin pérdida.
Sabemos que las estaciones de compresión se alimentan con algo del canal de
alimentación impulsado por las tuberías, y para sostener esta hipótesis se
considera el costo de este consumo como un costo adicional de nuestro
modelo llamado costo de oportunidad que representa la cantidad de dinero
que debe gastar si compramos la combustible de terceros.
• Cada ducto de la red tiene una dirección especificada previamente.
• Hay un número pre-determinado de compresores centrífugos idénticas

conectadas en paralelo en cada una de las estaciones de compresión.
MODELO
En este caso, los parámetros y datos se representan con letras mayúsculas,

mientras que las variables se representan en minúsculas.
Parámetros:
Vs: Conjunto de nodos de suplencia

Vd: Conjunto de nodos de demanda
V: Conjunto de todos los nodos en la red
Ap: Conjunto de arcos de tuberías
Ac: Conjunto de arcos de estaciones compresoras
A: Conjunto de todos los arcos en la red; A = Ap ∪ Ac
Uij: Capacidad de arco de la tubería (i,j); (i,j) ∈ Ap
Rij: Resistencia de la tubería (i,j); (i,j) ∈ Ap
Nij: Límite superior sobre el número de unidades de estaciones
compresoras (i,j); (i,j) ∈ Ac
PiL , PiU Límite de presión en cada nodo; L = límite inferior, U= límite superior; i
: ∈V
bi : Flujo de masa neta en cada nodo; bi > 0 si i ∈ Vs, bi < 0 si i ∈ Vd, bi = 0 de
otra manera.
37
Variables:
xij: Flujo de masa en el arco (i,j); (i,j) ∈ A
pi: Presión en el nodo i; i ∈ V
nij: Número de unidades de compresión trabajando en el nodo
(i,j); (i,j) ∈ Ac
Formulación.
Función objetivo
𝑚𝑖𝑛 ∑ 𝑔(𝑖,𝑗) (𝑥𝑖𝑗 , 𝑝𝑖 , 𝑝𝑗 )

{(𝑖,𝑗)∈𝐴}
(25)
Equilibrar la ecuación de flujo de masa en cada nodo, donde:
∑𝑖∈𝑉 𝑏𝑖 = 0 (26)
∑ 𝑥𝑖𝑗 − ∑ 𝑥𝑗𝑖 = 𝑏𝑖
{𝑗|(𝑖,𝑗)|∈𝐴} {𝑗|(𝑖,𝑗)|∈𝐴}
(27)
Capacidad de la tubería
𝑥𝑖𝑗 ≤ 𝑈𝑖𝑗 (𝑖, 𝑗) ∈ 𝐴𝑝 (28)
Dinámica de flujo de gas en cada tubería (estado estable)
𝑝𝑖2 − 𝑝𝑗2 = 𝑅𝑖𝑗 − 𝑥𝑖𝑗

2
(𝑖, 𝑗) ∈ 𝐴𝑝 (29)
Rango de presión
𝑃𝑖𝐿 ≤ 𝑝𝑖 ≤ 𝑃𝑖𝑈 𝑖∈𝑉 (30)
Límites operacionales en cada estación de compresión
𝑥𝑖𝑗
( , 𝑝𝑖 , 𝑝𝑗 ) ∈ 𝐷𝑗 (𝑖, 𝑗) ∈ 𝐴𝑐 (31)
𝑛𝑖𝑗
𝑥𝑖𝑗 , 𝑝𝑖 ≥ 0, 𝑛𝑖𝑗 ∈ {0,1,2, … . , 𝑁𝑖𝑗 } (32)
Es importante mencionar que una estación de compresión se compone de

varios compresores centrífugos, idénticos conectados en paralelo que pueden
ser encendidos o apagados, véase la figura 18.
38
Figura 18. Representación de una estación compresora.
3.3. ESTUDIOS ASOCIADOS CON EL DISEÑO ÓPTIMO DE NUEVAS REDES (MIN. COSTOS DE
INVERSIÓN.)
3.3.1. Localización óptima de las instalaciones de transmisión. [23]
Este modelo contiene 3 aspectos importantes en el problema de

localización, estos son:
a. Problema de topología.
b. Naturaleza de la demanda.
c. Características de la instalación (estación de
combustible).
En este problema, la estructura topológica se considera como una red de

que tiene demanda en varios lugares y las instalaciones son limitadas en
ductos de red y nodos. Los nodos en la red se consideran como puntos de
intersección de los ductos, y también áreas potenciales para las
instalaciones de estaciones compresoras; el tramo entre dos nodos de la
red se considera como ductos.
El segundo aspecto del problema es la naturaleza de la demanda, es decir,

los modelos tradicionales solo tratan de optimizar la demanda de consumo
final por cada nodo de la red, en este la demanda es considerada como
39
continua en los ductos, de hecho, el problema de flujo de demanda se
considera dentro de la red. Además, la demanda es neta y obtenida a partir
del volumen de consumo a través de los ductos de en la red.
El tercer aspecto del problema, son las características de la instalación que

se dan de la siguiente manera: el número de localizaciones potenciales de
las instalaciones es limitado, capacidad de las instalaciones, y el problema
es considerado como múlti-instalaciones. El objetivo del problema es
determinar los lugares óptimos para las instalaciones en los posibles
lugares considerados en la red, de tal manera que fluyan entre los ductos
de la red existente y cumplan plenamente reduciendo los costos al mínimo.
Este costo se considera como una función positiva y no descendente del
número de instalaciones puesto que, si aumentamos el número de
estaciones subirá el costo del consumo de combustible, entre otros.
También el problema mencionado tiene sus propias suposiciones y
limitaciones que se enfrentan a la función objetivo para la minimización de
costos. Los supuestos de problemas son los siguientes:
1- Tipos de clientes considerados en el sistema.

2- La red de ductos tienen jerarquía. Las instalaciones ubicadas (estación de
Gas Natural Comprimido), en los ductos con la prioridad más baja no puede
dar servicio a las demandas de los ductos con una prioridad más alta, pero
esto es posible a la inversa.
3- Asignación de la demanda de gas de los ductos se realiza en base a la
disposición lógica de camino más corto entre los nodos.
4- Propuesta de posibles nodos para la localización de la exposición
instalaciones basadas en la reducción de costos y la disponibilidad.
5- Las instalaciones están capacitadas teniendo en cuenta la máxima
capacidad y la cobertura total de la demanda.
6- Función objetivo es en la forma de número mínimo de instalaciones
teniendo en cuenta la cobertura total de la demanda.
El modelo de programación matemática para nuestro problema es el

siguiente:
Problema de cobertura de la demanda de tubos.
𝑀𝑖𝑛 (𝑓(𝑝)) (33)
40
Sujeto a:
𝐹
∑ 𝐿𝑓 = 𝑃
𝑓=1
(34)
𝑓
∑ 𝑑𝑖,𝑗 ≤ 𝐶𝑓 . 𝐿𝑓 , ∀𝑓 ∈ 𝐹
𝑖,𝑗|𝑖,𝑗∈𝐴𝑓
(35)
𝐹
𝑓
∑ 𝑑𝑖,𝑗 ≥ 𝐷𝑖𝑗 , ∀𝑖, 𝑗
𝑓=1
(36)
𝐶𝑓 ≤ 𝐶𝑚𝑎𝑥 . 𝐿𝑓 , ∀𝑓 ∈ 𝐹 (37)
𝑃 ∈ {𝐼𝑛𝑡𝑒𝑔𝑟𝑎𝑙 𝑛𝑜 𝑛𝑒𝑔𝑎𝑡𝑖𝑣𝑎} (38)
𝐿𝑓 ∈ {0,1}, ∀𝑓 ∈ 𝐹 (39)
𝑓
𝑑𝑖,𝑗 ≥ 0, ∀𝑖, 𝑗, 𝑓 (40)
𝐶𝑓 ≥ 0, ∀𝑓 ∈ 𝐹 (41)
Notación:
𝑃: Número de instalaciones necesarias
𝑓(𝑝): Función no decreciente de p, que indica los costos de construcción de

una instalación (estaciones de GNC) proporcional a su capacidad y distancia de
los principales nodos de la red.
𝐿𝑓 : Es igual a 1 si f es la apertura de la posición y de lo contrario es igual a 0.
F: Propuesta de nodos que se establecen para la localización de las

instalaciones.
𝑓
𝑑𝑖,𝑗 : Parte del ducto i-j que está cubierto por las instalaciones f.
𝐷𝑖𝑗 : La demanda total del ducto i-j de una manera tal que todos los puntos
entre los nodos i y j está cubierta por al menos una instalación.
𝐴𝑓 : El conjunto de todos los ductos, está sujeta a las instalaciones f

considerando la gama de vehículos.
𝐶𝑓 : La capacidad de las instalaciones de f.
𝐶𝑚𝑎𝑥 : Capacidad máxima de las instalaciones.
𝐴: Ductos existentes establecidos.
41
Las restricciones del modelo dicen:
 El costo mínimo de las estaciones compresoras a partir de la ubicación de las

estaciones compresoras.
 El número de estaciones compresoras que hay en el sistema es una restricción,
pues entre más haya más combustible consumirán.
3.4. ESTUDIOS INTEGRADOS GAS Y ELECTRICIDAD. [24]
La integración de los sectores de gas y electricidad naturales ha aumentado

considerablemente en la última década como consecuencia de las centrales térmicas
de gas natural. Varios países de Europa y Latinoamérica con generación a gas han sido
un factor importante en el crecimiento global del consumo de gas natural. Este
aumento del consumo de gas natural es producto de la cantidad de yacimientos de gas
no explorados. Este tipo de energía debe ser mayormente explotado puesto que tiene
varios beneficios como su bajo impacto ambiental y su competitividad económica en
comparación con otros combustibles fósiles.
El vínculo entre los dos sistemas es a través de las plantas eléctricas de gas (plantas de
generación de ciclo combinado o simple) que están conectados a la red de tuberías de
gas natural. Un caso de estudio se presenta en la integración del sistema de prueba
IEEE-14 y la red de gas Belga. El problema de flujo de potencia óptima de la
electricidad-gas integrado se resuelve utilizando un enfoque híbrido que combina el
algoritmo estrategia evolutiva con la de Newton y el método de punto interior. Este
enfoque híbrido toma las ventajas de ambos optimización estrategia evolutiva y
métodos clásicos.
Este modelo busca minimizar el costo total de la operación del sistema de gas y
energía eléctrica. Para este propósito se modela la red de tuberías y el efecto de los
compresores. El funcionamiento de las instalaciones de almacenamiento es tratado de
una manera similar como a un nodo de inyección de gas puesto que se considera un
problema el funcionamiento estático. Se dan las capacidades máxima y mínima de
salida. Esto no se mantiene para fines de planificación de la operación de múltiples
42
etapas en el que el funcionamiento dinámico del almacenamiento debe ser modelado
adecuadamente.
La formulación matemática de flujo de potencia eléctrica del gas integrado esta

formulada de la siguiente manera:
2
𝑀𝑖𝑛 ∑ (𝑎𝑖 + 𝑏𝑖 𝑃𝐺𝑖 + 𝑐𝑖 𝑃𝐺𝑖 ) + ∑ 𝑔𝑖 𝑠𝑖 (42)
𝑖∉Ωℯ 𝑖∈𝑁𝑠
Esta función objetivo está sujeta a las siguientes restricciones:
𝑃𝐺𝑖 − 𝑅𝑒(Ψ𝑖 (𝑉, 𝜃)) = 𝑃𝐿𝑖 ; ∀𝑖 ∈ 𝑁𝐵 (43)
𝑄𝐺𝑖 − 𝐼𝑚(Ψ𝑖 (𝑉, 𝜃)) = 𝑄𝐿𝑖 ; ∀𝑖 ∈ 𝑁𝐵 (44)
𝑉𝑖,𝑚𝑖𝑛 ≤ 𝑉𝑖 ≤ 𝑉𝑖,𝑚𝑎𝑥 ; ∀𝑖 ∈ 𝑁𝐵 (45)
|𝑃𝑖𝑗 (𝑉, 𝜃)| ≤ 𝑃𝑖𝑗,𝑚𝑎𝑥 ; ∀(𝑖, 𝑗) ∈ 𝑁𝐵 (46)
𝑃𝐺𝑖,𝑚𝑖𝑛 ≤ 𝑃𝐺𝑖 ≤ 𝑃𝐺𝑖,𝑚𝑎𝑥 ; ∀𝑖 ∈ 𝑁𝐵 (47)
𝑄𝐺 𝑚𝑖𝑛 ≤ 𝑃𝐺𝑖 ≤ 𝑃𝐺 𝑚𝑎𝑥 ; ∀𝑖 ∈ 𝑁𝐺 (48)
𝑆𝑖 + ∑ 𝑓𝑗𝑖 − ∑ 𝑓𝑖𝑗 − 𝑒𝑖 (𝑃𝐺,𝑖 ) = 𝑑𝑖 ; ∀(𝑖, 𝑗) ∈ 𝑁𝑁 (49)

𝑗/(𝑗,𝑖)∈∅ 𝑗
∈∅
𝑖,𝑗
𝑠𝑖𝑔𝑛(𝑓𝑖𝑗 )𝑓𝑖𝑗2 = 𝐶𝑖𝑗2 (𝑃𝑖2 − 𝑃𝑗2 ) ; ∀(𝑖, 𝑗) ∈ 𝐴𝑃 (50)
𝑠𝑖𝑔𝑛(𝑓𝑖𝑗 )𝑓𝑖𝑗2 ≥ 𝐶𝑖𝑗2 (𝑃𝑖2 − 𝑃𝑗2 ) ; ∀(𝑖, 𝑗) ∈ 𝐴𝑎 (51)
𝑠𝑖,𝑚𝑖𝑛 ≤ 𝑠𝑖 ≤ 𝑠𝑖,𝑚𝑎𝑥 ; ∀𝑖 ∈ 𝑁𝑠 (52)
𝑃𝑖𝑚𝑖𝑛 ≤ 𝑃𝐺𝑖 ≤ 𝑃𝑖𝑚𝑎𝑥 ; ∀𝑖 ∈ 𝑁𝑁 (53)
𝑓𝑖𝑗 ≥ 0 ; ∀(𝑖, 𝑗) ∈ 𝐴𝑎 (54)
Las restricciones dicen lo siguiente:
Donde la ecuación (42) (función objetivo) está compuesta por dos partes; una es el
costo de generación debido a la generación eléctrica y la otra es el costo del suministro
de gas natural. Dentro de las restricciones se encuentran principalmente limitantes
que buscan encontrar parámetros constantes, por ejemplo las ecuaciones (47, 48) son
los límites inferior y superior de potencia activa y reactiva de generación en cada nodo.
Lo anterior por la parte eléctrica y por la parte de gas que es lo que nos interesa
encontramos que las ecuaciones (52, 52) son los límites inferior y superior de
suministro de gas y las capacidades de las presiones en los nodos respectivamente.
43
Este modelo no puede ser aplicado al contexto colombiano puesto que este está
enfocado al sistema de trasmisión de Bélgica, puesto que este sistema es enmallado y
el sistema colombiano es ramal. Para lograr correr esta simulación con los datos
necesarios se requiere que el sistema sea tipo mallas.
3.5. OBSERVACIÓN SOBRE LOS MODELOS
 Para finalizar este capítulo dando un ejemplo sobre los modelos, primero se debe
ver de dónde vienen las ecuaciones o inecuaciones de los modelos de
optimización. Los modelos vienen divididos en dos partes; una función objetivo y
un conjunto de restricciones. La función objetivo está definida como aquella que es
objeto de maximización o minimización en un problema en nuestro caso de
transmisión de gas natural. Las restricciones vienen de las limitantes que tiene cada
uno de los aspectos que intervienen para el desarrollo del problema. Estas se
desprenden de los modelos matemáticos de los aspectos antes mencionados, en
este caso son las estructuras del sistema de transmisión de gas, estos son; la
estación de compresión que esta a su vez compuesta por el compresor de gas,
turbina de gas y enfriador de aire, tubería y nodo de entrega. Para entender más
este proceso veamos el siguiente ejemplo.
Las ecuaciones siguientes tratan la diferencia de presiones en un tramo de tubería,

donde estas generan un empuje o fuerza que lleve el gas del punto i al punto j. Esta
ecuación sale del modelo matemático de la tubería, la cual se muestra en el
capítulo 1; al remitirnos a este, la ecuación 34 muestra la explicación en estado
estacionario de una tubería, y se nota que estas son muy parecidas en forma.
𝑃𝑖2 − 𝑃𝑗2 = 𝐾𝑖𝑗 ∗ 𝑄𝑖𝑗

2
(5)
𝑝𝑖2 − 𝑝𝑗2 = 𝑐𝑖𝑗 ∗ 𝑢|𝑢|.𝑎 (11)
𝑝𝑖2 − 𝑝𝑗2 = 𝑅𝑖𝑗 − 𝑥𝑖𝑗

2
(29)
Esta restricción pertenece a la mayoría de modelos de optimización pues involucra

la presión que ejercen las estaciones compresoras en la tubería para el impulso del
gas, este fenómeno es el que lleva a cabo el transporte del gas.
 Con respecto a las siguientes sumatorias, estas toman un nodo cualquiera donde
la suma de los flujos debe ser cero en ese punto. Esto lleva a recordar las leyes de
Kirchhoff donde estas dicen lo mismo para los circuitos.
44
Las ecuaciones mostradas a continuación demuestran el equilibrio de masa que
debe haber en cada nodo.
∑ 𝑏𝑖 = 0
𝑖∈𝑉
(26)
∑ 𝑥𝑖𝑗 − ∑ 𝑥𝑗𝑖 = 𝑏𝑖
{𝑗|(𝑖,𝑗)|∈𝐴} {𝑗|(𝑖,𝑗)|∈𝐴}
(27)
Lo cual se relaciona directamente con el balance de flujo de masa donde la

diferencia entre el flujo de entrada y el flujo de salida debe ser igual a cero así como
fue mencionada en el capítulo 2:
𝑖=𝑡 𝑗=𝑢 𝑙=𝑤
∑ 𝑄𝑖 − ∑ 𝑞𝑗 − ∑ 𝐷𝑙 = 0
𝑖=1 𝑗=1 𝑙=1
(9)
 El modelo de gas y energía eléctrica integrado busca mejorar los costos de
combustible energía para el abastecimiento de gas y energía puesto al combinar
los ciclos de cada uno de ellos encontramos que la generación térmica a partir de
gas puede ser más barata que con otros compuestos como el petróleo. De este
modo las estaciones compresoras necesitan energía que es suministrada por la
generadora térmica lo que da un beneficio doble. De esta forma optimizar los
ciclos de cada red.
3.6. MODELOS DE PRUEBA PARA TRANSMISIÓN DE GAS NATURAL.
Los modelos de prueba son usados para comprobar el cumplimiento de los requisitos
necesarios para que una red de transmisión funcione adecuadamente y verifique con
los parámetros y restricciones los modelos de optimización que se expone.
Para este punto de tercer capítulo tendremos varios sistemas de prueba donde se
demostrara la funcionalidad de los modelos de optimización del gas natural.
Para los modelos de optimización de la transmisión del gas se han realizado varias
pruebas, donde se busca mejorar el transporte mediante una función objetivo y unas
restricciones. Para poner estos modelos a prueba debemos tener una red conocida y
los datos para poder correr estos datos en un programa Software especializado. Para
dar ejemplo a esta idea tenemos los siguientes modelos de prueba.
45
3.6.1. Modelo de Prueba 1. [24]
Se propone un modelo de optimización del flujo de potencia integrado de redes de gas

natural y electricidad. Un modelo matemático de este problema se formula como un
problema de optimización donde la función objetivo es minimizar el costo de la
operación del sistema de transporte de gas y energía sujeto a las limitaciones del
sistema eléctrico y de tuberías de gas natural. El vínculo entre ambos sistemas es a
través de las plantas eléctricas de gas (plantas de generación de ciclo combinado o
simple) que están conectados a la red de tuberías de gas natural. Para este caso se
presenta un sistema de prueba integrando de la IEEE-14 y de la red de gas de Bélgica.
Las tablas que se mostraran en este documento revelan que esta combinación es una
ventaja sobre otros métodos para resolver el problema de integración de gas y
electricidad.
Cabe aclarar que se considera que el coeficiente calorífico para el gas natural es 9,400
kcal/m3, y 1 BTU es equivalente a 0.252 kcal. Por lo tanto, el factor de conversión para
1 MBTU/mm3 del gas natural es 9,400 / 0,252.
Como se puede observar en la siguiente imagen, la unión de los dos sistemas de

transporte donde la parte eléctrica abastece a la parte de transmisión de gas para que
los equipos funcionen correctamente y sin interrupciones.
Figura 19. Red eléctrica y de gas natural combinada.
46
Vemos que hay dos generadoras que abastecen la red de gas natural para las
estaciones compresoras que son las encargadas de impulsar el gas a través de las
tuberías. A partir de este esquema de la unión de redes, tenemos los siguientes valores
necesarios para hacer simulaciones:
Si máx. Si di ei pi
Nodo
(mm3/día) (mm3/día) (mm3/día) (mm3/día) (mm3/día)
1 17,39 14,687 0,000 0,0000 65,0000

2 12,60 6,546 0,000 0,0000 64,9498
3 0,00 0,000 3,918 0,0000 64,7620
4 0,00 0,000 0,000 0,4739 64,1138
5 7,20 0,382 0,000 0,0000 58,0042
6 0,00 0,000 4,034 0,0000 57,8905
7 0,00 0,000 5,256 0,0000 58,3239
8 33,02 20,574 0,000 0,0000 53,9316
9 0,00 0,000 0,000 0,0000 65,0000
10 0,00 0,000 6,365 0,0000 64,6991
11 0,00 0,000 0,000 0,0000 64,4808
12 0,00 0,000 2,120 0,4142 63,7616
13 1,80 3,590 0,000 0,0000 63,4569
14 1,44 1,407 0,000 0,0000 63,2515
15 0,00 0,000 6,848 0,0000 61,5299
16 0,00 0,000 15,616 0,0000 59,8420
17 0,00 0,000 0,000 0,0000 53,7605
18 0,00 0,000 0,000 0,0000 59,0379
19 0,00 0,000 0,222 0,0000 57,6260
20 0,00 0,000 1,919 0,0000 53,486
Total: 73,45 47.186 46,298 0,8881
Tabla 5. Suministro de gas y presiones para estudio de optimización
Donde; 𝑆𝑖 es el suministro de gas que entra al nodo i que podría ser cualquiera dentro
de la red de transmisión, 𝑑𝑖 es la demanda no eléctrica de gas para la demanda del
nodo 𝑖, 𝑒𝑖 es la demanda de gas natural dirigida a la planta eléctrica de gas hacia el
nodo 𝑖 y por ultimo 𝑝𝑖 es la presión en el nodo 𝑖.
Estos datos son los valores a estudiar para mostrar un suministro óptimo a la hora de
transmitir el gas natural. Las presiones muestran que a través de las tuberías estas se
mantienen constantes a lo largo de la trasmisión, de esta forma llevar el gas de un
punto bajo parámetros constantes.
Por otro lado la siguiente tabla muestra los flujos de gas a través de las tuberías.
47
Capacidad fij
Tubería De Hasta
máx. (mm /día) (mm3/día)
3
1 1 2 125,08 7,344
2 1 2 125,08 7,344
3 2 3 102,13 10,617
4 2 3 102,13 10,617
5 3 4 49,06 17,315
6 5 6 13,15 0,382
7 6 7 16,01 -3,652
8 7 4 19,78 -8,908
9 4 14 33,73 7,933
10 8 9 111,88 18,191
11 8 9 13,65 1,957
12 9 10 55,94 13,425
13 9 10 6,11 6,723
14 10 11 50,03 8,467
15 10 11 6,11 5,315
16 11 12 38,60 12,068
17 12 13 39,56 9,534
18 13 14 111,88 13,124
19 14 15 79,11 22,464
20 15 16 50,03 15,616
21 11 17 9,42 1,714
22 17 18 3,33 1,714
23 18 19 1,71 1,714
24 19 20 6,93 1,919
Tabla 6. Flujo optimo por tuberías, donde fij es el flujo de gas a través de un tramo de tubería
Se ve en la tabla que estos son los flujos óptimos del sistema lo cual lleva a estudiarlos.
Analizando la tabla encontramos que hay varios casos donde los trayectos son los
mismos por ejemplo del nodo 9 al 10, sin embargo se observa que los valores son
completamente diferentes, esto lleva a concluir que las tuberías están paralelas pero la
cantidad de gas que transporta no es la misma. Podrían haber varias suposiciones, la
fricción que tiene las tuberías, pero por esta razón no habría tanta diferencia entre una
y otra, el diámetro de las tuberías podría ser la razón del flujo de gas.
Al ver la red completa de Bélgica encontramos que entre nodos si hay tuberías en
paralelo como se mostraba en la tabla 6. Y que lo que diferencia las cantidades es por
el diámetro de tuberías, pues la imagen demuestra que entre más gruesa la línea es
más grande la tubería.
48
Figura 20. Sistema de transmisión de Bélgica. [25]
3.6.2. Modelo de Prueba 2. [26]
Este es un ejemplo de un sistema de transmisión de gas natural unificado donde

cuenta con etapas, dos de estas son; el yacimiento donde el gas se procesa y
transporta a través de largas distancias de tubería, y la distribución del gas tratado
hacia el usuario industrial o los usuarios residenciales. Esto aplica para lograr la
optimización del rendimiento de la red de ductos de gas natural, de esta forma
maximizar recursos y mejorar los recursos de precio de venta para las diferentes
necesidades de los usuarios.
Para llevar a cabo el cálculo para la red de tuberías de gas natural a gran escala, la red
de tuberías consta de 100 tubos y 98 nodos, con nodos específicos encargados de
abastecer al usuario o inyección de gas a la red, como se muestra en la siguiente
figura:
49
Figura 21. Red de transmisión de gas natural.
La siguiente tabla muestra los datos para el modelo de prueba dando a conocer los
datos principales de los nodos.
Precio Presión, Mpa Flujo, 104m3/d

Nodo Tipo
Yuan/m3 Max Min Max Min
1 0,80 Inyección 4,0 3,0 600 300
3 1,10 Demanda 4,0 2,5 3,5 1,5
7 0,98 Demanda 4,0 2,5 1,5 0,0
8 0,80 Demanda 4,0 2,5 1,0 0,0
10 0,85 Inyección 4,0 2,5 100,0 40,0
29 1,00 Demanda 3,0 2,0 15,0 5,0
30 1,30 Demanda 3,0 2,0 10,0 5,0
32 0,90 Demanda 3,0 2,0 20,0 5,0
48 0,97 Demanda 4,0 2,0 150,0 50,0
49 1,10 Demanda 4,0 2,0 90,0 10,0
50
53 0,80 Inyección 3,0 2,0 10,0 0,0
62 0,99 Demanda 4,0 2,0 20,0 5,0
72 0,85 Inyección 4,0 2,5 320,0 220,0
94 1,00 Demanda 4,0 2,0 150,0 90,0
96 0,85 Inyección 4,0 2,5 370,0 280,0
98 0,80 Inyección 4,0 2,5 140,0 70,0
Tabla 7. Principales nodos de la red.
Con estos datos de gran importancia, define valores a estudiar para la red, tales como
el flujo y la presión ejercida en la red de transmisión. Al saber esto podemos definir
que los modelos analizados por esta tesis pueden ser usados en esta red.
3.6.3. Modelo de prueba 3. [27]
Este modelo de prueba es presentado para optimizar el suministro de gas natural que
abastece a una empresa generadora de energía eléctrica. Este problema tiene cabida
desde su yacimiento hasta su uso final. Los modelos de la red de gas, incluyendo el
efecto de los compresores para ampliar la capacidad de transporte de gas. Este
proceso de optimización se inserta en un modelo de mayor nivel que analiza la
fiabilidad de suministro del sistema eléctrico teniendo en cuenta las características del
sistema de gas natural.
El modelo se basa en la red de transmisión de Bélgica puesto que los lugares

mencionados en la tabla corresponden a las ciudades de este país. Por lo tanto la
Figura 19 también corresponde a este modelo, sin embargo, los datos y el análisis del
modelo aplicado es completamente diferente puesto que para este caso analizaremos
el modelado de la red de gas natural para estudios de confiabilidad de sistemas de
potencia, por lo tanto hay unas diferencias con el modelo de prueba anterior puesto
que se analizan los máximos y mínimos de cada variable. Para entender cuáles son los
datos necesarios para llevar a cabo la optimización debemos tener las siguientes tablas
que comprenden información vital para el sistema analizado.
pi pi
Si min Si máx. di min di máx. ei min ei máx.
Nodo Cuidad min máx.
(mm3/día) (mm3/día) (mm3/día) (mm3/día) (mm3/día) (mm3/día)
(bar) (bar)
1 Zeebrugge 8,87 17,39 0,00 0,00 0,00 0,00 0,00 77,00
2 Dudzele 0,00 12,60 0,00 0,00 0,00 0,00 0,00 77,00
3 Brugge 0,00 0,00 5,88 5,88 0,00 0,00 30,00 80,00
4 Zomergem 0,00 0,00 0,00 0,00 0,00 0,00 0,00 80,00
5 Loenhout 0,00 7,20 0,00 0,00 0,00 0,00 0,00 77,00
6 Antwerpen 0,00 0,00 6,05 6,05 0,00 0,00 30,00 80,00
7 Gent 0,00 0,00 7,88 7,88 0,00 0,00 30,00 80,00
8 Voeren 20,34 33,02 0,00 0,00 0,00 0,00 50,00 66,20
51
9 Bemeau 0,00 0,00 0,00 0,00 0,00 0,00 0,00 66,20
10 Liege 0,00 0,00 9,55 9,55 0,00 0,00 30,00 66,20
11 Wamand 0,00 0,00 0,00 0,00 0,00 0,00 0,00 66,20
12 Namur 0,00 0,00 0,00 0,00 0,00 1,71 0,00 66,20
13 Anderlues 0,00 1,80 0,00 0,00 0,00 0,00 0,00 66,20
14 Peronnes 0,00 1,44 0,00 0,00 0,00 0,00 0,00 66,20
15 Mons 0,00 0,00 10,27 0,00 0,00 0,00 0,00 66,20
16 Blaregnies 0,00 0,00 23,42 0,00 0,00 0,00 50,00 66,20
17 Wanze 0,00 0,00 0,00 0,00 0,00 0,00 0,00 66,20
18 Sinsin 0,00 0,00 0,00 0,00 0,00 0,00 0,00 63,00
19 Arlon 0,00 0,00 0,33 0,33 0,00 0,00 0,00 66,20
20 Petange 0,00 0,00 0,00 0,00 0,00 1,71 25,00 66,20
Tabla 8. Datos de los nodos de la red.
En la tabla vemos varios parámetros tales como los límites del suministro, los límites
muestran un rango para la inyección de gas natural a la red. El parámetro d expone los
límites de demanda que no son eléctricos como gasolina, entre otros. El parámetro e
expone los límites de demanda eléctricos. Y la p señala los límites de presión que son
generados en las tuberías a causa del pasa se gas. En la siguiente tabla veremos los
datos en los que interviene las tuberías, su origen y destino entre otros datos.
Tipo de
Tubería Desde Hasta Kij
tubería
1 Zeebrugge Dudzele 3,012 Pasiva
2 Zeebrugge Dudzele 3,012 Pasiva
3 Dudzele Brugge 2,459 Pasiva
4 Dudzele Brugge 2,459 Pasiva
5 Brugge Zomergem 1,181 Pasiva
6 Loenhout Antwerpen 0,317 Pasiva
7 Antwerpen Gent 0,386 Pasiva
8 Gent Zomergem 0,476 Pasiva
9 Zomergem Peronnes 0,812 Pasiva
10 Voeren Bemeau 2,694 Activa
11 Voeren Bemeau 0,329 Activa
12 Bemeau Liege 1,347 Pasiva
13 Bemeau Liege 0,164 Pasiva
14 Liege Wamand 1,204 Pasiva
15 Liege Wamand 0,147 Pasiva
16 Wamand Namur 0,929 Pasiva
17 Namur Anderlues 0,952 Pasiva
18 Anderlues Peronnes 2,694 Pasiva
19 Peronnes Mons 1,905 Pasiva
20 Mons Blaregnies 1,205 Pasiva
21 Wamand Wanze 0,227 Pasiva
22 Wanze Sinsin 0,080 Activa
52
23 Sinsin Arlon 0,041 Pasiva
24 Arlon Petange 0,167 Pasiva
Tabla 9. Datos de las tuberías de la red.
Donde la constante Kij es una constante que depende de las propiedades de la tubería
(longitud, diámetro y la rugosidad absoluta), en la composición del gas y de las
características del compresor. En la anterior tabla se encuentran los datos de tuberías
y tramos desde lugares específicos a otro punto. Esto con el fin de definir trayectos y
saber cuál es la constante de las tuberías. Vemos además de lo antes mencionado que
tiene tipos de arcos, pasivos y activos, donde los activos corresponden a un flujo
constante de gas y el pasivo es un flujo variable a causa de la dirección del mismo flujo.
53
INTERPRETACION Y CONCLUSIONES
 El sistema de transporte de gas en Colombia es radial por ende los modelos

presentados anteriormente pueden llegar a ser aplicados según sea el caso de
optimización presentado.
 Los modelos de optimización están compuestos por ecuaciones o inecuaciones,
desprendiéndose las restricciones y función objetivo; sin embargo, hay un enfoque
sobre las restricciones que caracterizan los diferentes tipos de modelos existentes
en los sistemas de transmisión de gas natural.
 Las restricciones en los modelos buscan limitar parámetros como la presión, flujo,
entre otras para encontrar cuantificaciones constantes. Consiguiendo así optimizar
las variables de la función objetivo.
 Caracterizando los modelos se hizo un enfoque al uso de parámetros que busquen
un beneficio con el uso eficiente de combustible, mejor ubicación de las
estaciones compresoras entre otros, de este forma enfocar este documento para
que sea usado en redes de transmisión en Colombia.
 El modelo integrado por gas y electricidad busca mejorar los costos de
combustible combinando los ciclos. Abasteciendo de energía a las estaciones
compresoras que vienen del suministro de gas a las generadoras térmicas. Una
solución muy astuta, que optimiza ambos procesos.
 Los modelos de prueba busca obtener los datos necesarios para realizar una
simulación sobre un sistema de transporte puntual así llévalos a un software
especializados para su futuro uso en la red estudiada.
 En Colombia no existe ningún modelo de prueba, por esta razón uno de los
principales objetivos de este documento es formar conciencia en la búsqueda del
mejoramiento de la transmisión de gas natural. Por esta razón la construcción de
este modelo de prueba con los datos necesarios para una futura simulación es una
tarea que se debe realizar.
 Con el análisis de los modelos de optimización de la transmisión de gas natural
encontramos que Colombia tiene que avanzar en la adquisición del conocimiento
sobre el tema, por lo que enfocar estudios sobre el gas natural podría llevar a
minimizar el uso de combustible para sus estaciones compresoras y mejorar el
flujo a través de las tuberías.
 La mayoría de modelos buscan minimizar el costo por el uso de combustible en las
estaciones compresoras, pues según varios modelos el uso del combustible genera
hasta un 8% del costo final de la transmisión, por lo tanto trabajar sobre este
parámetro puede ser muy lucrativo. Según nuestro análisis la integración de la
energía eléctrica con el transporte de gas es una solución productiva puesto que al
combinar los ciclos de estas dos ramas mejoraría la generación de energía
eléctrica y el transporte del gas.
54
 Con este documento buscamos ayudar al conocimiento sobre la optimización de
transmisión de gas natural, adquiriendo así propuestas sobre el tema llevando a
cabo un mayor análisis sobre el tema. Del mismo modo generar teorías de
optimización para el transporte de este elemento.
 Con este compendio de modelos generamos también valores necesarios para
hacer simulaciones que introducidos en un software especializados encontraría
soluciones para una red de transmisión. De este modo buscamos que este
documento sea usado en una futura tesis sobre las simulaciones y optimización de
redes de transmisión.
55
REFERENCIAS
[1] CREG. Comisión de regulación de energía y gas.
“www.creg.gov.co/html/i_portals/index.php?p_origin=internal&p_name=content&p_id=MI-
73&p_options=”
[2] Federal Energy Regulatory Comission, Energy Primer. “www.ferc.gov/market-

oversight/guide/energy-primer.pdf”
[3] Ecopetrol. Ecopetrol. “www.ecopetrol.com.co/contenido.aspx?catID=358&conID=42613”
[4] Agencia Internacional de la Energía. 2012. “ www.geografiainfo.es/tuberias”.
[5] Ecopetrol. Ecopetrol. “www.ecopetrol.com.co/contenido.aspx?catID=358&conID=42585”
[6]Transportadora de Gas Internacional. TGI S.A. E.S.P. “http://www.tgi.com.co/images/mapa-

int/index.html”
[7] Transportadora de Gas Internacional. TGI S.A. E.S.P. “www.tgi.com.co”
[8] Promigas. Promigas S.A. E.S.P. “www.promigas.com”
[9] Progasur. Promotora de Gases del Sur S.A. E.S.P. “www.progasur.com.co”
[10] Transgastol. Transgastol. “www.transgastol.com”
[11] Transoccidente. Transoccidente. “www.transoccidente.com.co”
[12] CREG. Comisión de regulación de energía y gas.

“www.creg.gov.co/html/i_portals/index.php?p_origin=internal&p_name=content&p_id=MI-
76&p_options=”
[13] NASA, 2008; “www.grc.nasa.gov/www/k-12/airplane/mflow.html”.
[14] Universidad de las Américas de Puebla;

“catarina.udlap.mx/u_dl_a/tales/documentos/leip/mateos_e_e/capitulo5.pdf”
[15] M. Chaczykowski⇑, A.J. Osiadacz, F.E. Uilhoorn. 2011. “Exergy-based analysis of gas
transmission system with application to Yamal-Europe pipeline”. Applied Enegy, 2.
[16] Compresores Centrífugos. Tipos. Componentes. Dimensionamiento. Performance.

Funcionamiento anormal, Universidad de Buenos Aires;
“materias.fi.uba.ar/6720/unidad3.PDF”
[17] Compresores axiales. Teoría alar. Efectos viscosos. Pérdidas. Performance, Universidad de
Buenos Aires; “materias.fi.uba.ar/6720/unidad4.PDF”
[18] Abraham Debebe Woldeyohannes, Mohd Amin Abd Majid; “Simulation model for natural
gas transmission pipeline network system”
56
[19] Librería de revista científicas de chile, Análisis en estado estacionario de tuberías,
”www.scielo.cl”
[20] Suming Wv Lic , R. Z. Rios-Mercado, E. A. Boyd Pros, L. R. Scott, 1999, “Model Relaxations
for the Fuel Cost Minimization of Steady-State Gas Pipeline Networks”
[21] Daniel De Wolf, Yves Smeers, 2000, The Gas Transmission Problem Solved by an Extension
of the Simplex Algorithm.
[22] Diana Cobos-Zaleta, Roger Z. Ríos-Mercado, 2002, “Model for a Minimizing Fuel
Consumption on Natural Gas Pipeline Networks”
[23] Abtin Boostani, Reza Ghodsi, Ali Kamali Miab, 2010, “Optimal location of Compressed
Natural Gas (CNG) refueling station using the Arc Demand Coverage Model”
[24] Clodomiro Unsihuay, J. W. Marangon Lima (SM), A. C. Zambroni de Souza; “Modeling the
Integrated Natural Gas and Electricity Optimal Power Flow”
[25] http://www.egssis.com/wp-content/uploads/2014/02/Carte-Fluxys-1-fr-450x400.png
[26] Zhang Huaixin, Yang Yi Beijing, Luo Ming Chong Qing; Optimization of Large-scale Natural
Gas Pipeline Network Base on the Combined GASA Algorithm
[27] Jorge Munoz, Noemi Jimenez-Redondo, Juan Perez and Julian Barquin; “Natural Gas
Network Modeling for Power Systems Reliability Studies”; IEEE.
57
ANEXOS
TABLA 1
Fuente: Tabla 2B. Factores de diseño asignados a los sistemas de tuberías en acero. NTC 3728
58
TABLA 2
Fuente: Tabla 3. Factor de eficiencia longitudinal de las juntas soldadas. NTC 3728
TABLA 3
Fuente: Tabla 4. Factor de degradación por temperatura. NTC 3728
59

Análisis de Modelos de Optimización de Redes de Transmisión Del G

Cargado por

Copyright:

Formatos disponibles

Análisis de Modelos de Optimización de Redes de Transmisión Del G

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Análisis de Modelos de Optimización de Redes de Transmisión Del G

Cargado por

Copyright:

Formatos disponibles

Universidad de La Salle

Ingeniería Eléctrica Facultad de Ingeniería

Análisis de modelos de optimización de redes de transmisión del

Camilo Andrés Moya Garzón

Follow this and additional works at: https://ciencia.lasalle.edu.co/ing_electrica

Part of the Electrical and Electronics Commons

Att: Camilo Andrés Moya Garzón

A mi compañero y amigo de Camilo Moya por toda la dedicación, esfuerzo y

Y por último pero no menos importante a mi director de Tesis el ingeniero Adrián

Att: Andrés Mauricio Macías Valencia.

1.2. Generalidades de la transmisión de Gas Natural………………………..……….6

1.2.1 Gasoductos entre departamentos…………………………………………….7

1.3. Usos del Gas Natural…………………………………………………………….………..…..7

1.4. Estructura de Gasoductos en Sudamérica y Centro América….……….….8

1.5. Empresas transportadoras de gas en Colombia…………………….……….…..9

1.5.1 Transportadora de Gas Internacional – TGI……………….……….…..9

1.6. Mapa de conexiones de SNT……………………………………………….…….…...15

1.7. Estructura del sector del gas……………………………………………….…….…...17

2.2. Sistema de Ecuaciones en Estado Estacionario…………………….…………..23

2.2.1. Formulación de ecuaciones de flujo de tubería………....…………23

MODELOS DE OPTIMIZACIÓN DE TRANSMISIÓN DE GAS NATURAL

3.1.1 Minimización de los costos de combustible al final de la red en las

3.2. ESTUDIOS ASOCIADOS CON LA EXPANSIÓN Y REFUERZO DE REDES

3.2.1 Minimizar el consumo de combustible en una red de tuberías…35

3.3. ESTUDIOS ASOCIADOS CON EL DISEÑO ÓPTIMO DE NUEVAS REDES

3.3.1 Localización óptima de las instalaciones de transmisión….38

3.4. ESTUDIOS INTEGRADOS GAS-ELECTRICIDAD………………………………41

3.6.1. Modelo de prueba 1………………………………………………………….44

En el primer capítulo se muestra lo que hay en el sistema colombiano con respecto a la

El segundo capítulo demuestra el origen de las restricciones y función objetivo a partir

El tercer capítulo muestra los modelos estudiados y analizados, donde muestra la

PLANTEAMIENTO DEL PROBLEMA

Este análisis se hace con la idea de entender los modelos de optimización en la

Se analizaran los modelos de optimización en transmisión de gas natural más

GENERALIDADES DEL MERCADO DEL GAS NATURAL

1.1. Historia del gas en Colombia

En Colombia el desarrollo de la industria del gas natural es reciente. Sin embargo en el

En 1993 el Gobierno Nacional definió que ECOPETROL liderara la interconexión

Entre 1997 y 1998 se otorgaron concesiones de áreas de distribución exclusiva de gas

El Gobierno Nacional interesado en promover el desarrollo de este combustible en

En el 2007 Ecopetrol, PDVSA (petrolera venezolana) y Chevron suscribieron un

1.2. Generalidades de la transmisión del Gas

Dentro de las principales generalidades de las redes de transmisión de gas natural

Hay tres tipos principales de tuberías a lo largo de la ruta de transporte: el sistema de

Los otros dos tipos de tuberías se pueden clasificar como departamentales y

1.2.1. Gasoductos entre departamentos

La red de gasoductos departamentales transporta gas natural procesado en las plantas

El gas natural que se transporta a través de gasoductos departamentales viaja a altas

El modelo matemático del sistema de transmisión de gas es representado por un

1.3. Usos del Gas Natural

El gas natural en el sector industrial y petroquímico puede ser utilizado como

Como combustible se emplea en varios tipos de equipos, por ejemplo: hornos,

Como se muestra en la Figura 1, el sistema de conexión en Sudamérica es poco denso

Figura 1. Mapa de Sudamérica con la estructura de gasoductos [4].

1.5. Empresas transportadoras de gas en Colombia

Las empresas transportadoras de gas natural en Colombia son:

 Transportadora de Gas Internacional – TGI.

A continuación encontraremos una descripción breve de algunas de las distribuidoras

1.5.1 Transportadora de Gas Internacional – TGI.