Ejercicio 1 (Basado en Ejemplo 9.2 y Ejemplo 9.3 Del Libro Análisis de Sistemas de Potencia de Grainger y Stevenson JR.)

UNIVERSIDAD DE EL SALVADOR
FACULTAD DE INGENIERÍA Y ARQUITECTURA

ESCUELA DE INGENIERÍA ELÉCTRICA
CICLO II - 2018
Autor: Cristian Antonio Aguilar Quintanilla.
Ejercicios de flujo de carga.
Ejercicio 1 (Basado en Ejemplo 9.2 y Ejemplo 9.3 del libro Análisis de Sistemas de Potencia
de Grainger y Stevenson Jr.). Considerar el siguiente sistema de cuatro barras, la potencia
trifásica base es de 100 MVA, y la tensión línea-línea base es de 230 kV.
Se dan las impedancias serie y las susceptancias en pu para los circuitos equivalentes
nominales π de los cálculos de las cuatro líneas identificadas que terminan en las barras.
1
𝑍̅𝑠𝑒𝑟𝑖𝑒 𝑌̅𝑠𝑒𝑟𝑖𝑒 = Admitancias en paralelo
𝑍̅𝑠𝑒𝑟𝑖𝑒
Línea, de barra a barra
𝑌̅
Rserie Xserie Gserie Bserie MVAR totales de carga
2
1-2 0.01008 0.05040 3.815629 -19.078144 10.25 0.05125
1-3 0.00744 0.03720 5.169561 -25.847809 7.75 0.03875
2-4 0.00744 0.03720 5.169561 -25.847809 7.75 0.03875
3-4 0.01272 0.06360 3.023705 -15.118528 12.75 0.06375
𝑌̅ 2
Nota: Para calcular se hizo así, 𝑄 = √3𝑉𝐿𝐿 𝐼𝑐𝑎𝑟𝑔𝑎 e 𝐼𝑐𝑎𝑟𝑔𝑎 = 𝑉𝐿𝐿 𝑌𝑐𝑎𝑟𝑔𝑎 /√3, se obtiene que 𝑄 = 𝑉𝐿𝐿 𝑌𝑐𝑎𝑟𝑔𝑎 .
2
2 −4
Para el primer dato 𝑌𝑐𝑎𝑟𝑔𝑎 = (10.25 MVAR)/(230 kV) = 1.9376 × 10 S, la impedancia base del sistema
2
será: 𝑍𝑏𝑎𝑠𝑒 = 𝑉𝐿𝐿 𝑏𝑎𝑠𝑒
/𝑆3𝜙𝑏𝑎𝑠𝑒 = (230 kV)2 /(100 MVA) = 529 Ω, por lo cual 𝑌𝑐𝑎𝑟𝑔𝑎 [pu] = 𝑌𝑐𝑎𝑟𝑔𝑎 × 𝑍𝑏𝑎𝑠𝑒 =
0.1025 pu, la admitancia en derivación será: 𝑌𝑑𝑒𝑟𝑖𝑣𝑎𝑐𝑖ó𝑛 [pu] = 𝑌𝑐𝑎𝑟𝑔𝑎 [pu]/2 = 0.05125 pu.
En la siguiente tabla se enlista los valores de P, Q y V para cada barra.
Generación Carga Tensión

Barra
P [MW] Q [MVAR] P [MW] Q [MVAR] V [pu] δ [grados]
1 ----- ----- 50 30.99 1.00 0°
2 0 0 170 105.35 1.00 0°
3 0 0 200 123.94 1.00 0°
4 318 ----- 80 49.58 1.02 0°
Nota: Las potencias reactivas de carga en las barras, están calculadas a partir del valor de las potencias activas,
considerando un factor de potencia de 0.85 en atraso.
Como las barras 2 y 3 no tienen generación, son barras de carga pura (barras PQ). La barra 4
posee carga y generación, pero no se especifica la potencia reactiva generada (es una
variable), así que es una barra de tensión controlada (barra PV, o sea tensión fija en magnitud,
pero con ángulo de tensión de barra ajustable). La barra 1, por no conocerse las potencias
generadas, se asigna como barra de referencia.
Realizar tres iteraciones para encontrara 𝑉̅2, 𝑉̅3, 𝛿4 , 𝑄𝐺4 , 𝑃𝐺1 y 𝑄𝐺1 , mediante:
a) Método de Gauss-Seidel.
b) Método de Gauss-Seidel con sobrerelajación de 1.6 (aceleración de convergencia).
Para el literal a), completar todo el algoritmo de Gauss-Seidel, para el literal b) se obtendrán
iguales resultados. Algunos aplican la aceleración de convergencia al terminar de calcular los
valores no acelerados de una iteración, esa forma no es matemáticamente correcta, aunque
da resultados satisfactorios, pero se requieren más iteraciones; lo mejor es acelerar cada valor
calculado y usarlo en la siguiente variable de estado a calcular.
Solución. Desde el paso 3, se establece la diferencia entre el método sin aceleración y con
aceleración. Los valores que no tengan unidades están en pu.
Paso 1. Establecer las inyecciones iniciales 𝑃𝑖 + 𝑗𝑄𝑖 en todas las barras del sistema:
No se puede determinar aun

