Informe Proyecto 2

El Problema de la Planificación de Unidades de Energı́a Térmica

Análisis de Algoritmos II
Escuela de Ingenierı́a de Sistemas y Computación
Johan Sebastián Tombe Campo - 2110051
César Alejandro Grijalba Zuñiga - 2110035
Marcelo Alejandro Garcı́a Millán - 1941427
Alejandro Chara Garcia - 2010017
November 2023

1 ¿Entendimos el problema?
- Según la primera solución ofrecida:

Efectivamente, cumple con todos los requisitos y su costo serı́a 191.

- Según la segundo solución ofrecida:

Efectivamente, cumple con todos los requisitos y su costo es 202.

- y por ultimo tenemos esta posible solución:

1
 No cumple con la restricción del lı́mite superior de ampliación, pues la UTPEE 1 del perı́odo 1 al 2 pasa
de una potencia 70 a 350, teniendo un aumento de 280, lo cual supera los 200 que tenia definido como lı́mite
de ampliación.

2 El Modelo
Definimos el siguiente modelo para el problema de la planificación de unidades de la energı́a térmica
(PUEnTe)

2.1 Parámetros
• J : representa el número de UTPEEs.
• K : representa el número de intervalos.
• Ej : representa el costo de encender la UTPEE j. (∀j ∈ {1, . . . , J})
• Aj : representa el costo de apagar la UTPEE j. (∀j ∈ {1, . . . , J})

• Gj : representa el indicador de si la UTPEE terminó generando energı́a o no en el ciclo anterior. j.

(∀j ∈ {1, . . . , J})
• Fj : representa el costo fijo de la UTPEE j. (∀j ∈ {1, . . . , J})
• Vj : representa el costo variable j. (∀j ∈ {1, . . . , J})

• P j : representa el lı́mite inferior de potencia de la UTPEE j. (∀j ∈ {1, . . . , J})

• P j : representa el lı́mite superior de potencia de la UTPEE j. (∀j ∈ {1, . . . , J})

• Supj : representa el lı́mite superior de ampliación de potencia de la UTPEE j. (∀j ∈ {1, . . . , J})
• Infj : representa el lı́mite inferior de la reducción de potencia de la UTPEE j. (∀j ∈ {1, . . . , J})

• Pj0 : representa la potencia generada por la UTPEE j justo antes del primer perı́odo. (∀j ∈ {1, . . . , J})
• Dk : representa la demanda global de energı́a para cada perı́odo k. (∀k ∈ {1, . . . , K})
• Rk : representa la reserva para cada perı́odo k. (∀k ∈ {1, . . . , K})

2.2 Variables
• Pj,k : representa la energı́a que generó la UTPEE j en el perı́odo k. (∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , K})
• Onj,k : representa si la UTPEE j se encendió en el perı́odo k. (∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , K})

• Of fj,k : representa si la UTPEE j se apagó en el perı́odo k. (∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , K})
• Genj,k : representa si la UTPEE j generó energı́a en el perı́odo k. (∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , K})

2
2.3 Restricciones
2.3.1 No Negatividad

(Pj,k ≥ 0), (∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , K}

2.3.2 Cuando una UTPEE está generando energı́a, la potencia generada no puede ser inferior
a un mı́nimo permitido

(Pj,k ≥ P j ), (∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , K}

2.3.3 Cuando una UTPEE está generando energı́a, la potencia generada no puede ser superior
a un máximo permitido

(Pj,k ≤ P j ), (∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , K}

2.3.4 Si la UTPEE j generó energı́a en perı́odos consecutivos, entonces el cambio de gen-

eración de potencia entre el perı́odo k-1 y el perı́odo k no puede ser superior a un
lı́mite de ampliación establecido

Require: ∀j ∈ {1, . . . , J}, ∀k ∈ {2, . . . , K}

if Pj,k−1 > 0 ∧ Pj,k > 0 then
(Pj,k − Pj,k−1 ) ≤ Supj
end if

Require: ∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , 1}

if Pj0 > 0 ∧ Pj,k > 0 then
(Pj,k − Pj0 ) ≤ Supj
end if

2.3.5 Si la UTPEE j generó energı́a en perı́odos consecutivos, entonces el cambio de gen-

eración de potencia entre el perı́odo k-1 y el perı́odo k no puede ser inferior a un lı́mite
de reducción establecido

Require: ∀j ∈ {1, . . . , J}, ∀k ∈ {2, . . . , K}

