一种圆柱体探测器对点源的源峰探测效率获取方法与流程

本发明涉及辐射环境监测技术领域，具体是指一种圆柱体探测器对点源的源峰探测效率获取方法。

背景技术：

环境中γ辐射监测必须采用相应的探测器，目前多采用圆柱体或同轴型(相当于大圆柱体探测器内挖出一个小的圆柱体探测器)探测器。为估算辐射环境待测样中指定放射性核素的活度，必须获得探测器的源峰探测效率(单位活度源能被探测器记录到的特定能量γ射线全能峰计数大小)。考虑到辐射安全性及各类标准源购置综合费用高的影响，以往多采用montecarlo(蒙特卡罗，后续简称mc)软件模拟4π角均匀发散的γ射线与位置固定的探测器相互作用过程并记录入射γ射线全部能量损失在上述探测器内的概率，以此反推探测器的源峰探测效率，但计算时间过长。为此根据γ场叠加原理，可将任意形状体源网格划分成多个“微体积源”，则任意形状体源的源峰探测效率值等价于网格划分后的所有“微体积源”的源峰探测效率值之和。在网格划分足够精细的情况下，“微体积源”可用点源替代。据此，陈立等人提出数值计算与mc相结合的方法，但存在以下两个问题：①所推导的公式仅能计算圆柱体探测器顶面正上方点源的源峰探测效率值；②公式中点源的本征效率参数(圆柱体探测器对指定位置点源以不同角度发射特定能量e的γ射线的源峰探测效率)阵列获取复杂。

技术实现要素：

本发明的目的在于克服上述背景技术的不足之处，提出一种圆柱体探测器对点源的源峰探测效率获取方法，核心在于推导圆柱体探测器对任意位置点源的源峰探测效率计算公式，结合mc模拟获得全能峰线衰减系数。

本发明的目的是通过如下技术方案实现的，一种圆柱体探测器对点源的源峰探测效率获取方法：

对于一个新的圆柱体探测器，在使用本方法前必须先标定其几何结构，包括圆柱体探测器底面半径和高度、各屏蔽体顶层及侧面厚度。

利用窄束γ射线在物质中的指数衰减规律及γ场叠加原理构建任意尺寸圆柱体探测器源峰探测效率的解析方程式。

将圆柱体探测器外所有屏蔽体用真空替代，采用mc软件模拟11个指定位置点源发射能量为e的γ射线在圆柱体探测器内源峰探测效率值。利用全能峰线衰减系数μe∈(0,μ)的特点，以其范围两端作为初始值，代入上述解析方程式计算上述mc模拟模型中各位置点源发射能量为e的γ射线在圆柱体探测器内源峰探测效率值和并计算其与模拟值间的加权偏差平方和ζa和ζb。当ζa≥ζb时，取a＝(a+b)/2,b＝b；反之，取a＝a,b＝(a+b)/2；重复步骤上述过程直至|a-b|≤μ/10⁵结束。此时如果ζa≥ζb，则μe＝b；反之，取μe＝a。

将探测器几何结构参数、点源位置、μe和其它屏蔽材料的总线衰减系数代入上述解析方程式，采用梯形法即可获得上述圆柱体探测器对任意位置点源发射的任意能量γ射线的源峰探测效率值。

与现有技术相比，本发明具有如下优点：

1、利用不同角度下γ射线在裸探测器内的计算值与mc模拟值间的加权偏差平方和最小的原则，获得γ射线在探测器介质内的全能峰衰减系数，解决了以往计算中需点源的本征效率参数阵列，方法更加简单，且获取参数费时更短。

2、推导了圆柱体探测器圆柱面外点源的探测效率计算方法，将该方法推广至能获得探测器外任意位置点源的源峰探测效率值，扩展了以往方案的应用范围。

附图说明

图1为本发明实施例1的处理流程图；

图2为scionix公司cebr3圆柱体探测器中心轴竖切面示意图及点源位置分区；

图3为mc模拟探测器结构及点源布点位置；

图4为加权偏差平方随迭代次数的变化；

图5为本发明方法计算值与实测值间相对偏差分布图。

具体实施方式

下面结合附图和具体实施例对本发明做进一步的描述，然而本发明的范围并不限于下述实施例。

实施例1：下面通过附图和具体实施方式对本发明作更为详细的描述。

实施例采用scionix公司cebr3圆柱体探测器，经过几何结构标定后(中心轴竖切面示意图如图2所示)参数分别为：圆柱体探测器底面半径r＝1.27cm、高l＝2.54cm；由内而外各屏蔽层的尺寸为：alo3反射层顶面垂直厚度d1＝0.05cm、侧面轴向长度c1＝0.05cm，聚乙烯减振层顶面垂直厚度d2＝0.83cm、侧面轴向长度c2＝0.355cm，al壳顶面垂直厚度d3＝0.08cm、侧面轴向长度c3＝0.125cm。cs-137源特征γ射线的能量e＝661.657kev，查找γ射线与物质相互作用截面库可知alo3反射层、聚乙烯减振层和al壳对该能量γ射线的总衰减系数分别为和点源外包围空气的总衰减系数分别为可以看出点源到cebr3晶体间的屏蔽层的个数m＝4。

