137csγ能谱的蒙特卡罗模拟(20201227170133).pdf

137Csγ能谱的蒙特卡罗模拟 李梦楠 物理科学与技术学院 2010 级核工程与核技术专业 摘要： 利用 MCNP5 程序模拟了多种情况下 137 Cs 的γ能谱，并利用所学知 识分析能谱特点及成因。 关键词： 137 Cs MCNP5 γ射线与物质相互作用 1. 实验测定的 137 Csγ能谱 1.1 137Cs 能谱及特点 137Cs 是实验室常用的γ放射源，半衰期约 30.17 年，衰变产物为 137m Ba，衰变产生 0.662Mev 的γ射线。 图 1: 137Cs 的衰变纲图 137Cs 的能谱（采用 NaI 探测器）如图 2 所示，由于 137Cs衰变产生的γ射线能量低于电子对 效应发生的能量下限， 因此γ射线与物质的相互作用方式主要有光电效应和康普顿散射。 这 些相互作用反映在能谱上为 0.662Mev 的全能峰和能量从 0~0.481Mev 的康普顿坪。考虑到 实际探测过程中， γ射线与闪烁体周围的物质包括放射源衬底材料及屏蔽材料等发生康普顿 效应时， 反散射光子返回闪烁体， 通过光电效应被记录下来， 同时根据入射光子和反散射光 子能量的关系，我们发现当入射光子能量变化很大时，反散射光子能量变化不大，总是在 200Kev 左右，因此能谱上大约 0.2Mev 的位置会出现反散射峰。 除此之外， 由于衰变产物 137Ba 的 K 层特征 X 射线的贡献，在能谱的低能端会出现 32Kev 的 X 射线峰。 图 2 ：137Cs 的γ能谱图 2. 137 Cs 的蒙特卡罗模拟 2.1 蒙特卡罗模拟的思想方法 中子和光子在物质中输运的宏观表现是大量粒子与原子核微观作用的平均结果， 蒙特卡 罗方法通过逐一模拟和记录单个粒子的历程来求解输运问题。 要得到比较合理的平均结果需 要跟踪大量的粒子， 至于单个粒子在其生命中的某一阶段如何度过， 可以在已知统计分布规 律的前提下通过抽取随机数来决定。 2.2MCNP5 的输入文件 2.2.1 探测器部分 本程序采用 3*3Inch 的 NaI 探头，外层采用 0.7cm 的铝膜包裹，后端为光电倍增管部分，利 用塑料近似光电倍增管的结构材料，最外层为 3cm 的铅板屏蔽。探测器的几何结构如图 3 所示： 图 3 ：NaI 探测器的几何结构 2.2.2 放射源 137 Cs 部分 放射源采用点源，能量为 0.661661Mev 2.2.3 其他 采用光子 - 电子联合运输方式； 能谱采用 500 道，最高能量 1.2Mev ； 模拟运算 107 次； 能谱展宽系数为 a=0.01688, b=0.06165, c=0.34599 。 2.3 蒙特卡罗模拟的能谱 2.3.1 γ能谱与吸收片厚度的关系 在上述探测器条件下， 设置点源距离探测器 5.5cm ，改变吸收铝片的厚度分别为 1cm，2.5cm ， 5cm ，观察能谱的变化。模拟结果如下图： 图 4 ：137Csγ能谱与吸收片厚度的关系 a. 随着吸收片厚度的增加，能谱计数率下降，峰位保持不变，与实验上的结果一致； b. 就 1cm 铝吸收片的能谱而言， 全能峰和康普顿坪