Real-Time Calculation Method for Temperature Distribution of Temperature-Controlled Radiofrequency Ablation

doi:10.1007/s12204-022-2481-y

Abstract

Abstract: Precise temperature control and temperature distribution prediction are of great significance forradiofrequency ablation. This research proposes a real-time calculation method for the temperature distribution of radiofrequency ablation combined with proportional-integral temperature control. The thermo-electricalcoupling was simplified into a linear relationship based on the study of the influence of temperature-dependentelectrical conductivity and thermal conductivity on the PI-controlled radiofrequency ablation temperature distribution, which increases the computational efficiency by 150 times. The average calculation time for radiofrequencyablation of 10 min is about 23 s, and the difference between the calculation results of this method and that fromCOMSOL Multiphysics is no more than 1 ?C. This method is not only used for single-probe, but also for doubleprobe radiofrequency ablation in this paper.

Key words: radiofrequency ablation, temperature distribution, proportional-integral (PI) controller

摘要： 精确的温度控制和温度场预测对射频消融具有重要意义，本研究提出了一种比例-积分温度控制射频消融过程中的温度场实时计算方法，基于电导率对温控射频消融影响的研究，将热电耦合的计算关系转换为线性的计算关系，使计算效率提升了150倍。对于10 min的射频消融，计算时间约为23 s，计算结果和相同情况下的COMSOL Multiphysics相比不超过1 ℃。本研究中该计算方法不仅被用于单针射频消融，还被用于双针射频消融中温度场的计算。

关键词: 射频消融，温度场分布，比例-积分控制

CLC Number:

TP 312

WANG Xuewei (王雪维)，WANG Yifei (王逸飞)，ZHANG Aili* (张爱丽). Real-Time Calculation Method for Temperature Distribution of Temperature-Controlled Radiofrequency Ablation[J]. J Shanghai Jiaotong Univ Sci, 2023, 28(4): 411-.

References

[1] SCHULLIAN P, PUTZER D, EBERLE G, et al.Simultaneous stereotactic radiofrequency ablation ofmultiple ( 4) liver tumors: Feasibility, safety, and ef-ficacy [J]. Journal of Vascular and Interventional Radiology, 2020, 31(6): 943-952.
[2] TANG R B, WANG C C. Numerical simulation oftissue-equivalent material experiments for radio frequency hyperthermia of tumors [J]. Journal of Tsinghua University (Science and Technology), 2002,42(5): 676-679 (in Chinese).
[3] SAMSET E. Temperature mapping of thermal ablation using MRI [J]. Minimally Invasive Therapy & Allied Technologies, 2006, 15(1): 36-41.
[4] CHEN R D, LU F, WU F, et al. An analytical solution for temperature distributions in hepatic radiofrequency ablation incorporating the heat-sink effect oflarge vessels [J]. Physics in Medicine and Biology,2018, 63(23): 235026.
[5] MOLINA J A L, RIVERA M J, BERJANO E. Analytical transient-time solution for temperature in nonperfused tissue during radiofrequency ablation [J]. Applied Mathematical Modelling, 2017, 42: 618-635.
[6] HAEMMERICH D, CHACHATI L, WRIGHT A S,et al. Hepatic radiofrequency ablation with internallycooled probes: Effect of coolant temperature on lesionsize [J]. IEEE Transactions on Bio-Medical Engineering, 2003, 50(4): 493-500.
[7] LU X W, KIKUCHI H, HIROOKA K, et al. Methodfor estimating the temperature distribution associatedwith the vessel cooling effect in radio frequency ablation [C]//2015 37th Annual International Conferenceof the IEEE Engineering in Medicine and Biology Society. Milan, Italy: IEEE, 2015: 4836-4839.
[8] OOI E H, LEE K W, YAP S, et al. The effects ofelectrical and thermal boundary condition on the simulation of radiofrequency ablation of liver cancer fortumours located near to the liver boundary [J]. Computers in Biology and Medicine, 2019, 106: 12-23.
[9] AUDIGIER C, MANSI T, DELINGETTE H, et al.Efficient lattice Boltzmann solver for patient-specificradiofrequency ablation of hepatic tumors [J]. IEEETransactions on Medical Imaging, 2015, 34(7): 1576-1589.
[10] HE Z Z, LIU J. An efficient thermal evolution modelfor cryoablation with arbitrary multi-cryoprobe configuration [J]. Cryobiology, 2015, 71(2): 318-328.
[11] ZORBAS G, SAMARAS T. Parametric study of radiofrequency ablation in the clinical practice with theuse of two-compartment numerical models [J]. Electromagnetic Biology and Medicine, 2013, 32(2): 236-243.
[12] DOSS J D. Calculation of electric fields in conductivemedia [J]. Medical Physics, 1982, 9(4): 566-573.
[13] PENNES H H. Analysis of tissue and arterial bloodtemperatures in the resting human forearm [J]. Journalof Applied Physiology, 1998, 85(1): 5-34.
[14] ZHENG Y B, ZHANG K W, ZOU J C, et al. Annoninvasive and impedance-ignored control strategyof the ablation zone in radiofrequency ablation therapy [C]//2019 41st Annual International Conferenceof the IEEE Engineering in Medicine and Biology Society. Berlin, Germany: IEEE, 2019: 5514-5517.
[15] TRUJILLO M, BERJANO E. Review of the mathematical functions used to model the temperature dependence of electrical and thermal conductivities of biological tissue in radiofrequency ablation [J]. International Journal of Hyperthermia, 2013, 29(6): 590-597.
[16] CHANG I. Finite element analysis of hepatic radiofrequency ablation probes using temperature-dependentelectrical conductivity [J]. Biomedical Engineering Online, 2003, 2: 12.