计算X射线通量的几种方法的比较研究.pdfVIP

236第三届全国信息与电子工程学术会议、四川省电子学会曙光分会第十四届学术年会暨院青年科协第八届学术年会论文集 计算X射线通量的几种方法的比较 杜川华‘许献国赵洪超 (中国工程物理研究院电子工程研究所 四川 绵阳 621900) 摘要综述了目前计算X射线在材料中通置的几种方法：利用质量衰减系数直接计算的理论 方程和描述粒子运动统计规律的蒙特卡罗方法．文中分别采用理论方程、蒙特卡罗软件MCNP4C和 EGSnrc三种途径计算了黑体谱X射线穿过几种常用复合材料的透射率和透射谱，根据计算过程和计 算结果。结合相关参考文献，分析比较了理论方程与蒙特卡罗软件之间、两种蒙特卡罗软件MCNP4C 和EGSnrc之间的差异． EGSnrc透射率透射谱 关键词X射线蒙特卡罗Mo岬4C 1 引言 计算X射线在材料中输运的通量问题，可以采用理论方程直接计算，或者采用蒙特卡罗(Monte Carlo)软件进行计算，目前国内外已经有大量相关的文献。本文简要介绍了这两种方法，并以计算 黑体谱x射线穿过几种常用复合材料的透射率、透射谱为例，根据计算过程和计算结果，结合相关 参考文献，分析比较了理论方程与蒙特卡罗软件之间、两种蒙特卡罗软件MCNP4C和EGSnrc之间的 差异。 2理论方程 数衰减【1】： ，=Ioexp[一(／a／p)m] (1) 式中，厶为入射X射线光子强度，Pg／cm3)为吸收材料的密度，[z／p(cm2／g)为质量衰减系数： P／P2tYmt／uA=(仃肛+叮∞曲+仃删+盯p咖+o帅)／uA (2) 式中， 口耐为光子与物质相互作用的总截面，通常情况下由以下5项组成：光电效应作用截面盯Ⅳ， 三产生截面叮帅。根据amm和Hubbcll等人的研究，在5MeV一--．40MeV的能量区间，还有一种光核截 面对总截面有大约5％～10％的贡献，但是因为作为稳定元素的各种同位素的质量A和原子序数Z相应 的核结构存在不规则变化的复杂关系，因此目前的x射线质量衰减系数数据库中都还没有考虑该项 因素。”为原子质量单位，A为吸收材料的相对原子质量。复合材料的质量衰减系数，按以下公式计 算： (∥／p)c=∑i(嘭·(∥／圳 (3) ’作者简介：杜川华， (1977．12)，女，助理研究员，硕士，毕业于电子科技大学，主研方向：辐射效应仿真计算及 试验研究．邮箱：绵15日919信箱522分箱：邮编：621900；电话：2484513E-mail：chhdu@sian．com． 第三届全国信息与电子工程学术会议、四川省电子学会曙光分会第十四届学术年会暨院青年科协第八届学术年会论文集 237 式中，(∥／p)，为复合材料中各组成元素的质量衰减系数，Wi为各组成元素的质量百分比。 当X射线具有一定能谱分布时，一般将能谱区间分为N个能量子区间，把每个子区间的X射线看 成具有中心能量的单能射线，并以每个子区间中心能量对应的质量衰减系数代替该子区间的质量衰 减系数。假定总入射光子强度为1．则透射率可写为： ，：．芒n1p一(声，户)“p．J ，=∑qP一‘p’％p。 Ill (4) 式中，q为每个子区间的光强占总辐射光强的比值。 3蒙特卡罗方法 蒙特卡罗方法也称为随机抽样技巧，是目前唯一能够对辐照目标内部三维空间核粒子碰撞轨迹 和能量分布进行精确仿真计算的工具，在无法用电子设备测量的材料界面或者区域，该方法也是一 个预测辐射剂量、注量以及能谱等物理参数的重要工具。用蒙特卡罗方法计算光子响应的问题，是 光子和电子的偶合输运问题【3】：光子由源发出后，与原子核发生作用时，如果发生康谱顿散射，把 光子存储起来，去模拟电