含氢介质内中子能谱测量

合集下载

中子探测的基本方法及13.5常用中子探测器

中子探测的基本方法及13.5常用中子探测器

13.6 中子注量率测量的主要指标
中子灵敏度
13.7 堆用探测器
13.6 中子注量率测量的主要指标 中子灵敏度
R 中子灵敏度定义:
0
反应的发生率 中子注量率
R N t ( E ) ( E )dE
Nt 为探测器灵敏体积 中辐射体的靶核数。
对能量低于30keV的中子: 30 keV v 0 0 R Nt ( E )dE 0 v 由 ( E ) n( E ) v n(E)为能量E处单位能量 间隔的中子密度。 v为中子速度。
反应截面与中子能量的关系:

100000
0 v0
v
1 1 v Tn
B-10 Li-6 He-3
capture cross section(barn)
10000
1000
100
10
1
0.1 1E-111E-10 1E-9 1E-8 1E-7 1E-6 1E-5 1E-4 1E-3 0.01
可选用微型裂变室,且电极涂235U+239Pu(可增殖, 总积分通量由1.7×1021提高到4.8×1021中子); 也可以用自给能探测器。
3) 功率量程:大于1010/cm2s;足够大,本底相 对较小;用电流型裂变室或硼电离室。
2. 堆芯探测器——堆芯内中子注量率的空 间分布。
要求体积小,寿命长; 典型工作条件:
~ 5 10 / cm s
8 2
8

本底 ~ 10 R / h
工作温度 ~ 300 C 2 ~ 2500 N / cm 工作压力
v0 为热中子的最可 其中 v 为中子的平均速度, 几速度。
对热中子,在T=20C时,v / v0

测井曲线的原理及应用 - 副本

测井曲线的原理及应用 - 副本
岩石电阻率与岩性,储油物性和含油性有着密切的关系。 电阻率测井的主要任务是根据测量的岩层电阻率,判断岩 性,划分油气水层,研究储集层的含油性、渗透性和孔隙性。
主要用途是: 1.确定地层的电阻率; 2.计算储层的含水饱合度; 3.判断油、气、水层。
电祖率测井主要作用
求解含油饱和度 例:Archie 公式(1942年)
(1)能准确地确定地层界面的深度,并能详细地划分薄地层。 (2)能判断地层的岩性和渗透性。 (3)能计算储集层的储集性和含油性参数。 (4)能划分和评价油层、气层和水层。
1、电阻率测井系列 提供地层真电阻率和侵入带电阻率以及泥浆侵入状况,确定 储层的含水饱和度。 2、岩性—孔隙度测井系列 用于识别岩性、计算地层孔隙度,判识油气、水层
地层因素 :
F

Ro Rw

a
m
电阻率增大系数: I Rt b
Ro
S
n w
含水饱和度:
Sw n
abRw
m Rt
泥浆侵入特征 在钻井过程中,井内泥浆柱的静压力通常大于地层压力,此压力差
使泥浆滤液进入渗透性地层,叫泥浆侵入。 泥浆中固体颗粒沉淀于井壁形成泥饼。泥饼的渗透性较差,因此形
岩石体积密度 ρb=(1-Vsh-)ρma+Vshρsh+ Sxoρf+ (1-Sxo)ρh
Ρh --岩石骨架、泥质、泥浆滤液和油气的密度; Vsh--泥质含量; --有效孔隙度; Sxo--冲洗带含水饱和度。
密度和岩性—密度测井的应用
1)确定岩性和孔隙度 (1)根据密度曲线和岩心分析资料回归“密度—孔隙度”经 验公式,或分析资料与密度、声波时差、中子等测井参数经多 元回归的经验公式,再计算新井的地层孔隙度。 (2)岩性单一时,也可以用以下公式计算孔隙度(φ)。

中子能谱测量中的解谱技术研究进展

中子能谱测量中的解谱技术研究进展

㊀第42卷㊀第4期2022年㊀7月㊀辐㊀射㊀防㊀护Radiation ProtectionVol.42㊀No.4㊀㊀July 2022㊃综㊀述㊃中子能谱测量中的解谱技术研究进展黄迁明,刘㊀斌,陆㊀婷,王㊀波,唐松乾,吕焕文,应栋川,翟梓安(中国核动力研究设计院核反应堆系统设计技术重点实验室,成都610213)㊀摘㊀要:中子能谱解谱技术为中子能谱测量系统必要的组成部分,近几十年来国内外开展了大量研究㊂本文首先介绍了中子能谱常规解谱流程,包括解谱模型㊁响应函数㊁解谱误差等内容;接着详细介绍了国内外中子能谱测量技术研究现状以及中子能谱解谱算法研究现状,包括比较成熟的最小二乘算法㊁最大熵算法等,也有新兴的神经网络算法㊁遗传算法等,总结了不同解谱算法的特点;接着介绍了根据不同解谱算法发展的解谱程序,对比了不同解谱算法及程序的优缺点,基于最小二乘算法开发的SAND 系列程序和基于最大熵算法开发的MAXED 程序是解谱功能强大㊁使用最广泛的程序;最后梳理了中子能谱解谱方法的发展脉络,总结了国内和国外研究的区别,未来开发包含多种解谱方法的综合性解谱程序具备较强的应用需求㊂关键词:中子能谱测量;中子能谱解谱方法;广义最小二乘算法;最大熵算法;蒙特卡罗算法;压缩感知理论中图分类号:TL8;O571.54文献标识码:A㊀㊀收稿日期:2021-07-21作者简介:黄迁明(1992 ),男,2015年毕业于四川大学核工程与核技术专业,2020年毕业于北京大学粒子物理与原子核物理专业,获博士学位,工程师㊂E -mail:qianming.huang@通讯作者:刘斌㊂E -mail:liubin871204@㊀㊀中子能谱测量是中子探测学的重要组成部分[1-4],主要用于获取准确的中子能谱信息,通过相应解谱算法的研究,可在一定程度上提高实验测量的精度,且解谱算法性能的提高有助于在有限测量条件下降低探测系统硬件的要求,在各类核装置的设计和运行㊁辐射防护计量学和屏蔽验证㊁核军备控制以及反恐等应用中发挥着关键作用㊂近年来,在核不扩散和核军备反恐安检的大应用背景下,基于中子能谱的核素识别技术对于核材料准确㊁高效的识别具有很强的优势[5-9],国内外在中子能谱的解谱方法上进行了广泛研究㊂中子能谱测量难度较大,主要面临中子能谱差异大㊁中子能量跨度区间大㊁伽马本底干扰㊁中子场的统计涨落和测量系统噪声等问题[10-12]㊂目前中子能谱测量的手段主要有:飞行时间方法[13]㊁有机闪烁体测量方法[14-19]㊁多球谱仪方法[20-23]以及多箔活化方法[24-26]等㊂除飞行时间法外,其余三种测量方法的中子能谱均不能直接获取,而是需要结合探测系统的响应函数从测量值中进行解谱获得,即通过探测器探测到一个输入能谱,根据探测器对不同中子的响应特性,反推出测量点的真实能谱㊂现阶段能谱测量的难点主要由探测器测量数据统计误差㊁测量系统的噪声干扰㊁解谱函数不具有唯一解以及先验谱的可靠性造成,通过优化探测系统和解谱算法能一定程度上提高中子能谱测量精度㊂1㊀国内外中子能谱解谱方法研究现状1.1㊀解谱模型㊀㊀在经典的解谱模型中,对于有机闪烁体㊁多球谱仪以及多箔活化测量系统,探测器测量结果均可用式(1)所示的第一类Fredholm 积分方程表征[27-28]:M i ʃεi =ʏE max E minR i (E )ϕ(E )d E (1)式中,M i 为探测器计数,对于有机闪烁体测量系统,其为第i 道闪烁光子脉冲高度;对于多球谱仪测量系统,其为第i 个慢化球对应的热中子探测器计数;对于多箔活化测量系统,其为第i 个活化箔辐照结束后的放射性活度;εi 为测量结果的不确定度,其通常与探测器计数的统计特性㊁中子场统计涨落以及测量系统电子学噪声等因素相关;R i㊃562㊃㊀辐射防护第42卷㊀第4期为第i 个探测器的中子能量响应函数,对于有机闪烁体测量系统,其为第i 道多道不同中子入射产生的闪烁光子脉冲高度分布;对于多球谱仪测量系统,其为第i 个慢化球对应的中子能量响应函数;而对于多箔活化测量系统,其为第i 个活化箔对应的活化反应的多群截面;ϕ(E )为所要求解的中子能谱㊂而求解式(1)面临诸多困难,包括:(1)方程组不具有唯一解;(2)误差引起解的剧烈震荡,从而造成求解的失败;(3)探测器计数的统计特性㊁测量系统的噪声及中子辐射场本底造成求解困难;(4)解的非负值约束㊁系统响应不确定性的传递等问题㊂因此,该类问题在测量领域被称为 逆问题 ,利用测量系统的响应函数对真实能谱进行倒推成为解谱的新思路㊂1.2㊀响应函数获取响应函数的获取是中子能谱解谱的前提,针对不同的中子能谱测量方式可以制定对应的响应函数测量方法,一般有实验测量方法和理论计算结合实验标定方法㊂1.2.1㊀实验测量㊀㊀在有机闪烁体测量实验中,通常的手段主要采用单色性较好的加速器中子源结合飞行时间法对中子能量进行甄别,比如Lawrence 等人[29]采用加速器中子源结合飞行时间窗对EJ -309和EJ -299-33有机闪烁体的响应函数进行测量得到图1的响应函数㊂而在多球谱实验中,主要的实现方法是采用加速器中子源结合慢化材料以实现不同能量中子的响应函数标定,同样采用飞行时间方法对中子能量进行甄别,图2为Pioch 等人[30]采用加速器中子源对多球谱仪系统进行能量标定的方案㊂图1㊀EJ -309和EJ -299-33有机闪烁体相关能点的系统响应函数Fig.1㊀Response functions for organic scintillator EJ -309and EJ -299-331.2.2㊀理论计算结合实验标定㊀㊀在有机闪烁体测量实验中,Dickens 等人[31-32]开发了基于蒙特卡罗算法的SCINFUL㊁NRESP7等程序进行闪烁体响应函数计算,但它们仅适用于圆柱几何,且无法描述复杂的粒子源分布㊂2007年,Pozzi 等人[33-34]基于MCNP 程序开发了MCNP -PoliMi㊁MCNP -PHOTRACK 等程序,其可对探测器的复杂几何结构进行描述且可描述复杂分布的中子源,其中MCNP -PoliMi 程序计算的系统响应函数如图3所示[33],由于其采用耦合计算的策略,后处理的工作量十分巨大㊂2014年,Hartwig 等人[35-37]建立了基于Geant4程序的有机闪烁体中子物理过程模拟方法,并成功应用于EJ -301有机闪烁体的系统能量响应函数及粒子甄别情况的模拟,模拟结果与实验结果吻合较好,如图4所示[37]㊂而在多球谱实验中,随着Monte-Carlo 计算方法的发展,采用MCNP 和Geant4程序进行多球谱㊃662㊃黄迁明等:中子能谱测量中的解谱技术研究进展㊀图2㊀采用加速器中子源标定多球谱仪能量响应方案Fig.2㊀Programme for using accelerator based neutron source to calibrate multi-spherespectrometer图3㊀采用MCNP-PoliMi 程序计算的系统响应函数Fig.3㊀System response functions usingMCNP-PoliMi图4㊀EJ -301有机闪烁体系统响应Geant4模拟Fig.4㊀Geant4based system response function for organic scintillator EJ -301仪系统响应函数计算的研究越来越多,图5(a)㊁5(b)为Mares 等人[22]采用MCNP 对不同热中子探测器6LiI(Eu)和3He 正比计数管多球谱仪能量响应函数的模拟结果㊂Geant4程序由于其具备对高能物理过程的模拟能力,近年来研究者们建立了基于Geant4程序进行多球谱仪系统能量响应函数的计算方法,图5(c )为Garny 等人[38]采用Geant4程序对10-11~100MeV 能量范围内多球谱仪能量响应函数的模拟㊂与有机闪烁体测量方法㊁多球谱仪测量方法㊃762㊃㊀辐射防护第42卷㊀第4期图5㊀多球谱仪系统响应函数MCNP 模拟以及Geant4模拟结果Fig.5㊀Multi-sphere spectrometer system response functions from MCNP and Geant4不同,多箔活化测量方法的能量响应函数为多箔活化材料的中子多群活化截面,能量响应函数的获取本质上为活化材料中子多群活化截面的制作㊂早在1969年,McElroy 等人[26]建立了用于SAND 程序解谱的640群活化截面㊂2008年,中国工程物理研究院的邓勇军等人[39]针对厚活化㊀㊀㊀㊀箔和活化箔包裹热中子吸收材料后多群截面的修正方法开展了相关研究,得到不同厚度(0.048mm㊁0.28mm)Au 活化箔640群多群截面和包裹0.5mm 厚度Cd 的Au 活化箔修正截面,理论和实验符合较好,如图6所示㊂图6㊀179Au 活化箔多群截面修正Fig.6㊀Correction for the multiple groups cross-section from179Au activated foil1.3㊀解谱不确定度㊀㊀中子能谱测量的误差主要来自两方面:一是输入数据的不确定度,二是用于描述物理过程的数学模型近似表示的不确定性以及算法本身缺少唯一解等,一般假设数学模型近似表示的不确定性可以忽略㊂Manfred Matzke [40]于2002年提出了针对最小二乘法和最大熵方法的中子能谱解谱误差传递方法,最小二乘算法可通过数学推导由先验谱㊁活化率㊁截面等信息的协方差数据精确给出中子能谱的不确定度,但该方法限制较多通用性不强,SAND 程序的适用性更强,但其中子能谱不确定度相对不完善,通常采用不考虑先验谱的影响和截面协方差的蒙特卡罗抽样方法,或将相关输入量的不确定度按相关性为零的假设进行处理,给出合成不确定度,其他解谱方法的误差处理方式类似㊂2009年,王松林等人[41]在采用多箔活化法测量Am-Be 中子源屏蔽辐照腔内的中子能谱时,考虑了初始谱引入的误差㊁解谱所用的截面误差㊂在活化箔片灵敏区不能覆盖或覆盖较弱的能区,最后的解谱对初始谱有很大的依赖,而在活化箔片灵敏区覆盖较好的能区,最后的解谱对初始谱依赖较小;不同的截面库之间存在微小差异,一般通过选取合适的截面库进行解谱,这部分误差可㊃862㊃黄迁明等:中子能谱测量中的解谱技术研究进展㊀忽略㊂2015年,陈晓亮等人[42]基于广义最小二乘法开发了NSAGLS程序并进行了误差分析,在考虑了输入谱㊁核反应截面及测量活度不确定度导致的误差后,解谱效果良好㊂2016年,李达等人[43]针对SAND-II程序解谱过程提出了一种基于先验谱㊁活化率和截面协方差的不确定度蒙特卡罗分析方法,首先建立基于线性变换的截面协方差抽样方法,然后利用MCNP计算误差,使用迭代方法估计先验谱的不确定度,最后结合活化率的测量不确定度,利用蒙特卡罗抽样方法计算中子能谱的不确定度,与传统方法计算的不确定度比较接近㊂2㊀主要解谱算法及程序㊀㊀中子能谱的解谱方法主要有:广义最小二乘算法㊁奇异值分解和正则化算法㊁最大熵算法㊁贝叶斯算法㊁蒙特卡罗算法㊁遗传算法㊁神经网络算法㊁MLEM算法以及基于压缩感知理论的中子能谱解谱方法,具体情况列于表1㊂表1㊀常用解谱程序介绍Tab.1㊀Neutron spectrum unfolding packages㊀㊀广义最小二乘算法由于发展成熟,目前已有大量解谱程序应用,其他解谱方法由于数学上的复杂性以及发展历史较短的缘故,基于其开发的解谱程序较少,而对于其他一些较新的方法,如神经网络算法以及压缩感知方法,目前研究者们仅对其建立了相应的算法,而未专门开展相关的程序开发㊂2.1㊀广义最小二乘算法㊀㊀早在1964年,Gold等人[44-45]就提出采用解谱计算值和实验测量值之间的最小二乘偏差作为求解目标,通过迭代策略的设置来保证解的非负性进行中子能谱解谱㊂2010 2015年,陈晓亮㊁孙征等人[46-47]基于广义最小二乘原理对迭代策略进行了修正并开发了程序NSAGLS和2NP㊂随后Chen等人[42]在有机闪烁体测量实验中采用GRAVEL算法进行了解谱研究,采用GRAVEL程序对测量结果进行了解谱处理,解谱结果与参考解基本吻合,如图7所示㊂Seghour等人[48-49]在2001年采用SAND-II程序对ANO例题进行了解谱研究,中国工程物理研究院在2014年开展了采用多箔活化方法对中物院某临界装置的中子能谱测量实验,并采用SAND-II程序对测量结果进行了解谱研究,得到不超过2%的偏差,证明了SAND-II解谱的准确度,如图8所示㊂广义最小二乘法通过将计算值和实验测量值之间的最小二乘偏差作为求解目标,从而将中子能谱解谱问题转化为围绕先验值的拟合问题,先验信息通常作为迭代的初值进行使用㊂目前应用最广的迭代策略为GRAVEL算法,以GRAVEL算法作为核心的解谱程序主要有SAND系列程序,主要包括SANDC㊁SAND㊁SAND-II㊁MSAND等㊂广义最小二乘解谱算法在有机闪烁体㊁多球谱仪和多箔活化中子能谱测量实验中进行了广泛的应用,其解谱结果与参考解吻合较好,对于闪烁体解谱的超定问题,该方法对与矩阵病态问题的适应能力尚有提高空间,解谱结果呈现出一定的震荡性;对于多球谱仪以及多箔活化解谱的欠定问题,求解准确程度对迭代的初值依赖性很大㊂㊃962㊃㊀辐射防护第42卷㊀第4期图7㊀GRAVEL 算法解谱结果Fig.7㊀Neutron spectrum unfolding results from GRAVELmethod图8㊀SAND-II 程序解谱结果Fig.8㊀Neutron spectrum unfolding results from SAND-II2.2㊀最大熵算法㊀㊀1989年,Weise 等人首次将 熵 的概念引入到解谱问题当中[50],并定义了解谱问题的熵函数,如(2)式所示:S =-ð{ϕjln(ϕj/ϕjDEF)+ϕj DEF -ϕj}(2)式中,ϕDEF为求解中子能谱的先验信息,在最大熵算法中称为预置谱;ϕj 为所求解的中子能谱,其中j 为中子能量分段数量;S 为所求的熵㊂1998年,Reginatto [51]采用基于最大熵算法的MAXED 程序对多球谱仪测量结果进行了解谱研究,多球谱仪由8个聚乙烯球组成,对普林斯顿大学等离子物理实验室的托卡马克聚变反应堆装置120m 处中子能谱进行了测量,如图9(a)所示㊂解谱结果靠近预置谱收敛,但由于测量值与理论值之间存在的偏差,造成解谱结果也与预置谱存在一定的偏差㊂Reginatto [52]于2002年针对大气层中~20km处宇宙射线激发中子的多球谱仪测量结果进行了解谱研究,Green 等人[53]于2008年采用MAXED程序对加速器驱动次临界装置的多箔活化测量结果进行了解谱研究并采用MCNP 程序对装置中子能谱进行了模拟,祝庆军等人[54]于2014年采用最大熵算法对多球谱仪测量结果开展了解谱研究,雄厚华等人[55]于2018年采用多箔活化方法对聚变包层DFLL-TBM 中子能谱进行测量,均证实MAXED 程序解谱效果与模拟值吻合较好,如图9(b)㊁9(c)㊁9(d)所示㊂最大熵算法通过求取约束条件下熵值最大的解来实现解谱,从而将解谱问题转化为优化问题,该方法可以给出非负㊁连续的中子能谱㊂最大熵方法在多球谱仪和多箔活化测量结果的解谱中应㊃072㊃黄迁明等:中子能谱测量中的解谱技术研究进展㊀图9㊀MAXED 程序解谱结果Fig.9㊀Neutron spectrum unfolding results from MAXED用较广,其对聚变㊁裂变以及放射性同位素中子源各种类型的中子能谱也具有较好的适应性㊂最大熵方法最大的特点为先验信息通过 预置谱 的方式在熵函数构造中进行使用,因此该方法又被认为是先验信息与实测信息相干性最小的解谱方法之一㊂2.3㊀贝叶斯算法㊀㊀2006年,Reginatto 等人[56]又提出基于贝耶斯分析方法进行中子能谱解谱,该方法以贝叶斯理论为基础,将中子能谱解谱问题转化为贝叶斯参数估计问题进行求解,并对多球谱仪测量结果进行解谱研究,解谱结果与Monte-Carlo 模拟结果吻合较好,相对偏差通常小于4%,如图10所示㊂2018年,Mazrou 等人[47]提出将中子注量率按热区㊁超热区以及快中子能量区间采用麦克斯韦分布㊁1/E 分布以及瓦特裂变谱的叠加:ϕ=ϕt +ϕe +ϕf(3)式中,ϕt=a tET 2e -ET 0,E ɤ0.1eV㊂图10㊀多球谱仪测量结果贝叶斯解谱结果Fig.10㊀Neutron spectrum unfolding results formulti-sphere spectrometer measurement usingBayes method ,vertical ordinate shows thenormalized neutron fluenceϕe=a e (1-e-E 2E2d)Eb -1e -Eβᶄ,0.1eV ɤE ɤ10keVϕf=a f E αe -Eβ,10keV ɤE ɤ20keV(4)式中,a t ㊁a e ㊁a f 分别代表热区㊁超热区㊁快中子能量区间能量峰的重要性;T 0为麦克斯韦分布最可几能量,取0.025eV;E d 为超热区最低能量,取㊃172㊃㊀辐射防护第42卷㊀第4期0.0707eV;b 和βᶄ分别为控制曲线上升和下降的斜率参数;α和β分别为描述快谱形状和峰值的参数㊂在他们的研究中采用贝叶斯算法对多球谱仪241Am -Be 源中子能谱测量结果进行解谱,如图11所示,采用贝叶斯算法解谱结果与GRAVEL㊁MAXED 等程序结果吻合较好,但在能量峰值上有细微差别,贝叶斯法中该值比MAXED 法㊁GRAVEL 法中稍低,这可能与是否设置预置谱和采用不同的解谱方法有关㊂采用贝叶斯算法解谱结果可以达到较高的精度,解谱结果中子总注量率和剂量率相对偏差小于1%㊂图11㊀多球谱仪测量241Am -Be 源中子能谱Fig.11㊀Neutron spectrum for241Am -Be source using multi-sphere spectrometer㊀㊀2017年,宋鸿鹄等人[57]采用贝叶斯算法对241Am -Be 源中子能谱有机闪烁体测量结果进行了解谱研究,解谱结果如图12所示㊂在形状上与ISO 的标准谱吻合较好,解谱不确定度与实际源的分布㊁ISO 选取有关㊂图12㊀241Am -Be 源中子能谱解谱Fig.12㊀Neutron spectrum unfolding results for241Am -Be measurement贝叶斯算法以贝叶斯理论为基础,通过将中子能谱进行参数化表征从而将解谱问题转化为基于贝叶斯理论的参数估计问题,目前的参数化表征方式主要有将中子注量率按热区㊁超热区以及快中子能量区间采用麦克斯韦分布㊁1/E 分布以及瓦特裂变谱的叠加表征方法㊂该方法经验性强,在先验较为充分的情况下可以取得比较精确的结果,但目前的参数化表征方式也限定了其使用范围㊂2.4㊀蒙特卡罗算法㊀㊀2007年,Bedogni 等人[58]提出了采用蒙特卡罗算法进行中子能谱解谱㊂在贝叶斯参数化方法的基础上,增加了蒸发谱和高斯谱模型对快中子能量分布进行表征,增加了蒸发谱对高能区中子能量分布进行表征,代表性的程序为FRUIT㊂Bedogni 等人[59]采用FRUIT 程序对241Am -Be 源中子能谱㊁252Cf 自发裂变中子能谱㊁252Cf(D 2O)中子能谱㊁12C 离子束碰撞中子源能谱㊁LINAC 放疗中子能谱以及CERF 装置中子能谱开展了解谱研究,如图13示㊂采用蒙特卡罗解谱可以取得较好的解谱效果,解谱结果与参考解吻合较好㊂对于241Am -Be 源中子能谱,解谱结果显得过于光滑,参考解中特征能峰部位有一定的偏差,如图13所示,表明蒙特卡罗算法可以取得较好的解谱精度,其解谱结果总注量率与参考解的最大偏差在5%以内㊂㊃272㊃黄迁明等:中子能谱测量中的解谱技术研究进展㊀图13㊀蒙特卡罗解谱方法求解效果Fig.13㊀Neutron spectrum unfolding results for Monte Carlo method蒙特卡罗中子能谱解谱算法在对中子能谱参数化的基础上,采用蒙特卡罗方法随机生成可表征中子能谱的参数集,然后依据该参数集表征的中子能谱和测量系统的响应函数计算探测器的计数并与实验值进行比较,重复该过程直至满足相应的收敛㊂该方法的使用具有很强的经验性,其优点为仅需指定中子场类型而无需提供初始能谱,其比较适用于无法通过理论计算获取初始谱的情况㊂由于采用蒙特卡罗方法作为迭代策略,该方法需要较长的计算时间,另外该方法解谱结果也有可能出现不符合物理意义的结果㊂2.5㊀其他解谱方法㊀㊀1999年,Freeman 等人[60]通过将探测器计数的计算值与实验值差值的平方和作为适应度函数,从而将解谱问题等价为可应用遗传算法进行求解的全局优化问题,然后采用标准遗传算法实现解谱㊂此后Mukherjee㊁王东等人[61-65]也进行了相关研究,但该算法始终存在精度较差和能谱不光滑特性㊂2002年,Braga 等人[66]首次利用神经网络技术对多球谱仪测量结果进行了中子能谱解谱,其采用三层神经网络结构,应用于中子能谱解谱的神经网络结构均比较简单且需要大量样本进行训练,在所求解中子能谱能群数目较少的情况下可以给出比较好的解,但也容易产生 过拟合 或 拟合不足 等学习现象㊂2013年,Pehlivanovic 等人[67]采用ML-EM 算法对特征峰中子能谱㊁252Cf 中子能谱的有机闪烁体测量结果进行了解谱研究,解谱结果与参考解基本吻合,但结果不够平滑,出现了较多不属于参考解的错误峰值,如图14所示㊂2017年,Molina 等人[68]使用ML-EM 方法对RECH -1反应堆中子能谱的多箔活化测量结果进行了解谱研究,如图15所示,结果表明在先验信息充足的情况下,ML-EM 方法可以取得较好的解谱结果㊂2004 2006年间由Candes㊁Romberg 和Tao等人[69-73]提出了压缩感知理论,其被认为是对经典采样理论Nyquist-Shannon 定理的突破和补充,已广泛应用于信号处理,医学㊁雷达图像重建,数据压缩㊁数据传输等领域㊂刘斌等人[74]于2019年建立了基于压缩感知理论的中子能谱解谱方法,对于核反应堆屏蔽结构典型中子能谱,解谱结果与参考解的总注量率相对偏差在4%以内,如图16所示,对于反应堆辐照监管中多处位置的中子能谱解谱结果与参考解相对偏差小于2%,说明其对高欠定程度的解谱问题也有较好㊃372㊃㊀辐射防护第42卷㊀第4期图14㊀闪烁体中子能谱Fig.14㊀Neutron spectrum fromscintillator图15㊀RECH -1反应堆Fig.15㊀Reactor RECH -1适应性㊂由于研究起步晚,遗传算法㊁神经网络算法等新兴算法的研究成果还比较少,在解谱稳定性或精度方面不及最小二乘算法等经典算法,但由于新兴算法具有较大的解谱精度和扩展潜力,若能解决现存问题并达到可应用的程度,遗传算法㊁神经网络算法等将有望在中子能谱解谱中发挥重要作用㊂3㊀总结和展望文章总结了中子能谱解谱模型㊁响应函数获取方法,以及解谱过程中的误差产生㊁处理方法,重点介绍了国内外中子能谱测量技术研究现状以及中子能谱解谱算法研究现状,总结了不同解谱算法的特点,接着介绍了根据不同解谱算法发展的解谱程序,对比了不同解谱算法及程序的优缺点,并推介了目前最适合用于中子能谱解谱的程序㊂针对解谱方法的原创性研究国外研究者占大多数,国内学者的研究主要集中在跟踪研究㊁已有算法应用改进以及依据已有算法开展解谱程序的研发等方面,对于解谱方法本身突破性的研究,尚有较大的提升空间㊂现有解谱程序主要依赖于发展较早㊁较成熟的算法开发,基于最小二乘算法的程序最多,包括SAND 系列程序㊁NSAGLS㊁ANNs 等,基于最大熵算法的有MAXED 程序,基于蒙特卡罗算法的有FRUIT 程序,功能较强㊁应用最广泛的还属SAND 系列程序和MAXED 程序㊂而对于较新的算法,如贝叶斯算法㊁正则化算法虽然已有较多研究,解谱效果也一般,尚无通用的程序公开,其他诸如遗传算法㊁神经网络算法等尚处于原理研究阶段,而遗传算法㊁神经网络算法具有功能强大㊁扩展性强的特点,若能解决现在遇到的一些问题,这些新兴算法将在中子能谱解谱中发挥重要作用,这方面还有很多工作可以做,包括解决遗传算法的解谱结果不光滑问题㊁优化算法提高求解速度等㊂经过梳理可发现,目前大多数解谱程序都是㊃472㊃图16㊀几种典型中子能谱解谱效果Fig.16㊀Unfolding results for several typical neutron spectrum采用单一解谱算法,尚无综合多种解谱方法的解谱程序研发,而就目前对解谱效果的调研来看,每种解谱方法对特定的中子能谱类型往往具有适用性,比如ML-EM方法对特征峰类型的能谱往往具有很好的适应性,正则化算法对于平滑的中子能谱具有很好的适应性,尚无对所有中子能谱㊀㊀㊀㊀㊀都有很好适应性的解谱方法㊂因此,开发包含多种解谱方法的综合性解谱程序具备较强的应用需求,相信随着中子能谱解谱技术的发展,必将在各类核装置的设计和运行㊁辐射防护剂量学和屏蔽验证以及核军备控制和反恐等领域发挥重要作用㊂参考文献:[1]㊀Ikuo Kondo,Kiyoshi Sakurai.Experimental evaluation of reactor neutron spectrum in irradiation field[J].Journal ofNuclear Science and Technology,1981,18(6):461-472.[2]㊀Griffin P J,Kelly J G,Luera T F,et al.Effect of new cross section evaluations on neutron spectrum determination[J].IEEE Transactions On Nuclear Science,1992,39(6):2078-2085.[3]㊀DavideChiesa,MassimilianoNastasi,Carlo Cazzaniga,et al.Measurement of the neutron flux at spallation sources usingmulti-foil activation[J].Nuclear Instruments and Methods in Physics Research A,2018:14-24.[4]㊀Klein H.Workplace radiation filed analysis[J].Radiation Protection Dosimetry,1997,70(1):225-234.[5]㊀Flaska M,Pozzi S A.Identification of shielded neutron sources with the liquid scintillator BC-501A using a digital pulseshape discrimination method[J].Nuclear Instruments and Methods in Physics Research A,2007,577:654-663. [6]㊀Pehlivanovic B,Avdic S,Marinkovic P,et parison of unfolding approaches for monenergetic and continuous fast-neutron energy spectra[J].Radiation Measurements,2013,49:109-114.[7]㊀Horst Klein,Sonja Neumann.Neutron and photon spectrometry with liquid scintillation detectors in mixed fields[J].Nuclear Instruments and Methods in Physics Research A,2002,476:132-142.[8]㊀Guillaume H V,Matthieu H,Normand S,et al.Pulse shape discrimination between(fast or thermal)neutrons and㊃572㊃。

脉冲中子全谱测井技术及其在冀东油田的应用

脉冲中子全谱测井技术及其在冀东油田的应用

脉冲中子全谱测井技术及其在冀东油田的应用吕俊涛【摘要】介绍了脉冲中子全谱测井仪(PNST)的结构、特点、测量原理及技术指标。

脉冲中子全谱测井能同时完成双源距碳氧比、中子寿命、脉冲中子-中子和能谱水流测井功能。

其测井资料不但能准确评价地层剩余油气饱和度、判断出水层位,还能识别低电阻率油层、气层及评价压裂效果等。

该仪器在冀东油田已应用27井次,依据PNST测井资料增油降水措施有效率达95%以上。

冀东油田现场应用的统计分析表明,该仪器具有测量剩余油饱和度精确性高,应用范围广等特点,能为高含水后期油藏监测,为堵水、补孔、压裂等方案的制订提供有效的资料。

【期刊名称】《石油管材与仪器》【年(卷),期】2016(000)003【总页数】5页(P53-57)【关键词】脉冲中子全谱测井仪;剩余油饱和度;压裂评价;气层识别;动态监测【作者】吕俊涛【作者单位】大庆油田有限责任公司测试技术服务分公司,黑龙江大庆163453【正文语种】中文【中图分类】P631.817冀东油田进入高含水开发后期,综合含水在90 %以上,测准地层的剩余油情况至关重要,影响到二次开发的成败。

