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

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

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


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

技术介绍

[0002]在射频消融手术中，常常使用仿真方法对消融区域及温度场分布进行估计，需要对电热耦合场进行迭代计算，目前最常用的有限元算法具有耗时高、内存占用高等缺点。
[0003]专利文献为CN111027225A的专利技术专利公开了一种基于数据拟合的消融仿真方法，该方法包括：计算同一消融功率、消融时间下的消融边界的短径、长径数据的实测平均值，并对靶区进行热场仿真，以短径数据为基准，将温度阈值所在等高线对应的长径为长径仿真值；在低消融功率下，建立温度阈值与消融功率关系模型；在高消融功率下，建立温度阈值与消融时间关系模型；建立修正系数与消融功率关系模型；输入消融功率、消融时间，根据消融功率高低选择温度阈值与消融功率关系模型或温度阈值与消融时间关系模型，得到温度阈值，根据温度阈值得到对应的消融仿真图像，利用修正系数与消融功率关系模型，对消融仿真图像的长径进行修正。但是上述方案基于经验数据统计，缺乏对个体化差异的适应；温度阈值难以确定，无法与真实消融过程建立直接联系。
[0004]唐茹冰发表的《射频加热治疗肿瘤体模实验的数值模拟》一文中提出的肿瘤体模的有限元稳态热仿真方法，虽然能够计算稳态导热，但是对于非稳态以及非射频加热的场景缺乏适应性。Mohammad H.N.Naraghi发表的《A BOUNDARY
‑
DISPATCH MONTE CARLO(EXODUS)METHOD FOR ANALYSIS OF CONDUCTIVE HEAT TRANSFER PROBLEMS》一文中提出的边界派发式的蒙特卡洛热仿真方法，对于边界不规则的情况适应程度不佳，无法处理内部热源点。

技术实现思路

[0005]针对现有技术中的缺陷，本专利技术的目的是提供一种基于蒙特卡洛法的射频消融快速仿真方法及系统。
[0006]根据本专利技术提供的一种基于蒙特卡洛法的射频消融快速仿真方法，包括如下步骤：
[0007]基础仿真模型构建步骤：输入CT/MR形式的影像数据，根据影像数据生成体素网格，构建基础仿真模型；
[0008]参数估计步骤：输入基础仿真模型，基于影像数据的分割结果，细化模型边界，对边界体素进行细分，对于基础仿真模型计算域中的每一处，根据影像数据的数值，估算差异化参数；
[0009]计算域约化步骤：基于多维度条件对计算域进行聚类约化；
[0010]电/热物理场耦合计算步骤：使用蒙特卡洛法对电物理场耦合计算进行简化，使用有限差分法对热物理场进行计算，迭代进行快速的热/电物理场耦合计算；
[0011]细胞凋亡计算步骤：计算细胞凋亡结果，判断是否满足停止条件，若判断结果为是，则输出消融区域；若判断结果为否，则跳转至计算域约化步骤。
[0012]优选地，每个体素单元之间使用Pennes生物导热方程的差分形式确定热流传导方向和关系。
[0013]优选地，所述计算域约化步骤包括：
[0014]基于影像数据的分割结果确定组织类型，根据组织边界调整计算域，将目标组织边界外的部位排除出计算域；
[0015]根据消融探针位置，将超出距离阈值以及不在手术范围内的区域排除出计算域；
[0016]依据仿真中间结果，对相变、凋亡等动态过程造成的参数变化进行响应；
[0017]对已经变性、烧焦的区域调整模型参数，根据对称性条件，对模型计算域中电/热属性以及距离构成的计算条件相似的位置进行聚类约化。
[0018]优选地，所述差异化参数包括差异化的导热率、电阻率、血流灌注率。
[0019]优选地，使用蒙特卡洛法对电物理场耦合计算包括：
[0020]步骤S1：确定射频发射位置为起始位置，根据预设参数生成数量为N的模拟粒子，有比例为T
b
/T
h
的粒子为高能粒子，其余为低能粒子，其中T
b
代表加热的目标温度，T
h
代表组织细胞的最高温度；
[0021]步骤S2：计算粒子的随机游走，每运行一个时间步长，模拟粒子根据模型的电传导可行方向游走M步，直到到达边界停止；
[0022]步骤S3：所有模拟粒子随机游走结束后，统计每个体素网格内的经过的高能粒子数和低能粒子数，网格电场强度估计为E
i
+N
h
/N
×
Δe，其中E
i
为原始场强数值，Δe为电场强度单位值；
[0023]步骤S4：根据焦耳热计算公式Q
j
＝J
·
E计算射频电场在一个时间步长内的产生的热能，根据pennes方程
[0024][0025]在当前网格模型中的差分形式，使用有限差分法求解传热方程，将结果反馈到温度场中。
[0026]优选地，采用CUDA对粒子游走过程进行了GPU加速。
[0027]优选地，使用CUDA对传热计算进行并行化硬件加速。
[0028]优选地，所述细胞凋亡计算步骤中，根据温度场分布，依照阿累尼乌斯方程计算细胞累积收到的热损伤
[0029][0030]其中A为消融损伤系数，ΔE为激活能量，R为通用气体常数，T为温度，当Ω超过阈值时，将细胞设置为凋亡状态，计算由消融热效应造成的组织细胞凋亡结果，从而推算消融区域大小。
[0031]优选地，使用CUDA对损伤计算进行并行化硬件加速。
[0032]根据本专利技术提供的一种基于蒙特卡洛法的射频消融快速仿真系统，包括如下模块：
[0033]基础仿真模型构建模块：输入CT/MR形式的影像数据，根据影像数据生成体素网格，构建基础仿真模型；
[0034]参数估计模块：输入基础仿真模型，基于影像数据的分割结果，细化模型边界，对边界体素进行细分，对于基础仿真模型计算域中的每一处，根据影像数据的数值，估算差异化参数；
[0035]计算域约化模块：基于多维度条件对计算域进行聚类约化；
[0036]电/热物理场耦合计算模块：使用蒙特卡洛法对电物理场耦合计算进行简化，使用有限差分法对热物理场进行计算，迭代进行快速的热/电物理场耦合计算；
[0037]细胞凋亡计算模块：计算细胞凋亡结果，判断是否满足停止条件，若判断结果为是，则输出消融区域；若判断结果为否，则跳转至计算域约化模块。
[0038]与现有技术相比，本专利技术具有如下的有益效果：
[0039]1、本专利技术通过参数估计，实现了高低精度混合的体素网格建模，缓解了传统建模方式导致的计算量冗余的问题。
[0040]2、本专利技术通过计算域约化，实现了对具有对称性、等价性位置的自动归约，自动化地将不需要计算的区域排除出计算域，减少了计算量。
[0041]3、本专利技术通过蒙特卡洛法，对电场计算进行了合理简化，结合GPU硬件加速，实现了快速的电场计算。
...

【技术保护点】

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

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