基于CLYC闪烁体的中子能谱测量及反演方法

doi:10.13700/j.bh.1001-5965.2019.0643

北京航空航天大学学报 ›› 2021, Vol. 47 ›› Issue (1): 106-114.doi: 10.13700/j.bh.1001-5965.2019.0643

基于CLYC闪烁体的中子能谱测量及反演方法

侯东辉^1,2, 张珅毅^1,2, 杨祎罡³, 王琦标⁴, 张斌全^1,5, 余庆龙^1,5

1. 中国科学院国家空间科学中心, 北京 100190;
2. 中国科学院大学, 北京 100049;
3. 清华大学工程物理系, 北京 100084;
4. 成都理工大学核技术与自动化工程学院, 成都 610051;
5. 天基空间环境探测北京市重点实验室, 北京 100190

收稿日期:2019-12-24 发布日期:2021-01-29
通讯作者: 张珅毅 E-mail:zsy@nssc.ac.cn
作者简介:侯东辉,女,博士研究生。主要研究方向:探测器设计仿真及数据处理;张珅毅,男,博士,研究员,博士生导师。主要研究方向:粒子探测器设计。
基金资助:
国家自然科学基金（41941001）

Neutron measurement and inversion based on CLYC scintillator

HOU Donghui^1,2, ZHANG Shenyi^1,2, YANG Yigang³, WANG Qibiao⁴, ZHANG Binquan^1,5, YU Qinglong^1,5

1. National Space Science Center, Chinese Academy of Sciences, Beijing 100190, China;
2. University of Chinese Academy of Sciences, Beijing 100049, China;
3. Department of Engineering Physics, Tsinghua University, Beijing 100084, China;
4. College of Nuclear Technology and Automation Engineering, Chengdu University of Technology, Chengdu 610051, China;
5. Beijing Key Laboratory of Space Environment Exploration, Beijing 100190, China

Received:2019-12-24 Published:2021-01-29

摘要/Abstract

摘要： 空间中子是影响航天器和航天员安全的重要辐射要素之一。优化中子探测器，提高测量效率，提升反演精度是中子测量的难点。中国空间站将搭载一种基于新型中子探测材料Cs2LiYCl6：Ce（CLYC）闪烁体的中子探测器，该探测器具有同时测量热中子和快中子，以及探测效率高等特点。针对该新型探测器的中子能谱反演，分析了不同能量中子在该探测器中的响应特点，分析了中子反演常用的概率迭代法和非负最小二乘（NNLS）法的优缺点，考虑到这2种方法在CLYC探测器反演应用中的不足，提出了基于增广矩阵的非负最小二乘（AM-NNLS）法。数值实验结果表明：AM-NNLS法具有反演运算效率高和反演相对误差小的特点，验证了所提方法的有效性。

关键词: 空间站, 蒙特卡罗仿真, 中子探测器, 最小二乘近似, 能谱反演, Cs2LiYCl6∶Ce(CLYC)

Abstract: Neutron is one of the important radiation factors that affect the safety of spacecraft and astronauts. It is difficult for neutron measurement to optimize the neutron detector, improve the measurement efficiency, and improve the inversion accuracy. A neutron detector based on a new type of neutron detection material Cs2LiYCl6:Ce (CLYC) scintillator will be installed on the space station in China. This detector has the characteristics of measuring thermal neutrons and fast neutrons simultaneously, and has high detection efficiency, etc. For the neutron spectrum inversion of this detector, the response characteristics of different energy neutrons in the detetor are analyzed, and the advantages and disadvantages of the probabilistic iterative method and Non-Negative Least Square (NNLS) method commonly used in neutron spectrum inversion are analyzed, considering the disadvantages of these two methods in the inversion application of CLYC detector, a Non-Negative Least Square method based on the Augmented Matrix (AM-NNLS) is proposed. The numerical experiment results show that the AM-NNLS method has the characteristics of high inversion operation efficiency and small inversion relative error, which verifies the effectiveness of the method.

Key words: space station, Monte Carlo simulation, neutron detector, least squares approximations, spectrum inversion, Cs2LiYCl6:Ce (CLYC)

中图分类号:

V447⁺.1

侯东辉, 张珅毅, 杨祎罡, 王琦标, 张斌全, 余庆龙. 基于CLYC闪烁体的中子能谱测量及反演方法[J]. 北京航空航天大学学报, 2021, 47(1): 106-114.

