Instituto Tecnológico de La Paz.

Electromecánica.

Profesora: Abigayl Zaret Chargoy Cortez.

Sistemas eléctricos de potencia.

Problemario U3
Soluciones de flujos de potencia.
La información principalmente que se obtiene de un estudio de flujos de potencia es la
magnitud y el ángulo de fase del voltaje en cada barra y las potencias real v reactiva que
fluyen en cada línea.

El problema de flujos de potencia.

Para resolver el problema de flujos de potencia, se pueden usar las admitancias propias y
mutuas que componen la matriz de admitancias de barra 𝑌𝐵𝑎𝑟𝑟𝑎 o las impedancias de punto
de operación y de transferencia que constituyen 𝑍𝐵𝑎𝑟𝑟𝑎. Se limitará el estudio a los métodos
que usan admitancias. Las líneas de trasmisión se representan por su equivalente
monofásico nominal 𝜋. Los valores numéricos para la impedancia serie 𝑍 y la admitancia
total de carga de la línea 𝑌 (generalmente en términos de los megavars de larga de la línea
a voltaje nominal del sistema) son necesarios para cada línea, de forma que la computadora
puede determinar todos los elementos de la matriz de admitancias de barra de 𝑁 𝑋 𝑁 de la
que un típico elemento 𝑌𝑦 tiene la forma

𝑌𝑦 = |𝑌𝑦|∠𝜃𝑦 = |𝑌𝑦| cos 𝜃𝑦 + 𝑗|𝑌𝑦| sin 𝜃𝑦 = 𝐺𝑦 + 𝑗𝐵𝑦

Para avanzar en el estudio de flujos de potencia a realizar, se deben dar ciertos voltajes de
barra y se deben conocer algunos de los valores de inyecciones de potencia, como se
analizará más adelante.
El voltaje en una barra típica (𝑖) del sistema está dado en coordenadas polares por

𝑉𝑖 = |𝑉𝑖|∠𝛿𝑖 = |𝑉𝑖|(cos𝛿𝑖 + 𝑗 sin 𝛿𝑖)

y el voltaje en otra barra (𝑗) se escribe de manera similar cambiando el subíndice 𝑖 por el
𝑗.

En cada barra (𝑖) se especifican dos de las cuatro cantidades siguientes: 𝛿𝑖, |𝑉𝑖|, 𝑃𝑖 y 𝑄𝑖, y se
calculan las dos restantes. Las cantidades especificadas se seleccionan de acuerdo con el
siguiente análisis:
1. Barras de carga. En cada barra que no tiene generación, llamada barra de carga,
𝑃𝑔𝑖, y 𝑄𝑔𝑖 son cero y la potencia real 𝑃𝑑𝑖 y la reactiva 𝑄𝑑𝑖, que son tomadas del sistema
por la carga (entradas negativas al sistema) se conocen de los registros históricos,
de la planeación de cargas o de mediciones.
2. Barras de voltaje controlado. Cualquier barra del sistema en la que se mantiene
constante la magnitud del voltaje se llama de voltaje controlado. En las barras en las
que hay un generador conectado se puede controlar la generación de mega watts
por medio del ajuste de la fuente de energía mecánica y la magnitud del voltaje
puede ser controlada al ajustar la excitación del generador.
3. Barra de compensación. Por conveniencia, la barra (1) será denominada barra de
compensación. El ángulo del voltaje en la barra de compensación sirve como
referencia para los ángulos de todos los demás voltajes de barra.

El método de Gauss-Seidel.
Se desarrollarán ecuaciones para un sistema de cuatro barras y después, se escribirán las
ecuaciones generales. Se denomina la barra de compensación con el número (1), y los
cálculos empiezan con la barra(2). Si 𝑃2,𝑝𝑟𝑜𝑔 y 𝑄2,𝑝𝑟𝑜𝑔 son las potencias real y reactiva
programadas, respectivamente, que entran a la red en la barra (2), se obtiene de la
ecuación con 𝑖 igual a 2 y 𝑁 igual a 4,

Al despejar el valor de 𝑉2 se tiene

Por ahora, suponga que las barras (3) y (4) son también barras de carga con potencias real
y reactiva especificadas. En la barra (3) se tiene

La solución de la ecuación

da el voltaje corregido
𝑉2(1)calculado de la ecuación

La igualdad se alcanzará con un buen grado de exactitud después de varias iteraciones y

