Flujos de potencia en una línea de transmisión

Estudio y modelado en powerworld de los flujos de potencia en las líneas de transmisión

Diego Yahir Ortega Figueroa

Universidad de Guadalajara
CUCEI
Guadalajara, México
diego.ortega5689@[Link]

Resumen—En el presente documento se estudiará acerca de Primeramente, hay que recordar los siguientes conceptos y
los flujos de potencia en los sistemas de transmisión, también se ecuaciones básicas para el desarrollo del problema:
hablará acerca de las formulas que intervienen en él, así como el ¿
modelado en un programa de computación para obtener los Si=V i I i (1)
resultados.
n

Términos clave—admitancia, bus, consumo, flujo de potencia, I i=∑ Y ij V j (2)

generador, impedancia, nodo i=1

V
I. INTRODUCCIÓN I= (3)
Z
Para tener una operación exitosa de los sistemas de potencia
en condiciones normales balanceadas de estado estable 1 (4)
Y=
trifásico, se requiere lo siguiente: Z
I. La generación abastece la demanda (carga) más las I =YV (5)
pérdidas.
Y ij =Gij +B ij j
II. Las magnitudes de voltaje en las barras permanecen
cercanas a sus valores nominales. n (6)
III. Los generadores operan dentro de límites
Si=S Gi+ S Ci=V i ∑ (Y ¿ ¿ij V ¿ ¿ j)¿ ¿ ¿
j=1
especificificados de potencia real y reactiva.
n (7)
IV. Las líneas de transmisión y los transformadores no Si=Pi +Qi j=V i ∑ [ ( Gij + Bij j ) V j ¿ ¿ ¿ ¿ ]
están sobrecargados. j=1

Los programas de computadora han sido de gran ayuda Finalmente, para poder trabajar de una manera más
como una herramienta básica para investigar los requerimientos cómoda, se expresarán las fórmulas de la siguiente manera:
de los sistemas de potencia. Se calcula la magnitud del voltaje n
y el ángulo en cada barra o bus en un sistema de potencia en
condiciones balanceadas de estado estable trifásico. También Pi=V i ∑ V j (Gij cos θ ij + Bij sen θij ) (8)
j=1
permite calcular los flujos de potencia real y activa para los
equipos que interconectan las barras, así como las pérdidas en n
los equipos. Son de interés tanto los elementos del sistema de Qi=V i ∑ V j (Gij sen θij −Bij cos θ ij ) (9)
potencia existente como los cambios propuestos, incluyendo la j=1
nueva generación y transmisión para satisfacer el crecimiento Del mismo modo, para conocer los componentes de la
de la carga proyectada. sumatoria G ij y Bij se realizará una matriz de admitancias de
El análisis nodal común, o de mallas, no es adecuado para nuestro sistema, también conocida como “Y de Bus”, para el
estudios de flujos de potencia porque los datos de entrada para cuál necesitamos conocer las impedancias de nuestras líneas de
las cargas por lo general se dan en términos de potencia, no de transmisión. Estas admitancias se calculan con la ecuación (4),
impedancia; así mismo, se considera a los generadores como dejando claro que la admitancia es el inverso de la impedancia.
fuentes de potencia, no fuentes de voltaje o corriente. Por lo Para esta matriz se realiza un análisis nodal como comúnmente
tanto, el caso de estudio se formula como un conjunto de se realizaría, tomando las admitancias que intervienen en cada
ecuaciones algebraicas no lineales adecuado para resolverlo por nodo y entre nodos.
computadora.
El análisis matemático que se va a realizar en este
documento aborda hasta el planteamiento de las ecuaciones a
solucionar para encontrar las incógnitas de nuestro sistema.

1
Fig. 1. Ejemplo de matriz de admitancias [1]

Por último, necesitamos conocer el tipo de nodos que

