ENTREGA SEMANA 5

INVESTIGACIÓN Y OPERACIONES

PROFESOR: FRANKLIN ROLANDO AGUILAR TRIANA

INTEGRANTES DEL GRUPO:

CARLOS ANDRES OSORIO ARANGO

INSTITUCIÓN UNIVERSITARIA POLITÉCNICO GRANCOLOMBIANO

2023
 TABLA DE CONTENIDO

1. LOGO................................................................................................................ 3
2. ANÁLISIS ESTADÍSTICO DEMANDA..............................................................3
3. ANÁLISIS ESTADÍSTICO PRODUCCIÓN.......................................................4
4. MODELO DE LA RED DE SUMINISTRO DE GAS (GNV)...............................5
5. MODELO ALGEBRAICO..................................................................................6
6. PROGRAMA EN GAMS....................................................................................9
 1. LOGO

2. ANÁLISIS ESTADÍSTICO DEMANDA

Demanda
1,400

1,200

1,000

800

600

400

200

-
1 9 17 25 33 41 49 57 65 73 81 89 97 105 113 121 129 137 145 153 161 169 177

De acuerdo con las gráfica, se hace notorio que los datos presentan una tendencia
creciente de forma débil, pero igual la demanda crece año tras año pues cada barra de
acuerdo al mes es superior a su homólogo del año anterior, también se evidencia el
comportamiento cíclico de los datos en un periodo de cada 12 meses, en la gráfica se
observa como cada año se compone de dos picos al alza y dos picos a la baja, para
aclarar, cada año se observa el mismo comportamiento teniendo como los meses de
mayor demanda noviembre y diciembre y los de menor demanda marzo y mayo, pero
además existe una estacionalidad que se evidencia en la baja de la demanda para los
primeros meses donde a mitad de año se marca aún más, pero además se evidencia que
los datos tienen sus picos en los mismos periodos estacionales donde su diferencia a la
línea de tendencia es similar. Ahora bien, teniendo en cuenta que existen estas tres
características en la serie de datos es importante determinar como método de pronóstico
holt-winters pues este considera nivel, tendencia y estacionalidad de la serie, pero lo
principal es que los valores de los factores λ, β y α se pueden optimizar, minimizando el
error del modelo y ajustando cada una de las características a las necesidades de la
serie.

3. ANÁLISIS ESTADÍSTICO PRODUCCIÓN

Producción
1,000

900

800

700

600

500

400

300

200

100

-
1 8 15 22 29 36 43 50 57 64 71 78 85 92 99 106 113 120 127 134 141 148 155 162 169 176

De acuerdo con la gráfica de la producción, también se observa que existe una tendencia
creciente leve para esta serie y cabe destacar que la tendencia creciente es menor a la de
la demanda, pero de acuerdo a la gráfica de la serie y al igual que los datos de la
demanda se presenta una ciclicidad en los datos por periodos de 1 año, donde los meses
de mayor producción son de octubre a diciembre y los de menor producción son junio,
mayo y septiembre, además, se presenta una estacionalidad para periodos de tres
meses, durante estos periodos primero se establece la más corta tendencia en los
primeros meses donde el comportamiento es estable, para después del cuarto mes se
genera una fuerte caída hasta el mes de julio que es seguida por un aumento en la
producción hasta final del año, es por lo anterior que para la producción se hace mucho
más notoria la estacionalidad de los datos pues el comportamiento en meses y su igual
comportamiento es más fácil de observar, por eso en razón a las mismas características
de la demanda se selecciona para el pronóstico el método Winter en razón a que permite
manejar estas tres características con mayor facilidad, como razón particular es
importante mencionar que para la producción se ajusta el componente lambda (ʎ) en
razón a la periodicidad y la estacionalidad observada.
De acuerdo con los métodos de pronósticos seleccionados se procede a realizar el
pronóstico del cual se obtienen los siguientes resultados.

Producción Demanda
2016 2016
700 366
707 366
698 216
699 341
601 211
701 399
692 478
690 531
414 807
896 957
897 1001
907 1197

4. MODELO DE LA RED DE SUMINISTRO DE GAS (GNV)