podría ser el valor correcto de 𝑉2 con los voltajes estimados, pero sin considerar la potencia
en las otras barras. Sin embargo, este valor podría no ser la solución para 𝑉2 en las
condiciones de flujo de potencia específicas, porque los voltajes sobre los que se basa el
cálculo de 𝑉2 son los valores estimados 𝑉2(0) y 𝑉4(0) en las otras barras y no se conocen
todavía los voltajes reales.
A medida que se encuentra el voltaje correcto en cada barra, su valor se va usando para
calcular el voltaje correcto en la siguiente barra. Por lo tanto, al sustituir en la

ecuación se obtiene, para el primer valor

calculado en la barra (3), la ecuación

El proceso se repite en la barra (4) y en cada barra de manera consecutiva a través de la

red (con la excepción de la barra de compensación) hasta completar la primera iteración en
la que se encontraron valores calculados para cada variable de estado. Entonces, se lleva
a cabo una y otra vez el proceso completo hasta que la cantidad por corregir en el voltaje
en cada barra es menor que algún índice de precisión determinado previamente. A este
proceso de solución de las ecuaciones de flujos de potencia se le conoce como el método
iterativo de Gauss-Seidel.

El método de Newton-Raphson.
La expansión en serie de Taylor para una función de dos o más variables es la base del
método de Newton-Raphson para resolver el problema de flujos de potencia. El estudio del
método se iniciará por el análisis de la solución de un problema en que intervienen
solamente dos ecuaciones y dos variables. Entonces, se verá cómo extender el análisis a
la solución de ecuaciones de flujos de potencia.

Considere la ecuación de una función ℎ1, de dos variables 𝑥1 y 𝑥2, que es igual a la constante
𝑏1 y que se expresa como

𝑔1(𝑥1,𝑥2, 𝑢) = ℎ1(𝑥1,𝑥2, 𝑢) − 𝑏1 = 0

y una segunda ecuación que contiene una función ℎ2 tal que

𝑔2(𝑥1,𝑥2,𝑢) = ℎ2(𝑥1,𝑥2, 𝑢) − 𝑏2 = 0

donde 𝑏2 es también una constante. El símbolo 𝑢 representa un control independiente que

se considera constante. Las funciones 𝑔1, y 𝑔2, se introducen por conveniencia para permitir
el análisis de las diferencias entre los valores calculados de ℎ1 y ℎ2 y sus valores
especificados respectivos 𝑏1 y 𝑏2.
Para un valor especificado de u se estimará que las soluciones de estas ecuaciones son
. Los superíndices cero indican que esos valores son estimados iniciales y no
son las soluciones reales . Se designarán las correcciones como los
valores que se tienen que sumar a para dar las soluciones correctas 𝑥1∗ y 𝑥2∗. Así,
se puede escribir

Ahora, el problema es encontrar la solución para que se hace al expandir las

ecuaciones en series de Taylor alrededor de la solución supuesta, para tener

donde las derivadas parciales de orden mayor que 1 en la serie de términos de la expansión

no han sido listadas. El término indica que la derivada parcial se evalúa para los valores
estimados . Términos como éstos se evalúan de manera similar.
Si se desprecian las derivadas parciales de orden mayor que 1, se pueden volver a escribir
las ecuaciones en forma matricial. Entonces, se tiene

donde la matriz cuadrada de derivadas parciales se llama jacobiana J o, en este caso,

(0) para indicar que se han usado los estimados iniciales para calcular los valores
numéricos de las derivadas parciales. Como se hizo anteriormente, se designará el valor
especificado de , menos el valor calculado de , como el error y se define de
manera similar el error . Se tiene entonces el siguiente sistema lineal de ecuaciones
de error

Se pueden determinar los valores de al resolver las ecuaciones de error, ya

sea por factorización triangular de la jacobiana o (para problemas muy pequeños)
invirtiendo la matriz. Sin embargo, como se truncó la expansión en serie, estos valores
añadidos a los iniciales no determinarán la solución correcta y nuevamente se hará un
intento suponiendo unos nuevos estimados , donde

Se repite el proceso hasta que la corrección es tan pequeña en magnitud que satisface el
índice de precisión seleccionado ; esto es, hasta que y sean ambas menores
que .

Estudios de flujos de potencia en el diseño y operación de sistemas.

Un estudio de flujos de potencia para un sistema que trabaja bajo condiciones de operación
reales o proyectadas como normales se llama caso base. Los resultados del caso base
constituyen un medio para la comparación de los cambios en los flujos de la red y de los
voltajes bajo condiciones anormales o de contingencia.
Un programa típico de flujos de potencia es capaz de manejar sistemas de más de
2 000 barras, 3 000 líneas y 500 transformadores. Los programas pueden, por supuesto,
ser expandidos a un tamaño aún más grande siempre que las instalaciones
computacionales disponibles sean lo suficientemente grandes.
Los datos suministrados a la computadora deben incluir los valores numéricos de los datos
de líneas y de barras y una indicación de si la barra es la de compensación, o una regulada
donde la magnitud del voltaje se mantiene constante por la generación de potencia reactiva
Q, o bien, una barra de carga donde P y Q están fijos. Donde no haya valores que se tengan
que mantener constantes, las cantidades dadas en las tablas se interpretan como
estimados iniciales. Generalmente, los límites de generación de P y Q se deben especificar,
así como también, los límites de los kilo volt amperes de las líneas.
Los mega vars totales de carga de la línea que se especifican para cada línea cuentan para
la capacitancia en paralelo o derivación y son iguales a y √3 por el voltaje de línea nominal
en kilo volts y por 𝐼𝑐𝑎𝑟𝑔𝑎, carga, divididos entre 103. Esto es,

(𝑀𝑣𝑎𝑟)𝑐𝑎𝑟𝑔𝑎 = √3|𝑉|𝐼𝑐𝑎𝑟𝑔𝑎𝑥10−3 = 𝜔𝐶𝑛|𝑉|2

donde |𝑉| es el voltaje línea a línea nominal en kilo volts, 𝐶𝑛 es la capacitancia línea a neutro
en farads para toda la longitud de la línea e 𝐼𝑐𝑎𝑟𝑔𝑎. El programa crea una representación
nominal 𝜋 de la línea, al dividir la capacitancia calculada del valor dado de los mega vars
de carga por partes iguales entre los extremos de la línea.
Los listados impresos de los resultados que provee la computadora consisten en una serie
de tablas. Por lo general, la información más importante a ser considerada es la tabla que
enlista cada número y nombre de barra, la magnitud del voltaje de barra en por unidad y su
ángulo de fase, la generación y la carga en cada barra en megawatts y megavars y, por
último, los megavars de los capacitores estáticos o reactancias en la barra. Junto con la
información de las barras está el flujo de megawatts y de megavars desde esa barra a cada
una de las líneas de trasmisión conectadas a ella. Los totales de la generación y de las
cargas del sistema se enlistan en megawatts y megavars.
En la operación de sistemas reales se registra el intercambio de potencia entre las áreas
con el fin de determinar si un área dada está produciendo la cantidad de potencia que
resultará en el intercambio deseado.

Transformadores regulantes.
Los transformadores regulantes se pueden usar para controlar los flujos de potencia real y
reactiva en un circuito.
La figura es una representación más detallada del transformador regulante de la. La
admitancia Y en por unidad es el recíproco de la impedancia en por unidad del
transformador que tiene la relación de transformación 1:t, como se muestra. La admitancia
Y se encuentra sobre el lado del transformador ideal más cercano al nodo (J), que es el
lado de cambio de derivaciones. Si se está considerando un transformador con una relación
de vueltas no nominal, la t debe ser real o imaginaria, como 1.02 para un 2% aproximado
de elevación en la magnitud del voltaje o 𝑗𝜋⁄60para un defasamiento aproximado de 3º por
fase.
En la figura se muestran las corrientes 𝐼𝑖 e 𝐼𝑗 asociadas con las dos barras y los voltajes con
respecto al nodo de referencia son 𝑉𝑖y 𝑉𝑗 Las expresiones complejas para la potencia dentro
del transformador ideal desde las barras (𝑖) y (𝑗) son, respectivamente,

Como se ha supuesto que se tiene un transformador ideal sin pérdidas, la potencia 𝑆𝑖 dentro
del transformador ideal desde la barra (𝑖) debe ser igual a la potencia −𝑆𝑗 que sale del
transformador ideal en el lado de la barra (𝑗); y así, se obtiene

𝐼𝑖 = −𝑡∗𝐼𝑗

La corriente 𝐼𝑗 se puede expresar a través de

𝐼𝑗 = (𝑉𝑗 − 𝑡𝑉𝑖)𝑌 = −𝑡𝑌𝑉𝑖 + 𝑌𝑉𝑖 Al

multiplicar por −𝑡∗ y al sustituir 𝐼𝑖 por −𝑡∗𝐼𝑗 se obtiene

𝐼𝑖 = 𝑡𝑡∗𝑌𝑉𝑖 − 𝑡∗𝑌𝑉𝑗

Al hacer 𝑡𝑡∗ = |𝑡|2 y al re arreglar las ecuaciones en la forma de 𝑌𝑏𝑎𝑟𝑟𝑎, se tiene

Si una línea de trasmisión en particular en un sistema está llevando potencia reactiva

demasiada pequeña o demasiado grande, un transformador regulante colocado en una de
las terminales de la línea puede hacer que la línea trasmita una potencia reactiva mayor o
menor. También, cualquier caída de voltaje apreciable en el lado primario de un
transformador, ocasionada por un cambio de la carga, puede que haga deseable cambiar
la selección de la derivación del transformador, siempre que se tengan derivaciones
ajustables con el fin de mantener el voltaje apropiado en la carga. Se puede investigar el
ajuste de la magnitud del voltaje en la barra por medio del cambio automático de derivación
en el programa de flujos de potencia. Por ejemplo, suponga que en el sistema de cuatro
barras se desea elevar el voltaje en la barra (3) al insertar entre la carga y la barra un
transformador regulante de la magnitud. Con el valor de t real, se selecciona 𝑖 = 𝑡 y se le
asigna el número 5 a la barra (𝑗) de la que ahora la carga se está alimentando. La
𝑌𝑏𝑎𝑟𝑟𝑎 de la red se expande en una fila y una columna para la barra (5), con el fin de
acomodar el regulador en las ecuaciones de flujos de potencia, y los elementos de las
barras (3) y (5) en la matriz anterior, se suman a los que previamente tenía la matriz de
admitancias de barras, con lo que se obtiene

Los elementos 𝑌𝑦 corresponden a los parámetros que ya estaban en la red antes de que se
añadiera el regulador. El vector de las variables de estado depende de cómo sea
considerada la barra (5) dentro del modelo de flujos de potencia. Hay dos alternativas:

La derivación t se puede considerar como un parámetro independiente con un valor

especificado antes de que se empiece la solución de los flujos de potencia.
Entonces, la barra (5) se trata como una barra de carga con un ángulo 𝛿5 y una
magnitud de voltaje 𝑉5 que se determina junto con las otras variables de estado. En
este caso, el vector de variables de estado es
𝑋 = [𝛿2, 𝛿3, 𝛿4, 𝛿5, |𝑉2|, |𝑉3|,|𝑉5|]𝑇
𝑜;
La magnitud del voltaje en la barra (5) puede ser previamente especificada.
Entonces, la derivación t reemplaza a la |𝑉5| como variable de estado que tiene que
determinarse junto con 𝛿5 en la barra (5) de voltaje controlado. En este
𝑋 = [𝛿2, 𝛿3, 𝛿4, 𝛿5, |𝑉2|,|𝑉3|,𝑡]𝑇 y la jacobiana cambia de acuerdo con esta modificación.

Método desacoplado de flujos de potencia.

Cuando se resuelven sistemas de trasmisión de potencia de gran escala, el método
desacoplado de flujos de potencia representa una alternativa para mejorar la eficiencia
computacional y reducir los requisitos de memoria. Este método hace uso de una versión
aproximada del procedimiento de Newton-Raphson. El principio sobre el que se basa el
enfoque de desacoplamiento se sustenta en dos observaciones:

• Un cambio en el ángulo de voltaje 𝛿 en una barra afecta principalmente al flujo de

potencia real P en las líneas de trasmisión y deja sin cambio, relativamente, a la
potencia reactiva 𝑄.
• Un cambio en la magnitud de voltaje |𝑉| en una barra afecta principalmente al flujo
de potencia reactiva 𝑄 en las líneas de trasmisión y deja al flujo de potencia real 𝑃,
sin cambiar, relativamente.
Se han hecho notar estos dos efectos en la sección pasada cuando se estudió el regulador
del deslizamiento de la fase y de la magnitud del voltaje. Esencialmente, la primera
observación establece que 𝜕𝑃𝑖⁄𝜕𝛿𝑗 es mucho mayor que 𝜕𝑄𝑖⁄𝜕𝛿𝑗 que por ahora se
considerará como cero. La segunda observación establece que 𝜕𝑄𝑖⁄𝜕|𝑉𝑗| es mucho mayor
que la 𝜕𝑃𝑖⁄𝜕|𝑉𝑗| que también se considerará aproximadamente cero.