目前,研究剩余油主要是通过完井的水淹层测井解释资料及生产井的各种动态监测资料来进行综合分析,但由于水淹层测井解释资料只能反映完井时的水淹状况,而动态监测资料虽能反映目前油层的动用状况,但大部分资料都是合层测试,并不能完全反映单个油层的动用状况,而且并不是所有的井或层都具有这种资料。

因此,只凭这些资料很难弄清每个小层剩余油分布状况,从而给高含水后期剩余油潜力分析带来了难点。

而地层参数测井的首要任务是判断油层水淹状况、高含水层、在老井中寻找潜力油层,确保了油田的稳产和提高油田的开发效益。

目前国内地层参数测井主要是过套管电阻率、PNN、C/O等,这些技术都有测量参数单一、测量精度不高和使用范围窄等问题,并且解释符合率不是很高。

脉冲中子全谱测井仪(PNST)通过合理的传感器结构设计和电路设计,优化综合时序发生器、自动稳谱等技术,一次下井能同时完成双源距碳氧比、中子寿命、脉冲中子-中子、能谱水流等4项测井功能。

应用于中子照相的热中子敏感微通道板

应用于中子照相的热中子敏感微通道板

应用于中子照相的热中子敏感微通道板潘京生;顾燕;孙建宁;苏德坦;陆强【摘要】介绍了中子成像探测技术的应用.由于直接在微通道板(MCP)玻璃中掺加中子灵敏核素可使MCP对热中子敏感,从而可将MCP事件计数成像探测器的优势成功地应用于以中子为探针的成像探测技术,本文开展了热中子敏感微通道板的研究.通过在MCP玻璃中掺摩尔百分比为3%的Gd2O3,并沿用MCP的制作方法,完成了直径为50 mm和106 mm的大面阵热中子敏感MCP的制作,并进行了基于这种大面阵热中子敏感MCP的中子事件计数成像探测器的中子成像实验.理论和实验结果都验证了掺摩尔百分比为3%的Gd2O3 MCP可取得对热、冷中子30%~50%的探测效率.最后,进一步介绍了目前开展的封装式中子敏感MCP增强管的研制工作.基于掺Gd2O3的MCP增强管并经光学耦合到CCD或CMOS相机的紧凑混合传感器结构是实现高时空分辨能力的中子照相无损检测技术的有效途径.【期刊名称】《光学精密工程》【年(卷),期】2015(023)003【总页数】7页(P660-666)【关键词】中子照相;中子相机;中子敏感微通道板;增强CCD;增强CMOS-APS 【作者】潘京生;顾燕;孙建宁;苏德坦;陆强【作者单位】北方夜视技术股份有限公司,江苏南京211102;北方夜视技术股份有限公司,江苏南京211102;北方夜视技术股份有限公司,江苏南京211102;北方夜视技术股份有限公司,江苏南京211102;北方夜视技术股份有限公司,江苏南京211102【正文语种】中文【中图分类】TL816.3;TN151.11 引言在凝聚体物理学、材料科学和工程、地球科学以及生物科学等多个基础科学领域,中子成像已成为一种有效的探测方式。

中子相互作用概率与核结构的相关性决定了不同材料对中子的吸收系数经常有着数量级的差别,因此中子成像可获得更高的分辨力和灵敏度,而通过中子计算机断层成像进一步保留中子吸收的空间分布信息,对其物理结构按照实际的空间相对位置和方向进行可视化成像,有助于理解材料缺陷发展和器件失效的基本机制[1-4]。

2012.10.18中子的测量方法和探测器1

2012.10.18中子的测量方法和探测器1

的中子灵敏度可达 :
2.3×10-14 A/(n.cm2.s) 电离室高压电极(负极) 的内壁 及收集极外壁涂硼,室内充1%的氦 和6%的氮,93%的氩气。 由于中子产生的电流较大,γ的影 响较小。
9
(2)正比计数器
10B
(n, α) 7Li
输出脉冲幅度为:
A是正比计数器的气体放大倍数通常可达到103 ~105 ,En是中子的动能;Q是反应能;e是电子电荷量;W 是 平均电离能;C是计数器的等效输出电容。 主要用于热中子的测量。
16
几种发射体特性数据表
中子截面 ( b ) 150 5.1 37 响应时间 T 1/2 42s 3.7 min 瞬时 反应 (n,γ)β ̄ (n,γ)β ̄ (n,γ)γ
发射体
铑(Rh) 钒(V) 钴(Co)
密度
g/cm3
相对灵敏 度 铑=100
100 17.5 1.5
1014通量 燃耗 %每年
32 1.6 11
能量为E n 的中子微观截面σ= σ0(1/υ)
2
2、核反冲法 中子与物质原子核发生弹性碰撞,原子核被反冲,且带一定正电荷,选用反冲核 弹性碰撞截面大的材料作为探测器灵敏物质,就可以简接测量中子的注量率。通常是 利用含氢物质作为灵敏体。 反冲核的反冲能表示为:
EA
4A 2 E cos n 2 (1 A)
各种闪烁体测量中子的性能对比
名 称 型 号 规 格 Φ32,50
Φ50×20 Φ50×1.5 Φ50×1.5 Φ40×10 Φ10~50 探测中子 效 率%
光衰时间 0.2μS 0.2μS 0.2μS 0.2μS 0.1μS 30 ns
光谱峰(1010m) 4500 4500 4500 4500 3950 4470

中子能谱测量方法

中子能谱测量方法

中子能谱测量方法
中子能谱测量可是个很有趣又有点复杂的事儿呢。

一种常见的方法是活化法。

这就像是给中子找个小助手,让它去和一些物质发生反应。

比如说,把某些特定的箔片放在有中子的环境里。

中子就会和箔片里的原子核相互作用,让原子核变得“兴奋”起来,也就是被活化啦。

然后呢,我们再通过测量活化后的原子核发射出的射线的能量和强度等信息,就可以倒推出中子的能谱啦。

这就像是根据一个人留下的脚印,去推测这个人的身高体重啥的一样,很神奇吧。

还有时间飞行法。

想象一下中子就像一群调皮的小飞鸟。

我们先让中子从一个地方出发,然后在另一个地方测量它们到达的时间。

因为不同能量的中子飞行速度不一样呀,就像不同体力的小鸟飞的快慢不同。

能量高的中子飞得快,能量低的中子飞得慢。

通过精确测量到达时间,我们就能把不同能量的中子区分开来,从而得到中子能谱。

这个方法就像是给中子们办了一场飞行比赛,根据到达终点的先后顺序来判断它们的能量高低。

多球谱仪法也很厉害哦。

就像是用好多不同的小笼子去捕捉中子。

这些小笼子就是不同大小、不同材料的球。

中子进入球里就会和里面的物质相互作用。

每个球对不同能量的中子反应不一样,就像每个小笼子对不同类型的小动物有不同的捕捉能力一样。

然后我们综合分析各个球的反应情况,就能知道中子的能谱啦。

这些测量方法都各有各的巧妙之处,科学家们就像聪明的侦探一样,通过这些方法去揭开中子能谱这个神秘的面纱呢。

不管是活化法的巧妙反应,时间飞行法的速度比拼,还是多球谱仪法的综合探测,都让我们对中子的能量分布有了更多的了解。

试管内样品中子吸收剂量的测定

