UNIVERSIDAD DEL

BIOBÍO
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA
ELÉCTRICA Y ELECTRÓNICA

Sistemas Eléctricos de Potencia II

Tarea #2
Fallas de Corto Circuito

2020-2
1
2
 Departamento de Ingeniería Eléctrica y Electrónica
Universidad del Biobío, Chile

TABLA DE CONTENIDOS

Contenido
LISTADO DE FIGURAS........................................................................................................................3
PROBLEMA DE FLUJO DE POTENCIAS EN SISTEMA ELÉCTRICO DE POTENCIA..................4
Parte a:.................................................................................................................................................4
Parte b:...............................................................................................................................................13
Parte d:...............................................................................................................................................19
Parte e:...............................................................................................................................................23
Parte f:...............................................................................................................................................25
Parte g:...............................................................................................................................................25
Parte h:...............................................................................................................................................26
Parte i:...............................................................................................................................................29
CONCLUSIONES.................................................................................................................................38

3
Departamento de Ingeniería Eléctrica y Electrónica
Universidad del Biobío, Chile

LISTADO DE FIGURAS
Fig. 1 - Diagrama Unilineal Sistema de Estudio.................................................................................5
Fig. 2 - Diagrama Modelo Transformador con Tap............................................................................10
Fig. 3 - Diagrama Admitancias con tap nominales.............................................................................13
Fig. 4 - Diagrama Admitancias con tap nominales y corrientes de rama............................................17
Fig. 5 - Diagrama Admitancias con tap..............................................................................................21
Fig. 6 - Flujo de Potencias con tap nominales....................................................................................28
Fig. 7 - Flujo de Potencias con ajustes al sistema...............................................................................37

4
 PROBLEMA DE FLUJO DE POTENCIAS EN SISTEMA ELÉCTRICO
DE POTENCIA
Antecedentes: El sistema eléctrico de potencia de la Figura 1 que opera a una frecuencia nominal
de 60 Hz, con parámetros definidos en la tarea se procede a desarrollar el cálculo necesario para
llevar a cabo el cálculo de cortocircuito del sistema en todas sus barras, fallas simétricas y
asimétricas. Tomando una potencia base de 100 MVA, voltaje base línea a línea de 230 kV.
Fig. 1 – Diagrama Unilineal Sistema de Estudio.

Parte 0:
Se procede a realizar el cálculo en por unidad en bases del sistema para obtener los
parámetros necesarios para construir las redes de secuencia del sistema eléctrico.
Se implementa un cálculo de admitancias de todos los elementos que componen el Sistema de Estudio
y luego se calculan sus admitancias en pu con los valores bases impuestos al sistema S base = 100 MVA y
Vbase = 230 kV.
Cálculo de Parámetros de Línea 4-5 con espaciamiento simétrico:
Primero, para calcular la inductancia L de la Línea 4-5 se tienen las siguientes ecuaciones
DMG=√3 5∗5∗5=5[m] (a1.1)

5
 −1
4
(a1.2)
r ' Eficaz=r Real∗e =0.019081[m]

DMG H (a1.3)
L4 −5 =2∗10−7 ln
( '
r Eficaz) [ ]
=1.113704∗10−6
m

H mH (a1.4)
L4 −5 =1.1137∗10−6 [ ] m [ ]
∗106=1.113704
km

Luego, para la capacitancia de la Línea 4-5:

2∗π∗8.854∗10−12 −11 F (a1.5)

C 4−5=
DMG
=1.045992∗10
m [ ]
ln
(
r Real )
F nF (a1.6)
C 4−5=1.0460∗10−11 [ ] m
∗1012=10.459925[ ]
km

Para el cálculo de la resistencia de la línea es necesario multiplicar el valor entregado por catálogo (dado
por el fabricante), por el largo del conductor en este caso 150 km:
Ω (a1.7)
R4 −5=0.252936 [ ]km
∗150[km]=37.9404[Ω]

Ahora ya calculados los parámetros R, L y C, calculamos los paramentos de secuencia (+) y (-) para el
modelo de Línea 4-5:

Z 4−5=R 4−5+ ( j∗2∗π∗60∗L4 −5 )∗150=37.940399+ j 62.978478[Ω]

(a1.8
)

Y 4−5 ( j∗2∗π∗60∗C 4−5 )∗150 −4

= = j 2.957474∗10 [ s ] (a1.9)
2 2

Finalmente procedemos a pasar estos valores a valores propios del sistema por lo que calculamos lo
siguiente:
( (230∗103 )2 ) (a1.10)
Z base = =529[Ω ]
100∗106
6
 Z 4 −5 (a1.11)
Z L 4−5= =0.071720+ j 0.1190519[ pu]
Z base

1 (a1.12)
Y 4 −5= =3.712792− j6.162981[Ω]
Z 4−5

Y 4−5 (a1.13)
Y L4 −5 = ∗Z base = j 0.156450[ pu]
2

Ahora ya calculados los parámetros R, L y C, calculamos los paramentos de secuencia (0) será 3 veces la
secuencia positiva para el modelo de Línea 4-5:

Z 4−5=3∗( R 4−5 + ( j 2 π∗60∗L4−5 )∗150)=113.821197+ j 188.935434[Ω]

