SEGUNDA UNIDAD DIDÁCTICA:

DESPACHO ECONÓMICO DE CARGA Y FLUJO DE

POTENCIA OPTMO
 DESPACHO ECONOMICO DE CARGA
1. SISTEMAS DE POTENCIA INTERCONECTADOS

INTRODUCCIÓN
La función de un sistema eléctrico de potencia (SEP), es satisfacer la demanda de los consumidores
respetando los siguientes requerimientos básicos:
 - confiabilidad y seguridad en el suministro
- mínimo costo
- Calidad de producto (frecuencia y niveles de voltaje)
- Optimizar el uso de los recursos energéticos disponibles
- Preservación del medio ambiente
El suministro de la energía demandada con mínimo costo es un requerimiento económico básico para
los sistemas eléctricos, que procura optimizar el uso de los recursos energéticos disponibles. Los
requerimientos de calidad, confiabilidad y seguridad en el suministro, así como la preservación del
medio ambiente, constituyen restricciones que en general encarecen el suministro de energía.
 INTRODUCCIÓN

• Para satisfacer en forma económica la energía y potencia requerida por los consumidores, se
requiere la realización de múltiples tareas de planificación. La tarea global comprende tanto la
planificación de la expansión de la generación, la red de transmisión y la demanda, como la
operación óptima del sistema.

Planificación de la expansión

Operación del sistema eléctrico

Fig, 1.1. Tarea de un SEP

• La planificación de la expansión del sistema, comprende la incorporación de nuevos componentes

como las modificaciones y/o reemplazo de componentes obsoletos. La potencia máxima que será
requerida por los consumidores representa un indicador básico para una adecuada planificación del
crecimiento del sistema. Debe preverse un sistema de generación que permita satisfacer la
demanda máxima, mas una cantidad de reserva suficiente para garantizar niveles mínimos de
confiabilidad
 INTRODUCCION
• La operación de un sistema eléctrico se realiza mediante los centros de despacho de carga y es
conveniente su análisis en dos etapas: programación de la operación, y control de la operación.
• En la etapa de programación, como su nombre indica, se planifica o programa la operación
futura del sistema eléctrico, mientras que, la etapa de control tiene que ver con la operación en
tiempo real del sistema, ambas requieren de la ayuda de computadoras y programas de
aplicación a sistemas eléctricos de potencia.

• Figura 1.2 Esquema de optimizaciones secuenciale

 INTRODUCCIÓN
• En la figura 1.2, se observa la relación entre las diferentes fases de la programación de la generación.
Dado el grado de acoplamiento que existe entre los diferentes niveles de solución, los resultados de
cada subproblema son utilizados como datos de entrada y/o restricciones para la solución del siguiente
subproblema.
• Los resultados más importantes que se esperan de la programación de mediano plazo (4 años) son:
estrategia de movimiento de agua en los embalses, previsiones de combustibles, plan de
mantenimientos preventivos del parque de generación, remuneración del transporte. Asimismo, brinda
el marco de referencia para el cálculo de la programación de corto plazo, en el cual se consideran con
mayor exactitud las restricciones impuestas al parque de generación y la red de transmisión.
• Como resultado de la programación de la operación de corto plazo, se define el predespacho y
el despacho de las unidades que participan en la operación, sujetas a restricciones técnicas y
operativas incluidas la red de transmisión, seguridad y otros usos del agua con nivel horario.
Se optimiza en forma conjunta el despacho semanal y el diario considerando la semana como
horizonte, para coordinar con los resultados de mediano plazo.
• En el despacho de tiempo real (cada hora, minuto) se define el despacho de activo y reactivo,
y la regulación de frecuencia y voltaje para el seguimiento de las variaciones temporales de la
demanda
 CENTROS DE CONTROL DE DESPACHO DE CARGA
Los términos Centros de Despacho de Carga y Centros de Control (CDC) han sido usados
indistintamente como sinónimos.
En general el control se realiza a través de una estructura jerárquica (Figura 1.3), en cuyo primer
nivel se ubica el centro de control nacional (CCN), el cual es responsable de la toma de todas las
decisiones relativas a la operación segura y económica del sistema, es responsable además del
control de la frecuencia, del control de las transferencias a través de las interconexiones y de la toma
de decisiones correctivas en casos de emergencia
En el segundo nivel se ubican los centros de control regionales (CCR), los cuales monitorean la
operación de las redes regionales bajo control, toman decisiones con influencia local y supervisan la
apertura y cierre de los interruptores. Además envían estimaciones de la demanda y la disponibilidad
de las unidades de generación al CCN para las tareas de planificación. En el tercer nivel se ubican las
subestaciones, centrales, etc.
Figura 1.3: Estructura jerárquica
de los centros de control
3. SISTEMAS DE ADMINISTRACION DE ENERGIA (EMS)
•La operación y control de un sistema eléctrico de potencia es una tarea compleja, que requiere de muchas
evaluaciones y toma rápida de decisiones. Algunas funciones operativas pueden basarse en simples
procedimientos, otras requieren del empleo de complejos modelos matemáticos, la mayoría de las funciones se
relacionan entre sí o son dependientes unas de otras,
•Un sistema EMS (Energy Manangement Systems) o Sistema de Administración de Energía tiene las siguientes
funciones de aplicación:
•Base de datos
•control supervisión y adquisición de datos (SCADA)
•pronostico de la demanda
•estimador de estados
•valoración de la seguridad
•análisis de la red en tiempo real
•despacho económico
•asignación de unidades (unit commitment)
•control automático de generación
 3. SISTEMAS DE ADMINISTRACION DE ENERGIA (EMS)
• Un sistema EMS es un sistema integrado muy complejo, incluye todos los niveles que intervienen
en su funcionamiento, desde la subestación, el centro de control regional hasta el centro de control
nacional. ORGANIZACION ESTRUCTURAL DE UN SISTEMA DE CONTROL EN LINEA O EMS

Mediciones y Telecomunicaciones (SCADA)

Estimador de Estado

Base de Datos

Analisis de Pronostico Despacho Control de la

Seguridad demanda Economico Generación

Cambios en Predespacho Distribución

el Sistema de Generacion entre Unidades

Operadores

Displays
4. SISTEMAS DE POTENCIA INTERCONECTADOS
En las primeras etapas de los sistemas de potencia, los sistemas de suministro eran pequeños y su
operación era simple, pero con el crecimiento de la demanda eléctrica, los sistemas de potencia se
fueron expandiendo. Muchos de los grandes sistemas de potencia se fueron interconectando
obteniendo así ventajas técnico-económicas inherentes de la operación coordinada de sistemas de
potencia. Los beneficios de la operación de sistemas de potencia interconectados son las siguientes:
4.1 REDUCCIÓN EN LA CAPACIDAD DE GENERACIÓN DEBIDO A LA DIVERSIDAD
DE LA DEMANDA
La demanda pico diaria (mensual/anual) de dos sistemas de potencia generalmente no ocurre durante
la misma hora del día. Como resultado de esta diversidad, la misma carga del sistema combinado en
un día podría ser menor que la suma de la máxima demanda de las máximas demandas de los
sistemas individuales en el mismo día.
4.2 REDUCCIÓN DE LA CAPACIDAD EN STAND BY
En orden a asegurar el suministro de potencia a los consumidores, la reserva de potencia tiene que ser
mantenida por encima de la capacidad requerida para satisfacer la demanda pico anticipada. Estas
reservas son requeridas para tomar en cuenta el mantenimiento programado de plantas de generación,
la salida de la unidad más grande de generación e incertidumbre del pronóstico de la demanda.
Compartiendo la reserva a través de las interconexiones, puede reducirse la reserva combinada ante salidas forzadas
de generación
 4. SISTEMAS DE POTENCIA INTERCONECTADOS
4.3 AUMENTO EN EL TAMAÑO DE LAS UNIDADES DE GENERACIÓN
Con la operación de sistemas interconectados, plantas de generación más grandes con unidades más grandes, con
costos de capital y operación más bajos, pueden ser localizados en sitios más económicos de una región como un
todo.
4.4 UTILIZACIÓN OPTIMA DE LA CAPACIDAD DISPONIBLE DE LAS PLANTAS Y
FACILIDADES DE TRANSMISIÓN
Con la operación interconectada, las plantas térmicas pueden ser usadas más efectivamente junto con las plantas
hidráulicas. Durante periodos lluviosos, las plantas hidráulicas pueden ser operadas con factores de carga más altos
debido a la disponibilidad extra de agua, así contribuyendo más energía al sistema, cuando alguna de las unidades
térmicas podría ser desconectada.
En periodos secos, las unidades térmicas podrían estar trabajando con un factor de carga tal que el agua en los
reservorios de las plantas hidráulicas podría ser conservada.
4.5 CONFIABILIDAD DE SUMINISTRO DE POTENCIA
Los sistemas contiguos pueden proporcionarse asistencia mutua en caso de emergencia. Emergencias de corta
duración ocurren en cada sistema debido a la salida de plantas de generación, equipos de subestaciones y líneas de
transmisión. Son numerosas las oportunidades de asistencia en emergencias proporcionadas por los sistemas
contiguos.
4.6 MEJORANDO LA FRECUENCIA
En un sistema que opera independientemente, la frecuencia fluctúa continuamente en pequeñas amplitudes; en
cambio cuando el mismo sistema esta Interconectado a un gran sistema, la frecuencia del sistema es más estable.
 
 5. PROGRAMACIÓN DE LA GENERACIÓN
• PLANIFICACIÓN DE LA OPERACIÓN
• Planificar la operación de los SEP tiene por finalidad determinar una política
de operación en los diferentes horizontes de estudio, con la finalidad de
abastecer la demanda de energía eléctrica con determinados niveles de
confiabilidad, calidad de servicio y bajo ciertas condiciones operativas, todo
ello al mínimo costo operativo posible.
•  Horizontes de Planificación de la Operación
• La planificación de la operación tiene como objetivo establecer un
programa de generación, abarcando diversos horizontes de estudio, que
mediante el uso de métodos de optimización dan como resultado procesos
de decisión secuenciales. Según Uribe (2000), los resultados de
planificación de mayor horizonte constituyen un marco referencial para
decisiones de más corto plazo.
• Usualmente para planificación operativa se define los siguientes horizontes
de estudio:
 PLANIFICACIÓN DE LA OPERACIÓN
Largo Plazo (LP): Corresponde a un horizonte multianual, en el que se decide la operación de
embalses con capacidad de regulación interanual y eventuales contratos de largo plazo de compra de
combustible y/o suministro de energía eléctrica de sistemas vecinos, la resolución usualmente es
semestral o anual. La demanda se representa de forma simplificada y mediante un análisis
probabilístico se establece la disponibilidad de unidades de generación e hidrologías (Uribe, 2000, p.
10). Horizonte de planificación utilizado en estudios de planificación de obras.
Mediano Plazo (MP): Comprende un horizonte multianual de menor duración que el horizonte de
largo plazo, con resolución usualmente mensual, en el cual se decide las estrategias de transferencias
estacionales entre embalses (Uribe, 2000, p. 11). Al igual que el largo plazo se realiza un análisis
probabilístico para la disponibilidad de unidades e hidrologías. Como resultado se obtiene costos
marginales y el valor marginal del agua, valor utilizado es estudios de menor horizonte.
Corto Plazo (CP): Comprende un horizonte semanal o mensual con resolución horaria, en este nivel
se considera conocidas las hidrologías, la disponibilidad de unidades y el detalle de demanda (Uribe,
2000, p. 11). Se busca determinar el programa de asignación de unidades que permita establecer la
seguridad de la operación mediante análisis de contingencia haciendo uso de reserva rotante y/o fría,
asimismo debe considerar el efecto de las incertidumbres de las decisiones de largo plazo.
Operación de Tiempo Real (OTR): Se establece el despacho económico realizando ajuste por
variaciones instantáneas de demanda y/o fallas imprevistas en la oferta, considerando limitaciones de
disponibilidad en el despacho de unidades y criterios de seguridad previamente establecidos (Uribe,
2000, p. 10).
 PLANIFICACIÓN DE LA OPERACIÓN

Figura 2.1 Horizontes de Planificación Operativa

Programa de Obras de nuevas centrales.

Manejo de Embalse de regulación anual.

Programación de la Operación
de Largo Plazo (Anual) Contratos de Combustible de Largo Plazo.

Mantenimiento Mayor de las centrales.

Despacho de centrales Hidráulicas.

Programación de la Operación
Manejo de reservorios de regulación diaria.
de Corto Plazo (Semanal)
Selección de Parque Térmico (Predespacho)

Programación de Operación Distribución de cargas entre las maquinas en

Diaria funcionamiento (Despacho)

Ajuste de la demanda ante variaciones de

Operación de Tiempo Real
demanda y/o fallas imprevistas
 5. PROGRAMACIÓN DE LA GENERACIÓN

INTRODUCCION
En este capítulo se realiza una introducción al problema de la programación de generación de las
centrales eléctricas. Con ese propósito, se revisan conceptos básicos asociados a los costos de
generación y se plantean dos preguntas importantes:
¿Qué unidades generadoras van a estar en servicio en cada momento?
¿Cuánto debe generar cada unidad para minimizar los costos?
Este problema, fundamental para el correcto funcionamiento de un sistema eléctrico teniendo en
cuenta que en cada momento se debe atender la demanda más las perdidas en el sistema, es
objeto de distintos planteamientos dependiendo del horizonte temporal y del objetivo concreto
perseguido:
• Largo Plazo
• Planificación de mantenimientos
• Planificación de los recursos de generación térmica e hidráulica
• Corto-medio plazo
• Planificación de arranques y paradas de unidades térmicas
• Coordinación Hidrotermica a corto plazo
• Despacho económico
 Actividades del Operador del Sistema

Las actividades del operador del sistema están divididas en:

 Pre-despacho
 Despacho
 Despacho instantáneo
Actividades del Operador del Sistema

Pre-despacho: actividad de planificación

 Puede abarcar un período de una semana a un día antes de la

operación actual

 El operador dispone de la predicción horaria de la demanda

agregada del sistema, base para definir el plan de operación de
las unidades en el período.

 Se incluye en el plan la disponibilidad del recurso hidráulico del

sistema, dependiendo de su nivel y otros factores.

 El plan incluye además los intercambios de potencia con otras

regiones
Actividades del Operador del Sistema

Despacho: programa de corto plazo

 El estado de despacho comienza alrededor de 30 minutos antes

de la operación actual.
 En este estado el operador lleva a cabo las actividades
programadas de corto plazo, tal como: flujo de potencia,
despacho económico de carga.
 En esta etapa el operador tiene una clara idea de la demanda
esperada en el sistema, incluyendo los intercambios de potencia
programados.
 El operador decide sobre los soportes de reactivos (switching) y
desprendimientos posibles de cargas.
 El operador es responsable de mantener la adecuada reserva en
el sistema
 Actividades del Operador del Sistema

Despacho Instantáneo: actividad en tiempo real

 Estado con alrededor de 5 minutos antes de la operación en tiempo real actual

 El operador lleva acabo los planes punteados en los estados de pre-despacho y
despacho, basados en los datos de tiempo real.
 Si es necesario inicia las actividades de control secundario de frecuencia.
 Actualiza el programa de generación basado en los factores de participación de
las unidades y la información de carga.
 Software de Apoyo

 Gran número de algoritmos y programas para la simulación y análisis están

disponibles para ayudar al operador en la toma de decisión y hacer la operación
del sistema más segura y económica.

 Dependiendo del grado de aproximación estos programas están disponible para

el pre-despacho
 5. PROGRAMACIÓN DE GENERACIÓN
En la anterior clasificación, el “despacho económico clásico” aparece como la herramienta
básica que permite introducir una serie de conceptos utilizados en el resto de los problemas.
Por ello, se estudiara con más detalle el “despacho económico clásico”, luego, a continuación
se realiza una breve introducción a los problemas de la programación a corto plazo de las
unidades térmicas (problema conocido como “unit commitment” o planificación de arranques y
paradas) y de la utilización de los recursos hidráulicos o coordinación Hidrotérmica.

DESPACHO ECONOMICO CLASICO

El problema del “despacho económico clásico” consiste en determinar la potencia que debe
suministrar cada unidad generadora en servicio para una demanda determinada PD, con el
objetivo de minimizar el costo total de generación. Para ello, es necesario conocer los costos
variables de los combustibles, los rendimientos térmicos de las unidades, la red de
transmisión, etc.
 Despacho Económico de Carga

 El DEC se ejecuta en la etapa de despacho

 Primariamente busca la repartición de la carga total del sistema
entre las unidades generadoras que están en línea, manteniendo el
costo total de operación en el mínimo.
 Generalmente, el DEC se ejecuta cada 5 minutos
 El modelo del DEC debe representar al sistema en el mayor detalle
posible
El Despacho Económico de Carga
 Considerar un sistema con N unidades de generación térmica,
cada unidad entrega su potencia de generación, Pi, al sistema
de transmisión, por el cual se abastece la carga total del
sistema, PD

 Fig. 2.1 Modelo de barra única
 Modelo de Barra Múltiple

En este modelo, se consideran que las barras, generadores y las cargas
del sistema, están
interconectados por líneas y transformadores. Se puede decir que éste 
modelo, es el más real, pero tiene
sus limitantes como el excesivo costo computacional por ser multidime
nsional y multietapa. En la Fig. se representa esquemáticamente el mod
elo multinodal
 El Despacho Económico de Carga

 Las pérdidas de potencia en el sistema de transmisión

son una función de las impedancias de las líneas y las
corrientes.

 Las corrientes pueden también expresarse en

términos de la variable independiente P. Así,
finalmente las perdidas pueden expresarse mediante
Pperd = f(P)
 DESPACHO ECONOMICO DE CARGA

CARACTERISTICAS DE LAS UNIDADES GENERADORAS

La descripción de una unidad térmica -generadora comienza con la especificación de la
cantidad de calor de entrada requerida para producir una cantidad de energía eléctrica
como salida.
Así, la característica Entrada–Salida de la unidad-generadora, tiene forma cuadrática -
convexa, como en la figura 2.1. En el eje de ordenadas esta la entrada de calor H
[Btu/h] y en el eje de abscisas, la potencia de salida P [kW].
Así, la función cantidad de calor H es igual a la siguiente expresión:
H = a + b Pg + c Pg2 [Btu/h]
Multiplicando la cantidad de calor H por el costo de combustible se obtiene la función
costo de combustible C [$UM/h]. El costo total de producción incluye el costo de
combustible, el consumo propio y el costo de operación - mantenimiento.
Se asume que esos costos son un valor o porcentaje fijo del costo de combustible y
generalmente se incluyen en la curva costo de combustible.
 DESPACHO ECONOMICO DE CARGA
• Esta información se obtiene, a partir de pruebas que se realizan al grupo turbina-generador, para
varios niveles de potencia de salida (100%, 75% y 50%). La tasa de calor o Heat Rate (HR), se
define como la relación entre la entrada de calor en Btu/h dividido por la potencia de salida en
kW. H [Btu/h]
F [$us/h]
HR = H/P [Btu/kWh] Costo
Ci(UM/h)
Fig. 2.1 Característica Entrada - Salida
F
P

Pmin Pmax SALIDA

P (KW)

• El Heat Rate ( el índice de calor) es el reciproco de la HR

[Btu/Kwh]
eficiencia o rendimiento. En la figura 2.2, (costo de calor),
se observa que la máxima eficiencia
HR min
de la unidad se obtiene en el mínimo de la función
HR, que se da para valores próximos a la potencia
máxima.
Pmin. P ef. Pmax. P (KW)
 DEC
• El Costo Incremental de Combustible (IC) es igual a la derivada de la función costo.
dH/dPg= IC = dC/dPg = b + 2c Pg [$UM/kWh]
• El Costo Medio de Producción es igual a la división de la función costo total de producción por la
potencia máxima de salida. Es decir: Costo Medio = CT/Pgmáx [$US/kWh]
CÁLCULO DEL COSTO DE CALOR
Una información importante, para el cálculo de las funciones costo es el dato del valor de calor
(Heat Rate) de la turbina, determinada en sitio, a partir de pruebas efectuadas al grupo turbinas a
gas-generador. En la figura 2.3, se observa, que los datos a ser tomados durante las pruebas son las
siguientes: temperatura del aire de entrada al filtro de aire de la turbina (temperatura ambiente),
presión atmosférica en el sitio, volumen de gas que ingresa a la cámara de combustión, potencia y
energía activa de salida del generador, medida en bornes, etc.

Fig. 2.3 Esquema de medición de la prueba

