IE-512, Análisis de Sistemas de Potencia, Universidad Nacional Autónoma de Honduras (UNAH) 1

Cálculo de Flujo óptimo de Potencia y

Aplicaciones
C.A. Canales, Estudiante IE, C.M. Romero, Estudiante IE, J.C. Mejía, Estudiante IE, L.A. Argueta,
Estudiante IE y L.D. Machuca, Estudiante IE

 La información que se obtiene de la resolución de los flujos de

Abstract – A continuación se encuentra planteado de manera potencia es la magnitud y ángulo de fase del voltaje de cada
concisa el cálculo de flujo óptimo de potencia utilizando varios nodo, y la potencia compleja que fluye en cada línea.
métodos numéricos entre los cuales están el algoritmo de flujo de
potencia óptimo y el método de ajuste de parábola, de estos Los métodos para resolver los problemas de flujos de potencia
métodos se decide cual es el adecuado para el flujo óptimo de son métodos numéricos como, Gauss-Seidel, Newton-
potencia dependiendo de las condiciones del sistema eléctrico de
Raphson, Newton-Raphson con desacoplado rápido, etc. Estos
potencia. Para cada uno de estos algoritmos se hace un análisis
profundo comparando las ventajas y desventajas de cada uno y así métodos difieren en la forma de sus algoritmos y son elegidos
poder concluir cuál de ellos es el indicado dependiendo las dependiendo los requerimientos y la capacidad de
características y condiciones que se presenten en el SEP. Además procesamiento de datos.
se toma en consideración factores económicos al momento de la
generación de potencia.
B. Flujo óptimo de potencia
Palabras Clave—Barra de generación, despacho, flujo óptimo El flujo óptimo de potencia, fue definido a comienzos de la
de potencia, gradiente, sistema eléctrico de potencia. década de los 60’s. En general, es un problema de optimización
para reducir las pérdidas en un SEP.
I. INTRODUCCIÓN
Existen dos objetivos básicos a cumplir en la operación de un
L a preparación, estructura y el estudio de la operación de los
sistemas de potencia requieren de análisis a fin de evaluar
el desempeño del sistema existente. Los estudios
SEP;

identifican y alertan las potenciales deficiencias en el sistema  Asegurar una operación segura.
de manera que sean corregibles o se puedan prevenir. Al  Encontrar un punto económico de operación.
análisis que da la solución de estado estacionario de un sistema
El punto económico de operación significa reducir los costos
de potencia bajo ciertas condiciones preestablecidas de
por uso de la energía eléctrica, incluyendo generación,
generación, carga y topología de la red se le denomina, flujo de
transmisión y distribución.
potencia.
El flujo de potencia óptimo consiste en despachar una serie
Las técnicas de técnicas de optimización de problemas de
de generadores con el fin de minimizar las pérdidas. Al
operación y planificación de sistemas eléctricos de potencia son
minimizar las pérdidas obtenemos un mejor funcionamiento de
un área en constante investigación debido a la importancia y el
manera económica, social y ambiental.
impacto económico que representan las pérdidas de potencia en
un SEP.
La solución que nos brinda el cálculo flujo de potencia
consiste en conocer los niveles de voltaje en todas las barras del
El primer método de solución propuesto fue el método de
sistema, el flujo de potencia por todos los elementos de la red y
gradiente reducido propuesto por Carpentier. Seguido, Dommel
las pérdidas. En cambio, el flujo óptimo de potencia considera
y Tinney abordaron el problema resolviendo las ecuaciones de
un problema de despacho económico para determinar los
Kuhn-Tucker usando una combinación del método del
valores de potencia en las barras de generación, sujeto a un
gradiente para un grupo conocido de variables independientes
conjunto de restricciones.
y funciones de penalización para violaciones en las
restricciones dependientes. En la actualidad existen muchos
métodos de resolución del problema de flujo de potencia
II. METODOLOGÍA
óptimo, se dividen en dos grandes categorías: los basados en
A. Flujo de potencia optimización matemática clásica y los basados en métodos
heurísticos. Dentro de la primera categoría está el método del
El flujo de potencia es una herramienta importante en los gradiente y en el segundo grupo está el método usando los
Sistemas Eléctricos de Potencia (SEP), para la operación, algoritmos evolutivos.
regulación y proyección del sistema.

Univerisdad Nacional Autónoma de Honduras

IE-512, Análisis de Sistemas de Potencia, Universidad Nacional Autónoma de Honduras (UNAH) 2

C. Algoritmo de flujo de potencia óptimo [Paso 4] Utilizar los valores de [𝜆] para calcular el
Para minimizar las pérdidas, la función objetivo utilizadas es gradiente. El gradiente mide la sensibilidad de la
la potencia entregada por la barra de referencia o nodo “slack”. función objetivo respecto a los cambios en las
En el cálculo de flujo de carga no se especifica la potencia variables de control [𝑢] sujeta a las restricciones de
generada por la barra de referencia. Esto se da debido a que esta igualdad [𝑔(𝑥, 𝑢)] = 0.
barra es la encargada de suplir las pérdidas del sistema.
Consecuentemente, al minimizar la potencia entregada por la
barra de referencia, se están minimizando a la vez las pérdidas
𝜕𝑓 𝜕𝑔 𝑇
[∇𝑓] = [ ] + [ ] [𝜆]
de potencia activas del sistema. Si las potencias netas son 𝜕𝑢 𝜕𝑢
mantenidas constantes, este procedimiento es llamada flujo
óptimo de reactivos. El problema de optimización, sujeto [Paso 5] Calcular una dirección factible [𝑑] usando (15).
solamente a las ecuaciones de flujo de carga del sistema, está Si la norma de la dirección factible es menor a una
dado por: tolerancia especificada, entonces pare.
min 𝑓(𝑥, 𝑢)
𝑠𝑎 [𝑔(𝑥, 𝑢)] = 0. (1) 𝑑 (𝑖 ) = 0 𝑠𝑖 ∆𝑓 (𝑖 ) < 0 𝑦 𝑢(𝑖 ) = 𝑢(𝑖 )𝑚𝑎𝑥
𝑑 (𝑖 ) = 0 𝑠𝑖 ∆𝑓 (𝑖 ) > 0 𝑦 𝑢(𝑖 ) = 𝑢(𝑖)𝑚𝑎𝑥
En ese caso se han incluido en la formulación las variables 𝑑(𝑖 ) = −∇𝑓(𝑖 )𝑜𝑡𝑟𝑜 𝑐𝑎𝑠𝑜 (15)
de control 𝑢. Estas son variables que pueden ser ajustadas sin
alterar las restricciones de igualdad. Típicamente estas
[Paso 6] Calcular un paso c y actualizar los valores de
variables representan las magnitudes de tensiones en las barras
[u] desplazándose a lo largo de la dirección factible
de generación, las cuales pueden ser ajustadas a un valor fijo.
El problema consiste entonces en ajustar estos parámetros de
[𝑢𝑘+1 ] = [𝑢𝑘 ] + [∆𝑢],
control de modo que se minimicen las pérdidas del sistema.
donde
Usando el método de optimización clásico de los
multiplicadores de Lagrange el mínimo de la función objetivo [∆𝑢] = −𝑐 [𝑑 ].
sujeta a las restricciones de igualdad se puede encontrar El valor de c indica el tamaño del paso que debe dar el
introduciendo una variable auxiliar 𝜆 para cada restricción y algoritmo en dirección contraria al gradiente. Si c es muy
minimizando la función lagrangeana irrestricta. pequeño asegura la convergencia, pero hará que el algoritmo
precise de muchas iteraciones. Si c es muy grande puede causar
𝑙(𝑥, 𝑢) = 𝑓(𝑥, 𝑢) + [𝜆]𝑇 [𝑔(𝑥, 𝑢]. (2) oscilaciones alrededor del mínimo. El valor óptimo de c
normalmente se encuentra utilizando un método de búsqueda
Al aplicar las condiciones de optimalidad de Karush Kuhn- unidimensional. Se proponen dos opciones diferentes: utilizar
Tucker se obtiene el siguiente sistema: un valor constante para c, y encontrar un valor óptimo para este
a partir de un método de ajuste de parábola.
𝜕𝑙 𝜕𝑓 𝜕𝑔 𝑇
[ ] = [ ] + [ ] [𝜆 ] = 0 (11)
𝜕𝑥 𝜕𝑥 𝜕𝑥 D. Método de ajuste de parábola.
𝜕𝑙 𝜕𝑓 𝜕𝑔 𝑇
[ ] = [ ] + [ ] [𝜆 ] = 0 (12) El método de búsqueda dimensional por ajuste de parábola
𝜕𝑢 𝜕𝑢 𝜕𝑢 consiste en calcular tres puntos diferentes de la función objetivo
𝜕𝑙
[ ] = [𝑔(𝑥, 𝑢)] = 0 (13) 𝑓0 , 𝑓1 𝑦 𝑓2 de modo que 𝑓0 > 𝑓1 y 𝑓2 > 𝑓1 . Estos puntos deben
𝜕𝜆 estar distanciados un valor igual a h como se muestra en la
figura (1).
Una forma de resolver este sistema de ecuaciones es
mediante el método del gradiente. Este método consiste en
partir de un punto factible en dirección contraria al gradiente y
llegar a un nuevo punto factible con un valor de función
objetivo menor. Repitiendo este tipo de movimiento,
eventualmente se llega al óptimo. El algoritmo de solución para
el método del gradiente es:

[Paso 1] Asumir los valores iniciales para las variables

de control [u].
Fig. 1. Ilustración del método de ajuste de parábola
[Paso 2] Encontrar una solución para el flujo de
potencia usando el método de Newton. El valor mínimo de la parábola que aproxima estos tres
[Paso 3] Resolver (14) para encontrar los valores de [𝜆]. puntos está dado por

𝜕𝑔 −𝑇 𝜕𝑓 ℎ(3𝑓0 − 4𝑓1 + 𝑓2 )
[𝜆 ] = − [ ] [ ] (14) 𝑐∗ = (16)
𝜕𝑥 𝜕𝑥 2(𝑓0 − 2𝑓1 + 𝑓2 )
IE-512, Análisis de Sistemas de Potencia, Universidad Nacional Autónoma de Honduras (UNAH) 3

La expresión (16) indica el valor de c que produce la mayor IV. REFERENCIAS

reducción de la función objetivo. El principal inconveniente de
esta metodología es encontrar tres valores de la función objetivo [1] Victor H. Hinojosa, "Flujo óptimo de potencia utilizando algoritmos
que sigan el patrón descrito por la figura (1). Para encontrar evolutivos programación en Digsilent", Universidad Técnica Federico
estos puntos se debe utilizar un proceso iterativo. Santa María, Valparaís, Chile, Tech, 30, 2015
[2] Jesús María López Lezama, "Flujo de potencia óptimo usando el método
del gradiente para reducción de pérdidas en sistemas de potencia",
Universidad EAFIT, Medellín, Colombia, Tech, ISSN 1794-9165, 2008
E. Planteamiento Del Flujo Óptimo De Potencia
[3] Julio Cesar Fernandez, “Modelo de flujo óptimo de potencia utilizando
técnicas de optimización”, Universidad Centroamericana, Antiguo
El estudio del flujo de potencia óptima nace por la necesidad Cuscatlán, El Salvador, 2011.
de mejorar las condiciones de operación en estado estacionario
de los sistemas eléctricos de potencia.
El método de flujo de potencia óptimo utiliza cantidades V. BIOGRAFÍA
controlables para optimizar una función objetivo y al mismo
tiempo se debe de satisfacer restricciones operativas.
Carlo Alexi Canales, nació en Minas De Oro,
Comayagua, Honduras, en 1996. Graduado de
La función objetivo contiene factores económicos, Ciencias y Técnicas en el Centro de Investigación e
medioambientales, de seguridad, entre otros y se resuelve Innovación Educativas en el año 2013. Recibió el
mediante técnicas de optimización adecuadas para cada caso curso de Básico de refrigeración en INFOP en el 2015
y cursos de redes de la academia Cisco-UNAH.
específico. Las restricciones son leyes de física las cuales rigen Actualmente es pasante de la carrera de ingeniería
los generadores, líneas de transmisión, equipos y estrategias eléctrica industrial en la Universidad Nacional
operativas. La función objetivo se expresa como una función no Autónoma de Honduras (UNAH).
lineal, mientras que las restricciones se representan como
funciones lineales y no lineales.

Se consideran varias funciones objetivos en un OPF, la más Luis Alonso Argueta, nació en Comayagua,
Comayagua, Honduras, en 1994. Graduado de
común son los costos de generación, la que refleja aspectos técnico en computación y bachiller en ciencias y
económicos del SEP. Entonces nos enfocamos en minimizar los letras en el Instituto La Inmaculada en el año 2012.
costos de generación de potencia activa. Matemáticamente Recibió el curso de instalaciones eléctricas
tenemos: residenciales en INFOP. Actualmente es estudiante
de ingeniería eléctrica industrial en la Universidad
Nacional Autónoma de Honduras (UNAH).

Carlos Mauricio Romero Reina, nació en Lajitas

Talgua Lempira, Honduras, el 01 de agosto de 1996.
Obtuvo su grado de Bachiller en ciencias y letras en
el Instituto Departamental Álvaro Contreras en Santa
Rosa De Copan. Actualmente cursa la carrera de
ingeniería eléctrica industrial en la Universidad
Nacional Autónoma de Honduras (UNAH), también
participa en el curso de nano PLC en el Instituto de
formación profesional (INFOP).

III. CONCLUSIONES
El modelo del flujo óptimo de potencia mediante el método
del gradiente es un método que se utiliza para reducir las al
momento de hacer los cálculos en los flujos de potencia para
diseñar un sistema eléctrico de potencia.
Este método tiene el inconveniente de que se debe tener al
menos tres puntos de la función la cual se puede aproximar a
una parábola.
Se hace un conjunto de pasos iterativos para poder encontrar
esos tres puntos, otro inconveniente es que la iteración siguiente
requiere un proceso computacional mayor al anterior.
IE-512, Análisis de Sistemas de Potencia, Universidad Nacional Autónoma de Honduras (UNAH) 4

José Carlos Mejía Martínez, nació en Distrito

central, Francisco Morazán, Honduras, en 1995.
Graduado de bachiller en ciencias y técnicas con
orientación en electromecánica en el Centro de
Investigación e Innovación Educativas (CIIE), en el
año 2013, recibió diversos cursos de formación en
INFOP, academia de redes CISCO-UNAH y CEETI.
Actualmente pasante de la carrera de ingeniería
eléctrica industrial en la Universidad Nacional
Autónoma de Honduras.

Luis Daniel Machuca Castillo, nació en La Paz, La

Paz, Honduras, el 11 de julio de 1986. Obtuvo su
grado de Bachiller Técnico Industrial en Electricidad
en el Instituto Polivalente Doroteo Varela Mejía.
Actualmente cursa la carrera de ingeniería eléctrica
industrial en la Universidad Nacional Autónoma de
Honduras (UNAH), y laborando en el proyecto de
Huawei MSP para Tigo en el área de Rollout B2B.
IE-512, Análisis de Sistemas de Potencia, Universidad Nacional Autónoma de Honduras (UNAH) 5