Distribución Óptima de Carga Entre Unidades Generadoras.

1
V. Cervantes Verdin, D. Gmez Martnez, Estudiantes de Maestra en Ciencias U. de G. CUCEI
Distribucin ptima de Carga entre Unidades

Generadoras dentro de una Planta y entre
Plantas
Resumen: El presente proyecto tiene como objetivo emplear

despacho econmico para determinar la distribucin ptima de
carga entre las unidades de una planta y entre las plantas de un
sistema de potencia. Se describen dos mtodos iterativos (Lambda
iterativo y Newton) que permiten encontrar la distribucin ptima
de carga entre las unidades de una planta con el mnimo costo
incremental. Para determinar la distribucin ptima entre plantas
se expone un caso de estudio y se resuelve empleando el software
PowerWorld.
Palabras Clave Despacho econmico, Operacin econmica
del SEP.
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
una planta con dos unidades de generacin cuyas funciones de
entrada salida son:
H 1=8 P1 +0.024 P21 +80
(1)
H 2=6 P2 +0.04 P22 +120
(2)
Las curvas entrada-salida para (1) y (2) son:
I. INTRODUCCIN
l despacho econmico ha sido ampliamente

estudiado en las ltimas dcadas debido a la
necesidad de encontrar la operacin ms econmica
del sistema elctrico de potencia SEP. Actualmente
la mayor parte de la generacin de energa elctrica proviene
de unidades trmicas (Carbn, combustleo etc.) cuyo costo
de sus insumos ha incrementado progresivamente. El empleo
de combustibles fsiles adems de significar un alto costo
implica la emisin de gases contaminantes, por lo que resulta
de suma importancia buscar un ptimo aprovechamiento de
los combustibles.
El despacho econmico adems obtener la distribucin
ms econmica de las cargas entre generadores contribuye a
reducir la emisin de agentes contaminantes. En este trabajo
se presentan los mtodos lambda iterativo y Newton para la
obtencin del punto ptimo de operacin econmica de las
unidades de una planta [1]. Tambin se emplear el software
PowerWorld para determinar el punto ptimo de operacin
entre las unidades de un SEP compuesto por 3 generadores y 3
cargas.
Fig. 1. Curva entrada-salida para (1)
Resulta de utilidad la obtencin de las curvas

caractersticas de las unidades de generacin para visualizar
con mayor facilidad su comportamiento de entrada-salida. Por
lo tanto un buen punto de partida es la obtencin de las curvas
de entrada salida, proporcin de calor y costo incremental de
combustible.
II. CURVAS CARACTERSTICAS PARA LAS UNIDADES DE
GENERACIN DE UNA PLANTA
Primeramente se obtendrn las curvas caractersticas para

Universidad de Guadalajara CUCEI, M.C en Ingeniera
Elctrica, Operacin y control de Sistemas Elctricos
(Victor_Cervantes91@outlook.com,
Daniel1790890@hotmail.com )
Fig. 2. Curva entrada-salida para (2)
En las figuras 1 y 2 se muestra la relacin entre la entrada
de potencia de combustible y la salida de la unidad de

generacin en MW.
En las figuras 3 y 4 se muestra la relacin que existe entre
la entrada de combustible en BTU y la salida de energa en
KWh en funcin de la potencia de salida en MW.
Fig. 5. Curva de costo incremental para (1)
Fig. 3. Curva de proporcin de calor para (1)
Fig. 6. Curva de costo incremental para (2)
III. DESPACHO ECONMICO SIN PRDIDAS PARA UNIDADES

GENERADORAS DENTRO DE UNA PLANTA
Fig. 4. Curva de proporcin de calor para (2)
Las figuras 5 y 6 muestran las curvas de costo incremental

obtenidas a partir de las ecuaciones (1) y (2) escaladas por el
costo de combustible y derivadas con respecto a P 1 y P2,
ecuaciones (3) y (4).
En este caso el despacho econmico busca obtener el

punto ptimo de operacin econmica, que se logra cuando
las unidades operan a un mismo costo incremental.
En este caso de estudio se tienen 3 unidades de generacin
dentro de una planta con una generacin total de 450 MW, sus
ecuaciones de entrada-salida son las siguientes:
2
H 1=0.072 P1+12
(3)
H 2=0.12 P2+ 9
(4)
H 1=( 8.4 P1+ 0.0025 P1 +225 )
MBtu
h
(5)
H 2=(6.3 P2 +0.0081 P2 +729)
MBtu
h
(6)
3
2
H 3=( 7.5 P3 +0.0025 P3 +400 )
MBtu
h
(7)
'
H 2=(6.3 P2 +0.0081 P2 +729)1.02
R
MBtu
(9)
Las curvas caractersticas correspondientes para las

unidades de generacin descritas por (5), (6) y (7) se muestran
a continuacin:
H '3=( 7.5 P3 +0.0025 P23 +400 )0.90
R
MBtu
(10)
Derivando las funciones de costo (8), (9) y
obtienen las funciones de costo incremental:
(10)
se
d H1
R
==( .004 P1 +6.72 )
d P1
MWh
(11)
d H2
R
==(.0165 P2 +6.426)
d P2
MWh
(12)
d H3
R
==( .0045 P3+ 6.75 )
d P2
MWh
(13)
En las ecuaciones (11), (12) y (13) se puede ver que

es el costo incremental comn de las unidades. El valor de
y las potencias de cada unidad se encuentran al hacer el
despacho econmico de la planta.
Fig. 7. Curvas de entrada-salida para las unidades (5), (6) y (7).
Las restricciones de generacin para cada unidad se

muestran en la tabla 1:
A. Mtodo lambda iterativo.
TABLA I
RESTRICCIONES DE GENERACIN DE (5), (6) Y (7)
Generado
r
Pmin
Pmax
H1
H2
H3
45
350
45
350
47.5
450
Es un mtodo iterativo que busca obtener un costo

incremental ptimo para las unidades de la planta. En
este mtodo se parte de un valor inicial propuesto y
mediante iteraciones se ajusta al valor ptimo.
El siguiente diagrama de flujo muestra el procedimiento
para obtener el valor ptimo de a partir de un valor
inicial propuesto.
La solucin al despacho econmico no deber violar los

lmites permisibles para las unidades de generacin mostrados
en la tabla 1.
Si se escalan las funciones de entrada-salida (5), (6) y (7)
por el costo de combustible en (R/MBtu) se obtienen las
funciones de costo como sigue:
H '1=( 8.4 P1+ 0.0025 P21 +225 )0.8
En este documento se utilizan los mtodos lambda iterativo

y Newton para hacer el despacho econmico de la planta, los
mtodos se describen a continuacin:
R
MBtu
(8)
La matriz Hessiana
H=
Lx
x
contiene las
segundas derivadas de la funcin de Lagrange
vector gradiente
Fig. 8. Diagrama de flujo del mtodo iterativo
L se define como sigue:
L . El
[]
(16)
El vector de incrementos es
x=xx 0 , siendo
L
P1
L
P2
L=
L
P3
L
Es aconsejable partir de un valor

coherente cercano
al valor ptimo para agilizar la convergencia del mtodo.
x 0 el vector de datos iniciales:
Para este caso de estudio se elabor un cdigo de Matlab

para determinar el valor ptimo de y las potencias de
cada unidad. El cdigo se basa en el mtodo iterativo .
La corrida del programa arroja los siguientes resultados:
[]
P1
P
x 0= 2
P3
0
TABLA II
SOLUCIN AL DESPACHO ECONMICO POR LAMBDA ITERATIVO
#It
P1
P2
P3
1
2
3
8
7.2
7.5437
320
120
205.925
95.25
46.84
67.6573
277.77
100
176.37
7
-243.033
183.15
0
(17)
En ste mtodo se puede proponer el valor de lambda

inicial que se desee, puesto que el mtodo converge muy
rpidamente.
Siguiendo el algoritmo propuesto en (14) se encuentra el
vector solucin x .
En la tabla 2 se puede observar que en la tercera iteracin

el error se vuelve 0, por lo tanto se llega a la solucin
en 3 iteraciones.
TABLA III
SOLUCIN AL DESPACHO ECONMICO POR NEWTON
B. Mtodo iterativo Newton.

El mtodo de Newton es una variante del mtodo del
gradiente, con la diferencia de que la matriz Jacobiana se
transforma en una matriz Hessiana. En este mtodo los
incrementos se calculan a partir de la siguiente expresin:
x=
En donde L
est dada por:
Lx
x
(14)
es la funcin de Lagrange sin prdidas y
i=1
i=1
L= F i ( Pi)+ ( PLOAD Pi )
(15)
# Ite.
1
2
0
7.5436
P1
150
205.905
2
P2
100
67.7346
P3
200
176.37
se vuelve [0],
por lo tanto el mtodo converge en dos iteraciones.
A partir de la segunda iteracin
5
SOLUCIN AL DESPACHO ECONMICO CASO A.
P1
P2
P3
7.29
143 MW
52 MW
120 MW
Las figuras 11, 12 y 13 muestran el costo incremental

correspondiente de cada unidad de generacin.
10,00
7,50
R/MWh
5,00
Fig. 9. Despacho econmico para la planta. Todas las unidades operan al

mismo costo incremental.
2,50
El mtodo de Newton requiere de un mayor nmero de

instrucciones por iteracin en comparacin con lambda
iterativo, sin embargo el nmero de iteraciones es menor. Por
lo tanto se puede concluir que ambos mtodos presentan
demandas de recursos computacionales similares.
0,00
0
175
Generator Power (MW)
263
350
Fig. 11. Curva de costo incremental para la unidad de generacin 1.
IV. DESPACHO ECONMICO ENTRE PLANTAS PARA UN SEP
16,0
CON Y SIN PRDIDAS
En la figura 10 se muestra el SEP que se compone de 3

plantas generadoras y 3 cargas, el despacho econmico
pretende encontrar la distribucin ptima de carga entre las
plantas. Se emplear el software PowerWorld para encontrar
la solucin ptima.
88
12,0
R/MWh
8,0
4,0
0,0
0
88
175
263
350
GEN 3
GEN 2
10,00
7,50
R/MWh
5,00
2,50
GEN
1
0,00
Fig. 10. Diagrama de impedancias para un sistema de nueve nodos: todas

las impedancias estn en pu con respecto a una potencia base de 100 MVA y
un voltaje base de 13.8, 16.5, 18 y 230 kV.
Las ecuaciones de los generadores 1, 2 y 3 corresponden a

(5), (6) y (7) y sus correspondientes curvas de entrada salida
son las mostradas en la figura 7.Se harn tres simulaciones
resolviendo despacho econmico para los siguientes tres
casos:
A. Sin considerar restricciones de generacin ni prdidas.
Al hacer la simulacin para el SEP se obtiene una
=7.29 para las tres mquinas y las potencias mostradas
en la tabla III.
TABLA IV
113
225
338
450
B. Considerando restricciones de generacin sin prdidas.

De la simulacin se obtiene una
=7.29 y las
potencias de cada unidad son las mismas que las mostradas en
la tabla IV. Esto se debe a que an con restricciones de
generacin ningn lmite de generacin (tabla I) es violado.
C. Considerando restricciones de generacin y prdidas.
Al hacer la simulacin para el SEP se obtiene una
=7.3 para las tres mquinas y las potencias mostradas
en la tabla V.
TABLA V
SOLUCIN AL DESPACHO ECONMICO CASO C.
P1
P2
P3
7.3
144 MW
53 MW
122 MW
Como se puede observar en la tabla V, al considerar las

prdidas aumenta la generacin y el costo incremental se eleva
ligeramente con respecto a los casos anteriores.
V. CONCLUSIONES
Las curvas caractersticas de las unidades de generacin
son una herramienta grfica til que permite conocer para una
determinada potencia de salida el costo de generacin.
Los mtodos iterativos expuestos en este documento son
slo algunos entre muchos mtodos que existen para hacer

despacho econmico por planta. Ambos mtodos presentan
sus particularidades tales como complejidad y nmero de
iteraciones. El mtodo de Newton implica menor nmero de
iteraciones que el mtodo lambda iterativo pero requiere un
mayor nmero de instrucciones por iteracin. En cuanto a
exactitud ambos convergen en la misma solucin.
Al encontrar la distribucin ptima entre plantas sin
considerar restricciones de transmisin, el costo incremental
es el mismo para las tres plantas. Al incluir las prdidas
de transmisin la generacin en las tres plantas aumenta, lo
que implica un aumento en el valor de .
VI.
[1]
REFERENCIAS
Allen J. Wood, Power Generation Operation and Control, Edit. WileyInterscience, Chapter 3 Economic Dispatch of Thermal Units and
Methods of Solution, 1996.

Distribución Óptima de Carga Entre Unidades Generadoras.

Cargado por

Copyright:

Formatos disponibles

Distribución Óptima de Carga Entre Unidades Generadoras.

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Distribución Óptima de Carga Entre Unidades Generadoras.

Cargado por

Copyright:

Formatos disponibles

1

V. Cervantes Verdin, D. Gmez Martnez, Estudiantes de Maestra en Ciencias U. de G. CUCEI

Distribucin ptima de Carga entre Unidades

Resumen: El presente proyecto tiene como objetivo emplear

H 1=8 P1 +0.024 P21 +80

H 2=6 P2 +0.04 P22 +120

Las curvas entrada-salida para (1) y (2) son:

l despacho econmico ha sido ampliamente

Fig. 1. Curva entrada-salida para (1)

Resulta de utilidad la obtencin de las curvas

Primeramente se obtendrn las curvas caractersticas para

Fig. 2. Curva entrada-salida para (2)

En las figuras 1 y 2 se muestra la relacin entre la entrada

de potencia de combustible y la salida de la unidad de

Fig. 5. Curva de costo incremental para (1)

Fig. 3. Curva de proporcin de calor para (1)

Fig. 6. Curva de costo incremental para (2)

III. DESPACHO ECONMICO SIN PRDIDAS PARA UNIDADES

Fig. 4. Curva de proporcin de calor para (2)

Las figuras 5 y 6 muestran las curvas de costo incremental

En este caso el despacho econmico busca obtener el

H 1=( 8.4 P1+ 0.0025 P1 +225 )

H 2=(6.3 P2 +0.0081 P2 +729)

H 3=( 7.5 P3 +0.0025 P3 +400 )

H 2=(6.3 P2 +0.0081 P2 +729)1.02

Las curvas caractersticas correspondientes para las

H '3=( 7.5 P3 +0.0025 P23 +400 )0.90

En las ecuaciones (11), (12) y (13) se puede ver que

Las restricciones de generacin para cada unidad se

A. Mtodo lambda iterativo.

Es un mtodo iterativo que busca obtener un costo

La solucin al despacho econmico no deber violar los

H '1=( 8.4 P1+ 0.0025 P21 +225 )0.8

En este documento se utilizan los mtodos lambda iterativo

segundas derivadas de la funcin de Lagrange

Fig. 8. Diagrama de flujo del mtodo iterativo

L se define como sigue:

Es aconsejable partir de un valor

x 0 el vector de datos iniciales:

Para este caso de estudio se elabor un cdigo de Matlab

En ste mtodo se puede proponer el valor de lambda

En la tabla 2 se puede observar que en la tercera iteracin

B. Mtodo iterativo Newton.

es la funcin de Lagrange sin prdidas y

Las figuras 11, 12 y 13 muestran el costo incremental

Fig. 9. Despacho econmico para la planta. Todas las unidades operan al

El mtodo de Newton requiere de un mayor nmero de

Fig. 11. Curva de costo incremental para la unidad de generacin 1.

IV. DESPACHO ECONMICO ENTRE PLANTAS PARA UN SEP

CON Y SIN PRDIDAS

En la figura 10 se muestra el SEP que se compone de 3

Fig. 12. Curva de costo incremental para la unidad de generacin 2.

Fig. 10. Diagrama de impedancias para un sistema de nueve nodos: todas

Las ecuaciones de los generadores 1, 2 y 3 corresponden a

Fig. 13. Curva de costo incremental para la unidad de generacin 3.

B. Considerando restricciones de generacin sin prdidas.

Como se puede observar en la tabla V, al considerar las

slo algunos entre muchos mtodos que existen para hacer