AEP Advances in Environmental Protection 2164-5485 Scientific Research Publishing 10.12677/AEP.2022.125130 AEP-57106 AEP20220500000_98767106.pdf 地球与环境 基于NaI(Tl)谱仪的3种γ能谱解谱方法比较 Comparative Study of Three γ Spectrum Solution Methods Based on NaI(Tl) Spectrometer 胡 有华 2 1 王 建超 2 1 新疆维吾尔自治区辐射环境监督站，新疆 乌鲁木齐 null 01 09 2022 12 05 1049 1055 © Copyright 2014 by authors and Scientific Research Publishing Inc. 2014 This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

目的：NaI(Tl)谱仪相对分辨率较低，通过比较3种解谱方法并选出最优解谱方法，提高γ能谱测量值的准确性和可靠性。方法：先用相对比较、逆矩阵、最小二乘法3种分析方法对已知活度的混合标准源进行活度计算，再比较计算活度与标准源活度的相对偏差。结果：使用相对比较法、逆矩阵法、最小二乘法计算出3种核素相对偏差分别为： 226Ra (13.62%、14.71%、9.12%)、 232Th (14.43%、20.59%、13.17%)、 40K (2.18%、4.91%、4.84%)。结论：使用最小二乘法对 226Ra和 232Th 2种核素的计算活度相对偏差最小，优于相对比较法与逆矩阵法；40K核素的活度计算采用相对比较法更合适。 Objective: The resolution of NaI(Tl) spectrometer is relatively low. By comparing the three spec-trum solution methods and selecting the best solution method, the accuracy and reliability of the measured values of the γ energy spectrum can be improved. Methods: First use relative compari-son, inverse matrix, and least square analysis methods to calculate the activity of the mixed standard source with known activity, and then compare the relative deviation between the calcu-lated activity and the standard source activity. Results: The relative deviations of the three nuclides calculated using the relative comparison method, the inverse matrix method, and the least square method: 226Ra (13.62%, 14.71%, 9.12%), 232Th (14.43%, 20.59%, 13.17%), 40K (2.18%, 4.91%, 4.84%). Conclusions: Using the least square method to calculate the relative deviation of the activities of 226Ra and 232Th nuclides is the smallest, which is better than the relative comparison method and the inverse matrix method. The activity calculation of 40K nuclides is more suitable for the relative comparison method.

γ能谱，相对比较，逆矩阵，最小二乘逆矩阵, γ Spectrum Relative Comparison Inverse Matrix Least Squares Inverse Matrix 摘要

目的：NaI(Tl)谱仪相对分辨率较低，通过比较3种解谱方法并选出最优解谱方法，提高γ能谱测量值的准确性和可靠性。方法：先用相对比较、逆矩阵、最小二乘法3种分析方法对已知活度的混合标准源进行活度计算，再比较计算活度与标准源活度的相对偏差。结果：使用相对比较法、逆矩阵法、最小二乘法计算出3种核素相对偏差分别为：226Ra (13.62%、14.71%、9.12%)、232Th (14.43%、20.59%、13.17%)、40K (2.18%、4.91%、4.84%)。结论：使用最小二乘法对226Ra和232Th 2种核素的计算活度相对偏差最小，优于相对比较法与逆矩阵法；40K核素的活度计算采用相对比较法更合适。

 关键词

γ能谱，相对比较，逆矩阵，最小二乘逆矩阵

 Comparative Study of Three γ Spectrum Solution Methods Based on NaI(Tl) Spectrometer<sup> </sup>

Youhua Hu*, Jianchao Wang#

Radiation Environmental Supervising Station of Xinjiang, Urumqi Xinjiang

Received: Sep. 19th, 2022; accepted: Oct. 18th, 2022; published: Oct. 25th, 2022

 ABSTRACT