CÁLCULO DE LA FUNCIÓN COSTO DE COMBUSTIBLE
La función costo de combustible (F), se determina a partir de las pruebas antes mencionadas, con la siguiente
información:
• Temperatura ambiente en [ºC]
• Presión atmosférica del sitio en [mbar]
• Poder calorífico inferior del gas [Btu/PC]
• Costo del combustible en [$UM/Btu]
• Potencia de salida del generador en [kWh]
• Cantidad de calor (Heat Rate) en [Btu/kWh] para tres estados de operación de la maquina, que son 100%, 75%
y 50% de carga.
La función consumo de combustible generalmente se representa como una función convexa cuadrática, de la
forma,
Hi = ai + bi PGi +ci PGi2
El consumo de calor o rendimiento térmico (Heat Rate), fue antes definido de la siguiente manera.
HRi = Hi / PGi
Luego, igualando con la expresión del consumo de combustible se obtiene,
Hi = HRi x PGi = ai + bi PGi + ci PGi2
En esta ecuación cuadrática, son conocidos los rendimientos térmicos para los tres estados de carga
mencionados y las potencias de salida respectivas, siendo solo incógnitas los coeficientes de la función (ai, bi, ci).
DEC: DESPACHO UNINODAL DE UNIDADES TÉRMICAS
• Dado un sistema uninodal de N unidades térmicas de generación, conectadas a
una barra simple y suministrando energía eléctrica a una carga PC.
• En este análisis se considera un sistema eléctrico sin pérdidas de transmisión y
en la figura 2.4, se observa el sistema uninodal, siendo F1=C1 la función costo y
PG1 la potencia eléctrica de salida del generador.

Figura 2.4 Despacho económico uninodal: Fi(Pi)

 Despacho uninodal
• El costo total de operación del sistema uninodal es igual a la suma de
los costos de operación de cada unidad de generación en [$US/h].
N
CT   Ci
i 1

• La restricción principal del sistema a ser cumplida es que la generación

sea igual a la demanda N
P
i 1
Gi  PC

• Entonces el despacho económico consistirá en encontrar el costo

mínimo de operación del sistema, resolviendo un problema de
optimización planteado de la siguiente manera:
N
min Ct   Ci Función objetivo
i 1

N s.a. Restricción de igualdad

Pc  P
i 1
i
 Despacho uninodal
Este es un problema de optimización con restricciones que para su mejor comprensión será resuelto por el “Método
de los operadores de Lagrange”. Definimos la función de Lagrange:

L  CT  
Siendo,  el operador de Lagrange y conocido también como el costo incremental de las unidades de generación, y la
función  puede ser definido de la siguiente manera: N
   Pi  Pc  0
i 1

dL dCi
   0
La condición para el mínimo de la función es que: 0 dPi
dPi

dCi

Luego se obtiene la “Ecuación de Coordinación” y es: dPi
Esta es la condición necesaria para la existencia de un punto de operación a mínimo costo para el sistema
térmico de generación, es decir, que los costos incrementales (lambda) de todas las unidades sean iguales.
Conclusión: El costo mínimo de operación se tendrá cuando todas las unidades funcionan con el mismo
costo incremental y se cumple el balance de potencia dado por la ecuación (2).
 DESPACHO UNINODAL
• Condiciones de Kuhn Tucker
Al problema de optimización debe agregarse los límites operativos de las unidades de generación. Es decir la
potencia de salida de cada unidad debe ser mayor o igual que la potencia mínima permitida (por ejemplo el
limite técnico) y menor o igual que la potencia máxima permitida en la unidad de generación. Luego el
problema de optimización del despacho económico, se plantea de la siguiente manera.
N
min Ct   Ci
i 1
Función objetivo
N s.a.
Restricción de igualdad
Pc   Pi
i 1 Restricción de desigualdad

P i
min
 Pi  P i
max

• Las condiciones de Kuhn Tucker complementan las condiciones del Lagrangiano para incluir las
restricciones de desigualdad como términos adicionales. Así, las condiciones necesarias para el
despacho económico sin perdidas son:
 Despacho uninodal
Así, las condiciones necesarias para el despacho económico sin perdidas
son:

d Ci
d Pi
 para
Pi
min
 Pi  P i
max

d Ci
d Pi
 
para Pi  P i
max

d Ci
d Pi
 para Pi  P i
min
• Operación económica de los Sistemas de Potencia:
Importante para recuperar y obtener beneficios del
capital que se invierte.
• Involucra la generación de potencia y el suministro
– Despacho económico: Costo mínimo de producción
– Suministro: Con pérdidas mínimas a las cargas
– Tarifas y conservación del combustible: Máxima
eficiencia
– Minimización del costo del kWh
• Aproximación clásica al despacho económico.
Despacho económico sin considerar pérdidas de la red
• Formulación del problema
Dado un sistema con n nudos y m generadores y todas las
potencias demandadas por las cargas SDi (i=1,2,...,n),
determinar la potencia activa que debe generar cada
generador Pi, cuyo costo de operación es Ci(Pi), con
i=1,2,......,m, para minimizar el costo total CT; es decir:
m
C T   C i (Pi )  C1 (P1 )  C 2 (P2 )  ........  C m ( Pm ) (1)
i 1

Sujeto a:
m
PD   Pi  P1 P2  .......  Pm (2)
i 1

Pi mín
 Pi  Pi máx
con i  1,2,..., m (3)
• Solución sin considerar límites de generación
– Costo incremental Marginal de la unidad i CIi: Costo
adicional por hora que tiene el aumentar la salida de la
máquina i en 1 MW. Representa la pendiente de la curva de
costo Ci de la Figura 1; es decir:
Entrada
Hi(Mbtu/h) CIi(Pi)
dCi (Pi )
Costo CI i  (4)
Ci(UM/h) dPi

