CN103116106A - 一种电力系统自由振荡和强迫振荡的辨别方法 - Google Patents

一种电力系统自由振荡和强迫振荡的辨别方法 Download PDF

Info

Publication number
CN103116106A
CN103116106A CN201310068235XA CN201310068235A CN103116106A CN 103116106 A CN103116106 A CN 103116106A CN 201310068235X A CN201310068235X A CN 201310068235XA CN 201310068235 A CN201310068235 A CN 201310068235A CN 103116106 A CN103116106 A CN 103116106A
Authority
CN
China
Prior art keywords
oscillation
peak
value
peak value
forced
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
CN201310068235XA
Other languages
English (en)
Other versions
CN103116106B (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.)
Tsinghua University
China Southern Power Grid Co Ltd
Original Assignee
Tsinghua University
China Southern Power Grid Co Ltd
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 Tsinghua University, China Southern Power Grid Co Ltd filed Critical Tsinghua University
Priority to CN201310068235.XA priority Critical patent/CN103116106B/zh
Publication of CN103116106A publication Critical patent/CN103116106A/zh
Application granted granted Critical
Publication of CN103116106B publication Critical patent/CN103116106B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种电力系统自由振荡和强迫振荡的辨别方法,属于电力系统分析技术领域,该方法包括:系统开始振荡后,监测发电机、变压器和线路的有功功率振幅,以振幅最大的有功功率数据进行分析;求取每个振荡周期内的有功功率的最大值和最小值,相减得到振荡的峰峰值,从起振时刻开始,选取峰峰值随时间增大的几个或十几个振荡周期进行分析;用各峰峰值减去前一个峰峰值,得到峰峰值的增量,对增量求取自然对数;对峰峰值增量的自然对数进行线性拟合,计算斜率指标,斜率指标大于零为自由振荡,斜率指标小于零为强迫振荡。本发明利用起振阶段波形特征辨别振荡类型,提高了振荡类型辨别的速度和准确性,帮助电网运行人员及时准确的辨别振荡类型,采取针对性的措施平息振荡。

Description

