顱內(nèi)溫控射頻消融的熱損傷有限元分析

洪焦，李鵬，高宏建，吳水才

北京工業(yè)大學(xué) 生命科學(xué)與生物工程學(xué)院，北京 100124

顱內(nèi)溫控射頻消融的熱損傷有限元分析

洪焦，李鵬，高宏建，吳水才

北京工業(yè)大學(xué) 生命科學(xué)與生物工程學(xué)院，北京 100124

目的 運用有限元方法模擬顱內(nèi)溫控射頻消融過程中的溫度場分布，建立臨床常用的激勵模型，以提高射頻消融治療顱內(nèi)病灶的效果。方法 基于2種激勵分布來建立3種有限元模型，并對其溫度變化、熱損傷范圍進行對比分析。結(jié)果 消融區(qū)域在消融開始60s內(nèi)增長迅速，60s后增長緩慢，且固態(tài)介質(zhì)模型的消融面積＞相應(yīng)液態(tài)模的消融面積。結(jié)論 恒溫激勵源的設(shè)計可用于模擬溫控射頻消融治療。

射頻消融儀；腦病灶；射頻消融；有限元法（FEM）；熱模型

1 研究背景

射頻消融術(shù)（Radio Frequency Ablation，RFA）已廣泛用于治療心律失常[1]、肝臟腫瘤[2]和腦部病變。顱內(nèi)病灶 （蒼白球（帕金森病）[3]、海馬-杏仁核（顳葉癲癇）、腫瘤[4]等）的射頻消融治療是將射頻消融針（頻率460 kHz）通過腦立體定向、顱骨鉆孔插入病灶部位作為熱源，在數(shù)分鐘內(nèi)產(chǎn)生局部高溫（＞50℃），使病灶組織凝固變性，達到原位毀損病灶的目的。合理有效地利用熱消融治療，可明顯提高對常規(guī)放、化療不敏感腫瘤的控制率，改善治愈率和生存率。同時，熱療還可在不降低現(xiàn)有治療手段和療效的基礎(chǔ)上，有效緩解癥狀，改善全身狀況，增加機體免疫功能[5]。目前，射頻消融治療中加熱功率和時間的確定還依賴于醫(yī)生的經(jīng)驗，人為因素影響很大。因此，對射頻消融熱損傷區(qū)域（組織溫度≥50℃區(qū)域）進行仿真分析、實現(xiàn)術(shù)前的手術(shù)規(guī)劃具有重要的臨床意義。

影響射頻消融溫度場模擬結(jié)果的因素：

（1） 水循環(huán)冷卻的影響。水冷循環(huán)的加入可增加有效消融時間，增大消融范圍[6]。

（2） 血液灌流的影響。血液灌流對模型電特性影響不大，但在一定程度上降低了模型消融溫度。

（3） 生物組織熱電參數(shù)隨溫度變化的影響。變化的熱電參數(shù)在一定程度上改變熱傳導(dǎo)率和射頻能量吸收率，從而改變消融溫度[7-9]。

目前，針對腦部病灶射頻消融的離體實驗中多使用體模或動物腦。由于動物腦組織存在體積小、消融區(qū)不均勻并且離體后容易變質(zhì)等因素，消融實驗效果不佳。因此，研究者更傾向于使用體模進行消融實驗。Eriksson等人[10-11]制造了一種透明、均勻、可重復(fù)使用的混合白蛋白體模，其熱物性參數(shù)與人腦組織相近，但其性質(zhì)為液態(tài)，實際消融過程中熱對流作用帶走較多熱量。針對這種情況J. Wren等人[12-13]在液態(tài)體模中混合明膠，使其在室溫下可保持膠狀，屏蔽了對流作用的影響。本文針對以上兩種不同情況設(shè)計了兩種能量分布，利用有限元方法計算恒溫激勵下蛋白質(zhì)體模和腦組織消融溫度分布規(guī)律，比較消融溫度場以及熱損傷區(qū)的變化。

