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

1、华中科技大学现代电力系统分析潮流计算作业现代电力系统分析作业-基于Matpower的电力系统潮流计算专 业：班 级：姓 名：学 号：基于Matpower的电力系统潮流计算1.本次潮流计算的目的及意义本次潮流计算的目的及意义主要是了解电力系统潮流计算及其意义，在此基础上，了解电力系统潮流计算的模型以及常规的潮流计算的方法，掌握并熟练使用电力系统潮流计算软件Matpower。最后，利用Matpower计算美国西部电网WSCC三机九节点系统的静态潮流及最优潮流并给出分析。2.电力系统潮流计算及其意义电力系统潮流计算是研究电力系统稳态运行情况的基本电气计算，电力系统潮流计算的任务是根据给定的网络结构及

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

3、方法有很多种。最常见的方法则是牛顿拉夫逊法、高斯-塞德尔法、快速解耦法及直流潮流法。牛顿拉夫逊法能有效地解决非线性方程组，有较好的收敛性；它利用了导纳矩阵的对称性、稀疏性以及节点标号顺序优化的技巧，已经成为电力系统潮流计算的最常用方法。因此，本次潮流计算将重点了解牛顿拉夫逊法。牛顿拉夫逊法实质上就是切线法，是一种逐步线性化的方法。设有非线性方程组 其近似解为，.，。设近似解与精确解分别相差，则可得到 将上式中的多元函数在初值附近分别展开成泰勒级数。则可得到 忽略上式中，的二次及以上阶次的各项，整理便可以得到 上式是对于修正量，的线性方程组，称为牛顿法的修正方程式。利用高斯消去法或三角分解法可以

4、求出修正量，。然后对初始解近似解进行修正 经过k次迭代即可将上式写成如下形式 如此反复的迭代，直至满足收敛判据 或 将牛顿-拉夫逊法用于潮流计算，要求潮流方程写成形如方程式的形式，由于节点电压可以采用不同的坐标表示，牛顿拉夫逊潮流计算也将相应地采用不同的计算公式。3.2 节点电压用直角坐标表示时的牛顿拉夫逊潮流计算采用直角坐标系时，节点电压可表示为 导纳矩阵元素表示为 则可得到节点相应的有功功率和无功功率为 假定系统中的第1，2，m号节点为PQ节点；假定系统中的第m+1，m+2，n-1号节点为PV节点，第i个节点的给定功率设为和，对节点可列写如下方程 由于第n号节点为平衡节点，其电压是给定的，

5、故不参加迭代。式总共包含了2(n-1)个变量，而由可知，待求的变量也是2(n-1)个。还注意到已经具备了方程组的形式。因此即可得到如下的修正方程式 式中而J则为雅克比矩阵。 有了这些，即可按照上面的牛顿拉夫逊法计算电力系统的潮流。由于极坐标下的牛顿拉夫逊潮流计算的推导过程与上类似，这里就不在详尽介绍。3.3 牛顿拉夫逊法潮流计算程序框图牛顿拉夫逊法潮流计算程序框图如下图所示。图 3-1牛顿拉夫逊法潮流计算程序框图4.选用的潮流计算的系统本次选用的潮流计算系统为美国西部电网WSCC三机九节点系统。该系统比较经典，经常被业内人士用来作为分析的典型例子。其接线图如下图所示图 4-1三机九节点系统连接

6、图统的基准值和频率取为：系统中设备和线路的参数如表14所示：表4-1 发电机参数标幺值参数名称VN/kVxdxdxqxqxlTd0Tq0HDG116.50.1460.06080.09690.09690.03368.96s0s23.64s0G2180.89580.11980.86450.19690.05216s0.535s6.4s0G313.81.31250.18131.25780.250.07425.89s0.6s3.01s0表2 变压器参数的标幺值变压器名称T1T2T3电抗X(pu)0.05760.06250.0586变比K111表3 线路参数的标幺值 线路名称line1line2line3

