劉曉宙 龔秀芬 盧 瑩

（近代聲學(xué)教育部重點(diǎn)實(shí)驗(yàn)室，南京大學(xué)聲學(xué)研究所，南京 210093）

引言

隨著醫(yī)學(xué)超聲的發(fā)展，超聲不僅廣泛地應(yīng)用于疾病的診斷，還應(yīng)用于腫瘤等的治療，特別是高強(qiáng)度聚焦超聲（HIFU）的出現(xiàn)，為治療超聲的應(yīng)用開辟了一個(gè)新的領(lǐng)域［1-2］。目前，HIFU的應(yīng)用雖已經(jīng)獲得了巨大的成功，但若干基礎(chǔ)問題亟待解決，如治療超聲中的非線性效應(yīng)、空化效應(yīng)、溫度場(chǎng)的估計(jì)及無損測(cè)溫等。吳君汝等理論分析了高斯聚焦聲束引起的溫度提升［3］，及高斯聚焦聲束在組織和骨界面引起的溫度提升［4］；Bacon實(shí)驗(yàn)研究了聚焦超聲引起的溫度提升［5］；吳君汝等實(shí)驗(yàn)測(cè)量了聚焦和非聚焦超聲引起的溫度變化［6］；薛洪惠等研究了聚焦超聲波通過層狀生物媒質(zhì)的二次諧波聲場(chǎng)［7］；錢盛友等使用Madsen的方法計(jì)算了聚焦超聲的聲場(chǎng)并研究了在熱療中由超聲作用下引起的溫度場(chǎng)分布［8］；但軸對(duì)稱的計(jì)算無法直接應(yīng)用于研究肋骨這樣的非軸對(duì)稱障礙物后方的非線性聲場(chǎng)，李俊倫等研究了超聲波通過肋骨后的三維聲場(chǎng)［9］；朱曉峰等用角譜近似的方法研究了聚焦超聲穿過層狀生物媒質(zhì)后的二次諧波聲場(chǎng)分布［10］。劉曉宙等理論和實(shí)驗(yàn)得到了超聲波通過肋骨后的溫度場(chǎng)［15］。目前的研究主要是利用Khokhlov-Zabolotkaya-Kuznetsov（KZK）方程［11］和生物熱傳導(dǎo)方程（Pennes方程）［12］來研究HIFU輻照后生物體的溫度場(chǎng)［13-14］。但由于活塞聚焦換能器聲場(chǎng)的復(fù)雜性，活塞聚焦超聲引起的溫度提升的理論研究往往是通過數(shù)值計(jì)算獲得，溫升與換能器參數(shù)的關(guān)系并不明確。本研究在非線性KZK方程基礎(chǔ)上，利用高斯聲束疊代法［16］描述活塞聚焦換能器的聲場(chǎng)，即用一系列的高斯聚焦聲束來描述活塞聚焦換能器的聲源。由于高斯聲源的特殊性，活塞聚焦換能器的基波聲場(chǎng)仍為一系列高斯聲束的疊加，而二階聲場(chǎng)則為這一系列的高斯聲束的自作用和相互作用之和［17-18］。在此基礎(chǔ)上，應(yīng)用Pennes提出的生物熱傳導(dǎo)方程，理論上得到聚焦活塞聲場(chǎng)引起的溫度提升的解析表達(dá)式。理論計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果進(jìn)行了對(duì)比，并討論了換能器參數(shù)和傳播媒質(zhì)特性對(duì)溫度變化的影響。

1 理論模型

要獲得聚焦活塞聲場(chǎng)引起的溫度場(chǎng)的分布，首先要得出聚焦活塞聲場(chǎng)的分布。在Parabolic近似及準(zhǔn)線性近似條件下，聚焦活塞換能器的聲場(chǎng)滿足KZK方程［11］：

如圖1活塞聚焦聲源的半徑為a，聚焦長度為R，聲源質(zhì)點(diǎn)速度為u0。利用高斯聲束疊代法，用一系列高斯聲束模擬活塞聲源［16］

圖1 活塞聚焦聲源結(jié)構(gòu)Fig.1 Structure of piston focused source

將式（3）代入式（1）和式（2）可得

由聲場(chǎng)到溫度場(chǎng)變換的依據(jù)為Pennes提出的生物熱傳導(dǎo)方程［12］：

式中，qv=2αI為媒質(zhì)的能量吸收，α為聲衰減系數(shù)，I=p2/2ρc為聲強(qiáng)，ρ為媒質(zhì)的密度，K為熱傳導(dǎo)系數(shù)，cv為熱容量，τ為灌注時(shí)間常數(shù)。當(dāng)時(shí)間足夠長即進(jìn)入穩(wěn)態(tài)后，式（6）簡(jiǎn)化為