2 原理與方法

2.1 恒溫激勵模型

目前，臨床腦病灶射頻消融儀大多為溫度控制型，即實際手術(shù)過程中醫(yī)生只需要規(guī)定消融溫度和時間，醫(yī)生只能根據(jù)經(jīng)驗來判定消融效果，因此消融能量不宜過高，但這可導(dǎo)致消融不徹底、易復(fù)發(fā)等結(jié)果。

研究結(jié)果表明：由于消融針的電導(dǎo)率比白蛋白體模以及腦組織的電導(dǎo)率大很多，使射頻能量大部分集中于消融針外部，而消融針內(nèi)部不會得到射頻能量，其升溫主要是通過熱傳導(dǎo)；消融針周圍組織獲得的射頻能量正比于1/r4，其中r為針半徑，這說明射頻能量主要集中于消融針周圍，且徑向上能量衰減非常快。對于溫控射頻消融儀，治療初始10 s左右即可達到設(shè)定溫度。由于射頻能量不恒定，很難通過射頻能量的有限元分析來模擬溫度變化。

本研究設(shè)計了兩種恒溫能量刺激模型，模擬針尖附近能量分布規(guī)律，以期找到合適的模型。其能量激勵分布區(qū)域，見圖1。圖1(a)用于模擬存在熱對流時的消融情況，其能量主要分布于針尖兩側(cè)；圖1(b)用于模擬不存在熱對流時的消融情況，其能量包裹在針尖四周，且絕緣桿與針尖連接處能量較強。根據(jù)不同材料，設(shè)計了3種有限元計算模型：白蛋白溶液（模型I，激勵分布a）；凝膠狀白蛋白（模型II，激勵分布b）；腦灰質(zhì)（模型III，激勵分布c）。

2.2 射頻針與消融區(qū)模型

腦病灶射頻消融中所用的射頻針結(jié)構(gòu)，見圖2。消融針總長44 mm，直徑2 mm；前段金屬發(fā)射部分為醫(yī)用合金，長4 mm；絕緣桿為聚氨酯，長40 mm。

消融區(qū)組織模型，見圖3。模型半徑30 mm，高度60 mm。組織邊界條件如下：組織初始溫度Tu=37 ℃，恒溫能量激勵溫度：Tc=90 ℃。射頻針和組織材料的熱電參數(shù)，見表1。腦灰質(zhì)密度和比熱來自Duck[14]的數(shù)據(jù)，導(dǎo)熱率來自Chato[15]的文獻數(shù)據(jù)，射頻針和絕緣桿的參數(shù)來自Louay Chachati等人[16]的文獻。

表1 射頻針和組織材料的熱電參數(shù)

2.3 模型邊界條件

考慮3種恒溫激勵源的有限元模型特點，運用有限元分析軟件ANSYS 13.0求解，模擬使用軸對稱分析和瞬態(tài)熱傳導(dǎo)模塊。使用的射頻針，發(fā)射針長4 mm，絕緣桿長30 mm，探針直徑2 mm。腦組織為30 mm×60 mm區(qū)域，記錄A（0.5 mm）、B（2.5 mm）、C（5 mm）三點的溫度變化，其中A點位于針內(nèi)部，距中軸線0.5 mm，B、C兩點位于針外部，距中軸線距離分別為2.5 mm、5 mm，見圖4。記錄消融開始后60 s和120 s時刻的溫度分布和熱損傷區(qū)域大小（以50 ℃為邊界的橢圓為參考，超過50 ℃后組織變性壞死）。模型初始和邊界條件如下：① 腦組織初始溫度為37 ℃。② 激勵源溫度設(shè)定為90 ℃。③ 假設(shè)消融時間為120s。④ 假設(shè)腦組織為腦灰質(zhì)，且均勻。⑤ 忽略血液灌流對消融效果的影響。⑥ 忽略組織代謝產(chǎn)熱對消融效果的影響。

3 結(jié)果與討論

