EJERCICIOS DEL CAPÍTULO 07

Instituto Especializado de Estudios Superiores Loyola

San Cristóbal, República Dominicana
Elizabet Benzant De la Rosa, 16-03-0024
[email protected]

Realizar los problemas 7.7, 7.8 7.9,7.10 y 7.11

7.7 Las funciones de costo de combustible en $ / h para dos plantas térmicas de 800 MW
están dadas por

C1 = 400 + 6.0P1 + 0.004P21

C2 = 500 + βP2 + γP22

donde P1 y P2 están en MW.

Entonces para comenzar con este ejercicio procedemos a derivas con respecto a la potencia este
será nuestro primer paso:

𝑑𝐶1
= 6 + 0.008P1
𝑑𝑃1

𝑑𝐶2
= β + 2γP2
𝑑𝑃2

(a) El costo incremental de la energía λ es de $ 8 / MWh cuando la demanda de energía total

es de 550 MW. Si se pierden las pérdidas, determine la generación óptima de cada planta.

λ=8
PD = 550MW

λ − βi
𝑃𝑖 =
γi

8−6 2
𝑃1 = = = 𝟐𝟓𝟎𝑴𝑾
0.008 0.008

𝑃2 = PD − 𝑃1 = 550 − 250 = 𝟑𝟎𝟎𝑴𝑾

(b) El costo incremental de la energía λ es de $ 10 / MWh cuando la demanda total de energía
es de 1300 MW. Despreciando las pérdidas, determine la generación óptima de cada planta.

λ = 10
PD = 1300MW

λ − βi
𝑃𝑖 =
γi

10 − 6 4
𝑃1 = = = 𝟓𝟎𝟎𝑴𝑾
0.008 0.008

𝑃2 = PD − 𝑃1 = 1300 − 500 = 𝟖𝟎𝟎𝑴𝑾

(c) De los resultados de (a) y (b) encuentre los coeficientes de costo de combustible β y γ de
la segunda planta.

β + 2γP2 = λ
β + 2γP2 = λ

Remplazamos los valores que tenemos:

𝛽 + 2γ(300) = 8
𝛽 + 2γ(800) = 10
Método de sustitución

𝛽 + 2γ(300) = 8 1
𝛽 = 10 − 2γ(800) 2

Sustituyo a 2 en 1

10 − 2γ(800) + 2γ(300) = 8
10 − 2γ(800 + 300) = 8
−2γ(1100) = 8 − 10
𝛄 = 𝟎. 𝟎𝟎𝟗𝟎𝟗
Sustituyendo a γ en la ecuación 01

𝛽 + 2γ(300) = 8
𝛽 + 2(0.00909)(300) = 8
𝛽 + 5.454 = 8
𝜷 = 𝟐. 𝟓𝟒𝟔
7.8 Las funciones de costo de combustible en $ / h para tres plantas térmicas están dadas
por

C1 = 350 + 7.20P1 + 0.0040P21

C2 = 500 + 7.30P2 + 0.0025P22
C3 = 600 + 6.74P3 + 0.0030P23

donde P1, P2 y P3 están en MW. Los gobernadores están configurados de tal manera que
los generadores comparten la carga por igual. Descuidando las pérdidas de línea y los límites
del generador, encuentre el costo total en $ / h cuando la carga total es

(i) PD = 450 MW

Como nos dice que los generadores comparten la carga por igual entonces asumimos que la PD
se divide entre las 3 cargas:

450
𝑃1 = 𝑃2 = 𝑃3 = = 150𝑀𝑊
3

𝐶𝑡 = 𝐶1 + 𝐶2 + 𝐶3

𝐶𝑡 = 350 + 7.20P1 + 0.0040P12 + 500 + 7.30P2 + 0.0025P22 + 600 + 6.74P3

+ 0.0030P32

𝐶𝑡 = 350 + 7.20(150) + 0.0040(150)2 + 500 + 7.30(150) + 0.0025(150)2 + 600

+ 6.74(150) + 0.0030(150)2

