钢板缺陷识别的Volterra-SVM模型研究_邓勇.pdf

年第 卷 月第 期机 械 科 学 与 技 术 ：收稿日期：基金项目：四川省科技支撑计划项目（）作者简介：邓勇（），高级工程师，硕士生导师，研究方向为模式识别和超声无损检测，邓勇，黄远伟，赖治屹钢板缺陷识别的 模型研究机械科学与技术，（）：钢板缺陷识别的 模型研究邓勇，黄远伟，赖治屹（西南石油大学 机电工程学院，成都；西南油气田分公司输气管理处，成都）摘要：针对钢板缺陷识别问题，结合超声波脉冲反射原理，提出一种基于 级数和支持向量机的钢板缺陷识别方法。首先，利用 级数模型建立起钢板缺陷的特征模型；其次，使用分数阶粒子群优化算法提取出原始信号中的特征参数，即 级数时域核；最后，将提取到的特征向量输入支持向量机模型进行训练与测试，完成对钢板缺陷的分类识别。设计实验得到多组数据样本，进行模型验证，实验结果表明：基于 级数和支持向量机的识别模型能够较好的完成对钢板缺陷的分类识别，识别准确率达。关键词：缺陷识别；级数；分数阶粒子群优化算法；支持向量机中图分类号：；文献标志码：文章编号：（），（，；，）：，（），（），：；无损检测技术因其不破坏被检件本身结构的特点成为工业发展必不可少的有效工具。超声检测作为常用无损检测技术之一，具有扫描速度快、分辨率高和缺陷量化准确等优点，被广泛应用于钢板的缺陷检测，近年来受到了越来越多的关注。随着机器学习技术的发展，利用信号分析方法处理分析超声信号，从超声信号波形图中提取特征，进而对缺陷进行识别已经成为解决缺陷识别问题的一个发展方向。近年来，国内外的学者对各类信号处理方法应用于缺陷识别的问题做了大量的研究，等使用经验模态分解方法得到原始信号本征模态函数作为信号特征，利用谱聚类和粗糙集属性约简进行特第 期 邓勇，等：钢板缺陷识别的 模型研究：征选择，输入神经网络完成对缺陷信号的分类；等提出一种基于深度学习的卷积神经网络来对缺陷进行分类热编码的方法；等使用小波包分解提取缺陷超声信号的特征，利用变尺度混沌遗传算法去除冗余和无关特征，再输入支持向量机模型完成对缺陷的分类。潘峰等采用支持向量机和主成分分析法消除特征向量中的冗余信息完成对腐蚀缺陷深度信号分类；车红昆等提出一种将支持向量机与贝叶斯推理相结合的多特征融合识别算法，对缺陷超声信号的识别有着较好的效果。目前，在超声波缺陷信号识别研究中，特征提取方法常见的有傅里叶变换（）、经验模态分解（）、小波包分解（）以及提取回波信号时频域无量纲参数作为特征等。上述方法大多只能完成对缺陷的定性识别，对于缺陷的大小、形状、深度的识别还需要进一步研究。除上述方法外，级数模型也广泛应用于信号特征提取研究中，它适合处理非线性、非平稳信号。超声检测钢板缺陷信号是典型的非线性、非平稳信号，使用 级数模型处理这种信号具有独特的优势。因此可以将 级数时域核辨识方法应用于钢板缺陷信号的特征提取。目前 级数模型辨识方法主要有最小二乘法，遗传算法，粒子群算法等。最小二乘法辨识时需要目标函数连续可导；遗传算法具有交叉、变异的特点，算法不够简洁。基于此，有必要寻求新的辨识算法。分数阶粒子群优化算法（）是基于群体的演化算法，将粒子群算法推广到分数阶，使用分数阶导数控制粒子群算法的收敛速度，克服传统粒子群迭代速度会爆炸式增长的缺陷。常用信号分类方法有神经网络、隐马尔科夫模型和支持向量机等。支持向量机对比其他方法，更适用于小样本条件下的分类问题，具有很好的泛化能力，能够解决实验获取的各缺陷的样本数量较少的问题，对于钢板缺陷信号的分类问题有着良好的适应性，并且能够很好地解决非线性、高维数和局部极小点等实际问题。基于对上述问题的分析，本文提出将 级数模型、分数阶粒子群优化算法（）和支持向量机（）相结合的识别方法，实现对钢板缺陷的分类识别。首先，利用 级数模型结合超声脉冲反射原理建立起缺陷信号的特征模型；其次，利用分数阶粒子群优化算法对建立起的 级数模型进行时域核辨识，提取出前三阶 级数时域核作为缺陷信号的特征向量；将提取出的特征向量输入训练好的多分类支持向量机模型，完成对钢板缺陷的量化识别。超声波脉冲反射法检测缺陷超声波脉冲反射法是检测精度较高的无损检测方法，利用超声波的反射特性，通过对需检测的对象发射超声波脉冲，根据有无缺陷回波或工作底面反射回波、回波幅度的大小、回波信号数量、回波在示波频时基线上的位置以及回波包络线形状变化来判断被检件是否存在缺陷，常采用收发共用的单探头。为满足进一步研究自动化检测需求，本实验所用超声换能器为超声探头轮，其原理如图 所示。超声换能器产生激发波，穿过耦合液和耦合固体到达工件表面，遇到缺陷产生缺陷回波，超声换能器接收到缺陷回波，得到超声信号，完成缺陷检测。图 脉冲反射原理示意图图 中，信号 为超声激励信号，、分别是界面、处超声回波信号，若钢板表面存在缺陷，则代表缺陷的超声回波信号，信号 为钢板底面回波信号。可从图 中得出，声波在多层非线性介质中传播，使得波形复杂，人工难以辨识，需要运用信号分析方法对原始信号进行处理，本文使用 级数核辨识的方法从原始信号中提取特征参数。原始信号特征提取 级数模型 级数是一种泛函级数，对于非线性系统而言，它能表征系统的本质特征，且不随输入输出信号改变。当系统内部发生变化，其时域核发生变化。超声缺陷检测系统可看作非线性系统，只有当系统内部结构发生变化时（出现缺陷），其 级数时域核发生变化，故可以将 级数各阶时域核作为超声波缺陷检测信号的特征。设超声激励信号为（），作为系统输入信号；回波信号为（），作为系统输出信号，这种非线性系统模型可以由 级数表示，即（）（）（）（）其离散形式可表示为机 械 科 学 与 技 术第 卷：（）（，）（）（）（）式中：（，）为 阶离散 时域级数核；为阶次；为 阶 级数记忆长度；（）为误差。当 和 取适当值时，（）可充分小。由于 级数存在维数灾难的难题，为了便于计算，分析实际系统常用离散系统的前 阶 级数进行近似描述。设系统输入向量（）（），（），（），为数据长度，定义输入矩阵为（），（），（），（）（）定义系统输出向量为（）（），（），（）（）定义 核向量为 （），（），（，），（，）（）根据 级数对称性原理，前 阶 核向量 的维数 可以表示为，（）（）（）（）由式（）式（）可得（）（）在已知系统输入、输出的基础上，使用分数阶粒子群优化算法求解式（），就可以得到 级数核向量，其核向量就作为缺陷信号的特征向量。分数阶粒子群分数微分概念揭示的一个重要性质是分数导数隐含着对所有过去事件的“记忆”。分数阶粒子群粒子速度由“记忆”中的粒子速度迭代得到，可通过分数阶次 加以控制。设粒子群规模为，粒子维数为，第 个粒子当前位置为，当前速度，第 个粒子找到的具有最佳适应度的位置为，称为局部最优点，其中最优者记为，称为全局最优点。，通过计算目标函数（）确定，粒子更新公式描述为：（）（）（）（）（）（）（）（）（）（）（）（）（）式中：，称为学习因子；，是两个独立的服从均匀分布的随机数。考虑到分数阶微积分观点，速度导数的阶数可以推广到 的实数，可使得粒子群速度迭代变换更平滑，记忆效应更长。设截断阶数 ，式（）可以被改写为（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）对于 核辨识来说，算法的目的就是寻求一组 维的向量，使得目标函数最小，目标函数为适应度值函数，即 （）（）通过对式（）进行更新迭代，搜索出适应度值最小的全局最优解，即得到 时域核向量。按照同样方法对多组缺陷超声信号求取前 阶 时域核，建立特征样本集，并输入支持向量机模型进行训练与测试。缺陷特征分类 模型是建立在统计学的 维理论和结构化风险最小的原理基础上的分类方法，能够很好地解决小样本、非线性、高维数和局部极小点等实际问题。根据钢板缺陷的大小、形状和深度对其编号，分别对应支持向量机的几个状态。支持向量机的分类就是寻找最优分类面 ，使不同类型样本之间分类间隔最大，利用 优化方法可以把上述最优分类面问题转化为其对偶问题，即：（），（）其中 是每个样本对应的 乘子。对于非线性问题，可以通过非线性变换将低维空间的非线性问题转换为高维空间的线性问题。在上诉对偶问题中，在高维空间实际上只需要进行内积计算，根据泛函理论，只要一种核函数 （，）满足 条件，它就对应某一变换空间的内积。目前研究最多的核函数主要有 类，线性核函数（，），多项式核函数（，）（），径向基核函数（，）。利用核函数将目标函数改写为（）（）（）第 期 邓勇，等：钢板缺陷识别的 模型研究：通过式（）建立样本决策函数：（）（，）（，）（）利用样本决策函数对特征样本集进行分类识别。实验结果与分析为获取缺陷超声信号原始波形数据，在 厚钢板设置如图 所示人工缺陷，并设计如图 所示实验平台对 厚带人工缺陷钢板进行超声检测。本实验中，检测的工件为钢板，其探测频率为，为了满足检测分辨率和检测灵敏度要求，探头选用中心频率.，收发一体的超声纵波直探头；为满足进一步研究自动化检测需求，设计如图 所示探头轮；设计有超声激发电路与接收电路，并设计 线接口便于超声信号的采集与显示；超声回波原始信号由 型数字示波器进行显示与采集并存储，其采样频率为 。图 部分人工缺陷示意图图 超声检测实验装置图实验流程如下：）超声激发电路发射脉冲信号，加载至超声换能器，产生激发波。）激发波穿过耦合液和耦合固体到达工件表面，若遇到缺陷则产生缺陷回波，超声换能器接收到缺陷回波信号。）激励信号与回波信号通过同轴线传输至数字示波器，进行采集、显示并存储。）利用分数阶粒子群算法对原始信号进行特征提取，提取出各信号的前 阶 级数时域核，用于建立缺陷特征样本集，缺陷参数如表 所示。表 不同种类缺陷序号大小 形状深度 圆形圆形矩形圆形圆形矩形不规则形不规则形无缺陷）对 模型进行训练，用训练好的 模型对测试集样本进行分类识别。实际工程应用中，钢板缺陷往往是缓变的，为了尽可能的满足实际检测需求，实验中设置的人工缺陷在形状、大小和深度上都按照一定规则均匀取值，并设置有一部分不规则缺陷。然而由于加工工艺限制，人工缺陷参数取值间隔较大，不能完全模拟实际缓变缺陷，减小缺陷参数取值间隔，提高缺陷识别精度，也是该研究进一步发展的方向。利用如图 所示实验装置对表 所示 种不同大小、不同形状、不同深度缺陷进行超声检测，采集原始超声信号，并引入无缺陷类型用于对照，对每个类型缺陷独立重复采样 次，共 组数据样本。部分超声回波数据原始信号如图 所示。机 械 科 学 与 技 术第 卷：图 原始信号图得到原始超声信号后，根据 级数模型理论确定每个缺陷特征向量的维数，由式（）可得，设：阶 时域核的记忆长度为，其长度为；阶 时域核的记忆长度为，其长度为；阶 时域核的记忆长度为，其长度为。缺陷特征向量，即前 阶 级数核，维数为，可表示为：（，）（），（），（）（，），（，），（，）（，），（，），（，）根据式（），利用分数阶粒子群算法对特征向量 进行求解。设粒子群规模为，粒子维数为，学习因子 .，.，迭代次数为 次，分数阶参数.，根据大量实验得出，取.时，算法收敛过程相对稳定，其适应度值曲线如图 所示。图 适应度值曲线由图 可看出，分数阶粒子群算法收敛速度快，迭代接近 次后完成对参数的寻优，且波动过程较稳定，不易陷入局部最优，其最大的优势在于其固有的记忆特性，它包含了更重要的过去的信息，对于解决非线性问题具有更好的适应性。利用分数阶粒子群算法提取出 组样本的 级数时域核参数，并在 中使用 函数绘制出部分缺陷类型特征向量的三维曲面图，其三维分别是：阶 级数时域核，阶 级数时域核，阶 级数时域核（维数不足以 补齐），如图 所示，不同缺陷特征有着较为明显的区分度。图 部分缺陷类型的三维曲面图第 期 邓勇，等：钢板缺陷识别的 模型研究：特征提取完成后，得到 个 维特征向量，根据缺陷类型对其编号，建立 的特征矩阵作为支持向量机模型分类识别的特征样本集。应用 工具箱进行编程，使用径向基核函数作为支持向量机的基本核函数。以样本总数的 （个）作为训练集，样本总数的（个）作为测试集，将支持向量机模型训练好后，