Despacho Economico

Despacho Econmico
Harold Salazar Isaza, Ph.D

Programa de Ingeniera Elctrica
Universidad Tecnolgica de Pereira
Abril 2009
Harold Salazar Isaza
Contenido

Introduccin a el despacho econmico

Modelaje de las curvas de costo
Formulacin matemtica del problema de despacho
econmico
Tcnicas de optimizacin (Ver presentacin Introduccin a
las tcnicas de optimizacin)

Solucin al problema de despacho econmico

Ejemplos
2
Introduccin
El despacho econmico es una tcnica ampliamente utilizada en
operacin de los sistemas de potencia. El objetivo es establecer
despacho optimo (despacho de menor costo) de un conjunto
generadores para una determinada demanda sujeto a un conjunto
restricciones operativas.
la
el
de
de
Esquematicamente:
Pg1=?
d1
Pgi=?
Red de
Transmision
Pgn=?
dk
Introduccin
De que depende el despacho econmico?
El costo de generacin de electricidad de las

distintas tecnologas
Se requiere el modelo
2.
La demanda
estos componentes para formular
3.
Las restricciones operativas
matemticamente el problema de
1.
matemtico de cada uno de
despacho econmico
Modelaje de la curva de costo

Los costos de generacin se deben principalmente a los siguientes
factores:
1. Costo de construccin
2. Costo operativos
Costo del combustible (energas primarias) necesarias para

producir electricidad
Costo de mano obra
Costo de mantenimiento
Estos costos son

lo que
primordialmente
determinan los
costos
operativos