𝐶𝑡 = 350 + 1080 + 90 + 500 + 1095 + 56.25 + 600 + 1011 + 67.50

𝑪𝒕 = 𝟒, 𝟖𝟒𝟗. 𝟕𝟓 $/𝒉

(ii) PD = 745 MW

745
𝑃1 = 𝑃2 = 𝑃3 = = 248.33𝑀𝑊
3

𝐶𝑡 = 𝐶1 + 𝐶2 + 𝐶3

𝐶𝑡 = 350 + 7.20P1 + 0.0040P12 + 500 + 7.30P2 + 0.0025P22 + 600 + 6.74P3

+ 0.0030P32

𝐶𝑡 = 350 + 7.20(248.33) + 0.0040(248.33)2 + 500 + 7.30(248.33)

+ 0.0025(248.33)2 + 600 + 6.74(248.33) + 0.0030(248.33)2
 𝐶𝑡 = 350 + 1787.98 + 246.67 + 500 + 1812.81 + 154.17 + 600 + 1673.74 + 185

𝑪𝒕 = 𝟕, 𝟑𝟏𝟎. 𝟒𝟔 $/𝒉

(iii) PD = 1335 MW

1335
𝑃1 = 𝑃2 = 𝑃3 = = 445𝑀𝑊
3

𝐶𝑡 = 𝐶1 + 𝐶2 + 𝐶3

𝐶𝑡 = 350 + 7.20P1 + 0.0040P12 + 500 + 7.30P2 + 0.0025P22 + 600 + 6.74P3

+ 0.0030P32

𝐶𝑡 = 350 + 7.20(445) + 0.0040(445)2 + 500 + 7.30(445) + 0.0025(445)2 + 600

+ 6.74(445) + 0.0030(445)2

𝐶𝑡 = 350 + 3204 + 792.1 + 500 + 3248.5 + 495.06 + 600 + 2999.3 + 594.08

𝑪𝒕 = 𝟏𝟐, 𝟕𝟖𝟑. 𝟎𝟒 $/𝒉

7.9 Descuidando las pérdidas de línea y los límites del generador, determine la programación
óptima de generación para cada condición de carga en el problema 7.8
(a) Mediante técnica analítica, utilizando (7.33) y (7.31).

𝐏𝐃 = 𝟒𝟓𝟎𝐌𝐖

𝑛𝑔 𝛽𝑖
𝑃𝐷 + ∑𝑖=1 2𝑦𝑖
λ=
1
2𝑦𝑖

7.2 7.3 6.74

450 + 0.008 + + 0.006 450 + 900 + 1460 + 1123.33
λ= 0.005 =
1 1 1 125 + 200 + 16.667
+ +
0.008 0.005 0.006

450 + 3483.33 3933.33

λ= = = 8 $/𝑀𝑊ℎ
491.667 491.667

λ − βi
𝑃𝑖 =
2γi

8 − 7.2 0.8
𝑃1 = = = 100 𝑀𝑊
2(0.004) 0.008
 8 − 7.3 0.7
𝑃2 = = = 140 𝑀𝑊
2(0.0025) 0.005

8 − 6.74 1.26
𝑃3 = = = 210 𝑀𝑊
2(0.003) 0.006

𝐏𝐃 = 𝟕𝟒𝟓𝐌𝐖

𝑛𝑔 𝛽𝑖
𝑃𝐷 + ∑𝑖=1 2𝑦𝑖
λ=
1
2𝑦𝑖

7.2 7.3 6.74

745 + 0.008 + + 0.006 745 + 900 + 1460 + 1123.33
λ= 0.005 =
1 1 1 125 + 200 + 16.667
+ +
0.008 0.005 0.006

745 + 3483.33 4228.33

λ= = = 8.6 $/𝑀𝑊ℎ
491.667 491.667

λ − βi
𝑃𝑖 =
2γi

8.6 − 7.2 1.4

𝑃1 = = = 175 𝑀𝑊
2(0.004) 0.008

8.6 − 7.3 1.3

𝑃2 = = = 260 𝑀𝑊
2(0.0025) 0.005

