用于γ全吸收装置的大体积BaF2探测单元的α/γ鉴别方法

上一篇

下一篇

邹翀, 张奇玮, 栾广源, 吴鸿毅, 罗淏天, 陈玄博, 王晓宇, 贺国珠, 任杰, 黄翰雄, 阮锡超, 鲍杰, 朱兴华. 用于γ全吸收装置的大体积BaF2探测单元的α/γ鉴别方法[J]. 物理学报, 2025, 74(10): 102801-1. doi: 10.7498/aps.74.20250017
引用本文: 邹翀, 张奇玮, 栾广源, 吴鸿毅, 罗淏天, 陈玄博, 王晓宇, 贺国珠, 任杰, 黄翰雄, 阮锡超, 鲍杰, 朱兴华. 用于γ全吸收装置的大体积BaF2探测单元的α/γ鉴别方法[J]. 物理学报, 2025, 74(10): 102801-1. doi: 10.7498/aps.74.20250017
Chong ZOU, Qiwei ZHANG, Guangyuan LUAN, Hongyi WU, Haotian LUO, Xuanbo CHEN, Xiaoyu WANG, Guozhu HE, Jie REN, Hanxiong HUANG, Xichao RUAN, Jie BAO, Xinghua ZHU. α/γ discrimination method for bulky BaF2 detector used in γ total absorption facility[J]. Acta Physica Sinica, 2025, 74(10): 102801-1. doi: 10.7498/aps.74.20250017
Citation: Chong ZOU, Qiwei ZHANG, Guangyuan LUAN, Hongyi WU, Haotian LUO, Xuanbo CHEN, Xiaoyu WANG, Guozhu HE, Jie REN, Hanxiong HUANG, Xichao RUAN, Jie BAO, Xinghua ZHU. α/γ discrimination method for bulky BaF2 detector used in γ total absorption facility[J]. Acta Physica Sinica, 2025, 74(10): 102801-1. doi: 10.7498/aps.74.20250017

用于γ全吸收装置的大体积BaF2探测单元的α/γ鉴别方法

    作者简介: 邹翀 E-mail: zouc@mail.chinaships.com .
    通讯作者: E-mail: zqwvictor@126.com.;  E-mail: lgyciae@hotmail.com.; 
  • 中图分类号: 28.20.Np, 07.85.Nc, 01.52.+r

