HIFU经颅聚焦脑肿瘤治疗温度场的数值仿真研究
本文关键词: HIFU 经颅 脑肿瘤 温度场 CT图像 建模 焦域 出处:《天津医科大学》2014年硕士论文 论文类型:学位论文
【摘要】:高强度聚焦超声(High intensity focused ultrasound,HIFU)是一种无创、无毒、无副作用的肿瘤治疗手术,通过超声波的机械效应、热效应、空化效应等作用机制使靶区的温度在数秒内上升到60℃以上,导致肿瘤细胞发生凝固性坏死,目前该技术已应用于实体肿瘤的临床治疗。而对于有颅骨遮挡的脑肿瘤HIFU经颅治疗而言,由于超声经颅传播时产生相位偏移和幅值衰减,在颅骨内不能形成焦域或形成的焦域偏离靶区;同时又由于颅骨与周边软组织的声阻抗相差较大,在HIFU治疗颅内浅表组织肿瘤过程中可能在颅骨处产生热量沉积,引起颅骨以及周围软组织的损伤等问题,限制了 HIFU治疗脑肿瘤在临床中的应用。研究目的参照人体CT头颅结构,建立三维HIFU经颅骨治疗模型,采用声波传播方程式和生物热传导方程数值仿真HIFU经颅骨治疗的温度场分布,利用时间反转法对焦域位置进行调控和颅骨温度较高处的热点消除,进而分析讨论输入声强、辐照时间、换能器与颅骨相对距离对颅内形成焦域的影响。其研究结果可为临床HIFU脑肿瘤治疗剂量及其治疗计划的制定提供理论参考数据和理论依据。研究方法本课题主要对HIFU经颅治疗进行仿真建模并数值仿真其温度场分布。.在圆柱坐标系下,参照人体头颅CT图像建立均质化头颅参数的结构与八元环换能器构成该模型,利用高强度超声传播的非线性理论和时域有限差分(Finite Difference Time Domain,FDTD)的数值仿真法,数值仿真计算不同时间和空间超声波非线性传播的声压,并结合Pennes生物热传导方程,在不同治疗条件下数值仿真HIFU治疗脑肿瘤,并对焦域位置进行调控和降低颅骨处温升以优化HIFU经颅治疗,仿真研究辐照因素对颅内形成焦域及其温度场分布的影响。而对于基于CT图像建立的人体头颅模型,在直角坐标系下采用64阵元换能器,在不同治疗条件下数值仿真HIFU治疗脑肿瘤过程中形成的温度场分布,将均质化头颅模型中分析讨论的相控方法应用在该模型中,进而分析辐照因素对颅内形成焦域及其温度场分布的影响。研究结果一、均质化头颅参数模型的数值仿真:1、采用时间反转法后形成的焦点位置与设定80mm的焦距仅相差0.4mm;2、采用热点消除法后颅骨处最高温升仅为2.4℃,达到了降低颅骨温升的目的;3、在辐照时间相同的条件下,当形成一定焦域体积时进行热点消除后的经颅HIFU治疗所需要的输入能量要大于无颅骨治疗时的输入能量;4、经过热点消除后,不同输入声强开始形成焦域所需要的时间不同,输入声强越大所需时间越短,且焦域体积随辐照时间的延长而非线性增加;输入声强越大,焦域体积随时间的增长率越大;5、声强越大,达到相同大小的焦域体积所需要的辐照时间越短;声强越小,达到相同大小的焦域体积所需要的辐照时间越长;6、在治疗深度为20mm的条件下,当换能器与颅骨之间距离不同时,形成的焦域体积随辐照时间变化率不同。二、人体头颅模型的数值仿真:1、采用时间反转法后形成的焦点位置与设定80mm的焦距仅相差0.75mm;2、采用热点消除法后颅骨处最高温升仅为9.7℃,达到了降低颅骨温升的目的;3、当输入声强一定时,颅骨温升以及焦点温升均随辐照时间延长呈非线性增加,并且经过热点消除后的焦点温升略有升高,在进行热点消除时对实际焦点处的温升随时间变化趋势影响不大;4、随着治疗深度由深到浅,采用时间反转相控聚焦法形成的焦域逐渐减小,而当热点消除后形成的焦域也随治疗深度的变浅而体积逐渐减小,当治疗深度为20mm时未形成焦域,因此根据肿瘤靶区位置不同,需要消除热点的案例也不同,形成的焦域体积随治疗深度也不同;5、当设定治疗深度一定为25mm时,换能器与颅骨之间距离不同时,经过热点消除后,能够形成焦域体积最大时的颅骨位置发生了改变,因此在制定治疗临床计划时,需要进行数值仿真计算来预测能够形成最大的焦域体积以选择最佳的治疗位置。综合上述仿真结果,得到如下结论:1、时间反转法可以调控经颅HIFU的焦域位置,使其位置能够到达目标靶区;2、在治疗肿瘤位于浅表组织位置时,利用热点消除法可降低颅骨的温升,同时焦点处的温升略有降低,且对焦点影响不大;3、HIFU经颅治疗所需的能量大于HIFU不经颅骨的治疗;4、在相同焦域体积条件下,辐照时间与辐照声强呈负相关;5、HIFU经颅治疗形成的焦域体积随输入声强以及辐照时间的增加而非线性增加;6、换能器与颅骨间的距离对焦域的形状、大小均有影响。
[Abstract]:High intensity focused ultrasound ( HIFU ) is a non - invasive , non - toxic and side - effect tumor treatment operation . In the cylindrical coordinate system , the structure of homogenized skull parameters and eight - membered ring transducer are established with reference to the human skull CT image . The influence of radiation factors on intracranial formation of focal area and temperature field distribution is studied by using the numerical simulation method of high intensity ultrasonic propagation . 3 . When the treatment depth is set to 25mm , the focal area of the skull can be decreased with the depth of treatment .
【学位授予单位】:天津医科大学
【学位级别】:硕士
【学位授予年份】:2014
【分类号】:TB559;R739.41
【相似文献】
相关期刊论文 前10条
1 常满,郭秀梅;海面电磁波散射的数值仿真[J];地球物理学报;1986年02期
2 王新慧;李旭东;;气液两相圆锥绕流的数值仿真[J];计算机仿真;2013年08期
3 ;院级重点学科 四川省高等学校数值仿真重点实验室简介[J];内江师范学院学报;2007年06期
4 杨冬梅;王晓鸣;;混凝土中爆炸数值仿真算法研究[J];爆炸与冲击;2005年06期
5 童启秀;王胜兵;;一阶常微分方程组加速算法及数值仿真[J];计算机与数字工程;2011年07期
6 刘文锋;半主动半被动结构控制(英文)[J];青岛建筑工程学院学报;1997年04期
7 孙学金;孙海洋;江志东;;不同大气条件雨滴下落速度的数值仿真[J];计算机仿真;2011年12期
8 杨纪华;;基于Matlab混沌系统的数值仿真[J];绵阳师范学院学报;2013年02期
9 徐晨,林祥都;FRM的非线性动力学模型及数值仿真(英文)[J];现代电子技术;1998年04期
10 王惠;张会生;许庆新;;自由液面下增阻器水动力特性的数值仿真[J];计算机仿真;2009年11期
相关会议论文 前10条
1 陈秉智;聂春戈;兆文忠;;列车大变形碰撞的数值仿真[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
2 赵振洋;杜家政;;公交车踏步结构数值仿真及设计改进[A];北京力学会第20届学术年会论文集[C];2014年
3 王欢;臧润清;;家用自然恒温保鲜箱的模拟研究[A];第五届全国冷冻冷藏产业(科技、管理、营销)创新发展年会论文集[C];2013年
4 黄英兰;徐元铭;;复合材料格栅壁板屈曲的数值仿真[A];北京力学会第14届学术年会论文集[C];2008年
5 陈虹;夏芳;沈卫超;宋磊;;高功率微波数值仿真通用数据格式及其I/O接口[A];2010年第16届全国信息存储技术大会(IST2010)论文集[C];2010年
6 陈虹;夏芳;宋磊;沈卫超;;高功率微波数值仿真通用数据格式及其I/O接口实现[A];中国工程物理研究院科技年报(2008年版)[C];2009年
7 姜晓婷;刘东顺;吴立顺;李文龙;菅喜岐;;高强度聚焦超声热疗焦区温度分布的时域有限差分法数值仿真[A];中国生物医学工程进展——2007中国生物医学工程联合学术年会论文集(上册)[C];2007年
8 陈睿;纪冬梅;沈忱;贡顶;;单粒子效应TCAD数值仿真的三维建模方法[A];第十六届全国核电子学与核探测技术学术年会论文集(上册)[C];2012年
9 李根国;;高性能计算在制造业和工程领域的若干应用[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
10 伍志辉;田怀璋;高原;陈林辉;梁俊杰;石毅登;;制冷剂在管内蒸发换热过程的数值仿真[A];第六届全国低温与制冷工程大会会议论文集[C];2003年
相关重要报纸文章 前1条
1 北京信息与控制研究所 王正中;仿真学科的研究与思考[N];大众科技报;2007年
相关博士学位论文 前6条
1 马艳琴;增程外弹道理论分析及数值仿真[D];南京理工大学;2003年
2 焦培刚;基于SPH的流体仿真数值算法及工程应用研究[D];山东大学;2010年
3 白兰;基于MEMS的无阀泵数值仿真与实验研究[D];中国科学院研究生院(长春光学精密机械与物理研究所);2005年
4 胡岳;基于数值仿真与PIV的楔形流量传感器流场研究[D];天津大学;2014年
5 李宝良;线接触磨损数值仿真及应用研究[D];大连交通大学;2010年
6 袁玲;不锈钢冶炼AOD炉烘烤新工艺过程数值仿真与优化研究[D];东北大学;2009年
相关硕士学位论文 前10条
1 王永飞;基于CFD的电子水泵数值仿真与流场分析[D];长安大学;2015年
2 顾闻阳;PIN二极管电特性的谱元法快速分析[D];南京理工大学;2015年
3 朱正象;履带底盘件应力数值仿真及实验研究[D];山东大学;2015年
4 刘俊雄;基于云计算的登峰变电所开放式数值仿真系统[D];东北石油大学;2015年
5 袁倩;单兵筒式武器液体平衡发射过程数值仿真及试验研究[D];南京理工大学;2016年
6 胡名雄;自旋型智能雷发射过程数值仿真及试验研究[D];南京理工大学;2016年
7 陆文山;航天器自旋人工重力场中航天员行走动力学与控制研究[D];南京航空航天大学;2016年
8 罗美千子;Trijet组合循环发动机排气系统数值仿真与分析[D];南京航空航天大学;2016年
9 周鑫;基于微量润滑磨削的喷嘴雾化仿真与实验研究[D];长沙理工大学;2015年
10 张千;HIFU经颅聚焦脑肿瘤治疗温度场的数值仿真研究[D];天津医科大学;2014年
,本文编号:1456070
本文链接:https://www.wllwen.com/guanlilunwen/gongchengguanli/1456070.html