CN105004792A - 一种用于微裂纹检测的非线性超声相控阵成像方法 - Google Patents

Abstract

一种用于微裂纹检测的非线性超声相控阵成像方法，本方法利用超声相控阵检测***，采集微缺陷在顺序聚焦与并行聚焦方式下的检测信号，首先确定合适的接收延时t_r与时域截取窗的宽度T，对截取的检测信号进行傅里叶变换；再分别计算顺序聚焦与并行聚焦方法的声动能E_s、E_p及两种聚焦方式的相对能量差γ；最后根据γ确定非线性源的位置，即微缺陷的位置。

Description

一种用于微裂纹检测的非线性超声相控阵成像方法

技术领域

本发明属无损检测领域，具体涉及一种非线性超声相控阵成像方法，该方法用于非线性源(如疲劳裂纹)的检测与识别。

背景技术

结构在长期服役过程中受到疲劳载荷、冲击及腐蚀、高温等环境因素的作用会逐渐老化，表现为出现不同程度的损伤。其中，疲劳裂纹是一种不稳定性缺陷，对结构的潜在危害性最大。在外界因素的作用下，疲劳裂纹可能发生迅速扩展使结构发生突然断裂。因此，采取有效的手段实现结构疲劳裂纹的早期检测、定位，对于确保重大基础设施的正常运行和人们生命财产安全至关重要。

非线性超声技术对材料微观组织结构、机械损伤以及组织病理的敏感性高，特别适合于结构中疲劳裂纹的检测。但目前的非线性超声测量结果或局限于空间平均值，或通过繁琐的成像实现隐式的局部化，缺乏有效的非线性超声成像手段。同时，这些隐式的非线性成像方法需要利用光学测量方法获得检测点信息，其检测范围因而仅局限于表面测量，使这些技术的应用受到了很大限制。

与之相反，线性超声成像因相控阵技术而变得极为容易。全矩阵成像是近年来广泛应用的一种阵列成像方法。与相控阵成像中阵元的并行激励不同，在全矩阵成像中，利用全矩阵数据采集模式，对阵元进行顺序激励，并在后处理中对全矩阵数据施加虚拟时延，以等价地仿效施加物理时延下的并行激励。将顺序激励并通过后处理实现的虚拟聚焦和阵元并行激励实现的真实聚焦分别称为顺序聚焦和并行聚焦。利用上述两种相控阵激励方式，将非线性与相控阵技术相结合，考虑两种激励模式下聚焦点处的非线性响应的差别表征聚焦点处是否存在疲劳裂纹，从而实现疲劳裂纹的空间定位。

本文创新点在于提出了一种基于并行激励、顺序激励两种相控阵激励方式下扩散声场能量差的非线性阵列成像方法，综合考虑超声波与微缺陷相互作用的多种非线性效应，能够较精确地实现微缺陷定位检测。

发明内容

本发明的目的在于发展一种基于两种相控阵聚焦方式下扩散声场能量差的非线性阵列成像方法，通过该方法可以较精确地确定微缺陷的位置。基于不同相控阵聚焦方式下扩散声场能量差的非线性成像方法利用顺序聚焦与并行聚焦方式中阵列成像区域中聚焦点的相对能量差表示该点的非线性效应，存在非线性源(微缺陷)时聚焦点的相对能量差较大，因而可以判断非线性源的有无。

本发明提出的基于两种相控阵聚焦方式下扩散声场能量差的非线性成像方法，其基本原理如下：

两种聚焦检测模式中超声波的线性传播特性相同，但其非线性传播特性不同。并行聚焦时，结构中聚焦位置处的绝对声压要高于顺序聚焦时任意单独激励下的绝对声压。因此，在并行聚焦检测时，若聚焦点处存在损伤等非线性源，其非线性效应将会导致更多的能量从基波频率向其它频率转移。因此，顺序聚焦和并行聚焦两种检测成像方法在聚焦点处激励的线性响应及成像结果完全相同，但两种方法在聚焦点处激发的非线性响应及激励带宽内(基波)的能量损失有很大的差别。

一般说来，介质的非线性响应不仅包含高次谐波，也包含分谐波以及激励带宽内频率分量的和频、差频的组合，而这些非线性效应的共性是有能量从激励带宽的转移。因此，利用激励带宽内(基波)的能量损失作为非线性表征指标，考虑到了疲劳裂纹可能产生的多种非线性效应，可望大大提高对疲劳裂纹的检测能力。

