李茂東， 林曉輝， 楊 波， 洪曉斌

(1. 廣州特種承壓設(shè)備檢測研究院，廣東 廣州 510663; 2. 華南理工大學機械與汽車工程學院，廣東 廣州 510640)

0 引 言

工業(yè)用金屬材料表面通常覆蓋有非金屬防腐層，如各種金屬板材以及石化行業(yè)中用到的內(nèi)襯防腐管道等，隨著使用時間增加，防腐層容易產(chǎn)生各種缺陷，其中界面分層是常見的一種缺陷?？傮w而言，分層缺陷的產(chǎn)生主要來自生產(chǎn)過程和使用過程，生產(chǎn)過程中，非金屬防腐層與金屬基材之間的結(jié)合通常要經(jīng)過加熱加壓保溫的步驟，若在此過程中出現(xiàn)加壓復合時壓力過小、防腐層厚度不均勻等，則可能導致防腐層材料不能充分產(chǎn)生塑性變形而殘存較大應(yīng)力，進而在冷卻時產(chǎn)生分層缺陷。而在使用過程中，可能存在反復的加載、卸載、升溫、降溫等復雜工況，會在應(yīng)力集中的區(qū)域形成很大的交變應(yīng)力，如不加以檢測防范，會使得非金屬防腐層產(chǎn)生分層缺陷，金屬基材將失去保護而遭到腐蝕甚至發(fā)生斷裂。

目前對防腐層界面損傷檢測可分為非超聲波防腐層界面損傷監(jiān)測方法和超聲波防腐層界面損傷監(jiān)測方法兩大類，其中非超聲波防腐層界面損傷監(jiān)測方法主要包括紅外熱波檢測方法、漏磁檢測方法、渦流檢測法等。如Han J S等[1]通過加熱被檢測區(qū)域，采用紅外熱像儀記錄試件表面溫度場，通過數(shù)據(jù)分析軟件對紅外圖像進行處理后得到涂層缺陷程度。J.Jun等[2]使用一個掃描式的漏磁探頭結(jié)合霍爾傳感器陣列來檢測復合材料中的裂紋情況，并對其檢測能力進行討論。JH.Rose等[3]提出了一種基于頻率變化阻抗的渦流檢測線圈，用來檢測金屬涂層的厚度和導電性。紅外熱波方法若要達到一定分辨率，對儀器設(shè)備要求較高；漏磁檢測需要在鐵磁材料上進行，不適用于非金屬材料的分層檢測；渦流檢測對金屬表面或近表面的缺陷有很高的檢出靈敏度，但是對于非金屬的熱障涂層損傷卻無法進行檢測。

超聲波檢測法是廣泛用于材料探傷的常用方法，也是最早用于復合材料無損評價的方法之一。它主要利用復合材料本身或其缺陷的聲學性質(zhì)對超聲波傳播的影響來檢測材料內(nèi)部和表面的缺陷，如氣泡、分層、裂紋、脫粘、貧膠等。超聲波應(yīng)用于防腐層界面損傷檢測主要有聲發(fā)射法、超聲脈沖回波法、超聲導波法等。如Gallego A等[4]基于聲發(fā)射技術(shù)，對在熱鍍鋅的鋼樣品表面防護層和本體金屬的結(jié)合程度進行評價。日本慶應(yīng)義塾大學的Y Suga等[5-6]利用超聲波方法研究了噴敷層脫粘現(xiàn)象，并探索了影響超聲波檢測效果相關(guān)因素，同時將這種方法應(yīng)用于檢測陶瓷噴敷層材料脫層現(xiàn)象實驗，發(fā)現(xiàn)處于熱沖擊試驗下產(chǎn)生的脫層現(xiàn)象可以利用這種方法有效檢測出來。韓國先進科技學院的H.Sohn等[7]探索了使用非接觸超聲導波檢測多層復合材料中隱藏的脫層現(xiàn)象可行性?，F(xiàn)有的超聲檢測方法大多為檢測復合材料的粘接質(zhì)量，而較少涉及對金屬和非金屬分層缺陷的研究。本文采用超聲導波檢測方法，搭建試驗平臺，利用超聲導波設(shè)備對分層缺陷進行檢測并提取有效數(shù)據(jù)，接著基于小波變換，提出利用小波能量值檢測分層缺陷損傷程度的評價方法。

1 非金屬防腐層分層缺陷測量原理

1.1 超聲導波在雙層介質(zhì)中的傳播規(guī)律

由于界面的存在和材料性質(zhì)的不同，超聲導波在雙層介質(zhì)中的傳播不同于傳播在無限大介質(zhì)或者層板中。建立雙層板模型如圖1所示。其中，ρ，λ和μ分別為介質(zhì)密度和Lame參數(shù)，h表示介質(zhì)層厚度，x軸、y軸、z軸如圖所示。在分層介質(zhì)中，對于按eiw簡諧變換的位移場，必定滿足Navier位移方程，即有：

圖1 雙層板模型

基于Navier位移方程結(jié)合Helmholtz分解可計算得到超聲導波在多層介質(zhì)中傳播的位移場勢函數(shù)通解[8]：

其中為任意參數(shù)；每一個勢函數(shù)可以看作向下和向上傳播的平面波兩項，相當于每一層板中包括兩個縱波和兩個剪切波的4個平面體積波。從上述理論分析可知，超聲導波在多層介質(zhì)中傳播存在多模態(tài)，多方向等較為復雜特點，難以利用特定模態(tài)對分層缺陷進行準確檢測。同時，隨著層數(shù)的增加，導波在介質(zhì)中傳播的模態(tài)大幅增加，從能量方面考慮，層數(shù)越多，在相同距離內(nèi)，導波發(fā)生的折射和模態(tài)轉(zhuǎn)換越多，即能量消耗越嚴重，因此，可根據(jù)接收信號能量大小對非金屬防腐層分層缺陷進行判斷。

1.2 基于小波包變換的信號能量值評估

小波變換是當前應(yīng)用數(shù)學中發(fā)展迅速的領(lǐng)域，它是以局部化函數(shù)形成的小波基作為基底展開的，具有許多特殊的性能和優(yōu)點，小波包變換是基于小波變換的進一步發(fā)展，能夠提供比小波變換更高的分辨率。對于信號函數(shù)x(t)使用小波Ψ(t)為基的小波變換[9]的定義為：

其中x(t)∈L2(R)，Ψ(t)是母小波，a，b分別為小波變換的尺度因子和時間因子，滿足條件：a>0。基于式（6）可以得到相對應(yīng)的小波重構(gòu)公式為

小波變換是一種積分變換，它的實質(zhì)在于將L2(R)空間中的任意函數(shù)表示為其在具有不同伸縮因子a和 平移因子b的 Ψ(t)之上的投影的疊加。通過調(diào)整伸縮因子a和 平移因子b可以得到具有不同時域?qū)挾鹊男〔ㄒ云ヅ湓夹盘柕娜我馕恢?，達到對信號的時頻局部化分析的目的。在這里，用一個信號的3層小波包分解樹說明小波分解的過程。如圖2所示，其中A表示低頻，D表示高頻，末尾的序號數(shù)表示小波分解的層數(shù)。由圖可知，分解級數(shù)越大，也就是選擇的小波包尺度越大，小波包系數(shù)對應(yīng)的空間分辨率就越低，利用這一點，可以在不同的空間分辨率上進行分析。

圖2 3層小波包分解樹

本文在基于小波能量值的分層缺陷評估中，將小波包變換引入對分層缺陷信號分析，通過n階小波包對接收信號進行分解，得到2n個信號小波子集分解后的子集可以表示如下：

其中j是頻率范圍，m指總的采集的數(shù)據(jù)，定義第i個經(jīng)過分解的能量信號如下：

因此，第i個能量集合定義為：

基于上面公式推導，定義第i個能量集合的小波能量值為：

非金屬防腐層分層缺陷損傷可以用信號的小波能量值Ew來評估。由于超聲導波在雙層介質(zhì)中的傳播損耗更大，因此理論上采用基于小波能量值的分層缺陷評估方法中，小波能量值越小，表明信號能量損耗越大，進而說明分層缺陷損傷越大。

2 基于小波能量值的分層缺陷評估實驗平臺

為檢測出非金屬防腐層分層缺陷，基于超聲波在雙層板中的傳播理論，本文首先搭建基于小波能量值分層缺陷損傷評估方法的實驗平臺，如圖3、圖4所示。實驗平臺由上位機、波形發(fā)生器、壓電放大器、試樣和PCI-20614數(shù)據(jù)采集卡組成。實驗過程中，首先在上位機編制發(fā)射信號，通過波形發(fā)生器生成信號，接著信號放大器將信號放大50倍之后通過發(fā)射探頭激勵試樣，經(jīng)過分層缺陷并被接收探頭接收，接收的信號通過采集卡采集并傳輸?shù)缴衔粰C進行進一步的信號處理。

圖3 基于小波能量值的分層缺陷評估實驗平臺

圖4 基于小波能量值的分層缺陷評估實驗流程圖

