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

1、广州学院Guangzhou College of South China University of Icchnology电力系统分析潮流计算程序设计报告题目：13节点配电网潮流计算学 院电气工程学院专业班级学生姓名学 号班内序号指导教师房大中提交日期 2015 年05月04日一、程序设计目的2二、程序设计要求4三、13节点配网潮流计算 43.1主要流程 错误!未定义书签。3.1.1第一步的前推公式如下（1-1）-（ 1-5）： 错误!未定义书签。3.1.2第二步的回代公式如下（1-6）（ 1-9）: 错误!未定义书签。3.2配网前推后代潮流计算的原理 73.3配网前推后代潮流计算迭代过程 8

2、3.3计算原理9四、计算框图流程10五、 确定前推回代支路次序 错误!未定义书签。六、 前推回代计算输入文件 11主程序：11输入文件清单： 12计算结果：13数据分析：13七、配电网潮流计算的要点 14八、自我总结14九、参考文献 15附录一 MATLAB的简介151华南理工大学广州学院电气工程学院课程设计报告18一、程序设计目的开式网络潮流计算：配电网的结构特点呈辐射状，在正常运行时是开环的；配电网 的潮流计算采用的方法是前推回代法，本程序利用前推回代法的基本原理、收敛性。 在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划网架, 选择无功补偿方案，满足规划水平年的大、小方式

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

4、S)123.3673.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系统负荷参数节点编号节点类型节点初始电（kV）Pi (MVA)Qi (MVA)1根节点10.4002中间节点10.00.03420.03013中间节点10.00.06930.06424中间节点10.00.08450.07635叶节

5、点10.00.02950.02616中间节点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. 看懂前推回代法计算程序；2. 报告叙述计算原理及计算流程；3. 绘制计算流程框图；4. 确定前推回代支路次序(广度优先，或深度优先)并编写前推回代计算输入文 件，然后进行潮流计算；5. 整理专利计算结果；6.

