电力系统短路潮流计算

1、2012年6月15日广州学院guangzhou college of south china university of technology电力系统分析课程报告题目：基于mat lab电力系统短路潮流计算系 别 电气工程系 专业班级 学生姓名 学生学号 指导教师 提交日期目录短路的原因及其后果1. 1短路的原因312 短路的后果3二. 短路类型2. 1对称短路-三相短路 32. 2不对称短路3三. 电力系统短路潮流计算的意义和冃的3. 1短路潮流计算的意义43.2短路潮流计算的目的5四. 短路潮流计算的基本原理五. 短路计算机计算的流程5.1 短路计算流程图65.2 短路计算源程序图75.3

2、 运行计算结果及分析12六总结七.参考文献 八.附录一、短路的原因及其后果1.1短路的原因1、电气设备及载流导体因绝缘老化，或遭受机械损伤，或因雷击、过电 压引起绝缘损坏；2、架空线路因大风或导线履冰引起电杆倒塌等，或因鸟兽跨接裸露导体 等；3、电气设备因设计、安装及维护不良所致的设备缺陷引发的短路；4、运行人员违反安全操作规程而误操作，如带负荷拉隔离开关，线路或 设备检修后未拆除接地线就加上电压等。1.2短路的后果1、强人的短路电流通过电气设备使发热急剧增加，短路持续吋间较长时, 足以使设备因过热而损坏甚至烧毁；2、巨大的短路电流将在电气设备的导体间产生很大的电动力，可能使导 体变形、扭曲或

3、损坏；3、短路将引起系统电压的突然大幅度下降，系统中主要负荷异步电动机 将因转矩下降而减速或停转，造成产品报废其至设备损坏；4、短路将引起系统中功率分布的突然变化，可能导致并列运行的发电厂 失去同步，破坏系统的稳定性，造成大面积停电。这是短路所导致的最严 重的后果。5、巨大的短路电流将在周围空间产牛很强的电磁场，尤其是不对称短路 吋，不平衡电流所产生的不平衡交变磁场，对周围的通信网络、信号系 统、晶闸管触发系统及自动控制系统产生干扰。二、短路类型2. 1、对称短路一一三相短路三相电流和电压仍是对称的2.2、不对称短路:两相短路:小1)单相接地短路:两相短路接地:相间短路与接地短路:相间短路：三

4、相短路、两相短路接地短路：单相接地短路、两相短路接地三、电力系统短路潮流计算的意义和目的3. 1短路潮流计算的意义短路潮流计算是电力系统的一项重要分析功能，是进行故障计算，继电保护 整定，安全分析的必要工具。matlab fl 1980年问世以来，它的强大的矩阵处理功能给电力系统的分析、 计算带来许多方便。在处理短路潮流计算吋，其计算机软件的速度已无法满足大 电网模拟和实时控制的仿真要求，而高效的短路问题相关软件的研究已成为大规 模电力系统仿真计算的关键。随着计算机技术的不断发展和成熟，对matlab短 路潮流计算的研究为快速、详细地解决大电网的计算问题开辟了新思路。短路潮流计算是电力系统分析

5、屮的一种基本的计算，其任务时要根据给定的 网络接线和其他已知条件，计算网络中的功率分布、功率损耗和未知的节点电压。 采用计算机计算，可以减轻计算人员的工作量，提高计算准确度和计算速度。短路潮流计算是电力系统分析中的一种基本的计算，其任务时要根据给定的 网络接线和其他已知条件，计算网络中的功率分布、功率损耗和未知的节点电压。 采用计算机计算，可以减轻计算人员的工作量，提高计算准确度和计算速度。32短路潮流计算的目的1、为选择和校验各种电气设备的机械稳定性和热稳定性提供依 据，为此，计算短路冲击电流以校验设备的机械稳定性，计算短路电 流的周期分量以校验设备的热稳定性；2、为设计和选择发电厂和变电所

6、的电气主接线提供必要的数据;3、为合理配置电力系统中各种继电保护和自动装置并正确整定 其参数提供可靠的依据。短路潮流计算的基本原理短路电流计算的关键是化简网络以求得各电源点相对短路点的转移阻抗。那 么根据上面介绍的数学模型，在程序设计时，关键也就在于如何化简节点导纳矩 阵，求得各电源节点相对于短路点的互导纳，进而求得转移电抗，查曲线数字表 求得短路电流。对总节点数为n,电源节点数为m的网络，设其节点导纳矩阵为 y,我们可以利用高斯消元法消去n-m- 1个中间节点和联络节点，将其化简 为仅含有m个电源节点和一个短路节点的等效网络。化简后的矩阵为转移导纳矩 阵，转移导纳矩阵仍然是对称矩阵，其阶数为

7、m+ 1阶，矩阵中各元素的意义与 化简前相同，b|j：对角线元素y y i i称为节点i的自导纳，其值等于接于节点 i的所有支路导纳之和。非对角线元素yy i j称为节点i、j间的互导纳，它 等于直接联接于节点i、j间的支路导纳的负值。所以各电源节点与电源节点之 间的转移阻抗就是转移导纳矩阵中对应的互导纳的倒数，于是便可求得各电源供 给的短路电流。各电源供给的短路电流之和就是短路点的总短路电流。这一系列 过程都可以由微机自动完成。其步骤如下。1 ）准备工作：将网络中的各节点用从1起的正整数编号，先电源节点后其 它节点。2）支路阻抗计算：各支路按其元件性质不同分别计算阻抗。3）节点导纳矩阵形成：

8、根据数学模型中介绍的节点导纳矩阵中各元素的意义，可由上一步中的各支路阻抗值形成节点导纳矩阵。4 ）节点导纳矩阵化简：利用高斯消元法进行化简。5 ）查曲线：运算曲线以曲线数字表的形式存放在磁盘中。应用时按电源的 性质（水轮发电机、汽轮发电机或系统）区分，调用不同的运算曲线。如果数值 介于两数之间，则用插值法求结果。实践证明，采用线性插值得到的结果已经能 够满足一般的计算要求。五、短路计算机计算的流程5. 1短路计算流程图v用公式（69）计算指定支路电流输出结果5. 2短路计算源程序图sbase_mva=100.fid=fbpen(,nodedata.txt,);n=textscan(fid, &

9、#39;%s %u %d %f %f %f %f %f %f);fclose(fid);busnumber=size(n 1 ,1)for i=l:busnumberbus(i).name=n 1 (i);bus(i).type=n2(i);bus(i).no=i;bus(i).base_kv=n3(i);bus(i).pg=n4(i);bus(i).qg=n5(i);bus(i).pl=n6(i);bus(i).ql=n7(i);bus(i).pb=n8(i);bus(i).v=1.0;bus(i).angle=o;endfid=fbpencaclinedata.txt)；a=textsca

10、n(fid, '%s %s %f%f%f%f); fclose(fid);aclinenumber=size(a 1, 1)for i=l:aclinenumberacline(i).fbname=a 1 (i);acline(i).tbname=a 2(i);aciine(i).base_kv=a 3 (i);acline(i).r=a 4(i);acline(i).x=a5(i);acline(i).hb=a6(i);for k=l:busnumberif strcmp(acline(i).fbname, bus(k).name) acline(i).fbno=bus(k).no;

11、endif strcmp(acline(i).tbname, bus(k).name)acline(i).tbno=bus(k).no;endendendfid=fopen(ft ransdata.txt*);t=textscan(fid, *%s %f%f%s %f%f%f%f);fclose(fid);tansnumber=size(t 1, 1)for i=l:tansnumbert rans(i). fbname=t 1 (i);trans(i).fbbase_kv=t 2 (i);t rans(i). fbratedk v=t 3(i);t rans(i). tbname=t 4(i

12、);trans(i).tbbase_kv=t5(i);trans(i).tbrated_kv=t 6 (i);trans(i).r=t7(i);trans(i).x=t 8 (i);for k=l:busnumberif strcrnp(trans(i).fbnarne, bus(k).name)trans(i).fbno=bus(k).no;endif strcmp(trans(i).tbname, bus(k).name)trans(i).tbno=bus(k).no;endendtrans(i).k=trans(i).tbrated_kv*trans(i).fbbase_kv/trans

13、 .fbi*ated_kv/trans(i).tb base_kv;tempx=trans(i).fbrated_kva2/trans(i).fbbase_kva2;trans(i).x=tempx*trans(i).x;t rans(i).r=tempx*trans(i) .r;end%n=0%trans(l)%trans(2)% for y=g+jb matrixg,b,b2=formymatrix(bus,busnumber,acline,aclinenumber,trans,tansnumber); %b:=b'b2:=b”dlmwrite(,gmatrix.txt, g, &

14、#39;delimiter', 't','precision1, 6);dlmwrite(,bmatrix.txt', b, 'delimiter*, t','precision1, 6);%g%b%b2yn=g+j*b%不包括发电机导纳和负荷导纳的节点导纳矩阵pausejp,jq=formjpqmatrix(bus,b,b2,busnumber);jp；ijp=-inv(jp);jq；ijq二 inv(jq);%pause%maxiteration=0for i=l:busnumbernodev(i)=bus(i).v;node

15、a(i)=bus(i).angle;vx(i)=bus(i).v*cos(bus(i).angle);vy(i)=bus(i).v*sin(bus(i).angle);dqgql(i)=bus(i).qg-bus(i).ql; dpgpl(i)=bus(i).pg-bus(i).pl;endnodev=nodev,;nodea=nodeaf;%vx=vx'%vy=vy'dqgql=dqgql：dpgpldpgpl*;%pause;%for nointer=l:10maxdp= 1.;maxdq= 1.;epsilon=0.000001;noiteration=0;while (

16、maxdp>epsilon)&(maxdq>epsilon)deltap,deltaq,maxdp,maxdq=formdpqvector(bus,nodev,nodea,dqgql,dpgpl,b,g,busnumber);deltap;deltaq;maxdp maxdqda=ijp*deltap;dv=ijq*deltaq;nodea=nodea+da;nodev=nodev+dv;noiteration=noiteration+1 pauseif noiteration>20breakendendfor i=l:busnumberbus(i).v=nodev(

17、i);nodev(i)=nodev(i)*bus(i).base_kv;bus(i).angle=nodea(i);n odea(i)=nodea(i) * 180/pi;endnoiteration%潮流计算结束nodea=nodea,nodev=nodev,%潮流计算节点电压模及相角y=1.0/(j*0.3)pausefor i=l:3yn(i,i)=yn(i,i)+y;%将发电机导纳加入节点导纳矩阵endfor i=4:9y=(bus(i).pl-j* bus(i).ql)/bus(i).va2yn(i,i)=yn(i,i)+y;%将负荷导纳加入节点导纳矩阵endyn %包括发电机导纳和

18、负荷导纳的节点导纳矩阵um=0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.01;%以第四个元素为1的单位向量z4=inv(yn)*um%z4:节点阻抗矩阵的第四列if=(bus(4).v*cos(bus(4).angle)+j*bus(4).v*sin(bus(4).angle)/z4(4) % 母线对 地短路电流标幺值ib230=100.0/sqrt(3)/230.0ifymzh=abs(ifl<ib230)%母线对地短路电流有名值1.6208kafor i=l:busnumbervo(i)=bus(i).v*cos(bus(i).angle)-fj*bus(i)

19、.v*sin(bus(i).angle);vf(i)=v0(i)-z4(i)*if;endvovf%由vo和vf计算各支路电流程序 补上for i=l:busnumberfor k=l:busnumberi(i,k)=(vf(k)-vf(i)*yn(i,j); %求出各支路电流endend5.3运行计算结果及分析sbase mva=100 busnumber =9 aclincnumbcr =6 tansnumber =3yn =0 -17.361h+ 17.0648i0 +17.361h-1.3652+11.604h -1.9422+10.5107i-16.0000i0+16.0000i0

20、-17.0648i3.3074-1.36522.5528 -17.3382i-1.1876+ 5.975h-1.94223.2242 -15.8409i+ 17.361h00000-39.3089i0+11.6041i0+ 10.5107i-1.2820 +5.5882i0 +16.0000i0 bo6944i-h7o648i0 +1736 二 i13652+11.604= 一9422+105107i+17.36 二 i'93333i0+160000i0 ,203982i33074,393089i'3652+二6041i39159,178835i11876 + 5.975 =5

21、500002i41700'6 一 562i+ 16.019422+ 105107i'282011876 + 5.975 =2.00047 ,354456i16171 +136980i97oc43i321539i6171 +136980i3oc016 ,236635i155 一0 +17064ooi2820 + 55882i11551 + 97843i2.4371'1876 + 5.975= 0 00 09.7x43i00'.2x20 +32539imaxdp=4.9-61£007 maxdq h7.6235c007 noiteration h9node

22、a h0 96687 4.7711nodevh 一oc一oc2.8047 ,354456i06171 +13.69oo0i 055g?2i 024066 43499 40173227220224229227'6171 +1369oo0i 00 0 2.7722 ,23.3032i' 1551 + 0 +17.064m 0 11551 + 97843i 2.43713.7991 0.6215 1.9256231z4 =0.0433 + 0.1207i0.0308 + 0.0657i0.0296 +0.067 li0.0517+ 0.1439i0.0547 +0.1143i0.05

23、24 + 0.1168i0.0372 + 0.0794i0.0392+ 0.077h0.0354 + 0.0802iif=1.9254 6.1630iib230 =0.2510ifymzh =1.6208v0 =1.0000 0.9858+0.1680i0.9965+0.0832i-0.0683i09940+0.0660i0.9856+0.0107i0.9861 0.0414i0.9549 - 0.0726i1.0028+ 00337i0.9731vf =0.1727+0.0348i 0.5216+0.2314i 0.5262+0.1367i o-o.ooooi 0450+0.0447i0.1

24、523 +0.0294ii 0.4331 +0427i0.4351 +0042i0.4406 + 0.0977i1 0000.6036 2.9983i000000000001.4195-1.4162i00000000000.6657 1,4607i-0.6036 + 2.9983i000-0.7165 + 1.6216i -0.6052 + 1.543 li0000000.7165 1.6216i000.9278+ 1.6048i000000.6052 1.5431i0000-0.7509 +1.5240i0-1.4195 + 1.4162i000.9278 -1.6048i0 00.5242 + 0.0900i0000000-0.5242 0.0900i00.0575+0.0616i00-0.6657 + 1.4607i000.7509 - 1.5240i0-0.0575 0.0616i0六、总结这次的程序实验让我明白了短路潮流计算如何通过mat lab软件进行计算， 也了解到电力系统屮出现短路故障时，系统功