matlab函数runpf,基于MATLAB的电力系统潮流计算.pptx

基于MATLAB的电力系统潮流计算;目录;一、潮流计算的定义及应用;二、电力系统潮流计算方法;(2)确定计算方法。非线性代数方程组的解法离不开迭代，因此，潮流计算方首先要求他是能可靠的收敛，并给出正确答案。现在比较常用的方法有牛顿-拉夫逊法(简称牛顿法)和P-Q分解法。牛顿法是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。P-Q分解法是从改进和简化牛顿法潮流程序的基础上提出来的，他的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓住主要矛盾，以有功功率误差作为修正电压向量角度的依据。;(3)制定计算流程并编制程序。无论采用什么编程语言，其基本流程都是绘制程序流程图，根据选定的计算方法进行循环迭代从而得到最终的计算结果。(4)上机计算及对计算结果进行分析。;三、 MATPOWER的功能与使用方法;MATPOWER数据格式; branch变量是一个矩阵，用来设置电网中各支路参数，其格式为：fbus、tbus、r、x、b、rateA、mteB、mteC、ratio、aJlgle、status。fbus和tbus用来设置该支路由起始节点(母线)编号和终止节点(母线)编号。r、x和b用来设置该支路的电阻、电抗和充电电纳。rateA、 rateB和rate(；分别用来设置该支路长期、短期和紧急允许功率。ratio用来设置该支路的变比，如果支路元件仅仅是导线为O；如果支路元件为变压器，则该变比为fl，us侧母线的基准电压与thus侧母线的基准电压之比。angle用来设置支路的相位角度，如果支路元件为变压器，就是变压器的转角；如果支路元件不是变压器，相位角度为0度。status用来设置支路工作状态，1表示投入运行，0表示退出运行。 gen变量也是一个矩阵，用来设置接入电网中的发电机(电源)参数，其格式为：bus、魄、Qg、Qmax、Qmin、Vg、mBase、status、Pmax、Pmin。bus用来设置接人发电机(电源)的母线编号。Pg和Qg用来设置接人发电机(电源)的有功和无功功率。Pmax和Pmin用来设置接入发电机(电源)的有功功率最大、最小允许值。9max和Qmin用来设置接入发电机(电源)的无功功率最大、最小允许值。Vg用来设置接人发电机(电源)的工作电压。mBase用来设置接入发电机(电源)的功率基准，如为缺省值，就是base～lVA变量的值。status用来设置发电机(电源)工作状态，l表示投入运行，0表示退出运行。;MATPOWER电力潮流和最优潮流;默认的潮流计算方法是标准的潮流法[12]，采用全雅克比矩阵，迭代求解。这种方法在 许多文教科书中都有提到。法则 2和法则 3是改进型快速解耦算法[10].MATPOWER采用 XB和 BX变换，参见文献[1].法则 4是标准的高斯－赛德尔方法[5],基于意大利 Bologna大学的 Alberto Borhetti的贡献的代码开发。要使用出默认的牛顿法之外的其他算法，PF_ALG选项必须正确的设置。比如说，要使用 XB快速解耦算法： >> mpopt= mpoption('PF_ALG',2); >> runpf(casename,mpopt);最后一种算法是直流潮流算法[13],它的使用是通过设置 PF_DC为 1，运行 runpf而进行的，或者直接使用 rundcpf。直流潮流的计算是通过直接的不迭代的的方法解母线电压相角和指定母线的有功注入获得，基于方程 2，3和 4。对于交流潮流解法，如果 ENFORCE_Q_LIMS选项被设为 true(默认为 false)，并且运行过程中有任何发电机组的无功越限，相应的母线被转换为 PQ母线(节点)，将无功出力设定在限制值，并且案例重新计算。该母线的电压幅值为满足无功限制的要求将偏离指定值。如果参考母线(节点)的有功出力达到限制值，该节点将自动转化为 PQ母线(节点)，在下一轮迭代中第一个依然存在的 PV母线(节点)将被当作松弛母线(节点)，这将导致该母线(节点)的机组有功出力稍微偏离指定值通常，没有 MATPOWER的潮流解法中不包含变压器分接头的改变或者操作，或者部分系统从网络中解列等。 潮流计算的解法，除了高斯－赛德尔法之外，都可以很好的解决甚至是大规模网络，因为这些算法和计算充分利用了 MATLAB的内部稀疏矩阵处理;最优潮流MATPOWER提供多种解算最优潮流问题(OPF)的方法，可以通过访问函数 runopf的方法实现。除了提供将计算结果输出到屏幕之外(默认)，runopf函数还可以通过设置以下的参数返回解到其他地方。>>[bas