介电特性(介电常数与有效电导率)、热物性参数(体积热容与热导率)、组织血流灌注率(blood perfusion rate)：何者是肺部微波消融(microwave ablation, MWA)模型生物物理参数变异性的主导来源？

《International Journal of Hyperthermia》：Dielectric properties, thermal properties, blood perfusion rate: which biophysical parameter is the leading source of variability of microwave ablation models in lung?

　　

摘要：肺部肿瘤微波消融(MWA)术后疗效存在较大个体差异，提示需建立考虑组织特异性生物物理特性的三维温度场量化方法来评估消融区。本研究旨在通过计算生物传热(bioheat transfer)模型，识别肿瘤-肺交界区各生物物理参数对消融区(ablation zone)影响程度的大小排序。研究人员构建2.45 GHz微波施加器置于直径2 cm肿瘤内的计算模型，采用60 W—10 min能量设置，分别取低/高值进行仿真：(1)相对介电常数(relative permittivity, εr) 28.3–50.4及有效电导率(effective conductivity, σ) 1.0–1.9 S·m?1；(2)体积热容(volumetric heat capacity, ρc) 1.3–3.5 MJ·m?3·K?1及热导率(thermal conductivity, k) 0.2–0.4 W·m?1·K?1；(3)肿瘤周边-肺实质血流灌注率(blood perfusion rate, Wbl) 5–9 kg·m?3·s?1。量化60 ℃等温线长度(L60)与直径(D60)，含5 mm边缘(TM5)及不含边缘(TM0)的肿瘤体积热覆盖率(volumetric thermal coverage)。结果显示：体积热容与肿瘤-肺实质边界血流灌注率的变异引起TM5差异5.7–11.4％、D60差异3.1–4.6 mm、L60差异3.0–4.4 mm；介电特性基线值变异仅引起TM5差异约3％、D60与L60差异约1.2 mm。结论：肺部MWA达到充分治疗边缘的过程中，热物性参数及肿瘤周围组织的血流灌注率是影响最显著的参数。该定量分析结果指出，后续实验应重点测定上述参数的组织特异性数值，以为肺部MWA治疗计划(treatment planning)提供循证依据。

论文解读：肺部微波消融模型中生物物理参数对消融区变异性的敏感性分析

