CN104298879A - 基于小波包分解的燃气轮机燃烧信号分析方法 - Google Patents

基于小波包分解的燃气轮机燃烧信号分析方法 Download PDF

Info

Publication number
CN104298879A
CN104298879A CN201410542171.7A CN201410542171A CN104298879A CN 104298879 A CN104298879 A CN 104298879A CN 201410542171 A CN201410542171 A CN 201410542171A CN 104298879 A CN104298879 A CN 104298879A
Authority
CN
China
Prior art keywords
combustion
wavelet packet
frequency
packet decomposition
gas turbine
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
CN201410542171.7A
Other languages
English (en)
Other versions
CN104298879B (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.)
Dongfang Electric Automatic Control Engineering Co ltd
Dongfang Electric Zhongneng Industrial Control Network Security Technology Chengdu Co Ltd
Original Assignee
SICHUAN ORIENT ELECTRIC AUTOMATIC CONTROL ENGINEERING 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 SICHUAN ORIENT ELECTRIC AUTOMATIC CONTROL ENGINEERING Co Ltd filed Critical SICHUAN ORIENT ELECTRIC AUTOMATIC CONTROL ENGINEERING Co Ltd
Priority to CN201410542171.7A priority Critical patent/CN104298879B/zh
Publication of CN104298879A publication Critical patent/CN104298879A/zh
Application granted granted Critical
Publication of CN104298879B publication Critical patent/CN104298879B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Of Engines (AREA)

Abstract

本发明公开了一种基于小波包分解的燃气轮机燃烧信号分析方法,使用传感器多通道采集燃烧信号,对采集的燃烧信号进行滤波、放大后经FPGA模块传送至CPU处理器,在CPU处理器内搭建小波包分解算法分频段地对燃烧信号进行分析,然后将分析结果分别传送至人机界面和控制系统。本发明能够准确采集燃气轮机燃烧压力脉动信号,并能根据采集的信号做出相应的指令以保证机组正常、稳定、安全的运行。

Description

基于小波包分解的燃气轮机燃烧信号分析方法
技术领域
本发明涉及燃气轮机数据采集分析领域,尤其涉及一种基于小波包分解的燃气轮机燃烧信号分析方法。
背景技术
随着能源可持续性发展战略的实施和环境保护意识的增强,燃气轮机技术引起了全世界的广泛重视。在燃气轮机关键技术中,预混燃烧技术的发展受到了燃烧不稳定性的限制,其表现为燃烧室出现的自激振荡现象,直接会对燃气轮机燃烧室火焰筒和过渡段的损坏,进而脱落的碎片会随气流进入透平段,造成二次破坏。燃气轮机燃烧室内的燃烧信号(压力波动)是表征燃气轮机燃烧稳定性的指标之一,对于燃气轮机燃烧信号稳定性的判定是通过对燃烧信号时-频转换、划分特征频带、确定频带上限、幅值判断的方式进行的。因此,如果准确高效地依据对燃烧室内壁各部位造成损伤的特征频点、特征频带和各特征频带上限,确定燃烧信号的最大频率解析范围                                                、最大采样频率、最小频率分辨率和频带上限,成为保护燃气轮机正常稳定运行的前提。其中,是自然数,表示燃烧信号的第采集通道;是自然数,表示在第通道下的第频带。
在现有技术中,对于燃烧信号的分析一般采用传统傅里叶变换(FFT)的方法进行,虽然FFT很容易将信号从时域转换到频域,但对于非稳定的燃烧信号,它存在如下缺陷:如果信号本身是非平稳的,FFT的过程可以理解它是在整个时间域上的一个积分,对于有限时间信号,类似做延拓后求取积分结果。如果频率是固定的,则延拓后该频率成分分布在整个时间轴上;如果某个频率成分随时间变化,那么延拓后,特定的频率成分显然是只占据了时间轴的一部分,如果某个信号持续时间很短,可以想象持续时间很短的信号有很大的幅度在频谱上可能反应不出来。而燃烧信号本身的非稳定性使得用FFT的方法对其分析存在局限性。
发明内容
本发明的目的在于克服现有技术中存在的上述问题,提供一种基于小波包分解的燃气轮机燃烧信号分析方法,本发明能够快速采集燃气轮机燃烧压力脉动信号,并能对采集的燃烧信号进行小波包分解,判断当前燃烧室内燃烧的稳定性,从而做出相应的指令以保证机组正常、稳定、安全的运行。
为实现上述目的,本发明采用的技术方案如下:
一种基于小波包分解的燃气轮机燃烧信号分析方法,其特征在于:使用传感器多通道采集燃烧信号,对采集的燃烧信号进行滤波、放大后经FPGA模块传送至CPU处理器,在CPU处理器内搭建小波包分解算法分频段地对燃烧信号进行分析,然后将分析结果分别传送至人机界面和控制系统。
所述小波包分解算法包括以下步骤:
(a)、根据燃烧信号中的采样频率和最小频率分辨率,通过公式确定小波包分解的层数n,小波包分解后在第n层对应的节点数为
(b)、根据燃烧信号确定的分解层数n,对燃烧信号进行小波包分解,并对每层分解的节点进行编号,将确定的最大频率解析范围作为起始节点(0,0),进而分解为(1,0)和(1,1),而后继续依次将其余节点进行分解,以此类推直到(n,0)、(n,1)(n,) ;
(c)、根据燃烧信号的特征频带和最小频率分辨率,通过小波函数对燃烧信号进行分解,确定特征频带所对应的频段数目和对应编号; 
(d)、对(c)步骤所确定的特征频带和隶属于该特征频带各节点的燃烧信号进行小波重构,从而获得特征频带中每一频率分量的小波重构结果,进而求得在此特征频带上,其重构均方根幅值的最大值所对应的特征频率点;
(e)、当特征频率点所对应的幅值超过该特征频率点所在特征频带的上限,则通过通讯接口分别将特征频率点及其对应幅值传送至人机界面报警和传送至控制系统进行燃烧调整,使这一特征频率点的脉动降低至上限范围内。
所述步骤(a)中,当对燃烧信号进行小波包分解层数为n时,解析到的最小频率分辨率是,解析到的最大频率分辨率为
所述步骤(b)中对燃烧信号的分解包括0—500HZ的频率解析和500—5000HZ的频率解析,对燃烧信号0—500HZ的频率解析在最小频率分辨率=1HZ下进行,其对应的小波包分解层数为12;对燃烧信号500—5000HZ的频率解析在频率分辨率 = 10HZ下进行,其对应的小波包分解层数为9,并调换小波包分解在第9层和第12层的各个频带节点的顺序,消除频谱混叠。
所述步骤(c)中的小波函数为db02小波函数。
采用本发明的优点在于:
一、本发明中,通过FPGA模块能够精确定时多通道燃烧信号采样频率,实现燃气轮机燃烧数据(压力波动)的多通道实时采集;通过分频段的实时小波包分解算法,对燃气轮机燃烧信号进行解析,求取其可能对燃烧室内壁各部位造成损伤的特征频点及其上限值,通过在多频率分辨率下解析燃烧信号,有效的提高了对于非稳定燃烧信号的检测精度。通过人机界面能够对燃气轮机的燃烧状态实时监视,通过控制系统能够根据采集的信号做出相应的指令以保证机组正常、稳定、安全的运行。另外,小波变换在时域和频域同时具有良好的局部化性质,而且由于对高频成分采用逐渐精细的时域或频域取样步长,从而可以聚焦到燃烧信号的任何细节,所以其在分析非稳定燃烧信号中显示了无可比拟的优越性。小波包变换建立在小波变换的基础上,实现频带均匀划分的同时避免了高频信号检测精度低的问题,提高了信号检测精度。
二、本发明中,所述步骤(b)中对燃烧信号的分解包括0—500HZ的频率解析和500—5000HZ的频率解析,对燃烧信号0—500HZ的频率解析在最小频率分辨率=1HZ下进行,其对应的小波包分解层数为12;对燃烧信号500—5000HZ的频率解析在频率分辨率 = 10HZ下进行,其对应的小波包分解层数为9,并调换小波包分解在第9层和第12层的各个频带节点的顺序,消除频谱混叠。这样既保证了对高频分量在低频率分辨率下的解析,又保证了低频分量在高频率分辨率下的解析,减小了数据运算量的同时,降低了CPU的负荷。
三、本发明中,由于小波分析中用到的小波函数没有唯一性,由此带来的问题是使用不同的小波基分析同一个问题会产生不同的结果,但目前没有一个选择最优小波基的同一方法,而主要是通过小波分析方法处理信号的结果与理论分析结果的误差来判定小波基的好坏,由此来选定小波基。通过对大量燃气轮机燃烧信号的分析,结果验证采用db02小波函数对燃烧信号的分析精度最高。
附图说明
附图1是燃气轮机模块化燃烧信号分析方法;
附图2是燃气轮机燃烧信号小波包分解流程图。
具体实施方式
图1是燃气轮机模块化燃烧信号分析方法实施方案,包括CPU处理器、输入/输出模块、FPGA模块、通讯模块、传感器、信号调理器、人机界面和控制系统,首先使用传感器多通道采集燃烧信号,采集到的燃烧信号经信号调理器进行滤波、放大后,由FPGA模块实时、调整采集,被采集到的燃烧数据通过PCI总线传送至CPU处理器,CPU处理器为实时操作系统,在其内部搭建小波包分解算法分频段地对燃烧信号进行分析,在分析结果中,其数据分为控制数据和显示数据,显示数据以以太网通讯方式传送至人机界面进行监视,控制数据以串口通讯的方式传送至控制系统进行综合处理,参与燃烧调整。通过对燃气轮机燃烧信号多通道采集,时—频分析,获取燃气轮机燃烧室内压力波动的上限值及其对应的特征频点,从而做出快速、有效的预测判断,进而进行燃烧调整,保障燃气轮机的稳定燃烧。其中,输入/输出模块与传感器配合实现对燃气轮机燃烧信号的多通道实时采集,FPGA模块通过精确定时控制燃烧信号的采样频率和采样点数N,通讯模块用于燃气轮机燃烧信号解析结果与外围设备的数据交互。
图2是燃气轮机燃烧数据小波包分解算法流程图,其具体步骤如下:
(a)、一方面根据燃气轮机燃烧信号的历史数据,确定用于解析燃气轮机燃烧信号的最大频率解析范围、最大采样频率、最小频率分辨率和上限阈值;另一方面,通过认知燃气轮机燃烧信号的最大频率解析范围,确定燃烧信号的采样频率、采样点数N、多频率分辨率)和),其中。根据燃烧信号中的采样频率和最小频率分辨率,通过公式确定小波包分解的层数n,小波包分解后在第n层对应的节点数为。低频分量的解析范围是0HZ—500HZ,高频分量解析范围是500HZ—5000HZ。
(b)、根据燃烧信号确定的分解层数n,对燃烧信号进行小波包分解,并对每层分解的节点进行编号,将确定的最大频率解析范围作为起始节点(0,0),进而分解为高频频段(1,0)和低频频段(1,1),而后继续依次将其余节点分高频频段和低频频段进行分解,以此类推直到(n,0)、(n,1)(n,) 。
(c)、根据燃烧信号的特征频带和最小频率分辨率,通过db02小波函数对燃烧信号进行分解,确定特征频带所对应的频段数目和对应编号;例如在采样频率,N = 1024的情况下,通过特定的db02小波函数对燃烧信号进行12层分解,同时确定对于燃烧信号的低频分量(0—500HZ)在最小频率分辨率=1HZ下进行解析,即位于小波包分解节点的12层,其节点编号从(12,0)—(12,499);高频分量(500—5000HZ)在频率分辨率 = 10HZ的下进行解析,即位于小波包分解节点的第9层,其节点编号从(9,50)—(9,499),调换小波包分解在第9层和第12层的各个频带节点的顺序,消除频谱混叠,达到正确的频带划分效果。从而实现对燃气轮机燃烧信号的多分辨率分析。如表1所示:
(d)、对确定的特征频带经过小波包分解、并对隶属于特征频带各节点的燃烧信号进行小波重构,别针对低频分量和高频分量进行,从而获得特征频带中每一频率分量的小波重构结果,进而求得在此特征频带上,其重构均方根幅值的最大值所对应的小波包节点,即特征频率点。
(e)、当特征频率点所对应的幅值超过该特征频率点所在特征频带的上限,则通过通讯接口分别将特征频率点及其对应幅值传送至人机界面报警和传送至控制系统进行燃烧调整,使这一特征频率点的脉动降低至上限范围内。
本发明中,所述步骤(a)中,当对燃烧信号进行小波包分解层数为n时,能够解析到的最小频率分辨率是,能够解析到的最大频率分辨率为,即可实现对燃烧信号的多分辨率分析
本发明中,小波分析的提出主要是针对非稳定信号的分析方法,在时间(空间)频率的局部化分析,它通过伸缩平移运算随信号(函数)逐步进行多尺度细化,最终达到高频处时间细分,低频处频率细分,能自动适应时频信号分析的要求,从而可以聚焦信号的任意细节。在时频平面上具有可变的时间和频率分辨率,即在高频区具有好的时间分辨率,在低频区具有好的频率分辨率。通过选择不同的小波函数可以得到不同的时频分布。往往对于非稳定信号的分析,是在多分辨率尺度下进行的,因而对于小波变换会进行延伸,这就是的小波包分解。小波变换这种独特的能力也使它成为一种分析非稳定信号的有力工具。

Claims (5)

1.一种基于小波包分解的燃气轮机燃烧信号分析方法,其特征在于:使用传感器多通道采集燃烧信号,对采集的燃烧信号进行滤波、放大后经FPGA模块传送至CPU处理器,在CPU处理器内搭建小波包分解算法分频段地对燃烧信号进行分析,然后将分析结果分别传送至人机界面和控制系统。
2.如权利要求1所述的基于小波包分解的燃气轮机燃烧信号分析方法,其特征在于所述小波包分解算法包括以下步骤:
(a)、根据燃烧信号中的采样频率                                               和最小频率分辨率,通过公式确定小波包分解的层数n,小波包分解后在第n层对应的节点数为
(b)、根据燃烧信号确定的分解层数n,对燃烧信号进行小波包分解,并对每层分解的节点进行编号,将确定的最大频率解析范围作为起始节点(0,0),进而分解为(1,0)和(1,1),而后继续依次将其余节点进行分解,以此类推直到(n,0)、(n,1)(n,) ;
(c)、根据燃烧信号的特征频带和最小频率分辨率,通过小波函数对燃烧信号进行分解,确定特征频带所对应的频段数目和对应编号; 
(d)、对(c)步骤所确定的特征频带和隶属于该特征频带各节点的燃烧信号进行小波重构,从而获得特征频带中每一频率分量的小波重构结果,进而求得在此特征频带上,其重构均方根幅值的最大值所对应的特征频率点;
(e)、当特征频率点所对应的幅值超过该特征频率点所在特征频带的上限,则通过通讯接口分别将特征频率点及其对应幅值传送至人机界面报警和传送至控制系统进行燃烧调整,使这一特征频率点的脉动降低至上限范围内。
3.如权利要求2所述的基于小波包分解的燃气轮机燃烧信号分析方法,其特征在于:所述步骤(a)中,当对燃烧信号进行小波包分解层数为n时,解析到的最小频率分辨率是,解析到的最大频率分辨率为
4.如权利要求2所述的基于小波包分解的燃气轮机燃烧信号分析方法,其特征在于:所述步骤(b)中对燃烧信号的分解包括0—500HZ的频率解析和500—5000HZ的频率解析,对燃烧信号0—500HZ的频率解析在最小频率分辨率=1HZ下进行,其对应的小波包分解层数为12;对燃烧信号500—5000HZ的频率解析在频率分辨率 = 10HZ下进行,其对应的小波包分解层数为9,并调换小波包分解在第9层和第12层的各个频带节点的顺序,消除频谱混叠。
5.如权利要求2所述的基于小波包分解的燃气轮机燃烧信号分析方法,其特征在于:所述步骤(c)中的小波函数为db02小波函数。
CN201410542171.7A 2014-10-15 2014-10-15 基于小波包分解的燃气轮机燃烧信号分析方法 Active CN104298879B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410542171.7A CN104298879B (zh) 2014-10-15 2014-10-15 基于小波包分解的燃气轮机燃烧信号分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410542171.7A CN104298879B (zh) 2014-10-15 2014-10-15 基于小波包分解的燃气轮机燃烧信号分析方法

Publications (2)

Publication Number Publication Date
CN104298879A true CN104298879A (zh) 2015-01-21
CN104298879B CN104298879B (zh) 2018-01-19

Family

ID=52318602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410542171.7A Active CN104298879B (zh) 2014-10-15 2014-10-15 基于小波包分解的燃气轮机燃烧信号分析方法

Country Status (1)

Country Link
CN (1) CN104298879B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106419898A (zh) * 2016-08-12 2017-02-22 武汉中旗生物医疗电子有限公司 一种去除心电信号基线漂移的方法
CN110598173A (zh) * 2019-08-31 2019-12-20 中国人民解放军陆军工程大学 基于运用效能的作战系统操作行为分析方法
CN110966100A (zh) * 2018-09-30 2020-04-07 中国航发商用航空发动机有限责任公司 一种燃烧振荡的监测装置及其方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030185247A1 (en) * 2002-03-26 2003-10-02 Hal-Wen Chen Method and system for wavelet packet transmission using a best base algorithm
CN102506444A (zh) * 2011-11-04 2012-06-20 国电南京自动化股份有限公司 基于智能控制的计算机视觉技术的炉膛火焰检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030185247A1 (en) * 2002-03-26 2003-10-02 Hal-Wen Chen Method and system for wavelet packet transmission using a best base algorithm
CN102506444A (zh) * 2011-11-04 2012-06-20 国电南京自动化股份有限公司 基于智能控制的计算机视觉技术的炉膛火焰检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张晗博 等: "基于小波变换的非平稳信号分析与处理", 《南京师范大学学报(工程技术版)》 *
王丽君 等: "小波包分解在氢发动机异常燃烧诊断中的应用", 《车用发动机》 *
赵凯 等: "基于小波包变换和功率谱估计的间谐波检测新方法", 《中国电力教育》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106419898A (zh) * 2016-08-12 2017-02-22 武汉中旗生物医疗电子有限公司 一种去除心电信号基线漂移的方法
CN110966100A (zh) * 2018-09-30 2020-04-07 中国航发商用航空发动机有限责任公司 一种燃烧振荡的监测装置及其方法
CN110966100B (zh) * 2018-09-30 2021-03-05 中国航发商用航空发动机有限责任公司 一种燃烧振荡的监测装置及其方法
CN110598173A (zh) * 2019-08-31 2019-12-20 中国人民解放军陆军工程大学 基于运用效能的作战系统操作行为分析方法
CN110598173B (zh) * 2019-08-31 2023-06-09 中国人民解放军陆军工程大学 基于运用效能的作战系统操作行为分析方法

Also Published As

Publication number Publication date
CN104298879B (zh) 2018-01-19

Similar Documents

Publication Publication Date Title
CN107065657B (zh) 机器状态无线监测设备、方法和系统
CN108168891B (zh) 滚动轴承微弱故障信号特征的提取方法及设备
CN105651376B (zh) 机械设备离线检测系统振动信号频谱分析和报警方法
CN103472008B (zh) 孵化前期近红外种蛋内鸡胚性别识别方法
CN106640722B (zh) 一种压气机气动稳定性诊断和控制的装置及方法
CN103743470B (zh) 一种汽车噪声频谱分析方法
WO2017143649A1 (zh) 基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统
CN104298879B (zh) 基于小波包分解的燃气轮机燃烧信号分析方法
CN102840907B (zh) 早期故障状态下滚动轴承振动信号特征提取和分析方法
JP2013542677A5 (zh)
CN109009087B (zh) 一种心电信号r波的快速检测方法
Shuuji et al. Low-speed bearing fault diagnosis based on improved statistical filtering and convolutional neural network
CN111160146A (zh) 基于时频转换的水电机组状态监测信号数字滤波方法、装置及系统
CN106769054A (zh) 一种基于声发射信号的水轮机组空蚀空化状态诊断方法
CN103852338A (zh) 空气净化器异常噪声故障检测方法
CN203414563U (zh) 一种局部放电ae位置检测系统
CN107610055A (zh) 傅里叶变换光谱仪干涉图噪声检测及抑制方法
CN203350268U (zh) 旋翼转速信号处理器
CN211477587U (zh) 在线安全预警装置和故障诊断系统
CN106303879A (zh) 一种基于时域分析的检测装置及检测方法
Liu et al. Research on weak signal detection for downhole acoustic telemetry system
CN103438983A (zh) 一种信号随机平均谱的数据处理方法
CN204629114U (zh) 一种应用于光纤预警系统中的自适应噪声抵消装置
CN109612573B (zh) 一种基于噪声频谱分析的树冠火和地面火检测方法
Zhang et al. Machine fault detection during transient operation using measurement denoising

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 618000 No. 18, Third Section of Lushan South Road, Deyang City, Sichuan Province

Patentee after: Dongfang Electric Automatic Control Engineering Co., Ltd.

Address before: 618099 No. 18, Section 3, Lushan South Road, Deyang City, Sichuan Province

Patentee before: Sichuan Orient Electric Automatic Control Engineering Co., Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220104

Address after: No.2828 Shuzhou Road, Zhengxing street, Tianfu New District, China (Sichuan) pilot Free Trade Zone, Chengdu, Sichuan 610000

Patentee after: Dongfang Electric Zhongneng industrial control network security technology (Chengdu) Co.,Ltd.

Address before: 618000 No. three, 18 South Mount Lu Road, Sichuan, Deyang

Patentee before: DONGFANG ELECTRIC AUTOMATIC CONTROL ENGINEERING Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220301

Address after: 618000 No. three, 18 South Mount Lu Road, Sichuan, Deyang

Patentee after: DONGFANG ELECTRIC AUTOMATIC CONTROL ENGINEERING Co.,Ltd.

Patentee after: Dongfang Electric Zhongneng industrial control network security technology (Chengdu) Co., Ltd

Address before: No.2828 Shuzhou Road, Zhengxing street, Tianfu New District, China (Sichuan) pilot Free Trade Zone, Chengdu, Sichuan 610000

Patentee before: Dongfang Electric Zhongneng industrial control network security technology (Chengdu) Co.,Ltd.