DESPACHO ECONÓMICO EN SISTEMAS

ELÉCTRICOS DE POTENCIA

1.- INTRODUCCIÓN

En la operación de los SEP's se combinan dos aspectos muy importantes, los

cuales están relacionados: la seguridad y la economía.
La seguridad en un sistema de potencia se puede definir como la capacidad del
sistema de potencia para suministrar la demanda ante la presencia de alguna
contingencia. Las contingencias se identifican en un estándar de seguridad para
verificar el funcionamiento del sistema ante la ocurrencia de las mismas. Para
operar el sistema de manera segura se deberán de tener márgenes operativos
adecuados, los cuales se determinan por medio de un análisis de contingencias.
El otro aspecto importante en la operación de un SEP es el económico, el cual
está relacionado principalmente con el combustible que se consume durante un
periodo específico de tiempo. Así, el objetivo principal en la operación de un sistema
de potencia es suministrar la demanda a un costo mínimo, cumpliendo requisitos de
calidad y seguridad operativa.
En mercados de energía, los factores físicos (la red de transmisión y su
operación) y los factores económicos (la estructura del mercado y sus reglas) son
elementos determinantes en su funcionamiento.
Además es fundamental tener en cuenta el marco legal regulatorio de cada país
que condicionará las actuaciones y medidas que será necesario emplear.
De lo anterior podemos entonces definir en la actualidad a un despacho
económico, como la acción de suministrar la demanda de un sistema eléctrico
de potencia a un costo mínimo cumpliendo requisitos de calidad y seguridad,
tomando en cuenta a la vez el marco legal regulatorio de cada país.

1.1.- ANÁLISIS DE COSTOS

La operación del sistema requiere una supervisión permanente del mismo,

considerando la continuidad del servicio al menor costo posible. El problema del
suministro de energía eléctrica a bajo costo está influenciado por:
 La eficiencia del parque generador
 El costo de su instalación
 El costo de combustible para las máquinas térmicas

Los costos involucrados en la producción de energía pueden clasificarse en:

 Costos fijos
 Costos variables
1.1.1.- COSTOS FIJOS

Incluyen las inversiones de capital, los intereses de los préstamos, los salarios,
los impuestos y otros gastos que son independientes de la demanda del sistema.
Los responsables por la operación directa del SEP, tienen un pequeño control sobre
estos costos.

1.1.2.- COSTOS VARIABLES

Estos costos pueden ser controlados por los operadores del SEP, y dependen,
entre otras cosas, de:

 La confiabilidad de todos los componentes del SEP(Máquinas Síncronas,

Líneas, Transformadores, etc,.)
 El tipo de combustible usado en su operación
 El control de las pérdidas en las Líneas de Transmisión
 La forma de operar en forma conjunta el parque hidrotérmico y las energías
renovables para atender los requisitos diarios de la demanda
 La compra y venta de energía

La demanda puede ser suministrada por un parque generador que puede estar
compuesto de diferentes tipos de centrales: termoeléctricas convencionales,
nucleares, hidráulicas, tecnologías basadas en energías renovables y la importación
de energía.

El problema fundamental es determinar:

1.- La combinación de dichas generaciones (Asignación de Unidades) para
atender la demanda pronosticada.
2.- La potencia que deberá suministrar cada unidad generadora de las elegidas
en el primer punto, para que el costo global de la operación sea mínimo
(Despacho Económico).

Debido a que el combustible utilizado por las máquinas térmicas puede ser de
diferentes tipos (gas natural, gasoil, material nuclear, carbón, etc), con costos
diferentes y variables entre sí, y además que la carga de un sistema varía
continuamente y en forma aleatoria con el tiempo, es necesario que el problema de
la operación económica deba ser replanteado frecuentemente, reprogramando la
distribución de la generación, en saltos discretos de tiempo.
Cuando se dispone de agua para la generación hidroeléctrica, dicha
disponibilidad puede tener diferentes valores a cada instante, y por lo tanto así su
“costo de agua” asociado. El uso de esta generación debe de ser integrado a la
operación del SEP para que este funcione al menor costo posible (Operación
Hidrotérmica).
 El intercambio de energía, entre sistemas interconectados, así como las energías
renovables, pueden ser aprovechados de manera ventajosa para minimizar los
costos de combustible cuando existan diferencias significativas en los costos de
generación de estos sistemas.

1.2.- CURVAS ENTRADA-SALIDA

En los estudios de despacho económico es común caracterizar a las unidades

termoélectricas con las denominadas curvas “Entrada-Salida”, que definen el calor
(Gcal/hr) necesario para generar en nivel de potencia (Mw).
Generalmente estas curvas son caracterizadas con polinomios de segundo
grado y un ejemplo de estas se puede observar en la figura 1, siendo la ecuación 1
la que las representa.
ENTRADA Gcal/Hr

SALIDA MW

Figura 1 Curva de Entrada-Salida

𝑄(𝑝) = 𝑎 + 𝑏𝑝 + 𝑐𝑝2 𝐸𝑐 − 1

Donde:

𝑄(𝑝) = 𝐹𝑢𝑛𝑐𝑖ó𝑛 𝑑𝑒 𝐶𝑜𝑛𝑠𝑢𝑚𝑜 (𝐺𝑐𝑎𝑙/𝐻𝑟)

