COMPENSACIÓN ÓPTIMA DE POTENCIA REACTIVA EN SISTEMAS

ELÉCTRICOS DE TRANSMISIÓN MEDIANTE BÚSQUEDA ARMÓNICA

 UNIVERSIDAD POLITÉCNICA SALESIANA
SEDE QUITO

CARRERA:
INGENIERÍA ELÉCTRICA

Trabajo de titulación previo a la obtención del título de

INGENIERO ELÉCTRICO

TEMA:
COMPENSACIÓN ÓPTIMA DE POTENCIA REACTIVA EN SISTEMAS
ELÉCTRICOS DE TRANSMISIÓN MEDIANTE BÚSQUEDA ARMÓNICA

AUTOR:
PABLO ANDRES ACHIG SANTAMARIA

TUTOR:
DIEGO FRANCISCO CARRIÓN GALARZA

Quito, enero 2019

Pablo Andres Achig Santamaria

COMPENSACIÓN ÓPTIMA DE POTENCIA REACTIVA EN SISTEMAS

ELÉCTRICOS DE TRANSMISIÓN MEDIANTE BÚSQUEDA ARMÓNICA.

Universidad Politécnica Salesiana

Ingeniería Eléctrica

Breve reseña historia e información de contacto:

Pablo Andres Achig Santamaria (Y'1993-M'11). Realizó sus

estudios secundarios en el colegio Técnico Salesiano Don Bosco,
se graduó de Bachiller Técnico Industrial en Electricidad-
Electrónica. Egresado de la Carrera de Ingeniería Eléctrica de la
Universidad Politécnica Salesiana. Su trabajo se basa en la
compensación óptima de potencia reactiva en sistemas eléctricos
de transmisión mediante búsqueda armónica.

[email protected]

Dirigido por:

Diego Francisco Carrión Galarza (Y'1981-M'12). Se graduó en

Ingeniería Eléctrica de la Universidad Politécnica Salesiana,
Ecuador en 2010 y en la actualidad está trabajando para lograr su
título de doctor en ingeniería de la Universidad Politécnica
Bolivariana – Medellín, Colombia. Es profesor e investigador de la
Universidad Politécnica Salesiana – Quito, Ecuador. En la
actualidad es miembro del grupo de investigación GIREI (Grupo
de Investigación en Redes Eléctricas Inteligentes – Smart Grid
Research Group).
[email protected]

Todos los derechos reservados:

Queda prohibida, salvo excepción prevista en la ley, cualquier forma de reproducción,
difusión de este texto con fines académicos o investigativos por cualquier medio, con la
debida notificación a los autores.
DERECHOS RESERVADOS ©2019
Universidad Politécnica Salesiana
QUITO-ECUADOR

I
 DECLARATORIA DE COAUTORÍA DEL DOCENTE TUTOR/A

Yo, Diego Francisco Carrión Galarza declaro que bajo mi dirección y asesoría fue
desarrollado el trabajo de titulación “Compensación Óptima de Potencia Reactiva en
Sistemas Eléctricos de Transmisión Mediante Búsqueda Armónica” realizado por Pablo
Andres Achig Santamaria, obteniendo un producto que cumple con todos los requisitos
estipulados por la Universidad Politécnica Salesiana para ser considerados como trabajo
final de titulación.

Quito, enero de 2019

Diego Francisco Carrión Galarza

Cédula de identidad: 1713703062

II
 CESIÓN DE DERECHOS DE AUTOR

Yo, Pablo Andres Achig Santamaria, con documento de identificación N°

1722976089, manifiesto mi voluntad y cedo a la Universidad Politécnica Salesiana
la titularidad sobre los derechos patrimoniales en virtud de que soy autor/es del trabajo
de grado intitulado: “Compensación Óptima de Potencia Reactiva en Sistemas Eléctricos
de Transmisión Mediante Búsqueda Armónica”, mismo que ha sido desarrollado para
optar por el título de: Ingeniero Eléctrico, en la Universidad Politécnica Salesiana,
quedando la Universidad facultada para ejercer plenamente los derechos cedidos
anteriormente.
En aplicación a lo determinado en la Ley de Propiedad Intelectual, en mi condición de
autor me reservo los derechos morales de la obra antes citada. En concordancia, suscribo
este documento en el momento que hago entrega del trabajo final en formato impreso y
digital a la Biblioteca de la Universidad Politécnica Salesiana.

Firma

Pablo Andres Achig Santamaria

Cédula: 1722976089
Quito, enero de 2019

III
1. ÍNDICE GENERAL

Resumen ............................................................................................................................. 1
Abstract .............................................................................................................................. 1
1. Introducción ................................................................................................................... 2
2. Compensacion de Potencia Reactiva en Sistemas Electricos de Transmisión ............... 4
2.1 Compensadores en Serie .......................................................................................... 5
2.2 Compensadores en Paralelo ..................................................................................... 5
2.3 Compensadores Combinados ................................................................................... 5
3. Modelacion de Dispositivos FACTS ............................................................................. 6
3.1 TSCS ........................................................................................................................ 6
3.2 SVC .......................................................................................................................... 6
3.3 UPFC ........................................................................................................................ 7
4. Metaheuristica de Búsqueda Armónica ......................................................................... 7
5. Formulación del problema ............................................................................................. 9
6. Análisis de resultados ................................................................................................... 12
7. Conclusiones ................................................................................................................ 18
8. Trabajos Futuros .......................................................................................................... 19
9. Referencias ................................................................................................................... 20
9.1 Anexos ................................................................................................................... 24
10. Estado del Arte ........................................................................................................... 26
2. ÍNDICE DE FIGURAS

Figura 1. Ubicación Óptima de FACTS en Sistemas Eléctricos de Transmisión Mediante

Tecnicas de Optimización .................................................................................................. 4
Figura 2. Modelo Matemático de Dispositivos TCSC ....................................................... 6
Figura 3. Modelo Matemático de Dispositivos SVC y UPFC ........................................... 6
Figura 4. Modelo de 14 barras de la IEEE ....................................................................... 12
Figura 5. Ubicación de Dispositivos SVC en el Modelo de 14 Barras de la IEEE .......... 13
Figura 6. Potencia Reactiva en el Modelo de 14 Barras de la IEEE ................................ 14
Figura 7. Perfil de Tensión en el Modelo de 14 Barras de la IEEE ................................. 14
Figura 8. Ángulo de Tensión en el Modelo de 14 Barras de la IEEE .............................. 15
Figura 9. Modelo de 30 barras de la IEEE ....................................................................... 15
Figura 10. Ubicación de Dispositivos SVC en el Modelo de 30 Barras de la IEEE ........ 16
Figura 11. Potencia Reactiva en el Modelo de 30 Barras de la IEEE .............................. 17
Figura 12. Perfil de Tensión en el Modelo de 30 Barras de la IEEE ............................... 17
Figura 13. Ángulo de Tensión en el Modelo de 30 Barras de la IEEE ............................ 18
Figura 14. Temática Estado del Arte................................................................................ 32
Figura 15. Formulación del Problema Estado del Arte .................................................... 32
Figura 16. Solución Estado del Arte ................................................................................ 32
3. ÍNDICE DE TABLAS

Tabla 1. Potencia Reactiva Flujo Óptimo de Potencia Sistema de 14 Barras GAMS ..... 13
Tabla 2. Capacidad de los Compensadores SVC en el Flujo Óptimo de Potencia Sistema
de 14 Barras GAMS ......................................................................................................... 14
Tabla 3. Porcentaje de Mejora de Potencia Reactiva Flujo Óptimo de Potencia Sistema de
14 Barras GAMS .............................................................................................................. 14
Tabla 4. Porcentaje de Mejora de Perfil de Tensión Flujo Óptimo de Potencia Sistema de
14 Barras GAMS .............................................................................................................. 15
Tabla 5. Potencia Reactiva Flujo Óptimo de Potencia Sistema de 30 Barras GAMS ..... 16
Tabla 6. Capacidad de los Compensadores SVC en el Flujo Óptimo de Potencia Sistema
de 30 Barras GAMS ......................................................................................................... 17
Tabla 7. Porcentaje de Mejora de Potencia Reactiva Flujo Óptimo de Potencia Sistema de
30 Barras GAMS .............................................................................................................. 17
Tabla 8. Porcentaje de Mejora de Perfil de Tensión Flujo Óptimo de Potencia Sistema de
30 Barras GAMS .............................................................................................................. 18
Tabla 9. Datos Líneas de Transmisión Modelo de 14 Bus IEEE .................................... 24
Tabla 10. Datos de los Límites de Generación Modelo de 14 Bus IEEE ........................ 24
Tabla 11. Datos de Carga en Nodos Modelo de 14 Bus IEEE ......................................... 24
Tabla 12. Datos Líneas de Transmisión Modelo de 30 Bus IEEE ................................... 24
Tabla 13. Datos de los Límites de Generación Modelo de 30 Bus IEEE ........................ 25
Tabla 14. Datos de Carga en Nodos Modelo de 30 Bus IEEE ......................................... 25
Tabla 15. Matriz Estado del Arte ..................................................................................... 31
COMPENSACIÓN ÓPTIMA DE POTENCIA REACTIVA EN
SISTEMAS ELÉCTRICOS DE TRANSMISIÓN MEDIANTE
BÚSQUEDA ARMÓNICA

Resumen Abstract
En este trabajo se presenta una This paper presents a methodology for the
metodología para la compensación de los compensation of electrical power systems
sistemas eléctricos de potencia mejorando by improving the voltage profile of the
el perfil de tensión del sistema de power system, for which an analysis of the
potencia, para lo cual se realiza un análisis behavior before and after the
del comportamiento antes y después de la compensation is carried out by means of
compensación mediante flujos óptimos de optimal flows of AC power in the models
potencia AC en los modelos de 14 y 30 of 14 and 30 bars of the IEEE. The optimal
barras del IEEE. La compensación óptima compensation of reactive power will be
de potencia reactiva se realizará mediante made by the location of FACTS devices in
la ubicación de dispositivos FACTS SVC the bars of a power system, allowing to
en las barras de un sistema de potencia, improve the reactive component, the
permitiendo mejorar la componente voltage profile in the nodes and the
reactiva, el perfil de voltaje en los nodos y voltage angle. The mathematical
el ángulo de tensión. La modelación modeling of the reactive compensation
matemática de la compensación reactiva will be done through the metaheuristic
se realiza mediante la metaheurística known as "Harmonic Search"; with which
conocida como “Búsqueda Armónica”; costs are minimized by the optimal
con el cual se minimiza los costos por la location of FACTS devices in the SEP,
ubicación optima de dispositivos FACTS with restrictions of active nodal power,
en el SEP, con restricciones de potencia reactive and magnitude of voltage.
nodal activa, reactiva y magnitud de Finally, the optimal compensation in the
voltaje. Finalmente, la compensación bars of a SEP improves the voltage profile
óptima en las barras de un SEP mejora el and the quality of the electric service.
perfil de tensión y la calidad del servicio
eléctrico.

Palabras Clave: Búsqueda Armónica, Keywords: Harmony Search, FACTS,

FACTS, Flujo de Potencia Óptimo, Perfil Optimal Power Flow, Voltage Profile,
de Tensión, Compensación de Potencia Optimal Location, Compensation of
Reactiva. Reactive Power.

