CN110263489A - 基于dicom数据的肝肿瘤微波消融三维温度场仿真方法 - Google Patents

基于dicom数据的肝肿瘤微波消融三维温度场仿真方法 Download PDF

Info

Publication number
CN110263489A
CN110263489A CN201910598064.9A CN201910598064A CN110263489A CN 110263489 A CN110263489 A CN 110263489A CN 201910598064 A CN201910598064 A CN 201910598064A CN 110263489 A CN110263489 A CN 110263489A
Authority
CN
China
Prior art keywords
tumour
liver
ablation
emulation
tissue
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201910598064.9A
Other languages
English (en)
Inventor
钱志余
王娟
赵金哲
沐勇杰
晋晓飞
冯宇
俞钦栋
刘伟清
李韪韬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201910598064.9A priority Critical patent/CN110263489A/zh
Publication of CN110263489A publication Critical patent/CN110263489A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Surgical Instruments (AREA)

Abstract

本发明公开一种基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法,包括如下步骤:步骤1,基于DICOM数据重建并修复肿瘤和血管模型,将其导入多物理场仿真模块中,选定工作平面画出消融针结构图,构建仿真几何模型;步骤2,对仿真几何模型中的计算域设定不同材料及参数,其具体材料种类分为:肝脏,肿瘤,PTFE;步骤3,构建耦合电磁波辐射模型及生物传热模型;步骤4,根据步骤3构建的模型,依据设定的微波频率对仿真几何模型划分网格,设定求解器求解方法;步骤5,对求解得到的仿真数据进行可视化处理,得到温度场分布数据。此种仿真方法可利用病人术前影像数据,建立三维有限元仿真模型,计算消融热场,指导手术方案的制定。

Description