其中試樣的構(gòu)成如圖5所示，由上層PC板和下層鍍鋅鋼板組成，兩層介質(zhì)之間通過水緊密連接。利用兩塊PC板之間的空隙模擬分層缺陷，通過控制移動兩塊PC板可以控制分層缺陷的位置和大小。PC板和鋼板的厚度均為3 mm。經(jīng)過測試不同頻率和入射角度的超聲斜探頭，發(fā)現(xiàn)不同的頻率和入射角度對實驗結(jié)果的總體趨勢影響不大，因此實驗選擇常用的2.5Z10×10A60的超聲斜探頭作為發(fā)射和接收探頭。發(fā)射探頭與接受探頭位置固定且距離為200 mm。分層缺陷產(chǎn)生時，在缺陷產(chǎn)生的區(qū)域，非金屬防腐層與金屬基材發(fā)生分離，該實驗平臺模擬的分層缺陷與實際的分層缺陷基本符合，同時本文中涉及的基于增益值的分層缺陷損傷評估方法和基于小波能量值的分層缺陷評估方法原理是利用接收信號能量在正常區(qū)域和缺陷區(qū)域的損耗不同來實現(xiàn)對分層缺陷的檢測，在這一方面模擬的分層缺陷與實際的分層缺陷對信號能量所造成的影響是一致的，因此模擬分層缺陷的實驗結(jié)果具有較高的可信度。

圖5 分層缺陷試驗實物圖

在超聲導波檢測領(lǐng)域中，由于經(jīng)過漢寧窗調(diào)制的正弦信號具有頻散特性弱、能量損失低、信號失真低等優(yōu)點，因此本實驗采用10周期漢寧窗正弦調(diào)制信號檢測分層缺陷，同時，基于探頭的諧振頻率，選擇調(diào)制信號的中心頻率為2 MHz。該漢寧窗正弦調(diào)制信號的具體時域頻域圖如圖6和圖7所示。

圖6 漢寧窗正弦調(diào)制信號時域圖

3 實驗結(jié)果及分析

3.1 基于增益值的分層缺陷損傷評估實驗

超聲檢測儀器的增益值是衡量回波幅值的重要指標，增益值越大說明回波幅值越小，即信號能量越低。利用增益值的這一特點，本文首先探索了一種基于增益值的分層缺陷損傷評估方法。實驗平臺中試樣和探頭與第2節(jié)中一致，發(fā)射和接收部分則由常規(guī)超聲儀器代替，如圖8所示。

基于增益值的分層缺陷損傷評估方法的實驗按照模擬的分層缺陷寬度大小分為21組進行，從分層缺陷寬度為0 mm開始到200 mm，每隔10 mm一組。對于每一組缺陷，通過移動發(fā)射探頭和接收探頭使得缺陷回波幅值達到最高，接著利用超聲檢測儀器的自動增益功能將回波幅值調(diào)到80%的高度，記錄下此刻儀器顯示的增益值。實驗過程中，始終保持分層缺陷位于發(fā)射探頭和接收探頭中間，實驗結(jié)果如圖9所示。

由圖可知，隨著分層缺陷寬度的增加，使接收信號達到80%高度的儀器增益值整體呈下降趨勢，結(jié)果表明利用接收信號的能量變化可以判斷分層缺陷的大小。然而，由于實際應(yīng)用這種基于增益值的分層缺陷損傷檢測方法時，需要通過移動發(fā)射探頭和接收探頭使得缺陷回波幅值達到最高，這個過程中產(chǎn)生的誤差影響其檢測精度。此方法只能在分層缺陷較大的情況下才能將其有效分辨出來，從圖中可以看出，當分層缺陷為0 ，10 ，20 mm時，增益值分別為47.2 ，47.3 ，47 dB，即在此范圍內(nèi)，增益值并沒有隨著分層缺陷的增大而單調(diào)減小，此時利用超聲儀器增益值的方法難以有效分辨出分層缺陷。下面利用小波能量法檢測非金屬防腐層分層缺陷并進行比較。

圖9 基于增益值的分層缺陷損傷檢測實驗結(jié)果

3.2 基于小波能量值的分層缺陷評估實驗

本文基于小波能量值的分層缺陷評估方法總共進行了41組實驗，為了更加細致地觀察不同分層缺陷寬度下接收信號的小波能量值的規(guī)律，從分層缺陷寬度為0 mm開始到200 mm每隔5 mm取一組缺陷進行實驗，每組實驗進行10次并取接收信號的平均值進行3層的小波包分解，如圖2所示。用A表示低頻，D表示高頻，末尾的序號數(shù)表示小波分解的層數(shù)，則得到小波包分解信號如圖10所示。

實驗過程中始終保持發(fā)射探頭和接收探頭位置不變。對接收信號利用Matlab程序提取小波能量值，記錄這41組數(shù)據(jù)，結(jié)果如圖11所示。

從圖中實驗結(jié)果可以發(fā)現(xiàn)，隨著分層缺陷寬度的增大，接收信號的小波能量值單調(diào)增大，表明分層缺陷寬度越大，超聲導波傳播過程中發(fā)生的反射、折射和損耗就越小，因此小波能量值越大。不同于基于增益值的評估方法，由于在此方法中，探頭始終保持固定，不需要每一組實驗都移動發(fā)射接收探頭，避免了移動過程中可能產(chǎn)生的誤差。在基于增益值的評估方法實驗結(jié)果中，增益值總體呈減少趨勢，但具體在相鄰的某幾組實驗中，并不一定是單調(diào)減少的。而在基于小波能量值的分層缺陷檢測實驗中，可發(fā)現(xiàn)接收信號的小波能量值具體到相鄰的實驗中，都是隨著缺陷寬度而單調(diào)增加。同時，對于分層缺陷寬度在0～20 mm的范圍內(nèi)，由于較多誤差存在使得基于增益值的分層缺陷評估方法并不能很好地反映出缺陷問題，對較小分層缺陷容易產(chǎn)生漏檢的情況。如表1所示，利用基于小波能量值的分層缺陷評估方法，即使在缺陷寬度為0～20 mm的情形下，也具有比較明顯的趨勢，雖然小波能量值的增長絕對值不大，但隨著分層缺陷寬度的增加仍有明顯的單調(diào)上升趨勢，即通過小波能量值可判斷出較小分層缺陷的存在。

圖10 分層缺陷寬度為5 mm時接收信號的三層小波包分解信號

圖11 基于小波能量值的分層缺陷評估實驗結(jié)果

表1 較小分層缺陷下的小波能量值

4 結(jié)束語

本文為了實現(xiàn)對非金屬防腐層的分層缺陷檢測，首先基于超聲導波在多層介質(zhì)中的傳播理論，建立了超聲導波在雙層板狀介質(zhì)中的傳播模型，提出了基于小波能量值的分層缺陷檢測方法，搭建了基于小波能量值的分層缺陷檢測平臺。實驗部分首先利用超聲儀器增益值的指標對分層缺陷寬度進行檢測，結(jié)果表明隨著分層缺陷寬度的增加，達到相同回波幅度的增益值逐漸下降，但由于存在操作誤差和儀器誤差，在局部范圍內(nèi)無法識別出分層缺陷。接著利用基于小波能量值的分層缺陷檢測方法對分層缺陷做進一步的檢測，實驗結(jié)果表明接收信號的小波能量值這一指標能更好地實現(xiàn)對分層缺陷的檢測，特別是在缺陷寬度較低的情況下，采用小波能量法能更好滿足分層缺陷檢測需求。

[1]HAN J S, PARK J H. Detection of corrosion steel under an organic coating by infrared photography[J]. Corrosion Science,2004, 46(4): 787–793.

[2]JUN J, LEE J, PARK D. NDT of a nickel coated inconel specimen using by the complex induced current–magnetic flux leakage method and linearly integrated hall sensor array[J].Journal of the Korean Society for Nondestructive Testing, 2007,27(5): 375–382.

[3]SETHURAMAN A, ROSE J H. Rapid inversion of eddy current data for conductivity and thickness of metal coatings[J].Journal of Nondestructive Evaluation, 1995, 14(1): 39–46.

[4]GALLEGO A, GIL J F, VICO J M, et al. Coating adherence in galvanized steel assessed by acoustic emission wavelet analysis[J]. Scripta Materialia, 2005, 52(10): 1069–1074.

[5]SUGA Y, LIAN D, KURIHARA S, et al. Effect of some measuring factors on detection of sprayed coating delamination in ultrasonic testing method[C]// Japan Welding Society, 1995:246–247.

[6]LIAN D, SUGA Y, SHOU G, et al. An ultrasonic testing method for detecting delamination of sprayed ceramic coating[J]. Journal of Thermal Spray Technology, 1996, 5(2):128–133.

[7]SOHN H, DUTTA D, YANG J Y, et al. Delamination detection in composites through guided wave field image processing[J].Composites Science & Technology, 2011, 71(9): 1250–1256.

[8]ROSE J L. Ultrasonic guided waves in solid media[J].Cambridge University Press , 2014.

[9]PUN C M, LEE M C. Log–polar wavelet energy signatures for rotation and scale invariant texture classification[J]. Hydrology& Earth System Sciences Discussions, 2003, 25(5): 590–603.