CN114755306A - 超声导波质量检测方法、系统及其应用 - Google Patents

Abstract

本发明公开了一种超声导波质量检测方法、系统及其应用，其中所述质量检测方法包括：采集检测中返回的超声波信号；通过改进的自适应噪声完备集合经验模态分解对该信号进行处理，获得信号分解后的固有模态函数，进而确定检测目标的质量情况。本发明的检测方法高效、精确，适于广泛应用。

Description

超声导波质量检测方法、系统及其应用

技术领域

本发明涉及超声导波用于质量检测的技术领域。

背景技术

锚杆是边坡、岩土深基坑等地表工程及隧道、采场等地下工程施工中广泛采用的一种加固支护工具，锚固是锚杆中伸入支座或基础中，用以进行受力传递的部分，一般情况下，锚杆中锚固段越长，其加固支护效果越佳。

在工程应用中，锚杆质量的好坏直接关系着项目的安全，一旦出现锚杆失效现象，会对岩土结构的稳定性产生巨大威胁，随之可能产生严重的现场事故和项目的巨大损失。引起锚杆失效的原因通常有多种，如：锚杆锚固段长度不足；锚固剂与杆体和围岩未能有效粘结；因地质条件变化或者爆破、重型机械和地震力引发的冲击使得锚杆被拉断；由于围岩的横向滑动，使锚杆的应力状态由轴向拉应力变为横向剪应力，导致锚杆被剪断；由于工作环境较恶劣，如水、潮湿空气、雨水等液体渗入到锚杆，致使锚杆发生腐蚀而引起应力集中，进而使得锚杆自应力集中处断裂等等。这些原因除锚杆本身的质量缺陷及强大外力破坏外，较多是在长期工作中自小的损伤产生的，需要通过对锚杆的及时有效的无损检测进行识别、预防和危险排除。

现有技术中开发了多种针对这类加固和/或支护工具的无损检测方法，如声频应力波法、光纤法、电磁波法和超声导波法等。其中，声频应力波法多采用小锤敲击震动从而产生应力波的形式，其检测结果容易因敲击力的不同而产生偏差；电磁波法利用了金属钢筋导电的特性，通过导体在阻抗变化界面会发生电势的突变的现象达到检测目的，但大多数的电磁波法，如电势差法、电驻波法自感应法等，只能检测锚杆的长度，对于锚杆密实度的检测能力很弱；光纤法大多用于锚固力的检测，其对焊接加工要求较高，需要预先埋置、不易操作，且检测使用的传感器、测力计等易受到环境的干扰。超声导波法相对于其他方法，对锚固质量更加敏感且仪器操作简单，但其信号处理难度相对于其他方法更高。如何在超声导波信号中快速、准确地获得被检测对象的质量信息是本领域一直期望解决的问题。

在信号的处理中，现有技术多依赖于傅里叶变换或快速傅里叶变换，部分技术手段使用了小波变换。其中，傅立叶变换主要基于线性、平稳的信号,对于非线性、非平稳的反射信号,容易出现信号失真；快速傅里叶变换和小波变换虽然都在一定程度上改善了傅里叶变换的缺点，但未能从根本上摆脱傅立叶分析方法的局限，且其中小波变换易出现因小波基的有限长导致能量泄露,从而使信号的能量频率时间的定量分布变的非常困难的情况。

发明内容

本发明的目的在于提出一种完全独立于傅里叶变换的，特别适用于非线性、非平稳信号处理的超声导波质量检测方法，其对于检测对象在检测过程中产生的弱反射和多次反射具有较强的提取能力。

本发明的目的还在于提出上述检测方法的应用。

本发明首先公开了如下的技术方案：

一种超声导波质量检测系统，其包括：发射低频超声导波的超声导波发射仪，与所述超声导波发射仪经传导介质相连的、固定于待检测目标端部的激励传感器，与所述激励传感器固定于同一端部且平行设置的接收传感器，与所述接收传感器传导介质相连的超声导波接收仪，与所述超声导波接收仪输出端相连的数据采集装置，与所述数据采集装置经数据传输通道相连的数据处理装置。

在一些优选实施方式中，所述低频超声导波选择20kHz-100kHz的纵向超声导波。

在一些优选实施方式中，所述超声导波发射仪和所述超声导波接收仪选自具有超声导波发射和接收功能的脉冲发生接收仪。

