4.第三章电力系统运行的灵敏度分析及应用

1、第三章电力系统运行的灵敏度分析及应用第一节灵敏度分析分析在给定的电力系统运行状态下，某些量发生变化时，会引起其他变量发生多大变化的问题。这一问题当然可通过潮流计算来解决，但计算工作量大。采用灵敏度分析法，计算量小，并可揭示各量之间的关系。但变化量大时，灵敏度分析法的精度不能保证。一、灵敏度分析的基本方法1、常规计算方法电力系统稳态运行的潮流方程一般性描述为：,f区u)二0(3-1)y二y(x,u)x为状态变量，如节点电压和相角；u为控制变量，如发电机输出功率或电压；y为依从变量，如线路上的功率。实际上，(3-1)中f(x,u)=0就是节点功率约束方程，y=y(x,U)是支路功率与节点电压的关系

2、式。设系统稳态运行点为(x,U)，受到扰动后系统的稳态运行点变为00(x+x,u+u)。为了求出控制量变化量与状态量变化量之间的关系，在(x,U)0000处将(3-1)按泰勒展开并取一次项，得：afaf3-2)f(x+x,u+u)=f(x,u)+x+u=00000axauy+y=y(x,u)+学x+学u000axau代入，有：f(x,u)=000y=y(x,u)000afaxu=0auau3-3)j“i-1OfuSuOuxuOyyx+uOxOux堆Oy3-4)空S+空uSu(OxxuOu丿yu其中S-xu,OfOfSyu(Ox丿OuOls(OxxuOu丿3-5)为u的变化量分别引起x和y变化量

3、的灵敏度矩阵。如果控制变量为各节点的有功、无功设定量，则竺diag11OuS就是潮流方程的雅可比矩阵的逆。u,x为两个不同状态间的变化量。xu.1，所以,2、准稳态灵敏度计算方法考虑到电力系统运行的实际：1)初始控制变量的改变量，与到达新稳态的最终改变量不同；2)一个控制量的变化可能使另一些控制量也发生变化。所以控制变量的初始改变量与最终改变量不同，表示为：uFu(0)u3-6)由此得到准稳态的灵敏度关系：xSuSFu(0)SRu(0)xuxuuxuIySuSFu(o)Sru(o)yuyuuyu第二节潮流灵敏度矩阵3-7)1、发电机母线电压改变量V与负荷母线电压改变量V之间的灵敏度关GD节点注

4、入无功的平衡量方程Q-VV(Gsin”-Bcos”)iijijijijijj“iQ+VB0ijij3-8)上式简化依据了电力系统结构和运行的特点。根据灵敏度分析的基本方法BB一,V一r，QDDBGDDGBGGD,VGD，QG将（3-8）在当前状态点泰勒展开舍去高次项，的受到扰动后各变量变化量之间的关系,Q+B,V0iijj,Q-B,Viijj写成矩阵形式，jwi并将负荷节点与发电机节点分开排列3-9)3-9）式与P-Q分解法V-Q迭代的修正方程式形式一致。但要注意在这里,Q、,Q是发电机和负荷的变化量。即（3-9）式表示了系统新稳态相对于旧DG稳态控制量的变化量与状态量的变化量之间的关系。假定

5、v调整后，负荷的无功功率不变化，即,Q0，则式（3-9）第一式为:GDB,V+B,V0DDDDGG变换得,V-B-1B,VS,V（3-10）DDDDGGDGG其中S-B-1B（3-11）DGDDDG为,v与,v之间的灵敏度矩阵。通过灵敏度矩阵可以知道哪些发电机对控制负DG荷母线电压最有效，从而实现对负荷电压的定量控制。几种情况讨论：（1）只调整部分发电机的电压，无功充足能维持电压不变（,v=0）的发G电机对（3-9）式没贡献，可从B中划去发电机电压能维持不变的节点对应的DG列。（2）被控量为部分负荷节点，即其它负荷节点的电压不关心，可从B、BDDDG中高斯消去不关心电压变化的负荷节点。3)无功

