非接觸式光聲成像技術及其應用

陳重江，周王婷，胡毅成，楊思華

華南師范大學生物光子學研究院 激光生命科學研究所、暨激光生命科學教育部重點實驗室，廣東 廣州 510631

引言

目前，在光聲成像技術中，光聲信號通常需要使用超聲換能器進行接觸式探測。由于超聲波（光聲信號）無法在空氣中傳播，在樣品和超聲換能器之間必須添加水或超聲凝膠作為超聲波的耦合介質，這使其在很多應用中受到限制。比如，在高溫、強輻射等的工業應用中，往往因為待檢測區域所處的環境惡劣，無法使用耦合劑進行耦合，或者超聲換能器無法置于被檢測區域；在燒傷檢測和潰瘍檢測等的生物醫學應用中，使用耦合劑會污染乃至感染待檢測組織區域，而且在手術導航中耦合劑的使用會對目標區域的檢測帶來不便，這就限制了其在臨床醫學上的應用。光學超聲檢測是一種很有前途的超聲波檢測方法，以非接觸的方式檢測超聲波信號，可以實現寬頻帶高靈敏度的光聲成像，此方法避免了傳統方法需要耦合劑而帶來的不便。因此，與傳統的光聲成像相比，非接觸式光聲成像技術更利于擴展到工業及生物醫學領域的應用。

過去幾年，越來越多的光學超聲檢測方法被提出，光學檢測超聲信號的方法已經被廣泛的應用于各種領域[1-21]。早在上個世紀90年代，Pouet等[22]在光折變晶體中混合雙光束的方法形成外差干涉儀，并將其用來探測物體粗糙表面的超聲信號。2004年，Carp等[23]利用改進的馬赫曾德爾干涉儀測量光聲效應引起的表面振動的位移，該系統具備4 ns的時間分辨率，0.3 nm位移靈敏度。亦有研究使用基于聚合物薄膜的法布里—珀羅干涉儀來檢測超聲波[24]。2010年，Berer等[2]利用雙波長混頻干涉儀實現了遠距離和非接觸式光聲成像。最近，有研究提出將相位敏感型光學相干層析（Optical Coherence Tomography，OCT）系統用于光聲信號的檢測[25]。在以上提到的各種光學檢測超聲的方法中，有的由于檢測靈敏度低并不能進行非接觸式測量，有的由于光學信號調制頻率受限導致系統檢測帶寬較窄，有的則是成像分辨率較低。為了實現非接觸式寬帶高分辨率的接收光聲信號，我們在本文中提出了一種基于低相干邁克爾遜干涉儀的光學檢測光聲信號的方法，并建立非接觸式光聲成像系統，驗證了該系統的成像能力及橫向分辨率，并通過模擬實驗和在體實驗驗證了該非接觸式光聲成像系統能夠高分辨率的成像生物組織的脈管系統，促進了光聲成像技術在生物醫學中的應用。

1 原理與方法

1.1 光學相干法測量光聲信號的原理

當吸收體被短脈沖激光照射時，吸收體會吸收光能產生熱膨脹，從而輻射熱彈波，也就是光聲信號。光聲信號本質上是機械波，它在傳播的過程中會引起吸收體表面的振動。因此，我們可以通過測量吸收體表面振動位移的大小來反演光聲信號的大小。邁克爾遜干涉儀是常用的精密測量設備，其精度可達到皮米量級甚至更高，在科研及工業領域中都有著廣泛的應用。本文中，我們利用邁克爾遜干涉儀對相位變化的高度敏感性來檢測光聲信號引起吸收體表面振動的位移，從而實現光聲信號的非接觸式檢測。

光學相干法測量光聲信號的原理，見圖1。光聲信號檢測光經分束器分為兩束光，一束光通過透鏡并聚焦在反射鏡上形成參考臂，另一束光依次通過二向色鏡和物鏡最后聚焦在樣品上形成樣品臂，從樣品臂和參考臂各自反射回去的光相遇后發生干涉，光電平衡探測器用來接收干涉信號，這就構成了邁克爾遜干涉儀。光聲信號檢測光與光聲信號激發光經過二向色鏡后變為同軸光束，一同經過物鏡聚焦于生物組織上，生物組織吸收激發光能量后產生光聲信號，光聲信號傳到組織表面時會引起表面的振動，振動幅度的大小就反映了光聲信號的大小。根據光聲信號在介質中的傳播性質，組織的振動位移Δx(t)與引起組織振動的光聲信號p(t)的一維簡化關系如下[26]：

