Soluciones de Flujo de

Potencia
INTRODUCCIÓN
• Los estudios de flujos de potencia son de gran importancia en la
planeación y diseño de a expansión futura de los sistemas de
potencia, así como también en la determinación de las mejores
condiciones de operación de los sistemas existentes.
• La información más importante contenida de un estudio de flujo de
potencia es:
o Magnitud y ángulo de fase del voltaje de cada barra.
o Corriente que fluye por cada línea.
o Potencia real y reactiva que fluye por 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 𝑍𝑏𝑎𝑟𝑟𝑎 .
• Los análisis se realizarán con las admitancias.
• Las líneas de transmisión se representan por su equivalente
monofásico nominal 𝜋 como se muestra en la siguiente figura.

• La matriz de admitancias de barra 𝑁𝑥𝑁 de la que un típico elemento

𝑌𝑖𝑗 tiene la forma:
𝑌𝑖𝑗 = |𝑌𝑖𝑗 |∡𝜃𝑖𝑗 = |𝑌𝑖𝑗 | cos 𝜃𝑖𝑗 + 𝑗|𝑌𝑖𝑗 | sin 𝜃𝑖𝑗 = 𝐺𝑖𝑗 + 𝑗𝐵𝑖𝑗
Fuente: Grainger, J. y Stevenson, W., Análisis de Sistemas de Potencia
EL PROBLEMA DE FLUJOS DE POTENCIA
• Otra información esencial:
o Valores nominales de los transformadores y sus impedancias.
o Capacidades de los capacitores en derivación.
o Taps de los transformadores que pueden ser usados.
• Para poder realizar un estudio de flujo de potencia se deben conocer
ciertos voltajes de barra y los valores de inyección de potencia de los
generadores.
• El voltaje en una barra típica 𝑖 del sistema esta dado en coordenadas
polares por:
𝑉𝑖 = |𝑉𝑖 |∡𝛿𝑖 = |𝑉𝑖 |(cos 𝛿𝑖 + 𝑗 sen 𝛿𝑖 )
• La corriente total que inyecta en la red a través de la barra 𝑖 en términos
𝑌𝑖𝑛 de 𝑌𝑏𝑎𝑟𝑟𝑎 , esta dada por la sumatoria:
𝑁

𝐼𝑖 = 𝑌𝑖1 𝑉1 + 𝑌𝑖2 𝑉2 + ⋯ + 𝑌𝑖𝑁 𝑉𝑁 = ෍ 𝑌𝑖𝑛 𝑉𝑛

𝑛=1
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
EL PROBLEMA DE FLUJOS DE POTENCIA
• Sean 𝑃𝑖 y 𝑄𝑖 las potencias real y reactiva totales que entran a la red a
través de la barra 𝑖. Entonces, el complejo conjugado de la potencia
que se inyecta a la barra 𝑖 es:
Retornar
Retornar
Slide 26

• Si sustituimos las ecuaciones de voltaje y admitancia, se obtiene:

Retornar
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
 • Sea 𝑃𝑔𝑖 la potencia programada que se esta generando en la barra 𝑖 y
𝑃𝑑𝑖 la potencia programada que demanda la carga en esta barra.
Entonces, la expresión 𝑃𝑖,𝑝𝑟𝑜𝑔 = 𝑃𝑔𝑖 − 𝑃𝑑𝑖 da la potencia programada
total que esta siendo inyectada dentro de la red en la barra 𝑖.

• Se nombra al valor calculado de 𝑃𝑖 como 𝑃𝑖,𝑐𝑎𝑙𝑐 y se llega a la

definición del error ∆𝑃𝑖 como:

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia

 • Si los valores calculados 𝑃𝑖,𝑐𝑎𝑙𝑐 y 𝑄𝑖,𝑐𝑎𝑙𝑐 igualan perfectamente a los
valores programados 𝑃𝑖,𝑝𝑟𝑜𝑔 y 𝑄𝑖,𝑝𝑟𝑜𝑔 se dice que los errores ∆𝑃𝑖 y
∆𝑄𝑖 son cero en la barra 𝑖 y se tienen las siguientes ecuaciones de
balance de potencia:

• Las cuatro cantidades potencialmente desconocidas que se asocian

con cada barra 𝑖 son: 𝑃𝑖 , 𝑄𝑖 , 𝛿𝑖 𝑦 |𝑉𝑖 |. A lo más hay dos ecuaciones
como las anteriores disponibles para cada nodo y así, se debe
considerar como se puede reducir el número de cantidades
desconocidas para que se tenga el mismo número de ecuaciones
disponibles antes de empezar a resolver el problema de flujos de
potencia.
Retorno
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
• En los estudios de flujos de potencia tenemos que identificar tres
tipos de barras en la red. En cada barra 𝑖 se especifican dos de las
cuatro cantidades siguientes: 𝑃𝑖 , 𝑄𝑖 , 𝛿𝑖 𝑦 |𝑉𝑖 | y se calculan las dos
restantes.
• Tipos de barras:
o Barras de carga: En cada barra que no tiene generación, 𝑃𝑔𝑖 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. Con frecuencia en la practica solo se
conoce la potencia real y la reactiva que se basa en factor de potencia
supuesto tal como 0.85 o mayor. Los valores desconocidos son 𝛿𝑖 𝑦 |𝑉𝑖 |.
o Barras de voltaje controlado: Cualquier barra del sistema en la que se
mantiene constante la magnitud del voltaje. En las barras en las que hay un
generador conectado se puede controlar la generación de megawatts 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. En este tipo de
barra se pueden especificar apropiadamente 𝑃𝑔𝑖 y |𝑉𝑖 |.
 o 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.
En la practica el ángulo de esta barra es 0°. No se definen errores para esta
barra, y así la magnitud del voltaje se especifica como la otra cantidad
conocida. Entonces no hay necesidad de incluir las ecuaciones del slide 7 en el
problema de flujos de potencia.

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia

RESUMEN DEL PROBLEMA DE FLUJOS DE POTENCIA

• Las funciones 𝑃𝑖 y 𝑄𝑖 de las ecuaciones del slide 5 son funciones no

lineales de las variables de estado 𝛿𝑖 𝑦 |𝑉𝑖 |. Por lo tanto, en general,
los cálculos de flujos de potencia se emplean técnicas iterativas
(como los procedimientos de Gauss-Seidel y de Newton-Raphson).
❖El método de Gauss-Seidel resuelve las ecuaciones del flujo de potencia en
coordenadas rectangulares (variable compleja) hasta que las diferencias en
los voltajes de barra de una iteración a otra son lo suficientemente pequeñas.
❖El método de Newton-Raphson resuelve la forma polar de las ecuaciones de
flujos de potencia hasta que los errores ∆𝑃 y ∆𝑄 en todas las barras caen
dentro de los limites especificados.

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia

Método de Gauss-Seidel
MÉTODO DE GAUSS-SEIDEL
• El método de Gauss-Seidel es conocido como el método de los
Desplazamientos Sucesivos. Para ilustrar la técnica, consideremos la
solución de la siguiente ecuación no lineal:
𝑓(𝑥) = 0
• La ecuación anterior se re-organiza de la siguiente manera:
𝑥 = 𝑔(𝑥)
• Si x(k) es un estimado inicial de la variable x, se forma la siguiente
secuencia iterativa:
𝑥 (𝑘+1) = 𝑔(𝑥 (𝑘) )
Se obtiene una solución cuando el valor absoluto de la diferencia entre dos
iteraciones sucesivas es menor que un Margen de Tolerancia establecido:
𝑥 (𝑘+1) − 𝑥 (𝑘) ≤ Ɛ
Donde Ɛ es el nivel de exactitud deseado.

Fuente: Saadat, H., Power System Analysis

MÉTODO DE GAUSS-SEIDEL
• Las soluciones digitales de los problemas de flujos de potencia siguen
un proceso iterativo al asignar valores estimados a los voltajes de
barra desconocidos y calcular nuevos valores para cada voltaje de
barra, a partir de los estimados en las otras barras y de las potencias
real y reactiva especificadas. Así, se obtiene un nuevo conjunto de
valores, para el voltaje en cada barra, que se usa para calcular otro
conjunto de voltajes de barra. A cada cálculo de un nuevo conjunto
de voltajes se le llama iteración.
• El proceso iterativo se repite hasta que los cambios en cada barra son
menores que un valor mínimo especificado.

