CN103176201A - 放射性物质巡检定位方法及设备 - Google Patents
放射性物质巡检定位方法及设备 Download PDFInfo
- Publication number
- CN103176201A CN103176201A CN2011104391994A CN201110439199A CN103176201A CN 103176201 A CN103176201 A CN 103176201A CN 2011104391994 A CN2011104391994 A CN 2011104391994A CN 201110439199 A CN201110439199 A CN 201110439199A CN 103176201 A CN103176201 A CN 103176201A
- Authority
- CN
- China
- Prior art keywords
- radiation intensity
- radiomaterial
- value
- radioactive radiation
- subregion
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000007689 inspection Methods 0.000 title abstract description 9
- 239000000941 radioactive substance Substances 0.000 title abstract 5
- 230000005855 radiation Effects 0.000 claims abstract description 114
- 230000002285 radioactive effect Effects 0.000 claims abstract description 85
- 238000001514 detection method Methods 0.000 claims abstract description 30
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 230000004807 localization Effects 0.000 claims description 27
- 238000000354 decomposition reaction Methods 0.000 claims description 7
- 230000007613 environmental effect Effects 0.000 abstract description 4
- 238000005070 sampling Methods 0.000 abstract description 2
- 239000000523 sample Substances 0.000 abstract 2
- 238000005259 measurement Methods 0.000 description 8
- 230000004888 barrier function Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 230000000903 blocking effect Effects 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
- G01V5/20—Detecting prohibited goods, e.g. weapons, explosives, hazardous substances, contraband or smuggled objects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
-
- 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/16—Measuring radiation intensity
- G01T1/169—Exploration, location of contaminated surface areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种放射性物质巡检定位方法及设备。所述方法包括步骤:提供环境的本底放射性辐射强度值;在巡检路径上的多个采样点用探测器采集来自探测区域的放射性辐射强度值;根据所采集的放射性辐射强度值和本底放射性辐射强度值计算探测区域的放射性辐射强度分布;和根据所述放射性辐射强度分布确定放射性物质所在位置。所述方法和设备能够基于探测器在巡检路线上的多点观测来得出探测区域的放射性物质的位置和放射性强度分布。
Description
技术领域
本发明涉及一种放射性物质巡检定位方法及用于实现该放射性物质巡检定位方法的设备。
背景技术
环保、安全等相关领域在监测和管理放射性物质时,往往需要测量某区域内的放射性分布情况或寻找丢失的放射源。移动式巡检谱仪是一种可以安装在汽车、直升机、船舶或其它载体上的放射性物质探测设备,自身带有定位系统,工作时可以实时获取当前的地理位置和放射性物质的伽马能谱或中子计数率。移动式巡检谱仪除了普通的谱仪所具有的探测器、电源、前放、电子学模块、多道和数据处理单元外,还包括GPS定位装置。另外,移动式巡检谱仪更注重能高效的对大量采集数据进行获取、保存、读取和分析。通过传统的移动式谱仪,可以测量行进路线上的放射性水平,但得不到路线所在局部区域的放射性分布,例如,当车载式的巡检谱仪沿着公路检测时,只能标记出公路上各个点的放射性水平,但不能得到公路两侧一定范围内的放射性分布的详细情况。如果是在寻找丢失的放射源,则不能知道放射源的具体位置。
也有一些专门用来确定放射源位置的设备。一类是通过成像来定位,例如伽马相机。一类是带有准直器和旋转装置的设备,通过准直器可以判断放射源所在的方向,利用多个位置的测量就可以确定放射源的位置。然而由于需要长时间定点测量和探测效率低等原因,这些设备都不适用于巡检,只能在巡检谱仪确定了可疑区域后,才使用它们进一步测量。
为此,需要一种能够迅速有效地探测巡检路径周围的探测区域的放射性辐射强度分布和确定放射性物质的所在位置的巡检方法和巡检设备。
发明内容
本发明的目的之一是提供一种能够基于探测器在巡检路线上的多点观测来得出探测区域的放射性物质的位置和放射性强度分布的巡检方法和巡检设备。
为了实现上述发明目的,本发明的技术方案通过以下方式来实现:
根据本发明的一个方面,提供一种放射性物质巡检定位方法,包括步骤:
(a)提供环境的本底放射性辐射强度值;
(b)在巡检路径上的多个采样点用探测器采集来自探测区域的放射性辐射强度值;
(c)根据所采集的放射性辐射强度值和本底放射性辐射强度值计算探测区域的放射性辐射强度分布;
(d)根据所述放射性辐射强度分布确定放射性物质所在位置。
优选地,所述方法在步骤(b)和步骤(c)之间还可以包括步骤(b1):
判定所采集到的放射性辐射强度值与本底放射性辐射强度值之差是否大于采集阈值,如果大于采集阈值,则执行步骤(c);如果不大于采集阈值,则通过将所采集到的放射性辐射强度值与本底放射性辐射强度值进行加权平均来对本底放射性辐射强度值进行更新,并返回步骤(b)重新采集放射性辐射强度值。
优选地,所述方法在步骤(b)和步骤(c)之间还可以包括步骤(b2):
判定在连续多个采样点处所采集到的放射性辐射强度值序列中是否存在先突降再突升的数据段,如果存在所述数据段,则删除所述数据段以消除障碍物的影响。
优选地,所述探测器可以包括朝向不同方向的探测区域的至少两个探测器,所述探测器之间由屏蔽部件隔开。
优选地,所述放射性辐射强度值可以由计数率或剂量率来表示。
进一步地,所述步骤(c)可以包括步骤(c1):
将探测区域分成多个子区域,提供所述探测器对于所述子区域的探测效率,并根据所述本底放射性辐射强度和所述探测效率确定每个子区域的放射性辐射强度。
进一步地,所述子区域的个数可以不超过采样点的个数。
进一步地,所述探测效率可以针对于所探测的辐射能量范围和放射性物质的种类进行标定。
进一步地,所述步骤(c)可以包括步骤(c2):
基于每个子区域的放射性辐射强度在探测区域上绘制放射性辐射强度分布图。
进一步地,所述子区域的划分可以在一维、二维或三维方向上进行。
进一步地,所述探测器设置在移动载体上,所采集到的放射性辐射强度值
其中m是子区域的标号,m=1,…,M,M为子区域的总数,bt是所述本底放射性辐射强度值,pm表示标号为m的子区域的位置,am是标号为m的子区域的放射性辐射强度值,qt是所述移动载体所处的位置,θt是所述移动载体的方向角,Φ(pm,qt,θt)是探测器在所述移动载体处于位置qt、方向角为θt时对pm位置的探测效率,t是采集到的放射性辐射强度值的标号,t=1,…,N,N为采集到的放射性辐射强度值的个数,εt是采集误差。
更进一步地,所述步骤(c1)可以通过向量C在向量集合D上的稀疏分解来实现,其中向量C=(c1-b1,…,cN-bN)T,向量集合D={Dm},
Dm=(Φ(pm,q1,θ1),…,Φ(pm,qN,θN))T。
再进一步地,所述确定每个子区域的放射性辐射强度的步骤还可以包括:
(c11)给定标号集合S和残差向量R,并设定S的初始值S0为空集,将残差向量R的初始值R0设定为向量C,将am的初始值设定成0;
(c12)将R0在向量集合D中的向量Dm上的投影定义为Pm,将Pm取最大值时对应的m值定义成n,将n添加到S0中以对标号集合S进行更新,并对残差向量R进行更新,令更新后的残差向量R为向量C与向量C在向量Dn上的投影之差;
(c13)判定更新后的残差向量R的范数是否小于残差阈值,如果小于残差阈值,则将向量C在向量Dm上的投影系数作为am的计算结果,其中仅当m的取值为所得的更新后的标号集合S中的元素时,am的计算结果不为零,而当m取其它值时,am的计算结果为零;如果更新后的残差向量R的范数不小于残差阈值,则将更新后的残差向量R和更新后的标号集合S分别替代R0和S0重新执行步骤(c12)以反复更新向量残差向量R和标号集合S直至更新后的残差向量R的范数小于残差阈值或重新执行步骤(c12)的循环次数达到了设定值为止。
再进一步地,所述步骤(d)可以包括步骤(d1):
将am的计算结果不为零的子区域确定为有源子区域,通过有源子区域的分布和am的计算结果来确定放射性物质所在的子区域。
再进一步地,所述步骤(d1)可以包括:
利用数值方法求出在误差函数取值最小时的有源子区域的位置,其中所述误差函数为
其中l=x1,…,xL,L为有源子区域的个数,x1,…,xL为有源子区域在所有子区域中的标号,所述在误差函数H取值最小时的pl即是放射性物质更准确的位置点。
根据本发明的另一方面,还提供了一种用于实现上述任一种放射性物质巡检定位方法的放射性物质巡检定位设备,包括:
两个或更多探测器,用于在巡检路径上的多个采样点中的每个采样点处采集来自巡检路径周围的探测区域的放射性辐射强度值;和
移动载体,用于承载所述探测器沿所述巡检路径移动经过所述采样点,
其中所述探测器朝向采样点周围的不同方向,彼此由屏蔽部件所隔离。
本发明的上述技术方案中的至少一个方面能够基于探测器在不同方向和距离上的探测效率的差异性通过在巡检路径上的多点观测来估计放射性物质的存在位置和分布情况。尤其为放射性物质的动态巡检提供了一种迅捷有效的途径。
附图说明
图1示出根据本发明的实施例的放射性物质巡检定位方法的示意性流程图;
图2示出根据本发明的实施例的放射性物质巡检定位设备的巡检路线和探测区域的示意图;
图3示出根据本发明的放射性物质巡检定位方法中采集到的放射性辐射强度分布的示意图;和
图4示出本发明的实施例的具有四个探测器的放射性物质巡检定位设备的示意图;和
图5示出本发明的实施例的具有两个探测器的放射性物质巡检定位设备的示意图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。在说明书中,相同或相似的附图标号表示相同或相似的部件。下述参照附图对本发明实施方式的说明旨在对本发明的总体发明构思进行解释,而不应当理解为对本发明的一种限制。
图1示意性地示出根据本发明的一实施例的一种放射性物质巡检定位方法。所述方法包括步骤:提供环境的本底放射性辐射强度值的步骤101、在巡检路径上的多个采样点用探测器采集来自探测区域的放射性辐射强度值的步骤102、根据所采集的放射性辐射强度值计算探测区域的放射性辐射强度分布的步骤105和根据所述放射性辐射强度分布确定放射性物质所在位置的步骤106。图1中虚线表示的流程步骤是可选的步骤。
在实际操作中,即使在环境中没有放射性物质的情况下,探测器也会显示出一定的探测信号值,即本底值。在探测开始前,需要首先提供本底放射性辐射强度值。该本底值的采集与有放射性物质存在的一个可能的差异是:本底值从各个方向测量都没有明显的变化;但由于放射性物质一般分布不均匀,如果探测设备各个方向的探测效率不相同,那么在各个方向上所得的测量结果会有明显的变化。可以利用这一特点判断采集的初始本底值是否有效。另一个方法是采集几个不同位置的本底值,如果差异不大,说明周围有放射性物质的可能性很小,从而也可以作为初始本底值。初始本底值是否有效也可以通过经验和以往的测量数据判断。
在本发明的实施例中,在巡检路径上可以设置有多个采样点。在提供环境的本底放射性辐射强度值之后,用探测器在所述多个采样点处采集来自探测区域的放射性辐射强度值。在巡检过程中可以以一定的时间间隔不断采样。
根据所采集的放射性辐射强度值和本底放射性辐射强度值来计算探测区域的放射性辐射强度分布。在该步骤105中,可以将探测区域分成多个子区域,提供所述探测器在各数据采集点对于所述子区域的探测效率,并根据所述本底放射性辐射强度和所述探测效率确定每个子区域的放射性辐射强度。所述子区域的划分可以在一维、二维或三维方向上进行。也就是说,每个子区域可以被构建成一维的线段、二维的平面区域或三维的空间体。图2示出了在二维平面中将探测区域划分成一个个方格形状的子区域的示例。从图2中可见,探测区域覆盖了巡检路线(即探测器的移动载体的行驶路线)及其两侧的大片区域。
探测效率是探测到的辐射粒子数与从辐射源发出的该粒子数之比。探测效率是放射性物质巡检设备自身的属性,但其与测量角度、距离,甚至放射性物质的能量范围和放射源类型都是相关的。因此,在探测或计算放射性物质的辐射强度分布之前,需要先标定或计算出探测效率,例如可以通过实际测量、模拟计算或两者相结合的方法来求出探测效率。
子区域的个数涉及探测区域的分辨率,受到探测器采样率以及探测器的移动载体的移动速度的限制。例如,可以将子区域的个数限定成不超过采样点的个数。
假定共有M个子区域,令子区域的标号为m,有m=1,…,M。令标号为m的子区域的放射性辐射强度值为am,标号为m的子区域的位置为pm。探测器所采集到的辐射强度值实际上是各个子区域的辐射强度的加权和。假定探测器共采集到N个放射性辐射强度值,令采集到的放射性辐射强度值的标号为t,有t=1,…,N,本底放射性辐射强度值为bt。假定探测器设置在移动载体上,移动载体所处的位置为qt,移动载体的方向角为θt,探测器在所述移动载体处于位置qt、方向角为θt时对pm位置的探测效率为Φ(pm,qt,θt)。则所采集到的放射性辐射强度值可以为
其中εt是采集误差。如上所述,式(1)中的探测效率Φ(pm,qt,θt)可以通过事先标定或插值计算等方式获得,在计算中为已知量。
从式(1)可以看出,计算每个子区域的放射性辐射强度实质上就是求出针对于各个子区域的am的值。探测器工作正常的情况下,采集误差εt相对较小。于是对am的计算实际上变成了求解线性方程组的问题。然而,在不另附加约束的情况下,am的解往往不唯一。考虑到巡检放射性物质的实际情况,放射性物质源一般不会太多,且位置相对集中,所以为式(1)的求解增加了“含有放射性物质的子区域最少”的约束条件。所谓“含有放射性物质的子区域”,即am的值大于零或大于某一阈值的子区域,这样的子区域可以被称为有源子区域。于是,问题变成了求出式(1)所表示的方程组中的am的解中零值最多的解。在测量点有限的情况下,这样的假定有利于发现放射性辐射强度最高的若干个子区域。
根据上述约束条件对式(1)进行求解的方法很多,但考虑到根据本发明的巡检定位方法和设备往往需要进行实时处理,因此,优选计算复杂度相对较低的方法来进行计算。例如,可以通过向量的稀疏分解来求解。假定向量C=(c1-b1,…,cN-bN)T,向量集合D={Dm},Dm=(Φ(pm,q1,θ1),…,Φ*pm,qN,θN))T,则可以通过向量C在向量集合D上的稀疏分解来求出所需的am的解。
同样,向量的稀疏分解的方法也有很多,例如Matching Pursuit方法等。在以下的实施例中,仅给出一种示例,但本发明的实施例不限于此。
在该实施例中,首先给定标号集合S和残差向量R,并设定S的初始值S0为空集,将残差向量R的初始值R0设定为向量C,将am的初始值设定成0。
接着,将R0在向量集合D中的向量Dm上的投影定义为Pm,将Pm取最大值时所对应的m值定义成n,将n添加到S0中以对标号集合S进行更新,并对残差向量R进行更新,令更新后的残差向量R为向量C与向量C在向量Dn上的投影之差。经过该步骤,标号集合S中的元素从零个增加到1个。
然后,判定更新后的残差向量R的范数是否小于残差阈值,如果小于残差阈值,则将向量C在向量Dm上的投影系数作为am的计算结果。在这种情况下,仅当m的取值为所得的更新后的标号集合S中的元素时,am的计算结果才不为零,而当m取其它值时,am的计算结果为零。反之,如果更新后的残差向量R的范数不小于残差阈值,则将更新后的残差R和更新后的标号集合S分别替代R0和S0重新执行上一步骤以反复更新残差向量R和标号集合S直至更新后的残差向量R的范数小于残差阈值为止。其中,残差阈值可以根据需要来指定,残差阈值越大,则运算的循环次数越少,得出的解的不为零的am值也越少,反之,残差阈值越小,运算的循环次数越多,得出的解的不为零的am值也越多,但是太小的残差阈值可能导致解的不稳定。也可以根据需要给出循环次数的设定值,当反复更新残差向量R和标号集合S进行循环运算的循环次数到达该设定值时,即终止该运算过程。
经过上述计算,假定最终获得的标号集合S包含L个元素,用x1,…,xL来表示,则意味着在集合{am}中,仅当m=x1,…,xL时,am值才不为零。为此,可以将这些am的值不为零的子区域确定为有源子区域。当然,也可以为am的值再设定一个大于零的阈值,而仅将am的值大于该阈值的子区域才确定成有源子区域。在上述步骤106中,可以基于所求出的每个子区域的放射性辐射强度值am在探测区域上绘制放射性辐射强度分布图来确定放射性辐射物质可能存在的子区域位置。
在上述步骤106中,在通过向量的稀疏分解求出子区域的放射性辐射强度值am之后,还可以进一步计算放射性物质所在的精确位置,尤其是在有源子区域较少且所巡检的放射性物质为点放射源时更是如此。该计算可以利用数值方法进行,给定误差函数
其中L为有源子区域的个数,x1,…,xL为有源子区域在所有子区域中的标号;t=1,…,N,是采集到的放射性辐射强度值的标号,利用数值方法求出在误差函数H取值最小时的pl即是放射性物质更准确的位置点。需要指出的是,在对式(2)进行计算时,不仅pl是未知量,有源子区域的放射性辐射强度值al也需要当作未知量来处理,因此,
不再是线性方程。
在上述利用向量的稀疏分解进行求解的过程中,向量集合D={Dm}中的向量对求解式(1)有较大影响,Dm之间的差异越大,即相关性越小,则式(1)所表示的方程组的病态性就越小,求得的结果的误差也越小,也越稳定。又因为Dm=(Φ(pm,q1,θ1),…,Φ(pm,qN,θN))T,即探测器在每个采样点处对同一子区域(同一位置)的探测效率。所以,如果探测器在各个方向上的探测效率的差异性越大,则Dm之间的相关性就越小,式(1)所表示的方程组的病态性就越小。
为了在不明显降低整体探测效率的基础上增加探测器在各个方向上的探测效率的差异性,可以安装多个探测器,并在彼此之间加装屏蔽材料,如图4和图5所示。例如,可以将探测器设置成包括朝向采样点四周的不同方向的两个、三个、四个或更多探测器。在本发明的实施例中,探测器可以为NaI探测器、中子探测器或其它探测器。
在一些实施例中,可能需要对不同能量区间或只对感兴趣的某种放射性物质进行定位,而对于不同能量或放射性物质,探测器的探测效率也有所不同。在此情况下,所述探测效率可以针对于所探测的辐射能量范围和放射性物质的种类进行标定。这样可以更加准确地估计不同能量或感兴趣的某种放射性物质在探测区域中的分布或放射源的位置。
如图1所示,在上述步骤102和步骤105之间还可以包含两个可选的步骤103和104。在步骤103中,判定所采集到的放射性辐射强度值与本底放射性辐射强度值之差是否大于采集阈值,如果大于采集阈值,则表明探测区域中明显存在放射性物质,需要对其位置和强度分布进行计算,因而向下执行步骤105或104;如果不大于采集阈值,则表明探测区域中没有明显的放射性物质的辐射,因此转到步骤107,在步骤107中,通过将所采集到的放射性辐射强度值与本底放射性辐射强度值进行加权平均来对本底放射性辐射强度值进行更新,更新后的本底放射性辐射强度值替代了最初在步骤101中提供的本底放射性辐射强度值,并返回步骤102重新对放射性辐射强度值进行采集。本领域技术人员可以根据诸如辐射的类型、强度、能量等因素设定合适的采集阈值,以满足实际巡检的需要。通过对本底放射性辐射强度值的持续更新,可以得出每个采样点附近的更加准确的本底放射性辐射强度值,从而能够防止由于不同的采样点之间的环境本底放射性的差异而带来的探测误差。
在一些实施例中,采样点的周围的环境可能是复杂的,例如存在能够使放射性辐射强度产生明显衰减的建筑物、山坡等障碍物。在这种情况下,为了消除由于障碍物的遮挡所导致的探测误差,可以对探测得到的数据进行预处理去除掉这些受到障碍物的遮挡的影响的数据。于是在上述步骤102和步骤105之间可以设置可选的步骤104。如果巡检设备经过被障碍物所遮挡的区域,则在进入该被遮挡的区域时,采集到的辐射强度值会突然下降,而在离开该被遮挡的区域时,采集到的辐射强度值又会突然上升。图3示出了存在障碍物遮挡情况的探测器采集的放射性辐射强度曲线,其中,放射性辐射强度由计数率来表示。在图3中,区域2和区域4即是受到障碍物遮挡的区域,其辐射强度相比于相邻的区域都显现出先突然下降而后突然上升的特点。因此在步骤104中,检查在连续多个采样点处所采集到的放射性辐射强度值序列中是否存在先突降再突升的数据段,并由此判定探测器是否被障碍物所阻挡。如果存在所述数据段,则删除所述数据段以消除障碍物的影响。
在图1中用虚线表示的上述两个可选的步骤103和104可以同时存在,也可以存在其中的任意一个,这两个可选的步骤的顺序也可以调换。
在本发明的一些实施例中,可能需要巡检多个探测区域或进行反复多次确定放射性物质的分布和位置估计,在这种情况下,可以在图1中的步骤106后面增加一判定步骤108,如果需要反复确定放射性物质的分布和位置估计,则返回步骤102重新采集检测数据继续流程,否则结束流程。当然,也可以不设置步骤108,而在执行步骤106以确定放射性物质所在位置之后直接结束整个流程。
在本发明的实施例中,放射性辐射强度值可以由计数率或剂量率等在本领域中已知的能够恒量辐射强度的参数来表示。本发明的上述实施例中的特征可以相互组合以形成落入本发明的保护范围的新的实施例,除非特征之间存在明显的相互矛盾。
除去上述的放射性物质巡检定位方法之外,本发明还涉及一种用于实现上述任一放射性物质巡检定位方法的放射性物质巡检定位设备,该巡检定位设备包括:探测器和移动载体。探测器用于在巡检路径上的多个采样点中的每个采样点处采集来自巡检路径周围的探测区域的放射性辐射强度值。而移动载体,例如汽车,用于承载所述探测器沿所述巡检路径移动经过所述采样点。所述探测器可以具有两个、三个、四个或更多个,且可以分别朝向采样点周围的不同方向,并彼此由屏蔽部件所隔离,例如由屏蔽板隔开。
虽然结合附图对本发明进行了说明,但是附图中公开的实施例旨在对本发明优选实施方式进行示例性说明,而不能理解为对本发明的一种限制。
虽然本发明总体构思的一些实施例已被显示和说明,本领域普通技术人员将理解,在不背离本总体发明构思的原则和精神的情况下,可对这些实施例做出改变,本发明的范围以权利要求和它们的等同物限定。
Claims (16)
1.一种放射性物质巡检定位方法,包括步骤:
(a)提供环境的本底放射性辐射强度值;
(b)在巡检路径上的多个采样点用探测器采集来自探测区域的放射性辐射强度值;
(c)根据所采集的放射性辐射强度值和本底放射性辐射强度值计算探测区域的放射性辐射强度分布;
(d)根据所述放射性辐射强度分布确定放射性物质所在位置。
2.根据权利要求1所述的放射性物质巡检定位方法,其特征在于,所述方法在步骤(b)和步骤(c)之间还包括步骤(b1):
判定所采集到的放射性辐射强度值与本底放射性辐射强度值之差是否大于采集阈值,如果大于采集阈值,则执行步骤(c);如果不大于采集阈值,则通过将所采集到的放射性辐射强度值与本底放射性辐射强度值进行加权平均来对本底放射性辐射强度值进行更新,并返回步骤(b)重新采集放射性辐射强度值。
3.根据权利要求1所述的放射性物质巡检定位方法,其特征在于,所述方法在步骤(b)和步骤(c)之间还包括步骤(b2):
判定在连续多个采样点处所采集到的放射性辐射强度值序列中是否存在先突降再突升的数据段,如果存在所述数据段,则删除所述数据段以消除障碍物的影响。
4.根据权利要求1所述的放射性物质巡检定位方法,其特征在于,所述探测器包括分别朝向不同方向的探测区域的至少两个探测器,所述探测器之间由屏蔽部件隔开。
5.根据权利要求1所述的放射性物质巡检定位方法,其特征在于,所述放射性辐射强度值由计数率或剂量率来表示。
6.根据权利要求1-5中任一项所述的放射性物质巡检定位方法,其特征在于,所述步骤(c)包括步骤(c1):
将探测区域分成多个子区域,提供所述探测器对于所述子区域的探测效率,并根据所述本底放射性辐射强度和所述探测效率确定每个子区域的放射性辐射强度。
7.根据权利要求6所述的放射性物质巡检定位方法,其特征在于,所述子区域的个数不超过采样点的个数。
8.根据权利要求6所述的放射性物质巡检定位方法,其特征在于,所述探测效率针对于所探测的辐射能量范围和放射性物质的种类进行标定。
9.根据权利要求6所述的放射性物质巡检定位方法,其特征在于,所述步骤(c)包括步骤(c2):
基于每个子区域的放射性辐射强度在探测区域上绘制放射性辐射强度分布图。
10.根据权利要求6所述的放射性物质巡检定位方法,其特征在于,所述子区域的划分在一维、二维或三维方向上进行。
11.根据权利要求6所述的放射性物质巡检定位方法,其特征在于,所述探测器设置在移动载体上,所采集到的放射性辐射强度值
其中m是子区域的标号,m=1,…,M,M为子区域的总数,bt是所述本底放射性辐射强度值,pm表示标号为m的子区域的位置,am是标号为m的子区域的放射性辐射强度值,qt是所述移动载体所处的位置,θt是所述移动载体的方向角,Φ(pm,qt,θt)是探测器在所述移动载体处于位置qt、方向角为θt时对pm位置的探测效率,t是采集到的放射性辐射强度值的标号,t=1,…,N,N为采集到的放射性辐射强度值的个数,εt是采集误差。
12.根据权利要求11所述的放射性物质巡检定位方法,其特征在于,所述步骤(c1)通过向量C在向量集合D上的稀疏分解来实现,其中向量C=(c1-b1,…,cN-bN)T,向量集合D={Dm},
Dm=(Φ(pm,q1,θ1),…,Φ(pm,qN,θN))T。
13.根据权利要求12所述的放射性物质巡检定位方法,其特征在于,所述确定每个子区域的放射性辐射强度的步骤还包括:
(c11)给定标号集合S和残差向量R,并设定S的初始值S0为空集,将残差向量R的初始值R0设定为向量C,将am的初始值设定成0;
(c12)将R0在向量集合D中的向量Dm上的投影定义为Pm,将Pm取最大值时对应的m值定义成n,将n添加到S0中以对标号集合S进行更新,并对残差向量R进行更新,令更新后的残差向量R为向量C与向量C在向量Dn上的投影之差;
(c13)判定更新后的残差向量R的范数是否小于残差阈值,如果小于残差阈值,则将向量C在向量Dm上的投影系数作为am的计算结果,其中仅当m的取值为所得的更新后的标号集合S中的元素时,am的计算结果不为零,而当m取其它值时,am的计算结果为零;如果更新后的残差向量R的范数不小于残差阈值,则将更新后的残差向量R和更新后的标号集合S分别替代R0和S0重新执行步骤(c12)以反复更新向量残差向量R和标号集合S直至更新后的残差向量R的范数小于残差阈值或重新执行步骤(c12)的循环次数达到了预定值为止。
14.根据权利要求13所述的放射性物质巡检定位方法,其特征在于,所述步骤(d)包括步骤:
(d1)将am的计算结果不为零的子区域确定为有源子区域,通过有源子区域的分布和am的计算结果来确定放射性物质所在的子区域。
15.根据权利要求14所述的放射性物质巡检定位方法,其特征在于,所述步骤(d1)包括:
利用数值方法求出在误差函数取值最小时的有源子区域的位置,其中所述误差函数为
其中l=x1,…,xL,L为有源子区域的个数,x1,…,xL为有源子区域在所有子区域中的标号,所述在误差函数H取值最小时的pl即是放射性物质更准确的位置点。
16.一种用于实现根据权利要求1-15中任一项所述的放射性物质巡检定位方法的放射性物质巡检定位设备,包括:
两个或更多探测器,用于在巡检路径上的多个采样点中的每个采样点处采集来自巡检路径周围的探测区域的放射性辐射强度值;和
移动载体,用于承载所述探测器沿所述巡检路径移动经过所述采样点,
其中所述探测器朝向采样点周围的不同方向,彼此由屏蔽部件所隔离。
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110439199.4A CN103176201B (zh) | 2011-12-23 | 2011-12-23 | 放射性物质巡检定位方法及设备 |
RU2014129589/28A RU2575582C2 (ru) | 2011-12-23 | 2012-12-20 | Способ и устройство для патрульной инспекции и локализации радиоактивного вещества |
PCT/CN2012/087012 WO2013091552A1 (zh) | 2011-12-23 | 2012-12-20 | 放射性物质巡检定位方法及设备 |
GB1412789.8A GB2512547B (en) | 2011-12-23 | 2012-12-20 | Method and device for patrol inspecting and locating a radioactive substance |
EP12858884.5A EP2796898B1 (en) | 2011-12-23 | 2012-12-20 | Method and device for inspecting and positioning radioactive substance |
US14/367,824 US9829601B2 (en) | 2011-12-23 | 2012-12-20 | Method and device for patrol inspecting and locating a radioactive substance |
HK13113202.2A HK1185950A1 (zh) | 2011-12-23 | 2013-11-26 | 放射性物質巡檢定位方法及設備 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110439199.4A CN103176201B (zh) | 2011-12-23 | 2011-12-23 | 放射性物质巡检定位方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103176201A true CN103176201A (zh) | 2013-06-26 |
CN103176201B CN103176201B (zh) | 2015-03-25 |
Family
ID=48636182
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110439199.4A Active CN103176201B (zh) | 2011-12-23 | 2011-12-23 | 放射性物质巡检定位方法及设备 |
Country Status (6)
Country | Link |
---|---|
US (1) | US9829601B2 (zh) |
EP (1) | EP2796898B1 (zh) |
CN (1) | CN103176201B (zh) |
GB (1) | GB2512547B (zh) |
HK (1) | HK1185950A1 (zh) |
WO (1) | WO2013091552A1 (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104237959A (zh) * | 2014-08-21 | 2014-12-24 | 北京辛耕普华医疗科技有限公司 | 对放射性物质实时动态追踪定位的方法和设备 |
CN107977543A (zh) * | 2016-10-19 | 2018-05-01 | 中国辐射防护研究院 | 放射性物品运输弥散事故中集体剂量的计算方法 |
CN108445527A (zh) * | 2018-01-30 | 2018-08-24 | 奕瑞新材料科技(太仓)有限公司 | 一种射线源位置坐标的获取方法 |
CN109490939A (zh) * | 2018-11-23 | 2019-03-19 | 南京航空航天大学 | 一种应用于小型无人机的放射源定位方法 |
CN110850463A (zh) * | 2019-10-09 | 2020-02-28 | 岭澳核电有限公司 | 一种核电站放射性监测方法、装置、控制终端及监测系统 |
CN112147163A (zh) * | 2020-08-20 | 2020-12-29 | 中国原子能科学研究院 | 一种用于多点测量脉冲萃取柱钚浓度的串联探测装置 |
CN112357606A (zh) * | 2020-11-12 | 2021-02-12 | 重庆建安仪器有限责任公司 | 一种基于辐射监测的辐射标识器自动投放控制系统 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6677568B2 (ja) * | 2016-04-18 | 2020-04-08 | 日立Geニュークリア・エナジー株式会社 | 放射能分布解析システムおよび放射能分布解析方法 |
ES2681094B1 (es) | 2017-02-10 | 2019-10-21 | Consejo Superior Investigacion | Sistema y metodo para la identificacion volumetrica e isotopica de distribuciones de escenas radioactivas |
CN108510827A (zh) * | 2018-04-12 | 2018-09-07 | 重庆建安仪器有限责任公司 | 一种模拟训练器 |
GB201807560D0 (en) * | 2018-05-09 | 2018-06-20 | Univ Cape Town | Neutron detector |
CN111914877B (zh) * | 2020-06-15 | 2021-09-28 | 国网安徽省电力有限公司电力科学研究院 | 省级移动巡检与地市级移动巡检的系统融合接入方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH116875A (ja) * | 1997-06-16 | 1999-01-12 | Mitsubishi Electric Corp | シンチレーションファイバを用いた放射線検出方法 |
CN1598553A (zh) * | 2003-09-18 | 2005-03-23 | 清华大学 | 一种对车辆放射性物质定位的监测方法及其装置 |
CN2690896Y (zh) * | 2003-09-03 | 2005-04-06 | 广州杰赛科技股份有限公司 | 数字化辐射巡测系统 |
CN1318841C (zh) * | 2002-11-06 | 2007-05-30 | 美国科技工程公司 | 移动式x射线反向散射检查车 |
CN101069091A (zh) * | 2004-08-12 | 2007-11-07 | 纳沃特克医药有限公司 | 对象身体内的放射性源定位 |
WO2007131536A1 (en) * | 2006-05-11 | 2007-11-22 | Advanced Accelerator Applications | Method and device of detecting, locating and/or analyzing a radioactive source(s) in a material, e.g. a biological tissue |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6054712A (en) * | 1998-01-23 | 2000-04-25 | Quanta Vision, Inc. | Inspection equipment using small-angle topography in determining an object's internal structure and composition |
JP2002006053A (ja) * | 2000-06-19 | 2002-01-09 | Japan Nuclear Cycle Development Inst States Of Projects | 指向型放射線検出装置 |
US7835853B2 (en) * | 2003-02-19 | 2010-11-16 | Sumitomo Electric Industries, Ltd. | Vehicle detection system |
DE10337877A1 (de) * | 2003-08-18 | 2005-03-17 | Basf Ag | Verfahren zur Detektion der durch einen Umwelteinfluss hervorgerufenen Eigenschaftsänderung einer Probe |
US20070034808A1 (en) * | 2005-04-27 | 2007-02-15 | Applied Research Associates, Inc. | Method and system for detection and identification of rapidly moving radioactive sources |
US7603138B2 (en) * | 2005-08-22 | 2009-10-13 | Toshiba American Research, Inc. | Environmental monitoring using mobile devices and network information server |
US7465924B1 (en) * | 2006-04-28 | 2008-12-16 | Uchicago Argonne Llc | Tracking of moving radioactive sources |
DE102007009639B4 (de) * | 2007-02-26 | 2010-04-08 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Verfahren zur Ermittlung einer hochgenauen Position von Fahrwegen und/oder Objekten |
WO2009006632A2 (en) * | 2007-07-05 | 2009-01-08 | Purdue Research Foundation | Nuclear detection via a system of widly distributed low cost detectors |
CN101470205B (zh) * | 2007-12-29 | 2011-07-27 | 同方威视技术股份有限公司 | 放射性物质方位探测设备和方法 |
GB0810638D0 (en) * | 2008-06-11 | 2008-07-16 | Rapiscan Security Products Inc | Photomultiplier and detection systems |
CN202522705U (zh) | 2011-12-23 | 2012-11-07 | 同方威视技术股份有限公司 | 放射性物质巡检定位设备 |
-
2011
- 2011-12-23 CN CN201110439199.4A patent/CN103176201B/zh active Active
-
2012
- 2012-12-20 GB GB1412789.8A patent/GB2512547B/en not_active Expired - Fee Related
- 2012-12-20 EP EP12858884.5A patent/EP2796898B1/en active Active
- 2012-12-20 WO PCT/CN2012/087012 patent/WO2013091552A1/zh active Application Filing
- 2012-12-20 US US14/367,824 patent/US9829601B2/en not_active Expired - Fee Related
-
2013
- 2013-11-26 HK HK13113202.2A patent/HK1185950A1/zh not_active IP Right Cessation
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH116875A (ja) * | 1997-06-16 | 1999-01-12 | Mitsubishi Electric Corp | シンチレーションファイバを用いた放射線検出方法 |
CN1318841C (zh) * | 2002-11-06 | 2007-05-30 | 美国科技工程公司 | 移动式x射线反向散射检查车 |
CN2690896Y (zh) * | 2003-09-03 | 2005-04-06 | 广州杰赛科技股份有限公司 | 数字化辐射巡测系统 |
CN1598553A (zh) * | 2003-09-18 | 2005-03-23 | 清华大学 | 一种对车辆放射性物质定位的监测方法及其装置 |
CN101069091A (zh) * | 2004-08-12 | 2007-11-07 | 纳沃特克医药有限公司 | 对象身体内的放射性源定位 |
WO2007131536A1 (en) * | 2006-05-11 | 2007-11-22 | Advanced Accelerator Applications | Method and device of detecting, locating and/or analyzing a radioactive source(s) in a material, e.g. a biological tissue |
Non-Patent Citations (1)
Title |
---|
张文 等: "基于GPS的γ射线车载巡测系统设计", 《核电子学与探测技术》, vol. 31, no. 4, 30 April 2011 (2011-04-30), pages 466 - 469 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104237959A (zh) * | 2014-08-21 | 2014-12-24 | 北京辛耕普华医疗科技有限公司 | 对放射性物质实时动态追踪定位的方法和设备 |
CN107977543A (zh) * | 2016-10-19 | 2018-05-01 | 中国辐射防护研究院 | 放射性物品运输弥散事故中集体剂量的计算方法 |
CN108445527A (zh) * | 2018-01-30 | 2018-08-24 | 奕瑞新材料科技(太仓)有限公司 | 一种射线源位置坐标的获取方法 |
CN109490939A (zh) * | 2018-11-23 | 2019-03-19 | 南京航空航天大学 | 一种应用于小型无人机的放射源定位方法 |
CN110850463A (zh) * | 2019-10-09 | 2020-02-28 | 岭澳核电有限公司 | 一种核电站放射性监测方法、装置、控制终端及监测系统 |
CN112147163A (zh) * | 2020-08-20 | 2020-12-29 | 中国原子能科学研究院 | 一种用于多点测量脉冲萃取柱钚浓度的串联探测装置 |
CN112357606A (zh) * | 2020-11-12 | 2021-02-12 | 重庆建安仪器有限责任公司 | 一种基于辐射监测的辐射标识器自动投放控制系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103176201B (zh) | 2015-03-25 |
US20150025841A1 (en) | 2015-01-22 |
GB2512547B (en) | 2019-04-24 |
EP2796898A1 (en) | 2014-10-29 |
US9829601B2 (en) | 2017-11-28 |
GB2512547A (en) | 2014-10-01 |
RU2014129589A (ru) | 2016-02-10 |
EP2796898B1 (en) | 2020-01-22 |
WO2013091552A1 (zh) | 2013-06-27 |
GB201412789D0 (en) | 2014-09-03 |
HK1185950A1 (zh) | 2014-02-28 |
EP2796898A4 (en) | 2015-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103176201B (zh) | 放射性物质巡检定位方法及设备 | |
CN108181635B (zh) | 一种用于输电线路交叉跨越分析的激光点云分类方法 | |
Gunatilaka et al. | On localisation of a radiological point source | |
CN105510952B (zh) | 飞行模式CdZnTe巡检系统和巡检方法 | |
Liu et al. | Sensor networks for the detection and tracking of radiation and other threats in cities | |
CN105203551A (zh) | 车载激光雷达隧道检测系统、基于隧道检测系统的自主定位方法及隧道灾害检测方法 | |
WO2008127442A9 (en) | High z material detection system and method | |
CN202522705U (zh) | 放射性物质巡检定位设备 | |
Sharma et al. | Three-dimensional localization of low activity gamma-ray sources in real-time scenarios | |
Bukartas et al. | A Bayesian method to localize lost gamma sources | |
CN104155677A (zh) | 用于pet探测器的符合响应线定位方法及装置 | |
Aucott et al. | Effects of background on gamma-ray detection for mobile spectroscopy and imaging systems | |
Mascarich et al. | Autonomous mapping and spectroscopic analysis of distributed radiation fields using aerial robots | |
KR20160081307A (ko) | 방사선원의 3차원 영상 처리 장치 및 방법 | |
US7863567B1 (en) | Multimodal radiation imager | |
CN107862747A (zh) | 一种γ放射源增强现实方法及系统 | |
WO2015057973A1 (en) | Systems and methods for high-z threat alarm resolution | |
Zabulonov et al. | Airborne gamma spectrometric survey in the Chernobyl exclusion zone based on oktokopter UAV type | |
CN113639970B (zh) | 一种评估卫星闪电成像仪探测地面定标能力的方法 | |
CN104485002B (zh) | 一种基于psd的交通流检测系统 | |
Tan et al. | A fast gamma-ray source localization method for mobile robots | |
US8878138B2 (en) | Multi-sensor neutron source location system | |
US20170146682A1 (en) | Systems and methods for identifying anomalous radiation measurements | |
Wang et al. | Adaptive noise filtering for single photon Lidar observations | |
CN112285758A (zh) | 一种核辐射场探头、探测器和探测系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1185950 Country of ref document: HK |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: GR Ref document number: 1185950 Country of ref document: HK |