其中va為超聲在介質中的傳播速度，ρ為介質的密度，t為激光脈沖的脈寬。組織表面的振動會改變干涉儀中樣品臂的相位，因此通過檢測相位變化就能得到光聲信號引起的振動位移的大小。為方便，將Δx(t)寫成Δx。光聲信號引起的組織表面位移與相位的變化關系如下：

其中λ0為探測光的中心波長。通過干涉儀最終探測到的干涉光強為：

A為光電平衡探測器的放大倍數，I1、I2分別為參考臂及樣品臂反射回去的光強。相位φ包含環境噪聲（光學平臺的振動及空氣擾動等）產生的相位φnoise及光聲信號引起的相位變化Δφ，即：φ=Δφ+φnoise。由(3)式可得，當時，即干涉信號經過差分放大后輸出為“0”時，系統對光聲信號的檢測靈敏度最高，此時作為光聲信號的檢測時間窗口。

圖1 光學相干法檢測光聲信號的原理圖

1.2 非接觸式光聲成像系統

非接觸式光聲成像系統原理圖，見圖2。光聲信號檢測系統由超發光二極管（中心波長1310 nm，帶寬45 nm）、光纖環形器、2×2光纖耦合器、二向色鏡、反射鏡光電平衡探測器等組成（圖2a）。從光纖耦合器輸出的檢測光通過光纖準直器和透鏡1聚焦在反射鏡上形成干涉儀的參考臂，另一路檢測光通過二向色鏡、擴束系統（L2焦距是30 mm，L3焦距是100 mm）及物鏡聚焦在樣品表面形成干涉儀的樣品臂，調整樣品臂與參考臂的長度使其相等，從而使干涉儀工作在零差模式。考慮到該系統主要用來進行血管成像，所以選用的脈沖激光器波長是532 nm，脈寬為 10 ns，重復頻率是 0.11～10 kHz（HLX-I-F005，Horus Laser），4倍物鏡用來聚焦光束。由前文我們知道只有在環境噪聲? =kπ ±的情況下系統的靈敏度最高，在這里我

noise們使用FPGA控制卡采集由光電平衡探測器輸出的差分信號，對采集到的差分信號進行實時分析，當差分信號在零點附近時可編程邏輯器件（Field-Programmable Gate Array，FPGA）控制卡會發出脈沖信號T1去觸發激光器，激光器就會發出激光脈沖來激發光聲信號，并同時輸出同步信號T2取觸發NI-5124采集卡采集光聲信號，此時完成對一個掃描點信號的采集，然后由PCI-1716控制卡發出位移信號T3，驅動二維電機掃描平臺移動到下一個掃描點，依次循環，實現對樣品的二維掃描，最后對各個采樣點將得到數據進行最大值投影，從而得到一幅二維光聲圖像。為了降低噪聲影響，將信號通過截止頻率為30 kHz的高通濾波器。通常，生物組織對檢測光的反射率較低，為進一步提高系統的檢測靈敏度，通常在組織表面覆蓋一層幾百微米厚的油層用來增強反射光。系統的工作時序，見圖2b。

首先，我們測試了這種方法檢測光聲信號的可行性。為此，我們將黑色膠帶作為吸收體，使用非接觸式光聲成像系統測試了黑膠帶的光聲信號，結果見圖2c，表示的是光聲信號引起吸收體表面振動的位移大小。根據第一節中等式(1)，將圖2c中的位移信號進行微分，即可得到光聲壓信號。這個信號與我們用超聲換能器測到的光聲信號是一致的，具有雙極性。因此，這種光學方法檢測光聲信號是可行的。

圖2 非接觸式光聲成像系統原理圖

2 結果

2.1 系統分辨率測試

為了驗證該非接觸式光聲成像系統的橫向分辨能力，我們對直徑為6 μm碳纖維絲進行光聲成像實驗，通過對所獲得的光聲圖像進行分析來量化上述系統的分辨率。實驗中，電機的掃描步距為2 μm，掃描速度為40 Hz，且對每個掃描點的數據平均10次。圖3a是直徑為碳纖維的光聲圖像。取其在x軸掃描方向上的一維信號，然后通過高斯擬合擬合出光滑的曲線，最后取其半峰全寬作為系統的橫向分辨率數值，結果見圖3b，可以看到該成像系統的橫向分辨率優于10.6 μm。

圖3 非接觸式光聲成像系統的橫向分辨率

2.2 模擬樣品實驗

接下來我們證明了該非接觸式光聲成像系統在散射介質中也有較好的成像能力。實驗對象為埋在濃度為1%的瓊脂凝膠中的菩提葉脈，葉脈在實驗前經過染料處理，使得它在532 nm激發光下可以激發出光聲信號，凝膠模擬組織對光的高散射特性。實驗中，電機掃描步距為10 μm，成像范圍為4 mm×4 mm。成像結果，見圖4。葉脈的紋理能夠被清晰地成像出來，且圖像具有很高的對比度。模擬實驗證明了該非接觸式光聲成像系統能夠以高的分辨率成像高散射介質中的吸收體。

圖4 菩提葉脈的光聲成像實驗

2.3 在體實驗

為了驗證該非接觸光聲成像系統的在體成像能力，我們用該系統對接種有黑色素瘤的小鼠耳朵的血管及瘤體進行了活體成像。根據實驗規范，在實驗前我們先用麻醉劑隊小鼠進行腹腔注射（40 mg/kg）麻醉。為了保證小鼠耳朵不被損害，實驗過程中對激光能量進行測量，將其控制在12 mJ/cm2以下。圖5為黑色素瘤小鼠耳朵光聲二維成像結果，其血管網絡清晰可見，成像范圍為3 mm×3 mm，綠色圓形虛線框為黑色素瘤實體部位，由圖可以看出腫瘤相關血管結構具有高度的紊亂性，其新生血管結構為點狀結構，見圖5中的綠色矩形框的箭頭所示。此在體實驗結果說明了該非接觸式光聲成像技術可以活體成像生物組織，在識別皮膚腫瘤疾病等相關血管結構及腫瘤實體方面具有潛在的應用前景，并且與傳統方法比較，該非接觸式光聲成像技術對于病灶有潰瘍的區域的成像及檢測更具優勢。

圖5 黑色素瘤小鼠耳朵血管及瘤體的光聲成像實驗

3 結論

本文我們首先介紹了非接觸式光聲成像的背景，分析了基于全光纖低相干邁克爾遜干涉儀的光學檢測光聲信號的原理，然后詳細描述了非接觸式光聲信號檢測的方法，最后，成功構建了非接觸式光聲成像系統，測試了系統的分辨率，通過模擬及在體實驗驗證了系統的成像能力，實驗表明非接觸式光聲成像系統具高分辨成像生物組織的能力。另外，由于該光聲成像系統的非接觸的特性，因而更加方便于實際應用，使得該系統在生物醫學研究中具有廣闊的應用前景。同時，本文中光聲信號檢測所使用的低相干邁克爾遜干涉儀是OCT系統中的關鍵部件，因而可以方便的實現光聲/OCT雙模態成像系統，能夠提供更豐富的生物組織信息。

[1] Rousseau G,Gauthier B,Blouin A,et al.Non-contact biomedical photoacoustic and ultrasound imaging[J].J Biomed Opt,2012,17(6):061217.

[2] Berer T,Hochreiner A,Zamiri S,et al.Remote photoacoustic imaging on solid material using a two-wave mixing interferometer[J].Opt Let,2010,35(24):4151-4153.

[3] Hochreiner A,Berer T,Grün H,et al.Photoacoustic imaging using an adaptive interferometer with a photorefractive crystal[J].J Biophoton,2012,5(7):508-517.

[4] Monchalin JP.Optical detection of ultrasound at a distance using a confocal Fabry-Perot interferometer[J].Appl Phys Lett,1985,47(1):14-16.

[5] Paltauf G,Nuster R,Haltmeier M,et al.Photoacoustic tomography using a Mach-Zehnder interferometer as an acoustic line detector[J].Appl Optics,2007,46(16):3352.

