文章快速检索    
  核技术  2018, Vol. 41 Issue (1): 010501   DOI: 10.11889/j.0253-3219.2018.hjs.41.010501
0

引用本文 [复制中英文]

刘文彪, 陈伟, 田自宁, 王雪梅, 苏川英, 冯天成. γ全能峰右高左低现象理论分析[J]. 核技术, 2018, 41(1): 010501. DOI: 10.11889/j.0253-3219.2018.hjs.41.010501. [复制中文]
LIU Wenbiao, CHEN Wei, TIAN Zining, WANG Xuemei, SU Chuanying, FENG Tiancheng. Theoretical study of the phenomenon that right is higher than left in full-energy peak of the γ spectrum[J]. Nuclear Techniques, 2018, 41(1): 010501. DOI: 10.11889/j.0253-3219.2018.hjs.41.010501.
[复制英文]

基金项目

国家自然科学基金(No.11405134)资助

第一作者

刘文彪, 男, 1990年出生, 2015年于兰州大学获硕士学位, 研究领域为核技术及其应用

文章历史

收稿日期: 2017-08-07
修回日期: 2017-10-24
γ全能峰右高左低现象理论分析
刘文彪, 陈伟, 田自宁, 王雪梅, 苏川英, 冯天成     
西北核技术研究所 西安 710024
摘要: 用HPGe γ谱仪对高活度放射性样品进行测量时,在能谱中出现了全能峰右侧明显高于左侧的现象,为准确定量分析γ放射性核素带来困难。通过理论分析,详细解释了该现象产生的原因,实验验证了理论分析结论是正确的。分析结果表明,全能峰右侧高于左侧的现象主要由偶然符合效应引起,并且总体上全能峰左右两侧高度的比值随死时间增大先增后减,当死时间小于等于1.0%时,全能峰左右两侧高度近似相等,当死时间位于1.0%-4.0%时,峰左右高度比值随死时间而增大,当死时间位于4.0%-56%时,峰左右高度比值随死时间增加而减小。并通过实验结果与LabSOCS计算结果对比,显示右高左低这种现象并不会对点源的能谱测量造成影响。
关键词: γ能谱    全能峰    偶然符合    
Theoretical study of the phenomenon that right is higher than left in full-energy peak of the γ spectrum
LIU Wenbiao , CHEN Wei , TIAN Zining , WANG Xuemei , SU Chuanying , FENG Tiancheng     
Northwest Institute of Nuclear Technology, Xi'an 710024, China
Received date: 2017-08-07; accepted date: 2017-10-24
Supported by National Natural Science Foundation of China (No.11405134)
First author: LIU Wenbiao, male, born in 1990, graduated from Lanzhou University with a master's degree in 2015, focusing on nuclear techniques and application
Abstract: Background: The phenomenon that right is higher than left in full-energy peak of the γ spectrum in the high active sample using HPGe gamma-ray spectrometry in laboratory, brings error into the gamma-radioactive nuclide analysis. Purpose: This paper aims to explain the phenomenon to realize accurate quantitative analysis. Methods: The theoretical formula of the right-to-left ratios was developed, and experiment was designed to prove it. Results: The experiment indicates that the theoretical calculation is right. The research achievement shows that the phenomenon is caused by accidental coincidence. Conclusion: The right-to-left ratios increase first, then decrease with the higher dead time. When the dead time is lower than 1.0%, the right is approximate equal to the left in full energy peak. With the dead time is higher than 1.0% and lower than 4.0%, the ratios increases with the dead time. As the dead time is higher than 4.0% and lower than 56%, the ratios decreases with the dead time. Comparing the experiment with LabSOCS, it indicated that the phenomenon has no influence on analysis of the point source.
Key Words: γ spectrometry    Full-energy peak    Accidental coincidence    

随着核技术的应用,环境放射性日益引起人们重视。采样后进行实验室测量分析是环境放射性调查的一种有效手段,γ能谱测量分析是关键技术之一。通过γ能谱的测量分析,可得到样品中发射γ射线的核素的种类和活度,相关研究文献报道较多[1-3]。在实际工作中,使用HPGe γ谱仪系统对高活度样品进行测量时,在γ能谱中发现了全能峰右侧明显高于左侧的现象,猜测这种现象可能会为核素活度的准确分析引入较大误差。

对于此现象,国内外尚未有详细分析的文献报道。为此,本课题对该现象的产生原因进行了理论研究,以求在γ能谱分析中活动更加准确地定量分析结果。

1 实验现象介绍

