朱先祥 ， 許軍才 ， 倪福海 ， 丁 偉

（1.無(wú)損檢測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室（南昌航空大學(xué)），南昌 330063；2.河海大學(xué) 力學(xué)與材料學(xué)院，南京 210098；3.安徽中盛建設(shè)工程試驗(yàn)檢測(cè)有限公司，合肥 230051）

0 引言

沖擊回波法是混凝土結(jié)構(gòu)質(zhì)量檢測(cè)中一種常見(jiàn)無(wú)損檢測(cè)方法。該方法是根據(jù)縱波在結(jié)構(gòu)中的傳播特性對(duì)結(jié)構(gòu)質(zhì)量進(jìn)行檢測(cè)，其原理如圖1 所示[1-2]。利用鋼錘敲擊待測(cè)結(jié)構(gòu)表面，沖擊產(chǎn)生的振動(dòng)波在結(jié)構(gòu)內(nèi)部傳播過(guò)程中，當(dāng)遇到一定尺度的缺陷或邊界，因2 種材料的波阻抗不同而發(fā)生反射和透射現(xiàn)象，結(jié)構(gòu)沖擊表面上的質(zhì)點(diǎn)就會(huì)產(chǎn)生振動(dòng)響應(yīng)，這種振動(dòng)響應(yīng)以加速度的信號(hào)形式被傳感器采集到，再通過(guò)對(duì)該信號(hào)進(jìn)行快速傅里葉變換（Fast Fourier Transform，F(xiàn)FT），得到反映結(jié)構(gòu)信息的頻譜圖，最后對(duì)頻譜圖分析來(lái)判斷結(jié)構(gòu)的內(nèi)部質(zhì)量情況。

圖1 沖擊回波法基本原理Fig.1 Principle of impact echo method

但是，在實(shí)際的檢測(cè)過(guò)程中很難得到理想的縱波反射信號(hào)，通常伴有干擾信號(hào)如面波干擾、隨機(jī)噪聲以及體波散射等[3-4]。為了解決這一問(wèn)題，學(xué)者提出了不同時(shí)頻分析方法，如短時(shí)傅立葉變換、連續(xù)小波變換及Hilbert-Huang 變換（HHT），來(lái)對(duì)沖擊回波信號(hào)進(jìn)行分析和研究[5-7]。這些分析方法中，HHT 相對(duì)于短時(shí)傅立葉變換和連續(xù)小波變換能獲取更高的頻率分辨率。集成經(jīng)驗(yàn)?zāi)B(tài)分解（Ensemble Empirical Mode Decomposition，EEMD）是HHT 過(guò)程中的關(guān)鍵步驟。但是，EEMD 通過(guò)加入高斯白噪聲消除模態(tài)混疊問(wèn)題，卻帶來(lái)了殘余噪聲的影響，導(dǎo)致分解出的IMF 不是唯一，而加總平均又增大了相應(yīng)的計(jì)算成本[8-11]。為消除EEMD方法的局限性，Torres 提出了一種對(duì)噪聲有自適應(yīng)分解特性的完備集成經(jīng)驗(yàn)?zāi)B(tài)分解（Complementary Ensemble Empirical Mode Decomposition，CEEMD），該方法效果明顯好于早期的EEMD方法[12-13]。

針對(duì)傳統(tǒng)方法對(duì)沖擊回波信號(hào)分析的局限性，本研究將CEEMD 引入沖擊回波信號(hào)分析，建立基于CEEMD 的沖擊回波分析方法，通過(guò)仿真與試驗(yàn)證明研究CEEMD 對(duì)沖擊回波信號(hào)分解后頻譜特征和有效信號(hào)提取性能，為實(shí)際的沖擊回波檢測(cè)方法提供工程參考。

1 方法原理

1.1 沖擊回波原理

沖擊回波法對(duì)混凝土結(jié)構(gòu)缺陷處反射信號(hào)的有無(wú)及強(qiáng)弱來(lái)識(shí)別缺陷的存在，入射彈性波的反射率越高檢測(cè)的分辨率就越高。

根據(jù)彈性波傳播理論，彈性波在遇到不同介質(zhì)的交界面時(shí)，會(huì)發(fā)生反射和折射現(xiàn)象。對(duì)于雙層介質(zhì)，當(dāng)彈性波從介質(zhì)1 傳播到介質(zhì)2 時(shí)，其反射率 R計(jì)算式為：

式中： z1、 z2分別為介質(zhì)1、介質(zhì)2 的波阻抗，對(duì)于彈性波而言， z =ρv（ ρ為 介質(zhì)密度， ν為波速）。

對(duì)于夾層介質(zhì)的情況，如圖2 所示，反射率計(jì)算式為：

式中： k2為 介質(zhì)2 中的波數(shù)， f 為頻率， d為中間層厚度。

