Matpower 学习总结和优化本人所编写的程序一、 Matpower 用法总结（潮流计算）应用Matpower计算潮流技巧的核心在于输入好三个矩阵和部分参数，清晰的知道输入 参数、矩阵中每一个元素的含义•下列以case（ ‘）为例子说明：★参数一%％ MATPOWER Case Format : Version 2mpc.version = '2'解释:目前普遍采用 2 形式的算法★ 参数二％% system MVA basempc.baseMVA = 100;解释：采用有铭值mpc.baseMVA = 100； （Matpower只能计算有铭值得网络）★ 矩阵一%% bus data％ bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 解释:bus data母线参数也就是我们所说的节点参数，下面逐条注释：1 bus number （positive integer） :第一列表示节点的编号（括号里面注释正整数）；2 bus type ：第二列表示节点的类型，一般只用得到1、 2、 3 三种节点类型,4 类型的节点目前没有接触到；PQ bus = 1PV bus = 2reference bus = 3isolated bus = 43 Pd, real power demand （MW） ：表示负荷所需要的有功功率（所有数据都是正数）;（有 铭值）4 Qd, reactive power demand （MVAr） :表示负荷所需要的无功功率（所有数据都是正数）； （有铭值）5 Gs, shunt conductance :表示和节点并联的电导，非线路上的电导，一般该列为0；6 Bs， shunt susceptance:表示和节点并联的电纳，非线路上的电纳，一般该列为0；7 area number, （positive integer） ：表示母线的断面号，一般设置为1；8 Vm, voltage magnitude （p.u。
：表示该节点电压的初始幅值（设置成标幺值）;9 Va， voltage angle （degrees） ：表示该节点电压的初始相位角度；10 baseKV， base voltage （kV） ：表示该节点的基准电压； （有铭值）11 zone, loss zone （positive integer） ：表示母线的省损耗区域，一般设置为1；12 maxVm, maximum voltage magnitude （pu.）:该节点所能接受的最大电压；（标幺值）13 minVm, minimum voltage magnitude （p u.） :该节点所能接受的最小电压；（标幺值）★ 矩阵二％％ generator data% 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 解释:表示发电机参数，下面逐条解释：1 bus number ：发电机节点的编号；2 Pg, real power output （MW） ：发电机节点输出的有功,如果为平衡节点则设置为0； （有铭值）3 Qg, reactive power output （MVAr ）：发电机节点输出无功，如果为平衡节点则设置为0； （有铭值）4 Qmax, maximum reactive power output （MVAr）:该节点能接受输出最大无功功率；（有 铭值）5 Qmin， minimum reactive power output （MVAr） ：该节点能接受输出最大有功功率； （有铭值）6 Vg， voltage magnitude setpoint（po u。
该节点电压的标幺值；7 mBase, total MVA base of this machine, defaults to baseMVA：该发电机的容量（有铭 值） ;8 status, > 0 — machine in service〈= 0 - machine out of service ：表示发电机的运行状态， 1表示投入， 0表示否9 Pmax， maximum real power output （MW） :允许输出的最大有功功率； （有铭值）10 Pmin， minimum real power output （MW） ：允许输出的最大无功功率；（有铭值）11 Pc1， lower real power output of PQ capability curve （ MW）12 Pc2, upper real power output of PQ capability curve （MW）13 Qc1min, minimum reactive power output at Pc1 （ MVAr）14 Qc1max， maximum reactive power output at Pc1 （ MVAr）15 Qc2min， minimum reactive power output at Pc2 （ MVAr）16 Qc2max, maximum reactive power output at Pc2 （MVAr）17 ramp rate for load following/AGC （MW/min）18 ramp rate for 10 minute reserves （MW）19 ramp rate for 30 minute reserves （MW）20 ramp rate for reactive power （2 sec timescale） （ MVAr/min）21 APF, area participation factor注释：红色区域参数均设置为0；矩阵%％ branch data％ fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax解释支路参数，下面逐条解释：1 f, from bus number :支路首端号；；2 t, to bus number: 支路末端号；3 r， resistance （p。
u ：支路电阻的标幺值；4 x, reactance （p.u ） ：支路电抗的标幺值；5 b， total line charging susceptance （p ） ：支路电纳的标幺值（注意:是整条支路的电 纳值）；6 rateA, MVA rating A （long term rating） :长距离输电支路所允许的容量（有铭值）；7 rateB， MVA rating B （short term rating） ：短距离输电支路所允许的容量（有铭值）；8 rateC, MVA rating C （emergency rating）:紧急输电支路所允许的容量（有铭值）；9 ratio， transformer off nominal turns ratio （ = 0 for lines ） （taps at 'from' bus， impedance at 'to' bus, i.e if r = x = 0, then ratio = Vf / Vt） ：支路变比,不含变压器设置为 0；含有变压 器变比为支路首端电压和末端电压之比:Matpower中变压器的模型，如下图所示：10 angle, transformer phase shift angle (degrees), positive =〉delay： 该参数设置为 0；11 initial branch status, 1 - in service, 0 一 out of service :该支路是否投入运行；12 minimum angle difference, angle(Vf)一 angle (Vt) (degrees):该支路所允许最小相位角 度13 maximum angle difference， angle(Vf) 一 angle (Vt) (degrees)： 该支路所允许最大相 位角度二、Matpower潮流计算结果和本人编写程序计算结果对比展示★各节点电压和相角| Bus Data |BusVoltageGun*crt ionLo^d辛JUg(pu? ATI?(deg)F (MV)Q (K^r)P (W 2! (JTVAr)1L000乩]?]L50. 00106.73D.1S0. 7320.9890,049--100. 0093. 6030.978-0-862--D. 007 9310.990-5. 941--100. 00SO. 00”0.931-6,978--20. 0015.0060.910-6.320--D. 070. OSt1.0007L-61102.55--lot al:22L-612D9.23220. 26L84. 43图一 matpower计算得出的结果输出各节点电压幅11 （与相角对应）E1.0000 0,3839 （1.S7S1 0,^04 O.S;11 0,9335 :・雜二各节齟曰眾；与嚨讨回昌俺5. 1709 0. 0492 -0.8622 -5. 9113 -6.9778 -6.3202 0图二本人程序计算得出的结果对比图一和图二可知，两种计算方法得出的节点电压和相位角度一致★ PV节点和平衡节点注入网络的功率輸岀各节瞬入的孵孤左往右檢节点黑腿)1.4991 + 1.05991 -L. DM0 - 0.93601 -0.0000 + 0.04981 -l.'JWfl - 0.80001 -3.2000 - 0.15001 -0.0008 - 140431 0. 7W1 + ：. 02551图三本人程序得出注入网络的功率Eus DataEus Val1 age ion LoadMag ipu)Angideg)P iW)Q（K臥t）F 'JW：'Q (mr)1].0005-171150.00106.730,190.7320. 9B90.血--100.0093.3030. 97S-D- B62--a.oo-du 9840. 9B0-5・ D41--10Q.Q0sa.50.93]-6- 970--20.00L5.Q060. 910-6- 320--a.a?0.087].000D. 000：＜7i.fi!102. S5--7Di al:2Z1.61209.20220.2618^.43图四matpower注入网络的功率对比图三和图四可知，两种计算方法得出PV节点和平衡节点注入网络的功率基本一致; 网络中的支路功率流动情况Bra。