牛拉法潮流计算程序(附3机9节点结果对比)VIP专享VIP免费

摘要电力系统潮流计算是研究电力系统稳态运行的一种重要方法,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态,包括各母线的电压、线路的功率分布以及功率损耗等等。潮流计算主要用于电网规划和静态安全分析，它可为扩建电力网络，以达到规划周期内所需要的输电能力提供依据；也可以对预想事故进行模拟和分析，校核预想事故下的电力系统安全性。本文简单介绍了牛顿-拉夫逊潮流计算的原理、模型与算法，然后用具体的实例，利用MATLAB对牛顿-拉夫逊法的算法进行了验证。关键词：电力系统潮流计算牛顿-拉夫逊法MATLAB一、牛拉法的数学模型对一个N节点的电力网路，列写节点电压方程，即˙I=Y˙V(1.1)式中，˙I为节点注入电流列相量，Y为节点导纳矩阵，˙V为节点电压列相量。由于异地测量的两个电流缺少时间同步信息，以注入功率替换注入电流作为已知量。即(1.2)其中,Yij=Gij+jBij，带入上式，得到有功功率和无功功率方程Pi=Vi∑j=1nVj(Gijcosθij+Bijsinθij)(1.3)Qi=Vi∑j=1nVj(Gijsinθij−Bijcosθij)(1.4)大部分情况下，已知PQ，求解Vθ。考虑到电网的功率平衡，至少选择一台发电机来平衡全网有功功率，即至少有一个平衡节点，常选择调频或出线较多的发电机作为平衡节点。具有无功补偿的母线能保持电压幅值恒定，这类节点可作为PV节点。潮流计算中节点分类总结如下：表1.1节点分类PQ节点负荷母线、变电站母线绝大部分是PQ节点已知PQ，求VθPV节点有无功储备的发电机母线很小部分是PV节点已知PV（V=1），求Qθ平衡节点选择一个发电机母线已知Vθ（V=1,θ=0）,求PQ电压相量作为参考轴θ=0已知电力系统有m个PQ节点，r个PV节点和1个平衡节点，则可以提取m+r个有功功率方程和m个无功功率方程，从而求解出m+r个θ和m个V，其余节点的有功和无功可通过式（1.3）、（1.4）求得，这样就完成了潮流计算。二、潮流计算流程用于潮流计算的有功和无功功率方程的阶数一般很高，常选择牛顿-拉夫逊迭代法进行求解高维的非线性方程组。牛顿-拉夫逊法的迭代过程如下：图2.1牛拉法迭代流程图将功率方程改成F(x)=0的形式ΔPi=Pis−Vi∑j=1nVj(Gijcosθij+Bijsinθij)=0(2.1)ΔQi=Qis−Vi∑j=1nVj(Gijsinθij−Bijcosθij)=0(2.2)可以得到[ΔθjΔVj]=−[∂ΔPi∂θj∂ΔPi∂Vj∂ΔQi∂θj∂ΔQi∂Vj]−1[ΔPiΔQi](2.3)X(k)ΔX(k)=−J−1(X(k))F(X(k))X(k+1)=X(k)+∆X(k)max{|ΔX(k)|}<ε？max{|F(X(k))|}<ε？结束为了便于计算，将上式改为[ΔθjΔVjVj]=−[∂ΔPi∂θj∂ΔPi∂VjVj∂ΔQi∂θj∂ΔQi∂VjVj]−1[ΔPiΔQi](2.4)求取功率方程的雅克比矩阵i≠j时，∂ΔPi∂θj=−ViVj(Gijsinθij−Bijcosθij)(2.5)∂ΔPi∂VjVj=−ViVj(Gijcosθij+Bijsinθij)(2.6)∂ΔQi∂θj=ViVj(Gijcosθij+Bijsinθij)(2.7)∂ΔQi∂VjVj=−ViVj(Gijsinθij−Bijcosθij)(2.8)i=j时，∂ΔPi∂θi=Vi2Bii+Qi(2.9)∂ΔPi∂ViVi=−Vi2Gii−Pi(2.10)∂ΔQi∂θi=Vi2Gii−Pi(2.11)∂ΔQi∂ViVi=Vi2Bii−Qi(2.12)牛拉法求解潮流的迭代过程如下图2.2牛拉法潮流迭代流程图三、MATLAB编程实例实例采用王锡凡主编的《现代电力系统分析》p326-327的例题。图3.1某小型电力网络图例题中，1号发电机为平衡节点，2号和3号发电机为PV节点，其余为PQ节点。为了便于编程，可将PV节点和平衡节点放在节点描述图的最后，故要重新对节点进行编号。。图3.2修改节点编号后的电力网络图支路表如下：%首节点末节点电阻电抗容纳之半变比120.0100.0850.0881130.0170.0920.0791240.0320.1610.1531360.0390.1700.1791450.00850.0720.07451560.01190.10080.104519100.0576017400.0625018600.058601节点表如下：%节点电压幅值电压相角节点有功节点无功节点类型/“1”为PQ节点，“2”为PV节点，“3”为平衡节点11.0000121.00-1.2500-0.5131.00-0.9-0.3141.0000151.00-1-0.35161.0000171.02501.63000281.02500.85000291.0400003程序清单：取ε=10-5a=load('zhilu.txt');b=load('jiedian.txt');N=size(a,1);%支路数Nbus=size(b,1);%节点数pq=0;fork=1:Nbus%PQ节点个数ifb(k,6)==1pq=pq+1;endendY=zeros(Nbus);%节点导纳矩阵fork=1:Nt...