FACULTAD DE ING.

ELCTRICA Y ELECTRNICA

UNIVERSIDAD NACIONAL DEL

CENTRO DEL PER

FACULTAD DE INGENIERA ELCTRICA Y

ELECTRNICA
TEMA:
FLUJO PTIMO DE POTENCIA
CTEDRA:
ANLISIS DE SISTEMAS DE POTENCIA II
CATEDRTICO:
Mg. TORRES MAYTA, Pedro
ESTUDIANTE:
ALANYA PARIONA, Jaime
CABRERA PALACIOS, Luis Alejandro
GALA FLORES, Lidia Cecilia
MATOS CAPCHA, Erik
MELGAR ASTETE, Milagros Nashely
ANLISIS DE SISTEMAS DE POTENCIA II

FACULTAD DE ING. ELCTRICA Y ELECTRNICA

37

HUANCAYO 2015

ANLISIS DE SISTEMAS DE POTENCIA II

DEDICATORIA:
A nuestros padres; Por su apoyo
incondicional que me brindan y por
estar siempre con nosotros.
A nuestros profesores: Quienes
son nuestros guas en el aprendizaje,
dndonos los ltimos conocimientos

para nuestro buen desenvolvimiento

en la sociedad.

Los Alumnos

AGRADECIMIENTO
Agradecemos en primer lugar, al ser Supremo, nico dueo de
todo saber y verdad, por iluminarnos durante este trabajo y
por permitirnos finalizarlo con xito; y en segundo lugar, pero
no menos importante, a nuestros queridos padres, por su
apoyo incondicional y el esfuerzo diario que realizan por
brindarnos una buena educacin.

NDICE GENERAL
1.1
1.2
1.3
1.4

I. Definiciones......... 1
Flujo ptimo de potencia..... 1
Despacho econmico.. 1
Planificacin de SEP....... 2
Operacin de SEP... 2
II. Importancia de flujo ptimo en la planificacin y operacin de SEP.. 2
III. Formulacin General....
3
a)
Variables...... 4
b)
Restriccin de Igualdad...... 5
3.1
Restriccin de Desigualdad..... 6
3.2
tiempos mnimos de encendido y apagado.......
7
3.3
Lmites de velocidad de cambio de generacin....
8
3.4
Lmites de flujos en las ramas..... 8
3.5
Perfiles de voltaje.....
9
3.6
Lmites de la posicin de TAPS...... 9
IV. Formulacin Matemtica.... 10
V. Mtodo de la Gradiente reducido para la Solucin del Flujo ptimo de
Potencia....... 14
VI. El problema de Flujos de Potencia ptimo con Restricciones de
Seguridad.... 17
VII. Restricciones y problemas de la Optimizacin.....
30
7.1
Restricciones..........
30
7.2
Restriccin de Desigualdad....... 31
7.3
Restriccin de Seguridad........
32
7.4
El problema de flujos de Potencia ptimo......
33
VIII. Descripcin General de una Red Elctrica..... 34

8.1
34
8.1.1
35
8.1.2
35
8.1.1
35
8.2
8.3
36

Generacin de energa elctrica.........

Generacin

hidroelctrica..........

Generacin trmica convencional..........

Generacin trmica de ciclo combinado............
Generador........... 36
Transformadores.............

Conclusiones............... 38
Bibliografa..... 39

INTRODUCCION
El presente trabajo enfoca el flujo ptimo de potencia en general con las que son variables
de estado, de la formulacin de las funciones objetivo, y la tcnica del gradiente reducido
a optimizacin.
El punto de operacin de un SEP se obtiene cuando las variables de control toman
valores reales que permiten la funcin objetivo planteado, tenga su valor mximo o ms
cercano a este.
Se presentan luego ejercicios desarrollados sobre las frmulas matemticas desarrolladas
y con el mtodo de la gradiente.

Los alumnos.

FLUJO PTIMO DE POTENCIA

1. DEFINICIONES:
1.1.

Flujo ptimo de potencia:

El estudio que permite la formulacin del flujo de potencia, optimizndolo en algn

sentido y cumplimiento al mismo tiempo un conjunto de restricciones, recibe el
nombre de flujo ptimo de potencia. En la formulacin del flujo ptimo se requiere
de alguna funcin objetivo, la cual se debe optimizar (maximizar o minimizar), y
de una tcnica de optimizacin. Alguna de las funciones objetivo son:

Minimizar los costos de generacin.

minimizar las perdidas des sistema.

minimizar la compensacin de potencia reactiva.

El proceso de optimizacin requiere la divisin de las variables del sistema, en

variables de control y variables de estado. La solucin ptima se obtiene
encontrando el valor de las variables de control que minimicen la funcin objetivo,
y al mismo tiempo satisfagan las estricciones del problema.

1.2. Despacho econmico: El proceso de repartir la potencia activa total de carga

entre las diferentes unidades de generacin, de un sistema de potencia, con el
objeto de realizar la mayor

economa en el funcionamiento del mismo, recibe

el nombre de despacho ptimo o despacho econmico. Dicho en otras palabras,

el despacho econmico es un flujo ptimo de potencia en el que la funcin
objetivo es minimizar los costos de generacin. Para realizar

despacho

econmico, se considera que las unidades de generacin son trmicas. Si las

unidades de generacin son hidrulicas, se considera que puede obtenerse
el equivalente en unidades trmicas o que previamente han sido optimizadas.

1.3. Planificacin

de SEP: La planificacin

de un sistema de potencia es el

estudio que se hace, de una posible expansin del sistema existente, debida
principalmente

al incremento de carga por la aparicin de nuevos centros de

consumo. Desde este punto de vista, la planificacin es un estudio a mediano o

