(完整)电力系统潮流计算方法分析

1、电力系统潮流分析基于牛拉法和保留非线性的随机潮流姓名：*学号:*1潮流算法简介1.1常规潮流计算常规的潮流计算是在确定的状态下.即：通过已知运行条件（比如节点功率或网络结构等）得到系统的运行状态（比如所有节点的电压值与相角、所有支路上的功率分布和损耗等）。常规潮流算法中的一种普遍采用的方法是牛顿-拉夫逊法.当初始值和方程的精确解足够接近时，该方法可以在很短时间内收敛.下面简要介绍该方法。1.1。1牛顿拉夫逊方法原理对于非线性代数方程组式（1T）,在待求量x初次的估计值x（0）附近，用泰勒级数（忽略二阶和以上的高阶项）表示它，可获得如式（1-2）的线性化变换后的方程组，该方程组被称为修正方程组。

2、f（x）是f（X）对于X的一阶偏导数矩阵，这个矩阵便是重要的雅可比矩阵J。TOC o 1-5 h zf（x,x,x）=0i=1,2,n（1-1）i12nf（x（0）+f（x（0）Ax（0）=0（12）由修正方程式可求出经过第一次迭代之后的修正量Ax（0），并用修正量Ax（0）与估计值x（0）之和，表示修正后的估计值x（1），表示如下（1一4）.Ax(0)=-f(x(0)-if(x(0)(13)x(1)=x(0)+Ax(0)(1-4)重复上述步骤.第k次的迭代公式为：f(x(k)Ax(k)=-f(x(k)(1一5)x(k+1)=x(k)+Ax(k)(1-6)=e+fiii=G+jBijijij假

3、设系统的网络中一共设有n个节点当采用直角坐标系解决潮流方程，此时待解电压和导纳如下式:（1-7）1-8)平衡节点的电压是已知的，平衡节点表示如下.V=e+jfnnn(完整)电力系统潮流计算方法分析除了平衡节点以外的所有2(n-1)个节点是需要求解的量。每个节点可列出两个方程式假定系统中前m个节点为P一Q节点，第m+1到n1个节点为P-V节点。对于PQ节点，P和Q的值是ii固定的，对于PV节点，P和V的值是固定的。iiaP=Pe工(GeB/)/工(G/+Be)=0,.iisijjjjjjjjj(1-9)Q=Q-f工(GeBf)+e乞(Gf+Be)=0i=1,2,mijjijjjijjijjji-

4、Bf)-f工(Gf+Be)=0(110)可Ji二m+1,m+2,n一1iisi*jaP=Pe工(GeiisiijjiijjaV2二V2(e2+/2)二oiisiiiij选定电压初始值，按泰勒级数展开，忽略Ae,Af二次方程及以后各项，得到修正方程如下:iiAW=JAU(1-11)其中：AW=APAQ】APAU二Ae1AfAeAfmAQmAem+1APm+1Afm+1AU2m+1Aen-1APn-1AfA，n-1AU2n-1-SAPbSe1SAQi-Se1SAPSf1SAQ1SAPSemSAQ1SemSAP-Sf+mSAQSfmSAPSem+1SAQSem+1SAPbSfm+1SAQi-Sfm+

5、1SAP1Sen-1SAQSen-1SAPSfn-1SAQi-Sfn-1SAPSe1SAQSe1SPm+1Se1SAU2Se1SAPmSf1SAQSf1SPm+1Sf1SAU2m+1Sf1SAPmSAQSemSPm+1SemSAU2SemSAPmSfmSAQSAPSAPn1Se1SAU2n1Se1SAPSf1SAU2Sf1SAPn-1SemSAU2n1Sem雅克比矩阵J各元素的计算公式如下:mSPm+1SfmSAU2m+1Sfmm1SAQSem1SPm+1Sem1SAU2m+1Sem+1SAPmSfm1SAQSfm1SPm+1Sfm1SAU2m+1Sfm+1SAPn-1SAQSen-1SPm+

6、1Sen-1SAU2Sen-1SAPmSfn-1SAQSfn-1SPm+1Sfn-1SAU2m+1Sfn-1SAPn1SfmSAU2SfmSAPSem+1SAU2Sem+1SAPn1Sfm1SAU2n1Sfm1SAPSen-1SAU2Sen-1SAPn-1Sfn-1SAU2Sfn-1（完整）电力系统潮流计算方法分析OAPOAQ丄)=-77=-(Ge+Bf)OAfijiijij:=Be-GfijiijijiOAeOAPOAQVi=iOAfOAejOAU2=OAU21-12)=0OfjOAPiOeiOAPiOfoaQ丄OeV1=-工(Ge-Bf)-Ge-Bfijjijjiiiiij=1=-工(Gf

7、+Be)-Gf+Beijjijjiiiiiiij=1(Gf+Be)-Gf+Beijjijjiiiiiij=1OAQr=工(GAe-Bf)+Ge+BfOfijjijjiiiiiijj=1OAU2i-=-2edei仝=-2fOfJii113)一般雅克比矩阵表示为：Hij=OAP=i=Oej-(Ge+Bf)ijiiji工(Ge-Bf)-Ge-Bf(j=i)j(j幻)ijjijjiiiiiiNijOAP=i=Ofj(Be-Gf)ijiiji(j幻)工(Gf+Be)+Be-Gf(j=i)ijjijjiiiiiiji=OAQiOej(Be-Gf)ijiiji(丰i)=、乙(Gf+Be)+Be-Gf(j=i

8、)ijjijjiiiiiijiLij=OAQiOfjRij=6AU2=idej=OAU2=iSijOfjI(Ge+Bf)ijiiji(j丰i)|-E(Ge-Bf)+Ge+Bf(j=i)Ijjjjiiiiiiiji=V0(jHi)=VI-2e(j=i)i=VI0(jHi)=VI-2f(j=i)i(1-14)牛顿拉夫逊方法求解框图如下：图1.1牛顿拉夫逊潮流计算法求解框图1.1。2保留非线性法求解过程与牛顿法的不同之处在于，第一是假设雅克比矩阵在迭代过程中不变，即取初值9和U形成的雅克比矩阵来迭代；第二是计算出来的修正量一直是初始值的修正量。由于保留非线性只对（完整）电力系统潮流计算方法分析直角坐

9、标形式的公式不存在截断误差，因此为了减小计算误差，本文以直角坐标形式的牛拉法为基础编写了保留非线性潮流计算方法的程序。迭代公式为：x（k+i）=J-iy（x（o）ys+y（Ax（k）（114）迭代过程和牛拉法相类似，流程图如下所示：图1。2保留非线性法求解框图1。2蒙特卡罗模拟法1。2。1蒙特卡罗模拟原理（完整）电力系统潮流计算方法分析蒙特卡罗模拟方法的思想是，是当求解问题是一不确定事件的平均值时，我们通过构建模型并采用某特定的“实验”，就可以实验中此事件发生的频率去估算概率.1.2。2蒙特卡罗模拟步骤1）根据不同新能源的特点建立新能源输出功率的样本，规模为N；2）将得到的N个样本值带入对应接

10、入新能源的各节点，得到接入光伏后的各节点的值。3）按照1。1所述的牛顿拉夫逊法进行确定性潮流计算，得到N组关于节点的电压，支路功率与网损的数据等。4）运用数学上的统计原理，可以求出输出变量的分布情况.1。3拉丁超立方采样法1。3。1拉丁超立方采样原理拉丁超立方采样由M.D.McKay、R。J.Beckman和W。J。Conover在1979年提出，它通过分层采样使采样点能够覆盖到整个随机变量的分布范围。该方法分成两步:1）采样：所有的输入变量可以通过分层采样，使得样本点更加准确均匀的分布；2）排列:改变初次采样得到的样本数据的顺序，令变量数据之间的关联程度最小,或者通过排序达到指定的相关系数。

11、1.3.2拉丁超立方采样优点1）可以使采样得到的数据较为全面地覆盖变量所分布的范围，同时分层使得采样时不会再采到一样或相似的数据，更准确地体现变量的总体情况，同时减小了样本规模。一些文献证明了拉丁超立方采样与简单随机采样在采样规模同是M时，两种方法抽取到的变量假设是独立的，那么它们的联合覆盖空间百分比平均值表示如下：M-12P=xlOO%116)1LM_M-112P=xlOO%m|_M+1_可以看出，当M大于等于2时，一式大于二式，表明拉丁超立方采样比随机采样覆盖的范围（完整）电力系统潮流计算方法分析(完整)电力系统潮流计算方法分析a一rand、x=F-1(z)=Ft(),a=1,2,Nkak

12、aA(1-19)(1-20)大比如当M=20时，按式(1T6)计算得：P=90.25%,P=81.86%。lm2)拉丁超立方采样的稳健性好。假设一输出随机变量Y满足下式Y=cX(117)iii=1c是常数，丫是输入随机变量X的线性函数。在相同采样规模下，进行一定次数的蒙特卡罗ii模拟，每一次都能获得一个关于丫的分布情况.由每个丫的分布的期望值可以得到一个新的分布。用方差O表示这个分布的离散程度。若O越大，表明不同仿真间的差异越大，算法的稳健性越ZZ不好。文献指出通过拉丁超立方采样法得到的方差要比随机采样得到的方差小1/N2。表明一共进行总数为N3的随机采样得到的方差a与只需进行N次拉丁超立方采

