DESPACHO ECONÓMICO DE GENERACIÓN TERMOELÉCTRICA USANDO

MÉTODOS DE SENSIBILIDAD AC
 UNIVERSIDAD POLITÉCNICA SALESIANA
SEDE QUITO

CARRERA:
INGENIERÍA ELÉCTRICA

Trabajo de titulación previo a la obtención del título de

INGENIERO ELÉCTRICO

TEMA:
DESPACHO ECONÓMICO DE GENERACIÓN TERMOELÉCTRICA USANDO
MÉTODOS DE SENSIBILIDAD AC

AUTOR:
CRISTIAN MARCELO DÍAZ GARCÍA

DIRECTOR:
CARLOS ANDRÉS BARRERA SINGAÑA

Quito, diciembre 2019

 . Datos de Catalogación Bibliográfica .
Cristian Marcelo Díaz García

DESPACHO ECONÓMICO DE GENERACIÓN TERMOELÉCTRICA USANDO

MÉTODOS DE SENSIBILIDAD AC

Universidad Politécnica Salesiana, Quito – Ecuador 2019

Ingeniería Eléctrica

Breve reseña histórica e información de contacto.

Cristian Marcelo Díaz García (Y’1990 – M’01). Realizó estudios

secundarios en el Instituto Tecnológico Superior Central Técnico
obteniendo el título de Bachiller Técnico industrial en Electrónica,
instituto en el cual posteriormente estudiaría hasta graduarse de
Tecnólogo Electrónico en el año 2012. cdiazg3@est.ups.edu.ec

Dirigido por:
Carlos Andrés Barrera Singaña (Y’1989). Realizó sus estudios
superiores en la Escuela Politécnica Nacional de Quito, donde se graduó
de Ingeniero Eléctrico en octubre de 2013. Además, cursó estudios de
posgrado en la Universidad de Zaragoza, España, donde obtuvo el título
en Máster en Energías Renovables y Eficiencia Energética con
Especialidad en Sistemas Eléctricos. Actualmente es profesor ocasional
a tiempo completo de la Universidad Politécnica Salesiana.
cbarrera@ups.edu.ec

Todos los derechos reservados:

Queda prohibida, salvo excepción prevista en la ley, cualquier forma de reproducción,

distribución, comunicación pública y transformación de esta obra para fines comerciales,
sin contar con la autorización de los titulares de propiedad intelectual. La infracción de
los derechos mencionados puede ser constitutiva de delito contra la propiedad intelectual.
Se permite la libre difusión de este texto con fines académicos o investigativos por
cualquier medio, con la debida notificación a los autores.

DERECHOS RESERVADOS
©2019 Universidad Politécnica Salesiana
QUITO - ECUADOR

I
 DECLARATORIA DE COAUTORÍA DEL DOCENTE TUTOR

Yo, Carlos Andrés Barrera Singaña declaro que bajo mi dirección y asesoría fue
desarrollado el trabajo de titulación “DESPACHO ECONÓMICO DE GENERACIÓN
TERMOELÉCTRICA USANDO MÉTODOS DE SENSIBILIDAD AC” realizado por
Cristian Marcelo Díaz García, obteniendo un producto que cumple con todos los requisitos
estipulados por la Universidad Politécnica Salesiana para ser considerados como trabajo
final de titulación.

Quito D.M., diciembre de 2019

………………………………….
Carlos Andrés Barrera Singaña
C.C.: 0503503336

II
 CESIÓN DE DERECHOS DE AUTOR

Yo, Cristian Marcelo Díaz García, con documento de identificación N° 1723523708,

manifiesto mi voluntad y cedo a la Universidad Politécnica Salesiana la titularidad sobre
los derechos patrimoniales en virtud de que soy autora del trabajo de grado/titulación
intitulado: “DESPACHO ECONÓMICO DE GENERACIÓN TERMOELÉCTRICA
USANDO MÉTODOS DE SENSIBILIDAD AC”, mismo que ha sido desarrollado para
optar por el título de: Ingeniero Eléctrico, en la Universidad Politécnica Salesiana,
quedando la Universidad facultada para ejercer plenamente los derechos cedidos
anteriormente.
En aplicación a lo determinado en la Ley de Propiedad Intelectual, en mi condición de
autor me reservo los derechos morales de la obra antes citada. En concordancia, suscribo
este documento en el momento que hago entrega del trabajo final en formato impreso y
digital a la Biblioteca de la Universidad Politécnica Salesiana.

Quito, D.M., diciembre de 2019

………………………………….
Cristian Marcelo Díaz García
C.C.: 1723523708

III
1. ÍNDICE GENERAL

1 INTRODUCCIÓN ................................................................................................................ 2
2 MARCO TEÓRICO .............................................................................................................. 3
2.1 Despacho económico de generación .................................................................................. 3
2.2 Generación termoeléctrica ................................................................................................. 3
2.2.1 Función de costos ..................................................................................................... 4
2.2.2 Costo incremental ..................................................................................................... 4
2.3 Despacho económico simplificado .................................................................................... 5
2.4 Despacho económico con pérdidas .................................................................................... 6
3 MÉTODOS DE SOLUCIÓN DEL DESPACHO ECONÓMICO CON PÉRDIDAS .......... 6
3.1 Método Lambda iterativo ................................................................................................... 6
3.1.1 Fórmula de pérdidas .................................................................................................. 7
3.2 Método del Jacobiano ...................................................................................................... 10
4 IMPLEMENTACIÓN DEL MODELO MATEMÁTICO .................................................. 12
5 ANÁLISIS DE RESULTADOS ......................................................................................... 14
6 CONCLUSIONES .............................................................................................................. 16
6.1 Trabajos futuros ............................................................................................................... 16
7 REFERENCIAS .................................................................................................................. 16
7.1 Matriz de Estado del Arte ................................................................................................ 19
7.2 Resumen de Indicadores .................................................................................................. 20

IV
2. ÍNDICE DE FIGURAS

Figura 1. Flujo de generación-costo ............................................................................................... 3

Figura 2. Esquema del Despacho Económico ................................................................................ 3
Figura 3. Esquema de una central termoeléctrica a gas ................................................................. 4
Figura 4. Curva entrada-salida para una unidad generadora térmica ............................................. 4
Figura 5. Ngen unidades termogeneradoras comprometidas a cubrir la carga Pload ......................... 6
Figura 6. Ngen unidades térmicas cubriendo una carga Pload a través de una red de transmisión. ... 6
Figura 7. Sistema de cuatro barras ................................................................................................. 7
Figura 8. Algoritmo del método lambda iterativo. ....................................................................... 13
Figura 9. Algoritmo del método Jacobiano. ................................................................................. 13
Figura 10. Sistema de prueba para DEP ...................................................................................... 14
Figura 11. Costo incremental resultante en cada método ............................................................ 15
Figura 12. Potencia despachada en cada método ......................................................................... 15
Figura 13. Sistema de prueba para DEP ...................................................................................... 15
Figura 14. Costo incremental resultante en cada método ............................................................ 16
Figura 15. Potencia despachada en cada método ......................................................................... 16
Figura 16. Resumen e indicador de la temática - Estado del arte ................................................ 20
Figura 17. Indicador de formulación del problema - Estado del arte........................................... 21
Figura 18. Indicador de solución - Estado del arte....................................................................... 21

3. ÍNDICE DE TABLAS

Tabla 1. Datos de líneas en [p.u.]................................................................................................. 14

Tabla 2. Datos de barras en [p.u.] ................................................................................................ 14
Tabla 3. Generación y voltajes resultantes del flujo de carga en [p.u.] ....................................... 14
Tabla 4. Comparación de resultados ............................................................................................ 14
Tabla 5. Datos de líneas en [p.u.], 5 barras. ................................................................................. 15
Tabla 6. Datos de barras en [p.u.] ................................................................................................ 15
Tabla 7. Comparación de resultados ............................................................................................ 15
Tabla 8. Matriz de estado del arte ................................................................................................ 19