试管内样品中子吸收剂量的测定
( 中国Z 程物理研究院核物理与化学研究所,四川绵阳, 621900; 国家核技术工业应用工程技术研究中心,四川绵阳, 621900
摘要: 根据吸收剂量与中子注量的关系。 采用活化法测量了试管内样品的中子吸收剂量。实验表明 采用具有反应阑的活化箔能较好地监测少量含氢材料中的中子注量,从而确定样品中的中子吸收剂量。
中子能谱改变不大, 吸收剂量与中子注量仍然呈
线形关系。
3
实验过程
中 子注量率采用2 9 裂变电 直接测量。 3 Pu 离室
由于辐照样品多, 注量点多, 采用电离室监测每 个位置每个注量点的做法不可行。实际上以电离 室的测量结果作为参考标准, 采用活化箔来监测
所有测点的中子注量。

2
吸收剂量的确定
吸收剂量D定义为: 任何电离辐射授予物质 的平均能量除以物质的质量。
关键词: CFBR-1 1堆; 吸收剂量; 中 子注量; 活化箔
中图分类号: TL375. 4 文献标识码: A
引 言 为了研究强辐射对生化制剂的辐照效应和 破坏机理, CFBR-II '决中子脉冲堆上进行了生 在 化制剂的中子辐照实验。生物制剂样品有短小芽 抱杆菌和枯草芽抱杆菌样品; 化学制剂样品为有 机磷酸醋 DMMP, TEP, TMP 几种样品。 样品大 部分为液体, 有少量固 体粉末, 均封装在直径 10 mm, 高度 40 mm 的玻璃或聚醋试管中。 生化制剂的辐照效应主要取决于吸收剂量。 当样品和中子能谱确定后,中子吸收剂量则取决 1
随距离的增加而迅速减小, 但是 Au 活化箔的比
计数率完全背离了这一趋势。这说明 Au 活化箔 对于低能中子的反应截面太大, 少量低能中子的 增加就会引起计数率的巨大改变, 所以不能应用 Au 活化箔来测量中子注量。 S 活化箔是闭探测器, 其阂值为3 MeV 左右, 低能中子的增加不会对测量结果产生影响。放置 在试管前的活化箔由于反射中子的存在,比计数 率应稍有增大, 但影响不大。试管后面的活化箔 由于中子的吸收、散射和反射的影响,比 计数率 要减小,此影响要大些。从表 2 可以看出,将比 计数率由试管前推到试管后, 所有外推数据均比 实验数据要大, 这符合实际情况。 所以中心中子 注量采用二者的平均值。 选择放置在液体试管前的 S 活化箔数据, 做

重水反应堆热柱中子能谱测量

重水反应堆热柱中子能谱测量

重水反应堆热柱中子能谱测量
包宗渝;陈军
【期刊名称】《青岛大学学报:自然科学版》
【年(卷),期】1997(010)002
【摘要】利用飞得时间法,用机械选择器测量了中国原子能科学研究院重水反应堆热柱的中子能谱。

结果表明,该热柱能谱和理论Maxewll谱有偏离。

【总页数】4页(P66-69)
【作者】包宗渝;陈军
【作者单位】中国原子能科学研究院;中国原子能科学研究院
【正文语种】中文
【中图分类】TL423.075
【相关文献】
1.核反应堆内中子能谱测量技术 [J], 王谷军
2.游泳池反应堆热柱中子束空间分布.能谱及剂量参数的测定 [J], 陈常茂;谢建伦
3.SPRR-300反应堆大热柱内中子注量率及能谱分布 [J], 窦海峰;代君龙
4.研究性重水反应堆厂房外环境中子、γ剂量当量率的测定 [J], 陈常茂;闻友勤
5.重水反应堆热柱中子能谱测量 [J], 包宗渝;陈军;岳骞;徐昆;
因版权原因,仅展示原文概要,查看原文内容请购买。

中子散射技术及其应用

中子散射技术及其应用
2. 3 中子散射的特点 中子散射技术中的许多实验方法 , 如中子衍
射 、中子小角散射 、中子反射技术等都是在相应的 X 射线实验方法的基础上发展起来的.
为什么有了 X 射线分析方法还需要发展中子 散射技术呢 ? 原因在于 , 中子散射有自己的特点 , 这些特点恰好弥补了 X射线分析固有的缺点 , 从而 使它和 X射线分析形成了互补的关系. 这些特点可 以归纳为 : (1) X射线对原子序数低的轻元素不灵 敏 ,但中子对轻 、重元素的灵敏度没有明显的差别 ; (2) X射线不能分辨原子序数相近的元素 , 而中子 通常可以分辨 ; ( 3) 中子可以区分同位素 ; ( 4 ) 中 子具有磁矩 ,因而可以研究磁性物质的磁结构和自 旋动力学. 常规 X 射线分析不能提供磁的信息. 近 年来虽然己经可以用同步辐射来研究物质的磁结 构 ,但中子作为微观磁结构的研究工具仍然是其他 方法无法代替的 ; ( 5) 中子对物质有较强的穿透能 力 ; (6) X射线只能研究物质的静态结构 ,不能研究 动力学问题 ,这是因为波长在 0. 1—1 nm 左右的 X 射线 ,其能量比原子 、分子的运动能量高几十万倍 , 所以不可能用它来研究物质的微观动力学特性.
35卷 (2006年 ) 11期 http: ΠΠwww. wuli. ac. cn
Hale Waihona Puke ·961·前沿进展
其中有些反应堆开始用于基础研究. 1946 年 ,美国 Oak R idge实验室的科学家 Wollan, Shull等首先在 Clinton反应堆上开展了中子衍射工作. 这个工作代表 了中子散射的一个方面 , 即用中子散射方法研究物 质的静态结构 ; 1950年前后 ,加拿大 Chalk R iver实验 室科学家 B rockhause开始用中子非弹性散射研究晶 格动力学 , 开创了中子散射技术的另一个方面 , 即利 用中子散射研究物质的微观动力学性质. 通常说的 “中子散射技术 ”是这两方面工作的总称.

中子能谱特征

中子能谱特征

中子能谱特征
中子能谱是研究核离子化过程中中子的能量分布的一种方法。

通常
情况下,中子能谱的测量是通过将中子束传输到样品中并收集散射反
应中的中子来实现的。

根据收集到的中子能量,可以绘制出中子能谱。

中子能谱具有一些特征:
1. 能峰
能峰是中子能谱中最明显的结构之一,表示中子在样品中被吸收的能量。

能峰的位置反映了特定核反应的阈值和强度。

2. 背景
背景是指中子计数与能峰之间的连续谱线。

它通常由样品对中子的散
射造成,并且不具有明显的峰值结构。

3. 反冲峰
反冲峰是一种在中子能谱中出现的有趣结构。

它是由于中子在样品中
被散射并与其他核反应重新交互后引起的。

反冲峰可以用于确定中子
在样品中的散射截面。

4. 能谱分辨率
能谱分辨率指能够分辨两个紧密而不同的能量间隔的程度。

它取决于测量系统的几何形状和探测器的性能。

能谱分辨率越高,能够分辨出的核反应中心的能量就越精确。

中子能谱是研究核反应、核能和核物理的重要方法之一。

通过测量中子能谱,可以了解核反应的能量和强度,为核能利用和核武器控制等领域提供重要参考。

同时,中子能谱的研究还有助于进一步了解中子的相互作用和能量分布规律。

飞行时间效应对含氢系统瞬发中子衰减常数测量的影响

飞行时间效应对含氢系统瞬发中子衰减常数测量的影响
tonsf r t r 1n uton s nstv e e t r r o he ma e r e ii e d t c o . Ke r s pr m pt n ton; s c iiaiy; fi httm e sow i own; t e m a n u y wo d : o eu r ub rtc lt lg - i ; l ng d hr l e— to r n;Ros i m e h d s— to
I fu n e o u r n Flg — m e Ef e to o ptNe t o n l e c f Ne t o i htTi f c n Pr m ur n
De a n t n e s r m e s f r Hy r g n u a t r c y Co s a tM a u e nt o d o e o sRe c o s
S ONG i g l L n —i H ENG u ,Z Ch n,JANG n I Yo g,L U a — o I Xio b
(nttt f cerP y k n kmi r Isiueo Nu la h s Sa d C e s y,C iaAcd myo g neig P y i ,M in a g 6 1 0 ,C ia t hn a e f En iern h s s c a y n 2 9 0 hn )
响; 对热 中子 灵 敏 的 .玻 璃 探 测 器 , 过减小 探测距 离 和屏蔽低 能 中子可 降低 中子飞 行时 间效应 的影响 。 I i 通 关 键 词 : 发 中 子 ; 临界 ; 瞬 次 飞行 时 间 ; 化 ; 中子 ; o s 口法 慢 热 R si 一 中 图分 类号 : 7 . O5 1 1 文献标志码 : A 文章 编 号 :0 0 6 3 ( 0 0 0 — 9 60 1 0 — 9 1 2 1 ) 80 0 —4

中子与物质的相互作用

中子与物质的相互作用

10-8 – 10-9
65Cu
31%
10-9 – 10-10
121Sb
57%
10-9 – 10-10
* 测量1小时 放射性计数100 探测效率10%
应用2:中子测井
中子—中子
中子—伽马
中子 测井
中子寿命
中子活化
利用中子与物质相互作 用的各种效应
研究钻井剖面岩层性质 的一组方法
中子-中子测井
中子源放出快中子,经过一系列碰撞而减弱到热能状态, 最后被吸收 中子的空间分布与物质的减速及吸收性质密切相关。
裂变碎片的总动能为150~170MeV,形成的脉 冲幅度比 本底脉冲幅度大得多,可用于强辐 射场内中子的测量。
热中子可引起的核裂变核:233U,235U,239Pu。 如235U的热中子截面为580b。对慢中子满足1/v 规律,仅适用于热中子的注量率测量。
一些重核只有当中子能量大于某一阈能才能 发生核裂变,可用此判断中子的能量。
2 中子的分类与性质
(1)中子的分类 1) 慢中子:0~1KeV。包括冷中子、热 中子、超热中子、共振中子。
热中子:与吸收物质原子处于热平衡状态, 能量为0.0253eV,中子速度~2.2×103m/s.
2) 中能中子:1KeV~0.5MeV。
3) 快中子:0.5MeV~10MeV。
4) 特快中子:>10MeV。
24He
01n
17.6MeV (Q-value)
Notes :
1) 1 mA D (accelerated potentiao 100~300keV) → 109n/s (D-D), 1011n/s(D-T). 2) Other neutron generators (higher energy): 9Be(d,n), 7Li(p,n), 3H(p,n)

活化法测量散裂靶中子能谱的实验验证

活化法测量散裂靶中子能谱的实验验证

活化法测量散裂靶中子能谱的实验验证曾丽娜;王强;艾自辉;郑春【摘要】散裂靶中子的能谱对加速器驱动次临界系统的倍增因数和嬗变率等影响很大,计算表明散裂靶中子谱在M eV能区与裂变中子谱相近。

本文利用活化法测量临界装置的泄漏中子谱和中子注量率,提出了用In、Al、Mg、Ti、Au、Zn、Ni、Rh、Fe和Co等活化箔测量散裂靶中子能谱和中子注量率的方案。

结果表明,将活化箔在散裂靶中子场中辐照5 h ,中子注量最高达5×1014 cm-2量级,辐照后1 h内取出活化箔,根据半衰期的长短安排测量顺序,可测量散裂靶的中子能谱和中子注量率。