𝑎 = 𝐶𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒𝑙 𝑡é𝑟𝑚𝑖𝑛𝑜 𝑖𝑛𝑑𝑒𝑝𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 (𝐺𝑐𝑎𝑙/𝐻𝑟)
𝑏 = 𝐶𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒𝑙 𝑡é𝑟𝑚𝑖𝑛𝑜 𝑙𝑖𝑛𝑒𝑎𝑙 (𝐺𝑐𝑎𝑙/𝑀𝑤𝐻𝑟)
𝑐 = 𝐶𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒𝑙 𝑡é𝑟𝑚𝑖𝑛𝑜 𝑐𝑢𝑎𝑑𝑟𝑎𝑡𝑖𝑐𝑜 (𝐺𝑐𝑎𝑙/𝑀𝑤 2 𝐻𝑟)
𝑝 = 𝑁𝑖𝑣𝑒𝑙 𝑑𝑒 𝑔𝑒𝑛𝑒𝑟𝑎𝑐𝑖ó𝑛 (𝑀𝑤)

La unidad típica de generación eléctrica comprende el sistema caldera-turbina-

generador, esto consiste en una caldera que genera vapor mediante algún tipo de
combustible y este vapor acciona una turbina acoplada mecánicamente a un
generador eléctrico, la salida es potencia eléctrica en MW’s que se entrega al
sistema.
La curva que describe el comportamiento de la caldera-turbina-generador es la
curva de entrada-salida.
Esta curva representa la cantidad de combustible o el costo del combustible que
la unidad consume por hora, para la potencia a la que la unidad está generando,
esta relación como se dijo al principio se representa generalmente por un polinomio
de segundo orden.

1.3.- CURVA DE CONSUMO ESPECÍFICO (HEAT RATE)

Para los diversos puntos de prueba de una unidad, la relación de calor expresada
en Kcal/KwHr o Gcal/MwHr contra la salida expresada en MW’s, es llamada curva
de consumo específico o heat rate figura 2, siendo la ecuación 2 la que la
representa.
$/MwHr
o
ENTRADA Gcal/MwHr

SALIDA MW

Figura 2 Curva de Consumo Específico

𝑄(𝑝) 𝑎
= + 𝑏 + 𝑐𝑝 𝐸𝑐 − 2
𝑝 𝑝

La curva de consumo específico o head rate (Régimen Térmico), representa el

costo por MwHr que tiene la unidad, para una determinada potencia de salida; esto
significa que las unidades serán las siguientes: en las ordenadas $/MwHr o bien
Gcal/MwHr y en las abscisas MW’s.
El método de prueba para la obtención de esta curva consiste básicamente en
tener una carga fija en la unidad y obtener el gasto de combustible consumido en
este lapso de tiempo, después se varía la carga de la unidad y se van tomando los
diferentes valores de consumo de combustible.
Periódicamente se realizan pruebas de consumo específico a las unidades,
generalmente cuando entran y salen de un mantenimiento, se comparan los valores
obtenidos con las curvas de los fabricantes de las unidades y se evalúan los
desplazamientos entre las curvas.
La curva de consumo específico empieza con un valor alto para generación
cercana a cero y va decreciendo hasta llegar a un mínimo, aproximadamente un
85% de la potencia máxima, y después comienza a ascender. A este punto mínimo
se le conoce como punto de máxima eficiencia 𝜶.
La curva de consumo específico también se puede encontrar si se conoce el
polinomio de la curva entra-salida y esta se divide entre la potencia de salida.
La curva de consumo específico representa la eficiencia en la conversión de la
energía suministrada por el combustible al sistema caldera-turbina-generador y
transformado en potencia eléctrica.

1.4.- COSTOS INCREMENTALES

La curva de costos incrementales representa la pendiente o derivada de la curva

entrada-salida, indica el incremento en costo o en consumo de combustible que se
requiere para un incremento en la potencia de salida. Conociendo la curva de
entrada-salida se puede encontrar la curva de costo incremental derivando la curva
de entrada-salida. La ecuación 3 representa a esta curva de costo incremental.

𝐶𝐼 = 𝑏 + 2𝑐𝑝 ($/𝑀𝑤𝐻𝑟) 𝐸𝑐 − 3

La característica de esta curva es ascendente, ya que para cada aumento de

producción corresponde un aumento de combustible.
En pruebas reales que se realizan a las unidades se ha encontrado que existen
algunos picos en esta curva los cuales son suavizados mediante técnicas de
segmentación. Las discontinuidades se presentan por el cierre y apertura de
válvulas en la turbina, presentando una restricción en la operación del generador.
La curva de costos incrementales es usada en estudios de despacho económico
de generación y en la planeación de sistemas de potencia. Esta curva también es
utilizada para fijar precios de venta de energía en algunas centrales.

“Cabe mencionar que existen procedimientos técnicos relacionados con

las energías renovables que hay que tomar en cuenta al momento de incluir
dichas generaciones dentro de un despacho económico. En estas notas solo
se cubrirá el aspecto relacionado a las Unidades Termoeléctricas”.
1.5.- Obtención de la curva entrada-salida por mediciones de la
planta

Una forma de obtener la curva de entrada-salida es tomando mediciones en

planta para el generador bajo análisis ante diferentes potencias de generación,
donde la medición para una potencia especifica del generador, representa la
cantidad de combustible que la unidad consume por hora.
A continuación se determinará la curva entrada-salida para un generador de 300
MW’s de capacidad y cuyas mediciones se presentan en la tabla 1:

Generación % Generación MW’s Consumo (Gcal/Hr)