對(duì)于小體積聲源dv來說，上述方程的解為式中，ξ=r/a，σ=z/r0，r0為瑞利長度，L=（Kτ/cv）1/2為灌注長度。

將聚焦活塞聲源的一階表達(dá)式代入式（8），經(jīng)過化簡(jiǎn)得到由基波引起的溫度提升為

同理可以得到由二次諧波引起的溫度改變?yōu)?/p>

式中，q2為二次諧波引起的媒質(zhì)的能量吸收。

2 數(shù)值計(jì)算及討論

在理論計(jì)算和實(shí)驗(yàn)中采用的換能器參數(shù)和樣品參數(shù)相同，活塞聚焦換能器的半徑a=12mm，焦距R=58mm，中心頻率為3MHz；傳播媒質(zhì)為脂肪和肝臟組織，脂肪組織的參數(shù)為：熱傳導(dǎo)系數(shù)K=0.6w/（m·℃），一階聲衰減系數(shù)α1=71.7Np/m，熱容量cv=3 800J/（kg·℃），密度ρ=950kg/m3，聲速c=1 445m/s，非線性系數(shù)β=6.5；肝臟組織的參數(shù)為：熱傳導(dǎo)系數(shù)K=0.5w/（m·℃），一階聲衰減系數(shù)α1=57.1Np/m，熱容量cv=3 113J/（kg·℃），密度ρ=1 224kg/m3，聲速c=1 614m/s，非線性系數(shù)β=4.65。二階聲衰減系數(shù)α2≈2α1。

計(jì)算中先假設(shè)媒質(zhì)為理想狀態(tài)，L=∞，即無熱耗散，由計(jì)算結(jié)果可知在初始聲強(qiáng)為3.0w/cm2下，脂肪組織中二階聲場(chǎng)引起的溫度提升和一階聲場(chǎng)引起的溫度相比，相差3個(gè)數(shù)量級(jí)，可見在聲源聲強(qiáng)不是很大且聚焦程度較弱的情況下，溫度的升高主要決定于基波聲場(chǎng)。圖2和圖3分別是由一、二階聲場(chǎng)共同作用引起的在脂肪和肝臟組織中的軸向的溫度穩(wěn)態(tài)分布，其中實(shí)線為初始聲強(qiáng)為3.0w/cm2的理論結(jié)果；長劃線為初始聲強(qiáng)為1.3w/cm2的理論結(jié)果；短劃線為初始聲強(qiáng)為0.2w/cm2的理論結(jié)果；圖中標(biāo)有在不同初始聲強(qiáng)下的實(shí)驗(yàn)結(jié)果?？梢婋S著聲源聲功率的增加，軸向溫度最高點(diǎn)提升的數(shù)值增加，但是最高點(diǎn)的位置基本保持不變。經(jīng)過數(shù)值計(jì)算，在聲強(qiáng)6.2w/cm2時(shí)，脂肪組織中二階聲場(chǎng)引起的溫度提升和一階聲場(chǎng)引起的溫度相比，相差2個(gè)數(shù)量級(jí)，二級(jí)聲場(chǎng)引起的溫度可以忽略，因此當(dāng)聲強(qiáng)低于6.2w/cm2時(shí)，溫度的提升由基波聲場(chǎng)決定。在聲強(qiáng)11.0w/cm2時(shí)，脂肪組織中二階聲場(chǎng)引起的溫度提升和一階聲場(chǎng)引起的溫度相比，只相差1個(gè)數(shù)量級(jí)，二級(jí)聲場(chǎng)引起的溫度不能忽略，因此當(dāng)聲強(qiáng)高于11.0w/cm2時(shí)，溫度的提升將由基波聲場(chǎng)和二次諧波聲場(chǎng)共同決定。

圖2 脂肪組織中聲源聲功率變化引起的溫度分布變化Fig.2 The distribution of temperature rise in fatty tissue induced by sound source

圖3 肝臟組織中聲源聲功率變化引起的溫度分布變化Fig.3 The distribution of temperature rise in liver tissue induced by sound source

在理論模型基礎(chǔ)上，計(jì)算換能器參數(shù)（含聲功率、頻率、半徑和焦距）和傳播媒質(zhì)的特性對(duì)溫度變化的依賴關(guān)系。

2.1 換能器參數(shù)對(duì)溫度變化的影響

