在侵入性激光热治疗过程中,生物组织中温度场的在线解耦估计

《Journal of Theoretical Biology》:Online decoupled estimation of temperature field in biological tissue during invasive laser thermal therapy

【字体: 时间:2026年02月11日 来源:Journal of Theoretical Biology 2

编辑推荐:

  针对侵入式激光诱导热疗(LITT)中温度场实时估计难题,提出基于解耦估计的在线优化模型,结合时空相关模型与最优估计框架,实现关键点温度解耦估计,有效克服传统状态空间模型的耦合问题,并通过实验与仿真验证其计算效率与准确性,为临床LITT提供可靠技术支持。

  
毛照辉|陈宏|王广军|李振彪
重庆大学能源与动力工程学院,中国重庆400044

摘要

侵入式激光诱导热疗(LITT)是一种用于治疗深部肿瘤的新兴微创技术。在线估计目标生物组织中的瞬态温度场对于LITT技术的研究与开发至关重要。现有的状态估计理论采用状态空间模型来描述系统,本质上属于耦合估计方法。这些方法无法实现系统状态的解耦估计,因此在将当前状态估计方法应用于LITT过程中的温度场估计时面临诸多挑战。本文重点研究了侵入式LITT过程中生物组织中温度场的在线解耦估计。通过构建最优估计模型和时空相关模型,实现了整个侵入式LITT过程中关键离散点的温度在线解耦估计。该方案有效克服了基于状态空间模型的现有温度场估计技术所固有的困难。利用现有实验数据验证了所构建生物模型的可靠性。通过数值模拟分析了影响时间域、测量噪声和模型不匹配的影响。与标准卡尔曼滤波器的比较表明,当只需要实时估计特定区域(如肿瘤区域)的温度时,本文提出的解耦估计方案在计算效率方面具有显著优势。这项研究为组织温度场的在线估计提供了可靠的技术支持。

引言

近年来,侵入式激光诱导热疗(LITT)作为一种微创技术,在去除深部恶性肿瘤方面受到了广泛关注(Vitorino等人,2023年;Wang等人,2023年)。其核心原理是利用激光产生的热效应加热深部肿瘤组织,诱导凝固性坏死,最终实现肿瘤消融的治疗目标(Yin等人,2020年)。
由于正常细胞和肿瘤细胞之间的致死温度阈值存在显著差异,LITT过程对温度场有严格的要求(Obonai等人,2024年)。一方面,肿瘤组织温度必须升高到43°C以上才能确保有效消融;另一方面,正常组织温度必须保持在41°C以下以避免不可逆损伤(Pang等人,2025年)。在治疗过程中实时监测温度场对于确保LITT的治疗效果和安全性至关重要(Nomura等人,2017年)。
目前,深部肿瘤的侵入式LITT过程中的温度监测主要采用侵入式或非侵入式方法(Liu等人,2025年;Pominova等人,2019年)。
侵入式方法使用热电偶或光纤布拉格光栅等点式传感器进行局部测量(Saccomandi等人,2013年)。尽管这些传感器体积小且响应速度快,但它们只能提供传感器附近的温度数据,无法捕捉肿瘤及其周围组织的完整三维温度分布。因此,侵入式单点监测无法完全满足临床对全面、实时热反馈的需求(Lo等人,2017年)。
非侵入式技术,如磁共振温度成像(MRTI)、光声成像和超声成像,可以监测较大区域,但存在局限性。MRTI成本高昂、速度较慢,且在检测温度信号方面存在挑战;而光声成像和超声成像受限于组织穿透深度,无法用于深部肿瘤的温度监测(Alaeian等人,2019年;Cui等人,2020年)。
为了解决这些局限性,通常会使用卡尔曼滤波器和Luenberger观测器等状态估计技术(Kalman,1960年;Luenberger,1971年)。这些方法利用可观测的温度数据实时预测未知的内部组织温度场,从而实现LITT过程中的在线估计。例如,Pacheco等人应用稳态卡尔曼滤波器对LITT过程中的深层组织进行在线二维温度估计(Pacheco等人,2020年)。Wang等人使用卡尔曼滤波器进行浅表肿瘤热疗的实时温度场估计,该方法后来被扩展到深层热疗场景(Li等人,2025年;Wang等人,2024年)。同样,Hamza等人利用Luenberger观测器对皮下有机组织进行实时温度场估计(El-Kebir等人,2022年)。
现有的状态估计方法使用状态空间模型来描述系统。这需要离散化生物组织中的热传递控制方程以构建相应的状态空间模型。这些方程包括经典的傅里叶和非傅里叶热传递模型,如热波模型和相位滞后模型(Azhdari等人,2023年;Ghanbari和Rezazadeh,2022年;Seyedpour等人,2024年)。通过引入状态观测器,可以实现内部状态的递归最优估计。然而,将这些方法应用于实际的侵入式LITT治疗时仍存在问题。
首先,大多数临床遇到的恶性肿瘤位于边界复杂、形状不规则的深部组织中。这种复杂性使得生物热传递方程的离散化和准确状态空间模型的构建变得复杂,直接影响了温度场估计的准确性(Pourkargar和Armaou,2015年;Singh和Kumar,2023年)。
其次,这些方法通过利用所有系统状态之间的联系以及输入(激光)与状态(温度)之间的关系来在线估计系统状态。它们本质上属于一种状态耦合估计方案。无论实际治疗是否需要监测整个区域,这些方案都需要估计整个生物组织的状态空间(温度场)。这不仅显著增加了不必要的计算成本,还可能严重影响实时性能。然而,在侵入式LITT过程中,实时性能是一个关键的安全指标(Pham等人,2017年)。
最后,由于需要同时估计生物组织的整个温度场,上述状态耦合估计方案通常需要足够丰富的测量信息才能获得满意的结果。然而,在临床实践中,深部肿瘤的LITT过程中的温度监测通常仅依赖于通过光纤传感器获得的单点温度测量数据。测量信息的匮乏严重限制了这些方案的可行性(Melnick等人,2021年)。
基于系统状态空间模型的状态估计及其固有的状态耦合估计特性是上述困难的根本原因。因此,研究侵入式热疗过程中生物组织温度场的解耦估计问题具有重要的理论和实际意义。
现有的在线解耦框架主要面向多变量控制领域。它们通常通过相对增益阵列分析或传递函数对角化等方法在有限维系统中实现静态或动态解耦,不适用于具有无限维分布参数特性的生物热疗温度场估计(Yang等人,2019年)。
Wang等人(Chen等人,2023年;Mao等人,2025年)在初步研究中发现,在热传递系统中的任意两个空间点之间的温度响应之间存在确定性相关机制。这种机制与内部热源无关,可以通过建立温度响应的时空相关模型直接重建系统的瞬态温度场。基于此,可以使用阶跃响应模型构建热传递系统的解耦状态参数估计模型。因此,可以利用可获取的温度信息解耦估计系统的瞬态温度场。
本研究基于解耦估计技术,解决了侵入式LITT过程中生物组织中温度场的在线估计问题。通过将Penes生物热方程与Beer–Lambert定律结合,开发了一个基础的热模型。基于该模型,建立了温度阶跃响应模型和时空相关模型。利用带有过程噪声的阶跃响应模型,最优估计框架可以过滤噪声测量数据,得到修正后的值。然后将这些值与时空相关模型结合,实现治疗过程中离散空间点的在线解耦温度估计。这项工作为侵入式热疗过程中的实时估计能力提供了科学验证。

