CN107015147B - 一种发电机组低频振荡检测方法及装置 - Google Patents
一种发电机组低频振荡检测方法及装置 Download PDFInfo
- Publication number
- CN107015147B CN107015147B CN201710214064.5A CN201710214064A CN107015147B CN 107015147 B CN107015147 B CN 107015147B CN 201710214064 A CN201710214064 A CN 201710214064A CN 107015147 B CN107015147 B CN 107015147B
- Authority
- CN
- China
- Prior art keywords
- frequency
- function
- generating set
- oscillation
- matrix
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/34—Testing dynamo-electric machines
- G01R31/343—Testing dynamo-electric machines in operation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种发电机组低频振荡检测方法及装置,包括以下步骤,步骤一,获得最小二乘滤波器;步骤二,确定函数采样样本密度及长度;步骤三,计算样本空间所包含的所有函数特征;步骤四,频率扫描和门槛值比较得到振荡特征;步骤五,通过趋势追踪确定振荡函数。本发明提出了一种全新的发电机组电气参数振荡测量方法;不仅适用于电网,而且同样适用于有周期性变化的过程的振荡特征的检测;对保证机组和电网安全具有重要意义;根据本发明方法制造的发电机组振荡监测装置,具有巨大的经济效益和社会效益。
Description
技术领域
本发明涉及发电机组低频振荡检测技术领域,尤其涉及一种发电机组低频振荡检测方法及装置。
背景技术
随着我国特高压技术的发展,电网进入了以特高压直流(±800kV)和交流(1000kV)为跨省连接骨干、以500kV为省内连接骨干的特大电网时代。电网安全性、稳定性问题比以往更加重要。特别是,以特高压交流连接的特大同步电网,送出端与负荷端在事故扰动时的同步趋势,直接关系着大电网的稳定。电网稳定有两个层面,暂态稳定和稳态稳定。暂态稳定是指在数个周波内的稳定,稳态稳定是指在十数个或数十个周波后的稳定。暂态稳定不一定代表稳态稳定,稳态稳定也不一定代表着暂态稳定。而二者任何一个不稳定都会破坏另一个稳定。因此,暂态和稳态不稳定都是竭力要避免的。暂态稳定控制的技术手段是靠励磁及一次调频技术,稳态稳定的控制手段是靠能量平衡管理(EMS)。前者无法人工干预,后者可以通过人工或人工智能实现自动管理。
电网低频振荡是一种电气参数呈现的周期性,其周期从数秒到数分钟不等的非正常振荡状态,如振荡持续,将带来电网稳定的破坏,严重时造成机组解列、电网瓦解,后果十分严重。前几年,在特高压电网尚未出现时,我国多个区域电网就已经发生过低频振荡,2010~2011年山东威海电厂600MW机组就发生过多次机网振荡和机组解列现象,山东电网也多次发生过500kV线路低频振荡现象。目前,低频振荡现象严重时,经验丰富的值班员可以发现,而振荡初期以及不甚明显的振荡靠人工难以发现,只有通过分析计算机系统的历史数据才会发现。
低频振荡是一种系统现象,但对发电机组来说,这种危害十分严重。电网发生振荡时,由于发电机组具有自动功率跟随(AGC)功能,它不仅受励磁和一次调频的暂态控制以同步于电网频率,而且也将在AGC控制下,跟随电网振荡功率指令而做相应跟随,使发电机组的转子受到更大阻尼冲击,这对机组转子寿命和安全性带来严重影响。2013年印度曾因电网频率波动发生试运中的300MW机组转子断裂重大事故。我国能源局于2014年下发了2014年161号文件,其中第5.1.8款规定,发电机组必须同步装设扭振保护装置。
实际上,电网和发电企业已经认识到低频振荡对电网和机组的危害。近几年,为降低低频振荡对机组的影响,发电机组普遍增加了PSS(电力系统稳定控制器)。PSS的作用从其功能模型上看是一个一阶惯性环节(阻尼环节),可抑制电网功率的突变。但其作用范围是有限的,而实际的低频振荡范围可能超出了这个区间。如果有一种装置能准确检测更大范围的低频振荡,不仅可以弥补PSS的不足,而且可以真正地对电网和机组安全运行提供指导和帮助。虽然PSS有利于电网的稳定,但它的副作用也很大,就是降低了机组的负荷响应速度,这对于快速调峰是不利的。有了低频振荡检测结果后,在正常情况下可以切除PSS,价值更大。
发明内容
本发明的目的就是为了解决上述问题,提供一种发电机组低频振荡检测方法及装置,将发电机定子电压和电流作为检测对象,通过对快速采样和无误差计算形成的区间历史数据的滚动扫描和最小二乘参数辨识,获得所有周期性函数的特征,最终通过幅值筛选,根据最大值原则确定振荡函数及其频率。
为了实现上述目的,本发明采用如下技术方案:
一种发电机组低频振荡检测方法,包括以下步骤,
步骤一,获得最小二乘滤波器;
步骤二,确定函数采样样本密度及长度;
步骤三,计算样本空间所包含的所有函数特征;
步骤四,频率扫描和门槛值比较得到振荡特征;
步骤五,通过趋势追踪确定振荡函数。
所述步骤一中最小二乘滤波器的获得方法为,
步骤1.1,将信号表示为一个成为带有直流分量、谐波成分的复合信号;
步骤1.2,将直流分量按泰勒基数展开并取前两项,将谐波成分按三角函数展开;
步骤1.3,根据最小二乘拟合原则得到:将采样时刻ti带入步骤1.2的函数得到的结果与采样样本u(ti)之间相等,不含ti的参数为需要辨识的未知变量;
步骤1.4,将N个采样方程按照步骤1.3的等式进行整理得到矩阵方程A·X=U;其中A是一个N行,2(M+1)列的定参数矩阵,M为谐波分量的最高次数,X是一个2(M+1)行1列的变量矩阵,U是一个N行1列的定参数矩阵;
步骤1.5,令2(M+1)=N,使得A为方阵,则A-1即为所选择的最小二乘滤波器。
所述最小二乘滤波器为:
其中τ为第i次采样时刻,ω为基波角频率,N为采样个数,M为谐波分量的最高次数。
所述步骤二中,采样角频率大于2ω;采样频率大于2kω,其中ω为采样数据中待识别的原信号的角频率,k为待识别的谐波次数。
所述步骤三中,样本空间所包含的所有函数特征为变量矩阵X所包含的函数特征,包括基波函数的峰-峰值,基波函数的有效值,相对于采样时间窗起始点的函数初始相角。
所述变量矩阵X为:
其中U0为直流分量,λ为直流分量衰减时间常数,UM为第M
次谐波分量峰值;则基波函数的峰-峰值为:基波函数的有效值为:相对于采样时间窗起始点的函数初始相角为:其
中x3、x4分别为矩阵X第3、4行的元素。
所述步骤四中,对采样样本按照从高频到低频原则逐个计算扫描,扫描时按ω的奇数倍逐一计算每个特征频率函数峰值,大于门槛值的为振荡特征,ω为基波角频率;门槛值为设定的振荡门槛值。
每个特征频率函数峰值按照公式进行计算,振荡函数角频率为:(2×i-1)ω,i为峰值表中的序号,i=1,...,N。
所述步骤五的方法为,通过步骤四发现振荡特征后,在下一周期再次出现时进行再次确认,再次确认的方法为重复所述步骤四的过程。
采用所述一种发电机组低频振荡检测方法的检测装置,包括相互通信的检测单元和监控站,所述检测单元包括主CPU模块,所述主CPU模块与振动相位测量模块、电气测量模块、AO模块及DO模块通过总线通信;所述电气测量模块的输入信号来自发电机定子电压和电流。
本发明的有益效果:
1.本发明提出了一种全新的发电机组电气参数振荡测量方法;
2.本发明提出的振荡测量不仅适用于电网,而且同样适用于有周期性变化的过程的振荡特征的检测;
3.本发明提出的发电机组振荡测量方法对保证机组和电网安全具有重要意义;
4.根据本发明方法制造的发电机组振荡监测装置,具有巨大的经济效益和社会效益。
附图说明
图1为本发明检测系统结构图。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
一种发电机组低频振荡检测方法,包括以下步骤:
步骤1:设计最小二乘滤波器;
步骤2:确定函数采样样本密度及长度;
步骤3:计算样本空间所包含的所有函数特征;
步骤4:频率扫描和门槛值比较发现振荡特征;
步骤5:通过趋势追踪确定振荡函数。
所述步骤1)设计最小二乘滤波器的过程如下:
任何信号都可表示为一个成为带有直流分量、特征分量及其谐波成分的复合信号,具体可表示为如下形式:
式中:U0—直流分量
λ—直流分量衰减时间常数
Uk—第k次谐波分量峰值,k=1,...,M
θk—第k次谐波的初始相位角
ω—基波角频率,可以是任意角频率。
U1和θ1是对应于基频ω下的相量幅值和初始相位角,是基波在采样时刻的主要特征,也是检测所需的关键特征。
首先将公式(1)做如下处理:
将公式(1)中的U0e-λt按泰勒基数展开并取前两项,Uksin(kωt+θk)按三角函数展开,则公式(1)可表示为如下形式:
根据最小二乘拟合原则,该函数在采样时刻ti,与采样样本u(ti)之间必然相等,即采样方程为:
ti—第i次采样时刻。
由于采样时刻ti已知,因此,公式(4)中sin(kωti),cos(kωti)也可认为是已知量,U0、λ、Ukcos(θk)、Uksin(θk),k=1,...,M,是需要辨识的未知参量。如果将N个采样方程按公式(4)逐一列出,并得到如下所示一个矩阵方程(为方程简洁起见,将ti用τ表示):
将系数矩阵用A,未知变量矩阵用X,采样矩阵用U表示,则矩阵方程(5)又可表示如下:
A·X=U (6)
因此,A是一个N行,2(M+1)列的定参数矩阵,X是一个2(M+1)行1列的变量矩阵,U是一个N行1列的定参数矩阵。如果令2(M+1)=N,则A就一个方阵,且最高可识别谐波倍数M=N/2-1。
根据矩阵理论,系数矩阵A存在逆矩阵A-1的条件为:|A|≠0,且A*A-1=I。I为单位矩阵,即对角线元素为1。
因此,待测的未知参数矩阵为:
X=A-1·U (7)
A-1即为所选择的最小二乘滤波器。
所述步骤2)确定函数样本密度及长度的过程如下:
确定函数样本密度及长度的原则是采样样本与所能辨识的函数周期必须符合申农采样定理。根据申农采样定理,数字采样频率必须为模拟信号频率的两倍以上。如果要从采样数据中识别角频率为ω的原信号,那么采样角频率必须大于2ω。推而广之,如果要辨识其中的k次谐波,则采样频率必须大于2kω。
因此,申农采样定理是采样样本密度和长度的选择原则。当然,如果既要识别低频信号又要识别高频信号,那么就需要保留足够多的采样样本。
所述步骤3)计算样本空间所包含的所有函数特征的过程如下:
公式(7)中的变量矩阵X即为公式(5)左侧的待辨识的未知参数列矩阵:
如果用xi表示矩阵X的第i行的元素(i=1,2M+2),则:
x1=U0,即直流分量
x3=U1cosθ1
x4=U1sinθ1
x3就是待测基波矢量复数的实部,x4为待测基波矢量复数的虚部。因此,基波函数的峰-峰值为:
基波函数的有效值为:
相对于采样时间窗起始点的函数初始相角为:
所述步骤4)频率扫描和门槛值比较发现振荡特征的过程如下:
由于被测函数是未知的,即,不能确定公式(5)左侧的系数矩阵中的基波角频率ω。因此,必须对采样样本按照从高频到低频原则逐个计算扫描,扫描时按ω的奇数倍按公式(8)逐一计算每个特征频率函数峰值。因此,将获得一个扫描频率计算表U,其元素个数为符合申农定理所覆盖的所有频率(ω奇数倍)。
如果只定义峰值大于某个门槛值的为振荡特征,则从特征峰值表里便可以找到振荡特征。即,发生振荡的条件为:
Ui≥UZ (10)
UZ——振荡门槛值
则,振荡函数角频率为:(2×i-1)ω (11)
其中,i为峰值表中的序号,i=1,...,N
所述步骤5)通过趋势追踪确定振荡函数的过程如下:
当通过步骤4)按公式(10)、(11)进行频率扫描发现振荡特征后,属于首次发现振荡特征,由于振荡是一种持续性现象,因此,在下一个周期必定再次出现,因此,从可靠性角度,必须再次确认。
再次确认的计算方法就是对新的采样样本重复步骤4)的计算过程。
因此,振荡特征需要第一次的发现和第二次的确认。
采用所述一种发电机组低频振荡检测方法的检测装置,该装置针对典型600MW机组设计,由两部分组成。一部分为上位监控和历史数据记录单元,另一部分为信号检测单元。信号检测单元内集承了一个主CPU模块(TCU),一个振动相位测量模块,一个电气测量模块,两个信号和报警输出模块。系统结构如图1所示。
其中,电气测量模块的输入信号来自发电机定子电压和电流,模块内采用了高性能数字处理器(DSP),DSP内部集成了A/D转换器、CPU、程序存储器等。DSP高速采样电气信号,并计算获得电气相量特征。计算结果通过内部TCU-32总线通信到TCU中,本发明的所有计算都在TCU中实现。
TCU计算结果由监控站监测和保存。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (7)
1.一种发电机组低频振荡检测方法,其特征是,包括以下步骤,
步骤一,获得最小二乘滤波器;
步骤1.1,将信号表示为一个成为带有直流分量、谐波成分的复合信号;
步骤1.2,将直流分量按泰勒基数展开并取前两项,将谐波成分按三角函数展开;
步骤1.3,根据最小二乘拟合原则得到:将采样时刻ti带入步骤1.2的函数得到的结果与采样样本u(ti)之间相等,不含ti的参数为需要辨识的未知变量;
步骤1.4,将N个采样方程按照步骤1.3的等式进行整理得到矩阵方程A·X=U;其中A是一个N行,2(M+1)列的定参数矩阵,M为谐波分量的最高次数,X是一个2(M+1)行1列的变量矩阵,U是一个N行1列的定参数矩阵;
步骤1.5,令2(M+1)=N,使得A为方阵,则A-1即为所选择的最小二乘滤波器;
所述最小二乘滤波器为:
其中τ为第i次采样时刻,ω为基波角频率,N为采样个数,M为谐波分量的最高次数;
步骤二,确定函数采样样本密度及长度;
步骤三,计算样本空间所包含的所有函数特征;
步骤四,频率扫描和门槛值比较得到振荡特征;
步骤五,通过趋势追踪确定振荡函数;
所述步骤五的方法为,通过步骤四发现振荡特征后,在下一周期再次出现时进行再次确认,再次确认的方法为重复所述步骤四的过程。
2.如权利要求1所述一种发电机组低频振荡检测方法,其特征是,所述步骤二中,采样角频率大于2ω;采样频率大于2kω,其中ω为采样数据中待识别的原信号的角频率,k为待识别的谐波次数。
3.如权利要求1所述一种发电机组低频振荡检测方法,其特征是,所述步骤三中,样本空间所包含的所有函数特征为变量矩阵X所包含的函数特征,包括基波函数的峰-峰值,基波函数的有效值,相对于采样时间窗起始点的函数初始相角。
4.如权利要求3所述一种发电机组低频振荡检测方法,其特征是,所述变量矩阵X为:
其中U0为直流分量,λ为直流分量衰减时间常数,UM为第M次谐波分量峰值;则基波函数的峰-峰值为:基波函数的有效值为:相对于采样时间窗起始点的函数初始相角为:其中x3、x4分别为矩阵X第3、4行的元素。
5.如权利要求4所述一种发电机组低频振荡检测方法,其特征是,所述步骤四中,对采样样本按照从高频到低频原则逐个计算扫描,扫描时按ω的奇数倍逐一计算每个特征频率函数峰值,大于门槛值的为振荡特征,ω为基波角频率;门槛值为设定的振荡门槛值。
6.如权利要求5所述一种发电机组低频振荡检测方法,其特征是,每个特征频率函数峰值按照公式进行计算,振荡函数角频率为:(2×i-1)ω,i为峰值表中的序号,i=1,...,N。
7.采用权利要求1所述一种发电机组低频振荡检测方法的检测装置,其特征是,包括相互通信的检测单元和监控站,所述检测单元包括主CPU模块,所述主CPU模块与振动相位测量模块、电气测量模块、AO模块及DO模块通过总线通信;所述电气测量模块的输入信号来自发电机定子电压和电流。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710214064.5A CN107015147B (zh) | 2017-04-01 | 2017-04-01 | 一种发电机组低频振荡检测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710214064.5A CN107015147B (zh) | 2017-04-01 | 2017-04-01 | 一种发电机组低频振荡检测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107015147A CN107015147A (zh) | 2017-08-04 |
CN107015147B true CN107015147B (zh) | 2019-08-06 |
Family
ID=59445090
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710214064.5A Active CN107015147B (zh) | 2017-04-01 | 2017-04-01 | 一种发电机组低频振荡检测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107015147B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732504A (zh) * | 2018-03-29 | 2018-11-02 | 上海理工大学 | 全电荷区域内锂电池等效电路模型参数辨识与优化方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113533861A (zh) * | 2021-06-10 | 2021-10-22 | 中国电建集团华东勘测设计研究院有限公司 | 一种基于衰减指数函数模型的动态谐波相量测量方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5291128A (en) * | 1992-05-19 | 1994-03-01 | General Motors Corporation | Motor testing apparatus utilizing inertial loading |
JP2890991B2 (ja) * | 1992-09-01 | 1999-05-17 | 日本電気株式会社 | 発振回路 |
CN101865963B (zh) * | 2010-07-08 | 2012-08-22 | 山东电力研究院 | 一种新的继电保护装置的快速检测方法 |
CN102157914B (zh) * | 2011-04-02 | 2014-03-12 | 山东电力研究院 | 一种新的继电保护装置的半波快速检测方法 |
CN104242325B (zh) * | 2014-09-18 | 2016-05-25 | 国家电网公司 | 一种电力系统低频振荡模式参数辨识方法 |
CN104931129B (zh) * | 2015-05-21 | 2018-05-08 | 山东泰景电力科技有限公司 | 一种发电机组转子扭振测量系统及测量方法 |
CN105067948B (zh) * | 2015-07-15 | 2018-06-12 | 山东泰景电力科技有限公司 | 一种小电流接地选线装置及单相接地检测方法 |
CN106202977B (zh) * | 2016-08-17 | 2018-09-14 | 华南理工大学 | 一种基于盲源分离算法的低频振荡模式分析方法 |
CN106526359B (zh) * | 2016-10-21 | 2019-08-06 | 国网新疆电力有限公司电力科学研究院 | 基于Prony算法和病态数据分析的电网低频振荡检测方法 |
-
2017
- 2017-04-01 CN CN201710214064.5A patent/CN107015147B/zh active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732504A (zh) * | 2018-03-29 | 2018-11-02 | 上海理工大学 | 全电荷区域内锂电池等效电路模型参数辨识与优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107015147A (zh) | 2017-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Brahma et al. | Real-time identification of dynamic events in power systems using PMU data, and potential applications—models, promises, and challenges | |
Aubert et al. | Kalman-filter-based indicator for online interturn short circuits detection in permanent-magnet synchronous generators | |
CN105470970B (zh) | 电力系统中次同步和超同步振荡的保护方法及系统 | |
CN108051739B (zh) | 用于监测永磁电机转子失磁故障的方法和系统 | |
Barzegaran et al. | Fault diagnosis of the asynchronous machines through magnetic signature analysis using finite-element method and neural networks | |
Lin et al. | Recursive algorithm for real-time measurement of electrical variables in power systems | |
CN109374270A (zh) | 一种gis异常振动分析及机械故障诊断装置及方法 | |
Concari et al. | Induction machine current space vector features to effectively discern and quantify rotor faults and external torque ripple | |
CN103198184A (zh) | 一种电力系统低频振荡特征类噪声辨识方法 | |
CN107015147B (zh) | 一种发电机组低频振荡检测方法及装置 | |
Usama et al. | Design and implementation of a wavelet analysis‐based shunt fault detection and identification module for transmission lines application | |
Avdakovic et al. | Identifications and monitoring of power system dynamics based on the PMUs and wavelet technique | |
Samal et al. | New signal subspace approach to estimate the inter‐area oscillatory modes in power system using TLS‐ESPRIT algorithm | |
Koochi et al. | Coherency detection and network partitioning supported by wide area measurement system | |
CN113285471A (zh) | 海上风电场次超同步振荡源感知与定位方法、装置及设备 | |
Shair et al. | Monitoring of subsynchronous oscillation in a series‐compensated wind power system using an adaptive extended Kalman filter | |
CN113075474A (zh) | 一种电能测量系统及方法 | |
Vanfretti et al. | A PMU-based fast real-time sub-synchronous oscillation detection application | |
CN108303609A (zh) | 一种次同步振荡在线监测及扰动源定位方法及系统 | |
CN109359822B (zh) | 电子式电压互感器测量状态评估方法及系统 | |
Yektaniroumand et al. | Optimal rotor fault detection in induction motor using particle-swarm optimization optimized neural network | |
CN115632410A (zh) | 一种针对新能源电力系统宽频振荡保护方法 | |
CN104931129B (zh) | 一种发电机组转子扭振测量系统及测量方法 | |
CN115313485A (zh) | 有源配电网孤岛检测方法、装置、电子设备及存储介质 | |
CN109557365A (zh) | 电网电压频率和相角检测方法及系统、单相电网检测系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |