MATLAB试验电力系统暂态稳定分析报告

1、标准文案实验三电力系统暂态稳定分析电力系统暂态稳定计算实际上就是求解发电机转子运动方程的初值问题，从而得出8 -t和co-t的关系曲线。每台发电机的转子运动方程是两个一阶非线性的常微分方程。因此，首先介绍常微分方程的初值问题的数值解法。-、常微分方程的初值问题（一）问题及求解公式的构造方法我们讨论形如式（3-1 ）的一阶微分方程的初值问题寸（x）=f（x, y）, axb金（3-1 ）y（xo） =y（）设初值问题（3-1）的解为y（x）,为了求其数值解而采取离散化方法，在求解区间 a,b上取一组节点a =&：:：xi ：- ： xi :二 xi 1 ： : xn =b称hi =xrxi （

2、i =0,1，n1 ）为步长。在等步长的情况下，步长为b _ ah 二n用yi表示在节点xi处解的准确值y（x）的近似值。设法构造序列 ）所满足的一个方程（称为差分方程）yi+ =y +h 阴（xi, yi ,h）（3-2 ）作为求解公式，这是一个递推公式，从（x, y）出发，采用步进方式，自左相右逐步算出y（x）在所有节点X上的近似值yi （i=1,2，n）。在公式（3-2 ）中，为求yi+只用到前面一步的值 yi ,这种方法称为单步法。 在公式（3-2 ）中的yi +由yi明显表示出，称为显式公式。而形如（3-3）y+=yi+h 甲3，yi, yi+, h）（3-3）的公式称为隐式公式，因

3、为其右端中中还包括y书。如果由公式求时，不止用到前一个节点的值，则称为多步法。(3-4)由式（3-1 ）可得dy = f (x, y)dx两边在x , xf上积分，得大全标准文案xi . 1-y（Xi+） =y（x。f （x, y（x） dx （3-5）Xi由此可以看出，如果想构造求解公式，就要对右端的积分项作某种数值处理。这种求解公式的构造方法叫做数值积分法。（二）一般的初值问题的解法.欧拉法和改进欧拉法对于初值问题（3-1）,采用数值积分法，从而得到（3-5）。对于（3-5）右端的积分用矩形公式（取左端点），则得到Xi 1f （x, y（x）dx : h f （Xi, y（xj Xi进而得

4、到（3-1 ）的求解公式（3-2 ） TOC o 1-5 h z yi + =yi +h f（Xi,yi）（ i =0, 1, 2, n-1 ）（3-6）此公式称为欧拉（Euler ）格式。如果对式（3-5）右端的积分用梯形公式X X1hf (x, y(x)dx : - (f(Xi, y(Xi) f(Xi 1, y(Xi 1)则可以得到初值问题（3-1 ）的梯形求解公式如式（3-7 ）yi+=yi +2 区）+f （为十，yq （ i =0, 1, 2, n-1 ）（3-7）式（3-7）是个隐式公式。可以采取先用欧拉格式求一个y（x+）的初步近似值，记作 与书，称之为预报值，然后用预报值yi+

5、替代式（3-7）右端的yi+,再计算得到y ,称之为校正值，这样建立起来的预报-校正方法称为改进欧拉格式(3-8)y 1 =yi - h f (Xi, yi)yi 1 =yi h f (Xi, yi) - f (x 1 ,yi 1) I22.龙格一库塔方法在单步法中，应用最广泛的是龙格库塔（ Runge-kutta ）法，简称R- K法。下面直接给出一种四阶的龙格一库塔法的计算公式（3-9 ）大全标准文案1Va =yi +(Ki +2K2 +2K3 +K4) 一 6(3-9) TOC o 1-5 h z Ki =h f (Xi ,yi) h1水2 =h f(Xi +-, yi +- Ki) 2

6、2. h 1K3 =h f (x +-,yi +-K2) 22K4 =h f(Xi - h,yi - K3)它也称为标准(古典)龙格-库塔法。例3-1研究下列微分方程的初值问题 HYPERLINK l bookmark36 o Current Document *1C 2jy =2y1 +x y(0) =o解：这是一个特殊的微分方程，其解的解析式可以给出，为应用龙格库塔法，取 h =0.25 ,根据式(3-9)编写一段程序，由零开始自左相右逐步算出y(x)在所有节点X上的近似值y。计算结果见表3-1。计算结果表明，四阶龙格库塔方法的精度是较高的。表3-1Xnyny(Xn ) -yn2.00.3

7、99956994.3e-54.00.235291592.5e-66.00.162161793.7e-78.00.123076839.2e-83-9)编写较复杂的实际上，MATLAB为常微分方程提供了很好的解题指令，使得求解常微分方程变得很容易，并且能将问题及解答表现在图形上。因此，我们可以不用根据式(程序，而只需应用MATLABt供的常微分方程解题器来解决问题。下面给出用MATLA面写的解题程序。首先编写描述常微分方程的OD或件，文彳名为 myfun，便于解题器调用它。function dy = myfun(x,y) dy = zeros(1,1);dy=1/(1+xA2)-2*yA2;再编写

8、利用解题器指令求解y的程序。clearx0=0;for i=1:4xm=2*i;大全标准文案y0=0;x,y = ode45( myfun ,x0 xm,y0);format longy(length(y) endplot(x,y,-)运行上述程序，在得到几个点的函数值的同时，也得到函数y的曲线，如图3-1所示。图3-1 根据运算结果画出y的曲线简单电力系统的暂态稳定性(一)物理过程分析某简单电力系统如图3-2(a)所示，正常运行时发电机经过变压器和双回线路向无限大系统供电。发电机用电势E作为其等值电势，则电势E 与无限大系统间的电抗为x t=x； +xT1+xT2(3-10 )2这时发电机发

9、出的电磁功率可表示为P=EUsin6 =Pp sin6(3-11)x如果突然在一回输电线路始端发生不对称短路，如图3-2(b)所示。故障期间发电机电势E ,与无限大系统之间的联系电抗为Xlx(XdXti)( Xt2)xq=(Xd +x)十(二+Xt2)+2 (3-12)2x在故障情况下发电机输出的电磁功率为P=EUsin6 =Pp sin 6(3-13)x在短路故障发生之后，线路继电保护装置将迅速断开故障线路两端的断路器，如图3-2(c)所示。此时发电机电势E与无限大系统间的联系电抗为大全标准文案x = Xd XT1 Xl XT2(3-14 )发电机输出的功率为sin :. pPmM sin、

10、.(3-15)Plll EU X IE jXd jXT1j jXLjXT2 UjXLjXT2jXLjX:(b)(c)图3-2简单电力系统及其等值电路（a）正常运行方式及其等值电路；（b）故障情况及其等值电路;(c)故障切除后及其等值电路如果正常时发电机向无限大系统输送的有功功率为P0,则原动机输出的机械功率Pt等P0。因此，可3-3所示。于R。假定不计故障后几秒种之内调速器的作用，即认为机械功率始终保持 以得到此简单电力系统正常运行、故障期间及故障切除后的功率特性曲线如图PrnpP.Pt =P0、0Qk:;c0m图3-3 简单系统正常运行、故障期间及故障切除后的功率特性曲线但是，实际工作需对于

11、上述简单电力系统， 我们可以根据等面积定则求得极限切除角。要知道在多少时间之内切除故障线路，也就是要知道与极限切除角对应的极限切除时间。解决这个问题，必须求解发电机的转子运动方程。（二）求解发电机的转子运动方程大全标准文案求解发电机转子运动方程可以得出8-t和co-t的关系曲线。其中8 -t曲线一般称为摇摆曲线。在上述简单电力系统中故障期间的转子运动方程为5：,、(1-1);，1.(3-16)d ,1 ,=(Pt -P M sin、.)dt Tj式中，二一一功率角，其单位为弧度；.一一转子角速度，标幺值；i 转子的同步角速度，即缶1= 2日=314.16 ,其单位为弧度/秒；Tj 发电机的惯性

12、时间常数，其单位为秒;R、p 分别为机械和电磁功率，标幺值。这是两个一阶的非线性常微分方程，它的起始条件是已知的，即1 Ptt =to =0,= = coo =1.0 . = = 60 = sin Pm故障切除后，由于系统参数改变，以致发电机功率特性发生变化，必须开始求解另一组 微分方程：d、.( -1),1.(3-17)d 1愉二T；(Pt -p m sin、)式中变量含义同前述，其中 P型也为标幺值。这组方程的起始条件为t=tc； 0 = 0c； CO = Oc其中tc为给定的切除时间；&、a为与tc时刻对应的6和6 ,它们可由故障期间的8 -t和CO -t的关系曲线求得(6和8都是不突变

13、的)。一般来说，在计算故障发生后几秒种的过程 中，如果8始终不超过 180o ,而且振荡幅值越来越小，则系统是暂态稳定的。当发电机与无限大系统之间发生振荡或失去同步时，在发电机的转子回路中， 特别是阻尼绕组中将有感应电流而形成阻尼转矩(也称为异步转矩)。当作微小振荡时，阻尼功率可表达为：Pd = D/ = DQ1)(3-18)式中，d称为阻尼功率系数；4为转子角速度的偏移量，标幺值；CO为转子角速度，标幺值。阻尼功率系数 d除了与发电机的参数有关外，还和原始功角、尴的振荡频率有关。在一般情况下它是正数。 在原始功角较小，或者定子回路中有串联电容使定子回路总电阻相对 于总电抗较大时，D可能为负数

14、。如果考虑阻尼功率的影响，则故障后的转子运动方程又可 表达为大全标准文案5：二(.：. 一1)，1(3-19 )dtPt DQ，_1)_P m sin c.dt Tj3-16 );另电力系统暂态稳定计算包括两类问题，一类是应用数值计算法得出故障期间的曲线后, 根据曲线找到与极限切除角对应的极限切除时间，此时只需要求解微分方程(类是已知故障切除时间，需要求出摇摆曲线来判断系统的稳定性，此时需要分段分别求解微 分方程(3-16)和(3-17)。如果考虑阻尼转矩的影响，则此时需要分段分别求解微分方程(3-16 )和(3-19 )。三、例题例3-2 某简单电力系统如图3-4所示，取基准值Sb=220M

15、VA Ub=209KV。换算后的参数 已经标在图中，其中一回线的电抗Xl =0.486, Tj =8.18秒。设电力线路某一回的始端发生两相接地短路。假定 E=常数。(1)计算保持暂态稳定而要求的极限切除角。(2)计算极限切除时间，并且作出在 0.15秒切除故障时的8 -t曲线。GT1Lx2 =0.432x =0.486XT1 =0.138LXd =0.295Xlo =4XL图3-4某简单电力系统的接线图解：计算系统正常运行方式，决定 E 和瓦。由3-3(a)的正序网络可得，此时系统的总电抗为x =0.295+0.138+0.243+0.122=0.798发电机的暂态电势为：0.2 0.798

16、 21.0 0.798 2 E = (1.0 .)() =1.41,1.01.0、.=tg一口798=34.53 o1.0 0.2 0.798(2)故障后的功率特性又由3-3(b)的负序、零序网络可得故障点的负序、零序等值电抗为(0.432 0.138) (0.243 0.122)X2 丁= =0.222一 (0.432 0.138) (0.243 0.122)大全标准文案_ 0.138(0.972 +0.122) =0 1230.138 (0.972 0.122)所以在正序网络故障点上的附加电抗为:=0.0790.222x0.123必一0.222 +0.123于是故障时等值电路如图3-3(c

17、)所示，则x =0.433 0.3650.433 0.365_0.079= 2.80因此，故障期间发电机的最大功率为:EUx .1.41 1.02.8= 0.504(3)故障切除后的功率特性故障切除后的等值电路如图3-3(d)所示x -0.295 0.138 0.486 0.122=1.041此时最大功率为P MEU1.41 1.0 =1.35XI1.0410. j, 1.0，、，h -180 -sin1.35=132.20E=1.41 j 0.295 j0.138 j 0.243 j 0.122U =1.0、0 =34.53f(1)(a)j 0.432 j 0.138 j 0.243F0 =

18、1.0Q0 =0.2 j 0.122f(2)j 0.138 j0.243 4 j 0.122f(0)U =1.0 j 0.122I(b)EV.41 j 0.295 j 0.138 j 0.243,j 0.079U = 1.0(c)El.41j 0.295 j 0.138 j 0.486 j 0.122q(d)图3-5 例题7-12的等值电路(a)正常运行等值电路；(b)负序和零序等值电路；(c)故障时等值电路；(d)故障切除后等值电路(4)计算极限切除角大全标准文案COS CcmPt ( ：0 -、)P，m COS.：,h -Pm COS%P |MP M1.0 (132.2 -34.53) 1

19、.35cos132.2 _0.504cos34.531801.35 -0.504=0.458、.cm =62.740(5)找出极限切除时间tcm根据(3-16),首先计算初值34.536。=* n=0.6027,切。=1.0180令y(1)= S, y(2)= 0。编写描述故障期间转子运动方程的OD政件，文彳名为 myequ，。function dy = myequ(t,y)dy = zeros(2,1);f=50;w1=2*pi*f;dy(1) = (y(2)-1)*w1;dy(2) = (1/8.18)*(1.0-0.504*sin(y(1);再编写利用解题器指令求解y的程序。cleart

20、0=0;tm=0.25;d0=(34.53/180)*pi;w0=1;T,Y = ode45( myequ ,t0 tm,d0 w0);plot(T,(Y(:,1)/pi)*180,- ,0.194,62.76,*)text(0.194,60, delta_cmax=62.76circ , FontSize ,10)text(0.194,56,t_cmax=0.194s, FontSize ,10)图3-6 例题7-12的S-t曲线图3-6给出短路发生后0秒到0.25秒期间的8 -t计算曲线，根据最大切除角6cm大全标准文案(=62.740)找到极限切除时间tcm为0.194秒。由图3-6可见

21、，如果故障切除时间大于0.194秒，则发电机的功角将不断地增大，最终失去暂态稳定。在极限切除时间之前切除故障，发电机的摇摆曲线的状况将在下面作计算、分析。(6)不考虑阻尼转矩的影响，当故障切除时间为 0.15秒时通过计算得出8 -t曲线首先编写描述故障期间转子运动方程的OD或件，文彳名为“ myfun01”。function dy = myfun01(t,y)f=50; w1=2*pi*f;TJ=8.18; Pt=1.0; P2m=0.504;dy = zeros(2,1);dy(1) = (y(2)-1)*w1;dy(2) = (1/TJ)*(Pt-P2m*sin(y(1);再编写描述故障切

22、除后转子运动方程的OD或件，文彳名为“ myfun02”。function dy = myfun02(t,y)f=50; w1=2*pi*f;TJ=8.18; Pt=1.0; P3m=1.35;dy = zeros(2,1);dy(1) = (y(2)-1)*w1;dy(2) = (1/TJ)*(Pt-P3m*sin(y(1);编写利用解题器指令求解y的小程序。cleart0=0; tc=0.15; tm=2.0;d0=(34.53/180)*pi; w0=1.0;T1,Y1 = ode45( myfun01 ,t0 tc,d0 w0);dc=Y1(length(Y1),1);wc=Y1(le

23、ngth(Y1),2);T2,Y2 = ode45( myfun02 ,tc tm,dc wc); TOC o 1-5 h z plot(T1,(Y1(:,1)/pi)*180,- ,T2,(Y2(:,1)/pi)*180,- ,tc,(dc/pi)*180,*)text(0.28,50,itt_c=0.15s, FontSize ,8)text(0.28,43,itdelta_c=51.71circ, FontSize ,8)xlabel( itt)ylabel( itdelta )计算结果表明，功角6沿着故障切除后的功角特性曲线根据等面积定则作等幅振荡，如图3-7所示。实际上，由于阻尼转矩的影响，振荡的幅度是逐渐衰减的，功角 6最终运行在6k =47.8 o 。因此，发电机能够保持暂态稳定。大全标准文案图3-7不考虑阻尼转矩