高聲阻層下指紋的高頻超聲成像實驗

吳 迪，韋小鶴，畢 超，張博南，王亞平，滕永平

(1.北京交通大學 理學院 物理系,北京 100044；2.中國鐵道科學研究院集團有限公司 金屬及化學研究所,北京 100081；3.北京信泰智合科技發展公司，北京 100098)

指紋識別技術廣泛應用于生物識別技術領域，超聲波指紋成像傳感器的穿透力和分辨率優于光學和電容指紋傳感器[1]. 大學物理超聲掃描成像實驗中，為了獲取高精度的超聲指紋圖像，需要分析超聲波在多層介質中的傳播特性和超聲波指紋傳感器的聲學成像原理.

中科院聲學所李明軒等[2-3]以鋼層下多層橡膠層的層狀均勻粘結結構為檢測對象，做理論推導，指出多層介質超聲回波是各界面信號的疊加，通過對實驗波形的信號處理成功檢測出鋼層下4個界面的空氣脫粘缺陷. 在超聲成像方面，日本本田電子公司K. Kobayashi等[4]提出了生物組織二維聲阻抗成像方法，將超聲反射波頻譜強度轉化為目標組織的局部聲阻抗，與參照物(水)進行對比，實現了小鼠的小腦皮層的聲阻抗成像.

超聲成像中空間分辨率是指圖像中能夠分辨的2個目標點之間的最小距離，體現成像系統區分微小結構的能力，分為橫向分辨率和軸向分辨率，根據被測物體不同位置的回波信號和幅值差異，可以對被測樣品進行三維成像[5]. 韓國Hanyang大學W. Y. Choi等[6]利用激光激發超聲波，用聚焦超聲換能器接收，實現模擬指紋成像，圖像可達到88 μm軸向分辨率和120 μm橫向分辨率. 美國加州大學B. E. Boser等[7]研制出壓電微機械超聲換能器，以150 μm的軸向分辨率在4.6 mm×3.2 mm區域內的模擬指紋實現成像，圖像可達到75 μm橫向分辨率. 美國北卡州立大學C. Peng等[8-9]研制出壓電材料為PZT-5A/H的高頻壓電超聲換能器，并利用該換能器對1.0 mm×2.0 mm區域內的模擬指紋成像，獲得了軸向分辨率為80 μm、橫向分辨率為70 μm的指紋圖像. 對于超聲回波信號的處理可以采用基于希爾伯特變換的數字正交解調算法，將回波信號取包絡進行計算，然后應用軟件通過C 掃描方式顯示出超聲圖像[10].

本文針對由水層、玻璃層、模擬皮膚表皮層、模擬皮膚真皮層和水層組成的5層介質系統，利用傳遞矩陣推導出超聲波在5層介質中界面0處的反射系數，并以此作為系統的響應函數，用Matlab進行解析解的數值仿真，并與COMSOL有限元數值模擬結果進行比對. 基于大學物理實驗中的水浸超聲掃描成像裝置，用高頻超聲換能器對玻璃層下模擬指紋進行C 掃描實驗測量，由所得的時域波形數據可得到較高精度的指紋圖像. 實驗數據符合理論結果，驗證了高頻超聲波穿透高聲阻抗玻璃層的傳播規律可用阻抗傳遞矩陣法獲取，為大學物理實驗中超聲精細掃描指紋成像實驗提供聲學傳播的理論基礎.

1 理論模型

人體的皮膚由最外面的表皮層、中間的真皮層和深層皮下組織層3部分構成,其中手指皮膚表皮層厚度約0.2 mm，真皮層厚度約1.4 mm. 當玻璃屏下有指紋，模擬實驗時高頻超聲波換能器向模擬皮膚發射超聲波并接收反射的回波信號. 構建如圖1所示的多層介質模型，dn(n=0，1，2，3，4)表示各層介質的厚度.

圖1 模擬5層介質的二維示意圖

5層介質系統信號傳播如圖2所示.由上層至下層，設層0～4的縱波聲速分別為v0～v4，各層聲阻抗分別為Z0～Z4，各層的縱波波數分別為k0～k4,各層衰減系數分別為α0～α4.

圖2 5層介質系統的信號傳播示意圖

層0入射聲壓的表達式可寫為

pt0=Ai0eik0xe-iω t,

(1)

層0反射聲壓的表達式可寫為

pr0=Ar0e-ik0xe-iω t,

(2)

層1至層3的入射聲壓和反射聲壓表達式與式(1)和式(2)類似，系數分別為At1，Ar1，At2，Ar2，At3，Ar3.層4入射波的表達式如式(1)所示，系數為At4，由于層4設為半無限大介質，沒有反射波.

多層介質界面邊界條件為：界面兩側的聲壓和法向質點速度相等，由此可以得到各界面兩側的聲壓和法向質點速度的關系[11].根據此邊界條件，由阻抗傳遞法可得：

R23exp [i(2k1d1+2k2d2)]+R01R23R34exp [i(2k3d3)]+

R01R12R34exp [i(2k2d2+2k3d3)]+R12R23R34exp [i(2k1d1+2k3d3)]+

R34exp [i(2k1d1+2k2d2+2k3d3)]}/{1+R01R12exp [i(2k1d1)]+R12R23exp [i(2k2d2)]+

R01R23exp [i(2k1d1+2k2d2)]+R23R34exp [i(2k3d3)]+R12R34exp [i(2k2d2+2k3d3)]+

R01R12R23R34exp [i(2k1d1+2k3d3)]+R01R34exp [i(2k1d1+2k2d2+2k3d3)]}.

(3)

其中，Rhs表示第s層相對第h層的反射系數，kp為第p層介質(p=0,1,2,3)中的波數，波數表達式里包含有衰減系數αp.式(3)為界面0處聲波透過多層介質后的總反射系數，各層界面結構的反射系數表達式，可以用系統響應函數對各層回波信號進行定量分析.

圖3為時域信號經過穩定線性時不變系統的示意圖.圖3中x(t)表示輸入的時域信號，y(t)表示輸出的時域信號，h(t)表示穩定系統的響應函數.三者在時域上的關系可以表示為輸出信號等于輸入信號與系統響應函數的卷積，即

圖3 時域信號經過穩定線性時不變系統示意圖

y(t)=x(t)*h(t).

(4)

將圖1所示的分層介質模型視為穩定線性時不變系統，利用傳遞矩陣得到的界面0總反射系數式(3)作為系統響應函數H(ω)，將入射信號做傅里葉變換后在頻域上與系統響應函數H(ω)相乘，再做傅里葉逆變換，可獲得時域上的反射回波信號.

2 多層介質模型的仿真

1)利用Matlab軟件模擬高頻窄脈沖超聲波，并通過仿真得出聲波分別經過無指紋的多層介質系統和有指紋的多層介質系統的反射波波形圖.

1)利用COMSOL軟件的聲學模塊進行高頻超聲波在多層介質中傳播的有限元仿真.

3)將有限元仿真與Matlab解析法的數值仿真結果進行比對，得出超聲波在多層介質系統中的總聲壓場分布圖.

2.1 多層介質模型的Matlab數值仿真

高頻換能器激發超聲波透過玻璃層入射到組織層. 模擬計算反射回波時，將模擬的大平底面反射回波作為輸入信號x(t)，由式(4)可得到時域響應信號. 模擬計算時各層介質參量如表1所示，考慮到聲波在各層介質中會發生衰減，式(3)中對應exp [i(2kndn)](n=0,1,2,3)處都乘以衰減因子exp [-(2αndn)]. 基于Matlab軟件進行數值模擬的結果如圖4所示.

表1 各層介質的模擬參量

為了顯示指紋組織層對聲波的影響，將玻璃上沒有組織層和存在組織層2種情況進行對比顯示. 因為指紋存在凸凹不平的紋線(脊和谷)，凹紋線(谷)和玻璃接觸類似于紅色線回波情況，凸紋線(脊)和玻璃接觸類似藍色線回波情況. 從圖4淺藍色方框標識區域可以看出，指紋組織層的界面反射會引起玻璃底面界面回波的變化，這是因為玻璃層和指紋界面存在聲阻抗差異，導致聲波進入指紋組織層后產生能量減弱. 用Matlab軟件可以得到圖4淺藍色方框中2個回波各自的幅度峰值并且相除，入射到指紋組織層后反射的回波幅值下降至無指紋組織層時反射回波幅值的81.6%.

圖4 20 MHz高頻超聲波透過玻璃層后有模擬指紋組織和無指紋組織的反射回波信號

2.2 多層介質模型的COMSOL有限元仿真

為了驗證上述解析法數值仿真的可靠性，應用COMSOL軟件的聲學模塊對超聲波在多層介質中傳播過程做有限元仿真.

1)添加物理場，選擇壓力聲學-瞬態求解，根據圖1的模型，構建多層介質系統的二維模型，如圖5所示. 水層的寬度為4 mm，厚度為0.1 mm；玻璃層的寬度為4 mm，厚度為0.5 mm；模擬表皮層的寬度為2 mm，厚度為0.2 mm；模擬真皮層的寬度為2 mm，厚度為1.4 mm；完美匹配層模擬半無限大邊界，該處聲波被完全吸收，無反射回波；空氣層模擬無指紋情況，與有指紋情況進行對比.

圖5 5層介質系統COMSOL模型示意圖

2)設置各層介質的材料. 其中選擇材料庫中的皮膚作為表皮層材料. 由于材料庫中可以選擇的材料有限，將與皮膚的真皮層性質最相近的導熱硅脂作為真皮層材料. 設置完美匹配層材料為水.

3)將整個多層介質的求解設置為壓力聲學瞬態，并將邊界設置為平面波輻射. 在平面波輻射中，設置入射壓力場，大小為1 kPa.

4)設置各層介質的聲阻抗. 網格設置中，網格大小選擇極細化，各層介質均采用自由三角形網格形式，完美匹配層則采用網格層數為8的四邊形映射網格的形式. 網格設置越小，獲得的圖像質量越高，但是隨著網格的變小，計算量會變得越大.

在瞬態求解中，設置時間步長為T/60，即超聲波振蕩周期的1/60，以便獲得平滑的曲線圖. 經過有限元計算，超聲波在無指紋和有指紋的多層介質系統傳播過程中的總聲壓場分布如圖6所示.

圖6 超聲波在無指紋和有指紋多層介質系統傳播過程的總聲壓場分布圖

在水層與玻璃層的交界面處設置2個觀察點，如圖6中點1和點2所示，用以求得超聲波在2種不同情況下的總聲壓場的變化，2個觀察點的總聲壓場變化曲線如圖7所示.

由圖7可知，在空氣層中幾乎觀察不到超聲波的總聲壓場分布，表明在玻璃層與空氣層的交界處，入射的超聲波幾乎完全被反射. 而在表皮層和真皮層中可以觀察到微弱的超聲波總聲壓場分布，表明有部分超聲波會透過玻璃層與表皮層、表皮層與真皮層的交界面，于是在水層與空氣層交界面處觀察到的總聲壓場會有所減小. 圖7中，總聲壓場經過約1 μs達到穩定狀態后，在水層與玻璃層交界面，有指紋處的超聲波總聲壓場小于無指紋處的總聲壓場，且該處的總聲壓場約為無指紋處的72%(達到穩定狀態后綠色線聲壓值和藍色線聲壓值之比). 也就是說，假如換能器在空氣層與水層的交界面處接收超聲回波信號，在玻璃屏下無指紋與有指紋2種情況下，接收到的反射回波信號峰值會有差異，且后者峰值比前者小，這與Matlab解析法數值仿真結果吻合. 后續模擬指紋超聲掃描實驗的成像方法即是基于界面回波幅值變化的數值仿真結果.