一种电力系统自由振荡和强迫振荡的辨别方法
技术领域
本发明属于电力系统分析技术领域,特别涉及一种电力系统自由振荡和强迫振荡的辨别方法。
背景技术
低频振荡是威胁电力系统安全运行的一个突出问题。当发生低频振荡时,电网中发电机的转子角、转速、输出功率,以及线路功率、母线电压等电气量发生持续的振荡,振荡频率一般在0.1~2.5Hz,严重威胁电网的安全。一旦发生振荡,电网运行部门需要及时采取有效的控制措施快速平息振荡。目前电网中发生的低频振荡,主要有两种类型,一种是由于系统阻尼不足引起的弱阻尼甚至负阻尼自由振荡,另一种是持续的周期性小扰动引起的强迫振荡。两种振荡的处理方法不同,对于自由振荡,主要措施是增加系统阻尼,当系统阻尼为正时,振荡就会衰减,但对于强迫振荡,增加系统阻尼并不能平息振荡,理论上只有将扰动源清除振荡才会完全平息。因此,振荡发生后,电网运行部门需要及时在线辨别振荡的类型,以采取针对性的措施平息振荡,但目前的振荡类型辨识方法还不成熟,不能满足电网实际运行的需要。
电力系统发生低频振荡时,一般为单一主导模式,此时系统的振荡可以用下面的单自由度系统振动方程的一般形式进行分析:
m x · · + c x · + kx = f ( t )
自由振荡时,f(t)=0,系统的振荡曲线为
Figure BDA00002881518900012
电网发生低频振荡事故时,都有一个振幅增加的过程,称为起振阶段。对于自由振荡,σ1>0,系统为负阻尼,起振阶段振幅指数增加,但振幅增大到一定程度后由于电网中的非线性环节,变为等幅振荡。
强迫振荡是在外施周期性扰动的作用下产生的振荡,主要关心共振的情况,即扰动源频率和系统的固有频率相等,此时f(t)=μsin(ωt+θ),系统的振荡曲线为
Figure BDA00002881518900013
其中
Figure BDA00002881518900014
为自由分量,σ2<0,该分量随时间衰减,
Figure BDA00002881518900015
为强迫分量,不随时间衰减。当振荡进入稳态时自由分量完全衰减,只剩下强迫分量,振荡表现为等幅振荡。但在起振阶段也存在振幅增大的过程。
因此,仅从表现上看,自由振荡和强迫振荡很接近,起振阶段振幅增加,增加到一定程度后基本维持等幅振荡,这就给振荡类型的辨识带来了困难。
发明内容
本发明的目的是为解决已有技术对振荡类型辨识无有效方法的问题,提供一种电力系统低频振荡的振荡类型辨别方法,在低频振荡发生后能快速准确的辨别振荡类型是自由振荡还是强迫振荡。
本发明提供一种电力系统自由振荡和强迫振荡的辨别方法,其特征在于,包括以下步骤:
A:系统开始振荡后,监测发电机、变压器和线路的有功功率振幅,以振幅最大的有功功率数据进行分析;
B:求取每个振荡周期内的有功功率的最大值和最小值,两者相减得到振荡的峰峰值,从起振时刻开始,选取峰峰值随时间增大的几个或十几个振荡周期,将各个周期的峰峰值记为Xi,i=0,1,2,...,N,N+1为振荡周期的个数;
C:从第二个振荡周期(i=1)开始,用该振荡周期的峰峰值减去前一个振荡周期的峰峰值,得到峰峰值的增量Yi=Xi-Xi-1,i=1,2,...,N,对Yi求取自然对数,得到Zi=lnYi
D:计算斜率指标如果S>0,振荡类型为自由振荡,如果S<0,振荡类型为强迫振荡。
本发明的有益效果
本发明的电力系统低频振荡类型辨别方法,利用起振阶段的波形特征,辨别振荡是自由振荡还是强迫振荡,帮助电网运行人员在开始振荡后快速准确的辨别振荡的类型,从而采取针对性的措施平息振荡。所述方法提高了电网的低频振荡在线分析和控制技术水平,一旦电网发生振荡,电网运行人员能够及时的辨别振荡的类型,从而采取针对性的控制措施平息振荡,保证电网安全性。
附图说明
图1是本发明所述电力系统自由振荡和强迫振荡辨别方法的流程图。
具体实施方式
下面结合附图,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
本发明所述的电力系统自由振荡和强迫振荡辨别方法的原理为:主要利用起振阶段的波形特征辨别振荡类型。电网发生持续时间较长的低频振荡时,在振荡稳态自由振荡和强迫振荡都表现为等幅或基本等幅的振荡形式,很难加以区分,但在振荡的开始阶段,也就是振幅增大的阶段,自由振荡和强迫振荡的特征是不同的。
对于自由振荡,根据式(1),起振阶段振幅随时间的变化规律为
A ( t ) = A 1 e σ 1 t
其中A1>0,σ1>0。
强迫振荡的振荡曲线为式(2),初始时x(0)=0,
Figure BDA00002881518900033
根据初始条件可得
Figure BDA00002881518900034
Figure BDA00002881518900035
两式平方相加得
Figure BDA00002881518900036
推导可得
Figure BDA00002881518900037
且A2+A3=0,由于A3>0,所以A2=-A3<0。
振荡的振幅
其中A3>0,σ2<0。
峰峰值等于振幅的两倍,因此随时间的变化规律完全相同。起振阶段的峰峰值可以表示为统一形式Aeσt+B,对于自由振荡A>0,B=0,σ>0,对于强迫振荡A<0,B>0,σ<0。
峰峰值的增量
Yi=Xi-Xi-1=AeiσT-Ae(i-1)σT=A(eσT-1)e(i-1)σT
其中T为振荡周期。
求取自然对数后
Zi=lnYi=ln(A(eσT-1))+(i-1)σT=(ln(A(eσT-1))-σT)+σT·i
对实测并计算得到的Zi进行线性拟合,得到斜率σT,对于自由振荡σT>0,对于强迫振荡σT<0,由此可以辨别振荡类型。
根据线性拟合的最小二乘算法,拟合出的斜率:
σT = Σ i = 1 N i Z i - N + 1 2 Σ i = 1 N Z i N ( N + 1 ) N - 1 12
称为斜率指标,S和σT符号相同。则对于自由振荡S>0,对于强迫振荡S<0。
本发明计算斜率指标然后根据S的符号辨别振荡类型。
本发明提供一种电力系统自由振荡和强迫振荡的辨别方法实施例,具体包括以下步骤:
A:系统开始振荡后,监测发电机、变压器和线路的有功功率振幅,以振幅最大的有功功率数据进行分析;
B:求取每个振荡周期内的有功功率的最大值和最小值,两者相减得到振荡的峰峰值,从起振时刻开始,选取峰峰值随时间增大的几个或十几个振荡周期,将各个周期的峰峰值记为Xi,i=0,1,2,...,N,N+1为振荡周期的个数;
C:从第二个振荡周期(i=1)开始,用该振荡周期的峰峰值减去前一个振荡周期的峰峰值,得到峰峰值的增量Yi=Xi-Xi-1,i=1,2,...,N,对Yi求取自然对数,得到Zi=lnYi
D:计算斜率指标如果S>0,振荡类型为自由振荡,如果S<0,振荡类型为强迫振荡。
所述步骤B中,从起振时刻开始,选取峰峰值随时间增大的几个或十几个振荡周期的具体作法为:从起振时刻开始,依次计算各个振荡周期的峰峰值,直到峰峰值开始下降,或者峰峰值的数量满足要求。一般的,峰峰值的数量一般达到15个,若峰峰值仍未下降,也能满足计算要求。
本发明的电力系统低频振荡类型辨别方法,利用起振阶段的波形特征,辨别振荡是自由振荡还是强迫振荡,帮助电网运行人员在开始振荡后快速准确的辨别振荡的类型,从而采取针对性的措施平息振荡,保证电网安全性。
以上实施方式仅用于说明本发明,而并非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型,因此所有等同的技术方案也属于本发明的范畴,本发明的专利保护范围应由权利要求限定。