Objective: The resolution of NaI(Tl) spectrometer is relatively low. By comparing the three spectrum solution methods and selecting the best solution method, the accuracy and reliability of the measured values of the γ energy spectrum can be improved. Methods: First use relative comparison, inverse matrix, and least square analysis methods to calculate the activity of the mixed standard source with known activity, and then compare the relative deviation between the calculated activity and the standard source activity. Results: The relative deviations of the three nuclides calculated using the relative comparison method, the inverse matrix method, and the least square method: 226Ra (13.62%, 14.71%, 9.12%), 232Th (14.43%, 20.59%, 13.17%), 40K (2.18%, 4.91%, 4.84%). Conclusions: Using the least square method to calculate the relative deviation of the activities of 226Ra and 232Th nuclides is the smallest, which is better than the relative comparison method and the inverse matrix method. The activity calculation of 40K nuclides is more suitable for the relative comparison method.

Keywords:γ Spectrum, Relative Comparison, Inverse Matrix, Least Squares Inverse Matrix

Copyright © 2022 by author(s) and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

 1. 引言

随着核能的广泛实践与不断发展，环境介质中存在的放射性核素也愈发增多，对自然环境与公众健康产生潜在的威胁。作为在辐射环境测量工作中最常见的仪器，γ能谱仪具有灵敏度高、快速便捷等特点，可对环境中的放射性核素进行精准的定量、定性分析 [ 1 ]。但考虑到环境放射性核素种类繁多且活度浓度往往较低，在实际测量过程中得到的γ能谱中包含了大量的核素信息，其谱形较为复杂。尤其是在使用碘化钠[NaI(Tl)] γ谱仪进行放射性核素活度浓度测量时，由于不同类型γ能谱仪灵敏区材料不同，与高纯锗(HPGe) γ谱仪相比，NaI(Tl) γ谱仪的能量分辨率相对较低，入射能量相近的γ射线容易在能谱上出现能峰重叠的情形。重叠能峰的剥离增加了后续能谱解析的难度，不利于放射性核素测量工作的开展。因此，使用合适的解谱方法进行分析计算以确保γ能谱测量分析结果的准确可靠就显得尤为重要。

本文利用NaI(Tl)闪烁体γ谱仪 [ 2 ] 对多个已知活度的标准放射源进行测量，使用相对比较法、逆矩阵法与最小二乘法分别对γ能谱进行解谱分析，比较3种不同解谱方法的活度计算结果，探讨3种解谱方法的适用性，旨在提高辐射环境监测工作的效率，增加γ能谱在环境介质放射性核素活度浓度测量工作中的可靠性与准确性。

 2. 材料与方法 2.1. γ能谱仪选择

本文所使用的探测系统是由北京核仪器厂设计制造的低本底环境γ能谱。该套γ能谱探测系统主要由BH1324F型探测器、BH1324型一体化多道分析器、NG401-261型铅室、Ums 38 γ能谱采集软件、谱分析软件等部分组成。能谱仪探头是尺寸为 φ 75 mm × H 75 mm的NaI(Tl) 低钾闪烁晶体，对于137Cs标准源在661.65 keV处的特征峰能量分辨率小于8.5%。

 2.2. 分析方法

首先选用由中国计量科学研究院校准的60Co与137Cs标准源对能谱仪进行能量刻度，其活度分别为3.82 μCi与2.31 μCi。其中，60Co选用能量为1173.20 keV与1332.50 keV的特征峰，而137Cs选用能量为661.65 keV的特征峰 [ 3 ]。其次，将中国计量科学研究院提供并校准的226Ra、232Th、40K单核素标准体源用于后续解谱分析 [ 4 ]，其校准核素活度分别为286.30 Bq、184.00 Bq、1053.90 Bq。接下来将中国计量科学研究院提供的混合标准源视为被测样品，混合标准源中包含放射性核素226Ra (74.50 Bq)、232Th (33.60 Bq)与40K (569.80 Bq)。

在对标准源进行测量前，先对完成能量刻度的γ能谱仪进行时长为43,200 s (12 h)的环境本底测量，以排除环境本底干扰。本底测量结束后，对226Ra、232Th、40K单核素标准源与混合标准源进行时长为21 600 s (6 h)的测量，使用谱分析软件对得到的能谱进行谱数据分析。

2.2.1. 相对比较法

相对比较法适用于现有且可用的标准源中已包含全部目标放射性核素的条件下对目标放射性核素的活度浓度进行分析，主要计算分析过程如下：

