ELT - 2761 SISTEMAS DE POTENCIA I

AUXILIAR: Univ. Paul Danny Alba Lozada “JHAMBAO”

Ejercicios De Flujos De Carga 2do EXAMEN PARCIAL
PRIMER METODO DE SOLUCION DE FLUJOS DE CARGA

1. METODO DE GAUSS-SEIDEL
Al finalizar el estudio del método, el capacitando tendrá la habilidad de obtener la solución de los Flujos de Carga,
para que a través de ésta pueda operar de una manera segura un sistema Eléctrico de Potencia.

1.1. INTRODUCCION
El primer método desarrollado para la solución de los Flujos de Carga, fue el método de Gauss-seidel. Este método
tiene la ventaja de que es muy sencillo de programar y de que los requerimientos de memoria de computadora necesarios
son mínimos. La gran desventaja de este método es que el número de iteraciones necesarias para obtener una solución es
muy grande y se incrementa a medida que el sistema procesado crece. Esto se debe a que el efecto de un ajuste en el voltaje
de un bus durante una iteración, se refleja únicamente en los buses adyacentes a éste, de tal manera que se requieran
varias iteraciones para que los ajustes se propaguen a través del sistema trayendo como consecuencia un mayor tiempo
de cómputo. En la actualidad esto se reduce en gran medida debido a la rapidez de los, procesadores actuales.
El método no funciona para sistemas donde la matriz Ybus no es diagonalmente dominante ( Sistemas Predominantemente
Radiales ), ni en sistemas donde se manejan impedancias negativas ( Susceptancias Capacitivas ).

1.2 DESARROLLO
Para poder deducir el algoritmo correspondiente al método de Gauss-Seidel utilizaremos la formulación matricial que
relaciona los voltajes y las corrientes nodales para una red de “ n” nodos, siendo esta:

I   Y V  (1)
Desarrollando la expresión matricial ec.( 1 )
 I 1   Y11 Y12 Y13  Y1P  Y1N   V1 
 I  Y   Y2 N  V2 
 2   21 Y22 Y23 Y2 P
 I 3   Y31 Y32 Y33  Y3 P  Y3 N   V3 
    
             (2)
 I P  YP1 YP 2 YP 3  YPP  YPN  VP 
    
             
 I  Y   YNN  V N 
 N   N 1 YN 2 YN 3 YNP
Desarrollando la ecuación para el nodo “P”, se obtiene:
I P  YP1V1  YP 2V2  YP3V3  .......  YPPVP  .......  YPN VN (3)
en forma compacta:
N
I P  YPPVP   Y pqVq (4)
q 1
q p

despejando V P en la ec.( 4 )
 N

1 
VP  I P   Y pqVq 
YPP  q 1

 q p 
(5)
P  1,2,3,........, N

Para los nodos de carga, la corriente de consumo se calcula mediante la formula

S* PP  jQ P
S *  VP* I P IP  IP  (6)
VP* VP*
 ELT - 2761 SISTEMAS DE POTENCIA I
AUXILIAR: Univ. Paul Danny Alba Lozada “JHAMBAO”
Ejercicios De Flujos De Carga 2do EXAMEN PARCIAL
Sustituyendo la ec.(6) en la ec.(5)
 
1  PP  jQ P N
VP    Y pqVq  (7)
YPP  VP*
q 1

 q p 
Para poder manejar de una manera más eficiente la expresión anterior se definen las siguientes constantes:
PP  jQ P
KLP  (8)
YPP
YPQ
YL PQ  (9)
YPP
Sustituyendo las constantes anteriores ec.(8), ec.(9) en la ec.(7) se obtiene:
KLP N
VP    YLPQVq
VP* q 1 (10)
q p

P  1,2,3,........, N
Para poder representar esta ecuación dentro de un proceso iterativo, puede representarse de la siguiente forma:

KLP P 1 N
VP( K 1) 
VP*
 
q 1
YL V
PQ q
( K 1)
 
P 1
YLPQVqK
(11)
P  1,2,3,..........., N ; P  S
Donde S es el nodo compensador. Como se conoce el voltaje y ángulo para este nodo, no se necesita su ecuación
para resolver el problema de flujos de carga.

Los pasos a seguir para resolver la ec.(11) aplicando el método de Gauss-Seidel se presentan a continuación:
 Se considera un conjunto de valores iniciales de voltajes. Usualmente son voltajes planos (10  ) para el caso
de una red eléctrica.
 Esta serie de valores se utiliza para resolver el conjunto de ecuaciones (2), de manera tal que el nuevo valor
calculado VP( K 1) remplaza inmediatamente a V PK y es usado en la evaluación de las subsecuentes ecuaciones.
 El proceso se repite hasta que las diferencias de todos los voltajes entre las iteraciones ( K+1 ) y ( K ) están
dentro de una tolerancia especificada.

La convergencia del método puede ser mejorada, si los cambios de voltajes obtenidos en las diferentes iteraciones se
multiplican por un factor de aceleración empírico ( ) de la siguiente forma:
VP( K 1)  VP( K )   (VP( K 1)  VP( K ) )
VP  VP( K 1)  VP( K ) (12)
siendo 1    2
el valor aceptado en la mayoría de los casos a resolver es de 1.4, valor deducido en base a la experiencia.

1.3 NODOS DE VOLTAJE CONTROLADO

En un sistema eléctrico de potencia los voltajes de los nodos de generación permanecen constantes ( Nodos P-V ), por
lo que el procedimiento de solución al problema de Flujos de Carga visto anteriormente debe modificarse. En los nodos
de voltaje controlado normalmente se especifica la magnitud del voltaje y la Potencia Activa, así como los límites de
potencia Reactiva del generador. En el método de Gauss-Seidel, la Potencia Reactiva en un nodo de voltaje controlado
debe de calcularse antes de proceder a evaluar la ecuación de voltaje respectivo.
La potencia Real y Reactiva nodal en cualquier red se calcula mediante la ecuación:
 ELT - 2761 SISTEMAS DE POTENCIA I
AUXILIAR: Univ. Paul Danny Alba Lozada “JHAMBAO”
Ejercicios De Flujos De Carga 2do EXAMEN PARCIAL
N
PP  jQ P  VP*  YPqVq (13)
q 1
tomando en cuenta las definiciones de voltaje y admitancia:
VP  e P  jf P
Vq  eq  jf q (14)
YPq  G Pq  jB Pq
sustituyendo las expresiones anteriores en la ec.(13) y separando las expresiones de potencia Real y Reactiva se obtiene:

 
N
QP  eP2 BPP  f P2 B   f P (eP GPq  f q BPq )  eP ( f q G Pq  eq BPq ) (15)
q 1
q P

los valores de e P y f P en el nodo de voltaje controlado debe de satisfacer la siguiente igualdad:

eP2  f P2  VP
2
especificado (16)
K
Antes de calcular la potencia Reactiva Q , para un nodo ( P-V ) es necesario ajustar e P y f P que satisfaga la ec.(16).
P
Por lo tanto:
f P((Kforzado)  VP especificado sen P
) (K )

eP( K( forzado)  VP especificado cos  P

) K
(17)

f P( K )
donde tg 1
eP( K )
Si al calcular la potencia reactiva QP( K ) , usando la ec.(15), se viola alguno de los límites de generación especificado,
entonces se realiza la siguiente sustitución:
QP( K )  QPMax QP( K )  QPMin (18)
según el límite que se viole, resolviendo la ec.(7).
Con la finalidad de aplicar los conceptos teóricos previamente desarrollados, se presenta un ejemplo de una red
eléctrica básica de 3 Nodos cuyos datos y características se definen a continuación.

Ejemplo # 1.- Aplicando el Método de Gauss-Seidel, resolver el problema de flujos de carga para el sistema que se
muestra en la figura #1. En las tablas 1 y 2 se presenta toda la información relacionada con el sistema. Los límites de
potencia reactiva mínima y máxima de la unidad de generación son 0 y 35 Mvar’s. Considere un factor de aceleración de
1.4, un error máximo de 0.0001 y una potencia base de 100 MVA.
G1 G2
50+j20

1 2

60+j25

Figura # 1.- Red Eléctrica de 3 Nodos para el ejemplo # 1.

 ELT - 2761 SISTEMAS DE POTENCIA I
AUXILIAR: Univ. Paul Danny Alba Lozada “JHAMBAO”
Ejercicios De Flujos De Carga 2do EXAMEN PARCIAL

P-Q IMPEDANCIA Zpq Y'pq Ypq