V
DESPACHO ECONÓMICO DE GENERACIÓN
TERMOELÉCTRICA USANDO MÉTODOS DE
SENSIBILIDAD AC
Resumen Abstract
La generación eléctrica está ligada a Electricity generation is linked to
factores técnicos propios de cada proceso technical factors characteristic of each
de transformación de la energía y también energy transformation process and also
a factores económicos que aquello genera, the economic factors that it generates, the
mientras menos costoso sea el producir less expensive it is to produce electricity,
electricidad, esta se puede generar a mayor this can be generated to a greater amount
cantidad satisfaciendo cada vez más a una satisfying more and more a greater
mayor cantidad de usuarios, es decir number of users, that is, covering the
cubriendo la demanda de potencia del power demand of the system with reduced
sistema con costos reducidos. Esta se costs. This becomes the great task to be
convierte en la gran tarea por realizar para done for operators of the electrical
operadores de los SEP, la cual es conocida systems, which is known as the Economic
como el Despacho Económico de Power Dispatch. This paper presents two
Generación. En este trabajo se presentan methods of resolution of the economic
dos métodos de resolución del despacho dispatch, each of these has its own
económico, cada uno de estos tiene su formulation to consider the power losses
propia formulación para considerar las that exist in the electrical transmission in
pérdidas de potencia que existen en la the electric networks. The implemented
transmisión eléctrica en las redes methods are iterative Lambda and the
eléctricas. Los métodos implementados Jacobian, which allow determine the
son Lambda iterativo y Jacobiano, los penalty applied to the generators after the
cuales permiten determinar la losses in the power system have increased;
penalización que se aplica a los in both cases it is necessary first to know
generadores tras incrementarse las the solution of the load flow of the system.
pérdidas en el sistema de potencia; en The results obtained for a test system
ambos casos es necesario primero conocer indicate that the iterative Lambda method
la solución del flujo de carga del sistema. solves a lower operating cost at a lower
Los resultados obtenidos para un sistema number of iterations, which represents
de prueba indican que el método Lambda economic and computational savings of
iterativo entrega un menor costo de calculation.
operación a un menor número de
iteraciones, lo cual representa ahorro
económico y computacional de cálculo.

Palabras Clave: Despacho Económico, Keywords: Economic Dispatch, Iterative

Lambda iterativo, Multiplicadores de Lambda, Lagrange Multipliers,
Lagrange, Lagrangiano, Flujo AC Lagrangian, AC Flows, Newton-Raphson,
Newton-Raphson, Costo Incremental, Incremental Cost, Thermoelectric
Generación Termoeléctrica, Jacobiano. Generation, Jacobian.

1
1 Introducción de un ciclo termodinámico, para
finalmente, transformar la energía
La energía eléctrica se ha convertido en mecánica en electricidad con el uso del
un factor esencial para el desarrollo de generador, principalmente sincrónico. El
cualquier país, siendo de uso ciclo termodinámico es básicamente un
indispensable en las aplicaciones proceso de baja eficiencia. Las
comerciales, industriales y residenciales. eficiencias más altas para las plantas
Este tipo de energía se genera en grandes generadoras modernas de gran tamaño
magnitudes en puntos específicos y luego está por encima del 40%, mientras que
se la transmite económicamente por las plantas más pequeñas pueden tener
largas distancias a los centros de eficiencias menores [1][3].
consumo, donde es transformada, En una central termoeléctrica, la
aprovechada y para aplicaciones [1]. turbina acoplada a un generador
Debido a su importancia, se dedica eléctrico, entrega energía mecánica para
muchos esfuerzos en el estudio de los que a su vez el generador eléctrico
métodos de generación eléctrica tratando transforme ese movimiento en energía
de que estos procesos sean cada vez más eléctrica con condiciones adecuadas para
eficientes y confiables; y en los últimos la integración de este generador al SEP.
años buscando también los métodos para El movimiento de la turbina es generado
producir electricidad de manera que se por inyección de vapor o gas, en cada
genere el menor impacto ambiental. caso con características constructivas
Aunque aún no se puede optimizar del distintas.
todo este último punto, pues sigue siendo En el caso de las centrales térmicas
indispensable el uso de fuentes primarias convencionales, el precio por la
de energía a partir de combustibles obtención de combustible es alto [4][5],
fósiles, generando por lo tanto emisiones además de la operación y mantenimiento
contaminantes, como por ejemplo las de dichas centrales. Es decir, el costo de
centrales térmicas. Sin embargo, el uso producción de electricidad en este tipo de
de este tipo de centrales es altamente centrales es considerable, por tanto, se
necesario para mejorar estabilidad y convierte en un factor altamente
robustez de los sistemas eléctricos de importante al momento de ser incluidas
potencia (SEP) debido a que las centrales como aporte para cubrir la demanda
termoeléctricas presentan mejor energética de un sistema.
respuesta inercial frente a otro tipo de El despacho económico de energía
generadoras, lo cual es de gran ayuda eléctrica se lo realiza en función del costo
para compensar las desviaciones de de cada central de generación, con el
frecuencia que ocurren en el sistema. [2]. propósito de abastecer la carga dentro de
Convencionalmente, la energía límites operativos normales en los
eléctrica en centrales térmicas se obtiene componentes del SEP.
mediante la conversión de combustibles En este documento se analiza los
fósiles (gas natural, carbón o petróleo) y métodos de despacho económico
fuentes de hidrocloruro de azufre. El considerando las pérdidas de potencia en
calor obtenido por la quema de el sistema. Los dos métodos que se
combustibles fósiles o por la fisión de analizan en el presente documento hacen
materiales nucleares da paso a la uso de los flujos de carga AC, los cuales
transformación de la energía calórica en parten de los datos conocidos de un flujo
energía mecánica en una turbina a través de carga para establecer la potencia

2
 óptima en cada generador en la cual el En la Figura 2 se muestra el esquema
costo de producción sea el mínimo. Esta de trabajo del despacho económico (DE),
idea se resume gráficamente en la Figura en donde se presenta una demanda
1; en la cual, la entrada del proceso es la determinada de potencia, el proceso de
inyección de combustible con su DE analiza los costos de cada generador,
respectivo costo, ingresan al proceso de y como respuesta ordena la cantidad de
generación eléctrica en cualquier tipo de potencia que cada uno de los generadores
planta generadora termoeléctrica, dentro ha de producir para que se cumpla la
de la planta de generación también condición de despacho al menor costo
existen costos de operación, para el satisfaciendo la demanda.
análisis en este trabajo se engloba a todos
los costosCAMARA
COMBUSTIBLE en una
DE sola función que UNIDADES
GENERADORAS
COMBUSTION
representa todos los costos de GEN 1 GESTION DE
DATO DEL
VALOR DE
producción incluyendo los valores de $ combustible DESPACHO
POTENCIA
ELECTRICIDAD
operación y combustible. Entonces a la GEN 2 DESPACHO
DEMANDADA
AIRE COMPRESOR TURBINA
salida del proceso se tiene electricidad $ combustible ECONÓMICO
GENERADOR
producida a un determinado costo de GEN N
$ combustible
producción. GEN1
[MW]

ENTRADA PROCESO SALIDA DEMANDA DE

GEN2
POTENCA
GENERACION
POTENCIA ASIGNACION DE [MW]
COMBUSTIBLE [MW]
ELECTRICA CARGA A CADA
TERMOELECTRICA
costo combustible costo de UNIDAD GEN
costo operación
producción N ENTREGA DE
[MW] POTENCIA
Figura 1. Flujo de generación-costo DEMANDADA AL
SISTEMA
2 Marco teórico Figura 2. Esquema del Despacho Económico

2.1 Despacho económico de 2.2 Generación termoeléctrica

generación Las centrales de generación
El despacho económico de generación es termoeléctrica usan la transformación de
una tarea importante en la gestión de energía calorífica para la generación de
sistemas modernos de energía, el cual electricidad, este tipo de energía se
tiene como objetivo, asignar para obtiene de variados combustibles como,
cualquier condición de carga en el por ejemplo; el gas natural, carbón,
sistema, la salida de potencia a cada diésel. Para cada tipo de combustible
central de generación (y a cada unidad de utilizado, el equipamiento y
generación dentro de una misma planta), configuración de la central es diferente.
haciendo énfasis en la coordinación y Las centrales termoeléctricas se pueden
disminución de los costos de generación configurar en diferentes ciclos
eléctrica en todas las plantas generadoras termodinámicos tales como el ciclo
operativas dentro del sistema [6][7][8]. Brayton, el cual corresponde al uso de la
Así, si se conoce el flujo de cargas en un turbina de gas; el ciclo Rankine conocido
SEP, es posible hallar la potencia activa como ciclo convencional de vapor; el
y reactiva que tendría que aportar cada ciclo Diésel, el cual usa combustibles
planta generadora para cumplir con la líquidos en motores de combustión
demanda de potencia [9][5][10], en este interna; y el ciclo Otto, el cual también
caso sin tomar en cuenta el costo de usa combustibles líquidos en motores
operación. alternativos de chispa o bujía. Un
ejemplo de central termoeléctrica es la

