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