6、总结配电网潮流计算的要点；三、设计内容1. 根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；2. 赋予各节点电压变量(直角坐标系形式)初值后，求解不平衡量；3. 形成雅可比矩阵；4. 求解修正量后，重新修改初值，从2开始重新循环计算；5. 求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和节点6. 上机编程调试；连调；7. 计算分析给定系统潮流分析并与手工计算结果作比较分析。8. 准备计算机演示答辩，书写该课程设计说明书(必须计算机打印)四、13节点配网潮流计算4.1牛拉法的原理及其基本方程牛顿迭代法(Newt ons method)又称为牛顿-拉夫逊(拉弗森)方法(N

7、ewto n-R aphson method)，它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。设r是f(x) = 0的根，选取x0作为r初始近似值，过点(xO,f(xO)做曲线y = f(x) 的切线L， L的方程为y = f(x0) f(xO)(x-xO)，求出L与x轴交点的横坐标x1 =xO-f(xO)/f(xO)，称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

8、的n+1次近似值，上式称为牛顿迭代公式。解非线性方程f(x)=O的牛顿法是把非线性方程线性化的一种近似方 法。把f(x)在xO点附近展开成泰勒级数f(x) = f(xO)+(x xO)f(xO)+(x xOF2*f(xO)/2! + 取其线性部分，作为非线性方程 f(x) = 0 的近似方程，即泰勒展开的前两项，则有 f(xO)+f(xO)(x xO)=f(x)=O 设f(xO)工0则其解为x仁xO f(xO)/f(xO)这样，得到 牛顿 法的一个 迭代序列：x(n+1)=x(n)f(x( n)/f(x( n)4.2 PQ分解法的原理及其基本方程PQ分解法的基本思想是根据电力系统实际运行特点：

9、通常网络上的电抗远大于电阻值，则系统母线电压幅值的微小变化U对母线有功功率的改变P影响很小。同样，母线电压相角的少许改变，也不会引起母线有功功率的明显改变Q ，因此，节点 功率方程在用极坐标形式表示时，它的修正方程式可简化为P H OQ O L U /U(2-1 )这就是把2 (n-1 )阶的线性方程组变成了两个n-1阶的线性方程组，将P和Q分开来进行迭代计算，因而大大地减少了计算工作量。但是H、L在迭代过程中仍然在不断的变化而且又都是不对称矩阵。对牛顿法的进一步简化，即把式(2-1)中的系数矩阵简化为在迭代过程中不变的对称矩阵。在一般情况下，线路两端的电压相角j是不大的(不超过1O2O)，因

10、此，COS ij 1G sin “B可以认为GjSin ijBij(2-2)此外，与系统各节点无功功率相应的导纳BLDi远小于该节点自导纳的虚部，即BiQiU2因而(2-3)考虑到以上关系，式(2-1 )的系数矩阵中的个元素可以表示为HjUiUjBj( i 戶,2,n-1)LjUjUjBj(i , j=1,2,m)而系数矩阵H和L则可以分别写成:U1B11U1U1B12U2U 1B1,n 1U n 1U2 B21U1U2B22UU 2 B2,n 1U n 1U1U2B21B2B22.Bin 1.B2,n 1U1U2Un 1Bn 1,1Bn 1,2. Bn 1n 1Un1UD1B UD 1(2-

11、4)U1B11U11U 1 B12U 2U 1 B1 mU mLU2B21U1U 2 B22U 2U 2B2mU mU m Bm1U1 U m Bm2U 2U m BmmU mU1B11B12.B1mU1U2B21B22.B2mU2U mBm1Bm2.BmmUmUD2B UD 2(2-5)YG“hjijjBj将式(2-4 )和(2-5)代人式UiUiej i Ui(cos ijsin J中，得到PUD1B Ud1QUd2B U用1UD1和1UD2分别左乘以上两式，便得U D11PB Ud1(2-6)U n 1 Bn 1,1U1U n 1Bn 1,2U 2U n 1 Bn 1,n 1U n11U

12、d2 Q B U(2-7)这就是简化了的修正方程式，它们也可以展开写成RU1B11B12B1,n 1U11P2B21b22B2,n 1U22U2Pn 1Bn 1,1Bn 1,2Bn 1,n 1U n 1n 1Un1Q1(2-8)U1B11B12B1, mU1Q2B21B22B2,mU2U2QmBm ,1Bm,2Bm ,mUmUm(2-9)因而系数矩阵是对(2-10)在这两个修正方程式中系数矩阵元素就是系统导纳矩阵的虚部, 称矩阵，且在迭代过程中保持不变。这就大大减少了计算工作量。用极坐标表示的节点功率增量为nP Rs Ui Uj(Gj cos ij Bsin j) 0j inQi Qis U

13、i U j (Gij cos ij Bsin ij) 0j i式（2-8 ）、（ 2-9 ）和（2-10）构成了 PQ分解法迭代过程的基本方程式。4.3配网前推后代潮流计算的原理前推回代法在配电网潮流计算中简单实用，所有的数据都是以矢量形式存储，因此节省了大量的计算机内存，对于任何种类的配电网只要有合理的R/X值，此方法均可 保证收敛。算法的稳定性也是评价配电网潮流算法的重要指标。一般情况下，算法的收 敛阶数越高，算法的稳定性越差，前推回代法的收敛阶数为一阶，因此它也具有较好的 稳定性。比较而言，前推回代法充分利用了网络呈辐射状的结构特点，数据处理简单， 计算效率高，具有较好的收敛性，被公认是

14、求解辐射状配电网潮流问题的最佳算法之一。4.4配网前推后代潮流计算迭代过程配网潮流前推回代潮流算法第n+1步的迭代公式如下：节点i的前推计算公式为由图1可知节点i+1的注入有功功率和无功功率分别为：NB NB-1P(i+1)= 刀 PLjE LP(j)Q(i+1)= 刀 QL+)刀 LQ(j)其中：i=1, 2,。NB-1，NB为节点数，PL(i)为第i节点所带负荷有功功率；QL(i) 为第i节点所带负荷无功功率；LP(i)为第i条线段上的有功功率损耗；LQ(i)为第i条线 段上的无功功率损耗。收敛条件以前后两次迭代的电压偏差作为迭代收敛条件，(8)式表明节点电压幅值最大偏差小 于设定阈值，即

15、认为迭代收敛，则迭代过程结束。max|Vk(i)-Vk- 1(i)|0.0001PQ2=PQ;PL=0.0;for i=1:NBkf=FT(i,1);kt=FT(i,2);%调用数据文件%节点数%支路数数%V初始电压相量%每一次迭代各节点的注入有功和无功相同%前推始节点号%前推终节点号x=(PQ2(kf,2)A2+PQ2(kf,3)A2)/(V(kf)A2); %PQ1(i,1)= RX(i,1)*x+PQ2(kf,2);PQ1(i,2)= RX(i,2)*x+PQ2(kf,3);PQ2(kt,2)= PQ2(kt,2)+PQ1(i,1);PQ2(kt,3)= PQ2(kt,3)+PQ1(i,

16、2); PL=PL+RX(i,1)*x;endan gle(1)=0.0;for i=NB:-1:1计算沿线电流平方 A%计算支路首端有功/MW RX(i,1)R %计算沿支路的无功损耗/Mvar RX(i,2)X%用PQ1去修正支路末端节点的有功%用PQ1去修正支路末端节点的有功%累积计算整个网络的功率损耗P单位MWQ单位Mvar%回代始节点号%回代终节点号kf=FT(i,2);kt=FT(i,1);dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2)/V(kf); dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1)/V(kf);V2(kt)

17、=sqrt(V(kf)-dv1)A2+(dv2)A2);an gle(kt)= ata nd(dv2/(V(kf)-dv1)+a ngle(kf);end%计算支路电压损耗的纵分量%计算支路电压损耗的横分量%计算支路末端电压/kV%计算支路电压相角(单位度)dv1dv2maxd=abs(V2 (2)-V(2);V2(1)=V(1);for i=3:1:NNif abs(V2(i)-V(i)maxd;maxd=abs(V2(i)-V(i);endend maxd k=k+1PQ1V=V2%潮流分布即支路首端潮流%节点电压模计算结果 kVMV Aan glePL end clear%节点电压角度计

18、算结果单位度%网损单位MW输入文件清单：fun ctio n PQ,FT,RX=case115()PQ=%节点电压有功无功10.40010.00.03420.030110.00.06930.064210.00.08450.076310.00.02950.026110.00.04740.040910.00.11760.095710.00.09460.085710.00.09160.085910.00.02710.022910.00.06960.064310.00.06760.057910.00.02980.0242;FT=%首端末端5 413412114310 311 37 698326 28

19、22 1;RX=% R X4.5245.043.5213.9662.4362.8661.1451.284.144.6961.3281.7632.7452.9652.2372.7562.3562.5410.8561.143.7434.2513.3673.685;计算结果:k =8PQ1 =0.02960.02990.06780.14440.02720.13790.11830.09200.38520.16610.18910.02620.02430.05810.12730.02300.12310.09640.08640.34450.13790.17500.72780.8116V =10.40009.

20、87959.69919.66529.63779.84929.78729.73269.68709.67649.65789.62359.6444an gle = 00.30110.40110.42500.44380.34310.39420.38990.42660.42090.44980.48230.4452PL =0.0488数据分析:经过8次迭代，各节点的电压误差均在0.0001kv以内，计算到此结束。七、配电网潮流计算的要点1. 需要注意系统中所有母线的电压是否在允许的范围以内。2. 系统中各种元件（线路、变压器等）是否会出现过负荷。3. 在进行电压和功率分布计算以前，先要对网络等值电路作些简

21、化处理。4. 确立前推回代的计算次序。5. 在开始编写程序之前，要了解 matlab仿真软件，包括其程序的编写，程序的运 行，我们都要清楚的知道6. 为了提高计算精度，重复计算，给出一个容许误差作为计算的收敛的判据。八、自我总结此次课程设计首先让我明白了要使电力系统运行的稳定，必须经过精密的设计和计算。在进行课题设计的过程中，加深了我对潮流计算的认识，尤其是对牛顿拉夫逊潮流 计算的求解思路有了比较透彻的理解。在此次程序设计中，我在 MATLAB程，13节点配电网潮流计算，WOR文档的编辑 方面均有所提高，之前学过的潮流计算仅限于书本，通过完成这个作业使我对这块知识 了解得更加全面，但也暴露出了

22、一些问题：理论知识还是不足，基本知识不充足，做起 这个大作业难度还是很大的，对 MATLAB勺性能和特点还不能有一个全面的把握，相信 通过以后的学习能弥补这些不足，从而达到一个新的层次。在过程设计中我发现了一些要点，要求得支路的首端功率，从与叶节点联接的支路 开始，该支路的末端功率即等于叶节点功率，利用这个功率和对应的节点电压计算支路 功率损耗，才能求首段功率。还有在迭代计算开始之前，先要确定支路的计算顺序问题， 即确定好是深度优先还是广度优先。 要求得各支路终结点的电压的话利用上述所得的支 路首端功率和本步骤刚算出的本支路始节点的电压，从电源点开始逐条支路进行计算。而且在此次课程设计中，我发

23、现了自己的基础知识有很多的不足。这些基础的缺乏 给我的设计计划造成了不小的障碍。在这个过程中，我明白了，只要用心去做，认直去 做，持之以恒，就会有新的发现，有意外的收获。虽然很多错误被克服了，但是系统中难免还有很多不足之处，有些地方感觉做的不 是很正确，希望在以后的学习中能够在同学们的帮助和老师的指导下做得更好。九、参考文献1 何仰赞温增银电力系统分析华中科技大学出版社，20022 蔡旭晖、刘卫国、蔡立燕 MATLAB基础与应用教程D.北京：人民邮电出版社2009附录一 MATLAB的简介MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式 程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及