3
que se muestra en el esquema de la [MBTU/h]. El ingreso del combustible
Figura 3, la cual corresponde a una está en función de la potencia de salida
central a gas de ciclo Brayton, en el cual (Pgi) dado en [𝑀𝑊]. Si se incluye en la
se comprime aire y en una cámara de curva de entrada - salida, el costo en
combustión se mezcla la inyección de dólares del combustible dado en [$/
combustible y el aire comprimido, esta 𝑀𝐵𝑇𝑈] u otras unidades equivalentes, la
mezcla al combustionarse aumenta su expresión queda de la siguiente manera
temperatura, la cual entra a la turbina de (𝐻(𝑃𝑔𝑖 ) [𝑀𝐵𝑇𝑈/ℎ])𝑥 [$/𝑀𝐵𝑇𝑈] y como
gas y se expande moviendo la turbina que resultado tenemos que la curva inicial se
está acoplada al generador eléctrico. En ha transformado en la curva de costos o
resumen, la Figura 3 muestra al conocida como la función de costos
combustible como componente de 𝐹(𝑃𝑔𝑖 ) [$/ℎ] que expresa el costo en [$/h]
entrada y como componente de salida la de producir una cantidad de [MW]. [4].
potencia eléctrica generada.

COMBUSTIBLE CAMARA DE
COMBUSTION

ELECTRICIDAD
AIRE COMPRESOR TURBINA

GENERADOR

ENTRADA PROCESO SALIDA

POTENCIA
GENERACION
COMBUSTIBLE ELECTRICA
TERMOELECTRICA
costo combustible costo de
costo operación
producción

Figura 4. Curva entrada-salida para una unidad

Figura 3. Esquema de una central termoeléctrica a gas
generadora térmica
2.2.1 Función de costos
El modelo económico de la central La función de costo puede ser una
termoeléctrica se expresa por el costo de expresión lineal, polinómica de segundo
la cantidad de combustible a emplearse o tercer orden, o una función constante a
para abastecer la potencia demandada, la trozos [4].
potencia eléctrica que demanda el
sistema es la entrada de la función y 𝐹(𝑃𝑔𝑖 ) = 𝑎𝑃𝑔𝑖 2 + 𝑏𝑃𝑔𝑖 + 𝑐 [$/ℎ] (1)
aquel costo corresponde a la salida de la
función, además se toma en cuanta 2.2.2 Costo incremental
restricciones como por ejemplo los El costo adicional en [$/h] necesarios
límites máximos y mínimos de para generar 1 [MW] adicional de
producción de los generadores y demás potencia, se lo conoce como costo
elementos de una central térmica [4]. incremental, denotado por λ.
En la Figura 4, se visualiza la curva El costo incremental se determina al
entrada-salida de una planta medir la pendiente de la curva de costos,
termoeléctrica representada por 𝐻(𝑃𝑔𝑖 ), es decir, la curva de costo incremental es
con unidades en [𝑀𝐵𝑇𝑈/ℎ]. Esta gráfica la derivada de la curva de costo con
muestra el ingreso de combustible a una respecto a la potencia de generación, de
generadora termoeléctrica, dado en tal manera que, si la curva de costos es

4
una función cuadrática, la curva de costo operar con un mismo valor de costo
incremental será una función lineal incremental, ver ecuación 4. A esto se lo
[4][8]. conoce como Despacho Económico [8],
el cual puede resolverse mediante un
Función de costo cuadrática: método de optimización.

𝐹(𝑃𝑔𝑖 ) = 𝑎𝑃𝑔𝑖 2 + 𝑏𝑃𝑔𝑖 + 𝑐 [$/ℎ] (2) 𝑑𝐹1

𝜆= = 𝑎1 𝑃𝑔1 + 𝑏1
𝑑𝑃𝑔1
Función de costo incremental: (4)
𝑑𝐹𝑖 𝑑𝐹2
𝜆𝑖 = = 2𝑎𝑃𝑔𝑖 + 𝑏 [$/𝑀𝑊ℎ] 𝜆= = 𝑎2 𝑃𝑔2 + 𝑏2
(3) 𝑑𝑃𝑔2
𝑑𝑃𝑔𝑖
En la Figura 5 se muestra un sistema
2.3 Despacho económico con una cantidad Ngen de generadores
simplificado termoeléctricos, como entrada en cada
Para comprender este concepto, se puede unidad se tiene la función de costo
suponer que existen dos unidades de correspondiente Fi, y a su salida se tiene
generación en una o más plantas de un la generación de potencia eléctrica Pgi de
sistema. La salida de potencia de la cada unidad conectadas directamente a
planta es provista por la suma de las una única barra en la cual también se
potencias que generan cada una de sus encuentra recibiendo potencia una carga
dos unidades, y cada unidad produce Pload. La suma de todas esas potencias
energía eléctrica a un costo incremental generadas alimenta la carga Pload, y el
específico para cada una, siendo uno de costo de generación (u operación) total es
estos costos más alto que el otro. la suma de los costos individuales.
Entonces, se puede suponer que la El despacho de generación en este
unidad que tiene el costo incremental caso se lo conoce como Despacho
mayor transfiere parte de su carga de Económico Simplificado (DES), en el
potencia a la unidad de menor costo, esa cual, se desprecia las pérdidas de
reducción de carga en la unidad más potencia en el SEP.
costosa da como resultado una La función objetivo (FO) en el DES es
disminución en el costo de producción, y minimizar el costo total (FT) para
este costo disminuido es mayor al costo alimentar a la carga teniendo como
que se incrementa en la segunda unidad restricción que la suma de las potencias
por la toma de esa carga. Ya que los generadas debe ser igual a la carga del
costos incrementales de las unidades son sistema. Esta restricción es conocida
funciones de las potencias generadas, el como Balance de Potencia [11], además
proceso de transferencia de carga de una de los límites máximos y mínimos de
unidad a otra puede continuar hasta generación dados en las curvas de
alcanzar el punto en el cual los costos capacidad de los generadores
incrementales de las unidades sean sincrónicos.
iguales. Esto parte del hecho de que los
Si se aplica el mismo razonamiento generadores termoeléctricos necesitan de
para sistemas con más de dos unidades, un determinado lapso de tiempo para su
se plantea el criterio en que, para tener un puesta en marcha, y también para su
despacho óptimo de carga en la parada y reanudación, esto hace difícil
generación eléctrica, las unidades deben

5
que se los pueda sacar de generación la generación debe cubrir la potencia de
fácilmente [2]. carga y las pérdidas de potencia.
La función objetivo para este
P1
despacho se mantiene igual que en el
P1 pero son las restricciones
anterior caso,
las que se modifican. En la restricción de
F1
F1 balance de potencia ahora incluye a las
P2
pérdidas deP2transmisión [11].
RED DE
F2 TRANSMISIÓN
P1 F2
P1PÉRDIDAS
CON
Pload Pload
F1 PN
F1 PN
P2
P2
FN
FN
RED DE
Figura 5. NF2
gen unidades termogeneradoras TRANSMISIÓN
F2
comprometidas a cubrir la carga Pload CON PÉRDIDAS
Pload Pload
La formulación matemática
PN de la
PN
solución del DES es la siguiente:
FN
𝐹𝑇 = 𝐹1 + 𝐹2 + 𝐹3 + ⋯ + 𝐹𝑁 (5) FN
Figura 6. Ngen unidades térmicas cubriendo una carga
𝑁 Pload a través de una red de transmisión.
𝐹𝑂: 𝑚𝑖𝑛 → 𝐹𝑇 = ∑ 𝐹𝑖 (𝑃𝑔𝑖 ) (6)
𝑖=1 Por tanto, se muestra las expresiones
S.A.: matemáticas de las restricciones
𝑁 complementadas para el DEP [4].
𝑃𝑙𝑜𝑎𝑑 − ∑ 𝑃𝑔𝑖 = 0 (7)
𝑖=1 S.A.:
𝑚𝑖𝑛 𝑚𝑎𝑥 𝑁
𝑃𝑔𝑖 ≤ 𝑃𝑔𝑖 ≤ 𝑃𝑔𝑖 (8)
𝑃𝑙𝑜𝑎𝑑 + 𝑃𝑙𝑜𝑠𝑠 − ∑ 𝑃𝑔𝑖 = 0 (9)
𝑖=1
2.4 Despacho económico con
𝑚𝑖𝑛 𝑚𝑎𝑥
pérdidas 𝑃𝑔𝑖 ≤ 𝑃𝑔𝑖 ≤ 𝑃𝑔𝑖 (10)
Una vez comprendido el DES como un
método rápido y simple de estimar la
distribución de generación, en la práctica 3 Métodos de solución del
es necesario considerar a los SEP con sus despacho económico con
pérdidas de potencia en el sistema
transmisión Ploss, como se representa en
pérdidas
la Figura 6, entonces el problema pasa a 3.1 Método Lambda iterativo
ser nombrado Despacho Económico con La función de Lagrange o conocida como
Perdidas (DEP). Lagrangiano de la ecuación (11), se
En el caso del DES, el criterio de construye al relacionar los costos totales
balance de potencia debe asegurar que la FT añadido al balance de potencia Ø con
generación cubra la carga, pero para el el costo incremental λ que es
caso de un sistema con pérdidas, el denominado también como multiplicador
criterio de balance de potencia dicta que