圖2 介質(zhì)含夾層時(shí)彈性波的傳播與反射Fig.2 Elastic wave propagation and reflection in the medium sandwich

沖擊回波測(cè)試結(jié)構(gòu)厚度是通過(guò)回波信號(hào)最大峰值的頻率來(lái)計(jì)算的，混凝土結(jié)構(gòu)底板厚度或缺陷位置有下列關(guān)系式：

式中：D 為底板厚度或缺陷位置，Vp為縱波波速，f 為最大峰值頻率，β=0.96 為經(jīng)驗(yàn)校正系數(shù)。

1.2 沖擊回波CEEMD 方法

CEEMD 方法在EEMD 的基礎(chǔ)上進(jìn)行改進(jìn)，是一種自適應(yīng)信號(hào)分解方法，有效地解決了EMD方法存在的分解不穩(wěn)定、模式混疊和端點(diǎn)效應(yīng)等問(wèn)題。CEEMD 方法基本原理利用添加一對(duì)振幅相同、相位相反的高斯噪聲來(lái)降低噪聲殘余，方法實(shí)現(xiàn)流程（圖3）如下：

首先，將振幅相等、相位相反的n 組輔助高斯噪聲加入原始沖擊回波信號(hào)中，從而產(chǎn)生2 個(gè)集合IMF：

式中：S 代表原始信號(hào)，N 為高斯噪聲，M1、M2加入噪聲后的正負(fù)對(duì)信號(hào)集合。

其次，對(duì)集合 M1和 M2中所有信號(hào)進(jìn)行EMD分解處理，分解后得到2 組集成的IMF 分量，即IMF1 和IMF2，IMF1 為正噪聲組集成的平均結(jié)果，IMF2 為負(fù)噪聲組集成的平均結(jié)果。

圖3 CEEMD 方法沖擊回波信號(hào)分解流程圖Fig.3 Decomposition flow chart of CEEMD

最后，算法分解為2 組IMF 分量，對(duì)IMF1、IMF2 集成平均，結(jié)果為：

式中：cj代表CEEMD 分離得到的第j 個(gè)IMF 分量，cij為i 信號(hào)分離所得的編號(hào)為j 的IMF 分量。

CEEMD 方法通過(guò)在沖擊回波原信號(hào)中加入正負(fù)對(duì)白噪聲，消除了后期進(jìn)行信號(hào)重構(gòu)時(shí)出現(xiàn)的殘余輔助噪聲，不僅繼承了EEMD 方法的優(yōu)點(diǎn)同時(shí)提高了計(jì)算效率，而且減小了添加輔助白噪聲帶來(lái)的重構(gòu)誤差[14]。

2 沖擊回波檢測(cè)信號(hào)分析

2.1 沖擊回波模擬信號(hào)分析

設(shè)計(jì)的仿真模型尺寸為200 mm×200 mm×20 mm 的混凝土板，在板內(nèi)部有一尺寸為30 mm×30 mm×16 mm 夾層空隙(圖4a)。模擬中計(jì)算參數(shù)混凝土密度為2500 kg/m3，彈性模量為37.8 GPa，泊松比為0.2[15]。相應(yīng)的夾層和底部的回波的頻率分別為12.3、9.8 kHz。利用有限元方法對(duì)該混凝土板沖擊回波檢測(cè)模擬，模擬出的沖擊回波信號(hào)如圖4b 所示。

根據(jù)圖4 中的沖擊回波信號(hào)，利用本研究提出的CEMMD 方法對(duì)沖擊回波信號(hào)進(jìn)行分析，將原始信號(hào)分解為相應(yīng)的IMF，原信號(hào)分解后的結(jié)果見(jiàn)圖5。

圖4 沖擊回波模擬試驗(yàn)Fig.4 Impact echo simulation testing

圖5 CEEMD 分解沖擊回波模擬信號(hào)的IMFFig.5 IMF decomposed on impact echo simulation signal by CEEMD

根據(jù)圖5 結(jié)果，CEEMD 方法能夠?qū)⒃紱_擊信號(hào)分解為不同的IMF，模態(tài)函數(shù)IMF1～I(xiàn)M10。從IMF1～I(xiàn)MF10 特征來(lái)看，模態(tài)函數(shù)振動(dòng)頻率排序逐漸減小，其中，IMF1～I(xiàn)MF4 主要以高頻為主的成份，IMF5～I(xiàn)MF10 反映了沖擊回波信號(hào)結(jié)構(gòu)振動(dòng)特征，IMF1～I(xiàn)M10 頻率特征如圖6 所示。

圖6 沖擊回波模擬信號(hào)的IMF1～I(xiàn)MF10 頻譜圖Fig.6 IMF1～I(xiàn)MF10 spectrogram on impact echo simulation

圖6 的結(jié)果顯示，IMF 能夠很好與沖擊回波信號(hào)特征相吻合。IMF1～I(xiàn)MF4 主要體現(xiàn)的是振動(dòng)沖擊回波信號(hào)高頻成份。而IMF5、IMF6 對(duì)應(yīng)的峰值頻率分別為12.3、9.8 kHz，該頻率是來(lái)自?shī)A層空隙和板底反射信號(hào)的頻率；由此可見(jiàn)，CEEMD方法能夠很好地將板底和夾層空隙處的缺陷很好分離出來(lái)，并能對(duì)缺陷位置進(jìn)行確定。

2.2 沖擊回波工程檢測(cè)信號(hào)分析

以美國(guó)弗吉尼亞州Gainesville 15 號(hào)公路的沖擊回波現(xiàn)場(chǎng)測(cè)試為例（圖7a），現(xiàn)場(chǎng)檢測(cè)過(guò)程中檢測(cè)位置進(jìn)行了鉆孔取芯，芯樣如圖7b 所示[15]。芯樣路面面層厚度為21.6 cm，面層縱波波速大約為4000 m/s，由厚度和頻率的關(guān)系式(3)，計(jì)算得到底部反射波頻率為9 kHz。但由于面波和噪聲干擾，圖7a 中的檢測(cè)信號(hào)的頻譜很難識(shí)別出底面的反射波頻率。

圖7 沖擊回波現(xiàn)場(chǎng)試驗(yàn)Fig.7 Impact echo on site testing

對(duì)該現(xiàn)場(chǎng)檢測(cè)的沖擊回波信號(hào)采用CEEMD方法進(jìn)行分析，將該檢測(cè)信號(hào)分解為相應(yīng)的IMF，原始信號(hào)分解后的結(jié)果見(jiàn)圖8。

根據(jù)圖8 結(jié)果，CEEMD 方法能夠?qū)⒃紱_擊信號(hào)分解為IMF1～I(xiàn)M10 特征模態(tài)函數(shù)。從IMF1～I(xiàn)MF10 特征來(lái)看，模態(tài)函數(shù)振動(dòng)頻率排序逐漸減小，其中IMF1～I(xiàn)MF4 主要以高頻為主的成份，主由噪聲信號(hào)和干擾信號(hào)組成。IMF5～I(xiàn)MF10 反映了沖擊回波信號(hào)結(jié)構(gòu)振動(dòng)特征，IMF7～I(xiàn)MF10 成分為單調(diào)函數(shù)為主，與上節(jié)仿真沖擊回波信號(hào)的分析結(jié)果一致。對(duì)IMF1～I(xiàn)M10 特征模態(tài)函數(shù)進(jìn)行頻譜變換，其頻率特征如圖9 所示。

圖9 顯示出IMF1～I(xiàn)M10 特征模態(tài)函數(shù)變化規(guī)律，IMF1～I(xiàn)MF4 主要體現(xiàn)的是振動(dòng)沖擊回波信號(hào)高頻成份，是由噪聲和面波干擾形成的高寬頻信號(hào)。而IMF5、IMF6 對(duì)應(yīng)的峰值頻率分別為17.5、9.0 kHz，可判斷出該頻率是來(lái)自?shī)A層空隙和底面反射信號(hào)形成。根據(jù)頻率與厚度的函數(shù)關(guān)系式(3)，脫空位置在10.5 cm 處，圖7b 中芯樣脫空位置大致在10.0 cm 附近；而底面反射波峰值頻率也為9.0 kHz，與實(shí)際情況相符。由此可見(jiàn)，沖擊回波信號(hào)受干擾波影響時(shí)，CEEMD 方法能正確識(shí)別出結(jié)構(gòu)中的缺陷位置和混凝土厚度。

3 結(jié)論

1）CEEMD 方法能夠有效地將沖擊回波信號(hào)從含有嚴(yán)重的面波和噪聲干擾中分解為相應(yīng)的IMF 分解為相應(yīng)的IMF，IMF 模態(tài)函數(shù)頻率由高到低按序排列。

2）IMF 對(duì)應(yīng)的頻譜顯示，高頻成份一般為噪聲干擾成份為主的面波信號(hào)，處于中間的IMF 則為夾層處缺陷和板底反射成份的信號(hào)。

3）分析缺陷和板底反射信號(hào)的對(duì)應(yīng)IMF，確定缺陷和板底反射信號(hào)頻率，再由頻率準(zhǔn)確計(jì)算內(nèi)部夾層缺陷位置和混凝土板厚。

圖8 CEEMD 對(duì)現(xiàn)場(chǎng)測(cè)試信號(hào)分解的IMFFig.8 IMF decomposed on site testing by CEEMD

圖9 現(xiàn)場(chǎng)檢測(cè)信號(hào)的IMF1～I(xiàn)MF10 頻譜圖Fig.9 IMF1～I(xiàn)MF10 spectrogram on site testing signal