13、样得到的方差a相同。ZZ1。3.3拉丁超立方采样步骤1)采样假设X,X，,X是随机潮流计算的N个输入变量。X的累积概率分布是：12NkZ=F(X),k=1,2，N(1-18)kkk取采样规模为A,采样步骤为：将Z的取值范围】0,1均匀分为A等份，即0,丄，丄,？,.,匕!,1;kAAAA从所有区间内依次抽取一个值作为一个采样值，区间内的抽取是随机的；c。由累积概率分布Zk的反函数变换后，便能得到输入变量Xk的样本数据.第a个区间z的采样值和X的第n个采样值如下:kka一randZ=,a=1,2,NkaA总共有N个输入变量，每个随机变量采样规模为A,假设将随机变量的数据以行为单位依次排列，那么最

14、终可以得到N*A阶的样本矩阵2）排序在求解随机潮流时，往往假设输入随机变量是独立的，但是按照上述方法得到的样本矩阵具有一定的相关性。我们需要分析和处理样本矩阵的关联性。使得变量数据值之间的关联性最小或者通过排序达到指定的相关系数。2系统模型建立光伏接入后的配电网系统主要由光伏发电系统、负荷和发电机三部分组成。太阳能光伏发电利用光伏电池可将光照转变为电动势的原理。在研究光伏并网后的随机潮流计算等有关问题时,首先要确定的是光伏发电的输出功率的随机特性,而此出力与太阳的光照强度密切相关，所以要想得到出力情况，必须先求出光照强度的随机分布30-3。4本次光伏发电，采用的是典型的Beta分布。此时我们可

15、以得到光照强度的概率密度函数为：r（a）r（p）Is、1max1S丿max(2-1)其中S是指光照强度统计时间内的实际值，s是指最大值.r是Gamma函数.a和p是形状max参数，将一段时间里太阳光照强度的期望值卩和方差5进行下式的变换便能得到形状参数35-36。卩（1-卩）1p=(1卩)卩（1-卩）(2-2)(2-3)假设光伏发电所用的电池方阵中有N个电池组，每个电池组的面积为A，光电转换效率为n耳n二1,2,N。n那么电池方阵总体的光电之间转化效率耳和方阵总的面积A分别是：A耳nn2-4)(2-5)门=-n=1AA亠nn=l此时这个电池方阵总的输出功率为：P二耳SA（26）N通过（24）-

16、（26）,在光照强度的概率密度函数基础上,便能推导出光伏输出功率的概率密度函数为：f(P)=r(a+p)、a1了DAp1r（a）r（p）IP丿ma1P丿max(2-7)其中，P二耳AS，为光伏出力的最大值.mama当a=0.8，p=2时，光照强度的概率分布曲线为:概率密度函数（IIIIIP0,10.20J04CL50.6OJ0,8C191S/Smax图2。1形状参数为0。8和2时光照强度的概率分布图配电网中可以将接入光伏的节点视为PQ节点，主要由于通过调节电容器可以使得功率因数恒定.3IEEE30节点算例3.1IEEE30节点系统介绍IEEE-30节点系统包括6台发电机，30个节点与41条支路

17、。选取系统的主要接线图如下:图3.1IEEE30节点系统接线图在计算时,为了简化计算对节点进行了重新编号。3.2两种常规潮流算法比较分别采用牛顿拉夫逊法和保留非线性法对IEEE30节点进行潮流计算，选取精度为108。牛拉法的迭代次数为6次，时间为0o031021s;保留非线性的迭代次数为12次，时间为0.022598s。保留非线性的迭代次数多但是总的计算速度快。牛拉法则是相反。以30个节点的电压为例，误差表示两值之差，计算的结果如表3.1所示。表3.1两种常规潮流算法对比保留非线性电压幅值/标幺牛拉、口误差相角/弧度保留非线性牛拉、口误差1o02991.02890.0010.097829-0.

18、09749-0。000341.02621.02470o0015-0o11691-0.116430.000481.02311.02030o0028-0o134780.13448-0。00031.01221.00550.0067-0.163360.16238-0.000981.03661.04390.00730o16185-0.163280。001431o02191.0408-0.0189-0o194560o196110.001551o04471o0460o0013-0.184610o181180.003431.02841o0321-0o0037-0.200440.197270.003171.02

19、251o0285-0.0060.201450.19954-0.001911o02781.03650.00870o19360o19242-0.001181.01821o0342-0.0160o19801-0o19870。000691.00991o0206-0o0107-0o21159-0.21063-0。000961o00561o01910o01350.214160o213830。00033(完整)电力系统潮流计算方法分析(a)保留非线性1.00891。02380.0149-0。210280.210420.000141。00981。0287-0。01890.20286-0。204310.0014

20、51.01051.0293-0。01880。202670.20420.001531。00951.0207-0。0112-0。208020。20776-0。000261.00081.0186-0.01780。21059-0.212590.0021。01141。02190.01050.212250。211630。000620.993611。0043-0。010690.219640。218860.000781.02681.03260。0058-0.20837-0。2064-0.001971。01731.01470.00260.14329-0。142930.000361。00711.0129-0.00

21、580。229690.22748-0.002210。995631。0015-0。005870。24498-0.24259-0。002391。03521。0340。00120.061255-0。060770。000491.01821。0060。01220.18095-0。17858-0。002371。02961.0230。00660。137560.136150。001411。09761。0910。00660.129120.130580.001461。09861.0880。01060。16399-0.160390。00361.051.050000在相同节点接入了相同的光伏发电，样本规模为500，采

22、用蒙特卡罗模拟法得到节点1电压的PDF与CDF如图3。1和3.2所示。可以看出两种算法还是存在差异的。(完整)电力系统潮流计算方法分析1(a)保留非线性000050211FDP-N5062OQ11.031.0321.0341.0361.0381.041.042U(b)牛顿拉夫逊图3.2两种算法下电压1的PDF图U1860000FDC-B（完整）电力系统潮流计算方法分析642000FDC-N1%25104TU（b）牛顿拉夫逊图3。3两种算法下电压1的CDF图3.3两种随机潮流算法的比较将以简单随机采样为基础的蒙特卡罗模拟法（MCSRS）和以拉丁超立方采样为基础的模拟法（MCLHS）得出的数据从准

23、确性和性能等方面做一个评估，全面比较两种随机潮流算法。3。3。1模型的准确性评估通过对输入随机变量的概率分布参数拟合，来分析所建立的模型的有效性和正确性。拟合的效果用相对误差指标来表示，表明分布情况的参数x的相对误差指标计算公式如下:x100%3-1)ccxfxbcxbc和c分别为参数x的样本拟合值和给定值.xfxb对光伏的输出功率采用Beta分布模型进行评估。Beta分布的两个形状参数的选取值为：a=0.9,卩二0.85.在一定规模下，根据光伏采样样本得到样本的平均值和方差，得到形状参数a,P的拟合值。并根据式（3-1）与实际的给定值0.9、0.85相比较得到误差.不同规模下分别采样50次后

24、，将平均值作为最终的相对误差指标来评估分布模型的准确性，以减小随机性对结果产生表3。2光伏形状参数相对误差指标对比表采样规模NMCSRSMCLHSaBaB1005.16582。88171。55131.54703004.54313.92500.52540。52676001。07972.10750.26200.260610000.89950。86650。15750。157830000.81560。58840。05240。052460000。36860。50690.02630。0263100000.21180.11610.01570.0157300000.19410。36040。00520.0052

25、的影响。由表可以看出，相同规模下，MCLHS比MCSRS的误差更小，用MCLHS生成的样本准确性更高。随着规模的增加，MCLHS和MCSRS生成的样本数据的正确性都有很大的提高。3。3.2性能评估通过算出的输出变量的平均值与标准差去评估MCLHS与MCSRS两种方法的计算精确度。计算公式如下：Euxcucuxfxbcuxbx100%Edxcdcdxfxbcdxbx100%32)33)上面两个式子式分别用来表示平均值与标准差的相对误差指标。采样规模为N时，一类输出变量便有N个数值，输出变量相对误差指标用这N个值的期望值表示.X分为mean、std、max和min四类。为减小随机性对结果产生的影响

26、，对两种方法在不同规模下分别采样50次，最后输出变量误差指标用50次误差的平均值mean表示，将这50次误差计算的标准差std、最大值max与最小值min用来评估上述方法收敛性与稳健性。cu和cd是误差计算的参考值.分别选取用xbxb20000次蒙特卡罗模拟得到的所选取的电压、功率和网损值来作为参考值。本次算例以节点18电压值、支路编号为3（3-4）的功率值与网损值作为研究对象.1）选取采样规模为500，以节点18电压值，支路3的功率值与网损值为研究对象，将得到的（完整）电力系统潮流计算方法分析平均值和标准差与参考值比较得到误差。两种方法均在此规模下进行50次仿真,得到50次计算结果的平均值、

27、标准差、最大值和最小值（单位）。表3。3两种方法在采样规模为500时的误差比较表仿真方法电压平均电压标准差平均值标准差最大值最小值平均值标准差最大值最小值MCLHS0。00310.00000。00310。00312。05380。13732.29921.6746MCSRS0。01550.00840.03190。00087。87095。076418.46160。1886仿真功率平均功率标准差方法平均值标准差最大值最小值平均值标准差最大值最小值MCLHS0.11590.00110。11880。11382。50380.23082.94651.9608MCSRS1。04450。59012.41200.03871