专利名称：故障预测模型建立方法及系统、风机监测预警系统及方法
技术领域：
本发明涉及一种故障预测模型建立方法及系统、风机监测预警系统及方法，特别是涉及一种以音频信号辅助振动信号的故障预测模型建立方法及系统、风机监测预警系统及方法。
背景技术：
风力发电厂的所有人越来越感到运营成本的压力，迫切希望降低运营成本；尤其在当今时代更是如此。要实现运营成本的降低，必须更有效的使用维护和维修服务的人员。当风力发电机需要维修和维护时，就需要聘请维修人员，租借大型吊车，以及采购更换的零部件。通过使用故障诊断系统，风力发电厂的所有人处于成本缩减的考虑，可以将多次的维修合并到一起；例如，一个风力发电机已经因故障停机，而另一台风力发电机(简称“风机”) 已经出现一些部件即将发生故障的信号；那么，就可以将两个风机的维修安排合并到一起来完成以此来降低维修成本。通过追踪风力发电机的实效模式，可以延长风机的维修间隔。通过追踪风机的状态，风力发电厂的所有人可以预计风机继续运行的时间，从而延长风机的服务和维修间隔。通过将故障诊断系统集成到风力发电机上，风力发电机的操作者可以实现对风力发电机的主动性维护。换言之，通过仔细的监控，他们可以预计何时对风力发电机实施必要的维护。这使得他们可以更准确的和更有效的安排维护工作。在风机故障中，滚动轴承、齿轮箱为容易出现故障的零部件，而滚动轴承、齿轮箱的高频故障都具有如下特点故障初期振动冲击明显，但随着故障程度的加深，振动冲击逐渐恢复至正常值，而振动能量值却明显增大。冲击信号虽然敏感但不够稳定；振动能量值不够敏感，但会随着故障的稳定而增长。这就要求在对这类故障进行诊断时，既要关注高频冲击值，又要关注低频有效值。当前，风机的在线监测和故障诊断系统都基于对风机轴承振动的监测，对振动的加速度传感器信号值进行时域波形分析、概率密度分析、自相关分析、趋势分析、幅值谱分析、功率谱分析、倒频谱分析、包络谱分析、细化谱分析、短时傅里叶分析、小波分析等等。然而，为了得到更精确、更灵敏的监测效果以对风机进行更好的维护并且在故障发生前及时维护、更换容易出现故障的零部件，期望找到一种灵敏度更高、故障监测更准确、预测数据更准确的监测预警系统及方法。

发明内容
本发明要解决的技术问题是为了克服现有技术中对风机故障的监测方式较为单一、监测结果灵敏度不足、故障预警准确度不高的缺陷，提供一种灵敏度更高、故障监测更准确、预测数据更准确的故障预测模型建立方法及系统、风机监测预警系统及方法。在实现风机的监测预警之前，需要先建立较为精准、全面的故障类型预测模型和故障概率预测模型，之后在根据建立的故障类型预测模型和故障概率预测模型再来实现风机的监测预警。在现有技术中，对于风机的在线监测都基于振动信号，从而导致了监测结果灵敏度不足、故障预警准确度不高的缺陷，本发明在对风机的多种故障深入研究、大量实验的基础上引入了音频信号监测的方式来实现风机的监测预警。本发明是通过下述技术方案来解决上述技术问题的一种故障类型预测模型的建立方法，其特点在于，包括以下步骤步骤S1、采集风机的振动信号，以及在采集该振动信号的同时采集风机的音频信号;步骤S2、预处理该振动信号并对预处理后的振动信号进行特征提取以获得振动特征参数，以及预处理该音频信号并对预处理后的音频信号进行特征提取以获得音频特征参数；步骤S3、将该振动特征参数代入一故障类型数据库中的故障类型振动模板以得到与该振动特征参数相对应的故障类型，该故障类型振动模板表示故障类型与振动特征参数 之间的对应关系；步骤S4、将该音频特征参数与该故障类型相关联以建立故障类型音频模板，该故障类型音频模板表示故障类型与音频特征参数之间的对应关系，多个故障类型音频模板构成该故障类型预测模型。鉴于音频信号具有携带信息丰富和非接触测量的特点，可对移动目标进行检测，可以采用音频信号来建立精准、全面的故障类型预测模型。在该故障类型预测模型建立之初，可以采用该故障类型数据库中根据历史振动数据建立的故障类型振动模板先将与该振动特征参数相对应的故障类型确定下来，由于音频信号和振动信号是同时采集的，它们之间有着唯一确定的关系，继而可以通过该唯一确定的关系将该音频特征参数与该故障类型相关联，并由此建立故障类型音频模板，多个故障类型音频模板构成该故障类型预测模型。这样，基于音频数据的故障类型预测模型就建立完成了。当然，故障类型预测模型建立之后还可以在采集了更丰富的振动数据、音频数据后不断对该故障类型预测模型进行完善。当该故障类型预测模型由足够多的数据完善到一定程度时，甚至可以完全基于音频特征参数来实现风机的监测预警。首先，对于振动信号来说，可以采用现有的预处理和特征提取方法来处理采集的振动信号，例如振动信号的采样按照现有的一维信号进行时频分析处理，提取信号处理的特征参数。其中，时频分析的方法，它的思路是设计时间和频率的联合函数，用它同时描述信号在不同时间和频率的能量密度或强度。利用时频分布来分析信号，能给出各个时刻的瞬时频率及其幅值，并且能够进行时频滤波和时变信号研究。时频分析的方法包括时域波形分析方法、概率密度分析方法、自相关分析方法、趋势分析方法、幅值谱分析方法、功率谱分析方法、倒频谱分析方法、包络谱分析方法、细化谱分析方法、短时傅里叶分析方法、小波分析方法。以上只是举例说明，对于振动特征参数的获得不局限于上述方法，还可以采用其他的现有手段实现。其次，对于音频信号来说优选地，步骤S2中预处理该音频信号还包括以下步骤步骤S21、对该音频信号分帧处理；步骤S22、预加重经分帧处理的该音频信号；步骤S23、加窗处理经预加重的该音频信号。
更优选地,步骤S21中将该首频信号分为每10ms-30ms —巾贞,且相邻两巾贞有的重叠；和/或，步骤S22中通过一预加重数字滤波器来实现预加重；和/或，步骤S23中通过汉明窗来实现加窗处理。即在对音频信号进行特征参数提取之前，先要对音频信号进行采样和预处理，预处理一般包括分帧处理、预加重和加窗处理三个部分。例如由一模数转换器以合适的采样频率对模拟的音频信号进行采样，再以适合的量化精度进行量化，得到计算机能够处理的离散的数字信号。为了防止混叠效应，采样频 率必须满足奈奎斯特定理。音频信号的时域波形显示出很强的时变特性，需要对音频信号进行分帧处理之后再进行特征提取。音频信号时一种典型的非平稳信号。因此音频信号可假定为短时平稳的，即在一时间段内，例如10-20ms的时间段内该音频信号是平稳的，其频谱特性和某些物理特性参量可近似地看作是不变的，因此在该时间段内可以采用平稳过程的分析处理方法来处理。也就是说，将音频信号进行分帧操作，将该音频信号分成多帧，例如分为每10ms-30ms —段，且相邻两帧之间有1/3到2/3的重叠。由于采样频率的差异，帧长和帧移所对应的实际采样点数也是不同的。假设所采用的音频信号的采样频率为22. 05kHz, 11. 6ms的帧长对应256个点，5. 8ms的帧移对应128个点。帧长的样点数=采样率*窗时间长，N= 22050*0.0116 = 255. 78，这个窗时间长取值为11. 6ms，N不是整数，一般分帧的时候取到点数为整数，用整段输入语音(整段输入语音即为采样的音频信号)除上N取整就是分巾贞的个数,最后一巾贞不够N大小的样点可以不做处理。Overlap-add (重叠相加法)就是解决语音分帧后产生不连续的问题。采集器采集的肯定是一帧一帧的。音频信号的高频端大约在800Hz以上按6dB/倍频程跌落，高频的音频信号频谱都很小，通常是频率越高幅值越小，为此要对语音信号进行预加重。预加重的目的是提升高频部分，使音频信号的频谱变得平坦，以便于进行频谱分析或声道参数分析。预加重通过两种渠道实现一种可以放在A/D (模拟/数字)变换前在反混叠滤波之前进行，这样不仅能够进行预加重，而且压缩了信号的动态范围，有效地提高了信噪比。另一种可以放在A/D变换之后进行，用具有倍频程的提升高频特性的预加重数字滤波器来实现。它一般是一阶的H(Z)=I-U z—1，其中ii值接近于I。在实际中选取U等于0.9375的滤波器来对音频信号进行预加重处理。音频信号也是一个一维数组，就是将音频信号代入上述公式中，得到一个新的信号，就是预加重处理了。为了保持语音信号(语音信号即为经分帧和预加重后的音频信号)的短时平稳性，利用窗函数来减少由截断处理导致的Gibbs效应(吉布斯现象Gibbs phenomenon,又叫吉布斯效应将具有不连续点的周期函数(如矩形脉冲)进行Fourier series (傅立叶级数)展开后，选取有限项进行合成。当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时，该峰起值趋于一个常数，大约等于总跳变值的9%。这种现象称为吉布斯现象)。理想的窗函数的频率响应要求主瓣无限狭窄且没有旁瓣(无频谱泄露)，但这种窗函数在实际工程中无法实现。所以，根据不同的应用我们可以采用不同的窗函数来逼近理想的频率响应。把音频信号看成在短时段上是一个平稳过程。可以用对平稳过程的分析方法进行分析，因此需要将语音信号划分为一个一个的短时段，每一短时段称为一巾贞，为从语音信号中切取含有N个样本的音频信号波形，需要用时间窗函数乘以原来的音频信号。而矩形窗的采用使得每一帧信号在起始点和终点处发生突变，不再连续，从而导致Gibbs现象的出现。为了减小起始点和终点处的不连续性，可以改变窗函数的时域特性，使其在起始点和终点处逐渐变为零，如Hamming窗(汉明窗)。优选地，步骤S2中对预处理后的音频信号进行特征提取的方法选自线性预测倒谱系数方法、美尔频率倒谱系数方法和小波包分析方法。具体来说,线性预测倒谱系数LPCC( linear forecast cepstrum coefficient)利用音频信号采样点之间的相关性，用过去的采样点的值来预测现在或者未来采样点的值。一个音频信号的抽样能够用过去若干个音频抽样或者他们的线性组合来逼近。通过使实际的音频信号抽样值与线性预测抽样值的均方误差达到最小，这样就能确定出唯一的一组线性预测系数。 而美尔频率倒谱系数MFCC (Mel Frequency Cepstrum Coeff icient)是基于人的听觉特性，在IOOOHz以下，感知能力与频率成线性关系；在1000Hz以上，感知能力与频率成对数关系。为了模拟这种人耳的感知特征，人们提出了频标的概念，意义为=IMel为1000Hz的音调感知程度的1/1000。小波包分析的音频信号特征提取包括(1)将采集到得音频信号进行小波包分解，音频信号经小波包分解可得到其在各频带内成分的特征，即各低频和高频系数。(2)重构小波包分解系数，求取各频带的信号能量，并进行归一化处理。即将各个小波包分解系数做选择性重构，抽取和故障相关的几个频带进行重构，提取各频段范围的信号，并计算各频段的信号能量。(3)构造特征向量。由于出现故障对各频带内的信号能量有较大的影响，因此能量变化突出的若干信号为特征向量。本系统对音频信号进行db3小波分解，对得到的第3层小波包分解系数进行重构。(4)对得到的各频带的系数求取能量值并进行归一化处理，得到各类音频信号的离散小波变换系数。同样地，以上只是举例说明，对于音频特征参数的获得不局限于上述方法，还可以采用其他的现有手段实现。优选地，步骤S3中该故障类型振动模板通过采用该故障类型数据库中的历史振动数据和BP算法获得。BP算法，即采用BP神经网络模型，BP神经网路模型把一组样本的输入/输出问题转化为一个非线性优化问题，使用了优化中最普遍的梯度下降法。用迭代运算求解全职的过程相应于学习记忆问题。这种网络训练过程是动态的，前向计算则是输入到输出的一种高度非线性映射。BP神经网络对于相同的输入一定得到相同的输出，对于相近的输入则可以得到相近的输出。这一特性进行推理，专家知识作为训练样本，经过学习后被存储在网络的权值中，推理过程对输入模式作前向计算，可得出推理结论。通过该故障类型数据库中的历史振动数据和BP算法建立起的该故障类型振动模板，就能先通过该振动特征参数确定该故障类型，以便将该音频特征参数与该故障类型相关联。优选地，步骤S1中采集的为该风机的机舱中的振动信号和音频信号以及该风机的传动系统的振动信号和音频信号。优选地，步骤S1中采集的该振动信号和该音频信号为模拟信号，步骤S1之后、步骤S2之前还包括以下步骤将该模拟信号转换为数字信号。
本发明还提供一种故障类型预测模型的建立系统，其特点在于，其包括一振动传感装置，用于采集风机的振动信号；一音频传感装置，用于在采集该振动信号的同时采集风机的音频信号；一第一数据处理装置，用于预处理 该振动信号并对预处理后的振动信号进行特征提取以获得振动特征参数；一第二数据处理装置，用于预处理该音频信号并对预处理后的音频信号进行特征提取以获得音频特征参数；一故障类型分析装置，用于将该振动特征参数代入一故障类型数据库中的故障类型振动模板以得到与该振动特征参数相对应的故障类型，该故障类型振动模板表示故障类型与振动特征参数之间的对应关系；一故障类型预测模型建立装置，用于将该音频特征参数与该故障类型相关联以建立故障类型音频模板，该故障类型音频模板表示故障类型与音频特征参数之间的对应关系，多个故障类型音频模板构成该故障类型预测模型。优选地，该故障类型振动模板是通过采用该故障类型数据库中的历史振动数据和BP算法获得的。优选地，该第一数据处理装置用于通过时频分析方法来预处理该振动信号；和/或，该第二数据处理装置还包括一分帧处理装置，用于对该音频信号分帧处理；—预加重装置，用于预加重经分帧处理的该音频信号；一加窗处理装置，用于加窗处理经预加重的该音频信号。优选地，时频分析方法选自时域波形分析方法、概率密度分析方法、自相关分析方法、趋势分析方法、幅值谱分析方法、功率谱分析方法、倒频谱分析方法、包络谱分析方法、细化谱分析方法、短时傅里叶分析方法和小波分析方法；和/或，该分巾贞处理装置用于将该音频信号分为每10ms-30ms —巾贞,且相邻两巾贞有了的
3 3
重叠；和/或，该预加重装置为一预加重数字滤波器；和/或，该加窗处理装置包括一汉明窗。优选地，第二数据处理装置用于通过以下方法中的一种或多种来对预处理后的音频信号进行特征提取该线性预测倒谱系数方法、美尔频率倒谱系数方法和小波包分析方法。优选地，该振动传感装置用于采集该风机的机舱中的振动信号和该风机的传动系统的振动信号，和/或，该音频传感装置用于采集该风机的机舱中的音频信号和该风机的传动系统和音频信号。优选地，该振动传感装置采集的该振动信号和该音频传感装置采集的该音频信号为模拟信号，该振动传感装置还包括一第一 A/D (模拟/数字)转换器，该音频传感装置还包括一第二 A/D转换器,该第一 A/D转换器和该第二 A/D转换器均用于将该模拟信号转换为数字信号。本发明还提供一种故障概率预测模型的建立方法，其特点在于，包括以下步骤
步骤Sn、采集风机的振动信号，以及在采集该振动信号的同时采集风机的音频信号;步骤S2t、预处理该振动信号并对预处理后的振动信号进行特征提取以获得振动特征参数，以及预处理该音频信号并对预处理后的音频信号进行特征提取以获得音频特征参数；步骤S3t、将该振动特征参数代入一故障类型数据库中的故障类型振动模板以得到与该振动特征参数相对应的故障类型，该故障类型振动模板表示故障类型与振动特征参数之间的对应关系；步骤S4t、将该音频特征参数与该故障类型相关联以建立故障类型音频模板，该故障类型音频模板表示故障类型与音频特征参数之间的对应关系，多个故障类型音频模板构 成故障类型预测模型；步骤S5t、采用该音频特征参数和HMM算法为每个故障类型建立故障概率预测模型。在确定了故障类型之后，针对每种故障类型，建立起HMM模型以便之后预测发生该种故障的概率，以及时提示用户实施相应的应对措施。隐马尔可夫模型HMM (HiddenMarkov Model)是一种描述随机过程统计特性的概率模型,它是在Markov链的基础上发展起来的，能够对多个观察样本进行有效融合而构成一个模型。HMM训练是指从同类故障的大量信号样本中提取统计信息，利用恰当的训练算法对模型参数反复修正直至收敛，最后得到模型的状态转移概率、观察概率分布和初始概率分布等参数。为了增加HMM故障诊断系统的稳健性和提高故障诊断的准确率，需要选取多个样本进行训练，以建立各类故障状态的HMM参数模型。记L各观察值序列(即L个样本)为Oa)、0(2)、…、0(1)，每个观察值序列的长度记为 T，则 Ow = IO11, O21,0/}，1=1，2, ...,L0优选地，步骤S3t中该故障类型振动模板通过采用该故障类型数据库中的历史振动数据和BP算法获得。优选地，步骤S2t中通过时频分析方法来预处理该振动信号；和/或，步骤S2t中预处理该音频信号还包括以下步骤步骤Sm、对该音频信号分帧处理；步骤S2t2、预加重经分帧处理的该音频信号；步骤S2t3、加窗处理经预加重的该音频信号。优选地，时频分析方法选自时域波形分析方法、概率密度分析方法、自相关分析方法、趋势分析方法、幅值谱分析方法、功率谱分析方法、倒频谱分析方法、包络谱分析方法、细化谱分析方法、短时傅里叶分析方法和小波分析方法；和/或，步骤Sm中将该音频信号分为每10ms-30ms —巾贞,且相邻两巾贞有"7的重叠；和丨
3 3
或，步骤S2t2中通过一预加重数字滤波器来实现预加重；和/或，步骤S2t3中通过汉明窗来实现加窗处理。优选地，步骤S2t中对预处理后的音频信号进行特征提取的方法选自线性预测倒谱系数方法、美尔频率倒谱系数方法和小波包分析方法。
优选地，步骤Sit中采集的为该风机的机舱中的振动信号和音频信号以及该风机的传动系统的振动信号和音频信号。优选地，步骤Sit中采集的该振动信号和该音频信号为模拟信号，步骤Sit之后、步骤S2t之前还包括以下步骤将该模拟信号转换为数字信号。优选地，步骤S5t中包括以下步骤采用该音频特征参数和HMM算法计算以下参数以为每个故障类型建立故障概率预测模型，该参数包括初始状态概率、状态转移概率矩阵、混合高斯函数权值、均值和协方差矩阵。具体来说，HMM模型训练是指从同类故障的大量音频信号样本中提取统计信息，利用恰当的训练算法对模型参数反复修正直至收敛，最后得到模型的初始状态概率n i、状态转移概
权利要求
1.一种故障类型预测模型的建立方法，其特征在于，包括以下步骤步骤S1、采集风机的振动信号，以及在采集该振动信号的同时采集风机的音频信号；步骤S2、预处理该振动信号并对预处理后的振动信号进行特征提取以获得振动特征参数，以及预处理该音频信号并对预处理后的音频信号进行特征提取以获得音频特征参数；步骤S3、将该振动特征参数代入一故障类型数据库中的故障类型振动模板以得到与该振动特征参数相对应的故障类型，该故障类型振动模板表示故障类型与振动特征参数之间的对应关系； 步骤S4、将该音频特征参数与该故障类型相关联以建立故障类型音频模板，该故障类型音频模板表示故障类型与音频特征参数之间的对应关系，多个故障类型音频模板构成该故障类型预测模型。
2.如权利要求I所述的故障类型预测模型的建立方法，其特征在于，步骤S3中该故障类型振动模板通过采用该故障类型数据库中的历史振动数据和BP算法获得。
3.如权利要求I所述的故障类型预测模型的建立方法，其特征在于，步骤S2中通过时频分析方法来预处理该振动信号；和/或， 步骤S2中预处理该音频信号还包括以下步骤 步骤S21、对该音频信号分帧处理； 步骤S22、预加重经分帧处理的该音频信号； 步骤S23、加窗处理经预加重的该音频信号。
4.如权利要求3所述的故障类型预测模型的建立方法，其特征在于，时频分析方法选自时域波形分析方法、概率密度分析方法、自相关分析方法、趋势分析方法、幅值谱分析方法、功率谱分析方法、倒频谱分析方法、包络谱分析方法、细化谱分析方法、短时傅里叶分析方法和小波分析方法； 步骤S21中将该音频信号分为每10ms-30ms —帧，且相邻两帧有的重叠；和/或， j 3 步骤S22中通过一预加重数字滤波器来实现预加重；和/或， 步骤S23中通过汉明窗来实现加窗处理。
5.如权利要求I所述的故障类型预测模型的建立方法，其特征在于，步骤S2中对预处理后的音频信号进行特征提取的方法选自线性预测倒谱系数方法、美尔频率倒谱系数方法和小波包分析方法。
6.如权利要求1-5中任意一项所述的故障类型预测模型的建立方法，其特征在于，步骤S1中采集的为该风机的机舱中的振动信号和音频信号以及该风机的传动系统的振动信号和音频信号。
7.如权利要求6所述的故障类型预测模型的建立方法，其特征在于，步骤S1中采集的该振动信号和该音频信号为模拟信号，步骤S1之后、步骤S2之前还包括以下步骤将该模拟信号转换为数字信号。
8.一种故障类型预测模型的建立系统，其特征在于，其包括 一振动传感装置，用于采集风机的振动信号； 一音频传感装置，用于在采集该振动信号的同时采集风机的音频信号； 一第一数据处理装置，用于预处理该振动信号并对预处理后的振动信号进行特征提取以获得振动特征参数； 一第二数据处理装置，用于预处理该音频信号并对预处理后的音频信号进行特征提取以获得音频特征参数； 一故障类型分析装置，用于将该振动特征参数代入一故障类型数据库中的故障类型振动模板以得到与该振动特征参数相对应的故障类型，该故障类型振动模板表示故障类型与振动特征参数之间的对应关系； 一故障类型预测模型建立装置，用于将该音频特征参数与该故障类型相关联以建立故障类型音频模板，该故障类型音频模板表示故障类型与音频特征参数之间的对应关系，多个故障类型音频模板构成该故障类型预测模型。
9.如权利要求8所述的故障类型预测模型的建立系统，其特征在于，该故障类型振动模板是通过采用该故障类型数据库中的历史振动数据和BP算法获得的。
10.如权利要求8所述的故障类型预测模型的建立系统，其特征在于，该第一数据处理装置用于通过时频分析方法来预处理该振动信号；和/或，该第二数据处理装置还包括 一分帧处理装置，用于对该音频信号分帧处理； 一预加重装置，用于预加重经分帧处理的该音频信号； 一加窗处理装置，用于加窗处理经预加重的该音频信号。
11.如权利要求10所述的故障类型预测模型的建立系统，其特征在于，时频分析方法选自时域波形分析方法、概率密度分析方法、自相关分析方法、趋势分析方法、幅值谱分析方法、功率谱分析方法、倒频谱分析方法、包络谱分析方法、细化谱分析方法、短时傅里叶分析方法和小波分析方法；和/或， 该分帧处理装置用于将该音频信号分为每10ms-30ms —帧，且相邻两帧有的重 J 3置；和/或， 该预加重装置为一预加重数字滤波器；和/或， 该加窗处理装置包括一汉明窗。
12.如权利要求8所述的故障类型预测模型的建立系统，其特征在于，第二数据处理装置用于通过以下方法中的一种或多种来对预处理后的音频信号进行特征提取该线性预测倒谱系数方法、美尔频率倒谱系数方法和小波包分析方法。
13.如权利要求8-12中任意一项所述的故障类型预测模型的建立系统，其特征在于，该振动传感装置用于采集该风机的机舱中的振动信号和该风机的传动系统的振动信号，和/或，该音频传感装置用于采集该风机的机舱中的音频信号和该风机的传动系统和音频信号。
14.如权利要求13所述的故障类型预测模型的建立系统，其特征在于，该振动传感装置采集的该振动信号和该音频传感装置采集的该音频信号为模拟信号，该振动传感装置还包括一第一 A/D转换器,该音频传感装置还包括一第二 A/D转换器,该第一 A/D转换器和该第二 A/D转换器均用于将该模拟信号转换为数字信号。
15.一种故障概率预测模型的建立方法，其特征在于，包括以下步骤 步骤Sit、采集风机的振动信号，以及在采集该振动信号的同时采集风机的音频信号； 步骤S2t、预处理该振动信号并对预处理后的振动信号进行特征提取以获得振动特征参数，以及预处理该音频信号并对预处理后的音频信号进行特征提取以获得音频特征参数； 步骤S3t、将该振动特征参数代入一故障类型数据库中的故障类型振动模板以得到与该振动特征参数相对应的故障类型，该故障类型振动模板表示故障类型与振动特征参数之间的对应关系； 步骤S4t、将该音频特征参数与该故障类型相关联以建立故障类型音频模板，该故障类型音频模板表示故障类型与音频特征参数之间的对应关系，多个故障类型音频模板构成故障类型预测模型； 步骤S5t、采用该音频特征参数和HMM算法为每个故障类型建立故障概率预测模型。
16.如权利要求15所述的故障概率预测模型的建立方法，其特征在于，步骤S3t中该故障类型振动模板通过采用该故障类型数据库中的历史振动数据和BP算法获得。
17.如权利要求15所述的故障概率预测模型的建立方法，其特征在于，步骤S2t中通过时频分析方法来预处理该振动信号；和/或， 步骤S2t中预处理该音频信号还包括以下步骤 步骤Sm、对该音频信号分帧处理； 步骤S2t2、预加重经分帧处理的该音频信号； 步骤S2t3、加窗处理经预加重的该音频信号。
18.如权利要求17所述的故障概率预测模型的建立方法，其特征在于，时频分析方法选自时域波形分析方法、概率密度分析方法、自相关分析方法、趋势分析方法、幅值谱分析方法、功率谱分析方法、倒频谱分析方法、包络谱分析方法、细化谱分析方法、短时傅里叶分析方法和小波分析方法；和/或，步骤Sm中将该音频信号分为每10ms-30ms —帧，且相邻两帧有的重叠；和/或， 3 J 步骤S2t2中通过一预加重数字滤波器来实现预加重；和/或， 步骤S2t3中通过汉明窗来实现加窗处理。
19.如权利要求15所述的故障概率预测模型的建立方法，其特征在于，步骤S2t中对预处理后的音频信号进行特征提取的方法选自线性预测倒谱系数方法、美尔频率倒谱系数方法和小波包分析方法。
20.如权利要求15-19中任意一项所述的故障概率预测模型的建立方法，其特征在于，步骤Sit中采集的为该风机的机舱中的振动信号和音频信号以及该风机的传动系统的振动信号和音频信号。
21.如权利要求20所述的故障概率预测模型的建立方法，其特征在于，步骤Sit中采集的该振动信号和该音频信号为模拟信号，步骤Sit之后、步骤S2t之前还包括以下步骤将该模拟信号转换为数字信号。
22.如权利要求15所述的故障概率预测模型的建立方法，其特征在于，步骤S5t中包括以下步骤采用该音频特征参数和HMM算法计算以下参数以为每个故障类型建立故障概率预测模型，该参数包括初始状态概率、状态转移概率矩阵、混合高斯函数权值、均值和协方差矩阵。
23.如权利要求15所述的故障概率预测模型的建立方法，其特征在于，步骤S5t中采用该音频特征参数、该振动特征参数和HMM算法为每个故障类型建立故障概率预测模型。
24.一种故障概率预测模型的建立系统，其特征在于，其包括 一振动传感装置，用于采集风机的振动信号； 一音频传感装置，用于在采集该振动信号的同时采集风机的音频信号； 一第一数据处理装置，用于预处理该振动信号并对预处理后的振动信号进行特征提取以获得振动特征参数； 一第二数据处理装置，用于预处理该音频信号并对预处理后的音频信号进行特征提取以获得音频特征参数； 一故障类型分析装置，用于将该振动特征参数代入一故障类型数据库中的故障类型振动模板以得到与该振动特征参数相对应的故障类型，该故障类型振动模板表示故障类型与振动特征参数之间的对应关系； 一故障类型预测模型建立装置，用于将该音频特征参数与该故障类型相关联以建立故障类型音频模板，该故障类型音频模板表示故障类型与音频特征参数之间的对应关系，多个故障类型音频模板构成该故障类型预测模型； 一故障概率预测模型建立装置，用于采用该音频特征参数和HMM算法为每个故障类型建立故障概率预测模型。
25.如权利要求24所述的故障概率预测模型的建立系统，其特征在于，该故障类型振动模板是通过采用该故障类型数据库中的历史振动数据和BP算法获得的。
26.如权利要求24所述的故障概率预测模型的建立系统，其特征在于，该第一数据处理装置用于通过时频分析方法来预处理该振动信号；和/或，该第二数据处理装置还包括 一分帧处理装置，用于对该音频信号分帧处理； 一预加重装置，用于预加重经分帧处理的该音频信号； 一加窗处理装置，用于加窗处理经预加重的该音频信号。
27.如权利要求26所述的故障概率预测模型的建立系统，其特征在于，时频分析方法选自时域波形分析方法、概率密度分析方法、自相关分析方法、趋势分析方法、幅值谱分析方法、功率谱分析方法、倒频谱分析方法、包络谱分析方法、细化谱分析方法、短时傅里叶分析方法和小波分析方法；和/或， 该分帧处理装置用于将该音频信号分为每10ms-30ms —帧，且相邻两帧有的重 3 3置；和/或， 该预加重装置为一预加重数字滤波器；和/或， 该加窗处理装置包括一汉明窗。
28.如权利要求24所述的故障概率预测模型的建立系统，其特征在于，第二数据处理装置用于通过以下方法中的一种或多种来对预处理后的音频信号进行特征提取该线性预测倒谱系数方法、美尔频率倒谱系数方法和小波包分析方法。
29.如权利要求24-28中任意一项所述的故障概率预测模型的建立系统，其特征在于，该振动传感装置用于采集该风机的机舱中的振动信号和该风机的传动系统的振动信号，和/或，该音频传感装置用于采集该风机的机舱中的音频信号和该风机的传动系统和音频信号。
30.如权利要求29所述的故障概率预测模型的建立系统，其特征在于，该振动传感装置采集的该振动信号和该音频传感装置采集的该音频信号为模拟信号，该振动传感装置还包括一第一 A/D转换器,该音频传感装置还包括一第二 A/D转换器,该第一 A/D转换器和该第二 A/D转换器均用于将该模拟信号转换为数字信号。
31.如权利要求24所述的故障概率预测模型的建立系统，其特征在于，该故障概率预测模型建立装置用于采用该音频特征参数和HMM算法计算以下参数以为每个故障类型建立故障概率预测模型，该参数包括初始状态概率、状态转移概率矩阵、混合高斯函数权值、均值和协方差矩阵。
32.如权利要求24所述的故障概率预测模型的建立系统，其特征在于，故障概率预测模型建立装置还用于采用该音频特征参数、该振动特征参数和HMM算法为每个故障类型建立故障概率预测模型。
33.一种风机监测预警方法，其特征在于，其包括以下步骤 步骤Sik、采集风机的振动信号，以及在采集该振动信号的同时采集风机的音频信号；步骤S2k、预处理该振动信号并对预处理后的振动信号进行特征提取以获得振动特征参数，以及预处理该音频信号并对预处理后的音频信号进行特征提取以获得音频特征参数；步骤S3k、将该振动特征参数代入一故障类型数据库中的故障类型振动模板以得到与该振动特征参数相对应的故障类型，和/或，将该音频特征参数代入一故障类型数据库中的故障类型音频模板以得到与该音频特征参数相对应的故障类型，其中，该故障类型振动模板表示故障类型与振动特征参数之间的对应关系，该故障类型音频模板表示故障类型与音频特征参数之间的对应关系； 步骤S4k、将该音频特征参数和/或该振动特征参数代入一故障概率预测模型中以计算该故障类型的概率，其中该故障概率预测模型采用如权利要求15-23中任意一项所述的故障概率预测模型的建立方法得到。
34.如权利要求33所述的风机监测预警方法，其特征在于，步骤S2k中通过时频分析方法来预处理该振动信号；和/或， 步骤S2k中预处理该音频信号还包括以下步骤 步骤S2ki、对该音频信号分帧处理； 步骤S2k2、预加重经分帧处理的该音频信号； 步骤S2k3、加窗处理经预加重的该音频信号。
35.如权利要求34所述的风机监测预警方法，其特征在于，时频分析方法选自时域波形分析方法、概率密度分析方法、自相关分析方法、趋势分析方法、幅值谱分析方法、功率谱分析方法、倒频谱分析方法、包络谱分析方法、细化谱分析方法、短时傅里叶分析方法和小波分析方法；和/或， 步骤S2ki中将该音频信号分为每10ms-30ms —帧，且相邻两帧有^■- 的重叠；和/或， 3 3 步骤S2k2中通过一预加重数字滤波器来实现预加重；和/或， 步骤S2k3中通过汉明窗来实现加窗处理。
36.如权利要求33所述的风机监测预警方法，其特征在于，步骤S2k中对预处理后的音频信号进行特征提取的方法选自线性预测倒谱系数方法、美尔频率倒谱系数方法和小波包分析方法。
37.如权利要求33-36中任意一项所述的风机监测预警方法，其特征在于，步骤Sik中采集的为该风机的机舱中的振动信号和音频信号以及该风机的传动系统的振动信号和音频信号。
38.如权利要求37所述的风机监测预警方法，其特征在于，步骤Sik中采集的该振动信号和该音频信号为模拟信号，步骤Sik之后、步骤S2k之前还包括以下步骤将该模拟信号转换为数字信号。
39.一种风机监测预警系统，其特征在于，其包括 一振动传感装置，用于采集风机的振动信号； 一音频传感装置，用于在采集该振动信号的同时采集风机的音频信号； 一第一数据处理装置，用于预处理该振动信号并对预处理后的振动信号进行特征提取以获得振动特征参数； 一第二数据处理装置，用于预处理该音频信号并对预处理后的音频信号进行特征提取以获得音频特征参数； 一故障类型分析装置，用于将该振动特征参数代入一故障类型数据库中的故障类型模板以得到与该振动特征参数相对应的故障类型，和/或，用于将该音频特征参数代入一故障类型数据库中的故障类型音频模板以得到与该音频特征参数相对应的故障类型，其中，该故障类型振动模板表示故障类型与振动特征参数之间的对应关系，该故障类型音频模板表示故障类型与音频特征参数之间的对应关系； 一故障概率分析装置，用于将该音频特征参数和/或该振动特征参数代入一故障概率数据库中的故障概率模型以计算该故障类型的概率，其中该故障概率预测模型采用如权利要求15-23中任意一项所述的故障概率预测模型的建立方法得到。
40.如权利要求39所述的风机监测预警系统，其特征在于，该第一数据处理装置用于通过时频分析方法来预处理该振动信号；和/或，该第二数据处理装置还包括 一分帧处理装置，用于对该音频信号分帧处理； 一预加重装置，用于预加重经分帧处理的该音频信号； 一加窗处理装置，用于加窗处理经预加重的该音频信号。
41.如权利要求40所述的风机监测预警系统，其特征在于，时频分析方法选自时域波形分析方法、概率密度分析方法、自相关分析方法、趋势分析方法、幅值谱分析方法、功率谱分析方法、倒频谱分析方法、包络谱分析方法、细化谱分析方法、短时傅里叶分析方法和小波分析方法；和/或， 该分帧处理装置用于将该音频信号分为每10ms-30ms —帧，且相邻两帧有的重 3 3置；和/或， 该预加重装置为一预加重数字滤波器；和/或， 该加窗处理装置包括一汉明窗。
42.如权利要求39所述的风机监测预警系统，其特征在于，第二数据处理装置用于通过以下方法中的一种或多种来对预处理后的音频信号进行特征提取该线性预测倒谱系数方法、美尔频率倒谱系数方法和小波包分析方法。
43.如权利要求39-42中任意一项所述的风机监测预警系统，其特征在于，该振动传感装置用于采集该风机的机舱中的振动信号和该风机的传动系统的振动信号，和/或，该音频传感装置用于采集该风机的机舱中的音频信号和该风机的传动系统和音频信号。
44.如权利要求43所述的风机监测预警系统，其特征在于，该振动传感装置采集的该振动信号和该音频传感装置采集的该音频信号为模拟信号，该振动传感装置还包括一第一A/D转换器，该音频传感装置还包括一第二 A/D转换器，该第一 A/D转换器和该第二 A/D转换器均用于将该模拟信号转换为数字信号。
全文摘要
本发明公开了一种故障类型预测模型的建立方法，包括采集风机的振动信号，以及在采集该振动信号的同时采集风机的音频信号；预处理该振动信号并对预处理后的振动信号进行特征提取以获得振动特征参数，以及预处理该音频信号并对预处理后的音频信号进行特征提取以获得音频特征参数；将该振动特征参数代入故障类型振动模板以得到与该振动特征参数相对应的故障类型；将该音频特征参数与该故障类型相关联以建立故障类型音频模板。本发明还公开了一种故障类型预测模型的建立系统、故障概率预测模型的建立方法及系统和风机监测预警系统及方法。本发明的风机监测预警系统及方法灵敏度更高、故障监测更准确。
文档编号G01M15/00GK102778358SQ20121018158
公开日2012年11月14日 申请日期2012年6月4日 优先权日2012年6月4日
发明者范兆骥 申请人:上海东锐风电技术有限公司