中国测试  2019, Vol. 45 Issue (6): 13-18

文章信息

滕忠斌, 田丽霞, 宋明哲, 倪宁, 张曦, 魏可新
TENG Zhongbin, TIAN Lixia, SONG Mingzhe, NI Ning, ZHANG Xi, WEI Kexin
$\gamma $ 射线空气比释动能空腔理论计算值准确度的蒙特卡罗研究
Study of the accuracy of γ-rays air kerma calculated by cavity theroy using Monte-Carlo method
中国测试, 2019, 45(6): 13-18
CHINA MEASUREMENT & TEST, 2019, 45(6): 13-18
http://dx.doi.org/10.11857/j.issn.1674-5124.2018090024

文章历史

收稿日期: 2018-09-10
收到修改稿日期: 2018-10-15
$\gamma $ 射线空气比释动能空腔理论计算值准确度的蒙特卡罗研究
滕忠斌1,2 , 田丽霞1 , 宋明哲2 , 倪宁2 , 张曦2 , 魏可新2     
1. 东华理工大学,江西 南昌 330013;
2. 中国原子能科学研究院,北京 102413
摘要:当使用基于空腔理论的石墨空腔电离室绝对测量 $\gamma $ 射线空气比释动能时,需要考虑空腔体积的大小变化对测量结果准确度的影响。利用EGSnrc程序计算得到10~5 000 cm3的球型石墨空腔电离室的石墨和空气的限制的阻止本领比 $(\overline{L}_ \Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$ 、入射光子与空腔气体直接作用沉积能量的份额Fair和Spencer-Attix空腔理论修正因子kSA。结果表明对于60Co和137Cs能量的光子,在计算 $(\overline{L}_ \Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$ 时选择合适的 $\Delta $ 值可以使kSA的值保持在0.998~1.000之间。此外kSA的大小随空腔体积的增大有逐渐减小的趋势。表明在设计和使用较大体积(≥50 cm3)的石墨空腔电离室时需要考虑合适大小的kSA值,以保证空腔理论计算值和绝对测量结果的准确。
关键词空腔理论    限制的阻止本领比    EGSnrc    电离室    
Study of the accuracy of γ-rays air kerma calculated by cavity theroy using Monte-Carlo method
TENG Zhongbin1,2 , TIAN Lixia1 , SONG Mingzhe2 , NI Ning2 , ZHANG Xi2 , WEI Kexin2     
1. East China University of Technology,Nanchang 330013, China;
2. China Institute of Atomic Energy, Beijing 102413, China
Abstract: When people absolutely measure the ${\rm{\gamma}} $ -rays air kerma using a cavity chamber based on cavity theory, it is necessary to consider the impact of the different sizes of the cavity volume on the accuracy of the measurement results. Based on the calculation of EGSnrc code,the restricted stopping-power ratios graphite to air, $(\overline{L}_ \Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$ ,the fraction of dose caused by photon interacting with cavity gas,Fair,and the Spencer–Attix correction factor, kSA,have been obtained,when the cavity volume changes from 10 cm3 to 5000 cm3.The results show that for photons of 60Co and 137Cs source, selecting the appropriate $\Delta $ value,during the calculation of $(\overline{L}_ \Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$ ,can keep the kSA between 0.998-1.000.And the kSA gradually decreases with the increase of cavity volume.Thus, to obtain accurate results of cavity theory calculation and absolute measurement, the appropriate value of kSA needs to be considered carefully when designing or using a large volume (≥50 cm3) graphite-walled cavity ionization chamber.
Key words: cavity theory     restricted mass stopping power ratios     EGSnrc     cavity chamber    
0 引 言

为了能够使用空气比释动能标准复现 ${\rm{\gamma}} $ 射线参考辐射场中的低空气比释动能率,实验室中通常使用替代法进行测量[1]。虽然该方法现实可行,但每次次级标准电离室的校准都会引入测量不确定度。由于基准(或最高标准)电离室对 ${\rm{\gamma}} $ 射线空气比释动能的绝对测量都是基于空腔理论,如果将电离室的空腔体积增大来得到更低的绝对测量下限从而减少校准次数,那么此时需要考虑空腔体积的变化对空腔理论计算结果准确度的影响。

MA[2]在1991年提出使用光子与空腔气体作用所沉积能量占空腔气体中总沉积能量的份额Fair来判断一个空腔是否可以作为理想条件下的Bragg-Gray空腔。但该因子只是考虑不同条件时的空腔是否可以作为满足空腔理论前提条件的理想空腔,并未说明当该空腔作为实际电离室空腔时是否可以准确测量。Borg[3]和Russa[4]在2000年和2009年先后提出,60Co和192Ir光子能量和较小的空腔体积时使用Spencer-Attix空腔理论修正因子kSA或注量扰动因子kfl来判断空腔理论计算值与真实值之间的差异。引入kSA的原因是实际空腔的存在打破了Bragg-Gray理论的前提条件[3~5]kSA可以用于判断基于Spencer-Attix空腔理论的石墨空腔电离室是否可以准确测量 ${\rm{\gamma}} $ 射线参考辐射场中某点处的空气比释动能。较小的kSA值对该点的空气比释动能的准确测量起到了关键作用。

研究目的是探究不同入射光子能量时,空腔电离室体积的变化对空腔理论计算值准确度的影响。使用EGSnrc计算不同光子能量和不同空腔体积时石墨和空气限制的阻止本领比(考虑了不同的密度效应修正和平均激发能推荐值)、必要的修正因子、FairkSA。计算结果可以用于判断较大空腔体积时空腔理论的计算结果是否准确,并且可以用于指导较大石墨空腔电离室(≥50 cm3)的制作,最终用于探究 ${\rm{\gamma}} $ 射线空气比释动能绝对测量下限和减小替代法得到的测量结果测量不确定度。

1 空腔理论

空腔理论可以将介质中的吸收剂量与探测器介质中的吸收剂量联系起来,使得可以用剂量计的读数来评估介质中的吸收剂量。实验室中通常使用基于Bragg-Gray空腔理论的石墨空腔电离室来直接测量 ${\rm{\gamma}} $ 辐射场中一点的空气比释动能[6]为:

$\begin{aligned} {K_{{\rm{air}}}} =& \left[(Q/{m_{{\rm{air}}}}){(\overline W /{{e}})_{{\rm{air}}}}{({\overline \mu _{{\rm{en}}}}/\rho )_{{\rm{air,c}}}} \times\right.\\ &\left.{(\overline s /\rho )_{{\rm{c,air}}}}/(1 - {\overline g _{{\rm{air}}}})\right]\prod\limits_{{i}} {{k_i}}\end{aligned} $ (1)

式中:Q——空腔内干空气中的电离电荷量;

$m_{\rm{air}} $ ——参考温度和压力条件下空腔中干空气的质量, ${m_{{\rm{air}}}} = {\rho _0}V$

ρ0——参考条件下干空气密度;

V——空腔的体积;

$(\overline W{{/}}e)_{\rm{air}}$ ——在干空气中生成一个离子对,电子所需要消耗的平均能量;

${\rm{(}}\overline \mu_ {\rm{en}}/\rho)_ {\rm{air,c}}$ ——空气与石墨的质能吸收系数之比;

${\rm{(}}\overline s {\rm{/}}\rho)_ {\rm{c,air}}$ ——石墨和空气的非限制的阻止本领之比;

${{\bar g}_{{\rm{air}}}}$ ——次级电子在空气中由于轫致辐射所消耗的能量占初始能量的份额;

$\prod\limits_i {{k_i}} $ ——因实际条件与参考条件的差异而对测量值进行的修正。

由于式(1)使用的是平行光子束,空腔内介质为干空气,所以在 $\prod {{k_i}} $ 中排除了空气湿度修正kh和轴向不均匀性修正kan,只保留壁效应修正kwall和Spencer-Attix理论修正kSAkwall可以修正室壁对初级光子的衰减和散射作用,用于重现电离室没有壁情况下空腔内的吸收剂量。kSA可以用来判断基于Spencer-Attix空腔理论的石墨空腔电离室在绝对测量时是否可以准确测量,也就是说kSA越接近于1,使用空腔理论计算的结果越准确。kSA的具体的表达式为:

$k_{\rm{SA}} = \frac{{k_{\rm{air}}(1 - \overline g_ {\rm{air}})}}{{D_{\rm{gas}}\left( {{{\overline S }}/{\rho }} \right)_{{\rm{gas}}}^{{\rm{wall}}}\left( {{{\overline\mu _{\rm{en}} }}/{\rho }} \right)_{{\rm{wall}}}^{{\rm{air}}}\prod {k_{{i}}} }}$ (2)

虽然式(1)使用了非限制的阻止本领比 ${\rm{(}}\overline s {\rm{/}}\rho) _{\rm{c,air}}$ ,但在国际计量局组织的 ${\rm{\gamma}} $ 射线空气比释动能国际比对中,各国基准使用的都是限制的阻止本领比[7-8]。那么此时使用的空腔原理亦称之为Spencer-Attix空腔理论。当空腔体积变化时,Spencer-Attix空腔理论较Bragg-Gray空腔理论更适用于空腔内沉积能量的计算。对于能量低于 $\Delta $ 的电子,认为它们不能穿越空腔,能量就地沉积[9]。在计算石墨中的吸收剂量与空腔气体中的吸收剂量时,针对能量低于 $\Delta $ 的电子使用了限制的阻止本领,二者之比为:

$ \begin{aligned} &\left( {\frac{{\overline {L_\Delta } }}{\rho }} \right)_{{\rm{air}}}^{\rm{c}} =\\& \frac{{\displaystyle\int_\Delta ^{E_{\max} } {{\varPhi_ {\rm{m}}}\left( E \right)\left( {\displaystyle\frac{{{{\overline L}_ {\Delta}} \left( E \right)}}{\rho }} \right)_{\rm{c}}{\rm d}E + {\varPhi_{\rm{c}}}\left( \Delta \right)\left( {\displaystyle\frac{{\overline S \left( \Delta \right)}}{\rho }} \right)_{\rm{c}}\Delta } }}{{\displaystyle\int_\Delta ^{E_{\max} } {\varPhi_ {\rm{air}}\left( E \right)\left( {\displaystyle\frac{{\overline L_\Delta \left( E \right)}}{\rho }} \right)_{\rm{air}}{\rm d}E + \varPhi_ {\rm{air}}\left( \Delta \right)\left( {\displaystyle\frac{{\overline S \left( \Delta \right)}}{\rho }} \right)_{\rm{air}}\Delta } }}\end{aligned}$ (3)

Spencer-Attix空腔理论与Bragg-Gray空腔理论有两个相同的前提条件:1)空腔足够小,空腔的存在不会改变介质中次级电子的注量和谱分布。2)光子在空腔中的作用可以忽略不计。但在实际 ${\rm{\gamma}} $ 射线辐射场中空腔内必然存在光子与空腔气体的作用,并且空腔的存在影响了次级电子的注量及其谱分布。当空腔体积逐渐增大时,实际空腔与理想空腔的差异将导致使用式(1)计算的准确度变低。因此在使用较大体积的石墨空腔电离室直接测量空气比释动能时,需要定性和定量地分析式(1)是否准确可靠。

2 方法和结果 2.1 计算模型

EGSnrc是被广泛应用于电离辐射剂量领域中的蒙特卡罗程序,其计算的电离室响应与真实值相差在0.1%之内[10]。使用EGSnrc中的用户程序计算式(2)中等号右边所有的物理量和修正因子,计算时的输运选项均使用缺省值。球型石墨空腔电离室被放置在真空中,空腔中充满干空气,石墨壁的厚度为3 mm,空腔体积大小在10~5 000 cm3之间变化。源项设置为0.662 MeV和1.25 MeV的单能光子以及两种参考辐射场的 ${\rm{\gamma}} $ 能谱(137Cs源加铅衰减后的 ${\rm{\gamma}} $ 能谱和60Co源加铅衰减后的 ${\rm{\gamma}} $ 能谱,分别记为Cs+Pb能谱和Co+Pb能谱),光子从电离室侧方平行入射,射束直径和室壁外径相等。计算模型如图1

图 1 EGSnrc计算使用的模型

在计算式(2)物理量和修正因子前,使用EGS_GUI内的用户程序创建计算时需要的PEGS4文件。它包含了计算所需要的材料的种类、密度、截面数据、密度效应修正、平均激发能和截止能量等参数。需要设置的截止能量包括:光子的输运截止能量PCUT和产生光子的阈值能量AP,电子输运截止能量ECUT和产生次级电子的阈值能量AE= $\Delta $ +511 keV,其中 $\Delta $ 与电离室的空腔大小有关。对于体积较小的电离室, $\Delta $ 通常取值为10 keV[5]。它等于电子在干空气中射程为L=4V/A时所对应的能量。LVA分别是空腔的平均弦长、体积和表面积。首先需要计算不同空腔体积对应的L值,然后利用电子在干空气中的CSDA射程和电子能量的关系插值得到相应的 $\Delta $ [9],结果如表1。在PEGS4文件中石墨密度为1.836 g·cm−3,使用ICRU37号报告[11]推荐的平均激发能Ig=78 eV,以及石墨在1.7 g·cm−3时的密度效应修正,空气的密度为1.025 kg·m−3

表 1 不同空腔体积的Δ值和不同条件时的石墨和空气的限制的阻止本领比
体积/
cm3
内径/
cm
外径/
cm
L/
cm
Δ/
keV
$(L_ \Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$
Cs+Pb 能谱 Ig=78 eV ρ=1.7 g·cm−3 Cs+Pb 能谱 Ig=81 eV ρ=2.265 g·cm−3 Co+Pb 能谱 Ig=78 eV ρ=1.7 g·cm−3 Co+Pb 能谱 Ig=81 eV ρ=2.265 g·cm−3
$\Delta $ =10 keV $\Delta $ ≠10 keV $\Delta $ =10 keV $\Delta $ ≠10 keV $\Delta $ =10 keV $\Delta $ ≠10 keV $\Delta $ =10 keV $\Delta $ ≠10 keV
10 2.677 3.277 3.57 31.3 1.010 9 1.009 4 1.003 7 1.002 7 1.001 4 0.999 9 0.993 2 0.992 3
30 3.857 4.457 5.14 38.4 1.010 9 1.009 2 1.003 7 1.002 6 1.001 4 0.999 8 0.993 2 0.992 2
50 4.610 5.21 6.15 42.5 1.010 9 1.009 0 1.003 7 1.002 5 1.001 4 0.999 7 0.993 2 0.992 1
100 5.759 6.359 7.68 48.3 1.010 9 1.008 9 1.003 7 1.002 4 1.001 4 0.999 5 0.993 2 0.992 0
300 8.306 8.906 11.07 59.4 1.010 9 1.008 6 1.003 7 1.002 1 1.001 4 0.999 3 0.993 2 0.991 8
500 9.847 10.447 13.13 65.4 1.010 9 1.008 5 1.003 7 1.002 0 1.001 4 0.999 2 0.993 2 0.991 8
1 000 12.407 13.007 16.54 74.6 1.010 9 1.008 2 1.003 7 1.001 8 1.001 4 0.999 0 0.993 2 0.991 7
3 000 17.894 18.494 23.86 91.8 1.010 9 1.007 8 1.003 7 1.001 5 1.001 4 0.998 8 0.993 2 0.991 5
5 000 21.216 21.816 28.29 101.1 1.010 9 1.007 6 1.003 7 1.001 3 1.001 4 0.998 6 0.993 2 0.991 3

2.2 光子在空腔中的作用

Fair是初级光子在空腔中沉积的能量占空腔中总沉积能量的份额,它可以用于判断空腔是否满足空腔理论的前提条件。利用CAVSPHnrc分别计算在Cs+Pb能谱、0.662 MeV、Co+Pb能谱和1.25 MeV4种不同入射光子能量时,不同体积的空腔内光子直接作用沉积的能量 $D_{{\rm{gas}}}^{\rm{p}}$ 以及相应条件时空腔中总的沉积能量 $D_{{\rm{gas}}}^{{\rm{total}}}$ 。具体方法为:创建PEGS4文件,其中石墨的截止能量AE大于入射光子能量、ECUT=512 keV、AP=PCUT=1 keV,空气截止能量为AE=ECUT=512 keV,AP=PCUT=1 keV。用CAVSPHnrc计算得到只有光子在空腔中沉积的能量 $D_{{\rm{gas}}}^{\rm{p}}$ 。然后创建新的PEGS4文件,将石墨和空气的截止能量均设置为AE=ECUT=512 keV,AP=PCUT=1 keV,计算空腔中总的沉积能量 $D_{{\rm{gas}}}^{{\rm{total}}}$ ,二者之比得到Fair,计算结果如图2,不同入射光子能量,随着空腔体积增加,Fair越来越大。可以得到,空腔体积的增大导致了空腔理论前提条件的失效。这种现象也说明了在式(1)中引入kSA的必要性。

图 2 空腔中光子作用沉积能量的份额Fair随空腔体积的变化

2.3 电离室内吸收剂量Dgas、壁效应修正kwall和空气比释动能kair

在用CAVSPHnrc使用计算Dgas时,源项有4种设置,分别是Cs+Pb能谱、单能0.662 MeV、Co+Pb能谱和单能1.25 MeV。在PEGS4文件中,石墨和空气的截止能量均为AE=ECUT=512 keV,AP=PCUT=1 keV。通常AE和ECUT越小,模拟计算的结果越准确。但此时计算时间会成倍增加。计算结果如表2。计算时设置光子输运历史数为9×108个,使计算结果的统计不确定度小于0.1%。在计算得到的输出文件中,可以得到相应条件时的壁效应修正因子kwall

表 2 不同光子能量、不同空腔体积和不同截止能量时的Dgas大小
体积/
cm3
Dgas/Gray
0.662 MeV
AE=ECUT=512 keV
0.662 MeV
AE=ECUT=521 keV
相对
偏差/%
1.25 MeV
AE=ECUT=512 keV
1.25 MeV
AE=ECUT=521 keV
相对
偏差/%
10 2.999×10−12 3.003×10−12 −0.13 5.229×10−12 5.235×10−12 −0.12
30 2.994×10−12 2.996×10−12 −0.07 5.224×10−12 5.229×10−12 −0.10
50 2.992×10−12 2.993×10−12 −0.03 5.219×10−12 5.223×10−12 −0.07
100 2.988×10−12 2.990×10−12 −0.07 5.215×10−12 5.217×10−12 −0.05
300 2.981×10−12 2.982×10−12 −0.03 5.204×10−12 5.209×10−12 −0.10
500 2.979×10−12 2.980×10−12 −0.03 5.201×10−12 5.206×10−12 −0.09
1 000 2.973×10−12 2.977×10−12 −0.13 5.194×10−12 5.198×10−12 −0.07
3 000 2.966×10−12 2.969×10−12 −0.10 5.185×10−12 5.188×10−12 −0.06
5 000 2.965×10−12 2.967×10−12 −0.07 5.179×10−12 5.184×10−12 −0.09

同样使用CAVSPHnrc计算空腔中的空气比释动能kair。设置室壁与空腔气体的材料都为空气,室壁的厚度为0.001 cm。在计算使用的PEGS4文件中,空气的AE=ECUT=3 MeV,AP=PCUT=1 keV。此时计算得到的空腔中的吸收剂量就等同于该空腔内的空气比释动能kair

2.4 质能吸收系数比 $({{\overline \mu}_ {\rm{en/}}}\rho )_{{\rm{air}},{\rm{c}}}$ ${{\overline g}_ {\rm{air}}}$

石墨和空气的质能吸收系数 ${(\overline \mu_ {\rm{en/}}}\rho )_{\rm{c}}$ ${{{(\overline \mu}_ {\rm{en/}}}\rho )_{\rm{air}}}$ 以及电子在空气中的辐射能量损失份额 $\overline g_ {\rm{air}}$ 都可以由g用户程序计算。在计算使用的PEGS4文件中,石墨和空气的截止能量均为AE=ECUT=521 keV,AP=PCUT=1 keV。设置计算的历史数为109,计算结果的统计不确定度小于0.1%。

2.5 石墨和空气的限制的阻止本领比 $(\overline{L}_ \Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$

SPRRZnrc是专门用于计算Spencer-Attix限制的阻止本领比的EGSnrc用户程序。在计算时开启光子再生功能,用于消除石墨对初级光子的衰减和散射作用,使石墨内达到带电粒子平衡[12]。在Bethe[11]公式中,材料的平均激发能I和密度效应修正的大小影响了其阻止本领的计算。ICRU37号报告中给出了Ig的推荐值为78 eV,而在更新的ICUR90[13]号报告中Ig的推荐值为81 eV。因此需要考虑不同Ig值对石墨和空气的限制的阻止本领比的影响。在计算 $(\overline{L}_ \Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$ 之前,需要创建两个PEGS4文件,一个使用石墨平均激发能Ig=78 eV,石墨密度效应修正ρ=1.7 g·cm−3,另一个使用Ig=81 eV,密度效应修正ρ=2.265 g·cm−3。此外, $\Delta $ 的选择会对阻止本领比的计算值产生影响[14]。利用表1中的 $\Delta $ 值,设置PEGS4文件中石墨和空气的截止能量均为AE=511 keV+ $\Delta $ ,ECUT=512 keV,AP=PCUT=512 keV,用于 $(\overline{L}_ \Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$ 的计算,计算结果在表1中。

2.6 Spencer-Attix理论修正因子kSA

根据上述方法,计算得到式(2)右侧的所有物理量和修正因子,利用式(2)便可以计算出在不同光子入射能量时,各空腔体积大小所对应的kSA大小。如图3图4所示,对于Co+Pb能谱和1.25 MeV能量以及Cs+Pb能谱和0.662 MeV能量的光子,随着空腔体积的增大,使用AE=511 keV+ $\Delta $ $\Delta $ 的取值与空腔体积有关,AE用于计算 $(\overline{L} _\Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$ )得到的kSA值的比使用AE=521 keV(其中 $\Delta $ =10 keV)得到的kSA值更接近于1。这说明对于不同大小的空腔,在计算 $(\overline{L}_ \Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$ 时需要考虑使用合适的 $\Delta $ 值,而不是仅将 $\Delta $ 设置为10 keV。并且从图3图4中可以看出,随着空腔体积的增大,kSA值有逐渐减小的趋势。这说明随着空腔的增大,空腔理论的计算值与真实值的偏差也逐渐增大。因此需要根据合适的kSA值以及石墨空腔电离室的实际加工难易程度来确定可以使用的最大空腔体积。

图 3 Co+Pb能谱和1.25 MeV光子能量下不同体积时的kSA

图 4 Cs+Pb能谱和0.662 MeV光子能量下不同体积时的kSA

2.7 密度效应修正和平均激发能对kSA的影响

不同的石墨的密度效应修正和石墨的平均激发能都会影响空腔气体沉积能量的计算结果,从而影响了kSA的大小。本次研究对比了使用石墨的平均激发能Ig=78 eV和石墨体密度效应修正计算的kSA与石墨的平均激发能Ig=81 eV和使用石墨颗粒密度效应修正计算的kSA。对比结果从表3中可以得到,不同的密度效应修正和平均激发能的选择对kSA的大小影响很小,相差小于0.1%。

表 3 不同入射光子能量、不同空腔体积、不同密度效应修正和平均激发能时kSA的大小
体积/
cm3
kSA
Cs+Pb 能谱
ρ=1.7 g·cm−3
Ig=78 eV
统计不
确定度/%
Cs+Pb 能谱
ρ=2.265 g·cm−3
Ig=81 eV
统计不
确定度/%
相对
偏差/%
Co+Pb 能谱
ρ=1.7 g·cm−3
Ig=78 eV
统计不
确定度/%
Co+Pb 能谱
ρ=2.265 g·cm−3
Ig=81 eV
统计不
确定度/%
相对
偏差/%
10 0.999 11 0.058 6 0.999 47 0.061 4 −0.036 0 0.999 84 0.036 6 0.999 67 0.037 0 0.000 17
30 0.998 87 0.057 6 0.999 7 0.060 4 −0.083 1 0.999 99 0.036 1 0.999 67 0.036 6 0.000 32
50 0.999 61 0.057 5 0.999 3 0.060 4 0.031 0 0.999 27 0.035 7 0.999 73 0.036 1 −0.000 46
100 0.998 77 0.057 1 0.999 37 0.060 0 −0.060 1 0.999 48 0.035 7 0.999 33 0.036 1 0.000 15
500 0.998 63 0.056 2 0.999 09 0.059 2 −0.046 1 0.999 14 0.035 2 0.999 45 0.035 7 −0.000 31
1 000 0.998 26 0.055 7 0.998 91 0.058 8 −0.065 1 0.998 86 0.035 2 0.999 49 0.035 7 −0.000 63
5 000 0.998 11 0.054 9 0.999 03 0.058 0 −0.092 2 0.998 62 0.034 8 0.999 0.035 3 −0.000 38

3 结束语

通过上述方法,探究了不同光子能量和不同空腔体积时的球型石墨电离室是否满足空腔理论的前提条件,计算对比了不同条件下 ${\rm{\gamma}} $ 射线空气比释动能的空腔理论计算值和真实值的差异,即Spencer-Attix空腔理论修正因子kSA的大小。且在相同空腔体积和入射光子能量时,本次计算的结果与Borg和Seltzer的计算结果相差小于0.2%。从图2中可以得到:当体积增加时Fair逐渐增大,此时的电离室空腔不能视为理想的Bragg-Gray空腔。但当使用合适的 $\Delta $ 值来计算 $(L _\Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$ 时,kSA的大小可以保持在0.998~1.000之间。从图3图4中可以得到,对于用于 ${\rm{\gamma}} $ 射线空气比释动能绝对测量的石墨空腔电离室, $\Delta $ 取值应该与空腔体积大小有关。并且在设计和使用较大体积的石墨空腔电离室时需要考虑合适的kSA值大小。从表3中可以得到:不同石墨的密度效应修正和平均激发能对最终kSA的影响很小,两次计算结果的相对偏差小于0.1%。因此在之后 $(L_ \Delta(E) / \rho)_{\mathrm{c}, \rm { air }}$ 的计算中可以考虑使用ICUR90号报告的Ig和密度效应修正的推荐值。本次研究方法与结果可以为较大体积的石墨空腔电离室在 ${\rm{\gamma}} $ 射线空气比释动能绝对测量中的设计和使用提供参考。

参考文献
[1]
γ射线空气比释动能计量器: JJG 2044-2010[S]. 北京: 中国标准出版社, 2011.
[2]
MA C M, NAHUM A E. Bragg-Gray theory and ion chamber dosimetry for photon beams[J]. Physics in Medicine & Biology, 1991, 36(4): 413-428.
[3]
BORG J, KAWRAKOW I, ROGERS D W, et al. Monte Carlo study of correction factors for Spencer-Attix cavity theory at photon energies at or above 100 keV[J]. Medical Physics, 2000, 27(8): 1804-1813. DOI:10.1118/1.1287054
[4]
RUSSA D J L, ROGERS D W O. Accuracy of Spencer-Attix cavity theory and calculations of fluence correction factors for the air kerma formalism[J]. Medical Physics, 2009, 36(9): 4173-4183.
[5]
COMMITTEE R T. A protocol for the determination of absorbed dose from high‐energy photon and electron beams[J]. Medical Physics, 1983, 10(6): 741-771. DOI:10.1118/1.595446
[6]
BÜERMANN L, BURNS D T. Air-kerma cavity standards[J]. Metrologia, 2009, 46(2): 24. DOI:10.1088/0026-1394/46/2/S03
[7]
ALLISYROBERTS P J, BURNS D T, KESSLER C. Summary of the BIPM.RI(I)-K1 comparison for air kerma in 60Co gamma radiation[J]. Metrologia, 2007, 44(1A): 06006. DOI:10.1088/0026-1394/44/1A/06006
[8]
ALLISYROBERTS P J, KESSLER C, BURNS D T. Summary of the BIPM.RI(I)-K5 comparison for air kerma in 137Cs gamma radiation[J]. Metrologia, 2013, 50(1A): 06001. DOI:10.1088/0026-1394/50/1A/06001
[9]
ATTIX F H. Introduction to Radiological Physics and Radiation Dosimetry[M].New Jersey:John Wiley and Sons, 1986.
[10]
ROGERS D W O, TREURNIET J. Monte Carlo calculated wall and axial non-uniformity corrections for primary standards of air kerma[R]. Nrcc Report Pirs, 1999.
[11]
ICRU. Stopping powers for electrons and positrons[R]. ICRU Report No.37. WashingtonDC: ICRU, 1984: 15-34.
[12]
RUSSA L, DANIEL J. The effect of low-energy electrons on the response of ion chambers to ionizing photon beams[J]. ProQuest Dissertations And Theses, 2009, 71(5): 3119.
[13]
ICRU. Key date for ionizing -radiation dosimetry: measurement stands and applications[R]. ICRU Report No.90. England: Oxford University Press, 2016: 15-34.
[14]
ROGERS D W, KAWRAKOW I. Monte Carlo calculated correction factors for primary standards of air kerma[J]. Medical Physics, 2003, 30(4): 521-532. DOI:10.1118/1.1563663
[15]
SELTZER S M, JR B P M. Changes in the U.S. primary standards for the air kerma from gamma-ray beams[J]. Journal of Research of the National Institute of Standards & Technology, 2003, 108(5): 359-381.