Advanced Search
Volume 41 Issue 1
Mar.  2024
Turn off MathJax
Article Contents

Long WEI, Lingqin KONG, Lei SHUAI, Xiuzuo LIANG, Xuanhou HU, Yiwen ZHANG, Xiaoming WANG. Radiation Imaging Methods and Intelligent Technology in Nuclear Safety[J]. Nuclear Physics Review, 2024, 41(1): 94-108. doi: 10.11804/NuclPhysRev.41.2023CNPC82
Citation: Long WEI, Lingqin KONG, Lei SHUAI, Xiuzuo LIANG, Xuanhou HU, Yiwen ZHANG, Xiaoming WANG. Radiation Imaging Methods and Intelligent Technology in Nuclear Safety[J]. Nuclear Physics Review, 2024, 41(1): 94-108. doi: 10.11804/NuclPhysRev.41.2023CNPC82

Radiation Imaging Methods and Intelligent Technology in Nuclear Safety

doi: 10.11804/NuclPhysRev.41.2023CNPC82
Funds:  National Natural Science Foundation of China (12375307, 12005234, 12105116)
More Information
  • Corresponding author: E-mail: shuail@ihep.ac.cn
  • Received Date: 2024-02-10
  • Rev Recd Date: 2024-03-01
  • Available Online: 2024-03-15
  • Publish Date: 2024-03-20
  • Radiation imaging is an important aspect of nuclear radiation measurement. In the field of radiation safety, typical radiation imaging methods mainly include pinhole imaging, (spatial) coded aperture imaging, time-coded imaging, and Compton scattering imaging, etc, each with its own advantages, disadvantages, and applicable scenarios. Based on the accumulation of existing imaging technology methods and the experience in the development of advanced imaging equipment, this article explores three feasible schemes for the three-dimensional positioning of radiation hotspots, the application of artificial intelligence in the dynamic imaging of radiation hotspots, nuclear identification, and the intelligentization of radiation imaging platforms. It addresses the actual needs of radiation detection (imaging) that is “faster, more accurate, and covers a larger range.” In the research on the three-dimensional positioning of radiation hotspots, a four-eye stereo imaging system has achieved a distance measurement accuracy of less than 20% for radioactive sources at a distance of 100 meters in long-distance imaging scenarios. In indoor scenes, three-dimensional positioning of four radioactive sources has been achieved through three-dimensional scene modeling and radiation field construction. In the application of artificial intelligence algorithms, real-time positioning of moving radiation hotspots has been realized, capable of dynamically tracking and positioning a 14 mCi 137Cs radioactive source moving at a speed of 36 km/h at a distance of 25 meters. In the research on neural network nuclear identification algorithms, a feature-enhanced convolutional neural network based on prior energy information has been proposed, which can quickly identify nuclear substances in environments with low counts and mixed nuclear species. In the development of intelligent radiation imaging platforms, an accelerator tunnel environmental inspection robot capable of detecting radiation hotspots has been developed, and a vehicle-mounted radiation emergency monitoring platform has been created. Equipped with multiple detection terminals such as unmanned aerial vehicle (UAV) mounted, vehicle-mounted, and portable gamma cameras, it can quickly mobilize measurements in urban areas during radiation emergencies. The various intelligent radiation imaging technologies presented in this paper enable people to obtain the three-dimensional distribution and nuclear information of radioactive sources more quickly and accurately, providing advanced and efficient radiation detection technical means for radiation environmental assessment, radiation protection, nuclear emergency drills, and decision-making.
  • [1] ANDERSON D N, STROMSWOLD D C, WUNSCHEL S C, et al. Technometrics, 2012, 48(2): 252. doi:  10.1198/004017005000000463
    [2] IAEA Releases Annual Data on Illicit Trafficking of Nuclear and other Radioactive Material | IAEA[EB/OL]. [2024-02-25] https://www.iaea.org/newscenter/pressreleases/iaea-releases-annual-data-on-illicit-trafficking-of-nuclear-and-other-radioactive-material.
    [3] 国家核安全局. 国家核安全局2022年年报[R]. 北京: 国家核安全局, 2023.

    National Nuclear Safety Administration. Annual Report 2022[R]. Beijing: National Nuclear Safety Administration, 2023. (in Chinese)
    [4] PRADEEP KUMAR K A, SHANMUGHA SUNDARAM G A, SHARMA B K, et al. Nuclear Engineering and Technology, 2020, 52(10): 2151. doi:  10.1016/j.net.2020.03.014
    [5] 陈伯显. 核辐射物理及探测学[M]. 2版. 哈尔滨: 哈尔滨工程大学出版社, 2021.

    CHEN Boxian. Nuclear Radiation Physics and Detection[M]. 2nd ed. Harbin: Harbin Engineering University Press, 2021. (in Chinese)
    [6] KIM J, BACK H, JOO K. Nuclear Engineering and Technology, 2020, 52: 456. doi:  10.1016/j.net.2019.07.027
    [7] KONG L, SHUAI L, LIANG X, et al. Review of Scientific Instruments, 2022, 93(11): 113103. doi:  10.1063/5.0097257
    [8] WEINFURTHER K, MATTINGLY J, BRUBAKER E, et al. Nucl Instr and Meth A, 2018, 883: 115. doi:  10.1016/j.nima.2017.11.025
    [9] SHAH N. Adaptive Imaging with a Cylindrical, Time-Encoded Imaging System[D]. Ann Arbor: University of Michigan, 2020.
    [10] MORTIMER R K, ANGER H O, TOBIAS C A. The Gamma-ray Pinhole Camera With Image Amplifier: UCRL-2524[R]. Berkeley: University of California, 1954.
    [11] AMGAROU K, HERRANZ M. Nuclear Engineering and Technology, 2021, 53(11): 3491. doi:  10.1016/j.net.2021.05.031
    [12] CIESLAK M. Development of Scintillator Based Coded-aperture Neutron Imager for Nuclear Decommissioning[D]. Lancaster: Lancaster University, 2019.
    [13] SUEOKA K, KATAOKA J, TAKABE M, et al. Nucl Instr and Meth A, 2018, 912: 115. doi:  10.1016/j.nima.2017.10.082
    [14] POTAPOV V N, IVANOV O P, SEMIN I A. Method to Improve Angular Resolution of Portable Gamma Camera with Pinhole Collimator[C]//2016 IEEE Nuclear Science Symposium, Medical Imaging Conference and Room-Temperature Semiconductor Detector Workshop (NSS/MIC/RTSD). 2016.
    [15] DICKE R H. Astrophysical Journal, 1968, 153: L101. doi:  10.1086/180230
    [16] ABLES J G. Publications of the Astronomical Society of Australia, 1968, 1(4): 172. doi:  10.1017/S1323358000011292
    [17] BRADY D J. Optical Imaging and Spectroscopy[M]. Washington D C: Wiley, 2009.
    [18] 帅磊. 编码孔径γ相机数据获取及图像重建系统研究[D]. 北京: 中国科学院大学, 2010.

    SHUAI L. Research on Data Acquisition and Image Reconstruction System of Coded Aperture Gamma Camera[D]. Beijing: University of Chinese Academy of Sciences, 2010. (in Chinese)
    [19] 孙世峰. 基于编码孔径的射线成像方法研究[D]. 北京: 中国科学院大学, 2016.

    SUN S F. Research on Radiographic Imaging Method Based on Coded Aperture [D]. Beijing: University of Chinese Academy of Sciences, 2016. (in Chinese)
    [20] BRENNAN J, BRUBAKER E, GERLING M, et al. Nucl Instr and Meth A, 2015, 802: 76. doi:  10.1016/j.nima.2015.08.076
    [21] ODA M. Applied Optics, 1965, 4(1): 143. doi:  10.1364/AO.4.000143
    [22] KIM H S, CHOI H Y, LEE G, et al. Journal of Radiological Protection, 2018, 38(1): 299. doi:  10.1088/1361-6498/aaa3c8
    [23] KORAL K F, ROGERS W L, KNOLL G F. Journal of Nuclear Medicine, 1975, 16(5): 402.
    [24] MARLEAU P, BRENNAN J, BRUBAKER E, et al. Results from the Coded Aperture Neutron Imaging System[C]//IEEE Nuclear Science Symposuim & Medical Imaging Conference. 2010: 1640.
    [25] BROWN S. Time-Encoded Thermal Neutron Imaging using Large-Volume Pixelated CdZnTe Detectors[D]. Ann Arbor: The University of Michigan, 2017.
    [26] KIM H S, YE S J, SHIN Y, et al. Journal of the Korean Physical Society, 2016, 69(11): 1644. doi:  10.3938/jkps.69.1644
    [27] MARLEAU P, BRUBAKER E. An Implementation of Zero Knowledge Confirmation Using a Two-dimensional Time-Encoded Imaging System[R]. Livermore: Sandia National Lab, 2016.
    [28] LIANG X Z, PANG X Y, CAO D Q, et al. Nucl Instr and Meth A, 2020, 951: 162964. doi:  10.1016/j.nima.2019.162964
    [29] BAO V L, KIM G. A Journal of Instrumentation, 2017, 12(03): C03013. doi:  10.1088/1748-0221/12/03/C03013
    [30] NA M, LEE S, KIM H S, et al. Journal of Instrumentation, 2017, 12(12): C12050. doi:  10.1088/1748-0221/12/12/C12050
    [31] KIM H S, NA M, OK J G, et al. IEEE Transactions on Nuclear Science, 2019, 66(9): 2114. doi:  10.1109/TNS.2019.2933442
    [32] OLESEN R J. Low-Information Radiation Imaging using Rotating Scatter Mask Systems and Neural Network Algorithms[D]. Montgomery: Air University, 2020.
    [33] OLESEN R J, HOLLAND D E, BRUBAKER E M, et al. Radiation Measurements, 2020, 137: 106441. doi:  10.1016/j.radmeas.2020.106441
    [34] SCHÖNFELDER V, HIRNER A, SCHNEIDER K. Nucl Instr and Meth, 1973, 107(2): 385. doi:  10.1016/0029-554x(73)90257-7
    [35] TODD R W, NIGHTINGALE J M, EVERETT D B. Nature, 1974, 251(5471): 132. doi:  10.1038/251132a0
    [36] 梁秀佐. 基于编码孔径与康普顿散射方法的关键成像技术研究[D]. 北京: 中国科学院高能物理研究所, 2022.

    LIANG X Z. Research on Key Imaging Techniques based on Coded Aperture and Compton Scattering Method[D]. Beijing: Institute of High Energy Physics, Chinese Academy of Sciences, 2022. (in Chinese)
    [37] 武传鹏, 李亮. 核技术, 2021, 44(5): 50403. doi:  10.11889/J.0253-3219.2021.hjs.44.050403

    WU C P, LI L. Nuclear Technology, 2021, 44(5): 50403. (in Chinese) doi:  10.11889/J.0253-3219.2021.hjs.44.050403
    [38] POITRASSON-RIVIÈRE A, MAESTAS B A, HAMEL M C, et al. Progress in Nuclear Energy, 2015, 81: 127. doi:  10.1016/j.pnucene.2015.01.009
    [39] ZHANG J, WEI L, SHUAI L, et al. IEEE Transactions on Nuclear Science, 2019, 67: 1691. doi:  10.1109/TNS.2019.2956180
    [40] FRANDES M, ZOGLAUER A, MAXIM V, et al. IEEE Transactions on Nuclear Science, 2010, 57: 144. doi:  10.1109/TNS.2009.2031679
    [41] YE Q, FAN P, WANG R, et al. Nucl Instr and Meth A, 2019, 937: 31. doi:  10.1016/j.nima.2019.05.022
    [42] HU Y, FAN P, LYU Z, et al. Nucl Instr and Meth A, 2022, 1023: 165971. doi:  10.1016/j.nima.2021.165971
    [43] PARADISO V, AMGAROU K, DE LANAUTE N B, et al. Nucl Instr and Meth A, 2018, 910: 194. doi:  10.1016/j.nima.2018.09.081
    [44] RATHNAYAKA P, BAEK S H, PARK S Y. Introducing a New Radiation Detection Device Calibration Method and Estimating 3D Distance to Radiation Sources[C]//KIM K J. Information Science and Applications. Berlin, Heidelberg: Springer, 2015: 335.
    [45] PARADISO V, CRIVELLARO A, AMGAROU K, et al. Nucl Instr and Meth A, 2018, 886: 125. doi:  10.1016/j.nima.2017.12.065
    [46] SATO Y. IApplied Radiation and Isotopes, 2023, 195: 110739. doi:  10.1016/j.apradiso.2023.110739
    [47] MIHAILESCU L, VETTER K, CHIVERS D. IEEE Transactions on Nuclear Science, 2009, 56(2): 479. doi:  10.1109/TNS.2009.2015304
    [48] CHIN J C, RAO N S, YAO D K Y, et al. ACM Transactions on Sensor Networks, 2010, 7(3): 50. doi:  10.1145/1807048.1807050
    [49] SUN S F, ZHANG Z M, SHUAI L, et al. Applied Radiation and Isotopes, 2016, 107: 177. doi:  10.1016/j.apradiso.2015.10.021
    [50] LIU Y T, SHUAI L, LI D, et al. AIP Advances, 2022, 12(8): 085317. doi:  10.1063/5.0096865
    [51] YU Y, LIU S, ZHANG Z, et al. IEEE Transactions on Nuclear Science, 2022, 69(8): 1931. doi:  10.1109/TNS.2022.3186433
    [52] TAKEUCHI K, KATAOKA J, NISHIYAMA T, et al. Nucl Instr and Meth A, 2014, 765: 187. doi:  10.1016/j.nima.2014.04.039
    [53] SATO Y, TERASAKA Y, OZAWA S, et al. Journal of Instrumentation, 2017, 12(11): C11007. doi:  10.1088/1748-0221/12/11/C11007
    [54] 邹艺. 基于编码孔径成像的移动放射源实时定位追踪算法研究[D]. 北京: 中国科学院大学, 2021.

    ZOU Y. Research on Real-time Location and Tracking Algorithm of Mobile Radioactive Source Based on Coded Aperture Imaging[D]. Beijing: University of Chinese Academy of Sciences, 2021. (in Chinese)
    [55] ZOU Y, LIU S Q, SUN X L, et al. Radiation Detection Technology and Methods, 2021, 5(2): 228. doi:  10.1007/s41605-020-00232-7
    [56] REDMON J, FARHADI A. YOLOv3: An Incremental Improvement. arXiv: 1804.02767, 2018. DOI: 10.48550/arXiv.1804.02767.
    [57] WOJKE N, BEWLEY A, PAULUS D. Simple online and realtime tracking with a deep association metric[C]//2017 IEEE International Conference on Image Processing (ICIP). 2017: 3645.
    [58] 张江梅, 任俊松, 李培培, 等. 核电子学与探测技术, 2016, 36(8): 856. doi:  10.3969/j.issn.0258-0934.2016.08.019

    ZHANG J M, REN J S, LI P P et al. Nuclear Electronics and Detection Technology, 2016, 36(8): 856. (in Chinese) doi:  10.3969/j.issn.0258-0934.2016.08.019
    [59] LI C, LIU S, WANG C, et al. Applied Radiation and Isotopes, 2022, 185: 110219. doi:  10.1016/j.apradiso.2022.110219
    [60] 李春苗. 基于计数型全方向探测的多放射源三维定位算法研究[D]. 北京: 中国科学院大学, 2022.

    LI C M. Research on Three-dimensional Location Algorithm of Multiple Radioactive Sources Based on count-type Omnidirectional Detection[D]. Beijing: University of Chinese Academy of Sciences, 2022. (in Chinese)
    [61] ZHANG R, GONG P, TANG X, et al. Nucl Instr and Meth A, 2019, 934: 41. doi:  10.1016/J.NIMA.2019.04.055
    [62] LIU Q, CHENG Y, TUO X, et al. Applied Radiation and Isotopes, 2021, 170: 109637. doi:  10.1016/j.apradiso.2021.109637
    [63] WEN H, WANG Y, XIUZUO L, et al. Review of Scientific Instruments, 2021, 92: 013106. doi:  10.1063/5.0035991
    [64] VETTER K, BARNOWSKI R, CATES J W, et al. Sensors, 2019, 19(11): 2541. doi:  10.3390/s19112541
    [65] CHIERICI A, MALIZIA A, DI GIOVANNI D, et al. The European Physical Journal Plus, 2021, 136(3): 314. doi:  10.1140/epjp/s13360-021-01276-4
    [66] WHITE S R, MEGSON-SMITH D A, ZHANG K, et al. Frontiers in Robotics and AI, 2020, 7. DOI: 10.3389/frobt.2020.499056.
    [67] CONNOR D T, WOOD K, MARTIN P G, et al. Frontiers in Robotics and AI, 2020, 6: 00149.
    [68] MARQUES L, VALE A, VAZ P. Sensors, 2021, 21(4): 1051. doi:  10.3390/s21041051
  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(24)  / Tables(2)

Article Metrics

Article views(73) PDF downloads(10) Cited by()

Proportional views

Radiation Imaging Methods and Intelligent Technology in Nuclear Safety

doi: 10.11804/NuclPhysRev.41.2023CNPC82
Funds:  National Natural Science Foundation of China (12375307, 12005234, 12105116)

Abstract: Radiation imaging is an important aspect of nuclear radiation measurement. In the field of radiation safety, typical radiation imaging methods mainly include pinhole imaging, (spatial) coded aperture imaging, time-coded imaging, and Compton scattering imaging, etc, each with its own advantages, disadvantages, and applicable scenarios. Based on the accumulation of existing imaging technology methods and the experience in the development of advanced imaging equipment, this article explores three feasible schemes for the three-dimensional positioning of radiation hotspots, the application of artificial intelligence in the dynamic imaging of radiation hotspots, nuclear identification, and the intelligentization of radiation imaging platforms. It addresses the actual needs of radiation detection (imaging) that is “faster, more accurate, and covers a larger range.” In the research on the three-dimensional positioning of radiation hotspots, a four-eye stereo imaging system has achieved a distance measurement accuracy of less than 20% for radioactive sources at a distance of 100 meters in long-distance imaging scenarios. In indoor scenes, three-dimensional positioning of four radioactive sources has been achieved through three-dimensional scene modeling and radiation field construction. In the application of artificial intelligence algorithms, real-time positioning of moving radiation hotspots has been realized, capable of dynamically tracking and positioning a 14 mCi 137Cs radioactive source moving at a speed of 36 km/h at a distance of 25 meters. In the research on neural network nuclear identification algorithms, a feature-enhanced convolutional neural network based on prior energy information has been proposed, which can quickly identify nuclear substances in environments with low counts and mixed nuclear species. In the development of intelligent radiation imaging platforms, an accelerator tunnel environmental inspection robot capable of detecting radiation hotspots has been developed, and a vehicle-mounted radiation emergency monitoring platform has been created. Equipped with multiple detection terminals such as unmanned aerial vehicle (UAV) mounted, vehicle-mounted, and portable gamma cameras, it can quickly mobilize measurements in urban areas during radiation emergencies. The various intelligent radiation imaging technologies presented in this paper enable people to obtain the three-dimensional distribution and nuclear information of radioactive sources more quickly and accurately, providing advanced and efficient radiation detection technical means for radiation environmental assessment, radiation protection, nuclear emergency drills, and decision-making.

Long WEI, Lingqin KONG, Lei SHUAI, Xiuzuo LIANG, Xuanhou HU, Yiwen ZHANG, Xiaoming WANG. Radiation Imaging Methods and Intelligent Technology in Nuclear Safety[J]. Nuclear Physics Review, 2024, 41(1): 94-108. doi: 10.11804/NuclPhysRev.41.2023CNPC82
Citation: Long WEI, Lingqin KONG, Lei SHUAI, Xiuzuo LIANG, Xuanhou HU, Yiwen ZHANG, Xiaoming WANG. Radiation Imaging Methods and Intelligent Technology in Nuclear Safety[J]. Nuclear Physics Review, 2024, 41(1): 94-108. doi: 10.11804/NuclPhysRev.41.2023CNPC82
    • 放射性物质探测,特别是核材料走私拦截和威胁探测,对国家安全至关重要[1]。根据国际原子能机构(International Atomic Energy Association, IAEA)的事件和贩运数据库(Incident and Trafficking Database, ITDB)报告,仅在2022年内,就发生了146起涉及核材料和其他放射性材料的非法或未经授权的活动事件。此外,ITDB表明,自1993年以来,在运输核材料或放射性材料期间发生的盗窃几乎占所有报告案件的52%,而这一比例在过去十年中则涨到了62%,凸显了加强安全运输措施的重要性[2]。我国是放射性同位素应用大国,据国家核安全局2022年年报,截至2022年底,国内目前在用放射源达164 028枚,其中Ⅰ类放射源为15 746枚[3],放射源的管理、运输、储存和使用是一项庞大且复杂的工作,而辐射监测的技术伴随着放射源从生产到最后回收处理的整个生命周期。

      辐射测量可以分为放射性污染控制(Radioactive Pollution Control)和紧急辐射监测(Emergency Radiation Monitoring, ERM)两类[4],如图1所示,对于核电厂等核设施的日常持续的监测属于放射性污染控制,而对于核事故造成的辐射污染评估等则属于ERM的范围。辐射测量按任务还可以细分为活度测量、辐射场量的测量、辐射能量(或能谱测量)、辐射剂量的测量、时间的测量以及辐射的空间分布测量与辐射成像等,辐射成像经过近几十年的迅猛发展,目前已成为独立学科[5]。核辐射测量仪器按技术特点可以分为报警仪、剂量仪、能谱仪、成像仪等。

      其中报警仪,如门式辐射监测系统、个人辐射剂量仪等,通常只能简单判断辐射剂量是否超标。剂量仪则能够以数值方式实时显示辐射剂量,但不能区分辐射来源与核素种类;能谱仪则能够显示辐射能谱,判断核素种类与剂量。图2(a)展示了Kim等[6]应用CsI(Tl)晶体配合SiPM(silicon-photomultiplier)开发的辐射探测背包所得到的多种放射源混合能谱,但其不能判断核辐射的来源。以上几种核辐射探测仪器,都需要通过近距离接触的方式进行探测,属于典型的“接近探测器(proximity detectors)”。辐射成像技术是随着探测器与成像算法的发展而产生的新技术,该技术采用“拍照”的方式,在短时间内远距离给出视野范围内放射性物质分布情况的二维图像。图2(b)展示了应用康普顿成像技术配合光学图像对放射源进行二维成像的结果[7]。图中直观指示了热点的位置,通过进一步的解析放射性物质能谱与图像重建算法优化,可以确定放射源的种类,实现核素识别的功能。

      辐射成像技术不仅可用于核事故应急响应时的“未知”放射源搜寻和排查,还可应用于核电站、核废料处理厂等的现场远距离放射性物质检测和分布评估,大大减少工作人员近距离接触放射性物质可能受到的各类辐射。如核电站例行检修时,将不同厂房按辐射剂量水平进行分类,对房间内的放射性热点进行定位以及有针对性的屏蔽,并对屏蔽效果进行评价。在做核设施的去污及退役的测量时,辐射成像技术能更加直观地对核设施退役时拆除的部件及核废料容器进行放射性残留测量,确定其辐射水平及放射性物质的分布情况,以最大限度地减少辐射暴露。

    • 在核安全领域中,对于所有需要γ射线(或中子)成像的核不扩散应用场景,一般有两类成像系统:一类是基于遮挡的成像系统(occlusion-based imaging systems),此类成像系统应用了射线沿直线传播的特点,通过准直器的遮挡实现对射线束流进行空间或时间调制成像,如针孔成像系统、编码孔径成像系统、时间编码成像系统等;另一类则是基于散射的成像系统(scatter-based imaging systems),也称为散射相机(scatter cameras),依赖于散射物理来定位辐射源。对于γ射线,散射相机依赖于康普顿散射原理;对于快中子,则依赖于弹性散射[89]

    • 针孔成像于20世纪中期运用于辐射成像领域,属于早期的核辐射成像技术[1011],也称为“暗室(dark chamber)”,其基本成像原理与“小孔成像”相同,如图3(a)所示,目标放射源发出的射线通过针孔(准直器),在像平面(位置灵敏探测器)上形成与物体点对点对应的倒立像。为了获得γ放射源图像,需要在探测器上包裹一层重金属材料屏蔽层,放射源产生的γ射线只能通过针孔到达探测器,因此针孔γ相机的灵敏度(探测效率)较低。针孔的尺寸还直接影响了针孔γ相机的角分辨率和信噪比(Signal-to-Noise Ratio, SNR),针孔直径越小,角度分辨率越好,信噪比越差[11],理想的单针孔相机需要一个无限小的针孔来获得最高的角度分辨率[12],因此,必须对这两个性能参数进行折中[11]。由于针孔γ相机一般具有厚实的屏蔽层,使得它可以适用于高剂量率环境。如在福岛第一核电站反应堆安全壳内,其辐射剂量率高达530 Sv/h,其他类型的γ相机可能会因信号堆积而无法工作[11, 13]图3(b)展示了Koki Sueoka等研发的针孔γ相机在680 Sv/h 剂量率环境下对60Co的成像结果。

    • 针孔成像是编码孔径成像的前身,多孔径编码成像方法于1968年被Dicke和Ables提出[1516],开发编码孔径成像的主要动机是在不降低系统分辨率的情况下增加射线的通量[17]。编码孔径成像的技术原理如图4所示,在视野张角内的放射源以球面光的方式发出γ射线或者中子,射线束流经过编码码板的空间调制,在探测器平面的晶体阵列上得到编码图案的投影(即投影矩阵),理想状况下,编码投影为射线源分布与编码函数(编码图案)的卷积。编码的本质其实是数学变换在成像中的物理表现形式,投影矩阵与解码矩阵经过直接反卷积(direct deconvolution)或迭代算法进而得到重建图像[1819]。编码孔径成像技术因其成像灵敏度高,图像信噪比高,近年来在核安全领域的诸多涉核场所中已经成为一种应用最广泛的辐射成像技术。

      在编码孔径成像中,编码图案通常被选择为最小化相关性(minimize correlations)[20],按照编码图案和系统点扩展函数的不同,编码图样可分为均匀冗余阵列编码(URA、MURA、PNP、M-P)和随机编码(RANDOM)这两种编码方式[18]。均匀冗余阵列编码需要特定排列的方形编码矩阵,编码维度数值必须为素数,它的系统点扩展函数为理想的δ函数,能够通过互关联解码实现快速图像重建。随机编码的编码维度数值可以为任意值,适应于对编码板的大小有特殊要求且不方便设计维度数值为特定素数的应用情况,因其系统点扩展函数有固有的本底噪声涨落,采用的图像重建算法多为极大似然期望最大化(Maximum Likelihood Estimation Maximization, MLEM)的迭代重建算法。

    • 时间编码孔径成像技术的概念最早于1965年由Minoru Oda提出并用于空间科学[2122],1974年应用于核医学成像[23]。但直到2010年,有人通过构建圆柱形时间编码成像(Cylindrical, Time-Encoded Imaging, c-TEI)系统对快中子进行了初步成像[9, 2425],时间编码成像技术才引起辐射安全研究者们的重视。在发射成像(emission imaging)领域,时间编码是指根据先验测量或已知的时变模式,调制从某一方向入射的射线的通量,在探测器上将得到随时间波动的计数率曲线,该计数率曲线与已知的时变模式相关,这种关系满足了辐射成像的基本要求:将测量与来自特定方向或空间位置的粒子强度联系起来[25]。典型的时间编码成像方式为旋转调制(rotational modulation),如图5所示,成像系统由一个探测器和一组栅条准直器构成,当两个准直器旋转时,栅条构成的“通孔”大小会发生变化,引起探测器内计数(率)的涨落(即时间调制),从而反推出放射源的分布信息。搭配具有脉冲形状辨别(pulse shape discrimination, PSD)功能的闪烁体探测器,可以构建γ射线/中子双粒子成像系统[26]

      基于时间编码的成像系统具有诸多优点,首先是成像系统只需要较低的通道计数,甚至是单通道计数,可以简化系统集成和校准过程,减少系统的不确定性,降低成本,提高系统稳定性。其次,重建图像的角分辨率主要取决于准直器的设计,角分辨率和能量分辨都可以独立优化,从而可以提高γ射线的探测效率和晶体选择的灵活性。此外,相对于空间编码孔径成像系统,时间编码成像系统的设计更加灵活多样,目前已有多种时间编码的准直器设计方案,如图6所示,包括具备360°视野的圆柱形旋转准直器(cylindrical, rotating collimator)[20, 2728]、半球旋转调制准直器 (hemispherical Rotational Modulation Collimator, H-RMC)[2931]、旋转散射码板(rotating scatter mask, RSM)[3233]等。

    • 康普顿成像技术利用康普顿散射效应来定位放射源位置,是辐射成像领域中较为复杂的成像技术,出现时间也最晚,于1973年由Schönfelder等[34]提出,用于伽马射线天文学的观测;1974年,Todd 等[35]提出了现代意义上康普顿相机,明确了散射-吸收的物理过程及反投影重建算法。康普顿成像的基本原理如图7所示,入射能量为E0的γ光子在散射点中发生康普顿散射,沉积能量为E1;散射光子在吸收点被吸收,沉积能量为E2,对于用于成像的“一次散射事例”,E1+E2即为入射光子能量E0mec2为电子的静止能量(其中me为电子的静止质量,c为真空光速),以散射点与吸收点连线的延长线为轴线,通过式(1)计算出散射角θγ

      由此得到“康普顿圆锥”,圆锥与成像空间相交即可得到一个圆环,通过多个事例的交叠,即可在成像空间中构建出放射源的位置,图中θγ为圆环的展宽,即康普顿散射角的不确定度,其主要的影响因素包括探测器能量分辨率、探测器空间分辨率、多普勒展宽效应[3637]。在成像系统中,一般需要预置多种核素的全能峰范围,成像时根据所测得的能谱情况,选择最大概率能量范围内的事例作为对应核素的成像事例,即在康普顿成像的重建算法中,入射γ射线的能量需要是已知的。

      典型康普顿成像系统的探测器为双层结构,如图8(a)所示,光子的散射、吸收过程分别发生在两层探测器内。在此基础上,为了提高成像系统的探测效率、扩大成像视野,研究人员开发了多层探测器结构、三维位置灵敏探测器结构的康普顿成像系统,如图8(b)所示,为通过激光内雕和双端读出技术构建的三维位置灵敏探测器,由此构建的康普顿相机具备4π全空间成像的功能[39]。此外,还有一种追踪反冲电子径迹(electron-tracking)的康普顿相机[40],如图8(c)所示,通过将“康普顿圆锥”的限制为“康普顿圆弧”,提高放射源定位精度,但技术复杂性也大大增加。康普顿成像则具有视野范围广的突出优点,但受限于系统的能量分辨与空间分辨,导致重建图像的角分辨率较差。

    • 表1对不同类型辐射成像系统的优缺点进行了总结[9, 11]

      类型优点缺点
      针孔成像角度分辨率高(1.9°~6.7°)设备重
      成像能量范围宽(241Am到60Co)灵敏度低
      具有良好的线性剂量率视场小(30°或者50°)
      信噪比高
      抗噪能力强
      编码孔
      径成像
      设备紧凑视场较小(45°~50°)
      灵敏度高能量分辨率一般
      角度分辨率高(2.5°~6°)
      剂量率线性可达10 Sv/h
      时间编
      码成像
      结构简单灵敏度低
      稳定性高
      角分辨率和能量分辨可以独立优化
      康普顿成像设备紧凑灵敏度低
      视场可达4π角度分辨率差(10°~30°)
      能量分辨率高

      上述成像技术手段在方法和原理上奠定了辐射成像的基础。除了上述几种典型的核辐射成像方法,近年来,清华大学提出了一种新型的自准直辐射成像技术,该技术不需要重金属准直器,通过优化探测器单元的空间排布,使得投影(光子事件密度分布)对入射伽马射线的方向敏感,目前已实现4π视场内对伽马射线进行成像[4142],具有极大的发展前景。然而,现有辐射成像应用场景的复杂性对成像方法提出了更多的需求,比如测量放射源的距离信息(三维定位)、辐射热点的动态跟踪成像、源项的核素分析以及各种智能化平台的应用等,这些均需要对现有技术方法做改进和提升。

    • 基于现有成像技术方法上的积累以及丰富的先进成像设备研制经验,作者及其研究团队在辐射成像方面做了大量的智能化应用实践,本小节在内容上分为三部分展开具体描述:1) 成像算法的拓展,如辐射热点三维定位的实现方法;2) 神经网络等人工智能的应用,包括移动放射源的动态成像、核素识别等;3) 辐射探测系统与智能化平台的联合应用。

    • 对辐射热点的三维定位是近几年辐射成像领域一个热门的话题[4345],由于辐射成像系统获得的图像通常为表示放射性相对强度的二维图像,但二维辐射热点图缺乏深度信息(即成像系统与热点之间的距离信息),无法确定放射源的具体位置[46]。而在核事故中,通过确定辐射热点的三维位置,估算辐射热点的活度,对于进一步减少个人辐射剂量、判断事故后续的影响、调整事故处理方案有着重要价值。放射源的三维定位有多种方法,如搭配深度相机法[45]、两点联合定位法[43]、多点定位探测法[4748]法等。图9展示了Paradiso等[45]通过γ相机搭配深度相机对辐射热点进行测距的实验结果,但此方法仅在放射源无遮挡情况下有效,尚未完全解决辐射热点的三维定位问题,以下就三种通过对辐射热点进行直接测量来获取深度信息的方法进行介绍。

      第一种方法基于编码投影放大效应[49]。如图10所示,在编码孔径成像中,辐射场中的放射性点源产生的射线,经过编码孔径准直器调制后,在探测器平面上得到的码板图样投影大小和该点源与探测器平面的距离有关,即不同的物距对应不同的放大系数。当选取的投影区域较大或较小时,重建图像的信噪比将会远差于选取合适大小投影区域时的重建图像,如图11所示,中间图像即为选取合适投影区域时的重建图像,该图像的伪影明显少于其他两幅。因此,可以通过对编码投影图像进行高分辨率精细采样,然后选取不同大小的投影区域进行图像重建;再对得到的重建图像序列进行图像质量评价来识别其放大系数,从而确定射线源距离成像系统的距离。图12为Sun等[49]分别用Rank 11和Rank 19的 MURA编码板,配合75×75 NaI(Tl)晶体阵列探测器,对99mTc点源(2.61 mCi)进行深度提取的测试结果。实验结果表明,在0.9~8.0 m范围内,通过投影放大效应能正确提取点源的深度信息,并能保持良好的线性度。

      辐射热点的三维定位的第二种实现方式则是构建多目相机,利用立体几何非线性逼近的最小面积来估计射线源到探测器的最大概率距离[5051]。基本原理如图13(a)所示,假设存在4个γ相机模块分别对放射源经行解码成像,4个成像模块的理论成像顶点分别为O1O2O3O4,连接顶点与辐射热点可以得到4条探测线。理论上,这4条线应在点源处相交;然而,在大多数情况下,由于重建投影的离散性,它们不会完全相交。此时,利用$ Z=Z^{\prime} $的垂直平面与四条探测线相交可以得到以Q1Q2Q3Q4为顶点的四边形,如图13(b)所示,横坐标为不同深度,纵坐标为围成的四边形面积,此四边形的面积A(Z') 反映了点源到探测器距离估计值Z'的接近程度。显然,A(Z')min(即四边形Q1Q2Q3Q4面积的最小值)处可以认为是点源到探测器的最大概率距离。

      图14(左)为Liu等[50]开发的四目立体成像系统实验装置,图14(右)为放射源三维定位能力的测试场景,实验中,1 mCi的137Cs点源到探测器的实际距离为10.0 m。成像系统分别在无遮挡、白板遮挡、铅玻璃遮挡三种情况下得到的估计距离为10.21, 10.03 和10.37 m。实验结果表明,即使在一定量辐射屏蔽下,利用多目相机及“最小面积算法”的测距能力依然能保持较好的测距精度。同时,该三维定位方法在多点源同时存在的复杂情况下,在更远距离(~100 m)下依旧具备较好的测距能力(例如,对30 m内的辐射热点测距精度小于12%,100 m测距精度小于20%)[50]。上述用到的多目立体成像系统是基于编码孔径成像实现的,同样,利用多个康普顿相机从多个角度测量也可以估算出放射源的距离信息[46, 5253]

      另一种放射源三维定位方法则是通过应用远程控制设备,携带剂量仪等辐射探测装置在辐射场中进行非成像的多点巡测,通过构建所探测空间的辐射场实现热点的三维定位。如图15所示,以无人车作为载体,搭载激光雷达、辐射探测器在辐射场中运动,激光雷达负责三维实物场景的点云采样和空间几何构建,获取辐射场实物的三维模型,根据探测场景的扫描结果,定制定点探测方案,规划出无人车的扫描路径,最后利用无人车辐射场景进行巡测和辐射场重构确定出热点位置。

      场景扫描及热点位置重建结果如图16所示,测试空间的范围为20.0 m×20.0 m×1.2 m,空间内存在4枚放射源,巡测过程中的探测点数量为 21 个。表2展示了全方向探测器放射源三维定位系统对四个点源辐射剂量场重建结果,可以看出,在XY平面的重建位置偏差均小于1.2 m。

      编号 放射源种类 放射源活度 实际位置/m 重构位置/m 偏差/m
      1 Cs-137 1 mCi (9.3, 11.8, 0) (9.6, 11.6, 0) 0.36
      2 Cs-137 1 mCi (16.5, −1.2, 0) (17.0, −0.2, 0) 1.12
      3 Tc-99m 1.5 mCi (8.2, 0.6, 0) (7.6, 0.7, 0) 0.61
      4 Tc-99m 1.5 mCi (12.8, 9.6, 0.3) (13.4, 9.8, 0) 0.70
    • 辐射测量是高标准差的泊松过程,因此辐射测量常出现高随机性。在对静态辐射热点成像的过程中,光子的累积使得探测器对空间放射性分布的测量趋于稳定,从而能定位静态辐射热点。辐射热点运动则会导致空间中放射性分布发生变化,探测器获取的投影不再固定,此时重建图像往往会出现定位偏移或运动模糊[54],此种情形也会出现在光子计数率较低时,若探测器采集信号的有效强度远低于本底辐射或放射源运动速度较快,探测器无法在一定时间内累积足够的光子从而获得较稳定的放射性分布,准确定位移动辐射热点将会变得十分困难。针对上述动态辐射热点问题,提高重建图像的信噪比与低活度移动辐射源的定位精度。Zou等[55]在成像过程中引入了卷积长短时记忆网络(Convolution Long Short-Term Memory network, CLSTM),对互相关解码法重建的移动放射源辐射图像进行自适应叠加与优化,减少移动辐射热点的运动模糊,提高了定位精度,在针对行人等以近似恒定速度移动的核安全探测场景中是可行的。

      移动辐射热点的实时定位追踪算法流程如图17所示,主要包括三个步骤:1) 使用CLSTM结合卡尔曼滤波方法(Kalman Filter, KF)优化互相关解码法重建辐射图像,卡尔曼滤波的作用是对重建图像所获得的放射源定位进行滤波和平滑,解决放射源定位随机性跳动与定位不准确的问题;2) 引入YOLOv3算法[56]与DeepSort算法[57]定位追踪光学目标;3) 使用热点-运动目标匹配方法确定放射源携带者,在公共场所对行人进行实时放射性检测的应用场景中,可认为放射源携带者与放射源的运动是一致的,据此可对光学图像中携带者的定位与辐射图像中热点定位进行相互匹配,从而确定放射源携带者,实现对放射源的追踪。图18展示了应用跟踪算法对注射了131I的病人进行实时跟踪的效果,成像过程中,行人以正常速度运动,探测成像系统保持静止不动。

      该跟踪算法除了应用于辐射热点移动的情况,也适用于成像系统运动而放射源相对静止的场景。刘彦韬等在行车过程中对放射源实现了对放射源的跟踪测量[50],测试场景如图19(a)所示,13.45 mCi 137Cs点源放置于O点,车辆携带着高灵敏度γ相机以36 km/h的速度从M点运动到N点,OH = 25 m,热点从视场右边缘出现[图19(b)]到左边缘消失[图19(d)]的过程中,始终与放射源实际位置符合得很好。

    • 典型核素识别过程为:探测器获取γ能谱后,将特征峰与标准核素库中的所有核素相匹配,然后利用峰面积、相邻峰间距等信息对匹配到的峰位进行排除和验证,从而确定核素种类。为了准确获取特征峰位信息,还需要对γ能谱进行平滑、寻峰和扣除本底等预处理。识别过程中运算量大、效率低,识别的准确率也不高。此外,放射源的活度、探测距离等也会对能谱的测量造成影响[58],特别是在低计数和多核素混合的环境下,不利于核素的快速识别。神经网络正适用于处理具有一定规律但不易用直观的数学模型表达的数据,在非线性数据中的信息提取有着广泛应用前景,自然地适用于能谱处理。

      Li等[59]开发的基于先验能量信息的特征增强卷积神经网络(Convolutional Neural Network, CNN)核素识别算法已经部署在多款高能所研制的核辐射设备中。该算法注重于低计数能谱与多核素混合场景下的核素识别问题,算法结构如图20(a)所示,它由特征增强器以及一维卷积神经网络两部分组成。其中特征增强器包括了特征提取和归一化两种操作,实现对放射性核素先验能谱的信息提取与特征增强;一维卷积神经网络主要包含分离器(Splitter)和分类器(Classifer),实现对场景中核素的种类进行确认。分离器的设计主要是对经过特征增强器的预处理数据进行进一步的特征过滤及信息提取,因为特征增强器中的归一化操作可能会放大一些通道中来自本底的计数信息。同时,通过分离器内的分离卷积和金字塔结构,可有效减少网络的参数数量,加快训练过程,并实现对能谱进行不同范围的信息提取。最后,分类器包含一个全连接结构,根据分离器的输出进行归纳,并给出最终的核素分类情况[5960]。在经过模拟和实验生成训练、测试能谱、神经网络参数训练等步骤后,实现了对多种不同活度放射源的核素识别,混合识别的准确率达90%以上,图20(b)展示了对多种不同活度的核素的能谱与识别结果。快速核素识别功能的实现,对于编码孔径成像、康普顿成像都是非常有利的;尤其是对后者,通过该核素识别算法的应用,理论上能够加快康普顿成像过程中对核素的判断与选择,提高图像质量(信噪比)与成像速度。

      除了在对动态目标成像与核素识别外,神经网络在图像重建、部分编码视场(Partially coded field-of-view, PCFOV) 放射源识别定位、连续晶体定位等方面均用应用[6163]。在成像分析方面,利用机器学习和深度学习算法,进行噪声抑制和图像增强、智能目标检测和分割、智能图像重建和恢复、智能数据分析和解释等,实现自动化的成像分析和定位是一个值得探索的方向。

    • 智能化平台应具备智能化控制、可跨平台和跨设备支持、多信使辐射成像、智能化分析等基本功能,具备综合性探测能力。如在中国散裂中子源的建设过程中,为了能将隧道内的辐射监测等信息实时传输至主控室,搭建了专门的轨道,并研制了隧道环境巡检系统,如图21(a)所示,巡检系统内部集成了轨道机器人、编码γ相机、NaI探测器、红外成像模块、高清光学视相机等,能对散裂中子源加速器隧道内的环境进行全方位的监测。如图21(b)展示了巡检机器人在隧道内的工作场景,图21(c)则为其中的伽玛相机对隧道内辐射热点的成像结果。

      近年来,使用无人设备进行大范围、区域性辐射探测也受到越来越多的关注,特别是日本福岛核事故发生后[6468]。沉降物绘制(mapping fallout)、定位(localization)、烟羽跟踪(plume tracking) 是针对大范围探测平台的三大主要任务[4],同时需要满足能够快速机动与部署、探测灵敏度高等特点,具备综合性探测能力的智能化平台是解决此类问题的关键工具。在辐射探测成像与定位方面,中国科学院高能物理研究所研制了辐射智能化应急平台,如图22所示。其内部搭载了无人机载、车载、便携式γ相机等探测终端,具备对几公里到百公里范围内的监测区域实施高灵敏度、高机动的空地一体化辐射监测能力。图23展示了无人机携带小型机载核辐射成像终端的测试场景与成像结果,图24则为无人机填图测量与热点定位结果。

      辐射成像智能化技术的发展,使得人们可以更快、更精确地获取放射源的位置与核素信息,帮助人们从二维到三维、从地面到空中全方位获取感兴趣空间内的辐射(热点)分布情况,同时提高探测结果的准确性,降低误报的可能,为辐射环境评估、辐射防护、核应急决策提供更有效的信息。此外,辐射成像监控的分布式动态网络的构建同样值得探讨和发展,结合便携式辐射成像仪、定点通道安检成像仪、小型无人机搭载成像终端、投放式小球探测器等多种探测手段,通过4G或者5G通信联合组网,既能进行区域性监测,又能精确定位辐射热点,实现应急情况下全方位的机动测量。

    • 核辐射成像技术在辐射安全监测领域的应用日趋广泛,随着技术的不断进步,涌现了大量的智能化成像技术手段。文章分析总结了核辐射成像技术中如空间编码成像、时间编码成像、康普顿成像等主流成像方法的技术特点,诸多成像技术能够以可视化的方式对放射源进行远距离成像定位,结合光学图像可直观地给出放射源的具体分布位置。然而,现有成像方法只能指示放射源的方位角信息即射线入射方向,缺乏源项的距离信息,属于二维成像定位的技术范畴,如何实现放射源的三维定位成为辐射成像领域的研究焦点。源项监测场景并非都是静态目标,比如注射放射性药物后未代谢完成的病人、携带放射性材料的运输车辆或者船只以及核应急演练中可能出现的移动放射源目标等等,针对动态放射源目标,需要开发专门的动态辐射成像算法,用以实时追踪携带放射源的目标。核素识别是对放射性材料做定性判断的关键技术,尤其涉及到特殊核材料筛查的应用场景,存在屏蔽的条件下,出射的伽玛射线量往往较少,而人工智能算法为高灵敏度核素识别技术带来了契机。辐射场景往往复杂多变,比如大科学装置中的加速器束流热点、核应急中面临的复杂地势地貌等,多终端、多信使探测技术的融合互补可实现辐射信息的全方位探测,辐射探测智能化平台便是通过智能化控制、跨平台和跨设备应用技术使得辐射探测与成像技术更加高效。

      针对上述问题,在智能化辐射探测与成像技术研究方面,本文系统地介绍了中国科学院高能物理研究所在远场放射源三维定位、动态辐射热点成像、核素识别与分析等方面开展的探测方法和优化研究,展示了辐射成像技术搭载于机器人、车载、无人机载等平台的智能化应用,实现了针对实际应用场景的辐射热点方测搜寻、针对一定区域的辐射态势感知及可视化、城市级辐射监测及地图绘制等功能,实现了诸多先进成像技术的可行性及其在核与辐射安全监测的行业应用。本文展示了辐射成像智能化应用中的多个亮点工作,在辐射热点三维定位的研究方面,通过四目立体成像系统实现了远距离成像场景中放射源100 m距离下小于20%的测距精度,在室内场景中通过三维场景建模和辐射场构建实现了4枚放射源的三维定位。在人工智能算法的应用方面,实现了移动辐射热点的实时定位,可对25 m远处以36 km/h速度移动的14 mCi 137Cs放射源进行动态定位追踪。在神经网络核素识别算法研究方面,提出了一种基于先验能量信息的特征增强卷积神经网络,可在低计数和多核素混合的环境下,实现核素的快速识别。在辐射成像智能化平台研制方面,研制了可实现辐射热点检测的加速器隧道环境巡检机器人,并开发了一套车载式辐射应急监测平台,搭载有无人机载、车载、便携式γ相机等多个探测终端,可在辐射应急情况下实现城市区域的快速机动测量。

      核辐射成像技术是与科技发展、实际需求紧密结合的学科,探测器、电子学、人工智能、无人化平台的发展极大促进了核辐射成像技术的发展,使得辐射成像技术更加高效、准确和智能化。未来的发展方向可能包括更高分辨率的探测器、更先进的数据处理算法、更强大的人工智能模型,以及更完善的三维可视化技术。从实际应用需求出发,辐射探测与成像领域中未来的研究热点和方向可能涉及到:1) 通过远距离测量实现放射性材料的量化,结合本文中的放射源三维定位技术,有望实现远距离测量条件下放射源活度的估算,从而减少近距离接触测量所带来的辐照损伤;2) 三维辐射场的动态实时构建,通过高灵敏度成像技术以及源项动态重建算法,可对核设施等场景做辐射水平的三维动态评估,从而规划辐照最优化工作路径,在核设施大修或者退役期间有效降低工作人员的受照剂量。这些进步将进一步提升核与辐射安全管理的水平,此外还可预期在海洋环境监测、天文学观测、医疗影像等领域发挥重要的作用。

Reference (68)

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return