基于相控阵聚焦原理的微裂纹定位方法的关键是需要获得两种聚焦检测模式中聚焦点处激励带宽内(基波)的能量损失。由于超声测量***无法直接获得聚焦点处的超声场，因此，需要对聚焦声场的传播及分布特性进行研究，以找到一种可用于聚焦点处基波声能损失计算的有效途径。

首先，分析两种检测模式下声能损失差异的时间、空间分布特点。在两种检测模式下，波的线性传播特性对声能的干涉作用是相同的，且仪器非线性和阵列与试件界面接触非线性对声能传输的影响也相同。同时，在非聚焦处，由于各阵元发射的超声波没有发生相互干涉作用，其非线性相互作用造成的激励带宽内(基波)的能量损失也是相同的。因此，激励频带内(基波)的能量差别仅仅是由于聚焦点处聚焦时刻的幅值差异造成的。因此，可以利用聚焦时刻以后任意瞬时两种聚焦检测方法得到的***总声能的差值来表征聚焦点处的相对非线性能量损失。为此，可以通过对后续的扩散声场分析来实现。

在理想扩散声场中，任意位置的声能密度均匀，且声能在各个方向传递的几率相等。实际上，由于超声波在晶界及试件边界的多次散射，初始相干激励在传播一段时间后可满足扩散声场的存在条件。在扩散声场中，任意点的声能与此时刻***的总能量成正比。因此，尽管聚焦点处的相对非线性能量损失仅出现在相干场传播的子集内，但它在扩散声场中是均匀分布的。因此，在并行聚焦和顺序聚焦的扩散声场中，任意位置处的能量差异为聚焦点处聚焦时刻的能量损失测量提供了一个近似值，相应地可用于该点处超声非线性的表征。

随着时间的增加，声场将收敛于扩散状态。但是，由于超声波传播的耗散性，信噪比会随着时间的增加而降低。同时，由于顺序激励时检测信号的幅值较低，其受非相干噪声的影响更大，导致两种聚焦方式下扩散声场能量比较的有效性随数据采集时延的增加而降低。因此，在使接收信号幅值最大和扩散声场收敛的情况下，应折中地选择接收信号的时延。由于扩散声场仅是统计意义上的静态，因此，可以利用一定时间窗口内(t_r,t_r+T)的积分来表示其能量值，其中t_r为信号的接收延时时间，T为时间窗宽。对于N阵元阵列，顺序聚焦检测中聚焦点处的扩散声动能E_S可表示为：

$E_S\left(\stackrel{\‾}{r}\right)=\underset{l=1}{\overset{N}{\Σ}}\left({\∫}_{{\mathrm{\ω}}_1}^{{\mathrm{\ω}}_2}{\mathrm{\ω}}^2{\left|\underset{k=1}{\overset{N}{\Σ}}{F}_{kl}\left(\mathrm{\ω}\right)e^{{\mathrm{i\ω\δ}}_k\left(\stackrel{\‾}{r}\right)}\right|}^{2}d\mathrm{\ω}\right)---\left(1\right)$

类似地，并行激励下声动能可以表示为：

$E_P\left(\stackrel{\‾}{r}\right)=\underset{l=1}{\overset{N}{\Σ}}\left({\∫}_{{\mathrm{\ω}}_1}^{{\mathrm{\ω}}_2}{\mathrm{\ω}}^2{\left|H_l(\stackrel{\‾}{r},\mathrm{\ω})\right|}^{2}d\mathrm{\ω}\right)---\left(2\right)$

利用顺序激励和并行激励下归一化的能量差，定义如下的非线性特征指标：

$\mathrm{\γ}\left(\stackrel{\‾}{r}\right)=\frac{E_S\left(\stackrel{\‾}{r}\right)-E_P\left(\stackrel{\‾}{r}\right)}{E_S\left(\stackrel{\‾}{r}\right)}---\left(3\right)$

其中，为非线性效应导致的激励带宽内(基波)的相对能量损失。f_kl(t)为k阵元顺序激励l阵元接收的时域信号。是聚焦在处时激励阵元k上施加的时延。为所有阵元在时延下并行激励，阵元l接收到的时域信号。F_kl(ω)和H_l(ω)为时间窗口(t_r,t_r+T)内，f_kl(t)和的频谱，和 $H_l\left(\mathrm{\ω}\right)={\∫}_{t_r}^{t_r+T}{h}_l\left(t\right)e^{-\mathrm{i\ωt}}\mathrm{dt},$ t_r为接收时延，(ω₁，ω₂)为基波频带范围。激励带宽对疲劳裂纹检测的敏感性以及聚焦点处激励能量的强度及转移量都有不同程度的影响，需要综合考虑。激励带宽选择的一般原则是无内部能量转移情况下的最大值。