6
de Lagrange, y corresponde al costo valor del costo incremental (λ) debe ser
incremental efectivo, para los sistemas igual para todas las unidades de
en los que se considera las pérdidas de generación. Con el uso de las ecuaciones
potencia en la transmisión [11]. (13) y (17), se construye el sistema de
ecuaciones en las cuales las variables
ℒ = 𝐹𝑇 + 𝜆Ø (11) serán las potencias Pgi y un único lambda
o conocido como lambda del sistema
Si se desea encontrar el valor mínimo [8][11].
de la función de Lagrange, cual es la
función objetivo del DEP, se toma la
3.1.1 Fórmula de pérdidas
derivada de la función de Lagrange con
Las pérdidas en la transmisión son una
respecto a cada uno de potencias
función de la potencia generada y tienen
generadas individuales y se iguala a cero
relación directa con las impedancias del
la ecuación, de modo que se obtiene una
sistema, para lo cual es necesario el
ecuación en la que están involucradas
cálculo de las corrientes de los
tanto la derivada de las pérdidas de
generadores como se muestra a
potencia como la derivada de la función
continuación.
de costos, las dos con respecto a la
Suponiendo un sistema de cuatro
potencia de cada unidad de generación y
barras, con la barra 1 y 2 con
también el término del costo incremental
generadores, y en la barra 3 y 4 con
[11].
cargas, y un nodo n definido como neutro
𝑑ℒ 𝑑 𝐹𝑖 𝑑𝑃𝑙𝑜𝑠𝑠 del sistema.
= − 𝜆 (1 − )
𝑑𝑃𝑔𝑖 𝑑𝑃𝑔𝑖 𝑑 𝑃𝑔𝑖 (12) GENERADOR 1
I3 3
=0 1 I1
CARGA
I4 ZBARRA I2
Si se ordena esta ecuación, da como CARGA
resultado la ecuación del costo 4 2
In GENERADOR 2
incremental para cada unidad de
n
generación, con respecto a su propia
Figura 7. Sistema de cuatro barras
función de costos, la potencia generada y
las pérdidas totales del sistema [11]. Entonces la carga del sistema ID se
compone de:
𝑑 𝐹𝑖 1 𝐼3 + 𝐼4 = 𝐼𝐷
( ) =𝜆 (13) (14)
𝑑 𝑃𝑔𝑖 1−𝑑𝑃𝑙𝑜𝑠𝑠
𝑑𝑃𝑔𝑖
Por tanto, cada carga es una porción
Donde: de la carga total:

1 𝐼3 = 𝑑3 𝐼𝐷 ; 𝐼4 = 𝑑4 𝐼𝐷 (15)
𝑑𝑃 : Factor de penalización
1− 𝑙𝑜𝑠𝑠
𝑑𝑃𝑔𝑖

𝑑𝑃𝑙𝑜𝑠𝑠 𝑑3 + 𝑑4 = 1 (16)
𝑑𝑃𝑔𝑖
: Pérdidas incrementales
Como cada una de las corrientes está
En el factor de penalización se basa el relacionada a las impedancias de las
DEP Como se había mencionado líneas en que circula, se plantea un
anteriormente, en el criterio DEP, el sistema de ecuaciones, y los voltajes se

7
miden con respecto a la barra de
referencia n. 𝐼1 1 . .
𝐼1
𝐼2 . 𝑍22 1
[ ]=[ ] [𝐼2 ]
𝑉1𝑛 𝑍11 𝑍12 𝑍13 𝑍14 𝐼1 𝐼3 −𝑑3 𝑡1 −𝑑3 𝑡2 −𝑑3 𝑡1 (25)
𝐼0
𝑉2𝑛 𝑍21 𝑍22 𝑍23 𝑍24 𝐼2 𝐼4 −𝑑4 𝑡1 −𝑑4 𝑡2 −𝑑4 𝑡1
[ ]=[
𝑉3𝑛 𝑍31 𝑍32 𝑍33
][ ]
𝑍34 𝐼3
(17)
𝑉4𝑛 𝑍41 𝑍42 𝑍43 𝑍44 𝐼4 𝐼1
𝐼1
𝐼2
[ ] = 𝐶 [𝐼2 ] (26)
Se toma la primera fila de la matriz 𝐼3
𝐼0
para reemplazar los términos de la 𝐼4
ecuación (15), y expresar la ecuación Entonces la ecuación de las pérdidas
(19) para ID. de potencia toma la siguiente expresión:
𝑉1𝑛 = 𝑍11 𝐼1 + 𝑍12 𝐼2 + 𝑍13 𝐼3 + 𝑍14 𝐼4 (18) 𝑃𝑙𝑜𝑠𝑠
𝐼1 ∗
= [𝐼1 𝐼2 𝐼𝑛0 ][𝐶 𝑇 𝑅𝑏𝑎𝑟𝑟𝑎 𝐶 ∗ ] [𝐼2 ] (27)
−𝑍11 −𝑍12
𝐼𝐷 = 𝐼 + 𝐼 + 𝐼0
𝑑3 𝑍13+ 𝑑4 𝑍14 1 𝑑3 𝑍13+𝑑4 𝑍14 2
−𝑍11 (19)
𝐼𝑛0
𝑑 𝑍 𝑑 𝑍
3 13+ 4 14 Donde:
Rbarra: es la componente real de Zbarra
El término In0 corresponde a la
corriente de carga nula: La ecuación (27) define las pérdidas
de potencia del sistema, y se puede notar
𝑉1𝑛 que la ecuación es una función de las
𝐼𝑛0 = − (20)
𝑍11 corrientes I1, I2, ya que la transformación
C es una expresión invariante y el
Se define los términos t1 y t2 para término In0 es un número complejo de
simplificar la ecuación (19) valor constante.
Como siguiente paso, para este
𝑍11 𝑍12
𝑡1 =
𝑑3 𝑍13+𝑑4 𝑍14
; 𝑡2 =
𝑑3 𝑍13+𝑑4 𝑍14
(21) análisis de las pérdidas de potencias
reales del sistema, se considera a la
potencia reactiva Qgi como una porción
𝐼𝐷 = −𝑡1 𝐼1 −𝑡2 𝐼2 −𝑡1 𝐼𝑛0 (22) s1 constante dentro de la potencia Pgi. Los
valores s1 y s2 son números reales.
Al reemplazar ID de la ecuación (22)
en (15)
𝑠1 = Q𝑔1 /𝑃𝑔1 ; 𝑠2 = Q𝑔2 /𝑃𝑔2 (28)
𝐼3 = −𝑑3 𝑡1 𝐼1 −𝑑3 𝑡2 𝐼2 −𝑑3 𝑡1 𝐼𝑛0 (23)

Por tanto, es lo mismo suponer que los

𝐼4 = −𝑑4 𝑡1 𝐼1 −𝑑4 𝑡2 𝐼2 −𝑑4 𝑡1 𝐼𝑛0 (24) generadores operan con fator de potencia
constante, en el mismo período de
Las ecuaciones (23) y (24), se definen análisis, de tal modo se puede expresar
como la transformación C de las las potencias de la siguiente forma:
corrientes iniciales I1, I2, I3, I4, expresadas
ahora como un conjunto de corrientes
𝑃𝑔1 + 𝑗Q𝑔1 = (1 + 𝑗s1 )𝑃𝑔1
nuevas I1, I2, In0. (29)
𝑃𝑔2 + 𝑗Q𝑔2 = (1 + 𝑗s2 )𝑃𝑔2

8
 Al usar estas nuevas expresiones en términos imaginarios fuera de la diagonal
los términos de las corrientes de los y se obtiene el duplo de los términos
generadores, se tiene las siguientes reales simétricos de Tα, que se expresan
ecuaciones: de esta manera.
𝐵10⁄
(1 − 𝑗s1 ) 𝐵11 𝐵12 2
𝐼1 = P𝑔1 = 𝛼1 P𝑔1 ∗
𝑉1∗ 𝐵20⁄ = 𝑇𝛼 + 𝑇𝛼 (33)
(30) 𝐵21 𝐵22 2
(1 − 𝑗s2 ) 2
𝐼2 = P𝑔2 = 𝛼2 P𝑔2 𝐵10 𝐵20⁄
𝑉2∗ [ ⁄2 2 𝐵00 ]

Expresadas en forma matricial: Si se añade este complejo conjugado a

la ecuación (32), se obtiene:
𝐼1 𝛼1 . . 𝑃𝑔1 𝑃𝑙𝑜𝑠𝑠
𝐼 𝛼2 𝐵10⁄
[ 2] = [ . .
] [ 𝑃2 ] (31) 𝐵11 𝐵12 2 𝑃𝑔1
𝐼𝑛0 . . 𝐼𝑛0 1 = [𝑃𝑔1 𝑃𝑔2 1 ] 𝐵21 𝐵20⁄ [𝑃 ]
𝐵22 2 𝑔2
𝐵10 𝐵20⁄ 1
[ ⁄2 2 𝐵00 ]
Se sustituye las corrientes en la
función de las pérdidas de potencia de la (34)
ecuación (27).
Donde B12 será igual a B21. Al realizar
la multiplicación filas por columnas
𝑃𝑙𝑜𝑠𝑠
𝑃𝑔1 𝑇 𝛼1 . . 𝛼1 . . ∗ 𝑃𝑔1 ∗ resulta.
= [ 𝑃2 ] [ . 𝛼2 . ] 𝐶 𝑇 𝑅𝑏𝑎𝑟𝑟𝑎 𝐶 ∗ [ . 𝛼2 . ] [𝑃𝑔2 ]
1 ⏟. . 𝐼𝑛0 . . 𝐼𝑛0 𝑃𝑙𝑜𝑠𝑠 = 𝐵11 𝑃𝑔1 2 + 2𝐵12 𝑃𝑔1 𝑃𝑔2 + 𝐵22 𝑃𝑔2 2
1
+ 𝐵10 𝑃𝑔1 + 𝐵20 𝑃𝑔2 + 𝐵00
Tα
2 2 2
(32) = ∑ ∑ 𝑃𝑔𝑖 𝐵𝑖𝑗 𝐵𝑔𝑗 + ∑ 𝐵𝑖0 𝑃𝑔𝑖 + 𝐵00
𝑖=1 𝑗=1 𝑖=1
Con el uso de la matriz de tipo (35)
Hermitiana, la cual es una propiedad de
Lo que equivale a decir:
las matrices que dice; la transpuesta de la
multiplicación de matrices también se 𝐵11 𝐵12 𝑃𝑔1
𝑃𝑙𝑜𝑠𝑠 = [𝑃𝑔1 𝑃𝑔2 ] [ ][ ]
puede expresar como la multiplicación 𝐵21 𝐵22 𝑃𝑔2
en orden contrario de sus transpuestas 𝐵
+ [𝑃𝑔1 𝑃𝑔2 ] [ 10 ] + 𝐵00
individuales (ABC)T=ATBTCT. Si se 𝐵20
incluye el complejo conjugado se obtiene
(ABC)T*=AT*BT*CT*; entonces se (36)
determina que el término Tα tomado de la Se puede expresar en forma general:
ecuación (32), que puede expresarse
𝑃𝑙𝑜𝑠𝑠 = 𝑃𝑔 𝑇 𝐵𝑃𝑔 + 𝑃𝑔 𝑇 𝐵0 + 𝐵00 (37)
como el complejo conjugado de la
transpuesta del mismo término. De esta Entonces para sistemas en los cuales
matriz Hermitiana se obtiene que los existen más de dos generadores, un
elementos de la diagonal son números número K de generadores, en la ecuación
reales y los elementos fuera de la de pérdidas el valor K define la
diagonal denominados mij son iguales al dimensión de las matrices y los vectores
complejo conjugado de su de tal ecuación, y esta queda expresada
correspondiente mji. De la suma de Tα de la siguiente forma:
con Tα* resulta la eliminación de los

9
 𝐾 𝐾 𝐾 Lo que en forma compacta es:
𝑃𝑙𝑜𝑠𝑠 = ∑ ∑ 𝑃𝑔𝑖 𝐵𝑖𝑗 𝑃𝑔𝑗 + ∑ 𝐵𝑖0 𝑃𝑔𝑖 + 𝐵00 ̅̅̅̅̅̅̅
∆𝑃, 𝑄 = |𝐽| + ̅̅̅̅̅̅̅
∆𝛳, 𝑉
𝑖=1 𝑗=1 𝑖=1 |J| = Jacobiano
(40)

(38) En donde las submatrices del

Jacobiano tienen la siguiente definición
El término B es una matriz simétrica [14].
conocida como la matriz de coeficientes
𝑑𝑃𝑖
de pérdidas. Los coeficientes llevan 𝐻𝑖𝑗 =
𝑑𝛳𝑗
como unidad [1/MW] cuando las 𝑑𝑃𝑖
unidades de las potencias Pg1 a PgK están 𝑁𝑖𝑗 =
𝑑𝑉𝑗
dadas en [MW], y solo el término B00 es 𝑑𝑄𝑖 (41)
adimensional [8][12]. 𝑀𝑖𝑗 =
𝑑𝛳𝑗
𝑑𝑃𝑖
3.2 Método del Jacobiano 𝐿𝑖𝑗 =
𝑑𝑉𝑗
La matriz Jacobiano es una matriz real
casi simétrica de los flujos de potencia Los cuales pueden ser determinados
AC [13]. Este método de desarrollo del con las potencias calculadas y términos
DEP, utiliza los flujos de potencia AC de conductancia y susceptancia de la
iterativamente para actualizar las matriz Ybarra [16]. Los términos de las
potencias de despacho en función de las submatrices H, N, M y L se indican a
correcciones de voltajes y ángulos de las continuación.
barras. En este caso se usa el método de
Newton-Raphson (NR), el cual realiza un
𝐻𝑖𝑖 = −𝑄𝑖 − 𝐵𝑖𝑖 |𝑉𝑖 |2 (42)
cálculo de las correcciones mientras
toma en cuenta todas las iteraciones k 𝑁𝑖𝑖 = 𝑃𝑖 + 𝐺𝑖𝑖 |𝑉𝑖 |2 (43)
[12][14]. A continuación, se resume cada
uno de los pasos usados para el desarrollo 𝑀𝑖 = 𝑃𝑖 − 𝐺𝑖𝑖 |𝑉𝑖 |2 (44)
del DEP con el uso del método del
Jacobiano. 𝐿𝑖𝑖 = 𝑄𝑖 − 𝐵𝑖𝑖 |𝑉𝑖 |2 (45)
Como punto de partida se realiza un 𝐻𝑖𝑗 = 𝐿𝑖𝑗 = |𝑉𝑖 ||𝑉𝑗 |(𝐺𝑖𝑗 sin 𝛳𝑖𝑗
flujo de potencia NR, dado cualquier − cos 𝛳𝑖𝑗 )
(46)
sistema con N número de barras, de las 𝑁𝑖𝑗 = −𝑀𝑖𝑗 = |𝑉𝑖 ||𝑉𝑗 |(𝐺𝑖𝑗 cos 𝛳𝑖𝑗
cuales; M número son barras de carga y + sin 𝛳𝑖𝑗 )
(47)
S son barras de voltaje controlado y una
barra oscilante. Se plantea 2M+S número Donde:
de ecuaciones para resolver 2M+S Vi, Vj: Magnitud de voltajes en las barras
número de incógnitas [15]. Se plantea la 𝛳𝑖 , 𝛳𝑗 : Ángulos de voltajes en las barras
ecuación matricial: 𝐺𝑖𝑗 : Conductancia de la matriz Ybarra en
∆𝑃𝑔1
ubicación ij
∆𝛳1
∆𝑃𝑔2 ∆𝛳2
𝐵𝑖𝑗 : Suceptancia de la matriz Ybarra en
. . ubicación ij
∆𝑃𝑔𝑁 = [ 𝐻 𝑁
] ∆𝛳𝑁 𝛳𝑖𝑗 = 𝛳𝑖 − 𝛳𝑗
𝑀 𝐿
∆𝑉1
(39)
∆𝑄𝑔1
. . Entonces el algoritmo de solución del
[∆𝑄𝑔𝑁 ] [ ∆𝑉𝑁 ]
método NR consiste en encontrar la
𝐽𝑎𝑐𝑜𝑏𝑖𝑎𝑛𝑜
actualización de magnitudes de voltajes y