Analisis comparativo de el costo promedio de diferentes energias primarias en los E.U
Fuente: http://www.eia.doe.gov/cneaf/electricity/epa/epat4p5.html
Observe que los costos de las energas primarias se pueden expresar en $/MMBtu
(Mayor informacin: http://www.eia.doe.gov/emeu/mer/append_a.html)
Los costos de las diapositivas anteriores reflejan los costos del

combustible. Estos costos no indican los costos de producir electricidad.
Como modelar matemticamente los costos de produccin?
Puesto que el costo de produccin esta definido principalmente por el costo

de combustible, la curva de costo solo tomara en consideracin este costo
para su construccin.

Como modelar el costo de generacin?
Curva de entrada-salida: Esta curva indica cuanto combustible por hora
es necesario para colocar a disposicin una cantidad determinada de MW
Entrada: Se puede expresar en

trminos de:
1. Pesos / hora
2. Toneladas de carbn / hora
3. Millones de m3 de gas / hora
4. Energa / hora
Salida: Potencia neta en MW

disponible para el sistema de
potencia

Curva de entrada-salida: Esta curva indica los requerimientos energticos
en MBTU o J (Joule) por unidad de tiempo (generalmente hora) para
colocar a disposicin una determinada potencia.
Entrada (MBTU/hr)
Entrada (kJ/hr)
Esta curva puede ser obtenida segn los datos de
diseo o por pruebas de laboratorio. Por ejemplo:
Para una planta de carbn, se puede establecer el
numero de toneladas por hora para generar un
valor dado en MW. Luego se puede multiplicar el
poder calorfico del carbn (MBTU/ton) para
establecer el valor de entrada (eje y)
Salida: Pg (MW)

Funcin de Produccin: Esta curva indica la potencia disponible como funcin de los
requerimientos energticos en MBTU o J (Joule) por hora.
Salida: Pg (MW)
Observe que a medida que se incrementa la

entrada no existe un incremento sustancial de
la salida. Este efecto es conocido, en
economa, como la ley de la disminucin de la
produccin marginal.
Esta disminucin al margen se explica por el
aumento de escape de energa de entrada a
medida que la temperatura aumenta.
Entrada (kJ/hr)
Entrada (MBTU/hr)
10

MWh se puede transformar
en MBTU para obtener la
eficiencia en MBTU/MBTU.
Ver ejemplo.
"Eficiencia" =
Energia de salida
MW
MWh
=
=
Energia de entrada MBTU / h MBTU
(MWh/MBTU)
Bajos niveles
de eficiencia
Nivel mas
eficiente de
generacin
Pg (MW)
Normalmente este es el valor nominal

de la unidad de generacin
11

Curva de proporcin de calor: (Heat rate curve) Esta curva es similar a la curva de
eficiencia excepto que los valores en el eje y se encuentran invertidos. Esto es, el eje y
contiene los valores que corresponden a el inverso de la eficiencia (curva de la diapositiva
anterior). La curva es denotada H(Pg). Esta curva es especialmente til para derivar la
curva de costo.
1/ = H(Pg)
(MBTU/MWh)
Esta curva es una caracterstica
Nivel mas eficiente

de generacin
Pg (MW)
Normalmente este es el valor nominal

12

Curva de costo: (Cost curve) La curva de costo se construye a partir de la curva de
proporcin de calor de la siguiente manera:
1.
Defina R como la razn a la cual la unidad de generacin utiliza combustible por hora,
esto es, la curva entrada-salida que es igual a la curva de proporcin de calor H(Pg)
multiplicada por Pg.
R ( Pg ) = Pg H ( Pg )
( MBTU / h )= ( MW ) x ( MBTU / MWh )

2.
Defina C como el costo de produccin establecido como el costo del combustible K

multiplicado por la razn a la cual la unidad de generacin utiliza combustible R
C ( Pg ) = K R( Pg )
( $ / h )= ( $ / MBTU ) x ( MBTU / h )
13

Curva de costo: (Cost curve) La curva de costo indica el costo por hora para una determinada
potencia de salida. Esta curva depende de las caractersticas de la unidad, la cual, estn
indicadas por la curva de proporcin de calor (heat reat curve) y el costo del combustible.
Dos caractersticas de la curva de costo:
C(Pg)
( $/ h)
Pg (MW)
1.
El costo incrementa con la generacin. Esta

es una caracterstica esperada pues a
mayor nivel de generacin se requiere mas
combustible incrementando de esa manera
el costo de generacin.
2.
La curva, por lo general y para efectos de

anlisis, se considera convexa. Esto es, el
segmento de lnea que une dos puntos
cualquiera de la curva SIEMPRE esta por
encima de la curva.
14

Curva de costo incremental: (Incremental cost curve) El costo incremental es el costo de
producir un kW adicional o kWh si la potencia se considera constante durante una hora. Se
define, matemticamente, como la derivada de la curva de costo, esto es:
CI ( Pg ) =
Note las unidades del
costo incremental $/MWh
dC ( Pg )
dPg
CI(Pg)
( $/ kWh)
Pg (kW)
Dos observaciones de la curva de costo

incremental:
1.
La curva de costo incremental es

creciente lo cual indica que el costo
de generar un MW adicional (o MWh
si la potencia es constante durante
una hora) incrementa con la
potencia.
2.
La curva de costo incremental es

lineal para curvas de costo
cuadrticas.
15

Ejemplo:
Tomado de http://powerlearn.ee.iastate.edu/asp/
Una unidad trmica de generacin de 100MW utiliza carbn cuyo contenido energtico es
12,000 BTU/lb. El costo del carbn es $1.5/MBTU (Note las unidades de costo del
carbn!!). Para la curva de demanda indicada en la tabla inferior, calcular:
1. La eficiencia para cada nivel de demanda.
2. La curva de proporcin de calor (Heat Rate Curve) en MBTU/MWh.
3. El costo por hora.
Hora
Pg (MW)
Toneladas de Carbon
empleada
12:00am-6:00am
40
105.0
6:00am-10:00am
70
94.5
10:00am-4:00pm
80
156.0
4:00pm-12:00am
100
270.0
16

Solucin
1. La eficiencia para cada nivel de demanda
Sea T el nmero de horas para producir P (MW) usando y toneladas de
carbn. Recuerde que 1W = 1J/s
Energia de salida
Energia de entrada
sec
watts
106
hr
MW
=
BTU
lb
joules
12000
2000
1054.85
ytons
lb
ton
BTU
P T 3600
17

2. La curva de proporcin de calor
La curva de proporcin de calor es la cantidad de MBTU usados en un
tiempo T divido por los MWh de salida en el mismo intervalo T
12,000
H=
BTU
lb
1MBTU
2000
ytons 6
lb
ton
10 BTU
P T
18

3. El costo por hora
Del punto 2 del problema
(diapositiva anterior)
R( Pg ) = Pg H ( Pg )
C ( Pg ) = K R( Pg ) = 1.5 R( Pg )
Costo del combustible en $/MBTU
19

Ecuacin
dispositiva 17
Solucin
Ecuacin
dispositiva 18
Ecuacin
dispositiva 19
T (hrs)
P (MW)
y (tons)
H (MBTU/MWh)
C ($/hr)
40
105.0
0.33
10.5
630
70
94.5
0.42
8.1
850
80
156.0
0.44
7.8
936
100
270.0
0.42
8.1
1215
20

La curva de costo se puede aproximar con una curva cuadrtica.
1300
fitted curve
1200
C(pg)=0.09033Pg2-2.955Pg+604.9
1100
Costo ($/h)
% Codigo en Matlab
Pg = [40 70 80 100]';
Costo = [630 850 936 1215]';
f = fittype('poly2');
c = fit(Pg,Costo,f);
plot(c);
xlim([0 120])
xlabel('Pg (MW)');ylabel('Costo ($/h)')
hold on
plot(Pg,Costo,'x')
1000
900
800
700
600
20
40
60
Pg (MW)
80
100
120
21

El costo incremental es la derivada de la curva de costo.
0.1807 x-2.955
16
14
C(pg)=0.09033Pg2-2.955Pg+604.9
CI(pg)=0.1807Pg-2.955
12
Costo Marginal ($/MWh)
% Codigo en Matlab
ezplot('0.1807*x-2.955',[20 100])
xlabel('Pg (MW)');ylabel('Costo Marginal ($/MWh)')
10
8
6
4
2
0
20
30
40
50
60
Pg (MW)
70
80
90
100
22

En general, la curva de costo de las generadores trmicos se
aproxima, matemticamente, a la siguiente expresin cuadrtica:
C ( Pg ) = aPg 2 + bPg + c
En donde:
a, b, y c son coeficientes que caracterizan al generador trmico
23
Formulacin matemtica del problema

El problema del despacho econmico se puede formular de la
siguiente manera:
Minimizar El costo TOTAL de la generacin
Sujeto a
1) Balance global: El despacho de todos los generadores debe ser igual a la
demanda total incluyendo las perdidas.
2) Limites de generacin: El valor despachado de cada generador debe estar
dentro de sus limites de operacin definido por sus valores mximo y
mnimo.
3) Capacidad de transmisin: El flujo de potencia por cada lnea del sistema
debe ser menor que su limite mximo de operacin.
24

