一种基于全聚焦成像算法改进边缘缺陷检测的方法与流程

本发明设计缺陷定位技术领域，特别是改进了一种常规超声不容易解决的小于半波长微小缺陷以及复杂结构试件的后处理算法。

背景技术：

超声相控阵技术逐渐应用于工业无损检测，特别是在核工业及航空工业等领域。如核电站主泵隔热板的检测；核废料罐电子束环焊缝的全自动检测及薄铝板摩擦焊缝热疲劳裂纹的检测。随着超声相控阵检测技术的不断普及，使得成像算法的研究越来越得到广泛的关注。其中HOMLES提出的全聚焦成像算法，使各阵元声束在检测区域内的划分的每一个点上进行聚焦，因此对于常规超声相控阵不容易解决的小于半波长的微小缺陷可以进行高精度识别。全聚焦成像提供了几种传播模式(直接或间接成像)。对于传统的全聚焦算法来说，一般直接采用直接成像，而实际检测中常常会由地面回波反射到缺陷上间接成像。

技术实现要素：

本发明的目的在于针对现有技术的不足，提供一种基于全聚焦成像算法改进边缘缺陷检测的方法。

本发明的目的是通过以下技术方案来实现的：一种基于全聚焦成像算法改进边缘缺陷检测的方法，包括以下步骤：

(1)采集超声阵列全矩阵数据：依次激发超声阵列探头的N个阵元发射超声波信号，每个阵元发射超声波信号时，全部阵元同时接收超声回波信号；将发射阵元i、接收阵元j采集的超声回波信号记为S_ij(i＝1,2,…,N；j＝1,2,…,N)；

(2)划分试块的检测区域，并将检测区域离散化；根据超声相控阵检测的检出率划分聚焦点类型，对于检出率高于检测标准的聚焦点，执行步骤3；对于检出率低于检测标准的聚焦点，执行步骤4；

(3)建立二维直角坐标系Oxz，坐标原点O设置在楔块下表面中心，利用延时规则将换能器中所有发射-接收阵元组合的超声回波信号在该聚焦点叠加，获得表征该聚焦点信号强度的幅值I(x,z)；分别获得被测区域内每一个目标点的幅值，完成整个被检测区域内的检测，I(x,z)的表达式如下：

其中t_ij(x,z)为得到该幅值的延迟时间，包括从阵元i传播到聚焦点P的时间，以及从聚焦点P传播到阵元j的时间，具体如下：

式中，(x_i,0)、(x_j,0)分别为发射阵元和接收阵元的坐标，c为超声波在试块中的传播速度；

(4)建立二维直角坐标系Oxz，坐标原点O设置在楔块下表面中心；根据扩散衰减系数重新定义幅值强度；利用延时规则将所有可能间接传播到聚焦点产生的幅值和直接传播到聚集点产生的幅值进行叠加，获得表征该聚焦点信号强度的幅值I(x,z)；分别获得被测区域内每一个目标点的幅值，完成整个被检测区域内的检测；

在直接传播过程中传播到聚焦点的扩散衰减系数k_d：

在间接传播过程中传播到聚集点的扩散衰减系数k_i：

其中，h为试块的检测厚度；

聚焦点(x,z)的幅值I(x,z)表示为：

其中，

本发明的有益效果是：本发明公开了一种基于全聚焦成像算法改进边缘缺陷检测的方法，现在的全聚焦成像算法仅仅考虑当阵列换能器的波之间发射到缺陷位置上，忽略了当阵列换能器的波传播到试块底部，反射到缺陷上产生的回波信号，本发明在考虑到多模型的全聚焦算法的基础上，提出了改进型的全聚焦成像算法，在考虑到波衰减的情况下，加进了底面反射回波对缺陷成像的影响，本发明方法能更加精确的仿真真实超声相控阵检测，提高边缘缺陷检测的精确度。

附图说明

图1为检测区域离散化结果示意图；

图2为检出率高于检测标准的聚焦点信号强度的幅值计算示意图；

图3为检出率低于检测标准的聚焦点信号强度的幅值计算示意图。

具体实施方式

下面结合附图和具体实施例对本发明作进一步详细说明。

本发明提出的一种基于全聚焦成像算法改进边缘缺陷检测的方法，该方法包括以下步骤：

(1)采集超声阵列全矩阵数据：依次激发超声阵列探头的N个阵元发射超声波信号，每个阵元发射超声波信号时，全部阵元同时接收超声回波信号；将发射阵元i、接收阵元j采集的超声回波信号记为S_ij(i＝1,2,…,N；j＝1,2,…,N)；S_ij是一组包含了每个时间采样点接收信号的幅值的数据。

(2)划分试块的检测区域，并将检测区域离散化，以规则的矩形试块为例，离散化结果如下图1所示；根据超声相控阵检测的检出率划分聚焦点类型，对于检出率高于检测标准的聚焦点(位置在没有用实线连接起来的中间区域)，执行步骤3进行更为精确的定位；对于检出率低于检测标准的聚焦点来说(除最外侧方框外，用方框框起来的位置)，受边界回波的影响，执行步骤4；检测标准可以选择90％。

(3)如图2所示，建立二维直角坐标系Oxz，坐标原点O设置在楔块下表面中心，利用延时规则将换能器中所有发射-接收阵元组合的超声回波信号在该聚焦点叠加，获得表征该聚焦点信号强度的幅值I(x,z)；分别获得被测区域内每一个目标点的幅值，完成整个被检测区域内的检测，I(x,z)的表达式如下：

其中t_ij(x,z)为得到该幅值的延迟时间，包括从阵元i传播到聚焦点P的时间，以及从聚焦点P传播到阵元j的时间，具体如下：

式中，(x_i,0)、(x_j,0)分别为发射阵元和接收阵元的坐标，c为超声波在试块中的传播速度；

(4)针对图1所示，用矩形框框起来的聚焦点来说，考虑到实际检测中常常会出现如图3所示的全聚焦模型，除了直接传播到P点的波，一部分波会直接传播到试块底部，再反射回来，所以对于传统的全聚焦来说，缺少了这一部分的值。但直接传播和间接传播，波衰减程度不一致，所以扩散衰减系数在区分这两部所得到的幅值，起到了重要作用。

如图3所示，建立二维直角坐标系Oxz，坐标原点O设置在楔块下表面中心；根据扩散衰减系数重新定义幅值强度；利用延时规则将所有可能间接传播到聚焦点产生的幅值和直接传播到聚集点产生的幅值进行叠加，获得表征该聚焦点信号强度的幅值I(x,z)；分别获得被测区域内每一个目标点的幅值，完成整个被检测区域内的检测。

在单一介质的二维模型中，使用距离平方根的反比来表示扩散衰减系数，故在直接传播过程中传播到聚焦点的扩散衰减系数k_d：

在间接传播过程中传播到聚集点的扩散衰减系数k_i：

其中，h为试块的检测厚度；

聚焦点(x,z)的幅值I(x,z)表示为：

其中，