BNCT治疗规划系统MCDB算法及测试.pdf

第29卷第5期 计 算 物 理 V01．29．No．5 CHINESE OFCOMPUTATIONALPHYSICS 2012年9月 JOURNAL Sep．，2012 文章编号：1001—246X(2012)05-0721-06 BNCT治疗规划系统MCDB算法及测试 李 刚1’2， 邓 力1， 陈朝斌1， 叶 涛1， 莫则尧1 (1．北京应用物理与计算数学研究所，北京100094；2．中国T-程物理研究院北京研究生部，北京100088) 擒 要：给出BNCT治疗计赳系统软件MCDB的算法及测试结果，通过使用材料矩阵描述BNCT的网格模型，配合 快速径迹算法和计数矩阵，进行粒子输运模拟及计数，模拟结果显示，MCDB的计算精度与MCNP相当，计算速度 较MCNP提高3．1—3。4倍．MCDB模拟l000万粒子的时间不足2个CPU小时，可保证大部分网格齐j量误差在5％ 以内，基本上达到了BNCT临床要求，可用于30kW的医院中子照射器的临床治疗． 关键词：BNCT；MCDB；快速径迹算法；蒙特卡罗 中图分类号：0571．51 文献标识码：A 0 引言 恶性肿瘤是危害人类的一大顽症，我国癌症患者现患病例310万余人，每年新增愈百万，其中头颈部肿 Neutron 瘤约10万人、皮下恶性黑色素瘤超过1万人．国际研究表明硼中子俘获治疗(BoronCaptureTherapy， 简称BNCT)作为一种崭新的二元靶向放射治疗手段，目前的治疗效果已经优于常规治疗手段¨1． BNCT治疗计划系统可以给出病人患部治疗剂量，根据病人患部硼浓度和照射中子注量率确定照射剂 用户定义区域，该格式中包含了几何和物理特征．所有的体素具有相同的长、宽和高，若干体素组合起来可以 配合Mesh计数．这两个软件通过多年的临床检验，技术已经趋于成熟． MCDB程序利用三维材料矩阵描述脑部模型的网格结构，在网格上采用快速径迹算法进行输运，同时配 规模多网格模型的所有网格进行计数，并且实现了对组合网格的模拟。．模拟结果显示MCDB的模拟速度 是MCNP的3．1～3．4倍，剂量计算结果可靠． 1 MCDB程序简介 Cado in MCDB(Monte Dosimetry 其中前处理用于读取和显示CT／MRI医学影像数据，配合中子照射条件产生粒子输运模拟输入文件；后处理 模块对粒子输运的结果进行处理和作图展示，以利于用户对计算结果的分析和掌握；粒子输运模块就是采用 蒙特卡罗方法对中子源产生的粒子进行模拟，它是整个治疗系统的核心部分，基于MCNP程序研制形成． 1．1三维材料矩阵 BNCT的单一网格模型结构简单，由完全相同的立方体网格(或长方体网格)组成，排列有序，邻域关系 清楚．随着网格尺寸的变小，网格数多达几十万，而材料和密度种类相对很少，最多几百种．MCNP程序在模 拟BNCT模型时计算粒子射线与网格面的交点要花费大量时间，除此以外，材料的质量密度、原子密度以及 收稿日期：2011—12—01；修回日期：2012—04—06 (2011A0103006)和国家磁约束聚变能研究专项(2010GBlll001)资助项目 万方数据 722 计 算 物 理 第29卷 网格的6个面号等都是网格编号的函数，这些数组很庞大，搜索时间相当长，特别是寻找粒子进入的下一个 网格，要先找到粒子所到达网格面的编号，再通过它找到另一侧的网格号．初期的MCDB程序只是在MCNP 的基础上发展了快速径迹算法，几何结构描述与其相同，模拟时间较MCNP节省了30％～36％，没有很好的 发挥出快速径迹算法的优势，并且不适于全部网格计数，因此考虑对几何结构采用其它描述方式，以求进一 步缩短模拟时间． 材料矩阵中进行输运，配合快速径迹算法，实现了全部网格计数．通过简单的拉伸映射把模型中的网格都映 度和原子密度改成材料编号的函数，同