高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于虚拟刻度原理对放射性体源效率刻度方法的发展

许丰 郭建友 於国兵 徐平 张建

许丰, 郭建友, 於国兵, 徐平, 张建. 基于虚拟刻度原理对放射性体源效率刻度方法的发展[J]. 原子核物理评论, 2022, 39(3): 352-358. doi: 10.11804/NuclPhysRev.39.2021087
引用本文: 许丰, 郭建友, 於国兵, 徐平, 张建. 基于虚拟刻度原理对放射性体源效率刻度方法的发展[J]. 原子核物理评论, 2022, 39(3): 352-358. doi: 10.11804/NuclPhysRev.39.2021087
Feng XU, Jianyou GUO, Guobing YU, Ping XU, Jian Zhang. Development of Efficiency Calibration Method for Radioactive Volume Sources Based on Virtual Calibration Principle[J]. Nuclear Physics Review, 2022, 39(3): 352-358. doi: 10.11804/NuclPhysRev.39.2021087
Citation: Feng XU, Jianyou GUO, Guobing YU, Ping XU, Jian Zhang. Development of Efficiency Calibration Method for Radioactive Volume Sources Based on Virtual Calibration Principle[J]. Nuclear Physics Review, 2022, 39(3): 352-358. doi: 10.11804/NuclPhysRev.39.2021087

基于虚拟刻度原理对放射性体源效率刻度方法的发展

doi: 10.11804/NuclPhysRev.39.2021087
基金项目: 国家自然科学基金资助项目(11935001)
详细信息
    作者简介:

    许丰(1996−),男,安徽铜陵人,硕士研究生,从事辐射监测技术研究;E-mail: xufengid@sina.com

    通讯作者: 郭建友,E-mail: Jianyou@ahu.edu.cn
  • 中图分类号: TL81

Development of Efficiency Calibration Method for Radioactive Volume Sources Based on Virtual Calibration Principle