[𝑃𝐺 − 𝑃𝐷 ] −1.7
[𝑃𝑖0 ] = 𝑖 𝑖
=[ ]
𝑆3𝜙𝑏𝑎𝑠𝑒 −2
2.38
[𝑄𝐺 − 𝑄𝐷 ] −1.054
[𝑄𝑖0 ] = 𝑖 𝑖
=[ ]
𝑆3𝜙𝑏𝑎𝑠𝑒 −1.239
Paso 2. Construcción de la matriz de admitancias de barra.
𝑦̅12 + 𝑦̅13 + 𝑦̅120 + 𝑦̅130 −𝑦̅12 −𝑦̅13 0

−𝑦̅21 𝑦̅21 + 𝑦̅24 + 𝑦̅210 + 𝑦̅240 0 −𝑦̅24
[𝑌̅𝑏𝑎𝑟𝑟𝑎 ] = [ ]
−𝑦̅31 0 𝑦̅31 + 𝑦̅34 + 𝑦̅310 + 𝑦̅340 −𝑦̅34
0 −𝑦̅42 −𝑦̅43 𝑦̅42 + 𝑦̅43 + 𝑦̅420 + 𝑦̅430
8.985 − 𝑗44.836 −3.816 + 𝑗19.078 −5.17 + 𝑗25.848 0
−3.816 + 𝑗19.078 8.985 − 𝑗44.836 0 −5.17 + 𝑗25.848
[𝑌̅𝑏𝑎𝑟𝑟𝑎 ] = [ ]
−5.17 + 𝑗25.848 0 8.193 − 𝑗40.864 −3.024 + 𝑗15.119
0 −5.17 + 𝑗25.848 −3.024 + 𝑗15.119 8.193 − 𝑗40.864
Paso 3. Literal a), Método Gauss-Seidel sin aceleración de convergencia.
Primera iteración.
Segunda iteración.
Tercera iteración.
Para seguir con los siguientes pasos, se realizó hasta la iteración 14, obteniendo:
Paso 4. Calculo de la potencia inyectada en la barra de referencia y de las potencias generadas
y demandadas de la misma y la generada en la barra de tensión controlada, luego de la
iteración 14.
Paso 5. Computo de los flujos de línea.

Las perdidas en las líneas serán:
Simulación en Power World
Paso 3. Literal b), Método Gauss-Seidel con aceleración de convergencia.
Primera iteración.
Segunda iteración.
Tercera iteración.
Ejercicio 2. Considerar el siguiente sistema de tres barras. Se presenta una tabla con la
potencia generada y/o demandada en cada barra. Aplicar el método de Gauss-Seidel con un
factor de aceleración de 1.4 y una tolerancia global de 0.02, encontrar las tensiones en las
barras 2 y 3 y resolver el flujo de carga, para las dos condiciones mostradas.
Líneas, de barra a barra ̅

𝑍𝑠𝑒𝑟𝑖𝑒 en pu 𝑌̅/2 en pu
1-2 j0.01 0
2-3 j0.02 0
1-3 j0.02 0
Condición de carga 1:
Barra Generación Carga Tensión

P [pu] Q [pu] P [pu] Q [pu] V [pu] δ [grados]
1 ----- ----- 0 0 1.0 0°
2 1.2 0.74 0 0 1.0 0°
3 0 0 3.64 2.5 1.0 0°
Barra Generación Carga Tensión Límites de potencia