Fuente: Saadat, H., Power System Analysis

EJEMPLO DEL MÉTODO DE GAUSS-SEIDEL
𝑓 𝑥 = 𝑥 3 − 6𝑥 2 + 9𝑥 − 4 = 0
Para encontrar a 𝑥, la expresión se escribe como:
1 3 6 2 4
𝑥=− 𝑥 + 𝑥 + =𝑔 𝑥
9 9 9

𝑥 (0) = 2
Primera iteración:
1 6 4
𝑥 (1) =𝑔 2 =− 2 3 + 2 2 + = 2.2222
9 9 9
Segunda iteración:
1 6 4
𝑥 (2) 3
= 𝑔 2.2222 = − 2.2222 + 2.2222 + = 2.5173 2
9 9 9
La secuencia de iteraciones resulta en 2.8966, 3.7398, 3.9568, 3.9988 y 4.0000. El
proceso es repetido hasta que el cambio en variable este dentro de la precisión
deseada. Se puede ver que este método necesita muchas iteraciones para llegar a
la precisión deseada, y no hay garantía para la convergencia.
Fuente: Saadat, H., Power System Analysis
EJEMPLO DEL MÉTODO DE GAUSS-SEIDEL

Fuente: Saadat, H., Power System Analysis

MÉTODO DE GAUSS-SEIDEL APLICADO A SISTEMA DE POTENCIA
• La complejidad de obtener una solución formal para el flujo de
potencia en un sistema eléctrico se debe a las diferencias en el tipo
de datos especificados para las diferentes clases de barras.
• A cada cálculo de un nuevo conjunto de voltajes se le llama iteración.
El proceso iterativo se repite hasta que los cambios en cada barra son
menores que un valor mínimo especificado.
• Primero 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 reales y reactiva
programadas, respectivamente, que entran a la red en la barra 2 , se
tiene de la ecuación del slide 5 con 𝑖 igual a 2 y N igual a 4.

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia

MÉTODO ITERATIVO DE GAUSS-SEIDEL
• Al despejar el valor de 𝑉2 se tiene:

• De igual forma se proceder para las barras 3 y 4.

• La solución se obtiene por la iteración que se basa en las potencias
real y reactiva programadas en las barras 2, 3 y 4, el voltaje en la
barra de compensación programado 𝑉1 = |𝑉1 |∡𝛿1 y las estimaciones
(0) (0) (0)
iniciales de voltaje 𝑉2 , 𝑉3 y 𝑉4 en las otras barras.
• Iterando la ecuación anterior tenemos:
(1) 1 𝑃2,𝑝𝑟𝑜𝑔 − 𝑗𝑄2,𝑝𝑟𝑜𝑔 (0) (0)
𝑉2 = 0 ∗
− 𝑌21 𝑉1 + 𝑌23 𝑉3 + 𝑌24 𝑉4
𝑌22 𝑉 2
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
MÉTODO ITERATIVO DE GAUSS-SEIDEL
(1) (0)
• El valor calculado de 𝑉2 y el valor estimado de 𝑉2 no serán iguales. Se puede
alcanzar una diferencia insignificante después de varias iteraciones y podrá ser el
valor correcto de 𝑉2 con los voltajes estimados, pero sin considerar la potencia en
las otras barras.
• A medida que se encuentra el voltaje correcto en cada barra, su valor se va
usando para (1)
calcular el voltaje correcto en la siguiente barra. Por lo tanto, al
sustituir 𝑉2 en la ecuación de la barra 3, se obtiene:

• 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 encuentran valores calculados para cada variable
de estado.
• Este proceso se lleva a cabo una y otra vez hasta que la cantidad por corregir en
el voltaje de cada barra es menor que algún índice de precisión determinado
previamente. A este proceso de solución de las ecuaciones de flujo de potencia
se le conoce como el método iterativo de Gauss-Seidel.
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
MÉTODO ITERATIVO DE GAUSS-SEIDEL
• La ecuación general para el voltaje calculado en cualquier barra i de
un sistema de barras, donde se programa P y Q, es:

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potenciac