Claims (1)

1.一种电力系统自由振荡和强迫振荡的辨别方法,其特征在于,包括步骤:
A:系统开始振荡后,监测发电机、变压器和线路的有功功率振幅,以振幅最大的有功功率数据进行分析;
B:求取每个振荡周期内的有功功率的最大值和最小值,两者相减得到振荡的峰峰值,从起振时刻开始,选取峰峰值随时间增大的几个或十几个振荡周期,将各个周期的峰峰值记为Xi,i=0,1,2,...,N,N+1为振荡周期的个数;
C:从第二个振荡周期(i=1)开始,用该振荡周期的峰峰值减去前一个振荡周期的峰峰值,得到峰峰值的增量Yi=Xi-Xi-1,i=1,2,...,N,对Yi求取自然对数,得到Zi=lnYi
D:计算斜率指标
Figure FDA00002881518800011
如果S>0,振荡类型为自由振荡,如果S<0,振荡类型为强迫振荡。
CN201310068235.XA 2013-03-04 2013-03-04 一种电力系统自由振荡和强迫振荡的辨别方法 Active CN103116106B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310068235.XA CN103116106B (zh) 2013-03-04 2013-03-04 一种电力系统自由振荡和强迫振荡的辨别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310068235.XA CN103116106B (zh) 2013-03-04 2013-03-04 一种电力系统自由振荡和强迫振荡的辨别方法

Publications (2)

Publication Number Publication Date
CN103116106A true CN103116106A (zh) 2013-05-22
CN103116106B CN103116106B (zh) 2015-05-20

Family

ID=48414532

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310068235.XA Active CN103116106B (zh) 2013-03-04 2013-03-04 一种电力系统自由振荡和强迫振荡的辨别方法

Country Status (1)