P [pu] Q [pu] P [pu] Q [pu] V [pu] δ [grados] QGmín [pu] QGmáx [pu]
1 ----- ----- 0 0 1.03 0.2° ----- -----
2 1.2 ------ 0 0 1.03 0° 0.87 1.6
3 0 0 3.64 2.5 1.0 0° ----- -----
Solución. Para la condición de carga 1, la barra 1 es la de referencia, la barra 3 es de carga y

la barra 2 al especificar sus potencias generadas (tanto activa como reactiva), su tensión de
barra en magnitud no es fija, se considera como una barra de carga (por condiciones de
clasificación). Para la condición de carga 2, ahora la barra 2 es una barra de tensión controlada
con restricciones en los niveles de potencia reactiva generada, si esos límites se violan, es
necesario restringir el valor de la potencia limitada a un valor dado por lo límites y tratar la
barra como PQ, y como no se tienen restricciones para tensión en esa barra no se puede
regresar a ser barra PV.
Paso 1. Inyecciones en barra.
Paso 2. Construcción de matriz de admitancia de barra.
Paso 3.
Primera iteración.
Segunda iteración.
Tercera iteración.
Cuarta iteración.
Se detiene, porque se cumple el requerimiento de tolerancia para ambas tensiones.

Paso 4. Determinación de potencia en la barra de referencia.
Simulación en Power World suponiendo niveles de tensión de 138 kV, 100 MVA.
Paso 1.
Paso 2. La misma matriz de la condición 1.

Paso 3.
Primera iteración.
Segunda iteración.
Tercera iteración.
Cuarta iteración. Se podría haber parado al final de la iteración anterior, pero para disminuir
el error se hizo una más.
Detener la iteración.
Paso 4. Calculo de la potencia inyectada en la barra de referencia.
Nota: En caso de violar la potencia reactiva, hacer 𝑄2 = 𝑄𝐺2𝑚á𝑥 − 𝑄𝐷2, y la barra será PQ
ahora. Ya que no se debe mantener la tensión constante en la barra 2 y como no se tienen
límites de tensión en dicha barra, esta se mantendrá como PQ, o sea no hay posibilidad que
vuelva a ser PV. Si se da el caso que la barra se vuelve PQ, se puede establecer 𝑉2 = 1 pu,
en lugar del valor dado cuando se consideraba PV, pero no es necesario, ya que el mismo
proceso iterativo lo corregirá.
Ejercicio 3. (Basado en el problema 6.36 de libro de Duncan). La matriz de admitancia de
barra en pu del sistema de potencia mostrada en la figura es la siguiente:
3 − 𝑗9 −2 + 𝑗6 −1 + 𝑗3 0
−2 + 𝑗6 3.666 − 𝑗11 −0.666 + 𝑗2 −1 + 𝑗3
[𝑌̅𝑏𝑎𝑟𝑟𝑎 ] = [ ]
−1 + 𝑗3 −0.666 + 𝑗2 3.666 − 𝑗11 −2 + 𝑗6
0 −1 + 𝑗3 −2 + 𝑗6 3 − 𝑗9
Las potencias complejas se muestran en las barras 2, 3 y 4 en la misma figura. Determinar

los valores de las tensiones en dichas barras producidas para la primera y segunda iteración,
usando los métodos:
a) Gauss-Seidel.
b) Newton-Raphson.
Asumir tensiones iniciales constantes de 1.0∠0° pu en las barras 2, 3 y 4.
Solución. Las barras 2, 3 y 4 son barras de carga pura (no tienen generadores conectados, o
en el caso de las barra 2 y 4 su potencia activa inyectada es positiva Pi=Pgi-Pdi, por lo que
podrían ser generadores que envían potencia activa constante y consumen potencia reactiva
constante por eso sus inyecciones son negativas, o sea son generadores síncronos
subexcitados, con los circuitos de campo y eje mecánico ajustados para esos valores
mostrados), la barra 3 puede ser una simple carga RC.
a) Método de Gauss-Seidel. Paso 2 se omite, porque ya dan la matriz de admitancia de
barra.
Paso 1.
Paso 3.
Primera iteración.
Segunda iteración.
El libro de Kothari, da un procedimiento para hacer los cálculos de una forma “mas simple”,
haciéndolos por partes:
Primera iteración.
Segunda iteración.
b) Método de Newton-Raphson.
Paso 1. Establecer las tensiones iniciales y calcular la matriz de admitancia de barra.

