“Ciencia y Tecnología al Servicio del País”

25-I

Facultad de Ingeniería Eléctrica y

Electrónica
Curso: Análisis de Sistemas de Potencia II

Docente: Mg. Ing. John Negreiros Terrones

2025-I
Despacho Económico
 Temario

• Despacho Económico de Unidades Térmicas y Métodos de

Solución

[Link] problema del despacho económico.

[Link] económico con funciones de coste lineales por partes
[Link]ón dinámica.
 Despacho económico

3.1 EL PROBLEMA DEL DESPACHO ECONÓMICO

𝐹𝑇 = 𝐹1 + 𝐹2 + 𝐹3 + ⋯ + 𝐹𝑁𝑔𝑒𝑛 (Ec 1)

𝑁𝑔𝑒𝑛

𝐹𝑇 = ෍ 𝐹𝑖 𝑃𝑖 (Ec 2)
𝑖=1

𝑁𝑔𝑒𝑛

∅ = 0 = 𝑃𝑙𝑜𝑎𝑑 − ෍ 𝑃𝑖 (Ec 3)
𝑖=1

Consumo de combustible: Mbtu/h

Tasa de costo: $/h

 Despacho económico

Para establecer las condiciones necesarias para un valor extremo de la función objetivo, agregue la función de
restricción a la función objetivo después de que la función de restricción haya sido multiplicada por un
multiplicador indeterminado
ℒ = 𝐹𝑇 + 𝜆 ∅ (Ec 4)
𝜕ℒ 𝑑𝐹𝑖 𝑃𝑖 𝑑𝐹𝑖
= −𝜆 o 0= −𝜆 (Ec 5)
𝜕𝑃𝑖 𝑑𝑃𝑖 𝑑𝑃𝑖
Además, hay dos desigualdades que deben satisfacerse para cada una de las unidades. Es decir, la potencia de
salida de cada unidad debe ser superior o igual a la potencia mínima permitida y también debe ser inferior o
igual a la potencia máxima permitida en esa unidad en particular. (Ec 6)
𝑑𝐹𝑖
−𝜆 𝑁𝑔𝑒𝑛 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛𝑠
𝑑𝑃𝑖
𝑃𝑖,𝑚𝑖𝑛 ≤ 𝑃𝑖 ≤ 𝑃𝑖,𝑚𝑎𝑥 2𝑁𝑔𝑒𝑛 𝑖𝑛𝑒𝑞𝑢𝑎𝑡𝑖𝑜𝑛𝑠
𝑁

෍ 𝑃𝑖 = 𝑃𝑙𝑜𝑎𝑑 1 𝑐𝑜𝑛𝑠𝑡𝑟𝑎𝑖𝑛𝑡
𝑖=1
 Despacho económico

Cuando reconocemos las restricciones de desigualdad, entonces las condiciones necesarias pueden
expandirse ligeramente como se muestra en el conjunto de ecuaciones que componen la Ecuación 7

𝑑𝐹𝑖
=𝜆 𝑓𝑜𝑟 𝑃𝑖,𝑚𝑖𝑛 ≤ 𝑃𝑖 ≤ 𝑃𝑖,𝑚𝑎𝑥
𝑑𝑃𝑖

𝑑𝐹𝑖
≤𝜆 𝑓𝑜𝑟 𝑃𝑖 = 𝑃𝑖,𝑚𝑎𝑥
𝑑𝑃𝑖

𝑑𝐹𝑖
≥𝜆 𝑓𝑜𝑟 𝑃𝑖 = 𝑃𝑖,𝑚𝑖𝑛
𝑑𝑃𝑖
 Despacho económico

Ejemplo
Unidad 1: Unidad de vapor : Unidad 2: Unidad de vapor :
Salida máxima = 600 MW Salida máxima = 400 MW
Salida mínima = 150 MW Salida mínima = 100 MW
Curva de entrada-salida: Curva de entrada-salida:
𝑀𝐵𝑡𝑢 𝑀𝐵𝑡𝑢
𝐻1 = 510.0 + 7.2𝑃1 + 0.00142𝑃12 𝐻2 = 310.0 + 7.85𝑃2 + 0.00194𝑃22
ℎ ℎ

Unidad 3: Unidad de vapor:

Salida máxima = 200 MW
Salida mínima = 50 MW
Curva de entrada-salida:
𝑀𝐵𝑡𝑢
𝐻3 = 78.0 + 7.97𝑃3 + 0.00482𝑃32
ℎ
 Despacho económico

Supongamos que deseamos determinar el punto económico de operación de estas tres

unidades al entregar un total de 850 MW. Antes de que se pueda resolver este problema,
se debe especificar el costo de combustible de cada unidad. Deje que los siguientes costos
de combustible estén en vigor.

Unit1 : fuel cost =1.1 $ /MBtu

Unit2 : fuel cost =1.0 $ /MBtu
Unit3 : fuel cost =1.0 $ /Mbtu

then
𝐹1 𝑃1 = 𝐻1 𝑃1 𝑥1.1 = 561 + 7.92𝑃1 + 0.001562𝑃12

