1. 引言
γ标准放射源活度的准确测量是γ放射性活度绝对测量计量传递的关键技术。60Co放射源被广泛应用于放射性临床诊断、治疗 [1] 、生物化学、动植物学及核科学等的基础研究,放射源活度测量的准确性直接影响到患者的治疗效果及科研结果的科学性。在不同领域的使用,对60Co放射源的测量精度有不同的要求;因此在进行相应的基础研究或诊断、治疗之前需要确定实验条件并选择合适的实验装置对放射源的活度进行测量 [2] 。符合法作为一种放射源活度的绝对测量方法在各个领域中得到广泛应用;开设核学科的高校都设有符合法测量源活度的实验,在进行符合法测源活度时,本科生并不能准确、快速的将源的活度测量出来;本文意在给出符合法测量源活度数据上的参考,并讨论出高校一般符合测量系统能够达到的最高测量精度。
2. 符合法原理
γ-γ符合法测量60Co放射源的活度,实际上就是利用β衰变产生,退激放出2个γ粒子γ1、γ2之间的符合;如图1(a),60Co的衰变为级联辐射,60Co经过β衰变到60Ni的激发态,通过放出一组级联γ1、γ2射线退激到60Ni的基态 [3] ;如图1(b),用两个探测器分别探测源放出的γ射线,输出的信号送入符合电路。
实验所用实验装置如图2(a)、图2(b),两个探头都为NaI(TI)闪烁探测器。
2.1. 放射源活度的计算公式推导
令60Co放射源的活度为A,实验装置第1道,第2道和符合道所接定标器测量折合的计数率分别为nγ1、nγ2、nγγ,设ε11、ε21分别为第1道对γ1、γ2射线的探测效率,ε12、ε22分别为第2道对γ1、γ2射线的探测效率,Ωγ1、Ωγ2分别为探测器1、探测器2对放射源所张的立体角;则有如下关系式:
第1道计数率:
(1)
第2道计数率:
(a)
(b)
Figure 1. (a) The decay plan of 60Co; (b) the illustration of coincidence method
图1. (a) 60Co衰变纲图;(b) 符合法示意图
(a)
(b)
Figure 2. (a) The block diagram of γ-γ coincides measurement system [4] ; (b) the physical map of γ-γ coincidence measurement system
图2. (a) γ-γ符合测量实验系统框图 [4] ; (b) γ-γ符合测量实验系统实物图
(2)
符合道计数率:
(3)
式(1)、(2)代入式(3),可得:
(4)
如果γ1、γ2的能量相近时有
,
,则
(5)
为了计算结果更加准确,在此基础上对nγ1、nγ2、nγγ进行本底修正、偶然符合计数修正后 [5] ,则
(6)
3. 最佳测量时间的测量
测量时间是符合测量过程中一个重要参数;设置的时间太长使得测量效率太低,太短则测量数据的统计涨落影响测量结果的准确度,因此需要设置一个恰当测量时间。
将测量时间依次设置为10 s、20 s、30 s、40 s、50 s、100 s、200 s,重复进行下述实验。
3.1. 探测器与仪器预热
按照图1连接实验装置图;将2个插件机箱的电源打开,仪器预热30分钟。测量探头的工作坪曲线,确定最佳工作电压,最佳工作电压取坪区1/3或1/2 [6] ,这里最佳工作电压取坪区1/3;计算得探头A的最佳工作电压为550 V,探头B的最佳工作电压为540 V,分别将高压插件 [7] 调至对应探头的最佳工作电压。
对各个道对应的定标器进行自检,定标器选择为“自检”,定时时间设为1 s (定器:k = 1,n = 0),若仪器显示“19,456”,则证明定标器计数正常 [8] ;
分别将137Cs源分别置于NaI (TI)探头1与NaI (TI)探头2之前,见图3(a);2探头之间放置铅砖,减小源与源之间的干扰;将2个NaI (TI)探头输出 [9] 、信号放大器的输出信号 [10] 、2个单道的输出信号 [11] 和符合插件 [12] 的输出信号分别接入示波器观察是否有信号输出、输出的信号是否失真。
(a)
(b)
Figure 3. (a) The position of 137Cs with probe; (b) the position of 60Co with probe
图3. (a) 137Cs与探头位置;(b) 60Co与探头位置
3.2. 60Co放射源活度测量
3.2.1. 脉冲宽度的设置
1) 单道的工作方式选择为“积分”;
2) 调节放大器增益的粗调和微调使放大器的输出脉冲的幅度为3 V;
137Cs输出的脉冲幅度为3 V时,在测量60Co源活度时,保持放大器的放大倍数不变,60Co源能量为1.137 MeV全能峰对应脉冲幅度大约6 V;由于单道的甄别阈为0~10 V,如若137Cs脉冲幅度太小,可能造成全能峰的重叠 [13] ;若137Cs脉冲幅度太大,60Co的脉冲幅度更大,使得单道的阈值对于60Co来说毫无意义;因此本实验将137Cs的脉冲幅度调为3 V。
3) 调节单道的甄别阈值为放大器输出脉冲幅度的75% (保证全能峰的计数都能够被记录下来;又去掉了X射线峰、反散射峰的计数贡献) (图4);
理想情况下阈值选择V3位置最好,但是由于探测器探测的全能峰会偏移 [14] ,全能峰可能产生的漏计数;如若阈值选择V1位置峰的偏移,可能会记录到反散射峰贡献的计数;彭朝华等 [15] 证实康普顿坪计数对活度测量影响为3%~5%;因此选择V2,V2靠近V3,既不会产生全能峰漏计数,也能减小康普顿坪计数对活度测量的影响;在上述工作电压、放大倍数的条件下,通过多道脉冲幅度分析器,得到137Cs该条件下的能谱;通过能谱测量得到V4 ≈ 2.8 V,V3 ≈ 2.45 V,V2 ≈ 2.1 V,V2/V4 ≈ 75%,峰偏移时,各个位置对应与脉冲幅度之间的比例不会改变,则V2的位置可近似看做脉冲幅度的75%。
4) 调节符合插件上符合成形时间,经过2017年殷志成等(待发表)大量实验证明,该符合测量系统的最佳脉冲个宽度为0.39 μs;因此将符合电路插件I和II路的“监测”输出脉冲宽度调到0.39 μs。
3.2.2. 偶然符合计数及符合道符合计数测量
1) 将两个137Cs源拿掉,并使60Co放射源夹于两个NaI (TI)闪烁探头之间,见图3(b);
2) 放大器增益的粗调和微调的倍数保持不变;单道工作方式仍为积分,通过示波器观察此时放大器输出脉冲幅度,调节对应单道的阈值为脉冲幅度的75%;
3) 固定将符合电路插件上I路延时固定某一位置,改变另一道的延时,延时从0.4 μs慢慢调节到最大处,若符合道的计数先是保持不变,在慢慢增大,增大后又慢慢减下,最后又保持不变,则最开始保持不变的计数即为偶然符合计数,该过程如图5所示;图中的1和4不产生符合计数,此时符合道的计数即为偶然符合计数。用定标器对符合道进行计数,计数时间为10秒,计数3次;
4) 再调节延迟使符合道脉冲计数为最大后,仪器预热30 s,用定标器对符合道进行计数,计数时间为10秒,计数3次
![](//html.hanspub.org/file/5-3150100x24_hanspub.png)
Figure 4. The 137Cs’γ spectrum and threshold selection
图4. 137cs的γ能谱及阈值选择示意图
5) 将源拿掉并置于远离探头的位置,对第1道、第二道和符合道的本底计数进行测量,用定标器对符合道进行计数,计数时间为10秒,计数3次。
6) 定标器的测量时间后,预热30秒,重复3.3的测量过程;得到的测量数据处理以后见表1。
γ-γ符合法测得的放射性活度A的物理量并不是直接得到的,它要通过诸如式(8)中各参量按着一定的函数关系计算得到。因此,实验中各个直接测量的物理量存在的误差会传递到间接测量的结果中去。所以对符合法测量误差进行分析可以考虑采用间接测量的误差传递公式。本符合测量中,物理量nγ1、nγ2中包含nγγ,所以这3个量是相关。为了选择独立的变量,可把nγ1分成两部分:一部分是符合计数nγγ,另一部分是未符合的计数用nx表示。对nγ2也是这样分成两部分即nγγ和ny。于是nγ1、nγ2可以写成。
(7)
式中,nx、ny、nc、nb1、nb2、nb、nrc均为互相独立的变量,根据误差传递公式可得A的方差 [5] :
(8)
其相对误差:
(9)
由表2可知:测量时间越长测量得到60Co放射源的相对误差越小,计数统计误差越小;测量时间100 s、200 s的相对误差分别为0.35%、0.24%,2者之间相对误差接近,但是200 s的测量时间要比100 s的测量时间多100 s,极大的减慢了测量的速度;为确保测量的效率,又使统计误差的影响降低;测量时间取100 s。
4. 60Co放射源活度测量
实验室构建两套符合测量实验系统。又偶然符合计数为式
,将式(6)计算放射源的公式变形成如下所示:
![](Images/Table_Tmp.jpg)
Table 1. The Calculated of nb1, nb2, nγ1, nγ2, nγγ, nrc and A
表1. nb1, nb2, nγ1, nγ2, nγγ, nrc及A计算值
![](Images/Table_Tmp.jpg)
Table 2. Measured A, σA and δA by different measurement times
表2. 不同测量时间测量得到的A, σA及δA
(10)
重复上述3.3.1的步骤,经测量计算得到2套符合测量系统探头最佳工作电压如表3所示。
4.1. 符合分辨时间τ的测量
1) 单道的工作方式选择为“积分”;调节放大器增益的粗调和微调使放大器的输出脉冲的幅度为3 V;调节单道的甄别阈值为放大器输出脉冲幅度的75%;
2) 将符合电路插件I和II路的监测输出接到示波器上,调节符合电路插件上的符合成形时间,使得Ⅰ和Ⅱ的监测输出脉冲的宽度为0.39 µs;
3) 分别调节三个定标器的定时时间设为100秒,并同时对第1道、第2道、符合道进行计数,计数3次,对源距探头5个不同距离重复测量;每次改变条件后,仪器预热100 s,再重新测量。
计数的相对误差要求小于1%,即计数率要大于100 s−1;改变源距探头距离的时候,要求符合道的计数有显著的改变,若改变的不显著,可能造成多点重合,增大拟合结果的误差(表4)。
根据公式(7),使用最小二乘法拟合得到结果见表5。
4.2. 符合计数及各道本底计数测量
1) 将两个137Cs源拿掉,并使60Co放射源夹于两个NaI (TI)闪烁探头之间,见图3(b)。
![](Images/Table_Tmp.jpg)
Table 3. The best working voltage of each probe
表3. 各探头最佳工作电压
![](Images/Table_Tmp.jpg)
Table 4. The calculate data of 2nγ1, nγ2 and nrc
表4. 2nγ1, nγ2, nrc计算数据
![](Images/Table_Tmp.jpg)
Table 5. Obtained τ by fitting 5 points
表5. 5点拟合得到的τ
2) 放大器增益的粗调和微调的倍数保持不变;单道工作方式仍为积分,通过示波器观察此时放大器输出脉冲幅度,调节对应单道的阈值为脉冲幅度的75%。
3) 固定将符合电路插件上Ⅰ路延时固定某一位置,改变另一道的延时,延时从0.4 μs慢慢调节到最大处,同时记录各个延迟时刻的符合计数,找到符合计数最大处。
4) 用定标器同时对符合道、第1道、第2道进行计数,计数时间为100秒,计数3次,改变条件后,仪器预热100 s,测量数据经计算后,结果见表6。
![](Images/Table_Tmp.jpg)
Table 6. The calculated of nb1, nb2, nγ1, nγ2 and nγγ
表6. nb1, nb2, nγ1, nγ2, nγγ计算值
5) 计算60Co放射源的活度A
将处理以后的数据带入式(7),可得60Co放射源的活度A1、A2、A3,相关计算数据计算值如表7所示。
本次实验室2套符合测量系统均对同一个60Co的放射源活度进行测量,且各次测量时间间隔较短,可认为在短时间60Co源的活度是不变的;测量误差都是由于仪器与统计误差造成。则引入塞贝尔公式 [16] :
(11)
其相对偏差为:
(12)
5. 结果分析
从表5可知,γ-γ符合测量实验系统在脉冲宽度设为0.39 µs,理论上符合测量系统的符合分辨时间应该为0.39 µs,但是2套符合测量装置测量数据拟合得出的符合分辨时间都小于0.39 µs,在0.37 µs上下波动,因此每次测量源活度之前都必须重新测量该符合测量的符合分辨时间。
该60Co放射源出厂时间为2015年1月,出厂活度为1.96 × 105 Bq;2017年6月进行了本次实验,计算得到实验时源的参考活度为141,086 Bq。由表8可知:第一套、第二套符合测量系统测量数据计算出的60Co放射源的活度分别为161,099 Bq、142,620 Bq,2套符合测量系统使用的插件型号均相同,实验过程中设置的参数均相同;理论上2套符合测量系统测出来的60Co放射源的活度应在误差范围内一致,但第一套符合测量系统测量的活度值要远远大于第二套符合测量系统;且第二套符合测量系统测量结果与参考源的活度相对偏差为1.09%在误差范围内一直,如若2套符合测量系统的第一道与第二道计数相同则第一套的符合测量系统真符合计数必定小于第二套符合测量系统。
符合插件其本身具有时间分辨本领,符合插件的符合分辨时间不能随意减小,本次实验2套符合测量系统设置的符合分辨时间都为0.39 µs,测量时间相同,第二套符合测量系统的真符合计数大于第一套符合测量系统,则说明第一套符合测量系统真符合计数产生了“漏计数”,即第二套符合测量系统的符合插件的性能好于第一套符合测量系统。
比较两个符合插件的性能好坏,可同时构建2套符合测量系统测量;先测量第一套符合测量系统的最佳脉冲宽度;再同时用2套符合测量系统测量同一个60Co放射源的活度,2套符合测量系统的符合分辨时间都设为第一套符合测量系统测出的最佳脉冲宽度且其他参数设置相同;通过比较两套符合测量系统测量出的放射源的活度;若第二套符合测量系统测量活度小于第一套系统,则第二套符合测量系统的符合插件性能好于第一套系统;反之第一套符合测量系统的符合插件的性能好于第二套。
![](Images/Table_Tmp.jpg)
Table 7. The calculated of A1, A2 and A3
表7. A1, A2, A3计算值
![](Images/Table_Tmp.jpg)
Table 8. The calculated of A ¯ , σ s and δ A ¯
表8.
,
及
计算值
文献报道中1990年许茂国 [17] 利用γ-γ符合法测量放射源60Co的活度,测量结果的相对偏差为1.4%;2002年程敏熙 [18] 使用γ-γ符合法对60Co放射源的活度进行测量,测量结果的相对偏差为3.93%;2003年马文彦等 [19] 通过β-γ符合法测量60Co放射源的活度,测量结果的相对偏差4.18%;2013年颜拥军等 [20] 利用数字符合法测量60Co放射射源的活度,测量结果的相对偏差为2.51%。
该实验测量出的60Co放射源的活度为142620Bq,相对偏差为0.217%,与其精度要比上述学者测量出的60Co放射源的精度都要高;与实验时放射源的参考活度的相对偏差为1.09%,在误差范围内一致;通过此实验步骤及插件的相关参数的设置,能够快速、准确的测量出60Co放射源的活度。
6. 结论
本文详细介绍了γ-γ符合法测量60Co放射源活度的原理,以及测量活度的实验步骤;阐述了为保证测量精度符合测量系统各插件相关参数设置方法;通过实验可知每套符合套测量系统都有其固有的符合分辨本领,实验之前必须找到其最佳脉冲宽度,且每次进行实验时都需对符合测量系统的符合分辨时间重新进行测量;本次实验测量出的60Co放射源活度为142,620 ± 924 Bq,测量精度为308 Bq,测量结果的相对偏差为0.217%;实验所得的源活度与实验时源的参考活度相对偏差为1.09%,在误差范围内一致,测量结果稳定可靠,实验精度要好于其他学者测量结果的精度。
基金项目
湖南省科技计划项目( 2015SK20063)。
NOTES
*通讯作者。