CN103368175A - 电力***动态稳定在线评估方法 - Google Patents

电力***动态稳定在线评估方法 Download PDF

Info

Publication number
CN103368175A
CN103368175A CN2013102804756A CN201310280475A CN103368175A CN 103368175 A CN103368175 A CN 103368175A CN 2013102804756 A CN2013102804756 A CN 2013102804756A CN 201310280475 A CN201310280475 A CN 201310280475A CN 103368175 A CN103368175 A CN 103368175A
Authority
CN
China
Prior art keywords
matrix
oscillation mode
mode
frequency
oscillation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2013102804756A
Other languages
English (en)
Other versions
CN103368175B (zh
Inventor
汪可友
李国杰
江秀臣
谢云至
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201310280475.6A priority Critical patent/CN103368175B/zh
Publication of CN103368175A publication Critical patent/CN103368175A/zh
Application granted granted Critical
Publication of CN103368175B publication Critical patent/CN103368175B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/22Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

一种电力***动态稳定在线评估方法,包括:步骤一,利用PMU进行高频采样,利用矩阵束辨识方法对采样数据提取振荡模式;步骤二,进行低频振荡检测,如果检测到的振荡模式阻尼比低于所设阈值,将发出低频振动警告;步骤三,通过频域分析方法对当前运行方式下***模型进行小扰动计算,获得特征值、特征向量、对应机组和参与因子信息;步骤四,根据矩阵束辨识方法确定的低频振荡模式,在小扰动计算结果中筛选出该振荡模式信息,并调整相关发电机出力,同时通过时域仿真得到相关区域间联络线潮流稳定极限。本发明有效避免噪声的影响,改善了辨识的精度,提高了计算速度和效率,明显改善了电力***在线动态安全评估效果。

Description

电力***动态稳定在线评估方法
技术领域
本发明涉及电力***在线安全稳定分析,特别是一种综合矩阵束辨识和频域特征值分析的电力***动态稳定在线评估方法。 
背景技术
随着我国电网建设继续跨越式大发展,截至2012年底,中国全国装机容量已超过美国跃居世界第一位,跨区互联电网规模也显著扩大,***稳定特性更趋复杂,表现为多种稳定问题相互交织,相互制约因素增多,从而使得增强电网结构加强稳定水平的效果不能简单评估,需要放在大***背景下通盘考虑。造成此现象的一个主要原因就是动态稳定水平在大规模交流互联电网中的各种稳定制约因素中的重要性不容忽视。比如总结近年来国调、各网调的大***仿真计算结论,小干扰计算结果基本上显示无弱或负阻尼机电振荡模式。但是***运行中还是难以杜绝的区域间功率振荡事件都表明运行中低频振荡风险始终存在。因此,动态稳定水平的评估就成为分析***稳定性、确定电网运行稳定控制限额的一个迫切需要解决的问题。 
目前国家电网公司***的调度部门已普遍在年度方式、日常方式等计算中开展了频域的小干扰稳定计算。其思想是对***的某一稳定运行点进行线性化处理,计算***状态矩阵的特征值。利用特征分析研究小扰动稳定问题有明显优势,可以求解特征根、特征向量、相关因子等大量与***动态稳定有关的重要信息,比如机组相关参与因子、相关比等。因此该法是复杂电力***动态稳定分析最有效的方法之一。然而从实践经验看,仅对某一运行方式进行小干扰计算是不完备的,由于低频振荡现象对发电机、励磁***及其它***参数敏感,并与***运行方式密切相关,使离线小干扰计算结果的可靠性和适应性不高。 
目前发展并逐步在电力***中应用起来的广域测量***(WAMS)能够利用相角测量装置(PMU)实时同步地测量电网不同地点的状态信息,并及时在调度中心数据平台汇总,对电网的动态过程进行监测,为实施区域间低频振荡模态信息的在线辨识和实时控制提供了有利条件。在线辨识方法是直接通过采样的数据序列进行计算后提取低频振荡模态信息。目前,在电力***中低频振荡模式辨识中应用比较广泛的算法是Prony算法,该算法是一种能够根据采样值直接估算出信号频率、衰减、幅值和初相位的分析方法,但是它对噪声尤为敏感。而矩阵束(Matrix Pencil)算法 可以去掉高次噪声的影响。因此,开展对基于WAMS的在线辨识算法研究,提高算法对噪声的抗干扰能力具有重要意义。 
电力***动态稳定在线评估要解决两个矛盾问题是:在线计算精度与速度的平衡、局部信息与全局信息的提取。目前在线辨识算法计算速度快,精度较好,但只能给出某个振荡模态的信息,没有全局信息,基于频域特征值分析的小干扰计算虽能给出全局信息,但受参数影响大,并且全部在线计算运算量过重。 
发明内容
为了克服上述现有技术的不足,本发明提供了一种电力***动态稳定在线评估方法。 
本发明所采用的技术方案是: 
一种电力***动态稳定在线评估方法,该方法包括如下步骤: 
步骤一,在电力***WAMS(广域测量***)中,利用PMU(同步相量测量单元)进行高频采样获得发电机功角信号或联络线功率曲线信号,然后利用矩阵束辨识算法对采样数据中的振荡模式进行提取; 
具体方法如下: 
y ( t ) = x ( t ) + n ( t ) ≈ Σ i = 1 M R i e ( s i t ) + n ( t ) - - - ( 1 )
其中,y(t)为通过PMU采样得到的发电机功角(或联络线功率曲线)信号,x(t)为***信号,n(t)为***中噪声,M为最大模态数,即振荡模式个数,第i个模式中,Ri为复参数,Sii+jωi为极点,αi为衰减因子,ωi为角频率; 
1)构造Hankel矩阵Y: 
将式(1)离散化获得采样时间序列: 
y ( kT s ) = Σ i = 1 M R i z i k + n ( kT s ) , k = 0,1 , . . . , N - 1 - - - ( 2 )
式中,
Figure BDA00003467243200023
为信号极点,Ts为采样周期,i为第i个采样点,k为第k个采样点,其中,是最大采样点,T为最大采样时间;根据采样值构造矩阵Y,L为矩阵束参数,通常设定为为N/4-N/3之间的整数, 
其中y(k)=y(kTs),k=0,1,…N-l 
2)求取最大模态数M, 
对矩阵Y进行奇异值分解:Y=UDVT
其中,U为酉矩阵,含有YYT的特征向量;V为酉矩阵,含有YTY的特征向量;D为对角阵,其主对角元素σc即为矩阵Y的奇异值,T表示转置; 
σ c σ max ≥ 10 - p , c = 1,2 , · · · L - 1 - - - ( 4 )
式中,σc即为矩阵Y的第c个对角元对应的奇异值,σmax为σc中最大的奇异值;p为十进制有效位数,可取3或4,此阀值可以有效防止噪声参与Y的重构。满足(4)式的奇异值的个数,即为M值。 
3)求取信号Z平面极点z的求取: 
构造以下两个矩阵: 
Figure 2013102804756100002DEST_PATH_IMAGE004
Figure 2013102804756100002DEST_PATH_IMAGE005
其中y(k)=y(kTs),k=0,1,…N-l 
计算矩阵 
G=Y1 +Y2   (7) 
式中,Y1 +为Y1的伪逆矩阵,矩阵G的M个非零特征值,即为矩阵Y2-λY1的广义特征值,也就是信号的极点zi,i=1,2,...,M; 
α i = - Re ( ln z i T s ) - - - ( 8 )
ω i = - Im ( ln z i T s ) - - - ( 9 )
式中αi和ωi分别为   
Figure 2013102804756100002DEST_PATH_IMAGE008
的实部与虚部。 
4)求取模式分量的参数Ri
Figure 2013102804756100002DEST_PATH_IMAGE009
利用最小二乘法求解Ri,到此Ri,αi,ωi,M均以得到,根据式(1)可以重构信号。 
步骤二,根据提取出来的振荡模式,进行低频振荡检测,根据检测到的振荡模式,如果该振荡模式阻尼比低于所设阈值,将发出低频振动警告: 
第i对特征根,一对共轭特征根代表一种振荡模式: 
λii±jωi,αi为衰减因子,ωi为角频率,计算阻尼比: 
ξ = - α α 2 + ω 2 - - - ( 11 )
若阻尼比ξ<0,负阻尼,***不可稳定运行; 
若阻尼比0<ξ<0.02,弱阻尼; 
若阻尼比0.02<ξ<0.03,较弱阻尼; 
若阻尼比0.03<ξ<O.04,适宜阻尼, 
若阻尼比ξ>0.05,正常阻尼,***运行的动态特性较好。 
步骤三,如步骤二发出低频振动警告,则启动频域分析算法对当前运行方式下***模型进行小扰动计算;获得特征值、特征向量、对应机组和参与因子信息; 
由***的微分代数方程 
x . = f ( x , y ) ; x ( 0 ) = x 0 - - - ( 13 )
0=g(x,y);y(0)=y0   (14) 
建立***模型的线性化状态方程   
Figure 2013102804756100002DEST_PATH_IMAGE012
式中Δx为***的状态变量增量的矢量形式,A为***矩阵: 
A = [ ∂ f ∂ x - ∂ f ∂ y ( ∂ g ∂ y ) - 1 ∂ g ∂ x ] - - - ( 15 )
根据***特征方程式|λI-A|=0对矩阵A特征根进行分析,得到A的特征值,左特征向量,右特征向量。再根据特征向量计算出参与因子矩阵: 
Figure 2013102804756100002DEST_PATH_IMAGE014
其中,pijijψji,φij为第j个右特征向量的第i个元素,ψji为第i个左特征 向量的第j个元素。 
步骤四,根据矩阵束辨识方法确定的低频振荡模式,在小扰动计算结果中筛选出该振荡频率信息,从而采取方式调整等措施提高该振荡模式的阻尼比,并同时通过时域仿真得到相关区域间联络线潮流稳定极限。 
根据矩阵束辨识算法确定的低频振荡模式,在小扰动计算结果中筛选出该振荡频率信息,根据参与因子Pij与特征根灵敏度的关系: 
p ij = ∂ λ j ∂ a ii - - - ( 17 )
推出对于所关心的低频振荡模式(假设第j个模式),对应的状态变量的参与程度高,再按图索骥由状态变量推出主导发电机或发电机群。调整主导发电机的有功出力,进行时域仿真,得到区域联络线的功率振荡曲线。根据功率振荡曲线判断稳定性,如果不稳定,继续调整主导发电机出力,直到得到相关区域联络线的潮流稳定极限。 
与现有技术相比,本发明的有益效果是 
本发明将矩阵束(Matrix Pencil)方法应用于低频振荡的在线辨识,不仅计算快速、辨识精度高、而且避免了Prony算法对噪声的敏感的缺点,大大增强了抗噪声干扰能力。 
本发明综合了低频振荡在线辨识和频域分析方法对电力***动态稳定性进行评估,很好克服了传统只利用单种方法的缺陷:无法达到精度与速度、局部信息获取与整体信息获取的平衡。本发明利用在线辨识算法计算速度快,精度较好的特点,先检测出低频振荡模态,然后再启动频域特征值分析计算,避免在线计算运算量过重,使得在线稳定评估可靠性高且不失效率。 
附图说明
图1为本发明的综合矩阵束辨识与频域分析的动态稳定在线评估流程。 
具体实施方式
下面结合附图对本发明进一步说明。 
步骤一,在电力***WAMS(广域测量***)中,利用PMU(同步相量测量单元)进行高频采样获得发电机功角信号或联络线功率曲线信号,然后利用矩阵束辨识算法对采样数据中的振荡模式进行提取: 
y ( t ) = x ( t ) + n ( t ) ≈ Σ i = 1 M R i e ( s i t ) + n ( t ) - - - ( 1 )
其中,y(t)为通过PMU采样得到的发电机功角(或联络线功率曲线)信号,x(t)为***信号,n(t)为***中噪声,M为最大模态数,即振荡模式个数,第i个模式中,Ri为复参数,Sii+jωi为S平面极点,αi为衰减因子,ωi为角频率; 
1)构造Hankel矩阵Y: 
将式(1)离散化获得采样时间序列: 
y ( kT s ) = Σ i = 1 M R i z i k + n ( kT s ) , k = 0,1 , . . . , N - 1 - - - ( 2 )
式中,   为Z平面极点,Ts为采样周期,i为第i个采样点,k为第k个采样点,其中,   
Figure DEST_PATH_IMAGE019
是最大采样点,T为最大采样时间; 
根据采样值构造矩阵Y,L为矩阵束参数,通常设定为为N/4-N/3之间的整数, 
Figure DEST_PATH_IMAGE020
其中y(k)=y(kTs),k=0,1,…N-l 
2)求取最大模态数M: 
对矩阵Y进行奇异值分解:Y=UDVT
其中,U为酉矩阵,含有YYT的特征向量;V为酉矩阵,含有YTY的特征向量;D为对角阵,其主对角元素σc即为矩阵Y的奇异值,T表示转置; 
σ c σ max ≥ 10 - p , c = 1,2 , · · · L - 1 - - - ( 4 )
式中,σc即为矩阵Y的第c个对角元对应的奇异值,σmax为σc中最大的奇异值;p为十进制有效位数,可取3或4,此阀值可以有效防止噪声参与Y的重构。满足(4)式的奇异值的个数,即为M值。 
3)求取Z平面极点z: 
构造以下两个矩阵: 
Figure DEST_PATH_IMAGE022
Figure DEST_PATH_IMAGE023
其中y(k)=y(kTs),k=,1,...N-T 
计算矩阵 
G=Y1 +Y2   (7) 
式中,Y1 +为Y1的伪逆矩阵,矩阵G的M个非零特征值,即为矩阵Y2-λY1的广义特征值,也就是信号的极点zi,i=1,2,…,M; 
α i = - Re ( ln z i T s ) - - - ( 8 )
ω i = - Im ( ln z i T s ) - - - ( 9 )
式中αi和ωi分别为   
Figure DEST_PATH_IMAGE026
的实部与虚部。 
4)求取模式分量的参数Ri
Figure BDA00003467243200093
利用最小二乘法求解Ri,到此Ri,αi,ωi,M均以得到,根据式(1)可以重构信号。 
步骤二,根据提取出来的振荡模式,进行低频振荡检测,根据检测到的振荡模式,如果该振荡模式阻尼比低于所设阈值,将发出低频振动警告; 
第i对特征根,一对共轭特征根代表一种振荡模式: 
λii±jωi,αi为衰减因子,ωi为角频率,计算阻尼比: 
ξ = - α α 2 + ω 2 - - - ( 11 )
若阻尼比ξ<0,负阻尼,***不可稳定运行; 
若阻尼比0<ξ<0.02,弱阻尼; 
若阻尼比0.02<ξ<0.03,较弱阻尼; 
若阻尼比0.03<ξ<0.04,适宜阻尼, 
若阻尼比ξ>0.05,正常阻尼,***运行的动态特性较好。 
步骤三,如步骤二发出低频振动警告,则启动频域分析算法对当前运行方式下***模型进行小扰动计算;获得特征值、特征向量、对应机组和参与因子等信息; 
由***的微分代数方程 
x . = f ( x , y ) ; x ( 0 ) = x 0 - - - ( 13 )
0=g(x,y);y(0)=y0   (14) 
建立***模型的线性化状态方程   
Figure DEST_PATH_IMAGE028
式中Δx为***的状态变量增量的矢量形式,A为***矩阵: 
A = [ ∂ f ∂ x - ∂ f ∂ y ( ∂ g ∂ y ) - 1 ∂ g ∂ x ] - - - ( 15 )
根据***特征方程式|λI-A|=0对矩阵A特征根进行分析,得到A的特征值,左特征向量,右特征向量。再根据特征向量计算出参与因子矩阵: 
Figure BDA00003467243200103
其中,pijijψji,φij为第j个右特征向量的第i个元素,ψji为第i个左特征向量的第j个元素。 
步骤四,根据矩阵束辨识方法确定的低频振荡模式,在小扰动计算结果中筛选出该振荡频率信息,从而采取方式调整等措施提高该振荡模式的阻尼比,并同时通过时域仿真得到相关区域间联络线潮流稳定极限。 
根据矩阵束辨识算法确定的低频振荡模式,在小扰动计算结果中筛选出该振荡频率信息,根据参与因子Pij与特征根灵敏度的关系: 
p ij = ∂ λ j ∂ a ii - - - ( 17 )
推出对于所关心的低频振荡模式(假设第j个模式),对应的状态变量的参与程度高,再按图索骥由状态变量推出主导发电机或发电机群。调整主导发电机的有功出力,进行时域仿真,得到区域联络线的功率曲线。如果功率曲线不稳定,将继续调整主导发电机出力,直到得到相关区域联络线的潮流稳定极限。 