𝐹2 𝑃2 = 𝐻2 𝑃2 𝑥1.0 = 310.0 + 7.85𝑃2 + 0.00194𝑃22

𝐹3 𝑃3 = 𝐻3 𝑃3 𝑥1.0 = 78.0 + 7.97𝑃3 + 0.00482𝑃32

 Despacho económico

Usando la Ecuación 5, las condiciones para un despacho óptimo son

𝑑𝐹1
= 7.92 + 0.003124𝑃1 = 𝜆
𝑑𝑃1

𝑑𝐹2
= 7.85 + 0.00388𝑃2 = 𝜆 𝑃1 + 𝑃2 + 𝑃3 = 850 𝑀𝑊
𝑑𝑃2

𝑑𝐹3
= 7.97 + 0.00964𝑃3 = 𝜆
𝑑𝑃3

𝜆 = 9.148 $Τ𝑀𝑊ℎ

𝑃1 = 704.6 𝑀𝑊 𝑃2 = 111.8 𝑀𝑊 𝑃3 = 32.6 𝑀𝑊

Esta solución cumple con la restricción que requiere que la generación total sea igual a 850 MW,
pero las unidades 1 y 3 no están dentro del límite.
 Despacho económico

Usando la ecuación 6
Suponga que la unidad 1 se establece en su salida máxima y la unidad 3 en su salida
mínima. El envío se convierte
𝑃1 = 600 𝑀𝑊

𝑃2 = 200 𝑀𝑊

𝑃3 = 50 𝑀𝑊

De la Ecuación 6, vemos que λ debe ser igual al costo incremental de la unidad 2 ya

que no está en ninguno de los dos límites. Entonces.

𝑑𝐹2 𝑑𝐹1 𝑑𝐹3

𝜆= ቤ = 8.626 $Τ𝑀𝑊ℎ ቤ = 8.016 $Τ𝑀𝑊ℎ ቤ = 8.452 $Τ𝑀𝑊ℎ
𝑑𝑃2 𝑃 𝑑𝑃1 𝑃 𝑑𝑃3 𝑃
2 =200 1 =600 3 =50
 Despacho económico

Tenga en cuenta que el costo incremental para la unidad 1 es menor que λ, por lo que la
unidad 1 debería estar en su máximo. Sin embargo, el costo incremental de la unidad 3 no
es mayor que λ, por lo que la unidad 3 no debe forzarse a su mínimo. Por lo tanto, para
encontrar el despacho óptimo, permita que el costo incremental en las unidades 2 y 3 sea
igual a λ de la siguiente manera:
𝑃1 = 600𝑀𝑊

𝑑𝐹2 𝑑𝐹3
= 7.85 + 0.00388𝑃2 = 𝜆 = 7.97 + 0.00964𝑃3 = 𝜆
𝑑𝑃2 𝑑𝑃3

entonces:
𝑃2 = 187.1𝑀𝑊
𝜆 = 8.576 $Τ𝑀𝑊ℎ
𝑃3 = 62.9𝑀𝑊
 Despacho económico

Nótese que este despacho cumple con las condiciones de la Ecuación 3.6 ya que

𝑑𝐹1
ቤ = 8.016 $Τ𝑀𝑊ℎ
𝑑𝑃1 𝑃
1 =600𝑀𝑊

que es menor que λ, mientras que (dF2/dP2) y (dF3/dP3) son iguales a λ.

 Despacho económico

3.2 Despacho económico con funciones de coste lineales por partes

Para todas las unidades en funcionamiento, comenzamos con todas ellas en Pmin, luego comenzamos a aumentar la
producción de la unidad con el segmento de costo incremental más bajo. Si esta unidad llega al extremo derecho de un
segmento, o si llega a Pmax, entonces encontramos la unidad con el siguiente segmento de costo incremental más bajo y
aumentamos su producción.
 Despacho económico

3.3.1 Funciones de costo lineales por tramos

LP es muy hábil en el manejo de restricciones de desigualdad siempre que el problema a resolver sea tal que
pueda linealizarse sin pérdida de precisión.
En la formulación que sigue, mostramos cómo el problema del despacho económico puede estructurarse
como un PL. Primero, abordamos el problema de expresar las funciones no lineales de insumo-producto o
costo como un conjunto de funciones lineales.
 Despacho económico

Comenzamos con una función de costo no lineal que se muestra en la figura 3.3. Podemos aproximar esta
función no lineal como una serie de segmentos de línea recta como se muestra en la Figura 3.4.

𝑚𝑖𝑛 + 𝑠 𝑃
𝐹𝑖 𝑃𝑔𝑒𝑛𝑖 = 𝐹𝑖 𝑃𝑔𝑒𝑛𝑖 𝑖1 𝑔𝑒𝑛𝑖1 + 𝑠𝑖2 𝑃𝑔𝑒𝑛𝑖2 + 𝑠𝑖3 𝑃𝑔𝑒𝑛𝑖3

𝑚𝑖𝑛
0 ≤ 𝑃𝑔𝑒𝑛𝑖𝑘 ≤ 𝑃𝑔𝑒𝑛𝑖𝑘
𝑓𝑜𝑟 𝑘 = 1,2,3

