基于噪声信号A计权分析的泵站机组故障诊断方法
技术领域
本发明涉及泵站机组故障诊断技术领域,通过监测泵站机组的噪声声压信号,对声压信号进行频谱分析,对各频率成分的声压级进行A计权网络修正,并通过声压级的能量叠加得到A计权噪声声压级,当A计权噪声声压级超过报警值时对声压信号进行1/3倍频程谱分析提取各倍频程频段能量特征,利用深度极限学习机,快速有效的进行特征学习,提取各个特征的隐含故障信息,从而实现水泵机组故障的智能诊断。
背景技术
泵站机组是泵站的主要组成部件,在城市排涝、引水、调水等工程中得到广泛应用,其运行状态对整个工程和泵站的安全可靠运行起着至关重要的影响。泵站机组稳定运行是工程效益正常发挥的重要保证。水泵机组投入运行后,随着工况的变化、运行时间的延长以及一些突发事件,一些安全隐患和设备故障将逐步显现,威胁设备的安全运行,甚至导致设备遭到破坏,严重时可能引发灾难性事故。因此,实时、全面地掌握机组运行状况,及时发现存在的安全隐患和设备故障,对安全运行和工程效益的发挥都起到十分重要的作用。
噪声一方面是作为一种有害的信息需要加以克服与控制,另一方面作为机械运行时发出的一种固有信号,噪声中必然携带着机械本身的结构信息和运行的状态信息,从理论上讲完全可以利用噪声信号来对设备的运行状态及故障进行在线监测和诊断。泵站机组在运行过程中会产生一定的噪声,当机组发生故障时,噪声的频率特性和能量分布会出现不同程度的变化,据此在不停机、不解体的情况下,就可以通过对噪声信号的分析推知设备的故障部位、故障原因及其严重程度。虽然在这方面已有大量的理论研究,但是由于泵站机组本身的复杂性以及不同声源之间的相互干扰,在具体应用上尚有许多问题亟待解决,其中对采集的噪声信号的有效处理和特征提取就是其中的关键。
泵站机组由于受水力、机械、电磁等诸多因素的影响,单一故障征兆无法准确的表征故障特征,而事实上,诊断结果往往取决于提取的特征的好坏,如果提取的特征不好,整个模型很难达到很好的分类效果。其次,传统的故障诊断往往利用支持向量机、决策树等浅层分类方法,带入到机器学习的分类算法中,这种浅层学习方法往往无法自动识别各个特征的权重,从而导致诊断精度不高。近年来,随着深度学习在各个领域如语音识别,图像识别等的发展,基于深度学习理论的智能故障诊断技术有了很大突破,通过泵站机组的噪声信号的A计权分析,快速计算得到A计权噪声声压级,当A计权声压级超过报警值时通过进行1/3倍频程分析得到各倍频程声压能量值,通过深度极限学习机,快速有效的进行特征学习,提取各个特征的隐含故障信息,从而实现水泵机组故障的智能诊断。
发明内容
本发明的目的在于提供一种通过监测泵站机组的噪声声压信号,对声压信号进行频谱分析,对各频率成分的声压级进行A计权网络修正,并通过声压级的能量叠加得到A计权噪声声压级,当A计权噪声声压级超过报警值时对声压信号进行1/3倍频程谱分析提取各倍频程频段能量特征,利用深度极限学习机,快速有效的进行特征学习,提取各个特征的隐含故障信息,从而实现水泵机组故障的智能诊断的方法。
为了实现上述目的,本发明采用的技术方案如下:
一种基于噪声监测的泵站机组故障分析方法,包括如下步骤:
S1,在泵站机组工作时的每种工况下,获取距离水泵叶轮外壳和电机外壳1米处的原始噪声声压信号;
S2,分别计算水泵和电机的噪声A计权声压级,采用FFT频谱分解加权修正的方法来计算A计权声压级;
S3,判断水泵和电机的噪声A计权声压级是否超标报警,如果是则启动S4-S6的机组故障诊断分析;
S4,将噪声信号在20Hz到20kHz人耳能听到的噪声频率范围内分为31段,每个段的划分原则利用了1/3倍频程的计算方式;对噪声声压原始信号进行FFT频谱分析,在每个频程范围内将声能量进行频带能量计算,得到31个频带声压能量值;
S5,建立深度极限学习机模型,将S4得到的特征矩阵代入极限学习机模型中,采用多层自编码结构进行无监督特征学习,得到具有故障特征的特征向量,其中,所采用的多层自编码结构具有对称的输入X和输出
所述具有故障特征的特征向量为多层自编码结构最终层的隐含节点;
S6,将S5得到的所述具有故障特征的特征向量代入单层极限学习机模型中进行计算,得到最终分类结果,完成故障诊断。
优选地,S1中采用自由场传声传感器采集水泵和电机的原始噪音声压信号。
优选地,所述原始噪音声压信号中包括水泵汽蚀、叶轮刮擦、叶轮不平衡、电机转子不平衡、电机基础松动、电机磁拉力不平衡故障信号。
优选地,S2中采用FFT频谱分解加权修正的方法来计算A计权声压级。
声压级是噪声总能量大小的反映,由于人耳对声音频率的敏感性是非线性的,声压相同的噪声对人耳产生的响应也可能不同,为了反映人耳的听觉特性,通常采用计权声级作为评价参数,计权网络有A、B、C、D四种,常用情况下我们都采用A计权进行测量,因为A计权与人类的主观感知有良好的相关性,针对人耳对低频信号不那么敏感的情况进行了补偿。A计权声压级的计算方法通常为:将噪声信号在20Hz到20kHz人耳能听到的噪声频率范围内分为31段,每个段的划分原则利用了1/3倍频程的计算方式,然后在每个段内分别利用带通滤波原理将噪声信号滤波后,利用能量的特点,将声能量进行积分计算,从而得到每个频程范围内的声能量,将这些声能量进行积分计算,从而得到每个频程范围内的声能量,将这些声能量按统一的标准转化为人们熟知的声压级之后,然后将频程的声压级值连成曲线即形成了标准的1/3倍频A声级,将各频带声压级按能量求和法即可算得A计权总声压级。该常规A计权声压级的计算过程要求采用带通滤波器组实现倍频程谱,低频带20Hz左右带通滤波器阶数会高达10万,同时分数倍频程又会使滤波器个数成倍增加,如1/3倍频程滤波器要30个,导致计算量巨大。
本发明采用FFT频谱分解加权修正的方法来计算A计权声压级,具体步骤如下:
A1,根据有效值定义,噪声声压信号的有效值按如下公式计算:
公式中pe是声压有效值,表示在一定时间间隔中,瞬时声压对时间取均方根值;其中T代表取平均的时间间隔,p(t)为时间间隔内某个时刻t的瞬时声压;
其中Lp为声压级,dB;p0=2×10-5Pa,为基准声压;pe为有效声压;
A2,根据傅里叶变换理论,任何一个函数都可以分解成一系列正弦函数的叠加,对噪声声压信号通过傅里叶变换分解为若干频率的纯音源,离散傅里叶变换的定义如下:
其中x(n)为N点有限长时域序列,e为自然常数,i为虚数单位,X(k)为N点有限长频域序列;
A3,假设声场中有n个不同噪声源,第i个噪声源的声压为pei,根据不相干声波叠加原理,合成声场的总声压与各声源的关系为
因此,多频噪声总声压级可表示为
A4,根据声压级定义公式(2),总声压级可表示为各分频声压级的能量叠加,公式为
A5,A计权网络是以40方等响曲线倒置,并使1000Hz值位于0时对噪声进行加权分析。在国标GB/T3785.1-2010中,给出了A计权网络曲线公式,如下:
其中f表示要计算的频率,f1=20.6Hz,f2=107.7Hz,f3=737.9Hz,f4=12194Hz;A1000是以分贝表示的归一化常数,相应于在1kHz提供0dB频率计权所需的增益,A1000=-2dB;
A6,对噪声信号的傅里叶变换频谱,对每个频谱幅值点,可按公式(7)计算该频率处的修正A计权分贝值,1000Hz处,修正分贝值为0dB,100Hz处,修正分贝值为-19.1dB;对各频谱幅值点计算声压级,对各频率幅值点的声压级进行A计权修正,然后合成总的A计权声压级,计算公式为:
其中Lpi是各频谱幅值所在频率处的声压级,取该频率处的有效值幅值,按公式(2)计算;Afi是该频率处的修正A计权分贝值,按公式(7)计算;
本发明的有益效果是:本发明实施例提供的一种基于噪声信号A计权分析的泵站机组故障诊断方法,通过监测泵站机组的噪声声压信号,对声压信号进行频谱分析,对各频率成分的声压级进行A计权网络修正,并通过声压级的能量叠加得到A计权噪声声压级,当A计权噪声声压级超过报警值时对声压信号进行1/3倍频程谱分析提取各倍频程频段能量特征,利用深度极限学习机,快速有效的进行特征学习,提取各个特征的隐含故障信息,从而实现水泵机组故障的智能诊断。
附图说明
图1是本发明提供的泵站机组故障诊断方法流程示意图;
图2是本发明提供的A计权网络曲线频率修正图;
图3是本发明提供的采用FFT频谱分解加权修正的方法来计算A计权声压级流程示意图;
图4是本发明提供的深度极限学习及诊断算法流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
如图1-4所示,本发明实施例提供了一种基于噪声信号A计权分析的泵站机组故障诊断方法,包括如下步骤:
S1,在泵站机组工作时的每种工况下,获取距离水泵叶轮外壳和电机外壳1米处的原始噪声声压信号;
S2,分别计算水泵和电机的噪声A计权声压级,采用FFT频谱分解加权修正的方法来计算A计权声压级;
S3,判断水泵和电机的噪声A计权声压级是否超标报警,如果是则启动S4-S6的机组故障诊断分析;
S4,将噪声信号在20Hz到20kHz人耳能听到的噪声频率范围内分为31段,每个段的划分原则利用了1/3倍频程的计算方式;对噪声声压原始信号进行FFT频谱分析,在每个频程范围内将声能量进行频带能量计算,得到31个频带声压能量值;
S5,建立深度极限学习机模型,将S4得到的特征矩阵代入极限学习机模型中,采用多层自编码结构进行无监督特征学习,得到具有故障特征的特征向量,其中,所采用的多层自编码结构具有对称的输入X和输出
所述具有故障特征的特征向量为多层自编码结构最终层的隐含节点;
S6,将S5得到的所述具有故障特征的特征向量代入单层极限学习机模型中进行计算,得到最终分类结果,完成故障诊断。
其中,S1中采用自由场传声传感器采集水泵和电机的原始噪声声压信号。
所述原始噪声声压信号中包括水泵汽蚀、叶轮刮擦、叶轮不平衡、电机转子不平衡、电机基础松动、电机磁拉力不平衡故障信号。
本发明的一个优选实施例中,S2中采用FFT频谱分解加权修正的方法来计算A计权声压级,可参见图3,具体步骤如下:
A1,根据有效值定义,噪声声压信号的有效值按如下公式计算:
公式中pe是声压有效值,表示在一定时间间隔中,瞬时声压对时间取均方根值;其中T代表取平均的时间间隔,p(t)为时间间隔内某个时刻t的瞬时声压;
其中Lp为声压级,dB;p0=2×10-5Pa,为基准声压;pe为有效声压;
A2,根据傅里叶变换理论,任何一个函数都可以分解成一系列正弦函数的叠加,对噪声声压信号通过傅里叶变换分解为若干频率的纯音源,离散傅里叶变换的定义如下:
其中x(n)为N点有限长时域序列,e为自然常数,i为虚数单位,X(k)为N点有限长频域序列;
A3,假设声场中有n个不同噪声源,第i个噪声源的声压为pei,根据不相干声波叠加原理,合成声场的总声压与各声源的关系为
因此,多频噪声总声压级可表示为
A4,根据声压级定义公式(2),总声压级可表示为各分频声压级的能量叠加,公式为
A5,A计权网络是以40方等响曲线倒置,并使1000Hz值位于0时对噪声进行加权分析。在国标GB/T3785.1-2010中,给出了A计权网络曲线公式,如下:
其中f表示要计算的频率,f1=20.6Hz,f2=107.7Hz,f3=737.9Hz,f4=12194Hz;A1000是以分贝表示的归一化常数,相应于在1kHz提供0dB频率计权所需的增益,A1000=-2dB;
A6,对噪声信号的傅里叶变换频谱,对每个频谱幅值点,可按公式(7)计算该频率处的修正A计权分贝值,1000Hz处,修正分贝值为0dB,100Hz处,修正分贝值为-19.1dB;对各频谱幅值点计算声压级,对各频率幅值点的声压级进行A计权修正,然后合成总的A计权声压级,计算公式为:
其中Lpi是各频谱幅值所在频率处的声压级,取该频率处的有效值幅值,按公式(2)计算;Afi是该频率处的修正A计权分贝值,按公式(7)计算;
本发明实施例提供的泵站机组故障诊断方法,至少具有以下优点:
1、本发明基于噪声信号A计权分析的泵站机组故障诊断方法中,通过监测泵站机组的噪声声压信号,对声压信号进行频谱分析,对各频率成分的声压级进行A计权网络修正,并通过声压级的能量叠加得到A计权噪声声压级。此方法摆脱了传统A计权声压级的获取需要对原始噪声声压信号的各倍频程频段进行计算量巨大的带通滤波后计算每个频程声压级最终按能量求和法得到总声压级的弊端,也避免了直接通过硬件A计权滤波网络获得A计权声压级而无法获取原始声压信号进行故障诊断的弊端,兼顾了获得噪声信号强度和对噪声作频谱分析的需求。
2、本发明基于噪声信号A计权分析的泵站机组故障诊断方法中,采用了对声压信号进行1/3倍频程谱分析提取各倍频程频段能量特征,将各倍频程频段能量作为输入,利用深度极限学习机,快速有效的进行特征学习,提取各个特征的隐含故障信息,从而实现水泵机组故障的智能诊断。泵站机组的原始噪声声压信号中包括水泵汽蚀、叶轮刮擦、叶轮不平衡、电机转子不平衡、电机基础松动、电机磁拉力不平衡等故障信号,各种故障信号会表现为各倍频程频段能量的变化,通过深度极限学习机对特征进行学习,可以很好地区分各种水泵机组故障。
3、与现有的故障诊断技术相比,本发明的泵站机组智能故障诊断方法能够直接从泵站机组噪声信号中获得噪声强度进行故障报警,并从噪声信号各倍频程频段能量中通过深度极限学习机提取特征,准确地识别泵站机组的故障类型,为解决泵站机组的故障诊断问题提供了一种新的有效途径,可广泛应用于电力、机械、冶金、化工等各重要领域的复杂系统当中。
具体实施例
本发明实施例提供的泵站机组故障诊断方法,可以采用如下方法进行实施:
步骤1、在泵站机组上的水泵和电机附近距离机组设备1米处安装自由场传声传感器,采集水泵和电机的原始噪声声压信号。
步骤2、分别计算水泵和电机的噪声A计权声压级,采用FFT频谱分解加权修正的方法来计算A计权声压级。
步骤3、当A计权声压级报警时,收集故障样本数据,即水泵汽蚀、叶轮刮擦、叶轮不平衡、电机转子不平衡、电机基础松动、电机磁拉力不平衡,对每种故障类型添加标签,建立训练样本集。
步骤4、特征提取阶段,将噪声信号在20Hz到20kHz人耳能听到的噪声频率范围内分为31段,每个段的划分原则利用了1/3倍频程的计算方式;对噪声声压原始信号进行FFT频谱分析,在每个频程范围内将声能量进行频带能量计算,得到31个频带声压能量值,作为极限学习机的输入特征向量。
步骤5、故障诊断阶段,将极限学习机分为6层,其中分别分为输入层、隐含层1、隐含层2、隐含层3,隐含层4、输出层,其中输入层数目为特征向量个数、输出层为故障类型数目,隐含层数目分别为100、200、400、800,用于稀疏特征。代入训练数据集,采用随机权重进行模型训练。
步骤6、上线运行阶段,模型训练好之后,保存模型结构和各个层的权重和偏置,部署在实际机组监测系统上,当水泵或电机的噪声A计权声压级报警时,启动诊断算法,重复特征提取阶段,并代入训练好的模型进行分类识别,得到的分类结果就是判断的故障类型。
通过采用本发明公开的上述技术方案,得到了如下有益的效果:本发明实施例提供的一种基于噪声信号A计权分析的泵站机组故障诊断方法,通过监测泵站机组的噪声声压信号,对声压信号进行频谱分析,对各频率成分的声压级进行A计权网络修正,并通过声压级的能量叠加得到A计权噪声声压级,当A计权噪声声压级超过报警值时对声压信号进行1/3倍频程谱分析提取各倍频程频段能量特征,利用深度极限学习机,快速有效的进行特征学习,提取各个特征的隐含故障信息,从而实现水泵机组故障的智能诊断。此方法摆脱了传统A计权声压级的获取需要对原始噪声声压信号的各倍频程频段进行计算量巨大的带通滤波后计算每个频程声压级最终按能量求和法得到总声压级的弊端,也避免了直接通过硬件A计权滤波网络获得A计权声压级而无法获取原始声压信号进行故障诊断的弊端,兼顾了获得噪声信号强度和对噪声作频谱分析的需求。本发明的泵站机组智能故障诊断方法能够直接从泵站机组噪声信号中获得噪声强度进行故障报警,并从噪声信号各倍频程频段能量中通过深度极限学习机提取特征,准确地识别泵站机组的故障类型,为解决泵站机组的故障诊断问题提供了一种新的有效途径。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。