一种基于多尺度形态分解能谱熵和支持向量机的泄漏声发射信号识别方法
【技术领域】
[0001] 本发明涉及一种基于多尺度形态分解能谱熵和支持向量机的泄漏声发射信号识 别方法，属于声发射信号模式识别领域。
【背景技术】
[0002] 在工业快速发展的今天，各类压力管道、高压锅炉随处可见，在使用过程中由于腐 蚀、磨损等原因可能会造成管道或炉壁材料破损导致泄漏，泄漏检测是石油、化工、天然气 和城市供水等领域中需要解决的重要问题之一。研究泄漏信号的检测理论、检测方法，实现 对泄漏信号的分类识别，对于维护管道、锅炉等的安全运行，避免资源浪费，有着重要的理 论价值和现实意义。声发射检测技术是一种利用局部材料快速释放能量所产生的瞬态弹性 波作为激励源的检测技术，在无损检测中占有重要的地位。目前声发射检测在刀具磨损、腐 蚀检测等方面已经取得了一定的研究成果，所以将其应用于工业管道泄漏声发射信号检测 中也是可行的。
[0003] 目前，对声发射信号的处理方法主要有时域参数法、时域波形法、频域法、小波分 析和经验模态分解等。这些方法从不同角度描述声发射信号的特征，提供了多种多样的声 发射信号检测方法。虽然声发射检测领域的文献很多，但是目前管道泄漏的声发射技术研 究主要是对有无泄漏和泄漏点位置的判定，根据泄漏声发射信号研究泄漏信号参数状态的 文献不多。
[0004] 判断泄漏声发射类型是泄漏声发射检测的主要目的之一，泄漏声发射类型的确定 有助于进一步对材料受损情况作出评估，以便提前采取相应补救措施。因此，研究一种能在 多干扰和低信噪比的复杂信号中，快速、准确地提取泄漏信息特征并识别泄漏类型的方法 具有重要意义。
[0005] 数学形态学非线性信号处理方法是一门建立在集合论和积分几何学上的一种方 法，数学形态学在考察信号时使用结构元素探针，通过结构元素探针在信号中不断移动来 提取有用信息从而进行特征分析和描述。多尺度形态学以数学形态学分析为基础，直接由 波形集合形态采用不同尺度结构元素对波形进行变换，通过形态学膨胀、腐蚀、开、闭运算， 获得各尺度形态谱及所研究对象的直观理解。支持向量机(SVM)是基于统计学理论下的一 种新型机器学习方法，能很好的解决实际应用中出现的小样本、过学习、高维数、局部最小 等问题，特别是在小样本、非线性情况下，具有较好的泛化能力。
[0006] 基于此，本发明首先用实验室设备数字声发射系统采集敲击、砂纸和断铅三种模 拟声发射信号;然后进行多尺度形态学分解，提取每一尺度上高低频段谱能量，然后计算能 谱熵;再计算每一尺度所占能谱熵的比例，构造特征向量;最后利用支持向量机对特征参数 进行训练和测试。结果表明，本文提出的方法可以有效地对泄漏声发射信号进行分类，能在 泄漏故障发生时及时准确判断出泄漏故障类型，为采取相应的补救措施提供参考。
【发明内容】
[0007] 本发明提供了一种基于多尺度形态分解能谱熵和支持向量机的泄漏声发射信号 识别方法，以用于解决泄漏声发射信号检测分类问题。
[0008] 本发明的技术方案是:一种基于多尺度形态分解能谱熵和支持向量机的泄漏声发 射信号识别方法，首先采用数字声发射系统进行实验数据采集;对采集的模拟泄漏声发射 信号进行多尺度形态分解，分别计算其在不同尺度上的谱能量，并计算能谱熵值;再计算每 一尺度所占能谱熵的比例，将其组成特征向量;最后利用支持向量机对特征向量进行训练 和测试。
[0009] 所述方法的具体步骤如下：
[0010] stepl、声发射信号采集:通过数字声发射系统采集模拟泄漏声发射信号N组，信号 记为f(X);
[0011] Step2、结构元素采用扁平结构元素 g= [000]，对信号f(x)进行多尺度形态分解； 其中，分解尺度为n，则f (X)的多尺度形态分解信号为di(x)，d2(x)，…，di(x)，…，dn(x)，di (x)，d2(x)，···，di(x)，…，dn-ι(χ)为第1~n-1尺度上求得的高频信号，dn(x)为第η尺度上求 得的低频信号；
[0012] Step3、计算各个尺度相应频段上的谱能量Ei及多尺度形态分解能谱熵HE;其中， ;M为数据窗;T为采样周期，x = kT，k为采样序数， λ一η-ivi
- ▲
良示信号在第i尺度的谱能量EiA总能量
的概率；
[0013] Step4、在求取形态分解能谱熵的基础上，分别求取每一尺度所占能谱熵的比例 每组信号经过公式计算后得到不同尺度下的比例Pl(i = 1，2,…，n)，将 所有的Pi组合起来，得到1 X η维的特征向量矩阵P = [Pi，P2，…，Pn]，按照这个方法，可以求 取每一组泄漏声发射信号的特征向量；
[0014] Step5、取每类信号的N1组特征向量作为训练集，剩余N2组特征向量作为测试集， 进行SVM测试，得到分类识别结果;其中，N1+N2 = N。
[0015] 所述模拟泄漏声发射信号选取敲击、砂纸和断铅模拟泄漏声发射信号中的任意两 种或以上。
[0016]本发明的工作原理是：
[0017]第一步:实验采用北京声华兴业科技有限公司生产的SAEU2S数字声发射系统进行 试验数据采集，实验时，保持传感器与事故点位置不变，在事故点处分别模拟敲击、砂纸和 断铅三种工况声源，并记录相关实验数据，为验证上述方法的有效性，在实验室条件下，按 照国际无损检测界的规定(美国ASTM976号文件），使用0.5mmHB铅笔芯进行断铅试验，可以 得到模拟断铅声源信号，并用砂纸和金属棒，在相同位置模拟出砂纸和敲击两种不同工况 声源信号，选取敲击、砂纸和断铅模拟泄漏声发射信号中的任意两种或以上作为模拟泄漏 声发射信号。
[0018]第二步:结构元素采用最简单的扁平结构元素 g= [000]，结构元素宽度L取决于分 解尺度n，即L = 2n+1，对原始信号f(x)进行多尺度形态分解，设分解尺度为n，则f(x)的多尺 度形态分解信号为di(x)，d2(x)，…，di(x)，…，dn(x)，其中，di(x)，d2(x)，…，di(x)，…，d n-1 (x)为第1~n-1尺度上求得的高频信号，dn(x)为第η尺度上求得的低频信号；则有：
，其中计算方法如下： /=]_
[001 ?
[002C
[0021
[002：
[0023] 其中，符号Θ表示膨胀运算，符号?表示腐蚀运算，符号?表示闭运算，符号〇表示 开运算；
[0024] 第三步:各个尺度相应频段上的谱能量的计算公式如下：
式 中，Μ为数据窗;Τ为采样周期，x = kT，k为采样序数。那么多尺度形态分解能谱熵HE的计算式 为：
表示信号在第i尺度的谱能量EiA总能量
的概率；
[0025] 第四步:在求取形态分解能谱熵的基础上，分别求取每一尺度所占能谱熵的比例， 具体公式为:Pi = _Pilnpi/HE，其中，pi的定义同上，HE为多尺度形态分解能谱熵，这样每组 信号经过公式计算后得到不同尺度下的比例&(1 = 1，2,…，n)，将所有的^组合起来，得到1 Xn维的特征向量矩阵···'"]，按照这个方法，可以求取每一组泄漏声发射信号 的特征向量；
[0026]第五步:取每类信号的N1组特征向量作为训练集，剩余N2组特征向量作为测试集， 选取最优C、γ值，进行SVM测试，得到分类识别结果；C代表惩罚因子参数，γ代表gamma函 数，N1+N2 = N。由上可知，选择合适的核函数和分类器参数可以得到识别正确率很高的分类 结果。
[0027]本发明的有益效果是：
[0028] 1、通过将泄漏声发射信号的类型识别问题转化为模式分类问题，在惩罚因子和γ 参数选择合适的情况下，SVM对测试样本包含的声发射信息具有较强的分辨能力，可以同时 区分敲击、砂纸、断铅信号，且识别正确率很高，可实现自动分类；
[0029] 2、支持向量机对小样本数据具有较好的统计学习能力；
[0030] 3、通过求声发射信号的多尺度分解能谱熵值，并分别求取每一尺度所占能谱熵的 比例，能较好的反映声发射信号状态信息，可以将其作为声发射信号的特征向量。
【附图说明】
[0031] 图1为本发明的流程图；
[0032] 图2为本发明三种类型模拟泄漏声发射信号的时域波形图；
[0033] 图3为本发明中敲击信号多尺度形态分解图；
[0034] 图4为本发明中砂纸信号多尺度形态分解图；
[0035] 图5为本发明中断铅信号多尺度形态分解图；
[0036] 图6为本发明支持向量机的体系结构图；
[0037]图7为本发明SVM训练样本可视化图像。
【具体实施方式】
[0038]实施例1:如图1-7所示，一种基于多尺度形态分解能谱熵和支持向量机的泄漏声 发射信号识别方法，首先采用数字声发射系统进行实验数据采集;对采集的模拟泄漏声发 射信号进行多尺度形态分解，分别计算其在不同尺度上的谱能量，并计算能谱熵值;再计算 每一尺度所占能谱熵的比例，将其组成特征向量;最后利用支持向量机对特征向量进行训 练和测试。
[0039]所述方法的具体步骤如下：
[0040] Stepl、声发射信号采集:通过数字声发射系统采集模拟泄漏声发射信号N组，信号 记为f(X);
[0041] Step2、结构元素采用扁平结构元素 g= [000]，对信号f(x)进行多尺度形态分解； 其中，分解尺度为n，则f (X)的多尺度形态分解信号为di(x)，d2(x)，…，di(x)，…，dn(x)，di (x)，d2(x)，···，di(x)，…，dn-ι(