1、电力系统分析潮流计算课程序设计方案及其MATLAB程序设计方案 电力系统分析潮流计算程序设计报告题目：13节点配电网潮流计算 学 院 电气工程学院 专业班级 学生姓名 学 号 班内序号 指导教师 房大中 提交日期 2015年05月04日 一、程序设计目的开式网络潮流计算：配电网的结构特点呈辐射状，在正常运行时是开环的；配电网的潮流计算采用的方法是前推回代法，本程序利用前推回代法的基本原理、收敛性。(1)在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划网架，选择无功补偿方案，满足规划水平年的大、小方式下潮流交换控制、调峰、调相、调压的要求。 (2)在编制年运行方式时，在预计负荷增
2、长及新设备投运基础上，选择典型方式进行潮流计算，发现电网中薄弱环节，供调度员日常调度控制参考，并对规划、基建部门提出改进网架结构，加快基建进度的建议。 (3)正常检修及特殊运行方式下的潮流计算，用于日运行方式的编制，指导发电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求。及电压质量要求。 (4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 图1 13节点配电网结构图表1 系统支路参数支路R（）X（）B/2(TK) （S）123.3673.6850.0232.3562.5410.0341.1451.280.0454.5245.040.0260
3、.8561.140.0672.7452.9650.0283.7434.2510.0892.2372.7560.03104.144.6960.03111.3281.7630.011122.4362.8660.04133.5213.9660.0表2 系统负荷参数节点编号节点类型节点初始电（kV）Pi（MVA）Qi（MVA）1根节点10.4002中间节点10.00.03420.03013中间节点10.00.06930.06424中间节点10.00.08450.07635叶节点10.00.02950.02616中间节点10.00.04740.04097叶节点10.00.11760.09578中间节点1
4、0.00.09460.08579叶节点10.00.09160.085910叶节点10.00.02710.022911中间节点10.00.06960.064312叶节点10.00.06760.057913叶节点10.00.02980.0242二、程序设计要求1. 看懂前推回代法计算程序；2. 报告叙述计算原理及计算流程；3. 绘制计算流程框图；4. 确定前推回代支路次序（广度优先，或深度优先）并编写前推回代计算输入文件，然后进行潮流计算；5. 整理专利计算结果；6. 总结配电网潮流计算的要点；三、设计内容1.根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；2.赋予各节点电压变量（直角坐标
5、系形式）初值后，求解不平衡量；3.形成雅可比矩阵；4.求解修正量后，重新修改初值，从2开始重新循环计算；5.求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和节点6.上机编程调试；连调；7.计算分析给定系统潮流分析并与手工计算结果作比较分析。8.准备计算机演示答辩，书写该课程设计说明书（必须计算机打印）四、13节点配网潮流计算4.1牛拉法的原理及其基本方程牛顿迭代法（Newtons method）又称为牛顿-拉夫逊（拉弗森）方法（Newton-Raphson method），它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。设r是f(x) = 0的根，选取x0作
6、为r初始近似值，过点（x0,f(x0)）做曲线y = f(x)的切线L，L的方程为y = f(x0) f(x0)(x-x0)，求出L与x轴交点的横坐标 x1 = x0-f(x0)/f(x0)，称x1为r的一次近似值。过点（x1,f(x1)）做曲线y = f(x)的切线，并求该切线与x轴的横坐标 x2 = x1-f(x1)/f(x1)，称x2为r的二次近似值。重复以上过程，得r的近似值序列，其中x(n+1)=x(n)f(x(n)/f(x(n)，称为r的n+1次近似值，上式称为牛顿迭代公式。解非线性方程f(x)=0的牛顿法是把非线性方程线性化的一种近似方法。把f(x)在x0点附近展开成泰勒级数 f
7、(x) = f(x0)+(xx0)f(x0)+(xx0)2*f(x0)/2! + 取其线性部分，作为非线性方程f(x) = 0的近似方程，即泰勒展开的前两项，则有f(x0)+f(x0)(xx0)=f(x)=0 设f(x0)0则其解为x1=x0f(x0)/f(x0) 这样，得到牛顿法的一个迭代序列：x(n+1)=x(n)f(x(n)/f(x(n)。4.2 PQ分解法的原理及其基本方程PQ分解法的基本思想是根据电力系统实际运行特点：通常网络上的电抗远大于电阻值，则系统母线电压幅值的微小变化对母线有功功率的改变影响很小。同样，母线电压相角的少许改变，也不会引起母线有功功率的明显改变，因此，节点功率方
8、程在用极坐标形式表示时，它的修正方程式可简化为 （2-1）这就是把2（n-1）阶的线性方程组变成了两个n-1阶的线性方程组，将P和Q分开来进行迭代计算，因而大大地减少了计算工作量。但是H、L在迭代过程中仍然在不断的变化而且又都是不对称矩阵。对牛顿法的进一步简化，即把式（2-1）中的系数矩阵简化为在迭代过程中不变的对称矩阵。在一般情况下，线路两端的电压相角是不大的（不超过1020），因此，可以认为 （2-2）此外，与系统各节点无功功率相应的导纳远小于该节点自导纳的虚部，即因而 （2-3）考虑到以上关系，式（2-1）的系数矩阵中的个元素可以表示为 ( i , j=1,2, n-1) ( i , j
9、=1,2, m)而系数矩阵H和L则可以分别写成： (2-4) (2-5)将式（2-4）和（2-5）代人式中，得到用 和 分别左乘以上两式，便得 (2-6) (2-7)这就是简化了的修正方程式，它们也可以展开写成 (2-8) (2-9)在这两个修正方程式中系数矩阵元素就是系统导纳矩阵的虚部，因而系数矩阵是对称矩阵，且在迭代过程中保持不变。这就大大减少了计算工作量。用极坐标表示的节点功率增量为 (2-10)式（2-8）、（2-9）和（2-10）构成了PQ分解法迭代过程的基本方程式。4.3配网前推后代潮流计算的原理前推回代法在配电网潮流计算中简单实用，所有的数据都是以矢量形式存储，因此节省了大量的计
10、算机内存，对于任何种类的配电网只要有合理的 R/X 值，此方法均可保证收敛。算法的稳定性也是评价配电网潮流算法的重要指标。一般情况下，算法的收敛阶数越高，算法的稳定性越差，前推回代法的收敛阶数为一阶，因此它也具有较好的稳定性。比较而言，前推回代法充分利用了网络呈辐射状的结构特点，数据处理简单，计算效率高，具有较好的收敛性，被公认是求解辐射状配电网潮流问题的最佳算法之一。4.4配网前推后代潮流计算迭代过程配网潮流前推回代潮流算法第n+1步的迭代公式如下：节点i的前推计算公式为由图1可知节点i+1的注入有功功率和无功功率分别为：NBNB-1P(i+1)= PL(j)+ LP(j)Q(i+1)= Q
11、L(j)+ LQ(j)其中：i=1，2，。NB-1，NB为节点数，PL(i)为第i节点所带负荷有功功率；QL(i)为第i节点所带负荷无功功率；LP(i)为第i条线段上的有功功率损耗；LQ(i)为第i条线段上的无功功率损耗。收敛条件以前后两次迭代的电压偏差作为迭代收敛条件，(8)式表明节点电压幅值最大偏差小于设定阈值，即认为迭代收敛，则迭代过程结束。max|Vk(i)-Vk-1(i)|0.0001 PQ2=PQ。 %每一次迭代各节点的注入有功和无功相同 PL=0.0。for i=1:NB kf=FT(i,1)。 %前推始节点号 kt=FT(i,2)。 %前推终节点号 x=(PQ2(kf,2)2+
12、PQ2(kf,3)2)/(V(kf)2)。 %计算沿线电流平方A PQ1(i,1)= RX(i,1)*x+PQ2(kf,2)。 %计算支路首端有功/MW RX(i,1)R PQ1(i,2)= RX(i,2)*x+PQ2(kf,3)。 %计算沿支路的无功损耗/Mvar RX(i,2)X PQ2(kt,2)= PQ2(kt,2)+PQ1(i,1)。 %用PQ1去修正支路末端节点的有功P单位MW PQ2(kt,3)= PQ2(kt,3)+PQ1(i,2)。 %用PQ1去修正支路末端节点的有功Q 单位Mvar PL=PL+RX(i,1)*x。 %累积计算整个网络的功率损耗endangle(1)=0.0
13、。for i=NB:-1:1 kf=FT(i,2)。 %回代始节点号 kt=FT(i,1)。 %回代终节点号 dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2)/V(kf)。 %计算支路电压损耗的纵分量dv1 dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1)/V(kf)。 %计算支路电压损耗的横分量dv2 V2(kt)=sqrt(V(kf)-dv1)2+(dv2)2)。 %计算支路末端电压/kV angle(kt)= atand(dv2/(V(kf)-dv1)+angle(kf)。 %计算支路电压相角（单位度）endmaxd=abs(V2(
14、2)-V(2)。V2(1)=V(1)。for i=3:1:NN if abs(V2(i)-V(i)maxd。 maxd=abs(V2(i)-V(i)。 end endmaxdk=k+1PQ1 %潮流分布 即支路首端潮流MVAV=V2 %节点电压模计算结果kVangle %节点电压角度计算结果单位度PL %网损单位MWendclear输入文件清单：function PQ,FT,RX=case115()PQ=%节点电压 有功 无功10.4 0 010.0 0.0342 0.030110.0 0.0693 0.064210.0 0.0845 0.076310.0 0.0295 0.026110.0
15、0.0474 0.040910.0 0.1176 0.095710.0 0.0946 0.085710.0 0.0916 0.085910.0 0.0271 0.022910.0 0.0696 0.064310.0 0.0676 0.057910.0 0.0298 0.0242。FT=%首端 末端 5 413 412 114 310 3 11 37 69 83 26 28 22 1。RX=% R X 4.524 5.043.521 3.9662.436 2.8661.145 1.284.14 4.6961.328 1.7632.745 2.9652.237 2.7562.356 2.5410.
16、856 1.143.743 4.2513.367 3.685。计算结果：k =8PQ1 = 0.0296 0.0262 0.0299 0.0243 0.0678 0.0581 0.1444 0.1273 0.0272 0.0230 0.1379 0.1231 0.1183 0.0964 0.0920 0.0864 0.3852 0.3445 0.1661 0.1379 0.1891 0.1750 0.8116 0.7278V =10.4000 9.8795 9.6991 9.6652 9.6377 9.8492 9.7872 9.7326 9.6870 9.6764 9.6578 9.6235
17、 9.6444angle = 0 0.3011 0.4011 0.4250 0.4438 0.3431 0.3942 0.3899 0.4266 0.4209 0.4498 0.4823 0.4452PL =0.0488数据分析：经过8次迭代，各节点的电压误差均在0.0001kv以内，计算到此结束。七、配电网潮流计算的要点1. 需要注意系统中所有母线的电压是否在允许的范围以内。2. 系统中各种元件(线路、变压器等)是否会出现过负荷。3. 在进行电压和功率分布计算以前，先要对网络等值电路作些简化处理。4. 确立前推回代的计算次序。5. 在开始编写程序之前，要了解matlab仿真软件，包括其程序的
18、编写，程序的运行，我们都要清楚的知道6. 为了提高计算精度，重复计算，给出一个容许误差作为计算的收敛的判据。八、自我总结此次课程设计首先让我明白了要使电力系统运行的稳定，必须经过精密的设计和计算。在进行课题设计的过程中，加深了我对潮流计算的认识，尤其是对牛顿拉夫逊潮流计算的求解思路有了比较透彻的理解。在此次程序设计中，我在MATLAB编程，13节点配电网潮流计算，WORD文档的编辑方面均有所提高，之前学过的潮流计算仅限于书本，通过完成这个作业使我对这块知识了解得更加全面，但也暴露出了一些问题：理论知识还是不足，基本知识不充足，做起这个大作业难度还是很大的，对MATLAB的性能和特点还不能有一个
19、全面的把握，相信通过以后的学习能弥补这些不足，从而达到一个新的层次。在过程设计中我发现了一些要点，要求得支路的首端功率，从与叶节点联接的支路开始，该支路的末端功率即等于叶节点功率，利用这个功率和对应的节点电压计算支路功率损耗，才能求首段功率。还有在迭代计算开始之前，先要确定支路的计算顺序问题，即确定好是深度优先还是广度优先。要求得各支路终结点的电压的话利用上述所得的支路首端功率和本步骤刚算出的本支路始节点的电压，从电源点开始逐条支路进行计算。而且在此次课程设计中，我发现了自己的基础知识有很多的不足。这些基础的缺乏给我的设计计划造成了不小的障碍。在这个过程中，我明白了，只要用心去做，认直去做，持
20、之以恒，就会有新的发现，有意外的收获。虽然很多错误被克服了，但是系统中难免还有很多不足之处，有些地方感觉做的不是很正确，希望在以后的学习中能够在同学们的帮助和老师的指导下做得更好。九、参考文献1何仰赞 温增银. 电力系统分析. 华中科技大学出版社，20022蔡旭晖、刘卫国、蔡立燕.MATLAB基础与应用教程D.北京：人民邮电出版社2009 附录一 MATLAB的简介MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究
21、、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解
22、算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且mathwork也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C+ ，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。MATLAB 产品族可以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理 技术 数字信号处理 技术 通讯系统设计与仿真 财务与金融工程 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特定类型问题。
