高纯锗γ能谱仪在放射性测量中具有分辨率强、灵敏度高、刻度方便、样品前处理简单等优点,广泛应用于日常和核事故应急状况下环境样品中放射性核素分析[1-2]。实验室间的比对可以提高γ能谱实验室放射性核素检测能力和水平,保证放射性核素γ能谱分析的质量[3-4]。本实验室常年承担北京市大气、水、土壤以及食品中放射性核素监测工作,为加强实验室质量控制,提高放射性核素实验室检测能力,本实验室于2014—2021年连续参加了由中国疾病预防控制中心组织的全国放射性核素γ能谱分析能力考核,对历年结果分析如下[5-6]。
1 材料与方法 1.1 样品来源比对组织机构为中国疾病预防控制中心,组织机构将均匀性检验合格、完成定值的样品分份,每份样品约500 g土壤样品粉末,装入相应的容器内密封后分组邮寄给各参加机构。γ能谱分析能力考核中使用的标准源为中国计量科学研究院提供的环境放射性标准物质,分别为河泥标准源(HN08095)和土壤监测效率校准源(7NTR/70-1201),具体情况见表1。
本实验室使用的仪器设备为ORTEC公司生产的GEM-MX7080P4型高纯锗 γ 谱仪测量系统,其对3″ × 3″ NaI (Tl) 探测器的相对探测效率为71%。谱仪对60Co的1332.5 keV γ 射线的能量分辨率约为2.2 keV,峰康比为88∶1,积分本底为3.01 s−1(20~2 400 keV)。配套使用 DSPEC jr2.0数字化多道分析器,谱测量分析软件采用Gamma Vision 32。谱仪测量时探测器置于壁厚118 mm、内腔278 mm (Φ) × 411 mm (H) 的复合屏蔽铅室内。谱仪在样品测量前,使用中国计量科学研究院提供的标准源对系统进行了能量刻度和效率刻度。
1.3 方法依据《土壤中放射性核素的γ能谱分析方法》GB/T 11743—2013[7]和《高纯锗γ能谱分析通用方法》GB/T 11713—2015[8]对样品进行制备、测量、计算及不确定度评定。
1.3.1 样品制备与测量收到考核样品后,将样品平铺在搪瓷盘中并放入电子恒温烤箱,温度设置为105℃,烘烤4 h后,将样品置于干燥器中冷却至室温。将样品封装于75 mm( Φ ) × 70 mm(H)样品盒中,称量土壤样品装样质量。按照本底、标准源、考核样品和本底的顺序放置在高纯锗γ能谱仪的探头上进行测量,测量时间均为72 h,所有样品测量条件保持一致。
1.3.2 计算分析方法2014—2021年考核的目标核素涉及的208Tl、228Ac、232Th、238U、40K和137Cs均包括在标准源中,其中232Th的放射性比活度可依据平衡条件下子体208Tl的特征能量583.19 keV或228Ac的特征能量911.07 keV来计算,238U的放射性活度浓度可依据平衡条件下子体234Th的特征能量63.29 keV来计算,40K 和 137Cs分别根据其特征能量1 460.75 keV和661.66 keV来计算。因此考核样品采用目标核素特征能量全能峰计数的相对比较法进行计算,目标核素比活度计算公式[9-10] :
$ {A}_{\mathrm{i}\mathrm{s}}=\frac{{N}_{\mathrm{i}\mathrm{s}}/{t}_{\mathrm{s}}-{N}_{\mathrm{i}\mathrm{b}}/{t}_{\mathrm{b}}}{{(N}_{\mathrm{i}\mathrm{c}}/{t}_{\mathrm{c}}-{N}_{\mathrm{i}\mathrm{b}}/{t}_{\mathrm{b}}){m}_{\mathrm{s}}{F}_{1}{F}_{2}}{A}_{\mathrm{i}\mathrm{c}} $ | (1) |
式中Ais 为考核样品第i种核素的比活度,Bq/kg;Aic是标准源第i种核素的活度,Bq;Nis , Nib和 Nic分别是考核样品、本底样品和标准源第i种核素的特征能量全能峰的净面积计数,counts;ts, tb 和 tc 分别是考核样品、本底样品和标准源测量活时间s;ms 为被测样品净重,kg;F1 自吸收校准因子;F2 为核素i校正到组织者规定日期的衰变校正系数。测量结果的相对不确定度采用公式:
$ {U}_{\mathrm{t}\mathrm{o}\mathrm{t}\mathrm{a}\mathrm{l}}=k\sqrt{{U}_{1}^{2} + {U}_{2}^{2} + {U}_{3}^{2} + {U}_{4}^{2}} $ | (2) |
式中,Utotal为合成相对不确定度;k为包含因子,取2;U1为计数率的相对不确定度;U2 标准源的相对不确定度;U3 样品称量相对不确定度;U4自吸收校正引入的相对不确定度。
1.4 比对结果判定方法参考值由组织方确定,以相对偏差RD、Z值,U值以及相对合成不确定度Urel 4个参数通过准确度和精确度来判定参加机构的考核结果,公式如下:
$ RD=\frac{{A}_{\mathrm{i}\mathrm{s}}-{A}_{\mathrm{i}\mathrm{r}}}{{A}_{\mathrm{i}\mathrm{r}}}\times 100\% $ | (3) |
$ Z=\frac{{A}_{\mathrm{i}\mathrm{s}}-{A}_{\mathrm{i}\mathrm{r}}}{\sigma } $ | (4) |
$ U=\frac{|{A}_{\mathrm{i}\mathrm{s}}-{A}_{\mathrm{i}\mathrm{r}}|}{2.58\sqrt{{U}_{\mathrm{i}\mathrm{s}}^{2} + {U}_{\mathrm{i}\mathrm{r}}^{2}}} $ | (5) |
$ {U}_{\mathrm{r}\mathrm{e}\mathrm{l}}=\sqrt{\Bigg({\frac{{U}_{\mathrm{i}\mathrm{s}}}{{A}_{\mathrm{i}\mathrm{s}}}\Bigg)}^{2} + \Bigg({\frac{{U}_{\mathrm{i}\mathrm{r}}}{{A}_{\mathrm{i}\mathrm{r}}}\Bigg)}^{2}}\times 100\% $ | (6) |
式中,RD为测量值与参考值的相对偏差,%;Ais为考核样品中核素i的测量值比活度,Bq/kg;Air为考核样品中核素i的参考值比活度,Bq/kg;σ为标准偏差的目标值,基于合理性原则,取值为0.20 × Air,Bq/kg;Uis为核素i的测量值的扩展不确定度(k = 2),Bq/kg;Uir为核素i的参考值的扩展不确定度(k = 2),Bq/kg。
1.5 判定标准依据以上4个参数的计算以及准确度和精确度的判定后,结合判定结果与提交结果报告资料的完整性打分, 其打分规则如下:
① 当比对样品中所有核素U ≤ 1且Urel ≤ 20%时考核结果为合格,打分60分。
② 当比对样品中所有核素U ≤ 0.4且Urel ≤ 15%时,加10分,打分70分。
③同时满足1、2条件后,对检测报告进行评分(分值10分),得80分方可参与优秀评比。
④再进行Q值评定,Q值为20分,总分大于等于95分为优秀。
2 结 果2014—2021年能力考核样品均为土壤样品,共涉及放射性核素6种,分别是天然放射性核素208Tl、228Ac、232Th、40K、238U和人工放射性核素137Cs。共检测上述放射性核素28次,分别检测208Tl、228Ac、232Th、40K、238U和137Cs为4、4、4、8、1和7次。2014—2021年放射性核素γ能谱分析能力考核结果见表2。RD、Z、U和Urel的分布范围分别为−11.82%~5.97%,−0.59~0.30,0.02~0.92和4.33%~10.49%,均值 ± 标准偏差分别为−2.50% ± 5.02%,−0.13 ± 0.25,0.28 ± 0.20和7.20% ± 1.50%。按照U ≤ 1且Urel ≤ 20%时考核结果为合格的标准,28次放射性核素测量结果均合格。按照U ≤ 0.4且Urel ≤ 15%时考核结果满足优秀的必要条件之一,28次放射性核素测量结果Urel均小于15%,U值 ≤ 0.4占比为79%(22/28),U值 > 0.4有6次,分别是2014、2015和2016年 40K(U = 0.46、0.92和0.63)3次,2016和2018年137Cs(U = 0.44、0.44)2次和2014年208Tl(U = 0.53)1次。2014—2021年考核结果均合格,其中2017年、2020年和2021年考核结果优秀,总体合格率为100%(8/8),总体优秀率为39%(3/8)。如果针对每次考核中单个核素而言,单个核素测量结果按照U ≤ 1且Urel ≤ 20%时为合格,U ≤ 0.4且Urel ≤ 15%时为优秀,核素测量结果合格率为100%(28/28),优秀率为75%(18/24),其中208Tl、228Ac和232Th测量结果的优秀率为92%(11/12)、40K测量结果优秀率为63%(5/8)、238U测量结果优秀率为100%(1/1)和137Cs测量结果优秀率为71%(5/7)。从图1和图2可以看出,随着逐年参加比对,U值逐渐降低至满足U ≤ 0.4考核优秀的条件,Urel虽有波动,但均满足Urel ≤ 15%考核优秀的条件。
本实验室在2014—2021年全国放射性核素γ能谱分析检测能力考核中合格率为100%,通过准确度(|RD|、|Z|和U)和精确度(Urel)来判定测量值与参考值的差异程度,|RD|均<20%,|Z|和U均<1,说明在99%的给定概率水平下测量值与参考值不存在显著差异,Urel均<15%表明测量值有优秀的精确度。说明实验室质量控制体系完善,相关技术人员熟练掌握了实验室γ能谱分析方法,能够准确测定样品中的放射性核素比活度。
比对考核中涉及天然放射性核素的分析,按照本底、标准源、考核样品和本底的顺序进行测量,在考核样品测量前后分别测量空样品盒作为本底样品,取均值作为本底结果,减少了实验室本底波动对测量结果的干扰。由于标准源中的核素包括历年考核中所涉及的核素,因此在测量中使用的方法是相对比较法,不考虑样品的自吸收和符合相加,此方法较效率曲线法误差更小。定期用标准源对高纯锗γ能谱仪进行刻度,保证峰位确定的准确性。以上因素均是取得良好考核结果的原因。
2014—2021年28次放射性核素测量结果Urel在5.0%~10.5%之间波动,但是由于U值6次 > 0.4而不满足优秀评定条件,其中 40K 3次,相对偏差RD分别为5.13%、−11.82%和−9.83,137Cs 2次,相对偏差RD分别为−7.79%和−8.06,208Tl 1次,相对偏差RD为−7.09%。由于在核素比活度计算中使用的是相对比较法,其基本假设样品源和标准源有相同的探测效率,这就要求测量样品与标准源有相同的体积、密度、形状、基质组成以及测量时与探测器相对位置相同。这在实践中是不可能完全匹配的,特别是样品与标准源基质成分相对含量会存在差别,样品与标准源密度不同引起不同的自吸收效应导致U值偏大。本实验室2019年前在考核中使用的标准源的表观密度(装样质量/样品盒体积)为0.97 g/cm3,而在考核中的土壤样品的表观密度均 > 0.97 g/cm 3,当样品密度大于标准源密度时,在使用相对比较法时,由于自吸收效应的不同导致测量值小于参考值,因此在6次U值 > 0.4的情况下,相对负偏差出现5次。
从2019年开始本实验室使用新购入的土壤监测效率校准源,其表观密度为 1.38 g/cm3,同时在装样过程中使用定制液压器制样,使样品和标准源的表观密度一致,减少自吸收效应对测量结果的影响,因此2019—2021年U值均≤0.4。如果样品与标准源密度不同应进行相应的自吸收校正,可以通过实验法或基于蒙特卡罗模拟的方法进行效率修正[11-13]。影响相对比较法测量准确性的另一个关键因素就是标准源的准确性,因此必须保证标准源定期经法定部门检定或校准。
在进行环境样品测量时,各种天然放射性核素γ射线会对考核样品中的待分析核素造成干扰,特别是对于单能放射性核素如40K和137Cs。232Th的子体228Ac的1459.2 keV γ射线会对40K的1460.8 keV γ射线产生干扰,226Ra的子体214Bi的665.45 keV γ射线会对137Cs的661.6 keV γ射线产生干扰。而环境样品尤其是土壤样品中均含有232Th和226Ra,尤其是在232Th和226Ra含量特别大时,由于射线能量非常接近,通过解谱软件进行重叠峰分离一般较难实现。在具体的分析中可以利用效率曲线法或混合体标准源实验法消除γ能谱中232Th对40K的干扰以及226Ra对137Cs的干扰[14]。
要在放射性核素γ能谱分析能力考核中取得良好的结果,也需要加强实验室γ能谱分析质量控制。检测报告和原始记录信息完整,名词术语准确表述、符号、单位和有效位数规范使用,原始记录还需公式表述正确,计算过程完整;实验设备进行定期维护保养、期间核查和检定或校准;相关技术人员参加国家组织的各种培训考核以及仪器厂家的技术培训;实验设备定期的能量刻度、效率刻度和本底测量等。
[1] |
孟庆华, 马永忠, 王欢, 等. 2009—2016年高纯锗γ能谱分析方法比对结果分析[J]. 中国辐射卫生, 2018, 27(2): 145-147,152. Meng QH, Ma YZ, Wang H, et al. Analysis of the intercomparison results on radioactivity measurements by HPGe γ spectrometer system during 2009-2016[J]. Chin J Radiol Health, 2018, 27(2): 145-147,152. DOI:10.13491/j.issn.1004-714X.2018.02.014 |
[2] |
Díaz-Asencio M, Sanchez-Cabeza JA, Ruiz-Fernández AC, et al. Calibration and use of well-type germanium detectors for low-level gamma-ray spectrometry of sediments using a semi-empirical method[J]. J Environ Radioact, 2020, 225: 106385. DOI:10.1016/j.jenvrad.2020.106385 |
[3] |
Glavič-Cindro D. Comparison of intercomparison results of gamma ray spectrometry of spiked and real samples[J]. Appl Radiat Isot, 2018, 134: 59-63. DOI:10.1016/j.apradiso.2017.07.009 |
[4] |
Mikami S, Ishikawa D, Matsuda H, et al. Guidance for in situ gamma spectrometry intercomparison based on the information obtained through five intercomparisons during the Fukushima mapping project
[J]. J Environ Radioact, 2019, 210: 105938. DOI:10.1016/j.jenvrad.2019.03.014 |
[5] |
王欢, 孟庆华, 马永忠, 等. 2017—2020年北京地区大气气溶胶中7Be和 210Pb放射性活度浓度监测与分析
[J]. 中华放射医学与防护杂志, 2021, 41(9): 690-694. Wang H, Meng QH, Ma YZ, et al. Monitoring and analysis of activity concentrations of 7Be and210Pb in atmospheric aerosol in Beijing from 2017 to 2020 [J]. Chin J Radiol Med Prot, 2021, 41(9): 690-694. DOI:10.3760/cma.j.issn.0254-5098.2021.09.009 |
[6] |
孟庆华, 马永忠, 王欢, 等. 2017-2018年北京地区土壤的放射性水平监测与分析[J]. 中华放射医学与防护杂志, 2020, 40(9): 702-706. Meng QH, Ma YZ, Wang H, et al. Monitoring and analysis of radioactivity in soils in Beijing during 2017-2018[J]. Chin J Radiol Med Prot, 2020, 40(9): 702-706. DOI:10.3760/cma.j.issn.0254-5098.2020.09.009 |
[7] |
国家质量监督检验检疫总局. GB/T 11743—2013 土壤中放射性核素的γ能谱分析方法[S]. 北京: 中国标准出版社, 2014. General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China. GB/T 11743—2013 Determination of radionuclides in soil by gamma spectrometry[S]. Beijing: Standards Press of China, 2014. |
[8] |
中华人民共和国国家质量监督检验检疫总局. GB/T 11713—2015 高纯锗γ能谱分析通用方法[S]. 北京: 中国标准出版社, 2016. General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China. GB/T 11713—2015 General analytical methods of high-purity germanium gamma spectrometer[S]. Beijing: Standards Press of China, 2016. |
[9] |
Wang S, Yang BL, Zhou Q, et al. Radionuclide content and risk analysis of edible mushrooms in northeast China[J]. Radiat Med Prot, 2021, 2(4): 165-170. DOI:10.1016/j.radmp.2021.10.001 |
[10] |
Tuo F, Xu CH, Zhang Q, et al. A review of nationwide radioactivity comparisons on gamma-ray spectrometry organized by the NIRP, China[J]. Appl Radiat Isot, 2014, 87: 435-438. DOI:10.1016/j.apradiso.2013.11.072 |
[11] |
李明生, 姜庆寰, 程金生. 应用Geant4模拟EBT3胶片的能量响应和质子射程测量[J]. 中国辐射卫生, 2022, 31(6): 669-673. Li MS, Jiang QH, Cheng JS. Energy response and proton range measurement with EBT3 film simulated by Geant4[J]. Chin J Radiol Health, 2022, 31(6): 669-673. DOI:10.13491/j.issn.1004-714X.2022.06.005 |
[12] |
赖忠俊, 宣志强, 俞顺飞, 等. 有源-无源γ能谱分析方法在分析环境样品中放射性核素的应用[J]. 中国辐射卫生, 2020, 29(4): 378-381,386. Lai ZJ, Xuan ZQ, Yu SF, et al. Application of sources and sources-less methods of γ spectrometry for analysis of radionuclides in environmental sample[J]. Chin J Radiol Health, 2020, 29(4): 378-381,386. DOI:10.13491/j.issn.1004-714X.2020.04.014 |
[13] |
孔淑颖, 杨宝路, 王仁生, 等. 基于Geant4模拟37Ar活度测量计数损失的修正方法
[J]. 中国辐射卫生, 2023, 32(2): 93-97,101. Kong SY, Yang BL, Wang RS, et al. Counting loss correction in 37Ar activity measurement based on Geant4 simulation [J]. Chin J Radiol Health, 2023, 32(2): 93-97,101. DOI:10.13491/j.issn.1004-714X.2023.02.001 |
[14] |
王崇杰, 张爱莲. 消除γ能谱中232Th对40K的干扰
[J]. 核技术, 2003, 26(11): 831-833. Wang CJ, Zhang AL. Eliminating interference of 232Th to 40K in γ energy spectra [J]. Nucl Tech, 2003, 26(11): 831-833. DOI:10.3321/j.issn:0253-3219.2003.11.006 |