8.6 − 6.74 1.86

𝑃3 = = = 310 𝑀𝑊
2(0.003) 0.006

𝐏𝐃 = 𝟏𝟑𝟑𝟓𝐌𝐖

𝑛𝑔 𝛽𝑖
𝑃𝐷 + ∑𝑖=1 2𝑦𝑖
λ=
1
2𝑦𝑖

7.2 7.3 6.74

1335 + 0.008 + + 0.006 1335 + 900 + 1460 + 1123.33
λ= 0.005 =
1 1 1 125 + 200 + 16.667
0.008 + 0.005 + 0.006
 1335 + 3483.33 4818.33
λ= = = 9.8 $/𝑀𝑊ℎ
491.667 491.667

λ − βi
𝑃𝑖 =
2γi

9.8 − 7.2 2.6

𝑃1 = = = 325 𝑀𝑊
2(0.004) 0.008

9.8 − 7.3 2.5

𝑃2 = = = 500 𝑀𝑊
2(0.0025) 0.005

9.8 − 6.74 3.06

𝑃3 = = = 510 𝑀𝑊
2(0.003) 0.006

(b) Utilizando el método iterativo. Comience con una estimación inicial de λ = 7.5 $ / MWh.

λ − βi
𝑃𝑖 =
2γi

7.5 − 7.2 0.3

𝑃1(1) = = = 37.5 𝑀𝑊
2(0.004) 0.008

7.5 − 7.3 0.2

𝑃2(1) = = = 40 𝑀𝑊
2(0.0025) 0.005

7.5 − 6.74 0.76

𝑃3(1) = = = 126.67 𝑀𝑊
2(0.003) 0.006

(i) ∆𝑃(1) = 450 − (37.5 + 40 + 126.67) = 245.83 𝑀𝑊

245.83 245.83
∆λ(1) = = = 0.5 $/𝑀𝑊ℎ
1 1 1 491.667
+ +
0.008 0.005 0.006

∆λ(2) = 7.5 + 0.5 = 8 $/𝑀𝑊ℎ

8 − 7.2 0.8
𝑃1(2) = = = 100 𝑀𝑊
2(0.004) 0.008

8 − 7.3 0.7
𝑃2(2) = = = 140 𝑀𝑊
2(0.0025) 0.005
 8 − 6.74 1.26
𝑃3(2) = = = 210 𝑀𝑊
2(0.003) 0.006

∆𝑃(2) = 450 − (100 + 140 + 210) = 0 𝑀𝑊

(ii) ∆𝑃(1) = 745 − (37.5 + 40 + 126.67) = 540.83 𝑀𝑊

540.83 540.83
∆λ(1) = = = 1.1 $/𝑀𝑊ℎ
1 1 1 491.67
0.008 + 0.005 + 0.006

∆λ(2) = 7.5 + 1.1 = 8.6 $/𝑀𝑊ℎ

8.6 − 7.2 1.4

𝑃1(2) = = = 175 𝑀𝑊
2(0.004) 0.008

8.6 − 7.3 1.3

𝑃2(2) = = = 260 𝑀𝑊
2(0.0025) 0.005

8.6 − 6.74 1.86

𝑃3(2) = = = 310 𝑀𝑊
2(0.003) 0.006

∆𝑃(2) = 745 − (175 + 260 + 130) = 0 𝑀𝑊

(iii) ∆𝑃(1) = 1335 − (37.5 + 40 + 126.67) = 1130.83 𝑀𝑊

1130
∆λ(1) = = 2.3 $/𝑀𝑊ℎ
491.67

∆λ(2) = 7.5 + 2.3 = 9.8 $/𝑀𝑊ℎ

9.8 − 7.2 2.6

𝑃1(2) = = = 325 𝑀𝑊
2(0.004) 0.008

9.8 − 7.3 2.5

𝑃2(2) = = = 500 𝑀𝑊
2(0.0025) 0.005

9.8 − 6.74 3.06

𝑃3(2) = = = 510 𝑀𝑊
2(0.003) 0.006

