法律状态公告日
法律状态信息
法律状态
2016-09-21
未缴年费专利权终止 IPC(主分类):H02J3/00 授权公告日:20130116 终止日期:20150729 申请日:20100729
专利权的终止
2013-01-16
授权
授权
2011-02-02
实质审查的生效 IPC(主分类):H02J3/00 申请日:20100729
实质审查的生效
2010-12-15
公开
公开
技术领域
本发明涉及电力系统动态稳定性分析方法,尤其涉及电力系统小干扰振荡稳定性分解式模态分析方法。
背景技术
电力系统中发生的低频振荡来源于发电机或发电机群转子之间的相对摇摆。它是当系统中出现局部的有功剩余或缺乏,相关的发电机或发电机群因此而加减速以求达到有功平衡时,由于阻尼缺乏而引起系统中发生的有功功率的持续振荡。振荡的频率范围一般在0.2~2.5Hz之间,故称为低频振荡,或机电振荡。近年来低频振荡在我国时有发生,严重影响了电网间的功率输送和安全稳定运行。在这种情况下,对低频振荡问题的研究备受关注。
电力系统低频振荡的机理分析和控制设计通常可以在系统的线性化模型上进行,然后通过非线性仿真来加以校验。所以,分析和设计面对的实际上是电力系统在小干扰下的低频振荡问题。在实践中,经常需要分析预测系统中某个指定的控制器或装置对系统小干扰下振荡稳定性的影响。而这一方面目前广泛使用的方法有二种:阻尼转矩分析方法(DTA:damping torque analysis)和模态分析方法(MA:modal analysis)。
基于经典控制理论的阻尼转矩分析方法建立在阻尼转矩的物理概念上,它是针对单机无穷大电力系统提出的,简单易懂。模态分析方法来源于现代控制理论的状态空间方法,是目前大规模复杂多机电力系统低频振荡分析和控制设计可以使用的常规方法,因为它适合于大型系统计算和分析,其编程简单,适用于软件的开发。目前流行的商业软件(包括电科院开发的商用软件PSASP和PSD-SSAP)中,低频振荡分析都采用模态分析法。模态分析方法的理论依据是线性系统的振荡模态,应而分析和控制设计都依赖于振荡模态的计算。其中用于分析预测电力系统中某个指定的控制器或装置对系统小干扰下振荡稳定性的影响主要是通过计算模态可控性指标,可观性指标和它们的乘积(残差)进行。
设在一个多机电力系统中,现需要分析预测某个指定的控制器或装置对系统的某个指定机电振荡模态阻尼的影响。该指定的控制器或装置可以是电力系统稳定器(PSS),某种灵活交流输电控制器(FACTS controller:flexible ac transmission systems controller),某种储能系统(ESS:energy storage system),某台电压调节器(AVR:automatic voltageregulator),或某条直流输电线(HVDC),等等。记该指定控制器或装置的输出信号为u(本文下文Δ都是指变量的偏差,变量上方加点都表示该变量的导数,变量下标加0都表示该变量在稳态的数值),输入信号为y,传递函数为T(s)(y=T(s)u)。指定分析的某个机电振荡模态为λi=-ξi±jωi。设含有指定控制器或装置的开环系统线性化模型为:
Δy=CΔX+DΔu
应有:
Avi=λivi,wiTA=wiTλi (2)
其中,λi为状态矩阵A的特征值(即指定分析的机电振荡模态),vi和wi分别为其相应的左右特征向量。根据模态控制理论,模态可控性指标bi,可观性指标ci和残差Ri定义为:
bi=wiTB,ci=Cvi,Ri=bici (3)
控制器或装置与模态的关系如图1所示。可见,模态可控性指标bi度量控制器或装置控制信号对模态能控性的大小;可观性指标ci度量模态在控制器或装置反馈信号中能被观察的程度。所以,残差Ri可以用来分析预测控制器或装置对模态的总影响。
显而易见,由残差Ri计算的模态分析给出的只是电力系统中控制器或装置与指定机电振荡模态的数学关系,它只是一个数学指标。这种模态分析不能表达指标的来源、分配和传递的物理过程,不能清晰地将分析计算与电力系统的物理量,特别是系统中的发电机这一与低频振荡的来源密切相关的具体装置联系起来,以加深对分析计算的理解。所以这种模态分析对应用于分析探索电力系统中的低频振荡机理,常常较为困难,因为它只给出一个数字计算结果。所以,对电力系统机电振荡的物理过程而言,它实际上是一种黑箱方法。
为使得模态分析与电力系统机电振荡的物理过程联系起来,电力研究和工作者一直在进行着不懈的努力。四十多年来,具有里程碑意义的进展是在八十年代由Perez-Arriaga,Verghese和Schweppe联合提出的选择模态分析方法(Selective modalanalysis with applications to electric power systems,part I and II,IEEE Transactions onpower apparatus and systems,No.9,Vol.PAS-101,pp3117-3134)。选择模态分析方法(由左右特征向量)通过计算参与性因子(participation factor)以揭示电力系统中发电机对机电振荡模态的灵敏度,从而成功地将模态分析与电力系统机电振荡的物理过程联系起来。如今,选择模态分析方法(计算参与性因子)已在大多数商业软件中采用。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种电力系统小干扰振荡稳定性分解式模态分析方法,将由残差Ri计算的模态分析及其计算结果与控制器或装置向电力系统中发电机提供的阻尼转矩和发电机对指定机电振荡模态的参与性(灵敏度)联系起来,使得残差Ri计算分析更加具有针对性,物理意义更清晰。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种电力系统小干扰振荡稳定性分解式模态分析方法,包括如下步骤:
(1)通过数据采集和监控系统SCADA系统、能量管理系统EMS获得电力系统稳态数据:发电机机端电压、机端有功、母线有功和母线无功;
(2)输入电力系统稳态数据:发电机机端电压、机端有功、母线有功和母线无功;
(3)输入电力系统静态数据:电力网络拓扑数据、线路阻抗导纳数据、变压器阻抗变比数据;
(4)发电机组固有数据:发电机内部电抗数据、励磁系统数据;
(5)利用潮流计算工具和初值计算工具,得到含指定控制器或装置的开环系统的线性化矩阵:
Δy=CΔX
以及:
其中,X为系统状态变量,A为系统线性化矩阵(不包含控制器自身),B为控制矩阵,C为输出矩阵,u为控制器输出信号,δ为发电机功角状态变量向量,ω为发电机转速状态变量向量,Z为系统其它状态变量向量,Δ为线性化算子,变量加点为微分算子,如为功角的导数,ω0I为对角阵,ω0为额定角速度;(4)式可如图2所示;
(6)计算指定控制器或装置向各台发电机提供的电磁转矩的通道因子Bj(λi),j=1,2,…N,即B(λi)的第j个分量,其中,N为系统中发电机的台数,λi=-ξi±jωi为系统中指定分析的某个机电的第i个振荡模态,由式(4)可得:
ΔZ=(sI-A33)-1(A31Δδ+A32Δω+B3Δu) (5)
所以
(6)式可如图3所示,根据图3可得针对第i个振荡模态的指定控制器或装置向各台发电机提供的电磁转矩向量为:
B(s)=B2+A23(sI-A33)-1B3
(7)计算指定控制器或装置的输入反馈信号用各台发电机转速状态变量重构的重构系数γj(λi),j=1,2,…N:
Δy=γj(λi)Δωj,j=1,2,…N (8)
因为在开环时,根据线性系统理论有
其中L为系统总的状态变量个数、λk为系统特征值、vk为相应的右特征向量、ak为系数。记
Δωj,j=1,2,…N也是系统的状态变量,所以应有
其中vk2j为vk中对应于Δωj,j=1,2,…N的分量。而输出变量为
所以
在上式中令s=λi可得
按式(14)即可计算指定控制器或装置的输入反馈信号用各台发电机转速状态变量重构的重构系数γj(λi),j=1,2,…N。
其中vi为相应于λi的右特征向量,vi2j为vi中对应于Δωj,j=1,2,…N的分量,j为第j台发电机,系统的输出方程为:
Δy=CX
(8)计算系统中各台发电机相应于指定振荡模态的参与因子Sij:
如果指定控制器或装置向各台发电机提供的阻尼转矩为DjΔωj,j=1,2,…N,针对第i个振荡模态,可以定义指定振荡模态对提供的阻尼转矩的灵敏度系数sij为:
显然,
上式表明:计算系统中各台发电机相应于指定振荡模态的参与因子即为式(15)定义的指定振荡模态对控制器或装置向各台发电机提供的阻尼转矩的灵敏度系数sij。
其中,DjΔωj,j=1,2,…N为各台发电机上获得的额外阻尼转矩,wi为相应于λi的左特征向量,wi2j为wi中对应于Δωj,j=1,2,…N的分量;
(9)通过对指定机电振荡模态的分解式模态分析方法计算残差:
(10)通过输出步骤(9)中得到的残差,安排指定控制器或装置在电力系统中的位置。
通过下面的证明过程可以反向证明式(17)的合理性:
记相应于λi的状态矩阵A的左特征向量wi为
wiT=[wi1Twi2Twi3T] (18)
由式(2)可得
由上式可得
wi3T=wi2TA23(λiI-A33) (20)
由式(3)、(4)、(5)和(7)得到可控性指标为:
由式(14)和(16)可有:
由式(3)和(21)可得:
证明完毕。
由式(7)可见,指定控制器或装置向各台发电机提供的电磁转矩为ΔTj=Bj(λi)Δu,j=1,2,…N。因为Δy=T(s)Δu,根据式(8),有:
ΔTj=Bj(λi)γj(λi)T(λi)Δωj,j=1,2,…N (24)
由式(16)和上式可见:在模态分析分解式(17)中,通道因子和重构系数的乘积Bj(λi)γj(λi)度量了指定控制器或装置向各台发电机提供的阻尼转矩;而参与因子度量了所提供的阻尼转矩是如何转化为对指定机电振荡模态的影响的。所以,模态分析分解式(17)的物理意义是:指定控制器或装置向每台发电机提供阻尼转矩,其大小由Bj(λi)γj(λi)度量;再乘以参与因子后,每台发电机获得的阻尼转矩就转化为指定控制器或装置通过各台发电机对指定机电振荡模态的影响,其大小由SijBj(λi)γj(λi)度量。残差Ri就是控制器或装置通过N台发电机对指定机电振荡模态的影响,为N项SijBj(λi)γj(λi)之和。这一模态分析分解式的物理意义如图4所示,指定控制器或装置首先向每台发电机G1、G2…GN提供阻尼转矩度量Bj(λi)γj(λi),然后每台发电机G1、G2…GN将获得的阻尼转矩转化对指定电机振荡模态影响的度量Sij。
有益效果:本发明的电力系统小干扰振荡稳定性分解式模态分析方法,能够清晰详细地给出指定控制器或装置是如何影响指定的电力系统机电振荡模态,从而为电力系统小干扰振荡稳定性分析和控制设计提供物理意义清晰的指导。
附图说明
图1为线性系统的模态分解示意图;
图2为系统的线性化状态方程示意图;
图3为系统的线性化状态方程压缩形式示意图;
图4为模态分析分解式的物理意义示意图;
图5为装有一台电池储能装置的简单四机二区域电力系统结构示意图;
图6为电池储能装置的数学模型示意图;
图7为电池储能装置装设阻尼控制器的阻尼控制效果非线性仿真结果示意图;
图8为无功阻尼控制器分解式模态分析结果示意图;
图9为无功阻尼控制器分解式模态分析结果示意图;
图10为一个实测大系统网络接线示意图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如附图5所示为简单四机二区域电力系统,通过图5所示系统说明本发明提出的模态分析的分解式计算方法的具体实施过程。模态分析的对象是安装在二区域联络线上的一台电池储能装置(BESS:battery energy storage systems)。模态分析的内容是预测在BESS上加设阻尼控制器的阻尼效果。四机二区域电力系统的参数为:
XT1=0.0027,XT2=0.0067,XT3=0.0027,XT4=0.0067,
X12=0.004,X26=0.22,X56=0.22,X45=0.01,X34=0.004,
L2=0.2,L5=4.0
Mi=8s.,Di=0,Xdi=0.8,Xdi′=0.2,Xqi=0.4,Td0i′=5,KAi=100,TAi=0.01
1.BESS线性化模型
根据已发表的文献(Chen Shen,Zhiping Yang,Crow,M.L,et al.,“Control ofSTATCOM with energy storage device”,Proceedings of IEEE PES Winter Meeting,Jan.2000,pp2722-2728),电池储能装置的数学模型如图6所示,为:
m=m0+KAC(s)(|Vs|-Vsref)+uq-pss (25)
φ=φ0+KDC(s)(Vdc-Vdcref)+up-pss
其中,Vs即为图6中电池储能装置在电力系统中接入点的电压V6,Vsx和Vsy分别为其在系统公共坐标下的x和y分量;Isx和Isy分别为电池储能装置注入电流在系统公共坐标下的x和y分量;Cdc和Vdc分别为AC/DC电压源转换器直流侧的电容和电压;VBESS和rBESS分别为电池的电压和等效电阻;m和φ分别为AC/DC电压源转换器脉宽调制算法(pulse width regulation)调制率和相位,m0和φ0为其相应的初值;KAC(s)和KDC(s)分别为由脉宽调制算法实现的交流和直流电压控制器的传递函数,uq-pss和up-pss分别为叠加在交流和直流电压控制器上的阻尼控制器输出控制信号;|Vs|为电池储能装置在电力系统中接入点电压的幅值;Vsref和Vdcref分别交流和直流电压控制的参考值;k为由AC/DC电压源转换器结构决定的常数。
记式(25)的最后一式的线性化为
Δγ=-Δφ+asTΔVsxy (26)
其中,由图6可有
Vsx+jVsy=jxs(Isx+jIsy)+Vc=jxs(Isx+jIsy)+mkVdccosγ+jmkVdccosγ(27)
由上式可得
利用式(26),可得上式的线性化为
ΔIsx=cx-dcΔVdc+cx-mΔm+cx-fΔφ+asxTΔVsxy (29)
ΔIsy=cy-dcΔVdc+cy-mΔm+cy-fΔφ+asyTΔVsxy
其中,
利用式(26)和(29)可得式(25)的第二式的线性化为
其中,
设:
交流电压控制器的传递函数KAC(s)的状态空间实现可以为:
Δm=CACXAC+DACΔVsxy+Δuq-pss
直流电压控制器的传递函数KDC(s)的状态空间实现可以为:
Δφ=CDCXDC+DDCΔVdc+Δup-pss
将式(31)和(32)代入(30)可得
由将式(31),(32)和(33)可得电池储能装置的状态方程为:
其中,
将式(31)和(32)代入(29)可得
ΔIsx=(cx-dc+cx-fDDC)ΔVdc+cx-mCACXAC+cx-fCDCXDC+(cx-mDAC+asxT)ΔVsxy+cx-mΔuq-pss+cx-fΔup-pss (35)
ΔIsy=(cy-dc+cy-fDDC)ΔVdc+cy-mCACXAC+cy-fCDCXDC+(cy-mDAC+asxT)ΔVsxy+cy-mΔup-pss+cy-fΔup-pss
记由式(34)和(35)可得:
ΔIsxy=CSXS+DSΔVsxy+Dq-pssΔuq-pss+Dp-pssΔup-pss (36)
其中,
式(34)和(36)即为BESS线性化模型。
2.全系统的线性化模型
设一个N机M节点电力系统中的节点电压和电流在系统坐标下分别表示为Vjx+jVjy和Ijx+jIjy,j=1,2,…M。不失一般性,设1至N号节点为发电机节点并记:
ΔVxy=[ΔV1x ΔV1y ΔV2x ΔV2y……ΔVNx ΔVNy]T
ΔIxy=[ΔI1x ΔI1y ΔI2x ΔI2y……ΔINx ΔINy]T
发电机的线性化模型是
ΔIxy=CgxyXg+DgxyΔVxy (37)
其中,Xg是发电机动态的状态变量。记除发电机和BESS接入节点(不失一般性,设为N+1号节点)外系统中的其它节点电压表示成的向量为:
ΔVoxy=[ΔV(N+1)x ΔV(N+1)y ΔV(N+2)x ΔV(N+2)y ……ΔVMx ΔVMy]T
将系统的网络方程写为:
其中gij+jbij,i,j=1,2,…M为系统中i和j号节点之间的导纳。将如上网络方程写为分块矩阵形式:
将式(36)和(37)代入(39)可得:
由上式可得:
其中,
由式(34)和(37)可得:
将式(41)代入(42)可得全系统的状态方程为:
其中,
3.阻尼控制器的输出方程
在图6所示电池储能装置上可以叠加二个阻尼控制器。一个叠加在交流电压控制器上,其输出控制信号是uq-pss,阻尼控制是通过调节电池储能装置和系统的无功功率交换实现的(以下称为无功阻尼控制器);另一个叠加在直流电压控制器上,其输出控制信号是up-pss,阻尼控制是通过调节电池储能装置和系统的有功功率交换实现的(以下称为有功阻尼控制器)。不失一般性,设二个阻尼控制器的反馈信号为电池储能装置安装地点线路上的有功功率偏差ΔPsj(即图5中所示ΔP65)。因为:
其中,*表示变量的共轭。由上式线性化可得:
利用式(41),由上式可得:
所以,关于阻尼控制器的输出方程为:
Δy=CX+DqΔuq-pss+DpΔup-pss (47)
其中,
4.模态分析分析结果及其分解
对图5所示的装有一台电池储能装置的简单四机二区域电力系统中的无功和有功阻尼控制器的阻尼效果预测的模态分析结果如表1所示。模态分析所关心的是系统中的区域振荡模态λa。在表1中也给出了阻尼控制器装设以后的特征值的计算结果。由表一结果可见:(1)在不同的联络线功率下,模态分析预测有功阻尼控制器的阻尼控制效果变化不大,而无功阻尼控制器的阻尼控制效果随联络线功率的增加而提高;(2)特征值的计算结果证实了模态分析预测的正确。非线性仿真的证实结果如图7所示,其中图7(a)和7(b)为装设无功阻尼控制器的阻尼控制效果非线性仿真结果;图7(c)和7(d)为装设有功阻尼控制器的阻尼控制效果非线性仿真结果;图7(a)和7(c)的P650=0.1p.u.,图7(b)和7(d)的P650=1.5p.u.。
表1模态分析预测结果和特征值计算验证:
可是,从表1的模态分析和特征值计算的结果无法得知:为什么无功和有功阻尼控制器的阻尼控制效果对联络线功率的改变有着不同的鲁棒性。所以,图8和图9给出了本发明分解式模态分析的计算结果。由图8和图9可见:(1)对无功(图8(a)和8(b))和有功(图9(a)和9(b))阻尼控制器,各台发电机的阻尼转矩灵敏度系数是一样的。这是因为发电机参与性与电池储能装置上装设的阻尼控制器的类型无关;(2)随着联络线输送功率的改变,各台机的阻尼转矩灵敏度系数(参与性)变化是不明显;(3)在不同的联络线功率下,有功阻尼控制器向各台机提供的阻尼转矩变化不大,所以阻尼控制的鲁棒性较好;(4)随着联络线输送功率的提高,无功阻尼控制器向各台机提供的阻尼转矩增加很多,所以阻尼控制的效果明显提高。因此,图8和图9的分解式模态分析分析给出了模态分析结果的物理解释。
从图8给出的分解式模态分析的计算结果还可见:无功阻尼控制器通过发电机G1和发电机G2对区域振荡模态λa的影响是正的(数值为负),但是通过发电机G3和发电机G4对区域振荡模态λa的影响却是负的(数值为正)。在联络线轻负荷时(如图8(a)所示P650=0.1p.u.时),通过发电机G3和发电机G4对区域振荡模态λa的影响更大,所以总的影响是负的。在联络线重负荷时(如图8(b)所示P650=1.5p.u.时),通过发电机G1和发电机G2对区域振荡模态λa的影响增加很多而大大超过通过发电机G3和发电机G4对区域振荡模态λa的影响,所以总的影响是正的。这也是为什么阻尼控制的效果随着联络线输送功率的提高而明显提高。
本发明提出的分解式模态分析方法在我国某一实际大系统上进行了测试。该实际大系统的简化接线图如图10所示。测试的对象是系统中将装设的一条高压直流线路(ZLC)的控制器(整流侧控制器和逆变侧控制器),目的是分析直流控制器对系统中的二个低频振荡模态(代码分别为FJM和AHM)的影响。
FJM模式主要表现为(区域1+区域2+区域3+区域4+区域5)对(区域6+区域7)的功率振荡。对FJM的模态分析计算的结果为:(1)ZLC整流侧控制器影响:Ri=0.4448-0.0509i;(2)ZLC逆变侧控制器影响:Ri=0.0016-0.0071i。AHM模式主要表现为(区域1)对(区域2+区域3+区域4+区域5区域6+区域7)的功率振荡。对AHM的模态分析计算的结果为:(1)ZLC整流侧控制器影响:Ri=-0.0630-0.0266i;(2)ZLC逆变侧控制器影响:Ri=-0.0001-0.0002i。这表明:
a.高压直流线路ZLC的整流侧控制器与其逆变侧控制器相比,对低频振荡模态的影响较大。这是因为整流侧作为功率送端,通过功率的调整,可以有效的影响功率振荡,而逆变侧为功率受端,是被动的功率变化,对功率振荡的影响有限。
b.对FJM而言,高压直流线路ZLC的控制对阻尼的影响是负的,即使得FJM向右半平面移动,阻尼变差;而对AHM而言,却是正的,但是影响与对FJM影响相比却小许多。对这一结果可以做进一步的分解式模态分析如下。
表2给出的是ZLC整流侧控制对FJM的模态分解式分析计算的结果。由表2中的结果可见:ZLC整流侧控制通过在区域6和7中的发电机对FJM的影响是使得FJM的阻尼增加,因为ZLC整流侧控制向区域6和7中的发电机提供正的阻尼转矩;但是,ZLC整流侧控制通过在其它区域中的发电机对FJM的影响却是使得FJM的阻尼减少许多,这样的总的结果是高压直流线路(ZLC)的控制使得FJM模态的阻尼减少。表3给出的是ZLC整流侧控制对AHM的模态分解式分析计算的结果。由表3中的结果可见:ZLC整流侧控制对AHM模态阻尼的影响较为复杂。它向在区域2和5中的大部分发电机(不是区域2和5中的全部发电机)提供正的阻尼转矩,从而使得AHM的阻尼增加较多。虽然ZLC整流侧控制向其它区域中的发电机提供负的阻尼转矩,但是总的影响是使得AHM的阻尼增加。所以,从分解式模态分析的结果可以清晰的看出整流侧控制器对FJM和AHM模态阻尼的影响是在全网如何分配,传递从而形成的。
表2:ZLC整流侧控制器对低频振荡模态FJM的分解式模态分析分析的详细结果
表3:ZLC整流侧控制器对低频振荡模态AHM的分解式模态分析分析的详细结果
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
机译: 基于多模态同步相量的电力系统振荡跟踪方法及系统
机译: 基于模态分析的电力系统实时振荡检测
机译: 基于模态分析的电力系统实时振荡检测