Matemticamente:
Por simplicidad no se consideran, inicialmente, las restricciones que impone la
red de transmisin.
min
Pgi
C ( Pg )
i =1
s.a:
Pg = Pd
i =1
Funcin Objetivo
j =1
+ PL
Pg imin Pgi Pg imax
Restricciones
i = 1 n
25

Definicin de las variables:
E n donde:
Ci
(P g i ):
C o s to d e l g e n e ra d o r i
P g i : V a lo r d e s p a c h a d o p a r a e l g e n e r a d o r i
Pd
: V a lo r d e la d e m a n d a j
P L : P e r d i d a s t o t a le s
Pg
m in
i
, Pg
m ax
i
: L im ite s d e g e n e ra c i n d e l g e n e ra d o r i
n : N m e ro to ta l d e g e n e ra d o re s
m : N m e ro to ta l d e c a rg a s o d e m a n d a s
26
El anterior problema se puede resolver empleando

tcnicas de programacin no-lineal. La presentacin
Introduccin a la Optimizacin presenta las tcnicas

comnmente empleadas en la solucin de este
problema.
27
Solucin al problema de despacho

econmico
El problema de despacho econmico se resolver de acuerdo a la siguiente
secuencia:
1. Despacho econmico ignorando perdidas y limites de generacin

2. Despacho econmico ignorando perdidas y considerando limites de
generacin
3. Despacho econmico considerando perdidas y limites de generacin
28

Ignorando prdidas y limites de generacin
Despacho econmico sin considerar perdidas e ignorando limites de
generacin.
Sistema el cual ignora la
red de transmisin
Pg1=?
Modelo matemtico
n
d1
min
Pgi
Pgi=?
s.a:
Pgn=?
dm
C ( Pg )
i =1
Pg = Pd
i =1
j =1
29

Modelo matemtico
Funcin objetivo
min
Pgi
C ( Pg )
i =1
s.a:
j =1
Restriccin
m
n
L( Pg1 , , Pg n , ) = Ci ( Pgi ) Pgi Pd j

i =1
j =1
i =1
Pg = Pd
i =1
Funcin de LaGrange
Multiplicador de LaGrange
30

m
n
L( Pg1 , , Pg n , ) = Ci ( Pgi ) Pgi Pd j

i =1
j =1
i =1
Funcin de LaGrange
CdPO
1. Derivadas con respecto a las potencias generadas:
L
=0
Pgi
=
Pgi Pgi
i = 1 n
m
n
n
C
Pg
Pg
Pd
i ( i )
i j = CI i ( Pgi ) = 0 i = 1 n
j =1
i =1
i =1

2. Derivada con respecto al multiplicador de LaGrange:
L
=0
Pg Pd
i =1
j =1
=0
31

El sistema de ecuaciones de la diapositiva anterior implica que el despacho
econmico es aquel en donde los costos marginales (o costos incrementales)
de todos los generadores son iguales.
CI i ( Pg i ) = 0 CI i ( Pg i ) = i = 1 n
n
Pg Pd
i =1
j =1
=0
Un sistema de n+1
ecuaciones lineales*
con n+1 incgnitas!!
* Si las funciones de
costo son cuadrticas
32

La interpretacin grfica para un sistema de tres generadores es la siguiente:

CI(Pg)
CI1
CI2
CI3
Igual costo incremental :
Pg1
Pg2
Pg3
Pg1+Pg2+Pg2=PD
Pg
33

Como solucionar el conjunto de ecuaciones de la dispositiva 32?
1. Si el sistema es pequeo (nmero reducido de generadores) y las curvas de
costo son cuadrticas entonces se puede resolver como un sistema de la forma
Ax=b.
Caso tres generadores:
CI1 ( Pg1 ) = 0 2a1 Pg1 + b1 = 0
CI 2 ( Pg 2 ) = 0 2a2 Pg 2 + b2 = 0
CI 3 ( Pg3 ) = 0 2a3 Pg3 + b3 = 0
Pg1 + Pg 2 + Pg3 = PD
34

Caso tres generadores:
CI1 ( Pg1 ) = 0 2a1 Pg1 + b1 = 0

Sistema de
ecuaciones
CI 2 ( Pg 2 ) = 0 2a2 Pg 2 + b2 = 0
CI 3 ( Pg3 ) = 0 2a3 Pg3 + b3 = 0
Pg1 + Pg 2 + Pg3 = PD
A
Sistema de
ecuaciones en
forma matricial
2a1
0
2a2
0
1
0
2a3
1
1 Pg1 b1
1 Pg 2 b2
=
1 Pg3 b3

0 PD
35

2. El problema se puede resolver en trminos generales con el mtodo iterativo
lambda. Este mtodo no requiere, como el caso anterior, que las curvas de costo sean
cuadrticas.
Mtodo iterativo:
1. Asumir un valor inicial de
2. Encontrar los valores Pgi a partir del asumido. Estos valores se encuentran con las curvas de
costo incremental.
3. Si se satisface la restriccin de generacin total igual a demanda total, entonces parar. De lo
contrario:
3.1) Si la suma de la generacin es mayor que la demanda total, disminuir e ir al paso 2
3.2) Si la suma de la generacin es menor que la demanda total, aumentar e ir al paso 2.
36

Una interpretacin grfica del mtodo iterativo lambda
CI(Pg)
CI1
CI2
CI3
2
1
o
Incrementando se
incrementa el despacho
total
Pg1
Pg2
Pg3
Pg1+Pg2+Pg2=PD
Pg
37

Otros dos mtodos (no cubiertos en este curso)

empleados para la solucin del despacho econmico son:
Mtodo del gradiente
Mtodo de Newton Raphson
38

Taller para ser resuelto en clase
Encuentre el despacho econmico para los siguientes tres generadores si la demanda total es de
550MW y las perdidas y los limites de los generadores son ignorados. Encuentre el despacho de la
siguiente manera:
1.
Por el mtodo iterativo de lambda asumiendo un despacho inicial de Pg1=195MW;

Pg2=193.3MW; y Pg3=78.94MW.
2.
Resolviendo el sistema de la forma Ax=b.

Curvas de entrada-salida
Costo del combustible
R ( Pg1 ) = 0.0025 Pg12 + 8.4 Pg1 + 225
Costo de combustible G1: $0.80/MBTU
R ( Pg 2 ) = 0.0035Pg 2 2 + 6.3Pg 2 + 729
R ( Pg3 ) = 0.0025Pg 32 + 7.5Pg 3 + 400
39

El significado de lambda
Considere los tres generadores del ejercicio anterior. El despacho ptimo esta dado
por:
PG1 229.1243
P 213.1920
G2 =
PG 3 107.6836

7.6365

Pg1=?
Pgi=?
d=550 MW
Pgn=?
Costo total de generacin: C1 ( PG1 ) + C2 ( PG 2 ) + C3 ( PG 3 ) = $5186.1/ hr
40

Considere un incremento de 1MW en la demanda total. El nuevo despacho, resolviendo
este nuevo problema, esta dado por:
Pd = 550MW
Pd = 551MW
Ctotal = $5186.1/ hr
Ctotal = $5193.7 / hr
PG1 229.1243
P 213.1920
G2 =
PG 3 107.6836

7.6365

PG1 229.5367
P 213.4325
G2 =
PG 3 108.0309

7.6381
Ctotal = $5193.7 / hr
Ctotal = 5193.7 = 5186.1 + 7.63
Costo para una
demanda de
550MW
para una
demanda de
550MW
41

El operador de LaGrange , en el despacho econmico, representa el

costo incremental del sistema. Es decir, el costo de suplir 1MW
adicional durante una hora.
Nota: Este resultado se puede demostrar, analticamente, para la formulacin
general (diapositiva 25) utilizando el teorema del sobre (envelope theorem). La
siguiente diapositiva es una demostracin para el caso en donde no existen
limites de generacin.
42

Demostracin:
dCi ( PGi )
CT =
PGi
dP
i =1
Gi
C T = C i ( PGi )
i =1
CTo + CT
= C i ( PGio ) + C i ( PGi )
n
i =1
i =1
= Ci ( P
n
i =1
o
Gi
)+
n
i =1
dC i ( PGi )
PGi
dPGi
i =1
i =1
= PGi = PGi
= PD
CT Cambio en el costo total $

=
PD
Cambio en la demanda MWh
= (continua en la parte derecha)
43

Ignorando prdidas y considerando limites de generacin
Considere ahora los limites de generacin
El despacho econmico considerando limites de generacin se puede resolver a travs
de varios programas de carcter comercial (Matlab, GAMS, etc). Sin embargo, para
sistemas pequeos (un nmero muy reducido de generadores) se puede utilizar
la funcin de LaGrange para resolver este problema. Una ventaja de resolver este tipo
de problemas empleando LaGrange es que permite un mejor entendimiento del
problema. Las siguientes diapositivas muestran como solucionar este problema
empleando el procedimiento de LaGrange.
44

Modelo matemtico
n
min
Pgi
s.a:
C ( Pg )
i
i=1
n
Pg = Pd
i =1
j =1
Pg imin Pgi Pg imax
i = 1 n
Funcin de LaGrange
L( Pg , ,
max
m
n
n
n max
max
) = Ci ( Pgi ) Pgi Pd j + i
Pgi Pg i
+ imin Pg imin Pgi
i =1
j =1
i =1
i =1
i =1
n
min
45

Funcin de LaGrange
L( Pg , ,
max
m
n
n
n max
max
Pgi Pg i
+ imin Pg imin Pgi
i =1
j =1
i =1
i =1
i =1
min
Para el caso de dos generadores:
L ( Pg1 , Pg 2 , , 1min , 1max , 2 min , 2 max ) =
C1 ( Pg1 ) + C2 ( Pg 2 ) ( Pg1 + Pg 2 PD ) +
( Pg
1 ( Pg
1
max
max
max
max
Pg
Pg
Pg
)
(
)+
1
1
2
2
2
min
min
Pg1 ) + 2 min ( Pg 2 min Pg 2 )
46

Observe que si ninguna restriccin de desigualdad esta activa (el despacho para un
generador no esta en su limite de generacin) los multiplicadores de LaGrange son
iguales a cero. (Ver presentacin introduccin a las tcnicas de optimizacin)
Funcin de LaGrange
L( Pg , ,
max
m
n
n
n max
max
Pgi Pg i
+ imin Pg imin Pgi
i =1
j =1
i =1
i =1
i =1
min
Cero si Pgi es diferente a Pgimax
Cero si Pgi es diferente a Pgimin
47

Con la consideracin de la diapositiva anterior, el problema del despacho se
puede resolver con el siguiente procedimiento:
Mtodo iterativo:
1. Resolver el despacho econmico ignorando los limites de generacin con el mtodo lamdba o
como un sistema lineal. Esta suposicin asume que ninguna restriccin esta activa, es decir,
que los operadores son cero.
2. Si la solucin no viola ninguna restriccin entonces parar. La solucin es la obtenida en el
punto 1. De lo contrario ir a 3 (Ver siguiente diapositiva)
48

Mtodo iterativo:
3.
Si un limite de generacin es violado, entonces identificar la ecuacin en las condiciones KKT

que corresponde a el limite violado. En esta ecuacin, el operador de LaGrange es diferente
de cero, por tanto, la restriccin de desigualdad esta activa, es decir, es igual a cero. El
sistema se resuelve nuevamente tomando en consideracin la desigualdad activa.
4.
Una vez resuelto el sistema anterior, verificar que el despacho se encuentre dentro de los
limites, si no se viola ninguna restriccin, entonces parar, de lo contrario ir nuevamente al
paso 3.
49

Ejemplo
Considere nuevamente los tres generadores de la diapositiva 39. Cuyos costos y limites de
generacin son los siguientes.
C1 ( Pg1 ) = 0.0020 Pg12 + 6.72 Pg1 + 180
50 Pg1 150
C2 ( Pg 2 ) = 0.0034 Pg 2 2 + 6.17 Pg 2 + 714
50 Pg 2 300
C2 ( Pg3 ) = 0.0024 Pg 32 + 7.12 Pg 3 + 380
100 Pg 2 300
Recuerde el modelo matemtico

3
min
Pgi
s.a:
C ( Pg )
i =1
Pg
i =1
= 550
Pg imin Pgi Pg imax

i = 1, 2,3
50

Ejemplo
La funcin de LaGrange y las condiciones KKT para este problema son:
L( Pg , ,
max
3
n
3 max
max
) = Ci ( Pg i ) Pgi 550 + i
Pg i Pg i
+ imin Pg imin Pg i
i =1
i =1
i =1
i =1
min
KKT
L
=0
Pgi
i = 1, 2,3
L
= 0.0040 Pg1 + 6.72 + 1max 1min = 0
Pg1
L
= 0.0068Pg 2 + 6.17 + 2max 2min = 0
Pg 2
L
= 0.0048 Pg3 + 7.12 + 3max 3min = 0
Pg 3
51

Ejemplo
KKT (Continuacin)
L
=0
gi ( x
*
i
Pg1 + Pg 2 + Pg3 550 = 0
)b 0
g ( x ) b = 0
* 0
i
i = 1, 2,3
Pg1 Pg1max 0
Pg 2 Pg 2 max 0
Pg 3 Pg 3 max 0
Pg1min Pg1 0
Pg 2 Pg 2 0
min
Pg 3 Pg 3 0
min
max
max
1
2
max
max
max
min
min
( Pg Pg ) = 0
( Pg Pg ) = 0
( Pg Pg ) = 0
( Pg Pg ) = 0
( Pg 2 Pg 2 ) = 0
( Pg 3 Pg 3 ) = 0
max
min
min
min
min
52

Ejemplo
Solucin:
La solucin de este problema sin restricciones corresponde a la solucin mostrada en la diapositiva
40, la cual es reproducida a continuacin.
PG1 229.1243
P 213.1920
G2 =
PG 3 107.6836

7.6365

Observe que el limite superior del generador 1 es violado, por tanto 1max0. El nuevo sistema de
ecuaciones, bajo esta consideracin es:
53

Ejemplo
Sistema de ecuaciones considerando que el limite superior de Pg1 es violado.
0.0040 Pg1 + 6.72 + 1max = 0

0.0068 Pg 2 + 6.17 = 0
0.0048 Pg 3 + 7.12 = 0
Pg1 + Pg 2 + Pg3 550 = 0
Pg1 150 = 0
Este sistema se puede simplificar de la siguiente manera:
+ 1max = 0.0040 150 6.72 = 7.32

0.0068Pg 2 + 6.17 = 0
0.0048Pg3 + 7.12 = 0
Pg 2 + Pg3 = 550 150 = 400
54

Ejemplo
El conjunto de ecuaciones anteriores es un sistema lineal (cuatro ecuaciones y cuatro incgnitas)
que se puede expresar de la forma Ax=b cuya solucin es la siguiente:
Sin restricciones de generacin
PG1 229.1243
P 213.1920
G2 =
PG 3 107.6836

7.6365

Con restricciones de generacin
Ctotal = $5203.9 / hr
PG1 150.00
P 247.41
G2
PG 3 = 152.58
7.8524

max
1 0.5324
Observe el
incremento
de estos dos
valores
respecto al
problema sin
restricciones.
Por qu?
55

Solucin
El significado de 1max
Con restricciones de generacin
Ctotal = $5203.9 / hr
Ctotal = $5203.4 / hr
50 Pg1 150
50 Pg1 151
PG1 150.00
P 247.41
G2
PG 3 = 152.58
7.8524

max
1 0.5324
PG1 151.00
P 247.00
G2
PG 3 = 152.00
7.8496

max
1 0.5256
Ctotal = $5203.9 / hr
Ctotal = 5203.9 0.5324 5203.4
Reduccin del costo si la
restriccin se aumenta 1MW
56

1. Solucione el ejercicio anterior, es decir, encuentre el despacho econmico, si el limite superior
del generador 2 se reduce a 150MW. Que ocurre con el costo total de generacin comparado
con el caso anterior?, Cuales son los valores y la interpretacin econmica de los
multiplicadores , 1max, 2max , 3max ?
Datos del problema:
C1 ( Pg1 ) = 0.0020 Pg12 + 6.72 Pg1 + 180
50 Pg1 150
C2 ( Pg 2 ) = 0.0034 Pg 2 2 + 6.17 Pg 2 + 714
50 Pg 2 150
C2 ( Pg3 ) = 0.0024 Pg 32 + 7.12 Pg 3 + 380
100 Pg 2 300
PD = 550MW
57

Solucin
Sin restricciones
de generacin
Con una restriccin

activa de generacin
Con dos restricciones

activas de generacin
Ctotal = $5203.9 / hr
Ctotal = $ 5259 / hr
PG1 229.1243
P 213.1920
G2 =
PG 3 107.6836

7.6365

PG1 150.00
P 247.41
G2
PG 3 = 152.58
7.8524

max
1 0.5324
PG1 150.00
P
150.00
G
2
PG 3 250.00
8.3200
max 1.0000
1
max
2 1.1300
58

Considerando prdidas e ignorando limites de generacin
Considere ahora las prdidas del sistema las cuales se denotan como PL
Modelo matemtico
n
min
Pgi
s.a:
C ( Pg )
i
i=1
n
Pg = Pd
i
i =1
j =1
+ PL
Funcin de LaGrange
m
n
L( Pg1 , , Pg n , ) = Ci ( Pgi ) Pgi Pd j PL

i =1
j =1
i =1
59

Una consideracin antes de encontrar las CdPO
El balance global del sistema se representa por la siguiente ecuacin:
n
Pg Pd
i
i =1
j =1
i =1
j =1
PL = 0
PL = Pg i Pd j
Observaciones:
1.
Recuerde que los Pdj son constantes. Por tanto, PL es funcin de las Pgi.
2.
Del flujo de potencia es sabido que Pg1 (la potencia neta inyectada) en el nodo slack (o nodo
de referencia) es funcin de las Pgi.
3.
Como consecuencia de las observaciones 1 y 2, PL es funcin de todas las potencias excepto la
potencia en el nodo slack, esto es PL PL ( Pg 2 , Pg3 , , Pg n )

60

m
n
Funcin de LaGrange L( Pg1 , , Pg n , ) = Ci ( Pgi ) Pgi Pd j PL

i =1
j =1
i =1
CdPO
1. Derivadas con respecto a las potencias generadas:
L
=0
Pgi
i = 1 n
m
n
L
n
=
C ( Pgi ) Pgi Pd j PL = CI1 ( Pg1 ) = 0
Pg1 Pg1 i =1 i
j =1
i =1

m
n
n
PL
L
=
Ci ( Pgi ) Pgi Pd j PL = CI i ( Pgi ) 1
= 0 i = 2 n
Pgi Pgi i =1
Pg
j =1
i
i =1

Recuerde que las prdidas esta

en funcin de las potencias
generadas, excepto Pg1
61

m
n
Funcin de LaGrange L( Pg1 , , Pg n , ) = Ci ( Pgi ) Pgi Pd j PL

i =1
j =1
i =1
CdPO
2. Derivada con respecto al multiplicador de LaGrange:
L
=0
Pg Pd
i =1
j =1
PL = 0
62

Las CdPO son, en resumen, las siguientes:
CI1 ( Pg1 ) = 0
PL
CI i ( Pgi ) 1
Pgi
n
Pg Pd
i =1
j =1
= 0 i = 2 n
PL = 0
Las ecuaciones del costo incremental se puede reescribir de esta manera (excepto la
asociada con el nodo slack):
1
PL
1
Pg i
CI i ( Pgi ) = i = 2 n
63

El trmino que multiplica a los costos incrementales en la diapositiva anterior se denomina el factor
de penalidad de prdidas del generador i y es denotado como Li, matemticamente:
1
i = 2 n
PL
1
Pg i
L1 = 1
Li =
Estos factores se pueden deducir

analticamente. Ver capitulo 11
del libro Power System Analysis
Arthur Bergen & V. Vittal
Las condiciones necesarias de optimalidad, teniendo en cuenta las perdidas, son las siguientes:
L1CI1 ( Pg1 ) = L2CI 2 ( Pg 2 ) = LnCI n ( Pg n ) =

n
Pg Pd
i =1
j =1
PL = 0
64

Comparacin entre el despacho econmico ignorando perdidas y considerando perdidas:
CdPO (ignorando prdidas)
CI1 ( Pg1 ) = CI 2 ( Pg 2 ) = CI n ( Pg n ) =
n
Pg Pd
i =1
j =1
=0
CdPO (considerando prdidas)
L1CI1 ( Pg1 ) = L2CI 2 ( Pg 2 ) = LnCI n ( Pg n ) =

n
Pg Pd
i =1
j =1
PL = 0
65

Ejemplo 1
(Tomado de Power System Analysis Arthur Bergen & V. Vittal)
Encuentre el despacho econmico para el siguiente sistema
Pd1=300 MW
Pd1=50 MW
CI1 ( Pg1 ) = 0.007 Pg1 + 4.1
CI 2 ( Pg 2 ) = 0.007 Pg 2 + 4.1
PL = 0.001( Pg 2 50 )
66

Ejemplo 1
Las CdPO para el despacho econmico considerando perdidas son:
L1CI1 ( Pg1 ) = L2CI 2 ( Pg 2 ) =

n
Pg Pd
i
i =1
j =1
PL = 0
Los factores de penalidad estn dados por:
L1 = 1
L2 =
1
1
=
PL 1.1 0.002 Pg 2
1
Pg 2
En el ejemplo 2 se muestra como

calcular esta expresion para un
sistema de dos nodos
67

Ejemplo 1
Las CdPO, por consiguiente, son:
L1CI1 ( Pg1 ) = 0.007 Pg1 + 4.1 =

L2CI 2 ( Pg 2 ) =
1
( 0.007 Pg 2 + 4.1) =
1.1 0.002 Pg 2
De forma equivalente
Pg1 =
4.1
0.007
1.1 4.1
Pg 2 =
0.007 + 0.002
68

Ejemplo 1
El sistema de ecuaciones es no lineal y requiere un mtodo iterativo.

Sistema de ecuaciones
Pg1 =
4.1
0.007
1.1 4.1
Pg 2 =
0.007 + 0.002
0.001( Pg 2 50 )
Pg1 + Pg 2 PL + 300 + 50 = 0
Sistema iterativo
Pg1
Pg1
PL
Balance
global
5.0
128.6
82.35
0.0786
-139.15
6.0
271.4
131.57
0.2214
52.78
5.5
200.0
108.33
0.1500
-41.81
5.75
231.7
120.27
0.1857
5.79
5.625
217.8
114.38
0.1679
-17.927
5.6875
226.78
117.34
0.1768
-6.0441
69

Ejemplo 2
Encuentre el despacho econmico para el siguiente sistema. Observe que es necesario

derivar una expresin para L2 pues no esta dada en trminos explcitos de Pg2.
V1 = 1e j 0
P1
V2 = 1e j2
P2
zl = 0.02 + j 0.08 pu
CI1 ( Pg1 ) = 0.7 Pg1 p.u + 4.1

CI 2 ( Pg 2 ) = 0.7 Pg 2 p.u + 4.1
PL ( Pg 2 ) = ?
Pd1=1 pu
70

Ejemplo 2
Para obtener PL en funcin de PG2 es necesario recurrir a las ecuaciones del flujo de
potencia, concretamente las ecuaciones de balance nodal. Recuerde que PL=P1+P2 y
estas potencias (P1+P2) se encuentran a partir de los balances nodales. Por lo tanto,
se necesita, como primer paso, el calculo de la Ybus.
v1 = 1e
v2 = 1e j2
j0
P1
P2
zl = 0.02 + j 0.08 pu
Ybus = Gbus + jBbus =

2.94 2.94
2.94 2.94 +
11.76 11.76
j
11.76 11.76
Pd1=3.0pu
71

Ejemplo 2
Las ecuaciones de balance nodal (ver notas de anlisis) son:
P/ i = Pgi Pdi = Vi V j ( Gij cos ij + Bij senij )

n
j =1
Pg1 Pd1 = P1 = V1V1 ( 2.94 cos 11 + 11.76 sen11 ) + V1V2 ( 2.94 cos 12 + 11.76 sen12 )
= 2.94 (1 cos 12 ) + 11.76sen12
Pg 2 Pd 2 = P2 = V2V1 ( 2.94 cos 21 + 11.76 sen 21 ) + V2V2 ( 2.94 cos 22 11.76sen 22 )

= 2.94 (1 cos 12 ) 11.76sen12
PL = P1 + P2 = 5.88 (1 cos 12 )
72

Ejemplo 2
El factor de penalidad de perdidas L2 se puede encontrar a travs de la siguiente

regla de la cadena:
dPL
dPL
dPL dPg 2
dPL
d12
=
=
d12 dPg 2 d12
dPg 2 dPg 2
d12
Derivando las expresiones de la diapositiva anterior
dPL
= 5.88sin 12
d12
dPg 2 dP2
=
= 2.94sin 12 11.76 cos 12
d12 d12
73

Ejemplo 2
El factor de penalidad de perdidas L2 se puede encontrar a travs de la siguiente

regla de la cadena:
L2 =
L2 =
1
1
=
5.88sen12
dPL
1
1
2.94sen12 11.76 cos 12
dPg 2
11.76 cos 12 2.94sen12 1 0.25 tan 12
=
11.76 cos 12 + 2.94s en12 1 + 0.25 tan 12
74

Ejemplo 2
Como encontrar el despacho econmico para este ejemplo?

1. Asumir un ngulo inicial 2
2. Calcular P1 y P2 basados en las ecuaciones de la diapositiva 72.
3. Calcular Pg1 y Pg2 con base en el punto 2 y las ecuaciones de la diapositivas 72.
4. Calcular el CI1 y CI2 con base en el punto 3 y los datos del problema
5. Calcular L2 con base en la suposicin 1. Recuerdo que L1=1
6. Se cumple la condicin de optimalidad?, esto es L1CI1=L2CI2
1. Si: Parar, la solucin es Pg1 y Pg2 encontradas en el punto 3.
2. No: Asumir un ngulo 2 diferente y repetir el procedimiento.
75

Finalizar el ejemplo anterior asumiendo un ngulo inicial 2=6.0o. Las ecuaciones que debe tener
presente, derivadas en el ejemplo anterior, son las siguientes:
Pg1 3.0 = P1 = 2.94 (1 cos 12 ) + 11.76 sen12
Pg 2 = P2 = 2.94 (1 cos 12 ) 11.76sen12

CI1 ( Pg1 ) = 0.7 Pg1 p.u + 4.1
CI 2 ( Pg 2 ) = 0.7 Pg 2 p.u + 4.1
L2 =
1 0.25 tan 12
1 + 0.25 tan 12
76
Taller
Taller:
Considere tres unidades de generacin cuyo costos de generacin estn
dados por la siguiente curvas de costo:
C1 ( Pg1 ) = 0.009 Pg12 + 5.00 Pg1 + 180 100 Pg1 300

C2 ( Pg 2 ) = 0.015Pg 2 2 + 5.50 Pg 2 + 714 100 Pg 2 300
C2 ( Pg3 ) = 0.018Pg 32 + 6.00 Pg 3 + 380 100 Pg3 300
77
Taller
Determine el despacho econmico, para una demanda de 700MW, en los siguientes
casos:
1.Ignorando limites de generacin.
2.Considerando los limites de generacin.
3.Considere que las unidades tiene la disposicin mostrada en la figura de
abajo. Encuentre el despacho econmico para ese sistema considerando las
perdidas del sistema e ignorando los limites de generacin.
V1 = 1.05e j 0
V2 = 0.98e j2
G2
G1
G3
zl = 0.01 + j 0.07 pu
Sb=100MVA
Vb=115KV
dT=700MW
78
Comentarios finales acerca del despacho

econmico
1.
El despacho econmico es una herramienta de anlisis bsica para el estudio

y anlisis de la operacin de sistemas de potencia.
2.
La solucin del despacho econmico considerando los efectos de red

requieren el uso de herramientas de optimizacin no cubiertas en estas
diapositivas. Sin embargo, las conclusiones que se obtiene de la solucin de
este problema con las herramientas estudiadas se extienden a problemas
que incorporan la red.
3.
La interpretacin de los operadores de LaGrange son de bastante utilidad en

el estudio y anlisis de los mercados elctricos. En futuras presentaciones se
hara uso de este importante concepto econmico.
79

Despacho Economico

Cargado por

Copyright:

Formatos disponibles

Despacho Economico

Cargado por

Información del documento

Descripción 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

Despacho Economico

Cargado por

Copyright:

Formatos disponibles

Despacho Econmico

Harold Salazar Isaza, Ph.D

Harold Salazar Isaza

Introduccin a el despacho econmico

Harold Salazar Isaza

Solucin al problema de despacho econmico

El costo de generacin de electricidad de las

estos componentes para formular

Las restricciones operativas

matemtico de cada uno de

Harold Salazar Isaza

Modelaje de la curva de costo

Harold Salazar Isaza

Costo del combustible (energas primarias) necesarias para

Costo de mano obra

Estos costos son

Modelaje de la curva de costo

Modelaje de la curva de costo

Los costos de las diapositivas anteriores reflejan los costos del

Harold Salazar Isaza

Puesto que el costo de produccin esta definido principalmente por el costo

Modelaje de la curva de costo

Entrada: Se puede expresar en

Harold Salazar Isaza

Salida: Potencia neta en MW

Modelaje de la curva de costo

Modelaje de la curva de costo

Observe que a medida que se incrementa la

Harold Salazar Isaza

Modelaje de la curva de costo

Como modelar el costo de generacin?

Normalmente este es el valor nominal

Modelaje de la curva de costo

Nivel mas eficiente

Normalmente este es el valor nominal

Modelaje de la curva de costo

( MBTU / h )= ( MW ) x ( MBTU / MWh )

Defina C como el costo de produccin establecido como el costo del combustible K

Modelaje de la curva de costo

El costo incrementa con la generacin. Esta

La curva, por lo general y para efectos de

Modelaje de la curva de costo

Harold Salazar Isaza

Dos observaciones de la curva de costo

La curva de costo incremental es

La curva de costo incremental es

Modelaje de la curva de costo

Harold Salazar Isaza

Modelaje de la curva de costo

Harold Salazar Isaza

Modelaje de la curva de costo

Harold Salazar Isaza

Modelaje de la curva de costo

Harold Salazar Isaza

Modelaje de la curva de costo

Harold Salazar Isaza

Modelaje de la curva de costo

Modelaje de la curva de costo

Harold Salazar Isaza