C j i = A j N e t j i (1)

其中， C j i 为标准源中第j种放射性核素的第i个特征峰的系数； A j 为标准源中第j种放射性核素的活度(Bq)； N e t j i 为标准源中第j种放射性核素的第i个特征峰的净计数率(计数/每秒)。

Q j = C j i ( A j i − A j i b ) W D j (2)

式中，混合标准源中第j种放射性核素的活度浓度为 Q j (Bq/kg)； A j i 为混合标准源中第j种放射性核素的第i个特征峰的净计数率(计数/每秒)； A j i b 为对应 A j i 的特征峰本底计数率(计数/每秒)；W为被测标准源的重量(kg)； D j 为放射性核素衰变修正系数，但考虑到在本研究中所测量的放射性核素半衰期均较长(226Ra的半衰期为1600 a，232Th的半衰期为1.40 × 1010 a，40K的半衰期为1.27 × 109 a)，所以此项可以被忽略 [ 5 ]。

 2.2.2. 逆矩阵法

逆矩阵法又被称作解线性方程组法，主要适用于被测样品中的放射性核素成分已知，但测量得到的γ能谱出现了能峰重叠的情况下进行被测样品放射性核素活度浓度分析。在扣除环境本底的贡献后，混合源的γ能谱在某一特征峰能道区内的计数率由该特征峰本身所对应的放射性核素、特征峰能量相近的放射性核素以及γ能量较高的放射性核素的康普顿散射等的叠加贡献 [ 6 ]。计算方法如下：

c i = ∑ j = 1 m a i j x j (3)

式中：j为混合源中第j种核素，m为混合源中全部的核素种类，i为混合源γ能谱在第i个特征峰道区内的计数率(计数/秒)， x j 为混合样品中第j种核素未知活度， a i j 为第i个特征峰道区对第j种核素的响应系数，如公式4)所示 [ 7 ]。

a i j = N e t j i A j (4)

式中， N e t j i 为第j种核素的标准谱在第i特征峰道区内的计数率； A j 为第j种核素的活度 [ 7 ]。

其矩阵方程为：

A m × m ⋅ X m × m = C m × m (5)

A m × m = ( X m × m − 1 ⋅ C m × m ) T (6)

式中： A m × m 为全部响应系数 a i j 组合而成的响应矩阵。

混合源中第j种核素的活度 x j 的计算如下：

X m × 1 = A m × m − 1 ⋅ C m × 1 (7)

x j = ∑ j = 1 m a j i − 1 c i (8)

 2.2.3. 最小二乘逆矩阵法

最小二乘逆矩阵法是逆矩阵法的延伸 [ 8 ]。逆矩阵法仅利用了每种放射性核素的一个特征峰道区的计数率对混合源中对应的核素活度进行计算，容易受到道址漂移的影响，并且未对γ能谱数据信息进行充分利用。而结合了最小二乘法统计学原理的最小二乘逆矩阵法对核素所具有的各个能量的γ射线都给予考虑，通过选择核素的多个主要的γ射线的能峰参与计算，即被选择的能峰道区数量n大于混合样品中的核素种类m。

C i = ∑ j = 1 m a i j X j ,     i = 1 , 2 , ⋯ , n (9)

根据最小二乘法的理论，其矩阵方程为：

X m × m ⋅ A m × n T = C m × n (10)

A n × m = ( X m × m − 1 ⋅ C m × n ) T (11)

混合样品中核素的活度的计算公式为：

X m × 1 = ( A T ⋅ A ) − 1 ⋅ A T ⋅ C (12)

 3. 结果 3.1. 3种方法计算结果 3.1.1. 相对比较法计算核素活度

在利用相对比较法对混合源中放射性核素活度进行计算时，先选择合适的核素特征峰，后在对应核素γ能谱中定位所选特征峰的道区，并计算出特征道区内的净计数率 [ 9 ]。特征峰与道区选取见表1。

Nuclide characteristic peaks and corresponding channel selected by relative comparison metho
放射性核素 226Ra 232Th 40K
特征峰能量/keV 609.32 911.32 1461.22
特征峰道址 305 447 730

表1. 相对比较法选取的核素特征峰与其对应的道址

放射性核素226Ra、232Th、40K在混合标准源中的活度通过公式(1)与(2)计算得到，结果见表2 [ 10 ]。混合标准源中的放射性核素226Ra、232Th、40K的活度分别为64.35 Bq、38.45 Bq、557.36 Bq。

Nuclide activity calculated by relative comparison metho
放射性核素 226Ra 232Th 40K
计算活度/Bq 64.35 38.45 557.36

表2. 相对比较法计算的核素活度

 3.1.2. 逆矩阵法计算核素活度

在使用逆矩阵法对放射性核素活度进行计算时 [ 11 ]，每种放射性核素选取1个特征峰能量，在226Ra、232Th、40K标准源谱、混合源谱与环境本底谱中标定3种核素对应特征峰的道区与计数率。在本研究中，逆矩阵法所选取的核素特征峰与相对比较法相同，见表1。由3个特征峰道区在全部能谱中的净计数率结合公式(3)~(6)计算得到的响应矩阵 A 3 × 3 为：

A = [ 0.004272 0.003597 0.0000916 0.0008 0.002415 0.000101 0.0005 0.000281 0.000295 ]

通过公式(7)、(8)的计算，得到了混合标准源中放射性核素226Ra、232Th、40K的活度，其计算活度值见表3。

Nuclide activity calculated by inverse matrix metho
放射性核素 226Ra 232Th 40K
计算活度/Bq 85.45 26.68 541.80

表3. 逆矩阵法计算的核素活度

 3.1.3. 最小二乘法计算核素活度

为了更加充分的利用γ能谱中的数据，在利用最小二乘法进行活度计算时，对3种放射性核素选取4个不同特征峰，特征峰能量及其道区见表4。

Nuclide characteristic peaks and their corresponding channel elected by least square metho
放射性核素 226Ra 232Th 40K
特征峰能量/keV 351.99 609.32 911.32 1461.22
特征峰道址 187 305 447 730

表4. 最小二乘法选取的核素特征峰与其对应的道址

通对226Ra、232Th、40K标准源、混合标准源与环境本底进行测量后，其相应的响应矩阵 A 4 × 3 为：

A = [ 0.005175 0.00403 0.000141 0.004272 0.003597 0.0000916 0.0008 0.002415 0.000101 0.0005 0.000281 0.000295 ]

通过公式(12)计算得到了混合标准源中放射性核素的活度，见表5。

Nuclide activity calculated by least square metho
放射性核素 226Ra 232Th 40K
计算活度/Bq 81.29 29.17 542.24

表5. 最小二乘法计算的核素活度

 3.2. 3种方法结果比较

使用相对比较法、逆矩阵法、最小二乘法3种解谱方法对混合标准源中的3种放射性核素226Ra、232Th、40K的活度计算结果见表6 [ 12 ]。根据表6计算得到3种解谱方法的放射性核素活度与活度标定值之间的相对偏差见表7。

Calculation of nuclide activity by different spectral analysis method
放射性核素 活度标定值/Bq 3种解谱方法计算活度(Bq)
相对比较法 逆矩阵法 最小二乘法
226Ra 74.5 64.35 85.45 81.29
232Th 33.6 38.45 26.68 29.17
40K 569.8 557.36 541.80 542.24

表6. 不同解谱方法的计算核素活度

Calculate the relative deviation between nuclide activity and calibration activit
放射性核素 3种计算方法的相对偏差(%)
相对比较法 逆矩阵法 最小二乘法
226Ra 13.62 14.71 9.12
232Th 14.43 20.59 13.17
40K 2.18 4.91 4.84

表7. 计算核素活度与标定活度之间的相对偏差

 4. 讨论