podemos encontrar en los sistemas, esto para conocer los datos
e incógnitas que nos proporcionan.
Fig. 2. Ejemplo de diagrama de un sistema de potencia a una línea [2]
Nodos de consumo: En estos nodos se conoce el consumo
de potencia activa y reactiva, siendo nula la potencia activa
generada. Quedando como incógnitas los voltajes nodales, así II. PLANTEAMIENTO DEL PROBLEMA
como su ángulo. Para este proyecto se analizará una red eléctrica en dos
Nodos de generación: Son los nodos donde se conoce la casos distintos, esta red se encuentra distribuida de la siguiente
magnitud de voltaje y la potencia activa generada. Quedando manera y conocemos los siguientes datos de esta misma:
como incógnitas la potencia reactiva y el ángulo del voltaje
nodal.
Nodo de referencia (slack): En este nodo se conocen tanto
el voltaje nodal como el ángulo del mismo. El ángulo se
específica como 0 ya que es una referencia para los demás
ángulos. Este nodo siempre es uno de generación.
La razón común de que sea difícil manejar una gran
transferencia de potencia es que hay pocos mecanismos para
controlar su ruta a través del sistema de transmisión desde el
generador hasta la carga distante.
A menudo esa ruta es indirecta, dependiendo de las
impedancias de las líneas y los lugares donde entra o sale la Fig. 3. Topología de la red eléctrica.
potencia del sistema.
Para determinar cómo se mueve la potencia a través de una
TABLE I. IMPEDANCIA DE LAS LÍNEAS
red de transmisión de generadores a cargas, es necesario
calcular el flujo de potencia real y reactiva en cada una de las Impedancia de las líneas
líneas de transmisión o transformadores, juntos con los voltajes Línea Resistencia Reactancia
de bus asociados.
1-10 0 0.0167
El problema de la visualización es hacer intuitivos estos
2-5 0 0.0167
conceptos. Una técnica simple pero efectiva para ilustrar el
flujo de potencia en una red eléctrica es usar el flujo animado 3-4 0 0.005
en líneas. 3-5 0.001 0.01
Para este proyecto se hará uso del programa computacional 3-6 (ambas) 0.011 0.11
Power-World, en este programa podemos interpretar de manera
10-5 0.0025 0.025
sencilla el comportamiento del flujo de potencia, el tamaño,
orientación y velocidad de las flechas indican la dirección del 11-9 0 0.0167
flujo de potencia en la línea. 12-8 0 0.0167
El programa nos proporciona gráficas de sectores 13-6 (ambas) 0.011 0.11
dimensionadas dinámicamente, la cual es otra forma útil de
13-7 0 0.005
visualización probada para para la rápida detección de
sobrecargas en una red grande. En el diagrama unifilar, el 13-8 0.001 0.01
porcentaje que ocupa cada sector en una gráfica indica qué tan 9-8 0.0025 0.025
cerca está cada línea de transmisión a su límite térmico.
TABLE II. DATOS DE LOS BUSES LOCALIZADOS EN LA RED

2
 Datos de los buses
Barra
No.
Tipo V (PU) Pg Qg Pc Qc
2 PV 1 7 - - -
3 PV 1 7.16 - - -
4 PQ 1 7 - 9.76 4.25
7 PQ 1 - - 17.65 7.2
1 SL 1 - - - -

III. CÁLCULOS Y PROCEDIMIENTO PARA LA SOLUCIÓN

Fig. 4. Flujo de potencia para el caso 1.
Primeramente, se calcula la Y de bus a partir de las
impedancias de las líneas en la red.
TABLE III. RESULTADOS PARA EL CASO 1
Por fines prácticos, no se realizará el cálculo de la matriz
manualmente, pero se obtiene mediante el procedimiento Buses
mencionado anteriormente en la introducción. Esta Y de bus
será calculada a partir de la simulación de la red en el programa No. de Pc Qc Pg Qg
bus
V (PU) θ (MVAR
powerworld. Sus resultados son presentados en el “anexo 1” al (MW)
)
(MW) (MVAR)
final del documento.
1 1 0 42.55 -298.96
Ahora, por medio de las ecuaciones (8) y (9) y el análisis de
2 1 -1.28 1 -741.84
la tabla de los datos de los buses, se empiezan a plantear las
ecuaciones para nuestra red. 3 1.2277 -1.94 1 1524.2

P2=7−0=7 (10) 4 1.2276 -1.96 9.76 4.2

(11) 5 1.1238 -1.29