La incorporación de estas aproximaciones en la jacobiana de la ecuación (9.45) hace que

los elementos de las sub matrices 𝐽12 y 𝐽21 sean cero. Entonces se tienen dos sistemas
separados de ecuaciones,

Estas ecuaciones están desacopladas en el sentido de que las correcciones del ángulo del
voltaje ∆𝛿 se calculan usando sólo los errores de la potencia real ∆𝑃, mientras las
correcciones de la magnitud del voltaje se calculan usando sólo los errores ∆𝑄.
Los dos conjuntos de ecuaciones podrían, por supuesto, resolverse alternadamente usando
en un conjunto las soluciones más recientes del otro conjunto. Pero este esquema todavía
requeriría la evaluación y factorización de las dos matrices de coeficientes en cada iteración.
Para evitar estos cálculos, se introducen más simplificaciones que se justifican a través de
la física de los flujos de potencia en líneas de trasmisión, como se explica en seguida.
En un sistema de trasmisión de potencia que está bien diseñado y apropiadamente operado:

• Las diferencias angulares (𝛿𝑖 − 𝛿𝑗) entre dos barras típicas del sistema son, por lo
general, tan pequeñas que
cos(𝛿𝑖 − 𝛿𝑗) = 1 ; sin(𝛿𝑖 − 𝛿𝑗) ≈ (𝛿𝑖 − 𝛿𝑗)
• Las susceptancias de las líneas 𝐵𝑦 son muchas veces más grandes que las
conductancias 𝐺𝑦 así que
𝐺𝑦 sin(𝛿𝑖 − 𝛿𝑗) ≪ 𝐵𝑦 cos(𝛿𝑖 − 𝛿𝑗)
• La potencia reactiva Q, que se inyecta a cualquier barra (𝑖) del sistema durante la
operación normal es mucho menor que la potencia reactiva que fluiría si todas las
líneas de la barra estuvieran en cortocircuito con la referencia. Esto es,
𝑄𝑖 ≪ |𝑉𝑖|2𝐵𝑖𝑖
Estas aproximaciones se pueden usar para simplificar los elementos de la jacobiana. Los
elementos fuera de la diagonal de 𝐽11 y 𝐽22 están dados por

Al aplicar en la ecuación pasada la identidad sin(𝛼 + 𝛽) = sin 𝛼 cos 𝛽 + cos 𝛼 sin𝛽, se obtiene

Las aproximaciones que se enlistaron anteriormente conducen a los elementos fuera de la

diagonal dados por

Se aplica la desigualdad dada por 𝑄𝑖 ≪ |𝑉𝑖|2𝐵𝑖𝑖 a esas expresiones, se llega a

Al sustituir las expresiones de las aproximaciones dadas en las ecuaciones en las matrices
de coeficientes𝐽11 y 𝐽12, se obtiene
Se multiplica la primera fila por el vector de corrección y entonces se divide la ecuación
resultante entre |𝑉2|con el fin de mostrar cómo se pueden quitar los voltajes de la matriz de
coeficientes dada por la ecuación pasada; el resultado es

También, se puede modificar la ecuación si se multiplica la primera fila por el vector de las
correcciones de ángulo y se re arregla el resultado para obtener

Al desarrollar todas las filas de la ecuación en una manera similar, se llega a dos sistemas
de ecuaciones desacoplados para las cuatro barras de la red

Generalmente, la matriz 𝐵 es simétrica y partida con elementos diferentes de cero que son
constantes y números reales, exactamente iguales al negativo de las susceptancías de
𝑌𝑏𝑎𝑟𝑟𝑎 En consecuencia, la matriz 𝐵 se forma fácilmente y sus factores triangulares
(una vez calculados al inicio del proceso de solución) no se tienen que volver a calcular,
con lo que se tienen iteraciones más rápidas. En las barras de voltaje controlado, la 𝑄 no
se especifica y ∆|𝑉| es cero; las filas y columnas que corresponden a tales barras se omiten
de la ecuación
Una estrategia de solución típica es:
1. Calcular los errores iniciales ∆𝑃⁄|𝑉|,
2. Resolver la matriz pasada para ∆𝛿,
3. Actualizar los ángulos 𝛿 y usarlos para calcular los errores ∆𝑄⁄|𝑉|,
4. Resolver la matriz pasada para ∆|𝑉|y actualizar las magnitudes |𝑉| y
5. Regresar a la matriz pasada para repetir la iteración hasta que todos los errores
estén dentro de las tolerancias especificadas.