6、已达界的发电机，不能作为控制变量，也不能维持节点电压不变，高斯消去这些发电机的节点。这些节点的AQ0。G高斯消去是等值变换，直接划去是不考虑它的影响2、发电机母线电压改变量AV，负荷母线电压改变量AV与发电机输出无G功的改变量AQ之间的灵敏度关系G将(3-9)变换为_AV_BB-1AQRRAQD-DDDGDDDDGDAVGBGDBGGAQGRGDRGGAQG负荷的无功功率不变，即AQD3-12)假定发电输出无功改变时，0，有3-13)3-14)AVRAQDDGGAVRAQGGGGR、R是灵敏度矩阵。DGGG几种情况讨论：(1)不是控制变量的PV节点，其电压可维持不变，可直接划去对应的行和列。(

7、2)不是控制变量的PQ节点，输出无功不变，当电压会发生变化，可将对应节点高斯消去。3)不关心的负荷节点，直接划去。3、负荷母线电压改变量AV与变压器变比改变量At之间的灵敏度关系D将节点无功平衡方程重写如下Q-VV(Gsin0-Bcos0)iijijijijijuQ+VB0ijij其中B是变压器变比的函数，不考虑节点注入无功的变化，将变压器变比ij作为控制变量，节点电压作为被控变量，写出灵敏度方程3-15),B乙BAV+乙V一jAt0ijjj,t1jijiIj上式中t为之路i,j的变压器变比。写成矩阵形式，包括所有负荷节点，并ij假定发电机母线电压不变，即认为发电机无功充足，可维持电压不变。A

8、t03-16),BBAV+1jVDI,tjijAV_IbliD,BjV,tjijAt3-17)为（3-15）式中第二项所组成的矩阵，行对应B仅包含负荷节点。|竺jVL,t.jij负荷节点，列对应可调变压器支路。每列中只有两个非零元素，分别在变压器支路的两个端点上。如果变压器支路有一个端点为PV节点，则由于PV节点电压不变，所以对应该变压器的支路只有一个非零元素。第三节分布因子分析节点注入有功功率变化、支路开断（结构变化）与支路潮流变化的灵敏度。1、支路开断分布因子分布因子：支路/基态有功潮流为p，支路/开断引起支路k功率变化量为lAPi，两者之间的关系表示为：kAPiDp（3-18）kk-ll

9、D为分布因子。k-l相似与无功平衡方程，由有功平衡方程可得节点有功注入变化量与节点电压相角变化量之间的灵敏度方程P二B0或0二XP（3-19）0B是以1为支路参数建立的导纳矩阵，X是B的逆。0 x0iPjPn考虑一条支路l（i,j）断开的情况。如图，假定支路开断不引起节点注入功率的变化，则支路开断后，新网络节点的注入功率变化量为P=b.P.-P.o（3-20）ll其中，节点i的改变量p（pp）二p，节点j的改变量p-（p+p）=-pO（3-20）iilljjll可表示为：P=b.1.-1.olp=Mp（3-21）lllM是节点-支路关联列矢量，行对应节点号，支路l离开节点元素为1，进入l节点元

10、素为-1，节点与支路无关元素为0o新网络的导纳矩阵变为B-Mx-1Mt，开端后节点电压相角的变化量由TOC o 1-5 h z0lll3-19）得0二(B-Mx-1Mt)-1P(3-22)0lll利用矩阵求逆辅助定理(B一Mx-1Mt)-1=B-1一B-1M(MtB-1Mx)-1MtB-10lll00ll0lll0二X-XM(MTXM-x)-1MTX(3-23)lllll=x一ncnTlll其中n=xmllc=(MtXM一x)-1=(X一x)-1lllll-llX=MTXM=lo.1.-1.oxlo.1.-1.0-Tl-lll=lx-x.x-x.x-x.x-x-b.1.-1i1j1iijiij