1-2 0.08 + j0.24 0 1.25 - j3.75
1-3 0.02 + j0.06 0 5 - j15
2-3 0.06 + j0.18 0 1.66 - j5
Tabla # 1.- Impedancias y Admitancias capacitivas para el ejemplo 1.

GENERACION CARGA
NODO VOLTAJE MW's MVAR's MW's MVAR's Tipo Nodo
1 1.05+J0 -------- -------- 0 0 Slack
2 1.03+J0 20 -------- 50 20 Controlado
3 1+J0 0 0 60 25 Carga
Tabla # 2.- Voltajes Iniciales, Potencias de Generación y de carga para el Ejemplo # 1.
La matriz de admitancias nodal YBUS  , para el ejemplo #1 se presenta a continuación:
(6.25  j18.75) (1.25  j3.75) (5  j15) 

YBUS    (2.916  j8.75) (1.666  j5) P.U
 (6.666  j 20) 
Los cálculos de los elementos YL Pq , y KLP se realizan utilizando las ecuaciones (8) y (9).
(1.25  j 3.75)
YL12   (0.2  j 0)
(6.25  j18.75)
(1.25  j 3.75)
YL 21   (0.4285714  j 0)
(2.916  j8.75)
(1.666  j 5)
YL 23   (0.5714285  j 0)
(2.916  j8.75)
(1.666  j 5)
YL32   (0.25  j 0)
(6.666  j 20)
(5  j15)
YL31   (0.75  j 0)
6.666  j 20)
P  jQ3 (0.6  j 0.25)
KL3  3   (0.02025  j 0.02325)
Y33 (6.666  j 20)

Se debe tener presente que la constante KL para el nodo 2 no se calcula ya que es un nodo de voltaje controlado ( PV ).
Las ecuaciones que se utilizan para la solución del problema de flujos de carga son las siguientes:
KL2
V2( K 1)  *(K )
 YL21V1  YL23V3( K ) (18)
V2
KL
V3( K 1)  *(K3)  YL31V1  YL32V2( K 1) (19)
V3
Por otra parte la ecuación del nodo 1 no se requiere, ya que es el nodo compensador y su voltaje y ángulo se conocen.
La primera ecuación a evaluar es la del nodo 2 que es un nodo generador, tomando un valor de 0  para el ángulo en la
primera iteración,  (1)  0  . Por lo tanto:
(e( 2) nuevo ) (1)  1.03 cos 0   1.03(1)  1.03
( f ( 2) nuevo ) (1)  1.03sen0   1.03(0)  0

La potencia reactiva se calcula mediante la ec.(6)

 ELT - 2761 SISTEMAS DE POTENCIA I
AUXILIAR: Univ. Paul Danny Alba Lozada “JHAMBAO”
Ejercicios De Flujos De Carga 2do EXAMEN PARCIAL
S 2*  V2* I 2  P2  jQ 2
P2  jQ 2  (1.03  j 0)(1.25  j 3.75)(1.05  j 0)  (2.916  j8.75)(1.03  j 0)  (1.666  j 5)(1  j 0)
 jQ 2   j 0.07725
o bien
Q2  0.07725
Este valor de potencia reactiva más la de la carga local se encuentra dentro del rango especificado, por lo que se utiliza
para evaluar la constante KL.
(0.2  0.5)  j (0.07725)
KL(21)   (0.002337  j 0.0335)
(2.91  j8.75)
A continuación se evalúa la ecuación del voltaje del nodo 2, para obtener una corrección del ángulo.
(0.002337  j 0.0335)
V2(1)   (0.4285  j 0.0)(1.05  j 0)  (05714  j 0.0)(1  j 0)
(1.03  j 0)
V2(1)  (1.01916  j 0.03253) P.U
V2  (1.011916  1.03)  j (0.03253  0)  (0.01084  j 0.03253) P.U
V2((1Acelerado)  (1.03  j 0)  1.4( 0.01084  j 0.03253)  (1.01482  j 0.04554)
)
P.U
En valor que se requiere utilizar para este caso es una mejor aproximación del ángulo, ya que el voltaje se debe de
mantener en 1.03 P. U. por lo tanto:
V2(1)  (1.02896  j 0.04618) P.U .
V2(1)  1.03  2.56972  P. U .
Para el voltaje en el nodo 3, se usa la ec.(19), obteniendo:
(0.02025  j 0.02325)
V3(1)   (0.75  j 0.0)(1.05  j 0)  (0.25  j 0.0)(1.0289  j 0.04618)
(1  j 0)
V3(1)  (1.02449  j 0.03479) P. U .
Para obtener un valor mejorado, se utiliza la ec.(12). Por lo tanto:
V3((1Acelerado)  (1  j 0)  1.4(1.02449  j 0.03479)  (1  j 0)
)

