1、课程设计(论文)题 目 名 称 电力系统潮流计算 课 程 名 称 电力系统稳态分析 学 生 姓 名 学 号 系 、 专 业 电气工程系电气工程及其自动化电力方向 指 导 教 师 2009 年 1 月 6 日前 言在如今的社会，电力已经成为人们必不可少的需求，而建立结构合理的大型电邵阳学院课程设计2力系统不仅便于电能生产与消费的集中管理、统一调度和分配，减少总装机容量，节省动力设施投资，且有利于地区能源资源的合理开发利用，更大限度地满足地区国民经济日益增长的用电需要。电力系统建设往往是国家及地区国民经济发展规划的重要组成部分。电力系统的出现，使高效、无污染、使用方便、易于调控的电能得到广泛应用，
2、推动了社会生产各个领域的变化，开创了电力时代，发生了第二次技术革命。电力系统的规模和技术水准已成为一个国家经济发展水平的标志之一。电力系统稳态分析包括潮流计算（或潮流分析)和静态安全分析。潮流计算针对电力革统各正常运行方式，而静态安全分析则要研究各种运行方式下个别系统元件退出运行后系统的状况。其目的是校验系统是否能安全运行，即是否有过负荷的元件或电压过低的母线等。原则上讲，静态安全分析也可 U 用潮流计算来代替。但是一般静态安全分析需要校验的状态数非常多，用严格的潮流计算来分析这些状态往往计算量过大，因此不得不寻求一些特殊的算法以满足要求。牛顿法是数学中解决非线性方程式的典型方法，有较好的收敛
3、性。解决电力系统潮流计算问题是以导纳距阵为基础的，因此，只要在迭代过程中尽可能保持方程式系数距阵的稀疏性，就可以大大提高牛顿法潮流程序的放率。自从 20 世纪 60 年代中期利用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、速度方面都超过了阻抗法，成为直到目前仍在广泛采用的优秀方法。邵阳学院课程设计3目 录第一章 系统概述 .41.1 设计目的与要求 41.1.1 设计目的 .41.1.2 设计要求 .41.2 设计题目 41.3 设计内容 4第二章 潮流计算设计题目 .52.1 潮流计算题目 52.2 对课题的分析及求解思路 6第三章 潮流计算算法及手工计算 .63.1 潮流计算算法 63
4、.2 关于电力系统潮流计算手工计算 83.2.1 节点导纳矩阵 .83.2.2 简化雅可比矩阵 .93.2.3 修正、迭代 10第四章 Matlab 概述 104.1 Matlab 简介 .104.2 矩阵的运算 114.2.1 四则运算 .124.2.2 与常数的运算 .124.2.3 基本数学运算 .124.2.4 逻辑关系运算 .12第五章 潮流计算流程图及源程序 .135.1 潮流计算流程图 135.2 潮流计算源程序 145.3 运行计算结果 19总 结 .20参考文献 .21邵阳学院课程设计4第一章 系统概述1.1 设计目的与要求1.1.1 设计目的1. 掌握电力系统潮流计算的基本
5、原理；2. 掌握并能熟练运用一门计算机语言（MATLAB语言或FORTRAN或C语言或C+语言）；3. 采用计算机语言对潮流计算进行计算机编程计算。 1.1.2 设计要求1. 程序源代码； 2. 给定题目的输入，输出文件； 3. 程序说明； 4. 给定系统的程序计算过程； 5. 给定系统的手算过程（至少迭代2次）。 1.2 设计题目电力系统潮流计算（牛顿-拉夫逊法、P-Q 分解法）1.3 设计内容1.根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵； 2.赋予各节点电压变量（直角坐标系形式）初值后，求解不平衡量；3.形成雅可比矩阵； 4.求解修正量后，重新修改初值，从2开始重新循环计算；
6、 5.求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和平衡节点功率; 6.上机编程调试；连调； 7.计算分析给定系统潮流分析并与手工计算结果作比较分析。 邵阳学院课程设计58.准备计算机演示答辩，书写该课程设计说明书（必须计算机打印） 。邵阳学院课程设计6第二章 潮流计算设计题目2.1 潮流计算题目图 2-1 电力系统接线图邵阳学院课程设计72.2 对课题的分析及求解思路此电力系统是一个 5 节点，4 支路的电力网络。其中包含 3 个 PQ 节点，一个 PV节点，和一个平衡节点。综合比较牛顿拉夫逊法（直角坐标、极坐标） 、PQ 分解法等多种求解方法的特点，最后确定采用牛顿拉夫
7、逊法（极坐标） 。因为此方法所需解的方程组最少。第三章 潮流计算算法及手工计算3.1 潮流计算算法本题采用了题目要求的牛顿拉夫逊潮流计算的方法。牛顿-拉夫逊法潮流计算的公式。把牛顿法用于潮流计算，采用直角坐标形式表示的如式（1-3）所示的形式。其中电压和支路导纳可表示为：（1-2）iijijjjjjijijijUefYGBf将上述表示式（1-2）代入（1-1）式的右端，展开并分出实部和虚部，便得：（1-3）11()()nniijijiijijjiiijijiijijj jPeGBfGfBeQffef按照以上的分类，PQ 节点的输出有功功率和无功功率是给定的，则第 i 节点的给定功率设为 和 （
8、称为注入功率） 。isPi假定系统中的第 1、2、m 节点为 PQ 节点，对其中每一个节点的 N-R 法表达式F(x)=0如 、 、 形式有些下列方程：0iSi 0iQ（1-4）11()()0nniisiisiijijiijijj jiisiisiijijiijijj jPeGBfGfBeQffef邵阳学院课程设计8=（1、2、m）iPV 节点的有功功率和节点电压幅值是给定的。假定系统中的第m+1、m+2、n-1 节点为 PV 节点，则对其中每一 PV 节点可以列写方程：（1-5）11222()()0nniisiisiijijiijijj jisiisiiPeGBfGfBeU=（m+1、m+2
9、、n-1）i(6)形成雅可比矩阵。N-R 法的思想是 ；本例 ；()0Fxx()PjQFx对 F(x)求偏导的式（1-6） 、式（1-7） ，即式（1-4） 、式（1-5）中的 、0i、 是多维变量的函数，对多维变量求偏导（ 、 、 、 、0iQUieijifijf、 、 、） ，并以矩阵的形式表达称为雅可比矩阵。ieijiPe当 j=i 时，对角元素为（1-6）112()()niijijiiijiijijiiijniijijiiijiijijiiijiiiiGeBfefNePfGfHfQfefLeGBffJfUeff 当 时,矩阵非对角元素为：ji邵阳学院课程设计9（1-7）22()0iii
10、jijijijjiiijijijijjjiijjPQGeBfNJeffHLfUef 由上式不难看出，雅可比矩阵有以下特点。 雅可比矩阵中的诸元素都是节点电压的函数，因此在迭代过程中，它们将随着节点电压的变化而不断的变化。 雅可比矩阵具有结构对称性，数据不对称。如非对角 ，ijjiH， 。ijijijHBeGfjiijijHBeGf 由式（1-7）可以看出，当导纳矩阵中非对角元素 为零时， 。雅可比矩阵ijY中相应的元素也为零，即矩阵是非常稀疏的。因此，修正方程的求解同样可以应用稀疏矩阵的求解技巧。正是由于这一点才使 N-R 法获得广泛的应用。3.2 关于电力系统潮流计算手工计算3.2.1 节点
11、导纳矩阵求得节点导纳矩阵 Y=iijyi0ikik各节点的导纳值如下：; ；10.8342.50Yj21.675.0Yj； ;267 938; ;13j 23j； ;450 40； ;1.1.Yj251.3.750Yj; ;31.675.0j 41j邵阳学院课程设计10; ;3210.3.0Yj420Yj; ;978531.53.70; ;34j 42j; ;50Y 50Y;51.15.0j;237Y;530j;4.56.218750Yj3.2.2 简化雅可比矩阵形成有功迭代和无功迭代的简化雅可比矩阵 B/和 B/B/= B/=将 B/ 和 B/进行三角分解：-0.234654 -0.5640
12、73 -0.443749-0.423068 2.357080-1.02902 -0.234654 -0.564073 -0.423068 邵阳学院课程设计113.2.3 修正、迭代给定 PQ 节点初值和各节点电压相角初值V1=1.050 。 ，V2 (0)=V3(0)=1.0，V4=1.12 (0)=3 (0)=0， 4 (0)=01 作第一次有功迭代，按公式计算节点有功功率不平衡量P2 (0)=-0.55-（-0.024037）=-0.525963P3 (0)=-0.30-（-0.022695）=-0.277305P4 (0)=0.500000P1 (0)/V1(0)=0.454545 P2
13、 (0)/ V2(0)=-0.525963P3 (0)/V3(0)=-0.2773092 做第一次无功迭代，按公式计算无功功率不平衡量，计算时电压相角最新的修正值。Q2 (0)=-0.13-(-0.001550)=-0.039594Q3 (0)=-0.18-(-0.14406)=-0.039588Q2 (0)/ V2(0)=-0.131553Q3 (0)/V3(0)=-0.039588解修正方程式，可得各节点电压幅值的修正量为V3 (0)） =-0.0148550524.)(2于是有：V2(1) = V2(0)+V2 (1)=0.964776V3(1) = V3(0)+V3 (1)=0.985
14、145到这里为止，第一轮有功迭代和无功迭代便做完了。3 按公式计算平衡节点功率，得：P1+jQ1=0.367885+j0.264696经过四轮迭代，节点不平衡功率也下降到 10-5以下，迭代到此结束。第四章 Matlab 概述4.1 Matlab 简介目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1）计算速度快（2）内存需要少（3）计算结果有良好的可靠性和可信性（4）适应性好，亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强（5）简单。MATLAB 是一种交互式、面向对象的程序设计语言，广泛
15、应用于工业界与学术邵阳学院课程设计12界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。MATLAB 程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB 与 C 语言和 FORTRAN 语言相比更容易被掌握。通过 M 语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。另外，MATLAB 提供了一种特殊的工具：工具箱（ TOOLBOXES）.这些工具箱主要包括：信号处
16、理（SIGNAL PROCESSING） 、控制系统（CONTROL SYSTEMS） 、神经网络（NEURAL NETWORKS） 、模糊逻辑(FUZZY LOGIC)、小波(WAVELETS)和模拟（SIMULATION）等等。不同领域、不同层次的用户通过相应工具的学习和应用，可以方便地进行计算、分析及设计工作。MATLAB 设计中，原始数据的填写格式是很关键的一个环节，它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计，主要应从使用的角度出发，原则是简单明了，便于修改。4.2.2 与常数的运算 常数与矩阵的运算即是同该矩阵的每一元素进行运算。但需注意进行数除时，常数通常只
17、能做除数。基本函数运算中，矩阵的函数运算是矩阵运算中最实用的部分，常用的主要有以下几个：det(a) 求矩阵 a 的行列式eig(a) 求矩阵 a 的特征值inv(a)或 a (-1) 求矩阵 a 的逆矩阵rank(a) 求矩阵 a 的秩trace(a) 求矩阵 a 的迹（对角线元素之和）我们在进行工程计算时常常遇到矩阵对应元素之间的运算。这种运算不同于前邵阳学院课程设计13面讲的数学运算，为有所区别，我们称之为数组运算。4.2.3 基本数学运算数组的加、减与矩阵的加、减运算完全相同。而乘除法运算有相当大的区别，数组的乘除法是指两同维数组对应元素之间的乘除法，它们的运算符为“.*”和“./”或
18、“.” 。前面讲过常数与矩阵的除法运算中常数只能做除数。在数组运算中有了“对应关系”的规定，数组与常数之间的除法运算没有任何限制。另外，矩阵的数组运算中还有幂运算（运算符为 . ） 、指数运算（exp） 、对数运算（log） 、和开方运算（sqrt）等。有了“对应元素”的规定，数组的运算实质上就是针对数组内部的每个元素进行的。矩阵的幂运算与数组的幂运算有很大的区别。4.2.4 逻辑关系运算 逻辑运算是 MATLAB 中数组运算所特有的一种运算形式，也是几乎所有的高级语言普遍适用的一种运算。邵阳学院课程设计14第五章 潮流计算流程图及源程序5.1 潮流计算流程图启动输入原始数据形成导纳矩阵 Y给
19、定电压初值 e 0 ， f 0置 K 0利用式 ( 4 - 4 6 ) ( 4 - 4 7 ) 计算d e l t P ( k ) , d e l t Q ( k ) , d e l t U ( k )I d e l t P ( k ) , d e l t Q ( k ) I 0.00001delt(5)=0;u(5)=1.06;for m=1:N1for n=1:N1+1pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n);qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(
20、n)-B(m,n)*cos(delt(m)-delt(n);endpp(m)=p(m)-sum(pt);qq(m)=q(m)-sum(qt);endfor m=1:N1for n=1:N1+1h0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n);n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n);j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(
21、m)-delt(n);L0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n);endH(m,m)=sum(h0)-u(m)2*(G(m,m)*sin(delt(m)-delt(m)-B(m,m)*cos(delt(m)-delt(m);N(m,m)=sum(n0)-2*u(m)2*G(m,m)+u(m)2*(G(m,m)*cos(delt(m)-delt(m)+B(m,m)*sin(delt(m)-delt(m);邵阳学院课程设计18J(m,m)=sum(j0)+u(m)2*(G(m,m)*cos(delt
22、(m)-delt(m)+B(m,m)*sin(delt(m)-delt(m);L(m,m)=sum(L0)+2*u(m)2*B(m,m)+u(m)2*(G(m,m)*sin(delt(m)-delt(m)-B(m,m)*cos(delt(m)-delt(m);endfor m=1:N1JJ(2*m-1,2*m-1)=H(m,m);JJ(2*m-1,2*m)=N(m,m);JJ(2*m,2*m-1)=J(m,m);JJ(2*m,2*m)=L(m,m);endfor m=1:N1for n=1:N1if m=nelseH(m,n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-de
23、lt(n)-B(m,n)*cos(delt(m)-delt(n);J(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n);N(m,n)=-J(m,n);L(m,n)=H(m,n);JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n);JJ(2*m,2*n-1)=J(m,n);JJ(2*m,2*n)=L(m,n);endendendfor m=1:N1PP(2*m-1)=pp(m);PP(2*m)=qq(m);enduu=-inv(JJ)*PP;precision=max(a
24、bs(uu);for n=1:N1delt(n)=delt(n)+uu(2*n-1);邵阳学院课程设计19u(n)=u(n)+uu(2*n);endk=k+1;endK=k-1,delt,u%the following program is used to calculate the S5 and S(m,n)for n=1:N1+1U(n)=u(n)*(cos(delt(n)+j*sin(delt(n);endfor m=1:N1+1I(m)=Y(5,m)*U(m);endS5=U(5)*sum(conj(I)for m=1:N1+1for n=1:N1+1S(m,n)=U(m)*(conj
25、(U(m)-conj(U(n)*conj(-Y(m,n);endendYS邵阳学院课程设计205.3 运行计算结果K =4；delt = -0.0461 -0.0839 -0.0896 -0.1044 0；U = 1.0365 1.0087 1.0073 1.0016 1.0600S5 =1.2982 + 0.2445iY =10.8340 -32.5000i -1.6670 + 5.0000i -1.6670 + 5.0000i -2.5000 + 7.5000i -5.0000 +15.0000i-1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +3
26、0.0000i 0 -1.2500 + 3.7500i-1.6670 + 5.0000i -10.0000 +30.0000i 12.9170 -38.7500i -1.2500 + 3.7500i 0 -2.5000 + 7.5000i 0 -1.2500 + 3.7500i 3.7500 -11.2500i 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 0 6.2500 -18.7500iS =0 0.2469 + 0.0815i 0.2793 + 0.0806i 0.5489 + 0.1333i -0.8751 - 0.0954i-0.2431 - 0.
27、0701i 0 0.1891 - 0.0121i 0 -0.3960 - 0.0677i-0.2746 - 0.0664i -0.1887 + 0.0132i 0 0.0633 + 0.0033i 0 -0.5370 - 0.0977i 0 -0.0630 - 0.0023i 0 0 0.8895 + 0.1387i 0.4087 + 0.1058i 0 0 0 邵阳学院课程设计21总 结此次课程设计首先让我明白了要使电力系统运行的稳定，必须经过精密的设计和计算。在进行课题设计的过程中，加深了我对潮流计算的认识，尤其是对牛顿拉夫逊潮流计算的求解思路有了比较透彻的理解。同时由于求解过程中用到求节
28、点导钠矩阵，求矩阵的逆等等，又使我对以前所学的知识有了一次很好的温习。同时也看到了研究性学习的效果，从研究中去学习，理论结合实际，将理论运用到实际，同时在实践中发现问题，然后解决问题，而且在此次课程设计中，我发现了自己的基础知识有很多的不足。这些基础的缺乏给我的设计计划造成了不小的障碍。在这个过程中，我明白了，只要用心去做，认真去做，持之以恒，就会有新的发现，有意外的收获。此次课程设计的完成，要感谢老师的帮助以及同学的协力合作，团结才能出成绩。 邵阳学院课程设计22参考文献1 何仰赞等.电力系统分析M. 武汉：华中理工大学出版社，2002.3 2 西安交通大学等.电力系统计算M.北京：水利电力出版社，1993.12 3 陈 衍.电力系统稳态分析M.北京：水利电力出版社，2004.1 4 李光琦.电力系统暂态分析M.北京： 水利电力出版社，2002.5 5 于永源，杨绮雯. 电力系统分析（第二版）M. 北京：中国电力出版社，2004.3
