powerpointpresentation-indico.pptVIP

中子活化瞬发伽马分析的MCNP模拟研究 报告人：朱传新 提要 引言 计算模型 结果与讨论 结论 引言 中子活化瞬发伽马分析在煤质成分检测、水泥成分分析等物料分析中，是一种重要的分析手段。 Carlos等科学家采用MCNP程序对于水泥生料密度和水分进行PGNAA系统的伽马计数率修正以及进行系统的刻度工作。 Robin和Charles对于NaI探测器的非线性问题在PGNAA装置上的应用进行了研究。 Zhang Wenchao 和Robin P. Gardner 对于PGNAA实验数据开展了近似线性插值处理方法的蒙特卡罗模拟研究。 Ali Asghar Mowlavi等人对于PGNAA系统伽马谱的脉冲堆积效应进行了蒙特卡罗模拟研究。 在本工作中，采用MCNP5程序和ENDF/B6数据库进行了锎源中子活化瞬发伽马分析的模拟计算，给出了所配制样品中Fe、S元素的特征谱，并与实验进行了对比分析。 计算模型 实验装置由252Cf中子源、聚乙烯、样品、NaI探测器等构成。 计算程序采用MCNP5程序，截面数据使用ENDF/B6库数据。装置的尺寸数据取用了实验中所用数据。 计算模型如图所示，通过对于实验装置的适当简化得到。 实验装置的计算模型 计算中，采用中子光子耦合输运方式，探测器上光子的记录采用了F4、F8等方式。F4记录方式计算的是探测器栅元平均光子通量，这种记录方式的计算效率比较高，因为它是直接反映伽马射线进入NaI探测器的总量情况。 F8记录方式计算的是由于探测器辐射所引起的脉冲能量分布，这种记录方式的优点是考虑了探测器对于伽马射线的效率因素，这样在与实验数据进行比较时，就可直接进行比较，但是缺点是计算的效率比较低，需要较长的计算时间（如3d～5d），才能得到不确定度较小的结果。 实验样品总质量约30kg，主要成分列于表1。主要由C、H、O、Si等元素构成。在这些主要成分基础上，按照一定质量比为0～2%掺入S元素和Fe元素进行计算，从而分别对于产生的S元素特征峰和Fe元素特征峰进行分析。 结果与讨论 样品中的硫元素计算分析 利用F4方式计算得到了252Cf中子源经过装置后与样品作用产生的瞬发伽马射线在NaI探测器上的体元平均通量。样品中的S元素的特征峰5.4MeV通量随着S元素含量的变化示于下图。 由图可以看出，当样品中的S含量增加时，特征峰计数近似成线性增加，这表明通过S元素特征峰的定量分析，根据这种成分和特征峰计数之间的线性关系，能够对应给出S含量数据。 样品中含2%的S元素时的γ能谱 结果表明：对于2MeV以下区域，计算值与实验值的差异在10%左右，吻合较好。 对全谱总数比较来看，计算值与实验值的差异在6%左右。计算值偏高。 4.4MeV至5.8MeV总数（三个峰的总量），比值为：0.94，二者是比较接近的。 通过谱形比较表明：实验与计算的峰位基本吻合。 Fe元素的计算结果与分析 对于掺入质量百分比为1%的Fe元素时，采用F8方式进行了计算，结果与对应的实验数据进行了绝对比较。掺入质量百分比为1%的Fe元素时，样品的光子能谱的计算结果与实验结果示于图4。 通过分析表明：从6.3MeV～8MeV铁峰总计数：计算值与实验值之比为1.04，表明二者是在不确定度范围内吻合的。 在描述252Cf中子源分布时，通常有两种方法，一种是采用Maxwell分布，一种是采用Watt谱，这两种谱不同描述，究竟对于计算有多大影响呢？ 针对这个问题，我们分别采用两种分布的相应参数对于所建立的模型进行了计算分析。 可以看出，在8MeV以下，二者基本重合，在10MeV以上有些差异，这也与高能部分计算的统计涨落较大有关，并且使用Watt分布谱中子源时，所得伽马能谱结果在高能部分离散比较大，相比与Maxwell分布条件下，与实验差异更大一些，这表明Watt谱条件下的计算，在8MeV以上的高能伽马谱部分收敛效果比采用Maxwell分布中子源时要差。 两种不同的中子能谱描述方式，对于中子活化瞬发伽马分析的MCNP理论计算而言，其影响是可以忽略的，并且采用Maxwell分布描述252Cf中子源时，计算的收敛效果较好。因此，针对这两种分布开展计算时，建议选用Maxwell分布方式。 结论 采用MCNP5程序和ENDF/B6数据库进行了锎源中子活化瞬发伽马分析的模拟计算，给出了所配制样品中分别掺入不同比例Fe、S元素后的特征谱。研究表明： （1）对于2MeV以下区域，计算值与实验值的差异在10%左右，吻合较好。 （2）对全谱总数比较来看，计算值与实验值的差异在6%左右。 （3）与实验对比，计算能起到校验作用，并且表明实验是可信的。 （4）通过采用Maxwell分布谱和Watt分布谱，这两种