在一些优选实施方式中，所述传导介质为同轴电缆。

在一些优选实施方式中，所述传输通道为USB接口。

在一些优选实施方式中，所述数据采集装置为双通道数据采集卡。

在一些优选实施方式中，所述数据处理装置为计算机。

在一些优选实施方式中，所述传感器使用压电式传感器。

本发明进一步公开了一种超声导波质量检测方法，其包括：

S1：采集发射至被检测目标并返回后的低频超声导波信号；

S2：对所述返回的超声导波信号进行改进的自适应噪声完备集合经验模态分解(ICEEMDAN)，获得固有模态函数(IMF)；

S3：根据所述固有模态函数(IMF)获得被检测目标的质量情况。

在一些优选实施方式中，步骤S1中发射至被检测目标中的超声波为20-100kHz的纵向超声导波。

在一些优选实施方式中，步骤S2包括：

S21：在所述返回的超声导波信号中加入多组高斯白噪声，生成相应的多组新信号；

S22：对加入高斯白噪声后的信号序列进行经验模态分解，获得多个包络线均值，对所述多个包络线均值进行加总平均，得到当前阶残差；

S23：自所述返回的超声导波信号中减去当前阶残差，获得当前阶固有模态函数，并对当前阶残差加入新的多组高斯白噪声，再次进行经验模态分解，获得多个新的包络线均值，对所述多个新的包络线均值进行加总平均，得到新的残差，自所述当前阶残差中减去新的残差，得到新的固有模态函数；

S24：重复步骤S23，至当前阶残差的极值不超过两个，完成分解。

在一些具体实施方式中，步骤S2可包括：

S21：在所述返回的超声导波信号的序列s中加入I组经过经验模态分解(EMD)、并相乘信噪比的高斯白噪声ω_i(t)，得到I组新的信号序列s_i，如下：

s_i＝s+β₀E₁(ω⁽ⁱ⁾)

其中，ω⁽ⁱ⁾表示均值和单位方差分别为0和1的高斯白噪声，且i＝1，2，...，I；其中，i表示添加的白噪声编号，I表示添加的白噪声的总个数，

表示白噪声系数，用于调节添加的噪声与添加噪声后的残差之间所需的信噪比，ε₀表示添加的噪声的振幅，通常设为0.2，std(·)为标准差算子；E₁表示进行EMD分解得到的第1个固有模态函数(IMF)的算子；

S22：对所述新的信号序列s_i进行I次EMD分解，获得满足IMF判定条件的包络线均值，将该多个包络线均值进行加总平均，得到所述返回的超声导波信号的第一阶残差r₁，自所述返回的超声导波信号的序列s中减去当前阶残差，得到第一阶固有模态函数IMF₁，如下：

r₁＝<M(s_i)>

IMF₁＝s-r₁

其中，M(·)表示满足IMF判定条件的包络线局部均值的算子；<·>代表进行平均值计算；

S23：继续在第一阶残差中加入I组高斯白噪声，构建新的待分解信号r₁+β₁E₂(ω⁽ⁱ⁾)，并通过EMD分解计算局部均值得到第二阶残差r₂和固有模态函数IMF₂，如下：

r₂＝<M(r₁+β₁E₂(ω⁽ⁱ⁾))>

IMF₂＝r₁-r₂，

其中，β₁表示第二次添加白噪声时设置的白噪声系数，E₂表示对添加的白噪声进行EMD分解得到的第2个固有模态函数。

S24：同理，获得第k阶残差r_k和固有模态函数IMF_k，k＝3，...，K，K表示获得IMF的总个数，如下：

r_k＝<M(r_k-1+β_k-1E_k(ω(i)))>

IMF_k＝r_k-1-r_k

其中，β_k＝ε₀std(r_k)，表示第k-1次添加白噪声时设置的白噪声系数，E_k(·)为进行EMD分解得到的第k个固有模态函数(IMF)的算子；

S26：重复步骤S24，当残差r_k为单调函数时，停止计算，记录获得的固有模态函数分量。

其中，所述经验模态分解包括：

a.通过三次插值法获得当前信号序列的局部极大值点与局部极小值点的包络线，分别为上包络线和下包络线；

b.将上包络线和下包络线的幅值相加并进行平均，获得两者包络线均值；