100 300 728
87.5 262.5 618.8
81.25 243.75 575.2
75 225 531.5
68.75 206.25 495.1
62.5 187.5 448.7
52.7 158.1 419.8

Tabla 1 Mediciones Generador 300 MW’s

La determinación de los puntos de medición requiere de controlar el proceso

térmico y permitir su estabilización. Para obtener una función analítica de la curva
entrada-salida con los puntos de medición de consumo del generador se asume un
polinomio que ajuste adecuadamente las mediciones obtenidas en las pruebas.
Típicamente se utiliza una función cuadrática para aproximar la curva de entrada-
salida. Existen varios métodos matemáticos para hacer esto dentro de los cuales se
puede mencionar al método de mínimos cuadrados. Para nuestro caso utilizaremos
la función polyfit del MatLab.

Ejemplo_1 Obtener el polinomio que caracteriza a la curva entrada-salida del

generador con las mediciones de la tabla 1.

X = [300 262.5 243.75 225 206.25 187.5 158.1]

Y = [728 618.8 575.2 531.5 495.1 448.7 419.8]

C = polyfit(X,Y,2)

C =

0.0059 -0.8643 385.0553

Por lo tanto

𝑄(𝑝) = 0.0059𝑝2 − 0.8643𝑝 + 385.0553

Para comprobar utilizamos el comando polyval de MatLab

Z = polyval(C,X)

Z =

728.7678 619.8516 572.4598 529.7788 491.8088 458.5496 415.8836

De esta manera queda determinada la ecuación de entrada-salida del generador

bajo estudio, la cual se repite a continuación.

𝑄(𝑝) = 0.0059𝑝2 − 0.8643𝑝 + 385.0553

2.- DESPACHO ECONÓMICO BÁSICO (sin límites y pérdidas)

El despacho económico consiste básicamente en utilizar los recursos energéticos

(térmicos, hidráulicos, solares, eólicos, etc.) disponibles para la generación de
energía eléctrica en una forma óptima de tal manera que se cubra la demanda
eléctrica a un mínimo costo, con un determinado grado de confiabilidad, calidad,
seguridad, tomando en cuenta el marco legal regulatorio de cada país.
Este consiste también en conocer la cantidad de potencia que debe suministrar
cada generador para satisfacer una condición de demanda de los consumidores
minimizando los costos de generación del sistema eléctrico sujeto a diferentes
restricciones operativas como son la de los generadores, la red de transmisión, así
como la seguridad del sistema tomando en cuenta el marco legal regulatorio. Lo
anterior representa una sintonización de todos los generadores cumpliendo
con las diferentes restricciones y operando a un mismo costo incremental.

2.1.- TÉCNICAS DE OPTIMIZACIÓN

Dado que el problema del despacho económico es un problema de optimización,

ya que se requiere despachar las unidades de generación de forma tal que estén
operando todas con un mismo costo (Costo Incremental), a continuación se
expondrá la teoría básica relacionada con las técnicas de optimización.

Para maximizar (max) o minimizar (Min) una función objetivo (FO) del tipo:

𝐹𝑂 = 𝑓(𝑋1 , 𝑋2 , … . , 𝑋𝑛 ) 𝐸𝑐 - 4

Se deben encontrar las primeras derivadas de dicha función respecto a cada

variable 𝑋, e igualarlas a cero. Luego sus segundas derivadas se usarán para
determinar si la condición encontrada corresponde a un máximo (max), un mínimo
(min) o a un punto de inflexión de la función FO.
 Como el objetivo es maximizar o minimizar una función matemática, se llamará
a esta función, Función Objetivo (FO), ecuación 4.

En los problemas reales, una función que se quiere maximizar o minimizar

generalmente tiene restricciones en los valores que pueden asumir las variables.
Estas restricciones pueden ser de igualdad o desigualdad. La región definida por
las restricciones generalmente se llama dominio de las variables independientes
(V.I). Si no existe ningún valor de las V.I que satisface todas las restricciones,
entonces el problema no tiene solución posible.

Para comprender lo dicho con anterioridad se desarrollará el siguiente ejemplo:

Minimizar la siguiente función objetivo:

𝑓(𝑋1 , 𝑋2 ) = 0.25𝑋1 2 + 𝑋2 2

Sujeta a las siguientes restricciones:

𝜔(𝑋1 , 𝑋2 ) = 0

𝜔(𝑋1 , 𝑋2 ) = 5 − 𝑋1 − 𝑋2

Si graficamos la función anterior 𝑓(𝑋1 , 𝑋2 ), junto con las restricciones se observa

que es necesario encontrar un punto mínimo donde una elipse toque a las
restricciones, en este caso 𝜔(𝑋1 , 𝑋2 ) = 5 − 𝑋1 − 𝑋2. Esto se puede apreciar en la
figura 3.

Figura 3 Gráfica de la función objetivo y restricciones

 En el punto (𝑋1 , 𝑋2 ) se calcula el vector gradiente de la función objetivo FO,
∇ 𝑓(𝑋1 , 𝑋2 ). Este vector es perpendicular a la FO en el punto (𝑋1 , 𝑋2 ), pero no lo es
a la restricción 𝜔, por lo tanto tiene una componente diferente de cero sobre 𝜔. Si
se realiza un movimiento por 𝜔 en el sentido de la componente diferente de cero
del ∇ 𝑓, se incrementará el valor de la FO. Por lo tanto para minimizar la FO, se
deberá ir en el sentido opuesto al sentido de dicha componente. En el punto óptimo
el ∇ 𝑓 (normal) es perpendicular tanto a la FO como a 𝜔. Entonces para garantizar
que ∇ 𝑓 es normal a 𝜔, se requiere que el ∇ 𝑓 y el gradiente ∇ 𝜔, sean vectores
linealmente dependientes, es decir vectores que tienen la misma dirección y no
necesariamente el mismo sentido y magnitud. Esto se puede apreciar en la figura
4.