為討論換能器半徑、焦距和聲源的頻率對(duì)溫度提升的影響，采用無量綱焦距σR=R/rl［3］，rl為修正瑞利長度：rl=［a2ω/（2c）］（1-c/cl），cl為換能器上固體透鏡中的聲速，如換能器不加固體透鏡，則此值為無窮大。隨著σR的減小，聚焦程度加強(qiáng)，從而溫度最高點(diǎn)的位置將向焦點(diǎn)靠攏，同時(shí)溫度最高點(diǎn)溫度值增加。圖4所示為軸向穩(wěn)態(tài)溫度分布與換能器半徑的關(guān)系。隨著半徑a的增大，溫度最高點(diǎn)的位置將向焦點(diǎn)靠攏，同時(shí)溫度最高點(diǎn)溫度值增加。在其他參數(shù)不變的情況下，換能器的半徑增加，rl增大，σR減小，提供能量的面積增大，溫度最高點(diǎn)的數(shù)值將增加。軸向穩(wěn)態(tài)溫度分布與焦距關(guān)系的計(jì)算結(jié)果示于圖5上。隨著焦距R的減小，溫度最高點(diǎn)的位置減小，同時(shí)溫度最高點(diǎn)溫度值增加。在其他參數(shù)不變的情況下，換能器的焦距減小，聚焦程度增強(qiáng)，σR減小。圖6為軸向穩(wěn)態(tài)溫度分布與頻率的關(guān)系的計(jì)算結(jié)果。因?yàn)樾拚鹄L度是一個(gè)與頻率有關(guān)的量，所以隨著發(fā)射頻率的增加，rl增大，σR減小，溫度最高點(diǎn)溫度值增加。

圖4 換能器的半徑變化引起的溫度分布變化Fig.4 The distribution of temperature rise caused by the change of the radius of the transducer

2.2 傳播媒質(zhì)的影響

隨著灌注長度L的增加，血液灌注熱量減小，當(dāng)L達(dá)到無窮大時(shí)達(dá)到理想狀況，血液不發(fā)生灌注，即血液不帶走能量，在同等狀況下，溫度提升達(dá)到最大值，如圖7所示。

圖5 換能器的焦距變化引起的溫度分布變化Fig.5 The distribution of temperature rise caused by the change of the focal length of the transducer

圖6 換能器的發(fā)射頻率變化引起的溫度分布變化Fig.6 The distribution of temperature rise caused by the change of the frequency of the transducer

圖7 灌注長度L變化引起的溫度分布變化Fig.7 The distribution of temperature rise caused by the change of the perfusion length of the medium

以上的研究表明，聲源的聲功率、頻率、換能器的半徑和焦距以及傳播媒質(zhì)的灌注長度將影響溫度的增加和溫度最高點(diǎn)的位置。在醫(yī)療應(yīng)用中，傳播媒質(zhì)的性質(zhì)一定，可以通過改變換能器的參數(shù)

（聲功率、頻率、半徑和焦距）來調(diào)整，得到適合于實(shí)際需要的溫度和溫度最高點(diǎn)的位置。本方法可預(yù)測(cè)由于超聲引起的最大溫升，為實(shí)際超聲治療腫瘤提供依據(jù)。

3 實(shí)驗(yàn)與結(jié)果

為驗(yàn)證理論分析的準(zhǔn)確性，建立了如下的實(shí)驗(yàn)系統(tǒng)，見圖8，函數(shù)發(fā)生器（Agilent 33250A，美國）發(fā)射頻率為3MHz的連續(xù)正弦信號(hào)，通過功率放大器（ENI A150，美國）加到聚焦活塞超聲換能器上，由此產(chǎn)生的超聲波作用于生物組織中，生物組織的溫升由MM100（Yokogawa，日本）熱電偶記錄，在樣品中插入四個(gè)探針，同時(shí)記錄樣品中4點(diǎn)的溫度提升。接收的超聲信號(hào)由一寬帶針狀水聽器（NP1000，美國）接收并經(jīng)前置放大器放大30dB后，由數(shù)字存貯示波器（Agilent 54810，美國）采樣，可得到空間聲場(chǎng)的分布。實(shí)驗(yàn)樣品為計(jì)算中所提到的脂肪和肝臟樣品。

圖8 實(shí)驗(yàn)系統(tǒng)方框圖Fig.8 Block diagram of the experimental system

實(shí)驗(yàn)結(jié)果見圖2和圖3，生物組織的最大溫升在軸向上的分布與理論的結(jié)果相比，理論計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果基本吻合。實(shí)驗(yàn)中，加熱足夠長時(shí)間后，溫升基本保持不變后，在脂肪組織內(nèi)溫度大約提升14.1℃，在肝臟組織內(nèi)的溫度大約提升為10.3℃，實(shí)驗(yàn)和理論計(jì)算結(jié)果均表明溫度最高點(diǎn)出現(xiàn)在幾何焦點(diǎn)前的微小偏移處，脂肪組織中的溫升大于在肝臟組織中的溫升，原因是脂肪組織中的聲衰減系數(shù)大于肝臟組織中的聲衰減系數(shù)，而且脂肪組織的熱容量大于肝臟組織中的熱容量，而它們的熱傳導(dǎo)系數(shù)相差不大。