P3=7.16−0=7.16
6 1.1190 -1.89
(12)
P4 =7−9.76=−2.76 7 1.0051 -1.59 17.65 7.2
(13)
Q4 =0−4.25=−4.25 (14)
8 0.9999 -2.32
9 1 -2.31
P7=0−17.65=−17.65 (15)
10 1.0499 -0.39
Q7=0−7.2=−7.2 (16) 11 1 -2.3 1 0

P2=|V 2|[|V 2|(G22)+|V 5|(G25 cos θ 25+ B25 sen θ 25)](17) 12 1 -2.31 1 0.15
13 1.0102 -1.82
P3=|V 3|[|V 3|(G33 )+|V 4|(G34 cos θ 34 +B 34 sen θ34 )+|V 5|(G35 cos θ35 +B 35 sen θ35)+|V 66|(G36 cos θ 36+ B36 sen θ36)]
TABLE IV. PÉRDIDAS EN LAS LÍNEAS PARA EL CASO 1
P4 =|V 4|[|V 3|(G43 cos θ 43+ B43 sen θ 43)+|V 4|(G 44)]
(18) Línea de transmisión
Q4 =|V 4|[|V 3|(G 43 sen θ 43−B43 cos θ43 )+|V 4|(−B 44 )] Pérdida de Q en
(19) Bus Pérdida de P (MW)
(MVAR)
P7=|V 7|[|V 7|(G77 )+|V 13|(G713 cos θ713 + B713 sen θ(20)
713 )] 1-10 0 15.23

Q7=|V 7|[|V 7|(−B77 )+|V 13|(G713 sen θ713 −B713 cos(21)

θ 713 )] 2-5 0 91.9
3-4 0 0
Una vez planteadas las ecuaciones, se podría continuar
manualmente mediante el método de Newton-Raphson, pero 5-3 10.86 108.65
por fines de precisión y de una visualización gráfica más 3-6 1.06 10.65
adecuada se realizará por medio de powerworld.
3-6 1.06 10.65
Para el primero de los casos, se realizará la simulación de 10-5 2.28 22.8
manera ordinaria, con los datos presentados anteriormente.
Enseguida se presenta el diagrama y el comportamiento del 6-13 1.06 10.65
flujo de potencia para el caso 1: 6-13 1.06 10.65
13-7 0 0.84
 9-8 0 0 12 1 -2.45 1 0.11
12-8 0 0 13 1.01028 -1.95
13-8 0.18 1.82
11-9 0 0 TABLE VI. PÉRDIDAS EN LAS LÍNEAS PARA EL CASO 2
Línea de transmisión
Pérdida de Q en
Para el segundo caso, se tiene un 50% más en los valores del Bus Pérdida de P (MW)
(MVAR)
nodo 4, dando los siguientes resultados en el programa:
1-10 0 15.33
2-5 0 91.91
3-4 0 0.01
5-3 10.88 108.77
3-6 1.06 10.65
3-6 1.06 10.65
10-5 2.29 22.95
6-13 1.06 10.65
6-13 1.06 10.65
13-7 0 0.84
9-8 0 0
Fig. 5. Flujo de potencia para el caso 1.
12-8 0 0
13-8 0.18 1.82
TABLE V. RESULTADOS PARA EL CASO 2
11-9 0 0
Buses

No. Pc Qc Pg Qg
de bus
V (PU) θ (MW)
(MVAR
(MW) (MVAR)
)
1 1 0 47.82 -299.27
2 1 -1.39 1 -742.12
3 1.22785 -2.08 1 1527.2 IV. CONCLUSIONES
4 1.22761 -2.11 15 6 El uso de programas computacionales es de gran utilidad
5 1.12394 -1.40 para el análisis de los flujos de potencia en las redes eléctricas,
se debería de implementar más el trabajar con estos softwares
6 1.11907 -2.02
ya que mejora en gran medida la comprensión de los temas con
7 1.00515 -1.73 17.65 7.2 los que se trabajan durante la carrera.
8 0.99999 -2.46
9 1.00001 -2.44
BIBLIOGRAFÍA
10 1.05001 -0.44
[1] A. Loaiza, “Analisis de redes”, (s.f).
11 1 -2.43 1 -0.02 [2] F. Pietryga, “Animation Of A Power System Using Powerworld
Simulator”, 2006.