CN114764538B - 一种设备声音信号模式识别方法 - Google Patents
一种设备声音信号模式识别方法 Download PDFInfo
- Publication number
- CN114764538B CN114764538B CN202011603599.XA CN202011603599A CN114764538B CN 114764538 B CN114764538 B CN 114764538B CN 202011603599 A CN202011603599 A CN 202011603599A CN 114764538 B CN114764538 B CN 114764538B
- Authority
- CN
- China
- Prior art keywords
- layer
- sound signal
- sound
- equipment
- feature
- 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
- 230000005236 sound signal Effects 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000013598 vector Substances 0.000 claims abstract description 40
- 238000012545 processing Methods 0.000 claims abstract description 30
- 238000000605 extraction Methods 0.000 claims abstract description 25
- 238000009432 framing Methods 0.000 claims abstract description 18
- 230000009467 reduction Effects 0.000 claims abstract description 17
- 238000012567 pattern recognition method Methods 0.000 claims abstract description 5
- 239000002245 particle Substances 0.000 claims description 29
- 238000011176 pooling Methods 0.000 claims description 25
- 238000003745 diagnosis Methods 0.000 claims description 21
- 238000001514 detection method Methods 0.000 claims description 19
- 238000013527 convolutional neural network Methods 0.000 claims description 18
- 238000004422 calculation algorithm Methods 0.000 claims description 15
- 238000005457 optimization Methods 0.000 claims description 13
- 238000012549 training Methods 0.000 claims description 13
- 210000002569 neuron Anatomy 0.000 claims description 12
- 230000004913 activation Effects 0.000 claims description 9
- 238000001994 activation Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000004891 communication Methods 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000013528 artificial neural network Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000003909 pattern recognition Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 29
- 230000006870 function Effects 0.000 abstract description 21
- 230000002159 abnormal effect Effects 0.000 abstract description 16
- 238000005516 engineering process Methods 0.000 abstract description 10
- 238000007689 inspection Methods 0.000 description 14
- 238000012423 maintenance Methods 0.000 description 12
- 230000005856 abnormality Effects 0.000 description 10
- 238000013135 deep learning Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 230000018109 developmental process Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 208000009119 Giant Axonal Neuropathy Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 201000003382 giant axonal neuropathy 1 Diseases 0.000 description 1
- 238000001093 holography Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 238000013526 transfer learning Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Software Systems (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- Biophysics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Medical Informatics (AREA)
- Testing And Monitoring For Control Systems (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种设备声音信号模式识别方法。所述设备声音信号模式识别方法包括分帧处理模块，所述分帧处理模块用于对所述原始声音信号进行加窗分帧处理；特征提取模块，所述特征提取模块用于对分帧处理后的声音信号进行特征提取，构造声音信号特征向量；PCA降维处理模块，所述PCA降维处理模块用于对所述声音信号特征向量进行降维处理，提取最能反应所述设备本体运行状态的声音信号的特征分量。本发明提供的设备声音信号模式识别方法基于物联网技术，在设备本体安装声音传感器，通过声音传感器采集设备的声音数据，从中提取声音特征值，利用设备声音信号异常识别模型，实现设备运行状态监测和设备故障报警等功能。
Description
技术领域
本发明涉及信号识别技术领域，尤其涉及一种设备声音信号模式识别方法。
背景技术
随着电厂设备的大型化、集中化、高速化、连续化和自动化，突发的故障停机造成的损失越来越大。由于大多数设备、备件的故障间隔离散性较大，定期维修不可避免地会造成较高的维修成本和较多的设备机器失效，因此做到预知维修非常重要，主动地发现问题和解决问题，把事故消灭在萌芽状态，能有效降低维修成本，降低事故停机率，因而具有很高的投资收益比。为了符合这一要求，设备状态监测和诊断已成为行业工具。状态监测是一种感知设备健康的技术，操作信息并分析此信息以量化设备状况。这样做是为了使潜在的问题能够在其发展的早期被检测和诊断，并且在它们变得严重之前通过适当的恢复措施进行校正。
随着物联网技术和深度学习技术的快速发展，基于模型驱动的技术可以提高到了一个新的高度：不再依赖于信号技术就可以提取更为精确的信号特征，特别是最近几年深度学习技术在图像识别领域的成功应用，其识别的能力几乎可以到达人的识别水平。在过去，电厂生产环境中有很大的噪音，不同设备的声音混杂在相对集中的空间中，由于针对声音信号没有很好的技术手段去处理，所以仅仅只能依靠具有丰富专业经验的巡检人员才能听出设备的异常声音，而现在有了深度学习技术就可以处理声音信号的识别问题了。在实际系统中，设备在不同工作条件下的运行状态不同，当工作条件发生改变时(如故障发生)，设备运行状态将发生改变。虽然状态改变过程中设备声学信号产生的机理比较模糊，但这种声学信号往往具有非平稳特性，故可以采用统计模型理论来分析处理。设备状态在改变过程中往往会引起声学信号结构的变化、出现不同的声学信号特征，可以通过设备状态声音的音色、音量的大小、频率的高低等声学信号特征的变化判断出设备的运行状态，甚至用以判别设备故障类型和发生部位。
国外对机械噪声研究起步比较早，英国人R.A.Collacott在二十世纪70年代就率先提出了利用噪声信号对机械设备进行运行状态监测以及故障诊断。但由于声音信号比振动信号容易受到环境的污染，故基于声音信号的设备状态监测研究非常少于基于振动和电流的信号分析。使用声音信号进行故障检测一直在吸引研究人员的兴趣。Y.Ono等人总结了用于故障检测的声学声音信号的分析，提出了一种利用F值对电机异常声音进行状态监测和轴承故障检测的方法。Akcay等在嘈杂环境下利用5个麦克风通过交叉谱算法对电机的噪声进行了精确识别。Y.Ono等提出提出了一种异常检测方法，用于在不使用异常信号的情况下观察运行中的电机的声音信号，有效地检测出一小部分特征中出现的异常。该方法最佳地估计各种特征对观察信号和正常信号分布之间的相异性得分的贡献率。所提出的方法强调了一小部分窄频率范围的声音，并且它实现了高达76％的误差减少率。
近年来随着深度学习的高速发展，国内外越来越多的研究人员利用深度学习方法对设备状态监测。特别是在振动监测领域，利用深度学习强大的特征自学习能力，将振动信号转为二维时频图像进行状态检测和故障诊断。振动信号作为一维时间序列信号，各时刻的数据点具有关联性，若直接将其转换为二维形式，则原始序列中的空间相关性会遭到破坏，有可能造成故障信息的丢失。目前绝大部分基于深度学习的故障诊断方法并没有直接从原始信号中获取数据信息。在基于声音信号对设备进行状态监测和故障诊断的研究更是少之又少。
发电行业的巡检人员利用听针来实现对设备异常的检查，但基于声音的设备故障预警与诊断系统鲜有建设。
WenboLu等提出了一种基于近场声全息(NAH)和灰度共生矩阵(GLCM)的基于声学的诊断轴承故障诊断新方案。在锅炉四管泄露监测方面也有应用。但这些应用基本上都是基于传统的机器学习对信号进行分析，适应性不强。
深度学习在设备健康监测方面有大量的研究，AE、CNN、LSTM、GANs和迁移学习等各种技术都有所涉及，但全部是基于振动监测信号的研究。
目前设备故障诊断的识别技术手段主要包括：基于人工巡视的现场定期检查、基于噪声传感器的声音监听系统、基于DCS或SIS的盘面数据人工监控三种。
其中基于人工的现场巡视是电厂制定了一系列的巡检制度，巡检人员通过手持设备可以实现对设备的振动、温度、声音等信号的收集，依靠人的专业素养来判断设备的运行状况，发现了不少设备运行的故障，但是这种制度过于依赖人的能力，再加上巡检人员与设备数据的不匹配造成巡检人员的工作量大大的增加，容易造成遗漏。
基于DCS的盘面数据人工监控主要是该类系统一般集成在DCS系统中，它关注的是系统的安全运行，一般通过设置多级阈值实现不同等级的告警，但没有对设备信号进行深入的分析，无法实现故障分类。
基于噪声传感器的声音监听系统是在现场安装噪声传感器来监听设备的声音，判断设备的运行状态，此类方法基于传感器监听现场的环境声音，通过声音的能量值来判断设别的运行状态，但真实环境是现场多种设备运转是发出的声音极大，通过能量值判断识别的故障非常少，由于现场一旦有其他高分贝的声音发出，则存在误报的可能。
因此，有必要提供一种新的设备声音信号模式识别方法解决上述技术问题。
发明内容
本发明为设备状态监测类系统的内部算法，用于实现基于物联网技术，在设备本体安装声音传感器，通过声音传感器采集设备的声音数据，从中提取声音特征值，利用设备声音信号异常识别模型，实现设备运行状态监测和设备故障报警等功能。
为解决上述技术问题，本发明提供的设备声音信号模式识别方法包括：
安装在设备本体上的声音传感器，所述声音传感器用于采集所述设备本体的原始声音信号；
无线通讯模块，所述无线通讯模块与所述声音传感器相连接；
服务器端，所述服务器端通过所述无线通讯模块与所述声音传感器相连接；
分帧处理模块，所述分帧处理模块用于对所述原始声音信号进行加窗分帧处理；
特征提取模块，所述特征提取模块用于对分帧处理后的声音信号进行特征提取，构造声音信号特征向量；
PCA降维处理模块，所述PCA降维处理模块用于对所述声音信号特征向量进行降维处理，提取最能反应所述设备本体运行状态的声音信号的特征分量；
OCSVM模型训练模块，所述OCSVM模型训练模块用于寻找一个最优超平面，确定正常样本的边界；
OCSVM异常检测模型，所述OCSVM异常检测模型用于对声音信号进行特征提取和PCA降维处理后的数据进行检测，实现对电力设备的模式识别；
粒子群优化算法模块，所述粒子群优化算法模块用于对所述OCSVM异常检测模型进行参数寻优与结构优化；
1D-CNN故障诊断模型，所述1D-CNN故障诊断模型用于模式识别后，对具体故障类型进行分类识别。
优选的，所述分帧处理模块对所述原始声音信号进行加窗分帧处理时，窗函数选用矩形窗。
优选的，所述特征提取模块包括时域特征提取模块和频域特征提取模块；所述时域特征提取模块用于提取声音信号在时域上的17种特征指标，分别为均值、绝对均值、平均均值、有效值、方根幅值、峰值、峰峰值、方差、标准差、偏斜度、峭度、波形指标、波峰指标、脉冲指标、裕度指标、偏斜指标、峭度指标，对每个声音信号提取时域统计特征，构造时域特征向量，记作α＝[α1,α2,α3,...,α16,α17]；所述频域特征提取模块用于提取原始声音信号的频域特征上的4种特征指标，分别为重心频率、均方频率、均方根频率和频率方差，对每个声音信号提取频域统计特征，构造频域特征向量，记作β＝[β1,β2,β3,β4]；将时域特征与频域特征结合组成21个声音信号特征指标，构造声音信号特征向量δ＝[δ1,δ2,δ3,...,δ20,δ21]。
优选的，所述OCSVM模型训练模块求解二次规划问题的函数为：
s.t.Φ(xi)ω...ρ-ξi,ξi...0
其中，xi为原空间中的样本，l为训练样本的数量，Φ为原始空间到特征空间的映射，ω和ρ分别为特征空间中所需超平面的法向量和补偿，可调参数ν∈(0,1)为控制总样本数中错误样本比例的上界，松弛变量ξi为某些训练样本被错误分类的程度；ω和ρ是这个二次规划问题的解，则其决策函数为：f(x)＝sgn(Φ(x)ω-ρ)。
优选的，所述粒子群优化算法模块的算法流程如下：
S1：设置最大迭代次数目标函数的自变量个数，粒子的最大速度，位置信息为整个搜索空间，在速度区间和搜索空间上随机初始化速度和位置，设置粒子群规模为M，每个粒子随机初始化一个飞翔速度；
S2：定义适应度函数，个体极值为每个粒子找到的最优解，从这些最优解找到一个全局值，叫做本次全局最优解；与历史全局最优比较，进行更新；
S3：更新速度和位置的公式：
Vid＝ωVid+C1random(0,1)(Pid-Xid)+C2random(0,1)(Pgd-Xid)
Xid＝Xid+Vid
其中ω称为惯性因子，调整ω大小，对全局寻优性能和局部寻优性能进行调整；C1和C2为加速常数；Pid表示第i个变量的个体极值的第d维，Pgd表示全局最优解的第d维；
S4：终止条件：达到设定迭代次数或者迭代数之间的差值满足最小界限。
优选的，所述1D-CNN故障诊断模型包括三部分：输入层、特征提取层、分类层；所述输入层为对原始声音数据分段处理后的直接输入；所述特征提取层包括3个卷积层和3个池化层，接收来自输入层的数据，对原始声音信号的特征提取，池化层选择最大池化算子实现对特征矢量的降维，同时提高非线性特征的鲁棒性；所述分类层由两个全连接层组成，第2个全连接层神经元个数与标签故障数目一致，利用Softmax回归分类器实现分类输出。
优选的，所述卷积层用于卷积核对前一层输出的特征矢量进行卷积操作，利用非线性激活函数构建输出特征矢量，每一层的输出均为对多输入特征的卷积结果，其数学模型可以描述为:
式中，Mj为输入特征矢量，l为第l层网络，k为卷积核，b为网络偏置，为第l层输出/>为第l层输入。
优选的，所述池化层用于池化，池化是非线性下采样的一种形式，通过减少网络的参数来减小计算量，并且能够在一定程度上控制过拟合，在所述卷积层的后面会加上一个所述池化层；最大池化是用不重叠的矩形框将输入层分成不同的区域，对于每个矩形框的数取最大值作为输出层，最大池化的变换函数如下所示：
式中：表示第l层第i个特征矢量中第t个神经元的值，t∈[(j-1)W+1,jW]；W为池化区域的宽度；/>表示第l+1层神经元对应的值。
优选的，所述全连接层是一种传统的前馈神经网络，之后在输出端使用Softmax函数作为激活函数，所述全连接层的所有神经元连接到前一层的所有激活；所述全连接层起到将学到的分布式特征表示映射到样本标记空间的作用，模型可表述为：
O＝f(bo+wofv)
式中：fv为特征矢量；wo、bo分别为偏差向量和权值矩阵。
与相关技术相比较，本发明提供的设备声音信号模式识别方法具有如下有益效果：
本发明提供一种设备声音信号模式识别方法，基于智能声学的设备状态监测系统程序的内部核心算法，基于该算法构建的设备状态监测系统能够复现运行管理人员及巡检人员的设备状态监控能力，实现基于声音识别的设备远程监控，为设备巡检及事故发展追溯提供分析支持与工具辅助，实现设备故障的早期识别，对于减小了设备故障的损伤规模，降低设备维修成本，妥善制定设备维修计划，缩短维修时间具有极为重要的现实意义；
经济效益分析上，实现设备的早期故障识别能够缩短检修时间，减少设备更换次数，避免因单一设备造成大范围的设备停机，挽回了巨大的经济损失；
本模型具备前期投入小，普适性强等特点，同时基于智能声学的设备状态监测系统作为电力生产运行部门的普遍需求，具备广阔的应用推广应用前景。
附图说明
图1为本发明提供的设备声音信号模式识别方法的一种较佳实施例的识别流程图；
图2为1D-CNN故障诊断模型结构示意图。
具体实施方式
下面结合附图和实施方式对本发明作进一步说明。
请结合参阅图1、图2，其中，图1为本发明提供的设备声音信号模式识别方法的一种较佳实施例的识别流程图；图2为1D-CNN故障诊断模型结构示意图。基于物联网技术，在设备本体安装声音传感器，通过声音传感器采集设备的声音数据，采用wifi的传输方式把声音数据传输到服务器端，在服务器端部署设备声音信号模式识别方法，建立一种“先异常检测后故障诊断”的设备模式识别方法。
1.提取声音传感器对电力设备采集到的原始声音信号,对原始声音信号进行加窗分帧处理，窗函数选用矩形窗。对分帧处理后的声音信号进行特征提取,首先提取声音信号在时域上的17种特征指标。指标名称及计算公式如表1-1所示，其中x(i)其中为分帧处理后的声音信号，K是声音信号的采样点数。
表1-1
对每个声音信号提取时域统计特征，构造时域特征向量，记作α＝[α1,α2,α3,...,α16,α17]。
提取原始声音信号的频域特征，选取4种特征指标，指标名称与计算公式如表1-2所示
表1-2
对每个声音信号提取频域统计特征，构造频域特征向量，记作β＝[β1,β2,β3,β4]；
将时域特征与频域特征结合组成21个声音信号特征指标，构造声音信号特征向量δ＝[δ1,δ2,δ3,...,δ20,δ21]。
2.对声音信号特征向量进行降维处理。PCA是一种数据分析方式，可以实现对高维数据的降维，可以提取数据中的主要特征分量。使用PCA方法对包含21个特征指标的样本进行降维，提取最能反应设备运行状态的声音信号的特征分量。
3.通过对声音信号进行特征提取和PCA降维处理后的数据，作为异常检测模型OCSVM的模型输入。异常样本标签为-1，正常样本标签为1，目标是寻找一个最优超平面，确定正常样本的边界。
OCSVM，其基本思路是寻找一个决策超平面，使大部分目标样本位于该超平面的一侧，而大部分非目标样本位于超平面的另一侧。OCSVM求解如下二次规划问题：
s.t.Φ(xi)ω...ρ-ξi,ξi…0
其中，xi为原空间中的样本，l为训练样本的数量，Φ为原始空间到特征空间的映射，ω和ρ分别为特征空间中所需超平面的法向量和补偿，可调参数ν∈(0,1)为控制总样本数中错误样本比例的上界，松弛变量ξi为某些训练样本被错误分类的程度。ω和ρ是这个二次规划问题的解，则其决策函数为：
f(x)＝sgn(Φ(x)ω-ρ)
4.利用粒子群优化算法(Particle Swarm Optimization，PSO)对OCSVM模型进行参数寻优与结构优化。粒子群算法通过设计一种无质量的粒子来模拟鸟群中的鸟，粒子仅具有两个属性：速度和位置，速度代表移动的快慢，位置代表移动的方向。每个粒子在搜索空间中单独的搜寻最优解，并将其记为当前个体极值，并将个体极值与整个粒子群里的其他粒子共享，找到最优的个体极值作为整个粒子群的当前全局最优解，粒子群中的所有粒子根据自己找到的当前个体极值和整个粒子群共享的当前全局最优解来调整自己的速度和位置。PSO算法流程如下：
步骤1.设置最大迭代次数目标函数的自变量个数，粒子的最大速度，位置信息为整个搜索空间，我们在速度区间和搜索空间上随机初始化速度和位置，设置粒子群规模为M，每个粒子随机初始化一个飞翔速度。
步骤2.定义适应度函数，个体极值为每个粒子找到的最优解，从这些最优解找到一个全局值，叫做本次全局最优解。与历史全局最优比较，进行更新。
步骤3.更新速度和位置的公式：
Vid＝ωVid+C1random(0,1)(Pid-Xid)+C2random(0,1)(Pgd-Xid)
Xid＝Xid+Vid
ω称为惯性因子，调整ω大小，对全局寻优性能和局部寻优性能进行调整。C1和C2为加速常数。Pid表示第i个变量的个体极值的第d维，Pgd表示全局最优解的第d维。
步骤4.终止条件：达到设定迭代次数或者迭代数之间的差值满足最小界限。
5.通过建立的OCSVM异常检测模型，对电力设备声音信号进行检测，实现对电力设备的模式识别。
6.模式识别后，对具体故障类型进行分类识别，建立一维卷积神经网络(1D-CNN)故障诊断模型。一维卷积神经网络包括：
卷积层
卷积核对前一层输出的特征矢量进行卷积操作，利用非线性激活函数构建输出特征矢量，每一层的输出均为对多输入特征的卷积结果，其数学模型可以描述为:
式中，Mj为输入特征矢量，l为第l层网络，k为卷积核，b为网络偏置，为第l层输出/>为第l层输入。
池化层
池化是非线性下采样的一种形式，通过减少网络的参数来减小计算量，并且能够在一定程度上控制过拟合。在卷积层的后面会加上一个池化层。最大池化是用不重叠的矩形框将输入层分成不同的区域，对于每个矩形框的数取最大值作为输出层，最大池化的变换函数如下所示：
式中：表示第l层第i个特征矢量中第t个神经元的值，t∈[(j-1)W+1,jW]；W为池化区域的宽度；/>表示第l+1层神经元对应的值。
全连接层
全连接层是一种传统的前馈神经网络，之后在输出端使用Softmax函数作为激活函数，这一层的所有神经元连接到前一层的所有激活。全连接层起到将学到的“分布式特征表示”映射到样本标记空间的作用，模型可表述为：
O＝f(bo+wofv)
式中：fv为特征矢量；wo、bo分别为偏差向量和权值矩阵。
7.现实工业生产中，正常样本远高于故障样本，1D-CNN故障诊断模型只通过故障数据进行训练。在OCSVM异常检测模型首先对原始声音信号进行异常识别后，将异常数据由人工进行类别标记，标记后输入到1D-CNN故障诊断模型进行训练。
8.通过1D-CNN故障诊断模型识别后，对具体故障类型进行报警。
基于声音信号的电力设备的模式模型主要包括：
1D-CNN故障诊断模型结构图如图2所示
本专利构建的一维卷积神经网络(1D-CNN)结构如图2所示，包括三部分：输入层、特征提取层、分类层。输入层为对原始声音数据分段处理后的直接输入。特征提取层包括3个卷积层和3个池化层，接收来自输入层的数据，对原始声音信号的特征提取。池化层选择最大池化算子实现对特征矢量的降维，同时提高非线性特征的鲁棒性。分类层由两个全连接层组成，第2个全连接层神经元个数与标签故障数目一致，利用Softmax回归分类器实现分类输出。
“先异常检测后故障识别”的设备模式识别具体实现过程如图1所示；
步骤01.采集电力设备的原始声音信号，信号包含正常状态数据和异常状态数据
步骤02.数据预处理，填写缺失值，识别或删除离群点解决不一致性来“清理”数据。对清理后的数据进行分帧加窗处理。
步骤03.特征提取，对原始声音信号进行时域和频域上的特征提取，共提取21个特征指标
步骤04.PCA数据降维，对21个特征指标进行降维后得到异常检测模型的训练集X0，X0为l×m的矩阵，其中l为声音信号的样本数量，m为经过降维后每个样本的维数，即每个样本所包含的特征数。
步骤05.降维后的声音数据作为单类支持向量机模型的输入
步骤06.首先训练基于正常数据的正常OCSVM异常检测模型，即识别正常状态数据和非正常状态数据。粒子群算法优化OCSVM参数，核函数选用径向基核函数。
步骤07.正常OCSVM异常检测模型识别出的正常样本数据输入到正常样本数据库，同时将数据反馈到正常OCSVM异常检测模型，更新模型，提升模型精度与泛化能力。被正常OCSVM异常检测模型识别出的异常数据，输入到1D-CNN故障诊断模型，该模型最初可以由历史故障数据训练或者由人为对异常数据进行准确的故障类型标定，标定后作为故障诊断模型的输入来训练模型。
步骤08.异常数据通过1D-CNN故障诊断模型识别后，在模型中加入的Softmax层会对模型判定某种类型故障给出判定概率。概率大于0.9，认定模型可以对故障做出准确判定；概率小于0.9，认定模型不能对故障做出准确判定。
步骤9.对不能准确判定故障类型的数据由人工进行标记，标记后作为输入更新模型，对故障诊断模型重新训练。
步骤10.所有给出准确判定的数据样本保存到正常和故障样本数据库。
与相关技术相比较，本发明提供的设备声音信号模式识别方法具有如下有益效果：
本发明提供一种设备声音信号模式识别方法，基于智能声学的设备状态监测系统程序的内部核心算法，基于该算法构建的设备状态监测系统能够复现运行管理人员及巡检人员的设备状态监控能力，实现基于声音识别的设备远程监控，为设备巡检及事故发展追溯提供分析支持与工具辅助，实现设备故障的早期识别，对于减小了设备故障的损伤规模，降低设备维修成本，妥善制定设备维修计划，缩短维修时间具有极为重要的现实意义；
经济效益分析上，实现设备的早期故障识别能够缩短检修时间，减少设备更换次数，避免因单一设备造成大范围的设备停机，挽回了巨大的经济损失；
本模型具备前期投入小，普适性强等特点，同时基于智能声学的设备状态监测系统作为电力生产运行部门的普遍需求，具备广阔的应用推广应用前景。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其它相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (5)
1.一种设备声音信号模式识别方法，其特征在于，包括：
安装在设备本体上的声音传感器，所述声音传感器用于采集所述设备本体的原始声音信号；
无线通讯模块，所述无线通讯模块与所述声音传感器相连接；
服务器端，所述服务器端通过所述无线通讯模块与所述声音传感器相连接；
分帧处理模块，所述分帧处理模块用于对所述原始声音信号进行加窗分帧处理；
特征提取模块，所述特征提取模块用于对分帧处理后的声音信号进行特征提取，构造声音信号特征向量；
PCA降维处理模块，所述PCA降维处理模块用于对所述声音信号特征向量进行降维处理，提取最能反应所述设备本体运行状态的声音信号的特征分量；
OCSVM模型训练模块，所述OCSVM模型训练模块用于寻找一个最优超平面，确定正常样本的边界；
OCSVM异常检测模型，所述OCSVM异常检测模型用于对声音信号进行特征提取和PCA降维处理后的数据进行检测，实现对电力设备的模式识别；
粒子群优化算法模块，所述粒子群优化算法模块用于对所述OCSVM异常检测模型进行参数寻优与结构优化；
1D-CNN故障诊断模型，所述1D-CNN故障诊断模型用于模式识别后，对具体故障类型进行分类识别；
所述1D-CNN故障诊断模型包括三部分：输入层、特征提取层、分类层；所述输入层为对原始声音数据分段处理后的直接输入；所述特征提取层包括3个卷积层和3个池化层，接收来自输入层输出的数据，对原始声音信号的特征提取，池化层选择最大池化算子实现对特征矢量的降维，同时提高非线性特征的鲁棒性；所述分类层由两个全连接层组成，第2个全连接层神经元个数与标签故障数目一致，利用Softmax回归分类器实现分类输出；
所述卷积层中的卷积核对前一层输出的特征矢量进行卷积操作，利用非线性激活函数构建输出特征矢量，每一层的输出均为对多输入特征的卷积结果，其数学模型可以描述为:
式中，Mj为输入特征矢量，l为第l层网络，k为卷积核，b为网络偏置，为第l层输出，/>为第l层输入；
所述池化层采用最大池化，池化是非线性下采样的一种形式，通过减少网络的参数来减小计算量，并且能够在一定程度上控制过拟合，在所述卷积层的后面会加上一个所述池化层；最大池化是用不重叠的矩形框将输入层分成不同的区域，对于每个矩形框的数取最大值作为输出层，最大池化的变换函数如下所示：
式中：表示第l层第i个特征矢量中第t个神经元的值，t∈[(j-1)W+1,jW]；W为池化区域的宽度；/>表示第l+1层神经元对应的值；
所述全连接层是一种传统的前馈神经网络，之后在输出端使用Softmax函数作为激活函数，所述全连接层的所有神经元连接到前一层的所有激活；所述全连接层起到将学到的分布式特征表示映射到样本标记空间的作用，模型可表述为：
O＝f(bo+wofv)
式中：fv为特征矢量；wo、bo分别为偏差向量和权值矩阵。
2.根据权利要求1所述的设备声音信号模式识别方法，其特征在于，所述分帧处理模块对所述原始声音信号进行加窗分帧处理时，窗函数选用矩形窗。
3.根据权利要求1所述的设备声音信号模式识别方法，其特征在于，所述特征提取模块包括时域特征提取模块和频域特征提取模块；所述时域特征提取模块用于提取声音信号在时域上的17种特征指标，分别为均值、绝对均值、平均均值、有效值、方根幅值、峰值、峰峰值、方差、标准差、偏斜度、峭度、波形指标、波峰指标、脉冲指标、裕度指标、偏斜指标、峭度指标，对每个声音信号提取时域统计特征，构造时域特征向量，记作α＝[α1,α2,α3,...,α16,α17]；所述频域特征提取模块用于提取原始声音信号的频域特征上的4种特征指标，分别为重心频率、均方频率、均方根频率和频率方差，对每个声音信号提取频域统计特征，构造频域特征向量，记作β＝[β1,β2,β3,β4]；将时域特征与频域特征结合组成21个声音信号特征指标，构造声音信号特征向量δ＝[δ1,δ2,δ3,...,δ20,δ21]。
4.根据权利要求1所述的设备声音信号模式识别方法，其特征在于，所述OCSVM模型训练模块求解二次规划问题的函数为：
subject to(ω·Φ(Xi))≥ρ-ξi,ξi≥0
其中，xi为原空间中的样本，l为训练样本的数量，Φ为原始空间到特征空间的映射，ω和ρ分别为特征空间中所需超平面的法向量和补偿，可调参数ν∈(0,1)为控制总样本数中错误样本比例的上界，松弛变量ξi为某些训练样本被错误分类的程度；ω和ρ是这个二次规划问题的解，则其决策函数为：f(x)＝sgn(Φ(x)ω-ρ)。
5.根据权利要求1所述的设备声音信号模式识别方法，其特征在于，所述粒子群优化算法模块的算法流程如下：
S1：设置最大迭代次数目标函数的自变量个数，粒子的最大速度，位置信息为整个搜索空间，在速度区间和搜索空间上随机初始化速度和位置，设置粒子群规模为M，每个粒子随机初始化一个飞翔速度；
S2：定义适应度函数，个体极值为每个粒子找到的最优解，从这些最优解找到一个全局值，叫做本次全局最优解；与历史全局最优比较，进行更新；
S3：更新速度和位置的公式：
Vid＝ωVid+C1random(0,1)(Pid-Xid)+C2random(0,1)(Pgd-Xid)
Xid＝Xid+Vid
其中ω称为惯性因子，调整ω大小，对全局寻优性能和局部寻优性能进行调整；C1和C2为加速常数；Pid表示第i个变量的个体极值的第d维，Pgd表示全局最优解的第d维；
S4：终止条件：达到设定迭代次数或者迭代数之间的差值满足最小界限。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011603599.XA CN114764538B (zh) | 2020-12-30 | 2020-12-30 | 一种设备声音信号模式识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011603599.XA CN114764538B (zh) | 2020-12-30 | 2020-12-30 | 一种设备声音信号模式识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114764538A CN114764538A (zh) | 2022-07-19 |
| CN114764538B true CN114764538B (zh) | 2024-04-26 |
Family
ID=82364120
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011603599.XA Active CN114764538B (zh) | 2020-12-30 | 2020-12-30 | 一种设备声音信号模式识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114764538B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116543789B (zh) * | 2023-07-06 | 2023-09-29 | 中国电信股份有限公司 | 设备异常识别方法、装置、设备及介质 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017192821A1 (en) * | 2016-05-06 | 2017-11-09 | Massachusetts Institute Of Technology | Method and apparatus for efficient use of cnc machine shaping tool including cessation of use no later than the onset of tool deterioration by monitoring audible sound during shaping |
| WO2019090878A1 (zh) * | 2017-11-09 | 2019-05-16 | 合肥工业大学 | 基于向量值正则核函数逼近的模拟电路故障诊断方法 |
| CN110785761A (zh) * | 2017-05-17 | 2020-02-11 | 美商新思科技有限公司 | 存在时钟动态重编程时仿真时间线的压缩方法 |
| CN110867196A (zh) * | 2019-12-03 | 2020-03-06 | 桂林理工大学 | 一种基于深度学习及声音识别的机器设备状态监测系统 |
| CN112037012A (zh) * | 2020-08-14 | 2020-12-04 | 百维金科(上海)信息科技有限公司 | 一种基于pso-bp神经网络的互联网金融信用评价方法 |
| WO2020255173A2 (en) * | 2019-06-20 | 2020-12-24 | Tata Consultancy Services Limited | Method and system for optimization of combined cycle gas turbine operation |
-
2020
- 2020-12-30 CN CN202011603599.XA patent/CN114764538B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017192821A1 (en) * | 2016-05-06 | 2017-11-09 | Massachusetts Institute Of Technology | Method and apparatus for efficient use of cnc machine shaping tool including cessation of use no later than the onset of tool deterioration by monitoring audible sound during shaping |
| CN110785761A (zh) * | 2017-05-17 | 2020-02-11 | 美商新思科技有限公司 | 存在时钟动态重编程时仿真时间线的压缩方法 |
| WO2019090878A1 (zh) * | 2017-11-09 | 2019-05-16 | 合肥工业大学 | 基于向量值正则核函数逼近的模拟电路故障诊断方法 |
| WO2020255173A2 (en) * | 2019-06-20 | 2020-12-24 | Tata Consultancy Services Limited | Method and system for optimization of combined cycle gas turbine operation |
| CN110867196A (zh) * | 2019-12-03 | 2020-03-06 | 桂林理工大学 | 一种基于深度学习及声音识别的机器设备状态监测系统 |
| CN112037012A (zh) * | 2020-08-14 | 2020-12-04 | 百维金科(上海)信息科技有限公司 | 一种基于pso-bp神经网络的互联网金融信用评价方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于粒子群优化SVM的面向对象软件缺陷预测模型;朱朝阳;陈相舟;王志宏;张信明;;计算机应用;20171220(第S2期);全文 * |
| 抗环境音干扰的设备声音故障监测方法;李兰村;廉东本;毛立爽;;计算机系统应用;20190615(第06期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114764538A (zh) | 2022-07-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109655259B (zh) | 基于深度解耦卷积神经网络的复合故障诊断方法及装置 | |
| Zan et al. | Control chart pattern recognition using the convolutional neural network | |
| Xia et al. | Fault diagnosis for rotating machinery using multiple sensors and convolutional neural networks | |
| CN110334764B (zh) | 基于集成深度自编码器的旋转机械智能故障诊断方法 | |
| CN106682814B (zh) | 一种基于故障知识库的风电机组故障智能诊断方法 | |
| CN112906473B (zh) | 一种旋转设备故障诊断方法 | |
| CN113255848B (zh) | 基于大数据学习的水轮机空化声信号辨识方法 | |
| CN111505424A (zh) | 一种基于深度卷积神经网络的大型实验装置电力设备故障诊断方法 | |
| CN107273924B (zh) | 基于模糊聚类分析的多数据融合的电厂故障诊断方法 | |
| CN110533167B (zh) | 一种电动阀门执行机构用故障诊断方法及诊断系统 | |
| CN113837000A (zh) | 一种基于任务排序元学习的小样本故障诊断方法 | |
| CN110647911A (zh) | 一种基于主成分分析和深度信念网络的轴承故障诊断方法 | |
| CN116380445B (zh) | 基于振动波形的设备状态诊断方法及相关装置 | |
| CN113657221A (zh) | 一种基于智能感知技术的电厂设备状态监测方法 | |
| CN114619292A (zh) | 一种基于小波降噪和注意力机制融合gru网络的铣削加工刀具磨损监测方法 | |
| CN111753891A (zh) | 一种无监督特征学习的滚动轴承故障诊断方法 | |
| CN115859077A (zh) | 一种变工况下多特征融合的电机小样本故障诊断方法 | |
| CN114118225A (zh) | 发电机剩余寿命预测方法、系统、电子设备和存储介质 | |
| CN114764538B (zh) | 一种设备声音信号模式识别方法 | |
| CN112763215B (zh) | 一种基于模块化联邦深度学习的多工况在线故障诊断方法 | |
| CN111783941B (zh) | 一种基于概率置信度卷积神经网络的机械设备诊断分类方法 | |
| CN113984389A (zh) | 一种基于多感受野与改进胶囊图神经网络的滚动轴承故障诊断方法 | |
| CN117056678A (zh) | 一种基于小样本的机泵设备运行故障诊断方法及装置 | |
| CN116821610A (zh) | 一种利用大数据优化风力发电效率的方法 | |
| CN117009916A (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 |