4 結(jié)論

基于Pennes方程及高斯聲束疊加法，本研究提出了計(jì)算聚焦活塞聲源引起的溫度穩(wěn)態(tài)分布的新方法。理論計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果基本相符，在理論分析的基礎(chǔ)上，本研究還計(jì)算了換能器參數(shù)和傳播媒質(zhì)特性對(duì)溫度提升和位置的影響，溫升隨換能器的功率、頻率和半徑增加而增加，隨焦距的增加而減小，隨媒質(zhì)的灌注長度的增加而增加。本研究的結(jié)果對(duì)超聲用于生物組織治療提供了有用的信息，具有參考價(jià)值。

［1］Wang Zhibiao，Bai Jing，Li Faqin，et al.Study of a”biological focalregion” of high-intensity focused ultrasound［J］.Ultrasound in Med＆ Biol，2003，29（5）：749-754.

［2］Hynynen K.Review of ultrasound therapy［J］.Proc IEEE International Ultrasonic Symposium，1997，2：1305-1313.

［3］Wu Junru，Du Gonghuan.Temperature elevation generated by a focused Gaussian beam of ultrasound［J］.Ultrasound in Med＆Biol，1990，16：489-498.

［4］Wu Junru，Du Gonghuan.Temperature elevation generated by a focused Gaussian ultrasonic beam at a tissue-bone interface［J］.J Acoust Soc Am，1990，87：2748-2755.

［5］Bacon DR，Carstensen EL.Increased heating by diagnostic ultrasound due to nonlinear propagation［J］.J.Acoust.Soc.Am，1990，88：26-34.

［6］Wu Junru，Chase JD，Zhu Zheminet al.Temperature rise in a tissue-mimicking material generated by unfocused and focused ultrasonic transducers［J］.Ultrasound in Med ＆Biol，1992，18：495-512.

［7］薛洪惠，劉曉宙，龔秀芬，等.生物媒質(zhì)中的二次諧波聲場(chǎng)的理論與實(shí)驗(yàn)研究［J］.物理學(xué)報(bào)，2005，54：5233-5238.

［8］錢盛友，王鴻樟.聚焦超聲源對(duì)生物媒質(zhì)加熱的理論研究［J］.物理學(xué)報(bào)，2000，50：501-506.

［9］李俊倫，劉曉宙，龔秀芬，等.條狀障礙物對(duì)超聲非線性聲場(chǎng)的影響研究［J］.物理學(xué)報(bào)，2006，55（6）：2809-2814.

［10］Zhu Xiaofen，Zhou Lin，Zhang Dong，etal.Nonlinear propagation of focused ultrasound in layered biological tissues based on angular spectrum approach［J］.Chin Phys，2005，8：1594-1599.

［11］Lucas BG，Muir TG.Field of a finite-amplitude focusing source［J］.J.Acoust Soc Am，1983，74：1522-1528.

［12］Pennes HH.Analysis of tissue and arterial blood temperature in the resting forearm［J］.J Appl Physiol，1948，1：93-122.

［13］Curra FP，Mourad PD，Khokhlova VA，et alNumerical simulations of heating patterns and tissue temperature response due to high-intensity focused ultrasound［J］.IEEE Transaction on Ultrasonics，F(xiàn)erroelectrics，and Frequency Control，2000，47：1077-1088.

［14］Meaney PM，Cahill MD，GR TH.The intensity dependence of lesion position shift during focused ultrasound surgery［J］.Ultrasound in Med ＆ Biol，2000，26：441-450.

［15］Liu Xiaozhou，Yin Chang，Gong Xiufen，et al.Theoretical and experimental study on temperature elevation behind ribs caused by weakly focused ultrasound［J］.Ultrasound in Med ＆ Biol，2010，36（10）：1704-1712.

［16］Wen JJ，Breazeale MA.Diffraction beam field expressed as the superposition of Gussian beam［J］.J Acoust Soc Am，1998，83（5）：1752-1756.

［17］Ding Desheng，Shui Yongan，Lin Jingbo，et al.A simple calculation approach for the second harmonic sound field generated by an arbitrary axial-symmetric source［J］.J.Acoust Soc Am，1996，100（2）：727-733.

［18］Zhang Dong，Gong Xiufen，Zhang Bo.Second harmonic sound field after insertion of a biological tissue sample［J］.J.Acoust Soc Am，2002，111：45-48.