华中科技大学现代电力系统分析潮流计算作业.docx

1、华中科技大学现代电力系统分析潮流计算作业现代电力系统分析作业基于Matpower的电力系统潮流计算专业:班级:姓名:学号:基于Matpower的电力系统潮流计算 1.本次潮流计算的目的及意义 2.电力系统潮流计算及其意义 3.电力系统潮流计算常规方法 牛顿-拉夫逊法 节点电压用直角坐标表示时的牛顿一拉夫逊潮流计算 牛顿一拉夫逊法潮流计算程序框图 4.选用的潮流计算的系统 5.利用软件matpower计算潮流 matpower 简介 基态潮流计算 基态潮流计算条件说明及数据输入 潮流计算结果 基态潮流结果计算分析 最优潮流计算 最优潮流计算条件说明及数据输入 最优潮流的理论结果 最优潮流仿真计算

2、结果 最优潮流结果计算分析 6.感想与小结 参考文献 附录一、基态潮流计算文件 附录二、最优潮流计算文件 基于 Matpower 的电力系统潮流计算1.本次潮流计算的目的及意义 本次潮流计算的目的及意义主要是了解电力系统潮流计算及其意义， 在此基础上，了解电力系统潮流计算的模型以及常规的潮流计算的方法， 掌握并熟练使用电力系统潮流计算软件 Matpower。最后，利用Matpower 计算美国西部电网 WSCC三机九节点系统的静态潮流及最优潮流并给出 分析。2.电力系统潮流计算及其意义电力系统潮流计算是研究电力系统稳态运行情况的基本电气计算， 电 力系统潮流计算的任务是根据给定的网络结构及运行

3、条件，求出电网的 运行状态，其中包括各母线的电压、 各支路的功率分布以及功率损耗等。潮流计算分为离线计算和在线计算两大类。 离线计算可以用于电力系统调度，并确定系统的运行方式；离线计算的结果还可以用于电力系统 规划方案的分析以及优化系统的运行状态；此外离线潮流计算可以作为 初值，用于配合系统的故障分析以及稳定性分析。在线计算可以用于能 量管理系统EMS中，主要用作静态分析；在线计算结果还可以帮助调度 员合理安排系统潮流以及进行系统的安全性分析；此外在线潮流计算还 可以帮助优化系统的潮流及运行状态。3.电力系统潮流计算常规方法牛顿- 拉夫逊法 电力系统潮流计算本质上是多元非线性方程组的求解问题，

4、而求解 的方法有很多种。最常见的方法则是牛顿拉夫逊法、高斯 - 塞德尔法、快速解耦法及直流潮流法牛顿一拉夫逊法能有效地解决非线性方程组，有较好的收敛性；它 利用了导纳矩阵的对称性、稀疏性以及节点标号顺序优化的技巧，已经 成为电力系统潮流计算的最常用方法。因此，本次潮流计算将重点了解 牛顿一拉夫逊法。牛顿一拉夫逊法实质上就是切线法，是一种逐步线性化的方法。设有非线性方程组其近似解为x10)，x20)，1 X1X 丄，Xn02 X1，X2，L，Xn0L Ln X1，X2，L，Xn0将上式(1.2)中的多元函数在初值附近分别展开成泰勒级数。 则可得Xn忽略上式(1.3)中乂0) , x20)，x,的

5、二次及以上阶次的各项，整理便可以得到节点电压用直角坐标表示时的牛顿一拉夫逊潮流计算导纳矩阵元素表示为(1.10)(1.11)Yij = Gij + jBij则可得到节点相应的有功功率和无功功率为nnPi ei Gij ej Bij f j fi Gij f j Bij ejj1 j 1nnQi fi Gij ej Bij f j ei Gij f j Bij ejj 1 j1假定系统中的第1,2,，m号节点为PQ节点；假定系统中的第n+1, n+2,，n-1号节点为PV节点，第i个节点的给定功率设为Rs和Q,s, 对节点可列写如下方程nnPi Pis Pi Pis ei Gij ej Bij

6、f j fi Gij fj Bij ej 0,i 1,2,., n 1j 1 j 1nnQi Qis Qi Qis fi Gij ej Bij f j ei Gij f j Bij ej 0,i 1,2,., m (1.12)j 1 j 12 2 2 2 Vi 2 Vis2 ei2 fi2 =0 ,i m 1,.,n 1由于第n号节点为平衡节点，其电压V&= &+ jfn是给定的，故不参加 迭代。式(1.12) 总共包含了 2(n-1) 个变量,而由 (1.9) 可知,待求的变量也 是 2(n-1) 个。还注意到 (1.12) 已经具备了方程组 (1.1) 的形式。因此即 可得到如下的修正方程