10
ángulos como se indican en las Se hace uso del siguiente criterio: en
ecuaciones (48) y (49), hasta cumplir con despacho económico, todos los
un criterio de convergencia. generadores se mantienen trabajando al
mismo costo de producción ante una
𝑉𝑖𝑘 = 𝑉𝑖𝑘−1 + ∆𝑉𝑖 (48) variación de potencia ∆P, es decir, que la
variación de costo de producción es cero
𝛳𝑖𝑘 = 𝛳𝑖𝑘−1 + ∆𝛳𝑖 (49) con un ∆P significativamente pequeño
[11][21], siendo:
Hasta que no se cumpla con la
convergencia, el proceso continúa de 𝐶𝑜𝑠𝑡𝑜 𝑇𝑜𝑡𝑎𝑙 𝑑𝑒 𝑃𝑟𝑜𝑑𝑢𝑐𝑐𝑖𝑜𝑛 = ∑ 𝐹𝑖 (𝑃𝑖 ) (53)
forma iterativa. Una vez que se cumpla la
convergencia se termina el proceso del
Se expresa la variación en el costo de
flujo NR. Entonces, se utiliza los últimos
producción al ocurrir un ∆Pi en la
valores actualizados para hallar un nuevo
generadora i.
Jacobiano, al cual se lo denomina como
Jacobiano del punto de solución
∆𝐶𝑜𝑠𝑡𝑜 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑐𝑖𝑜𝑛
[11][15]. Una vez obtenido este nuevo 𝑑 𝐹𝑖 (𝑃𝑖 ) 𝑑 𝐹𝑟𝑒𝑓 (𝑃𝑟𝑒𝑓 )
= ∆𝑃𝑖 + ∆𝑃𝑟𝑒𝑓 (54)
Jacobiano se puede encontrar los factores 𝑑𝑃𝑖 𝑑𝑃𝑟𝑒𝑓
de nodo directos de las barras. Estos
factores permiten conocer la relación de
Usando los términos de la ecuación
cambio de potencia en la barra de
(51), se tiene:
referencia Pref al existir una variación de
potencia ∆Pi. Si no existiese variación
∆𝐶𝑜𝑠𝑡𝑜 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑐𝑖𝑜𝑛
alguna, ∆Pref se mantiene como el 𝑑 𝐹𝑖 (𝑃𝑖 ) 𝑑 𝐹𝑟𝑒𝑓 (𝑃𝑟𝑒𝑓 ) (55)
negativo de ∆Pi. Sin embargo; los ajustes =
𝑑𝑃𝑖
∆𝑃𝑖 − 𝛽𝑖
𝑑𝑃𝑟𝑒𝑓
∆𝑃𝑖
en la generación causan cambios en los
flujos del sistema, estos cambios en los
En concordancia con el criterio
flujos producen que ∆Pref y ∆Pi no sean
mencionado, debe cumplirse que
necesariamente iguales [19]. Por lo tanto,
∆Costos de Producción=0, entonces debe
se expresa su relación como:
satisfacerse la igualdad:
∆𝑃𝑟𝑒𝑓 = −∆𝑃𝑖 +∆𝑃𝑙𝑜𝑠𝑠 (50)
𝑑 𝐹𝑖 (𝑃𝑖 ) 𝑑 𝐹𝑟𝑒𝑓 (𝑃𝑟𝑒𝑓 )
Entonces se define esta relación de = 𝛽𝑖 (56)
𝑑𝑃𝑖 𝑑𝑃𝑟𝑒𝑓
cambio con el termino β1, al cual se lo
conoce como factor de nodo o también
llamado factor de penalización [11]. Si se hace uso de la ecuación (13) del
costo incremental, y se expresa en
−∆𝑃𝑟𝑒𝑓 (∆𝑃𝑖 − ∆𝑃𝑙𝑜𝑠𝑠 ) términos del factor de nodo, se tiene la
𝛽𝑖 = = (51) siguiente ecuación:
∆𝑃𝑖 ∆𝑃𝑖

En otros términos: 𝑑 𝐹𝑖 1
( )=𝜆 (57)
𝑑 𝑃𝑔𝑖 𝛽𝑖
𝑑𝑃𝑙𝑜𝑠𝑠
𝛽𝑖 = 1 − (52)
𝑑𝑃𝑖 Ahora bien, se tiene el termino Pref, el
cual está en función de los ángulos de
El factor de penalización de la barra fase del sistema y de las magnitudes de
de referencia es βref =1.0 [20][21]. los voltajes.

11
 𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓
Con el Jacobiano calculado en el
∆𝑃𝑟𝑒𝑓 = ∑ ∆𝛳𝑖 + ∑ ∆|𝑉𝑖 | punto de solución y teniendo en cuenta
𝑑𝛳𝑖 𝑑|𝑉𝑖 |
𝑖 𝑖
𝑑𝑃𝑟𝑒𝑓 𝑑𝛳𝑖 𝑑𝑃𝑟𝑒𝑓 𝑑|𝑉𝑖 | (58) que las derivadas de la Pref con respecto a
=∑ ∆𝑃 + ∑ ∆𝑃𝑖 los ángulos y los voltajes corresponden a
𝑑𝛳𝑖 𝑑𝑃𝑖 𝑖 𝑑|𝑉𝑖 | 𝑑𝑃𝑖
𝑖 𝑖
los mismos valores de las submatrices H
y N de la matriz Jacobiano, de la
También expresada en términos de la operación de la ecuación (61) se obtiene
potencia reactiva: las derivadas de la potencia de referencia
con respecto a las potencias de
𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 generación Pref/Pi, que en otros términos
∆𝑃𝑟𝑒𝑓 = ∑ ∆𝛳𝑖 + ∑ ∆|𝑉𝑖 |
𝑑𝛳𝑖 𝑑|𝑉𝑖 | es el negativo del factor de nodo βi como
𝑖 𝑖
𝑑𝑃𝑟𝑒𝑓 𝑑𝛳𝑖 𝑑𝑃𝑟𝑒𝑓 𝑑|𝑉𝑖 | (59) se expresa en la ecuación (51).
=∑ ∆𝑄 + ∑ ∆𝑄𝑖
𝑑𝛳𝑖 𝑑𝑄𝑖 𝑖 𝑑|𝑉𝑖 | 𝑑𝑄𝑖 Al reemplazar los valores de los
𝑖 𝑖
factores de nodo en la ecuación (55) de
costo incremental para construir el
Donde los términos dPref/dϴi y sistema de ecuaciones convencional y
dPref/d|Vi| parten de la diferenciación de ejecutar el cálculo del despacho
las ecuaciones de flujo de potencia económico, se debe considerar que la
estándar para la barra de referencia, y los ecuación de costo incremental de la barra
términos dϴi/ dPi y d|Vi|/dP proviene del de referencia permanece sin alteración.
inverso de la matriz Jacobiano [11][21].
Entonces se puede expresar las
4 Implementación del
ecuaciones (58) y (59) en forma de
ecuaciones matriciales. modelo matemático
𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 Se ha realizado la construcción del
[ … ] modelo matemático en Matlab, en el cual
𝑑𝑃1 𝑑𝑄1 𝑑𝑃2 𝑑𝑄2 𝑑𝑃𝑁 𝑑𝑄𝑁
𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 −1 se aplican los dos métodos de resolución
=[ … ] [𝐽 ]
𝑑𝛳1 𝑑|𝑉1 | 𝑑𝛳2 𝑑|𝑉2 | 𝑑𝛳𝑁 𝑑|𝑉𝑁 | del DEP, tanto con el método de Lambda
Iterativo como también el método del
(60) Jacobiano, los respectivos algoritmos se
Si se transpone las ecuaciones indican en la Figura 8 y Figura 9. En el
matriciales: primer método se obtiene los resultados
del despacho económico con pérdidas; y,
𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 adicional a esto se obtiene los valores del
𝑑𝑃1 𝑑𝛳1 despacho económico sin considerar las
𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 pérdidas del sistema, estos últimos
𝑑𝑄1 𝑑|𝑉1 | resultados se los determina para efectos
𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 de comparación. Los dos métodos
𝑑𝑃2 𝑑𝛳2 conllevan procesos iterativos, los cuales
𝑑𝑃𝑟𝑒𝑓 = [𝐽𝑇−1 ] 𝑑𝑃𝑟𝑒𝑓 (61) buscan dar una respuesta cuasi óptima a
𝑑𝑄2 𝑑|𝑉2 |
. . un SEP. La evaluación de cada proceso
. . se la hace usando los dos métodos para
𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓 definir el despacho económico de un
𝑑𝑃𝑁 𝑑𝛳𝑁 mismo sistema, y comparar resultados
𝑑𝑃𝑟𝑒𝑓 𝑑𝑃𝑟𝑒𝑓
del despacho entre otros aspectos, como
[ 𝑑𝑄𝑁 ] [𝑑|𝑉𝑁 |] por ejemplo el número de iteraciones que
ocupa cada uno para su trabajo.