α/γ discrimination method for bulky BaF2 detector used in γ total absorption facility

    Corresponding authors: E-mail: zqwvictor@126.com.;  E-mail: lgyciae@hotmail.com.; 
  • MSC: 28.20.Np, 07.85.Nc, 01.52.+r

  • 摘要: γ全吸收型探测装置由40个BaF2探测单元组成, 用于在线测量中子辐射俘获反应截面数据, 填补国内实验数据的空白. 实验本底的一个重要来源是BaF2晶体自身包含的α粒子. 为扣除α粒子本底, 数据获取系统采用全波形采集的方式, 在线实验存储的数据量为118 MB/s, 产生了较长的死时间. 利用BaF2探测单元测量3种放射源(22Na, 137Cs及60Co)的实验数据, 确定了信号波形的积分长度为2000 ns时, 能够达到最佳能量分辨率; 使用快总成分比、脉冲宽度和时间衰减常数3种方法进行α粒子和γ射线的鉴别, 计算得到快总成分比(快成分5 ns总成分200 ns)方法的品质因子为1.19—1.41, 脉冲宽度(10%峰值)方法的品质因子为0.94—1.04, 时间衰减常数方法的品质因子为0.93—1.07. 通过品质因子的定量分析和能谱的比较, 确定快总成分比鉴别方法效果最好, 能够有效去除α粒子本底, 为下一步升级数据获取系统, 减少实验数据量, 降低截面数据不确定度奠定基础.
  • 加载中
  • 图 1  (a) BaF2探测单元; (b) GTAF

    Figure 1.  (a) BaF2 detector unit; (b) GTAF.

    图 2  BaF2探测单元的信号波形

    Figure 2.  Signal waveform of the BaF2 detection unit.

    图 3  BaF2探测单元的能谱

    Figure 3.  Energy spectrum of the BaF2 detection unit.

    图 4  不同积分长度能量分辨率的比较

    Figure 4.  Comparison of energy resolution of different integral length.

    图 5  品质因子计算分布图

    Figure 5.  Distribution of FOM calculation.

    图 6  波形的快成分与总成分

    Figure 6.  Fast and total components of waveform.

    图 7  不同快总成分比的品质因子 (a) 22Na; (b) 137Cs; (c) 60Co

    Figure 7.  FOM with different ratios of fast to total component: (a) 22Na; (b) 137Cs; (c) 60Co.

    图 8  能量与快总成分比的二维谱 (a) 22Na; (b) 137Cs; (c) 60Co

    Figure 8.  Energy versus ratio of fast to total component spectrum: (a) 22Na; (b) 137Cs; (c) 60Co.

    图 9  波形的脉冲宽度

    Figure 9.  Pulse width of waveform.

    图 10  不同波形脉冲宽度的品质因子

    Figure 10.  FOM with different pulse width of waveforms.

    图 11  能量与脉冲宽度的二维谱 (a) 22Na; (b) 137Cs; (c) 60Co

    Figure 11.  Energy versus pulse width spectrum: (a) 22Na; (b) 137Cs; (c) 60Co.

    图 12  波形的时间衰减常数拟合

    Figure 12.  Fitting of the time decay constant of waveform.

    图 13  能量与时间衰减常数的二维谱 (a) 22Na; (b) 137Cs; (c) 60Co

    Figure 13.  Energy versus time decay constant spectrum: (a) 22Na; (b) 137Cs; (c) 60Co.

    图 14  不同粒子鉴别方法处理后得能谱比较 (a) 22Na; (b) 137Cs; (c) 60Co

    Figure 14.  Comparison of energy spectra obtained by different particle identification methods: (a) 22Na; (b) 137Cs; (c) 60Co.

    表 1  不同粒子鉴别方法的品质因子比较

    Table 1.  Comparison of FOM of different particle identification methods.

    放射源 快总成分比(快成分
    5 ns/总成分200 ns)
    脉冲宽度
    (10%峰值)/ns
    时间衰减
    常数/s–1
    22Na 1.19 1.04 1.07
    137Cs 1.30 0.94 0.93
    60Co 1.41 1.01 0.96
    下载: 导出CSV
  • [1] 丁大钊 2001 中国基础科学 1 13 doi: 10.3969/j.issn.1009-2412.2001.01.003 Ding D Z 2001 China Basic Sci. 1 13 doi: 10.3969/j.issn.1009-2412.2001.01.003
    [2] Chandler D, Hartanto D, Bae J W, Burg K M, Robert Y, Sizemore C 2025 Ann. Nucl. Energy 211 110920 doi: 10.1016/j.anucene.2024.110920
    [3] Arnould M, Katsuma M 2008 International Conference on Nuclear Data for Science and Technology Nice, France, April 22–27, 2007 p7
    [4] Kompe D 1969 Nucl. Phys. 133 513 doi: 10.1016/0375-9474(69)90552-1
    [5] Wisshak K, Kappeler F, Reffo G 1984 Nucl. Sci. Eng. 88 594 doi: 10.13182/NSE84-A18376
    [6] Terada K, Katabuchi T, Mizumoto M, Arai T, Saito T, Igashira M, Hirose K, Nakamura S, Kimura A, Harada H, Hori J, Kino K, Kiyanagi Y 2015 Prog. Nucl. Energy 82 118 doi: 10.1016/j.pnucene.2014.07.017
    [7] Kobayashi K, Lee S, Yamamoto S 2004 Nucl. Sci. Eng. 146 209 doi: 10.13182/NSE04-A2404
    [8] Guber K H, Derrien H, Leal L C, Arbanas G, Wiarda D, Koehler P E, Harvey A 2010 Phys. Rev. C 82 057601 doi: 10.1103/PhysRevC.82.057601
    [9] Lee J, Hori J I, Nakajima K, Sano T, Lee S 2017 J. Nucl. Sci. Tech. 54 1046 doi: 10.1080/00223131.2017.1344575
    [10] Kim H I, Paradela C, Sirakov I, Becker B, Capote R, Gunsing F, KimG N, Kopecky S, Lampoudis C, LeeY O, Massarczyk R, Moens A, Moxon M, PronyaevV G, Schillebeeckx P, Wynants R 2016 Eur. Phys. J. A 52 170 doi: 10.1140/epja/i2016-16170-6
    [11] Mingrone F, Massimi C, Altstadt S, et al. 2014 Nuclear Data Sheets 119 18
    [12] Ren J, Ruan X C, Bao J, et al. 2019 Radiat. Detect. Tech. Methods 3 52 doi: 10.1007/s41605-019-0129-8
    [13] Wisshak K, Voss F, Kaeppeler F, Krticka M, Gallino R 2006 Phys. Rev. C 73 015802 doi: 10.1103/PhysRevC.73.015802
    [14] Endo S, Kimura A, Nakamura S, Iwamoto O, Iwamoto N, Rovira G, Terada K, Meigo S, Toh Y, Segawa M, Maeda M, Tsuneyama M 2022 J. Nucl. Sci. Techn. 59 318 doi: 10.1080/00223131.2021.1970040
    [15] Mosby S, Bredeweg T A, Couture A, Jandel M, Kawano T, Ullmann J, Henderson R A, Wu C Y 2018 Nucl. Data Sheets 148 312 doi: 10.1016/j.nds.2018.02.007
    [16] Laminack A, Blackmon J C, Couture A, Greene J P, Krticka M, Macon K T, Mosby S, Prokop C, Ullmann J L, Valenta S 2022 Phys. Rev. C 106 025802 doi: 10.1103/PhysRevC.106.025802
    [17] Zhong Q P, Zhou Z Y, Tang H Q, Chen X L, He G Z, Zhang Q W, Guo W X, Yuan J L, Ma X Y, Wang Q, Ruan X C, Li Z H 2008 Chin. Phys. C 32 102
    [18] 任杰, 阮锡超, 唐洪庆, 葛智刚, 黄翰雄, 敬罕涛, 唐靖宇, 黄蔚玲 2014 核技术 37 110521 doi: 10.11889/j.0253-3219.2014.hjs.37.100521 Ren J, Ruan X C, Tang H Q, Ge Z G, Huang H X, Jing H T, Tang J Y, Huang W L 2014 Nucl. Tech. 37 110521 doi: 10.11889/j.0253-3219.2014.hjs.37.100521
    [19] 张奇玮, 贺国珠, 黄兴, 阮锡超, 李志宏, 朱兴华 2014 原子能科学技术 48 612 doi: 10.7538/yzk.2014.48.S0.0612 Zhang Q W, He G Z, Huang X, Ruan X C, Li Z H, Zhu X H 2014 At. Energy Sci. Tech. 48 612 doi: 10.7538/yzk.2014.48.S0.0612
    [20] 张奇玮, 贺国珠, 黄兴, 程品晶, 阮锡超, 朱兴华 2016 原子能科学技术 50 536 doi: 10.7538/yzk.2016.50.03.0536 Zhang Q W, He G Z, Huang X, Cheng P J, Ruan X C, Zhu X H 2016 At. Energy Sci. Tech. 50 536 doi: 10.7538/yzk.2016.50.03.0536
    [21] 张奇玮, 栾广源, 任杰, 等 2021 物理学报 70 222801 doi: 10.7498/aps.70.20210742 Zhang Q W, Luan G Y, Ren J, et al. 2021 Acta Phys. Sin. 70 222801 doi: 10.7498/aps.70.20210742
    [22] 栾广源, 任杰, 张奇玮, 阮锡超, 贺国珠, 程品晶, 郭明伟 2022 同位素 35 273 doi: 10.7538/tws.2021.youxian.064 Luan G Y, Ren J, Zhang Q W, Ruan X C, He G Z, Cheng P J, Guo M W 2022 J. Isotopes 35 273 doi: 10.7538/tws.2021.youxian.064
    [23] 马霄云, 仲启平, 周祖英, 郭维新, 袁继龙, 王强, 张奇玮, 苏明, 阮锡超, 鲍杰, 黄翰雄, 蒋静, 聂阳波, 李霞, 刘刚, 兰长林 2009 原子能科学技术 43 180 doi: 10.7538/yzk.2009.43.02.0180 Ma X Y, Zhong Q P, Zhou Z Y, Guo W X, Yuan J L, Wang Q, Zhang Q W, Su M, Ruan X C, Bao J, Huang H X, Jiang J, Nie Y B, Li X, Liu G, Lan C L 2009 At. Energy Sci. Tech. 43 180 doi: 10.7538/yzk.2009.43.02.0180
    [24] 张奇玮, 栾广源, 贺国珠, 程品晶, 阮锡超, 朱兴华 2020 原子核物理评论 37 771 doi: 10.11804/NuclPhysRev.37.2019CNPC30 Zhang Q W, Luan G Y, He G Z, Cheng P J, Ruan X C, Zhu X H 2020 Nucl. Phys. Rev. 37 771 doi: 10.11804/NuclPhysRev.37.2019CNPC30
    [25] 张奇玮, 栾广源, 郭明伟, 贺国珠, 阮锡超, 邹翀, 朱兴华 2021 现代应用物理 12 040401 doi: 10.12061/j.issn.2095-6223.2021.040401 Zhang Q W, Luan G Y, Guo M W, He G Z, Ruan X C, Zou C, Zhu X H 2021 Mod. Appl. Phys. 12 040401 doi: 10.12061/j.issn.2095-6223.2021.040401
    [26] 彭猛, 贺国珠, 骆宏, 张奇玮, 石斌, 周祖英, 兰长林 2016 原子能科学技术 50 1866 doi: 10.7538/yzk.2016.50.10.1866 Peng M, He G Z, Luo H, Zhang Q W, Shi B, Zhou Z Y, Lan C L 2016 At. Energy Sci. Tech. 50 1866 doi: 10.7538/yzk.2016.50.10.1866
    [27] 石斌, 彭猛, 张奇玮, 贺国珠, 周祖英, 唐洪庆 2018 原子能科学技术 52 1537 doi: 10.7538/yzk.2017.youxian.0817 Shi B, Peng M, Zhang Q W, He G Z, Zhou Z Y, Tang H Q 2018 At. Energy Sci. Tech. 52 1537 doi: 10.7538/yzk.2017.youxian.0817
    [28] 张奇玮, 贺国珠, 栾广源, 程品晶, 阮锡超, 朱兴华 2021 强激光与粒子束 33 0440 doi: 10.11884/HPLPB202133.200220 Zhang Q W, He G Z, Luan G Y, Cheng P J, Ruan X C, Zhu X H 2021 Power Laser Part. Beams 33 0440 doi: 10.11884/HPLPB202133.200220
    [29] 王晓宇, 贺国珠, 张奇玮, 罗淏天, 徐阔之, 李倩, 解立坤, 栾广源, 阮锡超, 邹翀, 陈玄博, 吴鸿毅, 樊瑞睿, 蒋伟, 曹平, 余滔 2024 原子能科学技术 58 2262 doi: 10.7538/yzk.2023.youxian.0887 Wang X Y, He G Z, Zhang Q W, Luo H T, Xu K Z, Li Q, Xie L K, Luan G Y, Ruan X C, Zou C, Chen X B, Wu H Y, Fan R R, Jiang W, Cao P, Yu T 2024 At. Energy Sci. Tech. 58 2262 doi: 10.7538/yzk.2023.youxian.0887
    [30] 罗淏天, 张奇玮, 栾广源, 等 2024 物理学报 73 072801 doi: 10.7498/aps.73.20231957 Luo H T, Zhang Q W, Luan G Y, et al. 2024 Acta Phys. Sin. 73 072801 doi: 10.7498/aps.73.20231957
  • 加载中