V31( Acelerado)  (1.033429  j 0.04871) P. U .

Para esta iteración los valores de voltaje son:
 (1.05  j 0) 
V (1)
 (1.02896  j 0.044618)
 P. U .
 (1.03429  j 0.04871) 
El problema converge en 7 iteraciones y la solución digital se proporciona a continuación:
 (1.05  j 0) 
V (7)
 (1.028723  j 0.051262) 
 P. U .
(1.024195  j 0.034839)
Una vez que se conocen los voltajes nodales, se procede a calcular los flujos de potencia en las líneas de transmisión.
Por ejemplo supóngase que se desea calcular la inyección de potencia en la línea conectada en los nodos 1 y 2 así como
las pérdidas. Aplicando las siguientes ecuaciones tenemos:
YPq'
 V (VP  Vq )YPq  VP
* * 2
S Pq P (20)
2
similarmente la potencia inyectada por el nodo “q” es:
2 YPq'
S *
qP  V (Vq  VP )YPq  Vq
q
*
(21)
2
 ELT - 2761 SISTEMAS DE POTENCIA I
AUXILIAR: Univ. Paul Danny Alba Lozada “JHAMBAO”
Ejercicios De Flujos De Carga 2do EXAMEN PARCIAL
Para calcular las pérdidas a través de la línea usamos la siguiente ecuación:
S Pq  S qP  ( PP  Pq )  j (QP  Qq ) (22)
De lo anterior y aplicando la ec.(20), se obtiene:
*
S12  (1.05  j 0)(1.05  j 0)  (1.028723  j 0.051262)(1.25  j3.75)
S12  (0.22977  j 0.016494) P. U .
de igual forma aplicando la ec.(21), se obtiene:
*
S 21  (1.02419  j 0.034483)(1.024195  j 0.034839)  (1.05  j 0)(1.25  j3.75)
S 21  (0.22592  j 0.004942) P. U .
Las pérdidas se determinan aplicando la ec.(22), por lo tanto:
S pérdidas(1, 2)  (0.22977  j 0.016494)  (0.22592  j 0.0049423)
S pérdidas(1, 2)  (0.0038507  j 0.011552)
La potencia real y reactiva que suministra el generador conectado al nodo compensador se calcula sumando las
potencias inyectadas en las líneas que inciden en dicho nodo. Para nuestro ejemplo se obtiene:
S G1  S12  S13  (0.22977  j 0.016494)  (0.68419  j 0.2235047)
S G1  (0.913962  j 0.23999) P. U .
En la Tabla # 3 se proporcionan los flujos de potencia en cada una de las líneas. Si se suman las potencias en
columna de la Tabla # 3, se obtendrán las pérdidas totales de potencia activa y reactiva del sistema. Observando la
Tabla # 3 :
Ppérdidastotales  0.01373245 P. U . Q pérdidastotales  0.04119743 P. U .
P-Q P Q
1-2 0.22977 0.016494
1-3 0.68419 0.22350
2-1 -0.22592 -0.004942
2-3 -0.074147 0.055267
3-1 -0.67479 -0.195310
3-2 0.076442 -0.053816
SUMA 0.013732 0.041197
Tabla # 3.- Flujos de Potencia en las Líneas y Perdidas Totales.

G1 G2
24.01
50+j20
91.37 20 25.2

1 2

22.85 22.47
68.
51 1.67 .5 7.5 7
22
.34 5.

67 7
.57
19 7.5 3
.46 5.5

Potencia Real

Potencia Reactiva
60+j25
NOTA: ESPERO QUE LES SIRVA DE MUCHO YA QUE ESTA CLASE DE EJERCICIOS SON MUY LARGOS Y ESTE METODO
ESTA EXPLICADO A DETALLE Y BIEN RESUMIDO CUALQUIER DUDA NO DUDEN EN PREGUNTARME…… “GRACIAS”