圖5～7（a）記錄了3種模型的A、B、C 3點處的溫度時間變化曲線，從圖中可以看出，由于醫(yī)用金屬的熱傳導(dǎo)率遠大于腦組織，因此，A點可視為激勵源處溫度變化；B點較C點更靠近激勵源，初始時刻升溫迅速，隨著時間的推移升溫幅度減小直至基本不變；而C點處溫升緩慢，升溫速度基本保持恒定。圖5～7（b）和（c）分別為3種模型消融60 s和120 s時刻的溫度場分布，從中可以看出，消融區(qū)域近似為橢球形，并且消融前60 s熱損傷區(qū)域增長迅速，60～120 s時間段里損傷區(qū)的增長變得平緩。

比較模型I、II、III的B點和C點溫度隨時間的變化，見圖8。三種模型B、C點的溫升趨勢相似，但是模型II和III的溫度明顯高于模型I。說明白蛋白溶液體模存在的熱對流作用對消融結(jié)果影響較大。同時對比相同能量激勵下的模型II和III，可見模型II的溫度和消融區(qū)域要稍大于模型III，主要原因是其材料熱物性參數(shù)存在差異，明膠狀白蛋白體模對熱量的吸收能力（k/ρc）要大于腦灰質(zhì)。3種模型60 s和120 s時刻熱損傷區(qū)域大小（以50 ℃為消融邊界），見表2。其中X為消融區(qū)短半軸，Y為消融區(qū)長半軸，S為橢圓消融區(qū)的面積。可見模型II和III的消融面積分別較模型I大84.9%和68.8%，模型II的消融面積比模型III大9.2%。消融開始60 s內(nèi)消融區(qū)增長迅速，60 s之后增長緩慢。比較3種模型60 s和120 s兩時刻消融區(qū)域大小，其增長率分別為21.5%、26.8%、20.6%。

表2 3種模型熱損傷區(qū)域大小

4 結(jié)論

射頻消融過程中，為了有效控制熱損傷區(qū)域的大小，溫度監(jiān)控變得異常重要。實際的射頻消融儀只能通過探針內(nèi)部的熱電偶測量溫度，無法測量探針外部溫度，也無法檢測實際消融區(qū)的大小。同時隨著消融能量和介質(zhì)熱物性參數(shù)的變化，消融結(jié)果存在較大差異，臨床醫(yī)生多數(shù)情況下只能根據(jù)經(jīng)驗判斷消融能量和時間。因此，這就需要利用有限元分析工具對特定環(huán)境下的熱計量進行評估。

本研究對腦病灶射頻消融熱損傷區(qū)進行了有限元仿真。通過模擬臨床常用溫控射頻消融儀的工作機制設(shè)計出兩種能量分布，然后利用ANSYS軟件模擬液態(tài)白蛋白體模、凝膠白蛋白體模和腦灰質(zhì)的射頻消融情況，并比較消融溫度、溫度場分布和熱損傷區(qū)大小（邊界溫度超過50 ℃區(qū)域）。研究結(jié)果表明：液態(tài)體模（如蛋清和白蛋白溶液）內(nèi)部存在的熱對流作用對消融結(jié)果影響較大，實驗中要盡量避免使用液態(tài)體模；介質(zhì)材料的熱物性參數(shù)對熱劑量的計算至關(guān)重要；初始時刻介質(zhì)對能量的吸收能力強，熱損傷范圍迅速擴大，消融后1/2段吸收能力下降，熱損傷區(qū)域形成緩慢。

[1]Tungjitkusolmun S,Woo EJ,Cao H,et al.Finite element analyses of uniform current density electrodes for radio-frequency cardiac ablation[J].IEEE Trans Biomed Eng,2001,47(1):32-40.

[2]Goldberg SN.Radiofrequency tumor ablation:principles and techniques[J].Eur J Ultrasound,2001,13(2):129-147.

[3]Cosman E.Radiofrequency lesions.In:Gildenberg PL,Tasker R (eds) Textbook of stereotactic and functional neurosurgery[C]. Quebecor Printing, Kingsport,1996:973-985.

[4] 王啟弘,楊富明,楊世春.射頻毀損腦深部腫瘤的應(yīng)用研究[J].功能性和立體定向神經(jīng)外科雜志,1997,10(4):14-16.

[5]王梅.腫瘤熱療的研究進展[J].海南醫(yī)學(xué)院學(xué)報,2010, 16(2):245-247.

[6]徐曉菲,白景峰,陳亞珠.水冷式射頻消融過程中水冷溫度對熱損傷的影響[J].北京生物醫(yī)學(xué)工程,2008,27(4):393-395.

[7]Isaac Chang. Finite element analysis of hepatic radiofrequency ablation probes using temperature-dependent electrical conductivity[J].BioMed Central. Bio Medical Engineering OnLine,2003.

[8] 范曉筠,白景峰,陳亞珠.單針水冷式射頻消融過程中組織熱損傷區(qū)域的有限元分析[J].中國醫(yī)學(xué)物理學(xué)雜志,2007,24(3):197-199.

[9]侯惠亮,鄭小溪.熱療過程體溫監(jiān)測器的研制[J].中國醫(yī)療設(shè)備,2011,26(8):32-35.

[10] Eriksson O,Wardell K,Bylund N,et al.In vitro evaluation of brain lesioning electrodes using a computer-assisted video system[J].Neurological Res,1999,21(1):89-95.

[11] Eriksson O,Wren J,Loyd D,et al.A comparison between in vitro studies of protein lesions generated by brain electrodes and finite element model simulations[J].Med Biol Eng Comput,1999,37: 737-741.

[12]J Wren I O,Eriksson z K,W&rdell z D.Loyd 1.Analysis of temperature measurement for monitoring radio-frequency brain lesioning [J].Med.Biol.Eng.Comput,2001,39:255-262.

[13]J.Wren.On modelling and simulation of perfused tissue during thermal treatment-thermal analysis of blood perfusion and lesioning during ablative neurosurgery[M].Lic.Eng.Thesis, LinkGpings Univeristet,2000.

[14] Duck AF.Physical properties of tissue [M].The University Press, Cambridge,1990.

[15] Chato JC.Selected thermophysical properties of biological materials[M].In:Shitzer Eberhart eds.Heat Transfer in Medicine and Biology.Newyork;plenum press,1985,2:413-418.

[16] Louay Chachati,Andrew S.Wright,David M.Mahvi.Hepatic radiofrequency ablation with internally cooled probes:effect of coolant temperature on lesion size[J].IEEE Transactions on Biomedical Engineering,2003,50(4):493-500.

Finite Element Analysis of Encranial Thermal Injury Region with Temperature-controlled Radio Frequency Ablation

HONG Jiao, LI Peng,GAO Hong-jian, WU Shui-cai
College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China

Objective This paper introduces f nite element application of themerature-controlled radio frequency ablation(RFA) to simulate the temperature distribution of intracranial lesion and build the excitation models found in clinical practice, thus to improve the treatment effectiveness. Methods Three f nite element models with two different exitations were f rst bulit, and then the key parameters including temperatures range and thermal injury were compared. Results The ablation regions grew fast in the f rst 60 seconds, then it became slow. In addition, the ablation regions of solid medium models were much larger than those of liquid models. Conclusion The design of homeothermal excitation may be effective for simulating temperature-controlled radio frequency abltion.

radio frequency ablation device; brain lesions; radio frequency ablation(RFA); f nite element methods(FEM); thermal model

R318

A

10.3969/j.issn.1674-1633.2012.08.009

1674-1633(2012)08-0037-03

2012-03-01

北京市自然科學(xué)基金（3112005）資助。

本文作者：在讀碩士研究生。

吳水才，教授，博士生導(dǎo)師。

作者郵箱：15810948582@139.com