Subestaciones Eléctricas 1

Tema 3. Flujos de Potencia

Gómez Espinoza Macías, Daniel
219295737

Resumen—El estudio de los flujos de potencia nos permite II. FLUJOS DE POTENCIA
hacer aproximaciones precisas sobre el comportamiento de Como se mencionó en la introducción, el objetivo del análisis
la red, específicamente sobre la potencia y su relación entre de los flujos de potencia es obtener la magnitud y el ángulo de
cada uno de los nodos que componen al SEP a analizar. En fase del voltaje en cada barra, así como las potencias que fluyen
este trabajo se abordará la descripción física y matemática en cada línea. Para lograr esto “se pueden usar las admitancias
del cálculo para resolver problemas de flujos de potencia, propias y mutuas que componen la matriz de admitancias de
así como la deducción y aplicación de dos métodos barra 𝑌𝑏𝑎𝑟𝑟𝑎 o las impedancias de punto de operación y de
numéricos: Newton-Raphson y Gauss-Seidel, para transferencia que constituyen 𝑍𝑏𝑎𝑟𝑟𝑎 ” [1]. En esta ocasión se
aproximar respuestas de las variables de estado del sistema. realizará el cálculo a partir de las admitancias.
Palabras clave— Flujos de Potencia, Newton-Raphson, Gauss-
El primer paso es obtener el diagrama unifilar a partir de los
Seidel, Barras.
datos introducidos a la computadora, para conocer todos los
elementos de la matriz de admitancias. Un típico elemento 𝑌𝑖𝑗
I. INTRODUCCIÓN tiene la siguiente forma:

E l análisis de los flujos de potencia consiste principalmente

en la obtención de la magnitud y ángulo de fase del voltaje
𝒀𝒊𝒋 = |𝒀𝒊𝒋 |∠𝜽𝒊𝒋 = |𝒀𝒊𝒋 |𝒄𝒐𝒔𝜽𝒊𝒋 + 𝒋|𝒀𝒊𝒋 |𝒔𝒆𝒏𝜽𝒊𝒋 = 𝑮𝒊𝒋 + 𝒋𝑩𝒊𝒋 (1)

en cada barra y las potencias real y activa que fluyen en cada El voltaje en una barra típica “i” del sistema está dado en
línea. coordenadas polares por:

Los sistemas de potencia analizados en este trabajo se 𝑽𝒊 = |𝑽𝒊 |∠𝜹𝒊 (2)

encontrarán en condiciones normales balanceadas de estado
estable trifásico, para ello se debe cumplir lo siguiente: la La corriente total que se inyecta a la red a través de la barra
generación es suficiente para abastecer la demanda más las “i” en términos de la admitancia, se representa de la siguiente
pérdidas, las magnitudes de voltaje en las barras no sufren forma:
𝑰𝒊 = 𝒀𝒊𝟏 𝑽𝟏 + 𝒀𝒊𝟐 𝑽𝟐 + ⋯ + 𝒀𝒊𝑵 𝑽𝑵 = ∑𝑵
𝒏=𝟏 𝒀𝒊𝒏 𝑽𝒏 (3)
grandes variaciones con respecto a su valor nominal, los
generadores operan dentro de los límites especificados de El voltaje y la corriente “j” se describe de la misma manera
potencia real y reactiva, y los transformadores y las líneas de que en las ecuaciones (1) y (2), pero cambiando el subíndice “i”
transmisión no se encuentran sobrecargados. por el “j”. Sean 𝑃𝑖 y 𝑄𝑖 las potencias que totales que entran a la
red a través de la barra “i”, entonces el complejo conjugado de
El estudio de los flujos de potencia es sumamente complejo, la potencia es el siguiente:
debido a que se manejan grandes cantidades de potencia y
existen pocos mecanismos que nos permitan controlar la ruta 𝑷𝒊 − 𝒋𝑸𝒊 = 𝑽𝒊 ∗ ∑𝑵
𝒏=𝟏 𝒀𝒊𝒏 𝑽𝒏 (4)
del sistema, desde el generador hasta la última carga. Muchas
veces esa ruta es incluso indirecta, variando debido a las Si se sustituyen las ecuaciones (1) y (2), así como se
impedancias de las líneas de transmisión y a los lugares donde expanden las ecuaciones para obtener la potencia real y
entra o sale la potencia. Es por esto por lo que estos estudios se reactiva, obtenemos las siguientes expresiones.
realizan con equipos computacionales robustos, que entregan 𝑷𝒊 = ∑𝑵
𝒏=𝟏|𝒀𝒊𝒏 𝑽𝒊 𝑽𝒏 | 𝐜𝐨𝐬(𝜽𝒊𝒏 + 𝜹𝒏 − 𝜹𝒊 ) (5)
un panorama más completo de la red eléctrica. Este tipo de
estudios se hace principalmente con el objetivo de conocer y 𝑸𝒊 = − ∑𝑵
𝒏=𝟏|𝒀𝒊𝒏 𝑽𝒊 𝑽𝒏 | 𝐬𝐞𝐧(𝜽𝒊𝒏 + 𝜹𝒏 − 𝜹𝒊 ) (6)
determinar las condiciones de operación de los sistemas de
Estas últimas ecuaciones son la representación de las
potencia, para monitorear la red o en su caso, planear y diseñar
ecuaciones de flujo de potencia en su forma polar. Dichas
una futura expansión del Sistema Eléctrico de Potencia
ecuaciones dan valores calculados para la potencia real 𝑃𝑖 y la
potencia reactiva 𝑄𝑖 totales que entran a la red a través de una
barra típica i.

07 febrero 2023
Subestaciones Eléctricas 2

Sea la potencia 𝑃𝑔𝑖 la potencia programada que se está generador 𝑄𝑔𝑖 para mantener el voltaje controlado no se puede
generando en la barra i y 𝑃𝑑𝑖 la potencia programada que conocer, por lo que este parámetro junto con el ángulo de
demanda la carga en esa barra. Entonces la potencia voltaje 𝛿𝑖 serán las variables desconocidas para este tipo de
programada total se puede expresar con la siguiente ecuación: barras.

∆𝑷𝒊 = 𝑷𝒊,𝒑𝒓𝒐𝒈 − 𝑷𝒊,𝒄𝒂𝒍𝒄 (7) C) Barra de compensación

Esta barra será nuestra referencia para analizar al resto de
Dicha relación se representa con la Fig 1. barras del SEP. El ángulo y la magnitud particular del voltaje
que se le asigne a la barra de compensación no son relevantes
porque lo que nos interesa son las diferencias voltaje-ángulo
que determinarán los valores calculados de 𝑃𝑖 y 𝑄𝑖 , por ello es
común utilizar las siguientes referencias 𝛿𝑖 = 0 y |𝑉1 | = 1,
Fig 1. Notación para los flujos de potencia (activa y reactiva) en una conocidas como punto de arranque. No se definen errores para
barra la barra de compensación.
Donde se nombra al valor calculado de 𝑃𝑖 como 𝑃𝑖,𝑐𝑎𝑙𝑐 y se
llega a la definición del error ∆𝑃1 como el valor programado “Las magnitudes y ángulos de los voltajes de barra que no se
programaron en los datos de entrada del estudio de flujos de
𝑃𝑖,𝑝𝑟𝑜𝑔 menos el valor calculado 𝑃𝑖,𝑐𝑎𝑙𝑐 , tal como se describen
potencia se llaman variables de estado (porque sus variables
las fórmulas (8) y (9).
describen el estado del sistema). Por lo tanto, el problema de
flujos de potencia consiste en determinar los valores para todas
∆𝑷𝒊 = 𝑷𝒊,𝒑𝒓𝒐𝒈 − 𝑷𝒊,𝒄𝒂𝒍𝒄 = (𝑷𝒈𝒊 − 𝑷𝒅𝒊 ) − 𝑷𝒊,𝒄𝒂𝒍𝒄 (8)
las variables de estado, resolviendo un número igual de
∆𝑸𝒊 = 𝑸𝒊,𝒑𝒓𝒐𝒈 − 𝑸𝒊,𝒄𝒂𝒍𝒄 = (𝑸𝒈𝒊 − 𝑸𝒅𝒊 ) − 𝑸𝒊,𝒄𝒂𝒍𝒄 (9) ecuaciones de flujos de potencia que se basan en las
especificaciones de los datos de entrada” [1].
Si los valores calculados 𝑃𝑖,𝑐𝑎𝑙𝑐 y 𝑄𝑖,𝑐𝑎𝑙𝑐 igualan
perfectamente a los valores programados 𝑃𝑖,𝑝𝑟𝑜𝑔 y 𝑄𝑖,𝑝𝑟𝑜𝑔 , se “El análisis nodal común, o de mallas, no es adecuado para
dice que los errores ∆𝑃𝑖 y ∆𝑄𝑖 son cero en la barra i y se tienen estudios de flujo de potencia porque los datos de entrada para
las siguientes ecuaciones de balance de potencia: las cargas por lo general se dan en términos de potencia, no de
impedancias. Asimismo, se considera a los generadores como
𝒈𝒊 = 𝑷𝒊 − 𝑷𝒊,𝒑𝒓𝒐𝒈 = 𝑷𝒊 − (𝑷𝒈𝒊 − 𝑷𝒅𝒊 ) = 𝟎 (10) fuentes de potencia, no fuentes de voltaje o corriente. El
problema de flujos de potencia, por lo tanto, se formula como
𝒈𝒊 = 𝑸𝒊 − 𝑸𝒊,𝒑𝒓𝒐𝒈 = 𝑸𝒊 − (𝑸𝒈𝒊 − 𝑸𝒅𝒊 ) = 𝟎 (11) un conjunto de ecuaciones algebraicas no lineales adecuado
para resolverlo por computadora” [2].
La práctica general en los estudios de flujos de potencia es la
de identificar tres tipos de barras en la red. En cada barra i se
especifican dos de las cuatro cantidades siguientes: 𝛿𝑖 , |𝑉𝑖 |, 𝑃𝑖 III. MÉTODO NEWTON-RAPHSON
y 𝑄𝑖 , y se calculan las dos restantes. Para ello se hace una
Como se mencionó anteriormente, el cálculo de los flujos de
clasificación de las barras:
potencia no se puede obtener con análisis convencionales de
circuitos, sino que son necesarios métodos numéricos para
A) Barras de carga (P-Q)
obtener aproximaciones correctas. En esta ocasión se trabajará
Son las barras en donde no existe generación. Aquí 𝑃𝑔𝑖 y 𝑄𝑔𝑖
con el método Newton-Raphson multivariable, y en este
son cero y la potencia real 𝑃𝑑𝑖 y la reactiva 𝑄𝑑𝑖 que son tomadas método se obtienen los nuevos valores mediante linealizaciones
del sistema por la caga (entradas negativas del sistema) se sucesivas de las ecuaciones involucradas, expresadas como
conocen de los registros históricos de la planeación de cargas o f(x)=0.
de las mediciones. En la práctica sólo se suele conocer la Para este método las funciones o residuos a anular son la
potencia real y la potencia reactiva se basa en un FP supuesto diferencia entre los valores de potencia calculados y los
(0.85 o mayor). A este tipo de barra se le conoce también como especificados.
barra P-Q porque los valores programados 𝑃𝑖,𝑝𝑟𝑜𝑔 = −𝑃𝑑𝑖 y
𝑄𝑖,𝑝𝑟𝑜𝑔 = −𝑄𝑑𝑖 son conocidos y los errores ∆𝑃𝑖 y ∆𝑄𝑖 pueden En cada iteración debe resolverse el siguiente sistema lineal
definirse, por lo que las dos cantidades desconocidas serán 𝛿𝑖 y de ecuaciones:
|𝑉𝑖 |.
𝒌
𝑯 𝑵 𝒌 ∆𝜽 ∆𝑷 𝒌
[ ] [ ] =[ ] (12)
B) Barras de voltaje controlado o de generación (P-V) 𝑴 𝑳 ∆𝑽/𝑽 ∆𝑸
Es cualquier barra del sistema en la que se mantiene el voltaje
constante. En las barras en las que hay un generador conectado Cuya matriz de coeficientes es el jacobiano cambiando de
se puede controlar la generación de megawatts por medio de la signo, del sistema no lineal. Una vez resuelto el sistema, las
fuente de energía mecánica y la magnitud del voltaje puede se tensiones se actualizan mediante:
contralada al ajustar la excitación del generador. Por ello en
𝜽 𝒌+𝟏 𝜽𝒌 ∆𝜽 𝒌
estas barras se conocen 𝑃𝑔𝑖 y |𝑉𝑖 |, por esta razón se le conoce [ ]
𝑽
=[ ] +[ ]
𝑽 ∆𝑽
(13)
como barra P-V. Por otro lado, la potencia reactiva del
Subestaciones Eléctricas 3

Y sus valores se muestran en la siguiente tabla: diferenciales, entre otras representaciones complejas,
matemáticamente hablando. Y aunque la gran parte del análisis
de los flujos de potencia se realiza mediante sistemas
computacionales de alta capacidad, es importante entender en
qué se basan los programas que calculan este tipo de sistemas
complejos.

REFERENCIAS
[1] J.J. Grainger y W.D. Stevenson. “Soluciones de Flujos de Potencia”, en
Análisis de Sistemas Eléctricos de Potencia, 1ra ed, McGraw-Hill, 1996, pp.
310-331
“El proceso se detiene cuando todos los componentes del
vector de residuos [Δ𝑃/Δ𝑄] son menores que el error. Debe [2] J. Duncan Glover. M. S. Sarma. T. J. Overbye, “History of Electric Power
aclararse que las particiones que se han realizado en el sistema Systems”, en Power Systems, 5ta edición, Cengeage Learning, pp 17.
de ecuaciones pretenden simplemente visualizar mejor la
[3] A. Gómez Expósito, J.L. Martínez Ramos, J.A. Rosendo Macías, E. Romero
estructura del jacobiano. En la práctica, las dos ecuaciones de Ramos, J.M. Riquelme Santos. “Flujo de Cargas”, en Sistemas Eléctricos de
cada nodo de carga, y las variables respectivas, se ordenan Potencia, Pearson Education, 2005, pp. 55-58
consecutivamente” [3].