因此，基于并行聚焦和顺序聚焦的扩散声场分析，其一定时延下扩散场的声动能差可用于聚焦点处声能损失的表征。

本发明的技术方案如下：

本发明中使用了超声相控阵检测装置，如图1所示。该超声相控阵检测装置包括计算机1、超声相控阵***2、相控阵换能器3以及待检测试件4。其中，超声相控阵***2一端与计算机1连接，另一端与相控阵换能器3连接，相控阵换能器3与被测试件4通过耦合介质进行耦合。在计算机的控制下超声相控阵***激励/接收模块产生激励信号，通过相控阵换能器激励出超声波信号，沿待测试件传播，并通过相控阵传感器接收反射的超声波信号，然后通过相控阵***信号激励/接收模块进行接收传输到计算机中，通过计算机中相应采集软件即可获得检测的A扫波形。所述的相控阵换能器采用由32个阵元组成的线阵相控阵探头。

本发明提出的基于不同相控阵聚焦方式下扩散声场能量差的非线性成像方法，其流程图如图2所示。具体可以按照以下步骤实施检测：

步骤一：利用超声相控阵检测***采集顺序聚焦与并行聚焦方式下的检测信号。

1)将如图1所示的超声相控阵检测***设置为全矩阵采集模式，采集被测缺陷的全矩阵数据f_kl(t)(k＝1,2,3…N,l＝1,2,3…N)。其中，下标k表示阵列换能器中第k个阵元激励，l表示阵列换能器中第l个阵元接收；

2)将超声相控阵检测***设置为点聚焦模式，采集并行聚焦方式下每个聚焦点的检测信号h_l(t)(l＝1,2,3…N)。其中，下标l表示阵列换能器中第l个阵元接收；

步骤二：建立如图3所示的成像坐标系，定义相控阵探头中每个阵元在坐标系中的位置以及每个成像像素点在坐标系中的位置。坐标系原点定义在阵元序列的中心，其中(x_tx,z_tx)为激励阵元的坐标，(x_rx,z_rx)为接收阵元的坐标，(x,z)为成像点坐标。

步骤三：确定接收延时t_r与时域截取窗的宽度T，对截取的检测信号进行傅里叶变换：

$F_{kl}\left(\mathrm{\ω}\right)={\∫}_{t_r}^{t_r+T}{f}_{kl}\left(t\right)e^{-i\mathrm{\ω}t}dt---\left(4\right)$

$H_l\left(\mathrm{\ω}\right)={\∫}_{t_r}^{t_r+T}{h}_l\left(t\right)e^{-i\mathrm{\ω}t}dt---\left(5\right)$

步骤四：根据公式(1)、(2)和(3)，分别计算聚焦点处两种聚焦方法的声动能E_s、E_p及相对能量差γ；

步骤五：对两种聚焦方式下的的能量E_s、E_p进行成像，对两种聚焦方式的相对能量差γ进行成像，以定位非线性源。

附图说明

图1检测装置框图。

图2非线性超声阵列成像方法步骤。

图3被测试件。

图4成像坐标系。

具体实施方式

以下利用检测实验得到的数据，对本发明方法的有效性进行详细说明。

采用的试验***如图1所示，实验中采用的线性阵列相控阵探头，型号为5L-32A5，其参数如表1所示。被测试件如图4所示，材料为钢，尺寸为240mm×47mm×25mm，在其中部有一长20mm线切割缺口，利用疲劳试验机在线切割顶端处加工一长约9mm的疲劳裂纹，因此裂纹尖端距离检测表面20mm。激励信号中心频率为5MHz，采样频率为16.6MHz。

表1探头基本参数

步骤一：利用超声相控阵检测***采集顺序聚焦与并行聚焦方式下的检测信号；

(1)将如图1所示的超声相控阵检测***设置为全矩阵采集模式，采集被测缺陷的全矩阵数据f_kl(t)(k＝1,2,3…N,l＝1,2,3…N)，其中，其中，下标k表示阵列换能器中第k个阵元激励，l表示阵列换能器中第l个阵元接收；

(2)将超声相控阵检测***设置为单点聚焦模式，采集并行聚焦方式下每个聚焦点的检测信号h_l(t)(l＝1,2,3…N)。其中，下标l表示阵列换能器中第l个阵元接收；

步骤二：建立如图3所示的成像坐标系，定义相控阵探头中每个阵元在坐标系中的位置以及每个成像像素点在坐标系中的位置。坐标系原点定义在阵元序列的中心，其中(x_tx,z_tx)为激励阵元的坐标，(x_rx,z_rx)为接收阵元的坐标，(x,z)为成像点坐标。

步骤三：确定接收延时t_r＝0.9ms与时域截取窗的宽度T＝200μs，对截取的检测信号进行傅里叶变换；

步骤四：分别计算聚焦点处两种聚焦方法的声动能E_s、E_p及相对能量差γ；

步骤五：对两种聚焦方式下的的能量E_s、E_p进行成像，对两种聚焦方式的相对能量差γ进行成像，以定位非线性源。由图可知，γ值在疲劳裂纹尖端处附近较大，可以确定疲劳裂纹的位置。

Claims (3)

1.超声相控阵检测***，其特征在于：该***包括计算机(1)、超声相控阵***(2)、相控阵换能器(3)以及待检测试件(4)；其中，超声相控阵***(2)一端与计算机连接，另一端与相控阵换能器连接，相控阵换能器(3)与待检测试件(4)通过耦合介质进行耦合；在计算机的控制下超声相控阵***激励/接收模块产生激励信号，通过相控阵换能器激励出超声波信号，沿待测试件传播，并通过相控阵传感器接收反射的超声波信号，然后通过相控阵***信号激励/接收模块进行接收传输到计算机中，通过计算机中相应采集软件即可获得检测的A扫波形；所述的相控阵换能器采用由32个阵元组成的线阵相控阵探头。

2.依权利要求1所述的超声相控阵检测***，一种非线性超声相控阵成像方法，其特征在于：该方法具体识别步骤如下，

步骤一：利用超声相控阵检测***采集顺序聚焦与并行聚焦方式下的检测信号；

(1)将超声相控阵检测***设置为全矩阵采集模式，采集被测缺陷的全矩阵数据f_kl(t)(k＝1,2,3…N,l＝1,2,3…N)；其中，下标k表示阵列换能器中第k个阵元激励，l表示阵列换能器中第l个阵元接收；

(2)将超声相控阵检测***设置为点聚焦模式，采集并行聚焦方式下每个聚焦点的检测信号h_l(t)(l＝1,2,3…N)；其中，下标l表示阵列换能器中第l个阵元接收；

步骤二：建立成像坐标系，定义相控阵探头中每个阵元在坐标系中的位置以及每个成像像素点在坐标系中的位置；坐标系原点定义在阵元序列的中心，其中(x_tx,z_tx)为激励阵元的坐标，(x_rx,z_rx)为接收阵元的坐标，(x,z)为成像点坐标；

步骤三：确定接收延时t_r与时域截取窗的宽度T，对截取的检测信号进行傅里叶变换：

$F_{kl}\left(\mathrm{\ω}\right)={\∫}_{t_r}^{t_r+T}{f}_{kl}\left(t\right)e^{-i\mathrm{\ω}t}dt---\left(1\right)$

$H_l\left(\mathrm{\ω}\right)={\∫}_{t_r}^{t_r+T}{h}_l\left(t\right)e^{-i\mathrm{\ω}t}dt---\left(2\right)$

步骤四：分别计算聚焦点处两种聚焦方法的声动能E_s、E_p及相对能量差γ；

步骤五：对两种聚焦方式下的的能量E_s、E_p进行成像，对两种聚焦方式的相对能量差γ进行成像，以定位非线性源；

该方法利用聚焦时刻以后任意瞬时两种聚焦检测方法得到的***总声能的差值来表征聚焦点处的相对非线性能量损失；但是测量***无法直接获得聚焦点处的超声场，因此聚焦点的能量损失必须根据其后续特性来推断出；

在理想扩散声场中，任意位置的声能密度均匀，且声能在各个方向传递的几率相等；实际上，由于超声波在晶界及试件边界的多次散射，初始相干激励在传播一段时间后可满足扩散声场的存在条件；在扩散声场中，任意点的声能与此时刻***的总能量成正比；因此，尽管聚焦点处的相对非线性能量损失仅出现在相干场传播的子集内，但它在扩散声场中是均匀分布的；因此，在并行聚焦和顺序聚焦的扩散声场中，任意位置处的能量差异为聚焦点处聚焦时刻的能量损失测量提供了一个近似值，相应地可用于该点处超声非线性的表征；因此，基于并行聚焦和顺序聚焦的扩散声场分析，其一定时延下扩散场的声动能差可用于聚焦点处声能损失的表征；

由于扩散声场仅是统计意义上的静态，因此，利用一定时间窗口内(t_r,t_r+T)的积分来表示其能量值，其中t_r为信号的接收延时时间，T为时间窗宽；随着时间的增加，声场将收敛于扩散状态；但是，由于超声波传播的耗散性，信噪比会随着时间的增加而降低；同时，由于顺序激励时检测信号的幅值较低，其受非相干噪声的影响更大，导致两种聚焦方式下扩散声场能量比较的有效性随数据采集时延的增加而降低；因此，在使接收信号幅值最大和扩散声场收敛的情况下，应折中地选择接收信号的时延；

对于N个阵元组成的阵列，顺序聚焦检测中聚焦点处的扩散声动能可表示为：

$E_S\left(\stackrel{\‾}{r}\right)=\underset{l=1}{\overset{N}{\Σ}}\left({\∫}_{{\mathrm{\ω}}_1}^{{\mathrm{\ω}}_2}{\mathrm{\ω}}^2{\left|\underset{k=1}{\overset{N}{\Σ}}{F}_{kl}\left(\mathrm{\ω}\right)e^{{\mathrm{i\ω\δ}}_k\left(\stackrel{\‾}{r}\right)}\right|}^{2}d\mathrm{\ω}\right)---\left(3\right)$

其中(ω₁，ω₂)为基波频带范围，是聚焦在处时激励阵元k上施加的时延；

类似地，并行激励下声动能可以表示为：

$E_P\left(\stackrel{\‾}{r}\right)=\underset{l=1}{\overset{N}{\Σ}}\left({\∫}_{{\mathrm{\ω}}_1}^{{\mathrm{\ω}}_2}{\mathrm{\ω}}^2{\left|H_l(\stackrel{\‾}{r},\mathrm{\ω})\right|}^{2}d\mathrm{\ω}\right)---\left(4\right)$

利用顺序激励和并行激励下归一化的能量差，非线性效应导致的激励带宽内(基波)的相对能量损失可表示为：

$\mathrm{\γ}\left(\stackrel{\‾}{r}\right)=\frac{E_S\left(\stackrel{\‾}{r}\right)-E_P\left(\stackrel{\‾}{r}\right)}{E_S\left(\stackrel{\‾}{r}\right)}---\left(5\right)$

根据非线性超声相控阵成像结果可以确定非线性源的位置。

3.根据权利要求2所述的一种非线性超声相控阵成像方法，其特征在于：实验中采用的线性阵列相控阵探头，型号为5L-32A5，其参数如表1所示；被测试件的材料为钢，尺寸为240mm×47mm×25mm，在其中部有一长20mm线切割缺口，利用疲劳试验机在线切割顶端处加工一长约9mm的疲劳裂纹，因此裂纹尖端距离检测表面20mm；激励信号中心频率为5MHz，采样频率为16.6MHz；

表1  探头基本参数

步骤一：利用超声相控阵检测***采集顺序聚焦与并行聚焦方式下的检测信号；

(1)将超声相控阵检测***设置为全矩阵采集模式，采集被测缺陷的全矩阵数据f_kl(t)(k＝1,2,3…N,l＝1,2,3…N)，其中，下标k表示阵列换能器中第k个阵元激励，l表示阵列换能器中第l个阵元接收；

(2)将超声相控阵检测***设置为单点聚焦模式，采集并行聚焦方式下每个聚焦点的检测信号h_l(t)(l＝1,2,3…N)；其中，下标l表示阵列换能器中第l个阵元接收；

步骤二：建立成像坐标系，定义相控阵探头中每个阵元在坐标系中的位置以及每个成像像素点在坐标系中的位置；坐标系原点定义在阵元序列的中心，其中(x_tx,z_tx)为激励阵元的坐标，(x_rx,z_rx)为接收阵元的坐标，(x,z)为成像点坐标；

步骤三：确定接收延时t_r＝0.9ms与时域截取窗的宽度T＝200μs，对截取的检测信号进行傅里叶变换；

步骤四：分别计算聚焦点处两种聚焦方法的声动能E_s、E_p及相对能量差γ；

步骤五：对两种聚焦方式下的的能量E_s、E_p进行成像，对两种聚焦方式的相对能量差γ进行成像，以定位非线性源；由图可知，γ值在疲劳裂纹尖端处附近较大，可以确定疲劳裂纹的位置。