Pimín Pimáx Salida, Pi (MW)

Figura 1.- Característica E/S de un grupo T-G

En este caso, sólo se deben resolver las ecuaciones (1) y (2),
lo que se puede plantear de la siguiente forma:
– Valor mínimo de CT: Se tiene cuando dCT=0; es decir:
C T C T C T
dC T  dP1  dP2  .......... dPm  0
P1 P2 Pm
Ci sólo depende de la potencia Pi generada por ella:
dC1 dC 2 dC m
dC T  dP1  dP2  ..........  dPm  0 (5)
dP1 dP2 dPm
Si PD se supone ctte., su diferencial será cero:
dPD  dP1  dP2  .......  dPm  0 (6)
Multiplicando (6) por  (multiplicador de Lagrange) y
restando el resultado a (5) se obtiene:
 dC1   dC 2   dC m 
   dP1     dP2  ..........     dPm  0
 dP1   dP2   dPm 
 La ecuación anterior se satisface cuando:
dC1 dC 2 dC m
  ..........    (7)
dP1 dP2 dPm
– Conclusión: El costo mínimo de operación se tendrá
cuando todas las unidades funcionan con el mismo costo
incremental y se cumple el balance de potencia dado por la
ecuación (2).
Ejemplo 1: Las curvas de costo de dos generadores son:
C1 ( P1 )  900  45P1  0,01P12 ; C 2 ( P2 )  2.500  43P2  0,003P22
La carga total PD que debe ser suministrada es de 700 MW.
Determine la potencia que debe entregar cada máquina, el
costo incremental y el costo total.
Solución analítica: Se debe cumplir:
  CI1  45  0,02P1
  CI 2  43  0,006P2
P1  P2  700
Solución: P1=84,6 MW; P2=615,4 MW;
CI1=CI2== 46,69 UM/MWh; CT=34876,96 UM/h
C:\Archivos de programa\PowerWorld\Simulator\[Link]
Costos CI1(P1)
Incrementales
(UM/MWh)
CI2(P2)
46,69
45
43

84,6 615,4 P (MW)

Figura 2.- Solución gráfica Ejemplo 1

El Método Iteración-Lambda
Para entender como funciona este algoritmo, suponer
que se tiene un sistema con tres unidades, y se desea
encontrar el punto de operación óptimo.
Gráficamente la idea es dibujar las características de
CI de cada unidad i, y para cada valor de λ encontrar
la potencia generada por cada Unidad Generadora, y
la potencia total que debe equilibrar la demanda

Gráfica de un SEP con tres GENERADORES

La Figura siguiente muestra un diagrama en bloque del
algoritmo de cálculo Iteración-Lambda, que encuentra la
solución del despacho económico térmico sin considerar
pérdidas en las LT.
 Método de iteración en 
Paso 1: Elegir un valor inicial de 
Paso 2: Hallar las correspondientes potencias de los generadores, P1,
P2, ....,Pm
  m
 Pi  PD  0
Paso 3: Si i 1 aumentar  y volver al paso 2.
m
 Pi  PD  0
Si i 1 reducir  y volver al paso 2.
m
 Pi  PD  
Si i 1 proceso terminado.

El valor de  es un valor positivo que corresponde a la tolerancia

aceptable para la solución.
• Solución considerando límites de potencias generadas
Costos
Incrementales CI1(P1) CI2(P2)
3
(UM/MWh) 
2

1
CI3(P3)

P1mín P2mín P3mín P1máx P2máx P3máx P (MW)

Figura 3.- Curvas de CI con límites de Potencia generada

– Método de Iteración en : Suponer que todas las
unidades operan al mismo costo incremental. Si una o
varias unidades alcanzan uno de sus límites ya sea inferior
o superior, fijar la potencia a generar en el límite alcanzado
y continuar el proceso de ajuste de  con el resto de las
Ejemplo 2: Considere un rango de valores posibles para PD de
100 a 800 MW con las unidades del del Ejemplo 1, sujetas a los
límites: 50 MWP1200 MW; 50 MWP2600 MW. Represente
P1 y P2 en función de la potencia demandada PD.
Tabla 1.- Solución para el Ejemplo 2
PD (MW) P1 (MW) P2 (MW) CI (UM/MWh)  CT (UM/h)
100 50 50 43,30 dC2/dP2 7.832,5
200 50 150 43,90 dC2/dP2 12.192,5
300 50 250 44,50 dC2/dP2 16.612,5
400 50 350 45,10 dC2/dP2 21.092,5
500 50 450 45,70 dC2/dP2 25.632,5
550 50 500 46,00 1=2 27.925,0
600 61,54 538,46 46,23 1=2 30.230,8
680 80 600 46,60 1=2 33.944,0
700 100 600 47,00 dC1/dP1 34.880,0
 CI (UM/MWh)

50
49
48 CI1

46,6 47
46
CI2
45
44
43,3
43
42
0 50 100 200 300 400 500 600 P1 ; P 2 (MW)
80

Figura 4.- CI de generadores del Ejemplo 2

 P1; P2 (MW)

600

500

400 P2

300

200
100 P1
80
50

0 100 200 300 400 500 600 700 800 PD= P1+ P2 (MW)
550 680

Figura 5.- Reparto de carga entre los generadores

 Ejemplo 2
 

Dado un sistema uninodal de 3 unidades alimentando a una demanda de 52 MW,

determine a) la potencia de salida de cada unidad para obtener el despacho
económico, b) el costo de operación de cada unidad, c) el costo de operación total
del sistema y d) el Costo Incremental del sistema. Datos
Unidad H [Mbtu/h] Pmin [MW] Pmax [MW]

1 H1 = 0.234 P12 +2.112 P1 +112.8 12.63 21.05

2 H2 = 0.0022 P22 +8.71 P2 +69.91 11.29 19.50

3 H3 = 0.1032 P32 +6.119 P3 +79.61 12.16 20.23

Donde H es la entrada de calor a la turbina y el costo del combustible es la potencia

de salida de cada máquina es   1,76 ($US / MBtu )
 EJEMPLO 2
Solución
Multiplicando la entrada de calor por el costo del combustible obtenemos las funciones de costo de las
unidades en $US/h.
C1 = 0.412 P12 + 3.72 P1 +198.5 ($US/h)
C2 = 0.0039 P22 + 15.33 P2 +122 ($US/h)
C3 = 0.182 P32 + 10.77 P3 +140.11 ($US/h)
Derivando respecto a la potencia obtenemos el costo incremental de combustible, aplicando la
ecuación de coordinación y la ecuación de balance de potencia, resolvemos el sistema de ecuaciones
siguiente.
dC1/dP1 = 0.824 P1 + 3.72 = 
dC2/dP2 = 0.0078 P2 + 15.33 = 
dC3/dP3 = 0.364 P3 + 10.77 = 
P1 + P2 + P3 = 52 [MW]
La solución del sistema es:
a. La potencia de salida para el despacho económico de unidades
P1  15.92  MW 
P2  19.39  MW 
P3  16.98  MW 
 EJEMPLO 2
Se observa que todas las soluciones están dentro los rangos de
operación establecidos para cada unidad en la tabla de datos.
El costo de operación de cada unidad : CT= C1+ C2+ C3 = 1167.43 [$US/h]
C1 = 362.14 [$US/h]
C2 = 434.91 [$US/h]
C3 = 370.38 [$US/h]
El costo total de operación del sistema
CT  C1  C 2  C 3  1167 .43 [$US/h]

El costo incremental de combustible

 = 16.84 [$US/MWh]
 EJEMPLO 2
• Solución grafica del problema

[$ U S /M W h ] 17.80
17.40
17.00
16.60
16.20
15.80
15.40
15.00
11 12 13 14 15 16 17 18 19 20 21

P [MW]

dF1/dP1 dF2/dP2 dF3/dP3

Figura 2.5 Costos incrementales iguales

 P: 03 Ejemplo incluyendo límites en la generación :
El costo total de tres plantas térmicas en $/h está dada C 1  500  5.3 P1  0.004 P12
por:
C 2  400  5.5 P2  0.006 P22
C 3  200  5.8 P3  0.009 P32

Donde P1, P2 y P3 están en MW. La demanda total PD es 975 MW.

Los límites de generación son:

200  P1  450
150  P2  350
100  P3  225

Sin considerar pérdidas, encontrar el despacho óptimo.

5. 3 5.5 5. 8
975   
  0.008 0.012 0.018  9.16 $/MWh
1 1 1
 
0.008 0.012 0.018
9.16  5.3
P1   483 MW
2( 0.004)
9.16  5.5
P2   305 MW
2(0.006)
9.16  5.8
P3   187 MW
2( 0.009)
P1 viola el límite de Pmax, por lo que la “pego” al tope de 450MW y redespacho las otras dos
5 .5 5.8
525  
  0.012 0.018  9.4 $/MWh
1 1

0.012 0.018

9.4  5.5
P2   325 MW
2( 0.006 )

9.4  5.8
P3   200 MW
2( 0.009 )
 Despacho económico considerando pérdidas de la red
• Pérdidas del sistema Pp como función de Pi; es decir:
Pp  Pp (P1 , P2 ,......, Pm )
• Formulación del problema: Se trata de encontrar las
potencias Pi para minimizar el costo total CT (F)
m
C T   C i (Pi )  C1 (P1 )  C 2 (P2 )  ........  C m ( Pm ) (8)
i 1
m
 Pi  Pp (Pi )  PD 0 (9)
i 1

Pi mín
 Pi  Pi máx
con i  1,2,..., m (10)
• Solución sin considerar los límites de generación:
En este caso se obtiene:
 dC i Pp 
       0 i  1,2,.....,m
 dPi Pi 
O bien:
 
 1  dCi
   P  i  1,2,.....,m
 1  p  dPi
 Pi 
Donde:
Pp
Pi :Pérdidas incrementales de transmisión de la máquina i
1
Li  Pp
: Factor de Penalización de la máquina i
1  P
i
El sistema de ecuaciones a resolver queda:
  L i CI i i  1,2,.....m
m
(11)
 Pi  Pp (Pi )  PD  0
i 1

• Solución considerando los límites de generación