章节片段

生物组织热传递模型

在侵入式LITT中,光纤将红外激光能量传递到组织深处以热消融肿瘤(图1a)。对于热传导建模,本文采用了基于傅里叶定律的经典Penes生物热方程(Pennes,1948年)。尽管理论上也探讨了非傅里叶模型,但最近的实验研究表明,基于傅里叶的模型已经足以满足LITT等临床应用的精度要求(Azhdari等人,2025a,2026年;Seyedpour等人,2024年)。

生物组织热传递模型的验证

为了验证所提出的生物组织热传递模型的可靠性,将仿真结果与参考文献中的实验数据进行了对比(Bianchi等人,2021年)。
数值模拟使用商业软件进行,时间离散化采用隐式后向差分方案。模拟中使用的生物组织的光学和热参数详细列在表S1中。

局限性与实际应用

应当注意的是,所提出方案的性能取决于参数的适当选择,例如影响时间域D和正则化参数γ。在实际应用中,必须调整参数以平衡估计精度和计算负担。未来的研究应探索自适应参数选择算法。此外,如果实际组织参数与离线模型之间存在显著的时间变化偏差,当前的线性

结论

基于解耦估计模型,本文研究了侵入式LITT过程中组织温度场的在线估计。主要结论如下:
  • (1)
    基于Penes生物热方程建立了侵入式LITT的生物组织热传递模型。首次将分布式参数系统的解耦估计理论引入生物热疗,创建了一种依赖于单点测量的估计方案,无需复杂
  • CRediT作者贡献声明

    王广军:监督、资源管理、项目协调、资金获取。李振彪:撰写——审稿与编辑、可视化、验证、研究。毛照辉:撰写——初稿、方法论、研究、数据分析、概念化。陈宏:撰写——审稿与编辑、监督、资源协调

    利益冲突声明

    我们声明没有利益冲突。

    致谢

    本工作得到了国家自然科学基金(项目编号52276051)的支持。
    相关新闻
    生物通微信公众号
    微信
    新浪微博

    热点排行

      今日动态 | 人才市场 | 新技术专栏 | 中国科学人 | 云展台 | BioHot | 云讲堂直播 | 会展中心 | 特价专栏 | 技术快讯 | 免费试用

      版权所有 生物通

      Copyright© eBiotrade.com, All Rights Reserved

      联系信箱:

      粤ICP备09063491号