Figura 4 Vector ∇ 𝑓 y gradiente ∇ 𝜔

Matemáticamente se puede representar lo anterior de la siguiente manera:

∇𝑓 + 𝜆∇𝜔 = 0 𝐸𝑐 − 5

Esto es los dos gradientes se suman para que se anulen, cuando uno de ellos
esta multiplicado por un factor de escala. Este factor de escala variable 𝜆, es
llamado “Multiplicador de Lagrange”. En lugar de usar los gradientes se integra
la ecuación 4 y se obtiene:

𝐿(𝑋1 , 𝑋2 , … . . 𝑋𝑛 , 𝜆) = 𝑓(𝑋1 , … . . 𝑋𝑛 ) + 𝜆𝜔(𝑋1 , … . . 𝑋𝑛 ) 𝐸𝑐 − 6

Esta ecuación es llamada “Ecuación de Lagrange”, y consiste de (n+1) variables.

Para encontrar el óptimo de la FO, se calculan las derivadas parciales de L en
función de sus (n+1) variables desconocidas, es decir:
 ℓ𝐿
=0
ℓ𝑋1
•
•
ℓ𝐿
=0 en el punto óptimo 𝐸𝑐 − 7
ℓ𝑋𝑛

ℓ𝐿
=0
ℓ𝜆

Volviendo al ejemplo numérico y usando las ecuaciones 6 y 7, el problema se

resuelve:

𝐿(𝑋1 , 𝑋2 , … . . 𝑋𝑛 , 𝜆) = 𝑓(𝑋1 , … . . 𝑋𝑛 ) + 𝜆𝜔(𝑋1 , … . . 𝑋𝑛 )

𝐿(𝑋1 , 𝑋2 , … . . 𝑋𝑛 , 𝜆) = 0.25𝑋12 + 𝑋22 + 𝜆(5 − 𝑋1 − 𝑋2 )

ℓ𝐿
= 0.5𝑋1 − 𝜆 = 0
ℓ𝑋1

ℓ𝐿
= 2𝑋2 − 𝜆 = 0 Sistema de 3 ecuaciones con 3 incógnitas
ℓ𝑋2

ℓ𝐿
= 5 − 𝑋1 − 𝑋2 = 0
ℓ𝜆

Se debe notar de la expresión anterior, que la última ecuación es la restricción

original.
Resolviendo el sistema matricialmente se obtiene:

0.5 0 −1 𝑋1 0
[0 ] [
2 −1 𝑋2 ] = [ 0]
−1 −1 0 𝜆 −5

donde

𝑋1 = 4 𝑋2 = 1 𝜆=2

Cuando se tienen 𝑚 restricciones de igualdad, el problema es:

𝐹𝑂 = 𝑓(𝑋1 , … . . , 𝑋𝑛 ) → 𝐹𝑢𝑛𝑐𝑖ó𝑛 𝑂𝑏𝑗𝑒𝑡𝑖𝑣𝑜

𝜔1 (𝑋1 , … . , 𝑋𝑛 ) = 0

•
• Restricciones de igualdad

𝜔𝑚 (𝑋1 , … . , 𝑋𝑛 ) = 0

La suma de los vectores gradientes, y la ecuación de lagrange resultan:

∇𝑓 + 𝜆1 ∇𝜔1 + ⋯ + 𝜆𝑚 ∇𝜔𝑚 = 0 𝐸𝑐 − 8

𝐿(𝑋1 , 𝑋2 , . , 𝑋𝑛 , 𝜆) = 𝑓(𝑋1 , . , 𝑋𝑛 ) + 𝜆1 𝜔1 (𝑋1 , . , 𝑋𝑛 ) + 𝜆𝑚 𝜔𝑚 (𝑋1 , . , 𝑋𝑛 ) 𝐸𝑐 − 9

y las condiciones para encontrar el óptimo son:

ℓ𝐿
=0
ℓ𝑋1
•
•
ℓ𝐿
=0
ℓ𝑋𝑛
•
• en el punto óptimo 𝐸𝑐 − 10
ℓ𝐿
=0
ℓ𝜆1
•
•
ℓ𝐿
=0
ℓ𝜆𝑚

En el caso de tener un conjunto 𝒊 de restricciones de igualdad (𝜔𝑖 = 0), y uno

de desigualdad (𝑔𝑑 = 0), es decir:

𝐹𝑂 = 𝑓(𝑋1 , … . . , 𝑋𝑛 ) → 𝐹𝑢𝑛𝑐𝑖ó𝑛 𝑂𝑏𝑗𝑒𝑡𝑖𝑣𝑜

𝜔𝑖 (𝑋1 , … . , 𝑋𝑛 ) = 0
Restricciones de igualdad y desigualdad
𝑔𝑑 (𝑋1 , … . , 𝑋𝑛 ) ≤ 0
 Este sistema puede ser resuelto usando las condiciones de Kuhn-Tucker,