largo plazo. El estudio incluye la proyeccin de la demanda, lo cual exige de
tcnicas de proyeccin determinsticas, probabilsticas o estocsticas,
acuerdo con la cantidad, tipo e incertidumbre

de

de la informacin de que se

disponga.

1.4. Operacin de SEP: La operacin de sistemas de potencia es el estudio a corto

plazo, y en tiempo real, en el cual las fuentes de potencia del sistema existente,
deben ser programadas para satisfacer las condiciones de carga, de modo que el
sistema en conjunto opere en las mejores condiciones, tanto tcnicas como
econmicas.

2. IMPORTANCIA

DEL

PLANIFICACIN Y

FLUJO

PTIMO

EN

LA

OPERACIN DE SEP

Como hemos visto en el apartado anterior, siempre se ha tratado de con seguir un

modelo del sistema de potencia mediante el cual se optimice el funcionamiento del
mismo. Desde este punto de vista, el planificar un sistema de potencia tambin va

encaminado a conseguir las condiciones ptimas de operacin del sistema futuro.

Un estudio de

flujos de potencia ayudar

considerablemente

a conocer el

funcionamiento del sistema proyectado, y si este forma parte de un estudio de flujos

ptimos, se puede
transmisin

llegar

a conocer

distribucin,

las capacidades

que lgicamente

estarn

ptimas de
ligadas

generacin,
a una cierta

incertidumbre, ya que la carga proyectada tiene tambin una incertidumbre.

Dentro de la operacin de sistemas de potencia, el poder realizar

de pacho

econmico, minimizacin de prdidas, o ambos a la vez, convierte al flujo ptimo de

potencia en un modelo poderoso y til para

dicha actividad. Adems de poder

programar adecuadamente las centrales de generacin, se podr mantener

un

buen perfil de voltaje a lo largo de toda la red, mediante las restricciones que se
impongan al flujo ptimo,

lo cual permite que el sistema opere en las mejores

condiciones tcnicas y econmicas para una carga dada.

En general, podemos concluir que la operacin, actual o futura, de un sistema

elctrico de potencia, puede ser sintonizada en su punto ptimo, mediante el flujo
ptimo de potencia, de aqu que ste es una herramienta necesaria para la
planificacin y operacin de sistemas elctricos de potencia.

3. FORMULACION GENERAL
Se han considerado varias funciones objetivo en un FPO, pero la que ms
frecuentemente se usa toma en cuenta los costos de generacin, la que refleja
aspectos econmicos del sistema de potencia. De aqu que la formulacin matemtica
del FPO se enfoca en minimizar el costo de generacin de potencia activa por un

ajuste adecuado de las variables de control. De forma general, el FPO puede ser
formulado como un problema de optimizacin no lineal con restricciones, que
matemticamente se expresa como:

Min f( x,u)
h(x,u) = 0
y g(x,u) <0

(2.1)

Donde, xRn1 es un vector de n1 variables de estado, uRn2 es un vector de n2

variables de control, f(x,u) es la funcin objetivo, h(x,u)=0, representa r restricciones
de igualdad, g(x,u) representa m restricciones de desigualdad, que consideran los
lmites operativos de las variables de control y las restricciones operativas.

a) Variables:
Las variables que pueden ser ajustadas con el fin de encontrar la solucin ptima y
satisfacer las restricciones se denominan variables de control. Las variables de control
pueden incluir:

La potencia activa asignada a las unidades de generacin.

La potencia reactiva asignada a las unidades de generacin.

Magnitud de voltaje de los nodos de generacin.

Posicin del tap de los transformadores.

ngulo de los transformadores cambiadores de fase.

Potencia reactiva de los condensadores sncronos.

Potencia reactiva de un compensador de reactivos.

Magnitud de interrupcin de carga.

Parmetros de control de los dispositivos FACTS.

Las variables de control pueden tomar valores continuos dentro de sus lmites
operativos (como es el caso de la generacin de potencia activa de los generadores),
o valores discretos (por ejemplo la posicin del tap de un transformador).

Las variables que se calculan a partir de las variables de control se denominan

variables de estado (o dependientes). Las principales variables de estado son:

La potencia en el nodo slack,

Magnitud de voltaje en los nodos de carga.

Potencia reactiva asignada a los generadores.

ngulo de fase de todos los nodos (excepto el slack).

Flujos de potencia en la red de transmisin.

El costo de generacin de potencia activa se emplea ampliamente como funcin

objetivo en un FPO.

b) Restricciones de igualdad:
Las restricciones de igualdad son tpicamente las ecuaciones de balance de carga, las
que se obtienen al imponer una restriccin de balance de potencia activa y reactiva en
todos los nodos del sistema. En un punto de operacin en estado estable, la potencia
generada debe ser tal que sea suficiente para cubrir la demanda ms las prdidas en
la red. Las ecuaciones de balance consideran un punto de equilibrio de potencia activa
y reactiva que debe satisfacer cada una de los nodos:

Donde, es el nmero de nodos del sistema de potencia excepto el nodo slack; son
todos los nodos del sistema de potencia; PGi y QGi son la potencia activa y reactiva
de generacin programada para el i-simo generador, respectivamente; PDi y QDi son
la potencia activa y reactiva, respectivamente, de la carga conectada en el i-simo
nodo. y son la potencia activa y reactiva, respectivamente, inyectada en el i-simo
nodo. La Figura 2.1 muestra un diagrama unifilar de un nodo que incluye generacin,
carga, y lneas de transmisin.

3.1. Restricciones de desigualdad

Las restricciones de desigualdad consideran los lmites que deben satisfacer las
variables de control y estado. Estas restricciones reflejan los lmites operativos
impuestos a los dispositivos y al sistema elctrico de potencia. Las principales
restricciones de desigualdad consideradas en un FPO son:

Lmites de potencia activa y reactiva de generacin.- La potencia activa y reactiva

asignada a las unidades de generacin deben estar dentro de lmites de
operacin. Esta restriccin refleja los lmites operativos y trmicos que un
generador debe satisfacer; se les modela por (2.3).

Donde, PGi_MIN y PGi_MAX son los lmites de potencia activa mnima y mxima,
respectivamente, para el i-simo generador; QGi_MIN y QGi_MAX son los lmites
de potencia reactiva mnima y mxima, respectivamente, para el i-simo
generador.
3.2. Tiempos mnimos de encendido y apagado:
Cuando un generador cambia de estado, por motivos operativos y trmicos debe
satisfacer restricciones de un tiempo mnimo que la unidad debe permanecer en su
nuevo estado. De aqu: donde:

ON Timei : Representa el tiempo que la i-sima unidad ha permanecido encendida.

OFF Timei:: Representa el tiempo que la i-sima unidad ha permanecido apagada.

ON TimeMIN,i: Es el tiempo mnimo que la i-sima unidad debe permanecer

encendida antes de poder ser apagada nuevamente.

OFF TimeMIN,i: Es el tiempo mnimo que la i-sima unidad debe permanecer

apagada antes de poder ser encendida nuevamente.

3.3. Lmites de velocidad de cambio de generacin:

Los generadores deben ser capaces de pasar a un nuevo valor de generacin en un
perodo corto de tiempo. Sin embargo, si los generadores son unidades grandes, por
limitantes trmicas, stas no sern capaces de modificar la generacin de potencia
activa por encima de un valor denominado velocidad de cambio de generacin. Esta
caracterstica de los generadores requiere la inclusin de restricciones que consideren
este aspecto. Por lo tanto, cada unidad trmica debe obedecer a un lmite en la
velocidad de cambio de generacin de tal manera que:

Dnde:

3.4. Lmites de flujos en las ramas:

Con la finalidad de mantener la seguridad en los sistemas de potencia, los enlaces
(lneas o transformadores) no deben ser sobrecargados. Por lo tanto, es necesario
definir lmites para todas las ramas (o para un grupo de ellas). Estos lmites pueden
deberse a restricciones trmicas de los equipos o por consideraciones de seguridad
del sistema. Los lmites de los flujos se pueden formular como:

Donde, FLOW (m,n)MAX representa el mximo flujo de potencia activa permitido en la

rama que conecta los nodos m y n.
3.5. Perfiles de voltaje:
Debido a que el voltaje en los nodos es uno de los criterios de seguridad e ndice de
calidad de servicio ms importante, incluir una restriccin que mejore el perfil de
voltaje de los nodos de carga del sistema es un aspecto importante a ser considerado
en el problema de optimizacin. Los voltajes en los nodos de generacin son
constantes, mientras que el nivel de voltaje en los nodos de carga debe mantenerse
muy cercano a un voltaje de referencia. Matemticamente esta restriccin se puede
definir mediante:

Donde, Vi, representa la magnitud de voltaje en el i-simo nodo de carga, Vref,

representa la magnitud de voltaje de referencia, en general definido como Vref = 1.0
p.u.
3.6. Lmites de la posicin de taps:
Los transformadores con cambiador de derivaciones son utilizados como dispositivos
para el control de la magnitud de voltaje (y como resultado de esto el control del flujo
de potencia reactiva). Estos transformadores estn restringidos por valores del
cambiador de tap dentro de lmites mnimos y mximos.

Donde Tapi, es la posicin en la que se encuentra el i-simo tap del transformador;

Tapi_MIN y Tapi_MAX, son la posicin mnima y mxima, respectivamente, permitida
para el i-simo transformador.
Otras restricciones que pueden ser incluidas en un FPO son::

4.

Lmites del cambiador de fase de un transformador.

Lmites de las variables de control de los dispositivos FACTS.

Lmites de compensacin de potencia reactiva.

Requerimientos de reserva rodante.

Lmites de emisin de contaminantes al medio.

FORMULACION MATEMATICA:
El sistema supone en rgimen sinusoidal permanente, debido a la lenta evolucin de
la carga, y equilibrado (modelo monofsico). Para una red de N nudos el vector de
estado x contiene n=2N-1 elementos (uno de los ngulos se toma como origen de
fases). Considrese un vector Z compuesta m medidas las cuales se relacionen con el
vector de estado mediante la expresin:

Las ecuaciones para el clculo de flujo de carga estn dadas por (1) y (2)

En este caso, los valores calculados de las inyecciones de potencia activa y reactiva
estn dados por (3) y (4).

Dnde:

Las expresiones para los valores netos de potencias activa y reactiva en las barras
estn dadas por (5) y (6)

Dnde:

Cada nodo es caracterizado por cuatro variables: potencias activa y reactiva netas,
magnitud de voltaje y ngulo. En cada nodo, dos de estas variables son especificadas

y dos son desconocidas. Dependiendo de cuales variables son especificadas, las

barras se pueden dividir en tres tipos:
a. Barra de referencia: en esta barra la magnitud de voltaje y ngulo son
especificados, y se desconocen las potencias activas y reactivas netas.
b. Barras PQ: en estas barras los valores de las potencias activa y reactiva netas
son especificados, y se desconocen los valores de la magnitud de tensin y ngulo.
Normalmente son conocidas como barras de carga.
c. Barras PV: en estas barras el valor de la potencia activa neta y la magnitud del
voltaje son especificados, y se desconocen la potencia reactiva y el ngulo.
Normalmente son conocidas como barras de generacin.

El problema de flujo de potencia consiste en calcular las magnitudes de tensin y los

