CN103116106B - 一种电力***自由振荡和强迫振荡的辨别方法 - Google Patents
一种电力***自由振荡和强迫振荡的辨别方法 Download PDFInfo
- Publication number
- CN103116106B CN103116106B CN201310068235.XA CN201310068235A CN103116106B CN 103116106 B CN103116106 B CN 103116106B CN 201310068235 A CN201310068235 A CN 201310068235A CN 103116106 B CN103116106 B CN 103116106B
- 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.)
- Active
Links
- 230000010355 oscillation Effects 0.000 title claims abstract description 129
- 238000000034 method Methods 0.000 title abstract description 8
- 238000012544 monitoring process Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000013016 damping Methods 0.000 description 6
- 238000007405 data analysis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种电力***自由振荡和强迫振荡的辨别方法,属于电力***分析技术领域,该方法包括:***开始振荡后,监测发电机、变压器和线路的有功功率振幅,以振幅最大的有功功率数据进行分析;求取每个振荡周期内的有功功率的最大值和最小值,相减得到振荡的峰峰值,从起振时刻开始,选取峰峰值随时间增大的几个或十几个振荡周期进行分析;用各峰峰值减去前一个峰峰值,得到峰峰值的增量,对增量求取自然对数;对峰峰值增量的自然对数进行线性拟合,计算斜率指标,斜率指标大于零为自由振荡,斜率指标小于零为强迫振荡。本发明利用起振阶段波形特征辨别振荡类型,提高了振荡类型辨别的速度和准确性,帮助电网运行人员及时准确的辨别振荡类型,采取针对性的措施平息振荡。
Description
技术领域
本发明属于电力***分析技术领域,特别涉及一种电力***自由振荡和强迫振荡的辨别方法。
背景技术
低频振荡是威胁电力***安全运行的一个突出问题。当发生低频振荡时,电网中发电机的转子角、转速、输出功率,以及线路功率、母线电压等电气量发生持续的振荡,振荡频率一般在0.1~2.5Hz,严重威胁电网的安全。一旦发生振荡,电网运行部门需要及时采取有效的控制措施快速平息振荡。目前电网中发生的低频振荡,主要有两种类型,一种是由于***阻尼不足引起的弱阻尼甚至负阻尼自由振荡,另一种是持续的周期性小扰动引起的强迫振荡。两种振荡的处理方法不同,对于自由振荡,主要措施是增加***阻尼,当***阻尼为正时,振荡就会衰减,但对于强迫振荡,增加***阻尼并不能平息振荡,理论上只有将扰动源清除振荡才会完全平息。因此,振荡发生后,电网运行部门需要及时在线辨别振荡的类型,以采取针对性的措施平息振荡,但目前的振荡类型辨识方法还不成熟,不能满足电网实际运行的需要。
电力***发生低频振荡时,一般为单一主导模式,此时***的振荡可以用下面的单自由度***振动方程的一般形式进行分析:
自由振荡时,f(t)=0,***的振荡曲线为
电网发生低频振荡事故时,都有一个振幅增加的过程,称为起振阶段。对于自由振荡,σ1>0,***为负阻尼,起振阶段振幅指数增加,但振幅增大到一定程度后由于电网中的非线性环节,变为等幅振荡。
强迫振荡是在外施周期性扰动的作用下产生的振荡,主要关心共振的情况,即扰动源频率和***的固有频率相等,此时f(t)=μsin(ωt+θ),***的振荡曲线为
其中为自由分量,σ2<0,该分量随时间衰减,为强迫分量,不随时间衰减。当振荡进入稳态时自由分量完全衰减,只剩下强迫分量,振荡表现为等幅振荡。但在起振阶段也存在振幅增大的过程。
因此,仅从表现上看,自由振荡和强迫振荡很接近,起振阶段振幅增加,增加到一定程度后基本维持等幅振荡,这就给振荡类型的辨识带来了困难。
发明内容
本发明的目的是为解决已有技术对振荡类型辨识无有效方法的问题,提供一种电力***低频振荡的振荡类型辨别方法,在低频振荡发生后能快速准确的辨别振荡类型是自由振荡还是强迫振荡。
本发明提供一种电力***自由振荡和强迫振荡的辨别方法,其特征在于,包括以下步骤:
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),起振阶段振幅随时间的变化规律为
其中A1>0,σ1>0。
强迫振荡的振荡曲线为式(2),初始时x(0)=0,根据初始条件可得
两式平方相加得推导可得且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,由此可以辨别振荡类型。
根据线性拟合的最小二乘算法,拟合出的斜率:
记称为斜率指标,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,L,N,振荡周期的个数为N+1;
C:从第二个振荡周期开始,用该振荡周期的峰峰值减去前一个振荡周期的峰峰值,得到峰峰值的增量Yj=Xj-Xj-1,j=1,2,L,N,对Yj求取自然对数,得到Zj=lnYj;
D:计算斜率指标如果S>0,振荡类型为自由振荡,如果S<0,振荡类型为强迫振荡。
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 CN103116106A (zh) | 2013-05-22 |
CN103116106B true 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) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107102192B (zh) * | 2017-04-28 | 2019-07-23 | 许继集团有限公司 | 一种全相异步振荡判别方法及装置 |
CN108196146B (zh) * | 2017-12-26 | 2019-10-18 | 清华大学 | 电力***中低频振荡类型的判断方法 |
CN111626209A (zh) * | 2020-05-27 | 2020-09-04 | 江苏蔚能科技有限公司 | 基于高峰期电力***高低频振动实时监控方法 |
Citations (4)
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 | 武汉华中电力电网技术有限公司 | 一种负阻尼低频振荡与强迫功率振荡判别方法 |
-
2013
- 2013-03-04 CN CN201310068235.XA patent/CN103116106B/zh active Active
Patent Citations (4)
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)
Title |
---|
余一平等.基于能量函数的强迫功率振荡扰动源定位.《电力***自动化》.2010,第34卷(第5期),第1-6页. * |
梁志飞等.南方电网低频振荡控制策略探讨.《电力***自动化》.2011,第35卷(第16期),第54-58页. * |
陈磊等.基于振荡能量的低频振荡分析与振荡源定位(一)理论基础与能量流计算.《电力***自动化》.2012,第36卷(第3期),第22-27、86页. * |
陈磊等.基于振荡能量的低频振荡分析与振荡源定位(二)振荡源定位方法与算例.《电力***自动化》.2012,第36卷(第4期),第1-5、27页. * |
Also Published As
Publication number | Publication date |
---|---|
CN103116106A (zh) | 2013-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gargoom et al. | Automatic classification and characterization of power quality events | |
CN105651478A (zh) | 一种基于振动信号测试零部件疲劳寿命的分析方法 | |
CN103116106B (zh) | 一种电力***自由振荡和强迫振荡的辨别方法 | |
CN102928695A (zh) | 基于直线法判别负阻尼振荡与强迫振荡的方法 | |
CN109830972B (zh) | 一种新能源场站振荡源快速识别***及方法 | |
CN106772042B (zh) | 一种发电机原动***阻尼在线检测方法 | |
CN110176762A (zh) | 一种基于频域阻抗的次同步振荡风险在线评估方法及装置 | |
KR102161577B1 (ko) | 시간-주파수 이미지를 이용하여 회전체의 고장을 감지하는 장치 및 그 방법 | |
CN104620248A (zh) | 使用遗传编程估计剩余使用寿命 | |
Agrawal et al. | Electromechanical mode estimation in the presence of periodic forced oscillations | |
CN105092997B (zh) | 一种特高压输电线路雷电绕击与反击的识别方法 | |
CN104111405A (zh) | 一种基于阻尼转矩分析法的电力***低频振荡源定位方法 | |
CN103198184A (zh) | 一种电力***低频振荡特征类噪声辨识方法 | |
CN102928697A (zh) | 励磁调节器pss2a模型低频段阻尼检测方法及*** | |
CN102545245A (zh) | 基于端口供给能量的电力***振荡源定位方法 | |
Li et al. | Improved teager energy operator and improved chirp-Z transform for parameter estimation of voltage flicker | |
US20120303306A1 (en) | Method of subsynchronous resonance detection | |
CN106405210A (zh) | 一种采样数据异常检测方法及*** | |
CN105891625A (zh) | 基于能量流的电力***次同步振荡扰动源辨别方法 | |
CN104364623A (zh) | 驻极体型振动检测***、外部振动信息的生成方法、关于外部振动的传递函数信息的生成方法、外部振动信息的生成程序以及关于外部振动的传递函数信息的生成程序 | |
CN103493321A (zh) | 用于产生故障信号的方法 | |
CN110161344A (zh) | 基于尺度分解的变压器振声检测信号滤波方法和*** | |
CN103926494A (zh) | 一种干扰源的确定方法及装置 | |
CN106526359A (zh) | 基于Prony算法和病态数据分析的电网低频振荡在线检测算法 | |
CN109145773A (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 |