mostradas a continuación:
𝑑𝐹𝑖
= 𝜆 𝑃𝑎𝑟𝑎 𝑃𝑖,𝑚𝑖𝑛 < 𝑃𝑖 < 𝑃𝑖,𝑚𝑎𝑥
𝑑𝑃𝑖
𝑑𝐹𝑖
≤ 𝜆 𝑃𝑎𝑟𝑎 𝑃𝑖 = 𝑃𝑖,𝑚𝑎𝑥
𝑑𝑃𝑖
𝑑𝐹𝑖
≥ 𝜆 𝑃𝑎𝑟𝑎 𝑃𝑖 = 𝑃𝑖,𝑚𝑖𝑛
𝑑𝑃𝑖

Este grupo de ecuaciones se puede interpretar de la siguiente manera:

Los generadores que operan entre sus límites de potencia tienen costos marginales
idénticos y de valor 𝜆, los que operan a su mínima potencia tienen un costo marginal
igual o mayor que 𝜆, mientras los que operan a su límite superior tienen un costo
igual o menor que 𝜆.

2.2.- DESPACHO ECONÓMICO TÉRMICO DESPRECIANDO

PÉRDIDAS EN LA TRANSMISIÓN

Suponer que el sistema a estudiar se muestra en la figura 5, y como no se

consideran las pérdidas en la red, la generación de cada uno de los generadores se
concentra en un bus común donde se debe de cumplir que el total de la carga debe
de ser igual a la generación producida.

Figura 5 Generación total igual a carga total

El caso consiste de N unidades térmicas para generación de electricidad

conectadas a una barra que alimenta una carga 𝑃𝐷 . La entrada representa el costo
$
de combustible por hora de consumo [ℎ] de la unidad térmica 𝑖, 𝐶𝑖 (𝑃𝑖 ). La salida
eléctrica 𝑃𝑖 es la potencia eléctrica en [𝑀𝑊] generada por el generador 𝑖 asociado
a la unidad térmica 𝑖. El costo total es por supuesto la suma de los costos
individuales, con la restricción de que la suma de las potencias eléctricas generadas
individualmente sea igual a la potencia de la demanda.

El problema es: minimizar el costo total, con la restricción de igualdad entre

lo generado y lo consumido.

Matemáticamente se puede plantear fácilmente:

La función objetivo (la que se quiere optimizar, en este caso minimizar) es el

costo total CT, sujeto a la restricción que la suma de las generaciones individuales
debe de ser igual a la potencia total consumida.

𝐶𝑇 = ∑ 𝐶𝑖 (𝑃𝑖 ) 𝐹𝑢𝑛𝑐𝑖ó𝑛 𝑜𝑏𝑗𝑒𝑡𝑖𝑣𝑜 𝑎 𝑚𝑖𝑛𝑖𝑚𝑖𝑧𝑎𝑟 𝐸𝑐 − 11

𝑖=1

𝑃𝐷 − ∑ 𝑃𝑖 = 𝜙 𝑅𝑒𝑠𝑡𝑟𝑖𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑖𝑔𝑢𝑎𝑙𝑑𝑎𝑑 𝐸𝑐 − 12
𝑖=1

Se puede observar que en este problema de optimización, no se hace referencia

a pérdidas en la transmisión ni a ningún límite en la operación del sistema, por lo
que puede usarse el método de Lagrange visto en el punto anterior, para encontrar
su solución.

La ecuación de Lagrange está dada por:

𝐿 = 𝐶𝑇 + 𝜆𝜙

usando la ecuación 5, se obtiene:

ℓ𝐿 𝑑𝐶𝑇 𝑑𝜙 𝑑𝐶𝑇 𝑑𝐶𝑇

= + 𝜆 = −𝜆 =0 ⟶ = 𝜆 = 𝐶𝐼𝑖 [$/𝑀𝑊ℎ] 𝐸𝑐 − 13
ℓ𝑃𝑖 𝑑𝑃𝑖 𝑑𝑃𝑖 𝑑𝑃𝑖 𝑑𝑃𝑖

Es decir: La condición necesaria para la existencia de una condición de

operación al mínimo costo para un sistema de potencia térmico, es que el
costo incremental de todas las unidades sea igual a un valor a determinar,
llamado 𝝀.

Por supuesto que a esta condición necesaria, hay que agregarle la condición de
igualdad entre la potencia generada y la consumida.
Existen además dos condiciones de desigualdad que debe de cumplir cada
unidad, esto es la potencia generada debe de ser mayor o igual que una potencia
mínima determinada, y menor o igual que la máxima potencia especificada.
 Lo anterior puede resumirse de la siguiente manera:

𝑑𝐶𝐼
= 𝐶𝐼𝑖 = 𝜆 𝑁 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛𝑒𝑠 (𝑣𝑒𝑟 𝐸𝑐 − 3.6)
𝑑𝑃𝑖
𝑁

𝑃𝐷 − ∑ 𝑃𝑖 = 𝜙 1 𝑒𝑐𝑢𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑖𝑔𝑢𝑎𝑙𝑑𝑎𝑑
𝑖=1
𝑃𝑖𝑚𝑖𝑛 ≤ 𝑃𝑖 ≤ 𝑃𝑖𝑚𝑎𝑥 2𝑁 𝑟𝑒𝑠𝑡𝑟𝑖𝑐𝑐𝑖𝑜𝑛𝑒𝑠 𝑑𝑒 𝑑𝑒𝑠𝑖𝑔𝑢𝑎𝑙𝑑𝑎𝑑

de lo anterior existen (N + 1) incógnitas, formadas por:

 Las N potencias de los N generadores cuando trabajan al mismo valor del

