电力系统潮流分析—基于牛拉法和保留非线性的随机潮流姓名：**＊学号:＊**1潮流算法简介1.1常规潮流计算常规的潮流计算是在确定的状态下.即：通过已知运行条件（比如节点功率或网络结构等）得到系统的运行状态（比如所有节点的电压值与相角、所有支路上的功率分布和损耗等）常规潮流算法中的一种普遍采用的方法是牛顿-拉夫逊法.当初始值和方程的精确解足够接近时，该方法可以在很短时间内收敛.下面简要介绍该方法1牛顿拉夫逊方法原理对于非线性代数方程组式（1T）,在待求量x初次的估计值x（0）附近，用泰勒级数（忽略二阶和以上的高阶项）表示它，可获得如式（1-2）的线性化变换后的方程组，该方程组被称为修正方程组f'（x）是f（X）对于X的一阶偏导数矩阵，这个矩阵便是重要的雅可比矩阵Jf（x,x,x）=0i=1,2,…,n（1-1）i12nf（x（0））+f'（x（0））Ax（0）=0（1—2）由修正方程式可求出经过第一次迭代之后的修正量Ax（0），并用修正量Ax（0）与估计值x（0）之和，表示修正后的估计值x（1），表示如下（1一4）.Ax(0)=-[f'(x(0))]-if(x(0))(1—3)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)V =e+fiiiY =G+jBijijij假设系统的网络中一共设有n个节点当采用直角坐标系解决潮流方程，此时待解电压和导纳如下式:（1-7）1-8)平衡节点的电压是已知的，平衡节点表示如下.V=e+jfnnn(完整)电力系统潮流计算方法分析除了平衡节点以外的所有2(n-1)个节点是需要求解的量。
每个节点可列出两个方程式•假定系统中前m个节点为P一Q节点，第m+1到n—1个节点为P-V节点对于PQ节点，P和Q的值是ii固定的，对于PV节点，P和V的值是固定的iiaP=P—e工(Ge—B/)—/工(G/+Be)=0,.iisijjjjjjjjj(1-9)\Q=Q-f工(Ge—Bf)+e乞(Gf+Be)=0i=1,2,…'mijjijjjijjijjj©i-Bf)-f工(Gf+Be)=0(110)可Ji二m+1,m+2,…,n一1iisi*j©aP=P—e工(Geiisiijj©iijjaV2二V2—(e2+/2)二oiisiiiij□选定电压初始值，按泰勒级数展开，忽略Ae,Af二次方程及以后各项，得到修正方程如下:iiAW=—JAU(1-11)其中：AW=[APAQ】•…APAU二[Ae1Af…AeAfmAQmAem+1APm+1Afm+1AU2m+1…Aen-1APn-1AfA，n-1AU2n-1-SAPbSe1SAQi-Se1SAPSf1SAQ1SAPSemSAQ1SemSAP-Sf+mSAQSfmSAPSem+1SAQSem+1SAPbSfm+1SAQi-Sfm+1SAP1Sen-1SAQSen-1SAPSfn-1SAQi-Sfn-1SAPSe1SAQSe1SPm+1Se1SAU2Se1SAPmSf1SAQSf1SPm+1Sf1SAU2m+1Sf1SAPmSAQSemSPm+1SemSAU2SemSAPmSfmSAQSAPSAPn~1Se1SAU2n~1Se1SAPSf1SAU2Sf1SAPn-1SemSAU2n—1Sem雅克比矩阵J各元素的计算公式如下:mSPm+1SfmSAU2m+1Sfmm1SAQSem1SPm+1Sem1SAU2m+1Sem+1SAPmSfm1SAQSf]m1SPm+1Sfm1SAU2m+1Sfm+1SAPn-1SAQSen-1SPm+1Sen-1SAU2Sen-1SAPmSfn-1SAQSfn-1SPm+1Sfn-1SAU2m+1Sfn-1SAPn~1SfmSAU2SfmSAPSem+1SAU2Sem+1SAPn—1Sfm1SAU2n—1Sfm1SAPSen-1SAU2■Sen-1SAPn-1Sfn-1SAU2Sfn-1（完整）电力系统潮流计算方法分析OAPOAQ丄)=-^77^=-(Ge+Bf)OAfijiijij:=Be-GfijiijijiOAeOAPOAQVi=iOAfOAejOAU2=OAU21-12)=0OfjOAPiOeiOAPiOfoaQ丄OeV1=-工(Ge-Bf)-Ge-Bfijjijjiiiiij=1=-工(Gf+Be)-Gf+Beijjijjiiiiiiij=1(Gf+Be)-Gf+Beijjijjiiiiiij=1OAQr=—工(GAe-Bf)+Ge+BfOfijjijjiiiiiijj=1OAU2i-=-2edei仝=-2fOfJii1—13)一般雅克比矩阵表示为：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)ijjijjiiiiiij©i=OAQ—iOej(Be-Gf)ijiiji(］丰i)=<、'乙(Gf+Be)+Be-Gf(j=i)ijjijjiiiiiij©iLij=OAQ—iOfjRij=6AU2=idej=OAU2=iSijOfjI(Ge+Bf)ijiiji(j丰i)|-E(Ge-Bf)+Ge+Bf(j=i)Ijjjjiiiiiiij©i=V0(jHi)=VI-2e(j=i)i=VI0(jHi)=VI-2f(j=i)i(1-14)牛顿拉夫逊方法求解框图如下：图1.1牛顿拉夫逊潮流计算法求解框图1.1。
2保留非线性法求解过程与牛顿法的不同之处在于，第一是假设雅克比矩阵在迭代过程中不变，即取初值9和U形成的雅克比矩阵来迭代；第二是计算出来的修正量一直是初始值的修正量由于保留非线性只对（完整）电力系统潮流计算方法分析直角坐标形式的公式不存在截断误差，因此为了减小计算误差，本文以直角坐标形式的牛拉法为基础编写了保留非线性潮流计算方法的程序迭代公式为：△x（k+i）=—J-i[y（x（o））—ys+y（Ax（k））]（1—14）迭代过程和牛拉法相类似，流程图如下所示：图12保留非线性法求解框图12蒙特卡罗模拟法11蒙特卡罗模拟原理（完整）电力系统潮流计算方法分析蒙特卡罗模拟方法的思想是，是当求解问题是一不确定事件的平均值时，我们通过构建模型并采用某特定的“实验”，就可以实验中此事件发生的频率去估算概率.1.22蒙特卡罗模拟步骤1）根据不同新能源的特点建立新能源输出功率的样本，规模为N；2）将得到的N个样本值带入对应接入新能源的各节点，得到接入光伏后的各节点的值3）按照11所述的牛顿拉夫逊法进行确定性潮流计算，得到N组关于节点的电压，支路功率与网损的数据等4）运用数学上的统计原理，可以求出输出变量的分布情况.1。
3拉丁超立方采样法11拉丁超立方采样原理拉丁超立方采样由M.D.McKay、RJ.Beckman和WConover在1979年提出，它通过分层采样使采样点能够覆盖到整个随机变量的分布范围该方法分成两步:1）采样：所有的输入变量可以通过分层采样，使得样本点更加准确均匀的分布；2）排列:改变初次采样得到的样本数据的顺序，令变量数据之间的关联程度最小,或者通过排序达到指定的相关系数1.3.2拉丁超立方采样优点1）可以使采样得到的数据较为全面地覆盖变量所分布的范围，同时分层使得采样时不会再采到一样或相似的数据，更准确地体现变量的总体情况，同时减小了样本规模一些文献证明了拉丁超立方采样与简单随机采样在采样规模同是M时，两种方法抽取到的变量假设是独立的，那么它们的联合覆盖空间百分比平均值表示如下：M-12P=xlOO%1—16)1LM_「M-112P=xlOO%m|_M+1_可以看出，当M大于等于2时，一式大于二式，表明拉丁超立方采样比随机采样覆盖的范围（完整）电力系统潮流计算方法分析a一rand、x=F-1(z)=Ft(),a=1,2,…,NkakaA(1-19)(1-20)大•比如当M=20时，按式(1T6)计算得：P=90.25%,P=81.86%。
lm2)拉丁超立方采样的稳健性好假设一输出随机变量Y满足下式Y=£cX(1—17)iii=1c是常数，丫是输入随机变量X的线性函数在相同采样规模下，进行一定次数的蒙特卡罗ii模拟，每一次都能获得一个关于丫的分布情况.由每个丫的分布的期望值可以得到一个新的分布用方差O表示这个分布的离散程度若O越大，表明不同仿真间的差异越大，算法的稳健性越ZZ不好文献指出通过拉丁超立方采样法得到的方差要比随机采样得到的方差小1/N2表明一共进行总数为N3的随机采样得到的方差a与只需进行N次拉丁超立方采样得到的方差a相同3.3拉丁超立方采样步骤1)采样假设X,X，…,X是随机潮流计算的N个输入变量X的累积概率分布是：12NkZ=F(X),k=1,2，•…N(1-18)kkk取采样规模为A,采样步骤为：a. 将Z的取值范围】0,1］均匀分为A等份，即［0,丄］，［丄,？］,...,［匕!,1］;kAAAAb.。