if Pj,k−1 > 0 ∧ Pj,k > 0 then
(Pj,k−1 − Pj,k ) ≤ Infj
end if

3
Require: ∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , 1}
if Pj0 > 0 ∧ Pj,k > 0 then
(Pj0 − Pj,k ) ≤ Infj
end if

2.3.6 Si la UTPEE j NO generó energı́a en un perı́odo k-1 y SÍ generó energı́a en el perı́odo
k, significa que la UTPEE j fue encendida en el perı́odo k

Require: ∀j ∈ {1, . . . , J}, ∀k ∈ {2, . . . , K}

if Pj,k−1 = 0 ∧ Pj,k > 0 then
Onj,k = 1
else
Onj,k = 1
end if

Require: ∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , 1}

if Pj0 = 0 ∧ Pj,k > 0 then
Onj,k = 1
else
Onj,k = 1
end if

2.3.7 Si la UTPEE j SÍ generó energı́a en un perı́odo k-1 y NO generó energı́a en el perı́odo
k, significa que la UTPEE j fue apagada en el perı́odo k

4
Require: ∀j ∈ {1, . . . , J}, ∀k ∈ {2, . . . , K}
if Pj,k−1 > 0 ∧ Pj,k = 0 then
Of fj,k = 1
else
Of fj,k = 1
end if

Require: ∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , 1}

if Pj0 > 0 ∧ Pj,k = 0 then
Of fj,k = 1
else
Of fj,k = 1
end if

2.3.8 Indicador de si la UTPEE j generó energı́a en el perı́odo k

Require: ∀j ∈ {1, . . . , J}, ∀k ∈ {1, . . . , K}

if Pj,k > 0 then
Genj,k = 1
else
Genj,k = 1
end if

2.3.9 La suma de la potencia de las UTPEEs en un perı́odo k debe ser igual a la demanda
global de energı́a para ese perı́odo

J
X
(Pj,k ) = Dk , ∀k ∈ {1, . . . , K}
j=1

2.3.10 Indica la potencia total disponible en el perı́odo k

J
X
(P j × Genj,k ) ≥ (Dk + Rk ), ∀k ∈ {1, . . . , K}
j=1

2.4 Función Objetivo

Se busca minimizar el costo de generación de energı́a respetando las restricciones anteriormente mencionadas
J X
X K
(Fj × Genj,k + Vj × Pj,k + Ej × Onj,k + Aj × Of fj,k )
j=1 k=1

3 Pruebas
Se hizo uso de diferentes casos de prueba que definió el grupo con configuraciones que iban desde datos
pequeños a otro de dimensión algo más grande, esto con el fin de evaluar el comportamiento del modelo y
asegurar que los resultados arrojados eran correctos.

Entrada : J = 4, K = 10

5
 Entrada : J = 4, K = 8

Table 1: Instancia 1 (J = 4, K = 8)
U T P EE j 1 2 3 4
Ej 10 48 50 25
Aj 0.5 0.24 0.25 0.48
Gj 0 0 0 0
Fj 5 10 6 3
Vj 0.5 0.1 0.6 0.3
Pj 20 40 50 25
Pj 200 400 500 250
Supj 100 200 150 150
Infj 100 80 130 130
Pj0 0 0 0 0

P eríodo k 1 2 3 4 5 6 7 8
Dk 350 500 400 350 600 800 350 225
Rk 35 50 40 35 60 60 80 35

Salida : Costo = 821.73 (Solución Óptima)

Table 2: Instancia 2 (J = 4, K = 10)

U T P EE j 1 2 3 4
Ej 97 31 75 54
Aj 1.4 0.4 1.0 0.3
Gj 0 0 0 0
Fj 65 78 89 30
Vj 0.246 0.114 0.616 0.694
Pj 15 23 12 43
Pj 337 452 379 372
Supj 388 390 130 248
Infj 101 50 107 162
Pj0 0 0 0 0

P eríodo k 1 2 3 4 5 6 7 8 9 10
Dk 873 358 424 746 372 363 872 380 788 655
Rk 33 116 91 103 120 23 39 115 45 53

Salida : Costo = 2951.892 (Solución Óptima)

6
 Entrada : J = 5, K = 11

Table 3: Instancia 3 (J = 5, K = 11)

U T P EE j 1 2 3 4 5
Ej 60 120 38 95 50
Aj 1.8 0.6 0.8 1.9 0.8
Gj 0 0 0 0
Fj 79 90 61 68 88
Vj 0.714 0.741 0.193 0.233 0.832
Pj 41 57 91 46 99
Pj 588 565 256 396 270
Supj 72 290 383 128 188
Infj 222 104 180 142 159
Pj0 0 0 0 0

