昆明理工大学计算机辅助分析课程设计及matlab自动实现潮流调节与计算.doc

《昆明理工大学计算机辅助分析课程设计及matlab自动实现潮流调节与计算.doc》由会员分享，可在线阅读，更多相关《昆明理工大学计算机辅助分析课程设计及matlab自动实现潮流调节与计算.doc（65页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、重要声明： 1、本文只是提供一定参考，有错勿怪。2.、本文采用的潮流调节及计算只针对本实验，其中有很多固化的数据与参数，以及特殊的公式、请根据要求的不同进行一定的调整（潮流计算的程序适用任何情况、不同的只是系统参数）。3、本文潮流调节为自动调节生成，调节变量有发电机2、3电压，三台变压器分接头，没有进行无功补偿，请根据具体情况加入或删减变量。4、自动调节潮流程序中，控制变量为各点电压首次迭代值、最终各点电压。最终选择标准为网损最低的一组数据。5、本实验报告为草版，有错误，但懒得改。6、鉴于个人所有文档免费得原则，本文为免费文档，请尽量下载。7、最后！都免费了，要求不要太高。要求及原始数据电力系

2、统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。在电力系统中可能发生的各种故障中，危害最大且发生概率较高的首推短路故障。产生短路故障的主要原因是电力设备绝缘损坏。短路故障分为三相短路、

3、两相短路、单相接地短路及两相接地短路。其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路。短路故障大多数发生在架空输电线路。电力系统设计与运行时，要采取适当的措施降低短路故障的发生概率。短路计算可以为设备的选择提供原始数据。一、课程设计的目的与任务:该课程设计为电气工程及其自动化专业的主要专业课设计之一，是学习电力系统分析基础和电力系统计算机辅助分析课程后的一个重要的实践性教学环节。其目的在于巩固和加深对电力系统潮流和短路电流计算基本原理的理解，学习和掌握应用计算机进行电力系统设计和计算的方法，培养学生独立分析和解决问题的能力。二、课程设计的基本要求:掌握电力系统等值模型和参数计算，

4、以及潮流和短路计算的基本原理，学会应用计算机计算系统潮流分布和短路电流的方法。三、课程设计选题原则:该课程设计是根据电力系统分析基础和电力系统计算机辅助分析课程内容，结合实际工程和科研的电力系统网络进行系统的潮流和短路电流计算。四、课程设计内容:1、课程设计地点：机房和教室2、设计内容（1）电力系统潮流计算（2）电力系统短路电流计算具体内容包括：1）系统等值模型和参数计算。2）利用现有的潮流和短路电流计算程序进行计算，或编制电力系统潮流和短路计算程序。3）计算结果分析。4）编写说明书。3、时间安排1、系统等值模型和参数计算（2天）2、读懂现有的算法程序或视需要编制算法程序并进行计算（4天）3、

5、计算结果整理和分析（1天）4、编写说明书（2天）5、成绩考核（1天）五、课程设计主要参考资料：、李光琦.电力系统暂态分析.北京：水利电力出版社.第二版、陈珩.电力系统稳态分析.北京：水利电力出版社.第二版、孟祥萍.基于的电力系统设计与计算.长春：吉林人民出版社、孟祥萍.电力系统分析.北京：高等教育出版社、电力系统潮流和短路计算课程设计指导书.自编六、课程设计成绩考核：根据学生的平时表现，说明书，以及设计完成后的口试成绩等进行综合评定。课程结束时，要求学生有课程设计说明书。课程设计成绩以平时成绩（占50%）、设计说明书（占50%）进行综合评定。原始数据与思考题目标电网单线图系统参数表1. 线路参

6、数表线路编号线路型号线路长度（km）线路电阻/km线路正序电抗/km线路容纳之半S/km4-5LGJ-240/301040.0470.41.784-6LGJ-120/701100.0741.475-7LGJ-120/25980.0791.606-9LGJ-95/551630.0921.807-8LGJ-240/30880.0471.788-9LGJ-240/301220.0471.78说明：线路零序电抗为正序电抗3倍。表2. 变压器参数表线路编号变压器型号变压器变比（kV）短路电压百分数(%)2-7SSPL-24232.5%/20103-9SSPL-24232.5%/1561-4SSPL-24

7、232.5%/17.511说明：变压器零序电抗与正序电抗相等，且均为/Y0接法。表3. 发电机参数表1发电机额定功率MW额定电压kV额定功率因数120016.50.852180180.85310013.80.85 表4. 发电机参数表2发电机母线名（S）（S）（S）110.14600.06088.960.09690.0969047.28220.89580.11986.000.86450.19690.53512.80331.31250.18138.591.25780.25000.6006.02表5. 负荷数据表节点号有功负荷（MW）无功负荷（MVA）51255069030810035思考题1、

8、牛顿拉夫逊法与PQ分解法有哪些联系？有哪些区别？二者的计算性能如何？2、 为什么在用计算机对某网络初次进行潮流计算时往往是要调潮流，而并非任何情况下只一次送入初始值算出结果就行呢？要考虑什么条件？各变量是如何划分的？哪些可调？哪些不可调？3、 发电机节点的注入无功为负值说明了什么？4、 线路有功潮流最有可能的流向是？线路无功潮流最有可能的流向是？5、 负荷功率因数对系统潮流有什么影响？6、 潮流控制的主要手段有哪些？7、 如何降低系统网损？8、 如何校验主设备的短路容量？9、 利用节点阻抗矩阵计算短路电流和利用导纳矩阵计算短路电流的算法有何异同？10、减小短路电流的措施有哪些？设计主要内容：1

9、、完成调试潮流计算与短路计算程序调试；计算程序以电力系统计算机辅助分析所用教材中的例程为基础，要求掌握算法原理，建议进行适当改写以利于本次设计应用。2、根据设计资料准备计算数据；整理原始数据，使其符合计算程序的输入条件。要求比较采用变压器型模型前后两种情况的差异。3、为指定电网安排合理的稳态运行方式；合理设置目标电网中的控制变量，以1号机为平衡机，将潮流分布调整为可行的运行方式（要求所有节点的电压不得低于1.0 p.u.，也不能高于1.05p.u.）绘制潮流分布图。分析该运行方式的特点，提出进一步降低网损的改进的措施，并利用计算结果简单证明之。具体内容参考设计资料所附思考题。4、列表分析节点7

10、发生各种短路故障时的短路电流，以及运行方式对故障电流的影响；5、完成设计说明书。要求结构清晰，即从电网模型建立，到初态潮流分布计算，到问题分析，再到调整后的方案计算。故障电流计算也类似。不必逐一回答所有思考题。参数计算、潮流调节、与潮流计算matlab源程序%参数计算的clc; clear;L=104,110,98,163,88,122;r=0.047,0.074,0.079,0.092,0.047,0.047;b=1.78e-006,1.47e-006,1.60e-006,1.80e-006,1.78e-006,1.78e-006;disp(线路电抗X=);X1=100/2202*0.4*L

11、; %线路各段电抗（均为标幺值，取Ub=220KV,SB=100MVA，下同）disp(X1);for i=1:6 R(i)=100/2202*L(i)*r(i); %加分号避免出现过程项，下同end %线路各段电阻disp(线路电阻R=); disp(R); %输出R for i=1:6Z(i)=R(i)+j*X1(i);Z0(i)=R(i)+3*j*X1(i); enddisp(线路正序、负序阻抗Z1=Z2=); disp(Z); for i=1:6 Bl(i)=2202/100*b(i)*L(i); B(i)=0+j*Bl(i); end %线路各容纳之半disp(线路容纳之半B/2=)

12、; disp(B); %输出Bdisp(线路零序阻抗Z0:);disp(Z0);U0=11,10,6;SN=240,220,120;UN=242,242,242;for i=1:3; XT(i)=j*U0(i)*2422/SN(i)/2202; XT0(i)=XT(i);end;disp(变压器正、负序电抗值Xt：);disp(XT);disp(变压器零序电抗值Xt0：);disp(XT0);x=0.0608,0.1198,0.1813; %发电机暂态电抗有名值UGN=16.5,18,13.8;for i=1:3 K(i)=242/UGN(i);endfor i=1:3 Xdg(i)=j*10

13、0/2202*x(i)*K(i)*K(i);end %发电机暂态电抗disp(发电机暂态电抗Xdg=); disp(Xdg); %输出XdgPL=125,90,100;QL=50,30,35;for i=1:3 ZL(i)=100/(PL(i)2+QL(i)2)*(PL(i)+j*QL(i);end %负荷节点阻抗disp(负荷节点阻抗Zl=); disp(ZL); %输出ZL%支数参数形成的矩阵B1、支数参数形成的矩阵B2、节点号及其对地阻抗形成的矩阵X的形成B1=1,4,XT(1),0,1,0;2,7,XT(2),0,1,0;3,9,XT(3),0,1,0;4,5,R(1)+X1(1)*j

14、,B(1),1,0;4,6,R(2)+ X1(2)*j,B(2),1,0;5,7,R(3)+X1(3)*j,B(3),1,0;6,9,R(4)+X1(4)*j,B(4),1,0;7,8,R(5)+X1(5)*j,B(5),1,0;8,9,R(6)+X1(6)*j,B(6),1,0;disp(支数参数形成的矩阵：B1=); disp(B1); %输出支数参数形成的矩阵B1P=200,180,100;W=(P/100)+(P/0.85)*sqrt(1-0.85*0.85)/100)*j;B2=W(1),0,1.05,1,0,1; W(2),0,1,1,0,3; W(3),0,1,1,0,3;0,0

15、,1,0,0,2;0, 1.25+0.5i,1,0,0,2;0, 0.9+0.3i,1,0,0,2;0,0,1,0,0,2;0,1+0.35i,1,0,0,2;0,0,1,0,0,2;disp(支数参数形成的矩阵：B2=); disp(B2); %输出支数参数形成的矩阵B2%潮流自动调节的与自动计算的n=9;nl=9;isb=1;pr=0.0001;mm=0; X=1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0; 9,0;disp(*);disp(只考虑各点电压水平，得到的满足要求的数据如下：);for aa=0:5 B2(2,4)=1+0.01*aa;for ab=0:5B2

16、(3,4)=1+0.01*ab;for ba=0:5B1(1,5)=0.95+0.025*ba;for bb=0:5B1(2,5)=0.95+0.025*bb;for bc=0:5B1(3,5)=0.95+0.025*bc;Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);S1=zeros(nl);for i=1:n if X(i,2) =0; p=X(i,1); Y(p,p)=1./X(i,2); end endfor i=1:nl if B1(i,6)=0 p=B1(i,1); q=B1(i,2); else p=

17、B1(i,2); q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2; Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;endG=real(Y);B=imag(Y); for i=1:n e(i)=real(B2(i,3); f(i)=imag(B2(i,3); V(i)=B2(i,4); endfor i=1:n S(i)=B2(i,1)-B2(i,2); B(i,i)=B(i,i)+B2(i,5)

18、;endP=real(S);Q=imag(S);ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;while IT2 =0 IT2=0;a=a+1; for i=1:n if i=isb C(i)=0; D(i)=0; for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1); D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1); end P1=C(i)*e(i)+f(i)*D(i); Q1=f(i)*C(i)-D(i)*e(i); V2=e(i)2+f(i)2; if B2(i,6)=3 DP=P(i)-P1; DQ

19、=Q(i)-Q1; for j1=1:n if j1=isb & j1=i X1=-G(i,j1)*e(i)-B(i,j1)*f(i); X2=B(i,j1)*e(i)-G(i,j1)*f(i); X3=X2; X4=-X1; p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2; elseif j1=i & j1=isb X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i); X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i); X3=D(i

20、)+B(i,i)*e(i)-G(i,i)*f(i); X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2; end end else DP=P(i)-P1; DV=V(i)2-V2; for j1=1:n if j1=isb&j1=i X1=-G(i,j1)*e(i)-B(i,j1)*f(i); X2=B(i,j1)*e(i)-G(i,j1)*f(i); X5=0; X6=0; p=2*i-1;q=2

21、*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2; elseif j1=i&j1=isb X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i); X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i); X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2; end end end endend

22、for k=3:N0 k1=k+1;N1=N; for k2=k1:N1 J(k,k2)=J(k,k2)./J(k,k); end J(k,k)=1; if k=3 k4=k-1; for k3=3:k4 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; end if k=N0,break;end for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; endelse for k3=k1:N0 for k2=k1:N1 J(

23、k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; end endendfor k=3:2:N0-1 L=(k+1)./2; e(L)=e(L)-J(k,N); k1=k+1; f(L)=f(L)-J(k1,N);endfor k=3:N0 DET=abs(J(k,N); if DET=pr IT2=IT2+1; endendICT2(a)=IT2;ICT1=ICT1+1;for k=1:n dy(k)=sqrt(e(k)2+f(k)2);end for i=1:n Dy(ICT1,i)=dy(i); end for vb=1:n if Dy(1,v

24、b)=1.00; zy(vb)=1; else zy(vb)=0; endendend %用高斯消去法解w=-J*Vfor k=1:n V(k)=sqrt(e(k)2+f(k)2); O(k)=atan(f(k)./e(k)*180./pi;endE=e+f*j;for va=1:n if V(va)=1.000; yj(va)=1; else yj(va)=0; endendif yj(1)*yj(2)*yj(3)*yj(4)*yj(5)*yj(6)*yj(7)*yj(8)*yj(9)=1 if zy(1)*zy(2)*zy(3)*zy(4)*zy(5)*zy(6)*zy(7)*zy(8)*

25、zy(9)=1;disp(发电机2出口电压：);disp(B2(2,4);disp(发电机3出口电压：);disp(B2(3,4);disp(变压器1接头电压：);disp(B1(1,5);disp(变压器2接头电压：);disp(B1(2,5);disp(变压器3接头电压：);disp(B1(3,5);disp(各节点的电压大小V为(节点号从小到大排列)：);disp(V);for p=1:n C(p)=0; for q=1:n C(p)=C(p)+conj(Y(p,q)*conj(E(q); end S(p)=E(p)*C(p);endfor i=1:nl if B1(i,6)=0 p=B

26、1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); endSi(p,q)=E(p)*(conj(E(p)*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5)-conj(E(q)*conj(1./(B1(i,3)*B1(i,5);end for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); endSj(q,p)=E(q)*(conj(E(q)*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5)-conj(E(p)*conj(1./(

27、B1(i,3)*B1(i,5); endZDS=0;for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end DS(i)=Si(p,q)+Sj(q,p); ZDS=DS(i)+ZDS;endmm=mm+1;ZDSS(mm)=sqrt(real(ZDS)2+imag(ZDS)2);disp(*);for i=1:nZB224(mm)=B2(2,4);ZB234(mm)=B2(3,4);ZB115(mm)=B1(1,5);ZB125(mm)=B1(2,5);ZB135(mm)=B1(3,5);endende

28、ndendendendendenddisp(*);disp(满足各点电压要求的数据组共有：);disp(mm);LOL=min(ZDSS);for dota=1:mm if ZDSS(dota)=LOL mm=dota; endenddisp(当加入线路损耗时，得到其中线路损耗最小的一组组号及关键数据为：);disp(mm);disp(发电机2出口电压：);disp(ZB224(mm);disp(发电机3出口电压：);disp(ZB234(mm);disp(变压器1接头电压：);disp(ZB115(mm);disp(变压器2接头电压：);disp(ZB125(mm);disp(变压器3接头电

29、压：);disp(ZB115(mm);B2(2,4)=ZB224(mm);B2(3,4)=ZB234(mm);B1(1,5)=ZB115(mm);B1(2,5)=ZB125(mm);B1(3,5)=ZB135(mm);disp(用此组数据算的的潮流结果为：);disp(节点数：);disp(n);disp(节点支路数：);disp(nl);disp(平衡节点：);disp(isb);disp(精度：);disp(pr);disp(B1=); disp(B1); disp(B2=); disp(B2);disp(X=); disp(X);Y=zeros(n);e=zeros(1,n);f=zer

30、os(1,n);V=zeros(1,n);O=zeros(1,n);S1=zeros(nl);disp(节点数：);disp(n);disp(节点支路数：);disp(nl);disp(平衡节点：);disp(isb);disp(精度：);disp(pr);mm=0; X=1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0; 9,0;for i=1:n if X(i,2) =0; p=X(i,1); Y(p,p)=1./X(i,2); end endfor i=1:nl if B1(i,6)=0 p=B1(i,1); q=B1(i,2); else p=B1(i,2); q=B1(

31、i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2; Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;endG=real(Y);B=imag(Y); for i=1:n e(i)=real(B2(i,3); f(i)=imag(B2(i,3); V(i)=B2(i,4); endfor i=1:n S(i)=B2(i,1)-B2(i,2); B(i,i)=B(i,i)+B2(i,5);endP=real(S);

32、Q=imag(S);ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;while IT2 =0 IT2=0;a=a+1; for i=1:n if i=isb C(i)=0; D(i)=0; for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1); D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1); end P1=C(i)*e(i)+f(i)*D(i); Q1=f(i)*C(i)-D(i)*e(i); V2=e(i)2+f(i)2; if B2(i,6)=3 DP=P(i)-P1; DQ=Q(i)-Q1; for

33、j1=1:n if j1=isb & j1=i X1=-G(i,j1)*e(i)-B(i,j1)*f(i); X2=B(i,j1)*e(i)-G(i,j1)*f(i); X3=X2; X4=-X1; p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2; elseif j1=i & j1=isb X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i); X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i); X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2; end end else