第三章电力系统运行的灵敏度分析及应用第一节灵敏度分析分析在给定的电力系统运行状态下,某些量发生变化时,会引起其他变量发生多大变化的问题。这一问题当然可通过潮流计算来解决,但计算工作量大。采用灵敏度分析法,计算量小,并可揭示各量之间的关系。但变化量大时,灵敏度分析法的精度不能保证。一、灵敏度分析的基本方法1、常规计算方法电力系统稳态运行的潮流方程一般性描述为:),(0),(uxyyuxf(3-1)x为状态变量,如节点电压和相角;u为控制变量,如发电机输出功率或电压;y为依从变量,如线路上的功率。实际上,(3-1)中0),(uxf就是节点功率约束方程,),(uxyy是支路功率与节点电压的关系式。设系统稳态运行点为),(00ux,受到扰动后系统的稳态运行点变为),(00uuxx。为了求出控制量变化量与状态量变化量之间的关系,在),(00ux处将(3-1)按泰勒展开并取一次项,得:uuyxxyuxyyyuufxxfuxfuuxxf),(0),(),(0000000(3-2)将),(0),(00000uxyyuxf代入,有:uuyxxyyuufxxf0(3-3)uSuuySxyuuyxxyyuSuufxfxyuxuxu1(3-4)其中uySxySufxfSxuyuxu1(3-5)为u的变化量分别引起x和y变化量的灵敏度矩阵。如果控制变量为各节点的有功、无功设定量,则1...11diaguf,所以,xuS就是潮流方程的雅可比矩阵的逆。xu,为两个不同状态间的变化量。2、准稳态灵敏度计算方法考虑到电力系统运行的实际:(1)初始控制变量的改变量,与到达新稳态的最终改变量不同;(2)一个控制量的变化可能使另一些控制量也发生变化。所以控制变量的初始改变量与最终改变量不同,表示为:)0(uFuu(3-6)由此得到准稳态的灵敏度关系:)0()0()0()0(uSuFSuSyuSuFSuSxRyuuyuyuRxuuxuxu(3-7)第二节潮流灵敏度矩阵1、发电机母线电压改变量GV与负荷母线电压改变量DV之间的灵敏度关系节点注入无功的平衡量方程0Q)cossin(QiiijijjijijijijijjiBVBGVV(3-8)上式简化依据了电力系统结构和运行的特点。根据灵敏度分析的基本方法,将(3-8)在当前状态点泰勒展开舍去高次项,的受到扰动后各变量变化量之间的关系ijjijiijjijiVBQVBQ0写成矩阵形式,并将负荷节点与发电机节点分开排列GDGDGGGDDGDDQQVVBBBB(3-9)(3-9)式与P-Q分解法V-Q迭代的修正方程式形式一致。但要注意在这里DQ、GQ是发电机和负荷的变化量。即(3-9)式表示了系统新稳态相对于旧稳态控制量的变化量与状态量的变化量之间的关系。假定GV调整后,负荷的无功功率不变化,即0QD,则式(3-9)第一式为:0VBVBGDGDDD变换得GDGGDGDDDVSVBBV1(3-10)其中DGDDDGBBS1(3-11)为DV与GV之间的灵敏度矩阵。通过灵敏度矩阵可以知道哪些发电机对控制负荷母线电压最有效,从而实现对负荷电压的定量控制。几种情况讨论:(1)只调整部分发电机的电压,无功充足能维持电压不变(GV=0)的发电机对(3-9)式没贡献,可从DGB中划去发电机电压能维持不变的节点对应的列。(2)被控量为部分负荷节点,即其它负荷节点的电压不关心,可从DDB、DGB中高斯消去不关心电压变化的负荷节点。(3)无功已达界的发电机,不能作为控制变量,也不能维持节点电压不变,高斯消去这些发电机的节点。这些节点的0GQ。高斯消去是等值变换,直接划去是不考虑它的影响。2、发电机母线电压改变量GV,负荷母线电压改变量DV与发电机输出无功的改变量GQ之间的灵敏度关系将(3-9)变换为GDGGGDDGDDGDGGGDDGDDGDQQRRRRQQBBBBVV1(3-12)假定发电输出无功改变时,负荷的无功功率不变,即0QD,有GDGDQRV(3-13)GGGGQRV(3-14)DGR、GGR是灵敏度矩阵。几种情况讨论:(1)不是控制变量的PV节点,其电压可维持不变,可直接划去对应的行和列。(2)不是控制变量的PQ节点,输出无功不变,当电压会发生变化,可将对应节点高斯消去。(3)不关心的负荷节点,直接划去。3、负荷母线电压改变量DV与变压器变比改变量t之间的灵敏度关系将节点无功平衡方程重写如下0Q)cossin(QiiijijjijijijijijjiBVBGVV其中ijB是变压器变比的函数,不考虑节点注入无功的变化,将变压器变比作为控制变量,节点电压作为被控变量,写出灵敏度方程0ijljljijjijjijttBVVB(3-15)上式中ijt为之路ji,的变压器变比。写成矩阵形式,包括所有负荷节点,并假定发电机母线电压不变,即认为发电机无功充足,可维持电压不变。0tVBjijijDVtB(3-16)即tBVjijijDVtB1(3-17)B仅包含负荷节点。jijijVtB为(3-15)式中第二项所组成的矩阵,行对应负荷节点,列对应可调变压器支路。每列中只有两个非零元素,分别在变压器支路的两个端点上。如果变压器支路有一个端点为PV节点,则由于PV节点电压不变,所以对应该变压器的支路只有一个非零元素。第三节分布因子分析节点注入有功功率变化、支路开断(结构变化)与支路潮流变化的灵敏度。1、支路开断分布因子分布因子:支路l基态有功潮流为lP,支路l开断引起支路k功率变化量为lkP,两者之间的关系表示为:llklkPDP(3-18)lkD为分布因子。相似与无功平衡方程,由有功平衡方程可得节点有功注入变化量与节点电压相角变化量之间的灵敏度方程0PXθθBP或(3-19)0B是以x1为支路参数建立的导纳矩阵,X是0B的逆。考虑一条支路),(jil断开的情况。如图,假定支路开断不引起节点注入功率的变化,则支路开断后,新网络节点的注入功率变化量为0.........0llPPP(3-20)其中,节点i的改变量lliiPPPP)(,节点j的改变量lljjPPPP)(。(3-20)可表示为:lllTPPMP0...1...1...0(3-21)lM是节点-支路关联列矢量,行对应节点号,支路l离开节点元素为1,进入节点元素为-1,节点与支路无关元素为0。新网络的导纳矩阵变为TlllxMMB10,开端后节点电压相角的变化量由(3-19)得)(110PMMBθTlllx(3-22)利用矩阵求逆辅助定理TlllTlllTllTlllTllTlllcxxxηηXXMXMMXMXBMMBMMBBMMB1101101010110)()()((3-23)其中ijjjiijjijjiiiTjninjjijjiiijiTlTllllllllTllllXXXXXXXXXXXXXXXXxXxc20...1...1...0.........0...1...1...00...1...1...0)()(1111XXMMXMMXMη1P2PiPjPnPlPNllX为在原网络支路l两端节点i注入单位电流,节点j流出单位电流,其它节点注入电流为0的情况下,节点i与节点j的点位差,定义为端口ji的自阻抗。支路l开断后,支路),(nmk上有功潮流的变化量llklklTlllTkkTklkPDPxcxPMηηXMθM)((3-24)kM为支路k的节点-支路关联矢量。支路k与支路l之间的支路开断分布因子是lllklkkllllllklkkllllTllTklTkklTlllTklkxXxXxxXXXXxxXxcD/1/)/()/()(XMMXMMXMMMηηXM(3-25)其中mjninjmiTnjmjniminmTlTklkXXXXXXXXXXX0...1...1...0............0...1...1...00...1...1...011XXMM为在原网络支路l两端节点i注入单位电流,节点j流出单位电流,其它节点注入电流为0的情况下,支路k两端节点m与节点n的点位差,定义为端口ji与端口nm之间的互阻抗。若支路l开断后网络分解为互不连通的两部分,因这时lllxX使(3-25)无定义。推导:当有多条支路同时开断时,支路k的功率变化量mlllkkPDP1m为断开的支路数。2、发电机输出功率转移分布因子发电机输出功率转移分布因子定义为发电机输出功率变化引起支路潮流的变化量,表示为:iikikPGP(3-26)ikG为发电机输出功率转移分布因子。i为发电机号,k为支路号。假定发电机i输出功率变化后引起的功率不平衡完全由平衡节点吸收,其它节点的输出功率不变化,则节点电压相角的变化量:iiTiPPXXθ0......0(3-27)iX是阻抗矩阵X的第i个列矢量,X是直流潮流中0B矩阵的逆。支路k(两端节点号分别为m和n)上有功潮流的变化iknimikiiTkkTkikPxXXxPxPXMθM即knimiikiikikxXXGPGP(3-28)ikG为转移分布因子,miX、niX是X中的m行i列和n行i列的元素。3、准稳态发电机输出功率转移分布因子(1)准稳态发电机输出功率转移分布因子设Gn台发电机有功出力调整量为)0(GP,如果调整量之和不为零则产生功率不平衡,不平衡量为)0()0(GTGGiGiPP1如果大于零则为功率超额,否则为功率缺额。实际电网中功率不平衡由所有发电机按一定比例承担,Gn台发电机的承担系数矢量为Gα,并且01iGTGGiiα1(3-29)各发电机的实际调整量)0()0()0()0()0()0()(GuGTGGGGTGGGiGGGGGiPPFP1αIP1αPαPP即TGGGuGuG1αIFPFP)0((3-30)uF为GGnn的方阵,是准稳态响应的变换矩阵。将(3-28)推广到考虑多台发电机的情况,为knimiikniiikkxXXGPGP1写成矩阵的形式GTkkGkGGkkxPXeMGPG1(3-31)kM:支路k的关联矢量,X:包含所有节点在内的阻抗矩阵,Ge为GnN阶单位矩阵,每列都是一单位矢量,只在相应的发电机节点处有非零元1。GkG为一行矢量。knlmlknjmjknimiGkxXXxXXxXX...G将(3-30)代入(3-31)得uGkRGkGRGkkPFGGPG)0((3-32)为了使支路k的有功潮流变化kP,各发电机有功的调整量可由伪逆计算如下。kTRGkRGkTRGkGP1)0(])([)(GGGP(3-33)由kP反求(0)GΔP时,(3-32)式是一个不定方程,可由