P eríodo k 1 2 3 4 5 6 7 8 9 10 11
Dk 733 935 989 871 618 523 363 901 635 744 344
Rk 85 65 73 91 20 70 51 68 91 53 114

Salida : Costo = 4602.20 (Solución Óptima)

Entrada : J = 7, K = 15

Table 4: Instancia 4 (J = 7, K = 15)

U T P EE j 1 2 3 4 5 6 7
Ej 70 57 117 116 74 113 31
Aj 1.0 1.9 1.7 1.0 1.9 0.1 0.3
Gj 0 0 0 0 0 0 0
Fj 76 26 77 70 22 46 76
Vj 0.763 0.442 0.367 0.548 0.584 0.692 0.761
Pj 41 50 85 53 54 64 93
Pj 169 383 287 408 353 555 429
Supj 385 349 261 82 230 200 365
Infj 215 254 56 177 114 89 258
Pj0 0 0 0 0 0 0 0

7
P eríodo k 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Dk 840 839 850 445 761 960 987 366 823 726 852 759 407 901 711
Rk 76 31 45 53 30 47 53 115 35 41 23 59 48 100 58

Salida : Costo = 7240.983 (Solución Óptima)

Table 5: Instancia 5 (J = 18, K = 3)

U T P EE j 1 2 3 4 5 6 7 8 9
Ej 51 108 60 69 68 114 100 59 40
Aj 1.4 1.5 0.1 1.6 1.9 0.8 0.8 1.4 1.6
Gj 0 0 0 0 0 0 0 0 0
Fj 23 66 40 70 37 35 53 27 85
Vj 0.759 0.290 0.874 0.220 0.210 0.434 0.175 0.217 0.444
Pj 21 25 12 32 32 41 29 70 97
Pj 529 522 196 163 235 183 317 193 371
Supj 225 302 213 281 381 265 80 60 276
Infj 284 299 187 166 218 118 138 104 161
Pj0 0 0 0 0 0 0 0 0 0

U T P EE j 10 11 12 13 14 15 16 17 18
Ej 91 45 115 33 99 93 77 101 48
Aj 1.3 1.4 1.3 1.3 0.2 0.4 1.3 0.7 1.2
Gj 0 0 0 0 0 0 0 0 0
Fj 23 69 75 51 49 25 57 62 62
Vj 0.759 0.290 0.874 0.220 0.210 0.434 0.175 0.217 0.444
Pj 45 81 74 72 98 92 17 76 32
Pj 527 378 286 422 419 521 550 223 401
Supj 178 110 95 294 55 233 328 229 367
Infj 228 95 192 77 79 192 269 191 204
Pj0 0 0 0 0 0 0 0 0 0

P eríodo k 1 2 3
Dk 958 570 787
Rk 115 63 93

Salida : Costo = 1256.705 (Solución Óptima)

8
4 Análisis de Resultados
Las pruebas realizaran arrojaron resultados óptimos, lo que demuestra que el modelo ha sido bien imple-
mentado, al ser un problema de programación lineal con variables enteras, tiene la limitación de que con
entradas J y K muy grandes, no pueda encontrar solución en un tiempo factible.

5 Video Explicativo
El video explicativo se encuentra en el siguiente link: Link del Video

6 Conclusiones
En este proyecto hemos logrado desarrollar un modelo capaz de encontrar soluciones óptimas para difer-
entes instancias del problema, como se evidenció anteriormente en las salidas proporcionadas para diversas
configuraciones de las UTPEEs y perı́odos.
Es claro que el punto clave en este modelo es la correcta implementación de restricciones que garanticen
la viabilidad de las soluciones generadas y de acuerdo a las diferentes pruebas realizadas con los diferentes
conjuntos de datos proporcionados en la Baterı́a, logramos evidenciar que es posible garantizar soluciones
optimas y eficientes con nuestro modelo.
Tambien se vale destacar que se logra comprender un modelo lo suficientemente flexible para adaptarse a
instacias con amplia variacion en los parametros, que aunque en ocasiones pueda tomar tiempo esperar una
respuesta, se logra quedar satisfecho.
En conclusión, se logró cumplir el objetivo de diseñar un modelo robusto y eficiente para abordar el
problema de optimizacion en la generación de energı́a electrica, cabe aclarar que es importante continuar
evaluando su desempeño y quizá lograr realizar mejores ajustes que permitan encontrar soluciones en mejores
tiempos.