图( 14) 表( 1)
计量
  • 文章访问数:  34
  • HTML全文浏览数:  34
  • PDF下载数:  2
  • 施引文献:  0
出版历程
  • 收稿日期:  2025-01-05
  • 刊出日期:  2025-05-20

用于γ全吸收装置的大体积BaF2探测单元的α/γ鉴别方法

    通讯作者: E-mail: zqwvictor@126.com.; 
    通讯作者: E-mail: lgyciae@hotmail.com.; 
    作者简介: 邹翀 E-mail: zouc@mail.chinaships.com
  • 1. 中国原子能科学研究院, 核数据重点实验室, 北京 102413
  • 2. 中国船舶工业贸易有限公司, 北京 100044
  • 3. 华能山东石岛湾核电有限公司, 荣成 264312

摘要: γ全吸收型探测装置由40个BaF2探测单元组成, 用于在线测量中子辐射俘获反应截面数据, 填补国内实验数据的空白. 实验本底的一个重要来源是BaF2晶体自身包含的α粒子. 为扣除α粒子本底, 数据获取系统采用全波形采集的方式, 在线实验存储的数据量为118 MB/s, 产生了较长的死时间. 利用BaF2探测单元测量3种放射源(22Na, 137Cs及60Co)的实验数据, 确定了信号波形的积分长度为2000 ns时, 能够达到最佳能量分辨率; 使用快总成分比、脉冲宽度和时间衰减常数3种方法进行α粒子和γ射线的鉴别, 计算得到快总成分比(快成分5 ns总成分200 ns)方法的品质因子为1.19—1.41, 脉冲宽度(10%峰值)方法的品质因子为0.94—1.04, 时间衰减常数方法的品质因子为0.93—1.07. 通过品质因子的定量分析和能谱的比较, 确定快总成分比鉴别方法效果最好, 能够有效去除α粒子本底, 为下一步升级数据获取系统, 减少实验数据量, 降低截面数据不确定度奠定基础.