∆𝑃(2) = 745 − (325 + 500 + 510) = 0 𝑀𝑊

(c) Encuentre los ahorros en $ / h para cada caso en comparación con los costos en el
problema 7.8 cuando los generadores compartieron la carga por igual.
 (i) 𝑷𝟏 = 𝟏𝟎𝟎 𝑴𝑾 𝑷𝟐 = 𝟏𝟒𝟎 𝑴𝑾 𝑷𝟑 = 𝟐𝟏𝟎 𝑴𝑾

𝐶𝑡 = 350 + 7.20P1 + 0.0040P12 + 500 + 7.30P2 + 0.0025P22 + 600 + 6.74P3

+ 0.0030P32

𝐶𝑡 = 350 + 7.20(100) + 0.0040(100)2 + 500 + 7.30(140) + 0.0025(140)2 + 600

+ 6.74(210) + 0.0030(210)2

𝐶𝑡 = 350 + 720 + 40 + 500 + 1022 + 49 + 600 + 1415.40 + 132.3

𝑪𝒕 = 𝟒, 𝟖𝟐𝟖. 𝟕𝟎 $/𝒉

Comparación con el problema 7.8:

𝐶𝑡 = 4,849.75 $/ℎ − 4,828.70 $/ℎ

𝑪𝒕 = 𝟐𝟏. 𝟎𝟓 $/𝒉

(ii) 𝐏𝟏 = 𝟏𝟕𝟓 𝐌𝐖 𝐏𝟐 = 𝟐𝟔𝟎 𝐌𝐖 𝐏𝟑 = 𝟑𝟏𝟎 𝐌𝐖

𝐶𝑡 = 350 + 7.20P1 + 0.0040P12 + 500 + 7.30P2 + 0.0025P22 + 600 + 6.74P3

+ 0.0030P32

𝐶𝑡 = 350 + 7.20(175) + 0.0040(175)2 + 500 + 7.30(260) + 0.0025(260)2 + 600

+ 6.74(310) + 0.0030(310)2

𝐶𝑡 = 350 + 1260 + 122.5 + 500 + 1898 + 169 + 600 + 2089.4 + 288.3

𝑪𝒕 = 𝟕, 𝟐𝟕𝟕. 𝟐𝟎 $/𝒉

Comparación con el problema 7.8:

𝐶𝑡 = 7,783.04 $/ℎ − 7,277.20 $/ℎ

𝑪𝒕 = 𝟑𝟑. 𝟐𝟔 $/𝒉

(iii) 𝐏𝟏 = 𝟑𝟐𝟓 𝐌𝐖 𝐏𝟐 = 𝟓𝟎𝟎 𝐌𝐖 𝐏𝟑 = 𝟓𝟏𝟎 𝐌𝐖

𝐶𝑡 = 350 + 7.20P1 + 0.0040P12 + 500 + 7.30P2 + 0.0025P22 + 600 + 6.74P3

+ 0.0030P32
 𝐶𝑡 = 350 + 7.20(325) + 0.0040(325)2 + 500 + 7.30(500 ) + 0.0025(500 )2 + 600
+ 6.74(510) + 0.0030(510)2

𝐶𝑡 = 350 + 2340 + 422.5 + 500 + 3650 + 625 + 600 + 3437.4 + 780.3

𝑪𝒕 = 𝟏𝟐, 𝟕𝟎𝟓. 𝟐𝟎 $/𝒉

Comparación con el problema 7.8:

𝐶𝑡 = 12,783.04 $/ℎ − 12,705.20 $/ℎ

𝑪𝒕 = 𝟕𝟕. 𝟖𝟒 $/𝒉

Como podemos observar la diferencia de costo que existe entre el ejercicio 7.8 y 7.9 es muy
mínima.

7.10. Repita el problema 7.9 (a) y (b), pero esta vez considere los siguientes límites del
generador (en MW)

122 ≤ P1 ≤400
260 ≤ P2 ≤600
50 ≤ P3 ≤445