Primera iteración.
Paso 2. Calculo de los 𝛥𝑃𝑖 para las barras PQ y PV, y 𝛥𝑄𝑖 para las barras PQ, a partir de las
ecuaciones de balance de potencia (𝑔𝑖’ = 𝑃𝑖 − 𝑃𝑖𝑝𝑟𝑜𝑔 = 𝑃𝑖𝑐𝑎𝑙𝑐 − 𝑃𝑖𝑆𝑃 , 𝑔𝑖’’ = 𝑄𝑖 − 𝑄𝑖𝑝𝑟𝑜𝑔 =
𝑄𝑖𝑐𝑎𝑙𝑐 − 𝑄𝑖𝑆𝑃 ) 𝛥𝑃𝑖 = 𝑃𝑖𝑝𝑟𝑜𝑔 − 𝑃𝑖 , 𝛥𝑄𝑖 = 𝑄𝑖𝑝𝑟𝑜𝑔 − 𝑃𝑖 .
Las potencias programadas son las del diagrama, las mostradas, que vendrían a ser 𝑃𝑖𝑝𝑟𝑜𝑔 =
𝑃𝐺𝑖 − 𝑃𝐷𝑖 y 𝑄𝑖𝑝𝑟𝑜𝑔 = 𝑄𝐺𝑖 − 𝑄𝐷𝑖 .
Paso 3. Se plantea el vector de residuos y se calcula la matriz jacobiana.
El vector “r” solo se calculo para dejar indicado el orden los términos, la solución se trabaja
con el vector “R”.
Paso 4. Calculo de las correcciones.
Paso 5. Actualización de ángulos y magnitudes de tensión.

Se actualizan los vectores que contienen los ángulos, magnitudes y fasores de tensión en la
primera iteración.
Segunda iteración.
Paso 2.
Paso 3.
Paso 4.
Paso 5.

Ejercicio 1 (Basado en Ejemplo 9.2 y Ejemplo 9.3 Del Libro Análisis de Sistemas de Potencia de Grainger y Stevenson JR.)

Cargado por

Copyright:

Formatos disponibles

Ejercicio 1 (Basado en Ejemplo 9.2 y Ejemplo 9.3 Del Libro Análisis de Sistemas de Potencia de Grainger y Stevenson JR.)

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ejercicio 1 (Basado en Ejemplo 9.2 y Ejemplo 9.3 Del Libro Análisis de Sistemas de Potencia de Grainger y Stevenson JR.)

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD DE EL SALVADOR

FACULTAD DE INGENIERÍA Y ARQUITECTURA

Autor: Cristian Antonio Aguilar Quintanilla.

Ejercicios de flujo de carga.

En la siguiente tabla se enlista los valores de P, Q y V para cada barra.

Generación Carga Tensión

No se puede determinar aun

Paso 2. Construcción de la matriz de admitancias de barra.

𝑦̅12 + 𝑦̅13 + 𝑦̅120 + 𝑦̅130 −𝑦̅12 −𝑦̅13 0

Paso 5. Computo de los flujos de línea.

Líneas, de barra a barra ̅

Barra Generación Carga Tensión

Barra Generación Carga Tensión Límites de potencia

Solución. Para la condición de carga 1, la barra 1 es la de referencia, la barra 3 es de carga y

Paso 1. Inyecciones en barra.

Paso 2. Construcción de matriz de admitancia de barra.

Se detiene, porque se cumple el requerimiento de tolerancia para ambas tensiones.

Paso 5. Computo de los flujos de línea.

Paso 2. La misma matriz de la condición 1.

Paso 5. Computo de los flujos de línea.

Las potencias complejas se muestran en las barras 2, 3 y 4 en la misma figura. Determinar

Asumir tensiones iniciales constantes de 1.0∠0° pu en las barras 2, 3 y 4.

Paso 1. Establecer las tensiones iniciales y calcular la matriz de admitancia de barra.

Paso 5. Actualización de ángulos y magnitudes de tensión.

También podría gustarte

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.