Country Link
CN (1) CN103116106B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107102192A (zh) * 2017-04-28 2017-08-29 许继集团有限公司 一种全相异步振荡判别方法及装置
CN108196146A (zh) * 2017-12-26 2018-06-22 清华大学 电力系统中低频振荡类型的判断方法
CN111626209A (zh) * 2020-05-27 2020-09-04 江苏蔚能科技有限公司 基于高峰期电力系统高低频振动实时监控方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW200839255A (en) * 2007-03-22 2008-10-01 Joe-Air Jiang A multiple power quality events measurement system and method
CN102650657A (zh) * 2012-04-10 2012-08-29 中国电力科学研究院 基于二次差分法判断负阻尼振荡与强迫振荡的系统和方法
CN102914712A (zh) * 2012-10-18 2013-02-06 中国电力科学研究院 基于直线法判断负阻尼振荡与强迫振荡的系统和方法
CN102944798A (zh) * 2012-11-29 2013-02-27 武汉华中电力电网技术有限公司 一种负阻尼低频振荡与强迫功率振荡判别方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW200839255A (en) * 2007-03-22 2008-10-01 Joe-Air Jiang A multiple power quality events measurement system and method
CN102650657A (zh) * 2012-04-10 2012-08-29 中国电力科学研究院 基于二次差分法判断负阻尼振荡与强迫振荡的系统和方法
CN102914712A (zh) * 2012-10-18 2013-02-06 中国电力科学研究院 基于直线法判断负阻尼振荡与强迫振荡的系统和方法
CN102944798A (zh) * 2012-11-29 2013-02-27 武汉华中电力电网技术有限公司 一种负阻尼低频振荡与强迫功率振荡判别方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
余一平等: "基于能量函数的强迫功率振荡扰动源定位", 《电力系统自动化》 *
梁志飞等: "南方电网低频振荡控制策略探讨", 《电力系统自动化》 *
陈磊等: "基于振荡能量的低频振荡分析与振荡源定位(一)理论基础与能量流计算", 《电力系统自动化》 *
陈磊等: "基于振荡能量的低频振荡分析与振荡源定位(二)振荡源定位方法与算例", 《电力系统自动化》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107102192A (zh) * 2017-04-28 2017-08-29 许继集团有限公司 一种全相异步振荡判别方法及装置
CN108196146A (zh) * 2017-12-26 2018-06-22 清华大学 电力系统中低频振荡类型的判断方法
CN108196146B (zh) * 2017-12-26 2019-10-18 清华大学 电力系统中低频振荡类型的判断方法
CN111626209A (zh) * 2020-05-27 2020-09-04 江苏蔚能科技有限公司 基于高峰期电力系统高低频振动实时监控方法

Also Published As

Publication number Publication date
CN103116106B (zh) 2015-05-20

Similar Documents

Publication Publication Date Title
Gargoom et al. Automatic classification and characterization of power quality events
Follum et al. Detection of periodic forced oscillations in power systems
Rangel-Magdaleno et al. Hilbert spectrum analysis of induction motors for the detection of incipient broken rotor bars
Ashok et al. Experimental evaluation of cyber attacks on automatic generation control using a CPS security testbed
CN105606895A (zh) 电力系统次同步振荡成份的在线检测及滤除方法
CN110176762A (zh) 一种基于频域阻抗的次同步振荡风险在线评估方法及装置
CN103116106B (zh) 一种电力系统自由振荡和强迫振荡的辨别方法
CN109830972B (zh) 一种新能源场站振荡源快速识别系统及方法
CN106202977B (zh) 一种基于盲源分离算法的低频振荡模式分析方法
CN105203876A (zh) 一种利用支持向量机和相关分析的变压器在线监测状态评估方法
CN102928695A (zh) 基于直线法判别负阻尼振荡与强迫振荡的方法
CN101042426B (zh) 用于频率调整的方法和装置
CN102545245A (zh) 基于端口供给能量的电力系统振荡源定位方法
CN103198184A (zh) 一种电力系统低频振荡特征类噪声辨识方法
CN104620248A (zh) 使用遗传编程估计剩余使用寿命
CN104111405A (zh) 一种基于阻尼转矩分析法的电力系统低频振荡源定位方法
KR20180058113A (ko) 가전기기 고유값 추출을 이용한 가전기기 분류 시스템 및 방법
CN105891625A (zh) 基于能量流的电力系统次同步振荡扰动源辨别方法
CN105226637B (zh) 一种基于振荡分群辨识方法的低频振荡调度控制方法
CN105372568A (zh) 一种超声波发生系统及超声激励信号生成方法
CN104950843B (zh) 基于模型识别方法的流程工业厂级控制回路振荡源定位方法
Spyropoulos et al. Extraction of frequency information for the reliable screening of rotor electrical faults via torque monitoring in induction motors
CN103337866A (zh) 一种从随机响应数据中辨识电力系统低频振荡参数的方法
CN109145773A (zh) 一种多源趋势信息融合的故障预测方法
CN103493321A (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