Flujo de Potencia

Denotemos todas las variables en notación polar:

Vi = |Vi| < δi.
Yik = |Yik| < θik.
Además:
SGi : potencia generada en barra i. PGi + jQGi.
SLi : potencia consumida en barra i. PLi + jQLi.
SGi – SLi : potencia neta inyectada en la barra i.

EL4001 - Conversión de la Energía y Sistemas Eléctricos Eduardo Zamora D.

 Flujo de Potencia
n
Si = SGi − SLi = Vi ⋅ I = Vi ⋅ ∑ Yik* ⋅ Vk*
*
i
k =1
Conocidas
Si = ∑ ( Vi < δ i )⋅ ( Vk < −δ k )⋅ ( Yik < −θ ik )
n
las tensiones
k =1 en todas las
n barras es
Si = ∑ Vi ⋅ Vk ⋅ Yik < (δ i − δ k − θ ik ) posible
k =1
determinar
Entonces los flujos de
n potencia en
Pi = PGi − PLi = ∑ Vi ⋅ Vk ⋅ Yik ⋅ cos(δ i − δ k − θ ik ) todo el
k =1 sistema.
n
Qi = Q Gi − Q Li = ∑ Vi ⋅ Vk ⋅ Yik ⋅ sen (δ i − δ k − θ ik )
k =1
EL4001 - Conversión de la Energía y Sistemas Eléctricos Eduardo Zamora D.
 Clasificación de Barras
En la práctica no es posible fijar la tensión de todas las barras
del sistema.
Habrá consumos o generación que podrán modelarse como de
potencia aparente constante, independiente del voltaje.
En otras barras puede requerirse cierto nivel de tensión fijo para
alimentar determinada potencia activa.
Por último, se requiere de una barra de tensión y ángulo fijo,
que sirva de referencia para los ángulos de las tensiones.
Así, tendremos 3 tipos de barras, con 2 variables fijas y otras 2
a despejar según el modelo:

EL4001 - Conversión de la Energía y Sistemas Eléctricos Eduardo Zamora D.

 Clasificación de Barras
Barra Slack: Es la “barra infinita”. Su tensión y su ángulo (0°) son
fijos independiente de los flujos de potencia. Es única en el sistema y
se asocia a barras “grandes”. Se desconoce P y Q.
Barra PV: Aquella en la que se conoce la potencia inyectada P y
cuya tensión V es fija. Se puede tratar de barras cercanas a
generadores que puedan controlar tensión, o a barras que por
requerimiento de operación deban tener tensión fija. Se desconoce Q
y δ.
Barra PQ: Se conoce P y Q inyectados. En general se trata de barras
de consumo a potencia constante. También incluye barras de
“pasada” que tienen P = Q = 0.

EL4001 - Conversión de la Energía y Sistemas Eléctricos Eduardo Zamora D.

 Método de Gauss - Seidel
Queremos encontrar V, δ para todas las barras del sistema. Por cada
barra:
- 6 variables (P,Q,V,δ,I,φ)
- 2 ecuaciones derivadas de I = Y V (Re, Im)
- 2 ecuaciones derivadas de S = V I* (Re, Im)
- 2 variables conocidas
Como se trata de un sistema no lineal que involucra potencias, corrientes
y tensiones, se utilizan métodos iterativos para su resolución:
- Gauss – Seidel
- Newton - Raphson
EL4001 - Conversión de la Energía y Sistemas Eléctricos Eduardo Zamora D.
 Método de Gauss - Seidel
La “Receta” para Gauss – Seidel
1. Designar las barras del sistema: PQ, PV, Slack, y obtener los
valores de sus variables conocidas.
2. Definir un punto de origen para el método iterativo (voltajes
iniciales).
1. Si la barra es PQ, V0 = 1<0° [p.u.]
2. Si la barra es PV, V0 = V<0° [p.u.]
3. Si la barra es Slack, V = V<0° [p.u.] (fijo)
Esto corresponde a una recomendación. Salvo la barra Slack, las
tensiones de barra PQ y ángulo de barra PV son arbitrarios, pero
la convergencia del método requiere partir “cerca de la solución”.
EL4001 - Conversión de la Energía y Sistemas Eléctricos Eduardo Zamora D.
 Método de Gauss - Seidel
Resolvemos en orden para cada barra la iteración correspondiente:
Barras PQ:
PQ:
Estimamos la corriente Ii usando los voltajes de la iteración
anterior (m). *
 Pi + jQ i 
Ii =  m


 Vi 
Como: n
Ii = ∑ ik k
Y
k =1
⋅ V m

Despejamos Vim+1 para la iteración siguiente:

1  m
I i − ∑ Yik ⋅ Vk 
m +1 m +1 m +1
V i =V i <δ i =
Yii  k ≠i  Eduardo Zamora D.
EL4001 - Conversión de la Energía y Sistemas Eléctricos
 Método de Gauss - Seidel
Barras PV
Estimamos Qi con los datos de la iteración anterior.
 m n * m* 
Q i = ImVi ⋅ ∑ Yik ⋅ Vk 
 k =1 
Estimamos la corriente:
*
 Pi + jQ i 
Ii =  m

 Vi 
Despejamos δi para la siguiente iteración (Vi es dato):
1  m 
I i − ∑ Yik ⋅ Vk  
m +1 m +1
<V i =δ i =< 
 Yii  k ≠i 
EL4001 - Conversión de la Energía y Sistemas Eléctricos Eduardo Zamora D.
 Método de Gauss - Seidel
Barra Slack
No requiere entrar al proceso de iteración, ya que sus valores de
tensión y ángulo son fijos. Sin embargo, de todas formas se utiliza
este valor para los cálculos en las otras barras.

El método converge relativamente rápido (orden de 10

iteraciones), aunque claramente el tiempo dependerá de la
precisión que se quiera obtener.

EL4001 - Conversión de la Energía y Sistemas Eléctricos Eduardo Zamora D.

 Método de Gauss - Seidel
Variantes del método
Gauss – Seidel con actualización de variables
- En cada cálculo utilizo los valores más recientes de tensión, es
decir, para calcular Vim+1 utilizo V1m+1, ...Vi-1m+1, Vim, ...Vnm.

Gauss – Seidel con factores de aceleración

- Si llamamos Viom+1 al valor calculado con el método tradicional,
definimos ∆V = Viom+1 –Vim, entonces, Vim+1 = Vim + α∆V, donde α
se escoge usualmente entre 1.5 y 1.8

EL4001 - Conversión de la Energía y Sistemas Eléctricos Eduardo Zamora D.

 Método de Gauss - Seidel
Observaciones:
- Si en una barra PV asociada a un generador se obtiene Q superior
a Qmax dado por la carta de operación, esta barra debe cambiarse a
PQ, con Q = Qmax. Esto se debe a que la central pierde capacidad
de regular tensión al verse limitada en su entrega de reactivos.
- El flujo de potencia activa es desde barras de ángulos mayores
hacia barras de ángulos menores.

EL4001 - Conversión de la Energía y Sistemas Eléctricos Eduardo Zamora D.