%The spectrum of spallation neutron source is of great influence on the multi‐plication factor and transmutation rate of ADS .It is show n that the spectra of spallation target neutron and fission neutronare similar in MeV region .In this paper ,the leakage neutron spectrum and neutron fluence rate of the critical assembly w ere measured by activation method with In ,Al ,Mg ,Ti ,Au ,Zn ,Ni ,Rh ,Fe and Co foils as activationfoils .T he results show that the activation method can be used to measure the neutron spectrum in the spallation target neutron field .T he neutron fluence should be up to 5 × 1014 cm -2 in 5 h .The activation foils should be taken out in 1 h after irradiation .The measure order is according to the length of the half‐life .【期刊名称】《原子能科学技术》【年(卷),期】2015(000)004【总页数】5页(P700-704)【关键词】活化法;散裂源;中子能谱;ADS【作者】曾丽娜;王强;艾自辉;郑春【作者单位】中国工程物理研究院核物理与化学研究所,四川绵阳 621900; 中国工程物理研究院中子物理学重点实验室,四川绵阳 621900;中国工程物理研究院核物理与化学研究所,四川绵阳 621900; 中国工程物理研究院中子物理学重点实验室,四川绵阳 621900;中国工程物理研究院核物理与化学研究所,四川绵阳621900; 中国工程物理研究院中子物理学重点实验室,四川绵阳 621900;中国工程物理研究院核物理与化学研究所,四川绵阳 621900; 中国工程物理研究院中子物理学重点实验室,四川绵阳 621900【正文语种】中文【中图分类】TL326加速器驱动次临界系统(ADS)通过加速器产生高能质子轰击靶,产生的中子注入一次临界系统中使其倍增。

利用测热技术测量核反应堆中子通量密度

利用测热技术测量核反应堆中子通量密度

利用测热技术测量核反应堆中子通量密度杨铜锁;史永谦;朱庆福;杨维平;张建祥;王藩;王团部;鲁瑾;李来冬【摘要】一种新型中子探测器被研究,其原理是利用带电离子在矿物中沉积的能量退火时会以热量的方式释放出来,通过测量释放的热量而确定中子通量密度。

对新型中子探测器进行刻度,在反应堆内某位置测量的热中子通量密度为5.108×1011 cm -2・s-1,与标定的热中子通量密度(5.000×1011 cm -2・s-1)在2%内符合,说明该探测器可测量中子通量密度。

本文方法制作的探测器体积小,可制作成不同形状,便于反应堆不同环境下的中子通量密度测量。

选取相应中子能量反应截面较大的元素,该探测器还可测量不同中子能量的通量密度。

%One new neutron flux density detector was studied .Its principle was to mea‐sure annealing heat from mineral depositional energy to make sure neutron flux density . The new neutron detector was calibrated ,the measuring neutron flux density of certain location in reactor was 5.108 × 1011 cm -2 ・s-1 ,and the measuring value and calibrated neutron flux density (5.000 × 1011 cm -2 ・s-1 ) were coincidence in the 2% .The results show that the new neutron detector can use to measure neutron flux density . T his neutron detector has the advantages of small volume possession in nuclear reactor with required shape , and simplified multi‐point measurement . The detector can measure different energy neutron flux densities if select element is sensitive for corresponding energy neutron .【期刊名称】《原子能科学技术》【年(卷),期】2016(050)009【总页数】5页(P1629-1633)【关键词】中子探测器;退火热量;热中子通量密度【作者】杨铜锁;史永谦;朱庆福;杨维平;张建祥;王藩;王团部;鲁瑾;李来冬【作者单位】西京学院应用统计与理学系,陕西西安 710123;中国原子能科学研究院反应堆工程研究设计所,北京 102413;中国原子能科学研究院反应堆工程研究设计所,北京 102413;陕西师范大学民族教育学院,陕西西安 710062;西京学院应用统计与理学系,陕西西安 710123;中国原子能科学研究院反应堆工程研究设计所,北京 102413;西京学院应用统计与理学系,陕西西安 710123;中国原子能科学研究院反应堆工程研究设计所,北京 102413;中国原子能科学研究院反应堆工程研究设计所,北京 102413【正文语种】中文【中图分类】O571.53中子通量密度是核反应堆的重要参数,直接反映了核反应堆裂变反应速率的快慢。

中子检测方法

中子检测方法

中子检测方法中子检测的方法有很多种,主要包括核反冲法、核反应法、核裂变法、活化法等。

以下将简单介绍这些方法:1. 核反冲法:通过测量中子与原子核弹性散射后的反冲核在探测介质中引起的电离来反推原始中子的性质,这种方法只适用于探测快中子。

2. 核反应法:通过测量中子核反应产生的带电粒子来探测中子的方法,主要用于探测慢中子的强度,也可用来测定快中子的能谱。

3. 核裂变法:中子轰击重核时会引起核裂变,通过探测裂变碎片来探测中子的方法,常用235U,233U和239Pu作为裂变材料,主要用于热中子和慢中子的通量测量,强γ射线本底对测量也不会造成响。

4. 活化法:中子与原子核相互作用时生成复合核,一般复合核因处于激发态而不稳定,具有γ或β放射性,通过测量中子辐照后物质的放射性来推算中子强度。

5. 中子俘获法:利用中子与原子核相互作用的性质,当中子遇到原子核时,有时会被原子核捕获并变成一个新的核,同时释放出一些能量,利用这个方法可以观测中子在原子核中行进的路径和范围。

6. 热中子散射法:利用中子被物质散射的性质,将中子引入样品后,用探测器测量出中子的散射角度和能量分布,从而确定中子的性质。

7. 中子衍射法:通常使用晶体衍射技术,中子与晶体中原子核相互作用后,中子波的相位和振幅会发生变化,从而产生衍射图样,通过分析衍射图样可以了解中子的波长、动量和散射角度等信息。

8. 核磁共振法:利用了中子和原子核间的相互作用,当中子穿过样品时,其与样品中原子核的相互作用产生能量变化,从而激发了原子核中的自旋,进而观察到核磁共振信号。

9. 时间飞行法:通过测量中子到达探测器的时间差来确定它们通过样品的时间,从而推断中子的性质和路径。

以上内容仅供参考,如需获取更多信息,建议查阅相关文献或咨询专业人士。

核反应堆物理分析习题答案-第二章

核反应堆物理分析习题答案-第二章

1、 H 和O 在1000到1能量范围内的散射截面似为常数,分别为20b 和38b.计算2H O 的ξ以与在2H O 和中子从1000慢化到1所需要的碰撞次数。

解:不难得出,2H O 的散射截面与平均对数能降应有下列关系:222H O H O H H O O σξσξσξ⋅=⋅+⋅即2(2)2H O H O H H O O σσξσξσξ+⋅=⋅+⋅2(2)/(2)H O H H O O H O ξσξσξσσ=⋅+⋅+查附录3,可知平均对数能降: 1.000H ξ=,0.120O ξ=,代入计算得:2(220 1.000380.120)/(22038)0.571H O ξ=⨯⨯+⨯⨯+=可得平均碰撞次数:221ln()/ln(1.0001)/0.57112.0912.1C H ON E E ξ===≈2.设()f d υυυ''→表示L 系中速度速度υ的中子弹性散射后速度在υ'附近d υ'内的概率。

假定在C 系中散射是各向同性的,求()f d υυυ''→的表达式,并求一次碰撞后的平均速度。

解: 由: 212E m υ'='得:2dE m d υυ'=''()(1)dE f E E dE Eα'→''=-- E E E α≤'≤()f d υυυ''→=22,(1)d υυαυ''-- υυ≤'≤()f d υυυυυ='→''322(1)3(1)υαα=--6.在讨论中子热化时,认为热中子源项()Q E 是从某给定分解能cE 以上能区的中子,经过弹性散射慢化二来的。

设慢化能谱服从()E φ/E φ=分布,试求在氢介质内每秒每单位体积内由c E 以上能区,(1)散射到能量为()c E E E <的单位能量间隔内之中子数()Q E ;(2)散射到能量区间1g g g E E E -∆=-的中子数g Q 。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第38卷增刊原子能科学技术Vol.38,Suppl. 2004年7月Atomic Energy Science and TechnologyJ uly 2004含氢介质内中子能谱测量安力,陈渊,郭海萍,牟云峰,王新华,朱传新(中国工程物理研究院核物理与化学研究所,四川绵阳 621900)摘要:建立了直径34cm 的含氢慢化球和含氢慢化球与<24cm ×30cm 聚乙烯圆柱组合的2种基准装置,加速器的d 2T 中子入射到含氢慢化介质,用<18mm ×20mm 的 晶体闪烁探测器测量了2种实验装置内不同位置的1MeV 以上的中子能谱,并处理成不同能量阈值的中子数。

在0.95置信水平下,本测量方法的不确定度为4.8%。