(a1.14
)

Y 4−5 3∗( j2 π∗60∗C 4−5 )∗150 −4

= = j 8.872422∗10 [ s] (a1.15)
2 2

Para la línea L4-6 del sistema se procede de la misma manera por lo que se presenta una tabla
resumen con sus resultados obtenidos en el código de matlab Calculo_em_pu.m:
Tabla. 1 – Resultados Parámetros Línea 4-6 Secuencia (+), (-) y (0).
Parámetr Resultado
o
DMG 5 [m]
r’ eficaz 0.012305 [m]
L 1.201436 [mH/km]
C 9.662941 [nF/km]
R 0.089930 [Ω/km]
Z[Ω/km] 0.089930 + j*0.452931 [Ω/km]
Y/2[s/km] j*3.642843 [μs/km]
Z[Ω] 13.489500 + j*67.939642 [Ω]
Y/2[s] j*2.732132*10^-4 [s]
Z[pu] 0.0255 + j*0.1284 [pu]
Y[pu] 1.487349 – j*7.491008[pu]
Y/2[pu] j*0.144530 [pu]
Y(0) [pu] 4.462047-j*22.473024[pu]
Y/2(0) [pu] j*0.43359[pu]
Cálculo de Parámetros de Línea 5-7 con espaciamiento asimétrico:

7
 Primero, para calcular la inductancia L de la Línea 5-7 se tienen las siguientes ecuaciones

DMG=√3 5.8∗5.8∗10=6.9548[m] (a2.1)

−1 (a2.2)
r ' Eficaz=r Real∗e 4 =0.0167[m]

DMG H (a2.3)
L5−7=2∗10−7 ln
( '
)
r Eficaz
=1.2058∗10−6[ ]
m

H mH (a2.4)
L5−7=1.1137∗10−6 [ ]
m
∗106=1.2058
km[ ]
Luego, para la capacitancia de la Línea 5-7:

2∗π∗8.854∗10−12 −12 F (a2.5)

C 5−7=
DMG
=9.626235∗10
m [ ]
ln
(
r Real )
F μF (a2.6)
C 5−7=1.0460∗10−11 [ ]m
∗1012 =9.626235 [ ]
km

Para el cálculo de la resistencia de la línea es necesario multiplicar el valor entregado por catálogo (dado
por el fabricante), por el largo del conductor en este caso 300 km:
Ω (a2.7)
R5−7=0.169280 [ ]
km

Ahora ya calculados los parámetros R, L y C, calculamos los paramentos de secuencia (+) y negativa (-)
para el modelo de Línea 5-7 que corresponde al modelo de una línea larga por lo que se calcula de la
siguiente manera:

Z5 −7 =R5−7 + ( j∗2∗π∗60∗L 4−5 )=0.169280+ j 0.454586[Ω/km]

(a2.8)

Y 5−7=( j∗2∗π∗60∗C 5−7 )= j3.629005∗10−6 [ns /km] (a2.9)

γ 5−7= √ Z5 −7∗Y 5−7=0.000235+ j 0.001308[1 /km]

8
 (a2.10)

(γ ¿¿ 5−7∗300) (a2.11)
Z ' 5 −7 =Z 5−7∗300∗sin =48.297236+ j 133.483741[Ω ]¿
γ 5−7∗300

Y 5−7∗300 γ ∗300
Y'
5−7
=
2
∗tan 5−7
2 ( )
=2.584139 ¿ 10−6 + j5.511729∗10−4 [s ]
(a2.12
)
2 γ 5−7∗300
(2 )
Finalmente procedemos a pasar estos valores a valores propios del sistema por lo que calculamos lo
siguiente:

( (230∗103 )2 ) (a2.13)
Z base = =529[Ω ]
100∗106

Z 5−7 (a2.14)
Z ' L5−7 = =0.091299+ j 0.252332[ pu]
Z base

1 (a2.15)
Y ' L5−7= ∗Z =1.267919− j3.504270[ pu]
Y 5−7 base

Y ' L5−7 Y ' 5−7 (a2.16)

= ∗Z base=0.001367+ j 0.291570[ pu ]
2 2

Ahora ya calculados los parámetros R, L y C, calculamos los paramentos de secuencia (0) para el modelo
de Línea 5-7:

Z5 −7 =3∗(R5−7 + ( j∗2∗π∗60∗L4−5 ) )=0.50784 + j 1.363758[Ω /km]

(a2.14)

Y 5−7=3∗( j∗2∗π∗60∗C5−7 ) = j 10.887015∗10−6 [ns/km] (a2.15)

Z '5 −7 =3∗¿ (a2.16

)

9
 Y 5−7∗300 γ ∗300