Este modelo utilizado para el caso de estudio es un modelo multiperiodo y multinivel que
se característica por la capacidad de producción del nivel cero (0) para cada uno de los
diez periodos del análisis.
Los niveles que componen este modelo son:
Los campos de producción (Cusiana), el gasoducto Cusiana-Apiay, el gasoducto Apiay-
Bogotá y el Sistema de Distribución Bogotá (Estaciones de Servicio de GNV)
Y los periodos son los meses desde marzo a diciembre de 2016 obtenidos anteriormente.
Las variables del modelo que se presentan son:
Y= Capacidad de producción en un periodo especifico. (desde el mes 1 al 10)
X=Unidades que es envían de un nivel al siguiente. (desde el nivel 1 al 3)
I=Unidades que se envían desde un periodo a otro. (desde el periodo 1 al 9)
D=La demanda de cada periodo. (desde el mes 1 al 10)

El transporte de producto que se representa con la variable X se presenta en forma

vertical, mientras el almacenamiento y envio al siguiente periodo representado con la
variable I se establece de forma horizontal. Cabe destacar que los parámetros de
producción y demanda ya se encuentran presentados con sus valores respectivos.

5. MODELO ALGEBRAICO
Para modelar el sistema se deben encontrar e identificar todas sus partes.

Criterio de optimización
Como la idea de la empresa es reducir sus costos para atender la demanda de sus
clientes en el 2016 el criterio de optimización que se utilizara será de minimización de los
costos de producción, almacenamiento e inventario .

Conjuntos
Niveles (l): (0,1,2,3) Periodos (g):(1,2,3,4,5,6,7,8,9,10)

Variables
Las variables del modelo son:
Y(g) = La cantidad de gas en GigaBTU a producir en un periodo (g).
Donde g = 1,2,3,4,5,6,7,8,9,10
X(g,l) = La cantidad de gas en GigaBTU transportada desde el nivel (l + 1) en el periodo
(g).
Donde g = 1,2,3,4,5,6,7,8,9,10 y l = 1,2,3
I(g,l) = La cantidad de gas en GigaBTU a almacenar en el nivel (l) llegando en el periodo
(g).
Donde g = 1,2,3,4,5,6,7,8,9 y l = 1,2,3,4

Parámetros
D(g) = Demanda de gas en un periodo (g).
Donde g = 1,2,3,4,5,6,7,8,9,10
M(g) = Capacidad de producción en un periodo (g).
Donde g = 1,2,3,4,5,6,7,8,9,10
P(g) = Costo de una orden de producción en el periodo t en $US
Donde g = 1,2,3,4,5,6,7,8,9,10
T(g) = Costo de transportar una GigaBTU del nivel l hasta el nivel l + 1 en el periodo en
$US
Donde g = 1,2,3,4,5,6,7,8,9,10 y l = 1,2,3
A(g) = Costo de almacenar un GigaBTU en el periodo t y nivel l en $US
Donde g = 1,2,3,4,5,6,7,8,9 y l = 1,2,3,4

Función Objetivo
Minimización
Z =
1 2
1900 y 1+ 1300 y 2+1700 y 3 +1500 y 4 +1300 y 5 +2000 y 6 +1500 y 7 +2000 y 8 +1900 y 9+ 1700 y 10+3 x 1+ 4.5 x 1+ 3 x
Esta función se divide en costos de producción cuyas variables son Y, costos de
transportes cuyas variables con X, y los costos de almacenamiento o inventario que se
presentan con la variable I.
Restricciones
La familia inicial son las restricciones de inventario inicial, su utilidad radica en informarle
al sistema que al inicio de cada nivel el inventario inicial siempre será cero..
1 1 1 1
I 0=0 I 1=0 I 2=0 I 3=0
Ahora la presentación de las restricciones de entrada, esta familia de restricciones lo que
busca es que todo el producto que salga del arco de producción sea enviado al primer
nivel, con esto se alcanza el objetivo de que todo lo producido sea utilizado por el modelo,
se utiliza cada restricción para cada periodo.
1 1 0 1 1 1 1 1 1
X 1 + I 1=Y 1 + I 0 X 2 + I 2=Y 2 + I 1 X 3 + I 3=Y 3 + I 2
1 1 1 1 1 1 1 1 1 1 1 1
X 4 + I 4=Y 4 + I 3 X 5 + I 5=Y 5 + I 4 X 6 + I 6=Y 6+ I 5 X 7 + I 7=Y 7+ I 6
1 1 1 1 1 1 1 1 1
X 8+ I 8=Y 8+ I 7 X 9+ I 9=Y 9+ I 8 X 10 + I 10=Y 10+ I 9