costo incremental 𝜆.
 El valor del costo incremental 𝜆.

Cuando se trabajan con las 2N desigualdades, entonces las condiciones

necesarias pueden representarse de la siguiente manera (relaciones de Khun-
Tucker):

𝑑𝐶𝑇
= 𝜆 𝑃𝑎𝑟𝑎 𝑃𝑖𝑚𝑖𝑛 ≤ 𝑃𝑖 ≤ 𝑃𝑖𝑚𝑎𝑥 𝐸𝑐 − 14
𝑑𝑃𝑖
𝑑𝐶𝑇
≤ 𝜆 𝑃𝑖 = 𝑃𝑖𝑚𝑎𝑥 𝐸𝑐 − 15
𝑑𝑃𝑖
𝑑𝐶𝑇
≥ 𝜆 𝑃𝑖 = 𝑃𝑖𝑚𝑖𝑛 𝐸𝑐 − 16
𝑑𝑃𝑖

En la figura 6 se muestra la representación gráfica de dichas ecuaciones:

Figura 6 Representación gráfica del valor del costo incremental 𝜆

Ejemplo_1

Calcular la potencia a la que deben de operar los generadores 1, 2 y 3, para

minimizar los costos de operación y abastecer una demanda de 850 MWs. Los datos
se presentan a continuación en la tabla 1:

Unidad Pot Máx Pot Mín Curva Entrada-salida

No MWs MWs MBtu/h
1 600 150 𝐶𝐶1 = 510 + 7.2𝑃1 + 0.00142𝑃12
2 400 100 𝐶𝐶2 = 310 + 7.85𝑃2 + 0.00194𝑃22
3 200 50 𝐶𝐶3 = 78 + 7.97𝑃3 + 0.00482𝑃32

Tabla 1 Datos de generadores Ejemplo_1

Costo de combustible Unidad 1 𝐾1 = 1.1[$/𝑀𝐵𝑡𝑢]

Costo de combustible Unidad 2 𝐾2 = 1.0[$/𝑀𝐵𝑡𝑢]
Costo de combustible Unidad 3 𝐾3 = 1.0[$/𝑀𝐵𝑡𝑢]

Solución:

Como primer paso se calcula el costo de combustible de cada unidad

𝐶1 (𝑃1 ) = 𝐶𝐶1 (𝑃1 ) ∗ 1.1 = 561 + 7.92𝑃1 + 0.00156𝑃12 [$/ℎ]

𝐶2 (𝑃2 ) = 𝐶𝐶2 (𝑃2 ) ∗ 1.0 = 310 + 7.85𝑃2 + 0.00194𝑃22 [$/ℎ]

𝐶3 (𝑃3 ) = 𝐶𝐶3 (𝑃3 ) ∗ 1.0 = 78 + 7.97𝑃3 + 0.00482𝑃32 [$/ℎ]

𝐶𝑇 = 𝐶1 (𝑃1 ) + 𝐶2 (𝑃2 ) + 𝐶3 (𝑃3 ) [$/ℎ]

Usando la ecuación 6-17, las condiciones para el despacho son:

𝑑𝐶𝑇
= 7.92 + 0.00312𝑃1 = 𝜆 [$/𝑀𝑊ℎ]
𝑑𝑃1

𝑑𝐶𝑇
= 7.85 + 0.00388𝑃2 = 𝜆 [$/𝑀𝑊ℎ]
𝑑𝑃2

𝑑𝐶𝑇
= 7.97 + 0.00964𝑃3 = 𝜆 [$/𝑀𝑊ℎ]
𝑑𝑃3

𝑃1 + 𝑃2 + 𝑃3 = 850 [𝑀𝑊𝑠]
 Resolviendo el sistema de 4 ecuaciones con 4 incógnitas, se encuentra:

0.00312 0 0 −1.0
𝐴= [ 0 0.00388 0 −1.0]
0 0 0.00964 −1.0
1.0 1.0 1.0 0

−7.92
𝐵 = [−7.85]
−7.97
850.0

Por lo tanto

𝐴𝜆 = 𝐵

𝜆 = 𝐴−1 𝐵

393.1698
334.6038
𝜆=[ ]
122.2264
9.1483

De lo anterior se observa que se cumple con todas las restricciones, esto quiere decir
que cada unidad se encuentra trabajando entre sus límites de operación, y que la suma de
todas las generaciones es igual a la demanda de 850 MWs. Lo anterior se puede apreciar
en la figura 7.

393.2 334.6 122.2

9.148

Figura 7 Generadores operando al mismo costo incremental

 También podemos observar de la figura 7 que las máquinas están operando a diferentes
capacidades de acuerdo al costo incremental que debe de ser el mismo para todas. Es decir
la máquina más barata es el generador 1, le sigue el generador 2 y por último el generador
3.

Ejemplo_2

Igual al ejemplo 1, pero considerando que la unidad 1 tiene un costo especifico

de:

Costo de combustible Unidad 1 𝐾1 = 0.9[$/𝑀𝐵𝑡𝑢]

De donde

𝐶1 (𝑃1 ) = 𝐶𝐶1 (𝑃1 ) ∗ 0.9 = 459 + 6.48𝑃1 + 0.00128𝑃12 [$/ℎ]

𝑑𝐶𝑇
= 6.48 + 0.00256𝑃1 = 𝜆 [$/𝑀𝑊ℎ]
𝑑𝑃1