《International Journal of Hyperthermia》刊载的此项研究针对影像引导下肺部肿瘤微波消融(microwave ablation, MWA)术后消融体积存在显著个体间差异（文献报道约3–14 cm3），而现行临床治疗计划依赖离体动物组织测得的厂家预设消融区估算、未纳入患者特异性组织生物物理特性这一问题，开展生物传热传递建模(bioheat transfer modeling)的敏感性分析。研究目的为明确肺部MWA仿真模型中，介电特性(dielectric properties：相对介电常数relative permittivity ε_r、有效电导率effective conductivity σ)、热物性参数(thermal properties：体积热容volumetric heat capacity ρc、热导率thermal conductivity k)及肿瘤-肺实质边界血流灌注率(blood perfusion rate, W_bl)三组生物物理参数的基线值变异性，何者对消融区几何尺寸及带5 mm安全缘的体积热覆盖率(volumetric thermal coverage, TM5)影响最大，从而为后续人体肺组织生物物理特性实验测定方向及个体化MWA治疗计划提供量化依据。研究人员采用Pennes生物传热方程(Pennes bioheat equation)耦合亥姆霍兹方程(Helmholtz's equation)求解2.45 GHz时谐电场分布及电磁功率沉积，在轴对称二维旋转几何模型中置入直径20 mm球形肿瘤及外周肺实质域、中心放置6 mm辐射端水冷同轴单极天线，固定60 W输入功率持续10 min，分别单独改变三组参数的低/高文献报道值，提取60 ℃等温线包围的消融区长(L₆₀)与直径(D₆₀)及TM0（肿瘤完全被60 ℃覆盖）、TM5（肿瘤+5 mm边缘被60 ℃覆盖），比较各参数变异引起的指标变化幅度。研究得出结论：体积热容与热导率、肿瘤外围-肺实质血流灌注率差异是肺部MWA消融区尺寸与边缘覆盖率的主要影响因素；介电特性在生理至高温范围内的变异性对最终消融区影响可忽略。该结果提示未来肺部MWA个性化计划模型应优先获取或校准热物性参数与局部灌注信息，介电特性可采用已发表人群均值初始化。

主要关键技术方法：

研究人员建立含直径20 mm球形肿瘤及外围肺实质（半径50 mm、高100 mm）的轴对称几何模型，中心置入辐射段长6 mm的水冷同轴单极微波施加器（2.45 GHz，输入功率60 W，辐照10 min），肿瘤域设为三层以模拟瘤内到周边血流灌注梯度。采用Pennes生物传热方程（含温度依赖介电特性及灌注截断温区56–61 ℃）耦合频域亥姆霍兹方程进行电磁-热双向求解。分别取文献报道的人肺/肝组织低/高值代入三组参数单独做敏感性测试：(1) ε_r=28.3/50.4，σ=1.0/1.9 S·m^?1；(2) ρc=1.3×10⁶/3.5×10⁶J·m^?3·K^?1，k=0.2/0.4 W·m^?1·K^?1；(3)均匀灌注5 kg·m^?3·s^?1及分层模型中外周肺实质灌注9 kg·m^?3·s^?1、瘤周层(layer 3)灌注5 kg·m^?3·s^?1。评价指标为10 min时60 ℃等温面所围D₆₀、L₆₀及TM0、TM5。

研究结果

3.1. What is the influence of dielectric properties?

ε_r取极值差Δ=22.1、σ取Δ=0.9 S·m^?1时，L₆₀与D₆₀仅相差1.2 mm，TM5相差3％；TM0均为100％（2 cm肿瘤可被完全消融），但TM5均低于90％（5 mm边缘未全面达到）。由于升温过程中温度依赖性使ε_r与σ差值缩小（90 ℃时Δε_r≈12，Δσ≈0.6 S·m^?1），介电特性变异经短时初期后对后续热扩散主导的消融区生长几无作用。换用带扼流环天线及50 W—5 min设置验证，影响仍≤2.5 mm。研究人员由此判定介电特性基线值变异性对肺部MWA消融边缘达成的影响可忽略。

3.2. What is the influence of thermal properties?

体积热容ρc（Δ=2.2×10⁶J·m^?3·K^?1）与热导率k（Δ=0.2 W·m^?1·K^?1）取高低值时，D₆₀差5.2 mm，L₆₀差12 mm，TM5差7.9％（约为介电特性影响的2.6倍），TM0均为100％。低热容条件下升温快，组织电导率随温升骤降使电磁能更多沿天线纵向吸收而非径向扩散，且低k限制径向热传导，致消融区长轴比(aspect ratio, AR=D₆₀/L₆₀)由高热容时的~0.9变为~0.7（纵向增长43％，径向仅增17％）。体积热容通过改变升温速率间接调控温变介电特性及功率沉积剖面，是对温度场形态作用最显著的热物性参数。

3.3. What is the influence of blood perfusion rate?

当瘤周层(layer 3)与邻接肺实质灌注相等（均为5 kg·m^?3·s^?1），无论瘤内分层与否（layer 1/2高→低灌注梯度），L₆₀、D₆₀、TM0、TM5同均质灌注模型无显著差异——瘤中心区因快速升至＞56 ℃血管凝固，内部灌注异质性不影响消融边界。当肺实质灌注(9 kg·m^?3·s^?1)高于瘤周(layer 3, 5 kg·m^?3·s^?1)，D₆₀缩至25.7–29.2 mm（原28.8–33.8 mm），L₆₀缩至30.3–41.0 mm（原33.3–45.4 mm）；TM5在高/低热物性下分别下降11.4％/5.7％，即肿瘤-肺交界灌注差是边界冷却主因。温度-体积直方图显示此情形下TM5在60 ℃阈值处降至72％（高热物性）及85％（低热物性），提示较高肺实质灌注增加残存活灶风险。

讨论与结论（翻译自原文Conclusion）：

本研究证明：(1) 热物性参数及肿瘤边缘处的血流灌注率是肺部MWA消融区范围、形态及最终体积热覆盖率的首要影响生物物理参数；(2) 介电特性基线值的变异性对所估算含充分消融边缘的体积热覆盖率影响可忽略；(3) 仅依靠消融区长与直径评估肺部MWA成功与否可能不足，借助计算模型实现的体积热剂量学(volumetric thermal dosimetry)分析可为MWA治疗计划提供更有价值信息。