《137Cs能谱的蒙特卡罗模拟.docx》由会员分享，可在线阅读，更多相关《137Cs能谱的蒙特卡罗模拟.docx（8页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、137Cs能谱的蒙特卡罗模拟 李梦楠 物理科学与技术学院 2019级核工程与核技术专业摘要：利用MCNP5程序模拟了多种情况下137Cs的能谱，并利用所学知识分析能谱特点及成因。关键词：137Cs MCNP5 射线与物质相互作用 1. 实验测定的137Cs能谱1.1 137Cs能谱及特点137Cs是实验室常用的放射源，半衰期约30.17年，衰变产物为137mBa，衰变产生0.662Mev的射线。图1: 137Cs的衰变纲图137Cs的能谱（采用NaI探测器）如图2所示，由于137Cs衰变产生的射线能量低于电子对效应发生的能量下限，因此射线与物质的相互作用方式主要有光电效应和康普顿散射。这些相互
2、作用反映在能谱上为0.662Mev的全能峰和能量从00.481Mev的康普顿坪。考虑到实际探测过程中，射线与闪烁体周围的物质包括放射源衬底材料及屏蔽材料等发生康普顿效应时，反散射光子返回闪烁体，通过光电效应被记录下来，同时根据入射光子和反散射光子能量的关系，我们发现当入射光子能量变化很大时，反散射光子能量变化不大，总是在200Kev左右，因此能谱上大约0.2Mev的位置会出现反散射峰。除此之外，由于衰变产物137Ba的K层特征X射线的贡献，在能谱的低能端会出现32Kev的X射线峰。图2：137Cs的能谱图2. 137Cs的蒙特卡罗模拟2.1蒙特卡罗模拟的思想方法中子和光子在物质中输运的宏观表现
3、是大量粒子与原子核微观作用的平均结果，蒙特卡罗方法通过逐一模拟和记录单个粒子的历程来求解输运问题。要得到比较合理的平均结果需要跟踪大量的粒子，至于单个粒子在其生命中的某一阶段如何度过，可以在已知统计分布规律的前提下通过抽取随机数来决定。2.2MCNP5的输入文件2.2.1 探测器部分本程序采用3*3Inch的NaI探头，外层采用0.7cm的铝膜包裹，后端为光电倍增管部分，利用塑料近似光电倍增管的结构材料，最外层为3cm的铅板屏蔽。探测器的几何结构如图3所示：图3：NaI探测器的几何结构2.2.2 放射源137Cs部分放射源采用点源，能量为0.661661Mev2.2.3 其他 采用光子-电子联
4、合运输方式； 能谱采用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铝吸收片的能谱而言，全能峰和康普顿坪非常清晰，低能端与实验测量结果略有出入，可能是由于低能端的X射线峰能量较低,只有32Kev
5、，与程序中NaI的单一粒子能量截断相近，在模型较为简单，模拟次数较少的情况下，难以得到理想结果。c. 就1cm铝吸收片的能谱而言，反散射峰比较平坦，峰位并不突出，但可以观察到反散射峰的中心位置仍在200Kev左右，与实验上的结果相符。峰位不突出可能与能量展宽系数设置有关。d. 5cm铝吸收片的能谱与1cm，2.5cm时的能谱相比，康普顿坪较为平坦，峰康比明显减小，可能是由于吸收片厚度增加同时立体角变小而导致计数率较少的缘故，在后续模拟中我们将继续讨论这一问题。2.3.2能谱与放射源位置的关系在上述探测器条件下，采用1cm厚的铝吸收片，改变放射源与吸收片之间的距离分别为0.5cm，2cm，8cm
6、，观察能谱的变化，模拟结果如下：图5：137Cs能谱与点源位置的关系a. 随着点源与吸收片距离的增加，能谱计数率下降，峰位保持不变，与实验上的结果一致，原因是当距离增大时，点源对探测器所张的立体角变小，探测器能够接收的粒子数目减少，因此计数率下降；b. 当点源距离为2cm和8cm时，基本看不到突出的反散射峰。原因是点源距离越远，对于探测器所张的立体角越小，与探测器周围材料发生康普顿散射的机会也较小，康普顿散射产生的散射光子通过光电效应再次被NaI吸收的概率也相应的减小了，因此反散射峰不明显。2.3.3能谱与吸收片种类的关系在上述探测器条件下，保持点源位置不变，吸收片厚度不变，改变吸收片的种类，
7、观察能谱的变化，模拟结果如下：图6：137Cs能谱与吸收片种类的关系a 吸收片Al的全能峰远远高于Zn和Fe，原因是射线的光电吸收截面与Z5成正比，30Zn和26Fe都远大于13Al，因此对射线的吸收Zn和Fe也远大于Al，也就造成被探测器记录到的粒子数Al远大于Fe和Zn。b Zn和Fe的康普顿坪高于Al，原因是在能量较低的情况下，射线的康普顿散射截面与Z成正比，Zn和Fe的原子序数较大，因此在探测器中与NaI作用概率较大，计数较多。c Zn和Fe的反散射峰高于Al,原因是反散射峰叠加在康普顿坪上，因此相应的Zn和Fe的反散射峰要高于Al，从截面的角度分析，Zn和Fe的原子序数较大，不仅与N
8、aI康普顿散射作用截面大，与周围材料的康普顿散射截面也大，产生的散射光子相应多，于是重新进入探测器被记录到的概率也相对大，因此计数率较多。3. 存在的问题与分析a. 采用点源近似放射源与实验所使用的放射源存在差别，如果利用环形放射源，与实际情况更相符。b. 在时间允许的情况下，增加模拟运算的次数可以使得能谱更接近实际情况。c. 能谱低能端的模拟不是很理想，考虑增加运算次数和减小粒子的截断能量来改善能谱形状。d. 在分析过程中采用的数据较少，每一部分只有三组数据，进一步考虑增加模拟的组数。4. 附录：MCNP5输入文件3*3 inch NAI(Tl)-point source1 1 -13.6
9、(1 -2 -19):(2 -9 18 -19):(9 -13 17 -20):(13 -14 -20) $PB2 2 -3.67 5 -6 -15 $NaI3 3 -2.7 (3 -4 -17):(4 -11 16 -17):(11 -12 -17) $Al4 0 (2 -21 -18):(21 -9 17 -18):(4 -5 -16):(5 -11 15 -16): (12 -13 -17):(1 -9 19 -20)$air5 4 -1.18 6 -11 -15 $plastic6 0 -1:14:207 3 -2.7 (-3 21 -17) $Al1 PY 02 PY 3 $ 3cm
10、 Pb3 PY 19.5 $ 15.5cm air 4 PY 20.2 $ 0.7cm Al 5 PY 20.9 $ 0.7 cm Air 6 PY 28.52 $ 7.62cm NaI7 PY 29.52 $ 1cm light trans8 PY 39.52 $ 10cm 9 PY 40.2 $10 PY 40.411 PY 52.712 PY 52.913 PY 59.914 PY 64.915 CY 3.8116 CY 4.7 17 CY 5.018 CY 9.019 CY 14.020 CY 15.021 PY 14.5 $ 5cm AlMODE P EIMP:P 1 1 1 1 1 0 1M1 82000 1 $PB -13.6 M2 11000 1 53000 1 $NAI -3.67M3 13027 1 $Al -2.7M4 1001 -0.0806 6012 -0.5998 8016 -0.3196 $PLASTIC -1.18SDEF ERG=0.661661 POS=0 14 0F8:P 2E8 0 500i 1.2FT8 GEB 0.01688 0.06165 0.34599ELPT:E 1 0.01 1 1 1 1NPS 10000000PRINT第 8 页