FLUJOS DE POTENCIA EN SEP

ING. JESÚS E. SALAMANCA JAIMES

1. INTRODUCCIÓN
• Los estudios de flujo de potencia son de vital importancia en la
planeación y diseño de la expansión futura de los SEP, así como en la
determinación de las mejores condiciones de operación en SEP existentes
• Las variables obtenidas de un estudio de flujos de potencia es la
magnitud de la tensión (|V|) y la fase de la tensión en cada barra (δ),
así como la potencia real (P) y la potencia reactiva (Q) que fluyen en
cada barra del sistema
• Se realiza principalmente usando YBARRA
• Las líneas de transmisión se representan con su equivalente monofásico
en π
SISTEMA TIPO DE 4 BARRAS
2. ECUACIONES DE FLUJO DE POTENCIA
Pi y Qi son las potencias totales real y Pi _ prog = Pg _ i − Pd _ i
reactiva calculadas que entran a la red a
través de la barra i
N
Pi − jQi = Vi *
Y V
j =1
ij j

Pi =  Yij Vi V j cos (ij +  j −  i )

N

j =1

Qi = − Yij Vi V j sin (ij +  j −  i )

N

j =1
Qi _ prog = Qg _ i − Qd _ i
Pg_i es la potencia programada que se esta
generando en la barra i, Pd_i es la potencia
programada que demanda la carga
conectada a la barra i; de modo que
Pi _ prog = Pg _ i − Pd _ i da la potencia
programada total que esta siendo inyectada
dentro de a red a través de la barra i
MATRIZ DE IMPEDANCIA O ZBARRA
Los errores ΔP y ΔQ se obtienen de la diferencia entre los valores programados
(Pi_prog, Qi_prog), y los valores calculados (Pi y Qi), dichos errores aparecen cuando
durante el desarrollo de la solución de un problema de flujos de potencia los valores
programados no son iguales a los valores calculados
3. TIPOS DE BARRAS
1- Barra de compensación – Slack
• Es la barra 1 de la red
• La tensión en esta barra será la de referencia para la red por lo que δ1 será la
referencia para los ángulos de todas las demás tensiones de barra
• Por conveniencia se asigna δ1=0°,
• No se definen errores para la barra de compensación
• Las variables conocidas son |V1| y δ1
• Las variables desconocidas corresponden a las P y Q necesarias para mantener la
tensión de barra en los valores programados
• Suele corresponder a un generador de gran tamaño del sistema
3. TIPOS DE BARRAS
2- Barra de tensión controlada – Pv:
• Corresponden a cualquier barra del sistema que puede mantener la magnitud de la
tensión constante
• Normalmente corresponden a barras de generación
• También pueden ser barras con elementos compensadores de potencia reactiva a
base de electrónica de potencia que permiten controlar la tensión (FACTS, SVC), en
este caso Pg_i=0
• Las variables conocidas son Pg_i y |Vi|
• La potencia reactiva Qg_i que se requiere para mantener la tensión programada
|Vi| no se puede conocer de forma anticipada, después que se ha resuelto el
problema de flujo de potencia se puede calcular Qi
• Solo se dispone de el error ΔP y se tiene una variable desconocida: el ángulo de la
tensión δi
3. TIPOS DE BARRAS
3- Barra de carga – PQ:
• Corresponde a las barras del sistema que no tienen generación (Pg_i=Qg_i=0)
• Qi_prog= -Qd_i, Pi_prog= - Pd_i
• Los valores conocidos son las potencias P y Q demandadas en la barra que se
obtienen de registros históricos
• Los valores desconocidos corresponden a |Vi| y δi
• En la practica se conoce la potencia real y la potencia reactiva se calcula sobre un
factor de potencia supuesto, el cual debe ser mayor a 0.85
• Los errores ΔP y ΔQ se pueden definir con facilidad
Después de obtener los valores de |Vi| y δi en cada barra, usando la Ybarra, se
pueden calcular las corrientes por cada línea de a red y posteriormente el flujo de
potencia en todo el sistema; las perdidas de potencias entonces serán:
TIPOS DE BARRAS
4. MÉTODO DE GAUSS-SEIDEL
Aunque es posible formular un numero de ecuaciones para el sistema igual al numero
de variables de estado del sistema, la solución de forma directa por medio del sistema
de ecuaciones es poco conveniente, por lo cual se usan métodos iterativos que suponen
una solución inicial a las variables de estado y calculan nuevos valores para las
tensiones en cada barra, en el método GS, el proceso se realiza repetitivamente hasta
que las diferencias en las tensiones de barra de una iteración a otra son menores a un
error relativo determinado previamente
Debido a que la barra de compensación es la numero 1 y allí se conoce la magnitud
de tensión y su ángulo, los cálculos empiezan a partir de la barra 2
A medida que se va calculando el valor correcto de tensión en cada barra, este valor
se va usando para calcular el valor correcto de la siguiente barra hasta tener las
tensiones en cada barra del sistema y completar la iteración
N
Pi − jQi = Vi *
Y V
j =1
ij j
MÉTODO DE GAUSS-SEIDEL
Se evita la convergencia sobre una solución errónea si los valores iniciales supuestos
son de magnitud razonable y no difieren demasiado en fase, comúnmente a los valores
iniciales supuestos de las tensiones desconocidas en todas las barras de carga se les
asigna 1∟0°, tal inicio se conoce como inicio plano debido a la suposición de un
perfil uniforme de tensiones. De la ecuación general de flujo de potencia:
N
Pi − jQi = Vi *
Y V
j =1
ij j

