本发明涉及一种分步筛选有效信号的曲面结构缺陷全聚焦成像方法，其属于无损检测技术领域。
背景技术：
全聚焦方法(totalfocusingmethod，tfm)是一种利用全矩阵数据对待测区域逐点聚焦，实现缺陷成像的无损检测信号后处理方法，具有成像信噪比高，缺陷定位、定量准确等优点。其中，利用tfm对曲面结构工件内部缺陷实施检测时，一般采用水浸耦合方式，首先对曲面轮廓进行tfm成像，求解曲面轮廓位置信息，之后再对内部缺陷实施tfm成像，给出缺陷位置和尺寸信息。两次成像过程中，超声波在曲面和缺陷处散射现象复杂，全矩阵数据中存在大量冗余信号，导致计算量庞大，计算效率较低。因此，tfm在曲面结构缺陷成像检测应用方面受到限制。
为解决以上问题，研究人员根据信号幅值对进行筛选，通过减少两次tfm处理的信号数量来提高效率。考虑到发射阵元和接收阵元间距增大时，接收信号幅值降低，因此可以根据各阵元接收到信号的最高幅值设置一定阈值，从全矩阵数据中筛选对成像有显著作用的信号。阈值设置越高，计算效率越高，但可能会减少有用信号，降低曲面轮廓提取精度和缺陷成像质量。如何兼顾曲面结构的tfm计算效率和成像质量是一个亟待解决的问题。本发明提出了一种分步筛选有效信号的曲面结构缺陷全聚焦成像方法，可在减少冗余信号，提高计算效率的同时，保证成像质量。
技术实现要素：
本发明提供一种分步筛选有效信号的曲面结构缺陷全聚焦成像方法。其目的是针对曲面结构缺陷全聚焦成像存在冗余数据，成像质量与计算效率不能兼顾的问题，利用相控阵超声检测仪采集全矩阵数据，在对曲面轮廓和内部缺陷成像时分别设置合适阈值，以筛选有效信号。利用筛选后的有效信号，基于费马原理和全聚焦方法先后对曲面轮廓和内部缺陷成像，在提高计算效率的同时保证成像质量，实现曲面轮廓和缺陷特征信息提取。
本发明采用的技术方案是：利用相控阵超声检测仪和线性阵列相控阵探头，采取水浸耦合方式采集全矩阵数据；设置合理阈值进行第一次有效信号筛选，利用全聚焦方法实现曲面轮廓成像，并获取曲面轮廓位置信息；提高阈值范围，对第一次筛选出的有效信号进行第二次筛选，根据曲面轮廓信息和费马原理计算声束在水和待测工件中的传播路径，并利用全聚焦方法进行内部缺陷成像，实现缺陷定量；所述方法采用如下步骤：
（a）相控阵检测参数确定
根据待测工件的材料、几何尺寸及检测范围确定相控阵超声检测参数，包括相控阵探头频率、相控阵探头阵元数、水浸耦合距离、采样频率；
（b）全矩阵数据采集
采用步骤（a）中确定的相控阵检测参数，控制相控阵探头处于待测工件正上方，共获得n2个a扫描信号构成的数据矩阵，其中n为阵元数；
（c）坐标系建立
以相控阵探头主动轴方向为x轴正向，工件深度方向为y轴正向建立坐标系，并设探头首阵元位置为坐标原点；将被检区域划分为矩形网格，各网格节点作为图像重建点；任意图像重建点p的横、纵坐标分别定义为x和y；
（d）第一次有效信号筛选及曲面轮廓成像
在曲面轮廓成像前对全矩阵数据进行第一次有效信号筛选，考虑到曲面轮廓对声波散射作用，在大曲率位置处的散射波幅值较低，故可设置阈值，滤除这些幅值较低的信号；结合常规缺陷定量中采用的-12db法，即取每个阵元接收到的所有信号中最大幅值1/4为阈值，如式（1）所示
式中，smn为阵元m发射，并由阵元n接收的信号（1≤m，n≤n），thresholdm为发射阵元m对应阈值，只有高于此阈值的信号为有效信号，被用来计算轮廓成像；第一次筛选后，保留下来的由阵元m发射的有效信号数量为nm，各有效信号记为smn’（1≤n’≤nm），n’为有效信号接收阵元的重排序号；
利用第一次筛选后的有效信号进行延时叠加，获得p点信号幅值ap
式中，tmpn’为声波从阵元m发射，经过p点并被阵元n’接收时的传播时间，仅与耦合介质有关，cmn’为总校准系数，如式（3）所示
式中，dm、dn’分别为发射阵元m与接收阵元n’对应的指向性校准系数，bm、bn’分别为发射阵元m与接收阵元n’对应的扩散衰减系数；
利用式（2）对曲面轮廓待成像区域内各网格点进行逐点成像，得到第一次有效信号筛选后的曲面轮廓全聚焦重建图像，并拟合获得曲面轮廓位置信息；
（e）第二次有效信号筛选及内部缺陷成像
曲面结构对声波的散射作用会造成缺陷成像时产生伪像，需筛选去除干扰信号，因此要进一步提高筛选阈值，对经过第一次筛选后的有效信号进行第二次筛选；结合常规缺陷定量过程中采用-6db法，即取每个发射阵元对应所有接收信号中最大幅值的1/2作为新的阈值thresholdm’，如式（4）所示
只有高于该阈值的信号为有效信号，被用于内部缺陷成像；第二次筛选后，保留下来的由阵元m发射的有效信号数量为nm’，各有效信号记为smn’’（1≤n’’≤nm’），n’’为有效信号接收阵元的重排序号；根据步骤（d）得到的曲面轮廓位置信息，根据费马原理计算声束传播路径；
利用第二次筛选后的有效信号进行延时叠加，获得p点信号幅值ap’
式中，tmpn’’为声波从阵元m发射，经过p点并被阵元n’’接收时的传播时间，与耦合介质和待测工件有关，cmn’’为总校准系数，如式（6）所示
式中，dn’’和bn’’分别为接收阵元n’’对应的指向性校准系数和扩散衰减系数；
利用式（5）对内部缺陷待成像区域内各网格点进行逐点成像，得到第二次有效信号筛选后的内部缺陷全聚焦重建图像；
（f）缺陷定量检测
根据步骤（e）得到的曲面结构内部缺陷重建图像，读取成像区域内的峰值坐标点，利用-6db法对缺陷进行定量检测。
本发明的有益效果是：这种分步筛选有效信号的曲面结构缺陷全聚焦成像方法，在曲面轮廓重建和内部缺陷成像前分别对全矩阵数据进行筛选，根据不同成像对象调整阈值，减少冗余信号，提高计算效率；同时避免过度筛选有效信号，丢失特征信息，保证了成像质量，为解决曲面结构内部缺陷全聚焦成像质量与计算效率不能兼顾的问题提供新的思路。
附图说明
下面结合附图和实例对本发明做进一步说明。
图1是采用的相控阵超声检测系统示意图。
图2是含有横通孔的曲面结构碳钢工件示意图。
图3是第一次有效信号筛选结果。
图4是经过第一次有效信号筛选后的曲面轮廓全聚焦重建图像。
图5是第二次有效信号筛选结果。
图6是经过第二次有效信号筛选后的内部缺陷全聚焦重建图像。
图7是只进行一次-12db法有效信号筛选后的内部缺陷全聚焦重建图像。
图8是只进行一次-6db法有效信号筛选后的内部缺陷全聚焦重建图像。
具体实施方式
分步筛选有效信号的曲面结构缺陷全聚焦成像方法，采用的相控阵超声检测系统如图1所示，其中包括相控阵超声检测仪及配套笔记本电脑、64阵元线性阵列相控阵探头和水槽。具体检测及处理步骤如下：
（a）试验对象如图2所示碳钢工件，其高度为40.0mm，曲面结构部分高2.4mm，对应半径r为6.4mm，材料纵波声速为5890m/s，横波声速为3230m/s，试块内部圆弧正下方距试块底部20.0mm处加工一个φ2.0mm横通孔。
（b）采用标称频率5mhz的64阵元相控阵探头实施全矩阵数据采集，采样频率100mhz、探头距曲面结构圆弧最高处15.0mm。
（c）根据步骤（b）中确定的检测参数，采用水浸耦合方式采集获得64×64个a扫描信号构成的数据集，并以.txt文件形式导出。考虑曲面轮廓成像对数据的完整性要求较高，对全矩阵数据采用-12db法进行第一次有效信号筛选，筛选结果如图3所示，其中64×64个网格点分别代表smn，白色部分代表经第一次筛选后的有效信号smn’。计算可知，筛选后信号数量由4096减少为2247。利用第一次筛选后的信号进行全聚焦成像结果如图4所示，拟合得到曲面轮廓位置信息。
（d）对经过第一次筛选的有效信号进行第二次有效信号筛选。考虑到内部缺陷成像过程中存在伪像干扰，故设置较高阈值去除更多冗余信号。采用-6db法筛选方法进行第二次有效信号筛选，筛选结果如图5所示，白色部分代表经第二次筛选后的有效信号smn’’。计算可知，信号数量由2247减少为1508，需要处理的信号数量进一步减少。基于步骤（c）中得到的曲面轮廓信息，根据费马原理计算声束入射点与重建区域的最短声束路径，重建得到图6所示的内部缺陷图像。与图7给出的只进行一次-12db法有效信号筛选得到的内部缺陷成像结果相比，图6所用信号数量减少32%。与图8给出的只进行一次-6db法有效信号筛选得到的内部缺陷成像结果相比，图6中缺陷信号与伪缺陷信号分离度更好，缺陷信号幅值提高2.90db，图像信噪比提高1.28db。
（e）根据步骤（d）得到的内部缺陷重建图像，读取成像区域内的峰值坐标点，利用-6db法对缺陷进行定量检测。测量得到缺陷深度为距圆弧表面最高处22.2mm，定量误差0.2mm。
综上可见，分步筛选有效信号的方法针对曲面结构内部缺陷成像检测，能够根据不同测量对象设置不同阈值，大幅减少冗余信号，提高计算效率；同时能够保证成像质量，避免缺陷主要特征信息丢失，为解决曲面结构内部缺陷全聚焦成像质量与计算效率不能兼顾的问题提供新的思路。
1.一种分步筛选有效信号的曲面结构缺陷全聚焦成像方法，其特征在于：该方法利用相控阵超声检测仪和线性阵列相控阵探头，采取水浸耦合方式采集全矩阵数据；设置合理阈值进行第一次有效信号筛选，利用全聚焦方法实现曲面轮廓成像，并获取曲面轮廓位置信息；提高阈值范围对第一次筛选出的有效信号进行第二次筛选，根据曲面轮廓信息和费马原理计算声束在水和待测工件中的传播路径，并利用全聚焦方法进行内部缺陷成像，实现缺陷定量；所述方法采用如下步骤：
（a）相控阵检测参数确定
根据待测工件的材料、几何尺寸及检测范围确定相控阵超声检测参数，包括相控阵探头频率、相控阵探头阵元数、水浸耦合距离、采样频率；
（b）全矩阵数据采集
采用步骤（a）中确定的相控阵检测参数，控制相控阵探头处于待测工件正上方，共获得n2个a扫描信号构成的数据矩阵，其中n为阵元数；
（c）坐标系建立
以相控阵探头主动轴方向为x轴正向，工件深度方向为y轴正向建立坐标系，并设探头首阵元位置为坐标原点；将被检区域划分为矩形网格，各网格节点作为图像重建点；任意图像重建点p的横、纵坐标分别定义为x和y；
（d）第一次有效信号筛选及曲面轮廓成像
在曲面轮廓成像前对全矩阵数据进行第一次有效信号筛选，考虑到曲面轮廓对声波散射作用，在大曲率位置处的散射波幅值较低，故可设置阈值，滤除这些幅值较低的信号；结合常规缺陷定量中采用的-12db法，即取每个阵元接收到的所有信号中最大幅值1/4为阈值，如式（1）所示
式中，smn为阵元m发射，并由阵元n接收的信号（1≤m，n≤n），thresholdm为发射阵元m对应阈值，只有高于此阈值的信号为有效信号，被用来计算轮廓成像；第一次筛选后，保留下来的由阵元m发射的有效信号数量为nm，各有效信号记为smn’（1≤n’≤nm），n’为有效信号接收阵元的重排序号；
利用第一次筛选后的有效信号进行延时叠加，获得p点信号幅值ap
式中，tmpn’为声波从阵元m发射，经过p点并被阵元n’接收时的传播时间，仅与耦合介质有关，cmn’为总校准系数，如式（3）所示
式中，dm、dn’分别为发射阵元m与接收阵元n’对应的指向性校准系数，bm、bn’分别为发射阵元m与接收阵元n’对应的扩散衰减系数；
利用式（2）对曲面轮廓待成像区域内各网格点进行逐点成像，得到第一次有效信号筛选后的曲面轮廓全聚焦重建图像，并拟合获得曲面轮廓位置信息；
（e）第二次有效信号筛选及内部缺陷成像
曲面结构对声波的散射作用会造成缺陷成像时产生伪像，需筛选去除干扰信号，因此要进一步提高筛选阈值，对经过第一次筛选后的有效信号进行第二次筛选；结合常规缺陷定量过程中采用-6db法，即取每个发射阵元对应所有接收信号中最大幅值的1/2作为新的阈值thresholdm’，如式（4）所示
只有高于该阈值的信号为有效信号，被用于内部缺陷成像；第二次筛选后，保留下来的由阵元m发射的有效信号数量为nm’，各有效信号记为smn’’（1≤n’’≤nm’），n’’为有效信号接收阵元的重排序号；根据步骤（d）得到的曲面轮廓位置信息，根据费马原理计算声束传播路径；
利用第二次筛选后的有效信号进行延时叠加，获得p点信号幅值ap’
式中，tmpn’’为声波从阵元m发射，经过p点并被阵元n’’接收时的传播时间，与耦合介质和待测工件有关，cmn’’为总校准系数，如式（6）所示
式中，dn’’和bn’’分别为接收阵元n’’对应的指向性校准系数和扩散衰减系数；
利用式（5）对内部缺陷待成像区域内各网格点进行逐点成像，得到第二次有效信号筛选后的内部缺陷全聚焦重建图像；
（f）缺陷定量检测
根据步骤（e）得到的曲面结构内部缺陷重建图像，读取成像区域内的峰值坐标点，利用-6db法对缺陷进行定量检测。
技术总结