Las siguientes restricciones buscan satisfacer el balance de todos los nodos siguientes al
nivel 1, en estos se buscan que todo lo que ingresa al nodo salga, por medio de la
igualdad entre lo que ingresa al nodo por transporte y almacenamiento y lo que sale por
los mismos motivos.
1 1 1 2 2 2 1 3 3 3 2
X 2 + I 2=X 1+0 X 2 + I 2=X 1+ I 2 X 2 + I 2= X 1+ I 2
4 4 4 3 5 5 5 4 6 6 6 5
X 2 + I 2= X 1 + I 2 X 2 + I 2= X 1+ I 2 X 2 + I 2= X 1 + I 2
7 7 7 6 8 8 8 7 9 9 9 8
X 2 + I 2= X 1 + I 2 X 2 + I 2= X 1 + I 2 X 2 + I 2= X 1 + I 2
10 10 9 1 1 1 2 2 2 1
X 2 +0=X 1 + I 2 X 3 + I 3= X 2+ 0 X 3 + I 3= X 2+ I 3
3 3 3 2 4 4 4 3 5 5 5 4
X 3 + I 3= X 2+ I 3 X 3 + I3 = X 2 + I 3 X 3 + I 3= X 2+ I 3
6 6 6 5 7 7 7 6 8 8 8 7
X 3 + I 3= X 2 + I 3 X 3 + I 3= X 2 + I 3 X 3 + I 3= X 2 + I 3
9 9 9 8 10 10 10 9
X 3 + I 3= X 2 + I 3 X 3 + I 3 =X 2 + I 3

Ahora se muestra la familia de restricciones de salida en este grupo se verifica que la

demanda sea satisfecha por el modelo, se realiza una por cada periodo y verifica que la
demanda del periodo sumado al almacenamiento del periodo sea igual a loque se envía y
almacena desde el periodo anterior.
1 1 2 1 2 3 2 3
X 3 +0=D1 + I 4 X 3 + I 4 =D2+ I 4 X 3 + I 4 =D3 + I 4
4 3 4 5 4 5 6 5 6
X 3 + I 4=D 4 + I 4 X 3 + I 4 =D5 + I 4 X 3 + I 4 =D6 + I 4
7 6 7 8 7 8 9 8 9
X 3 + I 4 =D7 + I 4 X 3 + I 4 =D8 + I 4 X 3 + I 4 =D1 + I 4
10 9
X 3 + I 4= D 1+ 0

La siguiente familia de restricciones establece la capacidad de producción del arco en

todos los periodos con estas restricciones que se presentan para cada periodo se busca
que no se produzcan más unidades que las capaces de producción por el sistema en un
periodo determinado.

Y1≤ M1 Y2 ≤ M2 Y3 ≤ M3 Y4 ≤ M4 Y5≤ M5
Y6≤ M6 Y7 ≤ M7 Y8≤ M8 Y9 ≤ M9 Y 10 ≤ M 10
Por último, la restricción de no negatividad, en donde se obliga a que todas las variables
de decisión sean positivas.
I I
Yt≥0 X t ≥0 It ≥ 0

Proposiciones del sistema

Para el sistema se cumplen tres proposiciones especiales.
Proposición 1. La primera proposición define que una solución óptima del problema puede
ser descompuesta en pequeños subplanes, entendiendo un subplan como la satisfacción
de la demanda de un periodo inferior al total analizado por medio de la utilización de un
grupo especifico de periodos producidos, así las cosas, la sumatoria de subplanes origina
una solución optima del sistema.
Proposición 2. Esta proposición dice que en todo subplan existe como máximo un arco de
producción libre, es decir, cada subplan tiene como máximo un arco de producción que no
utiliza toda su capacidad de producción.
Proposición 3. Esta última proposición, argumenta que dentro del funcionamiento del
sistema siempre es posible encontrar un nodo donde la producción acumulada de
periodos anteriores sea igual al valor de producto enviado y almacenado del mismo nodo,
teniendo en cuenta que en algún momento el producto debe ser enviado o almacenado
para satisfacer la demanda de los diferentes periodos.

6. PROGRAMA EN GAMS

La formulación del problema en GAMS se realizó de forma directa en el programa, es

decir, se incluyó cada uno de los parámetros, variables, funciones y restricción de forma
manual, la función objetivo se presentó por medio de la sumatoria de los productos para la
producción, almacenamiento y transporte del sistema.
Presentación de variables

Parametros del modelo

Creacion de ecuaciones

Funcion objetivo

Y restricciones.
La solución que arroja el programa es la siguiente :