一种基于声纹成像的GIS故障诊断方法及系统与流程

一种基于声纹成像的gis故障诊断方法及系统
技术领域
1.本发明涉及电力设备噪声振动信号采集分析技术领域，具体地说是一种基于声纹成像的gis故障诊断方法及系统。


背景技术：

2.目前gis设备机械故障检测的主要方式为振动信号分析。通过在gis设备壳体表面放置振动传感器，分析传感器接收到的振动信号，根据振动信号频谱的变化即可判断gis是否出现机械故障，并对故障类型进行分类。但采用振动分析时，振动传感器只能安装在设备表面，只能获取局部位置的振动信号，难以反映gis设备整体的运行状况。
3.gis设备由于机械故障所造成的异常振动会辐射出声信号，根据声信号与振动信号的同理性，可用gis辐射出的声信号进行故障诊断。基于声学的故障诊断具有非接触式、操作简单快捷、不影响设备正常运行的优点。
4.但是，采用声学信号进行故障诊断最大的难点在于背景噪声的干扰，不同于振动信号，声信号的信噪比较振动信号低的多，特别是在环境干扰噪声较大的情况下，故障源辐射出的声信号极易被其他干扰源所覆盖，且当传声器阵列的测点位置选择不当会导致该处的声学特征采集效果变差，严重影响后续的故障诊断。


技术实现要素：

5.本发明所要解决的技术问题是克服上述现有技术存在的缺陷，提供一种基于声纹成像的gis故障诊断方法及系统，其适用于干扰噪声较强的声信号采集环境，利用声源可视化原理，提取gis辐射出的声信号并进行故障诊断。
6.为实现上述目的，本发明采用如下技术方案：一种基于声纹成像的gis故障诊断方法，其包括以下步骤：
7.s1、采集gis辐射出的声信号；
8.s2、对所述步骤s1中采集的gis辐射出的声信号进行时频分析和特征提取，对异常声源的位置进行定位；
9.s3、对采集到的声信号进行emd分解，并求取分解后imf分量的一维谱熵，作为训练样本代入svm模型模块中得到故障类型；
10.s4、将步骤s2和s3得到的故障诊断结果进行可视化处理，并同时进行存储和报警提醒。
11.进一步地，所述步骤s1中，采集gis辐射出的声信号所用的设备为传声器阵列和数据采集仪，两者之间采用数据线连接；传声器阵列获取gis辐射出的声信号采用的是远场声源信号模型：
[0012][0013]
式中，r为声源与传声器阵列中心的距离，l为各传声器之间的间距，λ为声信号的
波长；
[0014]
根据传声器的分布位置不同，不同位置上的传声器接收到的信号存在一定的时间延迟，以第一个传声器作为参考，则声信号到达其他传声器的时间相对于第一个传声器的时间延迟为：
[0015]
τm＝d cosθ/c，m＝1,2,λ,m，
[0016]
式中，c为声速，d为相邻两个传声器之间的距离，θ为远场声源的入射角度，τm为时间延迟，m为阵元数；
[0017]
阵列波束形成的输出为：
[0018][0019]
式中，wm为加权系数，pm为所对应阵元接收到的复声压信号，t为时间；
[0020]
用矢量表示上式可得：
[0021][0022]
式中，x(t)为阵列的m
×
1维快拍数据矢量，w(θ)为各传声器的方向矢量，p(t)为阵列接受的空间信号源矢量，为阵元m所对应的输出。
[0023]
进一步地，所述声信号的emd分解包括以下步骤：
[0024]
a、令原始信号为x(t)，为获取信号的上下包络线，首先找出所需的极大值点和极小值点，对这两类点用插值的方式分别拟合原始信号的波形，构造出上包络线和下包络线；
[0025]
b、计算上下包络线的平均值m1，并计算h1＝x(t)-m1，若h1满足imf分量的条件，则令h1为imf1，并进入下一个循环；
[0026]
c、若h1不满足条件，将h1替代原始信号x(t)，并重复步骤a和步骤b，得到h2，并判断新生成的h2是否满足上述条件，若还未满足条件，继续循环步骤a和步骤b，直到获得满足条件的imf1；
[0027]
d、若h1为imf1，则为计算下一个imf分量，需要将h1从原始信号x(t)中分离出来，得到可以继续分解的信号：
[0028]
r1＝x(t)-imf1；
[0029]
e、将r1替代需要分解的信号，重复循环步骤a-d，循环至满足结束条件，最终得到n个满足条件的imf分量和残差rn，信号x(t)表示为：
[0030][0031]
计算各个imf分量的一维谱熵：
[0032][0033]
式中，q(j)为第j个功率谱在整个imf中所占的百分比；n为imf包含的谱线条数；
[0034]
f、最后将这些imf分量的一维谱熵作为样本代入svm模型模块中训练样本，得到具体故障类型的诊断结果。
[0035]
更进一步地，步骤e中，所述的结束条件为rn的极值点数小于两个。
[0036]
更进一步地，分解得到的imf分量需要满足以下两个条件：
[0037]
条件一、在所需要分解的信号中，极点数和过零点数的差值需要控制在一个以内；
[0038]
条件二、在任意时刻，信号中由极大值点插值所构建的上包络线和极小值点所构建的下包络线均值为零。
[0039]
本发明采用的另一种技术方案为：一种基于声纹成像的gis故障诊断系统，其包括：
[0040]
采集模块，用于采集gis辐射出的声信号；
[0041]
成像定位模块，用于对采集的gis辐射出的声信号进行时频分析和特征提取，定位异常声源位置；
[0042]
故障诊断模块，用于对采集到的声信号进行emd分解，并求取分解后imf分量的一维谱熵，作为训练样本代入svm模型模块中得到诊断结果；
[0043]
pc模块，将成像定位模块和故障诊断模块得到的故障诊断结果进行可视化处理，并将故障诊断结果存储于存储器中，若故障诊断模块中判断gis设备出现故障，会触发pc模块中的报警器。
[0044]
进一步地，所述采集模块包括传声器阵列和数据采集仪，两者之间采用数据线连接。
[0045]
进一步地，所述成像定位模块包括时频分析模块和声源定位模块，声源定位模块用于定位异常声源位置，时频分析模块用于对采集的gis辐射出的声信号进行时频分析。
[0046]
进一步地，所述故障诊断模块包括特征提取模块和svm模型模块；所述的特征提取模块对采集到的声信号进行emd分解，并求取分解后imf分量的一维谱熵，所述的一维谱熵作为训练样本代入svm模型模块中，得到诊断结果。
[0047]
进一步地，所述pc模块包括显示屏、存储器和报警器。
[0048]
本发明具有以下有益效果：本发明通过声成像方法定位异常声信号，再对采集到的gis声信号做emd自适应分解，并计算各个imf分量的一维谱熵作为特征元素代入svm中，能够对gis的多种故障类型进行一次性的诊断判别，简单快捷，可以避免故障源辐射出的声信号被其他干扰源所覆盖，一定程度上可以提高传声器阵列对期望方向的声学特征采集效果，避免影响后续的故障诊断，且为判断gis的运行状态、评估gis的故障情况以及保障gis安全稳定运行提供了新方法。
附图说明
[0049]
图1为本发明一种基于声纹成像的gis故障诊断方法的流程示意图；
[0050]
图2为本发明一种基于声纹成像的gis故障诊断系统的结构示意图。
[0051]
图中：1、采集模块；11、传声器阵列；12、数据采集仪；2、成像定位模块；21、时频分析模块；22、声源定位模块；3、故障诊断模块；31、特征提取模块；32、svm模型模块；4、pc模块；41、显示屏；42、存储器；43、报警器。
具体实施方式
[0052]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
[0053]
实施例1
[0054]
参照图1，一种基于声纹成像的gis故障诊断方法，其包括以下步骤：
[0055]
s1、采集gis辐射出的声信号；
[0056]
s2、对步骤s1中采集的gis辐射出的声信号进行时频分析和特征提取，并对异常声源的位置进行定位；
[0057]
s3、对采集到的声信号进行emd分解，并求取分解后imf分量的一维谱熵作为训练样本代入svm模型模块中得到故障类型；
[0058]
s4、将步骤s2和s3得到的故障诊断结果进行可视化处理，并同时进行存储和报警提醒。
[0059]
svm(support vector machine，支持向量机)是机器学习算法中一种典型的监督学习方法，实现该方法的前提是需准备大量的训练样本，并对样本进行特征提取得到特征向量，将特征向量数据和样本的类型索引作为输入，训练得到svm模型模块。
[0060]
svm模型模块的训练步骤如下：
[0061]
步骤31、用分类器将所有样本集的特征分类；
[0062]
步骤32、根据选定的核函数分别计算每个分类器中每个特征向量的有关值；
[0063]
步骤33、根据这些特征相关值计算协方差矩阵空间；
[0064]
步骤34、对这个协方差矩阵空间进行householder变换；
[0065]
步骤35、计算特征系数，得到协方差矩阵以及其对应的超平面矩阵，根据这两个矩阵分别计算每个特征的特征系数，并将特征系数对协方差矩阵进行缩放；
[0066]
步骤36、获得模型参数；
[0067]
步骤37、根据训练得到的模型参数对某个测试集进行计算，验证该模型分类的准确性。
[0068]
在具体实施的过程中，选择正常运行的gis与已知故障类型的gis声信号作为样本集，并以3：1的比例划分训练样本和测试样本。将信号的imf分量的一维谱熵作为特征分量并与相关的分类器相对应，核函数采用rbf核函数进行分类，惩罚因子c＝1000，通过svm模型模块找到分类的最佳条件，即获取准确度较高(故障诊断正确率高)的超平面,从而在声纹诊断的过程中，根据此超平面去判断gis是否出现故障，提升了诊断的准确度和速度。
[0069]
s1中，采集gis辐射出的声信号所用的设备为传声器阵列和数据采集仪，两者之间采用数据线连接；传声器阵列获取gis辐射出的声信号采用的是远场声源信号模型：
[0070][0071]
式中，r为声源与传声器阵列中心的距离，l为各传声器之间的间距，λ为声信号的波长。
[0072]
根据传声器的分布位置不同，不同位置上的传声器接收到的信号存在一定的时间延迟，以第一个传声器作为参考，则声信号到达其他传声器的时间相对于第一个传声器的时间延迟为：
[0073]
τm＝d cosθ/c，m＝1,2,λ,m，
[0074]
式中，c为声速，d为相邻两个传声器之间的距离，θ为远场声源的入射角度，τm为时间延迟，m为阵元数；
[0075]
阵列波束形成的输出为：
[0076][0077]
式中，wm为加权系数，pm为所对应阵元接收到的复声压信号，t为时间；
[0078]
用矢量表示上式可得：
[0079][0080]
式中，x(t)为阵列的m
×
1维快拍数据矢量，w(θ)为各传声器的方向矢量，p(t)为阵列接受的空间信号源矢量，为阵元m所对应的输出。
[0081]
声信号的emd分解包括以下步骤：
[0082]
a、令原始信号为x(t)，为获取信号的上下包络线首先找出所需的极大值点和极小值点，对这两类点用插值的方式分别拟合原始信号的波形，构造出上包络线和下包络线；
[0083]
b、计算上下包络线的平均值m1，并计算h1＝x(t)-m1，若h1满足imf分量的条件，则令h1为imf1，并进入下一个循环；
[0084]
c、若h1不满足条件，将h1替代原始信号x(t)，并重复步骤a和步骤b，得到h2，并判断新生成的h2是否满足上述条件，若还未满足条件，继续循环步骤a和步骤b，直到获得满足条件的imf1；
[0085]
d、若h1为imf1，则为了计算下一个imf分量，需要将h1从原始信号x(t)中分离出来，得到可以继续分解的信号：
[0086]
r1＝x(t)-imf1；
[0087]
e、将r1替代需要分解的信号，重复循环步骤a-d，循环至满足结束条件(rn的极值点数小于两个)，最终得到n个满足条件的imf分量和残差rn，信号x(t)可以表示为：
[0088][0089]
首先对信号进行emd分解，分解得到n个imf分量，计算各个imf分量的一维谱熵：
[0090][0091]
式中，q(j)为第j个功率谱在整个imf中所占的百分比；n为imf包含的谱线条数。
[0092]
f、最后将这些imf分量的一维谱熵作为样本代入svm模型模块中训练样本，得到具体故障类型的诊断结果。
[0093]
步骤e中，所述的结束条件为rn的极值点数小于两个。
[0094]
分解得到的imf分量需要满足以下两个条件：
[0095]
条件一、在所需要分解的信号中，极点数和过零点数的差值需要控制在一个以内；
[0096]
条件二、在任意时刻，信号中由极大值点插值所构建的上包络线和极小值点所构建的下包络线均值为零。
[0097]
gis辐射出的声信号依次进行阵列信号处理、时频分析处理、特征提取处理和特征定位处理后，得到信号的故障定位。
[0098]
gis辐射出的声信号依次进行经过单通道声音信号处理、经验模式分解处理、一维
谱熵向量处理和svm模型模块故障诊断处理，得到信号的故障类型。
[0099]
实施例2
[0100]
参照图2，一种基于声纹成像的gis故障诊断的系统，包括：
[0101]
采集模块1，用于采集gis辐射出的声信号；
[0102]
成像定位模块2，用于对采集的gis辐射出的声信号进行时频分析和特征提取，定位异常声源位置；
[0103]
故障诊断模块3，用于对采集到的声信号进行emd分解，并求取分解后imf分量的一维谱熵，作为训练样本代入svm模型模块中得到故障类型；
[0104]
pc模块4包括显示屏41、存储器42和报警器43，将诊断结果进行可视化处理，并将成像定位模块2和故障诊断模块3得到的故障诊断结果存储于存储器中，若故障诊断模块中判断gis设备出现故障，会触发pc模块中的报警器。
[0105]
具体地，所述采集模块1包括传声器阵列11和数据采集仪12，两者之间采用数据线连接。
[0106]
具体地，所述成像定位模块2包括时频分析模块21和声源定位模块22，声源定位模块22用于定位异常声源位置，时频分析模块21用于对采集的gis辐射出的声信号进行时频分析。
[0107]
所述故障诊断模块3包括特征提取模块31和svm模型模块32；所述的特征提取模块31对采集到的声信号进行emd分解，并求取分解后imf分量的一维谱熵，所述的一维谱熵作为训练样本代入svm模型模块32中，得到诊断结果。
[0108]
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。