Se despeja y resuelve la sumatoria obteniendo la ecuación general para la tensión

calculada en cualquier barra i de un sistema de N barras donde se programan P y Q,
en barras PQ esta expresión es suficiente para obtener sus incógnitas |Vi| y δi
(el superíndice k indica el numero de iteración en la que se esta calculando la tensión,
el * indica valores complejos conjugados)
1  Pi , prog − jQi , prog i −1 N
( k −1) 
Vi ( k ) =  ( k −1)*
−  YijV j −  YijV j 
(k )

Yii  Vi j =1 j =i +1 
MÉTODO DE GAUSS-SEIDEL
Por ejemplo, para un sistema de 4 barras, donde la barra 1 es de compensación y las
barras 2, 3 y 4 son de carga las tensiones en cada barra se calculan por medio de:
1  P2, prog − jQ2, prog (0) 
V2 =
(1)
 (0)*
− Y21V1 − Y23V3 − Y24V4 
(0)

Y22  V2 
1  P3, prog − jQ3, prog (0) 
V 3
(1)
=  (0 )*
− Y31V1 − Y32V2 − Y34V4 
(1)

Y33  V3 
1  P4, prog − jQ4, prog (1) 
V 4
(1)
=  (0 )*
− Y41V1 − Y42V2 − Y43V3 
(1)

Y44  V4 
Se debe notar que a medida que se obtiene un resultado, este se usa inmediatamente
en la siguiente expresión; después de terminar todas las iteraciones se puede hallar la
P y Q de la barra Slack despejándolas de su respectiva expresión de V1
P1 − jQ1 = V1* Y11V1 + Y12V2 + Y13V3 + Y14V4 
MÉTODO DE GAUSS-SEIDEL

Se puede reducir el numero de iteraciones si se corrige cada resultado de cada tensión

usando un factor de aceleración α, el cual normalmente se encuentra entre 1.4 y 1.6,
de este modo para la tensión 2:
V2(1)_ ac = V2(0) +  (V2(1) − V2(0) )

En las barras PV, previo al calculo de la tensión de barra Vi, se debe calcular la
potencia reactiva necesaria en dicha barra para mantener la |Vi| en el valor
especificado
 ( k −1)*  i −1 N 
