CN112051601B - 一种基于虚拟源原理的源边界参数蒙特卡罗反演技术 - Google Patents
一种基于虚拟源原理的源边界参数蒙特卡罗反演技术 Download PDFInfo
- Publication number
- CN112051601B CN112051601B CN202010685178.XA CN202010685178A CN112051601B CN 112051601 B CN112051601 B CN 112051601B CN 202010685178 A CN202010685178 A CN 202010685178A CN 112051601 B CN112051601 B CN 112051601B
- Authority
- CN
- China
- Prior art keywords
- source
- virtual point
- hot zone
- steel frame
- peak
- 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.)
- Expired - Fee Related
Links
- 238000005516 engineering process Methods 0.000 title abstract description 19
- 238000001514 detection method Methods 0.000 claims abstract description 64
- 229910000831 Steel Inorganic materials 0.000 claims abstract description 48
- 239000010959 steel Substances 0.000 claims abstract description 48
- 238000004364 calculation method Methods 0.000 claims abstract description 42
- 229910001220 stainless steel Inorganic materials 0.000 claims abstract description 31
- 239000010935 stainless steel Substances 0.000 claims abstract description 31
- 239000004576 sand Substances 0.000 claims abstract description 27
- 238000000034 method Methods 0.000 claims description 29
- 239000002689 soil Substances 0.000 claims description 23
- 230000000694 effects Effects 0.000 claims description 22
- 238000001228 spectrum Methods 0.000 claims description 16
- 230000005855 radiation Effects 0.000 claims description 12
- 238000004088 simulation Methods 0.000 claims description 10
- 229910052778 Plutonium Inorganic materials 0.000 claims description 4
- OYEHPCDNVJXUIW-UHFFFAOYSA-N plutonium atom Chemical compound [Pu] OYEHPCDNVJXUIW-UHFFFAOYSA-N 0.000 claims description 4
- 239000010410 layer Substances 0.000 description 45
- 238000005259 measurement Methods 0.000 description 22
- 230000002285 radioactive effect Effects 0.000 description 15
- 239000013078 crystal Substances 0.000 description 13
- 238000012512 characterization method Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000000342 Monte Carlo simulation Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 5
- 238000011109 contamination Methods 0.000 description 4
- 239000002699 waste material Substances 0.000 description 4
- 229910052790 beryllium Inorganic materials 0.000 description 3
- ATBAMAFKBVZNFJ-UHFFFAOYSA-N beryllium atom Chemical compound [Be] ATBAMAFKBVZNFJ-UHFFFAOYSA-N 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011038 discontinuous diafiltration by volume reduction Methods 0.000 description 2
- 229910052732 germanium Inorganic materials 0.000 description 2
- GNPVGFCGXDBREM-UHFFFAOYSA-N germanium atom Chemical compound [Ge] GNPVGFCGXDBREM-UHFFFAOYSA-N 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- SAZUGELZHZOXHB-UHFFFAOYSA-N acecarbromal Chemical compound CCC(Br)(CC)C(=O)NC(=O)NC(C)=O SAZUGELZHZOXHB-UHFFFAOYSA-N 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005025 nuclear technology Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000003904 radioactive pollution Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000003900 soil pollution Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开一种基于虚拟源原理的源边界参数蒙特卡罗反演技术,包括不锈钢架平台、探测器、上层沙土衰减层、热区、深层干扰源和干扰源虚拟点源,所述热区上方设置有沙土衰减层,热区下方设置有深层干扰源,反演计算时需要将热区及深层干扰源虚拟成虚拟点源,虚拟点源位置在探测器对称轴线上,上层沙土衰减层和热区厚度及深层干扰源厚度可调,其最终厚度由建立峰谷比及虚拟点探测效率方程和均方偏差反演计算确定,小矩形钢架两条长边长度要长过大矩形钢架长边,且搭在大矩形钢架上,并固定,设计成滑轨结构,横向搭载两个平行不锈钢滑轨,滑轨与小矩形长边使用滑轨式接触,使两平行滑轨能够竖向滑动移动。
Description
技术领域
本发明专利涉及核技术应用领域,属于射性污染核素的测量及移动平台装置系统,具体涉及测量对象为沙土污染源,特别是放射性核素源边界参数蒙特卡罗反演技术。
背景技术
对于“热粒子”、“放射性汇集点”、“放射性汇集区”的测量,就地γ谱仪搜索扫描测量及伽玛相机成像技术,只能给出污染源的大概位置,而对于源边界参数的确定目前研究较少,尤其是无法给出核素深度方向的参数;通过出壳γ谱反推核弹头结构材料厚度是一件比较复杂和难度较大的工作,需要的参数多,解的不确定性大。本发明反向使用蒙特卡罗模拟表征晶体参数的技术,结合虚拟点源刻度方法等,根据测量能谱分析,就可以确定放射性污染区域深度方向分布参数,同时对反解核弹头惰层厚度具有参考意义。
发明内容
本发明的目的在于提供一种基于虚拟源原理的源边界参数蒙特卡罗反演技术,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:
一种基于虚拟源原理的源边界参数蒙特卡罗反演方法,包括不锈钢架平台、探测器、上层沙土衰减层、热区、热区虚拟点源、下层沙土衰减层、深层干扰源和干扰源虚拟点源,所述热区上方设置有沙土衰减层,热区下方设置有深层干扰源,热区和深层干扰源之间设置有下层沙土衰减层,上层沙土衰减层和下层沙土衰减层、热区及深层干扰源密度和成分一致,反演计算时需要将热区及深层干扰源虚拟成热区虚拟点源和干扰源虚拟点源,热区虚拟点源和干扰源虚拟点源位置在探测器对称轴线上,上层沙土衰减层和热区的厚度和深层干扰源厚度可调,其最终厚度由反演计算确定的虚拟点位置确定得到。
作为本发明进一步的方案:包括以下步骤:
S1:等效虚拟点源
为简化模型,将热区和深层干扰源分别等效虚拟成其虚拟点源,再将两个虚拟点源等效成一个等效虚拟点源,该点源包含所有源的辐射信息;对于钚污染土,241Am和239Pu的单位体积活度比一般为常数,根据这一比值可定出它们的能峰241Am的59.54keV(分支比0.359,26.4keV分支比0.024)和239Pu的129keV(分支比0.000062)射线等效虚拟点源探测效率(也即等效虚拟点在对称轴上的位置);
S2:建立峰谷比及虚拟点探测效率方程
根据等效虚拟点射线能峰探测效率和实测能谱的峰谷比参数,再结合MCNP程序建立测量对象源模型,进行理论反演计算源的边界参数,其反演主要依据如下三个方程:
59.54keV峰谷比方程:
59.54keV虚拟点探测效率方程:
ε1(h1)p1+ε2(h2)(1-p1)=ε0(h0) (2)
26.4keV或129keV虚拟点探测效率方程:
ε1(h1)*p1+ε2(h2)*(1-p1)=ε0(h0)* (3)
v表示谷区计数,n表示峰计数,p表示权重因子,h表示虚拟点位置(虚拟点位置一般在对称轴线上变化),0表示实测谱值,1表示虚拟点位置1,2表示虚拟点位置2;
S3:计算均方偏差σ(X)
等式(3)右边是定值E(Xi)(实验值,又叫数学期望值),其中式(2)和式(3)等式右边的虚拟点源位置h0处的探测效率值是通过核素活度比值定出,峰谷比通过能谱获得,左边是蒙卡模拟结果,使用蒙卡模拟等式左边探测效率及峰谷比等(不同衰减层厚、不同源厚和不同权重,两两组合),再使用最小二乘法计算均方偏差σ(X),如式(4)所示,找出均方偏差最小的组合即可,
作为本发明再进一步的方案:所述不锈钢架平台主要由钢管连接成网状结构,网状结构为一个大矩形钢架套一个小矩形钢架,对角线也用钢管连接,中间长度方向加一根钢管,不锈钢架平台大小设计为5.0×6.0m,小矩形钢架两条长边长度大于大矩形钢架长边,且搭在大矩形钢架上,并进一步固定设计成滑轨结构,横向搭载两个平行不锈钢滑轨,平行不锈钢滑轨与小矩形钢架长边使用滑轨式接触,使两平行不锈钢滑轨能够竖向滑动移动,平行不锈钢滑轨使用电机驱动滑动及遥控操作,滑轨小车也配有电机驱动装置和遥控装置,滑轨小车上承载准直器和探测器等测量装置,滑轨小车上开有准直器大小视野窗,在不锈钢架平台四个角加装四个电极驱动的轮子。
本发明专利涉及放射性核素源边界参数蒙特卡罗反演技术,包括测量模型和测量平台,所述的测量模型从上到下依次包括探测器、上层沙土衰减层、热区、热区虚拟点源、下层沙土衰减层、深层干扰源和干扰源虚拟点源,反演计算时需要将热区及深层干扰源虚拟成虚拟点源,虚拟点源位置在探测器对称轴线上,上层衰减层和热区厚度及深层干扰源厚度可调,其最终厚度由建立峰谷比及虚拟点探测效率方程和均方偏差反演计算确定。考虑到时效性,测量平台大小设计为5.0×6.0m,小矩形两条长边长度要长过大矩形长边,且搭在大矩形上,并固定,设计成滑轨结构,横向搭载两个平行不锈钢滑轨,滑轨与小矩形长边使用滑轨式接触,使两平行滑轨能够竖向滑动移动,滑轨车上承载准直器和探测器等测量装置,在平台四个角加装四个电极驱动的轮子,滑轨也假装电机驱动可实现远程遥控完成前后左右移动。
本发明使用虚拟技术刻度γ谱仪,验证其在就地γ谱仪搜索扫描测量细化污染源边界中的应用,主要将虚拟点源理论与蒙特卡罗参数表征技术两者相结合,建立源边界参数反演计算的理论模型及方法步骤,目前,理论计算和实验结果都验证所建立的计算模型和技术方法是正确的和实用的。目前该技术对于均匀分布的放射性核素已经能够精确确定放射性污染区域深度分布参数。
与现有技术相比,本发明的有益效果是:
该技术成果提高测量精度,大大缩小放射性可疑区的搜索范围,减少废物量,达到废物处置减容的目的,对于禁核试核查目标核弹头惰层厚度参数的确定具有重大的参考价值。
附图说明
图1为基于虚拟源原理的源边界参数蒙特卡罗反演技术的测量模型。
图2为基于虚拟源原理的源边界参数蒙特卡罗反演技术的测量滑轨车结构示意图。
图3为基于虚拟源原理的源边界参数蒙特卡罗反演技术中MCNP程序计算模型。
图中所示:1-探测器;2-上层沙土衰减层;3-热区;4-热区虚拟点源;5-下层沙土衰减层;6-深层干扰源;7-干扰源虚拟点源;8-滑轨小车;9-平行不锈钢滑轨;10-小矩形钢架;11-大矩形钢架;12-轮子。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1~3,本发明实施例中,一种基于虚拟源原理的源边界参数蒙特卡罗反演方法,包括不锈钢架平台、探测器1、上层沙土衰减层2、热区3、热区虚拟点源4、下层沙土衰减层5、深层干扰源6和干扰源虚拟点源7,所述热区3上方设置有沙土衰减层2,热区3下方设置有深层干扰源6,热区3和深层干扰源6之间设置有下层沙土衰减层5,上层沙土衰减层2和下层沙土衰减层5、热区3及深层干扰源6密度和成分一致,反演计算时需要将热区3及深层干扰源6虚拟成热区虚拟点源4和干扰源虚拟点源7,热区虚拟点源4和干扰源虚拟点源7位置在探测器1对称轴线上,上层沙土衰减层2和热区3的厚度和深层干扰源6厚度可调,其最终厚度由反演计算确定的虚拟点位置确定得到。
本发明包括以下步骤:
S1:等效虚拟点源
为简化模型,将热区3和深层干扰源6分别等效虚拟成其虚拟点源,再将两个虚拟点源等效成一个等效虚拟点源,该点源包含所有源的辐射信息。对于钚污染土,241Am和239Pu的单位体积活度比一般为常数,根据这一比值可定出它们的能峰241Am的59.54keV(分支比0.359,26.4keV分支比0.024)和239Pu的129keV(分支比0.000062)射线等效虚拟点源探测效率(也即等效虚拟点在对称轴上的位置);
S2:建立峰谷比及虚拟点探测效率方程
根据等效虚拟点射线能峰探测效率和实测能谱的峰谷比参数,再结合MCNP程序建立测量对象源模型,进行理论反演计算源的边界参数,其反演主要依据如下三个方程:
59.54keV峰谷比方程:
59.54keV虚拟点探测效率方程:
ε1(h1)p1+ε2(h2)(1-p1)=ε0(h0) (2)
26.4keV或129keV虚拟点探测效率方程:
ε1(h1)*p1+ε2(h2)*(1-p1)=ε0(h0)* (3)
v表示谷区计数,n表示峰计数,p表示权重因子,h表示虚拟点位置(虚拟点位置一般在对称轴线上变化),0表示实测谱值,1表示虚拟点位置1,2表示虚拟点位置2。
S3:计算均方偏差σ(X)
等式(3)右边是定值E(Xi)(实验值,又叫数学期望值),其中式(2)和式(3)等式右边的虚拟点源位置h0处的探测效率值是通过核素活度比值定出,峰谷比通过能谱获得,左边是蒙卡模拟结果,使用蒙卡模拟等式左边探测效率及峰谷比等(不同衰减层厚、不同源厚和不同权重,两两组合),再使用最小二乘法计算均方偏差σ(X),如式(4)所示,找出均方偏差最小的组合即可。
本发明包括从上到下依次为滑轨小车8、平行不锈钢滑轨9、小矩形钢架10、大矩形钢架11和轮子12,所述不锈钢架平台主要由钢管连接成网状结构,网状结构为一个大矩形钢架11套一个小矩形钢架10,对角线也用钢管连接,中间长度方向加一根钢管,考虑到时效性,不锈钢架平台大小设计为5.0×6.0m,小矩形钢架10两条长边长度大于大矩形钢架11长边,且搭在大矩形钢架11上固定,并进一步固定设计成滑轨结构,横向搭载两个平行不锈钢滑轨9,平行不锈钢滑轨9与小矩形钢架10长边使用滑轨式接触,使两平行不锈钢滑轨9能够竖向滑动移动,平行不锈钢滑轨9使用电机驱动滑动及遥控操作,滑轨小车8也配有电机驱动装置和遥控装置,满足整个车体在两个平行钢管上左右遥控滑动的功能,滑轨小车8上承载准直器和探测器等测量装置,滑轨小车8上开有准直器大小视野窗,为保证平台可以移动,在不锈钢架平台四个角加装四个电极驱动的轮子12,可实现远程遥控完成前后左右移动。
(一)实验研究
根据图1的理论计算模型,设计两种实验探测模式。探测模式1:75×25mm239Pu体源下方放置1.6cm厚沙土衰减层然后放置在距离平板型高纯锗探测器铍窗0.953cm位置处,源与衰减层密度都为1.34g/cm3;探测模式2:241Am点源通过支架放置在距离2.0cm处的0.5cm厚沙土衰减层上再将其放置在75×25mm厚241Am体源上,然后将它们整个放置在平板型高纯锗探测器铍窗进行放射性测量,源与沙土地密度都为1.43。谱获取时间足够长,保证使其统计涨落控制在1.0%以内,实验数据见表1所示。
表1实验谱能峰计数及处理结果
Table 1.The energy peak count of experimental spectrum and processresults
将计数通过计算转化为活度,活度计算结果见1中数据所示,其中探测效率数据是由蒙特卡罗模拟计算给出。对于探测模式1的情况,239Pu/241Am比值4.1(对于未知源项目标该比值一般是通过取样测量得到,或根据已有资料查询得到)。在探测模式2下,通过241Am两个能量分别计算得到各自总活度,其为点源活度加体源活度,且241Am(26.4keV)/241Am(59.54keV)比值理论上应该是1.0,但由于存在系统误差241Am(26.4keV)/241Am(59.54keV)比值计算结果为0.9,对于这个系统误差可以不用校正,因为后面计算都采用比值进行。241Am源谷区选择54-57keV能区,探测模式1下峰谷比为5.9,探测模式2下,峰谷比为6.9。
(二)理论计算研究
根据实验设计及探测器参数,建立如图3所示的蒙特卡罗计算模型,其中晶体的尺寸及位置坐标都已标出。探测器参数:晶体角死层被两个球及一个圆柱所切得到,球参数输入为sz 1.90 1.9,sz 3.4445 2.06。图中单位为cm。铍窗密度1.85,沙土成份在卡片中的输入为8016 -31.37 14028 -18.4 13000 -3.3 26056 -1.33 20040 -4.83 12024 -0.84。
对于探测器模式1,只有图1中沙土衰减层和体源,没有干扰源的存在,但模拟的过程中都是假设有干扰源存在的,其模拟计算过程如下:
1.等效虚拟点源
在对称轴上,计算不同位置处点源探测效率,其计算结果如表2中数据所示。
表2等效虚拟点源探测效率
Table 2.The detection efficiency of equivalent virtual point source
根据表2中探测效率值,结合表1中探测模式1的能峰数据计算不同位置坐标处的239Pu/241Am活度比值,可以看出,在坐标为-2.80处时239Pu/241Am比值为4.1,和表1中活度比值一致,因此,可以认为等效虚拟点源(如图2中)的位置坐标为(00 -2.80)。等效虚拟点源包含两部分的等效,一部分是测量对象体源的贡献,一部分是干扰源的贡献(干扰源可以是点源、面源和体源,我们统一将其转化成虚拟点源),将这两部分各自虚拟成点源,因此可以认为等效虚拟点源由体源虚拟点源和干扰源虚拟点源按不同的权重组合而成,因此下一步的目标就是找出这个权重和体源的虚拟点源位置坐标。
2.体源虚拟点源
等效虚拟点源的探测效率值7.0×10-3和9.45×10-3即为公式(2)和(3)右边的数学期望值,峰谷比5.9为公式(1)右边的数学期望值。根据表2中不同点位置处的探测效率,计算不同虚拟点位置和不同权重组合下的等效探测效率及等效峰谷比计算结果见表3所示,这些值作为公式(2)和(3)左边理论计算值,最后按理论计算公式(4)计算均方偏差,如表3所示。
表3不同组合下等效虚拟点的均方偏差
Table 3.The mean square deviation of equivalent virtual point atdifferent combination
省去中间过程,直接计算得到虚拟点坐标、权重及均方偏差的一一对应的表4。
表4不同组合下等效虚拟点的均方偏差
Table 4.The mean square deviation of equivalent virtual point atdifferent combination
表4续 不同组合下等效虚拟点的均方偏差
Table 4 follow.The mean square deviation of equivalent virtual pointat different combination
从表4中可以看出,均方偏差最小的值为0.167,其对应的体源虚拟点位置坐标为-2.40,干扰源对应的虚拟点位置坐标-3.20,而其对应的权重都为0.50,由于两个点位置很近权重又一样,因此,可以认为没有干扰源的存在,实际上等效虚拟点位置坐标就是体源的虚拟点位置坐标。
3.体源参数反演
如表2中坐标-2.80对应的探测效率及峰谷比作为公式(1)~(3)的右边的数学期望值,然后计算不同体源中心坐标下不同体源厚度下的探测效率及峰谷比值,将这些值作为理论计算值,并计算其与数学期望之间的均方偏差,计算结果如表5所示。
表5体源参数的反演计算结果
Table 5.The inversion result of volume source parameters
从表5中可以看出,均方偏差最小的为0.159,其对应的体源中心坐标为-2.45,体源厚度为2.5cm,和实际情况完全一致,充分说明该技术方法的可行性和准确性。对于探测模式2,探测器参数和探测模式1一致,但是沙土密度变为1.43,整个测量是贴近探测器进行的,129keV射线能量变成26.4keV。其计算过程如表6-8所示。
表6等效虚拟点探测效率
Table 6.The detection efficiency of equivalent virtual point source
表7不同组合下等效虚拟点的均方偏差
Table 7.The mean square deviation of equivalent virtual point atdifferent combination
表7续 不同组合下等效虚拟点的均方偏差
Table 7 follow.The mean square deviation of equivalent virtual pointat different combination
从表7中可以看出,均方偏差最小是0.112,其对应的上下虚拟点源权重都为0.5,上下虚拟点源坐标分别为0.0和-0.8cm,等效虚拟点源坐标为-0.20cm,比较靠近上虚拟点位置,因此可以认为,干扰源是存在的,且距上虚拟点位置较远,权重也占一半,因此,可以认为有一个强的干扰源的存在。干扰源只需证明它存在或者不存在即可,它的存在及其设置只是为更为准确地解出体源边界参数,在本发明中它的边界参数是不用求解出来的。
表8体源参数的反演计算结果
Table 8.The inversion result of volume source parameters
表8中,均方偏差的计算只用26.4keV射线能量探测效率和59.54keV射线峰谷比这两项,没有使用59.54keV射线能量探测效率。从表8中可以看出,最小均方偏差为0.0220,其对应的体源中心坐标为0.25cm,体源厚度为1.7cm,计算得沙土衰减层厚度为1.53mm,由于这一值相对体源厚度太小,在表层不可能存在这么薄的沙土衰减层,因此可以将其直接忽略,并入体源厚度,可认为体源厚度为2.0cm,相对真实值2.5cm,大概有20.0%左右的相对偏差,其偏差稍微偏大的原因是干扰点源使用低能26.4keV射线穿透3.0cm左右沙土衰减层后在探测器晶体内产生的有效计数太少,统计涨落过大,造成计算误差也大。因此,不能选用能量过低的低能射线,比如本发明选用的26.4keV能量射线就有些偏低,会造成结果一定的偏差。
本发明提供一种确定放射性可疑区范围的虚拟点源方法,其解决现有探测方法准确判定放射性污染区边界的技术问题。本发明可测量多种核素或测量发射多条特征能量射线的核素的源项目标,其过程是:首先根据核素活度浓度比值(一种核素该比值可认为是1.0,但要测量其两条及以上特征能量射线)找出等效虚拟点源,该虚拟点源特点是由该点算出的活度浓度比值和给出的一致,它包含目标体源及其他干扰源的信息;其次再将该等效虚拟点按一定的权重分成上下两个虚拟点源,若上下两个虚拟点源与等效虚拟点源距离过近且权重相等,则不需要将等效虚拟点源分成上下两个虚拟点源,它本身就是体源的虚拟点源,不存在干扰源,若上下两个虚拟点源距离等效点源距离较远或权重不一致,则可认为上虚拟点源就是目标体源的虚拟点源,该虚拟点源的探测效率及峰谷比和目标体源的一致;最后就是根据体源的虚拟点源的探测效率及峰谷比反演计算体源边界参数。该技术的优点是:1、利用虚拟点探测器理论及虚拟点源技术,成功地解决方程组法遇到深层干扰源情况下存在病态或不确定大的难题,使就地及航测γ谱仪发展为一项更为有效的探测手段。2、使用新原理及测量方式,不增加新的物理探测设备,大大减少经费开支,减少工作量及时间。
本发明包括从上到下依次为探测器1、上层沙土衰减层2、热区3、热区虚拟点源4、下层沙土衰减层5、深层干扰源6和干扰源虚拟点源7,源项又称热区3,上面一般会有沙土衰减层2,同时热区3下方有干扰源6的存在,热区3和干扰源6之间有衰减层5,一般热区3上面的沙土的衰减层2和干扰源6在一些源项简单测量中可能就不存在。探测器1晶体参数尤其是死层参数是经过点源实验表征过的,沙土衰减层2和5、热区3及深层干扰源6密度和成分一致,反演计算时需要将热区3及深层干扰源6虚拟成虚拟点源4和7,虚拟点源位置在探测器对称轴线上。上层衰减层2和热区3厚度及深层干扰源6厚度可调,其最终厚度由反演计算确定的虚拟点位置确定得到。
本发明包括以下步骤:
1】等效虚拟点源
为简化模型,将热区3和深层干扰源6分别等效虚拟成其虚拟点源,再将两个虚拟点源等效成一个等效虚拟点源,该点源包含所有源的辐射信息。对于钚污染土,241Am和239Pu的单位体积活度比一般为常数,根据这一比值可定出它们的能峰241Am的59.54keV(分支比0.359,26.4keV分支比0.024)和239Pu的129keV(分支比0.000062)射线等效虚拟点源探测效率(也即等效虚拟点在对称轴上的位置);
2】建立峰谷比及虚拟点探测效率方程
根据等效虚拟点射线能峰探测效率和实测能谱的峰谷比参数,再结合MCNP程序建立测量对象源模型,进行理论反演计算源的边界参数。其反演主要依据如下三个方程。
59.54keV峰谷比方程:
59.54keV虚拟点探测效率方程:
ε1(h1)p1+ε2(h2)(1-p1)=ε0(h0) (2)
26.4keV或129keV虚拟点探测效率方程:
ε1(h1)*p1+ε2(h2)*(1-p1)=ε0(h0)* (3)
v表示谷区计数,n表示峰计数,p表示权重因子,h表示虚拟点位置(虚拟点位置一般在对称轴线上变化),0表示实测谱值,1表示虚拟点位置1,2表示虚拟点位置2。
3】计算均方偏差σ(X)
等式(3)右边是定值E(Xi)(实验值,又叫数学期望值),其中式(2)和式(3)等式右边的虚拟点源位置h0处的探测效率值是通过核素活度比值定出,峰谷比通过能谱获得,左边是蒙卡模拟结果,使用蒙卡模拟等式左边探测效率及峰谷比等(不同衰减层厚、不同源厚和不同权重,两两组合),再使用最小二乘法计算均方偏差σ(X),如式(4)所示,找出均方偏差最小的组合即可。
本发明所需的装置包括从上到下依次为滑轨小车8、平行不锈钢滑轨9、小矩形钢架10、大矩形钢架11和轮子12。所述不锈钢架平台主要由钢管连接成网状结构,网的形状为一个大矩形11套一个小矩形10,对角线也用钢管连接,中间长度方向加一根钢管,考虑时效性,平台大小设计为5.0×6.0m,小矩形10两条长边长度要长过大矩形长边,且搭在大矩形11上,并固定,设计成滑轨结构,横向搭载两个平行不锈钢滑轨9,滑轨9与小矩形10长边使用滑轨式接触,使两平行滑轨9能够竖向滑动移动,平行滑轨9使用电机驱动滑动及遥控操作。滑轨小车8也配有电机驱动装置和遥控装置,满足整个车体在两个平行钢管上左右遥控滑动的功能,滑轨车8上承载准直器和探测器等测量装置,滑轨车上开有准直器大小视野窗。为保证平台可以移动,在平台四个角加装四个电极驱动的轮子12,可实现远程遥控完成前后左右移动。
本发明使用虚拟技术刻度γ谱仪,讨论其在就地γ谱仪搜索扫描测量细化污染源边界中的应用。主要将虚拟点源理论与蒙特卡罗参数表征技术两者相结合,建立源边界参数反演计算的理论模型及方法步骤,目前,理论计算和实验结果都验证所建立的计算模型和技术方法是正确的和实用的。目前该技术对于均匀分布的放射性核素已经能够精确确定放射性污染区域深度分布参数。
本发明的技术方案如下:HPGe探测器晶体一般为圆柱状,射线与晶体的作用可以认为是射线与晶体内一虚拟点探测器作用,即定义虚拟点探测器。有虚拟点探测器的理论,后面又出现虚拟点源概念,虚拟点源:对于辐射探测器测量对象为体源的情况下,在探测器中心对称轴线上存在有一个唯一的代表点位置,该位置的放射性点源全能峰探测效率与所对应的放射性体源的探测效率相等,即该点源叫做它们的虚拟点源。
在无源效率刻度中蒙卡模拟表征晶体参数的研究工作非常多技术方法也非常成熟,其基本原理是通过标准点源或体源(源的形状大小、材料和活度参数等完全已知)固定位置实验,获取探测器的实验刻度效率值,蒙特卡罗方法建立理论计算模型,调整理论计算模型中的晶体参数使其计算效率值和实验值一致,从而达到表征晶体参数的目的。在这一表征的过程中测量对象源的参数是已知的,测量仪器晶体参数是未知的,但在实际测量中会经常遇到大量测量对象未知的情况出现,针对这种情况,提出源边界参数表征技术,其基本原理和晶体参数表征技术类似。本发明采用逆向思维,晶体参数已知,源边界参数未知,通过蒙特卡罗方法建立计算模型,调整源边界参数使计算得到峰谷比及效率值与实验值一致。
一般测量时,源项又称热区上面一般会有沙土衰减层,同时热区下方有干扰源的存在,热区和干扰源之间有衰减层,一般热区上面的沙土的衰减层和干扰源在一些源项简单测量中可能就不存在。本发明的目标就是在有干扰源存在的情况下反演出热区位置及厚度参数。为简化模型,将热区和深层干扰源分别等效虚拟成其虚拟点源,再将两个虚拟点源等效成一个等效虚拟点源,该点源包含所有源的辐射信息。对于钚污染土,241Am和239Pu的单位体积活度比一般为常数,根据这一比值可定出它们的能峰241Am的59.54keV(分支比0.359,26.4keV分支比0.024)和239Pu的51.62keV(分支比0.000271,129keV分支比0.000062)射线等效虚拟点源探测效率(也即等效虚拟点在对称轴上的位置,该位置有且只有一个,通过该位置两个射线能量探测效率算得的239Pu/241Am活度比与已知条件一致);根据等效虚拟点射线能峰探测效率和实测能谱的峰谷比参数,再结合MCNP程序建立测量对象源模型,进行理论反演计算源的边界参数。
本发明的有益效果:该技术成果提高测量精度,大大缩小放射性可疑区的搜索范围,减少废物量,达到废物处置减容的目的,对于禁核试核查目标核弹头惰层厚度参数的确定具有重大的参考价值。
本发明为扩大测量面积和范围,设计矩形扫描测量支架,矩形对角线也用钢管连接,中间长度方向加一根钢管,考虑时效性,平台大小设计为5.0×6.0m,小矩形两条长边长度要长过大矩形长边,且搭在大矩形上,并固定,设计成滑轨结构,横向搭载两个平行不锈钢滑轨,滑轨与小矩形长边使用滑轨式接触,使两平行滑轨能够竖向滑动移动,平行滑轨使用电机驱动滑动及遥控操作。滑轨小车也配有电机驱动装置和遥控装置,满足整个车体在两个平行钢管上左右遥控滑动的功能,滑轨车上承载准直器和探测器等测量装置,滑轨车上开有准直器大小视野窗。为保证平台可以移动,在平台四个角加装四个电极驱动的轮子,可实现远程遥控完成前后左右移动。
尽管参照前述实施例对本发明进行详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于虚拟源原理的源边界参数蒙特卡罗反演方法,包括探测器(1)、上层沙土衰减层(2)、热区(3)、热区虚拟点源(4)、下层沙土衰减层(5)、深层干扰源(6)和干扰源虚拟点源(7),其特征在于:还包括不锈钢架平台,不锈钢架平台主要由钢管连接成网状结构,网状结构为一个大矩形钢架(11)套一个小矩形钢架(10),对角线也用钢管连接,中间长度方向加一根钢管,不锈钢架平台大小设计为5.0×6.0m,小矩形钢架(10)两条长边长度大于大矩形钢架(11)长边,且搭在大矩形钢架(11)上,并进一步固定设计成滑轨结构,横向搭载两个平行不锈钢滑轨(9),平行不锈钢滑轨(9)与小矩形钢架(10)长边使用滑轨式接触,使两平行不锈钢滑轨(9)能够竖向滑动移动,平行不锈钢滑轨(9)使用电机驱动滑动及遥控操作,滑轨小车(8)也配有电机驱动装置和遥控装置,滑轨小车(8)上承载准直器和探测器,滑轨小车(8)上开有准直器大小视野窗,在不锈钢架平台四个角加装四个电极驱动的轮子(12);所述热区(3)上方设置有沙土衰减层(2),热区(3)下方设置有深层干扰源(6),热区(3)和深层干扰源(6)之间设置有下层沙土衰减层(5),上层沙土衰减层(2)和下层沙土衰减层(5)、热区(3)及深层干扰源(6)密度和成分一致,反演计算时需要将热区(3)及深层干扰源(6)虚拟成热区虚拟点源(4)和干扰源虚拟点源(7),热区虚拟点源(4)和干扰源虚拟点源(7)位置在探测器(1)对称轴线上,上层沙土衰减层(2)和热区(3)的厚度和深层干扰源(6)厚度可调,其最终厚度由反演计算确定的虚拟点位置确定得到;反演计算包括以下步骤:
S1:等效虚拟点源
为简化模型,将热区(3)和深层干扰源(6)分别等效虚拟成其虚拟点源,再将两个虚拟点源等效成一个等效虚拟点源,该点源包含所有源的辐射信息;对于钚污染土,241Am和239Pu的单位体积活度比一般为常数,根据这一比值可定出它们的能峰241Am的59.54keV和239Pu的129keV射线等效虚拟点源探测效率,也即等效虚拟点在对称轴上的位置;注:241Am的能峰59.54keV分支比0.359,其26.4keV的能峰分支比0.024;239Pu的能峰129keV分支比0.000062;
S2:建立峰谷比及虚拟点探测效率方程
根据等效虚拟点射线能峰探测效率和实测能谱的峰谷比参数,再结合MCNP程序建立测量对象源模型,进行理论反演计算源的边界参数,其反演主要依据如下三个方程:
59.54keV峰谷比方程:
59.54keV虚拟点探测效率方程:
ε1(h1)p1+ε2(h2)(1-p1)=ε0(h0) (2)
26.4keV或129keV虚拟点探测效率方程:
ε1(h1)*p1+ε2(h2)*(1-p1)=ε0(h0)* (3)
v表示谷区计数,n表示峰计数,p表示权重因子,h表示虚拟点位置(虚拟点位置一般在对称轴线上变化),0表示实测谱值,1表示虚拟点位置1,2表示虚拟点位置2;
S3:计算均方偏差σ(X)
等式(3)右边是定值E(Xi),定值E(Xi)为实验值,又叫数学期望值,其中式(2)和式(3)等式右边的虚拟点源位置h0处的探测效率值是通过核素活度比值定出,峰谷比通过能谱获得,左边是蒙卡模拟结果,使用蒙卡模拟等式左边探测效率及峰谷比,即将不同衰减层厚、不同源厚和不同权重两两组合,再使用最小二乘法计算均方偏差σ(X),如式(4)所示,找出均方偏差最小的组合即可,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010685178.XA CN112051601B (zh) | 2020-07-16 | 2020-07-16 | 一种基于虚拟源原理的源边界参数蒙特卡罗反演技术 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010685178.XA CN112051601B (zh) | 2020-07-16 | 2020-07-16 | 一种基于虚拟源原理的源边界参数蒙特卡罗反演技术 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112051601A CN112051601A (zh) | 2020-12-08 |
CN112051601B true CN112051601B (zh) | 2023-03-14 |
Family
ID=73602316
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010685178.XA Expired - Fee Related CN112051601B (zh) | 2020-07-16 | 2020-07-16 | 一种基于虚拟源原理的源边界参数蒙特卡罗反演技术 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112051601B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113703034B (zh) * | 2021-07-14 | 2024-08-30 | 中国人民解放军63653部队 | 反演污染区虚拟点位置和边界的技术方法及设备 |
CN113885070A (zh) * | 2021-09-13 | 2022-01-04 | 中国人民解放军63653部队 | 确定均匀分布和指数分布体源虚拟点位置的方法和设备 |
CN115856994B (zh) * | 2023-02-02 | 2023-05-05 | 北京中科核安科技有限公司 | 伽玛探测器的效率校准方法、装置和系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010011844A1 (en) * | 2008-07-25 | 2010-01-28 | Tufts Medical Center | A system and method of clinical treatment planning of complex, monte carlo-based brachytherapy dose distributions |
CN104020483A (zh) * | 2014-06-13 | 2014-09-03 | 西北核技术研究所 | 采用虚拟源模拟刻度探测器探测效率的方法及装置 |
CN104483698A (zh) * | 2014-12-30 | 2015-04-01 | 西北核技术研究所 | 基于点源效率函数的氙气源虚拟点源刻度装置及方法 |
CN106291650A (zh) * | 2016-08-31 | 2017-01-04 | 广州市岱尼欣贸易有限公司 | 基于蒙特卡罗的剂量测量方法 |
DE102016225502A1 (de) * | 2016-12-19 | 2018-06-21 | Moba Mobile Automation Ag | Messsystem zur Schichtdickenerfassung |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999008066A1 (en) * | 1997-08-05 | 1999-02-18 | Medar, Inc. | Interference method and system for measuring the thickness of an optically-transmissive thin layer formed on a relatively planar, optically-reflective surface |
US9116100B2 (en) * | 2009-05-26 | 2015-08-25 | Kromek Limited | Method for the identification of materials in a container |
CN103853929B (zh) * | 2014-03-17 | 2016-06-15 | 东华理工大学 | 一种基于蒙卡响应矩阵的低分辨率γ能谱反演解析系统及方法 |
CN106896396A (zh) * | 2015-12-18 | 2017-06-27 | 中国辐射防护研究院 | 一种用于校准惰性气体流出物离线监测仪探测效率的方法 |
CN108287357B (zh) * | 2018-01-15 | 2020-12-25 | 东华理工大学 | 一种圆柱体探测器对点源的源峰探测效率获取方法 |
-
2020
- 2020-07-16 CN CN202010685178.XA patent/CN112051601B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010011844A1 (en) * | 2008-07-25 | 2010-01-28 | Tufts Medical Center | A system and method of clinical treatment planning of complex, monte carlo-based brachytherapy dose distributions |
CN104020483A (zh) * | 2014-06-13 | 2014-09-03 | 西北核技术研究所 | 采用虚拟源模拟刻度探测器探测效率的方法及装置 |
CN104483698A (zh) * | 2014-12-30 | 2015-04-01 | 西北核技术研究所 | 基于点源效率函数的氙气源虚拟点源刻度装置及方法 |
CN106291650A (zh) * | 2016-08-31 | 2017-01-04 | 广州市岱尼欣贸易有限公司 | 基于蒙特卡罗的剂量测量方法 |
DE102016225502A1 (de) * | 2016-12-19 | 2018-06-21 | Moba Mobile Automation Ag | Messsystem zur Schichtdickenerfassung |
Non-Patent Citations (2)
Title |
---|
大面积及大体积放射源HPGe探测器效率刻度;田自宁等;《现代应用物理》;20150326(第01期);全文 * |
蒙特卡罗方法计算高纯锗探测器的全能峰效率;卫晓峰等;《辐射防护》;20130520(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112051601A (zh) | 2020-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112051601B (zh) | 一种基于虚拟源原理的源边界参数蒙特卡罗反演技术 | |
Klem et al. | Measurement of asymmetries of inclusive pion production in proton-proton interactions at 6 and 11.8 GeV/c | |
CN104536029B (zh) | 基于多个NaI(T1)探测器的放射性物质二维定位方法 | |
CN103605149B (zh) | 一种刻度氙气样品HPGe探测效率的装置及方法 | |
CN104483693B (zh) | 一种非均匀分布源探测效率计算及模拟装置及方法 | |
CN104483698B (zh) | 基于点源效率函数的氙气源虚拟点源刻度装置及方法 | |
CN104020507B (zh) | 基于虚拟点探测器原理的放射性区域确定方法 | |
Mitrica et al. | A mobile detector for measurements of the atmospheric muon flux in underground sites | |
Yáñez et al. | A method to measure the integral vertical intensity and angular distribution of atmospheric muons with a stationary plastic scintillator bar detector | |
Serebrov et al. | Experiment neutrino-4 on searching for a sterile neutrino with multisection detector model | |
Prieto-Pena et al. | Impact of charge collection efficiency and electronic noise on the performance of solid-state 3D microdetectors | |
CN107884805B (zh) | 一种射源定位的方法及装置 | |
CN108287357A (zh) | 一种圆柱体探测器对点源的源峰探测效率获取方法 | |
Clark et al. | A Monte Carlo model for in situ prompt gamma-ray analysis probes | |
Oláh et al. | Improvement of cosmic-ray muography for Earth sciences and civil engineering | |
Aron et al. | Recovery of the Surface Density of Radionuclide Contamination of Soil from Aerogamma Survey Data | |
Myers et al. | Compton scattering from Li 6 at 86 MeV | |
Alwars et al. | Development of a high-precision gamma-ray source finder based on a single detector | |
Chase et al. | π− p Charge Exchange at Backward Angles | |
Guckes et al. | Development and deployment of the collimated directional radiation detection system | |
Shirakawa et al. | Remote sensing of nuclear accidents using a direction finding detector | |
Meckbach et al. | Unfolding of Ge-detector response by Monte Carlo simulation and applications in radiation protection | |
Mitrica et al. | A mobile detector for measurements of the atmospheric muon flux | |
Duarte | Muon tomography: application of image reconstruction algorithms within the LouMu project | |
Fedorov et al. | STUDY OF 14.1 MeV NEUTRONS INELASTIC SCATTERING ON IRON |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20230314 |
|
CF01 | Termination of patent right due to non-payment of annual fee |