Como el problema 7.9 (a) y (b) obtuvimos varias potencias ahora vamos a evaluar dentro de los
limites que nos acaban de poner. En la parte (a) (i) teníamos que la:

P1=100MW P2=140 MW P3= 210MW

Entonces podemos observar que la potencia P1 y P2 están por debajo de los límites, usaremos
para P1 y P2 los limites en este caso inferiores y buscaremos a P3. Pues entonces tendremos:

P3= PD-(P1+P2) = 450-(122+260) = 68MW

Seguimos evaluando en la parte (a) (ii) teníamos entonces:

P1=175MW P2=260 MW P3= 310MW

Entonces podemos observar que todas las potencias entran dentro de los limites establecidos. Ya
en la ultima parte para conseguir nuestro PD que es lo que andamos buscando en la parte (a) (iii)
teníamos entonces:

P1=325 MW P2=500 MW P3= 510 MW

Entonces podemos observar que en este caso solo la potencia P3 pasa de los límites establecidos,
entonces usaremos este caso el límite superior P3= 445 MW entonces tendríamos:
 PD= 1335-445= 890 MW

7.2 7.3
890 + 0.008 +
λ= 0.005
1 1
0.008 + 0.005
890 + 2360
λ=
325

𝛌 = 𝟏𝟎 $/𝑴𝑾𝒉

10 − 7.2
𝑃1 = = 350 𝑀𝑊
2(0.004)

10 − 7.3
𝑃2 = = 540 𝑀𝑊
2(0.0025)
Entonces para la parte (b) tenemos:

P1=325 MW P2=500 MW P3= 510 MW

Entonces podemos observar que en este caso solo la potencia P3 pasa de los límites establecidos,
entonces usaremos este caso el límite superior P3= 445 MW.

∆𝑃(2) = 1335 − (325 + 500 + 445) = 65 𝑀𝑊

65 65
λ(2) = =
1 1 325
+
2(0.004) 2(0.0025)

𝛌(𝟐) = 𝟎. 𝟐 $/𝑴𝑾𝒉

λ(3) = 0.2 + 9.8 = 10 $/𝑀𝑊ℎ

10 − 7.2 2.8
𝑃1(2) = = = 350 𝑀𝑊
2(0.004) 0.008

10 − 7.3 2.7
𝑃2(2) = = = 540 𝑀𝑊
2(0.0025) 0.005

𝑃3(2) = 445 𝑀𝑊

∆𝑃(2) = 1335 − (350 + 540 + 445) = 0 𝑀𝑊

 𝐶𝑡 = 350 + 7.20P1 + 0.0040P12 + 500 + 7.30P2 + 0.0025P22 + 600 + 6.74P3
+ 0.0030P32

𝐶𝑡 = 350 + 7.20(325) + 0.0040(325)2 + 500 + 7.30(500 ) + 0.0025(500 )2 + 600

+ 6.74(445) + 0.0030(445)2

𝐶𝑡 = 350 + 2340 + 422.5 + 500 + 3650 + 625 + 600 + 2999.3 + 594.075

𝑪𝒕 = 𝟏𝟐, 𝟕𝟐𝟒. 𝟑𝟖 $/𝒉

7.11. La función de costo de combustible en $ /h de dos plantas térmicas son:

C1 = 320 + 6.2P1 + 0.004P21

C2 = 200 + 6.0P2 + 0.003P22

donde P1 y P2 están en MW. Plantar fuera de lugar sujeto a los límites de guía (inMW)

50≤P1≤250
50≤P2≤350

La pérdida de potencia real del sistema por unidades con la generación expresada en la
unidad en una base de 100 MVA se da por:

PL (pu)= 0.0125P21(pu)+ 0.00625P22(pu)

La carga total es de 412.35 MW. Determinar el despacho óptimo de generación. Comience

con una estimación inicial de λ = 7 $ / MWh.

𝑃1 2 𝑃2 2
PL = [0.0125 ( ) + 0.00625 ( ) ] × 100𝑀𝑊
100 100

PL = 0.000125𝑃12 + 0.0000625𝑃22 𝑀𝑊

