王海濤,史 振,蔣 俊,梁 斌,甄 理,賈曉媛,殷少華

(1.南京航空航天大學(xué)自動化學(xué)院,南京 210016;2.南京鍋爐壓力容器監(jiān)督檢驗(yàn)研究院,南京 210002)

超聲波檢測作為五大常規(guī)無損檢測技術(shù)之一,由于它的穿透能力強(qiáng)、對人體無害,而且具有被測對象范圍廣、檢測深度大、缺陷定位準(zhǔn)確及成本低等諸多優(yōu)點(diǎn)[1],已廣泛應(yīng)用于工業(yè)及高技術(shù)產(chǎn)業(yè)中,并且成為國內(nèi)外應(yīng)用最廣泛、使用頻度最高且發(fā)展較快的一種無損檢測技術(shù)。超聲成像技術(shù)是現(xiàn)代超聲無損檢測技術(shù)中一種令人矚目的新技術(shù),已成為定量檢測的重要手段[2-3]。分辨率是衡量超聲成像系統(tǒng)優(yōu)劣的最基本和最重要的指標(biāo)。近年來人們一直都致力于改善和提高超聲成像系統(tǒng)空間分辨率的研究[4],采用的主要措施有：超聲數(shù)字形成技術(shù)、合成孔徑成像技術(shù)、動態(tài)多頻率掃描和超聲相控陣技術(shù)等。

合成孔徑成像技術(shù)最早出自于軍事上的合成孔徑雷達(dá)技術(shù)[5],一般分為機(jī)載和星載兩種。與直接成像技術(shù)相比較,合成孔徑成像技術(shù)的特點(diǎn)是可以用小孔徑的換能器以及低的工作頻率來獲得比較高的方位分辨率,可以在近場區(qū)(即菲涅爾區(qū))工作,在圖像的重建過程中,易于應(yīng)用各種圖像處理技術(shù),用于不同成像對象,以提高圖像清晰度及信噪比[6]。筆者在合成孔徑成像理論的基礎(chǔ)上,對合成孔徑成像的波束形成進(jìn)行了仿真研究,從單陣元合成孔徑聚焦技術(shù)和多陣元合成孔徑聚焦技術(shù)兩個(gè)方面分析了參數(shù)對波束主瓣寬度和副瓣高度的影響。

1 合成孔徑成像基本原理

假設(shè)有一小孔徑換能器,其水平孔徑大小為d。超聲換能器按照“掃描-采樣-掃描”[7]的方式工作。超聲換能器沿水平方向移動掃描,每隔一個(gè)恒定的間隔Δx發(fā)射并且接收其回波信號,如圖1。掃描結(jié)束后,對所有點(diǎn)的回波信號進(jìn)行處理及圖像的構(gòu)建。

圖1 合成孔徑聚焦成像原理

在進(jìn)行圖像重構(gòu)時(shí),根據(jù)重建點(diǎn)到各個(gè)采樣點(diǎn)(換能器陣列單元)的不同距離而造成各個(gè)采樣點(diǎn)的回波信號有不同的延時(shí)τn,即相位發(fā)生了變化,對每個(gè)采樣點(diǎn)的回波信號的相位進(jìn)行校正后再疊加,即為數(shù)據(jù)的聚焦過程。在有缺陷的地方,回波信號同相疊加,信號增強(qiáng);在無缺陷的地方,回波信號是無序隨機(jī)的,故回波信號較弱。

對于一個(gè)小孔徑換能器,水平孔徑為d,其換能器的半功率點(diǎn)波束角為0.88乘以波長與水平孔徑即陣列單元的長度之比[8],因此,對一個(gè)均勻照射的合成陣列來說,雙程3 dB波束寬度為θ3dB=0.88(λ/d)。如果目標(biāo)點(diǎn)Q到基元換能器的垂直距離為R0,則基元換能器輻射聲束在R0處的合成孔徑長度為：

在合成孔徑中,各個(gè)等效基元是依次順序發(fā)射并接收,因此,各等效基元間的相位差是由發(fā)射到接收的“雙程”距離差引起的。與實(shí)際同樣尺寸的直線陣列形成的“單程”指向特性相比,合成孔徑各等效基元間的相位差就大了一倍。因此,合成孔徑線陣的半功率波束角Bs為：

由式(1)和(2)可得：

對應(yīng)于合成孔徑的方位分辨率為：