𝑚𝑖𝑛 + 𝑃 + 𝑃
𝑃𝑔𝑒𝑛𝑖 = 𝑃𝑔𝑒𝑛𝑖 𝑖 𝑔𝑒𝑛𝑖1 + 𝑃𝑔𝑒𝑛𝑖2 + 𝑃𝑔𝑒𝑛𝑖3

𝑚𝑖𝑛 𝑚𝑖𝑛
𝐹𝑖 𝑃𝑔𝑒𝑛𝑖𝑘+1
− 𝐹𝑖 𝑃𝑔𝑒𝑛𝑖𝑘
𝑆𝑖𝑘 = 𝑚𝑖𝑛 𝑚𝑖𝑛
𝑃𝑔𝑒𝑛𝑖𝑘+1
− 𝑃𝑔𝑒𝑛𝑖𝑘
 Despacho económico

3.3.2 Despacho Económico con LP

Ahora la solución de PL del despacho económico se puede escribir como

Minimizar: 𝑁𝑔𝑒𝑛
𝑚𝑖𝑛 + 𝑠 𝑃
෍ 𝐹𝑖 𝑃𝑔𝑒𝑛𝑖 𝑖1 𝑔𝑒𝑛𝑖1 + 𝑠𝑖2 𝑃𝑔𝑒𝑛𝑖2 + 𝑠𝑖3 𝑃𝑔𝑒𝑛𝑖3
𝑖=1

Sujeto a:
𝑚𝑖𝑛 𝑓𝑜𝑟 𝑘 = 1,2,3 … 𝑓𝑜𝑟 𝑎𝑙𝑙 𝑔𝑒𝑛𝑒𝑟𝑎𝑡𝑜𝑟 𝑖 = 1 … 𝑁
0 ≤ 𝑃𝑔𝑒𝑛𝑖𝑘 ≤ 𝑃𝑔𝑒𝑛𝑖𝑘 𝑔𝑒𝑛

𝑃𝑖 = 𝑃𝑖𝑚𝑖𝑛 + 𝑃𝑔𝑒𝑛𝑖1 + 𝑃𝑔𝑒𝑛𝑖2 + 𝑃𝑔𝑒𝑛𝑖3 𝑓𝑜𝑟 𝑎𝑙𝑙 𝑔𝑒𝑛𝑒𝑟𝑎𝑡𝑜𝑟 𝑖 = 1 … 𝑁𝑔𝑒𝑛

𝑁𝑔𝑒𝑛

෍ 𝑃𝑖 = 𝑃𝑙𝑜𝑎𝑑
𝑖=1
 Despacho económico

Ahora mostraremos cómo LP resuelve el mismo despacho económico de tres generadores

dado en el Ejemplo 3A. Antes de investigar las soluciones que da el LP, es interesante ver
qué variables tratará el PL. Aquí hay una lista de todas las variables para los tres
generadores, problemas con tres segmentos usados para cada generador:

Variables del generador:

𝑃𝑔𝑒𝑛1 , 𝑃𝑔𝑒𝑛11 , 𝑃𝑔𝑒𝑛12 , 𝑃𝑔𝑒𝑛13 , 𝑃𝑔𝑒𝑛2 , 𝑃𝑔𝑒𝑛21 , 𝑃𝑔𝑒𝑛22 , 𝑃𝑔𝑒𝑛23 , 𝑃𝑔𝑒𝑛3 , 𝑃𝑔𝑒𝑛31 , 𝑃𝑔𝑒𝑛32 , 𝑃𝑔𝑒𝑛33
 Despacho económico

Ejemplo 3C:
Solución usando LP

2
𝐹1 𝑃𝑔𝑒𝑛1 = 561 + 7.92𝑃𝑔𝑒𝑛1 + 0.001562𝑃𝑔𝑒𝑛1

2
𝐹2 𝑃𝑔𝑒𝑛2 = 310.0 + 7.85𝑃𝑔𝑒𝑛2 + 0.00194𝑃𝑔𝑒𝑛2

2
𝐹3 𝑃𝑔𝑒𝑛3 = 78.0 + 7.97𝑃𝑔𝑒𝑛3 + 0.00482𝑃𝑔𝑒𝑛3

con límites mínimos y máximos de potencia del generador de

𝑚 𝑖𝑛 𝑚 𝑎𝑥
Generator 𝑃𝑔𝑒𝑛 𝑃𝑔𝑒𝑛
1 150 600
2 100 400
3 50 200
 Despacho económico

y con una carga total de 850 MW, podremos calcular el flujo de potencia óptimo (OPF)
usando LP. La única diferencia es que dependiendo del número de segmentos utilizados, la
solución será diferente a la obtenida por el método estándar.
Con K=1 segmento: para cada generador 𝑃𝑔𝑒𝑛1 , 𝑃𝑔𝑒𝑛11 , 𝑃𝑔𝑒𝑛2 , 𝑃𝑔𝑒𝑛21 , , 𝑃𝑔𝑒𝑛3 , 𝑃𝑔𝑒𝑛31
Usando el costo lineal por partes tenemos
𝑚𝑖𝑛
𝑃𝑔𝑒𝑛1 = 𝑃𝑔𝑒𝑛1 + 𝑃𝑔𝑒𝑛11 = 𝑃𝑔𝑒𝑛1 =400