( k −1)  
Qi( k ) = − Im Vi   YijV j +  YijV j  
(k )

  j =1 j =i  
MÉTODO DE GAUSS-SEIDEL
Luego de obtener Qi, esta se reemplaza en la expresión de calculo de Vi de la
respectiva barra incluyéndola como Qi_prog

1  Pi , prog − jQi , prog i −1 N

( k −1) 
Vi ( k ) =  ( k −1)*
−  YijV j −  YijV j 
(k )

Yii  Vi j =1 j =i +1 

Luego de calculado |Vi| en la barra PV, posiblemente no se obtenga la |V|

especificada, como es una restricción del sistema, esta se debe corregir para asegurar
esta magnitud de tensión en dicha barra, en la practica se hace por medio de la
corriente de campo del generador, la |V| del valor corregido calculado debe ser el
valor especificado en la barra
Vi k −1 *Vi ( k )
Vi (_kcorr
)
=
Vi ( k )
MÉTODO DE GAUSS-SEIDEL
Ejemplo 2: La figura muestra un sistema de potencia con generación en las barras 1 y
3, la magnitud del voltaje es 1,05 por unidad en la barra 1 y 1.04 en la barra 3
donde se inyecta un potencia real de 200 MW. La carga programada en la barra 2
es de 400MW y 250 MVAR, las impedancias se encuentran en por unidad sobre una
base de 100 MVA. Por medio del método Gauss-Seidel obtenga las P, Q y tensiones
en las barras del sistema
5. MÉTODO NEWTON - RAPHSON
• Es un procedimiento algorítmico que permite hallar las raíces de
funciones conocido un valor numérico cercano a la raíz con muy rápida
convergencia
• Su convergencia es cuadrática, lo que lo hace matemáticamente superior
y rápido comparado con GS que tiene convergencia lineal
• Es menos propenso a la divergencia en problemas mal condicionados
• En la actualidad este es el método iterativo más eficiente que se utiliza
para el análisis de flujos de potencia en redes eléctricas
• Se basa en la expansión de la serie de Taylor para una función de dos o
mas variables alrededor de un punto de solución
• Se usa la notación en coordenadas polares en vez de rectangulares
• Los ángulos de las variables deben estar en radianes
5. MÉTODO NEWTON - RAPHSON
La solución del método de NR parte de la ecuación de flujo de potencia
formulada en forma polar, la cual al ser decompuesta permite obtener las
expresiones de Pi y Qi N
Pi − jQi = Vi * I i = Vi *  YijV j
j =1

Pi =  Yij Vi V j cos (ij −  i +  j ) Qi = − Yij Vi V j sin (ij −  i +  j )

N N

j =1 j =1

Dichas expresiones permiten obtener un conjunto de ecuaciones

algebraicas no lineales en términos de las variables independientes: |V|
y δ (en radianes). En las barras de carga se tienen dos ecuaciones ya que
se conoce P y Q, en las barras de tensión controlada solo se conoce una
ecuación al conocer solo su P
5. MÉTODO NEWTON - RAPHSON
El uso de la matriz jacobiana permite la linealización de un sistema no
lineal alrededor de un punto de solución usando pequeños cambios en δ y
|V|
 P   J1 J 2    
 Q  =  J J    V 
   3 4 
Residuales de Incógnitas
Potencias Jacobiana

Los elementos de la matriz jacobiana son las derivadas parciales de las

ecuaciones de P y Q con respecto a Δδi y Δ|Vi|, J tiene m filas y m
columnas, donde m es el numero de incógnitas: la barra Slack no aporta
incógnitas, las barras PQ aportan 2 incógnitas y las barras PV 1 incógnita
5. MÉTODO NEWTON - RAPHSON
Solucionando las expresiones de potencia por medio de la serie de Taylor,
se obtiene:
 P2( k )   P2( k ) P2( k ) P2( k ) P2( k )    ( k ) 
     n  V2  Vn  
 2

   2

     

 (k )   (k ) (k ) 
 
 Pn   Pn Pn( k ) Pn( k ) Pn   ( k )
 n 
 P   J1 J 2          n  V2  Vn   
