CN109839657B - 一种辐射剂量仪能量响应补偿分析的方法 - Google Patents
一种辐射剂量仪能量响应补偿分析的方法 Download PDFInfo
- Publication number
- CN109839657B CN109839657B CN201711217984.9A CN201711217984A CN109839657B CN 109839657 B CN109839657 B CN 109839657B CN 201711217984 A CN201711217984 A CN 201711217984A CN 109839657 B CN109839657 B CN 109839657B
- Authority
- CN
- China
- Prior art keywords
- energy
- compensation
- radiation dosimeter
- value
- materials
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Measurement Of Radiation (AREA)
Abstract
本发明涉及一种辐射剂量仪能量响应补偿分析的方法,属于辐射防护领域。本发明所述的方法,包括以下步骤:1)通过实验测量剂量仪在不同能量E1、E2……En光子射线照射下的原始响应K1、K2……Kn;2)列举所有可能的补偿材料并对这些材料在不同能量Ei的光子照射下的线衰减系数μi值进行计算;3)列举出全部可能的补偿材料的空隙率p和厚度t的组合,并代入方程组求解值;4)用户利用计算机程序设置不同的选择条件,从计算结果中选择出符合用户需要的能量补偿方案。本发明可应用于辐射剂量仪的能量补偿,保证剂量监测准确性,具有实用性。
Description
技术领域
本发明属于辐射防护领域,具体涉及一种辐射剂量仪能量响应补偿分析的方法。
背景技术
辐射剂量仪被广泛应用于核电站、工业无损探伤、职业病预防、部队等存在放射性的工作场所的个人剂量监测,但是未经能量补偿的辐射剂量仪存在不同能量段响应不一致,严重影响测量准确度的问题。从上个世纪五六十年代开始,辐射剂量仪的能量响应补偿就是辐射剂量测量领域的一项重要内容,如何确定能量补偿方案是很多辐射测量设备研发人员需要花费很多心血解决的问题。已有的寻找优化的能量补偿方案的方法除通过大量的实验直接获取数据,从中挑选效果最好的方案外,还包括蒙特卡罗软件分析、数学理论计算等方法。
对于蒙特卡罗方法,首先人为给出一个能量补偿方案,然后利用蒙特卡罗方法模拟计算剂量仪的响应情况再根据模拟结果由人做出决策,调整参数后再进行模拟,不断重复上述过程,直到得到合适的结果为止。该方法不仅费时费力,而且人类的思维容易陷入局部最优化,如果实际上存在两种或两种以上的优化的能量补偿方案,该方法往往仅能找到其中一种。
除蒙特卡罗方法外,能量补偿方案的确定还可以使用如下数学理论计算方法:根据能量补偿理论,在射线能量为Ei时,剂量仪在无任何补偿措施时的原始响应记为Ki(Ei),选取某种补偿材料其线衰减系数为μi,设补偿材料的空隙率为p,材料厚度为t,则剂量仪有补偿材料部分的能量响应为空隙部分的能量响应为pKi(Ei)。经过补偿后,在射线能量为Ei时剂量仪能量响应可以表示为:
上述方程是在能量补偿材料确定以后,即ui确定以后,才可以求解,并不能解决材料选择的问题。
而且,在实际应用中,即便能量补偿材料已经确定,该方法也遇到一些问题。例如,由于建立数学方程时,“δi足够小”且“越小越好”的说法在数学上是不太严谨的。由于由i个δ值组成,利用数学公式不能严格描述这一具有模糊性的条件。现在通常是将条件设定为的值最小,认为此条件等同于“最小”,然后看此条件下是否所有的值均足够小。实际中,可能存在一种能量补偿方法,某一个能量点的较大,甚至不能满足辐射监测相关标准要求,但其它能量点的均很小,而且仍然是最小的;而另一种能量补偿方法中,各点的均满足标准要求但都不太低,这时的值较高。此时利用数学方程方法只能求得最小的补偿方法,可能会出现不能满足相关标准要求的“最佳”的能量补偿方法。而实际上,第二种情况是满足标准要求的,但利用数学方法却会将该方案遗漏。
发明内容
针对现有技术中存在的缺陷,本发明的目的是提供一种辐射剂量仪能量响应补偿分析的方法。该方法利用计算机软件和数学理论计算方法相结合,能够对所有可能的能量补偿方案的补偿结果进行计算,然后从中选择优化的能量补偿方案,保证剂量监测准确性。
为达到以上目的,本发明采用的技术方案是:一种辐射剂量仪能量响应补偿分析的方法,包括以下步骤:
1)测量剂量仪在不同能量E1、E2……En光子射线照射下的原始响应K1、K2……Kn;
2)列举所有可能的补偿材料并对这些材料在不同能量Ei的光子照射下的线衰减系数μi值进行计算;
3)列举出全部可能的补偿材料的空隙率p和厚度t的组合,并代入方程组(2)求解值;
方程组(2)中μ1、μ2……μn表示补偿材料在不同能量处的线衰减系数,δi指能量为Ei的射线入射时对应的δ值;
4)设置不同的选择条件,从计算结果中选择出符合用户需要的能量补偿方案。
进一步,所述步骤1)的方法包括:输入剂量计原始能量响应,剂量计探头未经补偿时,测量裸探头在不同能量E1、E2……En光子照射下的探测计数值,并归一化到某一能量点。
更进一步,所述的能量点是137Cs的γ射线能量662keV。
进一步,步骤2)中,所述的补偿材料有锡、铅、铜及其组合。
进一步,步骤2)中,以锡、铅、铜作为基础材料,每种材料以1%为变化步长,则有“1+2+3+……+101=5151”种选择,即第一种材料由0%至100%有101种选择,第一种材料确定后,第二种材料在不大于“1-第一种材料所占的百分数”的范围内选择,前两种材料及其比例确定后,第三种材料即严格确定。
进一步,步骤3)中,采用穷举法,列举出所有可能的空隙率p有n种,补偿材料的厚度t为m种,将全部可能的空隙率p和厚度t进行组合,则采用穷举法的组合数为:
A=C1 n·C1 m (3)
进一步,所述的补偿材料空隙率在0-1之间变化,空隙率的最小变化单位是0.1%,其中空隙率为1时补偿无意义,所以空隙率的变化有1000种可能。。
进一步,所述的补偿材料厚度的最大值设为5mm,以0.01mm为最小变化单位,其中厚度为0时补偿无意义,所以可列举500种厚度。
进一步,步骤4)中,对于每一种能量补偿方案,得到一组δi值,这一组δi值中绝对值最大的那个记为|δi|max,选择|δi|max最小的为该种能量补偿方案。
本发明的效果在于:采用本发明所述的方法,利用计算机软件和数学理论计算方法相结合,可以对所有可能的能量补偿方案的补偿结果进行计算,然后从中选择优化的能量补偿方案,保证剂量监测准确性。
附图说明
图1是本发明具体实施方式中所述一种辐射剂量仪能量响应补偿分析的方法的流程图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步描述。
本发明提出了一种新的方法,利用计算机程序采用穷举法对所有可能的方案进行计算,并列出所有满足标准要求的能量补偿方案,为剂量当量的准确测量提供支持。
一种辐射剂量仪能量响应补偿分析的方法,包括以下步骤:
1)通过实验测量剂量仪在不同能量E1、E2……En光子射线照射下的原始响应K1、K2……Kn;
2)列举所有可能的补偿材料并对这些材料在不同能量Ei的光子照射下的线衰减系数μi值进行计算;
3)列举出全部可能的补偿材料的空隙率p和厚度t的组合,并代入方程组(2)求解值;
方程组(2)中μ1、μ2……μn表示补偿材料在不同能量处的线衰减系数,δi指能量为Ei的射线入射时对应的δ值;
4)用户利用计算机程序设置不同的选择条件,从计算结果中选择出符合用户需要的能量补偿方案。
本实施例中,剂量仪能量补偿方案计算机程序实现方法的流程图如图1所示,该方法的具体步骤是:
步骤一、输入剂量计原始能量响应,剂量计探头未经补偿时,裸探头在不同能量光子照射下的探测计数值,并归一化到某一能量点,如本实施例中的137Cs的γ射线能量662keV。
步骤二、确定补偿材料及其线衰减系数,分析剂量计的归一化原始能量响应曲线,选择合适的补偿材料,主要采用的补偿材料有锡(Sn)、铅(Pb)、铜(Cu)等线衰减系数较大的材料,表1给出了三种材料锡、铅、铜的线衰减系数。确定补偿材料后即可计算相应的线衰减系数。
表1不同材料线衰减系数
步骤三、补偿材料的选择组合。以锡(Sn)、铅(Pb)、铜(Cu)作为基础材料,每种材料以1%为变化步长,则有“1+2+3+……+101=5151”种选择(第一种材料由0%至100%有101种选择,第一种材料确定后,第二种材料可在不大于“1-第一种材料所占的百分数”的范围内选择,前两种材料及其比例确定后,第三种材料即严格确定)
步骤四、补偿材料空隙率p和厚度t的选择组合。采用穷举法,列举出所有可能的空隙率p有n种,补偿材料的厚度t为m种。将全部可能的空隙率p和厚度t进行组合,则采用穷举法的组合数为:
A=C1 n·C1 m (3)
补偿材料空隙率显然只能在0~1之间变化,本方法中空隙率的最小变化单位是0.1%,其中空隙率为1时补偿无意义,所以空隙率的变化有1000种可能;考虑剂量仪内部空间限制,补偿材料厚度的最大可能值设为5mm,以0.01mm为最小变化单位,其中厚度为0时补偿无意义,所以可列举500种厚度。则公式(3)中,n=1000,m=500,则空隙率p和厚度t的组合数量为500000种。如果单一材料的补偿无法达到最佳效果,可能会采用两种材料,甚至多种材料,这样还涉及不同材料的比例,则空隙率p和厚度t的组合将会更多。将列出全部的空隙率p和厚度t可能组合,带入公式(4)求解δi。
步骤五、为选择最佳补偿材料、最佳空隙率和厚度,设置了两种计算方案,分别是:
1、|δi|max最小,即对于每一种能量补偿方案,可得到一组δi值(δi指能量为Ei的射线入射时对应的δ值),这一组δi值中绝对值最大的那个记为|δi|max,哪一种能量补偿方案的|δi|max最小,即选择该种能量补偿方案。
通过两种方案均可以给出相对应的补偿材料空隙率和厚度的最佳组合,而且可以给出满足相关标准的空隙率和厚度的范围,方便补偿实现。
将剂量仪未经能量补偿的原始归一化响应值带入程序进行计算,确定最佳的空隙率和厚度的组合。
例如:
某GM计数管在窄谱X射线和137Csγ射线照射下的原始归一化响应如表1所示.
表1剂量仪能量响应
能量 | 48keV | 65keV | 83keV | 100keV | 118keV | 164keV | 662keV |
原始归一化响应 | 3.70 | 3.93 | 3.36 | 2.90 | 2.36 | 1.56 | 1.00 |
利用以上计算过程,按照|δi|max值最小的方法,可搜索到|δi|max=0.07为|δi|max的最小值,此时的材料为锡:铅:铜=94:6:0,对应的空隙率p=0.265,厚度t=1.7mm。此时的能量响应计算结果如表2所示。
表2剂量仪补偿后理论计算的能量响应
能量 | 48keV | 65keV | 83keV | 100keV | 118keV | 164keV | 662keV |
原始归一化响应 | 3.70 | 3.93 | 3.36 | 2.90 | 2.36 | 1.56 | 1.00 |
理论计算响应 | 0.98 | 1.04 | 0.98 | 1.06 | 1.00 | 1.03 | 0.93 |
由于上述方案完全可以满足相关标准的要求(即|δi|max≤30%),因此选定采用锡:铅:铜=94:6:0的材料,空隙率p=0.265,厚度t=1.7mm作为能量补偿材料。如果两种方法筛选出来的结果不一致,也可根据实际需要具体选择采用何种方法。
本领域技术人员应该明白,本发明所述的方法并不限于具体实施方式中所述的实施例,上面的具体描述只是为了解释本发明的目的,并非用于限制本发明。本领域技术人员根据本发明的技术方案得出其他的实施方式,同样属于本发明的技术创新范围,本发明的保护范围由权利要求及其等同物限定。
Claims (10)
2.如权利要求1所述的一种辐射剂量仪能量响应补偿分析的方法,其特征是,所述步骤1)的方法包括:输入剂量计原始能量响应,剂量计探头未经补偿时,测量裸探头在不同能量E1、E2……En光子照射下的探测计数值,并归一化到某一能量点。
3.如权利要求2所述的一种辐射剂量仪能量响应补偿分析的方法,其特征是,所述的能量点是137Cs的γ射线能量662keV。
4.如权利要求1至3任一项所述的一种辐射剂量仪能量响应补偿分析的方法,其特征是:步骤2)中,所述的补偿材料有锡、铅、铜及其组合。
5.如权利要求4所述的一种辐射剂量仪能量响应补偿分析的方法,其特征是:步骤2)中,以锡、铅、铜作为补偿材料,每种材料以1%为变化步长,则有“1+2+3+……+101=5151”种选择,即第一种材料由0%至100%有101种选择,第一种材料确定后,第二种材料在不大于“1-第一种材料所占的百分数”的范围内选择,前两种材料及其比例确定后,第三种材料即严格确定。
6.如权利要求1所述的一种辐射剂量仪能量响应补偿分析的方法,其特征在于:步骤3)中,采用穷举法,列举出所有可能的空隙率p有n种,补偿材料的厚度t为m种,将全部可能的空隙率p和厚度t进行组合,则采用穷举法的组合数为:
A=C1 n·C1 m (3)。
7.如权利要求6所述的一种辐射剂量仪能量响应补偿分析的方法,其特征在于:所述的补偿材料空隙率在0-1之间变化,空隙率的最小变化单位是0.1%,其中空隙率为1时补偿无意义,所以空隙率的变化有1000种可能。
8.如权利要求6或7所述的一种辐射剂量仪能量响应补偿分析的方法,其特征在于:所述的补偿材料厚度的最大值设为5mm,以0.01mm为最小变化单位,其中厚度为0时补偿无意义,所以可列举500种厚度。
9.如权利要求8所述的一种辐射剂量仪能量响应补偿分析的方法,其特征在于:步骤4)中,对于每一种能量补偿方案,得到一组δi值,这一组δi值中绝对值最大的那个记为|δi|max,选择|δi|max最小的为该种能量补偿方案。
10.如权利要求8所述的一种辐射剂量仪能量响应补偿分析的方法,其特征在于:步骤4)中,或者计算每一种补偿方案得到的一组δi值的平方和∑δi 2,选择出其中的最小值作为补偿方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711217984.9A CN109839657B (zh) | 2017-11-28 | 2017-11-28 | 一种辐射剂量仪能量响应补偿分析的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711217984.9A CN109839657B (zh) | 2017-11-28 | 2017-11-28 | 一种辐射剂量仪能量响应补偿分析的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109839657A CN109839657A (zh) | 2019-06-04 |
CN109839657B true CN109839657B (zh) | 2022-05-20 |
Family
ID=66881149
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711217984.9A Active CN109839657B (zh) | 2017-11-28 | 2017-11-28 | 一种辐射剂量仪能量响应补偿分析的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109839657B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1204403A (zh) * | 1995-10-31 | 1999-01-06 | 比奥特雷斯公司 | 超低本底的多光子探测器 |
DE19827729A1 (de) * | 1998-06-22 | 2000-01-20 | Hans Schoefer | Energiekompensierte Festkörperdosismeßsonde zur gleichzeitigen Messung der Dosis, Flächendosis- und des Längendosisproduktes |
JP2008256630A (ja) * | 2007-04-09 | 2008-10-23 | Fuji Electric Systems Co Ltd | エネルギー補償型シンチレーション式光子線量計 |
CN101738626A (zh) * | 2008-11-14 | 2010-06-16 | 中国辐射防护研究院 | γ剂量计补偿过滤器的厚度和空隙率的确定方法 |
CN101796430A (zh) * | 2008-06-24 | 2010-08-04 | 富士电机系统株式会社 | 中子剂量计 |
CN101937090A (zh) * | 2010-08-12 | 2011-01-05 | 上海新漫传感技术研究发展有限公司 | 一种高灵敏宽量程X-γ周围剂量当量率仪探头 |
CN201945690U (zh) * | 2010-11-30 | 2011-08-24 | 中国辐射防护研究院 | 可对受照过程进行精细分析的个人剂量测量装置 |
CN102176047A (zh) * | 2011-03-04 | 2011-09-07 | 中国原子能科学研究院 | 一种具有能量补偿的中子个人剂量计 |
EP2420861A1 (en) * | 2010-08-20 | 2012-02-22 | Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO | A radiation dose meter for measuring radiation dose in an external magnetic field |
CN103675885A (zh) * | 2012-08-28 | 2014-03-26 | 柯尼卡美能达株式会社 | 闪烁体板和放射线检测面板 |
CN105841750A (zh) * | 2016-06-13 | 2016-08-10 | 成都新核泰科科技有限公司 | 放射源应急探测机器人 |
JPWO2016063391A1 (ja) * | 2014-10-23 | 2017-04-27 | 三菱電機株式会社 | 線量率測定装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080027974A1 (en) * | 2006-07-24 | 2008-01-31 | Collins Donald W | Intensity modulated radiation therapy filtration apparatus, system and method |
-
2017
- 2017-11-28 CN CN201711217984.9A patent/CN109839657B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1204403A (zh) * | 1995-10-31 | 1999-01-06 | 比奥特雷斯公司 | 超低本底的多光子探测器 |
DE19827729A1 (de) * | 1998-06-22 | 2000-01-20 | Hans Schoefer | Energiekompensierte Festkörperdosismeßsonde zur gleichzeitigen Messung der Dosis, Flächendosis- und des Längendosisproduktes |
JP2008256630A (ja) * | 2007-04-09 | 2008-10-23 | Fuji Electric Systems Co Ltd | エネルギー補償型シンチレーション式光子線量計 |
CN101796430A (zh) * | 2008-06-24 | 2010-08-04 | 富士电机系统株式会社 | 中子剂量计 |
CN101738626A (zh) * | 2008-11-14 | 2010-06-16 | 中国辐射防护研究院 | γ剂量计补偿过滤器的厚度和空隙率的确定方法 |
CN101937090A (zh) * | 2010-08-12 | 2011-01-05 | 上海新漫传感技术研究发展有限公司 | 一种高灵敏宽量程X-γ周围剂量当量率仪探头 |
EP2420861A1 (en) * | 2010-08-20 | 2012-02-22 | Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO | A radiation dose meter for measuring radiation dose in an external magnetic field |
CN201945690U (zh) * | 2010-11-30 | 2011-08-24 | 中国辐射防护研究院 | 可对受照过程进行精细分析的个人剂量测量装置 |
CN102176047A (zh) * | 2011-03-04 | 2011-09-07 | 中国原子能科学研究院 | 一种具有能量补偿的中子个人剂量计 |
CN103675885A (zh) * | 2012-08-28 | 2014-03-26 | 柯尼卡美能达株式会社 | 闪烁体板和放射线检测面板 |
JPWO2016063391A1 (ja) * | 2014-10-23 | 2017-04-27 | 三菱電機株式会社 | 線量率測定装置 |
CN105841750A (zh) * | 2016-06-13 | 2016-08-10 | 成都新核泰科科技有限公司 | 放射源应急探测机器人 |
Non-Patent Citations (1)
Title |
---|
张艳婷 等."环境级能量补偿实验方法的测试".《 第十八届全国核电子学与核探测技术学术年会》.2016,第288-291页. * |
Also Published As
Publication number | Publication date |
---|---|
CN109839657A (zh) | 2019-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kron et al. | Dosimetry of ionising radiation in modern radiation oncology | |
Ackermann et al. | Anisotropies in the diffuse gamma-ray background measured by the Fermi LAT | |
Ay et al. | Assessment of different computational models for generation of x‐ray spectra in diagnostic radiology and mammography | |
Shore et al. | Recent epidemiologic studies and the linear no-threshold model for radiation protection—considerations regarding NCRP commentary 27 | |
Hurier et al. | Combined analysis of galaxy cluster number count, thermal Sunyaev-Zel’dovich power spectrum, and bispectrum | |
Thierry-Chef et al. | Dose estimation for a study of nuclear workers in France, the United Kingdom and the United States of America: methods for the International Nuclear Workers Study (INWORKS) | |
Simon et al. | Estimating historical radiation doses to a cohort of US radiologic technologists | |
US8666711B2 (en) | Radiation analysis system and method | |
Fedon et al. | Internal breast dosimetry in mammography: experimental methods and Monte Carlo validation with a monoenergetic x‐ray beam | |
Granero et al. | Equivalent phantom sizes and shapes for brachytherapy dosimetric studies of and | |
CN114660096B (zh) | 一种材料热中子屏蔽性能的测试方法 | |
Romero‐Expósito et al. | Experimental evaluation of neutron dose in radiotherapy patients: Which dose? | |
CN109471999B (zh) | 一种非均匀源项分布的γ辐射场数据修正计算方法及系统 | |
Gressier et al. | Characterisation of the IRSN CANEL/T400 facility producing realistic neutron fields for calibration and test purposes | |
Salehi et al. | X-ray spectra and quality parameters from Monte Carlo simulation and analytical filters | |
CN109839657B (zh) | 一种辐射剂量仪能量响应补偿分析的方法 | |
Santos et al. | Evaluation of conversion coefficients relating air-kerma to H*(10) using primary and transmitted x-ray spectra in the diagnostic radiology energy range | |
Glushkov et al. | Estimation of primary cosmic ray energy registered at the EAS Yakutsk array | |
Takahara et al. | Statistical characterization of radiation doses from external exposures and relevant contributors in Fukushima prefecture | |
Park et al. | Counting efficiencies determined by Monte Carlo methods for in vivo measurement of 131I activity in thyroid | |
Van Den Elzen et al. | Alanine response to low energy synchrotron x-ray radiation | |
CN114740516B (zh) | 能谱-剂量的测量方法及装置 | |
Fonseca et al. | MaMP and FeMP: computational mesh phantoms applied for studying the variation of WBC efficiency using a NaI (Tl) detector | |
Kojo et al. | Radon exposure concentrations in Finnish workplaces | |
Bontempi et al. | Semiempirical simulation of x‐ray detectors for imaging applications |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |