电力系统分析潮流计算课程序设计及其MATLAB程序设计

1、PAGE 6-PAGE 6. z. - - . .可修编. . 电力系统分析潮流计算程序设计报告题目：13节点配电网潮流计算 学 院 电气工程学院 专业班级 学生* 学 号 班内序号 指导教师 房大中 提交日期 2015年05月04日 -. z. . .可修编. .目录TOC o 1-3 h u HYPERLINK l _Toc418380649一、程序设计目的 PAGEREF _Toc418380649 h 2HYPERLINK l _Toc418380650二、程序设计要求 PAGEREF _Toc418380650 h 4HYPERLINK l _Toc418380651三、13节点配网

2、潮流计算 PAGEREF _Toc418380651 h 4HYPERLINK l _Toc4183806523.1主要流程 PAGEREF _Toc418380652 h 4HYPERLINK l _Toc418380653第一步的前推公式如下1-1-1-5： PAGEREF _Toc418380653 h 4HYPERLINK l _Toc418380654第二步的回代公式如下1-61-9： PAGEREF _Toc418380654 h 4HYPERLINK l _Toc4183806553.2配网前推后代潮流计算的原理 PAGEREF _Toc418380655 h 5HYPERLIN

3、K l _Toc4183806563.3配网前推后代潮流计算迭代过程 PAGEREF _Toc418380656 h 5HYPERLINK l _Toc4183806573.3计算原理 PAGEREF _Toc418380657 h 6HYPERLINK l _Toc418380658四、计算框图流程 PAGEREF _Toc418380658 h 7HYPERLINK l _Toc418380659五、确定前推回代支路次序 PAGEREF _Toc418380659 h 7HYPERLINK l _Toc418380660六、前推回代计算输入文件 PAGEREF _Toc418380660

4、h 8HYPERLINK l _Toc418380661主程序： PAGEREF _Toc418380661 h 8HYPERLINK l _Toc418380662输入文件清单： PAGEREF _Toc418380662 h 9HYPERLINK l _Toc418380663计算结果： PAGEREF _Toc418380663 h 10HYPERLINK l _Toc418380664数据分析： PAGEREF _Toc418380664 h 10HYPERLINK l _Toc418380665七、配电网潮流计算的要点 PAGEREF _Toc418380665 h 11HYPERL

5、INK l _Toc418380666八、自我总结 PAGEREF _Toc418380666 h 11HYPERLINK l _Toc418380667九、参考文献 PAGEREF _Toc418380667 h 12HYPERLINK l _Toc418380668附录一 MATLAB的简介 PAGEREF _Toc418380668 h 12-. z.一、程序设计目的开式网络潮流计算：配电网的构造特点呈辐射状，在正常运行时是开环的；配电网的潮流计算采用的方法是前推回代法，本程序利用前推回代法的根本原理、收敛性。(1)在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划网架，选

6、择无功补偿方案，满足规划水平年的大、小方式下潮流交换控制、调峰、调相、调压的要求。 (2)在编制年运行方式时，在预计负荷增长及新设备投运根底上，选择典型方式进展潮流计算，发现电网中薄弱环节，供调度员日常调度控制参考，并对规划、基建部门提出改良网架构造，加快基建进度的建议。 (3)正常检修及特殊运行方式下的潮流计算，用于日运行方式的编制，指导发电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求。及电压质量要求。 (4)预想事故、设备退出运行对静态平安的影响分析及作出预想的运行方式调整方案。 图1 13节点配电网构造图表1 系统支路参数支路R*B/2(TK) S123.36

7、73.6850.0232.3562.5410.0341.1451.280.0454.5245.040.0260.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 系统负荷参数节点编号节点类型节点初始电kVPiMVAQiMVA1根节点10.4002中间节点10.00.03420.03013中间节点10.00.06930.06424中间节点10.00.08450.07635叶节点10.00.02950.02

8、616中间节点10.00.04740.04097叶节点10.00.11760.09578中间节点10.00.09460.08579叶节点10.00.09160.085910叶节点10.00.02710.022911中间节点10.00.06960.064312叶节点10.00.06760.057913叶节点10.00.02980.0242二、程序设计要求看懂前推回代法计算程序；报告表达计算原理及计算流程；绘制计算流程框图；确定前推回代支路次序广度优先，或深度优先并编写前推回代计算输入文件，然后进展潮流计算；整理专利计算结果；总结配电网潮流计算的要点；三、设计内容1.根据电力系统网络推导电力网络数

9、学模型，写出节点导纳矩阵；2.赋予各节点电压变量直角坐标系形式初值后，求解不平衡量；3.形成雅可比矩阵；4.求解修正量后，重新修改初值，从2开场重新循环计算；5.求解的电压变量到达所要求的精度时，再计算各支路功率分布、功率损耗和节点 6.上机编程调试；连调；7.计算分析给定系统潮流分析并与手工计算结果作比拟分析。8.准备计算机演示辩论，书写该课程设计说明书必须计算机打印四、13节点配网潮流计算4.1牛拉法的原理及其根本方程牛顿迭代法Newtons method又称为牛顿-拉夫逊拉弗森方法Newton-Raphson method，它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法

10、。设r是f(*) = 0的根，选取*0作为r初始近似值，过点*0,f(*0)做曲线y = f(*)的切线L，L的方程为y = f(*0) f(*0)(*-*0)，求出L与*轴交点的横坐标 *1 = *0-f(*0)/f(*0)，称*1为r的一次近似值。过点*1,f(*1)做曲线y = f(*)的切线，并求该切线与*轴的横坐标 *2 = *1-f(*1)/f(*1)，称*2为r的二次近似值。重复以上过程，得r的近似值序列，其中*(n+1)=*(n)f(*(n)/f(*(n)，称为r的n+1次近似值，上式称为牛顿迭代公式。解非线性方程f(*)=0的牛顿法是把非线性方程线性化的一种近似方法。把f(*

11、)在*0点附近展开成泰勒级数 f(*) = f(*0)+(*0)f(*0)+(*0)2*f(*0)/2! + 取其线性局部，作为非线性方程f(*) = 0的近似方程，即泰勒展开的前两项，则有f(*0)+f(*0)(*0)=f(*)=0 设f(*0)0则其解为*1=*0f(*0)/f(*0) 这样，得到牛顿法的一个迭代序列：*(n+1)=*(n)f(*(n)/f(*(n)。4.2 PQ分解法的原理及其根本方程PQ分解法的根本思想是根据电力系统实际运行特点：通常网络上的电抗远大于电阻值，则系统母线电压幅值的微小变化对母线有功功率的改变影响很小。同样，母线电压相角的少许改变 ，也不会引起母线有功功率

12、的明显改变 ，因此，节点功率方程在用极坐标形式表示时，它的修正方程式可简化为 2-1这就是把2n-1阶的线性方程组变成了两个n-1阶的线性方程组，将P和Q分开来进展迭代计算，因而大大地减少了计算工作量。但是H、L在迭代过程中仍然在不断的变化而且又都是不对称矩阵。对牛顿法的进一步简化，即把式2-1中的系数矩阵简化为在迭代过程中不变的对称矩阵。在一般情况下，线路两端的电压相角 是不大的不超过1020，因此，可以认为 2-2此外，与系统各节点无功功率相应的导纳远小于该节点自导纳的虚部，即因而 2-3考虑到以上关系，式2-1的系数矩阵中的个元素可以表示为 ( i , j=1,2, n-1) ( i ,

13、 j=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配网前推后代潮流计算的原理前推回代法在配电网潮流计算中简单实用，所有的数据都是以矢量形式存储，因此节省了大量的计算机内存，对于任

14、何种类的配电网只要有合理的 R/* 值，此方法均可保证收敛。算法的稳定性也是评价配电网潮流算法的重要指标。一般情况下，算法的收敛阶数越高，算法的稳定性越差，前推回代法的收敛阶数为一阶，因此它也具有较好的稳定性。比拟而言，前推回代法充分利用了网络呈辐射状的构造特点，数据处理简单，计算效率高，具有较好的收敛性，被公认是求解辐射状配电网潮流问题的最正确算法之一。4.4配网前推后代潮流计算迭代过程配网潮流前推回代潮流算法第n+1步的迭代公式如下：节点i的前推计算公式为由图1可知节点i+1的注入有功功率和无功功率分别为：NBNB-1P(i+1)= PL(j)+ LP(j)Q(i+1)= QL(j)+ L

15、Q(j)其中：i=1，2，。NB-1，NB为节点数，PL(i)为第i节点所带负荷有功功率；QL(i)为第i节点所带负荷无功功率；LP(i)为第i条线段上的有功功率损耗；LQ(i)为第i条线段上的无功功率损耗。收敛条件以前后两次迭代的电压偏差作为迭代收敛条件，(8)式说明节点电压幅值最大偏差小于设定阈值，即认为迭代收敛，则迭代过程完毕。ma*|Vk(i)-Vk-1(i)|0.0001 PQ2=PQ; %每一次迭代各节点的注入有功和无功一样 PL=0.0;for i=1:NB kf=FT(i,1); %前推始节点号 kt=FT(i,2); %前推终节点号 *=(PQ2(kf,2)2+PQ2(kf,

16、3)2)/(V(kf)2); %计算沿线电流平方A PQ1(i,1)= R*(i,1)*+PQ2(kf,2); %计算支路首端有功/MW R*(i,1)R PQ1(i,2)= R*(i,2)*+PQ2(kf,3); %计算沿支路的无功损耗/Mvar R*(i,2)* 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+R*(i,1)*; %累积计算整个网络的功率损耗endangle(1)=0.0;for i=NB:

17、-1:1 kf=FT(i,2); %回代始节点号 kt=FT(i,1); %回代终节点号 dv1=(PQ1(i,1)*R*(i,1)+PQ1(i,2)*R*(i,2)/V(kf); %计算支路电压损耗的纵分量dv1 dv2=(PQ1(i,1)*R*(i,2)-PQ1(i,2)*R*(i,1)/V(kf); %计算支路电压损耗的横分量dv2 V2(kt)=sqrt(V(kf)-dv1)2+(dv2)2); %计算支路末端电压/kV angle(kt)= atand(dv2/(V(kf)-dv1)+angle(kf); %计算支路电压相角单位度endma*d=abs(V2(2)-V(2);V2(1

18、)=V(1);for i=3:1:NN if abs(V2(i)-V(i)ma*d; ma*d=abs(V2(i)-V(i); end endma*dk=k+1PQ1 %潮流分布 即支路首端潮流MVAV=V2 %节点电压模计算结果kVangle %节点电压角度计算结果单位度PL %网损单位MWendclear输入文件清单：function PQ,FT,R*=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 0.0474 0.040

19、910.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;R*=% R * 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.856 1.143.74

20、3 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 9.6444angle

21、 = 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以内，计算到此完毕。七、配电网潮流计算的要点需要注意系统中所有母线的电压是否在允许的*围以内。系统中各种元件(线路、变压器等)是否会出现过负荷。在进展电压和功率分布计算以前，先要对网络等值电路作些简化处理。确立前推回代的计算次序。在开场编写程序之前，要了解matlab仿真软件，包括其程序的编写，程序的运行，我们都要清楚的知道为了提高计算精度，

22、重复计算，给出一个容许误差作为计算的收敛的判据。八、自我总结此次课程设计首先让我明白了要使电力系统运行的稳定，必须经过精细的设计和计算。在进展课题设计的过程中，加深了我对潮流计算的认识，尤其是对牛顿拉夫逊潮流计算的求解思路有了比拟透彻的理解。在此次程序设计中，我在MATLAB编程，13节点配电网潮流计算，WORD文档的编辑方面均有所提高，之前学过的潮流计算仅限于书本，通过完成这个作业使我对这块知识了解得更加全面，但也暴露出了一些问题：理论知识还是缺乏，根本知识不充足，做起这个大作业难度还是很大的，对MATLAB的性能和特点还不能有一个全面的把握，相信通过以后的学习能弥补这些缺乏，从而到达一个新

23、的层次。在过程设计中我发现了一些要点，要求得支路的首端功率，从与叶节点联接的支路开场，该支路的末端功率即等于叶节点功率，利用这个功率和对应的节点电压计算支路功率损耗，才能求首段功率。还有在迭代计算开场之前，先要确定支路的计算顺序问题，即确定好是深度优先还是广度优先。要求得各支路终结点的电压的话利用上述所得的支路首端功率和本步骤刚算出的本支路始节点的电压，从电源点开场逐条支路进展计算。而且在此次课程设计中，我发现了自己的根底知识有很多的缺乏。这些根底的缺乏给我的设计方案造成了不小的障碍。在这个过程中，我明白了，只要用心去做，认直去做，持之以恒，就会有新的发现，有意外的收获。虽然很多错误被克制了，但是系统中难免还有很多缺乏之处，有些地方感觉做的不是很正确，希望在以后的学习中能够在同学们的帮助和教师的指导下做得更好。九、参考文献1何仰赞 温增银. 电力系统分析. 华中科技大学，20022蔡旭晖、*卫国、蔡立燕.MATLAB根底与应用教程D.：人民邮电2009 附录一 MATLAB的简介MATLAB是由美国mathworks公司发布的主要面对科学计算、可视