复杂网络牛顿—拉夫逊法潮流分析.doc

《复杂网络牛顿—拉夫逊法潮流分析.doc》由会员分享，可在线阅读，更多相关《复杂网络牛顿—拉夫逊法潮流分析.doc（27页珍藏版）》请在三一办公上搜索。

1、山东交通学院电力系统分析课程设计报告书院(部)别 信息科学与电气工程学院 班 级 学 号 姓 名 指导教师 时 间 2014.06.9-2013.06.13 课 程 设 计 任 务 书题 目 复杂网络牛顿拉夫逊法潮流分析 院 (部) 信息科学与电气工程学院 专 业 电气工程及其自动化 班 级 学生姓名 学 号 6 月 9 日至 6 月 13 日 共 1 周指导教师(签字) 负责人(签字) 年 月 日一、设计内容及要求复杂网络牛顿拉夫逊法潮流分析与计算的设计电力系统潮流计算是电力系统中一项最基本的计算，设计内容为复杂网络潮流计算的计算机算法牛顿-拉夫逊法。首先，根据给定的电力系统简图，通过手算一

2、次迭代过程，从而加深对牛顿-拉夫逊法的理解，有助于计算机编程的应用。其次，利用计算机编程对电力系统稳态运行的各参数进行解析和计算；编程完成复杂网络的节点导纳矩阵的形成；电力系统支路改变、节点增减的程序变化；编程完成各元件的功率损耗、各段网络的电压损耗、各点电压、功率大小和方向的计算。二、设计原始资料给出一个46节点、环网、两电源和多引出的电力系统；参数给定，可以选用直角坐标表示的牛拉公式计算，也可以选用极坐标表示的牛拉公式计算。具体题目详见附录题单三、设计完成后提交的文件和图表1计算说明书部分设计报告和手算潮流的步骤及结果2图纸部分：电气接线图及等值电路；潮流计算的计算机算法，即程序；运算结果

3、等以图片的形式附在设计报告中。四、进程安排第一天上午：选题，查资料，制定设计方案；第一天下午第三天下午：手算完成潮流计算的要求；第四天上午第五天上午：编程完成潮流计算，并对照手算结果，分析误差第五天下午：答辩，交设计报告。五、主要参考资料电力系统分析（第三版） 于永源主编，中国电力出版社，2007年电力系统分析，何仰赞 温增银编著，华中科技大学出版社，2002年版；电力系统分析，韩桢祥主编，浙江大学出版社，2001年版；电力系统稳态分析，陈珩 编，水利电力出版社；成绩评定表指导教师成绩答辩小组成绩总评成绩目录摘要.21.1 潮流计算.31.1.1潮流计算要求.41.2 MATLAB简介.52.

4、1 实验设计要求.52.2 电力系统等值网络.5 2.2.1牛顿拉夫逊法程序流程图.52.2.2 实验设计步骤.72.2.3程序输出结果.10心得体会.13参考文献.14附录.15. 摘要 本课程设计通过对地区电网的设计，巩固和运用前面所学到的基础理论知识，掌握电力系统设计的一般原则和方法，培养分析问题和解决问题的能力。在此过程中，要考虑到各方面的相互关系和相互影响，综合地运用课程中所学到的知识，进行独立思考。 潮流计算，电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。潮流计算是电力系统分析中的一种最基本的计算，它的任务是对给

5、定的运行条件确定系统的运行状态，如各母线上的电压（幅值及其相角）、网络中的功率分布及功率损耗等。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。传运用 MATLAB软件进行仿真潮流计算，图形界面直观，运行稳定，计算准确，提高了计算速度，各个类的有效封装又使程序具有很好的模块性.可维护性和可重用性。因此，潮流计算在电力系统的规划计算、生产运行、调度管理及科学计算中都有着广泛的应用。也就是说，对于电气工程及其自动化专业的学生来说，掌握潮流计算是非常重要和必要的。 关键字：

6、 潮流计算 牛拉法 Matlab 1潮流计算1.1潮流计算概述 潮流计算是研究电力系统稳态运行情况的一种基本电气计算，常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。具体表现在以下方面：(1)在电网

7、规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。总结为在电力系统运行方式和规

8、划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和基础的计算。1.1.1潮流计算的要求电力系统运行必须满足一定技术和经济上的要求。这些要求够成了潮流问题中某些变量的约束条件，常用的约束条件如下：1.节点电压应满足 从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在

9、额定电压附近。PU节点电压幅值必须按上述条件给定。因此，这一约束条件对PQ节点而言。2.节点的有功功率和无功功率应满足 PQ节点的有功功率和无功功率，以及PU节点的有功功率，在给定是就必须满足上述条件，因此，对平衡节点的P和Q以及PU节点的Q应按上述条件进行检验。3.节点之间电压的相位差应满足 为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。这一约束的主要意义就在于此。 因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束条件。常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用约束条件进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至

10、修改系统的运行方式，重新进行计算。1 MATLAB简介1.2 MATLAB概述 MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 基本功能MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究必须进行有效数值计算的众多科学领域

11、提供了一种全面的解决方案。2.1实验设计要求 下图示网络中，变压器的变比、各支路阻抗和一半的对地电纳均以标么值标于图中，设5节点为平衡节点，电压为1，节点4为PV节点，电压为1,P=0.5。试求该网络的潮流分布，方法不限，求解精度为10e-5。2.2电力系统等值网络2.2.1牛顿拉夫逊法程序流程图2.2.2实验设计步骤第一步，节点编号。由已知条件可知G1为平衡节点，编号为5；G2为PV节点，编号为4；其余为PQ节点，编号为1，2，3；列出支路表与节点数据表如表1-1和1-2所示。表1-1 支路数据ijRXB/2(或k)1 1245233230.0250.030.02000.080.10.060

12、.19050.19050.070.090.051.051.05表1-2 节点数据iUPGQGPDQD类别12345待求待求待求1.01.00000.5待定000待定待定0.800.180000.530.12000PQPQPQPV平衡第二步，求节点导纳矩阵Y11=j0.09+1/(0.03+j0.1)+j0.07+1/(0.025+j0.08)=6.311-j20.4652Y12=Y21=-1/(0.025+j0.08)=-3.5587+j11.3879Y13=Y31=-1/(0.03+j0.1)=-2.7523+j9.1743Y22=1/(1.05*1.05*j0.1905)+j0.05+j0

13、.07+1/(0.025+j0.08)+1/(0.02+j0.06)=8.5587-j31.0292Y23=Y32=-1/(0.02+j0.06)=-5+j15Y24=Y42=-1/(1.05*j0.1905)=j4.9993Y33=1/(1.05*1.05*j0.1905)+j0.09+1/(0.03+j0.1)+j0.05+1/(0.02+j0.06)=7.7523-j28.7956Y35=Y53=-1/(1.05*j0.1905)=j4.9993Y44=1/(j0.1905)=-j5.2493Y55=1/(j0.1905)=-j5.2493Y14=Y41=Y15=Y51=Y25=Y52=

14、Y34=Y43=Y45=Y54=0则节点导纳矩阵为第三步，设定初值。给定个节点电压初始值第四步，根据节点功率方程 (1)式，求出修正方程式常数项向量 P1=-0.8055 Q1=-0.372;P2=-0.18 Q2=0.2475; P3=0 Q3=-0.3875 P4=0.5 U4=-2第五步：将电压初始值代入下式，求雅克比矩阵方程。(2)第六步：解修正方程式，求出修正量第七步：修正各节点电压第八步：将再代入（1）式，求出第九步：校验是否收敛，其收敛条件为式中，为向量中大分量的绝对值。这个收敛条件比较直观，它可以直接显示出最终结果的功率误差第十步：如果收敛，就进一步计算各段电力线路潮流和平衡节

15、点功率，并打印出计算结果；如果不收敛，转回3步进行下一次迭代计算，直到收敛为止。2.2.3程序输出结果注意：题目中节点编号改动，原来的5编为1，原来的4编为5，原来的1编为4B1 = Columns 1 through 5 4.0000 2.0000 0.0250 + 0.0800i 0 + 0.1400i 1.0000 4.0000 3.0000 0.0300 + 0.1000i 0 + 0.1800i 1.0000 2.0000 3.0000 0.0200 + 0.0600i 0 + 0.1000i 1.0000 1.0000 3.0000 0 + 0.1905i 0 1.0500 5.0

16、000 2.0000 0 + 0.1905i 0 1.0500 Column 6 0 0 0 0 0 B2 =Columns 1 through 5 0 0 1.0000 1.0000 0 0 0.1800 + 0.1200i 1.0000 0 0 0 0 1.0000 0 0 0 0.8000 + 0.5300i 1.0000 0 0 0.5000 0 1.0000 1.0000 0 Column 6 1.0000 2.0000 2.0000 2.0000 3.0000 导纳矩阵 Y=0 - 5.2493i 0 0 + 4.9994i 0 0 0 8.5587 -31.0292i -5.00

17、00 +15.0000i -3.5587 +11.3879i 0 + 4.9994i0 + 4.9994i -5.0000 +15.0000i 7.7523 -28.7956i -2.7523 + 9.1743i 0 0 -3.5587 +11.3879i -2.7523 + 9.1743i 6.3110 -20.4022i 0 0 0 + 4.9994i 0 0 0 - 5.2493i迭代次数： 4没有达到精度要求的个数： 7 7 3 0各节点的实际电压标幺值E为(节点号从小到大排列)： 1.0000 1.0104 - 0.1023i 1.0161 - 0.0982i 0.9814 - 0.

18、1275i 1.0000 - 0.0023i-各节点的电压大小V为(节点号从小到大排列)： 1.0000 1.0155 1.0208 0.9896 1.0000-各节点的电压相角sida为(节点号从小到大排列)： 0 -5.7840 -5.5225 -7.4010 -0.1321各节点的功率S为(节点号从小到大排列)：0.4911 + 0.1694i -0.1800 - 0.1200i 0.0000 + 0.0000i -0.8000 - 0.5300i 0.5000 + 0.1970i-各条支路的首端功率Si为(顺序同您输入B1时一致)： -0.4127 - 0.2547iS(4,2)=-0

19、.41269-0.25474i- -0.3873 - 0.2753iS(4,3)=-0.38731-0.27526i- -0.0979 - 0.1087iS(2,3)=-0.097919-0.10873i- 0.4911 + 0.1694iS(1,3)=0.49115+0.16945i- 0.5000 + 0.1970iS(5,2)=0.5+0.19702i-各条支路的末端功率Sj为(顺序同您输入B1时一致)： 0.4179 + 0.1307iS(2,4)=0.41792+0.13073i- 0.3930 + 0.1122iS(3,4)=0.39298+0.11221i- 0.0982 + 0

20、.0058iS(3,2)=0.098168+0.0058104i- -0.4911 - 0.1180iS(3,1)=-0.49115-0.11802i- -0.5000 - 0.1420iS(2,5)=-0.5-0.142i-各条支路的功率损耗DS为(顺序同您输入B1时一致)： 0.0052 - 0.1240iDS(4,2)=0.0052321-0.12401i- 0.0057 - 0.1630iDS(4,3)=0.0056674-0.16305i- 0.0002 - 0.1029iDS(2,3)=0.00024932-0.10292i- 0 + 0.0514iDS(1,3)=0+0.0514

21、24i- -0.0000 + 0.0550iDS(5,2)=-1.1102e-016+0.05502i-结果分析: 经分析，4节点既是无功分点，又是有功分点。总的功率损耗为0.01114882-0.283536i。 心得体会 在本次电力系统分析课程设计制作过程中，我投入了极大的热情和精力。通过对地区电网的设计，巩固和运用前面所学到的基础理论知识，掌握电力系统设计的一般原则和方法，培养了分析问题和解决问题的能力。通过对题目“地区电力网设计”的设计，对课上所学的知识有了更深刻的理解，学会了如何实际进行电力网络的潮流计算，既能考虑到细微之处又能整体统筹规划。在设计过程中，我也了解到了我国的各地区电网

22、的现状，了解了与先进技术的差异，这更激励我认真学习专业知识，为我国的电力事业尽一分力量。本次课程设计运用所学的基础知识，综合考虑了各个方面的相互关系，在设计的过程中培养了我独立思考的能力，使所学的理论知识和实践相结合，对电力网的基本知识有了初步的了解，为以后的工作打下了一定的基础。参考文献1.陈珩.电力系统稳态分析M,中国电力出版社，2007，第三版2.韩祯祥.电力系统分析M，浙江大学出版社，2005，第三版3.祝书萍.电力系统分析课程实际设计与综合实验M，中国电力出版社，2007，第一版4.电力系统稳态分析，百度文库5.电力系统分析（第三版） 于永源主编，中国电力出版社，2007年电力系统分

23、析，何仰赞 温增银编著，华中科技大学出版社，2002年版；电力系统分析，韩桢祥主编，浙江大学出版社，2001年版；电力系统稳态分析，陈珩 编，水利电力出版社；附录注意：题目中节点编号改动，原来的5编为1，原来的4编为5，原来的1编为4%本程序的功能是用牛顿拉夫逊法进行潮流计算% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳% 5、支路的变比；6、支路首端处于K侧为1，1侧为0不含变压器% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量% 6、节点分类标号：1为平衡节点（应为1号节点）

24、；2为PQ节点；% 3为PV节点；clear;n=5;%input(请输入节点数：n=);nl=5;%input(请输入支路数：nl=);isb=1;%input(请输入平衡母线节点号：isb=);pr=0.0001;%input(请输入误差精度：pr=);B1=4 2 0.025+0.08i 0.14i 1 0; 4 3 0.03+0.1i 0.18i 1 0; 2 3 0.02+0.06i 0.10i 1 0; 1 3 0.1905i 0 1.05 0; 5 2 0.1905i 0 1.05 0;B2=0 0 1 1 0 1; 0 0.18+0.12i 1 0 0 2; 0 0 1 0 0

25、 2; 0 0.8+0.53i 1 0 0 2; 0.5 0 1 1 0 3; %input(请输入各节点参数形成的矩阵： B2=);Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);% % %-for i=1:nl%支路数 if B1(i,6)=0%左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=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

26、); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2;%对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;%对角元1侧 end%求导纳矩阵disp(导纳矩阵 Y=);disp(Y)%-G=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);%PV节点电压给定模值 endfor i=1:n%给定各节点注入功率 S(i)=B2(i,1)-B2(

27、i,2); %i节点注入功率SG-SL B(i,i)=B(i,i)+B2(i,5);%i节点无功补偿量 end%=P=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率ICT1=0;IT2=1;N0=2*(n);N=N0+1;a=0; %迭代次数ICT1、a；不满足收敛要求的节点数IT2while IT2=0 % N0=2*n 雅可比矩阵的阶数；N=N0+1扩展列 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); %(

28、Gij*ej-Bij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1); %(Gij*fj+Bij*ej) end P1=C(i)*e(i)+f(i)*D(i); %节点功率P计算ei(Gij*ej-Bij*fj)+fi(Gij*fj+Bij*ej) Q1=C(i)*f(i)-e(i)*D(i); %节点功率Q计算fi(Gij*ej-Bij*fj)-ei(Gij*fj+Bij*ej)%求i节点有功和无功功率P,Q的计算值 V2=e(i)2+f(i)2;%电压模平方%= 以下针对非PV节点来求取功率差及Jacobi矩阵元素 = if B2(i,6)=3%非PV节

29、点 DP=P(i)-P1;%节点有功功率差 DQ=Q(i)-Q1; %节点无功功率差 %= 以上为除平衡节点外其它节点的功率计算 =%= 求取Jacobi矩阵 = for j1=1:n if j1=isb&j1=i%非平衡节点&非对角元 X1=-G(i,j1)*e(i)-B(i,j1)*f(i);% dP/de=-dQ/df X2=B(i,j1)*e(i)-G(i,j1)*f(i);% dP/df=dQ/de X3=X2; % X2=dp/df X3=dQ/de X4=-X1; % X1=dP/de X4=dQ/df p=2*i-1;q=2*j1-1; J(p,q)=X3;J(p,N)=DQ;m=p+1; % X3=dQ/de J(p,N)=DQ节点无功功率差 J(m,q)=