𝑚𝑖𝑛
𝑃𝑔𝑒𝑛2 = 𝑃𝑔𝑒𝑛2 + 𝑃𝑔𝑒𝑛21 = 𝑃𝑔𝑒𝑛2 = 400
𝑚𝑖𝑛
𝑃𝑔𝑒𝑛3 = 𝑃𝑔𝑒𝑛3 + 𝑃𝑔𝑒𝑛13 = 𝑃𝑔𝑒𝑛3 = 50

෍ 𝐹𝑖 𝑃𝑔𝑒𝑛𝑖
2 2 2
= 561 + 7.92 400 + 0.001562 400 + 310 + 7.85 400 + 0.00194 400 + 78.0 + 7.97 50 + 0.00482 50

෍ 𝐹𝑖 𝑃𝑔𝑒𝑛𝑖 = 3978.92 + 3760.4 + 488.55 = 8227.87

 Despacho económico

En la siguiente tabla, usamos 1, 3, 5, 10 y 50 segmentos para cada función de costo, y se

puede observar que la solución se acerca a la misma solución a medida que aumenta el
número de segmentos.

Number of Generator 1 Generator 2 Generator 3 Total Cost

Segments MW MW MW ($/h)
1 400 400 50 8227.870
2 375 350 125 8195.369
3 450 300 100 8204.105
5 400 340 110 8195.206
10 385 340 125 8194.554
50 393 335 122 8194.357
Standard 393.2 334.6 122.2 8194.356
solution with
lambda
search
 Despacho económico

3.4 EL MÉTODO DE ITERACIÓN LAMBDA

La figura 3.5 es un diagrama de bloques del

método de solución de iteración lambda para las
pérdidas totalmente térmicas que ignoran el
problema de despacho. Podemos acercarnos a la
solución de este problema considerando una
técnica gráfica para resolver el problema y luego
extendiéndola al área de los algoritmos
informáticos.
 Despacho económico

EL MÉTODO DE ITERACIÓN LAMBDA

Supongamos que tenemos un sistema de tres

máquinas y deseamos encontrar el punto
óptimo de funcionamiento económico. Un
enfoque sería trazar las características de costo
incremental para cada una de estas tres
unidades en el mismo gráfico, como se muestra
en la Figura 3.6.
 Despacho económico

EL MÉTODO DE ITERACIÓN LAMBDA

Por supuesto, nuestra primera estimación será
incorrecta. Si hemos asumido el valor del
costo incremental tal que la producción total
de energía es demasiado baja, debemos
aumentar el valor de l y probar otra solución.
Con dos soluciones, podemos extrapolar (o
interpolar) las dos soluciones para acercarnos
al valor deseado de potencia total recibida
(ver Figura 3.7).
 Despacho económico

Ejemplo 3D: suponga que desea utilizar funciones cúbicas para representar las características de
entrada y salida de las plantas generadoras de la siguiente manera
𝐻 𝑀𝐵𝑡𝑢 Τℎ = 𝐴 + 𝐵𝑃 + 𝐶𝑃2 + 𝐷𝑃3 𝑃 𝑖𝑛 𝑀𝑊

Para las tres unidades, encuentre el programa óptimo utilizando el método de iteración lambda.
A B C D
Unit 1 749.55 6.95 9.68𝑥10−4 1.27𝑥10−7
Unit 2 1285.0 7.051 7.375𝑥10−4 6.543𝑥10−8
Unit 3 1531.0 6.531 1.04𝑥10−3 9.98𝑥10−8

Suponga que el costo del combustible es de 1,0 $/MBtu por cada unidad y los límites de la unidad son
los siguientes:
320 𝑀𝑊 ≤ 𝑃1 ≤ 800 𝑀𝑊

300 𝑀𝑊 ≤ 𝑃2 ≤ 1200 𝑀𝑊

275 𝑀𝑊 ≤ 𝑃3 ≤ 1100 𝑀𝑊
 Despacho económico

Se muestran dos ejemplos de cálculos, ambos utilizando el diagrama de flujo de la Figura 3.5. En este
cálculo, el valor de l en la segunda iteración siempre se establece en un 10 % por encima o por
debajo del valor inicial según el signo del error; para las iteraciones restantes, lambda se proyecta
como en la Figura 3.7.
El primer ejemplo muestra la ventaja de comenzar 𝜆 cerca del valor óptimo.

𝑃𝑙𝑜𝑎𝑑 = 2500 𝑀𝑊
𝜆𝑠𝑡𝑎𝑟𝑡 = 8.00 $Τ𝑀 𝑊ℎ
Iteration 𝜆 Total Generation 𝑃1 𝑃2 𝑃3
(MW)

1 8.0000 1731.6 494.3 596.7 640.6

2 8.8000 2795.0 800.0 1043.0 952.0
3 8.5781 2526.0 734.7 923.4 867.9
4 8.5566 2497.5 726.1 911.7 859.7
5 8.5586 2500.0 726.9 912.7 860.4
 Despacho económico