ngulos para cada una de las barras. Una vez estos valores son conocidos, los
valores de potencias activa y reactiva desconocidas se pueden encontrar de forma
trivial mediante sustitucin en (3) y (4). Sea [x] el vector de las variables desconocidas
V y , el mtodo de Newton en forma polar para encontrar el valor del vector [x] est
dado por el algoritmo:
1. Seleccionar un nmero de ecuaciones de (1) y (2) igual al nmero de Incgnitas
para formar el vector [g(x)]. En este caso para las barras PQ se seleccionan ambas
ecuaciones dado que es posible calcular las Potencias activas y reactivas netas. Para
las barras PV solamente se selecciona la primera ecuacin, pues no se conoce la
potencia reactiva Neta.

2. Asignar un valor inicial para las magnitudes de tensin y los ngulos. Normalmente
se utiliza un inicio plano asignando un valor inicial de uno en por unidad para todas
las tensiones y cero para los ngulos.
3. Calcular los nuevos valores de potencias activas y reactivas inyectadas utilizando
(3) y (4), y compararlos con los valores netos en (1) y (2). Si estas ultimas ecuaciones
se cumplen dentro de determinada tolerancia especificada entonces pare, de lo
contrario contine al paso 4.
4. Encuentre un nuevo conjunto de voltajes y ngulos resolviendo (7) y (8). Regrese
al paso 3 utilizando los nuevos valores de voltajes y ngulos encontrados.

En (7) la derivada parcial de g(x) respecto a x es conocida como matriz Jacobiana y su

estructura est dada por:

Dnde:

5. MTODO

DEL

GRADIENTE

REDUCIDO

PARA

LA

SOLUCIN DEL FLUJO PTIMO DE POTENCIA

En general, cualquier punto de solucin factible, no necesariamente optimo, el flujo de
potencia puede ser resuelto. Una vez resuelto el flujo de potencia con a ayuda del
jacobiano de la ltima iteracin y con el vector de las derivadas f con respecto al
vector

XX , se puede encontrar el vector de los multiplicadores de lagrange

X
,

aplicando la ecuacin:

(1)

(2)

...(3)

Reemplazando este resultado en la ecuacion Xen la ecuacion 2 nos encontraremos en

general que

, ya que necesariamennte estaremos en el optimo . este vector

tiene un significado importante ; es el vector gradiente VXf u , el cual es ortogonal a

los contornos de los valores a los contornos de los valores constantes de la funcin de
objetivo.
Los contornos de los valores constantes de la funcin objetivo, en general pueden ser
curvas cncavas convexas. U na funcin es convexa si la interpolacin lineal entre los
valores de dos puntos cuales quiera, nos da resultado un valor no menor que el punto
que tiene igual abscisa que el punto de interpolacin, y que permanece en la curva
siguiente:

a)

En el caso de dimensiones tendramos:

Los contornos de igual valor de las funciones objetivos planteadas tienen , en el

caso de dos dimensiones , las siguientes formas :

En donde las curcas mas abiertas tienen mayor valor , asi ,f5>f4>>f1. Cada
una de las curvas indica el lugar geometrico en donde la funcion objetivo tiene
un valor constante asi por ejemplo se tendra un valor f4 de la funcion objetivo ,

cuando las variables de control tomen los valores (u11,u21), o cuando tomen los
valores (u12,u22).

Como vemos , para obtener un valor de la funcion objetivo, se puede obtener un

infinito numero de soluciones , de las cuales unas seran factibles , y otras no,
como se vera mas adelante cuando

se introduzcan

con la srestricciones de

desigualdad sobre los parametros de control.

Por otro , las ecuaciones (1),(2),(3) son ecualciones no lineales

y pueden

resolverse solo con metodos iteractivos .El esquema iterativo mas simple es el
metodo del descenso mas pronunciado, llamado tambien metodo de la
gradiente. La idea basica de este metodo es partir de una solucion factible
del flujo de potencia (un punto en la figura anterior ) y moverse alo largo de la
direccion del descenso mas pronunciado , osea en la direccion de la gradiente
negativo, para encontrar un nuevo punto factible de solucion, pero estara mas
cercano al punto optimo , osea que se tendra un menor valor en la funcion
objetivo , graficamente , en dos dimensiones , tendremos :

El grafico anterior, suponemos que A es el punto de solucin del flujo de potencia.

En este punto evaluamos el vector de la gradiente VXf u, y nos movemos al punto
B, a lo largo de la direccin de la gradiente negativa. En el punto B se tendr un
valor de la funcin objetivo que ser menor que el valor que tiene en el punto A,
por lo cual, estaremos cerca del ptimo de antes. Ahora en el punto B se deber
satisfacer la ecuacin (3), por lo que se deber resolver un flujo de potencia para
conocer el valor del valor de estado XX. Una vez satisfecha la ecuacin (3) se
tendr

que B es un nuevo punto de solucin

factible y se podr iniciar

nuevamente el proceso, llegando al punto C, luego al D, y as sucesivamente,

hasta satisfacer un criterio de convergencia y consecuentemente se encontrara el
valor mnimo de la funcin de objetivo.
Este es , a grandes rasgos , el mtodo del gradiente .el mtodo es llamado as
porque como se ver ms adelante no solo cambia el signo de la gradiente, sino
que ms adelante se le multiplica por la constante C , cuyo valor depender del
valor de la gradiente , as como el valor de la funcin objetivo , que se tenga en un
punto dado, o ms bien , grficamente , en el valor de C ser tal , en cada
iteracin , que permita un camino corto entre el punto inicial A inicial y el punto
ptimo.

6. EL PROBLEMA DE FLUJOS DE POTENCIA PTIMO CON

RESTRICCIONES DE SEGURIDAD
EJEMPLO 1.
En este caso se considera un sistema constituido por 5 nodos, 7 lneas y dos
generadores [1]. La susceptancia serie de las lneas de transmisin se muestran en el
diagrama unifilar (Figura. 4.1); la resistencia de las lneas no se considera.