7、式W J V (1.13)式 (1.13) 中而 J 则为雅克比矩阵。有了这些, 即可按照上面的牛顿拉夫逊法计算电力系统的潮流。 由 于极坐标下的牛顿拉夫逊潮流计算的推导过程与上类似,这里就不在 详尽介绍牛顿一拉夫逊法潮流计算程序框图牛顿一拉夫逊法潮流计算程序框图如下图所示。图3-1牛顿一拉夫逊法潮流计算程序框图4.选用的潮流计算的系统本次选用的潮流计算系统为美国西部电网 WSC三机九节点系统。该 系统比较经典，经常被业内人士用来作为分析的典型例子。其接线图如 下图所示图4-1三机九节点系统连接图统的基准值和频率取为：系统中设备和线路的参数如表14所示：表4-1 发电机参数标幺值参数名称VN/

8、kVXdX dXqX qXlT doTqOHDG10s0G2186s0G30表2变压器参数的标幺值变压器名称T1T2T3电抗X(pu)变比K111表3线路参数的标幺值线路名称line1lin e2lin e3lin e4Iine5lin e6阻抗RZ(pu)L导纳B(pu)表4负荷参数负荷名称LumpALumpBLumpC有功/MW12590100无功/Mvar5030355.利用软件 matpower 计算潮流 matpower 简介Matpower 更准确地讲，应该是 matlab 中的一个工具箱，而不是一个独立的软件。 它是由卡奈尔大学电气学院电力系统工程研究中心的 RAY ，CARLO

9、餉甘德强在ROBERT THOMA的指导下开发出来的。最初的基于 MATLA的电力潮流和最优潮流代码是为 POWERW项目的需要而编写的。Matpower是一个基于matlab m文件的组建包，用来解决电力潮流和 优化潮流的问题。它致力于为研究人员和教育从业者提供一种易于使用 和可更新的仿真工具。Matpower的设计理念是用尽可能简单、易懂，可 更新的代码来实现最优秀的功能。基态潮流计算基态潮流计算条件说明及数据输入本次潮流计算时选取Bus1为平衡节点，Bus2, Bus3为PV节点，其 余的节点为PQ节点；取发电机G2 G3的输出有功功率为其额定功率， 则依据matpower用户手册即可输

10、入以下数据（数据说明见注释）： function baseMVA,bus,gen,branch,areas,gencost = case9hbyq%定义 case9hbyq 为一个含变压器的三机九节点电力系统的潮流算例version =2;%定义系统的基值为 100MVAbaseMVA = 100;%bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin%母线参数，定义母线的序号、节点类型 （PV、PQ Swing节点）、注入的有功/无功、分区、%母线电压标幺值、基值、相角、电压标幺值的上下限等bus= 130000210 1220000

11、110 18 1320000110 14 1 0 0 0 0 2 1 0 230 15190 30 0 021 0 2301610 0 0 0 110 230171100 35 0 011 0 2301810 0 0 0 110 230191125 50 0 021 0 2301;%bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min %Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf %电机参数，定义电机的有功无功出力，无功输出上下限、端口电压标幺值、电机基值等

12、 gen = 10 0 300 300 1 100 1 250 10 0 0 0 0 0 0 0 0 0 0 0;20 300 300 1 100 1 300 10 0 0 0 0 0 0 0 0 0 0 0;30 300 300 1 100 1 270 10 0 0 0 0 0 0 0 0 0 0 0;%线路参数，定义两母线之间线路的阻抗和导纳值等%branch data%fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax branch = 1 4 0 0 250 250 250 1 0 1 -360 360;

13、45 250 250 250 0 0 1 -360 360;56 150 150 150 0 0 1 -360 360;3 6 0 0 300 300 300 1 0 1 -360 360;67 150 150 150 0 0 1 -360 360;78 250 250 250 0 0 1 -360 360;82 0 0 250 250 250 1 0 1 -360 360;8 9 250 250 250 0 0 1 -360 360;94 250 250 250 0 0 1 -360 360;% OPF Data % area data% area refbus areas = 1 5;%最

14、优潮流计算所需的条件% gen erator cost data% 1 startup shutdow n n x1 y1 . xn yn% 2 startup shutdow n n c(n-1) . c0gen cost =21500035 150;2200003600;23000031 335;；潮流计算结果将上述数据输入到文件中保存并运行后，在命令行窗口中输入run pf(case9hbyq) 指令即可得到基态潮流计算的结果。实际仿真结果见附录一。表5-1母线的潮流计算结果BusVoltageGen eratio nLoadMag(pu)An g(deg)P (MW)Q (MVAr)P

15、 (MW)Q (MVAr)1*-2-3-4-5-90306-7-8-9-Total :表5-2 线路的潮流计算结果BranchFromBusToBusFrom Bus Injectio nTo Bus Injectio nLoss (IA2 * Z)P (MW)Q (MVAr)P (MW)Q (MVAr)P (MW)Q (MVAr)114245346457569678789827939Total基态潮流结果计算分析由基态潮流计算结果可以得到，该三机九节点系统的潮流计算采用的 是牛顿-拉夫逊算法，总共迭代了 4次、共计秒就收敛了，且各母线电压 水平及有功无功功率在约束条件内。其次，可知发电机总共

16、发出的功率为Pg jQG (321.08 j41.31)MVA，而 负荷总共消耗了 Pl jQL (315.00 j115.00)MVA。显然 Pg Pl,Qg Ql。造成有功功率不等的主要原因是输电线路存在电阻，部分有功消耗在 了输电线路上，由潮流计算结果可得知Pg Pl Poss (315 6.08)MW 321.08MW ;三台发电机的有功输出分别为、和。三台发电机之所以有这样的有功输出，是因为的前提假设是与发电 机1相连的母线为平衡节点，与发电机 2、3相连的节点为PV节点，这 意味着发电机节点的有功出力可变，而发电机 2和3的有功出力不变。在给定初始数据时，已经假定发电机2和3按照额

17、定的有功功率出力和， 由有功功率平衡可知，此时的发电机 1的有功出力只能为。造成无功功率不等的主要原因是变压器和输电线路会产生或消耗部分无功功率。由基态潮流计算结果很容易验证发电机发出的无功功率加 上输电线电纳产生的无功功率之和与线路上消耗的无功功率加上负荷消 耗的无功功率之和。当然，改变初始条件，如发电机出力及部分母线电压时，所算出的潮 流会和此种情况有所不同，但是它们都有相同规律，那就是系统的有功 功率和无功功率要平衡，且满足极限运行条件。最优潮流计算在运用matpower计算本次潮流计算之前，首先要清楚什么是最优潮 流。这里所谓的最优潮流的指标就是电力系统的经济运行，即在满足极 限运行条

18、件下，电力系统用最少的成本，发出能满足负荷需求的功率。假定有n个火电厂，其燃料消耗特性分别为Fi(Pgi),F2(FG2), ,Fn(PGn),系 统的总负荷为Rd，暂不考虑网络中的功率损耗，假定各个发电厂的输出 功率不受限制，则系统负荷在n个发电厂间的经济分配问题可以表示为： 在满足下式(5.1)nBi Pld 0 (3.1)i 1的条件下，使目标函数nF Fi(FGi) (3.2)i 1为最小。这是多元函数求条件极值的问题。由拉格朗日乘数法可知，当其满足 以下(5.3)条件时，在发出满足负荷功率需求的条件下，所消耗的燃料最 少，即经济性最好。坐 (i 1,2,., n) (3.3)dFGi

19、当然，任一发电厂的有功功率和无功功率都不应该超出它的上、下限, 各节点电压也必须维持在其允许的变化范围内，及最优潮流计算条件说明及数据输入和基态潮流计算一样，最优潮流计算时依旧选取 Busi为平衡节点，Bus2, Bus3为PV节点，其余的节点为 PQ节点；取发电机G2 G3的输 出有功功率为其额定功率。需要在此说明的是 gen erator cost data 是 最优潮流计算的必须条件。最优潮流的理论结果由中输入的数据可得到三机九节点系统中三台发电机特性分别为F1 0.11PG12 5PG1 150F2 0.085PG22 1.2PG2 600 (3.5)2F3 0.1225FG3 Pg3

20、 335其中有巳1 Pg2 巳3 315MW (3.6)结合(3.3)、(3.5)与(3.6)即可求得巴1 86.564MW,Pg2 134.377MW,Pg3 94.057MW (3.7)最优潮流仿真计算结果在Matpower中，有专门的函数计算系统的最优潮流。在命令行窗口 中输入runopf(case9hbyq) 指令即可得到基态潮流计算的结果。实际 仿真结果见附录二。表5-1优化后母线的潮流计算结果BusVoltageGen eratio nLoadMag(pu)An g(deg)P (MW)Q (MVAr)P (MW)Q (MVAr)1*-2-3-4-5-90306-7-8-9-Tot

21、al :表5-2优化后线路的潮流计算结果BranchFromBusToBusFrom Bus Injectio nTo Bus Injectio nLoss (IA2 * Z)P (MW)Q (MVAr)P (MW)Q (MVAr)P (MW)Q (MVAr)114245346457569678789827939Total最优潮流结果计算分析除了可以得到与基态潮流计算结果类似的分析外， 还可以得到以下分 析结果。线路损耗明显比基态潮流计算结果的小。由最优潮流计算结果可以得到，三台发电机的有功功率输出分别为、 和，与理论计算结果(3.7)比较接近。之所以会有些差别是因为输电线路 存在电阻，而电阻

22、上会消耗有功功率，因此实际功率平衡条件(3.6)并没 有得到充分的满足。因此，理论计算结果和实际计算结果接近，但并不 完全一致。但这并不意味着理论计算结果毫无意义，当用软件计算电力 系最优统潮流时，完全可以利用理论上的结果来判断潮流计算的结果是 否合理或正确。6.感想与小结本次潮流计算让我收获颇多，我不仅学会了利用 matpower 计算电力 系统潮流和最优潮流的方法，还懂得了与人多交流的重要性，更领悟到 了一种学术研究的方式。通过本次潮流计算， 我了解到了电力系统潮流计算及其意义， 在此基 础上，了解电力系统潮流计算的模型以及常规的潮流计算的方法，掌握 并熟练使用电力系统潮流计算软件 Mat

23、power。最后，利用Matpower计 算美国西部电网 WSCC三机九节点系统的静态潮流及最优潮流并给出分 析。达到了本次潮流计算的最初目的。除此之外， 我还懂得了与人交流的重要性。 在本次潮流计算之前， 我 完全不知道该如何使用 matpower,因此不得不慢慢地详读 matpower手 册。然而手册上的某些参数的说明并不完善，且无例子可参靠，因此学 习进程相当缓慢。后来在一次和师兄师姐的聊天中,知道了他们用过 matpower，于是便向他们请教，最终解决了疑问，快速地完成了潮流计 算。从这次潮流计算中,我学会了与人交流,与人分享,交流和分享将 会推动集体的进步，会让自己受益颇多。最后，我

24、还领悟到了一种学术研究的方式。在开始做潮流计算之前， 我便有个疑问，我们为什么要做潮流计算？潮流计算该如何算？因此， 围绕这个问题，我又开始了自己的思考，并小有所获。这种收获在我的 报告中有所体现，我的报告的书写逻辑很明确，首先介绍什么是潮流计 算，潮流计算有何意义，常用的计算方法有哪些，最后才过渡到利用 matpower计算电力系统的潮流这一块。这意味着在干一件事之前，我们首先要搞清楚做这件事有何意义，因为做一件没有意义的事情只会浪费我们的时间。对于一个课题，我们认识到它是一个问题甚至比我们该如何去做更重要。总而言之， 我将好好利用这一次的收获， 帮助我在未来科研的道路上 越走越远，越走越顺

25、。参考文献1何仰赞 , 温增银 . 电力系统分析 : 上册 M. 华中科技大学出版社 ,2002.2丁 晓 莺 , 王 锡 凡 . 最 优 潮 流 在 电 力 市 场 环 境 下 最 新 发 展 J. 电 力 系 统 自 动 化,2002,26(13):1-7.3徐恒娇 , 王洪诚 , 胡江航 , 等. 基于 Matpower 的潮流计算方法 J. 物联网技术 , 2013, 3(1): 43-45.4肖宇, 董珉 . 应用 Matlab 计算电力系统的最优潮流 J. 机械工程与自动化 , 2006 (4):60-62.附录一、基态潮流计算文件 runpf(case9hbyq)MATPOWER

26、 Version , 14-Dec-2011 - AC Power Flow (Newton)Newtons method power flow converged in 4 iterations.Converged in seconds| System Summary |How many? How much? P (MW) Q (MVAr)Buses 9 Total Gen Capacity to Generators 3 On-line Capacity to Committed Gens 3 Generation (actual)Loads 3 LoadFixed 3 FixedBran

27、chesDispatchable 0 Dispatchable of Shunts 0 Shunt (inj)9 Losses (IA2 * Z)Tran sformers 3 Branch Charg ing (inj) -In ter-ties 2 Total In ter-tie FlowAreas 2Minimum Maximum. bus 6deg bus 2 MW line 8-9 MVAr line 8-2Voltage Mag nitude . bus 9Voltage Angle deg bus 9P Losses (IA2*R) -Q Losses (L2*X) -| Bus Data |Bus Vol