基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法
技术领域
本发明属于有限元仿真领域,特别涉及一种基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法。
背景技术
肝肿瘤是一种常见的腹部恶性肿瘤,微波消融作为针对该肿瘤的一种有效治疗方法已被纳入NCCN(美国国立综合癌症网络)治疗指南中。其术前方案制定涉及消融热剂量(消融功率、消融时间)设定以及消融针穿刺路径规划。此两点规划有赖于离体消融模型的体积实验数据和临床医生穿刺消融针的经验,缺乏针对不同病人的精准手术方案以及针对穿刺经验较少的临床医生的指导。
现有的肝肿瘤微波消融仿真方法多是基于二维轴对称几何模型,忽略在体消融时血管以及血流灌注的影响,计算得出的是完全对称的消融热场。因此如何在术前根据不同病人建立不同的消融模型,充分模拟出血管以及在体消融时血流灌注对消融结果的影响,调整热剂量以及进针方式,观察仿真消融热场分布,对手术计划的制定具有指导意义。
发明内容
本发明的目的,在于提供一种基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法,其可利用病人术前影像数据,建立三维有限元仿真模型,计算消融热场,指导手术方案的制定。
为了达成上述目的,本发明的解决方案是:
一种基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法,包括如下步骤:
步骤1,基于DICOM数据重建并修复肿瘤和血管模型,将其导入多物理场仿真模块中,选定工作平面画出消融针结构图,构建仿真几何模型;
步骤2,对仿真几何模型中的计算域设定不同材料及参数,其具体材料种类分为:肝脏,肿瘤,PTFE;
步骤3,构建耦合电磁波辐射模型及生物传热模型;
步骤4,根据步骤3构建的模型,依据设定的微波频率对仿真几何模型划分网格,设定求解器求解方法;
步骤5,对求解得到的仿真数据进行可视化处理,得到温度场分布数据。
采用上述方案后,本发明可以解决现有二维轴对称仿真方法中无法真实反映不同病人在相同热剂量下消融疗效差异的问题,能够对不同病人建立个性化三维有限元仿真模型,并采用动态生物材料参数计算得出温度场分布,得出与真实消融过程更为一致的预测模型,可仿真消融针不同穿刺路径的消融效果用以指导手术计划,有助于手术方案的制定。
附图说明
图1是本发明的流程图;
图2是本发明构建的几何模型图;
图3是本发明有限元模型图;
图4是本发明仿真结果,显示为50W,600s消融热剂量下肿瘤表面温度分布图;
图5(a)是本发明仿真结果,显示为50W,300s消融热剂量下,肿瘤切面温度场分布图;图5(b)是本发明仿真结果,显示为50W,300s消融热剂量下,肿瘤切面温度场及其60℃等温线范围图;
图6(a)是本发明仿真结果,显示为50W,600s消融热剂量下,肿瘤切面温度场分布图;图6(b)是本发明仿真结果,显示为50W,600s消融热剂量下,肿瘤切面温度场及其60℃等温线范围图;
图7是本发明在距离肿瘤及消融针不同位置处设定的4个测温点示意图;
图8是本发明在距离肿瘤及消融针不同位置处设定的4个测温点在50W消融功率下,600s内的温升曲线图。
具体实施方式
以下将结合附图,对本发明的技术方案及有益效果进行详细说明。
如图1所示,本发明提供一种基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法,以腹部CT,临近血管处肝肿瘤模型为例,包括以下步骤:
(一)在多物理场仿真模块中导入经DICOM数据重建并修复后的肿瘤和血管模型,选定工作平面画出消融针结构图,构建仿真几何模型;
1)基于DICOM数据重建肿瘤和血管的三维模型;
2)重建后的模型经过逆向建模模块Geomagic Studio修复,得出表面光滑且网格质量高的三维肿瘤和血管模型;
3)将肿瘤和血管模型导入多物理场仿真模块COMSOL Multiphsics中,设定相对简化容差为0.1,进一步修复网格并且生成实体;
4)在工作平面上建立二维轴对称消融针模型,旋转得出三维消融针几何模型,消融针内外导体简化为理想电导体边界,保留PTFE(聚四氟乙烯)介质为计算域;
5)以球形表征肝脏,本实施例设定球形半径为40mm,构建的几何模型如图2所示;
(二)对模型中计算域设定不同材料及参数,其具体材料种类分为:肝脏,肿瘤,PTFE(聚四氟乙烯);
1)设定肝脏/肿瘤的比热容CP,热导率k,密度ρ分别为:
生物组织相对介电常数ε,电导率σ随温度变化规律符合以下公式(记为f(T)),将肝脏和肿瘤的参数分别记为εliver,σliver,εtumor,σtumor,系数如表1所示:
表1
2)设定PTFE(聚四氟乙烯)介质参数为:εPTFE=2,σPTFE=0;
(三)耦合电磁波辐射模型与生物传热模型;
1)肝组织及肿瘤内平面横向磁场波动方程为:
其中,▽用于求解向量的梯度,μr为相对磁导率,E为电场强度,k0为自由空间波数,ε为组织的相对介电常数,σ为组织电导率,ω为角频率,ε0为真空介电常数。
2)将消融针内外导体简化为理想电导体边界;
3)消融针顶端设定为同轴端口,并定义输入能量大小:2450MHz微波频率下50W消融功率;
4)肝脏及血管外表面设定为电磁波散射边界;
5)肝组织及肿瘤内生物传热方程为:
其中,ρ为组织密度,Cp为组织比热容,ωb为血流灌注率,k为组织热导率,Q为外部产热,Qbio为组织新陈代谢产热,T代表温度。
6)生物组织对微波吸收由组织比吸收率SAR计算得出:
7)生物组织新陈代谢产热记为0,电磁产热记为:
8)定义消融针外表面为20℃,血管外表面为37℃,以分别模拟消融针水冷和血管热沉积效应;
9)对肝脏和肿瘤分别设定血流灌注率为0.016s-1,0.005s-1
(四)对几何模型划分网格,设定求解器求解方法;
1)耦合电磁波辐射模型与生物传热模型后,依据设定的微波频率划分有限单元网格,通常最大的单元网格不超过波长的0.2倍,可根据计算要求做出适当调整。本实施例中对所有域设定网格大小为Fine,最大单元尺寸与最小单元尺寸分别为6.8mm,0.85mm,观察网格划分结果,若存在低质量网络,可局部再细分,构建的有限元模型如图3所示;
2)添加时域研究步骤,600s内步长10s,以直接求解器的方法求解;
(五)对仿真数据可视化处理,得出温度场分布数据;
1)肿瘤外表面温度场分布;如图4所示,下方临近血管处,表面温度低于60℃,表明消融不完全,需要调整热剂量或者插针方式;
2)肿瘤多平面内观察温度场分布;图5(a)显示的是300s消融时长下一切面的温度分布图,图5(b)所示含有60℃等温线的温度分布图,可见等温线(外圈黑线)未超过肿瘤边界(内圈黑线)10mm,表明300s的消融时长不足以覆盖此肿瘤范围,需延长消融时间;图6(a)显示的是600s消融时长下一切面的温度分布图,图6(b)所示600s的消融时长下60℃等温线除临近血管处较为狭窄,其余均覆盖整个肿瘤边界,消融效果较好。可进一步调整切面位置观察温度场分布;
3)设定距离肿瘤和消融针不同距离的测温点观察消融效果,图7为设定的4个测温点示意图,图8为测温点温升曲线,可见2,3点距离消融针较近,最终温度超过60℃,达到消融效果,而1,4点由于收到血管热沉积效应的影响,温升曲线较慢,且均未超过60℃,表明手术方案需做调整。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。

Claims (5)

1.一种基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法,其特征在于包括如下步骤:
步骤1,基于DICOM数据重建并修复肿瘤和血管模型,将其导入多物理场仿真模块中,选定工作平面画出消融针结构图,构建仿真几何模型;
步骤2,对仿真几何模型中的计算域设定不同材料及参数,其具体材料种类分为:肝脏,肿瘤,PTFE;
步骤3,构建耦合电磁波辐射模型及生物传热模型;
步骤4,根据步骤3构建的模型,依据设定的微波频率对仿真几何模型划分网格,设定求解器求解方法;
步骤5,对求解得到的仿真数据进行可视化处理,得到温度场分布数据。
2.如权利要求1所述的基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法,其特征在于:所述步骤1的具体过程是:
步骤11,基于DICOM数据重建肿瘤和血管的三维模型;
步骤12,对步骤11重建的三维模型经过逆向建模模块Geomagic Studio修复,得出表面光滑且网格质量高的三维肿瘤和血管模型;
步骤13,将肿瘤和血管模型导入多物理场仿真模块COMSOL Multiphsics中,设定相对简化容差,进一步修复网格并生成实体;
步骤14,在工作平面上建立二维轴对称消融针模型,旋转得出三维消融针几何模型,消融针内外导体简化为理想电导体边界,保留PTFE介质为计算域;
步骤15,以球形表征肝脏,构建几何模型。
3.如权利要求1所述的基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法,其特征在于:所述步骤2中,各参数设定如下:
设定肝脏/肿瘤的比热容CP,热导率k,密度ρ分别为:
生物组织相对介电常数ε,电导率σ随温度变化规律符合公式f(T),将肝脏和肿瘤的参数分别记为εliver,σliver,εtumor,σtumor,系数如表1所示:
表1
设定PTFE介质参数为:εPTFE=2,σPTFE=0。
4.如权利要求1所述的基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法,其特征在于:所述步骤3的具体过程是:
步骤31,肝组织及肿瘤内平面横向磁场波动方程为:
其中,用于求解向量的梯度,μr为相对磁导率,E为电场强度,k0为自由空间波数,ε为组织的相对介电常数,σ为组织电导率,ω为角频率,ε0为真空介电常数;
步骤32,将消融针内外导体简化为理想电导体边界;消融针顶端设定为同轴端口,并定义输入能量大小;肝脏及血管外表面设定为电磁波散射边界;
步骤33,肝组织及肿瘤内生物传热方程为:
其中,ρ为组织密度,Cp为组织比热容,ωb为血流灌注率,k为组织热导率,Q为外部产热,Qbio为组织新陈代谢产热,T代表温度;
步骤34,生物组织对微波吸收由组织比吸收率SAR计算得出:
步骤35,生物组织新陈代谢产热记为0,电磁产热记为:
5.如权利要求1所述的基于DICOM数据的肝肿瘤微波消融三维温度场仿真方法,其特征在于:所述步骤4的具体过程是:
步骤41,依据耦合电磁波辐射模型与生物传热模型中设定的微波频率划分有限单元网格,观察网格划分结果,若存在低质量网络,则进行局部再细分,最终得到构建的有限元模型;
步骤42,进行时域研究,以直接求解器的方法求解。
CN201910598064.9A 2019-07-04 2019-07-04 基于dicom数据的肝肿瘤微波消融三维温度场仿真方法 Pending CN110263489A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910598064.9A CN110263489A (zh) 2019-07-04 2019-07-04 基于dicom数据的肝肿瘤微波消融三维温度场仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910598064.9A CN110263489A (zh) 2019-07-04 2019-07-04 基于dicom数据的肝肿瘤微波消融三维温度场仿真方法