11、jjinjn=X-X-X+Xiijiijjj=X+X-2XiijjijX为在原网络支路l两端节点i注入单位电流，节点j流出单位电流，其它ll节点注入电流为0的情况下，节点i与节点j的点位差，定义为端口ij的自阻抗。支路l开断后，支路k(m,n)上有功潮流的变化量“mt,0mt(xncnr)M厂小厂,Pi=k=kiii厶P=DPkxxlkllkk(3-24)M为支路k的节点-支路关联矢量。支路k与支路l之间的支路开断分布因子kMT(XncnT)MD=kiii厶klXkMTXMMTXMMTXM/(XX)=k1k111iiiXkXXX/(XX)=kk_iii1Xk3-25)其中X/x=klk1X/X

12、lllohb.1.1.0丄X=MtXM=t).1.1.klkl=lxX.XX.XXm1n1minimjnj=X+XXXminjnimj为在原网络支路1两端节点i注入单位电流，节点j流出单位电流，其它节点.1.1.0T注入电流为0的情况下，支路k两端节点m与节点n的点位差，定义为端口ij与端口m-n之间的互阻抗。若支路1开断后网络分解为互不连通的两部分，因这时x=X使(3-25)无lll定义。推导：当有多条支路同时开断时，支路k的功率变化量,p=Edpkkii1=1m为断开的支路数。ieG2、发电机输出功率转移分布因子发电机输出功率转移分布因子定义为发电机输出功率变化引起支路潮流的变化量，表示为

13、：Pi二GP(3-26)kk，iiG为发电机输出功率转移分布因子。i为发电机号，k为支路号。k，i假定发电机i输出功率变化后引起的功率不平衡完全由平衡节点吸收，其它节点的输出功率不变化，则节点电压相角的变化量：3-27)3-28)A0=xb.AP.ol=XAPiiiX是阻抗矩阵X的第i个列矢量，X是直流潮流中B矩阵的逆。io支路k(两端节点号分别为m和n)上有功潮流的变化5MTA0MTXPXX5Pi=k=kiJ=minAPkxxxiTOC o 1-5 h zkkkAPi=GAPkk，iiX-XG=minik，ixkG为转移分布因子,X、X是X中的m行i列和n行i列的元素。k，imini3、准稳

14、态发电机输出功率转移分布因子准稳态发电机输出功率转移分布因子设n台发电机有功出力调整量为AP，如果调整量之和不为零则产生功率GG不平衡，不平衡量为Yap=1apGGGiieG如果大于零则为功率超额，否则为功率缺额。实际电网中功率不平衡由所有发电机按一定比例承担，n台发电机的承担系数矢量为a，并且GG工a=1ta=1a0(3-29)iGGi各发电机的实际调整量P=P(o)-Ga工P(o)=P(o)，aItP(o)=(I-aIt)P=GGGGGGGGGGGiiGFP(o)uGK=FP（o）uGF=I一a1tuGGG的方阵，是准稳态响应的变换矩阵。3-30)将（3-28）推广到考虑多台发电机的情况，

15、为P=GPkk，iiVi=1X-XG=minik-Xk写成矩阵的形式P=GPkk，GGV1G=MtXek，Gxk3-31)kGM:支路k的关联矢量，X:k包含所有节点在内的阻抗矩阵，单位矩阵，每列都是一单位矢量，只在相应的发电机节点处有非零元1。G为k-G一行矢量。G=k，GXXminixkXXmjnjXkXXmlnlxk下。将（3-30）代入（3-31）=GrP(o)kG=GlPkk，GG）GR=GFk，Gk，Gu为了使支路k的有功潮流变化P，各发电机有功的调整量可由伪逆计算如k3-32)P=(Gr)tGr(Gr)t-iPGkGkGkGk3-33)由P反求AP（o）时，（3-32）式是一个不

