CN104020507A - 基于虚拟点探测器原理的放射性区域确定方法 - Google Patents

基于虚拟点探测器原理的放射性区域确定方法 Download PDF

Info

Publication number
CN104020507A
CN104020507A CN201410264383.3A CN201410264383A CN104020507A CN 104020507 A CN104020507 A CN 104020507A CN 201410264383 A CN201410264383 A CN 201410264383A CN 104020507 A CN104020507 A CN 104020507A
Authority
CN
China
Prior art keywords
detector
virtual point
active area
point source
formula
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.)
Granted
Application number
CN201410264383.3A
Other languages
English (en)
Other versions
CN104020507B (zh
Inventor
田自宁
欧阳晓平
王军
张显鹏
张忠兵
刘林月
阮金陆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwest Institute of Nuclear Technology
Original Assignee
Northwest Institute of Nuclear Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwest Institute of Nuclear Technology filed Critical Northwest Institute of Nuclear Technology
Priority to CN201410264383.3A priority Critical patent/CN104020507B/zh
Publication of CN104020507A publication Critical patent/CN104020507A/zh
Application granted granted Critical
Publication of CN104020507B publication Critical patent/CN104020507B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Radiation (AREA)

Abstract

本发明提供一种基于虚拟点探测器的放射性区域确定方法,包括拟合得到虚拟点探测器距离探测器下端面距离的拟合参数、求出地面放射性区域对应的虚拟点源位置、建立均匀分布地面放射性区域半径r的和其对应的虚拟点源位置x的函数关系式、地面放射性区域半径参数的反演计算等步骤,本发明利用虚拟点探测器理论及虚拟点源技术,成功地解决了以往高空无法确定放射性可疑区边界的问题,使空中发展为一项更为有效的探测手段。

Description

基于虚拟点探测器原理的放射性区域确定方法
技术领域
本发明涉及核技术应用领域中发现和确定可疑放射源位置及大小的探测技术。具体涉及测量对象为均匀分布放射性区域、点源、点探测器等装置。
背景技术
目前,在高空针对可疑放射性区域的测量,只能给出污染源的大概位置,而对于源的污染边界的确定是个很难的问题。对于放射性热点的确定主要是用GPS的定位加上航迹视频采图的判断,使得热点的纵向定位误差缩小到15m之内,横向误差达到132m,该技术根据视频采图判断存在很大的不确定性,误差较大的缺点。有研究者将卷积定理和反卷积定理应用到任意面状放射源产生的辐射场求解问题,该技术利用空中测量的某一高度平面上的照射量分布值,应用反卷积定理可求得地面上对应的放射性活度分布。该方法成立的前提条件必须做如下假设:不考虑空中悬浮存在的放射性物质贡献,忽略空气对γ射线的吸收和反散射作用,忽略α、β、n等其它射线的贡献,假设平均一次衰变放出一个γ射线。在做这些假设的时,由于与实际情况有差别,因此会产生很大误差,尤其忽略空气对γ射线的吸收和反散射作用,因为100keV射线经过100m空气几乎全部衰减。
发明内容
本发明目的是提供一种确定放射性可疑区范围的虚拟点探测器及点源的方法,其解决了现有探测方法准确判定放射性污染区边界的技术问题,该技术成果大大缩小了放射性可疑区的搜索范围,减少了工作量,提高了测量精度。
本发明的技术解决方案是:
基于虚拟点探测器的放射性区域确定方法,其特殊之处在于:包括以下步骤:
1)虚拟点探测器2在探测器1内或者在探测器外,按照下式进行最小二乘拟合得到虚拟点探测器2距离探测器1下端面距离的拟合参数h0
ϵ ( x 0 ) ϵ ( x ) - 1 = 1 ( x 0 + h 0 ) ( x - x 0 )
式中:
x0表示探测器1轴线上参考点位置,
x表示点源在探测器1轴线上任意点位置,h0表示虚拟点探测器的位置;
ε(x0)、ε(x)为已知的参考值,分别表示点源在参考点位置的和任意点位置时的探测效率;ε(x)是以x为自变量的应变量;
2)设地面放射性区域6探测效率为ε(x),则按照下式求出地面放射性区域6对应的虚拟点源8位置x1=x:
ϵ ( x 0 ) ϵ ( x ) e - μx - 1 = x - x 0 h 0 + x 0
3)建立均匀分布地面放射性区域6半径r的和其对应的虚拟点源8位置x的函数关系式r(x):
3.1)将标准点源3置于探测器轴线上不同位置,获取其探测效率,得到标准点源3探测效率与其高度的拟合函数:
εp(x)=a·eb·x
式中:
εp(x)表示点源在不同高度时的探测效率;
a,b表示拟合参数;
3.2)根据放射性区域6距探测器中的晶体一定距离的探测效率εd(r)等于点源效率可推出其虚拟点源8位置:
ϵ d ( r ) = ϵ p ( x ) = a · e b · x ⇒ x = ln ( ϵ d ( r ) a ) / b
3.3)对地面放射性区域6不同半径r对应的x,r(x)的函数关系式进行最小二乘法拟合,得到公式中的参数c,d,e的值:
r(x)=c·x2+d·x+e
式中:
c,d,e表示拟合参数;
4)地面放射性区域6半径参数的反演计算:
4.1)将步骤2中计算得到的虚拟点源8位置x=x1带入步骤3.3的公式中,计算得到地面放射性区域6半径r;
4.2)将求出的地面放射性区域6半径和其实际值进行比较,以验证该技术方法的正确性和有效性。
本发明专利的优点是:
1、本发明利用虚拟点探测器理论及虚拟点源技术,成功地解决了以往高空无法确定放射性可疑区边界的问题,使空中发展为一项更为有效的探测手段。
2、本发明由于使用了新的数据处理原理及测量方式,而不是增加新的物理探测设备,因此大大减少了经费开支,减少了工作量和工作时间。
3、本发明由于是对射线探测原理及方法的创新,因此,对于射线探测的其他相关领域都具有通用性,其应用范围非常广泛。
4、本专利将虚拟点源理论再结合国外的虚拟点探测器理论,可直接对测量能谱进行分析,由于是对射线探测原理的基础理论的发展创新不需要作任何假设就可以确定放射性污染区域的面积大小,因此具有非常高的准确性。
5、目前该技术对于均匀分布的放射性核素已经能够精确确定污染放射性区域半径大小。在核应急领域若使用航测技术准确探测污染源或丢失源的位置范围将会大大减少探测时间和及时找回丢失源,从而将危害降低到最低点。
6、本专利主要针对半径为100m的均匀分布137Cs(661.66keV射线)放射性区域进行了研究,测量高度为100m。
附图说明
图1是虚拟点探测器测量示意图;
图2是虚拟点探测器位置求解曲线示意图;图3是在对称中心点源高度和其效率的拟合曲线示意图;
图4是点源效率拟合曲线示意图;
图5是放射性区域半径和其虚拟点位置的拟合曲线示意图;
附图标记:1-探测器;2-虚拟点探测器;3-点源;4-对称轴;5-参考点源;6-地面放射性区域;7-探测区域;8-虚拟点源。
具体实施方式
探测器1位于探测区域7的正上方,地面放射性区域6在探测区域内,探测区域7的半径为R,探测器1相对探测区域7足够远时可认为它是一个点;
虚拟点探测器2在探测器1内或者在探测器外,虚拟点探测器2距离探测器1下端面距离为h0
参考点源5距离探测器1下端面距离为X。,点源3距离探测器1下端面距离为x,虚拟点源8距离探测器1下端面距离为x1
地面放射性区域6的半径为r,放射性核素均匀分布在放射性区域6内,探测器1探测区域7之间的吸收介质为空气。
对称轴4为探测器1的中心对称轴,虚拟点探测器2、参考点源5、点源3、虚拟点源8都在对称轴上;
虚拟点源8:
对于辐射探测器测量的对象地面放射性区域6,在探测器1中心对称轴线上存在有一个唯一的代表点位置,该位置的点源探测效率与地面放射性区域6效率相等,即该点源叫地面放射性区域6的虚拟点源8。
虚拟点探测器2:
对于辐射探测器1,在探测器1中心对称轴线上存在有一个唯一的代表点位置,其点探测器探测效率与探测器1的相等且满足平方反比定律,即该点探测器叫探测器1的虚拟点探测器2。
基于虚拟点探测器的放射性区域确定方法,包括以下步骤:
1】求解虚拟点探测器2的位置h0。虚拟点探测器2的计算模型如图1所示。根据虚拟点探测器理论参考点源5与点源3探测效率比或者计数率比与距离平方成反比定律:
ϵ ( x 0 ) ϵ ( x ) = C ( x 0 ) C ( x ) = ( x + h 0 ) 2 ( x 0 + h 0 ) 2 - - - ( 1 )
式(1)也可以改写成如下式(2)和式(3)。
C ( x 0 ) C ( x ) - 1 = 1 ( x 0 + h 0 ) ( x - x 0 ) - - - ( 2 )
ϵ ( x 0 ) ϵ ( x ) - 1 = 1 ( x 0 + h 0 ) ( x - x 0 ) - - - ( 3 )
式中:
x0表示探测器1轴线上参考点位置,x表示点源在探测器1轴线上任意点位置,h0表示虚拟点探测器的位置;
C(x0)、C(x)分别表示点源在参考点位置的和任意点位置时的峰计数;
ε(x0)、ε(x)分别表示点源在参考点位置的和任意点位置时的探测效率。
在求解公式(2)和(3)中:ε(x0)和x0为已知的参考值,ε(x)是以x为自变量的应变量,h0是拟合参数,因此对公式(3)进行最小二乘拟合可得拟合参数h0。由于实验测量时,一般得到的是测量源的峰计数,因此在实际应用中可以用峰计数率之比可以替代探测效率比,即主要使用公式(2)来进行计算,本专利主要使用公式(3)进行计算。
2】求出地面放射性区域6对应的虚拟点源8位置x1
设地面放射性区域6探测效率为ε(x),其中x1为放射性区域对应的虚拟点源位置,为未知参数。将地面放射性区域6探测效率等效为虚拟点源效率,其和参考位置点源5效率与它们的距离成平方反比的关系,同时考虑可得点源3与探测器1之间有吸收介质时,公式(3)平方反比定律变为下式。
ϵ ( x 0 ) ϵ ( x ) e - μx - 1 = x - x 0 h 0 + x 0 - - - ( 4 )
在公式(4)中,除了x外其他都是已知参数,因此可根据公式(5)求出虚拟点源8位置x=x1
3】建立均匀分布地面放射性区域6半径r的和其对应的虚拟点源8位置x的函数关系式r(x)。
将标准点源3置于探测器轴线上不同位置,获取其探测效率,得到点源3探测效率与其高度的拟合函数
εp(x)=a·eb·x (4)
式中:
εp(x)表示点源在不同高度时的探测效率;
a,b表示拟合参数。
可根据放射性区域6距晶体一定距离的探测效率εd(r)等于点源效率可推出其虚拟点源8位置,用公式可表示为:
ϵ d ( r ) = ϵ p ( x ) = a · e b · x ⇒ x = ln ( ϵ d ( r ) a ) / b - - - ( 5 )
得出地面放射性区域6不同半径r对应的x,r(x)的函数关系式如式(6)所示,可用最小二乘法拟合公式(6),得到公式中的参数c,d,e的值。
r(x)=c·x2+d·x+e (6)
式中:
c,d,e表示拟合参数。
4】地面放射性区域6半径参数的反演计算。
将步骤2中计算得到的虚拟点源8位置x=x1带入公式(6),计算可以得到地面放射性区域6半径r,求出的地面放射性区域6半径和实验值进行比较。
1.虚拟点探测器的位置h0
虚拟点探测器刻度方法首先用于点源探测效率的刻度,主要用参考点标准点源的探测效率值来计算不同位置处点源的探测效率。其计算模型如下图1所示。
HPGe探测器晶体一般为圆柱状体源,射线与晶体的作用可以认为是射线与晶体内一虚拟点探测器作用,即将晶体等效为一个虚拟的点,该虚拟点距离晶体下表面距离为h0,参考点源位置距离晶体下表面距离为x0,所求点源位置距离晶体上表面距离为x,如图1所示。则根据点源探测效率比或者计数率比与距离平方成反比定律有:
ϵ ( x 0 ) ϵ ( x ) = C ( x 0 ) C ( x ) = ( x + h 0 ) 2 ( x 0 + h 0 ) 2 - - - ( 1 )
C ( x 0 ) C ( x ) - 1 = 1 ( x 0 + h 0 ) ( x - x 0 ) - - - ( 2 )
ϵ ( x 0 ) ϵ ( x ) - 1 = 1 ( x 0 + h 0 ) ( x - x 0 ) - - - ( 3 )
从公式(1)中我们可以看出参考点源与所求点源的探测效率之比就等于它们之间的峰计数率之比,由于本文使用蒙卡模拟计算,因此在下文中探测效率比可以替代实验中的峰计数率之比,主要使用公式(3)来进行计算即可,而实际测量的过程中则可以使用公式(2)进行计算。当点源与探测器之间距离较远时,平方反比定律变为。
C ( x ) ( x + h 0 ) 2 = C 0 ( h 0 + x 0 ) 2 e - μx → C 0 C ( x ) e - μx - 1 = x - x 0 h 0 + x 0 - - - ( 4 )
C 0 C ( x ) e - μh - 1 = ϵ 0 ϵ ( x ) e - μx - 1 = x - x 0 h 0 + x 0 - - - ( 5 )
公式(5)中的峰计数率比也可以用探测效率比来替代,它可以用来确定虚拟点源地位置x1
因此针对本专利使用的BE3830型HPGe探测器,我们计算了其虚拟点探测器位置。其计算数据如表1中数据所示,x0=0.55cm,ε(x0)=7.84E-02,涨落为0.0015。
表1 点源在对称轴上不同高度探测效率数据
根据公式(3),将表1中的x-x0列对应的数据进行最小二乘拟合,如图2所示的结果。点源到晶体的距离必须为点源到晶体表面的距离,而不是到晶体碳窗的距离。
如图2中,根据公式(3),将一次函数的系数求倒数再减去x0,即可求出1/0.379≈2.6,h0=(2.6-0.55)cm=2.05cm。
2.放射性区域半径和其对应的虚拟点源位置的函数关系式
通过蒙卡程序我们可以根据准直射线的衰减情况计算出137Cs,661.66keV射线的线衰减系数为0.00924m-1。为了确定大放射性区域的半径和其虚拟点源的位置的函数关系式,我们需要做两个方面的模拟工作,一是在100m高空,模拟计算不同半径放射性区域的探测效率;二是建立大于100m高度时,对称轴上点源探测效率和其高度的函数关系式,由于放射性区域探测效率相对同高度的点源效率要小,因此其虚拟点源位置在100m以上,因此选择100m以上。首先,如表2所示为点源在对称中心的模拟实验数据,图3为它们的拟合结果。
表2 点源在对称中心的模拟实验数据
其次,如表3中点源在放射性区域上不同半径处的探测效率数据。放射性区域探测效率的求解公式如公式(6),(7)所示。
表3 100高度放射性区域模拟数据及计算结果
f s 1 S ∫ ∫ f p ( r ) ds = 2 π S ∫ 0 R r · f p ( r ) dr - - - ( 6 )
fp(r)=exp(a+br+cr2)(7)
表3中r表示点源距离对称轴的距离,fp(r)表示点源在不同半径位置处的探测效率,fs表示放射性区域探测效率,这三个参数含义和公式(6),(7)一样。公式(7)的最小二乘拟合函数如图4所示,将图4中的拟合曲线带入公式(6)进行不同半径0~100m,0~80m,0~60m,0~40m,0~20m的积分可以得到半径为100m,80m,60m,40m,20m的放射性区域探测效率,如表3种fs所对应的列所示。再将放射性区域探测效率fs的值代入到图3中的拟合关系式中就可以求出对应的虚拟点源的位置,如表3中“虚拟点源位置”所对应的数据所示。将放射性区域半径和虚拟点源位置进行最小二乘拟合处理即可得到放射性区域半径和其对应的虚拟点源位置的函数关系式如图5所示。
3.放射性区域半径的反演计算
假设表3中放射性区域探测效率fs为实验测量值,其对应的放射性区域半径为未知参数,则放射性区域探测效率fs对应的虚拟点源和虚拟点探测器的参考位置应该满足探测效率或峰计数率与它们的距离平方成反比,即公式(5),其中峰计数率换成探测效率比。根据公式(5)可以求解出fs对应的虚拟点源的位置(其中公式中的参数有:137Cs,661.66keV射线的线衰减系数为0.00924m-1,参考位置x0=0.55cm,h0=2.05cm,ε(x0)=7.84E-02,ε(x)=fs),通过公式(5)的计算可以得到虚拟点位置,如表3中“平方反比定律求出的虚拟点高度”所对应的列数据,最后再结合图5中的拟合关系式,我们可以最终反演求解得到放射性区域半径,如表3中“反演对应的放射性区域半径”那一列。求解出的放射性区域半径和实际计算使用的放射性区域半径进行比较,如表3中的相对偏差那一列所示,可以看出当污染区面积较小时偏差较大,这是因为放射性区域越小定位它的边界越困难,也符合实际情况,因为放射性区域越小,边界效率贡献越小;而放射性区域面积越大,定位其边界越准确,因此偏差也越小,这也和实际情况较为符合,因为半径大了边界效率贡献也越大,也就越好确定了。该技术成果只需在高空探测时测得污染核素的峰计数率就可以定出污染源的大概位置信息,在核应急及核查领域具有非常重大的核安全及军事意义。下一步工作重点是解决非均匀分布情况下的源的边界问题,若这个问题突破了,将会产生更为革命性的创新。由于未知放射性区域的探测效率是不知道的,因此根据步骤2中的放射性区域半径和其探测效率建立的函数关系式直接推出未知放射性区域的半径是行不通的,必须借助于虚拟点源、虚拟点探测器、平方反比定律和峰计数率比等条件才能求出实际的放射性区域半径。