Claims (5)

1.一种电力***动态稳定在线评估方法,其特征在于,该方法包括以下步骤: 
步骤一,在电力***WAMS中,利用PMU进行高频采样获得发电机功角信号或联络线功率曲线信号,然后利用矩阵束辨识方法对采样数据提取振荡模式; 
步骤二,根据提取出来的振荡模式,进行低频振荡检测,根据检测到的振荡模式,如果该振荡模式阻尼比低于所设阈值,将发出低频振动警告,并转入步骤三,否则,则表示电力***稳定; 
步骤三,通过频域分析方法对当前运行方式下***模型进行小扰动计算,获得特征值、特征向量、对应机组和参与因子信息; 
步骤四,根据矩阵束辨识方法确定的低频振荡模式,在小扰动计算结果中筛选出该振荡模式信息,并调整相关发电机出力,同时通过时域仿真得到相关区域间联络线潮流稳定极限。 
2.根据权利要求1所述的电力***动态稳定在线评估方法,其特征在于,所述的矩阵束辨识方法提取振荡模式的具体方法为: 
Figure 2013102804756100001DEST_PATH_IMAGE001
其中,y(t)为通过PMU采样得到的发电机功角信号或联络线功率曲线信号,x(t)为***信号,n(t)为***中噪声,M为最大模态数,即振荡模式个数,i为第i个振荡模式,Ri为复参数,sii+jωi为S平面极点,αi为衰减因子,ωi为角频率; 
2.1构造Hankel矩阵Y: 
将式(1)离散化获得采样时间序列: 
Figure 2013102804756100001DEST_PATH_IMAGE002
式中,
Figure FDA00003467243100012
为信号Z平面极点,Ts为采样周期,i为第i个采样点,k为第k个采样点,
Figure FDA00003467243100013
是最大采样点,T为最大采样时间,M为最大模态数; 
根据采样值构造矩阵Y,L为矩阵束参数,设定为N/4-N/3之间的整数, 
Figure 2013102804756100001DEST_PATH_IMAGE003
其中y(k)=y(kTs),k=0,1,…N-l 
2.2求取最大模态数M: 
对矩阵Y进行奇异值分解:Y=UDVT
其中,U为酉矩阵,含有YYT的特征向量;V为酉矩阵,含有YTY的特征向量;D为对角阵,其主对角元素σC即为矩阵Y的奇异值,T表示转置; 
Figure 2013102804756100001DEST_PATH_IMAGE004
式中,σC为矩阵Y的第c个对角元对应的奇异值,σmax为σc中最大的奇异值;p为十进制有效位数,取3或4,满足(4)式的奇异值的个数,即为M值。 
2.3求取信号Z平面极点z: 
构造以下两个矩阵: 
Figure 2013102804756100001DEST_PATH_IMAGE005
其中y(k)=y(kTs),k=0,1,…N-l 
计算矩阵 
G=Y1+Y2   (7) 
式中,Y1 +为Y1的伪逆矩阵,矩阵G的M个非零特征值,即为矩阵Y2-λY1的广义特征值,也就是信号Z平面极点zi,i=1,2,...,M; 
Figure DEST_PATH_IMAGE007
Figure DEST_PATH_IMAGE008
式中αi和ωi分别为   
Figure DEST_PATH_IMAGE009
的实部与虚部; 
2.4求取模式分量的参数Ri: 
Figure DEST_PATH_IMAGE010
利用最小二乘法求解Ri,到此Ri,αi,ωi,M均以得到。 
3.根据权利要求1所述的电力***动态稳定在线评估方法,其特征在于,所述的步骤二的具体方法为: 
根据步骤一提取信号的特征模式αii,推出***存在相应的振荡模式,第i对特征根,一对共轭特征根代表一种振荡模式: 
λii±jωi,αi为衰减因子,ωi为角频率, 
计算阻尼比: 
Figure DEST_PATH_IMAGE011
若阻尼比ξ<0,负阻尼,***不可稳定运行; 
若阻尼比0<ξ<0.02,弱阻尼; 
若阻尼比0.O2<ξ<O.03较弱阻尼; 
若阻尼比0.03<ξ<O.04,适宜阻尼, 
若阻尼比ξ>0.05,正常阻尼,***运行的动态特性较好。 
4.根据权利要求1所述的电力***动态稳定在线评估方法,其特征在于,所述的步骤三的的具体方法为: 
由***的微分代数方程 
Figure DEST_PATH_IMAGE012
0=g(x,y);y(0)=y0   (14) 
建立***模型的线性化状态方程   
Figure DEST_PATH_IMAGE013
式中Δx为***的状态变量增量的矢量形式,A为***矩阵: 
Figure DEST_PATH_IMAGE014
根据***特征方程式|λI-A|=0对矩阵A特征根进行分析,得到A的特征值,左特征向量,右特征向量; 
再根据特征向量计算出参与因子矩阵: 
Figure DEST_PATH_IMAGE015
其中,pijijψji,φij为第j个右特征向量的第i个元素,ψji为第i个左特征 向量的第j个元素。 
5.根据权利要求1所述的电力***动态稳定在线评估方法,其特征在于,所述的步骤四的的具体方法为: 
根据矩阵束辨识方法确定的低频振荡模式,在小扰动计算结果中筛选出该振荡频率信息,根据参与因子Pii与特征根灵敏度的关系: 
Figure DEST_PATH_IMAGE016
推出对于所关心的低频振荡模式,假设第j个模式,对应的状态变量的参与程度高,再按图索骥由状态变量推出主导发电机或发电机群;调整主导发电机的有功出力,进行时域仿真,得到区域联络线功率曲线;根据功率振荡曲线判断稳定性,如果不稳定,继续调整主导发电机出力,直到得到相关区域联络线的潮流稳定极限。 
CN201310280475.6A 2013-07-05 2013-07-05 电力***动态稳定在线评估方法 Active CN103368175B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310280475.6A CN103368175B (zh) 2013-07-05 2013-07-05 电力***动态稳定在线评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310280475.6A CN103368175B (zh) 2013-07-05 2013-07-05 电力***动态稳定在线评估方法

Publications (2)

Publication Number Publication Date
CN103368175A true CN103368175A (zh) 2013-10-23
CN103368175B CN103368175B (zh) 2015-09-09

Family

ID=49368960

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310280475.6A Active CN103368175B (zh) 2013-07-05 2013-07-05 电力***动态稳定在线评估方法

Country Status (1)

Country Link
CN (1) CN103368175B (zh)

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104037787A (zh) * 2014-05-13 2014-09-10 南方电网科学研究院有限责任公司 提高多小水电群送出能力的电力***稳定器协调优化方法
CN104201708A (zh) * 2014-08-01 2014-12-10 四川大学 基于改进矩阵束算法的多馈入直流***附加阻尼控制器的布点方法
CN104459409A (zh) * 2014-12-12 2015-03-25 国家电网公司 基于功率振荡次数评估***阻尼比的方法
CN105353305A (zh) * 2015-11-17 2016-02-24 浙江大学 基于完全自适应矩阵束的电机转子偏心故障诊断方法
CN105471110A (zh) * 2015-12-24 2016-04-06 国家电网公司 一种电力***类噪声信号的低频振荡模式分层检测方法
CN106410782A (zh) * 2015-07-31 2017-02-15 国家电网公司 一种基于阻尼贡献的互联电网动态稳定安控策略优化方法
CN106451436A (zh) * 2016-11-04 2017-02-22 广东电网有限责任公司电力科学研究院 一种综合量化评估发电机组动态特性的方法及***
CN106505587A (zh) * 2016-10-19 2017-03-15 福州大学 基于广义形态滤波与改进mp算法的低频振荡模态辨识方法
CN106786515A (zh) * 2016-12-01 2017-05-31 西南交通大学 一种电力***低频振荡模态分析方法
CN106776473A (zh) * 2016-12-16 2017-05-31 杭州电子科技大学信息工程学院 基于改进非线性变换的水电***频率非线性特性分析方法
CN108089095A (zh) * 2017-12-05 2018-05-29 广东电网有限责任公司电力科学研究院 一种电网低频振荡参数预测方法及装置
CN108110765A (zh) * 2018-01-20 2018-06-01 耿秋钰 一种新能源与电网配合供电方法及其校验装置
CN108270226A (zh) * 2018-01-20 2018-07-10 耿秋钰 一种快速恢复供电的方法
CN109001628A (zh) * 2018-06-05 2018-12-14 山东大学 基于同步相量数据的微型异步电动机组监控***、方法
CN109378837A (zh) * 2018-11-22 2019-02-22 中国南方电网有限责任公司 一种基于矩阵束法的超低频振荡辨识方法
CN109390957A (zh) * 2018-12-05 2019-02-26 江苏方天电力技术有限公司 一种风电功率波动诱发***强迫振荡的检测方法
CN109510217A (zh) * 2018-12-10 2019-03-22 东北电力大学 基于随机数据的电力***小干扰稳定分析在线识别方法
CN109507554A (zh) * 2018-12-10 2019-03-22 国网四川省电力公司电力科学研究院 一种电气设备绝缘状态评估方法
WO2019080393A1 (zh) * 2017-10-25 2019-05-02 东北电力大学 一种电力***动态稳定协同辨识方法
CN109711012A (zh) * 2018-12-14 2019-05-03 华北电力大学 一种基于奇异谱分析的pmu单通道丢失数据的恢复方法
CN110875599A (zh) * 2019-11-04 2020-03-10 云南电网有限责任公司 一种电网频率振荡的控制方法及***
CN110896217A (zh) * 2019-09-04 2020-03-20 中国电力科学研究院有限公司 一种提升受端电网电压稳定性的网源协调控制方法及装置
CN111426951A (zh) * 2020-06-09 2020-07-17 潍坊宏图环保设备有限公司 一种基于功率检测的发电机组电力输出稳定性检测***
CN112861328A (zh) * 2021-01-22 2021-05-28 东北电力大学 一种基于随机响应信号的发电机阻尼评估装置及方法
CN112861326A (zh) * 2021-01-21 2021-05-28 东北电力大学 一种基于量测的新能源电网发电机阻尼评估装置及方法
CN113036755A (zh) * 2021-03-09 2021-06-25 东北电力大学 基于随机数据驱动的区域联络线有功传输极限计算方法
CN113156247A (zh) * 2021-04-23 2021-07-23 北京建筑大学 一种电力***低频振荡的预警方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003018747A (ja) * 2001-06-29 2003-01-17 Mitsubishi Electric Corp 訓練用シミュレータおよび電力系統模擬方法
CN102170126A (zh) * 2011-04-06 2011-08-31 华中电网有限公司 一种电力***低频振荡在线分析及预警方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003018747A (ja) * 2001-06-29 2003-01-17 Mitsubishi Electric Corp 訓練用シミュレータおよび電力系統模擬方法
CN102170126A (zh) * 2011-04-06 2011-08-31 华中电网有限公司 一种电力***低频振荡在线分析及预警方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王宇静等: "《电力***振荡模态的矩阵束辨识法》", 《中国电机工程学报》 *
王官宏等: "《原动机调节***对电力***动态稳定的影响》", 《中国电机工程学报》 *