Funds: National Natural Science Foundation of China(11935001)
More Information
  • 摘要: 采用蒙特卡罗模拟计算的方法,标定能够适用于不同HPGe探测器上的体源虚拟点源位置,需要对点源、体源的探测效率进行模拟计算。通过 241Am、137Cs、60Co点源和体源研究了HPGe晶体尺寸、类型对它们的虚拟点源位置的影响,模拟结果表明241Am虚拟点源位置随着探测器尺寸、类型的不同呈现明显的差异性,说明虚拟刻度原理对于探测器表面的小体积样品测量在低能区间是不可取的,最后采用 137Cs、60Co源得出体源高度与虚拟点源位置的半经验公式。通过实验室两台HPGe探测系统对尺寸Φ70 mm×65 mm标准土壤体源的探测效率进行计算,与模拟探测效率值、虚拟点源效率表征结果对比分析,验证了拟合一组表征虚拟点源位置的半经验公式是可行的。实验结果表明,对$\gamma $能量区间在300~2 000 keV进行体源的虚拟刻度时,选择可靠的定标源即可建立体源和虚拟点源位置的关系。这为解决监测工作中样品与探测器之间重复进行探测效率校准的问题提供了新的途径。
  • 图  1  体源的虚拟点源刻度模型示意图(在线彩图)

    图  2  LabSOCS点源仿真计算示意图(在线彩图)

    (a) 为LabSOCS软件构建的点源模型;(b) 为点源计算模型的参数示意图。

    图  3  点源在探测器轴线上方不同位置的探测效率(虚线为指数函数拟合) (在线彩图)

    图  4  LabSOCS软件对圆柱体源构建的几何模型(在线彩图)

    (a)为LabSOCS软件构建的柱状体源模型;(b)为体源计算模型的参数示意图。

    图  5  不同尺寸及类型的探测器对体源虚拟点位置的比较(在线彩图)

    图  6  体源的虚拟源位置关系拟合(在线彩图)

    图  7  尺寸Ф70 mm×65 mm标准体源的实验测试能谱(在线彩图)

    图  8  虚拟源刻度方法对体源效率刻度实验验证结果对比(Exp:体源探测效率实验测量值;VPS:虚拟点源效率刻度值;LabSOCS:软件对体源模拟计算值) (在线彩图)

    表  1  点源位置(XP)与探测效率关系曲线拟合

    核素探测器编号c0αbR2
    241Am(59.54 keV)Coaxial-13.41×10−34.57×10−2−4.10×10−21.00
    Coaxial-26.46×10−49.31×10−2−3.11×10−21.00
    Coaxial-3−9.58×10−38.97×10−2−2.02×10−20.99
    Coaxial-48.02−7.977.09×10−50.98
    Planar-12.08×10−23.24×10−1−4.81×10−21.00
    Planar-21.82×10−23.65×10−1−3.94×10−21.00
    Planar-32.36×10−23.51×10−1−3.73×10−21.00
    137Cs(661.66 keV)Coaxial-14.65×10−31.65×10−2−3.22×10−21.00
    Coaxial-28.21×10−32.79×10−2−3.30×10−21.00
    Coaxial-31.29×10−24.26×10−2−3.17×10−21.00
    Coaxial-41.71×10−25.66×10−2−3.04×10−21.00
    Planar-15.09×10−31.89×10−2−3.51×10−21.00
    Planar-29.73×10−33.47×10−2−3.50×10−21.00
    Planar-31.27×10−24.30×10−2−3.28×10−21.00
    60Co(1 332.49 keV)Coaxial-12.12×10−31.90×10−2−6.24×10−21.00
    Coaxial-22.89×10−33.02×10−2−5.44×10−21.00
    Coaxial-34.75×10−34.64×10−2−4.84×10−21.00
    Coaxial-45.03×10−35.99×10−2−3.98×10−21.00
    Planar-11.39×10−31.83×10−2−6.00×10−21.00
    Planar-22.64×10−33.37×10−2−5.44×10−21.00
    Planar-33.82×10−33.94×10−2−4.82×10−21.00
    下载: 导出CSV

    表  2  不同尺寸及类型探测器对不同高度体源(D=70 mm)的探测效率${\varepsilon }^{}_{\mathrm{V}}(H,R)$计算值

    核素探测器编号H / mm
    5152535455565
    241Am(59.54 keV)Coaxial-11.63×10−21.22×10−29.43×10−37.55×10−36.22×10−35.24×10−34.51×10−3
    Coaxial-24.23×10−23.11×10−22.39×10−21.90×10−21.56×10−21.32×10−21.13×10−2
    Coaxial-34.42×10−23.27×10−22.53×10−22.02×10−21.66×10−21.40×10−21.21×10−2
    Coaxial-43.52×10−22.66×10−22.08×10−21.68×10−21.39×10−21.17×10−21.01×10−2
    Planar-11.52×10−11.05×10−17.83×10−26.14×10−24.99×10−24.19×10−23.59×10−2
    Planar-22.02×10−11.39×10−11.04×10−18.13×10−26.61×10−25.54×10−24.74×10−2
    Planar-32.17×10−11.51×10−11.13×10−18.91×10−27.27×10−26.11×10−25.23×10−2
    137Cs(661.66 keV)Coaxial-11.87×10−21.47×10−21.20×10−21.00×10−28.59×10−37.48×10−36.60×10−3
    Coaxial-23.20×10−22.50×10−22.04×10−21.71×10−21.47×10−21.28×10−21.13×10−2
    Coaxial-34.94×10−23.91×10−23.21×10−22.71×10−22.33×10−22.04×10−21.81×10−2
    Coaxial-46.59×10−25.26×10−24.35×10−23.68×10−23.18×10−22.78×10−22.47×10−2
    Planar-12.10×10−21.61×10−21.29×10−21.07×10−29.07×10−37.86×10−36.90×10−3
    Planar-23.90×10−23.00×10−22.41×10−22.01×10−21.71×10−21.48×10−21.31×10−2
    Planar-34.94×10−23.87×10−23.16×10−22.66×10−22.28×10−21.99×10−21.76×10−2
    60Co(1 332.49 keV)Coaxial-11.04×10−28.26×10−36.81×10−35.76×10−34.98×10−34.79×10−33.87×10−3
    Coaxial-21.82×10−21.45×10−21.19×10−21.01×10−28.69×10−37.63×10−36.77×10−3
    Coaxial-33.07×10−22.47×10−22.05×10−21.74×10−21.51×10−21.33×10−21.19×10−2
    Coaxial-44.29×10−23.47×10−22.89×10−22.47×10−22.15×10−21.90×10−21.69×10−2
    Planar-11.01×10−27.82×10−36.33×10−35.29×10−34.53×10−33.94×10−33.48×10−3
    Planar-22.05×10−21.59×10−21.29×10−21.09×10−29.31×10−38.12×10−37.18×10−3
    Planar-32.73×10−22.17×10−21.79×10−21.51×10−21.31×10−21.15×10−21.02×10−2
    下载: 导出CSV
  • [1] NOTEA A. Nucl Instr and Meth, 1971, 91(3): 513. doi:  10.1016/S0029-554X(71)80031-9
    [2] MAHLING S, ORION I, ALFASSI Z B. Nucl Instr and Meth A, 2006, 557(2): 544. doi:  10.1016/j.nima.2005.11.118
    [3] HUYNH D C, LE T N T, NGUYEN V H, et al. Applied Radiation and Isotopes, 2020, 162: 109179. doi:  10.1016/j.apradiso.2020.109179
    [4] SEDIGHEH S, ZAHRA M, MEHROOSH K, et al. Scientific Research Journal, 2020, 17(1): 15. doi:  10.24191/srj.v17i1.9345
    [5] 田自宁, 刘文彪, 陈伟, 等. 强粒子与激光束, 2018, 30(3): 0340011. doi:  10.11884/HPLPB201830.170278

    TIAN Zining, LIU Wenbiao, CHEN Wei, et al. High Power Laser And Particle Beams, 2018, 30(3): 0340011. (in Chinese) doi:  10.11884/HPLPB201830.170278
    [6] DEBERTIN K, HELMER R G. Nucl Instr and Meth, 1989, 44(2): 246. doi:  10.1016/0168-583X(89)90436-9
    [7] JULIO C A. Applied Radiation and Isotopes, 2008, 66(8): 1123. doi:  10.1016/j.apradiso.2007.12.007
    [8] BYUN J I, CHOI H Y, YUN J Y. Applied Radiation and Isotopes, 2010, 68(2): 370. doi:  10.1016/j.apradiso.2009.10.054
    [9] 聂鹏, 倪邦发, 田伟之, 等. 原子能科学与技术, 2010, 44(2): 206. doi:  10.7538/yzk.2010.44.02.0206

    NIE Peng, NI Bangfa, TIAN Weizi, et al. Atomic Energy Science and Technology, 2010, 44(2): 206. (in Chinese) doi:  10.7538/yzk.2010.44.02.0206
    [10] 熊文彬, 仇春华, 段天英, 等. 原子能科学与技术, 2011, 45(8): 999. doi:  10.7538/yzk.2011.45.08.0999

    XIONG Wenbin, QIU Chunhhua, DUAN Tianying, et al. Atomic Energy Science and Technology, 2011, 45(8): 999. (in Chinese) doi:  10.7538/yzk.2011.45.08.0999
    [11] 田自宁. 基于虚拟源原理的放射性体源无源效率刻度技术研究[D]. 北京: 清华大学, 2015.

    TIAN Zining. The Reasearch on Radioactive Volume Source-less Efficiency Calibration Technology Based on the Virtual Source Theory[D]. Beijing: Tsinghua University, 2015. (in Chinese)
    [12] TIAN Zining, OUYANG Xiaoping, LIU Linyue. Journal of Radioanalytical and Nuclear Chemistry, 2015, 304: 905. doi:  10.1007/s10967-014-3907-2
    [13] 田自宁, 陈伟, 韩斌, 等. 物理学报, 2016, 65(6): 062901. doi:  10.7498/aps.65.062901

    TIAN Zining, CHEN Wei, HAN Bin, et al. Acta Phys Sin, 2016, 65(6): 062901. (in Chinese) doi:  10.7498/aps.65.062901
    [14] 田自宁, 欧阳晓平, 陈伟, 等. 物理学报, 2019, 68(23): 232901. doi:  10.7498/aps.68.20191095

    TIAN Zining, OUYANG Xiaoping, CHEN Wei, et al. Acta Phys Sin, 2019, 68(23): 232901. (in Chinese) doi:  10.7498/aps.68.20191095
    [15] 董汝昆, 吴绍华, 王柯, 等. 红外技术, 2021, 43(5): 510.

    DONG Rukun, WU Shaohua, WANG Ke, et al. Infrared Technology, 2021, 43(5): 510. (in Chinese)
    [16] 白尔隽, 郑志鹏, 高德喜, 等. 原子核物理评论, 2016, 33(1): 52. doi:  10.11804/NuclPhysRev.33.01.052

    BAI Erjuan, ZHENG Zhipeng, GAO Dexi, et al. Nuclear Physics Review, 2016, 33(1): 52. (in Chinese) doi:  10.11804/NuclPhysRev.33.01.052
    [17] 罗奇, 胡世鹏, 孙慧斌, 等. 核技术, 2011, 34(12): 901.

    LUO Qi, HU Shiping, SUN Huibin, et al. Nuclear Techniques, 2011, 34(12): 901. (in Chinese)
    [18] CHEN Zhilin, SONG Guoyang, MU Long, et al. Chinese Physics C, 2010, 34(5): 560. (in Chinese) doi:  10.1088/1674-1137/34/5/009
    [19] 俞森标, 姜让荣. 辐射防护, 1998, 18(4): 314. doi:  10.3321/j.issn:1000-8187.1998.04.001

    YU Senbiao, JIANG Ranrong. Radiation Protection, 1998, 18(4): 314. (in Chinese) doi:  10.3321/j.issn:1000-8187.1998.04.001
    [20] VENKATARAMAN R, BRONSON F. Nucl Instr and Meth, 1999, 422(1-3): 450. doi:  10.1016/S0168-9002(98)01115-2
    [21] VENKATARAMAN R, BRONSON F, Atrashkevich V, et al. Journal of Radioanalytical and Nuclear Chemistry, 2005, 264(1): 213. doi:  10.1007/s10967-005-0696-7
    [22] HJ 1127-2020 应急监测中环境样品γ核素测量技术规范[S]. 北京: 中国环境科学出版社, 2020.

    HJ 1127-2020 Technical Specification for γ Nuclide Measurement of Environmental Samples in Emergency Monitoring[S]. Beijing: China Environmental Science Press, 2020.
  • 加载中