7、line4line5line6阻抗Z(pu)R0.01000.03200.01700.03900.00850.0119L0.0850i0.1610i0.0920i0.1700i0.0720i0.1008i导纳B(pu)0.1760i0.3060i0.1580i0.3580i0.1490i0.2090i表4 负荷参数负荷名称LumpALumpBLumpC有功/MW12590100无功/Mvar5030355.利用软件matpower计算潮流5.1 matpower简介Matpower更准确地讲，应该是matlab中的一个工具箱，而不是一个独立的软件。它是由卡奈尔大学电气学院电力系统工程研究中心的

8、RAY D.ZIMMENRman，CARLOS E.Murillo和甘德强在ROBERT THOMAS的指导下开发出来的。最初的基于MATLAB的电力潮流和最优潮流代码是为POWERWEB项目的需要而编写的。 Matpower是一个基于matlab m文件的组建包，用来解决电力潮流和优化潮流的问题。它致力于为研究人员和教育从业者提供一种易于使用和可更新的仿真工具。Matpower的设计理念是用尽可能简单、易懂，可更新的代码来实现最优秀的功能。5.2 基态潮流计算5.21 基态潮流计算条件说明及数据输入 本次潮流计算时选取Bus1为平衡节点，Bus2，Bus3为PV节点，其余的节点为PQ节点；取

9、发电机G2、G3的输出有功功率为其额定功率，则依据matpower用户手册即可输入以下数据（数据说明见注释）：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节点)、注入的有功/无功、分区、%母线电压标幺值、基值、

10、相角、电压标幺值的上下限等bus = 1 3 0 0 0 0 2 1 0 16.5 1 1.1 0.9; 2 2 0 0 0 0 1 1 0 18 1 1.1 0.9; 3 2 0 0 0 0 1 1 0 13.8 1 1.1 0.9; 4 1 0 0 0 0 2 1 0 230 1 1.1 0.9; 5 1 90 30 0 0 2 1 0 230 1 1.1 0.9; 6 1 0 0 0 0 1 1 0 230 1 1.1 0.9; 7 1 100 35 0 0 1 1 0 230 1 1.1 0.9; 8 1 0 0 0 0 1 1 0 230 1 1.1 0.9; 9 1 125 50

11、0 0 2 1 0 230 1 1.1 0.9;%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% 电机参数，定义电机的有功无功出力，无功输出上下限、端口电压标幺值、电机基值等gen = 1 0 0 300 300 1 100 1 250 10 0 0 0 0 0 0 0 0 0 0 0;2 163.2 0 300 300 1 100 1 300 10 0 0 0 0 0 0 0 0 0 0 0;3 108

12、.8 0 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 angmaxbranch = 1 4 0 0.0576 0 250 250 250 1 0 1 -360 360; 4 5 0.017 0.092 0.158 250 250 250 0 0 1 -360 360; 5 6 0.039 0.17 0.358 150 150 150 0 0 1 -360 3

13、60; 3 6 0 0.0586 0 300 300 300 1 0 1 -360 360; 6 7 0.0119 0.1008 0.209 150 150 150 0 0 1 -360 360; 7 8 0.0085 0.072 0.149 250 250 250 0 0 1 -360 360; 8 2 0 0.0625 0 250 250 250 1 0 1 -360 360; 8 9 0.032 0.161 0.306 250 250 250 0 0 1 -360 360; 9 4 0.01 0.085 0.176 250 250 250 0 0 1 -360 360;%- OPF Data -% area data% area refbusareas = 1 5;%最优潮流计算所需的条件% generator cost data% 1 startup shutdown n x1 y1 . xn yn% 2 startup shutdown n c(n-1) . c0gencost = 2 1500 0 3 0.11 5 150; 2 2000 0 3 0.085 1.2 600; 2 3000 0 3 0.1225 1 335;5.22 潮流计算结果 将上述数据输入到case9hbyq.m文件中保存并运行后，在命令行窗口中输入runpf(case9hbyq)指令