Los lmites mximos y mnimos de generacin, la velocidad de cambio de generacin

de las unidades trmicas, y los coeficientes de las funciones de costo se presentan en
la Tabla 4.1.
Los coeficientes de la curva de costos de generacin dados en la Tabla 4.1,
corresponden a una funcin cuadrtica definida por:

La curva de beneficio del consumidor se presenta como una ecuacin en funcin del
consumo de potencia activa, y se define mediante:

Donde

representan los i-simos coeficientes de la funcin de beneficio del

consumidor los cuales se presentan en la Tabla 4.1.

El costo por interrupcin de carga se define por la curva:

Donde

representa el i-simo coeficiente de interrupcin de carga el cul se

presenta en la Tabla 4.1

Los lmites de capacidad de transporte de las lneas de transmisin (en MW) se

exhiben en la Figura 4.1. El valor menor corresponde al lmite para el estado pre
contingencia, mientras el valor mayor (valor entre parntesis), es la capacidad de

transmisin de las lneas para el caso post-contingencia.

Antes de aplicar FPO-RS al sistema en estudio, se ejecuta un simple FPO, que no
considera contingencias. El objetivo de este problema es minimizar el costo total de
operacin sin el costo asociado por la ocurrencia de contingencias. De las ecuaciones
(2.9) y (2.10), la funcin objetivo para este ejemplo se define por:

El pronstico de demanda para el estado pre-contingencia es el que se ilustra en la

Figura. 4.1.

La solucin del problema de FPO, Fig. 4.2, indica que la mxima generacin de la
unidad ms barata (G1) est limitada por la mxima transferencia de potencia
permitida a travs de las lneas que conectan los nodos 1-2, y 2-3, que operan cerca
de sus lmites operativos. De aqu que la salida de estas lneas se consideran como
contingencias de mayor impacto para el problema de FPO-RS.
Para el caso de FPO-RS, se determina la asignacin de potencia activa para los
generadores, y se estudian dos casos. En el primero no se considera como variables
de control la demanda de potencia activa en los nodos de carga, y como segundo
caso se incluyen las demandas como variables de control.

CASO 1. La funcin objetivo se define por:

La Fig. 4.3 ilustra la estructura de la poblacin y una partcula para este caso. La
poblacin est constituida por tres subpoblaciones, una para el caso precontingencia y
una por cada una de las contingencias consideradas (salida de las lneas que unen los
nodos 1-2 y 2-3). Las dos variables de control que se consideran, corresponden a la
potencia activa de generacin asignada a las unidades trmicas disponibles en el
sistema de prueba.

El despacho de generacin resultante se despliega en la Tabla 4.2. De los resultado

obtenidos, comparando el despacho de generacin para la solucin del problema de
FPO, Figura. 4.2, con aquellos obtenidos por el despacho de generacin para el
problema de FPO-RS, se observa que el generador 1 disminuye su generacin de 188
MW a 146 MW (columna 3, Tabla 4.2). Adems, indica que cuando la salida de la lnea
que une los nodos 1-2 se considera como contingencia, el flujo a travs de la lnea 1-5
presenta un sobreflujo de 1.45 MW (columna 4 de la Tabla 4.2). Aunque el generador
1 puede disminuir su generacin para eliminar este sobreflujo, el generador 2 no
puede incrementar su generacin sin violar la restriccin de la velocidad de cambio de
generacin. Por lo tanto, una interrupcin de carga es necesaria para mantener la
seguridad del sistema. El valor ptimo de interrupcin de carga ser obtenido ms
adelante.

CASO 2. Para el segundo caso se adicionan los valores de demanda como variables
de control, lo que permite determinar los valores de stas para satisfacer la restriccin
que no pudo ser satisfecha en el caso previo.

La Figura. 4.4 esquematiza la estructura de la poblacin y una partcula.

Para este caso, la ecuacin dada por (4.7) se utiliza para calcular el costo para el
estado pre-contingencia. Para el estado post-contingencia, el trmino que considera el
costo por interrupcin de carga (4.3) debe ser agregado a la ecuacin (4.8).
La Tabla 4.3, despliega el despacho de generacin, carga y los flujos a travs de las
lneas.

Note que adems del re-despacho de generacin, la carga que se encuentra

conectada en el nodo 5 debe disminuir su consumo en el caso de que ocurra la salida
de la lnea 1-2; esto permite mantener los flujos de las lneas dentro de sus lmites.
Para este caso de estudio, los mejores resultados reportados en [1] se muestran en la
Tabla 4.4.
De la Tabla 4.3 y 4.4, se observa que la metodologa propuesta es capaz de encontrar
una mejor solucin, reduciendo la funcin objetivo de $ -4556.00, como mejor
resultado reportado en [1], a $ -4569.80 aplicando el OEP-OR.

EJEMPLO 2.
Se considera un sistema compuesto por 6 unidades de generacin, 26 nodos y 46
lneas de transmisin, Fig. 4.5. La demanda total del sistema es de 1263 MW, los
parmetros de las lneas y los datos de las cargas se toman de [2].
Los coeficientes de la curva de costo se presentan en la Tabla 4.5. Estos
corresponden a una curva de costos de generacin de potencia activa con una
componente senoidal que representa el efecto de apertura de las vlvulas de las
calderas de un generador, (4.9), [3].

El sistema tiene un total de 22 variables de control, 6 correspondientes a la generacin

de potencia activa de las unidades, 7 correspondientes a los taps de los
transformadores, y 9 para la potencia reactiva aportada por los dispositivos
conectados en derivacin. Los transformadores tienen un total de 20 pasos de ajuste,
en el rango de [0.9 a 1.1]. En lo que respecta al aporte de reactivos de los capacitores,
stos tiene un total de 10 pasos de ajuste en el rango de [0.0, 5.0 MVAR].