HOU Donghui, ZHANG Shenyi, YANG Yigang, WANG Qibiao, ZHANG Binquan, YU Qinglong. Neutron measurement and inversion based on CLYC scintillator[J]. Journal of Beijing University of Aeronautics and Astronautics, 2021, 47(1): 106-114.

参考文献

[1] MOUNTFORD P J,TEMPERTON D H.Recommendations of the international commission on radiological protection(ICRP)1990[J].European Journal of Nuclear Medicine,1992,19(2):77-79.
[2] MAZUR J E,CRAIN W R,LOOPER M D,et al.New measurements of total ionizing dose in the lunar environment[J].Space Weather-The International Journal of Research & Applications,2011,9(7):1-12.
[3] KOSHIISHI H,MATSUMOTO H,CHISHIKI A,et al.Evaluation of the neutron radiation environment inside the International Space Station based on the bonner ball neutron detector experiment[J].Radiation Measurements,2007,42(9):1510-1520.
[4] HASSLER D M,ZEITLIN C,WIMMER-SCHWEINGRUBER R F,et al.The radiation assessment detector (RAD) investigation[J].Space Science Reviews,2012,170(1-4):503-558.
[5] TANIGUCHI T,UEDA N,NAKAZAWA M,et al.Systematic study on spectral effects in the adjustment calculations using the NEUPAC-83 code[M].Berlin:Springer,1985:685-691.
[6] KOHLER J,EHRESMANN B,MARTIN C,et al.Inversion of neutron/gamma spectra from scintillator measurements[J].Nuclear Instruments and Methods in Physics Research B,2011,269:2641-2648.
[7] 王冬,何彬,张全虎.用遗传算法求解中子能谱[J].原子能科学技术,2010,44(10):1270-1275. WANG D,HE B,ZHANG Q H.Unfolding neutron spectrum using genertic algorithm[J].Atomic Energy Science and Technology,2010,44(10):1270-1275(in Chinese).
[8] 杨鑫,李润东,刘汉刚,等.基于概率迭代的NDP反演方法[J].计算物理,2012,29(6):891-900. YANG X,LI R D,LIU H G,et al.An unfolding method of NDP based on probability iteration[J].Chinese Journal of Computational Physics,2012,29(6):891-900(in Chinese).
[9] D'OLYMPIA N W,NATHAN W.Development of novel neutron and gamma-ray scintillators:Cs2LiYCl6:Ce and CeBr₃[D].Boston:University of Massachusetts,2013:8-40.
[10] 胡圣荣,戴纳新.病态线性方程组的新解法:增广方程组法[J].华南农业大学学报,2009,30(1):119-121. HU S R,DAI N X.A novel method for solving III-conditioned liner system:Augmented system method[J].Journal of South China Agricultural University,2009,30(1):119-121(in Chinese).
[11] 李肖.反符合探测杯在空间粒子探测中的应用[D].北京:中国科学院,2015:5-23. LI X.The application of anti-coincidence detective cup technology in space particle detection[D].Beijing:University of Chinese Academy of Science,2015:5-23(in Chinese).
[12] WHITNEY C M,SOUNDARA-PANDIAN L,JOHNSON E B,et al.Gamma-neutron imaging system utilizing pulse shape discrimination with CLYC[J].Nuclear Instruments and Methods in Physics Research Section A:Accelerators,Spectrometers,Detectors and Associated Equipment,2015,784:346-351.
[13] LEE D W,STONEHILL L C,KLIMENKO A,et al.Pulse-shape analysis of Cs2LiYCl6:Ce scintillator for neutron and gamma-ray discrimination[J].Nuclear Instruments and Methods in Physics Research Section A:Accelerators,Spectrometers,Detectors and Associated Equipment,2012,664(1):1-5.
[14] WANG Q B,TUO X G,DENG C,et al.Characterization of a Cs2LiYCl6:Ce³⁺ scintillator coupled with two silicon photomultiplier arrays of different sizes[J].Nuclear Instruments and Methods in Physics Research Section A:Accelerators,Spectrometers,Detectors and Associated Equipment,2019,942:1-5.
[15] D'OLYMPIA N W,CHOWDHURY P,LISTER C J,et al.Pulse-shape analysis of CLYC for thermal neutrons,fast neutrons,and gamma-rays[J].Nuclear Instruments and Methods in Physics Research Section A:Accereators,Spectrometers,Detectors and Associated Equipment,2013,714:121-127.
[16] BROWN D,CHADWICK M B,CAPOTE R,et al.ENDF/B-VIII.0:The 8th major release of the nuclear reaction data library with CIELO-project cross sections,new stantards and thermal scattering data[J].Nuclear Data Sheets,2018,148(2):1-142.
[17] AGOSTINELLI S,ALLISON J,AMAKO K,et al.GEANT4-A simulation toolkit[J].Nuclear Instruments and Methods in Physics Research Section A:Accelerators,Spectrometers,Detectors and Associated Equipment,2003,506(3):250-303.
[18] ZHANG S Y,ZHANG X G,WANG C Q,et al.The geometric factor of high energy protons detector on FY-3 satellite[J].Science China Earth Sciences,2014,57(10):2558-2566.
[19] KHARYTONOV A,BÖHM,E,WIMMER-SCHWEINGRUBER R F.Regularization methods used in error analysis of solar particle spectra measured on SOHO/EPHIN[J].Astronomy & Astrophysics,2009,495(2):663-675.
[20] GUO J N,BANJAC S,ROSTEL L,et al.Implementation and validation of the GEANT4/AtRIS code to model the radiation environment at Mars[J].Journal of Space Weather and Space Climate,2019,9(3):1-30.
[21] NOWAK R D,KOLACZYK E D.A statistical multiscale framework for Poisson inverse problems[J].IEEE Transactions on Information Theory,2000,46(5):1811-1825.
[22] WILLIAM H P,BRIAN P F,SAUL A,et al.Numerical recipes:The art of scientific computing[M].Cambridge:Cambridge University Press,2007:809-816.
[23] LAWSON C L,HANSON R J.Solving least squares problems[M].Philadelphia:Society for lndustrial and Applied Mathematics,1987:1-5.
[24] BÖHM E,KHARYTONOV A,WIMMER-SCHWEINGRUBER R F.Solar energetic particle spectra from the SOHO-EPHIN sensor by application of regularization methods[J].Astronomy & Astrophysics,2007,473(2):673-682.
[25] MORALE J L,NOCEDAL J.Remark on ‘algorithm 778:L-BFGS-B:Fortran subroutines for large-scale bound constrained optimization’[J].ACM Transactions on Mathematical Software,2011,38(1):1-4.