El segundo ejemplo muestra los problemas oscilatorios que se pueden encontrar con un enfoque de
iteración lambda..

𝑃𝑙𝑜𝑎𝑑 = 2500 𝑀𝑊
𝜆𝑠𝑡𝑎𝑟𝑡 = 10.00 $Τ𝑀 𝑊ℎ
 Despacho económico

3.5 DESPACHO ECONÓMICO MEDIANTE BÚSQUEDA BINARIA

Un algoritmo de iteración lambda muy útil, a veces denominado “búsqueda binaria”2, evita las
oscilaciones y siempre logra encontrar el despacho económico óptimo. Primero, el usuario debe
calcular el costo incremental en la producción máxima y mínima de generación para cada
generador.
Luego establezca 𝜆𝑚𝑖𝑛 en el valor más pequeño entre los costos incrementales en los valores
Pmin del generador y luego establezca 𝜆𝑚𝑎𝑥 como el valor más grande entre los costos
incrementales en los valores Pmax del generador. Si λ = λ_min, entonces el algoritmo de
búsqueda lambda establecería todos los generadores en P = Pmin, y si λ = λ_max, entonces el
algoritmo de búsqueda lambda establecería todos los generadores en P = Pmax.

𝜆𝑚𝑎𝑥 − 𝜆𝑚𝑖𝑛
∆𝜆 =
2
𝜆𝑖 = 𝜆𝑚𝑖𝑛 + ∆𝜆

(2 El nombre de “búsqueda binaria” y el desarrollo del método de despacho económico se atribuyen a Walter O. Stadlin.)
 Despacho económico

Ahora calcule la salida del generador que corresponde a cada generador que tiene este costo
incremental. Si el valor 𝜆 1 es menor que el costo incremental en Pmin, simplemente configure la
salida del generador en Pmin, y si el valor 𝜆 1 es mayor que el costo incremental en Pmax,
simplemente configure la salida del generador en Pmax; de lo contrario, calcule el valor P para el
generador a partir de la función de costo incremental. Ahora agregue todas las salidas del
generador:

𝑁𝑔𝑒𝑛
si
෍ 𝑃𝑖 > 𝑃𝑙𝑜𝑎𝑑
𝑖=1

debemos reducir lambda, entonces

∆𝜆
∆𝜆 =
2

𝜆𝑖+1 = 𝜆𝑖 − ∆𝜆
 Despacho económico

si
𝑁𝑔𝑒𝑛

෍ 𝑃𝑖 < 𝑃𝑙𝑜𝑎𝑑
𝑖=1

debemos aumentar lambda, entonces

∆𝜆
∆𝜆 =
2

𝜆𝑖+1 = 𝜆𝑖 + ∆𝜆

𝑁𝑔𝑒𝑛
Si
𝑎𝑏𝑠 ෍ 𝑃𝑖 − 𝑃𝑙𝑜𝑎𝑑 < 𝑡𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑒, 𝑎𝑙𝑔𝑜𝑟𝑖𝑡ℎ𝑚 𝑖𝑠 𝑑𝑜𝑛𝑒
𝑖=1

El valor delta lambda cambia a la mitad en cada iteración. La búsqueda binaria funcionará con
funciones de costo incremental lineal por partes e incluso funciones de costo incremental lineal
por partes que tienen saltos entre segmentos lineales.
 Despacho económico

3.6 DESPACHO ECONÓMICO MEDIANTE PROGRAMACIÓN

DINÁMICA

La solución DP al despacho económico se realiza como un problema de asignación. Con este

enfoque, no calculamos un único conjunto óptimo de salidas de MW del generador para una carga
total específica suministrada; más bien, generamos un conjunto de salidas, en puntos discretos,
para un conjunto completo de valores de carga.
Ejemplo 3E: suponga que hay tres unidades en el sistema; todos están en línea. Sus características
de entrada y salida no son suaves ni convexas. Los datos son los siguientes.
 Despacho económico
La demanda total es D = 310 MW. Esto no se ajusta exactamente a los datos, por lo que debemos
interpolar entre los valores más cercanos que están disponibles a partir de los datos, 300 y 325
MW.
Programando las unidades 1 y 2, encontramos el costo mínimo para la función
𝑓2 = 𝐹1 𝐷 − 𝑃2 + 𝐹2 𝑃2
sobre el rango permisible de P2 y para 100 ≤ D ≤ 350 MW. Los datos de búsqueda se dan en la
siguiente tabla. Necesitamos ahorrar el costo de servir cada valor de D que sea mínimo y el nivel
de carga en la unidad 2 para cada nivel de demanda
 Despacho económico
Nuevamente minimizamos
𝑓3 = 𝐹2 𝐷 − 𝑃3 + 𝐹3 𝑃3

para 50 ≤ P3 ≤ 175 MW y D = 300 y 325 MW. Programar la tercera unidad para los dos niveles de
demanda diferentes solo requiere dos filas de la siguiente tabla.