La Fig. 4.5 incluye las capacidades de transmisin de potencia activa en MW para los
estados pre- y post-contingencia.

Para este ejemplo se consideran 3 casos de estudio.

CASO 1. No incluye las restricciones de capacidad mxima de transporte de potencia

activa en el sistema de transmisin. Este caso, conocido como despacho econmico
(DE), se ejecuta para comparar y validar la tcnica de optimizacin propuesta con
resultados presentados en la literatura [3].

CASO 2. En ste se resuelve un simple FPO, el que permite, a travs del ndice de
desempeo PI, determinar las contingencias ms probables.

CASO 3. Se calculan las soluciones al problema de FPO-RS, considerando las

contingencias con los 2 valores ms altos del ndice de desempeo PI.

Para los tres casos de estudio, los valores de los parmetros del OEP-OR utilizados
son presentados en la Tabla 4.6, [4]-[5].

El caso de DE se plante con fines de validar la tcnica de optimizacin planteada. La

referencia [3] presenta resultados obtenidos para el sistema de prueba considerado,
empleando variantes del OEP y una variante de algoritmo gentico (AG). La Tabla 4.7
reproduce los resultados reportados en dicha referencia.

Para este caso, el trmino senoidal de la curva de costos de generacin no se

considera ya que la referencia indicada tampoco lo hace.
Las tcnicas de optimizacin heursticas, al manejar variables aleatorias, calculan
soluciones prximas al ptimo global que varan de una corrida a otra, por lo que se
requiere ejecutar varias veces el mismo problema de optimizacin para encontrar la
mejor solucin. Datos estadsticos, como la desviacin estndar, permiten comprobar
el comportamiento y estabilidad del algoritmo. Para el caso planteado, los datos
estadsticos presentados en la Tabla 4.8 son resultado de ejecutar el algoritmo de
optimizacin un total de 100 veces. Los valores de las variables de control para el
mejor caso obtenido se despliegan en la Tabla 4.9.
La Tabla 4.9 indica tambin el tiempo de cmputo empleado para cada uno de los

casos planteados. Este tiempo incluye el empleado por el algoritmo de optimizacin y

el algoritmo de flujos de potencia usando la tcnica de Newton Raphson. El trmino
senoidal se incluye en la funcin objetivo para los problemas de FPO y FPO-RS.

TABLA 4.9
ASIGNACIN PTIMA DE LAS VARIABLES DE CONTROL , SISTEMA DE 6 GENERADORES

Unidad

FPO-RS
pre
Lneas salidas
3 - 13
17 - 18
contingencia
448.2
493.3
449.4
199.6
200.0
200.0

DE

FPO

G,1
G,2

446.4
173.1

446.7
166.3

G,3

261.9

255.8

197.7

117.7

209.7

G,4
PG,5

136.5
171.1

150.0
166.6

148.4
164.6

150.0
199.4

119.1
200.0

G,6

86.1

90.3

120.0

120.0

110.2

Tap1

0.90

1.03

1.04

0.90

1.10

Tap2

0.97

1.10

0.95

0.91

1.09

Tap3

0.94

0.90

0.99

1.04

0.99

Tap4

0.98

0.98

0.90

1.07

0.90

0.95

1.10
1.10

0.93

Tap6

1.00
0.96

0.96

1.10

1.08
0.90

Tapy
Shunt-,

0.97
3.5

0.90
0.0

0.90
5.0

0.90
0.0

0.90
0.0

Shunt2

1.5

0.0

0.0

5.0

0.0

Shunt3

0.0

0.0

0.0

0.0

shunt4

5.0

5.0

4.5

5.0

shunt5

1.5
1

4.5
2.0

5.0
5.0

4.5
0.0

0.0
4.5

0.0
0.0

5.0

4.5

Shunt8

0
5

0.0

5.0

0.0
2.5

Shunt9

5.0

0.0

4.5

2.0

15662
15.5

15874
17.5

15770
25.3

P
P
P
P

Tap5

shunt6
Shunt7

Costo ($/h)
Prdidas
tiempo [min]

15436 15523
12.1
5.0

12.8
7.0

16.0

De la Tabla 4.9, la salida de la lnea que une los nodos 3-13 provoca que el

despacho obtenido en estado pre-contingencia para la unidad de generacin

conectada en el nodo 3 disminuya con respecto a los obtenidos para los casos de
DE y FPO, incrementado el costo de generacin. Este redespacho de generacin
permite que en el caso de ocurrir dicha contingencia, la unidad pueda modificar su
nivel de generacin satisfaciendo los cambios de velocidad de generacin y lmites
de flujos de potencia activa en las lneas de transporte.

7. RESTRICCIONES Y PROBLEMAS DE LA OPTIMIZACIN

7.1. Restricciones:
Las variables del sistema elctrico sujetas a lmites impuestos tanto por los propios
equipos que constituyen el sistema elctrico, como por los criterios de operacin del
sistema.
Cabe distinguir, por tanto, entre dos tipos de lmites:

Lmites sobre las variables del sistema que son objeto de control por parte de los
operadores, principalmente:
Potencias activa y reactiva de los generadores, esta ltima normalmente
asociada e intercambiable con la tensin en bornes del generador.
Relacin de transformacin de los transformadores con capacidad de
cambiar tomas en carga, teniendo en cuenta el carcter discreto de las
propias tomas de regulacin.
Bateras de condensadores y reactancias, equipos con capacidad de
generar

consumir

potencia

reactiva

en

incrementos

discretos

determinados por la conexin/desconexin de los distintos elementos.

