实验2+故障分析实验-学生版

1、PAGE PAGE 19l实验2 故障分析实验 一、实验目的1、熟悉利用Matlab/Simulink/PSB建立电力系统仿真模型方法。2、理解掌握同步发电机突然短路的暂态过程。二、实验设备Matlab/Simulink/PSB三、实验原理 同步发电机是电力系统中最重要和最复杂的元件，它由多个有磁耦合关系的绕组构成，定子绕组同转子绕组之间还有相对运动，同步电机突然短路的暂态过程要比稳态对称运行（包括稳态对称短路）时复杂得多。稳态对称运行时，电枢磁势的大小不随时间变化，而且在空间以同步速度旋转，它同转子没有相对运动，因此不会在转子绕组中感应电流。突然短路时，定子电流在数值上发生急剧变化，电枢反应

2、磁通也随着变化，并在转子绕组中产生感应电流，这种电流又反过来影响定子电流的变化。定子和转子绕组电流的互相影响是同步电机突然短路暂态过程的一个显著特点。1、同步发电机突然短路的暂态过程 下面的讨论均假设在暂态过程期间同步发电机保持同步转速以及在短路后励磁电压保持不变。 在分析同步发电机突然三相短路暂态过程时，可以利用替加原理，这样同步发电机机端突然短路相当于在发电机端口处突然加上了与电机短路前的端电压大小相等但方向相反的二相电压：在定子绕组上突然加以对称的相电压后，为了保持其无源闭合电路的磁链不变，在其定子绕组中将要引起相应的瞬变电流，而且这些瞬变电流还要按照一定的时间常数逐步衰减至稳态值。 当

3、发电机突然短路时，定子各绕组电流将包含基频分量、倍频分量和直流分量。到达稳态后，定子电流起始值中的直流分量和倍频分量将由其起始值衰减到零，而基频分量则由其起始值衰减为相应的稳态值。同样，在转子绕组中也包含直流分量和同频率交流分量。 引入衰减因子以后，定子电流的d轴和q轴分量分别为经过变换和整理，可得定子a相电流为转子绕组中的电流 在前面几式中，xd、xq为定子绕组纵轴、横轴的同步电抗；xf为纵轴绕组之间的电枢反映电抗；xad、xaq为发电机转子纵轴、横轴的电抗；xD、xQ表示D、Q阻尼绕组的电抗，xD为D阻尼绕组的漏抗；xd、xd分别为纵轴暂态电抗、次暂态电抗；xq为横轴次暂态电抗；Eq、Eq

4、分别为横轴暂态电势、次暂态电势；Ed为纵轴次暂态电势；Eq0、V0为短路前瞬间的空载电势、机端电压。2、对称三相短路 考虑同步电机在同步频率下运行，并且励磁不变。电枢末端处发生三相短路时，出现三个电枢电流和一个励磁电流。假定电机为空载则有励磁电流的初始值为当电枢末端处发生三相短路时，有有由于i0=0，当发生三相短路时，转子参照系下的电机电压方程式为上述方程为状态空间形式，写成下面的紧凑形式为或者如果转速为常数，那么状态空间方程式是线性的，可以通过拉氏变换来求其解析解。然而，Matlab具有强大的功能，可以模拟同步电机的矩阵形式的非线性微分方程式，若考虑转速变化，需要同步电机的动态方程式，这是一

5、个二阶微分方程即摆动方程。摆动方程也可以写成状态空间的形式，如写成两个一阶微分方程。由于转速的变化对故障后电流的瞬时值影响很小，所以，这里可以忽略。求出直轴和交轴电流之后，相电流就可以通过派克反变换求出由于i0=0可得相电流为Matlab中提供两个M文件，ode23和ode45用来对微分方程求数值解，用龙格库塔芬尔格数值积分方法。函数ode23求解中等精度的2阶和3阶的常微分方程组，函数ode45求较高精度的4阶和5阶常微分方程。3、不对称短路同步电机的大多数故障类型为相间短路故障和单相接地短路故障，分析这样的不平衡故障比较困难。d-q-0模型对不平衡故障不是很适用，需要更进一步的转换。这样，

6、最后的解析解会很复杂，而且是近似解。（1）相间短路当b相和c相发生直接短路时，vb=vc=0和ib=-ic由于a相没有短路，且最初a相空载，即ia=0，因此va=0直轴和交轴的电流的导数为发生相间故障时dq0-abc量的电压方程为4、接地短路设同步发电机的相对地短路va=0，短路前空载ib=ic=0，一种求发电机电压方程的简便方法是用dq0-abc量，结合短路条件和d-q-0量的电感，得到相接地短路的电压方程式为其中四、实验内容1、同步发电机突然三相短路暂态过程的数值计算 在已知发电机参数的情况下，可以利用Matlab对突然三相短路后的定子电流、转子电流暂态过程表达式进行数值计算分析这样将有助

7、于更好地理解短路的物理过程。 假设一台有阻尼绕组同步发电机，PN=200MW，UN=13.8kV，fN=50Hz，xd=1.0，xq=0.6，xd=0.30，xd=0.21，xq=0.31，r=0.005，xf=0.18，xaD=0.1，xQ=0.25，Td0=5s，TD=2s，Tq0=1.4s。若发电机空载，端电压为额定电压，端子突然发生三相短路，且0=0，利用Matlab对突然三相短路后的定子电流进行数值计算的基本步骤如下： 首先计算各衰减时间常数。 Ta=0.16s，Td=0.72s，Tq=0.34s，Td=1.64s。 由于空载时Eq0=Eq0=Eq0=V0=1，Ed0=0，0=0，则

8、可得a相定子电流表达式为 利用Matlab对上式进行数值计算并绘图的m文件程序清单如下：%发电机出口空载短路计算曲线、仿真图绘制N=48;t1=(0:0.02/N:1.00);m=size(t1);fai=0*pi/180;%空载短路全电流表达式Ia=(-cos(2*pi*50*t1+fai)-1.43*exp(-2.97*t1).*cos(2*pi*50*t1+fai)-2.34*exp(-0.608*t1).*cos(2*pi*50*t1+fai)+. 4*exp(-6.3*t1).*cos(-fai*pi/180)+0.77*exp(-6.3*t1).*cos(2*2*pi*50*t1+

9、fai);%基频分量Ia1=-cos(2*pi*50*t1+fai)-1.43*exp(-2.97*t1).*cos(2*pi*50*t1+fai)-2.34*exp(-0.608*t1).*cos(2*pi*50*t1+fai);%倍频分量Ia2=0.77*exp(-6.3*t1).*cos(2*2*pi*50*t1+fai);%非周期分量Iap=4*exp(-6.3*t1).*cos(-fai*pi/180);subplot(4,1,1);plot(t1,Ia);grid on;axis(0 1 -10 10);ylabel(Ia(p.u.);subplot(4,1,2);plot(t1,

10、Ia1);grid on;axis(0 1 -10 10);ylabel(Ia1(p.u.);subplot(4,1,3);plot(t1,Ia2);grid on;axis(0 1 -1 1);ylabel(Ia2(p.u.);subplot(4,1,4);plot(t1,Iap);grid on;axis(0 1 -10 10);ylabel(Iap(p.u.);xlabel(t/s);请根据以上程序绘制发电机端突然发生三相短路时的特征波形（A相定子电流，基频分量、倍频分量和非周期分量的波形）。注：当发电机端突然发生三相短路时，定子电流中的倍频分量是很小的(在实用计算中常忽略不计)，为了能

11、够在图中表示清楚，将其纵坐标的值取为-1，1，波形相应地被放大。 2、同步发电机突然三相短路暂态过程的仿真方法 针对以上的发电机参数，建立其Matlab/Simulink仿真模型如下图所示。 在图中，同步发电机采用p.u.标准同步电机模块，根据前面的计算，设置同步发电机采用p.u.标准同步电机模块参数。升压变压器T采用“Three-phase transformer(Two Windings)”模型，设置“Three-phase transformer(Two Windings)”模型参数。 由于同步发电机模块为电流源输出，因此在其端口并联了一个有功功率为5MW的负荷模块。 仿真开始前，要利用

12、Powergui模块对电机进行初始化设置。单击Powergui模块，打开潮流计算和电机初始化窗口，设置参数如图下所示。图中设定同步发电机为平衡节点“Swing bus”。初始化后，与同步发电机模块输入端口相连的两个常数模块Pm和Vf以及 “Init.Cond.”将会自动设置。设置故障时间，其他参数采用默认设置。问题：如何设置故障发生时间？ 选择Ode15s算法，仿真的结束时间取为1s。开始仿真，得到发电机端突然三相短路后的三相定子电流波形图（请绘制仿真波形）。%*绘制短路电流仿真图x=ScopeData1.time;subplot(3,1,1);plot(x,ScopeData1.signal

13、s.values(:,1);grid on;axis(0 1 -10 10);ylabel(Ia(p.u.)subplot(3,1,2);plot(x,ScopeData1.signals.values(:,2);grid on;axis(0 1 -10 10);ylabel(Ib(p.u.)subplot(3,1,3);plot(x,ScopeData1.signals.values(:,3);grid on;axis(0 1 -10 10);ylabel(Ic(p.u.)xlabel(t/s);问题：a相定子电流的冲击电流标幺值和理论计算值存在_0.91_的误差？%*绘制短路时id、iq、

14、励磁电流仿真图x=Pm3.time;subplot(3,1,1);plot(x,Pm3.signals.values);grid on;axis(0 1 -10 10);ylabel(id(p.u.)subplot(3,1,2);plot(x,Pm2.signals.values);grid on;axis(0 1 -10 10);ylabel(iq(p.u.)subplot(3,1,3);plot(x,Pm1.signals.values);grid on;axis(0 1 0 3);ylabel(if(p.u.)xlabel(t/s);请根据程序绘制ua波形、短路后定子电流的d轴和q轴分量

15、id、iq以及励磁电流if的仿真波形图。改变故障模块中的短路类型，就可以仿真同步发电机发生各种不对称短路时的故障情况。例如，设置在0.02025s时发生BC两相短路故障。请绘制得到发电机端突然两相短路后的三相定子电流仿真波形图。3、对称三相短路一台同步电机空载运行，额定容量为500MVA，额定工作电压为30kV，频率为60Hz，励磁电压为400V，在电枢末端处发生三相短路。假定短路时转子直轴与a相轴重合，并且假定转子转速不变，为同步转速。电机的参数如下：求相电流和励磁电流的暂态波形（根据程序绘制波形）。直流励磁电压为VF=400V，励磁电流的初始值为由于电机最初为空载，所以% symshort

16、.mfunction iprime =symfault(t,i)f=60.; w=2.*pi*f;d=0.; d=d*pi/180.; theta=w*t+d +pi/2;% Parameters of a 160 MVA, 15 kV Synchronous MachineLF = 2.500; LD = 0.0068; LQ = 0.0016; Ld = 0.0072; Lq = 0.0070;MF = 0.100; MD = 0.0054; MQ = 0.0026; MR = 0.1250;ra = 0.002; rF = 0.4000; rD = 0.015; rQ = 0.0150;

17、VF = 400; % DC field voltageV = 0; -VF; 0; 0; 0; % Voltage column vectork=sqrt(1.5);R= ra 0 0 w*Lq w*k*MQ 0 rF 0 0 0 0 0 rD 0 0 -w*Ld -w*k*MF -w*k*MD ra 0 0 0 0 0 rQ; L =Ld k*MF k*MD 0 0 k*MF LF MR 0 0 k*MD MR LD 0 0 0 0 0 Lq k*MQ 0 0 0 k*MQ LQ;Li=inv(L);iprime=- Li*V - Li*R*i;%*fprintf(n Wait simul

18、ation is in process. n)VF = 400; rF = 0.4; iF0 = VF/rF;f = 60; w=2.*pi*f;deg = 0; d=deg*pi/180;t0 = 0 ; tfinal = 0.80;i0 = 0; iF0; 0; 0; 0 ; % Initial currents%t,i = ode45(symshort,t0,tfinal,i0); % use for Matlab 4tspan = t0, tfinal; % use for Matlab 5t, i = ode45(symshort, tspan, i0); % use for Mat

19、lab 5theta = w*t+d +pi/2;ia = sqrt(2/3)*(cos(theta).*i(:,1)+sin(theta).*i(:,4);ib = sqrt(2/3)*(cos(theta-2*pi/3).*i(:,1)+sin(theta-2*pi/3).*i(:,4);ic = sqrt(2/3)*(cos(theta+2*pi/3).*i(:,1)+sin(theta+2*pi/3).*i(:,4);iF = i(:,2); % ia = Armature current iF=Field currentfigure(1), plot(t,ia), xlabel(t,

20、 s), ylabel(i_a, A), gridtitle(Three-phase short circuit i_a, ,delta = , num2str(deg),circ)figure(2), plot(t,ib), xlabel(t, s), ylabel(i_b, A), gridtitle(Three-phase short circuit i_b, ,delta = , num2str(deg),circ)figure(3), plot(t,ic), xlabel(t, s), ylabel(i_c, A), gridtitle(Three-phase short circu

21、it i_c, ,delta = , num2str(deg),circ)figure(4), plot(t,iF), xlabel(t, s), ylabel(i_F, A), gridtitle(Three-phase short circuit i_F, ,delta = ,num2str(deg),circ)4、不对称短路 实验内容3中同步发电机空载运行，常值励磁电压为400V。在电枢末端b相和c相之间发生相间短路。假定短路时刻转子直轴与a相轴重合，并且假定转子转速保持同步不变。利用ode45算法求b相电流和励磁电流的波形（根据程序绘制波形）。% llshort.mfunction i

22、prime =llshort(t,i)f=60.; w=2.*pi*f;d=0.; d=d*pi/180.; theta=w*t+d +pi/2;% Parameters of a 500 MVA, 30 kV Synchronous MachineLF = 2.500; LD = 0.0068; LQ = 0.0016; Ld = 0.0072; Lq = 0.0070;MF = 0.100; MD = 0.0054; MQ = 0.0026; MR = 0.1250;ra = 0.002; rF = 0.4000; rD = 0.015; rQ = 0.0150;VF = 400; % d

23、c field voltageV = -VF; 0; 0; 0; % Voltage column vectorK=sqrt(1.5); RT2=sqrt(2.0);R=RT2*K*w*MF*cos(theta) rF 0 0 RT2*K*w*MD*cos(theta) 0 rD 0 RT2*K*w*MQ*sin(theta) 0 0 rQ RT2*(ra+w*(Ld-Lq)*sin(2*theta) K*w*MF*cos(theta) K*w*MD*cos(theta) . K*w*MQ*sin(theta);L=RT2*K*MF*sin(theta) LF MR 0 RT2*K*MD*si

24、n(theta) MR LD 0 -RT2*K*MQ*cos(theta) 0 0 LQ RT2*(Ld*sin(theta)2+Lq*cos(theta)2) K*MF*sin(theta) K*MD*sin(theta) . -K*MQ*cos(theta);Li=inv(L);iprime=-Li*V - Li*R*i;%*fprintf(n Wait simulation is in process. n)VF = 400; rF = 0.4; iF0 = VF/rF;f = 60; w = 2.*pi*f;deg = 0; d = deg*pi/180;t0 = 0 ; tfinal

25、 = 0.80;i0 = 0; iF0; 0; 0; % Initial currents%t,i = ode45(llshort, t0, tfinal, i0); % use for Matlab 4tspan = t0, tfinal; % use for Matlab 5t,i = ode45(llshort, tspan, i0); % use for Matlab 5ib=i(:,1); iF=i(:,2); % ib = Armature current iF=Field currentfigure(1), plot(t,ib), xlabel(t, s), ylabel(i_b

26、, A), gridtitle(Line-line short circuit i_b, ,delta = , num2str(deg),circ)figure(2), plot(t,iF), xlabel(t, s), ylabel(i_F, A), gridtitle(Line-line short circuit i_F, ,delta = ,num2str(deg),circ)5、接地短路 实验内容4同步发电机空载运行，常值励磁电压为400V。在电枢末端发生a相接地短路。假定短路时转子直轴与a相轴重合，并且假定转子转速保持在同步转速不变。利用ode45算法求a相电流和励磁电流的波形（根

27、据程序绘制波形）。% lgshort.mfunction iprime =lgshort(t,i)f=60.; w=2.*pi*f;d=0.; d=d*pi/180.; theta=w*t+d +pi/2;% Parameters of a 500 MVA, 30 KV Synchronous MachineLF = 2.500; LD = 0.0068; LQ = 0.0016; Ld = 0.0072; Lq = 0.0070;MF = 0.100; MD = 0.0054; MQ = 0.0026; MR = 0.1250; L0 =0.002;ra = 0.002; rF = 0.4000; rD = 0.015; rQ = 0.0150;VF = 400; % DC field voltageV = 0; -VF; 0; 0; % Voltage column vectorK=sqrt(1.5);Ls=(L0+Ld+Lq)/3; Lm=(Ld-Lq)/3;R=ra-2*w*Lm*sin(2*theta) -w*MF*sin(theta) -w*MD*sin(theta) w*MQ*cos(theta) -w*MF*sin(theta) rF 0 0 -w*MD*sin(theta) 0 rD 0 w*MQ*cos(t