de modo que entre los niveles de demanda de 300 y 325 MW, la unidad marginal es la unidad 2.
(Es decir, está absorbiendo todo el aumento de demanda adicional entre 300 y 325 MW). Por lo
tanto, podemos interpolar para encontrar el costo a un nivel de carga de 310 MW o un nivel de
salida en la unidad 2 de 110 MW. Los resultados para un nivel de demanda de 310 MW son: P 1=
50, P2 = 110, and P3 = 150 for a total cost of 4286 $ / h
 Despacho económico
Dada una carga que debe suministrarse en incrementos de tiempo t = 1 … tmax, con niveles de
carga de Pt y generadores Ngen en línea para suministrar la carga,
𝑁𝑔𝑒𝑛

෍ 𝑃𝑖𝑡 = 𝑃𝑙𝑜𝑎𝑑
𝑡

𝑖=1
Cada unidad debe obedecer a un límite de velocidad tal que
𝑃𝑖𝑡+1 = 𝑃𝑖𝑡 + ∆𝑃𝑖

−∆𝑃𝑖𝑚𝑖𝑛 ≤ ∆𝑃𝑖 ≤ ∆𝑃𝑖𝑚𝑎𝑥

Entonces debemos programar las unidades para minimizar el costo de entregar energía durante el
período de tiempo como 𝑁
𝑇𝑚𝑎𝑥 𝑔𝑒𝑛

𝐹 𝑡𝑜𝑡𝑎𝑙 = ෍ ෍ 𝐹𝑖 𝑃𝑖𝑡
𝑡=1 𝑖=1
Sujeto a
𝑁𝑔𝑒𝑛
𝑃𝑖𝑡+1 = 𝑃𝑖𝑡 + ∆𝑃𝑖
෍ 𝑃𝑖𝑡 = 𝑃𝑙𝑜𝑎𝑑
𝑡
𝑓𝑜𝑟 𝑡 = 1 … 𝑡𝑚𝑎𝑥
𝑖=1
−∆𝑃𝑖𝑚𝑖𝑛 ≤ ∆𝑃𝑖 ≤ ∆𝑃𝑖𝑚𝑎𝑥
 Despacho económico

3.7 FUNCIÓN DE COSTE DE PRODUCCIÓN DE GENERACIÓN

COMPUESTA

Entonces se puede desarrollar una curva de costo compuesto para las unidades 1, 2, …, Ngen:

𝐹𝑠 𝑃𝑠 = 𝐹1 𝑃1 + ⋯ . 𝐹𝑁𝑔𝑒𝑛 𝑃𝑁𝑔𝑒𝑛
 Despacho económico

Donde
𝑃𝑠 = 𝑃1 + ⋯ + 𝑃𝑁𝑔𝑒𝑛

𝑑𝐹1 𝑑𝐹2 𝑑𝐹𝑁𝑔𝑒𝑛

= =⋯= =𝜆
𝑑𝑃1 𝑑𝑃2 𝑑𝑃𝑁𝑔𝑒𝑛

Si una de las unidades llega a un límite, su salida se mantiene constante.

Un procedimiento simple que permite generar Fs(Ps) consiste en ajustar λ de λ min a λ max en
incrementos específicos, donde
𝑑𝐹𝑖
𝜆𝑚𝑖𝑛 = 𝑚𝑖𝑛 𝑖 = 1 … 𝑁𝑔𝑒𝑛
𝑑𝑃𝑖

𝑑𝐹𝑖
𝜆𝑚𝑎𝑥 = 𝑚𝑎𝑥 𝑖 = 1 … 𝑁𝑔𝑒𝑛
𝑑𝑃𝑖
 Despacho económico

En cada incremento, calcule el consumo total de

combustible y la potencia de salida total para todas las
unidades. Estos puntos representan puntos en la curva
Fs(Ps). Los puntos se pueden usar directamente
asumiendo que Fs(Ps) consta de segmentos de línea
recta entre los puntos, o se puede ajustar una curva
suave a los puntos usando un programa de ajuste de
mínimos cuadrados. Tenga en cuenta, sin embargo, que
tales curvas suaves pueden tener propiedades
indeseables como la falta de convexidad (por ejemplo,
la primera derivada no es monótonamente creciente).
El procedimiento para generar los puntos sobre Fs(Ps)
se muestra en la Figura 3.9.
 Despacho económico

3.8 PUNTO BASE Y FACTORES DE PARTICIPACIÓN

Este método asume que el problema de despacho económico tiene que resolverse repetidamente
moviendo los generadores de un horario económicamente óptimo a otro a medida que la carga
cambia en una cantidad razonablemente pequeña. Partimos de un cronograma dado: el punto base.
Luego, el programador asume un cambio de carga e investiga cuánto necesita moverse cada unidad
generadora (es decir, “participar” en el cambio de carga) para que la nueva carga se atienda en el
punto de operación más económico.
 Despacho económico

Suponga que tanto la primera como la segunda derivada en la función de costo versus potencia de
salida están disponibles (es decir, existen tanto F′i como F″i). La curva de costo incremental de la i-
ésima unidad se muestra en la figura 3.12. A medida que la carga unitaria cambia en una cantidad
∆Pi, el costo incremental del sistema se mueve de λ0 a λ0 + ∆λ. Para un pequeño cambio en la
potencia de salida de esta sola unidad,