Cited By (44)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104037787B (zh) * 2014-05-13 2016-04-13 南方电网科学研究院有限责任公司 提高多小水电群送出能力的电力***稳定器协调优化方法
CN104037787A (zh) * 2014-05-13 2014-09-10 南方电网科学研究院有限责任公司 提高多小水电群送出能力的电力***稳定器协调优化方法
CN104201708A (zh) * 2014-08-01 2014-12-10 四川大学 基于改进矩阵束算法的多馈入直流***附加阻尼控制器的布点方法
CN104459409B (zh) * 2014-12-12 2017-04-12 国家电网公司 基于功率振荡次数评估***阻尼比的方法
CN104459409A (zh) * 2014-12-12 2015-03-25 国家电网公司 基于功率振荡次数评估***阻尼比的方法
CN106410782B (zh) * 2015-07-31 2019-03-08 国家电网公司 一种基于阻尼贡献的互联电网动态稳定安控策略优化方法
CN106410782A (zh) * 2015-07-31 2017-02-15 国家电网公司 一种基于阻尼贡献的互联电网动态稳定安控策略优化方法
CN105353305A (zh) * 2015-11-17 2016-02-24 浙江大学 基于完全自适应矩阵束的电机转子偏心故障诊断方法
CN105353305B (zh) * 2015-11-17 2018-01-12 浙江大学 基于完全自适应矩阵束的电机转子偏心故障诊断方法
CN105471110A (zh) * 2015-12-24 2016-04-06 国家电网公司 一种电力***类噪声信号的低频振荡模式分层检测方法
CN105471110B (zh) * 2015-12-24 2018-07-27 国家电网公司 一种电力***类噪声信号的低频振荡模式分层检测方法
CN106505587A (zh) * 2016-10-19 2017-03-15 福州大学 基于广义形态滤波与改进mp算法的低频振荡模态辨识方法
CN106505587B (zh) * 2016-10-19 2019-01-18 福州大学 基于广义形态滤波与改进mp算法的低频振荡模态辨识方法
CN106451436A (zh) * 2016-11-04 2017-02-22 广东电网有限责任公司电力科学研究院 一种综合量化评估发电机组动态特性的方法及***
CN106451436B (zh) * 2016-11-04 2019-03-12 广东电网有限责任公司电力科学研究院 一种综合量化评估发电机组动态特性的方法及***
CN106786515A (zh) * 2016-12-01 2017-05-31 西南交通大学 一种电力***低频振荡模态分析方法
CN106776473A (zh) * 2016-12-16 2017-05-31 杭州电子科技大学信息工程学院 基于改进非线性变换的水电***频率非线性特性分析方法
CN106776473B (zh) * 2016-12-16 2018-12-25 杭州电子科技大学信息工程学院 基于改进非线性变换的水电***频率非线性特性分析方法
US11271401B2 (en) 2017-10-25 2022-03-08 Northeast Electric Power University Synergic identification method for dynamic stability of power system
WO2019080393A1 (zh) * 2017-10-25 2019-05-02 东北电力大学 一种电力***动态稳定协同辨识方法
CN108089095B (zh) * 2017-12-05 2020-02-04 广东电网有限责任公司电力科学研究院 一种电网低频振荡参数预测方法及装置
CN108089095A (zh) * 2017-12-05 2018-05-29 广东电网有限责任公司电力科学研究院 一种电网低频振荡参数预测方法及装置
CN108270226A (zh) * 2018-01-20 2018-07-10 耿秋钰 一种快速恢复供电的方法
CN108270226B (zh) * 2018-01-20 2021-08-31 耿秋钰 一种快速恢复供电的方法
CN108110765A (zh) * 2018-01-20 2018-06-01 耿秋钰 一种新能源与电网配合供电方法及其校验装置
CN108110765B (zh) * 2018-01-20 2021-08-31 武汉西元智能科技有限公司 一种新能源与电网配合供电方法及其校验装置
CN109001628A (zh) * 2018-06-05 2018-12-14 山东大学 基于同步相量数据的微型异步电动机组监控***、方法
CN109378837A (zh) * 2018-11-22 2019-02-22 中国南方电网有限责任公司 一种基于矩阵束法的超低频振荡辨识方法
CN109390957A (zh) * 2018-12-05 2019-02-26 江苏方天电力技术有限公司 一种风电功率波动诱发***强迫振荡的检测方法
CN109390957B (zh) * 2018-12-05 2021-01-26 江苏方天电力技术有限公司 一种风电功率波动诱发***强迫振荡的检测方法
CN109507554A (zh) * 2018-12-10 2019-03-22 国网四川省电力公司电力科学研究院 一种电气设备绝缘状态评估方法
CN109510217A (zh) * 2018-12-10 2019-03-22 东北电力大学 基于随机数据的电力***小干扰稳定分析在线识别方法
CN109507554B (zh) * 2018-12-10 2020-11-24 国网四川省电力公司电力科学研究院 一种电气设备绝缘状态评估方法
CN109711012A (zh) * 2018-12-14 2019-05-03 华北电力大学 一种基于奇异谱分析的pmu单通道丢失数据的恢复方法
CN110896217A (zh) * 2019-09-04 2020-03-20 中国电力科学研究院有限公司 一种提升受端电网电压稳定性的网源协调控制方法及装置
CN110875599A (zh) * 2019-11-04 2020-03-10 云南电网有限责任公司 一种电网频率振荡的控制方法及***
CN111426951B (zh) * 2020-06-09 2020-09-11 潍坊宏图环保设备有限公司 一种基于功率检测的发电机组电力输出稳定性检测***
CN111426951A (zh) * 2020-06-09 2020-07-17 潍坊宏图环保设备有限公司 一种基于功率检测的发电机组电力输出稳定性检测***
CN112861326A (zh) * 2021-01-21 2021-05-28 东北电力大学 一种基于量测的新能源电网发电机阻尼评估装置及方法
CN112861328A (zh) * 2021-01-22 2021-05-28 东北电力大学 一种基于随机响应信号的发电机阻尼评估装置及方法
CN112861328B (zh) * 2021-01-22 2022-08-30 东北电力大学 一种基于随机响应信号的发电机阻尼评估装置及方法
CN113036755A (zh) * 2021-03-09 2021-06-25 东北电力大学 基于随机数据驱动的区域联络线有功传输极限计算方法
CN113156247A (zh) * 2021-04-23 2021-07-23 北京建筑大学 一种电力***低频振荡的预警方法及装置
CN113156247B (zh) * 2021-04-23 2024-04-23 北京建筑大学 一种电力***低频振荡的预警方法及装置

