Ingeniera Electromecnica.
Reporte:
Anlisis de fallas utilizando matriz de impedancias.
Materia:
Sistemas Elctricos de potencia de potencia II.
Alumno:
Delgado Martnez Daniel
Gonzlez Ponce Jess
Ortiz Martnez Julio Adrin
Numero de control:
IS12110165
Profesor:
Rafael Barajas Vidal
ndice
1.
Introduccin. ................................................................................................................... 3
2.
Objetivo. .......................................................................................................................... 4
3.
Marco terico. ................................................................................................................. 5
4.
Desarrollo. ....................................................................................................................... 8
5.
Caso de estudio. ............................................................................................................ 11
6.
Resultados. .................................................................................................................... 12
7.
Conclusin. .................................................................................................................... 14
8.
Bibliografa .................................................................................................................... 14
1. Introduccin.
Los estudios de fallas en redes elctricas permiten deducir la mejor forma para
enfrentar los problemas de ser que estos ocurriesen en la realidad, as mismo para
recrear situaciones ocurridas en la vida real. Cabe mencionar que a partir de los
anlisis de fallas se pueden escoger las protecciones adecuadas para las redes
elctricas ante las distintas contingencias. Existen dos mtodos para realizar el
anlisis de fallas, las cuales son:
Mtodo circuital: sistema por unidad, con la aplicacin de los, mtodos de la
F.E.M o Thevenin.
Mtodo de Zbus o impedancia de barra.
En el presente trabajo se mostrara un algoritmo en el software Matlab que realice
un anlisis de falla utilizando el mtodo de Zbus.
2. Objetivo.
Realizar un anlisis de falla mediante el mtodo de zbus por medio de un
algoritmo desarrollado en Matlab
3. Marco terico.
Anlisis sistemtico de fallas usando la matriz de impedancia de bus.
Un mtodo simple de calcular la corriente de una falla trifsica es utilizando el
teorema de Thvenin, por medio del cual se requiere conocer nicamente el
voltaje en el punto de falla antes de que esta ocurra y los parmetros de la red. El
circuito equivalente de Thvenin que se debe obtener consiste de un generador y
una impedancia conectada entre el generador y el punto del cortocircuito. El
voltaje del generador ( V f ) es el voltaje de fase en el punto de falla justo antes de
que esta ocurra por lo que se le llama voltaje prefalla. Para obtener la impedancia
equivalente de Thvenin ( Zth ) es necesario modelar los elementos de la red para
obtener un unifilar en el que se deben cortocircuitar todas las fuentes de voltaje
como generadores o motores sincrnicos, ya que estos contribuyen con la
corriente de falla. Luego se reduce la red entorno al punto de falla para obtener la
impedancia equivalente Zth . Una vez conocidos estos parmetros es posible
calcular la corriente subtransitoria I f " mediante la Ecuacin 3.1.
I "f
Vf
3.1
Zth
En este clculo no se considera la corriente de carga, es decir se desprecia la
corriente que fluye en las lneas antes de la falla. Esta consideracin implica que
los voltajes en todas las barras son iguales a V f porque de esta manera no habra
diferencia de potencial entre las lneas. Tambin es comn despreciar las
resistencias y capacitancias de carga, adems de los ajustes en los taps de
transformadores. En un sistema de n barras los voltajes en las barras se pueden
calcular por medio de la Ecuacin 3.2. Esta relacin se puede utilizar en el caso de
la falla trifsica si se forma la matriz Zbarra con las impedancias subtransitorias de
los generadores cortocircuitados.
V1 Z11 Z12
V Z
2 21 Z 22
Vn Z n1 Z n 2
Z n3
Z1n I1
Z 2 n I 2
Zn4 I n
3.2
Con el fin de simular una falla trifsica en una barra se cortocircuitan todas las
fuentes de voltaje y se coloca un generador en el punto de falla con un voltaje - V f
. La magnitud de este voltaje es igual al que haba en la barra antes de la falla,
esto con el fin de obtener que la tensin en ese punto es nula durante el
cortocircuito. Al suponer todas las condiciones anteriores la Ecuacin 3.2 se
convierte en la Ecuacin 3.3 que permite calcular los cambios de voltaje al restar
el voltaje durante la falla y voltaje prefalla. Esta ecuacin supone que ocurre un
cortocircuito trifsico en la k -sima barra de la red y que todas las dems tienen
un voltaje V f antes de la falla, por lo que no hay corrientes de carga o prefalla.
Z12
V1 V f V1 Z11
Vk V f V f Z k1
Vk 1 V f Vk 1 Z ( k 1)1 Z ( k 1)2
Vn V f Vn Zu1
Z1n 0
Z kn I "f
Z( k 1) n 0
Zun 0
Z1k
Z kk
Z k ( k 1)
Z( k 1)( k 1)
Zuk
Z u( k 1)
3.3
Durante la falla el voltaje en la k -sima barra ( Vk ) es nulo, entonces de la
ecuacin 3.3 es posible deducir la Ecuacin 3.4 para calcular la corriente de falla
trifsica en la k - sima barra del sistema.
I "f
Vf
3.4
Vkk
Adems de las Ecuaciones 3.3 y 3.5 se deduce la Ecuacin 3.6 que permite
calcular los voltajes en las otras barras. Esta frmula calcula el voltaje en la j sima barra de la red durante la falla.
V j V j V f
Z
V j V f Z jk I "f V f 1 jk ;
Zkk
3.5
jk
3.6
Al conocer los voltajes en las barras durante la falla es posible calcular las
corrientes subtransitorias entre dos barras i y j conectadas a travs de una lnea
de transmisin con impedancia Zij . Esto se logra al sustituir la Ecuacin 3.6 en la
3.7 para as obtener la Ecuacin 3.8.
I "f
Z Z jk
I "f ij I "f ik
Zij
Vi V j
Zij
V f Zik Z jk
Zij Zkk
3.7
3.8
4. Desarrollo.
A continuacin se muestra el algoritmo realizado en Matlab, para anlisis de fallas
mediante zbus.
Este algoritmo permite conocer la corriente de falla en el bus seleccionado, los
voltajes de falla y las corrientes de falla en las lneas.
El archivo (z_datos contiene) las caractersticas del sistema tales como nodo de
envi, nodo de recepcin, resistencia y reactancia. Adems se crea la matriz de
impedancia (z_bus).
clc
clear all;
%datos de impedancias
%
De
A
R
z
=[ 0
1
0
10
0
11
1
2
2
3
2
5
2
6
3
4
3
6
4
6
4
9
4
10
5
7
6
8
7
8
7
11
8
9
X
0.00
0.00
0.00
0.00
0.08
0.04
0.12
0.10
0.04
0.15
0.18
0.00
0.05
0.06
0.06
0.00
0.052
0.20
0.15
0.25
0.06
0.30
0.15
0.45
0.40
0.40
0.60
0.70
0.08
0.43
0.48
0.35
0.10
0.48];
Y = y_bus(z);
%matriz de admitancia del bus
Z_impedancia=inv(Y);
nl=z(:,1);
nr=z(:,2);
R=z(:,3);
X=z(:,4);
ZB = R + j*X;
nbr=length(z(:,1));
nbus=max(max(nl),max(nr));
En el archivo (y_bus) se realiza la matriz de admitancia.
function[Y] = y_bus(z_datos);
nl=z_datos(:,1); nr=z_datos(:,2); R=z_datos(:,3); X=z_datos(:,4);
nbr=length(z_datos(:,1)); nbus=max(max(nl), max(nr));
Z = R + j*X;
%impedancia de la rama
y = ones(nbr,1)./Z; %admitancia de la rama
Y = zeros(nbus,nbus); %matriz de ceros
%Formacion de las admitancias mutuas
for k=1:nbr;
if nl(k) > 0 & nr(k) > 0
Y(nl(k),nr(k)) = Y(nl(k),nr(k)) - y(k);
Y(nr(k),nl(k)) = Y(nl(k),nr(k));
end
end
%Formacion de las admitancias propias
for n = 1:nbus
for k = 1:nbr
if nl(k) == n | nr(k) == n
Y(n,n) = Y(n,n) + y(k);
else, end
end
end
Por ultimo en archivo (fallas) se selecciona el bus a fallar, la impedancia de falla,
se calcula la corriente de falla, los voltajes de falla y las corrientes entre los buses.
clc; clear all;
z_datos
bf=1;
Zf=2;
V0 = ones(nbus, 1);
angv=zeros(nbus,1);
%%
angv = zeros(nbus,1);
%----------------------DETERMINANDO EL NUNERO DE GENERADORES
ngen=0;
for i=1:nbus
if nl(i)==0
ngen=ngen+1;
else,end
end
%---------------------DETERMINANDO LA CORRIENTE DE FALLA
Ifa = V0(bf,1)/(Z_impedancia(bf,bf)+Zf)
%---------------------DETERMINANDO LOS VOLTAJES DURANTE LA FALLA
fprintf('
fprintf('
fprintf('\n
Voltajes de bus durante la falla (PU) \n\n');
Bus
Voltaje
Angulo ');
No.
Mag.
Grados
\n\n');
for i=1:nbus
angv(i) = angv(i) - Z_impedancia(i,bf)*Ifa;
%<<Saadat
%
delV(i) = delV(i) - Zbus(i,fll)/Zbus(fll,fll);
%<<Stevenson
V0(i) = V0(i) + angv(i);
fprintf(' %4.3g
%5.4g
%8.4f \n', i, abs(V0(i)),
angle(V0(i))*180/pi);
%
magV = abs(V);
%
angleV = angle(V)*180/pi;
end
%----------------------DETERMINANDO LAS CORRIENTES DURANTE LA FALLA
fprintf('\n')
fprintf('
Corrientes de linea durante la falla (PU) \n\n');
fprintf('
Bus
Bus
Corriente
Angulo ');
fprintf('\n
Envo
Recepcion
Mag.
Grados
\n\n');
IFL=zeros(nbus,nbus);
for k=1:nbr
if nl(k,1)~=0
if V0(nl(k,1)) > V0(nr(k,1))
IFL(nl(k,1),nr(k,1))=((V0(nl(k,1)))(V0(nr(k,1))))/(j*X(k)+R(k));
fprintf(' %4.3g
%5.3g
%13.4f
%8.4f \n',nl(k,1),
nr(k,1), abs(IFL(nl(k,1),nr(k,1))), angle(IFL(nl(k,1),nr(k,1)))*180/pi);
%
magIFL(k,1)=abs(IFL(nl(k,1),nr(k,1)));
%
angIFL(k,1)=angle(IFL(nl(k,1),nr(k,1)))*180/pi;
else
IFL(nr(k,1),nl(k,1))=((V0(nr(k,1)))(V0(nl(k,1))))/(j*X(k)+R(k));
fprintf(' %4.3g
%5.3g
%13.4f
%8.4f\n',nr(k,1),
nl(k,1), abs(IFL(nr(k,1),nl(k,1))), angle(IFL(nr(k,1),nl(k,1)))*180/pi);
%
magIFL(k,1)=abs(IFL(nl(k,1),nr(k,1)));
%
angIFL(k,1)=angle(IFL(nl(k,1),nr(k,1)))*180/pi;
end
end
end
5. Caso de estudio.
Realizar un anlisis de falla utilizando el algoritmo de z bus para el sistema
mostrado en la (fig. 5.1), fallando los bus 2 y 8 con una impedancia de falla de 0.
Sistema de 11 nodos. (fig. 5.1)
Datos del sistema. (Tabla. 5.1)
DATOS DEL SISTEMA
DE
0
0
0
1
2
2
2
3
3
4
4
4
5
6
7
7
8
A
1
10
11
2
3
5
6
4
6
6
9
10
7
8
8
11
9
R
0
0
0
0
0.08
0.04
0.12
0.1
0.04
0.15
0.18
0
0.05
0.06
0.06
0
0.052
X
0.2
0.15
0.25
0.06
0.3
0.15
0.45
0.4
0.4
0.6
0.7
0.08
0.43
0.48
0.35
0.1
0.48
6. Resultados.
Falla en bus 2.
Corriente de falla:
Ifa= 0.3588 - 6.5286i
Voltaje de bus durante la falla. Tabla. (6.1)
Voltajes de bus durante la falla (PU)
# Bus
1
2
3
4
5
6
7
8
9
10
11
Voltaje magnitud
0.2308
0
0.3025
0.6607
0.1541
0.3452
0.5825
0.5079
0.5693
0.778
0.7014
Angulo de voltaje
0
90
-5.8688
-4.9091
-10.0017
-5.6151
-3.8907
-4.1442
-3.9678
-2.7166
-2.3067
Corrientes de lnea durante la falla. Tabla. (6.2)
Corrientes de lnea durante la falla (PU)
De
1
2
2
2
4
3
4
4
10
5
6
7
11
9
A
2
3
5
6
3
6
6
9
4
7
8
8
7
8
I magnitud
3.8462
0.9744
0.9923
0.7413
0.869
0.1063
0.5702
0.1272
1.5058
0.9923
0.337
0.21
1.2021
0.1272
I angulo
-90
-80.9374
-85.0703
-80.6837
-80.0625
-88.1078
-80.1005
-86.3257
-80.6035
-85.0703
-83.9034
-82.4359
-84.6103
-86.3257
Falla en bus 8.
Corriente de falla:
Ifa= 0.3764 - 3.3105i
Voltaje de bus durante la falla. Tabla. (6.3)
Voltajes de bus durante la falla (PU)
Voltaje magnitud
Angulo de voltaje
0.8082
-1.818
0.7508
-2.5443
0.6882
-1.5987
0.7491
-2.4902
0.7007
-2.3762
0.5454
-1.0194
0.5618
-3.8128
0
0
0.3008
2.4499
0.8362
-1.4547
0.6866
-2.2272
# Bus
1
2
3
4
5
6
7
8
9
10
11
Corrientes de lnea durante la falla. Tabla. (6.4)
Corrientes de lnea durante la falla (PU)
De
1
A
2
I magnitud
0.9697
I angulo
-82.4034
2
2
3
5
0.2053
0.323
-87.8751
-79.9626
2
4
6
3
0.4427
0.1503
-81.6497
-88.4042
3
4
6
6
0.3556
0.3305
-88.0987
-82.3804
4
10
9
4
0.6229
1.1029
-81.3672
-82.6275
5
6
7
8
0.323
1.1274
-79.9626
-83.8944
7
11
8
7
1.582
1.2601
-84.0852
-85.141
0.6229
-81.3672
7. Conclusin.
Realizar un anlisis de falla en un sistema elctrico es de suma importancia
ya que de esto dependern los ajustes de las protecciones asociadas a
este sistema. El mtodo de z bus para el anlisis de falla es una muy
buena opcin ya que es muy rpido a diferencia de realizar un anlisis
mediante el mtodo cortocircuital.
Del presente reporte se puede concluir que realizando un programa en
Matlab que ejecute el procedimiento de anlisis de fallas por el mtodo de
matriz de impedancia de bus se pueden obtener los parmetros de voltaje y
corriente en cada nodo del sistema bajo una condicin de falla.
El mtodo de anlisis de falla utilizando la matriz z bus se utiliza con mayor
frecuencia en sistemas con gran nmero de nodos ya que es ms rpido
que el mtodo cortocircuital. Al realizar un anlisis de falla podemos
conocer la corriente de falla, corrientes de lnea y voltajes de falla en cada
nodo. Cuando se realiza un anlisis de falla por medio de z bus se
recomienda realizar un algoritmo que desempee esta actividad ya que de
manera analtica resulta muy complejo.
8. Bibliografa
Saadat, H. (1998). Analysis Power System . New York: McGraw-Hill.
