搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

多重库仑散射对小尺度物体缪子透射成像精度的影响

张建鸣 李志伟 刘芳 李景太 冒鑫 程雅苹 庞捷 冯鑫茁 倪四道 欧阳晓平 韩然

引用本文:
Citation:

多重库仑散射对小尺度物体缪子透射成像精度的影响

张建鸣, 李志伟, 刘芳, 李景太, 冒鑫, 程雅苹, 庞捷, 冯鑫茁, 倪四道, 欧阳晓平, 韩然

Influence of multiple Coulomb scattering on accuracy of muon transmission imaging of small-scale matter

Zhang Jian-Ming, Li Zhi-Wei, Liu Fang, Li Jing-Tai, Mao Xin, Cheng Ya-Ping, Pang Jie, Feng Xin-Zhuo, Ni Si-Dao, Ouyang Xiao-Ping, Han Ran
PDF
HTML
导出引用
  • 缪子透射成像方法是一种基于宇宙线缪子穿过目标物体前后的通量变化, 进而获得其内部密度结构的无损探测成像方法. 缪子透射成像方法假设缪子在低原子序数物质中沿直线运动, 但实际上多重库仑散射作用会使缪子路径一定程度上偏离直线, 有可能对成像精度造成影响. 为此, 本文使用Geant4软件包开展了缪子透射成像蒙特卡罗模拟, 针对数米尺度多种密度结构的模型, 定量分析了打开和关闭多重库仑散射物理过程时对成像精度的影响. 结果表明: 对于数米尺度标准岩石物质, 缪子透射成像方法能够很好地恢复内部密度异常几何特征;但多重库仑散射作用对目标物体内部区域近垂向缪子通量造成的偏差可达5%, 而在目标物体边界区域的偏差可达13%. 因此, 需要宇宙线缪子透射成像中考虑多重库仑散射作用的影响, 以获得数米尺度目标物体更准确的绝对密度值成像结果.
    The muon transmission imaging method is a non-destructive detection imaging method and can obtain the internal density structure of the target object by analyzing the flux change of cosmic ray muon before and after passing through the target object. This method assumes that muon travels along a straight line in a lowatomic number material. However, the multiple Coulomb scattering causes the muon deviate from the straight line to a certain extent when penetrating the material, which may have a certain influence on the accuracy of muon transmission imaging. This study uses the Geant4 software package to carry out Monte Carlo simulation of muon transmission imaging. Object models with multiple density structures and several meters are used to analyze the effect of multiple scattering on the accuracy of transmission imaging. In the experiment, we set up a model in which muons of different energy vertically pass through a rock with a certain thickness, we can intuitively see the influence of scattering on the penetration path of muons. By setting up rock-water block models with thickness in a range of 0.8, 2.4 and 4.0 m, the effect of Coulomb scattering on the transmission imaging of small-scale material muons is analyzed. The results suggest that the muon transmission imaging method can well restore the geometric shape and spatial distribution characteristics of density anomalies for objects with several-meter scale for standard rock materials with a scale of several meters. However, the flux deviation caused by multiple Coulomb scattering on the muons in the near-vertical direction can reach 5%, and up to 13% in the boundary areas of the standard rock and air. We limit the scattering angle of the muon, and select the muon with a scattering angle of less than or equal to 1° for imaging. The results of transmission imaging by using the selected muon have improved. The image does not have the illusion of an abnormally increased flux around the model caused by scattering, but the muon flux in the model area is reduced even more, thereby affecting the accuracy of restoring the absolute density of an object using flux differences. Therefore, the effects of multiple Coulomb scattering should be considered for recovering more accurate absolute density in small-scale muon transmission imaging study.
      通信作者: 李志伟, zwli@whigg.ac.cn ; 刘芳, liuf@ncepu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 42174076, 41974064)资助的课题.
      Corresponding author: Li Zhi-Wei, zwli@whigg.ac.cn ; Liu Fang, liuf@ncepu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 42174076, 41974064).
    [1]

    Gaisser T K 1990 Cosmic Rays and Particles Physics (Cambridge: Cambridge University Press) p126

    [2]

    Lesparre, Gibert N D, Marteau J, Komorowski J C, Nicollin F, Coutant O 2012 Geophys. J. Inter. 190 1008Google Scholar

    [3]

    Zhang Z X, Enqvist T, Holma M, Kuusiniemi P 2020 Rock Mech. Rock Eng. 53 4893Google Scholar

    [4]

    Hou L, Huo Y, Zuo W, Yao Q, Yang J, Zhang Q 2021 Nucl. Eng. Technol. 53 208Google Scholar

    [5]

    Gilboy W B, Jenneson P M, Simons S J R, Stanley S J, Rhodes D 2007 Nucl. Instrum. Meth. A 580 785Google Scholar

    [6]

    Borozdin K, Hogan G, Morris C, Priedhorsky W C, Saunders A, Schultz L J, Teasdale M E 2003 Nature 422 277

    [7]

    Schultz L J, Borozdin K N, Gomez J J, Hogan G E, McGill J A, Morris C L, Priedhorsky W C, Saunders A, Teasdale M E 2004 Nucl. Instrum. Meth. A 519 687Google Scholar

    [8]

    Liu Y Y, Chen Z Q, Zhao Z R, Zhang L, Wang Z T 2009 Tsinghua Sci. Technol. 14 313Google Scholar

    [9]

    Schultz L J, Blanpied G S, Borozdin K N, Fraser A M, Hengartner N W, Klimenko A V, Morris C L, Orum C, Sossong M J 2007 IEEE Trans. Image Process. 16 1985Google Scholar

    [10]

    Wang G, Schultz L J, Qi J 2009 IEEE Trans. Image Process. 18 1080Google Scholar

    [11]

    Gary B, Sankaran K, Dustin D, Craig M, Isabelle B, Michael S, Shawn M K, Beth N 2015 Nucl. Instrum. Meth. A 784 352Google Scholar

    [12]

    George E P 1955 Commonw. Eng. 1955 455

    [13]

    Alvarez L W, Anderson J A, Bedwei F E, Burkhard J, Fakhry A, Girgis A, Goneid A, Hassan F, Iverson D, Lynch G, Miligy Z, Moussa AH, Sharkawi M, Yazolino L 1970 Science 167 832Google Scholar

    [14]

    Morishima K, Kuno M, Nishio A, et al. 2017 Nature 552 386Google Scholar

    [15]

    Han R, Yu Q, Li Z W, Li J T, Cheng Y P, Liao B, Jiang L X, Ni S D, Liu T F, Wang Z 2020 J. Instrum. 15 06019Google Scholar

    [16]

    Li, J T, Li Z W, Han R, Cheng Y P, Mao X, Yu L, Feng X, Liu B, Jiang L, Ouyang X P 2022 J. Instrum. 17 0502Google Scholar

    [17]

    Cheng Y P, Han R, Li Z W, Li J T, Mao X, Dou W Q, Feng X Z, Ouyang X P, Liao B, Liu F, Huang L 2022 Nucl. Sci. Tech. 33 88Google Scholar

    [18]

    苏宁, 刘圆圆, 王力, 程建平 2022 物理学报 71 064201Google Scholar

    Su N, Liu Y Y, Wang L, Cheng J P 2022 Acta Phys. Sin. 71 064201Google Scholar

    [19]

    何伟波 2019 博士学位论文 (合肥: 中国科学技术大学)

    He W B 2019 Ph. D. Dissertation (Hefei: University of Science and Technology of China) (in Chinese)

    [20]

    Kume N, Miyadera H, Morris C L, Bacon J, Borozdin K N, Durham J M, Fuzita K, Guardincerri E, Izumi M, Nakayama K, Saltus M, Sugita T, Takakura K, Yoshioka K 2016 J. Instrum. 11 09008Google Scholar

    [21]

    Borozdin K, Greene S, Lukić Z, Milner E, Miyadera H, Morris C, Perry J 2012 Phys. Rev. Lett. 109 152501Google Scholar

    [22]

    Takamatsu K, Takegami H, Ito C, Suzuki K, Ohnuma H, Hino R, Okumura T 2015 Ann. Nucl. Energy 78 166Google Scholar

    [23]

    Bethe H A 1953 Phys Rev. 89 1256Google Scholar

    [24]

    Case C T, Battle E L 1968 Phys. Rev. 169 201Google Scholar

    [25]

    Hagmann C, Lange D, Wright D 2007 IEEE Nuclear Science Symposium Conference Record Hawaii, United States October 26–November 3, 2007 p1143

    [26]

    Groom D E, Mokhov N V, Striganov S I 2001 Atom. Data Nucl. Data 78 183Google Scholar

  • 图 1  缪子穿透物体的路径

    Fig. 1.  Path of the muon through the object.

    图 2  CRY随机抽样产生缪子的能量分布

    Fig. 2.  Energy distribution of muons generated by random sampling of CRY.

    图 3  散射模块验证模型(缪子为单能点源, 垂直穿过待测物体)

    Fig. 3.  Verification model of scattering module (A single energy muon at a point source, passing vertically through the object to be measured).

    图 4  单能缪子竖直穿透的标准岩石厚度分布 (a) 0.5 GeV; (b) 1.0 GeV; (c) 10.0 GeV

    Fig. 4.  Standard rock thickness distribution for vertical penetration of muons with different energies: (a) 0.5 GeV; (b) 1.0 GeV; (c) 10.0 GeV.

    图 5  不同能量缪子能竖直穿透的标准岩石的平均厚度 (误差棒给出1σ误差)

    Fig. 5.  Average thickness of standard rocks penetrated vertically by muons of different energies (Error bar is 1σ).

    图 6  多重散射特征分析模型 (100个1.5 GeV的缪子, 垂直穿过3 m厚度的标准岩石(红色线条), 蓝色为缪子径迹)

    Fig. 6.  Multiple scattering feature analysis model (100 muonsof 1.5 GeV vertically passes through a standard rock with a thickness of 3 m (red line), the blue is the muon track).

    图 7  岩石-空腔模型透射成像模拟

    Fig. 7.  Rock-cavity model transmission imaging simulation.

    图 8  点源缪子竖直穿过1 m厚标准岩石的出射位置分布 (a) 0.7 GeV; (b) 1.0 GeV; (c) 3.0 GeV

    Fig. 8.  Distribution of emission positions of point source muons vertically passing through 1 m thick standard rock: (a) 0.7 GeV; (b) 1.0 GeV; (c) 3.0 GeV

    图 9  单能点源缪子穿过不同厚度的标准岩石的多重库仑散射影响 (a) 路径平面平均偏移长度(误差棒给出0.5σ误差); (b) 累计散射角

    Fig. 9.  Multiple Coulomb scattering effects of single-energy point source muons passing through standard rocks with different thicknesses: (a) Average offset length of the path plane (the error bar is 0.5σ); (b) cumulative scattering angle.

    图 10  缪子穿过不同厚度模型前后各像素上计数的比值(俯视图) (a), (d), (g) 0.8 m; (b), (e), (h) 2.4 m; (c), (f), (i) 4.0 m (黑色实线为岩石模型的外边界, 黑色虚线为岩石与空腔的分界)

    Fig. 10.  The ratio of the counts on each pixel before and after the muons pass through the models with different thicknesses (top view): (a), (d), (g) 0.8 m; (b), (e), (h) 2.4 m; (c), (f), (i) 4.0 m (Solid black line is the outer boundary of the rock model, the dashed black line is the boundary between the rock and the cavity).

    图 11  不同厚度模型切线经过像素的计数比值分布 (a), (d) 0.8 m; (b), (e) 2.4 m; (c), (f) 4.0 m (A-A' 线在岩石和外边界附近, B-B' 切线横穿了岩石-空腔模型 (图10 (a)(c)))

    Fig. 11.  Distribution of count ratios of the tangent lines passing through pixels for different thickness models: (a), (d) 0.8 m; (b), (e) 2.4 m; (c), (f) 4.0 m (A-A' line is near the rock and the outer boundary, and B-B' tangent traverses the rock-cavity model (Fig. 10 (a)-(c))).

    图 12  经过缪子散射角筛选后不同厚度模型缪子透射成像结果(俯视图) (a) 0.8 m; (b) 2.4 m; (c) 4.0 m

    Fig. 12.  Transmission imaging results of the model of different thicknesses after the screening of the muon scattering angle (top view): (a) 0.8 m; (b) 2.4 m; (c) 4.0 m.

    表 1  3 GeV的缪子穿过10 cm厚不同材料的散射角

    Table 1.  Multiple scattering for 3 GeV muons passing through 10 cm of various materials.

    材料X0/cm理论θ0/mrad模拟θ0/mrad相对偏差/%
    Al8.994.804.43–7.8
    Fe1.8011.410.6–7.0
    Cu1.4812.612.60
    U0.3129.331.78.2
    下载: 导出CSV
  • [1]

    Gaisser T K 1990 Cosmic Rays and Particles Physics (Cambridge: Cambridge University Press) p126

    [2]

    Lesparre, Gibert N D, Marteau J, Komorowski J C, Nicollin F, Coutant O 2012 Geophys. J. Inter. 190 1008Google Scholar

    [3]

    Zhang Z X, Enqvist T, Holma M, Kuusiniemi P 2020 Rock Mech. Rock Eng. 53 4893Google Scholar

    [4]

    Hou L, Huo Y, Zuo W, Yao Q, Yang J, Zhang Q 2021 Nucl. Eng. Technol. 53 208Google Scholar

    [5]

    Gilboy W B, Jenneson P M, Simons S J R, Stanley S J, Rhodes D 2007 Nucl. Instrum. Meth. A 580 785Google Scholar

    [6]

    Borozdin K, Hogan G, Morris C, Priedhorsky W C, Saunders A, Schultz L J, Teasdale M E 2003 Nature 422 277

    [7]

    Schultz L J, Borozdin K N, Gomez J J, Hogan G E, McGill J A, Morris C L, Priedhorsky W C, Saunders A, Teasdale M E 2004 Nucl. Instrum. Meth. A 519 687Google Scholar

    [8]

    Liu Y Y, Chen Z Q, Zhao Z R, Zhang L, Wang Z T 2009 Tsinghua Sci. Technol. 14 313Google Scholar

    [9]

    Schultz L J, Blanpied G S, Borozdin K N, Fraser A M, Hengartner N W, Klimenko A V, Morris C L, Orum C, Sossong M J 2007 IEEE Trans. Image Process. 16 1985Google Scholar

    [10]

    Wang G, Schultz L J, Qi J 2009 IEEE Trans. Image Process. 18 1080Google Scholar

    [11]

    Gary B, Sankaran K, Dustin D, Craig M, Isabelle B, Michael S, Shawn M K, Beth N 2015 Nucl. Instrum. Meth. A 784 352Google Scholar

    [12]

    George E P 1955 Commonw. Eng. 1955 455

    [13]

    Alvarez L W, Anderson J A, Bedwei F E, Burkhard J, Fakhry A, Girgis A, Goneid A, Hassan F, Iverson D, Lynch G, Miligy Z, Moussa AH, Sharkawi M, Yazolino L 1970 Science 167 832Google Scholar

    [14]

    Morishima K, Kuno M, Nishio A, et al. 2017 Nature 552 386Google Scholar

    [15]

    Han R, Yu Q, Li Z W, Li J T, Cheng Y P, Liao B, Jiang L X, Ni S D, Liu T F, Wang Z 2020 J. Instrum. 15 06019Google Scholar

    [16]

    Li, J T, Li Z W, Han R, Cheng Y P, Mao X, Yu L, Feng X, Liu B, Jiang L, Ouyang X P 2022 J. Instrum. 17 0502Google Scholar

    [17]

    Cheng Y P, Han R, Li Z W, Li J T, Mao X, Dou W Q, Feng X Z, Ouyang X P, Liao B, Liu F, Huang L 2022 Nucl. Sci. Tech. 33 88Google Scholar

    [18]

    苏宁, 刘圆圆, 王力, 程建平 2022 物理学报 71 064201Google Scholar

    Su N, Liu Y Y, Wang L, Cheng J P 2022 Acta Phys. Sin. 71 064201Google Scholar

    [19]

    何伟波 2019 博士学位论文 (合肥: 中国科学技术大学)

    He W B 2019 Ph. D. Dissertation (Hefei: University of Science and Technology of China) (in Chinese)

    [20]

    Kume N, Miyadera H, Morris C L, Bacon J, Borozdin K N, Durham J M, Fuzita K, Guardincerri E, Izumi M, Nakayama K, Saltus M, Sugita T, Takakura K, Yoshioka K 2016 J. Instrum. 11 09008Google Scholar

    [21]

    Borozdin K, Greene S, Lukić Z, Milner E, Miyadera H, Morris C, Perry J 2012 Phys. Rev. Lett. 109 152501Google Scholar

    [22]

    Takamatsu K, Takegami H, Ito C, Suzuki K, Ohnuma H, Hino R, Okumura T 2015 Ann. Nucl. Energy 78 166Google Scholar

    [23]

    Bethe H A 1953 Phys Rev. 89 1256Google Scholar

    [24]

    Case C T, Battle E L 1968 Phys. Rev. 169 201Google Scholar

    [25]

    Hagmann C, Lange D, Wright D 2007 IEEE Nuclear Science Symposium Conference Record Hawaii, United States October 26–November 3, 2007 p1143

    [26]

    Groom D E, Mokhov N V, Striganov S I 2001 Atom. Data Nucl. Data 78 183Google Scholar

  • [1] 潘新宇, 毕筱雪, 董政, 耿直, 徐晗, 张一, 董宇辉, 张承龙. 叠层相干衍射成像算法发展综述. 物理学报, 2023, 72(5): 054202. doi: 10.7498/aps.72.20221889
    [2] 赵富, 胡渝曜, 王鹏, 刘军. 偏振复用散射成像. 物理学报, 2023, 72(15): 154201. doi: 10.7498/aps.72.20230551
    [3] 苏宁, 刘圆圆, 王力, 程建平. 秦始皇陵地宫宇宙射线缪子吸收成像模拟研究. 物理学报, 2022, 71(6): 064201. doi: 10.7498/aps.71.20211582
    [4] 霍勇刚, 严江余, 张全虎. 缪子多模态成像图像质量分析. 物理学报, 2022, 71(2): 021401. doi: 10.7498/aps.71.20211083
    [5] 许文慧, 宁守琮, 张福才. 部分相干衍射成像综述. 物理学报, 2021, 70(21): 214201. doi: 10.7498/aps.70.20211020
    [6] 孙雪莹, 刘飞, 段景博, 牛耕田, 邵晓鹏. 基于散斑光场偏振共模抑制性的宽谱散射成像技术. 物理学报, 2021, 70(22): 224203. doi: 10.7498/aps.70.20210703
    [7] 霍勇刚, 严江余, 张全虎. 缪子多模态成像图像质量分析. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211083
    [8] 严江余, 张全虎, 霍勇刚. 基于散射和次级诱发中子的缪子多模态成像. 物理学报, 2021, 70(19): 191401. doi: 10.7498/aps.70.20210804
    [9] 庄佳衍, 陈钱, 何伟基, 冒添逸. 基于压缩感知的动态散射成像. 物理学报, 2016, 65(4): 040501. doi: 10.7498/aps.65.040501
    [10] 张美, 张显鹏, 李奎念, 盛亮, 袁媛, 宋朝晖, 李阳. 中子散射成像探测角分辨研究. 物理学报, 2015, 64(4): 042801. doi: 10.7498/aps.64.042801
    [11] 古宇飞, 闫镔, 李磊, 魏峰, 韩玉, 陈健. 基于全变分最小化和交替方向法的康普顿散射成像重建算法. 物理学报, 2014, 63(1): 018701. doi: 10.7498/aps.63.018701
    [12] 相里斌, 张文喜, 伍洲, 吕笑宇, 李杨, 周志盛, 孔新新. 相干场成像技术接收镜精度对传递函数的影响. 物理学报, 2013, 62(22): 224201. doi: 10.7498/aps.62.224201
    [13] 陈友华, 王召巴, 王志斌, 张瑞, 王艳超, 王冠军. 弹光调制型成像光谱偏振仪中的高精度偏振信息探测研究. 物理学报, 2013, 62(6): 060702. doi: 10.7498/aps.62.060702
    [14] 赵吉祯, 欧阳晓平, 盛亮, 魏福利, 张美. 脉冲辐射成像绝对测量方法研究. 物理学报, 2013, 62(22): 225203. doi: 10.7498/aps.62.225203
    [15] 刘诚, 潘兴臣, 朱健强. 基于光栅分光法的相干衍射成像. 物理学报, 2013, 62(18): 184204. doi: 10.7498/aps.62.184204
    [16] 陈鹤, 于斌, 陈丹妮, 李恒, 牛憨笨. 超衍射成像中双螺旋点扩展函数的三维定位精度. 物理学报, 2013, 62(14): 144201. doi: 10.7498/aps.62.144201
    [17] 李晓春, 高俊丽, 刘绍娥, 周科朝, 黄伯云. 无序对二维声子晶体平板负折射成像的影响. 物理学报, 2010, 59(1): 376-380. doi: 10.7498/aps.59.376
    [18] 赵宗清, 丁永坤, 郝轶聃, 袁永腾, 蒲以康. 瞄准精度对中子半影成像的影响. 物理学报, 2008, 57(9): 5756-5760. doi: 10.7498/aps.57.5756
    [19] 乐贵明, 韩延本. 用银河宇宙线数据分析1991年3月24日CME的结构. 物理学报, 2005, 54(1): 467-470. doi: 10.7498/aps.54.467
    [20] 郏东耀, 丁天怀. 皮棉杂质透射检测及成像目标增强方法. 物理学报, 2005, 54(9): 4058-4064. doi: 10.7498/aps.54.4058
计量
  • 文章访问数:  2641
  • PDF下载量:  75
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-09-13
  • 修回日期:  2022-10-13
  • 上网日期:  2022-10-19
  • 刊出日期:  2023-01-20

/

返回文章
返回