关键词:d 2T 中子;二维装置;介质内能谱中图分类号:O571.54 文献标识码:A 文章编号:100026931(2004)S020089204Measurement of N eutron Spectrum in Medium Containing H ydrogenAN Li ,CHEN Yuan ,GUO Hai 2ping ,MOU Yun 2feng ,WAN G Xin 2hua ,ZHU Chuan 2xin(Institute of N uclear Physics and Chemist ry ,China Academy of Engineering Physics ,P.O.Box 9192213,Mianyang 621900,China )Abstract : Two benchmarks were established.One is moderation shell containing hydrogen which is 34cm in diameter.The other is composed of the shell and polyethylene cylinder of<24cm ×30cm which are combinatorial.Measurement of neutron energy spectrum above 1MeV changing with different positions of the experimental assemblies were carried out for in 2cident neutrons from outside using d 2T fusion source of accelerator.Meanwhile ,the spectra were transformed integral spectrum.The scintillation detector of stilbene crystal of <18mm×20mm was used to measure neutron spectrum.At 0.95level of confidence ,the uncertain 2ty of the measurement is 4.8%.K ey w ords :d 2T neutron ;two dimensional assembly ;neutron spectrum in medium 收稿日期:2004204215;修回日期:2004205215作者简介:安 力(1973-),男,四川仪陇人,助理研究员,硕士研究生,核物理专业 中子与介质相互作用会产生散射反应等各种过程,当介质是大块物质时,作用过程更加复杂。

中子与含C 、H 、O 、N 等低原子序数元素材料的作用不同于与含高原子序数元素材料的作用,它主要表现为弹性散射,能谱软化也较快。

在以往的基准实验中,中子源均置于球中心测量泄露出来的中子能谱,而对球、柱组合装置的基准实验开展较少。

鉴于此,本工作拟测量含氢慢化球、含氢慢化球与<24cm ×30cm 聚乙烯圆柱组合2种实验装置内部的中子能谱。

1 实验装置含氢慢化球外径34cm ,内径4cm ,厚度15cm ,由上、下两半组成,界面间有一条通球心的28mm ×28mm 方槽。

聚乙烯圆柱直径24cm ,两段各长15cm ,可以组合。

当在球和圆柱组合装置内测量时,含氢慢化球的孔道旋转90°,球柱组合实验装置示于图1。

图1 球柱组合实验装置示意图Fig.1 Schematic diagram of combination of the shell and polyethylene cylinder2 测量技术2.1 α粒子监测中子绝对产额通过测量d 2T反应产生的伴随α粒子得到。

金硅面垒半导体探测器中心与D +束成178.9°夹角。

为降低测量时伴随α粒子谱的本底,经过几次测量比较,对建立的2套靶室经过几次测量比较得知,将氚靶装在细管内、金硅面垒半导体探测器装在粗管内的靶室,其伴随α粒子谱的本底只有单一的相同尺寸的细管靶室的1/3~1/2,且本底仅占α粒子总数的2%~3%。

因此,本实验采用粗细两节管子构成的靶室进行相应的α粒子监测。

2.2 中子谱仪为减小空腔对测量结果的影响,建立了小型 晶体闪烁探测器,探测器的外径为19mm 。

晶体分为有外壳和没有外壳2种,对于无外壳的裸晶体尺寸为<18mm ×20mm ;带有外壳及光反射层的其晶体尺寸为<12mm ×16mm ,加上外壳后其尺寸也为<18mm ×20mm 。

对它们的性能进行了测试比较,有外壳及光反射层的 晶体的光输出比无外壳裸晶体的光输出大1/3左右;它们的(n 、γ)分辨性能均达到使用要求。

无外壳的裸晶体的能量分辨明显好于有外壳及光反射层的,对于14MeV 中子和137Cs γ射线源,无外壳的裸晶体的能量分辨分别为4.8%和11%。

而有外壳及光反射层的晶体对14MeV 中子能谱测量不能满足能量分辨要求。

因此实验选用没有外壳的<18mm ×20mm 的 晶体。

3 实验方法在0.75~14MeV 的中子能量范围内,相应的质子能量的光产额相差近80倍。

在如此宽的能量范围内一次测量中子能谱十分困难,它受脉冲形状甄别效果和探测器的线性范围的限制,也受多道分析器可用的线性道数范围的限制。

因此,在该能量范围内分2段测量,低能段为0.75~5MeV ,高能段为3~16MeV ,3.25~4.5MeV 是两段的衔接部分。

(n ,γ)分辨采用脉冲形状甄别实现,分辨效果与线路参数调节有很大关系,除了得到清晰的(n ,γ)分辨图形之外,更主要的是要保证在该过程中中子或γ数目不要丢失。

为此,将(n ,γ)分辨线路的常份额甄别器输出(相当于积分数)分一路输到定标器测量计数,同时也把分辨线路的时2幅变换输出分一路到另一个定标器监视计数,在正确调节分辨线路参数情况下,这2个定标器的计数非常接近,而且(n ,γ)分辨图形也清晰。

在测量低能段时,高于5MeV 质子能量的脉冲幅度都已经饱和,而线性门对含有饱和脉冲比没有饱和脉冲允许通过的能力要小得多。

因此,在测量低能段时,线性门的计数率控制在约600s -1,而对高能段,可以允许其控制在约800s -1左右。

为了定量测出线性门的门控损失,把线性门的门控输入、输出和多道分析器的“SCA ”输出同时输到三路计数器进行比较。

在上述计数率控制情况下,门控损失一般均小于09原子能科学技术 第38卷2%。

22Naγ射线源作为探测器的能量刻度用源,其2个能量分别为1.275和0.511MeV,对应的康普顿电子能量为1.062和0.341MeV。

利用质子和电子光产额关系的光表,即可得到脉冲幅度与质子能量的对应关系。

对于测量的高能段,用d2T反应的14MeV中子作为刻度标准,效果更好些。

4 实验测量含氢慢化球的直径与D+束重合,中子源距球面30cm,对于球柱组合该距离保持不变。

对于含氢慢化球,分别在球心、内表面、距球心9cm处和外表面4个位置进行了1MeV 以上中子能谱测量。

测量球表面处的中子能谱时,因为缺少样品方槽相应的塞块,样品球水平旋转90°。

对于球柱组合装置(球柱紧密接合,柱轴线在D束零度方向上重合),测量了球、柱交界面和距离球心31.8cm处的1MeV以上的中子能谱。

5 实验结果实验得到的大于1MeV的中子能谱是一个微分的结果。

2个实验装置的测量能谱随位置的变化关系示于图2和3。

由图2和3可以看出,介质内不同位置的能谱变化,反映了中子在介质内输运的宏观效果,越是远离源中子,谱的软化越明显;球柱交界面处的中子谱比没有聚乙烯柱时要硬一些,同时通量也高些,这是由于聚乙烯柱对中子反射的缘故。

为了进一步比较,将能谱处理成不同阈能时的积分中子数。

积分中子数就是对于一个源中子入射时,高于1MeV的不同能量阈值的总中子数。

单一球体和球柱组合2种情况下的积分中子数随中子能量的变化关系示于图4和5。

从图4和2可看出,在球心和内球面,散射中子的份额明显大于距球心9cm的结果;在外球面,不同能量阈值的积分中子数的分布曲线较平坦,这是因为探测器只能探测到泄漏的中子的缘故。

对于球柱组合,从图5和3可看出,增加聚乙烯柱时,明显增加了中子的慢化效果,低能段尤为明显。

6 不确定度分析实验布置、实验测量和中子源条件等是实验结果不确定度的主要来源。

由于中子源在实验装置外部,周围环境及各种物质散射的中子有可能进入样品,与入射的中子一起被探测器记录,从而引起误差,其不确定度估计为3%。

测量的2段反冲质子谱用伴随α粒子归一,能谱和不同阈能的积分中子数均归一到一个源中子,因此中子绝对产额监测误差对实验19增刊 安 力等:含氢介质内中子能谱测量结果的影响较大,从2.1节对不同结构靶室的α粒子监测结果的比较可看出,中子绝对源强监测的不确定度小于2.5%。

实验时使用的是新氚靶,加速器的D+束流最大不超过10μA,从监测的伴随α粒子谱看出,d2D反应的中子占源中子份额不足1%。

由于测量的中子能量范围较大,低能部分利用22Naγ源刻度,高能部分利用源中子刻度。

源中子的能量依赖于D+能量,而D+能量和加速器的高压及氚靶的厚度相关,因此源中子的能量存在误差。

同时由于系统能量刻度存在大约2%的不确定度,所以能量刻度总的不确定度为2.5%。

综合以上考虑,本实验结果的不确定度为4.8%。

7 结论在加速器的源条件下,得到了d2T中子入射到含氢慢化球、含氢慢化球与<24cm×30cm聚乙烯圆柱组合2种实验装置内部的能谱。

在实验的不确定度范围内,本结果可供二维数值模拟计算方法检验参考使用。

相关文档
最新文档