基于蒙特卡洛法的射频消融快速仿真方法及系统技术方案

本发明专利技术提供了一种基于蒙特卡洛法的射频消融快速仿真方法及系统，根据影像信息对计算域进行聚类预处理，构建约化仿真模型，对电物理场耦合计算进行合理简化，使用有限差分法对热物理场进行计算，迭代进行快速的热/电物理场耦合计算；基于CUDA使用GPU硬件加速仿真计算过程中的模型构建和数值计算过程。本发明专利技术实现了高低精度混合的体素网格建模，缓解了传统建模方式导致的计算量冗余的问题；实现了对具有对称性、等价性位置的自动归约，自动化地将不需要计算的区域排除出计算域，减少了计算量；对电场计算进行了合理简化，结合GPU硬件加速，实现了快速的电场计算；实现了电/热耦合场的加速运算，提供了实时的手术区域温度场计算。算。算。

【技术实现步骤摘要】
基于蒙特卡洛法的射频消融快速仿真方法及系统


[0001]本专利技术涉及医疗领域，具体地，涉及一种基于蒙特卡洛法的射频消融快速仿真方法及系统。

技术介绍

[0002]在射频消融手术中，常常使用仿真方法对消融区域及温度场分布进行估计，需要对电热耦合场进行迭代计算，目前最常用的有限元算法具有耗时高、内存占用高等缺点。
[0003]专利文献为CN111027225A的专利技术专利公开了一种基于数据拟合的消融仿真方法，该方法包括：计算同一消融功率、消融时间下的消融边界的短径、长径数据的实测平均值，并对靶区进行热场仿真，以短径数据为基准，将温度阈值所在等高线对应的长径为长径仿真值；在低消融功率下，建立温度阈值与消融功率关系模型；在高消融功率下，建立温度阈值与消融时间关系模型；建立修正系数与消融功率关系模型；输入消融功率、消融时间，根据消融功率高低选择温度阈值与消融功率关系模型或温度阈值与消融时间关系模型，得到温度阈值，根据温度阈值得到对应的消融仿真图像，利用修正系数与消融功率关系模型，对消融仿真图像的长径进行修正。但是上述方案基于经验数据统计，缺乏本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于蒙特卡洛法的射频消融快速仿真方法，其特征在于，包括如下步骤：基础仿真模型构建步骤：输入CT/MR形式的影像数据，根据影像数据生成体素网格，构建基础仿真模型；参数估计步骤：输入基础仿真模型，基于影像数据的分割结果，细化模型边界，对边界体素进行细分，对于基础仿真模型计算域中的每一处，根据影像数据的数值，估算差异化参数；计算域约化步骤：基于多维度条件对计算域进行聚类约化；电/热物理场耦合计算步骤：使用蒙特卡洛法对电物理场耦合计算进行简化，使用有限差分法对热物理场进行计算，迭代进行快速的热/电物理场耦合计算；细胞凋亡计算步骤：计算细胞凋亡结果，判断是否满足停止条件，若判断结果为是，则输出消融区域；若判断结果为否，则跳转至计算域约化步骤。2.根据权利要求1所述的基于蒙特卡洛法的射频消融快速仿真方法，其特征在于，每个体素单元之间使用Pennes生物导热方程的差分形式确定热流传导方向和关系。3.根据权利要求1所述的基于蒙特卡洛法的射频消融快速仿真方法，其特征在于，所述计算域约化步骤包括：基于影像数据的分割结果确定组织类型，根据组织边界调整计算域，将目标组织边界外的部位排除出计算域；根据消融探针位置，将超出距离阈值以及不在手术范围内的区域排除出计算域；依据仿真中间结果，对相变、凋亡等动态过程造成的参数变化进行响应；对已经变性、烧焦的区域调整模型参数，根据对称性条件，对模型计算域中电/热属性以及距离构成的计算条件相似的位置进行聚类约化。4.根据权利要求1所述的基于蒙特卡洛法的射频消融快速仿真方法，其特征在于，所述差异化参数包括差异化的导热率、电阻率、血流灌注率。5.根据权利要求1所述的基于蒙特卡洛法的射频消融快速仿真方法，其特征在于，使用蒙特卡洛法对电物理场耦合计算包括：步骤S1：确定射频发射位置为起始位置，根据预设参数生成数量为N的模拟粒子，有比例为T
b
/T
h
的粒子为高能粒子，其余为低能粒子，其中T
b
代表加热的目标温度，T
h
代表组织细胞的最高温度；步骤S2：计算粒子的随机游走，每运行一个时间步长，模拟粒子根据模型的电传导可行方向游走M步，直到到...

【专利技术属性】
技术研发人员：王皓宇，
申请(专利权)人：精劢医疗科技南通有限公司，
类型：发明
国别省市：