從式(4)可以看出,合成孔徑的方位分辨率只與換能器的孔徑有關(guān),與目標(biāo)點(diǎn)的位置和超聲波的波長無關(guān)。而直接成像系統(tǒng)方位分辨率ρa(bǔ)=0.84(λ/d)R,其中R為目標(biāo)離換能器的距離[6]。由此可知,要提高超聲成像系統(tǒng)的方位分辨率,通常只有兩種途徑：采用大孔徑換能器,或應(yīng)用高工作頻率的換能器,以得到短波長聲波。即使如此,方位分辨率還和目標(biāo)與換能器間的距離有關(guān),離得越遠(yuǎn),方位分辨率越差;而且在試件中,隨著頻率的升高,聲波的傳播損耗增大,其穿透深度越來越小;而且,大孔徑超聲換能器也很不方便操作。

由此可見,采用合成孔徑技術(shù)所帶來的好處是用小孔徑的實(shí)際基元換能器和較低的工作頻率,對目標(biāo)物體進(jìn)行高方位分辨率的觀察,且分辨率不隨目標(biāo)離換能器距離的增大而變差。

2 單陣元合成孔徑聚焦

聲場在成像場域的分布稱為波束形成。波束形成在整個(gè)超聲成像系統(tǒng)中處于核心的位置,它直接關(guān)系到成像的質(zhì)量。波束的主瓣寬度和副瓣高度理論上是判斷波束形成質(zhì)量的標(biāo)準(zhǔn)。主瓣越窄,成像的側(cè)向分辨率就越高;副瓣越低,成像的偽影就越低,對比度就越明顯。

對于單陣元合成孔徑,其每次發(fā)射和接收只有一個(gè)陣元是有效的,即一個(gè)陣元發(fā)射,同時(shí)該陣元進(jìn)行接收,每個(gè)陣元依次發(fā)射接收,其工作原理如圖2所示。

圖2 單陣元發(fā)射接收模式(共組成N個(gè)換能器陣列)

為了方便計(jì)算,假設(shè)上圖的換能器模擬陣列為基元運(yùn)動后的虛擬線性陣列,換能器初始位置為子孔徑1,發(fā)射并接收。接著移動到子孔徑2位置發(fā)射接收,一直移動到子孔徑N處。圖中的位置只是為了視覺上的便于分析,實(shí)際上它是按照水平方向移動的。合成后的虛擬線性陣列,通過對回波信號的處理,可以得到超聲波束[9]：

波束仿真參數(shù)的設(shè)置,由于K是一個(gè)固定值,對波束的形成無影響,可以不予考慮,D為陣元移動距離。仿真時(shí)假設(shè)K=1,陣元總數(shù)N=32,中心頻率f0=2.0 MHz,超聲聲速c=1.5 mm/μs,陣元間距分別取λ/2,λ/4和λ/8。不同陣元移動距離下波束如圖3所示。