Y '
5−7
=3∗(
2
∗tan 5−7
2 ( )
)=7.752417 ¿ 10−6 + j 1.653519∗10−3 [s ]
(a2.17
)
2 γ 5−7∗300
(2 )
Para las demás línea L6-9 del sistema se procede de la misma manera por lo que se presenta una tabla
resumen con sus resultados:
Tabla. 2 – Resultados Parámetros Línea 6-9 Secuencia (+), (-) y (0).
Parámetro Resultado
DMG 6.954811 [m]
r’ eficaz 0.011759 [m]
L 1.276498 [mH/km]
C 9.071566 [nF/km]
R 0.206310 [Ω/km]
Z[Ω/km] 0.206310 + j*0.481229 [Ω/km]
Y/2[ns/km j*3.419899*10^-6 [ns/km]
]
γ [1/km] 0.000269 +j0.001311 [1/km]
Z[Ω] 58.868917 + j*141.471574 [Ω]
Y/2[s] 2.796703*10^-6+j*5.193940*10^-4 [s]
Z[pu] 0.111283 + j*0.267432 [pu]
Y[pu] 1.326319 – j*3.187361 [pu]
Y/2[pu] 0.001479 + j0.274759 [pu]
Z(0)[pu] 0.333849+j*0.802296[pu]
Y(0)[pu] 3.978957-j*9.562083[pu]
Y/2(0)[pu] 0.004437+j*0.824277[pu]
Cálculo de Parámetros de Línea 7-8 con espaciamiento asimétrico:
Primero, para calcular la inductancia L de la Línea 7-8 se tienen las siguientes ecuaciones

DMG=√3 3∗3∗6=3.779763[m] (a3.1)

−1 (a3.2)
r ' Eficaz=r Real∗e 4 =0.024143[m]

DMG −6 H (a3.3)
L7−8=2∗10 ln
−7
( '
)
r Eficaz
=1.010686∗10
m [ ]
H mH (a3.4)
L7−8=1.1137∗10−6 [ ] m
∗106=1.010686
km [ ]
Luego, para la capacitancia de la Línea 7-8:

10
 2∗π∗8.854∗10−12 −12 F (a3.5)
C 7−8=
DMG
=11.581584∗10
m [ ]
ln
(
r Real )
F nF (a3.6)
C 7−8=1.0460∗10−11 [ ] m [ ]
∗1012 =11.581584
km

Para el cálculo de la resistencia de la línea es necesario multiplicar el valor entregado por catálogo (dado
por el fabricante), por el largo del conductor en este caso 500 km:
Ω (a3.7)
R5−7=0.044965 [ ]
km

Ahora ya calculados los parámetros R, L y C, calculamos los paramentos de secuencia (+) y negativa (-)
para el modelo de Línea 7-8 que corresponde al modelo de una línea larga por lo que se calcula de la
siguiente manera:

Z7 −8=R5−7 + ( j∗2∗π∗60∗L 4−5 )=0.044965+ j 0.381019[Ω /km]

(a3.8)

Y 7−8= ( j∗2∗π∗60∗C 5−7 )= j 4.366154∗10−6 [ns /km]

(a3.9)

γ 7−8= √ Z 5−7∗Y 5−7=0.000076+ j 0.001292[1 /km] (a3.10)