Claims (1)

1.基于虚拟点探测器的放射性区域确定方法,其特征在于:包括以下步骤:
1)虚拟点探测器(2)在探测器(1)内或者在探测器外,按照下式进行最小二乘拟合得到虚拟点探测器(2)距离探测器(1)下端面距离的拟合参数h0
ϵ ( X 0 ) ϵ ( x ) - 1 = 1 ( x 0 + h 0 ) ( x - x 0 )
式中:
x0表示探测器(1)轴线上参考点位置,
x表示点源在探测器(1)轴线上任意点位置,h0表示虚拟点探测器的位置;
ε(x0)、ε(x)为已知的参考值,分别表示点源在参考点位置的和任意点位置时的探测效率;ε(x)是以x为自变量的应变量;
2)设地面放射性区域(6)探测效率为ε(x),则按照下式求出地面放射性区域(6)对应的虚拟点源(8)位置x1=x:
ϵ ( x 0 ) ϵ ( x ) e - μx - 1 = x - x 0 h 0 + x 0
3)建立均匀分布地面放射性区域(6)半径r的和其对应的虚拟点源(8)位置x的函数关系式r(x):
3.1)将标准点源(3)置于探测器轴线上不同位置,获取其探测效率,得到标准点源(3)探测效率与其高度的拟合函数:
εp(x)=a·eb·x
式中:
εp(x)表示点源在不同高度时的探测效率;
a,b表示拟合参数;
3.2)根据放射性区域(6)距探测器中的晶体一定距离的探测效率εd(r)等于点源效率可推出其虚拟点源(8)位置:
ϵ d ( r ) = ϵ p ( x ) = a · e b · x ⇒ x = ln ( ϵ d ( r ) a ) / b
3.3)对地面放射性区域(6)不同半径r对应的x,r(x)的函数关系式进行最小二乘法拟合,得到公式中的参数c,d,e的值:
r(x)=c·x2+d·x+e
式中:
c,d,e表示拟合参数;
4)地面放射性区域(6)半径参数的反演计算:
4.1)将步骤2中计算得到的虚拟点源(8)位置x=x1带入步骤3.3的公式中,计算得到地面放射性区域(6)半径r;
4.2)将求出的地面放射性区域(6)半径和其实际值进行比较,以验证该技术方法的正确性和有效性。
CN201410264383.3A 2014-06-13 2014-06-13 基于虚拟点探测器原理的放射性区域确定方法 Active CN104020507B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410264383.3A CN104020507B (zh) 2014-06-13 2014-06-13 基于虚拟点探测器原理的放射性区域确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410264383.3A CN104020507B (zh) 2014-06-13 2014-06-13 基于虚拟点探测器原理的放射性区域确定方法