English Abstract

    • 中子辐射俘获反应即(n, γ)反应在热中子能区到keV能区是主要的中子消失反应, 其截面数据广泛的应用于先进核能研究[1]、反应堆设计[2]以及核天体物理研究[3]中. 国际上早期采用缓发γ法(活化法)测量(n, γ)反应截面数据, 但只适用于单能中子的实验条件, 也难以区分通过不同反应道得到相同产物核的数据. 20世纪60年代开始, 基于飞行时间技术(TOF)的瞬发γ法(在线法)被用 在大体积液闪[4]、Moxon-Rae[5]、NaI (Tl) [6]、BGO[7]等不同类型的探测器上测量(n, γ)反应截面数据, 但由于这些探测器自身的限制, 测得截面数据的不确定度较大. 20世纪90年代之后国际上几个重要的实验室如: 美国ORNL的ORELA[8]、欧洲JRC的GELINA[9], CERN的n_TOF[10]、日本KURRI-LINAC[11]以及中国散裂中子源(CSNS)的反角白光中子源(Back-n)[12]又相继建立了C6D6探测器, 用于测量稳定核素的(n, γ)反应截面数据. 20世纪80年代, 发现BaF2晶体具有探测效率高、时间分辨好的特点, 开始制作成闪烁体探测器用于γ射线的测量, 国际上先后建造了3套γ全吸收型BaF2探测装置, 分别是德国FZK Karlsruhe的4πBaF2 [13]、欧洲(瑞士) CERN 的TAC[14]和美国LANL的DANCE [15,16], 用于测量裂变核素、放射性核素以及样品量低截面较小核素的中子俘获反应截面. 目前A < 120和A > 180范围内, 裂变核素(n, γ)反应的高精度截面数据测量已经成为国际上核数据研究的热点.

      国内由于缺乏高质量的探测装置和测量方法, 再加上发达国家对关键核截面数据实行壁垒保护, 使得一些重要的裂变核素如: 235U, 239Pu, 241Am等缺少截面数据或者数据存在较大分歧, 无法满足工程建设和科学研究的需求. 为了填补这一空白, 中国原子能科学研究院核数据重点实验室建成了我国首套γ全吸收型探测装置(gamma total absorption facility, GTAF), 基于中国散裂中子源的反角白光中子源Back-n已经开展了一些稳定核素的(n, γ)反应截面在线测量实验, 得到了较好的实验结果[1719].

      GTAF由40个BaF2探测单元组成, 通过分析在线实验的数据, 发现BaF2探测单元自身包含的α粒子是本底的重要来源, 扣除α粒子本底以提高效应本底比, 是降低截面数据不确定度的关键[2022]. 目前GTAF的数据获取系统使用全波形采集的方式, 所有探测单元满足触发条件的信号波形都会被记录, 这种方法的优点是在离线数据处理的过程中, 可以通过波形分析方法扣除α粒子本底, 得到信号的能量信息和时间信息; 缺点是需要采集、传输、存储的数据量庞大. 在中国散裂中子源140 kW功率的束流条件下, GTAF在线实验的平均数据量为118 MB/s, 峰值数据量能达到160 MB/s, 对于电子学的触发判选和数据采集都有极高的要求, 数据获取系统的死时间达到了20%以上, 直接影响了截面数据的不确定度, 后续随着散裂中子源束流功率的提升, 在线实验的数据量会继续增加, 因此计划对数据获取系统进行升级, 目标是降低获取系统需要记录的数据量.

    • BaF2晶体是一种无机闪烁体[2326], 具有两种发光成分: 其中快成分的光衰减时间为0.6 ns, 波长峰值为225 nm, 这决定了BaF2晶体具有非常好的时间分辨; 慢成分的光衰减时间为620 ns, 波长峰值为325 nm; BaF2晶体的中子灵敏度低, 折射率接近玻璃, 使闪烁光更容易进入光电倍增管; 其溶解度小, 不潮解, 便于加工和使用. 这些优点使其适合用于基于飞行时间方法的(n, γ)反应截面在线测量实验.

      42个五棱台和六棱台可以组成一个内半径10 cm, 厚度15 cm的闭合球壳. BaF2晶体(端面直径14 cm, 高度 15 cm)与5 in光电倍增管耦合后, 封装成独立的BaF2探测单元(图1(a)). 为了方便中子束流的进出, GTAF束线进出的两个位置空缺. 其余40个位置由BaF2探测单元填充, 组成了GTAF(图1(b)). 在线实验时, 中子束流与GTAF中心的样品发生(n, γ)反应, 生成复合核. 处于激发态的复合核, 通过出射级联γ射线的方式退激. 由于GTAF几乎覆盖了4π立体角, 退激的各条级联γ射线会被BaF2探测单元探测到, 记录为(n, γ)反应事件. 同时, 记录源中子产生的时刻和级联γ射线到达探测器的时刻, 以计算发生 (n, γ)反应的中子飞行时间, 确定中子能量[2730].

      由于Ba与Ra是同族元素, 因此BaF2晶体的天然放射性本底主要是由226Ra衰变链上的4种α粒子组成. 分别为226Ra到222Rn衰变时发射的能量为4.7 MeV 的α粒子、从222Rn 到218Po衰变时发射的能量为5.5 MeV的α粒子、从218Po到214Pb衰变时产生能量为6.0 MeV的α粒子及从214Po到210Pb衰变时产生能量为7.7 MeV的α粒子.

    • 本文选取了GTAF的一个BaF2探测单元, 采用CAEN公司的数字化仪DT5751(10 bit, 1 GS/s)获取探测器的信号波形, 测量3种放射源(22Na, 137Cs及60Co)的γ射线. 如图2所示, BaF2探测单元主要有3类信号: 1) 电子学噪声, 有快成分, 几乎没有慢成分, 来源于光电倍增管的暗电流和其他电子学电路; 2) γ射线, 在线实验的有效信号, 有明显的快成分和慢成分; 3) BaF2晶体自身的α粒子, 几乎没有快成分, 只有慢成分.

      数据获取系统设置的参数为: 触发位置100 ns, 阈值20 mV, 基线0—50 ns. 通过积分信号波形的脉冲幅度, 得到能谱. 如图3所示, 可以明显看到22Na (0.511 MeV), 137Cs (0.662 MeV), 60Co (1.17/1.33 MeV)的全能峰以及4种能量的α粒子. 在扣除本底后, 拟合计算得到能量分辨率, 其中60Co的两个全能峰无法完全区分, 按照平均值1.25 MeV计算.

      图4为不同积分长度下(以触发位置为积分起点), BaF2探测单元测量3种放射源得到的能量分辨率的比较. 积分长度从200 ns开始逐渐增大, 探测单元的能量分辨率逐渐变好, 积分长度在1900—2000 ns之间时, 能量分辨率达到最佳, 之后积分长度继续增大, 能量分辨率又随之逐渐变差. 由于BaF2晶体慢成分的光衰减时间为620 ns, 波形脉冲在1900—2000 ns时基本回到基线, 此时刚好采集到γ射线信号的完整波形, 能量分辨率最好, 因此在线实验选择能谱的积分长度为2000 ns.

    • 现阶段的生产工艺中, 很难使用化学的方法去除BaF2晶体原料中的Ra杂质, 因此需要采用粒子鉴别的方法, 尽可能地排除BaF2探测单元中的α粒子对γ射线信号的干扰, 有效降低在线实验的本底.

      图5所示, 使用粒子鉴别方法可以得到α粒子和γ射线的波形特征值分布图, 根据统计学原理, 这两种类型粒子的波形特征值都满足高斯分布. 国际上通常采用鉴别品质因子FOM作为评估探测器粒子鉴别能力的标准, FOM越大, 表明鉴别能力越强, 鉴别效果越理想.

      式中, μαμγ是分别使用高斯拟合后得到的α峰和γ峰的均值; FWHMα和FWHMγ分别为α峰和γ峰的半高全宽. 本文使用的3种粒子鉴别方法, 都通过FOM值的计算评估鉴别效果.

    • 图6所示, BaF2探测单元探测到的γ射线信号波形有明显的快成分, 之后信号幅度迅速下降, 与慢成分重叠在一起; α粒子的信号波形则没有快成分, 只有慢成分, 信号幅度下降缓慢. 快总成分比鉴别法从信号触发位置开始, 分别设定不同积分长度得到每个信号的快成分与总成分, 计算两者的比值作为波形特征值得到FOM.

      本文使用BaF2探测单元分别测量了3种放射源, 计算得到波形信号快总成分比值的FOM. 图7的横坐标是总成分的积分长度(50—2000 ns), F5—F30分别表示快成分的积分长度为5—30 ns. 可以明显看到快成分积分越长, FOM越小; 总成分的积分长度在200 ns时, FOM最大, 因此选择快成分5 ns与总成分200 ns的比值用于αγ鉴别.

      在确定了快总成分的参数后, 处理得到图8展示的二维谱. 红框圈出的部分是α粒子(能量1—4 MeV, 快总成分比0.021—0.128), 能够明显与γ射线区分出来, 在线实验数据处理时, 会把此范围内的波形信号去除掉.

    • 通过分析γ射线和α粒子在BaF2探测单元中的信号波形发现, 两者的波形脉冲宽度有明显的不同, 可利用这一差异进行粒子鉴别. 图9是BaF2探测单元测得的一个波形信号, 首先通过波形分析得到脉冲的峰值, 再使用插值计算得到脉冲上升过程中到达10%峰值位置的时间T_start10%, 同样使用插值计算得到脉冲下降过程中到达10%峰值的时间T_end10%.

      式中, w1—w5分别代表信号波形10%—50%峰值的脉冲宽度.

      使用不同的波形脉冲宽度作为特征值, 得到3种放射源的FOM值, 如图10所示, 其中10%峰值的脉冲宽度粒子鉴别能力最好. 由于α粒子没有慢成分, 与γ射线相比较, 峰值比例越低, 对应位置的脉冲宽度变长得越明显, 也就越容易把两者区分出来. 图11是10%峰值脉冲宽度的二维谱, 红框圈出的部分是α粒子(能量1—4 MeV, 10%峰值脉冲宽度120—880 ns), 需在数据处理中去除掉.

    • BaF2探测单元闪烁光的特点是信号波形到达峰值位置后, 按照负e指数规律衰减, 波形幅度可以使用(3)式进行拟合:

      式中, V(t)是指t时刻的波形幅度, V0是指信号基线的幅度, A是指波形峰值的幅度, t0是指到达峰值的时刻, τ是指信号波形的时间衰减常数.

      γ射线的信号波形包含快成分与慢成分, 波形在快速上升到峰值位置后, 又会快速下降到慢成分的位置, 之后就缓慢下降到基线位置, 因此γ射线信号波形的时间衰减常数τ主要是由快成分决定的. α粒子的信号波形几乎没有快成分只有慢成分, 它的时间衰减常数τ主要是由慢成分决定的. BaF2探测单元的快成分光衰减时间为0.6 ns, 慢成分光衰减时间为620 ns, 两者差异明显, 利用这一差异就可以实现α粒子和γ射线的鉴别. 如图12所示, 通过波形分析得到脉冲的峰值, 选择峰值位置之后2000 ns的波形, 利用(3)式拟合得到时间衰减常数.

      图13展示了时间衰减常数的二维谱, α粒子是红框圈出的部分(能量1—4 MeV, 时间衰减常数是55—45 s–1), 可以看到在1—1.3 MeV范围内, γ射线与α粒子会有部分重叠, 不能明显区分出来.

    • 选取3种粒子鉴别方法中FOM最优的设置, 表1展示了对比结果. 同种方法应用在不同放射源上, 得到的FOM值变化不大, 证明了每种方法都适用于在线实验测量, 用于区分不同能量的γ射线与α粒子. 脉冲宽度鉴别法和时间衰减常数鉴别法的FOM值低于快总成分比鉴别法, 原因主要是由于大体积的BaF2探测单元的信号波形有比较明显的涨落(如图9图12所示出现很多毛刺), 在进行脉冲宽度计算和衰减常数拟合时, 会严重影响特征值计算结果的统计分布, 使得α峰和γ峰的FWHM值变大, 从而降低FOM值.

      使用3种粒子鉴别方法后, 从放射源的能谱中可以看到α粒子明显减少, 同时也显现出了被α粒子计数淹没的22Na (1.27 MeV)的全能峰(图14(a)), 以及60Co的两条级联γ射线被同一BaF2探测单元探测到产生的2.5 MeV的γ峰(图14(c)). 从结果来看, 快总成分比(快成分5 ns, 总成分200 ns)的α粒子和γ射线的鉴别效果最好, 脉冲宽度(10%峰值)与时间衰减常数的鉴别效果差别不大, 这也与FOM值的计算结果一致.

    • 基于大体积的BaF2探测单元对γ射线的测量结果, 研究信号波形的采样长度对能量分辨率的影响, 确定2000 ns的采样长度得到的能量分辨率最好. 针对BaF2探测单元的信号波形特点, 采用快总成分比、脉冲宽度和时间衰减常数3种方法, 进行α粒子和γ射线的鉴别能力研究, 并使用品质因子FOM作为评估值. 通过放射源能谱和FOM值的比较, 快总成分比的α/γ鉴别能力优于其他两种方法.

      目前, 由40个BaF2探测单元组成的GTAF开展(n, γ)反应截面在线测量实验, 数据获取系统采用信号全波形采集的方式, 每次实验记录、传输、存储的数据量非常庞大, 不利于数据的存储和处理, 同时也延长了数据获取系统的死时间, 影响截面数据的不确定度. 本文验证了2000 ns的采样长度和快总成分比鉴别方法的优点, 下一步计划升级在线实验的数据获取系统, 只记录每个信号波形的过阈时间(用于飞行时间方法)、过阈后5 ns幅度积分值(快成分)、200 ns幅度积分值(总成分)、2000 ns幅度积分值(能量)、探测单元编号等一系列参数, 不用记录全波形信息, 使用这些信息就可以完成α粒子本底去除和(n, γ)反应的在线实验数据处理. 通过估算, 升级后的获取系统存储的数据量从118 MB/s降低到24 MB/s, 可以明显减少数据获取系统的死时间, 从而提高截面数据的精度.

    参考文献 (30)

目录

/

返回文章
返回