La solución que se encuentra siguiendo el método anterior es:

0.00256 0 0 −1.0
𝐴= [ 0 0.00388 0 −1.0]
0 0 0.00964 −1.0
1.0 1.0 1.0 0

−6.48
𝐵 = [−7.85]
−7.97
850.0

Por lo tanto

𝐴𝜆 = 𝐵

𝜆 = 𝐴−1 𝐵

705.1467
112.1587
𝜆=[ ]
32.6946
8.2852
 De lo anterior se observa que la solución cumple con:

𝑃1 + 𝑃2 + 𝑃3 = 850 [𝑀𝑊𝑠]

Pero como las unidades 1 y 3 trabajan fuera de sus límites de operación, este problema se
resuelve usando las condiciones de Kuhn-Tucker.

Suponemos que la unidad 1 opera a potencia máxima 600 MWs, y la unidad 3 a potencia
mínima 50 MWs, por lo tanto el despacho se vuelve:

𝑃1 = 600 𝑀𝑊𝑠

𝑃2 = 200 𝑀𝑊𝑠

𝑃3 = 50 𝑀𝑊𝑠

De la ecuación 6-17, se observa que 𝜆 debe de ser igual al costo incremental de la

unidad 2, ya que ésta trabaja entre sus límites.

𝑑𝐶𝑇
| = 𝜆
𝑑𝑃2 𝑃2 =200

𝑑𝐶𝑇
| = 7.85 + 0.00388𝑃2 = 𝜆
𝑑𝑃2 𝑃2 =200

𝑑𝐶𝑇
| = 7.85 + 0.00388(200) = 𝜆
𝑑𝑃2 𝑃2 =200

𝑑𝐶𝑇
| = 𝟖. 𝟔𝟐𝟔 [$/𝑴𝑾𝒉]
𝑑𝑃2 𝑃2 =200

Calculando el costo incremental para las unidades 1 y 3

𝑑𝐶𝑇
| = 6.48 + 0.00256(600)
𝑑𝑃1 𝑃1 =600

𝑑𝐶𝑇
| = 𝟖. 𝟎𝟏𝟔 [$/𝑴𝑾𝒉]
𝑑𝑃1 𝑃1 =600

𝑑𝐶𝑇
| = 7.97 + 0.00964(50)
𝑑𝑃3 𝑃
1 =50
𝑑𝐶𝑇
| = 𝟖. 𝟒𝟓𝟐 [$/𝑴𝑾𝒉]
𝑑𝑃3 𝑃
1 =50

Podemos observar que el costo incremental (CI) de la unidad 1 es menor que el

𝝀, por lo tanto y usando la condición de Khun- Tucker 𝐸𝑐 − 14, la unidad 1 deberá
operar a su máxima potencia.
Sin embargo el costo incremental (CI) de la unidad 3 no es mayor que 𝝀, por lo
que la unidad 3 no está obligada a operar a su mínima potencia.
Por lo tanto para encontrar el despacho óptimo, se permite que el costo
incremental de las unidades 2 y 3 sea igual a 𝝀.

de donde:

𝑃1 = 600 𝑀𝑊𝑠

𝑑𝐶𝑇
= 7.85 + 0.00388𝑃2 = 𝜆
𝑑𝑃2

𝑑𝐶𝑇
= 7.97 + 0.00964𝑃3 = 𝜆
𝑑𝑃3

600 + 𝑃2 + 𝑃3 = 850 [𝑀𝑊𝑠]

𝑃2 + 𝑃3 = 250 [𝑀𝑊𝑠]

La solución que se encuentra siguiendo el método anterior es:

0.00388 0 −1
𝐴= [ 0 0.00964 −1]
1 1 0

−7.85
𝐵 = [−7.97]
250.0

Por lo tanto

𝐴𝜆 = 𝐵

𝜆 = 𝐴−1 𝐵
 187.1
𝜆 = [ 62.9 ]
8.576

Se observa que el despacho anterior cumple con las condiciones de Kuhn-Tucker,

ecuaciones 6-17,18 y 19.

𝑑𝐶𝑇
| = 𝟖. 𝟎𝟏𝟔 [$/𝑴𝑾𝒉]
𝑑𝑃1 𝑃1 =600

es menor que 𝝀, mientras que:

𝑑𝐶𝑇 𝑑𝐶𝑇
= = 𝝀
𝑑𝑃2 𝑑𝑃3

Lo anterior se puede apreciar en la figura 8.

187.1 62.9
600.0
8.576
8.016

Figura 8 Generadores operando a diferente costo incremental

2.2.1.- FACTORES DE PARTICIPACION

En un sistema eléctrico de potencia la demanda que este abastece, se considera

que está cambiando de una manera cíclica; es decir se ajusta a una curva de
demanda que se repite día con día (días festivos, lluviosos, normales, etc).
Esta demanda se hace más alta durante el día y parte de la tarde cuando las
cargas de las industrias son altas, sin embargo durante la noche y la madrugada
esta demanda será mucho menor.
 El despacho de generación óptimo necesita resolverse para diferentes
condiciones de carga, porque la demanda de carga 𝑃𝐷 es cambiante. Cuando los
cambios de carga son pequeños, es posible moverse de un despacho óptimo de
generación a otro sin resolver de nuevo el problema, usando FACTORES DE
PARTICIPACIÓN.

Podemos iniciar con un despacho óptimo de generación conocido,