Lmites de operacin, principalmente:

Tensin de los nudos.
Flujos de potencia por lneas y transformadores, cuantificada en trminos
de potencia aparente, activa, reactiva o incluso en forma de intensidad
mxima.

7.2. Restricciones de desigualdad:

Las restricciones de desigualdad consideran los lmites que deben satisfacer las
variables de control y estado. Estas restricciones reflejan los lmites operativos
impuestos a los dispositivos y al sistema elctrico de potencia. Las principales
restricciones de desigualdad consideradas en un FPO son:

Lmites de potencia activa y reactiva de generacin.- La potencia activa y reactiva

asignada a las unidades de generacin deben estar dentro de lmites de
operacin. Esta restriccin refleja los lmites operativos y trmicos que un
generador debe satisfacer; se les modela por (2.3).

Donde, PGi_MIN y PGi_MAX son los lmites de potencia activa mnima y mxima,
respectivamente, para el i-simo generador; QGi_MIN y QGi_MAX son los lmites
de potencia reactiva mnima y mxima, respectivamente, para el i-simo
generador.

7.3. Restricciones De Seguridad

Hasta ahora solamente se han considerado las restricciones de igualdad en el
algoritmo propuesto, sin embargo, existen lmites mnimos y mximos tanto en los
parmetros de control como en las variables independientes.

a) Restricciones en los parmetros de control:

Los parmetros de control pueden variar en un rango definido por
Estas restricciones de desigualdad pueden ser manejadas fcilmente haciendo un
ajuste en el algoritmo. En este caso al aplicar (15), si alguno de los nuevos
parmetros de control viola uno de sus lmites, este es ajustado a un valor igual
al lmite violado como se indica en (17).

b) Restricciones en las variables independientes:

Las restricciones en las variables independientes son de la forma

Para manejar este tipo de restricciones se ha utilizado un mtodo de penalidad, el

cual consiste en sumar a la funcin objetivo penalidades asociadas a la violacin
de los lmites en estas variables. La funcin objetivo penalizada queda de la forma:

Donde:

En este caso [s] es el vector de factores de penalidad. A estos factores se les

asigna un valor pequeo al inicio del proceso iterativo, el cual se incrementa en
cada una de las iteraciones. Esto hace que a medida que el algoritmo se acerque
al ptimo, las variables independientes tiendan a permanecer dentro de sus
lmites.
7.4. El Problema De Flujos De Potencia ptimo
Un estudio de flujos de potencia ptimo es utilizado ampliamente en la industria
elctrica para diferentes aplicaciones, que van desde estudios de planeacin hasta
operacin de los sistemas. El principal objetivo de un FPO es optimizar las
condiciones de operacin en estado estacionario de un sistema elctrico de potencia.
Un FPO ajusta las cantidades controlables para optimizar una funcin objetivo
mientras satisface un conjunto de restricciones operativas.
Una funcin objetivo puede incorporar aspectos econmicos, de seguridad o
medioambientales, que se resuelve utilizando tcnicas de optimizacin adecuadas.
Las restricciones son leyes fsicas que gobiernan a los generadores, el sistema de
transmisin, lmites constructivos de los equipos elctricos y estrategias operativas.
Esta clase de problema es expresado como un problema de programacin no lineal,
con la funcin objetivo expresada como una funcin no lineal, y las restricciones
expresadas como ecuaciones lineales y no lineales.

En el FPO convencional no se consideran restricciones de diferencias angulares. Sin

embargo, el sistema que est operando en un estado determinado por el FPO puede
no preservar la estabilidad transitoria cuando el sistema experimenta una
contingencia. Como los nuevos mercados elctricos estn inmersos en medios
competitivos, y el costo por prdidas de sincronismo como resultado de una
inestabilidad transitoria es muy alto, se ha vuelto un aspecto importante a considerar
en la operacin y planeacin de las redes elctricas, por lo que en los ltimos aos se
ha visto la necesidad de incluir tales aspectos dentro del problema de flujos de
potencia ptimo.

8. DESCRIPCION GENERAL DE UNA RED ELCTRICA

En esta seccin se describen las redes elctricas que representan a los SEP y se
detallan los elementos componentes de las redes que participan en la formulacin del
OPF, a travs de parmetros y variables de optimizacin.
Las redes elctricas se basan en la interconexin entre generadores y consumos, a
travs de lneas de transmisin y transformadores, con el fin de abastecer las cargas y
brindar una calidad de suministro adecuada; en sentido de proporcionar confiabilidad a
la red a travs de la conformacin de una red enmallada.

8.1. Generacin de energa Elctrica:

La generacin de energa elctrica, se realiza principalmente a travs de la
conversin de energa mecnica en elctrica. En el SIC coexisten tres tipos
principales de tecnologas para generar energa elctrica:

Generacin Hidroelctrica,
Generacin Trmica Convencional,
Generacin Trmica de Ciclo Combinado.

8.1.1. Generacin hidroelctrica:

La generacin hidroelctrica se basa predominantemente en la conversin de
la energa potencial que posee el agua acumulada en embalses. El agua,
conducida a travs de tuberas, es aprovechada para mover turbinas, en las
cuales se genera la energa elctrica.

8.1.2. Generacin trmica convencional:

La generacin trmica convencional se basa, en la utilizacin de algn
combustible fsil, para permitir el movimiento de las turbinas a vapor o gas y
la consecuente generacin de energa elctrica.

8.1.3. Generacin trmica de ciclo combinado:

La generacin trmica de ciclo combinado se basa en la combinacin de una
turbina a gas, cuyo calor excedente es aprovechado en una turbina a vapor,
incrementndose en forma importante la eficiencia del proceso. Existen otros
tipos de generacin, adems de los anteriores, como la generacin nuclear,
la elica, solar y celdas de combustible entre otras. Sin embargo, el SIC est
compuesto principalmente por generacin de los primeros tres tipos.