c.在原始信号序列中减去所述均值，得到新的数据序列；

d.根据判别条件判断该新的数据序列是否为固有模态函数，若是，则此均值为目标均值，若不是，重复步骤a-c，进行k次筛选，并再次进行判断至符合判别条件。

其中，所述判别条件为：

同时满足以下条件的数据序列h_1k(t)为固有模态函数：

a.该序列中，极值点数目和过零点数目相等或者最多相差1个；

b.该序列中的任意点处，由局部极大值点和局部极小值点构成的两条包络线平均值为0。

在一些优选实施方式中，条件b通过限制标准差I_SD的值来替换，

其中，

t表示信号序列的时间点，T表示信号序列的总长度。

当I_SD取0.2～0.3时可认为满足条件b。

在一些优选实施方式中，所述步骤S3包括：

S31：获得超声导波在钢筋中的传播波速；

S32：获得超声导波在锚固密实段中的传播波速；

S33：根据所述传播波速和固有模态函数的峰值点确定被检测目标质量。

在一些具体实施方式中，所述步骤S33包括：

找出固有模态函数的峰值点，将其与被检测目标的不同端面和/或内部异常处进行对应；

根据不同突变点的时间及所述传播波速，获得被检测目标不同端面间的长度，和/或内部异常产生的位置；

本发明进一步公开了上述检测系统和/或检测方法在长条形工具上的应用。

在一些具体实施方式中，所述工具为锚杆。

所述锚杆可含有锚固段。

本发明具备以下有益效果：

本发明采用了独立于傅里叶变换的，特别适用于非线性、非平稳信号处理的信号处理过程，对检测目标在检测过程中产生的弱反射和多次反射具有较强的提取能力，可快速准确地获得目标的质量信息。

本发明可快速准确地检测到含有不同组成部分的目标的质量情况，如含有锚固段的锚杆。

本发明可快速准确地检测到目标的多种质量参数，如锚杆的长度、其中锚固段的长度、锚杆的密实度以及内部缺陷的位置、缺陷的异常程度等，属于无损检测领域。

本发明通过超声导波对目标进行检测，且信号处理高效准确，相对于现有的质量检测方法，检测效果更佳。

附图说明

图1为本发明的检测系统的装配示意图；

图2为实施例1中某一锚杆的端头反射波形图；

图3为实施例1中锚杆端头反射信号经过ICEEMDAN分解后的所有IMF分量图；

图4为实施例1中缺陷锚杆端头反射信号经过ICEEMDAN分解后的第二阶IMF分量图。

具体实施方式

以下结合实施例和附图对本发明进行详细描述，但需要理解的是，所述实施例和附图仅用于对本发明进行示例性的描述，而并不能对本发明的保护范围构成任何限制。所有包含在本发明的发明宗旨范围内的合理的变换和组合均落入本发明的保护范围。

通过以下过程进行锚固质量检测：

S1：连接超声导波发射仪，激发低频超声导波，并通过超声导波接收仪采集所述低频超声导波自被检测锚杆的返回信号。

其中，所用系统如附图1所示，其包括可向锚杆试样6发射及接收低频超声导波的超声导波发射、接收仪如脉冲发生接收仪1，对所述脉冲发生接收仪1输出数据进行采集、存储的采集装置如双通道数据采集卡2，对采集数据进行处理的数据处理装置如计算机3，设置于脉冲发生接收仪1发射端至锚杆试样6的发射通路间且固定于所述锚杆试样6端部的激励传感器4，设置于脉冲发生接收仪1接收端至锚杆试样6的接收通路间且固定于所述锚杆试样6端部的接收传感器5；其中，所述激励传感器和/或接收传感器优选为压电式传感器。超声导波发射仪和接收仪可通过具有良好屏蔽效果的同轴电缆与锚杆相连，其激发的低频超声导波可通过经粘结剂固定于锚杆端部的激励传感器得到传播，该传感器将所述低频超声导波传播至锚杆锚固中后，由固定于锚杆端部的接收传感器接收所述超声导波在缺陷处与锚固界面产生的发射波，所述反射波信号再经接收传感器同轴电缆传输至超声导波接收仪，超声导波接收仪将反射波信号再输出至信号数据采集装置，如双通道数据采集卡，采集装置对信号数据进行存储，并通过传输通道如USB接口输入数据处理装置，如计算机，根据计算机处理结果获得锚杆相关参数准确地检测到目标的多种质量参数，如锚杆的长度、其中锚固段的长度、锚杆的密实度以及内部缺陷如图示脱粘缺陷x的位置、缺陷的异常程度等。

其中，所述低频超声导波选择20kHz-100kHz的纵向超声导波，其产生过程可使用如下过程：

通过汉宁窗调制产生20kHz的5周期正弦函数波，由超声发射仪激发该正弦函数波组成的波包作为激励信号，该激励信号经过中心频率为30kHz的传感器传播到锚固锚杆中，即产生频率为20kHz-100kHz的低频纵向超声导波。

S2：对返回信号通过改进的自适应噪声完备集合经验模态分解(ICEEMDAN)进行处理。

其可进一步选择如下过程：

S21：在所述返回信号的序列s中加入I组经过经验模态分解(EMD)、并相乘信噪比的高斯白噪声ω_i(^t)，得到I组新的信号s_i。

s_i＝s+β₀E₁(ω⁽ⁱ⁾)

其中，ω⁽ⁱ⁾表示均值和单位方差分别为0和1的高斯白噪声，且i＝1，2，...，I；

表示白噪声系数，用于调节添加的噪声与添加噪声后的残差之间所需的信噪比，ε₀通常设为0.2，std(·)为标准差算子；算子E_k(·)为进行EMD分解得到的第k个固有模态函数(IMF)。

S22：对所述s_i进行I次EMD分解，获得满足IMF判定条件的包络线均值，将该多个包络线均值进行加总平均，得到所述返回信号的第一阶残差r₁，将所述返回信号的序列s中减去当前阶残差，得到第一阶固有模态函数IMF₁，如下：