本研究结果显示：相对偏差最大的是逆矩阵法对232Th放射性核素活度的计算为20.59%；相对偏差最小的是相对比较法对40K放射性核素活度计算为2.18%；针对226Ra、232Th的计算活度相对偏差最小的是最小二乘法分别为9.12%、13.17%；3种解谱方法对放射性核素40K的计算活度相对偏差均低于5%；其中3种解谱方法对于核素232Th的活度计算结果相对偏差均较大，这可能是由于232Th所选择的特征峰能量主要集中在低能段，γ能谱低能段存在康普顿散射、韧致辐射等额外计数的干扰，且NaI(Tl)谱仪的能量分辨率较低所致。

综上，在计算232Th放射性核素活度时不适合选取逆矩阵法，其相对偏差已经超出一般质控偏差要求的15%；在计算混合源中226Ra和232Th的放射性活度时最小二乘法为最优解谱方法；在计算混合源中40K的放射性活度时相对比较法为最优解谱方法，但是相对比较法使用的前提是标准物质核素必须含有40K，因此在无标准物质不含有40K时，也可选用其他2种解谱方法。本文给出了三种常见核素在使用NaI(Tl)谱仪进行核素分析时，能针对不同的目标核素更快地选出使分析结果更加准确和可靠的最优解谱方法。

 文章引用

胡有华,王建超. 基于NaI(Tl)谱仪的3种γ能谱解谱方法比较Comparative Study of Three γ Spectrum Solution Methods Based on NaI(Tl) Spectrometer[J]. 环境保护前沿, 2022, 12(05): 1049-1055. https://doi.org/10.12677/AEP.2022.125130

 参考文献 References 杜金秋, 高会, 王震, 等. 土壤样品中放射性核素活度γ谱分析比对[J]. 分析实验室, 2016, 35(9): 1041-1045.
https://doi.org/10.13595/j.cnki.issn1000-0720.2016.0237 方登富, 韦应靖, 崔伟, 等. 基于NaI(Tl) γ能谱温漂修正方法研究[J]. 中国辐射卫生, 2021, 30(1): 19-23+33.
https://doi.org/10.13491/j.issn.1004-714X.2021.01.005 杨秀玉, 李庆光, 葛琴. 特征峰干扰判据K值法解析γ能谱[J]. 核电子学与探测技术, 2021, 41(1): 87-91. 李奇, 王世联, 樊元庆, 等. 环境比对样品中226Ra, 232Th和40K γ谱测量分析[J]. 核电子学与探测技术, 2011, 31(11): 1234-1236+1283. 卿云花, 刘盼, 郑琪珊, 等. 放射性核素γ谱分析方法探讨[J]. 核技术, 2017, 40(9): 57-62. 杨奎. 碘化钠能谱仪在测量核事故放射性物质活度中的应用[D]: [硕士学位论文]. 南昌: 东华理工大学, 2015. 赖忠俊, 宣志强, 俞顺飞, 等. 有源-无源γ能谱分析方法在分析环境样品中放射性核素的应用[J]. 中国辐射卫生, 2020, 29(4): 378-381+386.
https://doi.org/10.13491/j.issn.1004-714X.2020.04.014 房宗良, 曹剑锋, 王季红. 最小二乘逆矩阵法在建材γ谱分析中的应用[J]. 核电子学与探测技术, 2006, 26(6): 717-720. 张京, 周强, 拓飞, 等. 用标准样品刻度γ能谱仪效率值的变化比较[J]. 中国辐射卫生, 2016, 25(3): 361-363.
https://doi.org/10.13491/j.cnki.issn.1004-714x.2016.03.037 姚杰, 钱爱君, 高林峰. 2018年参加全国放射性核素γ能谱分析比对结果[J]. 中国辐射卫生, 2019, 28(3): 321-324.
https://doi.org/10.13491/j.issn.1004-714x.2019.03.028 王硕, 杨宝路, 周强, 等. 无源效率刻度在马林杯状气体源活度测量中的应用[J]. 中国辐射卫生, 2022, 31(1): 6-12.
https://doi.org/10.13491/j.issn.1004-714X.2022.01.002 许家昂, 宋钢, 苏静, 等. 某氯化稀土γ能谱测量结果分析[J]. 中国辐射卫生, 2021, 30(6): 693-696+700.
https://doi.org/10.13491/j.issn.1004-714X.2021.06.007
Baidu
map