12
 INICIO se presenta dos SEP diferentes a los
cuales se aplica los algoritmos
CALCULAR UN
DES planteados anteriormente.
INICIALIZA ?,Pgi
INICIO
CALCULAR
PL/Pgi
FLUJO DE CARGA
para cada i=1...N
CALCULA V, ?

GENERA
SISTEMA DE ALMACENA LOS VALORES
ECUACIONES, DE V, ?, CON RESULTADOS
INCLUYE EN EL PUNTO DE SOLUCION
PERDIDAS

RESUELVE ACTUALIZA EL VALOR

SISTEMA DE DE LA POTENCIA Pref
ECUACIONES
ACTUALIZA ?,Pgi
CALCULAR NUEVO
JACOBIANO
CALCULAR
€=Pgi(iter) – Pgi(iter-1)
CALCULAR COEFICIENTES
DE PERDIDAS dPgs/dPgi

PRIMERA INCLUIR LOS COEFICIENTES EN

ITERACION EL SISTEMA DE ECUACIONES
SI

NO RESOLVER SISTEMA DE
ECUACIONES
SI ACTUALIZAR ?, Pgi
|€|<TOLERANCIA

NO CALCULAR
€=Pgi(iter) – Pgi(iter-1)
GUARDAR
VOLVER A RESULTADOS
CALCULAR ? DEFINITIVOS
|€|<TOLERANCIA
NO
FIN SI
GUARDAR
RESULTADOS
Figura 8. Algoritmo del método lambda iterativo. DEFINITIVOS

5 Análisis de resultados
Los resultados se analizan en relación a FIN
cuál de los métodos muestra el despacho
económico más eficiente, en este caso es
Figura 9. Algoritmo del método Jacobiano.
el costo incremental el valor que dicta
cual método logró obtener la cantidad de El primer ejercicio consiste en un SEP
potencia para cada generador a la cual el de cuatro barras ilustrado en la Figura 10.
valor del costo fue menor que otro el Los costos incrementales de cada
método. Para el análisis de los métodos generador son:

13
 𝑑 𝐹1 Para el análisis de los resultados del
= 0.0080 ∙ 𝑃𝑔1 + 8 [$/𝑀𝑊ℎ]
𝑑𝑃𝑔1 sistema dado, y la evaluación de los
métodos utilizados para del DEP, se hace
𝑑 𝐹2 la comparación de los valores obtenidos
= 0.0096 ∙ 𝑃𝑔2 + 6.4 [$/𝑀𝑊ℎ]
𝑑𝑃𝑔2 en cada método en el proceso realizado
en la aplicación de Matlab. Al haber
GEN 1
SLACK 3 utilizado el mismo caso de estudio con
CARGA
los mismos datos de entrada, los
1 A resultados obtenidos se muestran en la
Tabla 4.
El resultado de estos métodos muestra
4 2 al método del Lambda iterativo como el
GEN 2
CARGA método más eficiente dando como
B
resultado un valor de λ menor que el
Figura 10. Sistema de prueba para DEP
método del Jacobiano, e incluso, menor
En las siguientes tablas se muestran que el despacho sin pérdidas. Además,
los parámetros eléctricos del sistema de este método muestra un reparto en la
prueba. potencia de generación casi igual para
cada generador.
Tabla 1. Datos de líneas en [p.u.]
Tabla 4. Comparación de resultados
Barra i Barra j R X B
1 4 0.00744 0.0372 0.0775 Potencias de Costo
1 3 0.01008 0.0504 0.1025 generación Incremental
Método
2 3 0.00744 0.0372 0.0775 [MW] [$/MWh]
2 4 0.01272 0.0636 0.1275 Gen 1 Gen 2 λ
Tabla 2. Datos de barras en [p.u.] DES 272.72 227.27 10.18
Lambda
Barra PG |V| ϴ [°] PC QC Tipo 268.08 241.67 10.14
Iterativo
1 - 1 0 - - 1 Jacobiano 369.14 140.18 13.91
2 3.18 1 - - - 2
3 - - - 2.2 1.3634 3
4 - - - 2.8 1.7352 3
El método del Jacobiano, se muestra
en este caso de estudio, como el menos
Los valores dados están dados en eficiente al dar como resultado el valor
[p.u.] con 230 [kV] y 100 [MVA] como más alto de λ. El reparto de potencia en
valores base. El tipo de barra tiene los generadores tiene una gran
asignación numérica y corresponde a: 1- diferencia, tanto así que al primer
Barra Slack, 2-Barra PV, 3-Barra PQ. generador se le asigna más del doble de
Como se mencionado, para la potencia que al segundo generador.
resolución del DEP se parte de un cálculo Los resultados en el DES son
previo de flujo de carga. similares al de fórmula de pérdidas, el
valor de λ difiere apenas por centésimas
Tabla 3. Generación y voltajes resultantes del flujo de y la distribución de potencia en los
carga en [p.u.] generadores lleva la misma tendencia, el
Barra P Q |V| ϴ [°] primer generador es el que despacha más
1 1.913 1.8722 1 0 potencia, pero se le asigna una cantidad
2 3.18 1.3254 1 2.44
3 0.96051 -1.08 mayor que en el caso del método del
4 0.94304 -2.63 Lambda iterativo. La Figura 11 y Figura
Total 5.093 3.1977 12 ilustran de forma gráfica las
afirmaciones realizadas en cuanto al

14
costo incremental y el reparto de Costos incrementales de generadores:
potencias de generación en cada método
implementado. 𝑑 𝐹1
= 0.016 ∙ 𝑃𝑔1 + 3.2 [$/𝑀𝑊ℎ]
𝑑𝑃𝑔1

𝑑 𝐹2
= 0.0092 ∙ 𝑃𝑔2 + 4.5 [$/𝑀𝑊ℎ]
𝑑𝑃𝑔2

Parámetros eléctricos del SEP.

Tabla 5. Datos de líneas en [p.u.], 5 barras.

Barra i Barra j R X B
1 2 0 0.2 0
1 3 0 0.15 0
1 4 0 0.1 0
2 3 0 0.1 0
2 5 0 0.15 0
Figura 11. Costo incremental resultante en cada
5 4 0 0.1 0
método

Tabla 6. Datos de barras en [p.u.]

Barra PG |V| ϴ [°] PC QC Tipo

1 5.5 1 0 1.5 0.5 1
2 3.5 0.96 4.41 1.5 0.4 2
3 - 1.1 -4.43 2.5 1 3
4 - 1.04 -5.11 2.5 1 3
5 - 0.97 -4.6 1 0.7 3

Los valores dados están dados en

[p.u.] con 230 [kV] y 100 [MVA] como
valores base. El tipo de barra tiene
asignación numérica y corresponde a: 1-
Barra Slack, 2-Barra PV, 3-Barra PQ.
Figura 12. Potencia despachada en cada método
Se aplica el mismo procedimiento que
Como segundo ejercicio se plantea un el ejercicio anterior y se obtiene los
SEP de cinco barras ilustrado en la resultados que se muestran en la Tabla 7.
Figura 13.
Tabla 7. Comparación de resultados
Potencias de Costo
generación Incremental
Método
SLACK GEN 2 [MW] [$/MWh]
Gen 1 Gen 2 λ
1 2 DES 380.15 519.84 9.28
Lambda
380.15 519.84 9.28
Iterativo
Jacobiano 410.94 489.05 6.48
3
4 CARGA 5
Con el uso de la Figura 14 y Figura 15
CARGA A CARGA se hace la comparación de los resultados.
B C Para este segundo caso los resultados
Figura 13. Sistema de prueba para DEP muestran una mejor eficiencia en el
método del Jacobiano, el cual determina

15
un costo mucho menor al de los primeros con muchos procesos y operaciones
métodos y un reparto de carga más dentro de los lazos de repetición,
equilibrado entre unidades de consumirán gran cantidad de memoria de
generación. las computadoras y el tiempo de
respuesta será más extenso.
Una vez ejecutado los algoritmos, se
han encontrado los valores de las
potencias que se debería generar en cada
unidad generadora para que la
producción de energía eléctrica sea lo
menos costosa posible, pero el resultado
de cada método tiene una pequeña
variación una de otra. Lo cual ha
permitido evaluar las eficiencias de cada
método y de esta manera proponerla
como mejor opción.
Figura 14. Costo incremental resultante en cada Por el análisis comparativo entre los
método métodos presentados en este trabajo, el
método del Lamba iterativo se muestra
como la mejor opción para realizar el
cálculo del despacho económico, además
que, este método es el que menos
cantidad de operaciones contiene dentro
del lazo principal iterativo en la
programación del Matlab.

6.1 Trabajos futuros

El análisis expuesto en este trabajo puede
servir como punto de partida para un
trabajo más extenso y minucioso en
Figura 15. Potencia despachada en cada método busca de los métodos más eficientes para
llevar a cabo la tarea del despacho
6 Conclusiones económico, con investigaciones como:
Al existir varios métodos matemáticos  Despacho económico con flujos de
con los cuales se pude estimar el potencia DC.
despacho de generación termoeléctrica,  Despacho económico considerando
es importante saber diferenciar cuál de los costos variables de operación.
los métodos es el más eficiente para el  Despacho económico en sistemas
sistema que se desea analizar. Y aunque que incluyan generación con energías
cada método presenta su dificultad en el renovables.
momento de formular el modelo  Despacho económico hidro-térmico.
matemático, se pude diferenciar los
métodos iterativos, observando los 7 Referencias
procesos que intervienen en el lazo
[1] D. P. Kothari and I. J. Nagrath,
principal iterativo, ya que, para sistemas
“Modern Power System
de mayor tamaño, el realizar el cálculo
Analysis,” Third Edit., New

16
 Delhi: The McGraw Hill, 2009, 2000.
pp. 1–5. [11] A. J. Wood, B. F. Wollenberg, and
[2] A. Sánchez and H. Aviña, “El G. B. Sheblé, Power Generation,
impacto de los combustibles en los Operation, and Control, Third
costos de generación edit. New Jersey: John Wiley &
termoeléctrica,” Univerisdad Sons, Inc., 2014.
Nacional Autónoma de México, [12] H. Saada, Power system analysis.
2016. TATA McGraw-Hill, 2002.
[3] R. J. Buist and P. G. Lau, [13] Z. Xu, Q. Jia, and L. Li, “The
“Thermoelectric Power Generator approximate expression of power
Design and Selection from TE flow Jacobian matrix and
Coolling Module Specifications,” analysis,” 1st Int. Conf. Sustain.
TETechnology,Inc., no. 0, 1997. Power Gener. Supply,
[4] H. A. Cardona, M. A. Burgos, J. SUPERGEN ’09, pp. 1–4, 2009.
W. González, I. A. Isaac, and G. J. [14] R. Jegatheesan, N. M. Nor, and M.
López, “Aplicación en Matlab F. Romlie, “Newton-raphson
para la Programación del power flow solution employing
Despacho Económico systematically constructed
Hidrotérmico,” vol. 6, no. 2, pp. jacobian matrix,” PECon 2008 -
42–53, 2012. 2008 IEEE 2nd Int. Power Energy
[5] D. Glover, M. Sarma, and T. Conf., no. PECon 08, pp. 180–
Ovebye, Power System Analysis 185, 2008.
and Desing. USA: Thomson, [15] F. J. Vinueza, “Flujo Optimo de
2008. Potencia,” Universidad
[6] S. Affijulla and S. Chauhan, “A Politécnica Nacional, 1984.
New Intelligence Solution for [16] Lof of A, G. Anderson, and T.
Power System Economic Load Smed, “Fast Calculation of a
Dispatch,” Dep. Electr. Eng. NIT Voltage Stability Index,” IEEE
Hamirpur , India, 2011. Trans Power Syst., p. 7(1):54-64,
[7] M. E. El-Hawary and G. S. 1992.
Christensen, Optimal Economic [17] W. R. Barcelo, “Dynamic
Operation of Electric Power Economic Dispatch Using The
System. New York: Academic, Extended Security Constrained
1979. Economic Dispatch Algorithm Fir
[8] J. J. Grainger and W. D. Stevenson +,” vol. 12, no. 2, pp. 961–967,
Jr., Análisis de Sistemas de 1997.
Potencia. México: McGRAW- [18] R. Jegatheesan, N. M. Nor, and M.
HILL, 2002. F. Romlie, “Power Flow Solution
[9] D. Cerjan, M., D, Marcic, M, Using Impedance Parameters,”
“Short term power system Second Int. Power Eng. Optim.
planning with water value and Conf. Malaysia, 2008.
energy trade optimization.,” in [19] W. F. Reinoso P and J. L. Méndez,
Energy Market (EEM), 8th “Despacho Económico de
International Conference on the Potencia Diaria Aplicando al
European, 2011, pp. 269–274. Sistema Nacional Ecuatoriano,”
[10] R. Mujal, Tecnología Eléctrica. Universidad Politécnica
Terrasa. España: Capitulo X, Salesiana, 2011.

17
[20] B. Stott, “Review of Load-flow
Calculation Methods,” IEEE
Proceedingsvol. 62, pp., pp. 916–
929, 1974.
[21] G. B. Salazar, “Factores de Nodo
del SNI, Métodos y Análisis de
Sensitividad,” Universidad
Politécnica Nacional, 2000.

18
7.1 Matriz de Estado del Arte
Tabla 8. Matriz de estado del arte

DESPACHO ECONÓMICO DE GENERACIÓN TERMOELÉCTRICA USANDO MÉTODOS DE SENSIBILIDAD AC

FORMULACION DEL RESTRICCIONES DEL ALGORITMO DE SOLUCIÓN

DATOS TEMÁTICA
PROBLEMA PROBLEMA OPTIMIZACIÓN PROPUESTA

Límites de operación de Centrales Generadoras

Reparto de Potencia en Generadores

Despacho Económico de Generación

Despacho Económico con Pérdidas

Despacho Económico Simplificado

Flujo de carga Newton-Raphson

Costos de Producción Eléctrica

Multiplicadores de Lagrange
TÍTULO DEL ARTÍCULO

Centrales Termoeléctricas
Generación Eléctrica

Fórmula de perdidas
Balance de potencia
Costo Incremental
Función de Costos

Lambda Iterativo

Factores de Nodo
Jacobiano

Jacobiano
CITAS
ITEM

AÑO

1 2009 Modern Power System Analysis 8111   

2 2016
El impacto de los combustibles en los costos de generación
termoeléctrica
0      
3 1997
Thermoelectric Power Generator Design and Selection from TE
Coolling Module Specifications
32 
4 2012
Aplicación en Matlab para la Programación del Despacho
Económico Hidrotérmico
1    
5 2008 Power System Analysis and Desing 1432   
6 2011
A New Intelligence Solution for Power System Economic Load
Dispatch
19   
7 1979 Optimal Economic Operation of Electric Power System 313   
8 2002 Análisis de Sistemas de Potencia. 488         
9 2011
Short term power system planning with water value and energy
trade optimization
7  
10 2014 Power Generation, Operation, and Control 9579       
The approximate expression of power flow Jacobian matrix and
11 2009
analysis
3  
Newton-raphson power flow solution employing systematically
12 2008
constructed jacobian matrix
14   
13 1992 Fast Calculation of a Voltage Stability Index 509  

19
 Dynamic Economic Dispatch Using The Extended Security
14 1997
Constrained Economic Dispatch Algorithm
69  
15 2008 Power Flow Solution Using Impedance Parameters 2  
16 2011
Despacho Económico de Potencia Diaria Aplicando al Sistema
Nacional Ecuatoriano
2    
17 1974 Review of Load-flow Calculation Methods 809 
18 2000 Factores de Nodo del SNI, Métodos y Análisis de Sensitividad 0 
CANTIDAD: 2 4 2 3 6 7 3 2 2 2 2 1 1 2 1 6 7 5

7.2 Resumen de Indicadores

TEMÁTICA
Generación Eléctrica
4
3,5
3
2,5
2
1,5
1
0,5
Reparto de Potencia en
0 Centrales Termoeléctricas
Generadores

Costos de Producción
Eléctrica

Figura 16. Resumen e indicador de la temática - Estado del arte

20
 FORMULACIÓN DEL PROBLEMA
8

0
Despacho Económico de Generación Despacho Económico Simplificado Despacho Económico con Pérdidas

Figura 17. Indicador de formulación del problema - Estado del arte

SOLUCIÓN DADA MEDIANTE

8
7
6
5
4
3
2
1
0
Multiplicadores de Lagrange Flujo de carga Newton- Jacobiano Factores de Nodo
Raphson

Figura 18. Indicador de solución - Estado del arte

21