7 − 6.2 0.8
𝑃1(1) = = = 82.0513 𝑀𝑊
2(0.004 + 7 × 0.000125) 0.00975

7−6 1
𝑃2(1) = = = 145.4545 𝑀𝑊
2(0.003 + 7 × 0.0000625) 0.006875

PL = [0.000125(82.0513 )2 + 0.0000625(145.4545)2 ] = 0.84155 + 1.3223

PL = 2.1639
 ∆𝑃1(1) = 412.35 + 2.1639 − (82.0513 + 145.4545) = 187.01
2
𝜕𝑃1 2 0.004 + 0.000125 × 6.2 0.003 + 0.000625 × 6
∑( ) = 2
+
𝜕λ 2(0.004 + 7 × 0.000125) 2(0.003 + 7 × 0.00625)2
𝑖=1

2
𝜕𝑃1 2 0.004775 0.00675
∑( ) = +
𝜕λ 0.000004753 0.004371
𝑖=1

2
𝜕𝑃1 2
∑( ) = 243.27
𝜕λ
𝑖=1

187.01
∆λ(1) = = 0.7687 $/𝑀𝑊ℎ
243.27

λ(2) = 0.7687 + 7 = 7.769 $/𝑀𝑊ℎ

7.769 − 6.2
𝑃1(2) = = 157.785 𝑀𝑊
2(0.004 + 7.769 × 0.000125)

7.769 − 6
𝑃2(2) = = 253.723 𝑀𝑊
2(0.003 + 7.769 × 0.000625)

PL(2) = [0.0125(157.785 )2 + 0.00625(253.723)2 ] = 7.136

∆𝑃1(2) = 412.35 + 7.136 − (157.785 + 253.723) = 7.978

2
𝜕𝑃1 2 0.004 + 0.000125 × 6.2 0.003 + 0.000625 × 6
∑( ) = 2
+
𝜕λ 2(0.004 + 7.769 × 0.000125) 2(0.003 + 7.769 × 0.000625)2
𝑖=1

2
𝜕𝑃1 2 0.004775 0.00675
∑( ) = +
𝜕λ 0.0000049424 0.001234
𝑖=1

2
𝜕𝑃1 2
∑( ) = 235.514
𝜕λ
𝑖=1

7.978
∆λ(2) = = 0.0339 $/𝑀𝑊ℎ
235.514
 λ(2) = 0.0339 + 7.769 = 7.803 $/𝑀𝑊ℎ

7.803 − 6.2
𝑃1(3) = = 161.055 𝑀𝑊
2(0.004 + 7.803 × 0.000125)

7.803 − 6
𝑃2(3) = = 258.425 𝑀𝑊
2(0.003 + 7.803 × 0.000625)

PL(3) = [0.0125(161.055 )2 + 0.00625(258.425)2 ] = 7.4163

∆𝑃1(3) = 412.35 + 7.4163 − (161.055 + 258.425 ) = 0.2863

2
𝜕𝑃1 2 0.004 + 0.000125 × 6.2 0.003 + 0.000625 × 6
∑( ) = 2
+
𝜕λ 2(0.004 + 7.803 × 0.000125) 2(0.003 + 7.803 × 0.000625)2
𝑖=1

2
𝜕𝑃1 2
∑( ) = 235.18
𝜕λ
𝑖=1

0.2863
∆λ(3) = = 0.0012 $/𝑀𝑊ℎ
235.18

λ(3) = 0.0012 + 7.803 = 7.804 $/𝑀𝑊ℎ

𝐶𝑡 = 320 + 6.2P1 + 0.004P12 + 200 + 6.0P2 + 0.003P22

𝐶𝑡 = 320 + 6.2(161.055) + 0.004(161.055)2 + 200 + 6 (258.425) + 0.003(258.425)2

𝐶𝑡 = 320 + 998.541 + 103.7549 + 200 + 1550.55 + 200.35

𝑪𝒕 = 𝟑, 𝟑𝟕𝟑. 𝟏𝟗𝟔 $/𝒉