Publications (2)

Publication Number Publication Date
CN104020507A true CN104020507A (zh) 2014-09-03
CN104020507B CN104020507B (zh) 2017-06-09

Family

ID=51437362

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410264383.3A Active CN104020507B (zh) 2014-06-13 2014-06-13 基于虚拟点探测器原理的放射性区域确定方法

Country Status (1)

Country Link
CN (1) CN104020507B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107195342A (zh) * 2017-04-05 2017-09-22 中广核研究院有限公司 一种测量核电厂核素活度分布的方法
CN109814147A (zh) * 2019-01-01 2019-05-28 中国人民解放军63653部队 使用探测器加衰减体测量小孔周围γ放射性核素分布的方法
CN110888155A (zh) * 2019-10-11 2020-03-17 中国人民解放军陆军防化学院 一种大体积液态放射源对γ谱仪有效贡献范围确定方法
CN113703034A (zh) * 2021-07-14 2021-11-26 中国人民解放军63653部队 反演污染区虚拟点位置和边界的技术方法及设备
CN113885070A (zh) * 2021-09-13 2022-01-04 中国人民解放军63653部队 确定均匀分布和指数分布体源虚拟点位置的方法和设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004245649A (ja) * 2003-02-12 2004-09-02 Central Res Inst Of Electric Power Ind バックグラウンド計数率の推定方法及びバックグラウンド計数率を補正する測定装置
CN201732166U (zh) * 2010-06-12 2011-02-02 中国人民解放军63653部队 基于虚拟源原理的航空放射性测量系统刻度装置
TW201115176A (en) * 2009-10-30 2011-05-01 Iner Aec Executive Yuan Method for determining penetration effect and detecting efficiency and calibrating detecting efficiency of crystal
CN103135122A (zh) * 2011-12-01 2013-06-05 中国辐射防护研究院 混合核素γ点源体积样品效率刻度方法
CN103605149A (zh) * 2013-11-26 2014-02-26 西北核技术研究所 一种刻度氙气样品HPGe探测效率的装置及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004245649A (ja) * 2003-02-12 2004-09-02 Central Res Inst Of Electric Power Ind バックグラウンド計数率の推定方法及びバックグラウンド計数率を補正する測定装置
TW201115176A (en) * 2009-10-30 2011-05-01 Iner Aec Executive Yuan Method for determining penetration effect and detecting efficiency and calibrating detecting efficiency of crystal
CN201732166U (zh) * 2010-06-12 2011-02-02 中国人民解放军63653部队 基于虚拟源原理的航空放射性测量系统刻度装置
CN103135122A (zh) * 2011-12-01 2013-06-05 中国辐射防护研究院 混合核素γ点源体积样品效率刻度方法
CN103605149A (zh) * 2013-11-26 2014-02-26 西北核技术研究所 一种刻度氙气样品HPGe探测效率的装置及方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
S. MAHLING 等: "The dependence of the virtual point-detector on the HPGe detector dimensions", 《NUCLEAR INSTRUMENTS AND METHODS IN PHYSICS RESEARCH》, vol. 557, no. 2, 15 February 2006 (2006-02-15), XP 025055020, DOI: doi:10.1016/j.nima.2005.11.118 *
景越峰 等: "减小电子涨落的指向概率虚拟源方法", 《强激光与粒子束》 *
熊文彬 等: "基于虚拟点探测器模型的HPGe探测器体源效率刻度", 《原子能科学技术》, vol. 45, no. 8, 20 August 2011 (2011-08-20) *
田自宁 等: "HPGe探测器体源效率面源模拟刻度技术", 《核电子学与探测技术》, vol. 33, no. 2, 20 February 2013 (2013-02-20) *
邹长贵: "面源模式中后退虚拟点源位置的固定化研究", 《辐射防护通讯》 *
陈晓秋: "用于计算较大尺度面源气载放射性释放的环境辐射剂量的一种套网格模型", 《辐射防护》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107195342A (zh) * 2017-04-05 2017-09-22 中广核研究院有限公司 一种测量核电厂核素活度分布的方法
CN107195342B (zh) * 2017-04-05 2019-07-30 中广核研究院有限公司 一种测量核电厂核素活度分布的方法
CN109814147A (zh) * 2019-01-01 2019-05-28 中国人民解放军63653部队 使用探测器加衰减体测量小孔周围γ放射性核素分布的方法
CN109814147B (zh) * 2019-01-01 2021-01-29 中国人民解放军63653部队 使用探测器加衰减体测量小孔周围γ放射性核素分布的方法
CN110888155A (zh) * 2019-10-11 2020-03-17 中国人民解放军陆军防化学院 一种大体积液态放射源对γ谱仪有效贡献范围确定方法
CN113703034A (zh) * 2021-07-14 2021-11-26 中国人民解放军63653部队 反演污染区虚拟点位置和边界的技术方法及设备
CN113885070A (zh) * 2021-09-13 2022-01-04 中国人民解放军63653部队 确定均匀分布和指数分布体源虚拟点位置的方法和设备