Publications (1)

Publication Number Publication Date
CN110263489A true CN110263489A (zh) 2019-09-20

Family

ID=67924340

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910598064.9A Pending CN110263489A (zh) 2019-07-04 2019-07-04 基于dicom数据的肝肿瘤微波消融三维温度场仿真方法

Country Status (1)

Country Link
CN (1) CN110263489A (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110882056A (zh) * 2019-12-11 2020-03-17 南京亿高微波系统工程有限公司 Ct下精准肿瘤微波消融系统
CN110942832A (zh) * 2019-12-11 2020-03-31 南京亿高微波系统工程有限公司 一种微波消融系统
CN110991014A (zh) * 2019-11-18 2020-04-10 中国计量大学 人工心脏无线供能系统的多物理场分析方法
CN111027225A (zh) * 2019-12-25 2020-04-17 南京亿高微波系统工程有限公司 一种基于数据拟合的消融仿真方法
CN111329454A (zh) * 2020-04-04 2020-06-26 吉林大学 基于人体三维温度场无损重构技术低温手术辅助诊疗系统
CN111528800A (zh) * 2020-04-30 2020-08-14 深圳开立生物医疗科技股份有限公司 一种肿瘤消融疗效预测方法、装置、设备及计算机介质
CN111653363A (zh) * 2020-06-23 2020-09-11 南京诺源医疗器械有限公司 一种基于仿真温度场的微波消融术中疗效评估方法
CN111685870A (zh) * 2020-06-23 2020-09-22 南京诺源医疗器械有限公司 一种用于肝脏的微波消融碳化调控方法
CN112562859A (zh) * 2020-10-20 2021-03-26 应葵 肿瘤热消融手术智能仿真模型训练系统、训练方法
CN113768484A (zh) * 2020-06-10 2021-12-10 上海美杰医疗科技有限公司 生物组织个性化血流灌注率评估方法及系统
CN113842210A (zh) * 2021-08-02 2021-12-28 应葵 椎骨肿瘤微波消融手术仿真方法及装置
CN116258024A (zh) * 2023-05-16 2023-06-13 赛福凯尔(绍兴)医疗科技有限公司 目标区域的电场可视化方法、装置和计算机设备
CN116705329A (zh) * 2023-06-08 2023-09-05 上海睿速创生医疗科技有限公司 一种基于COMSOL Muliphysics的消融针消融肝脏组织模拟方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107526928A (zh) * 2017-08-21 2017-12-29 北京工业大学 一种基于特性参数反馈的微波热消融温度场仿真方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107526928A (zh) * 2017-08-21 2017-12-29 北京工业大学 一种基于特性参数反馈的微波热消融温度场仿真方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
王娟 等: "基于有限元仿真的微波消融针设计方法的研究", 《生物医学工程研究》 *
赵金哲 等: "一种实用微波消融仿真模型的构建方法", 《中国医疗设备》 *
郑龙波: "射频消融技术在骨肿瘤治疗中的应用研究", 《中国优秀硕士学位论文全文数据库》 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110991014A (zh) * 2019-11-18 2020-04-10 中国计量大学 人工心脏无线供能系统的多物理场分析方法
CN110991014B (zh) * 2019-11-18 2023-06-16 中国计量大学 人工心脏无线供能系统的多物理场分析方法
CN110882056A (zh) * 2019-12-11 2020-03-17 南京亿高微波系统工程有限公司 Ct下精准肿瘤微波消融系统
CN110942832A (zh) * 2019-12-11 2020-03-31 南京亿高微波系统工程有限公司 一种微波消融系统
CN110942832B (zh) * 2019-12-11 2024-02-20 南京亿高医疗科技股份有限公司 一种微波消融系统
CN110882056B (zh) * 2019-12-11 2022-12-06 南京亿高微波系统工程有限公司 Ct下精准肿瘤微波消融系统
CN111027225A (zh) * 2019-12-25 2020-04-17 南京亿高微波系统工程有限公司 一种基于数据拟合的消融仿真方法
CN111027225B (zh) * 2019-12-25 2023-05-02 南京亿高医疗科技股份有限公司 一种基于数据拟合的消融仿真方法
CN111329454A (zh) * 2020-04-04 2020-06-26 吉林大学 基于人体三维温度场无损重构技术低温手术辅助诊疗系统
CN111528800A (zh) * 2020-04-30 2020-08-14 深圳开立生物医疗科技股份有限公司 一种肿瘤消融疗效预测方法、装置、设备及计算机介质
CN113768484B (zh) * 2020-06-10 2022-06-24 上海美杰医疗科技有限公司 生物组织个性化血流灌注率评估方法及系统
CN113768484A (zh) * 2020-06-10 2021-12-10 上海美杰医疗科技有限公司 生物组织个性化血流灌注率评估方法及系统
WO2021258639A1 (zh) * 2020-06-23 2021-12-30 南京诺源医疗器械有限公司 一种基于仿真温度场的微波消融术中疗效评估方法
CN111685870A (zh) * 2020-06-23 2020-09-22 南京诺源医疗器械有限公司 一种用于肝脏的微波消融碳化调控方法
CN111653363A (zh) * 2020-06-23 2020-09-11 南京诺源医疗器械有限公司 一种基于仿真温度场的微波消融术中疗效评估方法
CN112562859A (zh) * 2020-10-20 2021-03-26 应葵 肿瘤热消融手术智能仿真模型训练系统、训练方法
CN113842210A (zh) * 2021-08-02 2021-12-28 应葵 椎骨肿瘤微波消融手术仿真方法及装置
CN116258024B (zh) * 2023-05-16 2023-08-29 赛福凯尔(绍兴)医疗科技有限公司 目标区域的电场可视化方法、装置和计算机设备
CN116258024A (zh) * 2023-05-16 2023-06-13 赛福凯尔(绍兴)医疗科技有限公司 目标区域的电场可视化方法、装置和计算机设备
CN116705329A (zh) * 2023-06-08 2023-09-05 上海睿速创生医疗科技有限公司 一种基于COMSOL Muliphysics的消融针消融肝脏组织模拟方法

Similar Documents

Publication Publication Date Title
CN110263489A (zh) 基于dicom数据的肝肿瘤微波消融三维温度场仿真方法
Ekstrand et al. Influence of electrical and thermal properties on RF ablation of breast cancer: is the tumour preferentially heated?
Kok et al. Current state of the art of regional hyperthermia treatment planning: a review
CN103800075A (zh) 用于肝脏肿瘤消融术的患者特定建模的系统和方法
Ge et al. A multi-slot coaxial microwave antenna for liver tumor ablation
CN113693708B (zh) 射频消融设备的射频输出调整方法、装置及计算机存储介质
Singh et al. Coupled thermo-electro-mechanical models for thermal ablation of biological tissues and heat relaxation time effects
Gas Study on interstitial microwave hyperthermia with multi-slot coaxial antenna
CN112998849A (zh) 一种基于多针组合的微波适形消融方法
Wang et al. Numerical evaluation of ablation zone under different tip temperatures during radiofrequency ablation
Lodi et al. Towards the robust and effective design of hyperthermic devices: Improvement of a patch antenna for the case study of abdominal rhabdomyosarcoma with 3D perfusion
Baskaran et al. Design and experimental verification of 434 MHz phased array applicator for hyperthermia treatment of locally advanced breast cancer
Radosevic et al. Short pulsed microwave ablation: computer modeling and ex vivo experiments
Gao et al. Temperature simulation of microwave ablation based on improved specific absorption rate method compared to phantom measurements
Kok et al. Body conformal antennas for superficial hyperthermia: The impact of bending contact flexible microstrip applicators on their electromagnetic behavior
Wu et al. A simulation-experiment method to characterize the heat transfer in ex-vivo porcine hepatic tissue with a realistic microwave ablation system
Gas The s 11-parameter analysis of multi-slot coaxial antenna with periodic slots
Aghayan et al. An inverse problem of temperature optimization in hyperthermia by controlling the overall heat transfer coefficient
Sharma et al. A novel multiphysics optimization-driven methodology for the design of microwave ablation antennas
Surita et al. Finite element analysis for optimizing antenna for microwave coagulation therapy
Jessie et al. High‐Frequency Excitation and Surface Temperature Analysis of Breast Tissue for Detection of Anomaly
Fiser et al. Waveguide applicator system for head and neck hyperthermia treatment
Yamazaki et al. Development of a temperature distribution simulator for lung RFA based on air dependence of thermal and electrical properties
Puentes et al. Dual mode microwave tool for dielectric analysis and thermal ablation treatment of organic tissue
Beck et al. Experimental and computational evaluation of capacitive hyperthermia

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190920