[6] Nuster R,Holotta M,Kremser C,et al.Photoacoustic microtomography using optical interferometric detection[J].J Biomed Opt,2010,15(2):021307.

[7] Gratt S,Wurzinger G,Nuster R,et al.Free beam fabry-perotinterferometer as detector for photoacoustic tomography[J].SPIE,2013,8800(21):880002.

[8] Wild G,Hinckley S.Acousto-Ultrasonic Optical Fiber Sensors:Overview and State-of-the-Art[J].IEEE Sens J,2008,8(7):1184-1193.

[9] Grün H,Berer T,Burgholzer P,et al.Three-dimensional photoacoustic imaging using fiber-based line detectors[J].J Biomed Opt,2010,15(2):021306.

[10] Yang Q,Loock HP,Kozin I,et al.Fiber Bragg grating photoacoustic detector for liquid chromatography[J].Analyst,2008,133(11):1567.

[11] Yarai A.Fiber Bragg grating applied multi-functional sensor based on pulsed photoacoustic technique[J].Procedia Engin,2010,5:1180-1183.

[12] Draguioti E,Pérennès F,Beard PC,et al.Optical fiber photoacoustic–photothermal probe[J].Opt Lett,1998,23(15):1235-1237.

[13] Govindan V,Ashkenazi S.Bragg waveguide ultrasound detectors[J].IEEE Trans Ultrason Ferroelectr Freq Control,2012,59(10):2304-2311.

[14] Chao CY,Ashkenazi S,Huang SW,et al.High-frequency ultrasound sensors using polymer microring resonators[J].IEEE T Ultrason Ferr,2007,54(5):957-965.

[15] Paltauf G,Schmidtkloiber H,K?stli KP,et al.Optical method for two-dimensional ultrasonic detection[J].Appl Phys Lett,1999,75(8):1048-1050.

[16] Paltauf G,Schmidt-Kloiber H,Guss H.Light distribution measurements in absorbing materials by optical detection of laser nduced stress waves[J].Appl Phys Lett,1996,69(11):1526-1528.

[17] Wilkens V.Characterization of an optical multilayer hydrophone for use as broadband ultrasound reference receiver—comparison with PVDF membrane hydrophones[J].Proc IEEE Ultrason Symp,2002,1(1):773-776.

[18] Chow CM,Zhou Y,Guo Y,et al.Broadband optical ultrasound sensor with a unique open-cavity structure[J].J Biomedl Opt,2011,16(1):017001.

[19] Schilling A,Yava? O,Bischof J,et al.Absolute pressure measurements on a nanosecond time scale using surface plasmons[J].Appl Phy Lett,1996,69(27):4159-4161.

[20] Boneberg J,Briaudeau S,Demirplak Z,et al.Two-dimensional pressure measurements with nanosecond time resolution[J].Appl Phys A,1999,69(1):557-560.

[21] Yakovlev VV,Dickson W,Murphy A,et al.Ultrasensitive non-resonant detection of ultrasound with plasmonic metamaterials[J].Adv Mater,2013,25(16):2351.

[22] Pouet BF,Ing RK,Krishnaswamy S,et al.Heterodyne interferometer with two-wave mixing in photorefractive crystals for ultrasound detection on rough surfaces[J].Appl Phys Lett,1996,69(25):3782-3784.

[23] Carp SA,Venugopalan V.Optoacoustic imaging based on the interferometric measurement of surface displacement[J].J Biomed Opt,2007,12(6):064001.

[24] Zhang E,Laufer J,Beard P.Backward-mode multiwavelength photoacoustic scanner using a planar Fabry-Perot polymer film ultrasound sensor for high-resolution three-dimensional imaging of biological tissues[J].Appl Opt,2008,47(4):561-577.

[25] Li C,Wang RK,Wang Y.Noncontact photoacoustic imaging achieved by using a low-coherence interferometer as the acoustic detector[J].Opt Lett,2011,36(20):3975.

[26] Xu M,Wang L.Time-domain reconstruction for thermoacoustic tomography in a spherical geometry[J].IEEE T Med Imaging,2002,21(7):814-822.