Para efectos del este trabajo, se considera que una red elctrica est compuesta
principalmente por los siguientes seis elementos: barras o nodos de la red, cargas
o consumos, generadores, lneas de transmisin, transformadores y equipos
FACTS SVC.
Una barra permite la interconexin de componentes elctricas, representa un nodo
de un SEP. Su parmetro ms importante es el tensin, representado por su
magnitud |V| y ngulo , stas son variables de estado para OPF .Las cargas o

consumos, son las demandas por potencia de los usuarios de la red elctrica.
Estas demandas sern, en nuestro caso, datos para OPF. Son representadas por
PL y QL, siendo PL la potencia activa y QL la potencia reactiva demandada por
una carga. Las cargas se conectan a una barra especfica.

8.2. Generador:
Un generador corresponde a una mquina elctrica que es capaz de inyectar
potencia (activa P y reactiva Q) a la red, con el fin de abastecer los consumos. Sus
parmetros son PG y QG y son variables a optimizar por OPF. Este elemento
posee rangos de operacin determinados por caractersticas propias de cada
mquina y determinarn restricciones para OPF. Las lneas de transmisin y
transformadores, son los encargados de transmitir la potencia entre los nodos o
barras de la red elctrica. Sus parmetros estn resumidos en una matriz compleja
denominada matriz admitancia nodal [Y], la cual representa las admitancias entre
los nodos de la red, a partir de un modelo PI asociado a cada lnea de transmisin
o transformador.

8.3. Transformadores:
Los transformadores, son elementos que tienen la capacidad de transformar
tensiones alternas. Adems, pueden cambiar su razn de transformacin a travs
de los denominados taps, que dependiendo del tipo de transformador, pueden ser
manipulados de distintas formas por los operadores de la red. En el presente
estudio se modelan los taps de los transformadores como variables de
optimizacin continuas e inciden directamente sobre la matriz de admitancia.
El flujo de potencia ptimo consiste en despachar una serie de generadores para
satisfacer una demanda dada teniendo en cuenta las restricciones del sistema y

hacindolo al menor costo posible. La solucin general de un OPF consiste en

minimizar o maximizar una funcin objetivo sujeta a restricciones de igualdad y
desigualdad. La funcin objetivo puede tomar diferentes formas dependiendo del
propsito del OPF. Algunos ejemplos de funcin objetivo incluye la minimizacin de
un parmetro como las prdidas del sistema, el costo de operacin y emisiones.
Otras formas de funcin objetivo pueden contemplar la maximizacin de un
parmetro como el beneficio social neto.
El problema de flujo de potencia ptimo se puede abordar usando un modelo
simplificado de la red, lo cual da origen a un flujo de potencia ptimo dc.
El modelo dc de la red se obtiene ignorando las perdidas y linealizando las
ecuaciones de flujo de potencia. Este modelo simplificado naci a partir de los
trabajos de stott y alssac en OPF con restricciones de seguridad. Para resolver un
DCOPF se puede utilizar tcnicas de programacin lineal o programacin
cuadrtica.

CONCLUSIONES
El principal inconveniente al utilizar un paso constante es que no se conoce con
anticipacin un valor ptimo de ste. Esta situacin puede llevar al algoritmo a una
convergencia prematura, o realizar ms iteraciones de las necesarias. Al usar un paso
variable calculado a partir del mtodo de ajuste de parbola, se encontr que en general
se reduce el nmero de iteraciones.
El inconveniente de este mtodo es que se deben encontrar tres puntos de la funcin
objetivo cuya trayectoria pueda ser aproximada por una par abola. Para encontrar estos
puntos se debe utilizar un proceso iterativo y correr varios flujos de carga, lo cual hace
que el esfuerzo computacional de cada iteracin sea mayor.

Actualmente, el proceso de despacho de carga se realiza considerando un modelo de

costos marginalista. En este esquema las restricciones tcnicas del sistema se consideran
en forma simplificada, llevando a soluciones no necesariamente ptimas.
La utilizacin del OPF llevar a mejoras en este proceso, dado que generaliza el
despacho econmico de carga convencional utilizado. Esto se debe a que incluye
restricciones tcnicas, tales como: balances de flujo de potencia activa y reactiva, lmites
tcnicos de operacin de los equipos y limitaciones en la mantencin del factor de
potencia en la liberacin de carga. Esto se ver plasmado en una disminucin de costos
al determinar y fijar en forma ptima las variables de control de un SEP.

Se ha podido comprobar empricamente, que la introduccin adecuada de equipos SVC

entendidos como caso particular de los FACTS (Flexible AC Transmission Systems) en
una red de transmisin de energa elctrica, llevara a menores costos totales por
concepto de ahorro de generacin de potencia activa.

BIBLIOGRAFA
Anlisis y control de Sistemas elctricos de potencia. Gabriel Arguello Ros
Sistemas Elctricos de Potencia. Antonio Gmez Expsito
Tesis FLUJO PTIMO DE POTENCIA. Autor Francisco Javier Vinueza Gmez.
ESCUELA POLITCNIACNACIONAL - FACULTAD DE INGENIERA ELCTRICA
Tesis SOLUCIN DEL PROBLEMA DE GFLUJOS DE POTENCIA PTIMO CON
RESTRICCIONES DE SEGURIDAD POR UN OPTIMIZADOR DE PARTCULAS
MODIFICADO. Autor Pablo enrique Oate Yumbla. CENTRO DE INVESTIGACIN
Y DEESTUDIOS AVANZADOS DEL I.P.N. - UNIDAD GUADALAJARA
ESPECIALIDAD INGENIERA ELCTRICA