基于CLYC闪烁体的中子能谱测量及反演方法

Neutron measurement and inversion based on CLYC scintillator

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 8

编辑推荐

Metrics

本文评价

[1]	葛祥雨, 黄杰, 周前祥, 柳忠起. 空间站在轨维修操作复杂度评价及试验验证[J]. 北京航空航天大学学报, 2019, 45(11): 2228-2236.
[2]	程阳, 程林, 张庆振, 倪昆, 储培. 基于在线约束限制的飞行器预测校正制导[J]. 北京航空航天大学学报, 2017, 43(10): 2143-2153.
[3]	党庆庆, 金磊, 徐世杰. 基于极点配置的空间站角动量管理[J]. 北京航空航天大学学报, 2016, 42(8): 1715-1723.
[4]	杨东升, 刘猛, 庞丽萍, 余青霓, 黄勇. H₂O对空间站5A分子筛CO₂去除性能影响[J]. 北京航空航天大学学报, 2015, 41(8): 1485-1491.
[5]	王伟斌, 路辉, 郎荣玲. 基于数据拟合的雷达高度表告警信号测试方法[J]. 北京航空航天大学学报, 2010, 36(2): 180-183.
[6]	龙军, 康锐, 马麟, 程海龙. 任意寿命分布的多部件系统备件配置优化算法[J]. 北京航空航天大学学报, 2007, 33(06): 698-700.
[7]	付仕明, 徐小平, 李劲东, 潘增富. 空间站乘员睡眠区二氧化碳聚集现象[J]. 北京航空航天大学学报, 2007, 33(05): 523-526.
[8]	邢玉明, 崔海亭, 袁修干. 高温熔盐相变蓄热系统的数值模拟[J]. 北京航空航天大学学报, 2002, 28(3): 295-297.