Ejemplo de Gauss-Seidel
En la figura siguiente se muestra un diagrama unifilar de un sistema de
potencia sencillo. Los generadores están conectados en las barras 1 y 4,
mientras las cargas se indican en las cuatro barras. Los valores base para el
sistema de transmisión son 100 MVA y 230 kV. Los datos de las líneas están
en la tabla 1 donde están las impedancias series en por unidad y las
susceptancias de carga de la línea para los circuitos equivalentes nominales 𝜋
de las cuatro líneas identificadas por las barras en las que terminan. Los
datos de barras en la tabla 2 enlistan los valores para P, Q y V en cada barra.
Los valores de Q se calculan con los valores de P bajo el supuesto de un
factor de potencia de 0.85. Los valores programados totales, 𝑃𝑖,𝑝𝑟𝑜𝑔 y 𝑄𝑖,𝑝𝑟𝑜𝑔
son negativos en las barras de carga 2 y 3. No se especifica la 𝑄𝑔𝑖 generada
donde la magnitud del voltaje es constante. En la columna de voltajes, los
valores para las barras de carga son estimaciones de inicio plano. La
magnitud de |𝑉1 | y el ángulo 𝛿1 de la barra de compensación, así como la
magnitud de |𝑉4 | en la barra 4, se mantienen constante en los valores que se
enlistan. Se hace un estudio de flujo de potencia por el método de Gauss-
Seidel. Encuentre el valor de 𝑉2 para la primera iteración suponiendo que los
cálculos iterativos comienzan en la barra 2.
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
Factor de aceleración
• La experiencia con el método Gauss-Seidel para la solución de flujos
de potencia ha demostrado que se puede reducir el numero de
iteraciones requeridas si la corrección en el voltaje de cada barra se
multiplica por alguna constante que incremente la cantidad de
corrección para que el voltaje sea mas cercano al valor al que se esta
aproximando.
• La diferencia entre el valor de voltaje que recientemente se ha
calculado y el mejor que previamente se evaluó en la barra, se
multiplica por el factor de aceleración apropiado para obtener una
mejor corrección que se añadirá al valor previo.

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia

Valor acelerado
• De manera general, el valor acelerado para la barra i durante la iteración k esta
dado por:

• Si 𝛼 = 1, entonces el valor de 𝑉𝑖 calculado por Gauss-Seidel se almacena como el

valor actual.
• Si 0 < 𝛼 < 1, entonces el valor para ser almacenado es un promedio ponderado
del valor de Gauss-Seidel del valor almacenado en la iteración previa.
• Si 1 < 𝛼 < 2, entonces el valor a ser almacenado es esencialmente uno que esta
extrapolado.
• Por lo general en los estudios de flujo de potencia, 𝛼 tiene un valor de
aproximadamente 1.6 y no puede exceder a 2 si la convergencia esta por ocurrir.

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia

 • Se realiza un cálculo similar para la barra 3 y del siguiente valor:

• La barra 4 debe tratarse de manera diferente puesto que es de voltaje

controlado.

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia

Barras de voltaje controlado
• Cuando en una barra i se especifica la magnitud del voltaje en lugar
de la potencia reactiva, las componentes real e imaginaria del voltaje
para cada iteración se encuentra calculando primeramente un valor
para la potencia reactiva. De la ecuación del slide 5 se tiene:

• Ecuación con iteraciones k:

• Donde 𝐼𝑚 quiere decir “parte imaginaria de” y los subíndices indican

la iteración apropiada.
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
Barras de voltaje controlado
(𝑘)
• Luego que encontremos a 𝑄𝑖 lo sustituimos en la ecuación del slide
(𝑘)
18, para encontrar un nuevo valor de 𝑉𝑖 . La potencia reactiva de la
barra 4 se calcula con la siguiente formula:

Retornar
Slide 29
• Al sustituir la potencia reactiva se tiene:

Retornar
Slide 30

• Como 𝑉4 es una cantidad especificada, se corrige la magnitud de

(1)
𝑉4 de la siguiente forma:

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia

Barras de voltaje controlado
• Desde un punto de vista practico, la salida 𝑄𝑔 del generador debe caer
dentro de limites definidos dados por la desigualdad:

• Donde 𝑄𝑚𝑖𝑛 𝑦 𝑄𝑚𝑎𝑥 son los limites mínimos y máximos de diseño del
generador. Si durante el transcurso de la solución de los flujos de potencia
el valor calculado de 𝑄𝑔 esta fuera de cualquiera de los límites, entonces
𝑄𝑔 se hace igual al valor del límite que se violo, la magnitud del voltaje
originalmente especificado se expande a otros valores y la barra se trata
como barra P-Q para la que un nuevo voltaje se calcula por medio del
programa de computadora. En las iteraciones subsecuentes, el programa
intentara sostener el voltaje especificado en la barra mientras se asegura
que 𝑄𝑔 esta dentro del rango permitido de valores. Esto podría ser posible
por que otros cambios pueden ocurrir en otro lado del sistema para apoyar
la acción local de la excitación del generador, conforme se ajusta para
satisfacer el voltaje en terminales especificado.
Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia
Continuación del ejemplo
• Encuentre el voltaje en la barra 4 para completar la primera iteración
del procedimiento de Gauss-Seidel, considerando los valores
acelerados anteriormente indicados en lugar de los voltajes
estimados en las barras 2 y 3.
• En la tabla de admitancias se observa que 𝑌41 es igual a cero y así, la
ecuación del slide 27 da:

• Al sustituir los valores, se tiene:

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia

 • El valor encontrado de la potencia reactiva en la barra 4 se sustituye
en la ecuación del slide 27, para obtener:

(1)
• Por lo tanto, 𝑉4
es igual a 1.017929 y así, se debe corregir la
magnitud estimada en la barra 4.

Fuente: Grainger, J. y Stevenson, W. Análisis de Sistemas de Potencia

Ecuación de flujo de potencia
• Aplicando la LCK para un sistema de n barras, el resultado es:

(1)

(2)

• La potencia real y reactiva en la barra i es:

(3)

(4)

Fuente: Saadat, H., Power System Analysis

Ecuación de flujo de potencia
• Sustituyendo a 𝐼𝑖 en ecuación (2):
(5)

• Esta es una ecuación no lineal que puede ser resuelta mediante

iteraciones.

Fuente: Saadat, H., Power System Analysis

Método Gauss-Seidel
• En el estudio de flujo de potencia, es necesario resolver un conjunto
de ecuaciones no lineales representadas por la ecuación (5) para dos
variables desconocidas en cada nodo. Con el método Gauss-Seidel se
puede determinar 𝑉𝑖 con una secuencia de iteraciones:

(6)

• Donde 𝑦𝑖𝑗 en letras minúsculas es la admitancia actual en por unidad.

𝑃𝑖𝑠𝑐ℎ 𝑦 𝑄𝑖𝑠𝑐ℎ son la potencia real y reactiva programadas en por
unidad.
• En las barras de generación la potencia activa y reactiva tienen
valores positivos y en las barras de carga tienen valores negativos. Si
utilizamos la ecuación (5) para encontrar 𝑃𝑖 𝑦 𝑄𝑖 , tenemos:
(7)

(8)

• La ecuación de flujo de potencia es usualmente expresada en

términos de los elementos de la matriz de admitancias. Desde los
elementos fuera de la diagonal de la matriz de admitancia de barra
𝑌𝑏𝑢𝑠 , mostrada por letras mayúsculas, son 𝑌𝑖𝑗 = −𝑦𝑖𝑗 , y los elementos
de la diagonal son 𝑌𝑖𝑖 = σ 𝑦𝑖𝑗 , al sustituir en la ecuación (6) tenemos:
 • 𝑌𝑖𝑖 incluida la admitancia de
