CN101187593A - 电站汽轮机转子振动故障的柯尔莫果洛夫熵诊断方法 - Google Patents
电站汽轮机转子振动故障的柯尔莫果洛夫熵诊断方法 Download PDFInfo
- Publication number
- CN101187593A CN101187593A CNA2007100326612A CN200710032661A CN101187593A CN 101187593 A CN101187593 A CN 101187593A CN A2007100326612 A CNA2007100326612 A CN A2007100326612A CN 200710032661 A CN200710032661 A CN 200710032661A CN 101187593 A CN101187593 A CN 101187593A
- Authority
- CN
- China
- Prior art keywords
- entropy
- fault
- golog
- ermo
- husband
- 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.)
- Pending
Links
Images
Landscapes
- Control Of Turbines (AREA)
Abstract
本发明涉及一种电站汽轮机转子振动故障的柯尔莫果洛夫(Kolmogorov)熵诊断的方法,基于混沌对***初始条件的敏感依赖性和相空间轨道的稠密性,通过研究***故障所引起的吸引子结构的动态变化,依靠混沌特征量—Kolmogorov熵对***状态变化的敏感性,来实现故障准确定位。采用小波包滤波和Kolmogorov熵诊断相结合的技术,先对原始振动数据进行小波包滤波,再对滤波后的数据赋值给一个时间序列,然后提取该时间序列的混沌特征量—Kolmogorov熵来进行故障定位。本发明算法简单,容易实现,可以有效克服采样过程中出现的噪声干扰,计算精度比较高,这些良好的性能大大增加了它的应用范围。
Description
技术领域
本发明涉及一种电站汽轮机转子振动故障的柯尔莫果洛夫(Kolmogorov)熵诊断方法,是一种算法简单,容易实现特别适用于有噪声干扰的电站汽轮机转子振动故障的诊断方法。
背景技术
汽轮机发电机组的故障率较高,而且故障危害性也较大,所以汽轮机的故障预测及诊断问题历来受到有关研究机构、企业和管理部门的高度重视,是现代故障诊断技术应用的一个重要方面。汽轮机转子振动波形提供了丰富的故障征兆信息,如何准确、全面地提取征兆信息对于故障类型的确定、故障发展的预测以及汽轮发电机组的状态检修都具有十分重要的意义。
在目前研究的汽轮机转子振动故障诊断方法当中,有以下2种典型方法。
(1)快速傅立叶变换(FFT)法
目前,应用最广泛的故障信号分析处理方法是傅立叶(Fourier)分析和相应的FFT快速算法,借助于FFT算法实现的信号处理有频谱分析、相关分析、相干分析、传递函数分析、细化谱分析、时间序列分析、倒频谱分析、包络分析等。这些分析方法在故障诊断过程中起到了重要的作用,比如参考文献“汽轮发电机组振动”。
但是,这种方法存在下述问题:
FFT只适合于分析连续的、平稳的时域信号,但实际中的很多信号是非线性、非平稳的。
(2)分形法
近年来,国内有一些学者提出用分形的方法来进行故障诊断,通过相空间重构理论,对满足随机分形统计自相似性的振动信号序列进行相空间重构,计算其关联维数,从而再现动力学特性。比如参考文献“基于分形关联维的汽轮机转子的振动故障诊断”(梁平,龙新峰,樊福梅,华南理工大学学报(自然科学版)第34卷第4期),针对故障振动位移首先进行了一个测点两个方向上的关联维数的计算,结合对应故障的波形图和频谱图进行分析。
但是,这种方法存在下述问题:问题一,分形无标度区的确定存在着很大的人为因素,导致最后结果精确度不高。
问题二,不能消除噪声的干扰。
综上所述,对于一个性能良好的腐蚀深度预测方法,必须符合下面几点基本要求:
(1)算法简单;(2)计算速度快;(3)精确度较高;(4)适合工程实际应用,等等。而现有技术尚未能很好的解决这些问题。
发明内容
本发明的目的在于针对现有技术的存在的上述不足,提供一种电站汽轮机转子振动故障的柯尔莫果洛夫(Kolmogorov)熵诊断的方法。
本发明基于混沌对***初始条件的敏感依赖性和相空间轨道的稠密性,采用小波包滤波和混沌特征量-Kolmogorov熵相结合的技术,将原始振动数据先用小波包分解,提取有用的频段进行分析,以达到滤波的目的,然后利用混沌动力学理论对滤波后的信号提取混沌特征量-Kolmogorov熵进行故障定位。该方法原理清晰简单,算法简单,实现速度快,并且能有效的滤掉噪声的干扰,对于故障类型能够精确的定位。
本发明的方法包括如下具体步骤:
(1)利用Bently实验装置模拟汽轮机转子的振动过程,并使其转速达到所需要求范围内。通过设置预置故障,经传感器及放大滤波采样,采集到在一定转速下各种常见转子故障的振动位移数据,经过数据的格式转换保存在相应的目录下,便于实验结束后的数据处理与分析。
由于汽轮机转子的四种典型故障:不平衡,不对中,碰摩,轴承松动的主要频谱特征主要分布在:0~100Hz,所以将原始信号进行三层小波包分解,选取第一个频带:0~125Hz进行分析,可以滤掉噪声和不需要的频率成分。然后再进行小波包重构,可以保留该小波包的数据,使信号恢复到原来的时域分辨率。
(2)将小波包重构结果赋值给一个时间序列,运用C-C方法计算出该时间序列的时间延迟。
(3)设置嵌入维数初值,例如从m=3开始,进行相空间重构,则此相空间为M维,M=N-(m-1)τ,其中N为时间序列序列的长度,对于M维相空间中的一对相点:
Xi={xi,xi+1,...,xi+M-1},Xj={xj,xj+1,...,xj+M-1 },当满足|xi-xj|<r0时,只要计算出{|xi+1-xj+1|,|xi+2-xj+2|,...,|xi+M-1-xj+M-1|}中小于r0的连续数据的个数,即为所求的b值,当|xi+t-1-xj+t-1|≤r0,且|xi+t-xj+t|>r0,t即为所要求的b。经过b步演化后找到距离大于r0的概率为:p(b)=c(b-1)-c(b)=(ek-1)e-kb这一概率密度函数被称为几何概率密度函数。
可以把这一概率密度函数的归一化表示为
利用b的概率密度分布,可以通过最大似然法推导出熵的表达式 (其中长度标度: )
(4)根据K熵的表达式来计算各种故障小波包滤波后的时间序列的Kolmogorov熵,并判断K熵是否已经趋于饱和,如为是,则该饱和值即为所要求的K熵值,否则使m=m+1,一直到K熵趋于饱和为止。
(5)通过K熵值来进行故障定位,可以发现不同故障类型的Kolmogorov熵值明显不同。
本发明的基于小波包滤波的汽轮机转子振动故障Kolmogorov熵诊断算法,为电站汽轮机转子振动故障的诊断开辟了一条新的路径。通过小波包滤波,来滤掉噪声的干扰,再计算Kolmogorov熵来进行故障定位。该方法算法简单,容易实现,精确度高,这些良好的性能大大增加了它的应用范围。
本发明的效果可以通过以下性能分析进行验证:
(1)将分形关联维方法的计算结果和Kolmogorov熵计算结果相比较进行分析,其中表1为3号通道各故障类型初始与故障下的分形关联维数增量,表2为3号通道各种故障状态下的Kolmogorov熵值。
由表1数据可知,3号通道各故障类型初始与故障下的分维数增量中,碰摩与松动故障下的分维数较初始状况的增量均较大,但不平衡故障与其初始状况的差别不大。从表2可以看出3号通道各种故障的Kolmogorov熵明显不同,有较好的故障区分度。
表1
表2
故障类型 | 时间延迟 | 嵌入维数 | Kolmogorov熵(nats/s) | 平均值(nats/s) |
不平衡1不平衡2不平衡3不对中1不对中2不对中3轴承松动1轴承松动2碰摩1碰摩2 | 4444444444 | 998121111710910 | 0.28050.31110.30900.22520.27300.25540.45430.54630.66200.6419 | 0.30020.25120.50030.6520 |
附图说明
图1为本发明电站汽轮机转子振动故障的Kolmogorov熵诊断的方法的流程图;
图2a、图2b、图2c和图2d分别为实施例中不平衡故障、不对中故障、碰摩故障和轴承松动的原始振动信号波形图,其中不平衡、碰摩、松动故障的转速为1500rpm,不对中故障的转速为1000rpm;
图3为原始振动信号的小波包分解示意图;
图4a、图4b、图4c和图4d分别为滤波之后的不平衡故障、不对中故障、碰摩故障和轴承松动振动信号的重构示意图;
图5为实施例中碰摩故障的Kolmogorov熵和嵌入维数的关系图。
具体实施方式
下面结合附图和实施例对本发明的实施方式作进一步说明。
实施例
如图1所示:
(1)根据Bently实验台所采集的四种典型汽轮机转子振动故障3号通道(测点1的y方向,测点通道如表3所示)的数据,先运用小波包分析进行滤波,再引用混沌动力学理论提取故障特征量-Kolmogorov熵,实验中各种故障的采样频率为2000Hz,则小波包分解频率为1000Hz,其中不平衡、碰摩、松动故障的转速为1500rpm,不对中故障的转速为1000rpm。
表3
测点 | 键相信号 | 测点1处x方向 | 测点1处y方向 | 测点2处x方向 | 测点2处y方向 |
通道号 | 1#通道 | 2#通道 | 3#通道 | 4#通道 | 5#通道 |
下面以四种故障的一种状况为例对其进行小波包分析,选取5040个振动数据作为研究对象,四种故障分别为不平衡故障、不对中故障、碰摩故障和轴承松动,该四种故障的振动信号波形图分别如图2a、图2b、图2c和图2d所示,其中X轴为时间数据序列,Y轴为振幅。
把原始信号进行3级小波包分解,在尺度3上形成了8个频带,如图3所示,则小波包分解可表示为:
d(0,0)=d(3,0)+d(3,1)+d(3,2)+d(3,3)+d(3,4)+d(3,5)+d(3,6)+d(3,7)。
各频段小波分解系数对应的频带如表4所示:
表4
小波系数 | 频带f/Hz | 小波系数 | 频带f/Hz |
d(3,0)d(3,1)d(3,2)d(3,3) | 0~125125~250250~375375~500 | d(3,4)d(3,5)d(3,6)d(3,7) | 500~625625~750750~875875~1000 |
由于不平衡,不对中,碰摩,轴承松动的振动故障主要频谱特征分别为:1X;1X、2X;2X、4X、3X、5X;1X、2X、3X(其中X为基频),由表2可看出,四种故障的振动频谱特征均主要分布在第一个频带,所以选取第一个频带进行分析,可以滤去噪声和不需要的频率成分。选取第一个小波包后再进行重构,可以保留该小波包的数据,使信号恢复到原来的时域分辨率。通过小波包分解及重构,得到滤波之后的四种振动信号(不平衡故障、不对中故障、碰摩故障和轴承松动)的重构示意图分别如图4a、图4b。图4c和图4d所示,其中X轴为时间数据序列,Y轴为振幅。
通过比较四种故障的波形图和重构示意图可以很明显的看出,小波包分析可以滤除数据采集时所产生的噪声和不需要的频率成分,并且保持了信息的完整性。
(2)将小波包重构结果赋值给一个时间序列,记为:x1,x2,x3,...,xn(n=1,2,...),运用C-C方法计算出该时间序列的时间延迟τ。
(3)然后设置嵌入维数初值,从m=3开始计算各种故障的Kolmogorov熵的最大似然估计值。
(4)判断K熵是否已经趋于饱和,如为是,则该饱和值即为所要求的K熵值,否则使m=m+1,一直到K熵趋于饱和为止。图5为碰摩故障的Kolmogorov熵和嵌入维数的关系图。由图5可以看出,当m=9以后,随着m的增大,K熵已经不再有明显的变化,则m=9时所对应的Kolmogorov熵即为所求。四种故障的Kolmogorov熵值计算结果如表5所示。
(5)由表5可以看出,小波包滤波后的四种故障振动信号的Kolmogorov熵值均大于零,说明在这四种故障状态下***均处于混沌状态;并且四种故障诊断下所计算的Kolmogorov熵值明显不同,碰摩故障状态下的Kolmogorov熵值最大,平均值为0.6520nats/s,松动故障的Kolmogorov熵值次之,不平衡时的较小,不对中时计算的Kolmogorov熵值最小,平均值仅为0.2512nats/s。
结论:将经过小波包滤波后的信号进行二次处理,提取混沌特征量-Kolmogorov熵进行故障定位,分析结果表明:不同故障状况下所求得的Kolmogorov熵值是明显不同的,K熵的数值越大,表示***的信息损失速率越大,***的混沌程度越大,K熵值可以作为***的一个状态特征量来进行故障诊断。
表5
故障类型 | 时间延迟 | 嵌入维数 | Kolmogorov熵(nats/s) | 平均值(nats/s) |
不平衡1 | 4 | 9 | 0.2805 | 0.3002 |
不平衡2 | 4 | 9 | 0.3111 | |
不平衡3 | 4 | 8 | 0.3090 | |
不对中1 | 4 | 12 | 0.2252 | 0.2512 |
不对中2 | 4 | 11 | 0.2730 | |
不对中3 | 4 | 11 | 0.2554 | |
轴承松动1 | 4 | 7 | 0.4543 | 0.5003 |
轴承松动2 | 4 | 10 | 0.5463 | |
碰摩1 | 4 | 9 | 0.6620 | 0.6520 |
碰摩2 | 4 | 10 | 0.6419 |
Claims (1)
1.一种电站汽轮机转子振动故障的柯尔莫果洛夫熵诊断方法,其特征在于包括如下步骤:
(1)对原始振动数据进行小波包分解,提取各种故障频谱特征所在的信号频段再进行小波包重构,从而达到滤波的目的;
(2)对滤波后的结果赋值给一个时间序列,运用C-C方法计算出该时间序列的时间延迟;
(3)设置该时间序列的嵌入维数初值,对该时间序列进行相空间重构,考察该相空间吸引子上的两个初始点,该两个初始点之间的初始距离小于时间序列的平均绝对偏差,记下该两个初始点分离至间距大于所述平均绝对偏差所需的时间;分别记下吸引子上初始距离小于所述平均绝对偏差的各点对分离至所述点对间距大于所述平均绝对偏差所需的分离时间,由所有点对的分离时间的平均值计算柯尔莫果洛夫熵的大小;
(4)增加嵌入维数初值,再计算在该嵌入维数初值下的柯尔莫果洛夫熵值,一直到柯尔莫果洛夫熵值不再随嵌入维数初值变化,即柯尔莫果洛夫熵值达到饱和,记下该饱和时的柯尔莫果洛夫熵值;
(5)通过比较不同故障状况下所求得的步骤(4)所述饱和时的柯尔莫果洛夫熵值,进行故障定位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2007100326612A CN101187593A (zh) | 2007-12-18 | 2007-12-18 | 电站汽轮机转子振动故障的柯尔莫果洛夫熵诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2007100326612A CN101187593A (zh) | 2007-12-18 | 2007-12-18 | 电站汽轮机转子振动故障的柯尔莫果洛夫熵诊断方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101187593A true CN101187593A (zh) | 2008-05-28 |
Family
ID=39480078
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2007100326612A Pending CN101187593A (zh) | 2007-12-18 | 2007-12-18 | 电站汽轮机转子振动故障的柯尔莫果洛夫熵诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101187593A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101839805A (zh) * | 2010-03-19 | 2010-09-22 | 同济大学 | 发动机主动故障质检与智能故障诊断方法 |
CN102435436A (zh) * | 2011-11-24 | 2012-05-02 | 电子科技大学 | 风扇轴承状态退化评估方法 |
CN101995290B (zh) * | 2009-08-28 | 2013-04-24 | 威能极公司 | 风力发电机振动监测的方法和*** |
CN103913272A (zh) * | 2014-04-03 | 2014-07-09 | 华北电力大学 | 一种汽轮发电机组转子质量不平衡故障的定量诊断方法 |
CN103968937A (zh) * | 2014-05-09 | 2014-08-06 | 国网福建晋江市供电有限公司 | 一种基于emd样本熵和fcm的配电开关机械状态诊断方法 |
CN106644918A (zh) * | 2016-10-21 | 2017-05-10 | 天津大学 | 基于混沌相空间重构理论的电化学噪声数据解析方法 |
CN109409341A (zh) * | 2018-12-10 | 2019-03-01 | 中国航发四川燃气涡轮研究院 | 一种近场的航空发动机噪声源辨识方法 |
CN110672319A (zh) * | 2019-10-11 | 2020-01-10 | 大连海事大学 | 基于吸引子的大型船舶轴系校中安全状态评估方法 |
CN111397896A (zh) * | 2020-03-08 | 2020-07-10 | 华中科技大学 | 一种旋转机械故障诊断方法、***及存储介质 |
CN112557004A (zh) * | 2020-11-23 | 2021-03-26 | 广西电网有限责任公司电力科学研究院 | 一种高压断路器机械故障检测方法 |
CN116046385A (zh) * | 2023-03-07 | 2023-05-02 | 天津汉云工业互联网有限公司 | 基于齿形图的齿轮故障识别方法、装置、介质和设备 |
-
2007
- 2007-12-18 CN CNA2007100326612A patent/CN101187593A/zh active Pending
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101995290B (zh) * | 2009-08-28 | 2013-04-24 | 威能极公司 | 风力发电机振动监测的方法和*** |
CN101839805A (zh) * | 2010-03-19 | 2010-09-22 | 同济大学 | 发动机主动故障质检与智能故障诊断方法 |
CN102435436B (zh) * | 2011-11-24 | 2015-05-27 | 电子科技大学 | 风扇轴承状态退化评估方法 |
CN102435436A (zh) * | 2011-11-24 | 2012-05-02 | 电子科技大学 | 风扇轴承状态退化评估方法 |
CN103913272B (zh) * | 2014-04-03 | 2016-05-18 | 华北电力大学 | 一种汽轮发电机组转子质量不平衡故障的定量诊断方法 |
CN103913272A (zh) * | 2014-04-03 | 2014-07-09 | 华北电力大学 | 一种汽轮发电机组转子质量不平衡故障的定量诊断方法 |
CN103968937A (zh) * | 2014-05-09 | 2014-08-06 | 国网福建晋江市供电有限公司 | 一种基于emd样本熵和fcm的配电开关机械状态诊断方法 |
CN106644918A (zh) * | 2016-10-21 | 2017-05-10 | 天津大学 | 基于混沌相空间重构理论的电化学噪声数据解析方法 |
CN109409341A (zh) * | 2018-12-10 | 2019-03-01 | 中国航发四川燃气涡轮研究院 | 一种近场的航空发动机噪声源辨识方法 |
CN110672319A (zh) * | 2019-10-11 | 2020-01-10 | 大连海事大学 | 基于吸引子的大型船舶轴系校中安全状态评估方法 |
CN110672319B (zh) * | 2019-10-11 | 2021-04-23 | 大连海事大学 | 基于吸引子的大型船舶轴系校中安全状态评估方法 |
CN111397896A (zh) * | 2020-03-08 | 2020-07-10 | 华中科技大学 | 一种旋转机械故障诊断方法、***及存储介质 |
CN112557004A (zh) * | 2020-11-23 | 2021-03-26 | 广西电网有限责任公司电力科学研究院 | 一种高压断路器机械故障检测方法 |
CN112557004B (zh) * | 2020-11-23 | 2022-06-14 | 广西电网有限责任公司电力科学研究院 | 一种高压断路器机械故障检测方法 |
CN116046385A (zh) * | 2023-03-07 | 2023-05-02 | 天津汉云工业互联网有限公司 | 基于齿形图的齿轮故障识别方法、装置、介质和设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101187593A (zh) | 电站汽轮机转子振动故障的柯尔莫果洛夫熵诊断方法 | |
Lu et al. | Effects of underdamped step-varying second-order stochastic resonance for weak signal detection | |
Feng et al. | Vibration signal models for fault diagnosis of planetary gearboxes | |
CN101221066B (zh) | 工程非线性振动检测方法 | |
Ou et al. | A new rolling bearing fault diagnosis method based on GFT impulse component extraction | |
CN105651504A (zh) | 基于自适应小波能量的旋转机械故障特征提取方法 | |
CN105784366A (zh) | 一种变转速下的风电机组轴承故障诊断方法 | |
CN103149029B (zh) | 利用倒频谱识别航空发动机转静碰摩部位的方法 | |
CN103575523A (zh) | 基于FastICA-谱峭度-包络谱分析的旋转机械故障诊断方法 | |
CN103091096A (zh) | 基于eemd和小波包变换的早期故障敏感特征提取方法 | |
CN102269644B (zh) | 基于最优自适应小波滤波器的滚动轴承与齿轮冲击性故障诊断方法 | |
Li et al. | Time–frequency analysis of typhoon effects on a 79-storey tall building | |
CN110926594B (zh) | 一种旋转机械信号时变频率特征提取方法 | |
Dahiya | Condition monitoring of wind turbine for rotor fault detection under non stationary conditions | |
Wang et al. | Bearing fault diagnosis of direct-drive wind turbines using multiscale filtering spectrum | |
Zhao et al. | Vibration health monitoring of rolling bearings under variable speed conditions by novel demodulation technique | |
Lin et al. | Instantaneous frequency estimation for wheelset bearings weak fault signals using second-order synchrosqueezing S-transform with optimally weighted sliding window | |
CN106053080A (zh) | 基于能量切片小波变换的滚动轴承故障特征提取方法 | |
CN104596766A (zh) | 一种轴承早期故障确定方法 | |
Wang et al. | Doppler Effect removal based on instantaneous frequency estimation and time domain re-sampling for wayside acoustic defective bearing detector system | |
Tian et al. | Rolling element bearing fault diagnosis using simulated annealing optimized spectral kurtosis | |
Zhang et al. | Doppler shift removal based on instantaneous frequency estimation for wayside fault diagnosis of train bearings | |
Di Lorenzo et al. | Dynamic characterization of wind turbine gearboxes using Order-Based Modal Analysis | |
Zhao et al. | Rolling element bearing fault diagnosis based on adaptive local iterative filtering decomposition and Teager–Kaiser energy operator | |
Cui et al. | Fault diagnosis of offshore wind turbines based on component separable synchroextracting transform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20080528 |