Also Published As

Publication number Publication date
CN103368175B (zh) 2015-09-09

Similar Documents

Publication Publication Date Title
CN103368175B (zh) 电力***动态稳定在线评估方法
CN103838959A (zh) 偏最小二乘回归应用于配电网谐波源定位与检测中的方法
CN102135571B (zh) 超高压/特高压多回路输电线路零序阻抗抗干扰测量方法
CN102819646B (zh) 一种线路舞动的电力***运行模拟方法
CN105653764A (zh) 电网安全事故风险等级在线评估与预警方法
CN104360235B (zh) 含dg复杂配电网的电压暂降源定位方法
CN103713217A (zh) 一种直流偏磁条件下电力变压器运行状态的在线监测方法
CN103474992A (zh) 电力***节点电压稳定状态的实时在线辨识判据
CN106127047B (zh) 一种基于延森-香农距离的电力***恶意数据检测方法
CN103245826A (zh) 配电线路雷电感应电压确定方法及***
CN104297628B (zh) 含dg的配电网的区段故障检测与定位方法
CN107919666B (zh) 一种基于广域响应的地区电网暂态稳定在线综合预判方法
CN104242267A (zh) 一种风力发电送出输电线路距离保护方法
CN104111405A (zh) 一种基于阻尼转矩分析法的电力***低频振荡源定位方法
CN100554976C (zh) 基于同步相量测量的区域电压稳定性监视方法
CN105067948A (zh) 一种小电流接地选线装置及单相接地检测方法
CN102033177B (zh) 一种配电网电力线路功率角测量方法及测量***
CN104808088A (zh) 一种基于雷电定位***记录与线路行波数据的雷电绕击与反击识别方法
Jin et al. WAMS light and its deployment in China
CN112018784B (zh) 一种基于同步相量测量数据的次同步谐振溯源方法
CN105676157A (zh) Wams低频振荡识别功能的测试***及方法
CN108802570A (zh) 一种用于交直流混联微电网的故障检测***及检测方法
Hosur et al. Subspace-driven output-only based change-point detection in power systems
CN106872861A (zh) 超特高压输电工程线路直击导线的雷电流反演方法及***
CN108828332A (zh) 一种计算雷电定位***探测效率的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant