DETERMINACIN DE LA POTENCIA FIRME
REMUNERABLE A TRAVS DEL FLUJO DE
POTENCIA PTIMO
Optimizacin de Sistemas Elctricos de Potencia
Despeje del Mercado de Potencia
Formulacin matemtica del problema
El problema del despeje simultneo del mercado de potencia con las restricciones de la
red de transmisin, se puede formular como un problema de programacin matemtica
no lineal, equivalente al flujo de potencia ptimo.
El problema general del flujo de potencia ptimo corresponde a la determinacin de un
punto de operacin del sistema de potencia que optimiza una funcin objetivo y satisface
un conjunto de restricciones fsicas y de operacin. Matemticamente el problema puede
ser formulado como:
()
Sujeto a
() = 0
() 0
Dnde:
1()
() = [ ]
()
1()
() [ ]
()
Las restricciones de igualdad y desigualdad corresponden al modelaje esttico de la red y
a las limitaciones en la operacin de los equipos utilizados, respectivamente.
Funcin Objetivo
En este caso, la funcin objetivo corresponde a minimizar el precio total de la potencia, y
se puede expresar como:
( )
=1
Optimizacin de Sistemas Elctricos de Potencia
Restricciones de Igualdad
Las restricciones de igualdad corresponden a las ecuaciones de balance de potencia
activa en cada barra:
= 1,
Dnde:
= ( )2 ( cos + sin )
Y las ecuaciones de balance de potencia reactiva en cada barra:
+ 2 =
= 1,
Dnde:
= ( )2 ( + ) ( sin + cos )
, :
: .
1.
: "shunt" del circuito km
: "shunt"
Optimizacin de Sistemas Elctricos de Potencia
Restricciones de desigualdad
Las restricciones de desigualdad corresponden a las ecuaciones de los lmites fsicos.
As, se tiene los siguientes:
(, )
(, )
= 1,
= 1,
= 1,
: tap
, :
, :
, :
, :
Las variables en este problema son: , , .
De la formulacin se puede observar que se trata de un problema de optimizacin de
gran tamao con funcin objetivo no lineal y restricciones de igualdad y desigualdad no
lineales; lo cual exige la aplicacin de metodologas de optimizacin no lineal para su
solucin.
Optimizacin de Sistemas Elctricos de Potencia
El precio de la Potencia
La determinacin del precio de la potencia requiere encontrar la solucin ptima del
problema descrito en la seccin anterior.
Se sabe por teora de dualidad que asociado a cada restriccin no lineal existe un
Multiplicador de agrange. En el caso de las restricciones de balance de potencia activa,
existe un multiplicador en cada barra k; el cual mide la variacin del precio total de la
potencia (funcin objetivo) con respecto a la variacin de la demanda (lado derecho de la
restriccin). Este multiplicador corresponde al valor de la potencia en cada barra de la res
elctrica.
= 1,
La determinacin de los multiplicadores asociados a cada barra, se obtiene directamente
de la solucin del problema de optimizacin y no requiere ningn esfuerzo computacional
adicional.
Optimizacin de Sistemas Elctricos de Potencia
Ejemplo de Aplicacin
En esta seccin se muestra la formulacin matemtica para un sistema elctrico de
prueba. Una vez formulado el problema, el despeje del mercado de potencia se
efectuar con la asistencia de un solver de uso comercial.
Sistema de Prueba de 6 barras
La figura siguiente muestra la configuracin del sistema elctrico de 6 barras.
Figura 1.1. Sistema ejemplo de 6 Barras
Datos de Barras:
La tabla siguiente muestra los datos de las barras:
Nmero Carga Carga Potencia Firme Capacidad Max Capacidad Min Precio Limite Inf Limite Sup
Barra Activa Reactiva Activa Generacin Reactiva Generacin Reactiva Potencia Tensin Tensin
(MW) (MVAR) (MW) (MVAR) (MVAR) (US$/kW) (pu) (pu)
1 0 0 50 100 -100 70 0.90 1.10
2 0 0 75 100 -100 50 0.90 1.10
3 50 30 0 0 0 0 0.90 1.10
4 0 0 100 100 -100 10 0.90 1.10
5 60 36 0 0 0 0 0.90 1.10
6 80 48 0 0 0 0 0.90 1.10
Datos de Lneas:
La tabla siguiente muestra los datos de las lneas:
Nmero Barra Barra Resistencia Reactancia Conductancia Susceptancia Lmte Superior Limite Inferior
Lnea Salida Llegada Lnea Lnea Lnea Lnea Potencia Aparente Potencia Aparente
(Ohm) (Ohm) (Mho) (MVA) (MVA)
1 1 2 0.05 0.1 4.00 -8.00 1 -1
2 1 6 0.05 0.1 4.00 -8.00 1 -1
3 2 3 0.05 0.1 4.00 -8.00 1 -1
4 2 6 0.05 0.1 4.00 -8.00 1 -1
5 3 4 0.1 0.2 2.00 -4.00 2 -2
6 5 6 0.05 0.1 4.00 -8.00 1 -1
Optimizacin de Sistemas Elctricos de Potencia
Formulacin Matemtica:
Funcin Objetivo:
La funcin objetivo est dada por:
Min 700001 + 500002 + 100004
Restricciones del Flujo de Potencia Activa y Reactiva:
Las ecuaciones del flujo de potencia activa en los circuitos son:
= ( )2 ( cos + sin )
Circuito 1-2:
12 = 4(1 )2 1 2 (4 cos 12 8 sin 12 )
Circuito 1-6:
16 = 4(1 )2 1 6 (4 cos 16 8 sin 16 )
Circuito 2-3:
23 = 4(2 )2 2 3 (4 cos 23 8 sin 23 )
Circuito 2-6:
26 = 4(2 )2 2 6 (4 cos 26 8 sin 26 )
Circuito 3-4:
34 = 2(3 )2 3 4 (2 cos 34 4 sin 34 )
Circuito 5-6:
56 = 4(5 )2 5 6 (4 cos 56 8 sin 56 )
Las ecuaciones del flujo de potencia reactiva en los circuitos son:
= ( )2 ( + ) ( sin + cos )
Circuito 1-2:
12 = 8(1 )2 + 1 2 (8 cos 12 4 sin 12 )
Optimizacin de Sistemas Elctricos de Potencia
Circuito 1-6:
16 = 8(1 )2 + 1 6 (8 cos 16 4 sin 16 )
Circuito 2-3:
23 = 8(2 )2 + 2 3 (8 cos 23 4 sin 23 )
Circuito 2-6:
26 = 8(2 )2 + 2 6 (8 cos 26 4 sin 26 )
Circuito 3-4:
34 = 4(3 )2 + 3 4 (4 cos 34 2 sin 34 )
Circuito 5-6:
56 = 8(5 )2 + 5 6 (8 cos 56 4 sin 56 )
Restricciones de Balance de Potencia
Barra 1:
1 8(1 )2 + 1 2 (4 cos 12 8 sin 12 ) + 1 6 (4 cos 16 8 sin 16 ) = 0
1 16(1 )2 1 2 (8 cos 12 4 sin 12 ) 1 6 (8 cos 16 4 sin 16 ) = 0
Barra 2:
2 12(1 )2 + 2 1 (4 cos 21 8 sin 21 ) + 2 3 (4 cos 23 8 sin 23 )
+ 2 6 (4 cos 26 8 sin 26 ) = 0
2 24(1 )2 1 2 (8 cos 21 4 sin 21 ) 2 3 (8 cos 23 4 sin 23 )
2 6 (8 cos 26 4 sin 26 ) = 0
Barra 3:
6(3 )2 + 3 2 (4 cos 32 8 sin 32 ) + 3 4 (2 cos 34 4 sin 34 ) 0.5 = 0
12(3 )2 3 2 (8 cos 32 4 sin 32 ) 3 4 (4 cos 34 2 sin 34 ) 0.3 = 0
Optimizacin de Sistemas Elctricos de Potencia
Barra 4:
4 2(4 )2 + 3 4 (2 cos 43 8 sin 43 ) = 0
4 4(4 )2 3 4 (4 cos 43 2 sin 43 ) = 0
Barra 5:
4(5 )2 + 5 6 (2 cos 56 8 sin 56 ) 0.6 = 0
8(5 )2 5 6 (8 cos 56 4 sin 56 ) 0.36 = 0
Barra 6:
12(6 )2 + 1 6 (4 cos 61 8 sin 61 ) + 2 6 (4 cos 62 8 sin 62 ) + 5 6 (4 cos 65 8 sin 65 ) 0.8 = 0
24(6 )2 5 6 (8 cos 65 4 sin 65 ) 2 6 (8 cos 62 4 sin 62 ) 1 6 (8 cos 61 4 sin 61 ) 0.48 = 0
Restricciones de lmites
Circuito 1-2:
(4(1 )2 1 2 (4 cos 12 8 sin 12 )) + (8(1 )2 + 1 2 (8 cos 12 4 sin 12 ))
(4(1 )2 1 2 (4 cos 12 8 sin 12 )) + (8(1 )2 + 1 2 (8 cos 12 4 sin 12 ))
Circuito 1-6:
(4(1 )2 1 6 (4 cos 16 8 sin 16 )) + (8(1 )2 + 1 6 (8 cos 16 4 sin 16 ))
(4(1 )2 1 6 (4 cos 16 8 sin 16 )) + (8(1 )2 + 1 6 (8 cos 16 4 sin 16 ))
Circuito 2-3:
(4(2 )2 2 3 (4 cos 23 8 sin 23 )) + (8(2 )2 + 2 3 (8 cos 23 4 sin 23 ))
(4(2 )2 2 3 (4 cos 23 8 sin 23 )) + (8(2 )2 + 2 3 (8 cos 23 4 sin 23 ))
Circuito 2-6:
(4(2 )2 2 6 (4 cos 26 8 sin 26 )) + (8(2 )2 + 2 6 (8 cos 26 4 sin 26 ))
(4(2 )2 2 6 (4 cos 26 8 sin 26 )) + (8(2 )2 + 2 6 (8 cos 26 4 sin 26 ))
Optimizacin de Sistemas Elctricos de Potencia
Circuito 3-4:
(2(3 )2 3 4 (2 cos 34 8 sin 34 )) + (4(3 )2 + 3 4 (4 cos 34 2 sin 34 ))
(2(3 )2 3 4 (2 cos 34 8 sin 34 )) + (4(3 )2 + 3 4 (4 cos 34 2 sin 34 ))
Circuito 5-6:
(4(5 )2 5 6 (4 cos 56 8 sin 56 )) + (8(5 )2 + 5 6 (8 cos 56 4 sin 56 ))
(4(5 )2 5 6 (4 cos 56 8 sin 56 )) + (8(5 )2 + 5 6 (8 cos 56 4 sin 56 ))
Restricciones del Lmite de la Tensin:
Barra 1:
V1 1.1
V1 0.9
Barra 2:
V2 1.1
V2 0.9
Barra 3:
V3 1.1
V3 0.9
Barra 4:
V4 1.1
V4 0.9
Barra 5:
V5 1.1
V5 0.9
Barra 6:
V6 1.1
V6 0.9
Optimizacin de Sistemas Elctricos de Potencia
Restricciones del Lmite de la generacin Activa y Reactiva
Generador 1:
1 0.5
1 0
1 1
1 1
Generador 2:
2 50.
2 0.75
2 1
2 1
Generador 4:
4 1
4 0
4 1
4 1
Restriccin del Angulo
Circuito 1-2:
12 0.15
12 0.15
Circuito 1-6:
16 0.15
16 0.15
Optimizacin de Sistemas Elctricos de Potencia
Circuito 2-3:
23 0.15
23 0.15
Circuito 2-6:
26 0.15
26 0.15
Circuito 3-4:
34 0.15
34 0.15
Circuito 5-6:
56 0.15
56 0.15
Optimizacin de Sistemas Elctricos de Potencia
Resultados obtenidos
La solucin del problema formulado en la seccin anterior se efectu con el solver del
Excel, la solucin alcanzada se resumen a continuacin:
Resultados de Barras
Nm ero Nom bre Generacin Generacin Carga Carga Capacidad
Barra Barra Activa Reactiva Activa Reactiva Activa
(MW) (MVAR) (MW) (MVAR) (MW)
1 Barra 1 0.4738 0.9239 0 0 0.5
2 Barra 2 0.7500 0.5513 0 0 0.75
3 Barra 3 0.0000 0.0000 0.5 0.3 0
4 Barra 4 0.8618 0.0361 0 0 1
5 Barra 5 0.0000 0.0000 0.6 0.36 0
6 Barra 6 0.0000 0.0000 0.8 0.48 0
2.0856 1.5112 1.90 1.14
Capacidad MaxCapacidad Min Costo Tensin Lim ite Inf Lim ite Sup Angulo
Reactiva Reactiva Genreacin Barra Tensin Tensin Barra
(MVAR) (MVAR) (US$/MW) (pu) (pu) (pu) (Grados)
1 -1 70 1.067 0.90 1.10 -32.16
1 -1 50 1.051 0.90 1.10 -30.47
0 0 0 1.027 0.90 1.10 -27.85
1 -1 10 1.100 0.90 1.10 -19.25
0 0 0 0.900 0.90 1.10 -36.79
0 0 0 0.974 0.90 1.10 -34.04
Potencia Potencia Costo Angulo Costos Factor
Activa Reactiva Operacin Barra Marginales Penalizacin
(MW) (MVAR) (US$) (Radianes) (US$/MW)
0.00 0.00 33.17 -0.5613 70.00 1.00
0.00 0.00 37.50 -0.5318 68.20 0.97
0.00 0.00 -0.4860 64.89 0.93
0.00 0.00 8.62 -0.3360 52.10 0.74
0.00 0.00 -0.6421 80.17 1.15
0.00 0.00 -0.5942 74.36 1.06
0.00 0.00 79.283
Optimizacin de Sistemas Elctricos de Potencia
Resultados de Lneas
Nm ero Barra Barra Resistencia Reactancia Conductancia Susceptancia
Lnea Salida Llegada (Ohm ) (Ohm )
1 1 2 0.05 0.1 4.00 -8.00
2 1 6 0.05 0.1 4.00 -8.00
3 2 3 0.05 0.1 4.00 -8.00
4 2 6 0.05 0.1 4.00 -8.00
5 3 4 0.1 0.2 2.00 -4.00
6 5 6 0.05 0.1 4.00 -8.00
Flujo Flujo Flujo km Lim te Sup Lim te Inf Pridas Flujo m k
Activa Reactiva Aparente Aparente Aparente Potencia Aparente
(MW) (MVAR) (MVA) (MVA) (MVA) (MW) (MVA)
-0.1948 0.2704 0.33 1 -1 0.0049 0.33
0.6686 0.6535 0.93 1 -1 0.0384 0.85
-0.2890 0.4096 0.50 1 -1 0.0114 0.49
0.8392 0.4023 0.93 1 -1 0.0392 0.86
-0.8003 0.0869 0.81 2 -2 0.0615 0.86
-0.6000 -0.3600 0.70 1 -1 0.0302 0.76
Barra Barra Variacin
Angular Absoluta
Salida Llegada (Rad)
1 2 0.03
1 6 0.03
2 3 0.05
2 6 0.06
3 4 0.15
5 6 0.05
Punto de operacin
El punto de operacin optimo encontrado se muestra en la figura 1.2. Es importante
destacar que la potencia activa inyectada por cada generador corresponde a su potencia
firme remunerable.
Precio de Potencia
El precio de la potencia obtenido en cada barra del sistema se muestra en la figura 1.3.
Optimizacin de Sistemas Elctricos de Potencia
47.38 75.00 50.00
92.39 55.13 30.00
-19.48 0.00
27.04 0.00
1 2
-28.90 0.00
3
66.86 40.96 0.00
65.35 83.92 -80.03
40.23 8.69
0.00 0.00
0.00 0.00
0.00
0.00
6 0.00 -60.00 5 4
0.00 -36.00
80.00 60.00 86.18
48.00 36.00 3.61
# Potencia Activa # Flujos de Potencia Activa
# Potencia Reactiva # Flujos de Potencia Reactiva
Figura 1.2. Punto de Operacin Optimo
68.20
70.00 64.89
1 2 3
74.36
6 5 80.17
4 52.10
Direccin del Flujo de Potencia Activa
# Costos Marginales
Figura 1.3. Precio de la Potencia en cada Barra
Optimizacin de Sistemas Elctricos de Potencia
