有关中子多球谱仪的计算问题.pdfVIP

有关中子多球谱仪的计算问题 陈军，杨小芸，容超凡 (中国原子能科学研究院，北京，102413) 摘要：为解决辐射场所中子能谱和剂量当量率的现场测量，中国原子能科学研究院放射性计 量测试部建立了两种类型的多球谱仪：分别基于主动式和被动式探测器。主动式谱仪的中心 探测器为球形3He正比计数器，被动式谱仪的中心探测器为球形金箔。利用MCNP程序计算 M 系统的响应函数，解谱程序有HEPRO程序包和UG程序包。 关键词： 中子能谱，中子剂量当量率，多球谱仪，球形3He正比计数管，球形金箔 1概述 八十年代初以来，美国、法国、德国、日本、英国、加拿大等许多国家先后开展了对反 应堆安全壳内外、轨道空间站内、航天航空飞行器内、高能加速器周围等辐射环境的中子能 谱和剂量当量率的测量工作，采用的主要方法包括：多球谱仪、3He谱仪、组织等效正比计 数器、雷姆仪、CR-39固体径迹探测器、热释光探测器、半导体探测器、气泡探测器等。国 内早期主要是采用雷姆仪对现场中子剂量当量率进行监测，随着我国核电、航空航天等事业 的迅速发展，对辐射场所的中子辐射剂量防护水平的评估提出了更高的要求，各种方法的研 究逐渐开展。中国原子能科学研究院放射性计量测试部建立了一套中子多球谱仪，并用其开 展辐射场所中子能谱和剂量当量率现场测量方法的研究。本文针对多球谱仪的有关计算问题 进行一些说明和讨论。 2多球谱仪 同直径的聚乙烯球、中心置有热中子探测器(如：3He正比计数器、BF3正比计数器、6LiI(Eu) 晶体探测器或活化箔等)组成。尽管其能量分辨率低，但由于可测中子能量范围宽(热能～ GeV)、适于现场应用、响应几乎各向同性及易于操作等优点，仍广泛地被各实验室采用。 2．1基本原理 从物理的角度，不同能量的单色快中子在不同大小的聚乙烯球内被慢化的程度不同，因 而响应也不同，这可以反映出中子能量的变化情况。对于小球，热中子大部分被慢化体(聚 乙烯)俘获吸收(H(n，丫)反应)，低能中子慢化成热中子后(与H和C发生弹性散射)有较大 的几率到达中心的热中子探测器并被探测到，而高能中子(与H和C相互作用后仍具有较高 能量)趋于逃逸，因而其响应函数的峰值位于低能区；对于大球，热中子基本被俘获吸收， 低能中子慢化成热中子后也被大量俘获吸收，而高能中子慢化成熟中子后有更大的几率到达 中心的热中子探测器并被探测到，因而其响应函数峰值位于高能区，如图l所示。从数学的 角度，假设有m个不同直径的球，第i个球的响应函数为Ri(Z)，被测中子能谱在不同能量的 注量为烈囝，则第i个球的探测器计数M(或饱和活度Ai)应为： (1) M(4)=卜(E)R，(E)dE，i=I，2…，m 考虑离散能量情况，则式(1)可以转化为： (2) M(4)=∑①(q)尺，(g)-i=l…2．，m；j=l…2．．，n J 370 其中，妒(岛)为能量区间蝇局+l】的入射中子注量，蜀喝)为能量区间蝎，局+l】的平均响应函数。 式(2)是式(1)的近似，近似程度随能量区间划分个数n的增加而增高。一般来讲，mn， 无法用解析的方法求得真解，但可以通过一些数学方法(如：迭代法、M．C方法)求得近似 解。 图l 多球谱仪响应函数分布图 of Distributionsofthe functions Fig．1 response multi-spherespectrometer 2．2谱仪的结构 本工作中，主动式多球谱仪包括直径分别为2．5、3、3．5、4、5、6、8、10和12英寸九 和120kPa的氪气。被动式多球