在其他條件不變的情況下,對比D對波束形成的影響,當(dāng)D=a/2時(shí),發(fā)現(xiàn)在±90°處出現(xiàn)了很大的副瓣,這將造成偽影,甚至?xí)θ毕菰斐烧`判。這是由于D間距過大,達(dá)到了那奎斯特采樣定理的臨界值;當(dāng)D=λ/4時(shí),發(fā)現(xiàn)副瓣消失了;當(dāng)繼續(xù)減小D時(shí),可以發(fā)現(xiàn),雖然副瓣消失了,但隨之而來的主瓣寬度將變得很寬,這使得空間分辨率將減小。因此對于陣元移動間距的選擇并不是越小越好。

由合成孔徑的方位分辨率ρa(bǔ)=d/2可知,陣元孔徑d越小,分辨率應(yīng)當(dāng)越高。而事實(shí)上得出式(4)時(shí)做了菲涅爾近似,如果超聲換能器孔徑d無限制減小,波束開角過大,菲涅爾近軸條件將被破壞,此時(shí)式(4)的結(jié)果不再成立。理論分析表明,方位分辨率的極限值是λ/4[10],λ是工作波長,此時(shí)可得d=λ/2。仿真時(shí)用線性虛擬陣列,忽略了基元的尺寸,移動的距離D非常小,甚至小于陣元的尺寸,這對探頭部分的移動裝置有很高的要求。

3 多陣元合成孔徑聚焦

在合成孔徑中,陣元的尺寸較小,其輻射的聲功率不大,波束角很寬。因而回波信號較弱。為了提高合成孔徑聚焦的信噪比,人們提出了多陣元合成孔徑聚焦技術(shù),原理如圖4所示。

在這種工作情況下,發(fā)射時(shí)同時(shí)有幾個(gè)陣元同時(shí)發(fā)射,然后接收。將這幾個(gè)陣元移動一段距離,再次發(fā)射接收。假設(shè)虛擬陣元N,子孔徑為L個(gè)陣元,子孔徑的個(gè)數(shù)(移動的次數(shù))為K,則K=N-L+1。在這里為了簡化計(jì)算,假設(shè)子陣元的間距和移動距離相同,并且忽略子孔徑陣元間發(fā)射信號的相互影響。通過計(jì)算可得,其合成波束為：

圖3 不同陣元移動距離不單陣元孔徑的聚集波束

圖4 多陣元發(fā)射接收模式(共組成N 個(gè)換能器陣列)

圖5 不同子孔徑陣元合成孔徑聚集波束

為了和單陣元合成孔徑聚焦比較,參數(shù)設(shè)置為：陣元總數(shù)N=32,中心頻率f0=2.0 MHz,超聲聲速c=1.5 mm/μs,D=λ/2,L分別取4,8,16和32,不同子孔徑陣元數(shù)下多陣元合成孔徑聚焦波束見圖5。

比較圖3(a)和5(b)可以發(fā)現(xiàn),多陣元合成孔徑聚焦對旁瓣的抑制效果很明顯,可使圖像對比度提高明顯。當(dāng)子孔徑陣元個(gè)數(shù)小于總的虛擬陣元個(gè)數(shù)的一半時(shí),波束對旁瓣的抑制效果較好,且主瓣寬度較窄。從公式上也可以看出,當(dāng)L=32時(shí),此時(shí)的多陣元合成波束函數(shù)較單陣元合成波束的函數(shù)幾乎是一樣的,因此波束也應(yīng)該基本相同。這一點(diǎn)從圖3(b)和5(d)得到了驗(yàn)證。因此對陣元子孔徑的個(gè)數(shù)選取并不是越多越好,這和系統(tǒng)的總體陣元個(gè)數(shù)有關(guān)。

4 結(jié)語

合成孔徑成像技術(shù)是一種先進(jìn)的成像技術(shù),在雷達(dá)上廣泛應(yīng)用,在超聲無損檢測方面也有著巨大的應(yīng)用潛力。介紹了合成孔徑成像的基本原理,在此基礎(chǔ)上分別從單陣元合成孔徑和多陣元合成孔徑兩個(gè)方面做了波束仿真分析。對比可知道,多陣元合成孔徑比單陣元合成孔徑對旁瓣有較好的抑制效果。

[1]沈建中.超聲無損檢測的進(jìn)展[J].無損檢測,1998,20(2)：31-35.

[2]沈建中.超聲成像技術(shù)及其在無損檢測中的應(yīng)用[J].無損檢測,1994,16(7)：202-206.

[3]龐勇,韓焱.超聲成像方法綜述[J].華北工學(xué)院測試技術(shù)學(xué)報(bào),2001,15(4)：280-284.

[4]Bilgutay N M.Enhanced ultrasonic imaging with splitspectrum processing[J].IEEE Trans,1991,37(10)：1590-1592.

[5]劉永坦.雷達(dá)成像技術(shù)[M].哈爾濱：哈爾濱工業(yè)大學(xué)出版社,1999.

[6]孫寶申,沈建中.合成孔徑聚焦聲成像(一).應(yīng)用聲學(xué),1993,12(3)：43-48.

[7]王鄧志,羅斌,羅宏建.合成孔徑聚焦超聲成像系統(tǒng)研究[J].東北電力技術(shù),2002(12)：42-44.

[8]張忠敏,范錄宏,皮亦鳴,等.噪聲干擾對逆合成孔徑雷達(dá)成像的影響分析[J].電子對抗技術(shù),2005,20(5)：39-42,45.

[9]彭虎.超聲成像算法導(dǎo)論[M].合肥：中國科學(xué)技術(shù)大學(xué)出版社,2008.

[10]孫寶申,沈建中.合成孔徑聚焦聲成像(二)[J].應(yīng)用聲學(xué),1993,12(5)：39-45.