16、定方程，可由伪逆法得到一个解。定kG3-34)当以上i，j)节点对多于一个时，定义a是上调节点集,P是下调节点集。义A是mn阶矩阵，At(AAt)-1是A伪逆矩阵。功率传输转移分布因子节点之间传输功率变化时引起的支路潮流的变化量。电力市场节点间购销合同的变化。当节点(i,j)之间的传输功率变化了,P时，节点i，ijj注入功率的变化量分别为+,P，-,P。由于没有产生功TOC o 1-5 h zijij率不平衡，可用常规法计算。任一支路k(m,n)上功率的变化量：,Pij,Pi+,PjkkkX-X人XXmi*,Pm/j,PxijxijkkXXX+Xmimjnini,PxijkAPyG,Pkkij

17、ijXXX+XGjmimjninjkxk当上调量为,P+，下调量为,P，且上调总量与下调总量相等，支路k(m,n)上潮a卩流的变化量,Pa卩。k曲ap丄(”(MTX,P)-”(MTX曲)kxkiikjj3-35)ki“aj“aMtX(e,P+e,P-)PXkaa卩卩kxkk-ak-p,P-p厂XMTXeXX.XXkakam1n1miniXMTXeXX.XX1k-pkpm1n1mjnjIKXXi1,2,3,“aj1,2,3,.“卩3-36)第四节中枢点电压及联络线功率控制的潮流计算1、中枢点电压控制的潮流计算系统中枢点电压设定值为Vsp，实际运行值为V，其差值为3-37)ii,V二VspViii

18、为了将中枢点电压限制在Vsp，i需要改变（部分）发电机的输出电压（或无功功率）。按照灵敏度分析的原理，参见（3-9）式写出灵敏度方程_BBB|_,V|-0-DDDiDGDBBB,V=0iDiiiGiBBB,V，QGDGiGGGG3-38)上式中包含了电网中的所有节点，下标G为其电压与节点i的电压有强关联的发电机（无功补偿装置）节点的集合，在潮流计算中作为PV节点；下标D为除节点i外的所有节点的集合，在潮流计算中作为PQ节点。并且假定当调整PV节点电压时，PQ节点与i节点的无功功率不发生变化。在（3-38）中对应于节点集D的常数项为0用高斯消去法消去D中的节点后，节点i和节点集G的常数项不发生变

19、化，得BBitiGBBGiGG,Vi,VG0，QG3-39)如果已知（3-38）系数矩阵的因子表，由下述方法结算可节省计算量。根据因子表的形成方法知（3-39）式系数矩阵因子表可以从（3-38）式系数矩阵的因子表中直接取出和节点i、节点集G相关的部分而得，并且有：BiiBGiBiGBGG10|巴01LGLL0D0LtGiGGGGGGddLT|iiiiGidLdLLt+LLtiiGiiiGiGiGGGG3-40)式中L、L、D、d是GiGGGGii3-38）系数矩阵因子表中的元素。由（3-40）得3-41)3-43)3-44)3-45)3-46)3-47)3-48)3-49)BdTOC o 1-

20、5 h z,iiiiBdLtiGiiGi将（3-41）的关系代入（3-39）第一方程式有AV-B-iBAV-LtAV（3-42）iiiiGGGiG式中，AV是由（3-37）确定的常数，Lt是一个行矢量。AV是待求的参与iGiG调节的所有发电机节点电压的改变量，满足（3-42）式的AV可以有无穷多组解。G为了得到一组定解将求解（3-42）转化为一个优化问题，如取控制量的变化量最小作为目标函数，构成如下包含等约束的优化问题：1minAVtAV，2GGs.t.AV+LtAV0iGiG为求解此优化问题，建立该优化问题的拉格朗日函数L-AVtAV+九（AV+LtAV）2GGiGiGL达极值的必要条件为AV+九L0QAVgGiG归AV+LtAV0视iGiG由(3-45)得TOC o 1-5 h zAV九LGGi将(3-47)代入(3-46)有(LtL)-1AVGiGii将(3-48)代入(3-47)得AV-(LTL)-1LAVGGiGiGii用（3-49）修正发电