𝑃10 , 𝑃20 , 𝑃30 , … , 𝑃𝑁0 ,
para una carga particular 𝑃𝐷 . Este despacho es tomado como punto
base y su correspondiente costo incremental es 𝜆0. Dejemos que se de un pequeño
incremento en la carga ∆𝑃𝐷 . Para satisfacer a este incremento de carga, las
generaciones van a ser incrementadas en ∆𝑃1 , ∆𝑃2 , ∆𝑃3 , ∆𝑃𝑁 , correspondiendo a un
incremento en el costo incremental ∆𝜆.

Como se sabe la curva Entrada-Salida para un generador 𝑖, se define en forma

general como:

𝐶𝑖 = 𝛼𝑖 𝑃𝑖2 + 𝛽𝑖 𝑃𝑖 + 𝛾𝑖 𝐸𝑐 − 17

Calculando el costo incremental 𝐶𝐼𝑖 , se obtiene:

𝐶𝐼𝑖 = 2𝛼𝑖 𝑃𝑖 + 𝛽𝑖 = 𝜆 𝐸𝑐 − 18

Cambios pequeños en el costo incremental producen cambios pequeños en la

generación, por lo tanto, haciendo 𝛽𝑖 en la ecuación 17:

∆2𝛼𝑖 𝑃𝑖 = ∆𝜆

∆𝜆
∆𝑃𝑖 = 𝑝𝑎𝑟𝑎 𝑖 = 1, … . , 𝑛 𝐸𝑐 − 19
2𝛼𝑖

El total de cambios en la generación es igual al total de cambios en la carga, por

lo tanto:

∑ ∆𝑃𝑖 = ∆𝑃𝐷 𝐸𝑐 − 20
𝑖=1

Sustituyendo la Ec-18 en la Ec-19

𝑁
1
∆𝑃𝐷 = ∆𝜆 ∑ 𝐸𝑐 − 21
2𝛼𝑖
𝑖=1

Obteniendo la relación entre las ecuaciones 18 y 20.

 ∆𝜆 1
∆𝑃𝑖 2𝛼𝑖 2𝛼𝑖
= = = 𝐾𝑖 𝑖 = 1, … … 𝑁 𝐸𝑐 − 22
∆𝑃𝐷 1 1
∆𝜆 ∑𝑁
𝑖=1 2𝛼 ∑ 𝑁
𝑖=1 2𝛼
𝑖 𝑖

∆𝑃𝑖
A la relación representada por el parámetro 𝐾𝑖 , se le llama Factor de
∆𝑃𝐷
Participación.

De lo anterior podemos afirmar que los cambios en las generaciones se pueden

dar como:

∆𝑃𝑖 = 𝐾𝑖 ∆𝑃𝐷 𝑖 = 1, … … , 𝑁 𝐸𝑐 − 23

Ejemplo_3

Igual al ejemplo 1, pero considerando un incremento en la demanda de 900

MWs.

7.92 + 0.00312𝑃1 = 𝜆 [$/𝑀𝑊ℎ]

7.85 + 0.00388𝑃2 = 𝜆 [$/𝑀𝑊ℎ]

7.97 + 0.00964𝑃3 = 𝜆 [$/𝑀𝑊ℎ]

Calculando los coeficientes de participación para cada generador:

1
𝐾1 = 7.92 = 0.3330
1 1 1
7.92 + 7.85 + 7.97

1
𝐾1 = 7.85 = 0.3360
1 1 1
7.92 + 7.85 + 7.97

1
𝐾1 = 7.97 = 0.3309
1 1 1
7.92 + 7.85 + 7.97
 Calculando los porcentajes para cada generador:

El porcentaje de incremento de la demanda fue de 50 Mws, por lo tanto:

∆𝑃1 = 0.3330 ∗ 50 = 16.65

∆𝑃2 = 0.3360 ∗ 50 = 16.80
∆𝑃3 = 0.3309 ∗ 50 = 16.54

De donde:

𝑃1 = 393.1698 + 16.65 = 𝟒𝟎𝟗. 𝟖𝟏𝟗𝟖 𝑴𝒘𝒔

𝑃2 = 334.6038 + 16.80 = 𝟑𝟓𝟏. 𝟒𝟎𝟑𝟖 𝑴𝒘𝒔
𝑃3 = 122.2264 + 16.54 = 𝟏𝟑𝟖. 𝟕𝟔𝟔𝟒 𝑴𝒘𝒔

Ejemplo_4

Igual al ejemplo 2, considerando un incremento en la demanda de 900 MWs.

7.85 + 0.00388𝑃2 = 𝜆

7.97 + 0.00964𝑃3 = 𝜆

Calculando los coeficientes de participación para cada generador:

1
𝐾2 = 7.85 = 0.5038
1 1
+
7.85 7.97

1
𝐾3 = 7.97 = 0.4962
1 1
+
7.85 7.97

Calculando los porcentajes para cada generador:

El porcentaje de incremento de la demanda fue de 50 Mws, por lo tanto:

∆𝑃2 = 0.5038 ∗ 50 = 25.19

∆𝑃3 = 0.4962 ∗ 50 = 24.81
De donde:

𝑃2 = 187.1 + 25.19 = 𝟐𝟏𝟐. 𝟐𝟗 𝑴𝒘𝒔

𝑃3 = 62.9 + 24.81 = 𝟖𝟕. 𝟕𝟏 𝑴𝒘𝒔

2.3.- DESPACHO ECONÓMICO TÉRMICO CON PÉRDIDAS EN LA

TRANSMISIÓN