一种基于分形的超声相控阵无损检测方法
【技术领域】
[0001] 本发明属于现代超声信号处理领域，特别涉及一种基于分形的超声相控阵无损检 测方法。
【背景技术】
[0002] 超声相控阵技术在上世纪八十年代被提出，初期主要是应用于医学超声成像诊 断。近年来国内外超声相控阵技术发展很快，在医学诊断和工业检测等方面的研宄很热门。 随着超声相控阵系统设计、系统模拟以及实际应用的不断发展，相控阵技术慢慢走上数字 化。随着超声相控阵技术开发和应用的不断发展，相控阵技术已被用于复合材料等的无损 检测。然而，由于复合材结构多种多样，使用要求也不尽相同。由于超声相控阵各阵元得到 的回波信息具有不同的信号特征，单个信息不能够全面描述被测体；超声波在复合材料中 幅值衰减大，所以超声阵列得到的回波幅度很小，后续信息不容易处理；当多个阵元同时发 射时，各阵元的发射波、回波之间会发生干涉作用，当各波相位不同时，信号强度会减弱。所 以存在以下难点：仅仅利用超声相控阵方法还难以胜任其质量的检测与评定，实际检测工 作中往往需要针对不同检测对象和要求，采用不同的检测技术和方法。
[0003] 根据上述技术难点，中国发明专利申请号201310152221. 6通过只采用一种线性 探头，运用超声相控阵技术可对不同形状及尺寸的轴类进行无损检测，提高了检测效率，解 决了轴类传统超声波检测中检测不同的轴需要不同型号的探头，并配备不同的轴类对比 试块进行对比试验的问题，大大减少了人力和物力的投入。然而，此发明只是针对轴类对 象，而对于平面的复合材料并不适用。
[0004] 中国发明专利申请号201310166560.X提出一种钢箱梁U肋角焊缝的超声相控阵 检测方法，本发明方法解决了钢箱梁U肋角焊缝熔深和内部缺陷检测的难题，使钢箱梁U 肋角焊缝的焊接质量得到有效控制。可是，此发明不仅需要超声相控阵探头还需要对比试 块，而且灵敏度的调节在试块上，不好操作，灵敏度误差大。
[0005] 国家发明专利申请号201310209297. 8提出一种基于超声相控阵和支持向量机的 焊缝缺陷检测识别方法，解决了厚板电子束焊缝缺陷无损检测与识别问题。然而，此发明需 要构造三个基于支持向量机的子模型，只要一个子模型出现错误，就会影响整个识别结果。
[0006] 中国发明专利申请号200510124212. 1提出使用相控阵超声自动检测系统对焊缝 横向裂纹进行检测的方法，解决了检测焊缝横向裂纹的难题，且检测结果准确可靠。可是， 此发明要使用四支超声相控阵探头一焊缝为对称轴呈十字对称排列，所以进行检测前必须 要知道横向焊缝的位置，而对于一般要检测的对象，我们事先并不知道焊缝在哪里。要实现 对未知缺陷情况的复合材料进行无损检测，通过硬件实现成本较高，实现比较困难。
[0007] 在近一个世纪内，信号处理研宄领域不断发展，人们已开始把研宄的焦点从传统 的信号处理到现代信号处理上来。现代信号主要是非线性和非平稳性等复杂信号。分形理 论就是处理现代信号的一种手段，并为人们解决非线性问题开辟了新的途径。小波变换与 分形算法具有深刻的内在联系。分形是一种几何语言，小波变换是一种分析工具，小波变换 的平移不变性与分形变换的尺度不变性相一致。分形的定义是不规则几何体在动力学演 化过程中，在一定的标度尺度范围内，相应的测度不会随尺度改变而发生变化，也就是系统 或整体中每一个元素或局部都能在一定程度上体现整个系统的特性和信息，这称为自相似 性。所以，通过从大到小不同尺度的变换，可以在越来越小的尺度上观察到越来越丰富的细 T。
[0008] 把分形算法和贝叶斯信息融合算法结合起来用于超声相控阵的无损检测，既可以 充分利用相控阵各阵元得到的信息，又可以快速准确的判别多种被测对象的未知缺陷。而 目前尚未有文献记载此类技术。
【发明内容】
[0009] 为解决上述技术难点，本发明的目的是提供一种基于分形的超声相控阵无损检测 的方法。该方法不仅能够提高各超声相控阵阵元所采集信息的利用率，得到对现实环境更 为准确、可靠的描述，还可以快速准确的判别复合材料的未知缺陷。
[0010] 本发明是通过下述技术方案加以实现的，一种基于分形的超声相控阵无损检测方 法，包括以下步骤：
[0011] (一）采用超声相控阵设备，设计设备参数让超声相控阵探头的各个阵元依次单 独发射，1号阵元单独发射时，2号阵元采集回波信号；2号阵元单独发射时，3号阵元采集回 波信号，以此类推，最后将采集到的所有信号组成一个数组，得到多个数组；
[0012] (二）对得到的数组采用分位图法，得到参加信息融合的最少融合数据，组成最佳 融合新数组；
[0013] (三）得到最佳融合新数组后进行最后的信息融合，从而得到最有结果，贝叶斯估 计的建立：
[0014] 1)求出每个最佳融合新数组m个元素彼此之间的置信距离测度duaj= 1，2,. . .，n)，组成置信矩阵D，
【主权项】
1. 一种基于分形的超声相控阵无损检测方法，包括以下步骤： (一） 采用超声相控阵设备，设计设备参数让超声相控阵探头的各个阵元依次单独发 射，1号阵元单独发射时，2号阵元采集回波信号；2号阵元单独发射时，3号阵元采集回波信 号，以此类推，最后将采集到的所有信号组成一个数组，得到多个数组； (二） 对得到的数组采用分位图法，得到参加信息融合的最少融合数据，组成最佳融合 新数组； 得到最佳融合新数组后进行最后的信息融合，从而得到最有结果，贝叶斯估计的建 立： 1) 经过分位图法可以得到新的数组，求出每个新数组的m个元素彼此之间的置信距离 测度(Iij (i, j = 1，2,. . .，η)，组成置信矩阵D，
PiUIxi)是在第i个样本值取Xi的条件下的概率密度，σ i是第i个样本值的方差；
4 经过数学推导得， 其中，Φ (b)是标准正态分布的概率，b表 > 示变量t的上限； 2) 求出关系矩阵R，R是由r。组成，
其中，ru是第i个样本与第j个样本的相关系数，β u是置信距离测度du的界线值； 3) 由关系矩阵R，得到每个数组参加信息融合的最佳融合数据集； 4) 依新数组各元素的pdf (概率密度）曲线作为各特征函数，记为p(X); 5) 求出决策函数(Kx1, X2, .·.，xn)，其中X1, X2, ...，xn是来自总体X的样本值； 6) 由决策函数求出风险函数R(u|d) =Eu{L(u，d(x1，x2，...，xn))}，其中u是参数，即所 要求的融合结果，L(u，d(X 1, x2, ...，xn))是损失函数的二次式； 7) 对风险函数求平均，得到决策函数(1(11，12，...，1)的贝叶斯风险艮其中，
h (u)是参数u的先验分布密度； 8) 若有 cf (X1, x2,. . .，xn)使得 B (cf) = min {B (d)}，d e Φ，则称 cf为参数 u 的贝叶斯 估计量，也称为最优估计； 9) 损失函数取二次式L (u, d) = [U-CKx1, X2, ...，xn)]2,则u的贝叶斯估计量为
要想求出估计量，只要求出P (U I (X1, X2, . . .，Xn))即可；其中，
其中，α是正则化因子， 贝ljp(u| H-.'XJ)
其中〇k是第k个样本的标准差，σ ^是所有样本 9 的整体标准差，Utl是样本均值；由于我们采集的数据波形并不是标准的正态分布，特征函数 P(X)很难直接求出，所以我们采用以下步骤推导出所要求的最优估计； 10) 假设P (u I (Xl，x2,...，Xn))服从正态分布"(〃、.σ()，其中 Un是正态分布的均值， 2是正态分布的方差；所以
由以上两式 比较，得
11) 通过以上推导，得到参数U的贝叶斯估计为也 就是超声相控阵阵元得到的最优息融合结果；
其特征在于，还包括： (四）分形过程的建立
1) 得到贝叶斯信息融合算法的最优融合结果，共有M个点，画出M个最有融合点的波形 图； 2) 选取标度网格长度δ (只能取整数，且不小于1)，横坐标被分成 .个网格； 3) 在横坐标的第i区间内，波形最高值减去波形最低值，再除以δ，即使该区间内波形 与网格相交的盒子数目，若余数不等于零，则盒子数加 1 ; 4) 将所有区间的盒子数目相加，得到盒子总数化（F); 5) 取不同的δ，重复上述2到4,得到相应的盒子总数Ns (F); 6) 画出Ns (F)和δ的双对数曲线； 7) 取出曲线线性关系好的一段，用最小二乘法做一元线性回归拟合出直线的斜率，作 为盒维数的估计值Dim; 8) 根据盒维数的估计值Dim，判断复合材料有无缺陷以及缺陷类型。
【专利摘要】本发明公开了一种基于分形的超声相控阵无损检测的方法，方法包括贝叶斯信息融合算法和信号分形处理两部分，其中，贝叶斯信息融合算法有以下步骤：取出超声相控阵探头各阵元依次单独采集信号的中位数和上、下四分位数，得四分位离散度，求淘汰点，得到有效数据融合集，求各有效数据的特征函数和决策函数，得到风险函数，提取决策函数的贝叶斯风险，得出参数的最优估计值；信号分型处理有以下步骤：通过小波分解系数图分析图像分形特性，确定无标度区，计算分形维数，根据分形维数判断有无缺陷以及缺陷类型，在使用贝叶斯信息融合算法来提高信息利用率的基础上，本发明使用分形技术处理非线性的现代信号，从而判断出复合材料有无缺陷以及缺陷的类型。
【IPC分类】G01N29-44
【公开号】CN104535663
【申请号】CN201410632519
【发明人】詹湘琳, 蔡玉杰, 刘岱, 刘涛
【申请人】中国民航大学
【公开日】2015年4月22日
【申请日】2014年11月11日