r₁＝(M(s_i)>

IMF₁＝s-r₁

其中，M(·)表示满足IMF判定条件的包络线局部均值的算子；<·>代表进行平均值计算。

S23：继续在第一阶残差中加入I组高斯白噪声，构建新的待分解信号r₁+β₁E₂(ω⁽ⁱ⁾)，并通过EMD分解计算局部均值得到第二阶残差和IMF：

r₂＝<M(r₁+β₁E₂(ω⁽ⁱ⁾))>

IMF₂＝r₁-r₂

S24：同理，对于k＝3，...，K，，可以得到第k阶残差：

r_k＝<M(r_k-1+β_k-1E_k(ω⁽ⁱ⁾))>

式中，β_k＝ε₀std(r_k)。

S25：计算第k阶IMF：

IMF_k＝r_k-1-r_k

S26：重复步骤S24和S25，当残差r_k为单调函数时，停止计算，记录获得的IMF分量。

其中，所述经验模态分解包括：

a.通过三次插值法获得当前信号序列的局部极大值点与局部极小值点的包络线，分别为上包络线和下包络线；

b.将上包络线和下包络线的幅值相加并进行平均，获得两者包络线均值；

c.在原始信号序列中减去所述均值，得到新的数据序列；

d.根据判别条件判断该新的数据序列是否为固有模态函数，若是，则此均值为目标均值，若不是，重复步骤a-c，并再次进行判断至符合判别条件。

其中，所述判别条件为：

同时满足以下条件的数据序列h_1k(t)为IMF：

c.该序列中，极值点数目和过零点数目相等或者最多相差1个；

d.该序列中的任意点处，由局部极大值点和局部极小值点构成的两条包络线平均值为0。

其中，条件b在实际的处理过程中，是难以完全满足的，可通过限制标准差I_SD的值来替换，

其中，

当I_SD取0.2～0.3时可认为满足条件b。

S3：根据处理后的信号，结合导波在钢筋中和锚固密实段中的波速，获得锚杆的质量参数，如锚杆长度、锚固段长度、锚杆密实度和内部缺陷位置等。

其具体过程可采用如：

S31：获得超声导波在钢筋中的传播速度C_f。

具体可采用：

S310：对已知长度L的自由锚杆(即不含有锚固段的锚杆)根据步骤S1的过程进行超声导波检测，在其返回信号产生的时间-振幅谱上识别出其中激发波信号与锚杆底部的反射波包信号。

S311：根据峰值法，将时间-振幅谱上激发波的峰值与锚杆底部的反射波包的峰值对应的时间分别作为激发时间和反射时间，二者之差为双程旅行时间Δt_p，由此可得到导波在自由锚杆和上的单程传播时间

S312：通过自由锚杆的长度L与导波在自由锚杆上的单程传播时间的比值，得到普遍的锚杆自由段波速C_r。

S32：获得超声导波在锚固密实段中的传播速度C_b。

具体可采用：

S320：对已知长度L的全长锚固密实锚杆根据步骤S1的过程进行超声导波检测，在其返回信号产生的时间-振幅谱上识别出其中激发波信号与锚杆底部的反射波包信号。

S321：根据峰值法，将时间-振幅谱上激发波的峰值与锚杆底部的反射波包的峰值对应的时间分别作为激发时间和反射时间，二者之差为双程旅行时间Δt_b，由此可得到导波锚固密实锚杆上的单程传播时间

S322：通过锚杆的长度L与导波在锚固密实锚杆上的单程传播时间的比值，得到普遍的锚杆锚固密实段波速C_b。

S33：根据步骤S2获得的所述固有模态函数的峰值，确定锚杆质量。

具体如：

找出固有模态函数的峰值点，分析出所述峰值点与锚杆底部信号突变、锚固端面信号突变及内部异常信号突变(即内部缺陷导致的信号突变)的对应，确定其各自对应的反射时间t，将导波在钢筋中和锚固密实段中的波速分别与不同段的传播时间相乘，即可分别得到锚杆的长度、锚固段的长度和内部缺陷的位置，并可进一步获得锚杆的锚固段长度与总长度之比，即锚杆的密实度。

具体如：

对于锚固密实的锚杆，锚杆的长度可以下式计算：

式中，L表示锚杆长度，L₁为自由段长度，L₂为锚固体长度，C_f为超声导波在钢筋的波速，C_b为超声导波在锚固密实段的波速，t_b为锚固体上界面的反射时间，t_r为锚杆底部的反射时间，t_g为激发波激发时间。

对于锚固缺陷锚杆，缺陷处可看作钢筋，锚杆的长度按下式计算：

式中，L₃为缺陷的长度，L₂₁和L₂₂表示灌浆密实的锚固段长度，t_dd为锚固缺陷下界面的反射时间，t_ud为锚固缺陷上界面的反射时间。

实施例1

通过上述实施方式进行仿真实验：

检测目标为总长3m、直径为2cm的锚杆，锚固体长度为2.7m。

其根据上述实施方式产生的返回信号如图2所示，图中可以观察到两个信号特征，第一个为激励信号，第二个为锚杆底部的反射信号，其对应的峰值点时间为2.6ms。两个信号之间没有其他明显的反射信号，说明锚固体内无缺陷，锚固质量较好。

经过ICEEMDAN分解，得到了10阶IMF，如图3所示。

其中的第一阶IMF所表征的是激励信号的波形图，从第二阶IMF中可以观察到锚固体上界面的特征反射信号，点A表示反射信号的峰值，其对应的时间为0.29ms。

另外，获得超声导波在钢筋和锚固密实段的波速分别为5172m/s和2285m/s。

根据峰值点的时间和导波波速可计算得到锚杆长度为3.03m，其中锚固段长度为2.66m，锚杆密实度为87.8％。

锚杆长度测试相对误差为1％，锚固段测试相对误差为1.5％，锚固密实度测试相对误差为2.5％。

以上实施例仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例。凡属于本发明思路下的技术方案均属于本发明的保护范围。应该指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下的改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种超声导波质量检测系统，其特征在于：包括：发射低频超声导波的超声导波发射仪，与所述超声导波发射仪经传导介质相连的、固定于待检测目标端部的激励传感器，与所述激励传感器固定于同一端部且平行设置的接收传感器，与所述接收传感器经传导介质相连的超声导波接收仪，与所述超声导波接收仪输出端相连的数据采集装置，与所述数据采集装置经数据传输通道相连的数据处理装置。

2.根据权利要求1所述的质量检测系统，其特征在于：所述传导介质为同轴电缆，和/或所述传输通道为USB接口，和/或所述数据采集装置为双通道数据采集卡，和/或所述数据处理装置为计算机，和/或所述超声导波发射仪和所述超声导波接收仪选自具有超声导波发射和接收功能的脉冲发生接收仪，和/或所述传感器为压电传感器。

3.根据权利要求1所述的质量检测系统，其特征在于：所述低频超声导波为20kHz-100kHz的纵向超声导波。

4.一种超声波质量检测方法，其特征在于：包括：

S1：采集发射至被检测目标并返回后的低频超声导波信号；

S2：对所述返回的超声波信号进行改进的自适应噪声完备集合经验模态分解，获得固有模态函数；

S3：根据所述固有模态函数获得被检测目标的质量情况；

其中，发射至被检测目标的超声波为20-100kHz的纵向超声导波。

5.根据权利要求4所述的质量检测方法，其特征在于：步骤S2包括：

S21：在所述返回的超声导波信号中加入多组高斯白噪声，生成相应的多组新信号；

S22：对所述多组新信号进行经验模态分解，获得多个包络线均值，对所述多个包络线均值进行加总平均，得到当前阶残差；

S23：自所述返回的超声导波信号中减去当前阶残差，获得当前阶固有模态函数，并对当前阶残差加入新的多组高斯白噪声，再次进行经验模态分解，获得多个新的包络线均值，对所述多个新的包络线均值进行加总平均，得到新的残差，自所述当前阶残差中减去新的残差，得到新的固有模态函数；

S24：重复步骤S23，至当前阶残差的极值不超过两个，完成分解。

6.根据权利要求5所述的质量检测方法，其特征在于：其中，S21所述新信号通过以下计算模型得到：

s_i＝s+β₀E₁(ω⁽ⁱ⁾)，i＝1，2，...，I

其中，s_i表示该新信号，s表示所述返回的超声导波信号，ω⁽ⁱ⁾表示均值和单位方差分别为0和1的高斯白噪声，i表示添加的白噪声编号，I表示添加的白噪声的总个数，β₀表示白噪声系数，用于调节添加的噪声与添加噪声后的残差之间所需的信噪比，ε₀表示添加的噪声的振幅，通常设为0.2，std(·)为标准差算子；E₁表示进行EMD分解得到的第1个固有模态函数的算子；

和/或，S22-S23中获得的第k阶残差r_k和第k个固有模态函数IMF_k通过以下计算模型得到：

r_k＝<M(r_k-1+β_k-1E_k(ω⁽ⁱ⁾))>，k＝3，...，K，

IMF_k＝r_k-1-r_k

β_k＝ε_kstd(r_k)，

其中，M(·)表示满足固有模态函数判定条件的包络线局部均值算子，<·>表示求平均值，E_k(·)为进行EMD分解得到的第k个固有模态函数的算子，β_k＝ε₀std(r_k)，表示第k-1次添加白噪声时设置的白噪声系数。

7.根据权利要求5所述的质量检测方法，其特征在于：所述经验模态分解包括：

通过三次插值法获得当前信号序列的局部极大值点与局部极小值点的包络线，分别为上包络线和下包络线；

将所得上包络线和下包络线的幅值相加并进行平均，获得两者包络线均值；

在原始信号序列中减去所述均值，得到新的数据序列；

根据判别条件判断该新的数据序列是否为固有模态函数，若是，则此均值为目标均值，若不是，重复步骤a-c，并再次进行判断至符合判别条件；

其中，所述判别条件为：同时满足以下条件的数据序列h_1k(t)为固有模态函数：

a.该序列中，极值点数目和过零点数目相等或者最多相差1个；

b.该序列中的任意点处,由局部极大值点和局部极小值点构成的两条包络线平均值为0。

8.根据权利要求7所述的质量检测方法，其特征在于：所述条件b的判定为：若限制标准差I_SD为0.2～0.3，可认为满足条件b，其中：

t表示信号序列的时间点，T表示信号序列的总长度。

9.根据权利要求4所述的质量检测方法，其特征在于：所述步骤S3包括：

S31：获得超声导波在钢筋中的传播波速；

S32：获得超声导波在锚固密实段中的传播波速

S33：根据所述传播波速和固有模态函数的峰值点确定被检测目标质量；

优选的，所述步骤S33包括：

找出固有模态函数的峰值点，将其与被检测目标的不同端面和/或内部异常处进行对应；

根据不同突变点的时间及所述传播波速，获得被检测目标不同端面间的长度，和/或内部异常产生的位置。

10.权利要求1-3所述的检测系统和/或权利要求4-9所述的检测方法在锚杆检测上的应用。

Images