línea a tierra.
(9)
• En la barra slack son dados
ambos componentes del voltaje.
• Hay 2(𝑛 − 1) ecuaciones para
(10) ser resultas por un método
iterativo.
(11) • La magnitud y ángulo de fase del
voltaje en las barras de cargas
está por debajo al valor de la
barra slack.
• Para las barras 𝑃 − 𝑄, la potencia real y reactiva 𝑃𝑖𝑠𝑐ℎ 𝑦 𝑄𝑖𝑠𝑐 son
conocidas. Iniciando con un estimado inicial. La ecuación (9) es
resuelta para encontrar las componentes real e imaginaria del voltaje.
• Para las barras de control de volteje (barras P-V) donde 𝑃𝑖𝑠𝑐ℎ 𝑦 |𝑉𝑖 |
son especificados, primero resuelve con la ecuación (11) para 𝑄𝑖𝑠𝑐ℎ y
(𝑘+1)
luego es usada la ecuación (9) para resolver 𝑉𝑖 . Sin embargo,
(𝑘+1)
como |𝑉𝑖 | es especificado, solo la parte imaginaria de 𝑉𝑖 es
retenida,
(12)

(13)
Factor de aceleración

• Donde ∝ es el factor de aceleración. Su valor depende del sistema. Se

encuentra en el rango de 1.3 a 1.7 para sistemas eléctricos de
potencia típicos.
• Se realiza un proceso de iteración hasta que los cambios de las
componentes real e imaginaria del voltaje de barra estén dentro de
una exactitud especificada:
Un error de exactitud típico es
0.001 pu
Flujo de líneas y perdidas
Método de Newton-Raphson
Método de Newton-Raphson
• El método más utilizado para resolver de forma simultaneas
ecuaciones algebraica no lineales es el método Newthon-Raphson.
• El método de Newton es un procedimiento de aproximación sucesiva
basado en una estimación inicial de lo desconocido y usa la serie de
expansión de Taylor. Considerar la solución de la ecuación
unidimensional dada por:
𝑓 𝑥 =𝑐
• Si 𝑥 (0) es un valor estimado inicial de la solución, y ∆𝑥 (0) es una
pequeña desviación de la solución correcta, obtenemos:
𝑓 𝑥 0 + ∆𝑥 0 = 𝑐

Fuente: Saadat, H., Power System Analysis

 • Expandiendo a mano izquierda la ecuación anterior en serie de Taylor
sobre 𝑥 (0) :

• Asumiendo que el error ∆𝑥 (0) es muy pequeño, los términos de

orden superior pueden ser despreciados, el cual resulta en:

• Donde:

Fuente: Saadat, H., Power System Analysis

 • Añadiendo ∆𝑥 (0) para un estimado inicial resultara en la segunda
aproximación:

• El uso sucesivo de este procedimiento produce el algoritmo de

Newton-Raphson:

Fuente: Saadat, H., Power System Analysis

 • Se puede arreglar como:

• Donde

• La relación de la ecuación anterior demuestra que la ecuación no

lineal 𝑓 𝑥 − 𝑐 = 0 es aproximadamente por la tangente en la curva
de 𝑥 (𝑘) . Por lo tanto, una ecuación lineal es obtenida en términos de
los pequeños cambios en la variable. La intersección de la línea
tangente con el eje x da como resultado 𝑥 (𝑘+1) .
Fuente: Saadat, H., Power System Analysis
Ejemplo del Método Newton-Raphson
• Use el método Newton-Raphson encontrar una raíz de la ecuación
dada en el ejemplo que utilizamos el método de Gauss-Seidel.
Asumiendo un estimado inicial de 𝑥 (0) = 6.
𝑓 𝑥 = 𝑥 3 − 6𝑥 2 + 9𝑥 − 4 = 0

Fuente: Saadat, H., Power System Analysis

 • La solución analítica dada por el algoritmo de Newton-Raphson es:

• Por lo tanto, el resultado en el fin de la primera iteración es:

Fuente: Saadat, H., Power System Analysis

 • Las subsecuentes iteraciones resultan en:

Fuente: Saadat, H., Power System Analysis

Realizado en MATLAB

Fuente: Saadat, H., Power System Analysis

Ecuaciones n-dimensional
• Ahora vamos a considerar la ecuación n-dimensional dada por:

• Expandiendo a mano izquierda las ecuaciones con la serie de Taylor

Fuente: Saadat, H., Power System Analysis

En forma de Matriz

En la forma corto:

Fuente: Saadat, H., Power System Analysis

