|作者：臧雨宸 林偉軍 蘇 暢
(1 中國科學院聲學研究所)
(2 中國科學院大學）
本文選自《實體》2021年第11期
摘要聲波是在彈性媒體中傳播的經典機械波。當聲波在物體表面發生反射、折射、散射等效應時，會與物體産生動量與能量的交換。這種交換在宏觀上表現為聲波對物體施加力的作用，稱為聲輻射力，瑞利和朗之萬兩位科學家正是聲輻射力的最早發現者。經過一百多年的探索，聲輻射力的基本理論和研究方法已經得到極大的豐富和發展，成為聲學的前沿與熱點問題。與此同時，在生物醫學、材料科學、納米科學等領域，基于聲輻射力的無損精準操控技術正受到越來越多的關注，展現出巨大的應用前景。
關鍵詞聲輻射力，聲懸浮，聲操控，聲輻射力彈性成像，聲輻射力天平
1什麼是聲輻射力
19世紀中葉，電磁學的集大成者麥克斯韋在完善電磁場理論後預言：光是一種電磁波，當光入射到物體表面時，會對物體施加壓力，稱為光壓。1901年，俄國實體學家列别捷夫首次用實驗測得了光壓的存在，有力地證明了電磁場的物質性。在此基礎上，科學家們利用光壓理論成功解釋了彗星出現彗尾的現象。和電磁波不同，聲波是在彈性媒體中傳播的經典機械波，反映了由于聲源振動而引起的彈性媒體周期性的膨脹與壓縮。與其他波動一樣，聲波也能攜帶動量與能量。當聲波在傳播過程中遇到障礙物時，會在物體表面産生反射、折射、散射等實體效應，并與物體發生動量與能量的交換，這種交換在宏觀上表現為聲波對物體的力的作用，稱為聲輻射力。
根據聲場的動量守恒定律可以推導得到聲輻射力的一般表達式為
其中，d = d，為物體表面的法向機關矢量，為物體表面，稱為平均聲輻射應力張量，其具體表達式為
其中，和分别表示聲壓和質點速度，表示二階機關張量，和分别是流體的密度與縱波聲速，表示對實體量求時間平均。正因為聲輻射力是一個時間平均量，故而又被稱為穩态聲輻射力。
從(2)式不難看出，聲輻射力由兩部分組成。其一是聲場平方項的貢獻，對應着(2)式的第一項。線上性聲學的範圍内，聲壓在一個周期内正負值恰好完全相抵。然而我們知道，描述聲波的三個基本方程——運動方程、連續性方程和媒體的本構方程都是非線性的，如果将這些非線性項考慮進來，聲壓的正負分量便不會完全抵消，在時間平均後會出現不為零的“直流”分量。圖1形象地顯示了随着傳播距離的增加，空氣中聲壓級為140 dB、頻率為1000 Hz的聲波産生的非線性波形畸變，從上至下分别對應着馬赫數(質點振速與聲速之比)為0.002、0.003和0.004的情況。從圖中可以看出，随着馬赫數的增加，聲波的間斷距離明顯減小，即非線性畸變更加顯著。其二是流體動量流的貢獻，對應着(2)式的第二項，這一項與聲場的非線性無關。值得一提的是，電磁場是由嚴格線性的麥克斯韋方程組描述，其輻射力僅僅源于電磁場動量流的存在，與非線性無關。從這點來看，聲輻射力的産生機理要更加複雜一些。
圖1 空氣中聲壓級為 140 dB、頻率為 1000 Hz 的聲波随着傳播距離的增加發生了非線性畸變
聲輻射力是聲場的二階量，因而一般情況下産生的聲輻射力很小。例如對130dB(相當于火箭發射點附近的聲強)的聲波，其産生的聲輻射壓還不到0.1Pa，但當聲壓級增加到174 dB時，聲輻射壓可以達到1000 Pa。由此看來，聲輻射力主要在高聲強作用下較為顯著。
2聲輻射力的研究脈絡與方法
早在20世紀初，聲學界的兩位集大成者瑞利和朗之萬就分别提出了聲輻射力的概念。瑞利對聲輻射力的描述可以總結為：在平面波傳播的情況下，随流體質點一起運動的表面上受到一個不為零的平均壓力；假設該流體具有相同的平均密度但處于靜止狀态，後者也将承受一個平均壓力，前者與後者之差定義為瑞利聲輻射力。可以看出，瑞利聲輻射力是針對無限大流體媒體中的一維平面行波而言的，實際狀态下，這樣的條件過于苛刻，很難得到滿足。但從線性聲學的角度出發，在一個理想剛性壁的管中，如果管子的直徑小于四分之一波長，則其中隻有平面波傳播，進而為瑞利聲輻射力提供了一個可以實作的實體圖像。後來，朗之萬取消了這一限制，重新定義了可以适用于有限寬波束的聲輻射力：聲場中随流體質點一起運動的物體表面上受到的時間平均力與該物體後面未被擾動流體中的壓力之差稱為朗之萬聲輻射力。
然而不可否認的是，上述兩位科學家提出的聲輻射力理論都比較抽象，與實際的聲操控還相去甚遠。1934年，加拿大實體學家King首次利用嚴格的散射理論計算了理想流體中剛性球受到的聲輻射力，為之後的聲輻射力研究打下了重要的理論基礎。後來，日大學學家Awatani、Yoisoka、Hasegawa等人用同樣的方法研究了彈性球、彈性柱等常見模型在行波和駐波場中的聲輻射力，給出了若幹富有價值的理論結果。進入20世紀90年代，随着超聲換能器技術的不斷發展，使用各種新型聲波進行微粒的聲操控與聲捕獲已經成為可能。于是，以高斯波束為代表的聚焦超聲下的粒子聲輻射力成為一時的研究熱點。1990年，佛蒙特大學的吳君汝教授計算得到了聚焦波束作用下球形和柱形粒子的聲輻射力，并進行了實驗驗證。21世紀以來，Marston、Mitri等人将聲輻射力的研究拓展到貝塞爾波束領域，以期充分利用貝塞爾波束的非衍射特性。其中，負向聲輻射力的發現尤為值得關注，它從理論上論證了利用單聲束制成的聲鑷子進行粒子捕獲的可能。随後，Silva等人将聲輻射力的計算推廣到任意波束。近年來，随着聲操控技術的不斷發展，複雜模型和複雜邊界下的聲輻射力研究成為了新的熱點，本文作者所在團隊在此方面也有一些進展。
綜合來看，目前聲輻射力的研究方法大緻可以分為散射法和數值法。其中，散射法基于嚴格的波動理論，可以得到有利于參數化分析的解析解，但當模型較為複雜時計算比較繁複，是以散射法主要用于高斯聲束、貝塞爾聲束等作用下圓柱、球等規則物體的聲輻射力計算。利用散射法可以計算得到平面行波場對剛性球的聲輻射力随無量綱頻率的變化關系，是聲波的波數，剛性球的半徑，結果表明，剛性球的聲輻射力随着的增加先增大而後趨于穩定。當聲場比較複雜、物體形狀不規則時，利用散射法推導聲輻射力的解析解往往是很困難的，此時可以考慮借助純數值的方法來計算，如有限差分法、有限元方法、邊界元方法等，其中應用最廣泛的是有限差分法。該方法直接從聲波的三個基本方程出發，不需要任何形式的導出方程，對所研究對象的形狀和尺寸沒有任何限制，但計算所消耗的資源較大。2010年，中國科學院深圳先進技術研究院蔡飛燕等人利用有限差分法計算了彈性柱的聲輻射力特性，其結果與解析解符合得很好。2015年，陝西師範大學孫秀娜等人運用該方法分析了離軸情況下剛性球的聲輻射力，結果表明剛性球在高斯聲場中的聲輻射力随着粒子離軸距離的增大而減小。
此外，在高頻和低頻情況下，還可以分别利用聲線法和梯度法進行簡化處理。高頻條件下，物體的尺寸遠大于聲波波長，此時聲波的波動性可以忽略，隻需考慮其粒子性即可。是以，我們可以利用聲線的疏密來表征聲能量的強弱，進而用射線理論來計算此時的聲輻射力。2005年，Lee等人最早利用聲線理論研究了流體媒體中任意位置粒子的聲輻射力特性，從理論上證明了聲捕獲的可能。低頻條件下，物體的尺寸遠小于聲波波長，此時物體對聲場的影響可以忽略，可以認為聲能密度在空間是連續分布的，因而聲輻射力可以近似表示為某種勢函數的空間梯度，如同電場力可以表示為電勢能的空間梯度一樣。前蘇聯科學家Gor′kov最早計算得到了低頻情況下聲輻射力所對應的時間平均勢函數，進而大大簡化了聲輻射力的運算。
3對抗重力的聲懸浮
從古至今，人們一直渴望能夠擺脫地面重力的束縛，像小鳥一樣在空中自由翺翔。明代的萬戶是世界上第一個想到借助火箭推力升空并付諸實踐的人，但卻為此付出了生命的代價。近代的工業革命促進了熱氣球、飛機等飛行器的相繼問世，無疑讓人們實作了亘古的飛天夢想。然而，人們與重力的對抗才剛剛揭開序幕。1901年，英國科幻小說家赫伯特·喬治·威爾斯在《月球上最早的人類》一書中最早提出了反重力裝置的概念，他設想有一種能夠屏蔽引力的金屬叫“卡弗質”，隻要把它放在飛船下面，就可以屏蔽地球引力。他制造了一艘遊艇，在下面屏蔽了地球引力，上面向月球打開引力，遊艇就把他送上了月球。遺憾的是，迄今為止這種夢幻式的宇宙航行還沒有實作，人們還是在地球上費盡九牛二虎之力設法對抗地球引力，通過加大能源的消耗來實作“反重力”，讓自己飛向茫茫太空。
“卡弗質”的出現也許還在很遙遠的未來，但擺脫重力的懸浮技術已經變成觸手可及的現實。在航天工業和材料科學等領域，無容器環境對研究凝聚态實體是非常重要的實驗條件。目前地面上模拟空間無容器環境的方法主要有自由落體方法和懸浮方法。與自由落體方法相比，懸浮方法可以獲得持續的微重力和無容器狀态，因而應用更加廣泛。常見的懸浮方法主要有氣動懸浮、電磁懸浮和聲懸浮。氣動懸浮的懸浮力來源于氣體掠過表面時的動量減少，其穩定性較差。電磁懸浮的穩定性很好，但隻适用于導電體與磁體。而聲懸浮技術對材料的實體限制較少，且較為穩定。
圖2 聲波使Kundt管中的塵埃微粒有序懸浮跳動
有趣的是，曆史上的聲懸浮完全誕生于一次意外的偶然。1866年，年僅27歲的德國實體學家Kundt在進行固體和氣體的聲速測量實驗時意外發現：諧振管中的聲波能夠讓塵埃顆粒有序地懸浮和舞動，這些灰塵微粒有序的懸浮跳動正是聲懸浮。圖2是Kundt當時所用諧振管的示意圖，管内兩列振幅相同、傳播方向相反的行波疊加形成駐波場，波節處聲壓為零，微粒始終保持靜止，其餘位置處的微粒将來回跳動，其中波腹處的聲壓最強，微粒跳動最劇烈。Kundt也是以成為發現聲懸浮現象的第一人，但當時人們并不十厘清楚這一現象的産生機理。直到1934年King給出理想流體中剛性球聲輻射力的計算結果，揭示了聲懸浮來源于高聲強作用下的非線性現象。1964年，美國明尼蘇達州立大學的Hanson等人基于King的理論制造了第一台用于研究單個液滴動力學行為的聲懸浮器。1975年，美國科學家Whymark利用聲懸浮實驗研究了鋁、玻璃等聚合物在無容器條件下的熔化和凝固過程。1992年，科學家在航天飛機上利用超聲駐波場産生輻射壓力使材料浮起不接觸容器壁溶解，合成了用于制造非可視域低損耗光纖的高純度玻璃材料。經曆了一個多世紀的發展，聲懸浮從一項偶然的發現逐漸發展成為了一種成熟的技術手段，為各種材料的無容器處理和液滴動力學等領域的研究開拓了廣闊的新天地。
圖3 小球穩定懸浮在駐波場的波節下方
聲懸浮技術中最關鍵的兩點是懸浮位置和懸浮能力，我們以最簡單也是最常見的單軸式聲懸浮系統為例作一簡單說明。在穩定的一維駐波場中，聲壓的波節位置對應着時間平均勢函數的極小值點，也是聲輻射力取零值的位置。當不考慮微小粒子的重力時，粒子将穩定懸浮于駐波場的波節處。在考慮粒子重力的情況下，必須借助一定的聲輻射力來克服自身的重力，此時的穩定平衡位置将會稍稍偏離駐波場的波節，如圖3所示。有必要指出，聲輻射力是一個瞬時變化的實體量，這裡的聲輻射力取零是時間平均的效果。由于聲波一個周期的時長很短，我們肉眼根本觀察不到小球瞬時的快速振動，而隻能看到其運動的平均效應，進而認為其在聲場中保持穩定的平衡。懸浮能力主要展現為該系統可以懸浮物體的最大密度。為了提升懸浮能力，往往需要增加發射換能器的聲壓振幅，并對發射端和反射端進行相應的優化設計。此外，當諧振模式較低，聲壓波節較少時，懸浮能力也會得到一定的增強。2002年，西北工業大學解文軍等人成功實作了對地面上密度最大的固體銥和密度最大的液體汞的穩定聲懸浮(圖4)，進而證明了原則上可以對地面上的任意物質進行聲懸浮。該課題組還利用聲懸浮裝置進行了液滴物性參數的非接觸反演測量研究，包括密度、粘性系數、表面張力和比熱等，并對懸浮液滴的動力學行為進行了理論計算與實驗觀測。
圖4 地面上密度最大的固體銥(a)和密度最大的液體汞(b)的聲懸浮
4聲鑷子和聲馬達
4.1 聲鑷子
日常生活中，我們經常會用鑷子來抓取那些肉眼可見但用手卻抓不住的小物體，但對于細胞、分子等微小粒子的操控，普通的鑷子就無能為力了。20世紀70年代，美國科學家阿什金(A. Ashkin)最早提出可以利用光的輻射壓力對微小物體施加力的作用。1986年，他成功将這一設想變成了現實，制成了世界上第一台用雷射來實作粒子操控的裝置，稱為光鑷子。現如今，光鑷子已經被廣泛運用到病毒、細菌、細胞等微粒的無創操控中，大大促進了生物技術的發展，阿什金教授也是以榮獲2018年諾貝爾實體學獎。
盡管如此，光鑷子仍然具有自身的局限性。首先，光鑷子對樣本的透明性要求很高；其次，聚焦雷射的高能量往往會引起粒子的局部升溫進而造成對樣品的損傷。在阿什金制成光鑷子僅僅5年後，吳君汝教授便成功進行了聲鑷子的設想與實驗，他采用兩個平行對立放置的3.5 MHz聚焦超聲換能器形成駐波場，實作了對直徑為270μm的乳膠粒子和蛙卵的捕獲和操控，首次證明了聲鑷子的可行性。誠然，聲波的能量遠不能與雷射相比拟，但聲波可以在固體、流體等任意媒體中傳播，不受媒體電磁特性和透明性的影響。另一方面，其能量與工作頻率和醫學超聲成像系統的參數相當，可以實作對單個細胞或顆粒的操控，并確定生物組織和目标粒子有足夠的安全性。
圖5 利用聲表面波實作對微粒的排列和操控，其中(a)，(b)是放置在晶片上的聲鑷系統，(c)，(d)和(e)，(f)分别是駐波場作用下對微粒的一維和二維操控
駐波場具有較強的空間能量梯度，因而更容易産生明顯的聲輻射力來實作對特定微粒的捕獲。當微粒落入這樣的聲場中，其在聲輻射力的作用下将被推到并鎖定在時間平均勢函數的極小值位置，如同被鑷子牢牢夾住一樣。在此基礎上，通過調節聲場改變波腹和波節的分布，進而将微粒移動到理想的位置。從操控的次元來看，聲鑷子既可以實作一維和二維粒子的空間排列，也可以實作粒子在三維空間的移動變換。此外，由于聲波波長的尺度跨度很大，原則上聲鑷子可以操控從納米到厘米尺寸的所有微粒，甚至還可以操控流體媒體，形成特定的流場環境。2009年，美國賓州州立大學的Shi等人基于聲表面波的原理，利用兩個相對放置的叉指換能器産生駐波聲場，實作了對微小粒子的一維和二維排列與操控，如圖5所示。随後，該團隊又利用叉指換能器和單層微流控通道實作了晶片上液滴的分類和排序，在晶片上實作了液滴微流控技術。
圖6 利用高頻超聲實作對細胞的捕獲。在換能器激發聲場之前，細胞位于距離勢阱中心 50 μm 處(a)，當激發聲場後，細胞受到輻射力的牽引作用向勢阱中心移動(b)，(c)，最終在勢阱中心處再次保持靜止(d)
駐波聲場至少需要一組換能器，或一個換能器與一個反射面，在實際的聲操控中這樣的條件往往很難得到滿足。與駐波場聲操控相比，單波束聲鑷子的裝置更加簡單，受到了研究人員的廣泛青睐。2009年，南加州大學的Lee等人用正弦信号驅動中心頻率為30 MHz的铌酸锂壓電換能器來捕獲軸線上的微粒，發現當水中的球形粒子偏離軸線時，聲輻射力會将其拉回到軸線上。這一回複力現象為單波束聲鑷子提供了理論依據。随後，他們在聲線理論的基礎上，使用30 MHz的高頻聚焦超聲成功捕獲了直徑為126 μm的油酸微粒，又利用100 MHz的高頻超聲成功實作了對細胞量級粒子的聲捕獲，如圖6所示。2016年，Ma等人基于微粒在行波場中的共振特性，在聲表面波微流控腔體内實作了對粒子的精确篩選。2017年，Chen等人通過改變激勵頻率制成“變尺度”聲鑷子，進而可以操控不同尺度大小的微粒。2018年，Zhao等人在PZT圓片上利用不同扇區的聲波幹涉效應産生近似的貝塞爾聲束，通過其負向聲輻射力形成的勢阱捕獲了亞毫米尺度的微粒。同年，英國布裡斯托爾大學Drinkwater教授和西班牙納瓦拉公共大學Marzo博士在美國國家科學院院刊(PNAS)上首次展示了多粒子的獨立懸浮與操控，該裝置被稱為全息聲鑷子(HAT)。HAT包含一個由256個直徑均為1cm的揚聲器組成的陣列，陣列單元之間有約23 cm的間隔，并且工作在40 kHz的超聲頻段。通過調節聲波的實體特征，成功實作了二維和三維空間的聲波牽引技術，使25顆直徑在1—3 mm的聚苯乙烯泡沫塑膠小球執行了一系列令人驚歎的空中動作(圖7)。Marzo博士在接受西班牙科學新聞周刊的采訪時稱，他希望未來能通過全息聲鑷子在微觀尺度上操縱3D細胞并搭建出二維培養皿中無法存在的立體生物學結構。倘若這一設想變為現實，無疑将成為結構生物學領域的又一次偉大革命。
圖7 利用全息聲鑷子實作三維空間的聲波牽引
總體來看，與技術相對成熟的光鑷子相比，目前的聲鑷子在靈敏度方面還存在一定的劣勢，但我們有理由相信，随着科學家們從原理、裝置和應用等方面對聲鑷子技術進行改良與拓展，聲鑷子必将朝着更高精度、更成熟、實用性更強的方向發展。
4.2 聲馬達
聲波不僅可以攜帶動量，還可以攜帶角動量，我們将攜帶角動量的聲場稱為渦旋聲場。渦旋聲場在傳播過程中會産生波陣面繞軸的扭轉，進而使其傳播相位關于聲束中心呈現螺旋變化的關系。定義拓撲荷數或聲場階數為一個波長的傳播距離内波陣面發生扭轉的次數，圖8顯示了拓撲荷數取不同值時聲場的波陣面形态圖。顯然，當=0時，渦旋聲場将退化為普通的平面波場。此外，拓撲荷數的正負表示波陣面的扭轉方向。當拓撲荷數非零時，我們稱聲場攜帶軌道角動量(OAM)資訊。與電磁波不同，流體中的聲波是不存在偏振現象的，因而不存在自旋角動量(SAM)，其軌道角動量就是最終的總角動量。
圖8 拓撲荷數取不同值時聲場的波陣面形态
當聲波與物體之間發生動量的傳遞時，宏觀上表現為聲輻射力的作用。我們自然地想到，聲波能否與物體發生角動量的傳遞，進而表現為力矩的作用？答案是肯定的，這種力矩稱為聲輻射力矩或聲輻射轉矩，最早由美國學者Maidanik于1958年提出。與聲輻射力有所不同，聲輻射力矩的産生必須具備兩個條件：其一，聲場是渦旋聲場；其二，物體能夠産生一定的聲吸收。華盛頓州立大學的張黎坤教授和Marston教授系統總結了聲輻射力矩的理論架構，并從理論上證明了對在軸對稱目标而言，聲輻射力矩正比于聲吸收功率、聲場的拓撲荷數，反比于聲波的角頻率。2012年，鄧迪大學的Demore等人利用1000個陣元的換能器陣列産生不同階的貝塞爾渦旋聲場來操控水中的吸聲圓片，首次從實驗上驗證了聲輻射力矩和聲吸收功率的正比關系。同年，Silva和Mitri等人成功推導出了任意渦旋聲束對在軸球形粒子聲輻射力矩的表達式，并以貝塞爾波這一最常見的渦旋聲束為例進行了仿真。近年來，張黎坤教授課題組又詳細分析了貝塞爾聲束對離軸物體的聲輻射力矩特性和力矩反轉的機理。本文作者所在團隊首次在三維聲輻射力矩的問題中引入阻抗邊界，分析了渦旋聲場對邊界附近粘彈性球殼的力矩。
圖9 在銅闆上刻蝕的阿基米德螺旋線，用其實作了聲渦旋
利用聲輻射力來操控微粒運動的技術稱為聲鑷子，利用聲輻射力矩來驅動微粒轉動的技術則被形象地稱為聲馬達。聲馬達的出現擴充了現實中聲操控的自由度，大大豐富了微粒運動的類型與方式。聲馬達技術的關鍵是渦旋聲場的産生，目前來看主要有三種方式。第一種方法是通過表面螺旋結構的聲源實作聲渦旋。2016年，武漢大學劉正猷教授課題組在銅闆上刻蝕了阿基米德螺旋線，成功獲得了貝塞爾渦旋聲場(圖9)，并驅動水中的聚苯乙烯顆粒呈環狀分布。2020年，該課題組基于同樣的原理利用矽基片制作聲渦旋透鏡，成功操控了蝦卵在聲場中的轉動，其觀測到的動力學特性與理論符合得很好。該方法具有技術簡單和成本低的優勢，但存在難以調節聲場參數的困難，且在低頻情況下龐大的結構尺寸限制了其在實際中的應用。第二種方法是通過有源相控陣技術實作聲渦旋。2016年，西北工業大學洪振宇課題組利用64陣元換能器産生二維貝塞爾渦旋聲場實作了對聚苯乙烯顆粒和水組成的兩相混合物的操控(圖10)。當顆粒尺寸較小時，會随流體一起旋轉進而實作聲渦旋的可視化。相控陣技術可以靈活地調節聲場和改變參數，但對每個換能器的電路均需要獨立控制，成本較為昂貴。第三種方法是通過聲人工結構和超材料實作聲渦旋。2016年，南京大學程建春教授課題組利用聲學共振原理，通過在厚度和半徑都遠小于聲波波長的平面狀共振超表面中有效地調控由聲學雜化共振所引起的等效聲波數分布，将入射的平面聲波轉換為攜帶軌道角動量的渦旋聲場(圖11)。該方法基于聲學超構材料的特殊聲學特性和亞波長尺度的幾何優勢，有利于産生更穩定、寬頻帶且可靈活調控的渦旋聲場。
圖10 利用64陣元換能器陣列産生貝塞爾渦旋聲場，實作了微流控的精确操控
圖11 利用聲學共振結構産生聲渦旋
基于聲渦旋的聲馬達技術開辟了新的操控自由度，因而已經逐漸成為聲操控領域的熱點問題。順便指出，聲渦旋除了具備在聲操控、聲懸浮等領域向物體傳遞角動量的力學特性外，在聲學通信領域也具有巨大的應用潛力。攜帶不同OAM的渦旋聲場之間彼此正交，恰好對應着不會發生混疊的不同信道。倘若OAM能夠與傳統的多路複用傳輸系統相容，則有望實作水下通信技術領域的巨大變革。
5聲輻射力彈性成像
人們很早就發現了生物組織的彈性資訊與其病理狀況緊密相關。漢代醫聖張仲景的《傷寒雜病論》中就對胸腹部的按診進行了詳細的論述，現已成為診斷和治療疾病的重要依據。盡管如此，單純依靠手指的觸摸去定性地感覺組織的彈性無疑對醫生的經驗有較高的依賴，甚至不同的醫生會得出截然不同的診斷結果。此外，對于體内較深部位的器官，人工按診也無法準确獲知其彈性和硬度的資訊。為了彌補傳統按診法的缺陷，人們緻力于發展多種對組織的成像方法，希望能夠精确有效地反演出組織的彈性資訊，在這一需求背景下，超聲彈性成像方法應運而生。
根據激勵源的不同，我們可以将目前在臨床上廣泛使用的超聲彈性成像方法分為兩大類。第一種是準靜态壓縮彈性成像，最早由美國教授Ophir和日本學者Yamakoshi于20世紀90年代提出。該方法通過在組織外表面施加一定的壓力使人體組織發生形變，再利用B超波束對其進行掃描監測産生組織的應變圖。這種方法的局限性在于它仍然隻能探測到淺表的組織資訊并且容易受到邊界條件的影響。此外，它無法給出組織硬度的具體數值，隻能反映和背景組織相對比的一個相對大小。第二種是基于聲輻射力的動态彈性成像技術，該技術使用超音波對組織施壓，很好地克服了傳統準靜态壓縮彈性成像中無法對深部組織成像的問題以及易受邊界條件影響的不足，受到了國内外的廣泛關注。聲輻射力彈性成像主要可以細分為諧波運動成像(HMI)、剪切波彈性成像(SWEI)和聲輻射力脈沖成像(ARFI)等手段，這裡分别對它們作簡要介紹。
圖12 慢性肝病患者的肝髒剪切波彈性成像及定量分析圖。下圖是原始的成像圖，上圖的藍色區域則是在成像區域内選取直徑為3 mm的圓形定量檢測區域進行彈性模量的檢測
諧波運動成像最早由Maleke等人提出。通過探頭産生的超音波對組織施加周期性的聲輻射力，使組織産生諧波運動。通過擷取組織被激發前後的B超圖像并對其進行運動估計，可以間接反映組織的彈性狀況。該方法可以準确反映深部組織的局部彈性狀況，但由于需要在激發信号的同時監測組織的運動，因而采集到的B超信号容易受到激發信号的嚴重幹擾，這是該成像方法迄今為止尚未得到完美解決的主要問題。當聲輻射力作用在組織上時會激發出剪切波，剪切波彈性成像正是通過檢測這一剪切波的傳播來達到成像的目的，最早由美國學者Sarvazyan和Emelianov提出。他們利用高強度聚焦超聲探頭産生聲輻射力，在人體組織中産生了明顯的剪切波，并運用磁共振成像方法監測該剪切波的傳播情況。由于剪切波速與人體組織的粘彈特性緊密相關，是以通過檢測剪切波的傳播速度等參數可以擷取組織的彈性資訊。圖12顯示了某一慢性肝病患者的肝髒彈性成像圖。聲輻射力脈沖成像技術最早由杜克大學的學者Nightingale等人提出，該方法通過高能聚焦超聲探頭在局部組織中産生短暫的高強度聲輻射力，進而推動組織産生一定的應變，再通過追蹤這些微小應變和剪切波的傳播來判斷組織的粘彈特性。研究人員發現，組織對于脈沖聲輻射力的瞬态響應直接相關于局部區域的硬度大小，而通過求解亥姆霍茲方程的反問題方法可以量化組織中剪切波的波速。Nightingale等通過詳細的數值仿真和體模研究證明了聲輻射力脈沖成像的可行性，并得到了與傳統B超圖像相比更高的對比度和分辨力。如今的聲輻射力脈沖成像已在臨床上具備多種髒器檢查的能力，圖13顯示了某位良性纖維瘤患者的乳腺聲輻射力脈沖成像圖。
圖13 乳腺良性纖維瘤患者的聲輻射力脈沖成像圖
基于聲輻射力的彈性成像技術是對傳統準靜态彈性成像技術的重大突破，使得彈性成像對較深部位或有局部遮擋的組織器官也有了一定的用武之地，在離體成像實驗和臨床診斷治療中均顯示了良好的應用前景。當然，如何提升成像速度、改善成像品質、確定劑量安全性也是亟待解決的重要問題。随着技術的不斷發展，相信聲輻射力彈性成像技術将會在不久的将來更好地造福人類。
6聲輻射力天平
聲學中，我們将聲源在機關時間内向外輻射的聲能量定義為聲功率。超聲功率的測量是聲學研究中的熱點問題，對評估超聲裝置的安全風險、控制超聲裝置的輻射劑量、確定超聲診斷和治療過程的可控性與安全性有着重要意義。目前，關于聲功率的測量方法主要有聲輻射力法、聲光法、量熱法等，其中聲輻射力法是國際電工委員會(IEC)推薦的首選測量方法，也是目前唯一國際公認的标準方法。早在1976年，波蘭科學院的Zieniuk和英國薩裡大學的Chivers就對聲輻射力法和量熱法進行了詳細的比較分析，認為前者的适用範圍更廣、測量更加簡潔有效。
圖14 利用全吸收靶(左)和全反射靶(右)測量聲功率的實驗裝置示意圖
利用聲輻射力測量聲功率的裝置被稱為聲輻射力天平。對于小振幅聲波而言，在橫向尺寸不受限制的行波傳播方向上，兩種媒體的分界面或聲場不均勻處存在單向的穩态聲輻射力。事實上，這一壓力正是前文所述的朗之萬聲輻射力，機關面積上的聲輻射力稱為聲輻射壓。對理想的平面行波而言，聲輻射壓和聲強、波速存在關系：
其中，是聲波的能量密度。在測量中，将全反射靶或吸收靶放置于水中聲場裡感覺聲能流，并用天平來精确測量此時靶受到的聲輻射力。最後，将聲輻射力換算成聲功率即得到最終的測量結果。圖14給出了利用全吸收靶和全反射靶測量聲功率的實驗裝置示意圖。全反射靶和全吸收靶的主要差別在于靶的結構特性不同。全反射靶通常是用薄金屬皮制成的中空凸圓錐體，其錐角為90°，它能将接收到的超音波能量幾乎完全反射回去。全吸收靶的表面則是由許多鋸齒形吸聲尖劈排列而成，進而能夠吸收絕大部分的聲能量。圖15顯示了兩種靶的實物照片。對于全吸收靶測量而言，聲功率和輻射天平的讀數滿足關系= ，其中是測量當地的重力加速度；對于全反射靶測量而言，上述關系變為= /(2cos)，其中是聲束在圓錐形全反射靶上的入射角。當然，在聲吸收不可忽略時，我們還應當計及聲波在水中的吸收聲衰減。
圖15 全吸收靶(左)和全反射靶(右)的實物照片
2012年，上海交通大學壽文德等人提出利用聲輻射力天平測量電聲器件參數的方法，并給出了平面活塞型、球面聚焦型和圓柱面聚焦型換能器的實驗測量結果。近年來，他們又提出一種基于凸球面吸收靶的聲輻射力天平測量方法，成功測量了相控陣換能器的發射功率。2013年，北京理工大學沈洋等人首次針對平闆傾斜放置時的聲輻射力天平測量裝置進行了詳細的分析，給出了理論計算結果和相應的實驗驗證，大大簡化了對傳統聲輻射力天平測量的實驗要求。2014年，黃鴻鑫等人利用聲輻射力天平法制成了聲功率測量儀。
7結束語
聲輻射力既是聲場傳遞動量的反映，也是聲波非線性的重要展現。曆經了一個多世紀的探索與發展，關于聲輻射力的研究仍然方興未艾。一方面，不少學者已經給出了各種模型下聲輻射力的大量計算結果，為實際的聲操控提供了必不可少的理論基礎。另一方面，寬頻聲場的動态聲輻射力、多粒子的聲輻射力和血管壁等複雜邊界下的聲輻射力等成為了新的熱點，使得聲輻射力的基礎理論研究與實際操控結合得更加緊密。同時，随着聲渦旋理論的發展，基于軌道角動量的聲輻射力矩開辟了新的聲操控自由度，讓這一領域更加呈現出生機盎然的發展前景。
在生命科學、材料科學、納米科學等領域，對微觀粒子的無損精準操控已經成為不可或缺的核心技術。随着這一技術的不斷發展，人們已經不僅僅滿足于根據固定的聲場資訊來操控特定微粒，而是試圖面向特定應用場景、根據不同的操控需求反演所需要的聲場結構和操控系統，實作對操控對象的實時精準操控。在圖像回報技術和控制技術的輔助下，聲操控正逐漸走向微型化、自動化與智能化。另一方面，物體的複雜動力學特性暗含了其本身的實體參數資訊，利用聲輻射力反演其物性參數對精準操控和非接觸測量有重要的理論指導意義。作為聲操控的實體基礎，聲輻射力這隻看不見的“手”也将會在衆多領域中有更多用武之地。
參考文獻
[1] Hasegawa T，Saka K，Inoue N et al. J. Acoust. Soc. Am.，1988，83：1770
[2] King L V. Proc. Roya. Soc. London Ser. A，1934，147：212
[3] Awatani J. J. Acoust. Soc. Jpn.，1953，9：140
[4] Yosioka K. Acta Acust. United Ac.，1955，5：167
[5] Yosioka K，Hasegawa T，Omura A. Acustica，1969，22：145
[6] Hasegawa T，Yosioka K. J. Acoust. Soc. Am.，1969，46：1139
[7] Hasegawa T，Watanabe Y. J. Acoust. Soc. Am.，1978，46：1733
[8] Hasegawa T. J. Acoust. Soc. Am.，1979，65：32
[9] Wu J R，Du G H. J. Acoust. Soc，Am.，1990，87：997
[10] Marston P L. J. Acoust. Soc. Am.，2006，120：3518
[11] Marston P L. J. Acoust. Soc. Am.，2009，125：3539
[12] Mitri F G. Ann. Phys.，2008，323：2840
[13] Mitri F G. Ultrasonics，2009，49：794
[14] Mitri F G. J. Phys. A：Math. Theor.，2009，42：245202
[15] Mitri F G. Eur. Phys. J. E，2009，28(4)：469
[16] Silva G T. J. Acoust. Soc. Am.，2011，130：3541
[17] Baresch D，Thomas J，Marchiano R. J. Acoust. Soc. Am.，2013，133：25
[18] Sapozhnikov O A，Bailey M R. J. Acoust. Soc. Am.，2012，133：661
[19] Silva G T. IEEE Trans. Ultrason. Ferroelectr. Freq. Control，2015，62：576
[20] Wang J，Dual J. Ultrasonics，2012，52：325
[21] Qiao Y P，Zhang X F，Zhang G B. J. Acoust. Soc. Am.，2017，141：4633
[22] Qiao Y P，Zhang X F，Zhang G B. Wave Motion，2017，74：182
[23] Qiao Y P，Shi T Y，Zhang X F. Wave Motion，2018，83：111
[24] Hasheminejad S M，Sanaei R. J. Comput. Acoust.，2007，15：377
[25] Mitri F G. Ultrasonics，2016，66：27
[26] Mitri F G. J. Appl. Phys.，2015，118：214903
[27] Mitri F G. Wave Motion，2015，57：231
[28] Mitri F G. EPL，2015，112：34002
[29] Mitri F G. Ultrasonics，2017，74：62
[30] 王明升，歐陽傑，劉浩 等. 應用聲學，2020，39：550
[31] Zang Y C，Lin W J，Su C et al. Chin. Phys. B，2021，30：044301
[32] Zang Y C，Lin W J. Results Phys.，2020，16：102847
[33] Cai F Y，Meng L，Jiang C X et al. J. Acoust. Soc. Am.，2010，128：1617
[34] 孫秀娜，張小鳳，張光斌 等. 陝西師範大學學報(自然科學版)，2015，43：28
[35] Lee J W，Ha K L，Shung K K. J. Acoust. Soc. Am.，2005，117：3273
[36] Gor′kov L P. Sov. Phys. Doklady，1962，6：773
[37] Kundt A. Ann. Phys. Chem.，1866，127：497
[38] Hanson A R. Rev. Sei. Instrum.，1964，35：1031
[39] Whymark R R. Ultrasonics，1975，13：251
[40] Hatano H. IEEE Trans. Ultrason. Ferroelectr. Freq. Control，1994，41：772
[41] Jackson D P，Chang M H. Am. J. Phys.，2021，89：383
[42] Xie W J，Cao C D，Lu Y J et al. Phys. Rev. Lett.，2002，89：104304
[43] Ashkin A. Phys. Rev. Lett.，1970，24：156
[44] Ashkin A，Dziedzic J M. Appl. Phys. Lett.，1971，19：283
[45] Ashkin A，Dziedzic J M. Science，1975，187：1073
[46] Ashkin A，Dziedzic J M. Science，1987，235：385
[47] Ashkin A. Proc. Natl. Acad. Sci. USA，1997，94：4853
[48] Wu J R. J. Acoust. Soc. Am.，1991，89：2140
[49] Shi J J，Daniel A. Lab Chip，2009，9：2861
[50] Li S. Anal. Chem.，2013，85：5468
[51] Lee J，Teh S Y，Lee A et al. Appl. Phys. Lett.，2009，95：073701
[52] Lee J W，Lee C Y，Shung K K. IEEE Trans. Ultrason. Ferroelectr. Freq. Control，2010，57：2305
[53] Lee J W. Biotechnol. Bioeng.，2011，108：1643
[54] Ma Z，Collins D J，Guo J et al. Anal. Chem.，2016，88：11844
[55] Chen X Y，Lam K H，Chen R M et al. Biotechnol. Bioeng.，2017，114：2637
[56] Zhao L R，Kim E S. Focused ultrasound transducer with electrically controllable focal length. Proceedings IEEE Micro Electro Mechanical Systems，2018：1088
[57] Marzo A，Drinkwater B W. PNAS，2019，116：84
[58] 梁彬，程建春. 實體，2017，46：658
[59] Maidanik G. J. Acoust. Soc. Am.，1958，30：620
[60] Zhang L K，Marston P L. J. Acoust. Soc. Am.，2011，129：1679
[61] Zhang L K，Marston P L. Phys. Rev. E，2011，84：065601
[62] Zhang L K，Marston P L. J. Acoust. Soc. Am.，2011，129：2381
[63] Demore C E M，Yang Z Y，Volovick A et al. Phys. Rev.Lett.，2012，108：194301
[64] Silva G T，Lobo T P，Mitri F G. EPL，2012，97：54003
[65] Mitri F G. Phys. Rev. E，2012，85：026602
[66] Mitri F G，Lobo T P，Silva G T. Phys. Rev. E，2012，86：059902
[67] Zhang L K. Phys. Rev. Appl.，2018，10(3)：034039
[68] Gong Z X，Marston P L. Phys. Rev. Appl.，2019，11：064022
[69] 臧雨宸，林偉軍，蘇暢 等. 聲學學報，2021，46：92
[70] Zang Y C，Lin W J，Zheng Y F et al. J. Sound Vib.，2021，509：116261
[71] Wang T，Ke M，Li W et al. Appl. Phys. Lett.，2016，109：123506
[72] Zhang R Q，Guo H L，Deng W Y et al. Appl. Phys. Lett.，2020，116：123503
[73] Hong Z，Zhang J，Drinkwater B W. Phys. Rev. Lett.，2015，114：214301
[74] Jiang X，Zhao J J，Liu S L et al. App. Phys. Lett.，2016，108：203501
[75] Ophir J，Cespedes I，Ponnekanti H et al. Ultrasonic Imaging，1991，13：111
[76] Maleke C，Konofagou E E. Phys. Med. Bio.，2008，53：1773
[77] Sarvazyan A P，Rudenko O V，Swanson S D et al. Ultrasound Med. Bio.，1998，24：1419
[78] 曾婕，吳莉莉，鄭榮琴 等. 中華超聲醫學雜志，2012，9：781
[79] Nightingale K，Soo M S，Nightingale R et al. Ultrasound Med. Bio.，2002，28：227
[80] 歐冰，羅葆明，智慧 等. 中國醫療器械資訊，2012，6：18
[81] Zieniuk J，Chivers R C. Ultrasonics，1976，14：161
[82] 楊德俊，陳沈理，何卓斌 等. 電子品質，2018，7：90
[83] 壽文德，餘立立，胡濟民 等. 聲學技術，2012，31：107
[84] 壽文德，賈龍洋，吉祥 等. 聲學技術，2019，38：629
[85] 沈洋，胡亞琴，胡毅. 科協論壇，2013，12：201
[86] 黃鴻新，徐紅蕾，胡昌明 等. 實驗技術與管理，2014，31：79
轉載内容僅代表作者觀點
不代表中科院實體所立場
來源：中國實體學會期刊網
編輯：fiufiu