1
1. Introducción capacidad de transferencia de energía
eléctrica [3].
Con el transcurso del tiempo el consume
de energía eléctrica ha incrementado Los dispositivos FACTS cuentan con la
sustancialmente, por el contrario, la ventaja de variar el control de potencia
expansión de la generación y transmisión reactiva y el control potencia activa y así
de energía ha sido gravemente poder ajustarse al mismo tiempo al
condicionada, este condicionamiento ha control de magnitud de voltaje debido a
sido debido a impedimentos ambientales su control rápido y su gran elasticidad.
y limitaciones económicas. Como La ubicación óptima y adecuada de estos
resultado de estas limitaciones muchas dispositivos no permite contar con los
de las líneas de transmisión se beneficios propios de los dispositivos
encuentran en sus límites máximos de FACTS [4].
operación, cargabilidad y estabilidad lo
cual convierte en un factor limitante al Los dispositivos FACTS también
sistema de transferencia de potencia [1]. conocidos como Sistemas Flexibles de
Transmisión AC, tienen la característica
La condición del sistema eléctrico de de brindar una gran mejoría al
transmisión para trasladar energía incrementar la capacidad de transmisión
eléctrica se ve severamente afectado por eléctrica del sistema y la elasticidad del
las diferentes limitaciones de estado control del flujo de potencia eléctrica [5].
dinámico o limitaciones de estado Los dispositivos FACTS son
estacionario tales como: estabilidad convertidores mediante la
dinámica, estabilidad transitoria, implementación de un sistema
estabilidad angular, límites térmicos y electrónico de potencia que tiene las
magnitud de tensión. Estas restricciones condiciones de controlar una gran
y limitaciones son las que determinan los variedad de parámetros eléctricos en
valores máximos de potencia eléctrica sistemas de transmisión eléctrica e
que se puede transmitir sin causar ningún instalaciones, tanto en el flujo de
desperfecto a los equipos eléctricos y a potencia en estado dinámico en estado
las líneas de transmisión [2]. estacionario [6]. Los tipos de FACTS
son: Compensador Estático
En el principio, las restricciones en los (STATCOM), Compensador Estático de
medios de transferencia de potencia Potencia Reactiva (SVC), Compensador
eléctrica se aligeraban por medio del de Serie Controlado por Tiristor (TCSC),
reacondicionamiento e instalación de Controlador de Flujo de Potencia
nuevas líneas de transmisión, generación Unificado (UPFC) y otros [7].
y distribución eléctrica. Opcionalmente,
los dispositivos (FACTS) o Sistemas El dispositivo FACT más utilizado en las
Flexibles de Transmisión AC, puede redes de los sistemas eléctrico de
llegar a cumplir las mismas metas sin potencia es el SVC, esto es debido a su
alteraciones mayores en la red eléctrica. excelente desempeño mejorando la
Los beneficios más importantes de los estabilidad de la red. Este dispositivo es
dispositivos FACTS son: menor costo de un compensador estático de potencia
inversión en la operación, reactiva conectado en paralelo con una
mantenimiento y transmisión, mayor salida variable, la cual os permite el
confiabilidad y seguridad el sistema intercambio de corriente tanto capacitiva
eléctrico de potencia, y aumento de la como inductiva para así proporcionar
con estabilidad de tensión en el sistema
eléctrico de potencia. Al instalar en una

2
ubicación óptimo el dispositivo SVC nos  Análisis de costos utilizando el
permite disminuir las pérdidas de energía flujo óptimo de potencia.
en el sistema eléctrico [8].  Análisis de estabilidad de
tensión.
El dispositivo STATCOM es un  Análisis de sensibilidad de
controlador estático de derivación y pérdidas.
además es uno de los compensadores
más fundamentales de los FACTS, ya Las metodologías de optimización
que dichos dispositivos se utilizan cada heurística han sido extensamente
vez con ms frecuencia en líneas de aplicadas en la solución del problema
transmisión largas en sistemas eléctricos de la colocación y ubicación optima de
de potencia actuales [9]. Los dispositivos SVC y STATCOM. En esta
STATCOM pueden contar con varias investigación se presenta un análisis
funcionalidades dependiendo la detallado de las metodologías de
aplicación en el funcionamiento, optimización heurística para la
operación y control del SEP, como en la colocación óptima de STATCOMs y
disminución de la cantidad de SVCs. Este estudio adjunta importantes
componentes asimétricos que metodologías de optimización
amortiguan las variaciones y heurística, tales como: Búsqueda
oscilaciones de potencia eléctrica Armónica (HS), Optimización de
mejorando así la estabilidad transitoria Enjambre de Partículas (PSO),
del sistema, también ayudan con la Recocido Simulado (SA), Algoritmos
programación del flujo óptimo de Genéticos (GA) y algoritmos de
potencia [10]. estrategias de evolución (ES) [14].
Últimamente las metodologías de
Los mayores beneficios de la optimización heurística se han
compensación de potencia reactiva van a convertido en la opción más óptima
de depender de la ubicación optima y la para resolver problemas con gran
capacidad de los compensadores. La complejidad que dado de otro modo
instalación y ubicación de dispositivos sería muy dificultoso resolverlos por
FACTS en todas las barras del sistema otras metodologías tradicionales [15].
eléctrico de potencia es indebido e En los últimos tiempos, se ha venido
innecesario debido a sus consideraciones evolucionando varias metodologías de
económicas [11]. La determinación de la optimización a nivel mundial, las
mejor ubicación de los controladores cuales están basadas en analogías
estáticos de potencia reactiva involucra inspiradas en la naturaleza. La mayoría
el cálculo de las circunstancias de estado de estas metodologías heurísticas están
estacionario para la red del sistema basadas en la población, las cuales son
eléctrico de potencia. No obstante, la conocidas también como algoritmos de
complejidad del problema aumenta propósito general debido a su gran y
debido a la no linealidad de las variada aplicación en una distinta
ecuaciones del flujo de carga [12]. cantidad de problemas [16-25].

Se utilizan una gran variedad de A continuación se observa en la figura

metodologías para optimizar la 1 la ubicación óptima de FACTS en el
colocación de los Sistemas Flexibles de sistema eléctrico de potencia mediante
Transmisión AC en la red. Estas técnicas de optimización.
metodologías pueden ser [13]:

3
 Figura 1. Ubicación óptima de FACTS en el sistema eléctrico de potencia mediante técnicas de optimización.

En esta investigación se propone un A continuación esta investigación está

procedimiento que permite comprobar y desarrollada de la siguiente forma. En la
validar el proceder del flujo óptimo en parte II se trata acerca de la
sistemas eléctricos de potencia, compensación óptima de potencia
proponiendo un caso de estudio de reactiva en un SEP, el cual es el objeto de
optimización para mejorar la estudio del presente documento. En la
compensación reactiva en el SEP parte III se detalla el formulamiento del
mediante la ubicación optima de FACTS, caso de estudio mediante el algoritmo de
mediante la cual es posible mejorar la optimización planteado. En la parte IV se
componente reactiva en un SEP, el perfil examina el resultado de la modelación
de tensión en los nodos de un SEP y el obtenida en la presente investigación. Se
ángulo de la tensión. Para ello se utilizará concluye esta investigación en la parte V
herramientas computacionales como con las conclusiones obtenidas del
Matlab y GAMS, que nos ayudaran a estudio.
presentar un análisis técnico del sistema
eléctrico de transmisión planteado 2. Compensación de Potencia
Se pretende realizar la compensación
Reactiva en Sistemas
reactiva en un flujo óptimo AC en los Eléctricos de Transmisión.
modelos propuesto de 14 y 30 barras del
IEEE, mediante la utilización de la En la actualidad existen varias clases de
metaheurística Búsqueda Armónica la controladores, los cuales se clasifican
cual coloca óptimamente los FACTS. según su manera de instalación y
Para ello se hará uso de un Flujo de conexión a la red del sistema eléctrico de
Potencia Optimo, para calcular los flujos potencia, tales como: Compensadores
reactivos del sistema. Los cuales servirán Paralelo, Compensadores Serie y
para realizar una optimización basada en Compensadores Combinados es decir
costos de instalación de los dispositivos serie y paralelo. A continuación se
FACTS. muestran algunos tipos de
compensadores.

4
2.1 Compensadores en serie. del compensador paralelo se base en
entregar potencia reactiva a la red de
Los controladores serie están transmisión, para incrementar la potencia
representados como una impedancia activa al sistema eléctrico de potencia por
cambiable o como un reactor. El medio de la cual se mejorará el factor de
compensador serie introduce un voltaje a potencia, manteniendo los parámetros de
la línea de transmisión dado a su tensión en los límites de estabilidad y
comportamiento como una impedancia operatividad [27].
cambiable, incrementada por la corriente
introducida a la línea de transmisión. De Los compensadores paralelo son:
la misma forma la tensión está alineada
con la corriente, el compensador - Reactor Controlado por Tiristores
introduce potencia reactiva a la línea. (TCR)
Existen otras combinaciones posibles las - Condensador Conmutado por
cuales dependerán de la transferencia de Tiristores (TSC)
potencia activa en el sistema eléctrico de - Compensador de VAR Estático (SVC)
potencia. De la misma manera el efecto - Compensador Estático de VAR
principal de esta clase de dispositivos es (STATCOM)
controlar el flujo de corriente y al - Generador Síncrono Estático (SVG)
cargabilidad en el SEP [26].
2.3 Compensadores combinados.
Los compensadores serie son:
El controlador combinado une las
aplicaciones de los compensadores serie
- Compensador Serie Conmutado por y de los compensadores paralelos
Tiristores (TSSC) obteniendo de esta manera mejores
- Compensador Serie Controlado por cualidad técnicas y económicas en sus
Tiristores (TCSC) diversas aplicaciones en un sistema
- Compensador Serie de VAR Estático eléctrico de potencia. El controlador más
(SSVC) común es: El Controlador Unificado de
- Compensador Serie Sincrónico Flujos de Potencia. El controlador de
Estático (SSSC) Interlinea y el controlador Unificado
Generalizado tienen la capacidad de
- Reactor serie controlado por monitorear la potencia de dos o más redes
tiristores (TCSR) de transmisión [28].
- Compensadores por ángulo de fase
(TCPAR) Los compensadores combinados son:

2.2 Compensadores en paralelo. - Controlador Interlínea de Flujos de

Potencia, (IPFC)
El controlador en paralelo introduce - Controlador Unificado de Flujos de
corriente reactiva el sistema eléctrico de Potencia Generalizado, GUPFC
potencia permitiendo de esta forma se (GUPFC)
verifica la tensión en la barra, muy - Controlador Unificado de Flujos de
independiente de la red de transmisión a Potencia, UPFC
la que están conectadas al nodo, a
divergencia de un compensador serie, El compensador combinado tiene la
requiere de compensadores unitarios, más característica de realizar una
aun si una línea de trasmisión tiene un compensación que cambia el flujo de la
fallo. De esta manera el funcionamiento potencia de una red de transmisión. No

5
obstante, cada una presenta ciertas
restricciones y limitantes relacionadas a
los reactivos y la conmutación de la
potencia activa con la red de trasmisión.

3 Modelación de Dispositivos
FACTS Figura 2. Modelamiento matemático del dispositivo
TCSC
Los dispositivos FACTS se conforman de
equipos estáticos utilizados para 3.2 Compensador de VAR Estático
aumentar la capacidad de transferir la (SVC).
potencia de la red del SEP. Son también
utilizados para disminuir las pérdidas de Los dispositivos SVCs tienen el
transmisión y de esta forma mejorar el comportamiento de una reactancia
perfil de voltaje del sistema eléctrico de cambiable la cual está conectada en
potencia. paralelo, la misma que se encarga de
generar y absorber potencia reactiva
3.1 Compensador Serie Controlado siendo el objetivo primordial ajustar la
por Tiristores (TSCS). magnitud de la tensión en la conexión al
nodo de la red del sistema eléctrico de
Al cambiar el valor de la reactancia de la potencia. Su uso principal es administrar
línea de transmisión, el TCSC actúa como de forma instantánea potencia reactiva,
un compensador inductivo o capacitivo. dando apoyo a la variación del voltaje,
En este trabajo, el valor de la reactancia controlando las sobrecargas, mejorando
de la línea de transmisión se ajusta directo la confiabilidad del sistema de corriente
por el dispositivo TCSC. La agrupación alterna dando un mejor amortiguamiento
de TCSC va a depender del valor de la ante las variaciones de potencia. Los
reactancia de la línea de transmisión controladores operar de manera
donde se encuentra ubicado el TCSC: coordinada por medio de un monitoreo de
control teniendo una respuesta inmediata.
𝑿𝒊𝒋 = 𝑿𝑳í𝒏𝒆𝒂𝒊𝒋 + 𝑿𝑻𝑪𝑺𝑪𝒊𝒋 (1) Se modelan con una administración de
potencia reactiva en el bus del nodo i, a
𝑿𝑻𝑪𝑺𝑪𝒊𝒋 = 𝒓𝒕 𝒄𝒔𝒄 + 𝑿𝑳í𝒏𝒆𝒂𝒊𝒋 (2) continuación se observa en la figura 3
[30].
Donde:
𝑋𝐿í𝑛𝑒𝑎 =Es la reactancia de la línea de
transmisión.
𝑟𝑡 𝑐𝑠𝑐 = Es el coeficiente que representa
el grado de compensación de TCSC.