Algoritmo de Newton-Raphson (n-dimensional)
𝐽(𝑘) es llamada la matriz Jacobiana. Los
elementos de esta matriz son las derivadas
Donde: parciales evaluadas en 𝑋 (𝑘) . Esto asumiendo que
𝐽(𝑘) tiene un inverso durante cada iteración. El
método de Newthon, al aplicar un conjunto de
ecuaciones no lineales, reduce el problema a
resolver un conjunto de ecuaciones lineales para
determinar los valores que mejoran la precisión
de las estimaciones.

La solución por inversión de la matriz es muy

ineficiente. Eso no es necesario para obtener el
inverso de 𝐽(𝑘) . En cambio, se obtiene una
solución directa mediante una factorización
Fuente: Saadat, H., Power System Analysis
triangular óptimamente ordenada
Ejemplo aplicando el método Newthon-Raphson
Gráficamente, la solución para este sistema esta
representado por la intersección de el circulo con la
curva. Los puntos intersecados son (1, -1.7) y (-1.8,
0.8).

Valores iniciales: 𝑥1 = 0.5 𝑦 𝑥2 = −1

Fuente: Saadat, H., Power System Analysis

SOLUCIÓN DE FLUJO DE POTENCIA CON
NEWTON-RAPHSON
SOLUCIÓN DE FLUJO DE POTENCIA CON NEWTON-RAPHSON
Barra de carga

Fuente: Saadat, H., Power System Analysis

Reducción de la matriz
• En lo anterior, se supone que el bus 1 es el bus inactivo. La matriz
jacobiana da la relación linealizada entre pequeños cambios en el
(𝑘)
ángulo de voltaje ∆𝛿 (𝑘) y la magnitud del voltaje ∆ 𝑉𝑖 con los
(𝑘) (𝑘)
pequeños cambios en la potencia real y reactiva ∆𝑃𝑖 𝑦 ∆𝑄𝑖 .
• Los elementos de la matriz Jacobiana son las derivadas parciales de
las ecuaciones de las potencias real y reactiva, evaluada en ∆𝛿 (𝑘) y
(𝑘)
∆ 𝑉𝑖 . En forma corta, se puede escribir:

−1
∆𝛿 𝐽1 𝐽2 ∆𝑃
=
∆|𝑉| 𝐽3 𝐽4 ∆𝑄
Fuente: Saadat, H., Power System Analysis
Barras de voltaje controlado
• En este tipo de barras, las magnitudes de los voltajes son conocidas.
Sin embargo, si m barras del sistema son de voltaje controlado, m
ecuaciones involucran ∆𝑄 𝑦 ∆𝑉 y las columnas correspondientes de la
matriz Jacobiana son eliminadas. Hay n-1 restricciones de potencia
real y n-1-m restricciones potencia reactiva, respectivamente, y la
matriz Jacobiana es del orden 2𝑛 − 2 − 𝑚 ∗ 2𝑛 − 2 − 𝑚 .
• 𝑱𝟏 es del orden 𝑛 − 1 ∗ 𝑛 − 1 , 𝑱𝟐 es del orden 𝑛 − 1 ∗
𝑛 − 1 − 𝑚 , 𝑱𝟑 es del orden 𝑛 − 1 − 𝑚 ∗ 𝑛 − 1 , y 𝑱𝟒 es del
orden 𝑛 − 1 − 𝑚 ∗ (𝑛 − 1 − 𝑚)

Fuente: Saadat, H., Power System Analysis

 • Los elementos de la diagonal y fuera de la
diagonal de 𝐽1 son: Ecuaciones principales

• Los elementos de la diagonal y fuera de la

diagonal de 𝐽2 son:

Fuente: Saadat, H., Power System Analysis

 • Los elementos de la diagonal y fuera de la
diagonal de 𝐽4 son: Ecuaciones principales

• Los elementos de la diagonal y fuera de la

diagonal de 𝐽4 son:
❖ Los valores estimados para la
barra de voltaje contralado
son:

Fuente: Saadat, H., Power System Analysis

Ejemplo del método Newton-Raphson
• Obtiene la solución del flujo de potencia con el método Newton-
Raphson del siguiente circuito.

Fuente: Saadat, H., Power System Analysis

 • Las impedancias de líneas convertidas en admitancias son:
𝑦12 = 10 − 𝑗20
𝑦13 = 10 − 𝑗30
𝑦23 = 16 − 𝑗32
• Estos resultados en la matriz de admitancias de la barra:

20 − 𝑗50 −10 + 𝑗20 −10 + 𝑗30

𝑌𝑏𝑎𝑟𝑟𝑎 = −10 + 𝑗20 26 − 𝑗52 −16 + 𝑗32
−10 + 𝑗30 −16 + 𝑗32 26 − 𝑗62
• Llevamos la matriz de admitancias a valores polares:

53.85∡ − 1.9029 22.36∡2.0344 31.62∡1.89

𝑌𝑏𝑎𝑟𝑟𝑎 = 22.36∡2.0344 58.14∡ − 1.107 35.77∡2.0344
31.62∡1.89 35.77∡2.0344 67.231∡ − 1.174

Fuente: Saadat, H., Power System Analysis

 • Las expresiones para la potencia real en las barras 2 y 3, y la potencia
reactiva e la barra 2 son:

• Los elementos de la matriz Jacobiana son obtenidos realizando las

derivadas de las ecuaciones anteriores con respecto a 𝛿2 , 𝛿3 𝑦 |𝑉2 |.

Fuente: Saadat, H., Power System Analysis

Fuente: Saadat, H., Power System Analysis
 • La carga y generación expresada en por unidad son:

• El voltaje de la barra slack es 𝑉1 = 1.05∡0° 𝑝𝑢, y la magnitud de

(0) (0) (0)
voltaje es 𝑉2 = 1.0 , 𝛿2 = 0.0 , y 𝛿3 = 0.0 , la potencia
residuales son computados:

Fuente: Saadat, H., Power System Analysis

 • Evaluado los elementos de la matriz Jacobiana con el estimado inicial,
el conjunto de ecuaciones lineales en la primera iteración:

Fuente: Saadat, H., Power System Analysis

 • Segunda iteración:

• Tercera iteración:

Fuente: Saadat, H., Power System Analysis

 • La solución converge en 3 iteraciones con un máximo error en
potencia de 2.5 𝑥 10−4 . Resultados:
𝑉2 = 0.97168∡ − 2.696°
𝑉3 = 1.04∡ − 0.4988°
• Luego se procede a obtener la generación de potencia real y reactiva
de la barra 3 y la barra slack (1):

Fuente: Saadat, H., Power System Analysis

Transformadores Regulantes
Transformadores Regulantes
• Como ya sabemos el flujo de potencia real a través de una línea de
transmisión es determinado por la diferencia de ángulo de los voltajes
de las terminales o barras. La potencia real y reactiva puede ser
controlada usando un transformador con cambiador de tap y
transformadores regulantes.
• Cuando la relación está en valor nominal, el transformador es
representado por una admitancia en serie 𝑦𝑡 en por unidad.
• Cuando la relación no está con el valor nominal, la admitancia en p.u
es diferente en ambos lados del transformado, y la admitancia puede
ser modificada para incluir el efecto de la relación fuera del nominal.

Fuente: Saadat, H., Power System Analysis

 • Considerar un transformador con una admitancia 𝑦𝑡 en serie y un
transformador ideal representando una relación de tap fuera del
nominal 1:a como se muestra.

• 𝑦𝑡 es la admitancia en serie en por unidad basada en la relación

nominal y “a” es la posición del tap fuera del nominal en p.u, con el
cual se realizan pequeños ajustes en voltaje de un ±10% usualmente.
• En el caso de transformadores cambiadores de fase, “a” es un numero
complejo.

Fuente: Saadat, H., Power System Analysis

• Considerar una barra x ficticia entre • También se tiene:
la relación de giro y la admitancia
del transformador.
• Asumiendo la dirección de las
corrientes, tenemos: • Sustituyendo a 𝐼𝑖 :

• Resultado en forma de matriz:

• Sustituyendo a 𝑉𝑥 :
Fuente: Saadat, H., Power System Analysis
Ejemplo

• Encontrar la matriz de admitancia.

Modelo 𝜋:
Matriz de admitancia
Grupo 2
Grupo 3