Also Published As

Publication number Publication date
CN104020507B (zh) 2017-06-09

Similar Documents

Publication Publication Date Title
CN106199676B (zh) 一种伽玛探测器无源效率刻度方法
CN104020507A (zh) 基于虚拟点探测器原理的放射性区域确定方法
CN103163548B (zh) 基于伽马相机的放射性物质探测方法及其装置和系统
CN104020483B (zh) 采用虚拟源模拟刻度探测器探测效率的方法及装置
CN104536029B (zh) 基于多个NaI(T1)探测器的放射性物质二维定位方法
US10838078B2 (en) Method for searching for and detecting gamma radiation sources
CN104483698B (zh) 基于点源效率函数的氙气源虚拟点源刻度装置及方法
Anderson et al. Mobile robotic radiation surveying with recursive Bayesian estimation and attenuation modeling
CN112051601B (zh) 一种基于虚拟源原理的源边界参数蒙特卡罗反演技术
US9423361B2 (en) Inner image generating apparatus and method thereof
CN104483693A (zh) 一种非均匀分布源探测效率计算及模拟装置及方法
CN104570035B (zh) 一种放射性气体核素β射线自吸收校正方法
CN109521456B (zh) 一种基于正则化最小二乘法的γ辐射源项反演方法及系统
KR20160081307A (ko) 방사선원의 3차원 영상 처리 장치 및 방법
CN103727962B (zh) 大视场红外光电经纬仪精度标定方法
JP6677568B2 (ja) 放射能分布解析システムおよび放射能分布解析方法
CN107783174B (zh) 一种核设施液态流出物在线监测仪探测效率的校准方法
CN203025355U (zh) 一种车载伽玛能谱大气氡消除装置
CN105300385A (zh) 一种线阵红外地球敏感器的四通道姿态解算方法
CN104501785B (zh) 基于电子经纬仪的精测系统的磁试验室磁场基准建立方法
Kowash A rotating modulation imager for the orphan source search problem
CN113447974A (zh) 一种放射源强度三维分布的估计方法
Aron et al. Recovery of the Surface Density of Radionuclide Contamination of Soil from Aerogamma Survey Data
Falkner et al. Modeling minimum detectable activity as a function of detector speed
Arinc et al. Defined solid angle alpha counting at NPL

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant