lOMoARcPSD|42593754

Aplicación EN Matlab PARA LA Programación

Catedra Integradora: Políticas Públicas SEPTIMO CICLO y Proyectos Sociales

(Universidad Nacional de Loja)

Escanea para abrir en Studocu

Studocu no está patrocinado ni avalado por ningún colegio o universidad.

Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

Revista Investigaciones Aplicadas | ISSN 2011-0413 | Medellín - Colombia

Vol. 6, No.2 (2012) Julio – Diciembre | PP.42-53
http://revistas.upb.edu.co/index.php/investigacionesaplicadas

APLICACIÓN EN MATLAB PARA LA PROGRAMACIÓN DEL DESPACHO ECONÓMICO

HIDROTÉRMICO

Hugo A. Cardona*†, Miguel A. Burgos **, Jorge W. González*, Idi A. Isaac*, Gabriel J. López*
*
Universidad Pontificia Bolivariana, Cir. 1 #70-01, Blq. 11, Medellín, Colombia.
**
PC Mejia S. A., Cra. 52 # 85-42, Itagüí, Colombia.

Recibido 15 Julio 2012; aceptado 29 Noviembre 2012

Disponible en línea: 21 Diciembre de 2012

Resumen: En este artículo se presenta una aplicación computacional para llevar a

cabo el despacho económico en un sistema de energía eléctrica utilizando
MATLAB ®. El proceso consiste en reducir al mínimo el costo de generación
hidro-térmica utilizando la técnica de los multiplicadores de Lagrange; e
identificando los parámetros de cada generador, tales como el consumo de energía
máximo y mínimo, el volumen máximo y el mínimo, la ecuación de costos, la
ecuación de flujo, la eficiencia, y el consumo total de potencia por día. La caja de
herramientas de interfaz gráfica de usuario de MATLAB se utiliza con el fin de
mostrar el proceso y los resultados de una manera didáctica. La aplicación también
se puede utilizar como una herramienta pedagógica para la enseñanza y el
aprendizaje del proceso y los términos utilizados en el Despacho Económico.
Además, se desarrolló un manual del usuario para lograr un adecuado de esta
aplicación.

Palabras clave: Despacho Económico, MATLAB, Multiplicadores de Lagrange,

Generación Hidro-térmica.

Abstract: In this article it is presented an application for performing the Economic

Dispatch in an electrical power system using MATLAB®. The process consist of
minimizing the cost of hydro-thermal generation using the technique of Lagrange’s
multipliers; and identifying parameters of each generator like maximum and
minimum power consumption, maximum and minimum volume, cost equation, flux
equation, efficiency; and the total power consumption per day. The MATLAB®
toolbox GUI is used in order to show the process and results in a didactic way. The
application can also be used like a pedagogic tool for teaching and learning the
process and terms used in the Economic Dispatch. Also, a User Manual was made
for a good use to this application.

Keywords: Economic Dispatch, MATLAB, Lagrange Multipliers, Hydro-Thermal

Generation.

†
Autor al que se le dirige la correspondencia:
Tel. (+574) 4488388.
E-mail: [email protected] (Hugo A. Cardona).

Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])

 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

1. INTRODUCCIÓN 2.2. Control de un Sistema de Potencia

El estudio de sistemas eléctricos por medio de La demanda de potencia a lo largo de un día será
herramientas computacionales ha avanzado en variable dependiendo del día de la semana
gran medida en las últimas décadas por lo cual es considerado, de la estación del año, e incluso de
prioritario en las universidades y en este entorno la ubicación geográfica del sistema. A pesar de
hacer uso de ellas para la resolución de problemas estas condiciones, la curva de carga no variará
que atañen a diario el campo de acción de la mucho, presentando picos y valles. En la Fig. 1 se
ingeniería eléctrica (Hosseinzadeh et al., 2011). puede distinguir un valor mínimo, la carga base y
En este caso concreto, como es el Despacho un valor máximo, la punta de carga.
Económico, es de vital importancia disponer de
métodos de análisis y herramientas
computacionales para el correcto funcionamiento
del sistema energético de cualquier país (Agrawal
y Porate, 2010), (Mariano et al., 2009).

Debido a la gran cantidad de variables

involucradas, a la complejidad de las relaciones
matemáticas de optimización requeridas y a la
inclusión del componente de costos de
generación, un programa para el estudio del Figura 1. Curva genérica de demanda de potencia
despacho hidrotérmico es conveniente en pro del diaria
mejoramiento en la calidad de enseñanza en
ingeniería eléctrica. Con el presente proyecto se 2.3. Gestión de Sistemas Térmicos
usa la herramienta MATLAB® para la
elaboración de un sistema capaz de realizar el En el caso de la gestión a corto plazo se deben
despacho hidrotérmico de varios generadores en tener presente ciertas restricciones que
función de una demanda establecida. contemplen las características de operación de
cada central:
Se hace necesario conocer las aplicaciones del
software para el cálculo de las diferentes - Mínimos tiempo de arranque y parada
variables como potencia, carga y tiempos,
teniendo en cuenta las múltiples restricciones que - Costos de arranque en función del tiempo de
los generadores presentan, y la interfaz gráfica de parada
usuario.
- Comportamiento de la subida de carga
2. MARCO TEÓRICO - Mínimos parámetros técnicos

2.1. Despacho Económico El objetivo de la optimización a corto plazo

consiste en decidir los arranques y paradas de
El conocimiento del flujo de cargas en un sistema cada uno de los generadores térmicos. Dadas las
eléctrico de potencia permite hallar la potencia características de la curva de demanda a cubrir, el
activa y reactiva que debe entregar cada unidad ciclo básico de operación resulta ser semanal, por
generadora para atender una demanda de potencia tanto, el problema del corto plazo debe ser
determinada. El reparto de cargas o potencias resuelto con horizonte semanal.
entre generadores depende de las condiciones de
operación que se impongan. El problema que se 2.4. Modelo de sistemas térmicos
plantea es el siguiente: conocida la demanda de
potencia total en el sistema, averiguar la potencia En una central térmica el generador eléctrico
que debe entregar cada unidad para que el costo convierte en energía eléctrica la energía mecánica
total de la generación sea mínimo. A esto se le entregada por la turbina. El aporte de vapor a la
conoce como despacho económico (Mujal, 2000), turbina es suministrado de diferentes formas
(Glover et al., 2008), (Cerjan et al. 2011). según se trate de una central térmica
convencional (en tal caso el vapor será generado
en una caldera), o de una central térmica nuclear
rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 43
Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

(el vapor será generado en el reactor de fisión). modo que la relación puede representarse con la
En las centrales térmicas convencionales se debe siguiente ecuación:
pagar un alto precio por la obtención de
combustible (Glover et al., 2008). ∂C
λ= = AP + B
El modelo económico de una central térmica está ∂P (2)
dado por una entrada que puede ser el costo de la
cantidad de combustible necesario para generar la Donde A representa los costos variables por la
potencia solicitada, que es la característica de producción que dependen de la potencia de salida,
salida, y teniendo en cuenta ciertas restricciones, y B son los costos fijos.
como el valor máximo y mínimo de generación
de cierta turbina, entre otros.
2.6. . Criterio para optimizar la distribución
2.5. Distribución de cargas entre diversas
unidades térmicas Utilizando el método de multiplicadores de
Lagrange se puede encontrar el criterio a tener en
La Fig. 2 muestra una curva típica de entrada- cuenta para una distribución óptima de la carga
salida en una central térmica: entre diferentes unidades.

Se supone que la carga total va a ser repartida

entre k unidades. El costo total y la potencia total
de la central serán las sumas de los costos o de las
potencias de todas las unidades en
funcionamiento:

k
CT = C 1 + C 2 + … + C k = ∑ C n ,
n=1 (3)
Fig. 2. Curva típica de una unidad generadora k
térmica. Tomado de Mujal (2000). PT = P1 + P2 + … + Pk = ∑ Pn ,
n=1 (4)

El costo del combustible se calcula fácilmente. El

Donde C T : costo total de la central.
criterio para distribución de carga entre dos
C n : costo de cada unidad en funcionamiento.
unidades se basa en, que si se aumenta la carga en
la primera unidad, disminuyendo en la misma P T : Potencia total consumida
proporción en la segunda unidad, resulta un P n : Potencia de cada unidad en funcionamiento.
aumento o disminución en el costo total. Así, se
está tratando la variación del costo que se
determina por las pendientes de las curvas El objetivo es obtener un valor mínimo de C T
características. para un P T dado, lo que exige que el diferencial
total sea dF T =0. Puesto que el costo del
C n : Entrada a la unidad n (costo por hora) combustible depende de la salida de potencia de
P n : Salida a la unidad n (potencia) cada unidad.

La variación del costo de combustible de la El criterio para el reparto económico de la carga

unidad será: entre diferentes unidades de una central, es que
∂C todas ellas deben funcionar con la misma
λ= variación del costo de combustible (Mujal,
∂P (1)
2000).
Normalmente, la variación del costo de
Se puede calcular el valor de λ de la siguiente
combustible puede aproximarse a través de una
manera:
función lineal, sin generar mayores errores, de tal

rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 44

Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

Se sabe que para cada unidad: salto neto y es por tanto igual al salto bruto
menos las pérdidas de carga durante las
∂C n conducciones.
λn = = An Pn + Bn ,
∂Pn (5)

por tanto:

λn − B n
Pn = ,
An (6)

entonces:

k k
λn − Bn
PT = ∑ Pn = ∑ ,
n=1 n=1 An (7) Fig. 3. Modelo de Central Hidroeléctrica.
Tomado de Garrido (1986).
Si λ1 = λ2 = …= λk = λ entonces:
La potencia que se obtiene en el alternador queda
P expresada como:
λ= k T ,
1 − Bn
∑n=1 An
(8) P = KQH n ηa ηt , (10)

Como An y Bn son parámetros conocidos de cada Donde:

unidad la potencia a suministrar por cada una se
puede calcular utilizando: P : potencia en el alternador.
K : factor de conversión.
λ − Bn Q : caudal turbinado.
Pn = , Hn : Salto neto.
An (9) nt : rendimiento de turbina.
na : rendimiento del alternador.
Esta relación se puede definir como P = f ( Q,
Hn )
2.7. Modelo de sistemas hidráulicos

Un modelo sencillo de central hidroeléctrica 2.8. Gestión hidrotérmica a corto plazo

consiste en un embalse limitado por una presa y
una central en la que se produce la energía como La gestión hidrotérmica a corto plazo cubre la
muestra la Fig. 3. La primera característica que se programación con detalle horario y horizonte
debe definir para cada central es la función de semanal y en algunos casos mensual. Se
rendimiento, es decir, la relación existente entre considera que con este horizonte son conocidas la
la potencia generada y el caudal a través de las demanda, las aportaciones y la disponibilidad de
turbinas. La energía disponible para la conversión los generadores. Las previsiones de la demanda y
eléctrica procedente del agua contenida en el de las aportaciones en diversos puntos de las
embalse es una función del salto, es decir, la cuencas se deben hacer previamente basadas en
diferencia de alturas entre la superficie de agua en datos históricos y de predicción de variables
el embalse y el nivel del agua en el desagüe meteorológicas (Correa, 2004).
central (CREG, 2012), (Garrido, 1986), (Rio et
al., 2012). Este salto se denomina salto bruto. El También son conocidos el estado inicial de los
salto disponible en la propia turbina es embalses así como los niveles deseados en cada
ligeramente menor que el salto bruto debido a las uno de ellos, al final del período en estudio. Esta
pérdidas de carga ocasionadas por fricción en la información debe provenir de estudios a medio y
toma, tubería forzada, tubo de aspiración y galería largo plazo.
del desagüe. El salto así determinado se denomina
rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 45
Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

- Problema matemático básico

Se plantea un primer modelo teniendo en cuenta

las siguientes condiciones:

- Se tiene una sola central térmica y una

hidráulica
- No hay pérdidas de caudal (Vertimiento)
- La cantidad de agua a gastar es fija.
- El salto neto constante.
Fig. 5. Acople de centrales para el cumplimiento
En la Fig. 4 se indican las variables a tener en de la demanda. . Tomado de Garrido (1986).
cuenta en el modelo de la hidroeléctrica:

El problema consiste en minimizar el costo total

de operación térmica:

jmax
minFT = ∑n F ,
j=1
j j
(11)
Fig. 4. Variables en el modelo de generación
hidráulica. Tomado de Garrido (1986). Sujeto a las siguientes dos condiciones:

jmax
Donde:
j : etapa.
∑n q
j=1
j j = qTOT , (12)

rj : aporte durante la etapa j.

Vj : volumen al final de la etapa j. (12) es el gasto total de agua.
qj : caudal turbinado durante la etapa j.
sj : vertido durante la etapa j. PLj − PHj − PTj = 0, (13)
En la Fig. 5 se observa el modelo circuital de la
potencia generada y entregada, además se definen (13) es el balance de potencia para cada etapa.
las siguientes variables para plantear la solución:
Se supone el funcionamiento con una diferencia
Tmax : número total de horas del período. de niveles constante (salto constante) y por tanto
nj : número de horas hasta la etapa j. que se dispone de la relación caudal turbinado –
Fj : Costo horario de operación de la central potencia q = q (PH).
térmica en la etapa j.
PHj : potencia generada por la central hidráulica Para hallar el mínimo de la función de costo se
durante la etapa j. construye la función de Lagrange:
PTj : potencia generada por la central térmica

[∑ ]
jmax jmax
durante la etapa j.
L= ∑ [n j F (P Tj )+λ j (P Lj− P Hj− PTj ) ]+ γ n j q j (P Hj )− q TOT
PLj : demanda a cubrir en la etapa j. j=1 j= 1
(14)

Si se deriva (14) con respecto a las dos variables

independientes y se iguala a cero para cumplir la
condición de mínimo se obtendrá para una cierta
etapa j = k:

rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 46

Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

∂L dF (PTk )
= 0 → nk = λk ,
∂PTk dPTk (15)
∂L dq(PHk )
= 0 → γnk = λk ,
∂PHk dPHk (16)

A estas ecuaciones se les conoce como

ecuaciones de coordinación. Se resuelven a través
de métodos numéricos, utilizando software de
cálculo como MATLAB en el caso de este
trabajo.

3. METODOLOGÍA

El despacho económico a implementar tiene las Fig. 6. Tabla de consumo

siguientes características:

- Entradas: Otra forma de hacerlo es a través de la generación

o Número indeterminado de unidades de una curva polinómica, indicando el consumo
térmicas o hidráulicas. base, las horas valle y las horas pico de consumo.
o Curva de consumo diario. Esto se hace utilizando la función polyfit. Se
o Curvas de costo para las unidades obtiene una curva como la que muestra la Fig. 7:
térmicas. 1000

o Curvas de flujo para embalses. 900

o Restricciones de potencia máxima y 800

mínima para generadores térmicos. 700

o Restricciones de volúmenes iniciales y 600

finales en embalses. 500

400
o Eficiencia de generadores. 300

200

- Salidas 100

o Distribución horaria de la carga total para 0

0 5 10 15 20 25

cada unidad.
o Energía consumida y costo para las
Fig. 7 Curva de consumo creada por el usuario
unidades térmicas.
o Simulación animada del funcionamiento
del despacho, y datos con resolución de
cinco minutos 3.2. Distribución en Unidades Térmicas

Dada las condiciones reales y consideraciones

especiales como tiempos de sincronización, las
3.1. Curva de consumo diario
unidades térmicas son difíciles de variar, de modo
que se utilizan para manejar la carga mínima
Para la curva de consumo diario se almacena en
requerida por la central durante el período del
un vector los datos ingresados por el usuario. En
despacho, mientras que la carga variable se deja
la Fig. 6 se muestra un ejemplo de la tabla que
para las unidades hidráulicas que son más
debe llenar el usuario.
controlables. De esta manera, el valor mínimo del
consumo se distribuye entre las diferentes
unidades térmicas de acuerdo con el criterio de
minimización de costos. Para esto es necesario
que el usuario ingrese los parámetros necesarios
para determinar la curva de costo, y la respectiva
variación.

rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 47

Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

El usuario ingresa los parámetros a, b y c en la 3.3. Distribución entre las unidades

siguiente ecuación cuadrática para cada unidad. hidráulicas

C = aP 2 + bP + c, (17) La carga restante caracterizada por ser variable y

menor, puede ser atendida por las unidades
los datos son almacenados en una matriz, donde hidráulicas, teniendo en cuenta que el volumen de
las filas representan cada unidad térmica y las agua utilizado no exceda un valor máximo. Uno
columnas los valores a, b y c respectivos. Al de los mayores problemas de esta distribución,
realizar la derivada: como se mostraba en la sección anterior, es el
desconocimiento del flujo de agua del afluente,
que este caso será ingresado por el usuario como
dC un valor constante durante el período. Sin
= 2aP + b,
dP (18) embargo, será adicionado a éste una variable
estocástica para simular un proceso aleatorio, lo
de donde se puede concluir que A = 2a, y B = b. cual será más real.
De esta manera se puede calcular el valor de λ y
la potencia suministrada por cada unidad: En este caso no existe una restricción de costo.
Por tanto, se supone que no existe un problema de
Lambda = Potencia / (cumsum((1-b)/2a)) optimización con respecto a esta variable. Ahora
potenciastermicas = (lambda-b)/2a con respecto al volumen, se utiliza la siguiente
relación cuadrática del flujo del embalse con
En esta distribución es importante revisar que la respecto a la potencia suministrada:
potencia asignada a las unidades no sobrepase el
máximo permitido o no alcance el valor mínimo Q = aP 2 + bP + c, (19)
requerido para su funcionamiento. En MATLAB,
se realiza primero la distribución y luego se
verifica esta condición. En caso que se exceda la teniendo en cuenta que el volumen se puede
calcular como el flujo multiplicado por el tiempo,
potencia máxima, esta unidad se deja trabajando
en su valor máximo y el excedente de potencia se y el tiempo se mide en horas (que es la unidad del
distribuye de acuerdo con el criterio del menor período), el volumen puede expresarse de la
costo (no variación). Por el contrario, si la siguiente forma:
potencia no alcanza el valor mínimo, entonces ese
generador no se conecta a la red y se redistribuye V = aP 2 + bP + c, (20)
la potencia entre las unidades restantes.
Vi sería el volumen de agua gastado por un
En la Fig. 8 se muestra un diagrama de flujo de la embalse en una hora, entregando una potencia
distribución de carga entre las unidades térmicas. determinada. Con esta relación, que es similar a la
del costo en el caso de la distribución térmica, se
puede calcular la potencia entregada por cada
unidad teniendo en cuenta los valores para los
cuales la variación del volumen es igual para cada
unidad:

dV
= 2aP + b,
dP (21)

La gran diferencia en esta distribución radica en

que mientras en las unidades térmicas el consumo
total no varía durante el período, para las
unidades hidráulicas, este cálculo debe realizarse
para cada hora, y además teniendo en cuenta que
Fig. 8. Diagrama de flujo distribución de carga el volumen entregado al final del período no
entre unidades térmicas puede exceder la diferencia entre el volumen final

rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 48

Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

e inicial ingresados por el usuario, además que el 4. RESULTADOS Y ANÁLISIS

volumen actual del embalse también dependerá
del flujo del afluente. Como se mencionó
anteriormente, como se trata de la unidad de Debido a restricciones de información, por su
tiempo la cantidad de volumen es equivalente a la carácter confidencial, no fue posible disponer de
magnitud del flujo: modelos matemáticos y parámetros de las
unidades térmicas e hidráulicas en Colombia.
Dada esta imposibilidad, sólo se tratan casos
Vi = aP 2 + bP + c + flujoentrada, (22) hipotéticos que permitan verificar el buen
funcionamiento de los algoritmos propuestos.
en este caso la programación debe ser dinámica
debido a que el valor del volumen de restricción Para la obtención de resultados y para el
irá disminuyendo o aumentando dependiendo de respectivo análisis, se realizará una simulación
la potencia utilizada en la hora anterior. teniendo en cuenta los datos de las Tablas 1 y 2.
Se consideran tres unidades térmicas y tres
Vi = Vi −1 + flujoentrada − flujosalida, (23) unidades hidráulicas. Se aclara que las ecuaciones
asignadas no corresponden a un sistema físico
para cada hora se realiza la distribución de cargas, específico sino que están adecuadas para la
además se debe comprobar que el volumen obtención de resultados que demuestren la
requerido no exceda el valor máximo. En caso de capacidad de la aplicación.
ser así, es necesario disminuir la potencia para esa
unidad, y ser asignada a otras unidades, Tabla 1. Especificaciones unidades térmicas
garantizando que el volumen de las demás no sea
PotMáx PotMín
excedido. Esto se realiza bajo el criterio del Unidad Tipo
(MW) (MW)
Ecuación %Ef
menor gasto de agua.
UT1 T 500 20 P2+0.2P+100 0.86
Para el caso de la distribución hidráulica se 2
UT2 T 600 30 3P +2P+30 0.91
realiza el diagrama de flujo de la Fig. 9:
UT3 T 750 45 2P2+0.5P+60 0.9

Tabla 2. Especificaciones unidades hidráulicas

Ini Fin
Unidad Tipo Ecuación Flujo %Ef
Hm3 Hm3
UT1 H 5000 300 3P2+4.5P+3 60 0.95

UT2 H 4000 500 0.2P2+3P+5 85 0.93

UT3 H 5000 500 0.3P2+3.2P+0 80 0.86

Se agregan las unidades y el siguiente paso es

ingresar la curva de consumo. Para este caso se
genera una curva automática teniendo en cuenta
un consumo base de 400MW y una tendencia a
aumentar a las 12:00 horas. La Fig. 10, muestra la
curva y los generadores agregados.

Una vez realizada la curva de consumo, el

siguiente paso es realizar el despacho. Ver Fig.
11.

Hasta aquí ya se ha cumplido el objetivo de la

aplicación. Se realizó la respectiva distribución de
Fig. 9. Diagrama de Flujo para distribución unidades manera que para las unidades
hidráulica hidráulicas la carga es variable y en la mayoría de
los casos de menor valor que la carga manejada
rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 49
Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

por las unidades térmicas (que además se La distribución de la carga base entre los
mantienen constantes durante todo el período). generadores térmicos ha sido realizada de manera
Ver Fig. 11. óptima de tal manera que esa combinación genera
un costo mínimo en el uso de combustibles. Este
cálculo permite la comparación con otras
combinaciones y comprobar que efectivamente es
la más óptima, tal como se observa en la Tabla 3.
Nótese que incluso al asignar toda la carga al
generador que produce menos costos, al final el
costo total sigue siendo mayor que el valor
obtenido con los multiplicadores de Lagrange.

Tabla 3. Comprobación del costo mínimo

UT1 UT2 UT3 Costo 1 Costo 2 Costo 3 Costo Total

193 63 96 37.387 12.063 18.540 67.990

100 200 52 10.120 120.430 5.494 136.044

52 100 200 2.814 30.230 80.160 113.204

300 25 25 90.160 1.955 1.322 93.437

352 0 0 124.074 30 60,00 124.164

0 352 0 100 372.446 60,00 372.606

0 0 352 100 30 248.044 248.174

Luego de realizar el despacho, vienen las

opciones adicionales de la aplicación, como lo es
la simulación, en la cual se aumenta la resolución
del tiempo a 5 minutos extrapolando los datos
Fig. 10. Curva de Consumo horarios; y se simulan fluctuaciones en la carga
debido a la aleatoriedad de los consumidores y a
la eficiencia en el funcionamiento de los
generadores. En la Fig. 12 se muestra la
simulación para el consumo total y la generación
de la planta UH1.

Fig. 12. Simulación

Fig. 11. Tabla de Despacho

rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 50

Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

Una vez terminada la simulación se pueden

observar todos los datos obtenidos, en la
ventana Resultados, en donde se puede
comparar gráficamente el comportamiento de
todos los consumos durante el período. Por
ejemplo, para el consumo de UT2, se observa
en la Fig. 13 cómo arranca desde cero (como
era de esperarse) y se mantiene oscilando
(línea verde) alrededor del valor óptimo
calculado (línea blanca). Las fluctuaciones se
deben a la variabilidad en el consumo y
generación.
Figura 14. Consumo de UH2 durante la
simulación

Fig. 13. Consumo UT2 durante Simulación

Por otro lado es posible observar también las Figura 15. Consumo de UH1 durante la
gráficas para la generación hidráulica. Para UH2 simulación
por ejemplo, existe mucha variabilidad aunque
con potencia menor a la de los generadores
A trav de varias tablas, la aplicación muestra la
térmicos. Nótese en la Fig. 14 el efecto de la
totalidad de datos utilizados durante el
restricción de volumen de agua alrededor de la
funcionamiento de la aplicación (programación
hora 9, en donde la potencia y por ende el flujo
del despacho y simulación) y es posible su
de salida es disminuido. Sin embargo, debido al
exportación a hojas de cálculo para análisis
afluente para las horas siguientes, es posible que
posteriores. Además que se cuenta con las
retorne a generar mayor potencia. Al observar la
posibilidad de guardar en el disco duro (o
Fig. 15, esa disminución de potencia entregada
cualquier otro dispositivo) los parámetros de las
por parte de la unidad UH2 es asignada a la
unidades creadas durante una simulación, para
unidad UH1, que presenta un aumento en ese
luego ser utilizadas nuevamente sin necesidad de
mismo tiempo.
realizar el ingreso nuevamente. Se dispone de la
presentación gráfica del comportamiento de la
potencia entregada por cada una de las diferentes
unidades durante y después de la simulación. El
uso de barras de herramientas permite opciones
de acercamiento en las gráficas, selección de
muestras de coordenadas e incluso impresión.

rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 51

Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

5. CONCLUSIONES Correa, C.M. (2004). Desarrollo de una