圖7 無指紋和有指紋情況下水層與玻璃層交界面處總聲壓場變化

3 玻璃層下指紋的高頻超聲成像

為驗證上述多層介質理論模型的模擬結果，使用水浸超聲波掃描成像實驗儀檢測玻璃層下指紋的高頻超聲成像. 實驗裝置各模塊的功能框圖如圖8所示，實驗采用JSR DPR300脈沖發生器作為激勵裝置，高頻針式換能器的激發中心頻率為20 MHz，壓電晶片尺寸為0.8 mm×0.8 mm，橫向聲束分辨率約為100 μm. 數字示波器(RIGOL DS1074B)接收并且存儲信號，脈沖發生器產生的同步信號接入示波器的外同步輸入口. 電腦控制DPR300脈沖發生器產生激發高電壓和同步信號，同時控制水槽掃描電機以步長0.05 mm的步進精度運動. 在每個步進后的掃描位置，示波器通過以太網自動將波形以csv文件格式保存到控制電腦.

圖8 超聲實驗系統的各模塊

在測量過程中，高頻超聲探頭按如圖9所示的運動方式掃查. 超聲探頭在每個測量點處都會采集到1個超聲回波信號xq r(t)，q=1,2，…，m，r=1，2，…，n，其中m表示探頭沿x方向進行逐行掃查的點數，n表示沿著y方向進行逐列掃查的點數.

圖9 高頻超聲探頭運動掃查方式

掃查完畢后，利用希爾伯特變換求取每個掃描點回波信號的頻譜包絡峰值，可獲得信號矩陣P，將矩陣P中的每個幅值數據映射成彩色圖像值，即可得到C掃描成像結果，如圖10所示.

圖10 指紋掃描圖像處理效果圖

圖11紅色框內為在模擬指紋試樣實際掃描區域用光學顯微鏡觀察試樣后拍攝的圖像. 將圖10和圖11中數字(1～8)標注的圖像區域進行對比，圖10顏色偏紅的位置對應圖11試樣中存在空氣(即指紋中谷所在位置)，可以看出成像結果和實際掃描區域對應很好，掃描圖能夠顯示出模擬指紋的走向和形貌. 在實際高頻超聲掃描過程中，C掃描圖像雖然清晰且目標指紋形態完整，但仍存在空氣邊界不銳利以及對比度不高等問題[12]. 在大學物理實驗教學過程中，學生可使用不同參量進行掃描，使峰值測量更準確，保證較高成像分辨率，通過對比不同掃描成像結果，從而得到最佳的掃描參量.

圖11 顯微鏡下觀察的實際樣品圖

4 結 論

利用傳遞矩陣推導了超聲波在5層介質系統界面0處的反射系數,并以此作為系統響應函數. 利用Matlab數值模擬,得到了包含不同介質層信息的超聲回波波形，解析法給出了高聲阻抗玻璃層入射到指紋組織層中脊和谷的反射回波幅值存在反射波幅差異，并且與COMSOL有限元模擬結果進行比對，結果一致，因此可用于模擬指紋試樣的C掃描成像. 由實驗結果可見，超聲精細掃描成像技術應用于大學物理教學實驗具有可行性. 利用理論模型數值求解結果，學生能夠較好地理解超聲波透過高聲阻抗介質層在指紋組織層中傳播的全過程. 通過實驗學生可以直接觀察到超聲C 掃描圖像的掃描過程，完成數據采集和數據處理，可以實時分析圖像的處理結果. 這一過程耗時較短，實驗簡單，容易操作，可進行多組對比實驗. 同時通過超聲成像拓展實驗，學生可對時下熱門的屏下指紋解鎖聲學原理及超聲精細掃描成像技術有更直觀、更清晰的認識，實驗教學效果良好.