∆𝜆𝑖 = ∆𝜆 ≡ 𝐹𝑖′′ ∆𝑃𝑖

Esto es cierto para cada una de las N unidades del sistema, de modo que
∆𝜆
∆𝑃1 = ′′
𝐹1
∆𝜆
∆𝑃2 =
𝐹2′′

∆𝜆
∆𝑃𝑁 = ′′
𝐹𝑁
 Despacho económico

El cambio total en la generación (= cambio en la demanda total del sistema) es, por supuesto, la
suma de los cambios de unidades individuales. Sea PD la demanda total de los generadores (donde
PD = Pcarga + Ploss), entonces
∆𝑃𝐷 = ∆𝑃1 + ∆𝑃2 + ⋯ ∆𝑃𝑁

Entonces el factor de participación para cada unidad es

∆𝑃𝑖 1Τ𝐹𝑖′′
=
∆𝑃𝐷 σ𝑖 1Τ𝐹𝑖′′
 Despacho económico

Ejemplo 3G: Partiendo de la solución económica óptima encontrada en el Ejemplo 3A, utilice el
método del factor de participación para calcular el despacho para una carga total de 900 MW.
Usando la Ecuación 3.8,
−1
∆𝑃1 0.003124 320.10
= −1 + 0.00388 −1 −1 = = 0.47
∆𝑃𝐷 0.003124 + 0.00964 681.57

Similarmente
∆𝑃2 0.00388 −1
= = 0.38
∆𝑃𝐷 681.57
∆𝑃3 103.73
= = 0.15
∆𝑃𝐷 681.57

∆𝑃𝐷 = 900 − 850 = 50

 Despacho económico

El nuevo valor de la generación se calcula utilizando

∆𝑃𝑖
𝑃𝑛𝑒𝑤𝑖 = 𝑃𝑏𝑎𝑠𝑒𝑖 + ∆𝑃𝐷 𝑓𝑜𝑟 𝑖 = 1,2,3
∆𝑃𝐷
Entonces para cada unidad

𝑃𝑛𝑒𝑤1 = 392.2 + 0.47 50 = 416.7

𝑃𝑛𝑒𝑤2 = 334.6 + 0.38 50 = 353.6

𝑃𝑛𝑒𝑤3 = 122.2 + 0.15 50 = 129.7

 Despacho económico

3.9 DESPACHO DEL SISTEMA TÉRMICO CON CONSIDERACIÓN DE

PÉRDIDAS DE RED
La figura 3.13 muestra simbólicamente un sistema de generación de energía completamente
térmico conectado a una barra de carga equivalente a través de una red de transmisión. El problema
de despacho económico asociado con esta configuración particular es un poco más complicado de
establecer que el caso anterior.

𝑁𝑔𝑒𝑛

𝑃𝑙𝑜𝑎𝑑 + 𝑃𝑙𝑜𝑠𝑠 − ෍ 𝑃𝑖 = ∅ =0
𝑖=1
 Despacho económico

El mismo procedimiento se sigue en el sentido formal para establecer las condiciones necesarias para
una solución de operación de mínimo costo. La función de Lagrange se muestra en la Ecuación 3.10. Al
tomar la derivada de la función de Lagrange con respecto a cada una de las salidas de potencia
individuales, Pi, debe reconocerse que la pérdida en la red de transmisión, Ploss, es una función de las
impedancias de la red y las corrientes que fluyen en la red.

𝜕ℒ 𝑑𝐹𝑖 𝜕𝑃𝑙𝑜𝑠𝑠
= − 𝜆 1− =0
𝜕𝑃𝑖 𝑑𝑃𝑖 𝜕𝑥𝑃𝑖

𝑑𝐹𝑖 𝜕𝑃𝑙𝑜𝑠𝑠
+ 𝜆 =𝜆
𝑑𝑃𝑖 𝜕𝑃𝑖

𝑁𝑔𝑒𝑛

𝑃𝑙𝑜𝑎𝑑 + 𝑃𝑙𝑜𝑠𝑠 − ෍ 𝑃𝑖 =0
𝑖=1
 Despacho económico

El otro enfoque básico para la solución de este problema es incorporar las ecuaciones de flujo de
potencia como restricciones esenciales en el establecimiento formal del problema de optimización. Este
enfoque general se conoce como el flujo de potencia óptimo.

Ejemplo 3H: Comenzando con las mismas unidades y costos de combustible que en el Ejemplo 3A,
incluiremos una expresión de pérdida simplificada:

𝑃𝑙𝑜𝑎𝑑 = 0.00003𝑃12 + 0.00009𝑃22 + 0.00012𝑃32

Esta fórmula de pérdida simplificada será suficiente para mostrar las dificultades en el cálculo de un
envío para el que se contabilizan las pérdidas. Tenga en cuenta que las fórmulas de pérdida del mundo
real son más complicadas que la utilizada en este ejemplo.
Aplicando las ecuaciones de coordinación,
𝑑𝐹𝑖 𝜕𝑃𝑙𝑜𝑠𝑠
= 𝜆 1−
𝑑𝑃𝑖 𝜕𝑥𝑃𝑖
 Despacho económico

