Elt 3861 Operación Económica Y Planificación de Sistemas Electricos

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 7

Facultad Nacional de Ingeniería Ingeniería Eléctrica e Ingeniería Electrónica

ELT 3861 OPERACIÓN ECONÓMICA Y PLANIFICACIÓN DE SISTEMAS ELECTRICOS


Primer Examen Parcial. Oruro, 06 de noviembre de 2020.

1.- Explicar el concepto de economía de escala en sistemas eléctricos. (5%)

2.- Resolver el Problema de Optimización No Lineal mediante el método del gradiente del ascenso más
pronunciado: (25%)

max 2𝑥1 + 6𝑥2 + 10


𝑠. 𝑎:
3𝑥22 + 2𝑥1 + 𝑥12 − 10 = 0
𝑥2 variable independiente,
𝑥1 variable dependiente.

3.- Sean los datos de los generadores de un sistema: (70%)

𝐹1 = 250 + 8.5𝑃1 + 0.003𝑃12 35 ≤ 𝑃1 ≤ 300 𝑀𝑊


𝐹2 = 380 + 8.0𝑃2 + 0.0042𝑃22 30 ≤ 𝑃2 ≤ 250 𝑀𝑊
𝐹3 = 245 + 6.1𝑃3 + 0.0040𝑃32 45 ≤ 𝑃3 ≤ 400 𝑀𝑊
Pérdidas:
𝑃𝑝𝑒𝑟𝑑 = 0.00008𝑃12 + 0.00006𝑃22 + 0.00009𝑃32 𝑀𝑊

Para una demanda de 𝑃𝐷 = 200 𝑀𝑊

a) Realizar el despacho de carga por el método iterativo lambda.


b) Realizar el despacho de carga por el método del gradiente de primer orden.
c) Realizar el despacho de carga por el método del gradiente de segundo orden.
d) Para una demanda de 𝑃𝐷 = 230 𝑀𝑊, aplicar el método de factores de participación

Resolución

1.- Resolver el Problema de Optimización No Lineal mediante el método del gradiente del ascenso más
pronunciado: (30%)

max 2𝑥1 + 6𝑥2 + 10


𝑠. 𝑎:
3𝑥22 + 2𝑥1 + 𝑥12 − 10 = 0
𝑥2 variable independiente,
𝑥1 variable dependiente.

La función de Lagrange, es:


ℒ(𝑥1 , 𝑥2 , 𝜆) = 2𝑥1 + 6𝑥2 + 10 + 𝜆(3𝑥22 + 2𝑥1 + 𝑥12 − 10)

Las condiciones de optimalidad, son:


𝜕ℒ (𝑥1 , 𝑥2 , 𝜆)
= 2 + 𝜆(2 + 2𝑥1 ) = 0 (𝛼)
𝜕𝑥1
Facultad Nacional de Ingeniería Ingeniería Eléctrica e Ingeniería Electrónica

𝜕ℒ (𝑥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 )

𝑥12 + 2𝑥1 + (3𝑥22 − 10) = 0


−𝑏±√𝑏2 −4𝑎𝑐
Considerando la ecuación general: 𝑎𝑥12 + 𝑏𝑥1 + 𝑐 = 0, sus raíces son: 𝑥1 = , tomando el
2𝑎
signo positivo, se tiene:

−2 + √22 − 4(3𝑥22 − 10) −2 + 2√1 − (3𝑥22 − 10)


𝑥1 = = = −1 + √11 − 3𝑥22
2 2
Los pasos del algoritmo para resolver mediante el gradiente de ascenso más pronunciado, son:
Paso 1.- 𝑥1 = −1 + √11 − 3𝑥22 de (γ)
1
Paso 2.- 𝜆 = − 1+𝑥 de (α)
1
Paso 3.- ∇ℒ = 6 + 6𝜆𝑥2 de (β)
Paso 4.- 𝑥1𝑖+1 = 𝑥1𝑖 + 𝑐∇ℒ

Considerando 𝑐 = 0.08 y 𝑥20 = 1.5


Iteración 𝑖 = 0
Paso 1.- 𝑥1 = −1 + √11 − 3 ∙ 1.52 = 1.061553
1 1
Paso 2.- 𝜆 = − 1+𝑥 = − 1+1.061553 = −0.485071
1
Paso 3.- ∇ℒ = 6 + 6𝜆𝑥2 = 6 + 6(−0.485071) ∙ 1.5 = 1.634359
Paso 4.- 𝑥21 = 𝑥20 + 𝑐∇ℒ = 1.5 + 0.08 ∙ 1.634359 = 1.630749

En la Tabla siguiente, se tiene las iteraciones:

𝑖 𝑥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

Para un 𝜖 = 0.00001, son necesarios 7 iteraciones. 𝑥1 = 0,658313, 𝑥2 = 1,658312,


𝜆 = −0,603022 y la función objetivo es: 𝑓(𝑥1 , 𝑥2 ) = 2𝑥1 + 6𝑥2 + 10 =21,26649

3.- Sean los datos de los generadores de un sistema: (70%)


Facultad Nacional de Ingeniería Ingeniería Eléctrica e Ingeniería Electrónica

𝐹1 = 250 + 8.5𝑃1 + 0.003𝑃12 35 ≤ 𝑃1 ≤ 300 𝑀𝑊