Para evitar una sobrecompensación, el

rango de trabajo del TCSC se escoge
entre: Figura 3. Modelo matemático de dispositivos SVC
y UPFC
−𝟎. 𝟕 ≤ 𝑿𝑳í𝒏𝒆𝒂 ≤ 𝟎. 𝟐 (3)

La potencia administrada en la barra del

A continuación se observa en la figura 2 nodo i es:
[29].
∆𝑸𝒊𝒔 = ∆𝑸𝑺𝑽𝑪 (4)
6
Donde: improvisan sus notas musicales en sus
∆𝑄𝑖𝑠 = Es la variación de potencia instrumentos para de esta manera
reactiva en el nodo i en MVAR. encontrar una mejor armonía [32]. De la
∆𝑄𝑆𝑉𝐶 = Es el valor de potencia reactiva misma forma que los instrumentos
del dispositivo SVC en MVAR. musicales pueden ser interpretados con
acordes musicales que están basados en
3.3 Controlador Unificado de Flujos un proceso de improvisación aleatorio,
de Potencia (UPFC). como resultado se puede obtener
variables optimas con ciertos valores
El modelamiento de un UPFC fue basados en procesos aleatorios mediante
propuesto en 1991 por Gyugyi. Este inteligencia artificial o computacional
dispositivo tiene la capacidad de [33].
controlar de una forma selectiva y
simultánea, las variables que afecten el Los intérpretes de la música mejoran sus
flujo de potencia en la red es decir acordes tomando como referencia ciertos
parámetros como: ángulo de fase, estándares, mientras que las variables de
impedancia y voltaje. El dispositivo diseño en la memoria computacional
UPFC está integrado con dos pueden mejorarse tomando como
convertidores de conmutación el cual esta referencia una función objetivo. Las
operado desde un enlace común de funciones del algoritmo de Búsqueda
corriente continua. Armónica se encuentran de forma
continua y discontinua, ya que estas
El rango de trabajo del ángulo UPFC está ecuaciones no requieren gradientes
entre -180° y +180° [31]. diferenciales, tampoco no requieren de na
variación del valor inicial para las
̅ 𝑭𝑨𝑪𝑻𝑺 = ∆𝑼
∆𝑼 ̅ 𝑼𝑷𝑭𝑪 (5) variables, las cuales están libres de
divergencia [34].
Las corrientes inyectadas en el nodo i y el
nodo j se pueden expresar de la siguiente El algoritmo de Búsqueda Armónica se
manera: ha utilizado para determinar la instalación
y ubicación óptima de los dispositivos
̅ 𝑼𝑷𝑭𝑪
∆𝑼 (6) FACTS tales como: SVC, TCSC y UPFC
∆𝑰̅𝒊𝒔 = en la red de un sistema eléctrico de
𝒁𝒊𝒋
potencia para de esta forma mejorar la
confiabilidad, estabilidad y seguridad del
̅ 𝑼𝑷𝑭𝑪
∆𝑼 (7)
∆𝑰̅𝒋𝒔 = − sistema eléctrico de potencia [35]. La
𝒁𝒊𝒋 cargabilidad, sobrecarga de la línea de
transmisión y el sobre voltaje en la barra
Donde: del nodo se resuelven controlando la
∆𝑈̅𝐹𝐴𝐶𝑇𝑆 = Es la variación de potencia potencia activa y reactiva del controlador
reactiva en el nodo i en MVAR. en serie.
∆𝑈̅𝑈𝑃𝐹𝐶 = Es el valor de potencia reactiva
del dispositivo UPFC en MVAR. El algoritmo de búsqueda armónica se ha
implementado actualmente en una
4 Metaheurística de Búsqueda metodología con el proceso de
Armónica improvisación de las notas musicales,
donde cada intérprete de música
El algoritmo de Búsqueda Armónica es improvisa los acordes de sus
una metodología de optimización que instrumentos para así obtener una mejor
está basada en los músicos que armonía musical. [36].

7
 En el algoritmo de búsqueda armónica Paso 3: Asignación.
cada posible resultado es considerado
una “armonía”, la cual es representada Imprimir (“Ingrese el número de
mediante un vector de n posibles iteraciones”)
dimensiones. Las posibles soluciones Leer (NI)
iniciales candidatas son agrupadas e
inicializadas de una forma aleatoria Imprimir (“Ingrese el tamaño de la
dentro de una memoria. De esta manera memoria”)
una nueva posible solución es generada Leer (HMS)
producto de uno de los componentes
que contiene en la memoria armónica Imprimir (“Ingrese el rango de
del algoritmo, por medio de una consideración de la memoria”)
operación aleatoria o mediante una Leer (HMCR)
operación de variación del “tono” de un
vector que este contenido dentro de la Imprimir (“Ingrese el ajuste del tono”)
misma. Finalmente, la memoria Leer (PAR)
armónica se actualiza mediante las
comparaciones del nuevo resultado Imprimir (“Ingrese el ancho de banda”)
candidato y el peor de todos los ventores Leer (BW)
que están contenidos dentro de la
memoria armónica, si el resultado es Imprimir (“Ingrese la función objetivo”)
mejor procederá a remplazar el lugar del Leer (X)
vector dentro de la memoria armónica,
caso contrario se seguirá manteniendo Paso 4: Proceso.
la misma solución resultado candidata y
no abra ninguna variación ni ningún Inicialización de la memoria armónica,
cambio. Este procedimiento se realizará construcción del vector HM y
hasta que se cumpla los parámetros reconocimiento del peor vector en HM.
establecidos de parada. La estructura
principal del algoritmo de búsqueda 𝑿𝟏
armónica se fundamenta en tres etapas: 𝑿 (8)
𝑯𝑴 = [ 𝟐 ]
inicialización de las variables, ⋮
improvisación de las nuevas armonías y 𝑿𝑯𝑴𝑺
la actualización de la memoria armónica
[37]. 𝑿𝟏
𝑿 (9)
𝑿𝒑𝒆𝒐𝒓 ∈ [ 𝟐 ]
⋮
Algoritmo de búsqueda armónica.
𝑿𝑯𝑴𝑺
Búsqueda Armónica Improvisar una nueva armonía.
Paso 1: Inicio del algoritmo. 𝑿′ = ∅ (10)
for 𝒊 = 1, … , 𝑁 do (11)
Paso 2: Declaración. if (𝑼(0,1) ≤ 𝑯𝑴𝑪𝑹) then (12)
begin
HMS, HMCR → numéricos
PAR, BW → numéricos 𝑿𝟏𝒊
𝟐 (13)
NI, X → numéricos 𝑿′𝒊 ∈ 𝑿𝒊
⋮
𝑯𝑴𝑺
[𝑿𝒊 ]
if (𝑼(0,1) ≤ 𝑷𝑨𝑹) then (14)

8
 end Función de Costo
else
𝑿′𝒊 ∈ 𝑿𝒊 (15) De acuerdo con [38], las funciones de
end if costos para TCSC, SVC y UPFC se
end for desarrollan de la siguiente manera:

Actualizar la memoria armónica. TCSC

𝑪𝑻𝑪𝑺𝑪 = 𝟎. 𝟎𝟎𝟏𝟓𝒔𝟐 − 𝟎. 𝟕𝟏𝟑𝟎𝒔
if (𝒇(𝑿′ ) < 𝒇(𝑿 𝒑𝒆𝒐𝒓
) then (16) + 𝟏𝟓𝟑. 𝟕𝟓 (17)

Incluir 𝑿′ a la HM SVC
Excluir 𝑿𝒑𝒆𝒐𝒓 a la HM 𝑪𝑺𝑽𝑪 = 𝟎. 𝟎𝟎𝟎𝟑𝒔𝟐 − 𝟎. 𝟑𝟎𝟓𝟏𝒔
+ 𝟏𝟐𝟕. 𝟑𝟖 (18)
Comprobar las condiciones y criterios de
parada. UPFC
𝑪𝑼𝑷𝑭𝑪 = 𝟎. 𝟎𝟎𝟎𝟑𝒔𝟐 − 𝟎. 𝟐𝟔𝟗𝟏𝒔
+ 𝟏𝟖𝟖. 𝟐𝟐 (19)
while (el criterio de terminación no es
especificado por NI)
Donde:
Regresar a improvisar una nueva armonía 𝐶𝑇𝐶𝑆𝐶 , 𝐶𝑆𝑉𝐶 𝑦 𝐶𝑈𝑃𝐹𝐶 = Están en 𝑈𝑆 $ /
y actualizar la memoria armónica. 𝑘𝑉𝑎𝑟
𝑠 = Es el rango de operación del
end if dispositivo FACT en MVAR.

Acepta los nuevos tonos armónicos, es Restricciones

decir si las soluciones son las mejores.
Para las restricciones planteada se han
Paso 5: Visualizar. tomado en cuenta los valores de los
límites de tensión tanto superior como
Imprimir (“Los nodos de ubicación de inferior sabiendo que siempre los
FACTS son:”, 𝑿′ ) parámetros deben estar dentro de los
márgenes del ±10% de la tensión nominal
Paso 6: Fin del algoritmo. del SEP, de la misma manera la potencia
reactiva del SEP debe encontrarse dentro
5. Formulación del problema de los parámetros establecidos tanto
inferiores como superiores en el cual se
encuentra programado nuestro sistema
En los últimos de tiempos se ha realizado
como valores adecuados para una
la utilización de una gran cantidad de
correcta optimización, de la misma forma
metodologías y algoritmos para dar
la potencia reactiva debe encontrase en
solución a interrogantes y problema de
los límites adecuados del sistema, de esta
ingeniería. Pero los grandes e
forma se logrará una optimización con
inadecuados problemas computacionales
mayor precisión de los valores, a
generados por los métodos numéricos han
continuación se detalla la ecuación (20)
exigido a los investigadores a buscar
con las respectivas restricciones y
nuevas metodologías de investigación y
parámetros de tensión máxima y mínima,
nuevos algoritmos como lo son las
en la ecuación (21) se presenta las
Metaheuristicas, las cuales combinan
restricciones y limitaciones de potencia
aleatoriedad y parámetros en sus procesos
reactiva con un valor máximo y mínimo,
para buscar una solución a un fenómeno
y por último en la ecuación (22) se
natural.

9
presenta las limitaciones y restricciones Función Objetivo
de potencia activa con los valores
máximos y mínimos. El problema objetivo se expresa de la
siguiente manera:
𝑽𝒋 𝒎𝒊𝒏 ≤ 𝑽𝒋 ≤ 𝑽𝒋 𝒎𝒂𝒙 (20)
𝑵

𝑸𝒋 𝒎𝒊𝒏 ≤ 𝑸𝒋 ≤ 𝑸𝒋 𝒎𝒂𝒙 (21) 𝐦𝐢𝐧 ∑ 𝒁𝒊 ∗ 𝑿𝒊 (23)

𝒊=𝟏
𝑷𝒋 𝒎𝒊𝒏 ≤ 𝑷𝒋 ≤ 𝑷𝒋 𝒎𝒂𝒙 (22) Donde:
𝒁𝒊 = Es el costo de la instalación,
Donde: operación, mantenimiento y pérdidas de
Vj = Es el voltaje en el nodo j. los dispositivos FACTS en el nodo i.
Qj = Es la potencia reactiva en el nodo j. 𝑿𝒊 = Es la variable binaria del problema.
Pj = E la potencia activa en el nodo j.
En esta ecuación descrita se toma en
Los superíndices min, max. Son los cuenta la minimización de los
límites mínimos y máximos de las dispositivos FACTS, tomando el valor de
variables. la potencia reactiva la cual es el producto
del análisis del flujo óptimo de potencia,
Ubicación Óptima de Dispositivos con la cual se procederá a calcular el
FACTS. detalle del valor del costo de la ubicación
de los dispositivos FACTS.
Los FACTS podrían instalarse y ubicarse
óptimamente en cualquier parte del ̂
𝒇(𝒙) ≥ 𝟏 (24)
sistema eléctrico de potencia de forma
aleatoria, pero para incrementar y contar Donde:
con una mayor optimización de la 𝒇(𝒙) = Es un vector de función que
ubicación e instalación de los dispositivos representa a las restricciones del
FACTS se debe hacer una selección problema.
adecuada del lugar del sistema eléctrico ̂ = Es un vector que posee entradas
𝟏
de potencia donde deberán ubicarse estos iguales a 1.
dispositivos. Para poder reconocer el
lugar más idóneo y óptimo donde deben En esta ecuación descrita se estima un
ser ubicados los FACTS se toman en valor de 1 si el dispositivo FACTS se
cuenta los siguientes parámetros: instalara en un nodo i y se considerara un
valor de 0 si no se lo hace, la ecuación
- Nivel de Tensión. propuesta es una ecuación únicamente de
- Valor de Potencia Reactiva que variables binarias es decir solo
circula en cada nodo del SEP reconocerá valores de 0 y 1 según
aplique.
Estos valores se los obtiene después de
hacer un análisis y correr un flujo óptimo Por medio de la matriz de conectividad es
de potencia, otro parámetro importante de donde se obtienen las limitaciones y
con el cual debemos contar es la matriz de restricciones, de esta manera se
conectividad del sistema eléctrico de determina las barras del sistema eléctrico
potencia planteado. Con estos valores de potencia en donde podrá ser ubicado
considerados se realizará la optimización un dispositivo FACT, es por tal motivo
del problema descrito. por el que se procede a realizar una
sumatoria de los nodos interconectados
por cada línea de transmisión.

10
 transmisión, así como cada uno de los
𝑮(𝒙) ≤ 𝟎 (25) transformadores modelados en forma de
un reactor. La operación y el computo a
En esta ecuación se asume un valor de 0 detalle de todos los componentes del
si en el nodo i existe conectado algún vector de admitancia (𝑌𝑏𝑢𝑠), divididos
generador. en las siguientes categorías, elementos
pertenecientes a la diagonal principal
Flujos Óptimos de Potencia (𝑌𝑖𝑖) de la misma forma los elementos
que están fuera de la diagonal secundaria
El análisis de un flujo óptimo de potencia (𝑌𝑖𝑗); El elemento de la diagonal es
es de mucha importancia para de esta obtenido por medio del cálculo de la
manera determinar las condiciones y sumatoria de la admitancia de los
limitaciones del sistema eléctrico de componentes que conforman el nodo
potencia existente, al mismo tiempo tiene [39]. El elemento que se encuentra en la
la capacidad de entregar información parte exterior de la diagonal son valores
valiosa para analizar las mejores desfavorables [40]. El vector de
estrategias que conlleven como resultad admitancia es simétrico, como se puede
final al mejoramiento de la operación y observar en siguiente formulación:
eficiencia del sistema eléctrico de
potencia; para el caso del presente estudio 𝑌𝑖𝑗=𝑌𝑗𝑖 (26)
este análisis facilita la ubicación e
instalación de los dispositivos FACTS En la siguiente formulación se observa la
para la compensación de potencia matriz de admitancias [41].
reactiva, indicando de alguna forma con
los valores entregados las ubicaciones 𝐘𝐢𝐢 𝐘𝐢𝐣 ⋯ 𝐘𝐢𝐧
más optimas y adecuadas para la 𝐘𝐣𝐢 𝐘𝐣𝐣 ⋯ 𝐘𝐣𝐧 (27)
colocación de los dispositivos FACTS 𝒀𝒃𝒖𝒔 = [ ]
⋮ ⋮ ⋱ ⋮
para de esta forma mejorar el rendimiento 𝐘𝐧𝐢 𝐘𝐧𝐣 … 𝐘𝐧𝐧
y la confiabilidad del sistema eléctrico de
potencia, mejorando el perfil de tensión y
disminuyendo las pérdidas en la En las ecuaciones (28) y (29) se observa
transmisión de energía eléctrica. Los los componentes de la matriz de
análisis de flujos de potencia óptimos admitancias, los cuales representan a la
entregan como resultado el ángulo y la conductancia (𝐺𝑏𝑢𝑠) y susceptancia
magnitud de fase de la tensión en cada (𝐵𝑏𝑢𝑠) [42].
uno de los nodos de la red del sistema
eléctrico de potencia, las pérdidas de 𝐆𝐢𝐢 𝐆𝐢𝐣 ⋯ 𝐆𝐢𝐧
transmisión en el sistema eléctrico de 𝐆𝐣𝐢 𝐆𝐣𝐣 ⋯ 𝐆𝐣𝐧 (28)
𝑮𝒃𝒖𝒔 = [ ]
potencia y el flujo de potencia reactiva y ⋮ ⋮ ⋱ ⋮
activa en el sistema eléctrico de potencia 𝐆𝐧𝐢 𝐆𝐧𝐣 … 𝐆𝐧𝐧
[38].
𝑩𝒃𝒖𝒔
El parámetro más fundamental a tomar en 𝐁𝐢𝐢 𝐁𝐢𝐣 ⋯ 𝐁𝐢𝐧 (29)
cuenta para la formulación de los flujos 𝐁𝐣𝐢 𝐁𝐣𝐣 ⋯ 𝐁𝐣𝐧
óptimos de potencia es la matriz de =[ ]
⋮ ⋮ ⋱ ⋮
admitancias, la misma que está formada 𝐁𝐧𝐢 𝐁𝐧𝐣 … 𝐁𝐧𝐧
por todos los componentes existentes
entre cada una de las barras y los nodos
de conexión, se consideran los valores de Una vez que se obtiene la matriz de
capacitancia e impedancia de las líneas de admitancia se puede calcular las
potencias activa y reactiva, las cuales son

11
expresadas en las siguientes ecuaciones Modelo de la IEEE de 14 barras.
(30) y (31) respectivamente [43]:
Con el objetivo de simular se lo ha
𝑵
realizado en el modelo establecido de la
𝑷𝒊 = ∑ 𝑽𝒋 ∗ 𝑽𝒊 ∗ (𝑮𝒊𝒋 𝒄𝒐𝒔𝜽 (30) IEEE de 14 barras que es un sistema
𝒊=𝟏
+ 𝑩𝒊𝒋 𝒔𝒆𝒏𝜽) normalizado que tiene las características
de contar con 5 generadores síncronos, 5
transformadores, 14 barras, por último se
𝑵
encuentra a la barra 1 como barra slack
𝑸𝒋 = ∑ 𝑽𝒋 ∗ 𝑽𝒊 ∗ (𝑮𝒊𝒋 𝒔𝒆𝒏𝜽 (31) usando como potencia base 100 MVA, en
𝒋=𝟏
− 𝑩𝒊𝒋 𝒄𝒐𝒔𝜽) la colocación óptima de FACTS en un
SEP, se realizará el respectivo análisis en
el modelo de 14 barras de la IEEE que se
Donde: observa en la figura 4.
𝑃𝑖 = Potencia activa que ingresa al nodo.
𝑄𝑖 = Potencia reactiva que ingresa al
nodo.
𝑉𝑖 = Voltaje en el nodo.
𝐺𝑖𝑗 = Conductancia en los nodos.
𝐵𝑖𝑗 = Susceptancia en los nodos.
𝜃𝑖𝑗 = Diferencia de ángulos del nodo.

6. Análisis de resultados
Para verificar la eficacia del algoritmo de
Búsqueda Armónica para la colocación
de FACTS en redes eléctricas de
transmisión, se hicieron uso de las redes
de prueba de la IEEE como el modelo de
14 y 30 barras. Se evaluaron aspectos
como el valor económico de la
colocación de los FACTS, la disminución
de energía eléctrica en la red y si había
presencia de variaciones de voltaje,
cuidando los parámetros de voltaje en Figura 4. Modelamiento de 14 barras de la IEEE
cada una de las barras del SEP.
El flujo óptimo de potencia es calculado
Los parámetros de inicialización tales mediante el optimizador GAMS, Para lo
como la velocidad de ajuste de banda, el cual se utilizan los datos presentados de
ancho de banda, el rango de variables y la potencia reactiva en [44]. La potencia
tasa de consideración de armonía se reactiva determina el costo de instalación
mantienen constantes para cada caso. Los de los FACTS.
parámetros que varían son el número de
variables, el número de restricciones, el
número de memorias armónicas y el
número de iteraciones ya que eso
depende directamente de cada uno de los
modelos.

12
Los valores observados en la tabla 1 son
producto del resultado del flujo óptimo de
potencia los cuales son resultado de
utilizar la herramienta computacional
GAMS.
Tabla 1. Potencia reactiva flujo óptimo de potencia
sistema de 14 barras GAMS.

POTENCIA
Nodos
REACTIVA [P.U]
N1 -9.114
N2 -9.470
N3 -0.820
N4 -20.990
N5 -21.768
N6 -10.866
N7 -11.538
N8 0
N9 -21.160
N10 -3.950
N11 -0.018
N12 -1.753
N13 -5.217 Figura 5. Ubicación de dispositivos SVC en el
N14 -1.749 modelo de 14 barras de la IEEE

Para validar la heurística propuesta se

El modelo de la IEEE de 14 barras posee compara el resultado obtenido en
14 variables binarias una para cada barra GAMS con el que se obtuvo en Matlab.
del SEP y 15 restricciones. La ubicación que determino GAMS de
los dispositivos FACTS en las barras del
Para validar los resultados de la heurística modelo de 14 barras propuesto por la
propuesta se compara el resultado de la IEEE es igual a la ubicación que
ubicación óptima de FACTS obtenida en determino Matlab. También se validó
Matlab con el resultado de la ubicación esta heurística mediante la comparación
óptima de FACTS obtenida en GAMS. de resultados que se obtiene en el
modelo de 14 barras de la IEEE
Ubicación Óptima de SVC en el propuesto por [45]. La capacidad de los
modelo de la IEEE de 14 barras. compensadores que se deberán instalar
en las barras 4, 5, 7, 11 y 13 del modelo
En la figura 5 se muestra la ubicación de 14 bus de la IEEE.
optima resultante de la metaheurística
Búsqueda Armónica, se observa la Los datos presentados en la tabla 2
ubicación de los SVC en color rojo, en las hacen referencia a la capacidad de los
barras 4, 5, 7, 11 y 13 del modelo de 14 compensadores SVC en las barras 4, 5,
bus de la IEEE. 7, 11 y 13 de modelo de 14 bus de la
IEEE.

13
Tabla 2. Capacidad de los compensadores SVC en el Tabla 3. Porcentaje de mejora de potencia reactiva
flujo óptimo de potencia sistema de 14 barras GAMS con la ubicación de SVC en el flujo óptimo de
. potencia sistema de 14 barras GAMS.

POTENCIA REACTIVA
Nodos
[P.U] MVAR POTENCIA POTENCIA
N4 -20.990 REACTIVA REACTIVA %
Nodos
N5 -21.768 [P.U] [P.U] MEJORA
N7 -11.538 SIN SVC CON SVC
N11 -0.018 N1 -9.114 -8.851 2.88 %
N13 -5.217 N2 -9.470 -9.466 0.04 %
N3 -0.820 0.210 74,39 %
N4 -20.990 -20.990 0%
En la figura 6 se muestra que al tener N5 -21.768 -21.768 0%
dispositivos SVC instalados en el SEP, N6 -10.866 -7.939 26.93 %
existe una mejor distribución de la N7 -11.538 -11.538 0%
potencia reactiva en las barras de cada N8 0 0.156 84.4 %
uno de los nodos, aliviando así el flujo N9 -21.160 -18.400 13.04 %
reactivo de los generadores instalados. N10 -3.950 -9.038 0%
N11 -0.018 -0.018 0%
N12 -1.753 -1.753 0%
N13 -5.217 -5.217 0%
N14 -1.749 -1.749 0%

En la figura 7 se muestra que al instalar

dispositivos FACTS en el SEP el voltaje
tiende a 1 (p.u), como se puede observar
las barras 10 y 11 presentan una mejora
de estabilidad de voltaje al momento de
compensar los nodos del SEP con
Figura 6. Modelamiento de la potencia reactiva con
SVC en el modelo a prueba de 14 barras propuesto dispositivos SVC.
por la IEEE

Los valores observados en la tabla 3 son

producto del resultado del flujo óptimo
de potencia, después de haber colocado
un dispositivo SVC en las barras 4, 5, 7,
11 y 13 de modelo de 14 bus de la IEEE.
Se observa que al colocar un SVC en las
barras mencionadas la cargabilidad de
potencia reactiva en los nodos del
sistema se distribuye de mejor forma.
Figura 7. Comportamiento del perfil de tensión con
SVC en el modelo de 14 barras de la IEEE

Los valores observados en la tabla 4 son

producto del resultado del flujo óptimo
de potencia, después de haber colocado
un dispositivo SVC en las barras 4, 5, 7,
11 y 13 del modelo de 14 bus de la
IEEE. Se observa que al colocar un SVC
en las barras mencionadas el perfil de

14
tensión en cada uno de los nodos del Modelo de la IEEE de 30 barras.
sistema tiende a 1 (p.u).
Con el objetivo de simular se lo ha
Tabla 4. Porcentaje de mejora del perfil de tensión realizado en el modelo establecido de la
con la ubicación de SVC en el flujo óptimo de
potencia sistema de 14 barras GAMS. IEEE de 30 barras que es un sistema
normalizado que tiene las características
de contar con 6 generadores síncronos, 7
PERFIL PERFIL transformadores, 30 barras, por último se
DE DE
% encuentra a la barra 1 como barra slack
Nodos TENSIÓN TENSIÓN
MEJORA usando como potencia base 100 MVA,
[P.U] [P.U]
SIN SVC CON SVC
para la colocación optima de FACTS en
N1 0.950 0.950 0% un SEP, se realizará el respectivo análisis
N2 0.950 0.950 0% en el modelo de 30 barras de la IEEE el
N3 0.950 0.950 0% cual se muestra en la figura 9.
N4 0.950 0.950 0%
N5 0.950 0.994 4.42 %
N6 0.950 0.950 0%
N7 0.950 0.950 0%
N8 0.950 0.978 2.86 %
N9 0.950 0.950 0%
N10 1.100 0.950 13.63 %
N11 1.100 0.950 13.63 %
N12 0.950 0.950 0%
N13 0.950 0.950 0%
N14 0.950 0.950 0%

En la figura 8 se muestra que al instalar

dispositivos SVC en el sistema el
ángulo de la tensión se mejora.

Figura 9. Modelamiento de 30 barras propuesto por

la IEEE
Figura 8. Comportamiento del Angulo de Tensión
con SVC en el modelo de 14 barras de la IEEE
El flujo óptimo de potencia es calculado
mediante el optimizador GAMS, Para lo
cual se utilizan los datos presentados de
potencia reactiva. El parámetro
necesario para realizar la optimización
es la potencia reactiva por cada nodo, la
misma que determina el costo de la
instalación de los dispositivos FACTS.

15
Los valores observados en la tabla 5 son Ubicación Óptima de SVC en el
producto del resultado del flujo óptimo modelo de la IEEE de 30 barras.
de potencia. De la potencia reactiva
necesaria para realizar la optimización. En la figura 10 se muestra la ubicación
optima resultante de la metaheurística
Tabla 5. Potencia reactiva flujo óptimo de potencia
sistema de 30 barras GAMS.
Búsqueda Armónica, se observa la
ubicación de los SVC en color rojo, en las
POTENCIA barras 3, 6, 7, 9, 12, 16, 19, 21, 23, 26 y
Nodos
REACTIVA [P.U] 29 del modelo de 30 bus de la IEEE.
N1 8.009
N2 -0.327
N3 5.279
N4 -0.016
N5 3.095
N6 -4.138
N7 -0.109
N8 -0.450
N9 2.220
N10 -14.811
N11 1.885
N12 1.824
N13 0.749
N14 19.533
N15 17.198
N16 -0.018
N17 -0.058
N18 15.119
N19 1.138
N20 1.261
N21 -2.234
N22 4.441
N23 -19.464
N24 -18.578
N25 -2.212
N26 -0.023
N27 -8.700
N28 -7.805
N29 -7.352
N30 -10.055

Figura 10. Ubicación de dispositivos SVC en el

El modelo de la IEEE de 30 barras posee modelo de 30 barras de la IEEE
30 variables binarias una para cada
barra del SEP y 31 restricciones. Para validar la heurística propuesta se
Para validar los resultados de la compara el resultado obtenido en
heurística propuesta se compara el GAMS con el que se obtuvo en Matlab.
resultado de la ubicación óptima de La ubicación que determino GAMS de
FACTS obtenida en Matlab con el los dispositivos FACTS en cada uno de
resultado de la ubicación óptima de los nodos del modelo a prueba de 30
FACTS obtenida en GAMS. barras propuesto por la IEEE es igual a
la ubicación que determino Matlab.

La capacidad de los compensadores que

se deberán instalar en las barras 3, 6, 7,

16
9, 12, 16, 19, 21, 23, 26, y 29 del modelo Tabla 7. Porcentaje de mejora de potencia reactiva
con la ubicación de SVC en el flujo óptimo de
de 30 bus de la IEEE. potencia sistema de 30 barras GAMS.

Los datos presentados en la tabla 6 POTENCIA POTENCIA

hacen referencia a la capacidad de los REACTIVA REACTIVA %
Nodos
compensadores SVC en las barras 3, 6, [P.U] [P.U] MEJORA
SIN SVC CON SVC
7, 9, 12, 16, 19, 21, 23, 26, 29 de
N1 8.009 7.997 0.14 %
modelo de 30 bus de la IEEE. N2 -0.327 -0.327 0%
N3 5.279 5.343 0%
Tabla 6. Capacidad de los compensadores SVC en el
flujo óptimo de potencia sistema de 30 barras GAMS.
N4 -0.016 -0.016 0%
N5 3.095 3.317 0%
POTENCIA REACTIVA N6 -4.138 -3.589 13.26 %
Nodos N7 -0.109 -0.109 0%
[P.U] MVAR
N3 5.279 N8 -0.450 -0.450 0%
N6 -4.138 N9 2.220 -3.89 75.22 %
N7 -0.109 N10 -14.811 -14.811 0%
N9 2.220 N11 1.885 1.591 15.59 %
N12 1.824 N12 1.824 -1.824 0%
N16 -0.018 N13 0.749 0.766 0%
N19 1.138 N14 19.533 18.055 7.56 %
N21 -2.234 N15 17.198 16.068 6.57 %
N23 -19.464 N16 -0.018 -0.018 0%
N26 -0.023 N17 -0.058 -1.675 0%
N29 -7.352 N18 15.119 11.018 27.12 %
N19 1.138 6.665 0%
N20 1.261 3.353 0%
En la figura 11 se muestra que al tener N21 -2.234 -2.234 0%
dispositivos SVC instalados en el SEP, N22 4.441 1.77 60.14 %
la potencia reactiva en las barras de cada N23 -19.464 -19.464 0%
uno de los nodos se distribuye de mejor N24 -18.578 -18.671 0%
N25 -2.212 -2.256 0%
forma, aliviando así el flujo reactivo de N26 -0.023 -0.023 0%
los generadores instalados. N27 -8.700 -8.700 0%
N28 -7.805 -8.798 0%
N29 -7.352 -7.352 0%
N30 -10.055 -10.055 0%

En la figura 12 se muestra que al instalar

dispositivos FACTS en el SEP el voltaje
tiende a 1 (p.u).

Figura 11. Comportamiento de la potencia reactiva

con SVC en el modelo de 30 barras de la IEEE

Los valores observados en la tabla 7 son

producto del resultados de flujo óptimo
de potencia, después de haber colocado Figura 12. Comportamiento del perfil de tensión
un dispositivo SVC en las barras 3, 6, 7, con SVC en el modelo de 30 barras de la IEEE
9, 12, 16, 19, 21, 23, 26 y 29 del modelo
de 30 bus de la IEEE.

17
Los valores observados en la tabla 8 son En la figura 13 se muestra que al instalar
producto del resultado del flujo óptimo dispositivos SVC en el sistema el
de potencia, después de haber colocado ángulo de la tensión se mejora.
un dispositivo SVC en las barras 3, 6, 7,
9, 12, 16, 19, 21, 23, 26 y 29 del modelo
de 30 bus de la IEEE. Se observa que al
colocar un SVC en las barras
mencionadas el perfil de tensión en cada
uno de los nodos del sistema tiende a 1
(p.u).

Tabla 8. Porcentaje de mejora del perfil de tensión

con la ubicación de SVC en el flujo óptimo de
potencia sistema de 30 barras GAMS.
Figura 13. Comportamiento del Angulo de Tensión
PERFIL DE PERFIL DE con SVC en el modelo de 30 barras de la IEEE
TENSIÓN TENSIÓN %
Nodos
[P.U] [P.U] MEJORA
7. Conclusiones
SIN SVC CON SVC
N1 8.009 7.997 0.14 %
N2 -0.327 -0.327 0%
En la presente investigación se presenta
N3 5.279 5.343 0% un procedimiento el cual mejora la
N4 -0.016 -0.016 0% compensación optima de potencia
N5 3.095 3.317 0% reactiva mediante la colocación óptima
N6 -4.138 -3.589 13.26 % de FACTS en un SEP minimizando
N7 -0.109 -0.109 0% perdidas de potencia en el sistema y
N8 -0.450 -0.450 0% optimizando los costos, con la ayuda del
N9 2.220 -3.89 75.22 %
optimizador GAMS y la metaheurística
N10 -14.811 -14.811 0%
N11 1.885 1.591 15.59 %
de Búsqueda Armónica tomando en
N12 1.824 -1.824 0% cuenta los límites de voltaje y de potencia
N13 0.749 0.766 0% reactiva.
N14 19.533 18.055 7.56 %
N15 17.198 16.068 6.57 % Se han establecido dos escenarios para el
N16 -0.018 -0.018 0% análisis del comportamiento antes y
N17 -0.058 -1.675 0% después de la compensación reactiva
N18 15.119 11.018 27.12 %
mediante flujos óptimos de potencia AC
N19 1.138 6.665 0%
N20 1.261 3.353 0%
y la colocación optima de los FACTS en
N21 -2.234 -2.234 0% el modelo a prueba de 14 barras
N22 4.441 1.77 60.14 % propuesto por la IEEE, se han colocado 5
N23 -19.464 -19.464 0% dispositivos FACTS en las barras
N24 -18.578 -18.671 0% 4,5,7,11 y 13 de diferentes capacidades
N25 -2.212 -2.256 0% logrando con esto tener una mejora en la
N26 -0.023 -0.023 0% compensación de potencia reactiva en
N27 -8.700 -8.700 0%
todo el SEP, mejorando el perfil de
N28 -7.805 -8.798 0%
N29 -7.352 -7.352 0%
voltaje de todas las barras del sistema
N30 -10.055 -10.055 0% pudiendo tener una estabilidad de voltaje
entre 0.95 a 1.05 p.u, por último se tiene
al sistema compensando con 5
dispositivos FACTS logrando tener las
menores perdidas de potencia con
respecto al escenario inicial, ubicando los
dispositivos SVC en los nodos óptimos

18
para lograr una compensación óptima de
potencia reactiva con mayor beneficio al
sistema eléctrico y entregando mayor
confiabilidad para el usuario.

Por medio de la metaheurística de

Búsqueda Armónica es posible realizar
una óptima ubicación de los FACTS en
un SEP, minimizando la cantidad de
FACTS instalados en la red y por ende
disminuyendo el precio de operación,
instalación y mantenimiento de los
mismos. Las simulaciones muestran que
la colocación óptima de FACTS en el
SEP ayuda a mejorar la componente
reactiva, el perfil de tensión, y el ángulo
de tensión de algunas barras del sistema
eléctrico de potencia acercándose al
1(p.u).

Los resultados obtenidos de la

modelación nos permiten redistribuir de
mejor forma los flujos reactivos en las
líneas de transmisión del sistema
eléctrico de potencia, ayudando así a
aliviar los reactivos en los generadores de
las centrales de generación eléctrica.

8. Trabajos Futuros
- Dimensionamiento de los tipos de
FACTS para mejorar la compensación
reactiva en las líneas de transmisión.

- Selección del tipo de FACTS que se

debe usar para disminuir pérdidas de
energía en la red del sistema eléctrico de
transmisión.

- Minimizar los costos de instalación de

FACTS utilizando dispositivos UPFC y
TCSC en las líneas de transmisión.

19
 9. Referencias [8] Gyugyi M.L.,Shauder
C.D.,SenK.K.,Static
[1] Saravanan M.,Slochanal Synchronous Series
S.M.R.,Venkatesh P.,Abraham Compensator A solid state
J.P.S., Application of particle Approach to the series
swarm optimization technique compensation of transmission
for optimal location of FACTS line, IEEE Transactions on
devices considering cost of Power Delivery,12(1997),No.
installation and system 3,406- 417.
loadability, Electric Power [9] Thukaram D.,Lomi A., Selection
System Research,77 (2007),276– of static VAR compensator
283. location and size for system
[2] Hingorani N.G., High Power voltage stability improvement,
Electronics and Flexible AC Electr Power Syst
Transmission System, IEEE Res,54(2000):139–150.
Power Engineering Review ,7 [10] Mithulanantham N.,Cañizares
(1988),3-4. C.A.,Reeve J.,Comparison of
[3] Galiana G.D. , Assessment and PSS, SVC and STATCOM
Control of the Impact of FACTS controllers for damping power
devices on Power System system oscillation, IEEE Trans
Performance, IEEE Transactions Power Syst, 18 (2003),786–792.
on Power Systems, 11 (1996), [11] Sode-Yome A.,Mithulananthan
No. 4,1931- 1936. N.,Comparison of shunt
[4] Sirjani R., Mohamed A., Shareef capacitor, SVC and STATCOM
H., Optimal placement and sizing in static voltage stability margin
of Static Var Compensators in enhancement.Int J Elec Eng Educ
power systems using Improved , 41(2004),No. 2,158–171.
Harmony Search Algorithm, [12] Khajehzadeh M., Taha M.R.,El-
PRZEGLĄD Shafie A. ,Eslami M., A Survey
ELEKTROTECHNICZNY,87 on Meta-Heuristic Global
(2011),No. 7 ,214-218. Optimization Algorithms,
[5] Azadani E.N. , Hosseinian S.H. , Research Journal of Applied
Janati M. , Hasanpor P.,Optimal Sciences, Engineering and
Placement of Multiple Technology,3 (2011), No. 6,569-
STATCOM, 12th International 578.
Middle-East Power System [13] Rechenberg I., Cybernetic
Conference,(2008),523-528. solution path of an experimental
[6] Zhang, W., Li F. , Tolbert L. M. , problem, Royal Aircraft
Optimal Allocation of Shunt Establishment, brary Translation
Dynamic Var Source SVC and ,(1965), 1122.
STATCOM: A Survey, [14] Holland J., An introduction with
Department of Electrical and application to biology, control
Computer Engineering, the and artificial intelligence
University of Tennessee, Adaptation in Natural and
Knoxville, (2006),1-6. Artificial System. MIT Press,
[7] Saranjeet K., Evolutionary Cambridge, MA, (1975).
Algorithm Assisted Optimal [15] Holland J.H., Adaptation in
Placement of FACTS Controllers natural and artificial systems.
in Power system, Master Thesis, The University of Michigan
Thapar University,(2009). Press,AnnArbor, Michigan,

20
 (1975). Evolutionary Strategies,
[16] Kirkpatrick S.,Gelatt Proceedings of IEEE
C.D.,Vecchi M.P., Optimization international conference on
by simulated annealing, Transmission and Distribution
Science,220 (1983), No. Evolutionary
4598,671-680. Computation,(2006),1-6.
[17] Farmer J.D., Packard N.H. [25] Dominguez-Navarro
,Perelson A.S. , The immune J.A.,Bernal-Agustin J.L., Diaz
system, adaptation, and machine A.,Requena D.,Vargas E.P.,
learning. Physica D., Optimal parameters of FACTS
22(1986),No. 1-3, 187-20413. devices in electric power systems
[18] Dorigo M. , Maniezzo V. , applying evolutionary
Optimization, Learning and strategies,Electrical Power and
Natural Algorithms. Ph.D. Energy Systems,29 (2007),83–
Thesis,Politecnico di Milano, 90.
Italy,(1992). [26] Talebi M.A.,Kazemi A.,Gholami
[19] Kennedy J. , Eberhart R. , A.,Rajabi M.,Ootimal Placement
Particle swarm optimization. of Static Var Compensators in
IEEE International Conference Distribution Feeders For Load
on Neural Networks Perth, Balancing By Genetic
Australia, (1995) ,1942- 1948. Algorithm, 18th International
[20] Geem Z.W., Kim J.H., Conference and Exhibition on
Loganathan G. ,A new heuristic Electricity Distribution,(2005),1-
optimization algorithm: 6.
Harmony search, Simulation, 76 [27] Zhang Y.,Milanovic
(2001), No.2, 136-701. J.V.,Voltage Sag Cost Reduction
[21] Nakrani S. , Tovey C. , On honey With Optimally Placed FACTS
bees and dynamic server Devices, 9th Internatinal
allocation in internet hosting Conference on Electrical
centers.Adapt. Behav., 12(2004), Power,Quality and Utilisation,
No. 3-4, 223. (2007),1-6.
[22] Rezaee-Jordehi A. ,Joorabian M. [28] Farsabgi M.M., Nezamabadi-
, Optimal Placement of pour H., Song Y.H., Lee
MultiType FACTS Devices in K.Y.,Placement of SVCs and
Power Systems Using Evolution Selection of Stabilizing Signals
Strategies,The 5th International in Power Systems,IEEE
Power Engineering and TRANSACTIONS ON POWER
Optimization SYSTEMS,22 (2007), NO.3
Conference,(2011),352-357. ,1061-1071.
[23] Mendes J.C.,Saavedra [29] Baghaee H.R., Jannati M.,
O.R.,Pessanha J.O., Power Vahidi B.,Hosseinian
system restoration restoration S.H.,Rastegar H., Improvement
with priority loads using an of Voltage Stability and Reduce
evolutionary strategy, Power System Losses by
Proceedings of 34th North Optimal GA-based Allocation of
American Power Symposium, Multi-type FACTS Devices. 11th
Arizona,(2002),254- 260. International Conference on
[24] Santiago M., Maldonado R., Optimization of Electrical and
Optimal Placement of FACTS Electronic
Controllers in Power systems via Equipment,(2008),209214.

21
[30] Huang J.S. ,Negnevitsky M. , A Management via Optimal
Messy Genetic Algorithm Based Location of FACTS Devices in
Optimization Scheme for SVC Deregulated Power Systems,
Placement of Power Systems Third International Conference
under Critical Operation on Electric Utility Deregulation
Contingence, International and Restructuring and Power
Conference on Computer Technologies,(2008),1592-1957.
Science and Software [37] Chang C.S.,Huang J.S., Optimal
Engineering,(2008),467-472. SVC placement for voltage
[31] Pisica ,I., Bulac C.,Toma L. , stability reinforcement, Electric
Eremia M., Optimal SVC Power Systems Research, 42
Placement in Electric Power (1997),165 172.
Systems Using a Genetic [38] Pirayesh ,A.,Parastar A.,
Algorithms Based Method , Nikoukar J.,Optimal location of
IEEE Bucharest Power Tech FACTS devices in a power
Conference (2009),1-6. system using modified particle
[32] Sundareswaran K.,Bharathram swarm optimization, 42nd
P.,Siddharth M.,Vaishnavi International Universities Power
G.,Shrivastava N.A.,Sarma Engineering
H.,Voltage Profile enhancement Conference,(2007),1122-1128.
through optimal placement of [39] Carrión D., González W., López
FACTS devices using G., Optimal fault location in
QueenBee-Assisted GA,Third transmission lines using hybrid
International Conference on method, Innovative Smart Grid
Power Systems,(2009),1-5. Technologies Conference - Latin
[33] Berizzi A., Bovo C. , Ilea V., America (ISGT Latin America),
Optimal placement of FACTS to 2017 IEEE PES.
mitigate congestions and inter- [40] Carrión D., Inga E., Gonzalez
area oscillations,IEEE W., Huncapié R., Optimal
Trondheim PowerTech, 2011,1- Geographical Placement of
8. Phasor Measurement Units based
[34] Tavakoli Bina M.,Rezaei on Clustering Techniques, Power
Siahbidi J.,Kanzi K.,Application Engineering Conference
of Averaging Technique to the (UPEC), 2016.
Power System Optimum [41] Caisapanta W., Carrión D.,
Placement and Sizing of Static Expansión de sistemas de
Compensators, The 7th transmisión electrica usando
International Power Engineering criterios de óptima potencia AC,
Conference,(2005),1-6. Universidad Politécnica
[35] Samimi A.,Golkar M.A. , A Salesiana, 2016..
Novel Method for Optimal [42] Sanchez A., Carrión D.,
Placement of STATCOM in Modeling of the Behavior Power
Distribution Networks Using Flow on Transmission Lines
Sensitivity Analysis by Based on Voltage and Current
DIgSILENT Software, Asia- Synchronopasors, IEEE Latin
Pacific Power and Energy America Transsactions, 2018.
Engineering Conference [43] Sánchez A., Carrión D.,
,(2011),1-5. Modelación del comportameinto
[36] Gitizadeh M.,Kalantar M., A de flujo de potencia en sistemas
New Approach for Congestion eléctricos de transmisión

22
 mediante sincrofasores e tensión
y corriente. Universidad
Politécnica Salesiana, 2018.
[44] Escudero P., Carrión D., Modelo
de expansión de un sistema de
transmisión basado en
linealización de flujos óptimos
de potencia AC, Universidad
Politécnica Salesiana, 2018.
[45] Pérez D., Carrión D., Ubicación
óptima de dispositivos FACTS
en redes electricas de transmisión
mediante búsqueda inteligente,
Universidad Politécncia
Salesiana, 2017.
[46] Yánez C., Carrión D.,
Caracterización de la carga
eléctrica por medio de fasores de
voltaje y corriente, Universidad
Politécnica Salesiana, 2017.
[47] Pillajo J., Carrión D., Estabilidad
de sistemas eléctricos de potencia
mediante técnicas de
conmutación en sistemas de
transmisión, Universidad
Politécnica Salesiana, 2017.
[48] Ayala C., Carrión D.,
Modelación de la propagación de
ármonicos en sistemas de
transmisión eléctrica,
Universidad Politécnica
Salesiana, 2017.
[49] Yaguana G., Carrión D.,
Optimización de la expansión de
los sistemas de transmisión
usando Gams, Universidad
Politécncia Salesiana, 2016.
[50] Jacome D., Aguila A.,
Localizacion y
mimensionamiento optimo de
dispositivos SVC para la
compensacion de potenca
reactiva en SEP, Universidad
Politécncia Salesiana, 2017.

23
9.1 Anexos Datos del modelo de prueba de la
Datos del modelo de prueba de la IEEE IEEE DE 30 barras.
DE 14 barras.
Tabla 12. Datos Líneas de Transmisión Modelo de
Tabla 9. Datos Líneas de Transmisión Modelo de 30 bus IEEE
14 bus IEEE.
Linea DESDE HASTA IMPEDA IMPEDA
Linea DESDE HASTA IMPEDA IMPEDA No BARRA BARRA R X
No BARRA BARRA R X 1 1 2 0.0192 0.0575
1 1 2 0.01938 0.05917 2 1 3 0.0452 0.1852
2 1 5 0.05403 0.22304 3 2 4 0.057 0.1737
3 2 3 0.04699 0.19797 4 3 4 0.0132 0.0379
4 2 4 0.05811 0.17632 5 2 5 0.0472 0.1983
5 2 5 0.05695 0.17388 6 2 6 0.0581 0.1763
6 3 4 0.06701 0.17103 7 4 6 0.0119 0.0414
7 4 5 0.01335 0.04211 8 5 7 0.046 0.116
8 4 7 0.00000 0.20912 9 6 7 0.0267 0.082
9 4 9 0.00000 0.55618 10 6 6 0.012 0.042
10 5 6 0.00000 0.25202 11 6 9 0 0.208
11 6 11 0.09498 0.19890 12 6 10 0 0.556
12 6 12 0.12291 0.25581 13 9 11 0 0.208
13 6 13 0.06615 0.13027 14 9 10 0 0.11
14 7 8 0.00000 0.17615 15 4 12 0 0.256
15 7 9 0.00000 0.11001 16 12 13 0 0.14
16 9 10 0.03181 0.0845 17 12 14 0.1231 0.2559
17 9 14 0.12711 0.27038 18 12 15 0.0662 0.1304
18 10 11 0.08205 0.19207 19 12 16 0.0945 0.1987
19 12 13 0.22092 0.19988 20 14 15 0.221 0.1997
20 13 14 0.17093 0.34802 21 16 17 0.0824 0.1932
22 15 18 0.107 0.2185
Tabla 10. Datos de los Límites de Generación 23 18 19 0.0639 0.1292
Modelo de 14 bus IEEE. 24 19 20 0.034 0.068
25 10 20 0.0936 0.209
Gen P. MIN P. MAX Q.MIN Q.MAX 26 10 17 0.0324 0.0845
No (MW) (MW) (MVAR) (MVAR) 27 10 21 0.0348 0.0749
G1 10 262 -40 75 28 10 22 0.0727 0.1499
G2 0 80 -40 50 29 21 22 0.0116 0.0236
G3 0 0 0 40 30 15 23 0.1 0.202
G4 0 0 -6 24 31 22 24 0.115 0.179
G5 0 0 -6 24 32 23 24 0.132 0.27
33 24 25 0.1885 0.3292
Tabla 11. Datos de Carga en Nodos Modelo de 14 34 25 26 0.2544 0.38
bus IEEE
35 25 27 0.1093 0.2087
36 28 27 0 0.369
Linea P Q
37 27 29 0.2198 0.4153
No (MW) (MVAR)
38 27 30 0.3202 0.6027
1 0.0000 0.0000
39 29 30 0.2399 0.4533
2 0.2170 0.1270
40 8 28 0.0636 0.2
3 0.9420 0.1900
41 6 28 0.0169 0.0599
4 0.4780 -0.0390
5 0.0760 0.0160
6 0.1120 0.0750
7 0.0000 0.0000
8 0.0000 0.0000
9 0.2950 -0.0240
10 0.0900 0.0580
11 0.0350 0.0180
12 0.0610 0.0160
13 0.1380 0.0580 24
14 0.1490 0.0500
 Tabla 13. Datos de los Límites de Generación
Modelo de 30 bus IEEE.

Gen P. MIN P. MAX Q.MIN Q.MAX

No (MW) (MW) (MVAR) (MVAR)
G1 5 200 -30 70
G2 20 80 -20 100
G3 15 50 -15 880
G4 10 35 -15 60
G5 10 30 -10 50
G6 12 40 -15 60

Tabla 14. Datos de Carga en Nodos Modelo de 30

bus IEEE.

Linea P Q
No (MW) (MVAR)
1 0.00 0.00
2 21.70 12.70
3 2.40 1.20
4 7.60 1.60
5 94.20 19.00
6 0.00 0.00
7 22.80 10.90
8 30.00 30.00
9 0.00 0.00
10 5.80 -17.00
11 0.00 0.00
12 11.20 7.0000
13 0.1380 7.50
14 6.20 1.60
15 8.20 2.50
16 3.50 1.80
17 9.00 5.80
18 3.20 9.00
19 9.50 3.40
20 2.20 0.70
21 17.50 11.20
22 0.00 0.00
23 3.20 1.60
24 8.70 2.70
25 0.00 0.00
26 3.50 2.30
27 0.00 0.00
28 0.00 0.00
29 2.40 0.90
30 10.60 1.90

25
10 Estado del Arte

PROPUESTAS
RESTRICCIONES
FORMULACIÓN DEL PROBLEMA PARA SOLUCIÓN
DATOS TEMÁTICA DEL
FUNCIONES OBJETIVO RESOLVER EL PROPUESTA
PROBLEMA
PROBLEMA

LOCALIZACION OPTIMA DE

LIMITACIONES DINAMICAS
EN SISTEMAS ELÉCTRICOS

MINIMIZACIÓN DE COSTOS

REDUCCION DE PERDIDAS

RENDIMIENTO DEL FLUJO

FIABILIDAD DEL SISTEMA
CALIDAD DE LA ENERGIA
FACTS - (POWER SYSTEM)
COMPENSACIÓN ÓPTIMA
DE POTENCIA REACTIVA

CONTROL DE POTENCIA
CONTROL DE POTENCIA
CAPACIDAD DE CARGA
MEDIANTE BÚSQUEDA

PERFIL DE VOLTAJE
GENETIC ALGORITM
TRANSFERENCIA DE
ACTIVO Y REACTIVO

OPTIMIZAR COSTOS
METAHEURÍSTICAS
HARMONY SEARCH

PARTICLE SWARM
DE TRANSMISIÓN

CONFIABILIDAD
OPTIMIZACION

OPTIMIZACION
OPTIMIZATION

RENDIMIENTO
FLEXIBILIDAD
OPTMIZACION

HEURÍSTICA o
ALGORITMOS,
DE POTENICA
ESTABILIDAD

ESTABILIDAD

ESTABILIDAD
SEGURIDAD

SEGURIDAD
CAPACIDAD

FIABILIDAD
ELECTRICA
ARMÓNICA

POTENCIA
STATCOM

STATCOM
COSTOS

COSTO
CITAS
ITEM

AÑO

Application of particle swarm optimization

technique for optimal location of FACTS
1 2015
devices considering cost of installation and
1                              
system loadability, Electric Power System.
High Power Electronics and Flexible AC
2 2015 Transmission System, IEEE Power 1                                
Engineering Review.
Assessment and Control of the Impact of
FACTS devices on Power System
3 2015
Performance, IEEE Transactions on Power
1                                 
Systems.
Optimal placement and sizing of Static Var
4 2014 Compensators in power systems using 4                              
Improved Harmony Search Algorithm.
Optimal Placement of Multiple STATCOM,
5 2014 12th International Middle-East Power 1                                
System Conference.
Optimal Allocation of Shunt Dynamic Var
Source SVC and STATCOM: A Survey,
6 2014
Department of Electrical and Computer
1                                
Engineering, the University of Tennessee.
Evolutionary Algorithm Assisted Optimal
7 2014 Placement of FACTS Controllers in Power 2                               
system.
Synchronous Series Compensator A solid
state Approach to the series compensation of
8 2014
transmission line, IEEE Transactions on
2                                
Power Delivery.
Selection of static VAR compensator location
9 2014 and size for system voltage stability 3                             
improvement, Electr Power System.

26
 PROPUESTAS
RESTRICCIONES
FORMULACIÓN DEL PROBLEMA PARA SOLUCIÓN
DATOS TEMÁTICA  FUNCIONES OBJETIVO  DEL
RESOLVER EL PROPUESTA 
PROBLEMA  PROBLEMA 



TRANSMISIÓN
POTENCIA REACTIVA EN SISTEMAS
DE

LOCALIZACION OPTIMA DE FACTS -

HEURÍSTICA o METAHEURÍSTICAS
CONTROL DE POTENCIA ACTIVO Y


PARTICLE SWARM OPTIMIZATION






TRANSFERENCIA DE POTENCIA
MEDIANTE BÚSQUEDA ARMÓNICA


REDUCCION DE PERDIDAS DE


LIMITACIONES DINAMICAS
MINIMIZACIÓN DE COSTOS

RENDIMIENTO DEL FLUJO


GENETIC ALGORITM

FIABILIDAD DEL SISTEMA

OPTIMIZAR COSTOS
CALIDAD DE LA ENERGIA

CONTROL DE POTENCIA


CAPACIDAD DE CARGA

CONFIABILIDAD
OPTIMIZACION 
ÓPTIMA

RENDIMIENTO
PERFIL DE VOLTAJE
FLEXIBILIDAD



HARMONY SEARCH


SEGURIDAD

FIABILIDAD
CAPACIDAD
(POWER SYSTEM)






OPTIMIZACION
OPTMIZACION

ALGORITMOS,

ESTABILIDAD

ESTABILIDAD

ESTABILIDAD

SEGURIDAD
ELECTRICA

POTENICA
REACTIVO

STATCOM

STATCOM
COSTOS

COSTO
CITAS
ITEM

AÑO

DE
COMPENSACIÓN

ELÉCTRICOS

Comparison of PSS, SVC and STATCOM

10 2013 controllers for damping power system 1                               
oscillation, IEEE Trans Power Syst.
Comparison of shunt capacitor, SVC and
11 2013 STATCOM in static voltage stability margin 1                               
enhancement.Int J Elec Eng Educ.
A Survey on Meta-Heuristic Global
Optimization Algorithms, Research Journal
12 2013
of Applied Sciences, Engineering and
1                              
Technology.
Cybernetic solution path of an experimental
13 2013 problem, Royal Aircraft Establishment, 1                          
brary Translation.
An introduction with application to biology,
control and artificial intelligence Adaptation
14 2013
in Natural and Artificial System. MIT Press,
3                                  
Cambridge.
Adaptation in natural and artificial systems.
15 2013
The University of Michigan Press,AnnArbor.
3                        
16 2013 Optimization by simulated annealing. 75                              
The immune system, adaptation, and
17 2012
machine learning. Physica.
4                                
Optimization, Learning and Natural
18 2011 Algorithms. Ph.D. Thesis,Politecnico di 13                                 
Milano.
Particle swarm optimization. IEEE
19 2011 International Conference on Neural 1                              
Networks Perth, Australia.
A new heuristic optimization algorithm:
20 2010
Harmony search, Simulation.
8                                
On honey bees and dynamic server allocation
21 2010
in internet hosting centers.Adapt.
1                                
Optimal Placement of MultiType FACTS
Devices in Power Systems Using Evolution
22 2009
Strategies,The 5th International Power
3                               
Engineering and Optimization .

27
 PROPUESTAS
RESTRICCIONES
FORMULACIÓN DEL PROBLEMA PARA SOLUCIÓN
DATOS TEMÁTICA  FUNCIONES OBJETIVO  DEL
RESOLVER EL PROPUESTA 
PROBLEMA  PROBLEMA 



DE

TRANSMISIÓN
POTENCIA REACTIVA EN SISTEMAS

LOCALIZACION OPTIMA DE FACTS -

HEURÍSTICA o METAHEURÍSTICAS
CONTROL DE POTENCIA ACTIVO Y


PARTICLE SWARM OPTIMIZATION




MEDIANTE BÚSQUEDA ARMÓNICA

TRANSFERENCIA DE POTENCIA


REDUCCION DE PERDIDAS DE


LIMITACIONES DINAMICAS
MINIMIZACIÓN DE COSTOS


RENDIMIENTO DEL FLUJO

GENETIC ALGORITM

FIABILIDAD DEL SISTEMA

OPTIMIZAR COSTOS
CALIDAD DE LA ENERGIA

CONTROL DE POTENCIA


CAPACIDAD DE CARGA

CONFIABILIDAD
OPTIMIZACION 
ÓPTIMA

RENDIMIENTO
PERFIL DE VOLTAJE
FLEXIBILIDAD



HARMONY SEARCH


SEGURIDAD

FIABILIDAD
CAPACIDAD
(POWER SYSTEM)





OPTIMIZACION
OPTMIZACION

ALGORITMOS,

ESTABILIDAD

ESTABILIDAD

ESTABILIDAD

SEGURIDAD
ELECTRICA

POTENICA
REACTIVO

STATCOM

STATCOM
COSTOS

COSTO
CITAS
ITEM

AÑO

DE
COMPENSACIÓN

ELÉCTRICOS

Power system restoration restoration with

priority loads using an evolutionary strategy,
23 2009
Proceedings of 34th North American Power
31                                
Symposium.
Optimal Placement of FACTS Controllers in
Power systems via Evolutionary Strategies,
24 2009 Proceedings of IEEE international 10                             
conference on Transmission and Distribution
Evolutionary.
Optimal parameters of FACTS devices in
electric power systems applying evolutionary
25 2008
strategies,Electrical Power and Energy
1                                
Systems.
Optimal Placement of Static Var
Compensators in Distribution Feeders For
26 2008 Load Balancing By Genetic Algorithm, 18th 35                               
International Conference and Exhibition on
Electricity.
Voltage Sag Cost Reduction With Optimally
Placed FACTS Devices, 9th Internatinal
27 2008
Conference on Electrical Power,Quality and
61                              
Utilisation.
Placement of SVCs and Selection of
28 2006
Stabilizing Signals in Power Systems.
6                          
Improvement of Voltage Stability and
Reduce Power System Losses by Optimal GA-
based Allocation of Multi-type FACTS
29 2006
Devices. 11th International Conference on
6                                  
Optimization of Electrical and Electronic
Equipment
Optimization Scheme for SVC Placement of
Power Systems under Critical Operation
30 2015 Contingence, International Conference on 18                         
Computer Science and Software
Engineering.

28
 PROPUESTAS
RESTRICCIONES
FORMULACIÓN DEL PROBLEMA PARA SOLUCIÓN
DATOS TEMÁTICA  FUNCIONES OBJETIVO  DEL
RESOLVER EL PROPUESTA 
PROBLEMA  PROBLEMA 



DE

TRANSMISIÓN
POTENCIA REACTIVA EN SISTEMAS

LOCALIZACION OPTIMA DE FACTS -

HEURÍSTICA o METAHEURÍSTICAS
CONTROL DE POTENCIA ACTIVO Y


PARTICLE SWARM OPTIMIZATION




MEDIANTE BÚSQUEDA ARMÓNICA

TRANSFERENCIA DE POTENCIA


REDUCCION DE PERDIDAS DE


LIMITACIONES DINAMICAS
MINIMIZACIÓN DE COSTOS


RENDIMIENTO DEL FLUJO

GENETIC ALGORITM

FIABILIDAD DEL SISTEMA

OPTIMIZAR COSTOS
CALIDAD DE LA ENERGIA

CONTROL DE POTENCIA


CAPACIDAD DE CARGA

CONFIABILIDAD
OPTIMIZACION 
ÓPTIMA

RENDIMIENTO
PERFIL DE VOLTAJE
FLEXIBILIDAD



HARMONY SEARCH


SEGURIDAD

FIABILIDAD
CAPACIDAD
(POWER SYSTEM)





OPTIMIZACION
OPTMIZACION

ALGORITMOS,

ESTABILIDAD

ESTABILIDAD

ESTABILIDAD

SEGURIDAD
ELECTRICA

POTENICA
REACTIVO

STATCOM

STATCOM
COSTOS

COSTO
CITAS
ITEM

AÑO

DE
COMPENSACIÓN

ELÉCTRICOS

Optimal SVC Placement in Electric Power

Systems Using a Genetic Algorithms Based
31 2015
Method , IEEE Bucharest Power Tech
1                              
Conference.
Voltage Profile enhancement through
optimal placement of FACTS devices using
32 2015
QueenBee-Assisted GA,Third International
1                                
Conference on Power Systems.
Optimal placement of FACTS to mitigate
33 2014 congestions and inter-area oscillations,IEEE 1                                 
Trondheim PowerTech.
Application of Averaging Technique to the
Power System Optimum Placement and
34 2014 Sizing of Static Compensators, The 7th 4                              
International Power Engineering
Conference.
Novel Method for Optimal Placement of
STATCOM in Distribution Networks Using
35 2014 Sensitivity Analysis by DIgSILENT 1                                
Software, Asia-Pacific Power and Energy
Engineering Conference.
A New Approach for Congestion
Management via Optimal Location of FACTS
Devices in Deregulated Power Systems,
36 2014
Third International Conference on Electric
1                                
Utility Deregulation and Restructuring and
Power Technologies.
Optimal SVC placement for voltage stability
37 2014 reinforcement, Electric Power Systems 2                               
Research.

Optimal location of FACTS devices in a

power system using modified particle swarm
38 2014
optimization, 42nd International
2                                
Universities Power Engineering.

29
 PROPUESTAS
RESTRICCIONES
FORMULACIÓN DEL PROBLEMA PARA SOLUCIÓN
DATOS TEMÁTICA  FUNCIONES OBJETIVO  DEL
RESOLVER EL PROPUESTA 
PROBLEMA  PROBLEMA 



DE

TRANSMISIÓN
POTENCIA REACTIVA EN SISTEMAS

LOCALIZACION OPTIMA DE FACTS -

HEURÍSTICA o METAHEURÍSTICAS
CONTROL DE POTENCIA ACTIVO Y


PARTICLE SWARM OPTIMIZATION




MEDIANTE BÚSQUEDA ARMÓNICA

TRANSFERENCIA DE POTENCIA


REDUCCION DE PERDIDAS DE


LIMITACIONES DINAMICAS
MINIMIZACIÓN DE COSTOS

RENDIMIENTO DEL FLUJO


GENETIC ALGORITM

FIABILIDAD DEL SISTEMA

OPTIMIZAR COSTOS
CALIDAD DE LA ENERGIA

CONTROL DE POTENCIA


CAPACIDAD DE CARGA

CONFIABILIDAD
OPTIMIZACION 
ÓPTIMA

RENDIMIENTO
PERFIL DE VOLTAJE
FLEXIBILIDAD



HARMONY SEARCH


SEGURIDAD

FIABILIDAD
CAPACIDAD
(POWER SYSTEM)






OPTIMIZACION
OPTMIZACION

ALGORITMOS,

ESTABILIDAD

ESTABILIDAD

ESTABILIDAD

SEGURIDAD
ELECTRICA

POTENICA
REACTIVO

STATCOM

STATCOM
COSTOS

COSTO
CITAS
ITEM

AÑO

DE
COMPENSACIÓN

ELÉCTRICOS

Optimal fault location in transmission lines

39 2017 using hybrid method, Innovative Smart Grid 3                             
Technologies Conference - Latin America.
Optimal Geographical Placement of Phasor
40 2018 Measurement Units based on Clustering 1                               
Techniques, Power Engineering Conference.
Expansión de sistemas de transmisión
41 2017 electrica usando criterios de óptima potencia 1                               
AC, Universidad Politécnica Salesiana.
Modeling of the Behavior Power Flow on
Transmission Lines Based on Voltage and
42 2018
Current Synchronopasors, IEEE Latin
1                              
America Transsactions.
Ubicación óptima de dispositivos FACTS en
redes electricas de transmisión mediante
43 2017
búsqueda inteligente, Universidad
1                          
Politécncia Salesiana.
Modelo de expansión de un sistema de
transmisión basado en linealización de flujos
44 2018
óptimos de potencia AC, Universidad
3                                  
Politécnica Salesiana.
Modelación del comportameinto de flujo de
potencia en sistemas eléctricos de
45 2017 transmisión mediante sincrofasores e tensión 3                        
y corriente. Universidad Politécnica
Salesiana.
Caracterización de la carga eléctrica por
46 2018 medio de fasores de voltaje y corriente, 3                              
Universidad Politécnica Salesiana.
Estabilidad de sistemas eléctricos de
potencia mediante técnicas de conmutación
47 2017
en sistemas de transmisión, Universidad
4                                
Politécnica Salesiana.
Modelación de la propagación de ármonicos
48 2018 en sistemas de transmisión eléctrica, 2                                 
Universidad Politécnica Salesiana.

30
 50
49
ITEM


AÑO

2018
2017
COMPENSACIÓN ÓPTIMA DE

Politécncia Salesiana.

Politécncia Salesiana.
POTENCIA REACTIVA EN SISTEMAS
ELÉCTRICOS DE TRANSMISIÓN
MEDIANTE BÚSQUEDA ARMÓNICA
DATOS

dispositivos SVC para la compensación de

Localización y dimensionamiento óptimo de
Optimización de la expansión de los sistemas

potencia reactiva en un SEP, Universidad

de transmisión usando Gams, Universidad

2
1
CITAS

CANTIDAD:
LOCALIZACION OPTIMA DE FACTS -

43
(POWER SYSTEM) 
OPTMIZACION 
23
CALIDAD DE LA ENERGIA
19

ELECTRICA 
ESTABILIDAD
16


CAPACIDAD DE CARGA
35


 TEMÁTICA

SEGURIDAD
22


FIABILIDAD DEL SISTEMA
19


COSTOS
11


6

CAPACIDAD
15

OPTIMIZACION 

MINIMIZACIÓN DE COSTOS
21


CONTROL DE POTENCIA ACTIVO Y
29

REACTIVO 
STATCOM
7


Tabla 15. Matriz Estado del Arte
21

FLEXIBILIDAD
FUNCIONES OBJETIVO

38

REDUCCION DE PERDIDAS DE


POTENICA 
FORMULACIÓN DEL PROBLEMA

RENDIMIENTO DEL FLUJO

26


COSTO 
6

ESTABILIDAD 
26
DEL

TRANSFERENCIA DE POTENCIA 
13
PROBLEMA

LIMITACIONES DINAMICAS 
17
RESTRICCIONES

HARMONY SEARCH
13


26

GENETIC ALGORITM

PARTICLE SWARM OPTIMIZATION

7


PARA

PROBLEMA

ALGORITMOS,
PROPUESTAS

RESOLVER EL


HEURÍSTICA o METAHEURÍSTICAS 
7

STATCOM
23


OPTIMIZACION
32


19

SEGURIDAD

CONTROL DE POTENCIA
20


ESTABILIDAD
29


PERFIL DE VOLTAJE
23


SOLUCIÓN
PROPUESTA

26

RENDIMIENTO


14

OPTIMIZAR COSTOS
31
23

FIABILIDAD
16

CONFIABILIDAD
                               
                             
 TEMÁTICA
LOCALIZACION OPTIMA
DE FACTS - (POWER…
50
COSTOS 40 OPTMIZACION
30
20
FIABILIDAD DEL 10 CALIDAD DE LA
0
SISTEMA ENERGIA ELECTRICA

SEGURIDAD ESTABILIDAD

CAPACIDAD DE CARGA

Figura 14. Temática Estado del Arte

FORMULACIÓN DEL PROBLEMA

40
35

30

25
20

15

10
5

0
CAPACIDAD OPTIMIZACION MINIMIZACIÓN CONTROL DE STATCOM FLEXIBILIDAD REDUCCION DE RENDIMIENTO COSTO
DE COSTOS POTENCIA PERDIDAS DE DEL FLUJO
ACTIVO Y POTENICA
REACTIVO

Figura 15. Formulación del Problema Estado de Arte

SOLUCIÓN DADA MEDIANTE

35

30

25

20

15

10

0
SEGURIDAD CONTROL DE ESTABILIDAD PERFIL DE RENDIMIENTO OPTIMIZAR FIABILIDAD CONFIABILIDAD
POTENCIA VOLTAJE COSTOS

Figura 16. Solución Estado de Arte

32