IV. MÉTODO GAUSS-SEIDEL

Este método, a diferencia del Newton-Raphson, nos permite
trabajar con las componentes rectangulares de nuestra red, para
de esta manera converger una solución. Este método es poco
competido en la actualidad debido a que su convergencia es
bastante pobre, pero a pesar de ello aún tiene cierto interés. El
método aplica recursivamente la siguiente fórmula:

𝟏 𝑷𝒊 𝒆𝒔𝒑 −𝒋𝑸𝒊 𝒆𝒔𝒑

𝒖𝒊 𝒌+𝟏 = [ ∗ − ∑𝒊−𝟏
𝒋=𝟏 𝒚𝒊𝒋 𝒖𝒋
𝒌+𝟏
− ∑𝒏𝒋=𝒊+𝟏 𝒚𝒊𝒋 𝒖𝒋 𝒌] 𝒊 = 𝟏, 𝟐 … 𝒏 − 𝟏
𝒚𝒊𝒊 (𝒖𝒊 𝒌)
(14)

Y se detiene cuando:

𝒎𝒂𝒙𝒊 |𝒖𝒊 𝒌+𝟏 − 𝒖𝒊 𝒌 | ≤ 𝜺 (15)

Donde 𝜀 es un valor suficientemente pequeño.

Como se mencionó anteriormente, este método tiene una
convergencia pobre (el número de iteraciones es del orden n),
por lo que una mejora consiste en utilizar un factor de
aceleración 𝛼:

𝒖𝒊𝒂𝒄𝒆𝒍 𝒌+𝟏 = 𝒖𝒊 𝒌 + 𝜶(𝒖𝒊 𝒌+𝟏 − 𝒖𝒊 𝒌 ) (16)

V. CONCLUSIONES
Con estos resultados se pueden concluir varias cosas,
primeramente, se demuestra la importancia del uso de los
métodos numéricos en el cálculo de los flujos de potencia, así
como la capacidad de estos métodos de arrojar resultados con
alto grado de precisión con ecuaciones complejas.

Cuando se está trabajando con sistemas tan complejos, como

lo son los sistemas eléctricos de potencia, los cuales involucran
variables y consideraciones que quedan fuera del rango de
aplicación de los métodos comunes de análisis de circuitos, es
necesario hacer uso de métodos numéricos, o de otro tipo de
herramientas que nos permitan trabajar con matrices,
ecuaciones simultáneas de segundo o mayor grado, ecuaciones