(1)利用窄束γ射线在物质中的指数衰减规律及γ场叠加原理构建任意尺寸圆柱体探测器源峰探测效率的解析方程式：

①当点源p位于圆柱体探测器顶面正上方时(如图2中x区)

②当点源p位于圆柱体探测器顶面侧上方时(如图2中y区)

对积分区间内的任意值，如果那么

反之，

③当点源p位于圆柱体探测器圆柱面外(如图2中z区)

(2)全能峰线衰减系数μe的获取

①将圆柱体探测器外所有屏蔽体用真空替代，采用mc软件模拟图3中t＝11个位置点源发射能量为e的γ射线在圆柱体探测器内源峰探测效率值(k＝1,2,…,t)；

②因μe的最大值为能量e的γ射线在圆柱体探测器内的总线衰减系数μ(本实施例取值为0.387011cm^-1)，则假设a＝0,b＝μ；

③取μe＝a与μe＝b代入(1)中解析方程式计算上述mc模拟模型中各位置点源发射能量为e的γ射线在圆柱体探测器内源峰探测效率值和计算其与模拟值间的加权偏差平方和

④当ζa≥ζb时，取a＝(a+b)/2,b＝b；反之，取a＝a,b＝(a+b)/2；重复步骤③直至|a-b|≤μ/10⁵结束。此时如果ζa≥ζb，则μe＝b；反之，取μe＝a。

上述③中采用mc软件模拟，其mc软件可能采用通用模拟软件也可采用自编软件，本实施例中mc模拟的方法具体如下：

基于vc++平台编写，软件能模拟1kev-3mev的γ射线在圆柱体探测器内的全能峰探测效率值。截面库选用xcom，包括γ射线与物质发生光电效应、康普顿散射和电子对效应。计算过程中采用的方差减小技术为：当γ射线与探测器介质外其它物质发生相互作用时则γ射线历史终止；利用指数变换提高γ射线在探测器介质中的权重。并利用gpu并行计算技术进一步缩减模拟时间。

图4给出每次迭代min(ζa,ζb)随迭代次数的变化规律，本实施例计算结果μe＝9.675266×10^-2cm^-1。

(3)将探测器几何结构参数、点源位置、μe和其它屏蔽材料的总线衰减系数代入(1)中解析方程式，采用梯形法即可获得上述圆柱体探测器对任意位置点源发射的任意能量γ射线的源峰探测效率值。

在上述cebr3圆柱体探测器外布置110个验证点(以探测器中心轴为z轴间角度夹角0°-90°，以10°为间隔；测点与cebr3圆柱体中心直线距离分别为5cm、10cm、15cm、20cm、30cm、50cm、70cm、100cm、135cm、170cm和200cm)测量cs-137点源源峰探测效率值，本发明方法计算值与实测值间相对偏差分布图如图5所示，可以结果在±8％范围内符合。

本发明所涉及的所有参数符号的具体含义为：

ηe：指定位置点源发射能量为e的γ射线在圆柱体探测器内的源峰探测效率值；

ε：点源发射能量为e的γ射线的几率；

r：圆柱体探测器底面半径；

l：圆柱体探测器高度；

μe：能量e的γ射线在圆柱体探测器内的全能峰线衰减系数；

μ：能量e的γ射线在圆柱体探测器内的总线衰减系数，查找截面库获得；

di：从圆柱体探测器向外第i屏蔽体顶层垂直厚度；

ci：从圆柱体探测器向外第i屏蔽体侧面径向厚度；

能量e的γ射线在第i屏蔽体内的总线衰减系数，查找截面库获得；

m：圆柱体探测器前屏蔽体总数目；

h：点源p与其在圆柱体探测器顶面上的投影q之间直线距离；

w：点源p在圆柱体探测器顶面上的投影q与圆柱体轴线间的垂直距离；

θ：点源p出射γ射线与pq线间的夹角；

φ：点源p出射γ射线在圆柱体探测器顶面上投影线与oq线间的夹角；

能量e的γ射线从顶面入射进入圆柱体探测器；

能量e的γ射线从侧面入射进入圆柱体探测器；

积分中间变量；

将圆柱体探测器外所有屏蔽体用真空替代时，mc软件模拟图5中第k个位置点源发射能量为e的γ射线在圆柱体探测器内源峰探测效率值；

ζi：将圆柱体探测器外所有屏蔽体用真空替代时，第i种μe取值计算值与mc模拟值间偏差平方和；

t：真空无屏蔽环境下点源模拟个数；

a,b：μe不同取值；

xcom：一个γ射线与任意物质相互作用截面库，由美国国家标准与技术研究院m.j.berger等人编制，详见https://physics.nist.gov/physrefdata/xcom/text/intro.html；

gpu：全称graphicsprocessingunit，图形处理器。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本实用新型的保护范围之内。