(γ ¿¿ 7−8∗300) (a3.
Z ' 7 −8=Z 7−8∗300∗sin =19.461206+ j 177.748936[Ω]¿ 11)
γ 7 −8∗300

Y 7−8∗300 γ ∗300 (a3.12)

Y '
7 −8
=
2
∗tan 7−8
2 ( )
=0.000005+ j 0.001131[ s]
2 γ 7−8∗300
(
2 )
Finalmente procedemos a pasar estos valores a valores propios del sistema por lo que calculamos lo
siguiente:
( (230∗103 )2 ) (a3.13)
Z base = =529[Ω ]
100∗106
11
 Z 7−8 (a3.14)
Z ' L7 −8 = =0.036789+ j 0.336009[ pu]
Z base

1 (a3.15)
Y ' L7−8= ∗Z base=0.321935− j 2.940855[ pu]
Y 7−8

Y ' L7−8 Y ' 7−8 (a3.16)

= ∗Z base=0.002571+ j 0.598292[ pu]
2 2

Ahora ya calculados los parámetros R, L y C, calculamos los paramentos de secuencia (0) para el modelo
de Línea 7-8 que se calcula de la siguiente manera:

Z7 −8=3∗( R ¿ ¿ 5−7+ ( j2 π∗60∗L4−5 ) )=0.134895+ j1.143057 [Ω /km]¿

(a3.17)

Y 7−8=3∗( j∗2∗π∗60∗C5 −7 ) = j 13.098462∗10−6 [ns/km]

(a3.18)

Z ' 7 −8=3∗¿
(a3.19)

γ 7−8∗300 (a3.20)
Y '
7 −8
=3∗(
Y 7−8∗300
tan ( 2 )
)=0.000015+ j 0.003393[s ]
2 2 γ 7−8∗300
( 2 )
Para la línea L8-9 del sistema se procede de la misma manera por lo que se presenta una tabla resumen
con sus resultados:
Tabla. 3 – Resultados Parámetros Línea 8-9 Secuencia (+), (-) y (0).
Parámetr Resultado
o
DMG 3.779763 [m]
r’ eficaz 0.008099 [m]
L 1.229122 [mH/km]
C 9.436057 [nF/km]
R 0.062951 [Ω/km]
Z[Ω/km] 0.062951 + j*0.463368 [Ω/km]
Y/2[s/km] j*3.557310*10^-6 [ns/km]

12
 γ [1/km] 0.000087 +j*0.001286 [1/km]
Z[Ω] 27.283061 + j*216.372325 [Ω]
Y/2[s] 4.513452*10^-6+j*9.211543*10^-4 [s]
Z[pu] 0.051575 + j*0.409021 [pu]
Y[pu] 0.303455 – j*2.406596 [pu]
Y/2[pu] 0.002388 + j*0.487290 [pu]
Z(0)[pu] 0.154725+j*1.227063[pu]
Y(0)[pu] 0.910365-j*7.219788[pu]
Y/2(0)[pu] 0.007164+j*1.46187[pu]
Luego de terminar con las admitancias de las líneas, procedemos a realizar el cálculo de admitancias
de los transformadores que están conectados al sistema de estudio, los parámetros de los
transformadores se encuentran en valores propios por lo que es necesarios evaluarlos en valores
propios del sistema, considerando que los transformadores cuentan con cambiador de tap es
fundamental tomar en cuenta el modelo con tap de la Figura 2. A continuación se procede a realizar
el desarrollo:
Fig. 2 – Diagrama Modelo Transformador con Tap.

Transformador 1 (T1): Primero determinamos la razón de transferencia del transformador:

13.8∗103 (a4
a T 1= = 0.06 .1)
230∗103

Luego se calcula la impedancia propia del transformador referido al secundario con la siguiente
ecuación:
R 1+( j∗2∗pi∗f∗L1 ) (a4.2)
ZT 1 = +( R2 + j∗2∗pi∗f∗L2 )[Ω ]
a2

Reemplazando valores obtenemos:

2.722∗10−3 +( j 2 π f ∗0.1321∗10−3 )
ZT 1 = 2
+(5.2905∗10−3+ j2 π f ∗0.4704 ¿ 10−3 ) (a4.3
a )

Luego pasamos a valores propios del sistema, dividiendo por la Z base obtenemos impedancia de
secuencia (+) y (-):
0.761541+ j 14.010816 (a4
Z ' T 1= =0.001439+ j 0.026485[ pu] .4)
Z base

13
Una vez obtenida la impedancia del transformador en pu, procedemos a calcular la admitancia del
transformador siguiendo los parámetros de un transformador con cambiador de tap, se considera en
este caso en la posición t T1 = 0.95 para el transformador T 1 por lo que se procede de la siguiente
manera:

1 (a4
∗1 .5)
Y T 1 Z 'T 1
= =2.153858− j39.626660[ pu]
tT1 tT 1

Para las admitancias paralelas del transformador tenemos:

(a4
Y T 1∗ ( t−1t )= Z1' ∗( t−1t )=−0.107693+ j 1.981333[ pu ]
T1
.6)

1−t 1 1−t (a4

Y T 1∗
( )
t 2
ZT1 t ( )
= ' ∗ 2 =0.113361− j 2.085614 [ pu] .7)

Impedancia de secuencia (0):

0.761541+ j 14.010816 (a4
Z ' 0T 1=3∗( )=0.004317+ j 0.079455[ pu ] .4)
Z base

Una vez obtenida la impedancia del transformador en pu, procedemos a calcular la admitancia del
transformador de secuencia (0) con cambiador de tap, se considera en este caso en la posición t T1 =
0.95 para el transformador T1 por lo que se procede de la siguiente manera:

1 (a4
∗1 .5)
Y T1 Z 'T 1
=3∗( )=6.461574− j 118.87998 [ pu]
tT1 tT1

Para las admitancias paralelas del transformador tenemos:

(a4
( t−1t )=3∗( Z1 ∗( t−1t ))=−0.323079+ j5.943999[ pu]
Y 0 T 1∗ '
T1
.6)

1−t 1 1−t (a4

Y 0 T 1∗
( )t 2
ZT1 t ( )
=3∗( ' ∗ 2 )=0.340083− j 6.256842[ pu] .7)

Impedancia de Aterrizaje:

14
 0.015+ j 0.075 (a4
Z ' AT 1= ( Z base )
=2.835539∗10−5+ j 1.417769∗10−4 [ pu ] .8)

Para los demás transformadores se procede exactamente de la misma manera por lo que se adjunta la
Tabla 4 y Tabla 5 para el transformador 2 y 3 respectivamente.
Tabla. 4 – Resultados Transformador 2 Secuencia (+), (-) y (0).2.
Parámetro Resultado T2
aT 2 0.06
ZT2[Ω] 0.978419 + j*29.051128 [Ω]
t2 1.03
Z’T2[pu] 0.001849 + j*0.054917 [pu]
Y’T2/t2 [pu] 0.594737 - j*17.658878 [pu]
0.017842 - j*0.529766 [pu]
Y T 2∗ ( t−1t )
1−t -0.017322 – j*0.514336 [pu]
Y T 2∗
( ) t2
Z’T2(0)[pu] 0.005547+j*0.16475[pu]
Y’T2/t2(0) [pu] 1.784211-j*52.976634[pu]
0.053526-j*1.589298[pu]
Y T 2 (0)∗( t−1t )[pu
]
1−t -0.051966-j*1.543008[pu]
Y T 2 (0)∗
( ) t2
[pu

]
Tabla. 5 – Resultados Transformador 3 Secuencia (+), (-) y (0).3.
Parámetro Resultado T3
aT 3 0.06
ZT3[Ω] 0.555608 + j*32.720796 [Ω]
T3 1.07
Z’T3[pu] 0.001050+ j*0.061854[pu]
Y’T3/t3 [pu] 0.256489 - j*15.105113 [pu]
0.017954 - j*1.057358 [pu]
Y T 3∗ ( t−1t )
1−t -0.016779 + j*0.988185 [pu]
Y T 3∗
( ) t2
Z’T3(0)[pu] 0.00315+j*0.185562[pu]
Y’T3/t3(0) [pu] 0.769467-j*45.315339[pu]
-0.050337+j*2.964555[pu]
Y T 3 (0) ( t−1t )
1−t -0.155898-j*4.629024[pu]
Y T 3 (0)
( ) t2
Z’T3AT[pu] 3.591682*10^-5+1.550095*10^-4[pu]
A continuación se procede a realizar el desarrollo del cálculo de parámetros de secuencia (+), (-) y

15
(0) de los Generadores del sistema:
Generador 1 (G1):
ZG 1¿ ¿ (a5.1)

Reemplazando valores obtenemos:

ZG 1¿ ¿
(a5.2
)
Para secuencia (-):
ZG 1¿ ¿ (a5.3)

Reemplazando valores obtenemos:

ZG 1¿ ¿
(a5.4
)
Para secuencia (0):
100 MVA (a5.5)

ZG 1(0)=
( 0.0058+ j0.1790 )∗ ( 247.5 MVA )
∗1

Z base

Reemplazando valores obtenemos:

ZG 1(0)=0.004848+ j 0.229494[ pu]
(a5.6
)
Luego impedancia de aterrizaje del Generador 1:
0.0103+ j 0.0971 (a5
Z ' ATG 1 = =1.947069∗10−5 + j1.835539∗10−4 [ pu] .7)
Z base

Para los demás Generadores se procede exactamente de la misma manera por lo que se adjunta la
Tabla 6 y Tabla7 para el generador 2 y 3 respectivamente.
Tabla.6 – Resultados Generador 2 Secuencia (+), (-) y (0).
Parámetro Resultado G2
Z(+)G2[pu] 0.0050 + j*0.097396 [pu]
Z(-)G2[pu] 0.0050+j*0.096354 [pu]
Z(0)G2[pu] 0.0050 + j*0.248438 [pu]
ZATG2[pu] 1.417769*10^-5 + j*9.924386*10^-4 [pu]
Tabla. 7 – Resultados Generador 3 Secuencia (+), (-) y (0).
Parámetro Resultado G3
Z(+)G3[pu] 0.003867+j*0.119333 [pu]
Z(-)G3[pu] 0.003867 + j*0.110000 [pu]
Z(0)G3[pu] 0.003867+j*0.291333 [pu]

16
ZATG3[pu] 7.939509*10^-6 + j*6.427221*10^-5 [pu]

17
Parte a:

Dibuje los diagramas de admitancias (o impedancias) de redes de secuencia y determine las matrices de
impedancia de secuencia positiva, negativa y cero. Para ello, especifique todas las consideraciones
pertinentes para efectos de cálculos de corrientes de cortocircuito simétrica inicial.
A continuación en las próximas figuras se pueden distinguir los diagramas de impedancia de secuencia
positiva, negativa y cero. Esta parte se divide en dos secciones, donde la primera sección corresponde a
diagramas con Tap nominales de los transformadores del sistema en estudio y la otra sección se
consideran las impedancias paralelas de modelo de transformador con cambiador de Tap.

Fig. 1 – Diagrama Admitancias Secuencia Positiva con Tap Nominales.

18
En la Figura 2 se observa en diagrama de admitancias de secuencia negativa con tap nominal con
valores en por unidad del sistema de estudio, se destaca que es exactamente idéntico al diagrama de la
Figura 1 debido a que los parámetros de secuencia positiva y negativa son idénticos.
Fig. 2 – Diagrama Admitancias Secuencia Negativa con Tap Nominales.

En la Figura 3 se realiza el diagrama de secuencia cero del sistema con tap nominales se destaca que
los valores de resistencia, admitancia y capacitancia de todos los elementos del sistema exceptuando
generadores han sido multiplicados por tres veces sus valores de secuencia positiva para obtener los
valores de secuencia cero.

19
 Fig. 3 – Diagrama Admitancias Secuencia Cero con Tap Nominales.

Mediante la construcción por bloques, se determina la matriz de admitancias de barra del sistema de
estudio considerando Tap de transformadores en su posición nominal.
Los diagramas de admitancias en por unidad realizados anteriormente son la base para poder
construir la matriz, los puntos de conexión de los elementos en las barras del sistema están
representados por la siguiente distribución por bloque:

20
Considerando los puntos de conexión en común en cada bloque, para el caso de la matriz Y bus considerando los taps nominales de los transformadores y
despreciando admitancias paralelas, en tendremos la siguiente matriz Y bus sec (+):
B1 B2 B3 B4 B5 B6 B7 B8 B9
B1 −Y T 1 0 0
Y T 1 +Y G 1 0 0 0 0 0

B2 0 0
0 Y T 2 +Y G 2 0 0 0 −Y T 2 0

B3 0 −Y T 3
0 0 Y T 3 +Y G 3 0 0 0 0

B4
−Y T 1 0 0 Y T 1 +Y L45 +Y L46 −Y L45 −Y L46 0 0 0

B5
Ybus 0 0 0 −Y L45 Y L45 +Y L57 0 −Y L57 0 0
sec(+)
= B6 −Y L69
0 0 0 −Y L46 0 Y L46 +Y L69 0 0

B7
0 −Y T 2 0 −Y L57 0 Y T 2 +Y L57 +Y L78 −Y L78 0
0
B8

0 0 0 0 0 −Y L78 Y L78+ Y L 89 −Y L89

0
B9

0 0 −Y T 3 0 −Y L69 0 −Y L89 Y T 3 +Y L69 +Y L89

0

21
 Reemplazando los valores en p.u. de cada elemento, la matriz Ybus será la siguiente:

2,582566 – -2,046165 +
0 0 0 0 0 0 0
j48,149851 j37,645327
1.138289 – -0,612579 +
0 0 0 0 0 0 0
j28,429004 j18,188644
0,545686 – -0,274444 +
0 0 0 0 0 0 0
j24.533570 j16,162471
– 2,046165 + 7,246306 – -3,712792 + -1,487349 +
0 0 0 0 0
j37,645327 j51,299317 j6,162981 j7,491008
Ybus
-3,712792 + 4,9807106 – -1,267919 +
sec(+) 0 0 0 0 0 0
= j6,162981 j9,667251 j3,504270
-1,487349 + 2,813668 – -1,326319 +
0 0 0 0 0
0 j7,491008 j10,678369 j3,187361
-0,612579 + -1,267919 + 2,206421 – -0,321985 +
0 0 0 0 0
j18,188644 j3,504270 j24,633769 j2,940855
-0,321985 + 0,625441 – -0,303455 +
0 0 0 0 0 0
j2,940855 j5,347451 j2,406596
-0,274444 + -1,326319 + -0,303455 + 1,904219 –
0 0 0 0 0
j16,162471 j3,187361 j2,406596 j21,756428

Destacar que la matriz obtenida es idéntica para la matriz de secuencia negativa del sistema eléctrico en estudio por lo que su construcción se omitirá y se
procederá a mostrar la matriz de secuencia cero que es distinta a las mencionadas anteriormente debido a la topología y tipo de conexionado de algunos
elementos del sistema.

22
 B1 B2 B3 B4 B5 B6 B7 B8 B9
B1 0
Y G 1 +3 Y ATG 1 0 0 0 0 0 0 0

B2
0 Y G 2 +3 Y ATG 1 0 0 0 0 0 0 0
B3 −Y T 3 −3 Y ATT 3
0 0 Y T 3 +3 Y ATT 3 +Y G 3+3 Y 0ATG 3 0 0 0 0

B4
0 0 0 Y T 1 +3 Y ATT 1 +Y L45 −Y
+Y L46
L45 −Y L46 0 0 0

B5
0 0 0 −Y L45 Y L45 +Y L57 0 −Y L57 0 0
Ybus sec(0)
=
B6
0 0 0 −Y L46 0 Y L46 +Y L69 0 0 −Y L69

B7
0 0 0 0 −Y L57 0 Y L57 +Y L 78 −Y L78 0

B8
0 0 0 0 0 0 −Y L78 Y L78+ Y L 89 −Y L89

B9

0 0 −Y T 3 0 −Y L69 0 −Y L89 Y T 3 +Y ATT 3+Y L69 +Y L

0

A continuación se muestran los resultados obtenidos de la matriz de secuencia cero en p.u.

23
 0,092679 –
0 0 0 0 0 0 0 0
j4,344987
0,081469 –
0 0 0 0 0 0 0 0
j4,018687
0,252749 – -0,206980 +
0 0 0 0 0 0 0
j9,475201 j6,045587
2,421425 – -1,237597 + -0,495783 +
0 0 0 0 0 0
j17,031901 j2,054327 j2,497003
Ybus
sec(0)
-1,237597 + 1,660237 – -0,422639 +
0 0 0 0 0 0
= j2,054327 j3,222417 j1,168089
-1,487349 + 0,937889 – -0,442106 +
0 0 0 0 0
0 j7,491008 j3,559456 j1,062454
-0,422639 + 0,529968 – -0,107328 +
0 0 0 0 0 0
j1,168089 j2,148374 j0,980285
-0,107328 + 0,208480 – -0,101152 +
0 0 0 0 0 0
j0,980285 j1,782484 j0,802199
-0,206980 + -0,442106 + -0,101152 + 0,750238 –
0 0 0 0 0
j6,055587 j1,062454 j0,802199 j7,910240

24
 Departamento de Ingeniería Eléctrica y Electrónica
Universidad del Biobío, Chile

En la Figura 4 podemos observar la red de secuencia positiva correspondiente al sistema eléctrico en

estudio, considerando las admitancias paralelas correspondientes a los Tap de los transformadores y
admitancias paralelas de las líneas que conectan el sistema, esto se realiza para poder demostrar cómo se
comportan las niveles de cortocircuito del sistema para así ver que tan influyentes son estas admitancias
en la contribución del cortocircuito.
Fig. 4 – Diagrama Admitancias Secuencia Positiva con Admitancias Paralelas.

25
Departamento de Ingeniería Eléctrica y Electrónica
Universidad del Biobío, Chile

Fig. 5 – Diagrama Admitancias Secuencia Negativa con Admitancias Paralelas.

26
 Departamento de Ingeniería Eléctrica y Electrónica
Universidad del Biobío, Chile

Fig. 6 – Diagrama Admitancias Secuencia Cero con Admitancias Paralelas.

En consecuencia, considerando lo propuesto por la Figura 4, se procede a construir la matriz Y BUS de

secuencia positiva:

27
 B1 B2 B3 B4 B5 B6 B7 B8 B9
B1 −Y T 1 0 0
Y T 1 +Y G 1 +Y TAPBT 1 0 0 0 0 0

B2 0 0
0 Y T 2 +Y G 2+Y TAPBT 2 0 0 0 −Y T 2 0

B3 0 −Y T 3
0 0 Y T 3 +Y G 3+Y TAPBT 3 0 0 0 0

B4
−Y T 1 0 0 Y T 1 +Y TAPAT 1 +Y L45 +Y−Y +Y PL 45+Y−Y
L46L45 PL46
L46 0 0 0

B5
Ybus 0 0 0 −Y L45 Y L45 +Y L57+ Y PL 45+Y
0 PL57 −Y L57 0 0
sec(+)
= B6 −Y L69
0 0 0 −Y L46 0 Y L46 +Y L69 +Y PL 46+Y PL69
0 0

B7
0 −Y T 2 0 −Y L57 0 Y T 2 +Y L57 +Y L78+ +Y PL57−Y PL78 + Y TAPAT 2 0
+Y L78
0
B8

0 0 0 0 0 −Y L78 Y L78+ Y L 89+Y PL78 +Y−Y

PL 89
L89

0
B9

0 0 −Y T 3 0 −Y L69 0 −Y L89 Y T 3 +Y TAPAT 3 +Y L69 +Y

0

A continuación se muestran los resultados obtenidos de la matriz de secuencia positiva en p.u.

28
 2,576898 – -2,153858+
0 0 0 0 0 0 0
j48,045571 j39,626660
1,137769 – -0,594737+
0 0 0 0 0 0 0
j28,413606 j17,658878
0,544511 – -0,256489 +
0 0 0 0 0 0 0
j24,464397 j15,105113
-2,153857+ 7,473326 – -3,712792 + -1,487349 +
0 0 0 0 0
j39,626660 j55,175052 j6,162981 j7,490081
Ybus
sec(+)
-3,712792 + 4,982078 – -1,267919 +
0 0 0 0 0 0
= j6,162981 j9,219230 j3,504269
-1,487349 + 2,815148 – -1,326319 +
0 0 0 0 0
0 j7,491008 j10,259080 j3,187361
-0,594737+ -1,267919 + 2,171762 – -0,321985 +
0 0 0 0 0
j17,658878 j3,504269 j22,714784 j2,940855
-0,321985 + 0,630399 – -0,303455 +
0 0 0 0 0 0
j2,940855 j4,261868 j2,406596
-0,256489 + -1,326319 + -0,303455 + 1,874449 –
0 0 0 0 0
j15,105113 j3,187361 j2,406596 j19,013483

29
 Departamento de Ingeniería Eléctrica y Electrónica
Universidad del Biobío, Chile

B1 B2 B3 B4 B5 B6 B7 B8 B9
B1 Y G 1 +3 Y ATG 1 +Y TAPBT01 0 0 0 0 0 0 0
B2 0 Y G 3 +3 Y ATG 3 +Y TAPBT02 0 0 0 0 0 0
B3 0 0 Y G 3 +3 Y ATG 3 +Y TAPBT 3 +3
0 Y ATT 3 +Y T 3 0 0 0 0 3 Y ATT 3 +Y T 3
B4 0 0 0 Y L45 +Y L46 +Y PL 45+ Y −Y
PL 46
L45+Y TAPAT−Y
1 +3 Y ATT 1 +Y T 1 0
L46 0 0
B5 0 0 0 −Y L45 Y L45 +Y L57+ Y PL 450 +Y PL57 −Y L57 0 0
Ybus
sec(0)=
B6
0 0 0 −Y L46 0 Y L46 +Y L69 +Y PL 46+Y
0 PL69 0
−Y L69
B7 0 0 0 0 −Y L57 0 Y L57 +Y L 78+Y PL57−Y
+Y PL
L7878 +Y TAPAT 2 0
B8 0 0 0 0 0 0 −Y L78 Y L78+ Y L 89+Y PL78 +Y−YPLL89
89
B9 Y L69+ Y L 89+Y PL69 +Y PL 89+Y TAPAT 3
0 0 3 Y ATT 3 +Y T 3 0 0 −Y L69 0 −Y L89

30
Reemplazando los valores en p.u. de cada elemento, la matriz Ybus para este caso será la siguiente:

-0,247414+
0 0 0 0 0 0 0 0
j1,911854
0,133437 –
0 0 0 0 0 0 0 0
j5,561696
-0,200873
0,296982 –
0 0 0 0 0 0 0 +
j12,264459
j5,869991
– 2,153858 -1,237597 -0,495783
2,815947 –
+ 0 0 + + 0 0 0
j23,368242
j39,626660 j2,054327 j2,497003
-0,422639
-1,237697 + 1,664338-
0 0 0 0 + 0 0
Ybus j2,054327 j1,878355
j1,168089
=
-0,442106
-1,487349 + 0,942328 –
0 0 0 0 0 +
0 j7,491008 j2,301589
j1,062454
-0,107328
-0,422639+ 0,491329–
0 0 0 0 0 + 0
j1.168089 j2,019279
j0,980285
-0,107328 -0,101152
0,223357 –
0 0 0 0 0 0 + +
j1,474264
j0,980285 j0,802199
-0,442106
-0,200873 + -0,101152+ 0,708687 –
0 0 0 0 + 0
j5,869991 j0,802199 j2,677881
j1,062454

31
Departamento de Ingeniería Eléctrica y Electrónica
Universidad del Biobío, Chile

V B 1−V B 4 (g1.1)
I G 1= =1.683101+ j 0.435949[ p .u ]
ZT 1

SG 1=V B 1∗I ¿B 1 =168.3101 MW − j 43.5949 Mvar (g1.2)

32
Departamento de Ingeniería Eléctrica y Electrónica
Universidad del Biobío, Chile

Parte h:
Efectúe un mapeo de los flujos de potencia activa y reactiva a través de todos los elementos del
sistema, mostrando la potencia activa y reactiva inyectada por cada generador y consumida
por cada carga, así también la cargabilidad de cada elemento del SEP.
Fig. 6 – Flujo de Potencias con tap nominales.

La cargabilidad de un elemento en este caso un dispositivo eléctrico como un transformador

corresponde a relación entre la potencia eléctrica con la que está trabajando este elemento versus la
potencia eléctrica de diseño del transformador. En el caso de los generadores corresponde algo
similar pero estos equipos cuentan con una curva de operación donde el fabricante establece ciertos
límites de trabajo de operación que no pueden ser sobrepasados, debido a lo delicado y valioso que
es este dispositivo. Para las líneas de transmisión el fabricante del cable establece un límite de
corriente de trabajo del conductor por lo tanto para determinar la cargabilidad de la línea se utiliza la
corriente de operación versus la corriente máxima de operación así podemos obtener en que
porcentaje de carga está operando la línea. A continuación procedemos a desarrollar el cálculo de
cargabilidad de cada elemento del sistema de estudio.
Cargabilidad Generadores:
Generador 1.-

S¿G 1 173.8643 MVA (h1.1)

ηG 1 = ∗100= ∗100=70.25 %
SG 1 247.4969 MVA

S¿Gi: Potencia aparente con la que está operando el generador.

SGi: Potencia aparente máxima de operación del generador.
Generador 2.-

S¿G 2 142.9232 MVA (h1.2)

ηG 2 = ∗100= ∗100=74.44 %
SG 2 192 MVA

Generador 3.-

S ¿G 3 130.11 MVA (h1.3)

ηG 3 = ∗100= ∗100=86.74 %
SG3 150 MVA

Cargabilidad Transformadores:
Transformador 1.-

33
Departamento de Ingeniería Eléctrica y Electrónica
Universidad del Biobío, Chile

S ¿T 1 175.63 MVA (h1.4)

ηT 1 = ∗100= ∗100=70.25 %
ST 1 250 MVA

S¿Ti: Potencia aparente con la que está operando el Transformador.

STi: Potencia aparente máxima de operación del Transformador.
Transformador 2.-

S ¿T 2 148.76 MVA (h1.5)

ηT 2 = ∗100= ∗100=74.38 %
ST 2 200 MVA

Transformador 3.-

S ¿T 3 136.77 MVA (h1.6)

ηT 3 = ∗100= ∗100=88.24 %
ST 3 155 MVA

Cargabilidad Líneas:
Línea 45.-

I ¿L 45 250.56 A (h1.7)
η L45= ∗100= ∗100=83.52%
I L 45 300 A

I ¿Li: Corriente con la que está operando Línea.

I Li: Corriente máxima de operación de la Línea.
Línea 46.-

I ¿L 46 187.34 A (h1.8)
η L46 = ∗100= ∗100=62.45 %
I L 46 300 A

Línea 57.-

I ¿L57 166.48 A (h1.9)

ηT 57= ∗100= ∗100=41.62 %
I L57 400 A

Línea 69.-

34
Departamento de Ingeniería Eléctrica y Electrónica
Universidad del Biobío, Chile

I ¿L69 105.69 A (h1.10)

η L69= ∗100= ∗100=26.42 %
I L69 400 A

Línea 78.-

I ¿L78 287.24 A (h1.11)

η L78= ∗100= ∗100=57.45 %
I L78 500 A

Línea 89.-

I ¿L89 247.07 A (h1.12)

ηT 89= ∗100= ∗100=49.41 %
I L89 500 A

35