entonces
7.92 + 0.003124𝑃1 = 𝜆 1 − 2 0.00003 𝑃1

Similarmente para P2 y P3

7.85 + 0.00388𝑃2 = 𝜆 1 − 2 0.00009 𝑃2

7.97 + 0.00964𝑃3 = 𝜆 1 − 2 0.00012 𝑃3

𝑃1 + 𝑃2 + 𝑃3 − 850 − 𝑃𝑙𝑜𝑠𝑠 = 0
 Despacho económico

Ya no tenemos un conjunto de ecuaciones lineales como en el Ejemplo 3A. Esto requiere un

procedimiento de solución más complejo como sigue.

Paso 1 Elija un conjunto de valores iniciales para P1, P2 y P3 que sumen la carga.

Paso 2 Calcular las pérdidas incrementales ∂ Ploss/∂ Pi así como las pérdidas totales Ploss. Las
pérdidas incrementales y las pérdidas totales se considerarán constantes hasta que volvamos al paso 2.

Paso 3 Calcule el valor de 𝜆 que hace que P1, P2 y P3 sumen la carga total más las pérdidas. Ahora esto
es tan simple como los cálculos del Ejemplo 3A ya que las ecuaciones son nuevamente lineales.

Paso 4 Compare P1, P2 y P3 del paso 3 con los valores utilizados al comienzo del paso 2. Si no hay un
cambio significativo en alguno de los valores, vaya al paso 5; de lo contrario, vuelva al paso 2.

Paso 5 Listo.
 Despacho económico

Usando este procedimiento, obtenemos:

Paso 1 Seleccione los valores iniciales de P1, P2 y P3 como

𝑃1 = 400.0 𝑀𝑊
𝑃2 = 300.0 𝑀𝑊
𝑃3 = 150.0 𝑀𝑊
Paso 2 Las pérdidas incrementales son
𝜕𝑃𝑙𝑜𝑠𝑠
= 2 0.00003 400 = 0.0240
𝜕𝑃1

𝜕𝑃𝑙𝑜𝑠𝑠
= 2 0.00009 300 = 0.0540
𝜕𝑃2

𝜕𝑃𝑙𝑜𝑠𝑠
= 2 0.00012 150 = 0.0360
𝜕𝑃3
 Despacho económico

Total losses are 15.6 MW.

7.92 + 0.003124𝑃1 = 𝜆 1 − 0.0240 = 𝜆 0.9760

7.85 + 0.00388𝑃2 = 𝜆 1 − 0.0540 = 𝜆 0.9460

7.97 + 0.000964𝑃3 = 𝜆 1 − 0.0360 = 𝜆 0.9640

𝑃1 + 𝑃2 + 𝑃3 − 850 − 15.6 = 𝑃1 + 𝑃2 + 𝑃3 − 865.6 = 0

Estas ecuaciones ahora son lineales, por lo que podemos resolver para λ directamente. los resultados
son
𝜆 = 9.5252 $Τ𝑀 𝑊ℎ
y las salidas del generador resultantes son

𝑃1 = 440.68 𝑀𝑊 𝑃2 = 299.12 𝑀𝑊 𝑃3 = 125.77 𝑀𝑊

 Despacho económico

Paso 4 Dado que estos valores para P1, P2 y P3 son bastante diferentes de los valores iniciales,
regresaremos al paso 2.

Paso 2 Las pérdidas incrementales se recalculan con los nuevos valores de generación:
𝜕𝑃𝑙𝑜𝑠𝑠
= 2 0.00003 440.68 = 0.0264
𝜕𝑃1

𝜕𝑃𝑙𝑜𝑠𝑠
= 2 0.00009 299.12 = 0.0538
𝜕𝑃2

𝜕𝑃𝑙𝑜𝑠𝑠
= 2 0.00012 125.77 = 0.0301
𝜕𝑃3

Las pérdidas totales son de 15,78 MW.

Paso 3 Las nuevas pérdidas incrementales y las pérdidas totales se incorporan a las ecuaciones y se
resuelve un nuevo valor de l y P1, P2 y P3 para
 Despacho económico

Paso 3 Las nuevas pérdidas incrementales y las pérdidas totales se incorporan a las ecuaciones y se
resuelve un nuevo valor de l y P1, P2 y P3 para

7.92 + 0.003124𝑃1 = 𝜆 1 − 0.0264 = 𝜆 0.9736

7.85 + 0.00388𝑃2 = 𝜆 1 − 0.0538 = 𝜆 0.9462

7.97 + 0.000964𝑃3 = 𝜆 1 − 0.0301 = 𝜆 0.9699

𝑃1 + 𝑃2 + 𝑃3 − 850 − 15.78 = 𝑃1 + 𝑃2 + 𝑃3 − 865.78 = 0

resulting in λ = 9.5275 $/MWh and

𝑃1 = 433.94 𝑀𝑊 𝑃2 = 300.11 𝑀𝑊 𝑃3 = 131.74 𝑀𝑊

 Despacho económico

La Tabla 3.2 resume el proceso iterativo utilizado para resolver este problema.
¡MUCHAS GRACIAS!