herramienta para ejecutar el despacho
Se presentó una aplicación computacional para económico Colombiano en condiciones de
llevar a cabo el despacho económico en un contingencia. Tesis de Grado UPB.
sistema de energía eléctrica utilizando MATLAB.
El proceso consiste en reducir al mínimo el costo CREG Comisión de Regulación de Energía y Gas
de generación hidro-térmica utilizando la técnica (2007). En línea, en [http://www.creg.gov.co],
de los multiplicadores de Lagrange; e consultado en 2012-02-01.
identificando los parámetros de cada generador,
tales como el consumo de energía máximo y Garrido, J.A. (1986). Coordinación Hidrotérmica,
mínimo, el volumen máximo y el mínimo, la Fundamentos Básicos y Métodos para la
ecuación de costos, la ecuación de flujo, la solución del Problema. Universidad
eficiencia, y el consumo total de potencia por día. Intemacional Menéndez Pelayo. Santander.
Como herramienta didáctica, la aplicación
desarrollada fue diseñada de tal modo que tuviera Rio, J.C., L. Bayon, E. Fernandez (2012). El
una forma visualmente agradable y de fácil problema de la coordinación hidrotérmica en el
acceso a los datos. La aplicación incluye además nuevo mercado eléctrico. En línea, en
el uso de los conceptos del despacho económico [http://www.revistadyo.com/index.php/dyo/arti
como estrategia para la motivación y fácil cle/viewFile/263/263]. Consultado en 2012-11-
comprensión del proceso como tal y el 15.
significado de términos como curvas de consumo,
Glover, D., M. Sarma, T. Overbye (2008). Power
distribución de cargas, variación del costo, curvas
Systems Analysis and Design. Thomson, USA.
de costos, restricciones de potencia, restricciones
de volumen. Hosseinzadeh, N. M.R.Hesamzadeh y M. Korki
(2011). Role of laboratory in the education of
Dentro de los problemas encontrados en el
modern power systems. Power and Energy
trabajo, estuvo la dificultad de realizar
Society General Meeting, 2011 IEEE: 1-6.
simulaciones de distribuciones de carga de
situaciones reales debido al escaso acceso a Mariano, S.J.P.S. M.R.A Calado, L.A.F.M
información sobre los modelos matemáticos y Ferreira (2009). Dispatch of Head Dependent
parámetros de las unidades térmicas e hidráulicas Hydro Units: Modeling for optimal generation
en Colombia, datos absolutamente necesarios in electricity market. PowerTech, 2009 IEEE
para la programación del despacho en la Bucharest: 1- 6.
aplicación. Motivo por el cual se trataron casos
hipotéticos que permitieran verificar el buen Mujal, R. (2000). Tecnología Eléctrica. Terrasa.
funcionamiento de los algoritmos propuestos. Capítulo X. España.
El toolbox de interfaz gráfica de usuario de
MATLAB, resultó ser una herramienta útil para el
manejo de resultados y gráficos. SOBRE LOS AUTORES

Hugo A. Cardona
REFERENCIAS
Ingeniero Electricista, MSc.y estudiante PhD. Es
Agrawal, S.P., K.B Porate (2010). Economic Profesor Titular de tiempo completo e
Dispatch of Thermal Units with the Impact of Investigador de la Universidad Pontificia
Wind Power Plant. Emerging Trends in Bolivariana, Colombia. Ha sido pasante
Engineering and Technology (ICETET), 2010: investigador en XM ESP.
48- 53.
Miguel A. Burgos.
Cerjan, M., D. Marcic, M. Delimar (2011). Short
term power system planning with water value Egresado de Ingeniería Eléctrica. Ha apoyado al
and energy trade optimization. Energy Market grupo de investigación en TyD de la UPB. Desde
(EEM), 8th International Conference on the el 2010 labora en PC Mejía S.A. dedicado al
European: 269- 274. diseño y montaje eléctrico.

rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 52

Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])
 lOMoARcPSD|42593754

H. Cardona, M. Burgos, J. González, I. Isaac, G. López | Aplicación en Matlab para la programación…

Jorge Wilson González

Ingeniero Electricista, MSc - PhD. Es Profesor

Titular de tiempo completo e Investigador de la
Universidad Pontificia Bolivariana, Colombia.
Laboró en HMV Ingenieros y en Siemens,
Erlangen en PTD (FACTS y HVDC).
Investigador invitado en la Univ. Kempten,
Alemania y en el Power Systems Institute en
Erlangen Univ. Alemania en el 2005.

Idi A. Isaac

Ingeniero Electricista, MSc - PhD. Es Profesor

Titular de tiempo completo e Investigador de la
Universidad Pontificia Bolivariana, Colombia.
Investigador invitado en la Univ. Kempten,
Alemania y Univ. Comillas, España.

Gabriel López

Ingeniero Electricista, MSc.y estudiante PhD. Es

Profesor Asociado de tiempo completo e
Investigador de la Universidad Pontificia
Bolivariana, Colombia. Investigador invitado en
la Univ. Kempten, Alemania.

rev.investig.apl | Vol. 6, No.2 (2012) Julio – Diciembre | p. 53

Descargado por DEYBIT JAIR QUEVEDO HERRERA ([email protected])