= 
 Q  =  J   2
  V  Q2( k )   Q2( k ) Q2( k ) Q2( k ) (k )  
Q2 

   3 J4          V2
(k ) 

    2  n  V2  Vn 
     

 (k )   (k ) (k ) 
 
 Qn   Qn Qn( k ) Qn( k ) Qn   
      V (k )

   2  n  V2  Vn   n 
5. MÉTODO NEWTON - RAPHSON
Los valores residuales de potencia se obtienen como:

Pi (k )
= Pi _ prog − Pi (k )

Qi( k ) = Qi _ prog − Qi( k )

Los nuevos valores estimados para las tensiones y ángulos de las barras
son:
 i( k +1) =  i( k ) +  i( k )
Vi ( k +1) = Vi ( k ) +  Vi ( k )
5. ALGORITMO NR
1. En las barras PQ establecer como valor inicial Vi=1∟0°,en barras PV
asignar la tensión especificada y δi=0°
2. En barras de carga y compensación calcular los valores que apliquen
de Pi(k) , Qi(k), ΔPi(k) y ΔQi(k), para después calcular la matriz de
residuales de potencia
3. Obtener los elementos de la matriz jacobiana
4. Resolver el sistema de ecuaciones matricial para obtener los Δ|Vi(k)| y
Δδi(k)
5. Comprobar que los nuevos valores de ΔPi(k) y ΔQi(k) sean menores al
error permisible establecido, sino continuar iterando
MÉTODO DE NR
Ejemplo 3: La figura muestra un sistema de potencia con generación en las barras 1 y
3, la magnitud del voltaje es 1,05 por unidad en la barra 1 y 1.04 en la barra 3
donde se inyecta un potencia real de 200 MW. La carga programada en la barra 2
es de 400MW y 250 MVAR, las impedancias se encuentran en por unidad sobre una
base de 100 MVA. Por medio del método Newton-Raphson obtenga las P, Q y
tensiones en las barras del sistema. Considere un error de 2.5x10-4
6. FLUJOS DE POTENCIA EN LAS LÍNEAS Y PERDIDAS
Luego de la solución iterativa de las tensiones de barra, el siguiente paso es el calculo
de los flujos de potencia en las líneas y sus perdidas, considere la línea que conecta
las barras i con j, Iij es la corriente de línea medida en la barra i, Iji es la corriente de
línea medida en la barra j

I ij = I L + I i 0 = Yij (Vi − V j ) + Yi 0Vi

I ji = − I L + I j 0 = Yij (V j − Vi ) + Y j 0V j

Las potencias complejas que fluyen desde el bus i hacia el j y desde el bus j hacia el i
son: Sij = Vi I ij* S ji = V j I *ji

Las perdidas en la línea entonces serán:

S L _ ij = Sij + S ji
PROCEDIMIENTO DE SOLUCION DE FLUJOS DE POTENCIA

1. Obtener la matriz de admitancias

2. Definir los valores programados en cada barra del sistema
3. Aplicar GS o NR
a) *Desarrollar las iteraciones hasta obtener el error deseado
b) *Actualizar valores de Q en las barras PV teniendo en cuenta las ultimas tensiones
obtenidas por el proceso iterativo
c) *Calcular P y Q en la barra 1
4. Calcular las corrientes en las líneas de transmisión
5. Obtener el flujo de potencias en las líneas de transmisión
6. Calcular las perdidas en las líneas de transmisión
FLUJOS DE POTENCIA EN LAS LÍNEAS Y PERDIDAS
Ejemplo 4: Determine el flujo de potencia y las perdidas en las líneas del
ejemplo 2, construya un diagrama indicando la dirección de los flujos