徐琰鋒，胡文祥，王浩

?

缺陷形貌的超聲成像方法

徐琰鋒，胡文祥，王浩

(同濟大學聲學研究所，上海 200092)

傳統(tǒng)超聲成像方法如合成孔徑聚焦技術、超聲衍射時差法，只能確定缺陷的位置，無法對整個缺陷形貌進行成像。逆時偏移方法基于波動方程、有限元或有限差分等方法計算聲源的正延拓聲場分布和接收回波的逆時延拓聲場分布，并對其采用互相關等成像條件，可以成像缺陷的形貌。針對平底狹縫鋁塊試樣，采用多陣元超聲換能器掃描，將接收數據用于成像，并與合成孔徑聚焦技術的成像比較，結果很好驗證了逆時偏移成像方法的優(yōu)越性。

形貌成像；逆時偏移；多元陣

0 引言

傳統(tǒng)超聲成像方法[1,2]僅僅能夠確定缺陷的位置及橫向尺寸，無法獲得缺陷的形貌信息。無法準確評估缺陷，為結構健康監(jiān)控帶來了極大的困難。比如，檢測存在平底豎直狹縫的試樣和內部有圓孔缺陷的試樣，只要兩者缺陷的上表面形狀尺寸相當，合成孔徑聚焦技術(Synthetic Aperture Focus Technique, SAFT)[3]及其他傳統(tǒng)成像方法得到的成像結果幾乎完全相同。其原因是，傳統(tǒng)的超聲成像方法實現方式較為簡單，并未利用聲場中的眾多信息，如多次反射波。相反，為了減小偽像，提高成像信噪比，傳統(tǒng)成像方法往往采用各種濾波方法壓制多次反射波。然而，缺陷的下表面及側面結構信息恰恰隱藏在多次反射波中。

有效利用多次反射波信息的成像方法需要滿足兩個條件。首先，成像方法必須基于波動方程。大多數傳統(tǒng)成像方法基于射線理論，如時域SAFT，通過計算表面接收換能器到達成像點的傳播時間，延時疊加獲得成像。使用射線理論計算多次反射波的傳播路徑是相當不切實際的。其次，成像方法必須能夠處理界面的反射與折射。頻域SAFT[4]是以波動方程為基礎，為了解決時域SAFT在多層介質成像的困難而被提出的。但是，其計算方法無法使聲波在界面上反射，多次反射波不能夠聚焦到缺陷上，卻可能在下層介質中形成偽像[5]。

源于地震成像領域的逆時偏移方法(Reverse Time Migration, RTM)[6]基于波動方程，使用有限元或有限差分等數值方法計算聲波在假設的聲場模型中傳播，因此可以處理多次反射波。上世紀80年代，該方法由Whitmore[7]，Baysal[8]，Mcmechan[9]等學者相繼提出，可以理解為將表面接收數據沿著時間減小的方向反向傳播的過程，由反射點散射的聲波重新聚焦回到散射點位置。

疊前逆時偏移的成像條件都是以Claerbout的時間一致性理論[10]為基礎，即反射界面存在于震源波場和接收波場在時間和空間上一致的地方，成像條件主要有三種[11]：激發(fā)時間成像條件，振幅比成像條件和互相關成像條件。激發(fā)時間成像條件，利用射線追蹤確定聲波到達成像點的時間，使用有限差分逆推聲場到這一時刻，將該位置在此時刻的位移作為成像值。激發(fā)時間成像條件難以處理多次反射波，成像質量較低，無法獲得缺陷形貌信息。對于復雜結構，激發(fā)時間成像條件較難得到準確的傳播時間。振幅比成像條件利用接收回波逆時延拓聲場，將其和聲源沿時間正向延拓聲場的比值表征反射系數，但是穩(wěn)定性不高。相對而言，互相關成像條件[12]應用廣泛，由Claebout[10]提出的零延遲互相關成像條件發(fā)展而來，使用聲源沿時間正向延拓聲場和接收回波沿時間逆向延拓聲場的互相關表征界面成像值。根據時間一致性理論，聲場互相關系數越大的位置存在缺陷或界面。

本文針對平底豎直狹縫的鋁塊試樣，采用線陣換能器掃描。并且對實驗數據進行逆時偏移，采用互相關條件成像。本文使用商業(yè)有限元軟件COMSOL Multiphysics計算聲場。聲場的計算模型的尺寸及參數與試樣均相同，但不存在任何缺陷。

1 逆時偏移技術

逆時偏移基于波動方程，分別計算聲波沿時間增大方向的正向傳播和沿時間減小方向的逆時傳播。根據時間一致性原理，應用成像條件表征成像值。一般，逆時偏移的聲場計算基于聲學波動方程[13]

采用多陣元換能器對缺陷試樣進行掃描，記錄下所有陣元的接收回波數據。接收數據應用于逆時偏移，其成像方法可以分為以下三個步驟：

(1) 假設合適的聲場模型，基于上述波動方程計算聲源沿時間正向傳播聲場。