En este caso se tiene una solución análoga a la del caso
sin pérdidas; es decir, se deben hacer funcionar todos los
generadores que estén dentro de sus límites de tal forma
que se cumpla que LiCIi=. Si en este proceso una unidad
de generación alcanza uno de sus límites, fijar la potencia a
generar por la unidad en ese límite (máximo o mínimo) y
continuar el proceso de ajuste de  con el resto de las
unidades.
Ejemplo 4.3: Considere las características de costo de los
generadores del Ejemplo 1. La carga total PD que debe ser
suministrada es de 700 MW. La expresión simplificada de las
pérdidas es de la forma: Pp  ( 0,00003P1
2
 0,00009P2
2
) MW
Determine, utilizando el método de iteración en , la potencia
que debe entregar cada máquina, las pérdidas en el sistema, el
costo incremental y el costo total.
Resolución: Las ecuaciones a considerar son:
1
  L 1 CI 1   45  0,02P1 
1  0,00006P1
1
  L 2 CI 2   43  0,006P2 
1  0,00018P2
Pp  (0,00003P12  0,00009P22 )
P1  P2  Pp  700  0
Para utilizar el método de iteración en , conviene escribir
las ecuaciones de la siguiente forma:
  45
P1 
0,02  0,00006
  43
P2 
0,006  0,00018
Pp  (0,00003P12  0,00009P22 )
P  P1  P2  Pp  700  0
Considerando como valor inicial para , el determinado en
ejemplo 1; es decir =46,69, se obtiene:
=50,378 UM/MWh; P1=233,6 MW; P2=489,65 MW;
CT=36.231,5 UM/h
–Comparación con resultados de Ejemplo 1
* Los Costos incremental y total son ahora mayores que en
el caso anterior.
* La potencia entregada por la máquina 1 aumentó y la
entregada por la máquina 2 disminuyó, a pesar de que el
costo incremental de esta última es menor. Ello se debe a
que su Factor de Penalización es mayor.

El proceso iterativo se muestra en la Tabla 2

 Tabla 2.- Iteración en , Ejemplo 3
n (UM/MWh) P1 (MW) P2 (MW) Pp (MW) P (MW ) CT (UM/h)
1 46,690 74,12 256,18 6,07 -375,78 18.002,67
2 50,000 217,39 466,67 21,02 -36,96 34.375,20
3 50,500 238,82 497,02 23,94 11,89 36.830,05
4 50,300 230,25 484,92 22,75 -7,58 35.848,67
5 50,350 232,40 487,95 23,05 -2,70 36.094,09
6 50,370 233,25 489,16 23,17 -0,75 36.192,24
7 50,375 233,47 489,46 23,20 -0,27 36.216,78
8 50,380 233,68 489,77 23,23 0,22 36.241,31
9 50,378 233,60 489,65 23,21 0,03 36.231,5
10 50,377 233,55 489,59 23,21 -0,07 36.226,59

..\[Link]
– Iteración en la potencia de los generadores
* Elegir valores de Pi tales que m P  P  0 ; es decir, en este
 i D
i 1
paso, no se consideran las pérdidas.
* Con los valores de Pi, calcular Li y Pp
* Resolver el sistema de ecuaciones (12), para determinar 
y nuevos valores para Pi
  L i CI i i  1,2,.....m
m (12)
 Pi  Pp (Pi )  PD  0
i 1

* Comparar los nuevos valores de Pi con los anteriores. Si

hay una diferencia importante, volver al paso 2. En caso
contrario, terminar el proceso.
– Otros Métodos: Gradiente de primer y segundo orden.
• Cálculo de las pérdidas en la red
– Método de los Coeficientes de pérdidas B
m m
Pp    B ij Pi Pj (13)
i 1 j1
 B11 B12 .. B1m   P1 
B B 22 .. B 2 m  P 
Pp   P1 P2 .. Pm   21   2    P  t  B  P 
 : : .. :   : 
B Bm2 .. B mm  P 
 m1  m
Los coeficientes Bij no son cttes., varían con la carga del
sistema. Se obtienen a partir de un CFP (caso base) y sus
valores dependen de las condicionbes de carga del sistema.
A partir de (13), con Bij=Bji se puede escribir:
Pp m
 2 B ij Pj (14)
Pi j1
– Formulación para un caso simple (Figura 6)
I1 I2
G1 G2
a b
1 3 2
c I1 + I 2
4 Carga
Figura 6.- Sistema simple con dos generadores
Ra, Rb y Rc : Resistencias de líneas de transmisión a, b y c
Pp: Potencia activa total perdida en el sistema de transmisión
2 2 2
Pp  3 I1 R a  3 I 2 R b  3 I1  I 2 R c
   

Si se supone que las corrientes están en fase:

I1  I 2  I1  I 2
Pp  3I (R a  R c )  6I1 I 2 R c  3I (R b  R c )
2
1
2
2

Sean:
P1; P2: Potencias trifásicas de salida de las máquinas
fp1; fp2: Factores de potencia de salida de las máquinas
V1; V2: Voltajes de línea en barras de los Generadores
I1; I2: Corrientes en la líneas de transmisión
P1 P2
I1  I2 
3 V1 fp1 3 V2 fp 2
Por lo que las pérdidas quedan:
Pp  P B11  2P1 P2 B12  P B 22
1
2
2
2
 Donde: Ra  Rc
B11  2 2
V1 (fp1 )
Rc
B12 
V1 V2 (fp1 )(fp 2 )
Rb  Rc
B 22  2 2
V2 (fp 2 )
Si los voltajes se expresan en kV, las resistencias en )fase,
las unidaes de los coeficientes B quedan expresadas en
MW-1 y Pp en MW.
El cálculo también pueden hacerse en por unidad