𝐹2 = 380 + 8.0𝑃2 + 0.0042𝑃22 30 ≤ 𝑃2 ≤ 250 𝑀𝑊
𝐹3 = 245 + 6.1𝑃3 + 0.0040𝑃32 45 ≤ 𝑃3 ≤ 400 𝑀𝑊
Pérdidas:
𝑃𝑝𝑒𝑟𝑑 = 0.00008𝑃12 + 0.00006𝑃22 + 0.00009𝑃32 𝑀𝑊

Para una demanda de 𝑃𝐷 = 200 𝑀𝑊

a) Realizar el despacho de carga por el método iterativo lambda.

La función de Lagrange, es:


3

ℒ (𝑃1 , 𝑃2 , 𝑃3 , 𝜆) = ∑ 𝐹𝑖 + 𝜆(𝑃𝐷 + 𝑃𝑝𝑒𝑟𝑑 − 𝑃1 − 𝑃2 − 𝑃3 )


𝑖=1

Las condiciones de optimalidad, son:

𝜕𝐹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

Las pérdidas marginales, son:


𝜕𝑃𝑝𝑒𝑟𝑑
= 0.00016𝑃1
𝜕𝑃1
𝜕𝑃𝑝𝑒𝑟𝑑
= 0.00012𝑃2
𝜕𝑃2
𝜕𝑃𝑝𝑒𝑟𝑑
= 0.00018𝑃3
𝜕𝑃3

Reemplazando las pérdidas marginales, se tiene:

𝜕𝐹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

Las potencias en función de la 𝜆0 , son:

𝜆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

En la iteración 𝑖 = 1 y con la proyección de Lambda, 𝑖 = 2, la unidades 𝑃1 y 𝑃2 violan las restricciones


de desigualdad, sus cotas mínimas son: 𝑃1𝑀𝑖𝑛 = 35 y 𝑃2𝑀𝑖𝑛 = 30. Significa que son unidades costosas,
pero tienen que trabajar por estar programadas y que entreguen su mínima generación. Por tanto:

𝑃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

𝑃𝑝𝑒𝑟𝑑 = 0.00008𝑃12 + 0.00006𝑃22 + 0.00009𝑃32 𝑀𝑊

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 𝑀𝑊

c) Realizar el despacho de carga por el método del gradiente de segundo orden.

La unidad 3, es la unidad de dependiente 𝑥. No se considera las pérdidas del sistema. Considerando la


siguiente solución factible:
𝑃1 = 50 𝑀𝑊
𝑃2 = 50 𝑀𝑊
𝑃3 = 100 𝑀𝑊
Para 𝑃𝐷 = 200 𝑀𝑊

Los costos marginales, son:


𝑑𝐹1
𝜆1 = = 8.5 + 0.0060𝑃1 = 8.80
𝑑𝑃1
𝑑𝐹2
𝜆2 = = 8.0 + 0.0084𝑃2 = 8.42
𝑑𝑃2
𝑑𝐹3
𝜆3 = = 6.1 + 0.0080𝑃3 = 6.90
𝑑𝑃3

Las segundas derivadas de la función costo, son:


Facultad Nacional de Ingeniería Ingeniería Eléctrica e Ingeniería Electrónica

𝜕 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

0.0060 + 0.0080 0.0080 ∆𝑃 8.80 − 6.90


[ ] [ 1] = − [ ]
0.0080 0.0084 + 0.0080 ∆𝑃2 8.42 − 7.90

0.0140 0.0080 ∆𝑃1 1.90


[ ][ ] = −[ ]
0.0080 0.0164 ∆𝑃2 1.52

∆𝑃1 0.0140 0.0080 −1 1.90


[ ] = −[ ] [ ]
∆𝑃2 0.0080 0.0164 1.52

∆𝑃1 99.033816 −48.309178 1.90


[ ] = −[ ][ ]
∆𝑃2 −48.309178 84.541062 1.52

∆𝑃1 −114.734300
[ ]=[ ]
∆𝑃2 −36.714976

∆𝑃3 = −∆𝑃1 − ∆𝑃2 = 151.449276


Las nuevas potencias:
𝑃1 = 50 − 114.734300 = −64.7343 𝑀𝑊
𝑃2 = 50 − 36.714976 = 13.285924 𝑀𝑊
𝑃3 = 100 + 151.449276 = 251.449276 𝑀𝑊

Las unidades 1 y 2, hay que ajustar a sus generaciones mínimas:

𝑃1 = 35 𝑀𝑊
𝑃2 = 30 𝑀𝑊
𝑃3 = 135 𝑀𝑊

d) Para una demanda de 𝑃𝐷 = 230 𝑀𝑊, aplicar el método de factores de participación

Al incrementar la demanda en ∆𝑃𝐷 = 30 𝑀𝑊 y como las unidades 1 y 2, no pueden aumentar su


generación la única posibilidad es que la unidad 3 aumente su generación en ∆𝑃3 = 30 𝑀𝑊. Las
generaciones son:

𝑃1 = 35 𝑀𝑊
𝑃2 = 30 𝑀𝑊
𝑃3 = 165 𝑀𝑊
Facultad Nacional de Ingeniería Ingeniería Eléctrica e Ingeniería Electrónica

También podría gustarte