图(8) / 表 (2)
计量
  • 文章访问数:  380
  • HTML全文浏览量:  118
  • PDF下载量:  21
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-11-25
  • 修回日期:  2022-01-25
  • 刊出日期:  2022-09-20

基于虚拟刻度原理对放射性体源效率刻度方法的发展

doi: 10.11804/NuclPhysRev.39.2021087
    基金项目:  国家自然科学基金资助项目(11935001)
    作者简介:

    许丰(1996−),男,安徽铜陵人,硕士研究生,从事辐射监测技术研究;E-mail: xufengid@sina.com

    通讯作者: 郭建友,E-mail: Jianyou@ahu.edu.cn
  • 中图分类号: TL81

摘要: 采用蒙特卡罗模拟计算的方法,标定能够适用于不同HPGe探测器上的体源虚拟点源位置,需要对点源、体源的探测效率进行模拟计算。通过 241Am、137Cs、60Co点源和体源研究了HPGe晶体尺寸、类型对它们的虚拟点源位置的影响,模拟结果表明241Am虚拟点源位置随着探测器尺寸、类型的不同呈现明显的差异性,说明虚拟刻度原理对于探测器表面的小体积样品测量在低能区间是不可取的,最后采用 137Cs、60Co源得出体源高度与虚拟点源位置的半经验公式。通过实验室两台HPGe探测系统对尺寸Φ70 mm×65 mm标准土壤体源的探测效率进行计算,与模拟探测效率值、虚拟点源效率表征结果对比分析,验证了拟合一组表征虚拟点源位置的半经验公式是可行的。实验结果表明,对$\gamma $能量区间在300~2 000 keV进行体源的虚拟刻度时,选择可靠的定标源即可建立体源和虚拟点源位置的关系。这为解决监测工作中样品与探测器之间重复进行探测效率校准的问题提供了新的途径。

English Abstract

许丰, 郭建友, 於国兵, 徐平, 张建. 基于虚拟刻度原理对放射性体源效率刻度方法的发展[J]. 原子核物理评论, 2022, 39(3): 352-358. doi: 10.11804/NuclPhysRev.39.2021087
引用本文: 许丰, 郭建友, 於国兵, 徐平, 张建. 基于虚拟刻度原理对放射性体源效率刻度方法的发展[J]. 原子核物理评论, 2022, 39(3): 352-358. doi: 10.11804/NuclPhysRev.39.2021087
Feng XU, Jianyou GUO, Guobing YU, Ping XU, Jian Zhang. Development of Efficiency Calibration Method for Radioactive Volume Sources Based on Virtual Calibration Principle[J]. Nuclear Physics Review, 2022, 39(3): 352-358. doi: 10.11804/NuclPhysRev.39.2021087
Citation: Feng XU, Jianyou GUO, Guobing YU, Ping XU, Jian Zhang. Development of Efficiency Calibration Method for Radioactive Volume Sources Based on Virtual Calibration Principle[J]. Nuclear Physics Review, 2022, 39(3): 352-358. doi: 10.11804/NuclPhysRev.39.2021087
    • $\gamma $辐射探测仪器在环境放射性监测、出入境检验检疫、核应急等诸多领域有着广泛的应用。日常监测分析中也会碰到大量的体源样品,由于体源的几何形状不一以及介质本身自吸收系数的影响,致使探测效率随体源与探测器之间的空间位置的变化而有所差异。这就意味着在实际工作中需要对各类样品与探测器之间重复进行探测效率校准,成本高昂且耗时费力。随着蒙特卡罗方法的发展,基于无源效率刻度方法对放射源的探测效率计算逐步发展起来,相对于制作标准源模型,无源效率刻度方法的优势在于方便快捷,可以根据体源形状计算出探测效率,然而其缺点在于对探测器内部参数、样品的组成成分等信息仿真不够精确。

      针对上述问题,1971年Notea等[1]提出虚拟点探测器模型的概念,该模型在处理$\gamma $探测器效率刻度的过程中,被国内外有关学者广泛使用。通过蒙特卡罗模拟的方法[2-3]或者神经网络算法[4]研究不同探测器的有效作用深度[5](Effective Interaction Depth, EID),该模型可实现对均匀面源、体源的效率刻度计算[6-7],此外也有学者利用该方法研究放射源的搜寻、定位等[8],这都极大地推动了该技术的进展。近年来,国内学者在该领域也取得了不错的进展[9-10],其中西北核技术所Tian[11]将虚拟点探测器刻度技术推广到现场测量工作中,并由此提出了虚拟刻度原理。虚拟刻度原理在面源及体源的虚拟点源位置计算 [12-13]、放射性热点以及污染区域范围估算等方面[14]取得了较为成功的应用,展现出其在就地测量中不依赖于探测器进行重复表征的一定优势。考虑到在实验室环境样品监测工作中,常常会遇到由于$\gamma $辐射探测器种类及晶体尺寸多样[15-16],得到的效率曲线经验公式往往仅针对于所研究的特定探测器可以使用[17-18],即使遇到此类测量场景,也常会由于实验采样不足而需要进行几何修正[19]。综上所述,由于$\gamma $辐射探测器种类及晶体尺寸多样,关于单台探测器的虚拟点源测量方法及经验公式是否适用多台仪器就显得尤为重要。

      本文基于虚拟刻度原理,将分别建立多台HPGe同轴型及平面型探测器仿真模型,对体源中$\gamma $放射性核素241Am、137Cs、60Co的虚拟点源位置进行表征,研究体源虚拟源位置的表征结果对探测器类型、能量的依赖性。

    • 使用$\gamma $辐射探测器测量时,可以认为是射线与探测器等效点之间的作用,因此对于不同能量的$\gamma $射线,在探测器中心轴线上存在唯一代表点位置,该位置为探测器的有效作用深度,即EID。然后通过反比平方定律去计算任意空间位置的点源探测效率,这一概念被称为虚拟点探测器模型。而对于单个放射性面源或放射性体源样品采用γ辐射探测器进行辐射测量时,根据虚拟点探测器原理同样也可以对放射源样品虚拟化,因此在探测器中心对称轴线上,也应该有一个唯一的代表点位置(Xp),该位置的放射性点源全能峰探测效率与所对应的放射性体源或面源的探测效率相等,将这一虚拟化的点源叫做放射性面源或放射性体源的虚拟点源(Virtual Point Source, VPS)[11]。因此利用虚拟刻度原理对面源或体源探测效率刻度的关键是虚拟点源位置的确定,原理示意图如图1所示,XVPS为半径为R、高度H体源的虚拟点源位置,${\varepsilon }^{}_{\mathrm{P}}\left({X}_{\mathrm{V}\mathrm{P}\mathrm{S}}\right)$为虚拟点源(VPS)的效率,${\varepsilon }^{}_{\mathrm{V}}(H,R)$为体源效率,EID为探测器的有效作用深度,探测器端面中心为坐标轴原点。

      图  1  体源的虚拟点源刻度模型示意图(在线彩图)

    • 本文以探测器端面中心处为坐标轴原点,将标准点源置于探测器表面对称轴上不同位置,获取其探测效率,然后采用指数函数拟合点源位置与探测效率的关系,具体函数形式如式(1):

      $$ {\varepsilon }^{}_{\mathrm{P}}\left(x\right)=a{{\rm e}}^{b \boldsymbol \cdot{ X}_{\mathrm{P}}}+{c}_{0} , $$ (1)

      式中:${\varepsilon }^{}_{\mathrm{P}}\left(x\right)$表示点源在不同位置x时的探测效率;a、b表示拟合参数;c0为指数函数常数项;XP表示点源位置。然后将实验或理论计算得到的体源的探测效率${\varepsilon }^{}_{\mathrm{V}}(H,R)$等于式(1)中应变量$ {\varepsilon }_{\mathrm{P}}\left(x\right) $,即可推出它的虚拟点源位置($ {X}_{\mathrm{V}\mathrm{P}\mathrm{S}} $)。虚拟点源位置($ {X}_{\mathrm{V}\mathrm{P}\mathrm{S}} $)计算公式如式(2):

      $$ {X}_{\mathrm{V}\mathrm{P}\mathrm{S}}=\frac{1}{b}\mathrm{ln}\left(\frac{{\varepsilon }_{{\rm{V}}}\left(H,R\right)-{c}_{0}}{a}\right) \text{,} $$ (2)

      式中:下标V表示体源;${\varepsilon }^{}_{\mathrm{V}}(H,R)$表示为体源探测效率,可以通过实验或模拟方法得到体源的探测效率。值得注意的是,日常辐射监测中由于标准体源的尺寸一般为标称,因此R为定值,而对于实际采样过程中,常常会遇到采样量不足导致几何高度不一致的情况,因此本文仅考虑体源高度H即可。最后计算出体源不同几何高度时虚拟点源位置,通过指数函数拟后,分析虚拟点源位置与体源高度的关系:

      $$ {X}_{\mathrm{V}\mathrm{P}\mathrm{S}}={A\boldsymbol\cdot {\rm e}}^{B\boldsymbol\cdot H}+C , $$ (3)

      式中:$ {X}_{\mathrm{V}\mathrm{P}\mathrm{S}} $表示不同尺寸的体源对应的虚拟点源位置;A、B、C为拟合参数;H为体源高度。

    • 探测效率的模拟计算采用无源效率刻度软件LabSOCS,该软件是由法国Canberra公司与美国Los Alamos国家实验室共同开发的产品,国内外学者或机构将该软件与传统的有源效率刻度技术进行对比研究,结论表明与传统方法的分析结果一致[20],LabSOCS主要原理是通过建立MCNP模型,用一系列可溯源的多能量点源在多个位置对探测器进行验证以及仿真模型的优化[21]。考虑到测量工作中的实际需求,本工作分别模拟了四种同轴型探测器,晶体直径分别为Coaxial-1: Φ50.00 mm×48.05 mm、Coaxial-2: Φ60.06 mm×45.28 mm、Coaxial-3: Φ70.00 mm×60.00 mm、Coaxial-4: Φ81.06 mm×65.94 mm;三种平面型HPGe探测器,晶体直径分别为Planar-1: Φ59.94 mm×19.58 mm、Planar-2: Φ69.98 mm×25.18 mm、Planar-3: Φ78.44 mm×30.66 mm,然后分别模拟探测器轴线上方不同位置处241Am(59 keV)、137Cs(661.66 keV)、60Co(1 332.49 keV)点源的探测效率。LabSOCS软件对点源仿真模型如图2所示,吸收介质层厚度d4.1=0,d5.1=0,源到探测器表面距离d6.1=5, 15, ..., 65 mm,模拟计算过程中统计误差小于2%。

      图  2  LabSOCS点源仿真计算示意图(在线彩图)

      结合多种尺寸及类型的无源表征文件,我们模拟计算了点源在不同位置处(点源与探测器端面中心处距离为XP,单位为mm)探测器对241Am、137 Cs、60Co的探测效率,结果图3所示。

      图  3  点源在探测器轴线上方不同位置的探测效率(虚线为指数函数拟合) (在线彩图)

      通过图3可以看出,点源探测效率随着轴向位置变化呈现明显的指数变化规律,同时结合文献[11],本实验选取式(1)对241Am、137Cs、60Co点源探测效率进行拟合,拟合的详细参数abc0以及拟合优度R2表1所列。

      表 1  点源位置(XP)与探测效率关系曲线拟合

      核素探测器编号c0αbR2
      241Am(59.54 keV)Coaxial-13.41×10−34.57×10−2−4.10×10−21.00
      Coaxial-26.46×10−49.31×10−2−3.11×10−21.00
      Coaxial-3−9.58×10−38.97×10−2−2.02×10−20.99
      Coaxial-48.02−7.977.09×10−50.98
      Planar-12.08×10−23.24×10−1−4.81×10−21.00
      Planar-21.82×10−23.65×10−1−3.94×10−21.00
      Planar-32.36×10−23.51×10−1−3.73×10−21.00
      137Cs(661.66 keV)Coaxial-14.65×10−31.65×10−2−3.22×10−21.00
      Coaxial-28.21×10−32.79×10−2−3.30×10−21.00
      Coaxial-31.29×10−24.26×10−2−3.17×10−21.00
      Coaxial-41.71×10−25.66×10−2−3.04×10−21.00
      Planar-15.09×10−31.89×10−2−3.51×10−21.00
      Planar-29.73×10−33.47×10−2−3.50×10−21.00
      Planar-31.27×10−24.30×10−2−3.28×10−21.00
      60Co(1 332.49 keV)Coaxial-12.12×10−31.90×10−2−6.24×10−21.00
      Coaxial-22.89×10−33.02×10−2−5.44×10−21.00
      Coaxial-34.75×10−34.64×10−2−4.84×10−21.00
      Coaxial-45.03×10−35.99×10−2−3.98×10−21.00
      Planar-11.39×10−31.83×10−2−6.00×10−21.00
      Planar-22.64×10−33.37×10−2−5.44×10−21.00
      Planar-33.82×10−33.94×10−2−4.82×10−21.00
    • 体源仿真模型采用实验室常见的核素均匀分布的固体样品,体源介质为标准土壤,元素占比分别为Si: 32.7 %、Al: 13.2 %、Fe: 3.4 %、O: 50.6 %,密度1.43 g/cm3。根据应急监测标准中对环境样品$\gamma $核素测量的要求[22],建议样品尺寸为Ф70 mm×65 mm,因此实验研究的固体放射性体源直径D=70 mm;样品盒材质为聚乙烯,厚度2.5 mm,密度0.91 g/cm3。计算样品探测效率时,体源直径不变,改变体源高度,高度调整范围为5~65 mm,变化梯度10 mm,样品贴近探测器表面,与探测器同轴放置,通过LabSOCS软件计算上述样品源的探测效率,分别得到不同体源高度下$\gamma $核素241Am、137Cs、60Co的探测效率值。LabSOCS软件对圆柱体源的几何配置图如图4,吸收介质层厚度d4.1=0,d5.1=0,d1.1=2.5 mm,d2.1=2.5 mm,材质为聚乙烯,d1.2=70 mm,d3.1=5, 15…, 65 mm,放置于探测器表面(d6.1=0),样品材质为上述描述的标准土壤,模拟计算过程中统计误差小于2%。

      图  4  LabSOCS软件对圆柱体源构建的几何模型(在线彩图)

      通过上述参数设置,结合四种同轴型探测器、三种平面探测器的仿真计算模型,通过LabSOCS模拟计算得到体源样品中241Am、137Cs、60Co的探测效率[$ {\varepsilon }_{\rm V}(H,R) $],体源高度区间为5~65 mm,直径D=70 mm。具体计算数据如表2所列。

      表 2  不同尺寸及类型探测器对不同高度体源(D=70 mm)的探测效率${\varepsilon }^{}_{\mathrm{V}}(H,R)$计算值

      核素探测器编号H / mm
      5152535455565
      241Am(59.54 keV)Coaxial-11.63×10−21.22×10−29.43×10−37.55×10−36.22×10−35.24×10−34.51×10−3
      Coaxial-24.23×10−23.11×10−22.39×10−21.90×10−21.56×10−21.32×10−21.13×10−2
      Coaxial-34.42×10−23.27×10−22.53×10−22.02×10−21.66×10−21.40×10−21.21×10−2
      Coaxial-43.52×10−22.66×10−22.08×10−21.68×10−21.39×10−21.17×10−21.01×10−2
      Planar-11.52×10−11.05×10−17.83×10−26.14×10−24.99×10−24.19×10−23.59×10−2
      Planar-22.02×10−11.39×10−11.04×10−18.13×10−26.61×10−25.54×10−24.74×10−2
      Planar-32.17×10−11.51×10−11.13×10−18.91×10−27.27×10−26.11×10−25.23×10−2
      137Cs(661.66 keV)Coaxial-11.87×10−21.47×10−21.20×10−21.00×10−28.59×10−37.48×10−36.60×10−3
      Coaxial-23.20×10−22.50×10−22.04×10−21.71×10−21.47×10−21.28×10−21.13×10−2
      Coaxial-34.94×10−23.91×10−23.21×10−22.71×10−22.33×10−22.04×10−21.81×10−2
      Coaxial-46.59×10−25.26×10−24.35×10−23.68×10−23.18×10−22.78×10−22.47×10−2
      Planar-12.10×10−21.61×10−21.29×10−21.07×10−29.07×10−37.86×10−36.90×10−3
      Planar-23.90×10−23.00×10−22.41×10−22.01×10−21.71×10−21.48×10−21.31×10−2
      Planar-34.94×10−23.87×10−23.16×10−22.66×10−22.28×10−21.99×10−21.76×10−2
      60Co(1 332.49 keV)Coaxial-11.04×10−28.26×10−36.81×10−35.76×10−34.98×10−34.79×10−33.87×10−3
      Coaxial-21.82×10−21.45×10−21.19×10−21.01×10−28.69×10−37.63×10−36.77×10−3
      Coaxial-33.07×10−22.47×10−22.05×10−21.74×10−21.51×10−21.33×10−21.19×10−2
      Coaxial-44.29×10−23.47×10−22.89×10−22.47×10−22.15×10−21.90×10−21.69×10−2
      Planar-11.01×10−27.82×10−36.33×10−35.29×10−34.53×10−33.94×10−33.48×10−3
      Planar-22.05×10−21.59×10−21.29×10−21.09×10−29.31×10−38.12×10−37.18×10−3
      Planar-32.73×10−22.17×10−21.79×10−21.51×10−21.31×10−21.15×10−21.02×10−2
    • 表1中点源位置效率函数关系的拟合参数代入式(1),结合表2中标准土壤体源的探测效率模拟计算值,可以联立求解该类体源的虚拟点源位置,具体过程由式(2)分析得到。最后针对上述常见晶体类型及尺寸的探测器,绘制其体源高度与虚拟点源位置的关系图,结果如图5所示。

      图  5  不同尺寸及类型的探测器对体源虚拟点位置的比较(在线彩图)

      图5可以看出,对241Am、137Cs、60Co$\gamma $射线进行虚拟点位置计算,对于低能$\gamma $射线,241Am(59.54 keV)的虚拟点位置与探测器尺寸、类型相关性最大,且随体源深度的延长,探测器表征的虚拟点位置差异逐渐增加,具体表现为Hv=5 mm时,标准偏差σ=6.16;Hv=65 mm时,σ=11.43。对于中高能$\gamma $射线,137Cs(661.66 keV)、60Co(1 332.49 keV)的虚拟点位置则与探测器相关性不大,且随着深度的变化,其位置表征结果在各探测器上仍趋于一致,具体表现为Hv=5 mm时,标准偏差$\sigma $=0.97, 0.82;Hv=65 mm时,$\sigma $=1.50, 1.45。这一差异,主要是由于低能区射线线衰减系数较大,受射线在探测晶体中衰减距离的影响明显,若适当增加样品距离会降低这种差异性。而对于实验室日常测量工作,环境样品一般贴近于探测器表面的情况,低能射线使用虚拟刻度原理需要根据探测器的类型、尺寸进行再次表征,而对于中高能射线是不再需要的。最后对图5中不同能量、晶体尺寸、体源高度得到的虚拟点源计算结果取平均值,拟合得到体源高度与虚拟点源位置(XVPS-H)的半经验公式,如图6所示,图中误差棒为上述谱仪系统对体源的虚拟点源位置计算结果的均值与各探测器上的计算结果之间的均方差。

      图  6  体源的虚拟源位置关系拟合(在线彩图)

    • 为了验证上述方法得到的经验公式在快速效率刻度中的可行性,我们通过两套Canberra公司生产的HPGe谱仪测量系统开展实验验证工作,一套为同轴P型HPGe探测器系统,型号GC4019,对60Co的1 332.49 keV的能量分辨率为1.78 keV,相对探测效率为41.80%;另一套为低本底(Ultra-Low Background)宽能型HPGe探测器系统,型号BE5030,能量分辨率为1.619 keV,相对探测效率为47.3%。模拟计算中选取上述两套测量系统配套的蒙特卡罗方法表征过的探头文件,将表征探头导入商用无源效率刻度软件LabSOCS中构建仿真计算模型。实验测量以及模拟计算均采用国标规定尺寸Ф70 mm×65 mm[22](活度信息:49K 1 681.9 Bq、232Th 641.0 Bq、238U 978.8 Bq、226Ra 961.5 Bq)的标准源,标准源介质为模拟土壤,为CuO(150 g)、SiO2(110 g)混合搅拌制备,密度为1.74 g/cm3。在300~1 000 keV区间内采用图6中的公式(a)标定,对1 000~2 000 keV区间采用图6中的公式(b)标定。模拟及实验过程中的分析核素为:228Ac(338.32 keV)、214Bi(609.32 keV)、228Ac(911.20 keV)、234mPa(1 001.10 keV)、40K(1 460.75 keV)、214Bi(1 764.49 keV)。图7为Ф70 mm×65 mm在两套HPGe谱仪系统(GC4019、BE5030)上的测试能谱。

      图  7  尺寸Ф70 mm×65 mm标准体源的实验测试能谱(在线彩图)

      通过图7的谱图分析计算,得到了Ф70 mm×65 mm标准体源在探测器GC4019、BE5030上的效率实验值,探测效率实验值的计算方法为$\varepsilon $=n/(ηAk),n$\gamma $特征能量计数率,η为射线发射强度,k为衰变修正因子,A为点源活度,$\varepsilon $为探测效率实验计算值。同时采用模拟方法,获得上述体源的效率计算值,模型建构方法参照图4所示。虚拟刻度方法中虚拟源位置结合图6可以给出:300 ~1 000 keV区间内,尺寸Ф70 mm×65 mm的虚拟源位置为40.9 mm,1 000 ~2 000 keV区间内虚拟源位置为37.8 mm,将实验与模拟计算的结果整理如图8所示。

      图  8  虚拟源刻度方法对体源效率刻度实验验证结果对比(Exp:体源探测效率实验测量值;VPS:虚拟点源效率刻度值;LabSOCS:软件对体源模拟计算值) (在线彩图)

      图8可以看出,采用137Cs、60Co点源分别标定300~1 000 keV、1 000~2 000 keV能区的虚拟点源位置,VPS刻度方法得到的探测效率与体源模拟效率和实验值基本一致。同时,图8中对228Ac的338.32 keV实验结果偏差较大,BE5030探测器上,体源实验计算结果与虚拟源计算值在338 keV偏差达到11.2%;GC4019探测器上,在338 keV两者偏差达到18.4%。低能区出现较大偏差的原因主要有如下两个方面:(1) 选取137Cs作为300~1 000 keV能区定标源对于低能段效果稍差,今后可对该能段再进行细分,将有利于改善效率刻度精度;(2) 虚拟点源位置的点源效率由模拟计算得到,然后对多套系统上的计算结果进行均值化处理,而该方法的前提是仿真计算模型的准确性,但实际情况是探测器内部死层等参数未知(值得注意的是,型号GC4019探测器购置于2005年3月,型号BE5030探测器购置于2020年1月,表征文件探头与GC4019的差异性会更大)。在利用VPS刻度方法对效率计算的过程中,分别采用137Cs(661.66 keV)、60Co(1 332.49 keV)标定300 ~1 000 keV、1 000 ~2 000 keV能区的虚拟点源,这从理论上来说是存在误差的,但图8中VPS结果与LabSOCS模拟和实验结果的一致性也说明该误差是比较小的。综上所述,通过对虚拟点源位置的比较、对能区的效率刻度等研究,表明虚拟源的位置具有不依赖于特定的探测器进行重复表征的优势,选取合适的点源用于表征一定能区范围的探测效率也是可行的,可以适用于放射性监测中的放射源测量工作。由于本文着重论述实验室测量工作,样品一般均是贴合探测器表面测量,今后如果能通过实际混合点源对源项虚拟点源位置进行表征效果将更好,有利于提高虚拟源原理对源项效率刻度的精确度。

    • 本文通过将137Cs、60Co点源分别作为300 ~1 000 keV能区、1 000~2 000 keV能区虚拟点源位置的标定源,提出了适用于300 ~2 000 keV能区的虚拟点源位置半经验公式。对于体源的虚拟点位置进行多方面研究,通过模拟计算点源及体源探测效率得到的数据为计算基础,比较了多种晶体尺寸的HPGe探测器上体源虚拟点位置关系,得出探测器类型的改变,对于体源中高能$ \gamma $射线的虚拟点位置表征结果影响不大。最后拟合给出了体源高度与其虚拟点位置的半经验关系式,模拟和实验结果表明该计算方法在300 keV以上能区可广泛适用于多种类型的探测器,为辐射环境应急监测的准确、快速效率刻度提供一条新的途径。

      在今后的样品测量工作中,可以使用合适的标准点源,如本文137Cs、60Co表征出一组虚拟点源的位置经验公式,用点源替代标准面源或者标准体源进行探测效率刻度,从而避免了体源批量制作、放射性废物处理困难、运输及携带不便等诸多问题,同时对日常辐射环境监测工作具有参考价值,并具有一定应用前景。

参考文献 (22)

目录

    /

    返回文章
    返回