(2) 將實驗測量或仿真得到的記錄回波數據作為聲源，在步驟(1)假設的聲場模型中計算沿時間逆向傳播聲場。

(3) 應用互相關成像條件，計算空間各點成像值。

逆時偏移方法的步驟(2)，聲場沿時間反向傳播類似于Fink的時間反轉概念[14]。時間反轉方法通過將多陣元相控陣探頭的接收信號反轉，并入射到原試樣，實現聲波在缺陷處聚焦的目的，從而增強缺陷回波的幅度。聲場的逆時延拓則是在無缺陷的假設模型下，模擬各陣元的接收回波沿時間減小的方向傳播的過程。本文的處理方法首先將接收信號時間反轉，然后將反轉后的信號輸入到各個陣元，計算聲波沿時間正向傳播的聲場。這里需要說明，將接收信號時間反轉后，計算其沿時間正向傳播的聲場，其計算結果和接收信號沿時間反向傳播的聲場是相同的。時間反轉是通過實驗的方式達到增強缺陷信號的目的，而逆時偏移則是通過數值計算方法模擬聲波的逆時傳播場，并最終獲得成像的后處理算法。兩者的思想相似，但實現方式及目的完全不同。

圖1 多陣元超聲測試示意圖

2 實驗

實驗系統(tǒng)采用Ultrasonix公司開放性多陣元相控陣研究平臺，并使用128通道的信號采集器SonixDAQ采集信號，采樣頻率為40 MHz。多陣元探頭128陣元，陣元寬度為0.2698 mm，間隔為0.3048 mm。圖2是平底裂縫鋁塊和多陣元探頭，鋁塊厚度為23 mm，裂縫寬度為0.5 mm，高為12 mm，相對于線陣孔徑尺度，鋁塊長度很大，因此左右邊界的反射波可以不予考慮。發(fā)射頻率為3 MHz，接收延遲時間為6.7 μs，接收信號的時間長度要能夠包含多次反射波，但是不同模式成分的多次反射波也是產生偽像的主要原因。因此，需要選擇合理的接收時間長度，這里接收總時間設置為16.5 μs，恰好能接收到鋁塊底部的二次反射波。掃描過程是第4+1個陣元發(fā)射聲信號(=0, 1, 2,…,31)，所有陣元同時接收回波。

圖2 實驗樣品及多陣元換能器

3 結果及討論

圖3是B型掃描圖，存在激發(fā)電干擾信號(標記3)。槽表面的回波(標記1)到達時間約為3.5 μs，鋁塊下界面的第一次回波(標記2)到達時間是7.3μs，7.3 μs后接收到的回波都是多次反射波。雖然多次反射波的幅度很小，但是包含缺陷的側面形貌信息。

圖4是第61個陣元(=15)發(fā)射，所有128個陣元接收到的回波信號。從圖4可以看到表面波的幅度很大，且與缺陷回波混雜在一起。表面波的存在會對成像結果造成極大的影響，因此在計算成像前需要先濾去表面波。表面等間隔陣元的表面波到達時間間隔相同，從圖4可以看到其視速度不變，因此表面波成分在頻率波數域呈傾斜直線分布(如圖5箭頭標記)。利用這一特點，采用FK濾波就可以有效地壓制接收信號中的表面波成分。FK濾波后的波形如圖6所示。從圖6可見，中表面波被抑制的同時，完好地保留了缺陷及界面的回波。

圖3 平底狹縫鋁塊的B型掃描圖。1-裂縫缺陷回波；2-鋁塊底部回波；3-電干擾信號。

圖4 接收信號。1-缺陷回波；2-底部界面回波；3-表面波

圖5 接收信號(圖4)的頻率波數域分布(箭頭標記的斜線表示表面波能量)

圖6 FK濾波。1-缺陷回波；2-底部界面回波

圖7 RTM 成像(矩形框標記裂縫缺陷)

合成孔徑聚焦技術(SAFT)廣泛應用于超聲成像領域，可以提高成像分辨率。這里，為了驗證逆時偏移方法獨特的優(yōu)點，將逆時偏移結果和SAFT的成像結果進行比較。如圖8是對圖3的B型掃描數據采用頻域SAFT方法得到的成像結果，箭頭標記的是SAFT成像結果的缺陷像，圖8中的橫線是由于過濾電干擾信號(圖3中的標記3)引起的信號跳變。SAFT雖然能夠準確地確定狹縫上端的位置，但無法獲得整個缺陷的形貌信息。顯然，傳統(tǒng)方法無法成像缺陷的形貌，逆時偏移成像技術可以彌補傳統(tǒng)成像方法在形貌成像上的劣勢。

圖8 SAFT成像(箭頭標記出裂縫位置)

目前，信噪比不佳是逆時偏移成像結果中存在的不足之處。圖8中，狹縫缺陷周圍分布著較多偽像。聲波在固體中傳播時，換能器接收到的信號存在縱波和橫波等模式轉換波。式(1)的標量波動方程忽略模式轉換波的影響，只能考慮一種模式波的傳播。因此，接收回波中的其他成分會以與實際不相等的傳播速度在模型中傳播。這些模式轉換波無法在缺陷位置聚焦，反而會形成偽像。在上述計算模型中，采用鋁塊的縱波速度作為聲傳播速度。那么回波中的橫波成分以及轉換波成分，被以縱波速度重新發(fā)射出去，這些回波無法在散射點位置發(fā)生聚焦，就形成了偽像，降低整個成像質量。也就是說，當以縱波速度作為聲傳播速度時，只有未發(fā)生過任何模式轉換的純縱波能夠準確地聚焦在狹縫位置，而其他模式波是產生偽像的主要原因。

4 結論

針對固體塊中的縱向分布裂縫缺陷，采用逆時偏移技術進行了成像方法研究。實驗采用多陣元換能器掃描平底狹縫的鋁塊，將接收數據應用于逆時偏移，并使用互相關成像條件表征成像值，得到的成像結果能夠清晰地分辨槽的位置及形貌，獲得了良好效果。

實驗的成像結果表明逆時偏移能夠較好成像缺陷的形貌，相對于傳統(tǒng)的合成孔徑成像方法有明顯的優(yōu)越性。但成像結果也存在偽像。造成偽像的主要原因是聲場計算均以縱波速度為聲波速度，接收回波中的轉換波模式被以與其不相同的聲波速度傳播，無法在缺陷位置聚焦。分離縱橫波模式是解決偽像的最直接方法，也是目前正在開展的工作。相關結果將另文討論。

當然，逆時偏移成像方法也存在其他缺點，如計算量巨大，極其耗時，高效率的聲場快速算法也將是需要考慮的課題。

[1] 彭虎. 超聲成像算法導論[M]. 合肥: 中國科學技術大學出版社, 2008.

[2] Stepinski T. Ultrasonic nondestructive inspection of solid objects[C]// AIP Conference Proceedings. 2012, 1433: 11.

[3] Spies M, Rieder H, Dillh?fer A, et al. Synthetic aperture focusing and time-of-flight diffraction ultrasonic imaging—past and present[J]. Journal of Nondestructive Evaluation, 2012, 31(4): 310-323.

[4] Busse L J. Three-dimensional imaging using a frequency-domain synthetic aperture focusing technique[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 1992, 39(2): 174-179.

[5] Olofsson T. Phase shift migration for imaging layered objects andobjects immersed in water[J]. IEEE Transactions on Ultrasonics, Ferroelectrics andFrequency Control, 2010, 57(11): 2522-2530.

[6] 牟永光. 地震數據處理方法[M]. 北京: 石油工業(yè)出版社, 2007.

[7] Whitmore N D. Iterative depth migration by backward time propagation[C]// 1983 SEG Annual Meeting. 1983.

[8] Baysal E, Kosloff D D, Sherwood J W C. Reverse time migration[J]. Geophysics, 1983, 48(11): 1514-1524.

[9] McMechan G A. Migration by extrapolation of time-dependent boundary values[J]. Geophysical Prospecting, 1983, 31(3): 413-420.

[10] Claerbout J F. Toward a unified theory of reflector mapping[J]. Geophysics, 1971, 36(3): 467-481.

[11] Chattopadhyay S, McMechan G A. Imaging conditions for prestackreverse-time migration[J]. Geophysics, 2008, 73(3): S81-S89.

[12] Kaelin B, Guitton A. Imaging condition for reverse time migration[C]// 2006 SEG Annual Meeting. 2006.

[13] Sava P, Hill S J. Overview and classification of wavefield seismic imaging methods[J]. The Leading Edge, 2009, 28(2): 170-183.

[14] Fink M. Time reversal of ultrasonic fields. I. Basic principles[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 1992, 39(5): 555-566.

[15] Müller S, Niederleithinger E, Bohlen T. Reverse time migration: a seismic imaging technique applied to synthetic ultrasonic data[J]. International Journal of Geophysics, 2012, Article ID: 128465.

Ultrasonic method for imaging appearance of defect

XU Yan-feng, HU Wen-xiang, WANG Hao

(Tongji University, Instituteof Acoustics, Shanghai 200092, China)

Traditional ultrasonic imaging methods, such as the synthetic aperture focusing technique and the time of flight diffraction, can find the position of a defect. However, it cannot image shape and appearance of the defect. The reverse time migration (RTM) method based on wave equation and implying the cross-correlation imaging condition between source propagation field and back-propagation field, which are calculated from measuring data by using finite element or finite differential method, can image shapes and appearances ofdefects. A multi-elements transducer is used to scan aluminum sample with a slot in our experiment. Comparing with synthetic aperture focusing technique, RTM imaging results show that it is a method superior to traditional one.

imaging of defect appearance; reverse time migration; multi-element array

O426.9

A

1000-3630(2014)-03-0222-05

10.3969/j.issn1000-3630.2014.03.008

2013-11-30;

2014-02-20

國家自然科學基金資助項目(11374230, 11104205)。

徐琰鋒(1988－), 男, 上海人, 博士研究生, 研究方向為超聲無損檢測。

胡文祥, E-mail: wxhu@#edu.cn