B102 - Clase 6 - Fallas Balanceadas

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 51

Universidad Nacional del Callao

Facultad de Ingeniería Eléctrica y Electrónica


Sección de Postgrado

Maestría en Ingeniería Eléctrica


Mención en Gestión de Sistemas de Energía Eléctrica

Curso B102 – Modelación de Sistemas de


Energía Eléctrica

Expositor : Manfred F. Bedriñana Aronés


Doctor en Ingeniería Eléctrica
Ingeniero Electricista C.I.P. Nº 95644

2013
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

Capítulo 8

Análisis de Fallas Balanceadas

27/10/2013 [email protected] 2
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

A. Análisis de Fallas
„ Tipos de Análisis
„ Flujo de Potencia – evaluar las condiciones de operación normal.
„ Análisis de Fallas – evaluar las condiciones de operación anormal.
anormal
„ En que consiste el análisis de fallas?
„ Determinar las tensiones en barra y corrientes en las líneas durante varios tipos
de fallas.
„ Uso de los resultados del análisis de fallas
„ Especificar las características técnicas de interruptores y fusibles.
„ Ajuste de relés de protección.
„ Especificar
ifi la l impedancia
i d i de d transformadores
f d y generadores.
d

27/10/2013 [email protected] 3
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

A. Análisis de Fallas
„ Ti
Tipos de
d Fallas:
F ll
„ Fallas desbalanceadas
ƒ Falla fase a tierra (A)
( )
ƒ Falla fase a fase (B)
ƒ Falla dos fases a tierra (C)
„ Fallas balanceadas
ƒ Falla trifásica (D)
ƒ Falla trifásica a tierra (E)
„ Otras en instalaciones
subterráneas
ƒ Falla piloto a fase (F)
ƒ Falla piloto a tierra (G)

27/10/2013 [email protected] 4
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

A. Análisis de Fallas
„ Ej
Ejemplos:
l
„ Falla línea a tierra „ Falla línea a línea

27/10/2013 [email protected] 5
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

A. Análisis de Fallas
„ La magnitud de las corrientes de falla dependerá de:
„ La impedancia de la red.
„ La impedancia interna de los generadores
generadores.
„ La resistencia de la falla (resistencia del arco).
„ Las impedancias de la red están gobernadas por:
„ Impedancia de la línea de transmisión.
„ Conexiones e impedancia
p del transformador.
„ Conexiones y resistencias de tierra.
„ El comportamiento del generador está dividido en tres periodos:
„ Periodo sub transitorio: abarca los primeros ciclos.
„ Periodo transitorio: abarca un periodo de tiempo relativamente largo (30 ciclos).
„ Periodo de estado estacionario

27/10/2013 [email protected] 6
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

B. Estudios de Fallas
„ P i d sub
Periodo b transitorio,
i i xg = xd”:
„ Determinar la capacidad de interrupción de los interruptores de muy alta tensión.
„ Determinar el tiempo de operación del sistema de relés de protección para las
redes de muy alta tensión

„ Periodo it i xg = xd’:
P i d ttransitorio,
„ Determinar la capacidad de interrupción de los interruptores de alta tensión.
„ Determinar el tiempop de operación
p del sistema de relés de protección
p ppara las
redes de alta tensión.
„ Estudios de estabilidad transitoria.

27/10/2013 [email protected] 7
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

C. Representación de Fallas Balanceadas


„ Una falla representa un cambio en la estructura de la red
„ Equivale a la adición de una impedancia en el lugar de la falla.
„ Si la impedancia de falla es cero,
cero la falla se denomina falla sólida o franca
franca.

„ Fallas trifásicas o balanceadas


„ Definido como un cortocircuito simultáneo que cruza las tres fases.
„ Considerado como el más severo tipo
p de falla encontrada.
„ Debido que la red es balanceada, los cortocircuitos trifásico pueden ser resueltos
usando modelos equivalentes monofásicos.
„ P d usarse ell método
Puede ét d de
d Thevenin
Th i para ell análisis
áli i de
d falla.
f ll

27/10/2013 [email protected] 8
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

D. Método de Thevenin
„ Premisas iniciales
„ Las resistencias de la red se desprecian.
„ Los generadores se modelan como una fem detrás de la reactancia sub transitoria
o transitoria.
„ Las capacitancias shunt se desprecian.
„ Se considera al sistema sin cargas.

27/10/2013 [email protected] 9
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

D. Método de Thevenin
Ej
Ejemplo
l aplicativo
li i
Se tiene un diagrama unifilar de impedancias como el mostrado. Los
ggeneradores están en fase y mantienen una tensión constante igual
g a 1 pp.u.
Determinar las corrientes en líneas, tensiones en barras durante una falla
trifásica en la barra 3 con valor de Zf = 0.16 p.u.
„ La falla se simula por la adición (switcheo) de una impedancia de falla en la
barra de aplicación de la falla.

Diagrama de impedancias
27/10/2013 [email protected] 10
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

D. Método de Thevenin
„ Las variaciones
L i i de
d las
l tensiones
i de
d lal redd es equivalente
i l a adicionar
di i una
tensión de barra pre falla en la barra con todas las fuentes cortocircuitadas.

Red en el estado prefalla.

Red en estado de falla.


Red considerando las variaciones
de tensiones por corriente de falla.
27/10/2013 [email protected] 11
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

D. Método de Thevenin
„ Para ell cálculo
P ál l ded la
l impedancia
i d i Thevenin
Th i en la
l barra
b 3 (Z33) convertimos
i ell
sistema de delta a estrella y procedemos a reducir la red.

27/10/2013 [email protected] 12
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

D. Método de Thevenin
„ Para ell cálculo
P ál l ded la
l corriente
i de
d falla
f ll usamos lal tensión
ió Thevenin
Th i Vth,
V h que
es igual en este caso a 1 p.u. (las fuentes de tensión de los generadores están
en fase y con tensión 1 p.u.).

27/10/2013 [email protected] 13
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

D. Método de Thevenin
„ L corrientes
Las i de
d falla
f ll en los
l generadores
d serán:
á

„ Las variaciones de tensión en la falla serán:

„ Las tensiones en la falla serán:

27/10/2013 [email protected] 14
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

D. Método de Thevenin
„ P soluciones
Para l i más
á exactas
„ Uso de las tensiones pre falla en barra las cuales se pueden obtener a partir de los
resultados de una solución de flujo de potencia.
„ Incluir las cargas – para preservar la linealidad, convertir la carga usando el
modelo de impedancia constante.
„ El Teorema de Thevenin p permite obtener las variaciones de tensiones en barra.
„ Las tensiones de barra se obtienen por la superposición de las tensiones de pre
falla y las variaciones de tensiones en la barra.
„ La corriente en cada rama o línea puede calcularse.
calcularse

27/10/2013 [email protected] 15
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

E. Capacidad de Cortocircuito (SCC)


„ Mide lla robustez
Mid b eléctrica
lé i de d la
l barra.
b
„ Se determina en MVA.
„ Determina la dimensión de barras y la capacidad de interrupción de los
interruptores.
„ Definición (MVA):

„ En por unidad:
id d

„ SCC (MVA), SB (MVA) y Xkk (p.u.).


Asume la tensión en barra igual a 1 p.u.

27/10/2013 [email protected] 16
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

E. Capacidad de Cortocircuito (SCC)


„ E
Encontrar
t ell SCC para la
l barra
b 3.
3

27/10/2013 [email protected] 17
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

F. Análisis de Fallas usando la Matriz de Impedancia


„ R d ió de
Reducción d la
l redd por ell método
ét d dde Th
Thevenin
i no es eficiente
fi i t
„ Difícil de aplicar a sistemas de grande porte (miles de barras).
„ Formación mediante álgebra de matrices
„ Buscar una matriz donde los elementos diagonales representan la impedancia del
generador para las barras.
„ Considere el siguiente sistema:

+ Opera bajo condiciones balanceadas.


+ Cada generador es representado por
una fem constante detrás de una
reactancia adecuada (xd, xd’ o xd”).
+ Líneas representadas por su
equivalente
i l t en modelo
d l π.

27/10/2013 [email protected] 18
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

F. Análisis de Fallas usando la Matriz de Impedancia

„ Coloque las tensiones pre falla en un vector.

„ Reemplace
p la carga
g ppor un
modelo de impedancia
constante con las tensiones
pre-falla de barra.

„ Las variaciones de las tensiones


causados por la falla son equivalentes
a colocar una tensión de falla en la
barra de falla con todas las fuentes
de tensión cortocircuitadas.

27/10/2013 [email protected] 19
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

F. Análisis de Fallas usando la Matriz de Impedancia


„ Usando superposición, las tensiones de falla
se calculan como la suma de tensiones pre falla
con las variaciones de tensión en falla.

„ Las variaciones de tensiones en la barra


se pueden calcular a partir de la
matriz de la red.

27/10/2013 [email protected] 20
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

F. Análisis de Fallas usando la Matriz de Impedancia

Será necesario
invertir la matriz
admitancia de barra

Variaciones de tensiones en barra


(estado en falla)

27/10/2013 [email protected] 21
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

F. Análisis de Fallas usando la Matriz de Impedancia

Es mejor usar la
matriz impedancia
de barra

Corriente de falla

27/10/2013 [email protected] 22
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

G. Ejemplo aplicativo
„ Diagrama de admitancias

Será necesario
invertir la matriz
admitancia de barra

27/10/2013 [email protected] 23
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

G. Ejemplo aplicativo
„ Calculo de tensiones en barra en estado de falla

Todos los cálculos usan elementos de la matriz


impedancia de barra

27/10/2013 [email protected] 24
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

H. La Matriz Impedancia de Barra


„ Definición:

„ Formación directa de la matriz


„ Inversión de la matriz de admitancias de barra presenta un esfuerzo de n3.
„ Para redes de pequeño y mediano tamaño, el desarrollo directo de la matriz es de
menor esfuerzo.
„ Para redes de gran tamaño,
tamaño es preferible la programación dispersa de la matriz
con la técnica de eliminación gaussiana.

27/10/2013 [email protected] 25
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

I. Formación de la Matriz de Impedancia de Barra


„ Técnicas de teoría de grafos ayudan a explicar el proceso de construcción.
„ Grafos de la red antes de la aplicación de la falla:

Nodo de barra = Vértice de grafo Rama o línea = Borde

Grafo conectado

27/10/2013 [email protected] 26
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

I. Formación de la Matriz de Impedancia de Barra


„ Árbol de un grafo conectado es un subgrafo conectando todos los nodos sin
Á
formar un lazo.
„ Número de ramas de un árbol: b = n – 1 (n número de nodos incluyendo el
nodo de referencia 0).
„ En cada árbol existen elementos remanentes (restantes) que son llamados de
enlaces El número de enlaces es igual a: l = e – b (e número de elementos
enlaces.
del grafo)

Rama de Enlace o rama


árbol asociado a un
Árbol y co-árbol
co árbol lazo básico
27/10/2013 [email protected] 27
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

I. Formación de la Matriz de Impedancia de Barra


„ Construcción básica de la red y la matriz. Se asume inicialmente que existe
una matriz Zbus.

Red parcial con m barras y


un nodo de referencia 0

27/10/2013 [email protected] 28
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

J. Adición de una Línea


„ Al adicionar una línea desde una barra existente (barra p o referencia 0) a
una nueva barra q, una nueva barra es adicionada a la red parcial creando
una nueva fila y columna

27/10/2013 [email protected] 29
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

K. Adición de una Línea a una Línea Existente

27/10/2013 [email protected] 30
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

L. Adición de una Línea desde la Referencia

27/10/2013 [email protected] 31
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

M. Cerrar un Lazo
„ Cuando el elemento adicionado es un enlace de un coárbol entre p y q
ninguna nueva barra es creada.

27/10/2013 [email protected] 32
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

N. Reducción de Kron
„ Primera referencia:

„ Eliminación de un nodo del sistema

27/10/2013 [email protected] 33
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

O. Adición de una Línea entre dos Líneas

„ A continuación, ejecute la reducción de Kron en Zll

27/10/2013 [email protected] 34
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

P. Adición de una Línea desde una Línea y Referencia

„ A continuación, ejecute la reducción de Kron en Z11

27/10/2013 [email protected] 35
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

Q. Reducción de Kron
„ La reducción de Kron elimina un eje (fija y columna) de una matriz sin
perder la influencia numérica del eje.

27/10/2013 [email protected] 36
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

R. Reglas de Construcción de la Matriz Zbus


„ Regla 1: Adición de una rama a la referencia
„ Comenzar con la red existente de la matriz [m x m]
„ Crear una nueva matriz de red [(m + 1) x (m + 1)] con
„ Una fila y columna con elementos fuera de la diagonal iguales a 0.
„ El elementos de la diagonal (m + 1),(m + 1) debe ser llenado con el elemento de
impedancia zq00

„ Regla
g 2: Adición de una rama a una barra existente
„ Conexión de la barra existente p
„ Comenzar con la matriz de la red existente [m x m]
„ C
Crear una nueva matriz
i de
d redd [(m
[( + 1) x ((m + 1)] con
„ La fila y columna nueva con elementos fuera de la diagonal copiado de la fila p y
columna p.
„ El elemento
l t de
d la
l diagonal
di l (m
( + 1)
1),(m
( + 1) ddebe
b ser llenado
ll d con ell elemento
l t de
d
impedancia zpq más la diagonal de impedancia Zpp

27/10/2013 [email protected] 37
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

R. Reglas de Construcción de la Matriz Zbus


„ Regla 3: Adición de una rama de enlace entre barras que no son referencia
„ Conexión entre barras existentes p y q
„ Comenzar con la matriz de la red existente [m x m]
„ Crear una nueva matriz de red [(m + 1) x (m + 1)] con
„ Una fila y columna nueva con elementos fuera de la diagonal llenados con una copia
de la fila q menos la fila p y columna q menos la columna p.
p
„ El elemento de la diagonal (m + 1),(m + 1) llenado con zpq + Zpp + Zqq – 2 Zpq
„ Realizar la reducción de Kron en la fila y columna m + 1.

„ Regla 4: Adición de una rama de enlace entre barra existente y la referencia


„ Conexión de la barra existente p y la referencia
„ Comenzar con la red existente de la matriz [m x m]
„ Crear una nueva matriz de red [(m + 1) x (m + 1)] con
„ Una fila y columna nueva con elementos fuera de la diagonal llenados con una copia
de la fila q menos la fila p y columna q menos la columna p.
„ El elemento de la diagonal (m+1),(m+1) llenado con zp0 + Zpp
„ Realizar la reducción de Kron en la fila y columna m + 1.
1

27/10/2013 [email protected] 38
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

S. Ejemplo aplicativo
„ Construir la matriz impedancia del ejemplo aplicativo.

„ La adición de las líneas en el siguiente orden: 1-0, 2-0, 1-3, 1-2 luego 2-3.

27/10/2013 [email protected] 39
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

S. Ejemplo aplicativo

Adición de elemento 1-0

Adición de elemento 2-0 (existente 0, nuevo 2)

Adición de elemento 11-3


3 (existente 1, nuevo 3)

27/10/2013 [email protected] 40
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

S. Ejemplo aplicativo

Adición de elemento 1-2


(existente 1 y 2)

Reducción de Kron de barra 4

27/10/2013 [email protected] 41
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

S. Ejemplo aplicativo

Adición de elemento 2-3


(existente 2 y 3)

Reducción de Kron de barra 4

27/10/2013 [email protected] 42
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

T. Ejemplo MATLAB
„ Se tiene
S i un diagrama
di unifilar
ifil de
d impedancias
i d i como ell mostrado. d Los
L
generadores están en fase y mantienen una tensión constante igual a 1 p.u.
Determinar la matriz impedancia de barra.
„ Usar el programa zbuild.m. Ingresar los datos de la siguiente forma:

Barra Barra R X
Ini Fin (p.u.) (p.u.)

27/10/2013 [email protected] 43
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

T. Ejemplo MATLAB
„ I
Ingresar ddatos ddell sistema:
i
>> zdata = [ 0 1 0 0.2
0 2 0 0.4
1 2 0 0.8
1 3 0 0
0.4
4
2 3 0 0.4];

„ Obtener matriz Z:
>> Zbus = zbuild(zdata)
This program forms the complex bus impedance matrix by the method
of building algorithm. Bus zero is taken as reference.

++ The program requires:


a. Impedance data: zdata

--> Adding a branch from a new bus to reference bus 0


--> Adding a branch from new bus to an existing bus
--> Adding a link between two old buses
Zbus =
0 + 0.1600i 0 + 0.0800i 0 + 0.1200i
0 + 0.0800i 0 + 0.2400i 0 + 0.1600i
0 + 0
0.1200i
1200i 0 + 0.1600i
0 1600i 0 + 0.3400i
0 3400i

27/10/2013 [email protected] 44
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

U. Ejemplo MATLAB
„ Se tiene
S i un diagrama
di unifilar
ifil ded impedancias
i d i como ell mostrado.d Los
L
generadores están en fase y mantienen una tensión constante igual a 1 p.u.
Determinar las corrientes en líneas, tensiones en barras durante una falla
trifásica en la barra 3 con valor de Zf = j*0.16
j*0 16 p.u.
pu
„ Usar los programas zbuild.m y symfault.m. Ingresar los datos de la
siguiente forma:

Nota: No se consideran los modelos de capacitores shunt de línea y cargas.

27/10/2013 [email protected] 45
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

U. Ejemplo MATLAB
„ I
Ingresar ddatos ddell falla
f ll y calcular
l l cortocircuito
i i trifásico:
ifá i
>> symfault(zdata, Zbus)
The program symfault is designed for the balanced three-phase
fault analysis of a power system network

--> Three-phase balanced fault analysis

Enter Faulted Bus No. -> 3

Enter Fault Impedance Zf = R + j*X in complex form (for bolted fault enter 0). Zf = j*0.16

Balanced three-phase fault at bus No. 3


Total fault current = 2.0000 per unit

Bus Voltages during fault in per unit

Bus Voltage Angle


No. Magnitude degrees
1 0.7600 0.0000
2 0.6800 0.0000
3 0.3200 0.0000

Line currents for fault at bus No. 3 Contribuciones de la corriente


From To Current Angle
de cortocircuito.
Bus Bus Magnitude degrees
G 1 1.2000 -90.0000 Suma iiguall a
S
1 2 0.1000 -90.0000
1 3 1.1000 -90.0000 If bus 3 = 2.0 p.u.
G 2 0.8000 -90.0000
2 3 0.9000 -90.0000
3 F 2
2.0000
0000 -90.0000
90 0000

27/10/2013 [email protected] 46
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

V. Ejemplo MATLAB
„ Se tiene
S i un diagrama
di unifilar
ifil
mostrado. Desprecie las
capacitancias shunt de las líneas y
las cargas
cargas. Las tensiones pre falla en
las barras son mantenidas en 1∠0°.
Determinar las corrientes en líneas,
tensiones en barras durante una falla
trifásica sólida en la barra 8.
„ Usar los programas zbuild.m y
symfault m Ingresar los datos
symfault.m.
como:

27/10/2013 [email protected] 47
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

V. Ejemplo MATLAB
„ I
Ingresar ddatos
t del
d l sistema
it y obtener
bt matriz
t i Z:
Z
>> zdata = [ 1 0 0.00 0.20
10 0 0.00 0.15
11 0 0.00 0.25
1 2 0
0.00
00 0
0.06
06
2 3 0.08 0.30
2 5 0.04 0.15
2 6 0.12 0.45
3 4 0.
0.10
0 0.
0.400
3 6 0.04 0.40
4 6 0.15 0.60
4 9 0.18 0.70
4 10 0.00 0.08
5 7 0.05 0.43
6 8 0.06 0.48
7 8 0.06 0.35
7 11 0.00 0.10
8 9 0
0.052
052 0
0.48];
48]
Zbus=zbuild(zdata)
This program forms the complex bus impedance matrix by the method
of building algorithm. Bus zero is taken as reference.

++ The program requires:


a. Impedance data: zdata

--> Adding a branch from a new bus to reference bus 0


--> Adding a branch from new bus to an existing bus
--> Adding a link between two old buses
27/10/2013 [email protected] 48
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

V. Ejemplo MATLAB
„ I
Ingresar ddatos ddell falla
f ll y calcular
l l cortocircuito
i i trifásico:
ifá i
>> symfault(zdata, Zbus)
The program symfault is designed for the balanced three-phase
fault analysis of a power system network

--> Three-phase balanced fault analysis

Enter Faulted Bus No. -> 8

E t
Enter F
Fault
lt I
Impedance
d Zf = R + j*X i
in complex
l f
form (f
(for b
bolted
lt d f
fault
lt enter
t 0)
0). Zf = 0

Balanced three-phase fault at bus No. 8


Total fault current = 3.3319 per unit

Bus Voltages during fault in per unit

Bus Voltage Angle


No. Magnitude degrees
1 0.8082 -1.8180
2 0.7508 -2.5443
3 0.6882 -1.5987
4 0.7491 -2.4902
5 0.7007 -2.3762
6 0
0.5454
5454 -1
1.0194
0194
7 0.5618 -3.8128
8 0.0000 0.0000
9 0.3008 2.4499
10 0.8362 -1.4547
11 0
0.6866
6866 -2.2272
2 2272

27/10/2013 [email protected] 49
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

V. Ejemplo MATLAB
„ I
Ingresar ddatos ddell falla
f ll y calcular
l l cortocircuito
i i trifásico:
ifá i
>> symfault(zdata, Zbus)

Line currents for fault at bus No. 8

From To Current Angle


Bus Bus Magnitude degrees
G 1 0.9697 -82.4034
1 2 0.9697 -82.4034
2 3 0
0.2053
2053 -87.8751
87 8751
2 5 0.3230 -79.9626
2 6 0.4427 -81.6497
3 6 0.3556 -88.0987 Contribuciones de la
4 3 0.1503 -88.4042
4 6 0.3305 -82.3804 corriente de
4 9 0.6229 -81.3672 cortocircuito.
5 7 0.3230 -79.9626
6 8 1.1274 -83.8944
7 8 1.5820 -84.0852 Suma igual a
8 F 3.3319 -83.5126 If bus 8 = 3.3319 p.u.
9 8 0.6229 -81.3672
G 10 1.1029 -82.6275
10 4 1.1029 -82.6275
G 11 1
1.2601
2601 -85
85.1410
1410
11 7 1.2601 -85.1410

27/10/2013 [email protected] 50
Maestría en Ingeniería Eléctrica
Curso B102 – Modelación de Sistemas de Energía Eléctrica

W. Ejemplo MATLAB
„ Se tiene
S i un diagrama
di unifilar
ifil
mostrado. Considerar las
capacitancias shunt de las líneas, pero
no las cargas
cargas. Las tensiones pre falla
en las barras son mantenidas en
1∠0°. Determinar las corrientes en
líneas, tensiones en barras durante
una falla trifásica sólida en la barra 8.
„ Usar los programas zbuildpi.m y
symfault m Ingresar los datos como:
symfault.m.

27/10/2013 [email protected] 51

También podría gustarte