Elt 3861 Operación Económica Y Planificación de Sistemas Electricos
Elt 3861 Operación Económica Y Planificación de Sistemas Electricos
Elt 3861 Operación Económica Y Planificación de Sistemas Electricos
2.- Resolver el Problema de Optimización No Lineal mediante el método del gradiente del ascenso más
pronunciado: (25%)
Resolución
1.- Resolver el Problema de Optimización No Lineal mediante el método del gradiente del ascenso más
pronunciado: (30%)
𝜕ℒ (𝑥1 , 𝑥2 , 𝜆)
= 6 + 6𝜆𝑥2 = 0 (𝛽)
𝜕𝑥2
𝜕ℒ (𝑥1 , 𝑥2 , 𝜆)
= 3𝑥22 + 2𝑥1 + 𝑥12 − 10 = 0 (𝛾)
𝜕𝜆
Reescribiendo (𝛾), como 𝑥1 es una variable dependiente:
𝑥1 = 𝑓(𝑥2 )
𝑖 𝑥2 𝑥1 𝜆 ∇ℒ
0 1,500000 1,061553 -0,485071 1,634359
1 1,630749 0,738383 -0,575247 0,371499
2 1,660469 0,651827 -0,605390 -0,031389
3 1,657957 0,659377 -0,602636 0,005132
4 1,658368 0,658146 -0,603083 -0,000805
5 1,658304 0,658339 -0,603013 0,000127
6 1,658314 0,658308 -0,603024 -0,000020
7 1,658312 0,658313 -0,603022 0,000003
𝜕𝐹1 𝜕𝑃𝑝𝑒𝑟𝑑
= 8.5 + 0.0060𝑃1 + 𝜆 ( − 1) = 0
𝜕𝑃1 𝜕𝑃1
𝜕𝐹2 𝜕𝑃𝑝𝑒𝑟𝑑
= 8.0 + 0.0084𝑃2 + 𝜆 ( − 1) = 0
𝜕𝑃2 𝜕𝑃2
𝜕𝐹3 𝜕𝑃𝑝𝑒𝑟𝑑
= 6.1 + 0.0080𝑃3 + 𝜆 ( − 1) = 0
𝜕𝑃3 𝜕𝑃3
𝜕𝐹1
= 8.5 + 0.0060𝑃1 + 𝜆(0.00016𝑃1 − 1) = 0
𝜕𝑃1
𝜕𝐹2
= 8.0 + 0.0084𝑃2 + 𝜆(0.00012𝑃2 − 1) = 0
𝜕𝑃2
𝜕𝐹3
= 6.1 + 0.0080𝑃3 + 𝜆(0.00018𝑃3 − 1) = 0
𝜕𝑃3
Reescribiendo, se tiene:
8.5 + 0.0060𝑃1 + 0.00016𝜆𝑃1 − 𝜆 = 0
8.0 + 0.0084𝑃2 + 0.00012𝜆𝑃2 − 𝜆 = 0
6.1 + 0.0080𝑃3 + 0.00018𝜆𝑃3 − 𝜆 = 0
Facultad Nacional de Ingeniería Ingeniería Eléctrica e Ingeniería Electrónica
Iteración 𝑖 = 0
𝜆0 − 8.5
𝑃1 =
0.0060 + 0.00016𝜆0
𝜆0 − 8.0
𝑃2 =
0.0084 + 0.00012𝜆0
𝜆0 − 6.1
𝑃3 =
0.0080 + 0.00018𝜆0
Para 𝜆0 = 9
𝑃1 = 67.2043
𝑃2 = 105.4852
𝑃3 = 301.4553
𝑃𝑝𝑒𝑟𝑑 = 9.2077
Y el error 𝑒0 = 𝑃𝑑 + 𝑃𝑝𝑒𝑟𝑑 − ∑3𝑖=1 𝑃𝑖 = − 264.9371, como el error es negativo significa que hay
exceso de generación, es necesario reducir 𝜆0 en 10 %
𝜆1 = 0.90𝜆0 = 8.10
Iteración 𝑖 = 1
𝑃1 = −54.8246
𝑃2 = 10.6701
𝑃3 = 211.4612
𝑃𝑝𝑒𝑟𝑑 = 4.2717
𝑒1 = 36.965
La proyección Lambda
𝜆1 − 𝜆0 8.1 − 9
𝜆2 = 𝜆1 − 𝑒1 = 8.1— 36.965 = 8.2102
𝑒1 − 𝑒0 36.965 − (−264.9371)
Iteración 𝑖 = 2
𝑃1 = −39.6251
𝑃2 = 22.3965
𝑃3 = 222.64540
𝑃𝑝𝑒𝑟𝑑 = 4.6171
𝑒2 = −0.7997
𝑃1 = 35 𝑀𝑊
𝑃2 = 30 𝑀𝑊
Quedando por determinar la potencia de generación de la unidad 3 que debe cubrir parte de la carga
como las pérdidas. Se debe satisfacer la ecuación del balance de potencia:
Facultad Nacional de Ingeniería Ingeniería Eléctrica e Ingeniería Electrónica
𝑃𝐷 + 𝑃𝑝𝑒𝑟𝑑 − 𝑃1 − 𝑃2 − 𝑃3 = 0
Con 𝑃𝐷 = 200, 𝑃1 = 35 𝑀𝑊 y 𝑃2 = 30 𝑀𝑊
200 + 0.00008 ∙ 352 + 0.00006 ∙ 302 + 0.00009𝑃32 − 35 − 30 − 𝑃3 = 0
0.00009𝑃32 − 𝑃3 + 135.152 = 0
Las soluciones son:
−(−1) ± √(−1)2 − 4 ∙ 0.00009 ∙ 135.152
𝑃3 =
2 ∙ 0.00009
𝑃31 = 10974.2739 𝑀𝑊
𝑃32 = 136.8372 𝑀𝑊
Descartándose, la primera por ser una solución imposible.
𝑃1 = 35 𝑀𝑊
𝑃2 = 30 𝑀𝑊
𝑃3 = 136.8372 𝑀𝑊
$
𝐹𝑇 = 2329.5595
ℎ
b) Realizar el despacho de carga por el método del gradiente de primer orden.
En este caso, no es posible aplicar el método del gradiente de primer orden. Debido a que solamente
queda un solo generador, la unidad 3:
𝑃3 = 200 − 35 − 30 = 135 𝑀𝑊
𝜕 2 𝐹1
𝐹1′′
= = 0.0060
𝜕𝑃12
′′
𝜕 2 𝐹2
𝐹2 = = 0.0084
𝜕𝑃22
𝜕 2 𝐹3
𝐹3′′ = = 0.0080
𝜕𝑃32
Utilizando la ecuación:
𝐹1′′ + 𝐹3′′ 𝐹3′′ ∆𝑃1 𝜆1 − 𝜆3
[ ′′ ] [∆𝑃 ] = − [𝜆 − 𝜆 ]
𝐹3′′ ′′
𝐹2 + 𝐹3 2 2 3
∆𝑃1 −114.734300
[ ]=[ ]
∆𝑃2 −36.714976
𝑃1 = 35 𝑀𝑊
𝑃2 = 30 𝑀𝑊
𝑃3 = 135 𝑀𝑊
𝑃1 = 35 𝑀𝑊
𝑃2 = 30 𝑀𝑊
𝑃3 = 165 𝑀𝑊
Facultad Nacional de Ingeniería Ingeniería Eléctrica e Ingeniería Electrónica