所用的HPGe γ谱仪系统由法国堪培拉公司(Canberra Industries)生产的BE3830 P型HPGe探测器、Inspect 2000多道脉冲幅度分析器、777A型铅屏蔽室、计算机组成,谱分析系统软件“Genie-2000”。某样品的γ能谱全能峰右侧高于左侧现象(图 1(a)),这与平常工作中认识到的全能峰左侧高于右侧的现象(图 1(b))明显不同。

图 1 全能峰右侧高于左侧(a)、左侧高于右侧(b)的γ能谱 Figure 1 The energy spectrum which right is higher than left (a), left is higher than right (b).
2 理论分析

γ射线与物质相互作用时,单次散射的几率远大于多次散射[4],因此在考虑康普顿散射时,只考虑单次散射。由Klein-Nishina公式得到入射探测器的γ射线在E'min-E'max能区的康普顿散射总截面σe, c,见式(1)。

$ \begin{array}{c} {\sigma _{{\rm{e}}, c}} = {\rm{ \mathsf{ π} }}r_0^2\frac{1}{{{a^2}}}\frac{1}{{{E_0}}}\int_{{{E'}_{{\rm{min}}}}}^{{{E'}_{{\rm{max}}}}} {[\frac{2}{a} + \frac{1}{{{a^2}}} + (a-\frac{2}{a}-2)\frac{{{E_0}}}{{E'}} + } \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\frac{{E_0^2}}{{{{E'}^2}}} + \frac{1}{a}\frac{{E'}}{{{E_0}}}]{\rm{d}}E' \end{array} $ (1)

式中:r0为经典电子半径,r0=e2/(mec2)=2.8×10-15m,其中me为电子静止质量;a=E/(mec2),其中E为入射γ射线的能量;E0为电子静止能量,E0=mec2E′为散射光子的能量。

假设关注全能峰对应的能量为Em,在其右侧还有能量为Em+1Em+2,…,Em+KK是能量大于Em的全能峰的个数)的全能峰,则由式(1)可知,能量为Ei的γ射线在Em左侧和右侧某段康普顿坪上的散射截面σiLσiR可分别用式(2)、(3)描述。

$ \begin{array}{l} {\sigma _{i{\rm{L}}}}{\rm{ = \mathsf{ π} }}r_0^2\frac{1}{{a_i^2}}\frac{1}{{{E_0}}}\int_{{E_{m{\rm{LL}}}}}^{{E_{m{\rm{LR}}}}} {[\frac{2}{{{a_i}}} + \frac{1}{{a_i^2}} + ({a_i}-\frac{2}{{{a_i}}}-2)\frac{{{E_0}}}{{E'}} + } \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\frac{{E_0^2}}{{{{E'}^2}}} + \frac{1}{{{a_i}}}\frac{{E'}}{{{E_0}}}]{\rm{d}}E' \end{array} $ (2)
$ \begin{array}{l} {\sigma _{i{\rm{R}}}}{\rm{ = \mathsf{ π} }}r_0^2\frac{1}{{a_i^2}}\frac{1}{{{E_0}}}\int_{{E_{m{\rm{RL}}}}}^{{E_{m{\rm{RR}}}}} {[\frac{2}{{{a_i}}} + \frac{1}{{a_i^2}} + ({a_i}-\frac{2}{{{a_i}}}-2)\frac{{{E_0}}}{{E'}} + } \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\frac{{E_0^2}}{{{{E'}^2}}} + \frac{1}{{{a_i}}}\frac{{E'}}{{{E_0}}}]{\rm{d}}E' \end{array} $ (3)

式中:EmLLEmLREmRLEmRR是能量为Ei的γ射线发生康普顿散射时,γ能谱中能量为Em的全能峰左侧和右侧的康普顿坪上的关注能量区间取值,左侧的取值为EmLL-EmLR,右侧的取值为EmRL-EmRR

EmLR-EmLLEmEmRR-EmRLEm时,E′Em,式(2)、(3)可以分别写为式(4)、(5)。

$ \begin{align} &{{\sigma }_{i\rm{L}}}\approx {\rm{ = \mathsf{ π} }}r_{0}^{2}\frac{1}{a_{i}^{2}}\frac{1}{{{E}_{0}}}[\frac{2}{{{a}_{i}}}+\frac{1}{a_{i}^{2}}+({{a}_{i}}-\frac{2}{{{a}_{i}}}-2)\frac{{{E}_{0}}}{{{E}_{m}}}+ \\ &\ \ \ \ \ \ \ \ \ \ \ \ \frac{E_{0}^{2}}{E_{m}^{2}}+\frac{1}{{{a}_{i}}}\frac{{{E}_{m}}}{{{E}_{0}}}]\left( {{E}_{m\rm{LR}}}-{{E}_{m\rm{LL}}} \right) \end{align} $ (4)
$ \begin{align} &{{\sigma }_{i\rm{R}}}\approx {\rm{ = \mathsf{ π} }}r_{0}^{2}\frac{1}{a_{i}^{2}}\frac{1}{{{E}_{0}}}[\frac{2}{{{a}_{i}}}+\frac{1}{a_{i}^{2}}+({{a}_{i}}-\frac{2}{{{a}_{i}}}-2)\frac{{{E}_{0}}}{{{E}_{m}}}+ \\ &\ \ \ \ \ \ \ \ \ \ \ \ \frac{E_{0}^{2}}{E_{m}^{2}}+\frac{1}{{{a}_{i}}}\frac{{{E}_{m}}}{{{E}_{0}}}]\left( {{E}_{m\rm{RR}}}-{{E}_{m\rm{RL}}} \right) \end{align} $ (5)

则能量为Ei的γ射线平均每个γ光子在能量为Em的全能峰的左侧EmLR-EmLL能区和右侧EmRR-EmRL能区产生康普顿散射光子的概率PiLPiR分别可用式(6)、(7)描述。

$ {{P}_{i\rm{L}}}=1-{{\rm{e}}^{-{{\sigma }_{i\rm{L}}}{{N}_{0}}d}} $ (6)
$ {{P}_{i\rm{R}}}=1-{{\rm{e}}^{-{{\sigma }_{i\rm{R}}}{{N}_{0}}d}} $ (7)

式中:N0为探测器晶体的单位体积内的原子个数,m-3d为晶体厚度,m。因为σiLN0d≪ 1、σiRN0d ≪1,故有PiLσiLN0dPiRσiRN0d

假设N为测量时长t内进入探测器的光子数,则能量为Ei的γ射线在能量为Em的全能峰左侧和右侧对应能区上产生的康普顿散射光子计数率niLniR分别为:

$ {{n}_{i\rm{L}}}=\frac{1}{t}N{{P}_{i\rm{L}}} $ (8)
$ {{n}_{i\rm{R}}}=\frac{1}{t}N{{P}_{i\rm{R}}} $ (9)

其中,N的计算如下:

$ N=\int_{0}^{t}{{{A}_{i}}{{\rm{e}}^{-{{\lambda }_{i}}t}}\varepsilon {{P}_{{{\gamma }_{i}}}}{\rm{d}}t={{A}_{i}}\varepsilon {{P}_{{{\gamma }_{i}}}}}\int_{0}^{t}{{{\rm{e}}^{-{{\lambda }_{i}}t}}{\rm{d}}t=\frac{{{A}_{i}}\varepsilon {{P}_{{{\gamma }_{i}}}}}{{{\lambda }_{i}}}\left( 1-{{\rm{e}}^{-{{\lambda }_{i}}t}} \right)} $ (10)

式中:Ai为能量Ei的γ射线对应的核素的活度,Bq;Pγi为能量Ei的γ射线的发射概率;λi为该核素的衰变常数;ε为光子入射到探测器表面的概率,对具体的探测器和测量位置条件,可看作常数。将式(10)代入式(8)、(9)可得:

$ {{n}_{i\rm{L}}}=\frac{1}{t}{{A}_{i}}\frac{{{P}_{{{\gamma }_{i}}}}\varepsilon \left( 1-{{\rm{e}}^{-{{\lambda }_{i}}t}} \right)}{{{\lambda }_{i}}}{{P}_{i{\rm{L}}}}=\frac{{{N}_{0}}d}{t}{{A}_{i}}{{k}_{i}}{{\sigma }_{i\rm{L}}} $ (11)
$ {{n}_{i\rm{R}}}=\frac{1}{t}{{A}_{i}}\frac{{{P}_{{{\gamma }_{i}}}}\varepsilon \left( 1-{{\rm{e}}^{-{{\lambda }_{i}}t}} \right)}{{{\lambda }_{i}}}{{P}_{i{\rm{R}}}}=\frac{{{N}_{0}}d}{t}{{A}_{i}}{{k}_{i}}{{\sigma }_{i\rm{R}}} $ (12)

式中:ki为核素发射的平均每个γ光子进入探测器的概率,ki=Pγiε(1-e-λt)/λi。能量为Em的全能峰左侧对应能区的康普顿散射光子计数率nL是能量为EmEm+1Em+2,…,Em+K的γ射线的总贡献,右侧对应能区的康普顿散射光子计数率nR是能量为Em+1Em+2,…,Em+K的γ射线的总贡献,且由式(4)、(5)可知,当能量为Em的全能峰左右两侧关注的康普顿散射能区取相同值,即EmLR-EmLL=EmRR-EmRL时,有σiLσiR,因此,这两个能区中的康普顿散射光子计数率比值可用式(13)表示。

$ \begin{align} &\frac{{{n}_{\rm{L}}}}{{{n}_{\rm{R}}}}=\frac{\sum\nolimits_{m}^{m+K}{{{n}_{i\rm{L}}}}}{\sum\nolimits_{m+1}^{m+K}{{{n}_{i\rm{L}}}}}=\frac{\frac{{{N}_{0}}d}{t}\sum\nolimits_{m}^{m+K}{{{k}_{i}}{{A}_{i}}{{\sigma }_{i\rm{L}}}}}{\frac{{{N}_{0}}d}{t}\sum\nolimits_{m+1}^{m+K}{{{k}_{i}}{{A}_{i}}{{\sigma }_{i\rm{R}}}}} \\ &\ \ \ \ \ =\frac{{{k}_{m}}{{A}_{m}}{{\sigma }_{m\rm{L}}}+\sum\nolimits_{m+1}^{m+K}{{{k}_{i}}{{A}_{i}}{{\sigma }_{i\rm{L}}}}}{\sum\nolimits_{m+1}^{m+K}{{{k}_{i}}A{}_{i}{{\sigma }_{i\rm{R}}}}}\approx 1+\frac{{{k}_{m}}{{A}_{m}}{{\sigma }_{m\rm{L}}}}{\sum\nolimits_{m+1}^{m+K}{{{k}_{i}}{{A}_{i}}{{\sigma }_{i\rm{R}}}}} \\ \end{align} $ (13)

Em的全能峰对应核素的活度Am很小时,${{k}_{m}}{{A}_{m}}{{\sigma }_{m\rm{L}}} \ll \sum\nolimits_{m+1}^{m+K}{{{k}_{i}}{{A}_{i}}{{\sigma }_{i\rm{R}}}}$。由式(13)可知,此时nL/nR≈1,即关注核素活度很低时,其对应全能峰左右两侧高度近似相等;随着Am增大,nL/nR逐渐增大,全能峰左侧逐渐高于右侧(图 1(b))。因此,如果没有其他效应的影响,全能峰左侧高度将随核素活度增大而越来越大于其右侧。根据偶然符合计数计算公式[5]M=2τnM为符合率;τ为谱仪分辨时间;n为产生偶然符合的两种能量γ光子计数的乘积),可将全能峰左右两侧关注能区的康普顿计数比值关系用式(14)表示。

$ \begin{align} &\frac{{{n}_{\rm{L}}}}{{{n}_{\rm{R}}}}=\frac{\frac{{{N}_{0}}d}{t}\sum\nolimits_{i=m}^{m+K}{{{k}_{i}}{{A}_{i}}{{\sigma }_{i\rm{L}}}}-\frac{1}{t}2\tau {n}'{{k}_{m}}{{A}_{m}}}{\frac{{{N}_{0}}d}{t}\sum\nolimits_{i=m+1}^{m+K}{{{k}_{i}}{{A}_{i}}{{\sigma }_{i\rm{R}}}}+\frac{1}{t}2\tau {n}'{{k}_{m}}{{A}_{m}}} \\ &\ \ \ \ \approx 1-\frac{2{{k}_{m}}\left( 2\tau {n}'-{{N}_{0}}d{{\sigma }_{m\rm{L}}} \right)}{\frac{{{N}_{0}}d}{{{A}_{m}}}\sum\nolimits_{i=m+1}^{m+K}{{{k}_{i}}{{A}_{i}}{{\sigma }_{i\rm{R}}}}+2\tau {n}'{{k}_{m}}} \\ \end{align} $ (14)

式中:n′为γ能谱低能段计数。由式(14)可知,当2τn'-N0mL < 0时,nL/nR的值随着关注核素活度Am而增大,全能峰左侧会比右侧越来越高;当2τn'-N0mL > 0时,nL/nR的值随着关注核素活度Am增大而减小,全能峰右侧会比左侧越来越高,出现图 1(b)所示现象。

3 实验验证

实验布局见图 2,将241Am/60Co(活度715Bq/186 Bq)标准混合点源放置于探测器轴线探窗上,将137Cs(活度3.7×105Bq)标准点源放置于混合点源上方,中心点与混合点源中心点对齐,调整137Cs点源与探窗距离,分别测得不同距离对应的γ能谱,并在每个γ能谱的661.7 keV全能峰左右两侧康普顿坪上选取不同能区宽度计算nL/nR。同一个γ能谱中,能区宽度EmLR-EmLL=EmRR-EmRL,结果见表 1

图 2 实验示意图 Figure 2 Schematic of the experiment.
表 1 nL/nR计算结果 Table 1 The calculated results of nL/nR.

表 1可知,在干扰核素活度不变时,对于不同的能区宽度,均有在死时间小于等于1.0%时,nL/nR≈1;死时间位于1.0%-4.0%时,nL/nR随死时间而增大;当死时间位于4.0%-56.0%时,此时关注核素的活度远大于干扰核素,nL/nR随死时间增大而减小。该结果与理论分析结果吻合,表明理论分析结果是正确的。并从表 1可以看出,全能峰左侧高于右侧的现象出现在死时间小于等于20%时;当死时间大于等于20%时,全能峰右侧开始高于左侧。

为初步估计这种现象对分析峰面积和最终活度定值的影响,分别使用本次实验测得的探测效率与无源效率刻度软件LabSOCS的计算值进行对比,相对偏差随死时间变化曲线见图 3

图 3 探测效率相对偏差随死时间变化 Figure 3 The difference of efficiency between LabSOCS and experiment with time.

图 3可见,对于点源,软件LabSOCS计算的探测效率与实验测得的探测效率的相对偏差随着死时间的增大,在7%上下浮动。所以认为,右高左低这种现象并不会对点源的能谱测量造成影响。

4 结语

针对γ能谱分析中出现的γ全能峰右侧高于左侧的现象,进行了详细的理论分析,并通过实验验证了理论分析的正确性。分析表明:全能峰右侧高于左侧的现象由偶然符合现象引起。实验验证结果表明,在高纯锗探测器测量中,对于点源γ能谱,全能峰左侧高于右侧的现象主要出现在死时间小于等于20%时,并且高度比值先增后减;当死时间大于等于20%时,全能峰右侧开始高于左侧,结合理论分析结论可知,右侧高度将越来越大于左侧高度。而对于低能射线(比如241Am),由于进入探测器的前端散射比较多,所以不容易出现全能峰右高左低现象。随着死时间增长,当出现全能峰右侧高于左侧现象时,并不会对点源的γ能谱分析造成影响。

参考文献
[1]
吴伦强, 韦孟伏, 张连平, 等. γ能谱技术鉴定放射性核素的应用研究[J]. 核电子学与探测技术, 2009, 29(5): 931-934.
WU Lunqiang, WEI Mengfu, ZHANG Lianping, et al. Application research of γ spectrometry on the authentication of radionuclide[J]. Nuclear Electronics & Detection Technology, 2009, 29(5): 931-934.
[2]
伍怀龙, 刘素萍, 刘晓亚, 等. γ能谱解谱技术研究[J]. 核技术, 2005, 28(6): 430-434.
WU Huailong, LIU Suping, LIU Xiaoya, et al. Investigation of gamma spectra analysis[J]. Nuclear Techniques, 2005, 28(6): 430-434.
[3]
冯天成, 吴睿, 苏川英, 等. 就地γ能谱中天然核素峰谷比随砂土容重的变化[J]. 核电子学与探测技术, 2013, 33(4): 417-424.
FENG Tiancheng, WU Rui, SU Chuanying, et al. PVR's variety with soil-bulk density in in-situ γ spectra of natural radionuclides[J]. Nuclear Electronics & Detection Technology, 2013, 33(4): 417-424.
[4]
冯天成, 贾明雁, 冯元举, 等. 就地γ谱仪峰谷比法测量137Cs深度分布[J]. 原子能科学技术, 2011, 45(6): 722-729.
FENG Tiancheng, JIA Mingyan, FENG Yuanju, et al. Peak-to-valley ratio method to determine depth distribution of 137Cs in soil by in-situ γ spectrometry[J]. Atomic Energy Science and Technology, 2011, 45(6): 722-729.
[5]
李德, 过雅南. 随机信号偶然符合计数率公式的推导[J]. 核电子学与探测技术, 2000, 20(1): 1-4.
LI De, GUO Yanan. Formula to compute the accidental coincidence rate for random signals[J]. Nuclear Electronics & Detection Technology, 2000, 20(1): 1-4.