搜索

x

留言板

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

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

高能质子照相中基于角度准直器设计的理论研究

陈锋 许海波 郑娜 贾清刚 佘若谷 李兴娥

引用本文:
Citation:

高能质子照相中基于角度准直器设计的理论研究

陈锋, 许海波, 郑娜, 贾清刚, 佘若谷, 李兴娥

Theoretical study of angle-cut collimator based design in high-energy proton radiography

Chen Feng, Xu Hai-Bo, Zheng Na, Jia Qing-Gang, She Ruo-Gu, Li Xing-E
PDF
HTML
导出引用
  • 角度准直器在高能质子照相中有着重要作用, 既可以利用准直器提高图像对比度, 又能通过二次成像实现材料诊断及密度重建, 因此减小通过准直器后通量值的误差具有重要意义. 本文通过理论分析, 提出了一种高能质子照相中准直器设计的方法, 通过Geant4程序建立了1.6 GeV的质子成像系统, 该系统分别使用理想准直器、拉伸型准直器和利用该方法设计的准直器, 并对比通过客体后的通量分布. 结果表明, 在使用理想准直器和该方法设计的角度准直器时, 二者得到的客体的通量分布符合较好, 而使用拉伸型准直器时, 与使用理想准直器得到的结果相差较大. 因此利用理想准直器方法设计的准直器可以很好地减小通量误差.
    The angle-cut collimator plays an important role in high-energy proton radiography. By using the collimator, the image contrast can be improved, and the material diagnosis and density reconstruction can be realized through secondary imaging. As all these techniques depend on the flux value, it is of great significance to reduce the error of the detected flux value. The ideal collimator is a much thin surface, but thick enough to block protons outside the collimation region. It is designed by stretching the aperture of the collimation plane. The shape is cylindrical, and it will increase the error of the flux value with the angle truncation. The initial bunch is defined and the phase diagram of the bunch within the angle-cut is ideal in the theoretical model. The equation of designing the collimator is given by theoretical analysis. It is given by the transfer matrix, the radius of the object and the angle-cut. The pore structure is oval-shaped by calculating and simulating. The proton imaging system of 1.6 GeV is established by Geant4 program, and the detector is ideal. The round copper plate and the concentric spheres are chosen as objects respectively. The parameters of the designed collimator is given by this method. The ideal collimator, tensile collimator and designed collimator are used in simulation, the radius of object is 5 cm and the angle-cut is 2 mrad and 3.5 mrad. The results show that when using the ideal and the designed angle-cut collimator, the flux distributions are in good agreement, while when using the tensile collimator, the result is quite different from that obtained by using the ideal collimator. Therefore, the collimator designed by this method can effectively reduce the error of the detected flux value.
      通信作者: 许海波, 13641017929@163.com
    • 基金项目: 国家自然科学基金(批准号: 11675021)和国家自然科学青年科学基金(批准号: 11805018)资助的课题
      Corresponding author: Xu Hai-Bo, 13641017929@163.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11675021) and the Young Scientists Fund of the National Natural Science Foundation of China (Grant No. 11805018)
    [1]

    Gavton A, Morris C L, Ziock H J, et al. 1996 Los Alamos National Report 96 420

    [2]

    Mottershead C T, Zumbro J D 1997 Proceedings of the 1997 Particle Accelerator Conference Vancouver B C, Canada, May 12–16, 1997 p1397

    [3]

    Jason A J, Barlow D B, Blind B, Kelley J P, Lysenko W P, Neri F, Walstrom P L, Waynert J, Schulze M 2001 Proceedings of the 2001 Particle Accelerator Conference Chicago, USA, June 18–22, 2001 p3374

    [4]

    King N S P, Ables E, Adams K, et al. 1999 Nucl. Instrum. Meth. Phys. Res., Sect. A 424 84

    [5]

    Rigg P A, Schwartz C L, Hixson R S, Hogan G E, Kwiatkowski K K, Mariam F G, Marr-Lyon M, Merrill F E, Morris C L, Rightly P, Sauders A, Tupa D 2008 Phys. Rev. B 77 220101Google Scholar

    [6]

    Morris C L, Ables E, Alrick K R, et al. 2011 J. Appl. Phys. 109 104905Google Scholar

    [7]

    Antipov Y M, Afonin A G, Vasilevskii A V, et al. 2010 Instrum. Exp. Tech. 53 319Google Scholar

    [8]

    Golubev A A, Demidov V S, Demidova E V, Dudin S V, Kantsyrev A V, Kolesnikov S A, Mintsev V B, Smirnov G N, Turtikov V I, Utikin A V, Fortov V E, Sharkov B Y 2010 Tech. Phys. Lett. 36 177Google Scholar

    [9]

    Varentsov D, Antonov O, Bakhmutova A, et al. 2016 Rev. Sci. Instrum. 87 023303Google Scholar

    [10]

    Yang J J, Zhen X, Wei S M, Lv Y L, Wang F, Zhang Y W, Wen L P, Liu J Y, Cai H R, Ge T, Zhang S P, Cao L, Zhang T J, Li Z G 2016 CYC2016 Proceedings of the 21 st International Conference on Cyclotrons and their Applications Zurich, Switzerland, September 11–16, 2016 p401

    [11]

    Merrill F E 2015 Rev. Accel. Sci. Technol. 8 165Google Scholar

    [12]

    Sheng L N, Zhao Y T, Yang G J, Wei T 2014 Laser Part. Beams 32 651Google Scholar

    [13]

    Wei T, Yang G J, Li Y D, et al. 2014 Chin. Phys. C 38 087003Google Scholar

    [14]

    Aufderheide M B, Park H, Hartouni E P 1999 AIP Conference Proceedings Sydney, Australia, June 28–July 2, 1999 p497

    [15]

    Zumbro J D, Acuff A, Bull J S, et al. 2005 Radiat. Prot. Dosim. 117 447

    [16]

    Fesseha G, Mariam, John P 2011 2011 High-Energy-Proton Microscopy Workshop Summary Report New Mexico, USA, October 27–28, 2011 p54

    [17]

    Antipov Y M, Afonin A G, Gusev I A, et al. 2013 At. Energ. 114 359Google Scholar

    [18]

    Varentsov D, Bogdanov A, Demidov V S 2013 Phys. Medica 29 208Google Scholar

    [19]

    Yan Y, Sheng L N, Huang Z W, et al. 2015 Laser Part. Beams 33 439Google Scholar

    [20]

    Kantsyrev A V, Skoblyakov A V, Bogdanov A V 2018 J. Phys.: Conf. Ser. 946 012019Google Scholar

    [21]

    Agostinelli S, Allison J, Amako K A, et al. 2003 Nucl. Instrum. Meth. Phys. Res. Sect. A 506 250Google Scholar

    [22]

    Allison J, Amako K, Apostolakis J, et al. 2006 IEEE Trans. Nucl. Sci. 53 270Google Scholar

  • 图 1  质子成像系统示意图

    Fig. 1.  Schematic diagram of proton radiography system.

    图 2  一定截断角以内的质子束团在通过准直空间时的形状变化 (a) z = 0 m; (b) z = 0.6 m; (c) z = 1.2 m; (d) z = 1.8 m; (e) z = 2.4 m

    Fig. 2.  Shape changed of proton bunch within certain angle-cuts as it passes through the collimation space: (a) z = 0 m; (b) z = 0.6 m; (c) z = 1.2 m; (d) z = 1.8 m; (e) z = 2.4 m.

    图 3  质子成像系统参数示意图

    Fig. 3.  Diagram of parameters of proton radiography system.

    图 4  截断角为2 mrad、客体尺寸为5 cm时的目标束团边界线

    Fig. 4.  Boundary lines of the target bunch when the angle-cut is 2 mrad and the object size is 5 cm.

    图 5  端口处的孔径值随截断角的变化 (a) $z = 20\;{\rm{cm}}$; (b) $z = 120\;{\rm{cm}}$

    Fig. 5.  Aperture size varies with the angle-cut at the ports: (a) $z = 20\;{\rm{cm}}$; (b) $z = 120\;{\rm{cm}}$.

    图 6  准直器孔径的形状 (a) x-y平面; (b) y-z平面

    Fig. 6.  Shape of aperture of the collimator: (a) x-y plane; (b) y-z plane.

    图 7  客体示意图 (a)铜板; (b)同心球体

    Fig. 7.  Diagram of the object: (a) The round copper plate; (b) the concentric spheres.

    图 8  通过铜板的通量分布 (a)截断角为2 mrad; (b)截断角为3.5 mrad

    Fig. 8.  Flux distribution after passing the round copper plate: (a) Angle-cut of 2 mrad; (b) angle-cut of 3.5 mrad.

    图 9  通过同心球的通量分布 (a)截断角为2 mrad; (b)截断角为3.5 mrad

    Fig. 9.  Flux distribution after passing the concentric spheres: (a) Angle-cut of 2 mrad; (b) angle-cut of 3.5 mrad.

    表 1  1.6 GeV质子成像系统参数

    Table 1.  Parameters of the proton radiography system of 1.6 GeV.

    s/ml/mG/T·m–1t/m
    1.20.88.090.5
    下载: 导出CSV

    表 2  准直器的孔径参数

    Table 2.  Aperture parameters of the collimator.

    截断角/mrad准直器类型前端/cm 后端/cm厚度/m外半径/m材料
    xy xy
    2理想型 0.6110–73Al
    拉伸型0.61 0.611 & 0.43W
    设计型2.042.46 0.6113W
    3.5理想型 1.0710–73Al
    拉伸型1.07 1.071 & 0.43W
    设计型2.343.08 1.0713W
    下载: 导出CSV
  • [1]

    Gavton A, Morris C L, Ziock H J, et al. 1996 Los Alamos National Report 96 420

    [2]

    Mottershead C T, Zumbro J D 1997 Proceedings of the 1997 Particle Accelerator Conference Vancouver B C, Canada, May 12–16, 1997 p1397

    [3]

    Jason A J, Barlow D B, Blind B, Kelley J P, Lysenko W P, Neri F, Walstrom P L, Waynert J, Schulze M 2001 Proceedings of the 2001 Particle Accelerator Conference Chicago, USA, June 18–22, 2001 p3374

    [4]

    King N S P, Ables E, Adams K, et al. 1999 Nucl. Instrum. Meth. Phys. Res., Sect. A 424 84

    [5]

    Rigg P A, Schwartz C L, Hixson R S, Hogan G E, Kwiatkowski K K, Mariam F G, Marr-Lyon M, Merrill F E, Morris C L, Rightly P, Sauders A, Tupa D 2008 Phys. Rev. B 77 220101Google Scholar

    [6]

    Morris C L, Ables E, Alrick K R, et al. 2011 J. Appl. Phys. 109 104905Google Scholar

    [7]

    Antipov Y M, Afonin A G, Vasilevskii A V, et al. 2010 Instrum. Exp. Tech. 53 319Google Scholar

    [8]

    Golubev A A, Demidov V S, Demidova E V, Dudin S V, Kantsyrev A V, Kolesnikov S A, Mintsev V B, Smirnov G N, Turtikov V I, Utikin A V, Fortov V E, Sharkov B Y 2010 Tech. Phys. Lett. 36 177Google Scholar

    [9]

    Varentsov D, Antonov O, Bakhmutova A, et al. 2016 Rev. Sci. Instrum. 87 023303Google Scholar

    [10]

    Yang J J, Zhen X, Wei S M, Lv Y L, Wang F, Zhang Y W, Wen L P, Liu J Y, Cai H R, Ge T, Zhang S P, Cao L, Zhang T J, Li Z G 2016 CYC2016 Proceedings of the 21 st International Conference on Cyclotrons and their Applications Zurich, Switzerland, September 11–16, 2016 p401

    [11]

    Merrill F E 2015 Rev. Accel. Sci. Technol. 8 165Google Scholar

    [12]

    Sheng L N, Zhao Y T, Yang G J, Wei T 2014 Laser Part. Beams 32 651Google Scholar

    [13]

    Wei T, Yang G J, Li Y D, et al. 2014 Chin. Phys. C 38 087003Google Scholar

    [14]

    Aufderheide M B, Park H, Hartouni E P 1999 AIP Conference Proceedings Sydney, Australia, June 28–July 2, 1999 p497

    [15]

    Zumbro J D, Acuff A, Bull J S, et al. 2005 Radiat. Prot. Dosim. 117 447

    [16]

    Fesseha G, Mariam, John P 2011 2011 High-Energy-Proton Microscopy Workshop Summary Report New Mexico, USA, October 27–28, 2011 p54

    [17]

    Antipov Y M, Afonin A G, Gusev I A, et al. 2013 At. Energ. 114 359Google Scholar

    [18]

    Varentsov D, Bogdanov A, Demidov V S 2013 Phys. Medica 29 208Google Scholar

    [19]

    Yan Y, Sheng L N, Huang Z W, et al. 2015 Laser Part. Beams 33 439Google Scholar

    [20]

    Kantsyrev A V, Skoblyakov A V, Bogdanov A V 2018 J. Phys.: Conf. Ser. 946 012019Google Scholar

    [21]

    Agostinelli S, Allison J, Amako K A, et al. 2003 Nucl. Instrum. Meth. Phys. Res. Sect. A 506 250Google Scholar

    [22]

    Allison J, Amako K, Apostolakis J, et al. 2006 IEEE Trans. Nucl. Sci. 53 270Google Scholar

  • [1] 向鹏程, 蔡聪波, 王杰超, 蔡淑惠, 陈忠. 基于深度神经网络的时空编码磁共振成像超分辨率重建方法. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211754
    [2] 游志明, 王洁, 高勇, 范佳锟, 张静, 胡耀程, 王盛, 许章炼, 张琦. 超级质子-质子对撞机束屏内气体密度演化规律研究. 物理学报, 2021, 70(16): 166802. doi: 10.7498/aps.70.20201594
    [3] 陈锋, 郝建红, 许海波. 考虑磁透镜边缘场的质子成像系统优化设计. 物理学报, 2021, 70(2): 022901. doi: 10.7498/aps.70.20201141
    [4] 李曜均, 岳东宁, 邓彦卿, 赵旭, 魏文青, 葛绪雷, 远晓辉, 刘峰, 陈黎明. 相对论强激光与近临界密度等离子体相互作用的质子成像. 物理学报, 2019, 68(15): 155201. doi: 10.7498/aps.68.20190610
    [5] 朱炳辉, 杨爱香, 牛书通, 陈熙萌, 周旺, 邵剑雄. 100 keV质子与低高能质子在绝缘微孔中输运特性的对比分析. 物理学报, 2018, 67(1): 013401. doi: 10.7498/aps.67.20171701
    [6] 琚安安, 郭红霞, 张凤祁, 郭维新, 欧阳晓平, 魏佳男, 罗尹虹, 钟向丽, 李波, 秦丽. 铁电存储器中高能质子引发的单粒子功能中断效应实验研究. 物理学报, 2018, 67(23): 237803. doi: 10.7498/aps.67.20181225
    [7] 陈锋, 郑娜, 许海波. 质子照相中基于能量损失的密度重建. 物理学报, 2018, 67(20): 206101. doi: 10.7498/aps.67.20181039
    [8] 戚俊成, 陈荣昌, 刘宾, 陈平, 杜国浩, 肖体乔. 基于迭代重建算法的X射线光栅相位CT成像. 物理学报, 2017, 66(5): 054202. doi: 10.7498/aps.66.054202
    [9] 张天奎, 于明海, 董克攻, 吴玉迟, 杨靖, 陈佳, 卢峰, 李纲, 朱斌, 谭放, 王少义, 闫永宏, 谷渝秋. 激光高能X射线成像中探测器表征与电子影响研究. 物理学报, 2017, 66(24): 245201. doi: 10.7498/aps.66.245201
    [10] 王心怡, 范全平, 魏来, 杨祖华, 张强强, 陈勇, 彭倩, 晏卓阳, 肖沙里, 曹磊峰. Fresnel波带片编码成像的高分辨重建. 物理学报, 2017, 66(5): 054203. doi: 10.7498/aps.66.054203
    [11] 强帆, 朱京平, 张云尧, 张宁, 李浩, 宗康, 曹莹瑜. 通道调制型偏振成像系统的偏振参量重建. 物理学报, 2016, 65(13): 130202. doi: 10.7498/aps.65.130202
    [12] 罗尹虹, 张凤祁, 郭红霞, 郭晓强, 赵雯, 丁李利, 王园明. 纳米静态随机存储器质子单粒子多位翻转角度相关性研究. 物理学报, 2015, 64(21): 216103. doi: 10.7498/aps.64.216103
    [13] 何小亮, 刘诚, 王继成, 王跃科, 高淑梅, 朱健强. PIE成像中周期性重建误差的研究. 物理学报, 2014, 63(3): 034208. doi: 10.7498/aps.63.034208
    [14] 聂永发, 朱海潮. 利用源强密度声辐射模态重建声场. 物理学报, 2014, 63(10): 104303. doi: 10.7498/aps.63.104303
    [15] 王林元, 刘宏奎, 李磊, 闫镔, 张瀚铭, 蔡爱龙, 陈建林, 胡国恩. 基于稀疏优化的计算机断层成像图像不完全角度重建综述. 物理学报, 2014, 63(20): 208702. doi: 10.7498/aps.63.208702
    [16] 晏骥, 郑建华, 陈黎, 林稚伟, 江少恩. X射线相衬成像技术应用于高能量密度物理条件下内爆靶丸诊断. 物理学报, 2012, 61(14): 148701. doi: 10.7498/aps.61.148701
    [17] 袁志林, 杨睿, 杨柳, 宋丽丹, 孙莉萍, 马雨虹, 王猛, 陈定康, 郭金平, 唐丽红. 基于单准直透镜的阵列准直器研究. 物理学报, 2012, 61(18): 184217. doi: 10.7498/aps.61.184217
    [18] 余波, 应阳君, 许海波. 中子半影成像的两种非线性重建方法研究. 物理学报, 2010, 59(8): 5351-5357. doi: 10.7498/aps.59.5351
    [19] 秦晓刚, 贺德衍, 王骥. 基于Geant 4的介质深层充电电场计算. 物理学报, 2009, 58(1): 684-689. doi: 10.7498/aps.58.684
    [20] 张百钢, 姚建铨, 路 洋, 纪 峰, 张铁犁, 徐德刚, 王 鹏, 徐可欣. 抽运光角度调谐准相位匹配光学参量振荡器的研究. 物理学报, 2006, 55(3): 1231-1236. doi: 10.7498/aps.55.1231
计量
  • 文章访问数:  2799
  • PDF下载量:  46
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-11-04
  • 修回日期:  2019-11-28
  • 刊出日期:  2020-02-05

高能质子照相中基于角度准直器设计的理论研究

  • 1. 中国工程物理研究院研究生院, 北京 100088
  • 2. 北京应用物理与计算数学研究所, 北京 100094
  • 通信作者: 许海波, 13641017929@163.com
    基金项目: 国家自然科学基金(批准号: 11675021)和国家自然科学青年科学基金(批准号: 11805018)资助的课题

摘要: 角度准直器在高能质子照相中有着重要作用, 既可以利用准直器提高图像对比度, 又能通过二次成像实现材料诊断及密度重建, 因此减小通过准直器后通量值的误差具有重要意义. 本文通过理论分析, 提出了一种高能质子照相中准直器设计的方法, 通过Geant4程序建立了1.6 GeV的质子成像系统, 该系统分别使用理想准直器、拉伸型准直器和利用该方法设计的准直器, 并对比通过客体后的通量分布. 结果表明, 在使用理想准直器和该方法设计的角度准直器时, 二者得到的客体的通量分布符合较好, 而使用拉伸型准直器时, 与使用理想准直器得到的结果相差较大. 因此利用理想准直器方法设计的准直器可以很好地减小通量误差.

English Abstract

    • 高能质子照相是1995年由美国洛斯·阿拉莫斯国家实验室(LANL)的科学家Chris Morris 提出的辐射成像方法[1], 与X射线成像相比, 质子照相在穿透能力、空间分辨率、密度分辨率和多次成像等方面有着明显优势[2,3]. 最近十年, 美国LANL通过洛斯·阿拉莫斯中子科学中心(LANSCE) 800 MeV的质子加速器和布鲁克·海文国家实验室(BNL) 24 GeV的质子加速器验证了质子照相技术的连续成像能力及空间分辨能力[4-6]. 俄罗斯科学院利用70 GeV的U70加速器和800 MeV的TWAC加速器设计的质子照相系统, 展开了一系列流体动力学诊断研究[7,8]. 德国重离子研究中心(GSI)利用反质子和离子研究装置(FAIR)的4.5 GeV质子加速器进行了质子照相研究[9]. 中国原子能研究院设计了100 MeV的质子成像系统[10], 兰州近代物理研究所在重离子加速器(HIRFL)设计了600 MeV/u的碳离子照相装置[11,12], 中国工程物理研究院流体物理研究所设计了11 MeV的质子成像装置[13].

      高能质子成像装置由两组四极磁透镜对和准直器组成, 如图1所示. 该系统的中心平面是傅里叶平面(准直平面), 到达该平面的质子位置仅与库仑散射角有关, 而与初始位置无关, 库仑散射角越大, 离轴越远. 这样可以在准直平面处放置角度准直器来阻挡库仑散射角较大的质子. 因此可以通过角度准直调节图像的对比度, 而且, 如果串联两组成像系统, 并使用不同的角度准直器, 则可实现二次成像进行材料诊断[14]. 由此可见, 角度准直器在质子成像中有着重要作用.

      图  1  质子成像系统示意图

      Figure 1.  Schematic diagram of proton radiography system.

      理想的角度准直器是没有厚度的, 但是在实际中要有效阻挡高能质子, 必须有一定的厚度. 目前, 实际的角度准直器一般以给定截断角时, 得到准直平面处的半径($R = {m_{12}}{\theta _{\rm{c}}}$, 其中${m_{12}}$为传输矩阵元, ${\theta _{\rm{c}}}$为截断角值), 并以此半径拉伸一定厚度, 将其设计成圆柱形或椭圆柱形的(拉伸型), 有极少数设计成针尖孔状的. 2005年, LANL利用MCNP5模拟质子成像系统, 使用了厚度为24英尺(1 in = 2.54 cm)、材料为钨的圆柱形准直器[15]. 2011年, LANL设计了800 MeV的放大型质子成像系统, 使用的准直器形状为椭圆柱形[16], 设计的24 GeV成像系统中使用了厚度为1.2 m的钨材料圆柱形准直器[6]. 2013年, 俄罗斯高能所在50 GeV的质子成像系统中使用了圆柱形准直器[17], Varentsov等[18]首次将高能质子系统用于生物学成像, 其中使用了圆柱形准直器. 2014年, 中国工程物理研究院设计的11 MeV质子成像系统获得视场40 mm, 并使用圆柱孔准直器[13]. 2015年, 兰州近代物理研究所利用Geant4模拟2.6 GeV的质子成像系统时, 使用了理想准直器[19]. 2016年, 中国原子能物理研究院设计100 MeV的质子成像系统时, 将准直器设计成针尖孔状的[10]. 2018年, Kantsyrev等[20]利用Geant4模拟247 MeV的放大型质子成像系统时, 准直器设计为椭圆柱形.

      若将准直器设计成圆柱形, 则截断角以内的部分质子也将被阻挡, 这将造成通量误差, 从而影响密度重建和材料诊断. 本文通过理论分析, 建立了高能质子照相中设计角准直器的模型, 并基于能量为1.6 GeV的质子成像系统进行了数值模拟, 得出利用该模型设计的准直器和使用理想准直器时的结果相吻合, 准直器的厚度将不影响截断角以内的质子通量.

    • 通过数值模拟可以看出截断角以内质子束团的形状在放置角度准直器的漂移空间(准直空间)内是沿漂移位置变化的椭圆, 如图2所示. 如果准直器孔径的形状与该质子束团的形状相匹配, 此时准直器的厚度将不影响截断角内质子束团的通量, 因此, 设计角度准直器即为求该截断角内质子束团在准直空间的半轴(边界线)分布. 做以下假设: 1) 恰好离开客体的质子束团近似为圆(初始束团); 2) 一定截断角以内的初始束团(目标束团)相图为理想相图.

      图  2  一定截断角以内的质子束团在通过准直空间时的形状变化 (a) z = 0 m; (b) z = 0.6 m; (c) z = 1.2 m; (d) z = 1.8 m; (e) z = 2.4 m

      Figure 2.  Shape changed of proton bunch within certain angle-cuts as it passes through the collimation space: (a) z = 0 m; (b) z = 0.6 m; (c) z = 1.2 m; (d) z = 1.8 m; (e) z = 2.4 m.

      根据质子成像系统结构, xy方向传输矩阵形式一致, 因此将相空间$\left( {x, x'} \right)$$\left( {y, y'} \right)$记为$\left( {u, u'} \right)$, 在相空间内质子初始坐标用$\left( {{u_{\rm{o}}}, u{'_{\rm{o}}}} \right)$表示, 多次库仑散射角以$\theta $表示, 根据点对点成像要求坐标应满足:

      ${u'_{\rm{o}}} = {w_u}{u_{\rm{o}}} + \theta , $

      其中${w_u}$u方向成像系统的匹配参数, 由传输矩阵决定.

      束团通过传输矩阵R后的坐标记为$\left( {u, u'} \right)$, 则有

      $\left( {\begin{aligned} u \\ {u'} \end{aligned}} \right) = \left( {\begin{array}{*{20}{c}} {{R_{11}}}&{{R_{12}}} \\ {{R_{21}}}&{{R_{22}}} \end{array}} \right)\left( {\begin{aligned} {{u_{\rm{o}}}} \\ {{u'_{\rm{o}}}} \end{aligned}} \right), \; {\text{且}}\det ({R}) = 1.$

      将(1)式代入(2)式, 并整理可得

      $ \left( {\begin{aligned} u \\ {u'} \end{aligned}} \right) = \left( {\begin{array}{*{20}{c}} {\left( {{R_{11}} + w_u{R_{12}}} \right){u_{\rm{o}}} + {R_{12}}\theta } \\ {\left( {{R_{21}} + w_u{R_{22}}} \right){u_{\rm{o}}} + {R_{22}}\theta } \end{array}} \right), $

      设初始束团尺寸为${r_{\rm{c}}}$, 截断角为${\theta _{\rm{c}}}$, 则有

      ${u_{\rm{o}}} \in \left[ { - {r_{{\rm{c}}}},\;{r_{\rm{c}}}} \right],~~\theta \in \left[ { - {\theta _{\rm{c}}},\;{\theta _{\rm{c}}}} \right].$

      求目标束团在准直空间的边界线即求u的最值, 求一阶导数:

      $ \left\{ {\begin{aligned} & {\dfrac{{\partial u}}{{\partial {u_{\rm{o}}}}} = {R_{11}} + w_u{R_{12}}}, \\ & {\dfrac{{\partial u}}{{\partial \theta }} = {R_{12}}}, \end{aligned}} \right. $

      $\left\{ {\begin{aligned} & {\dfrac{{\partial u}}{{\partial {u_{\rm{o}}}}} = 0}, \\ & {\dfrac{{\partial u}}{{\partial \theta }} = 0} , \end{aligned}} \right. $

      根据(5)和(6)式可得

      $\left\{ {\begin{aligned} & {{R_{11}} = 0} \\ & {{R_{12}} = 0} \end{aligned}} \right.\quad \Rightarrow \quad \det ({R}) = 0.$

      由(2)式可知, 该结果不符合传输矩阵条件, 因而驻点不存在, 所以u的最值在边界上求得

      $\begin{split} u_{\rm m} =& \big\{u(r_{\rm c}, \theta _{\rm c}),u(- r_{\rm c}, \theta _{\rm c}),u(r_{\rm c}, - \theta _{\rm c}), \\ \, &\; u(- r_{\rm c}, - \theta _{\rm c}) \big\}_{\max (\min )}. \end{split}$

      通过上述分析, 如果已知成像系统传输矩阵、客体或视场半径和截断角值即可得到准直器参数. xy方向传输矩阵的具体形式如下:

      $\begin{split} {{{R}}_x} =& \left( {\begin{array}{*{20}{c}} 1&z \\ 0&1 \end{array}} \right)\left( {\begin{array}{*{20}{c}} {\cosh kl}&{\dfrac{1}{k}\sinh kl} \\ {k\sinh kl}&{\cosh kl} \end{array}} \right)\left( {\begin{array}{*{20}{c}} 1&t \\ 0&1 \end{array}} \right) \\ & \times\left( {\begin{array}{*{20}{c}} {\cos kl}&{\dfrac{1}{k}\sin kl} \\ { - k\sin kl}&{\cos kl} \end{array}} \right)\left( {\begin{array}{*{20}{c}} 1&s \\ 0&1 \end{array}} \right), \\ {{{R}}_y} =& \left( {\begin{array}{*{20}{c}} 1&z \\ 0&1 \end{array}} \right)\left( {\begin{array}{*{20}{c}} {\cos kl}&{\dfrac{1}{k}\sin kl} \\ { - k\sin kl}&{\cos kl} \end{array}} \right)\left( {\begin{array}{*{20}{c}} 1&t \\ 0&1 \end{array}} \right)\\ & \times \left( {\begin{array}{*{20}{c}} {\cosh kl}&{\dfrac{1}{k}\sinh kl} \\ {k\sinh kl}&{\cosh kl} \end{array}} \right)\left( {\begin{array}{*{20}{c}} 1&s \\ 0&1 \end{array}} \right), \\[-20pt] \end{split} $

      其中z是准直空间的漂移距离, l是磁透镜厚度, t是磁透镜对之间的内漂移距离, s是客体平面到第一个磁透镜的外漂移距离, 如图3所示. k = $ \sqrt {{G}/({{B\rho }})} $, G是磁场梯度, $B\rho $为粒子磁刚度. 通过(8)和(9)式可以得到目标束团在准直空间的边界线, 以1.6 GeV质子成像系统进行数值计算, 表1列出该系统的具体参数.

      s/ml/mG/T·m–1t/m
      1.20.88.090.5

      表 1  1.6 GeV质子成像系统参数

      Table 1.  Parameters of the proton radiography system of 1.6 GeV.

      图  3  质子成像系统参数示意图

      Figure 3.  Diagram of parameters of proton radiography system.

    • 图4是截断角为2 mrad、客体尺寸为5 cm时的目标束团边界线(半轴分布). 中间位置即为准直平面, 此时xy方向的半轴相等, 即束团形状在x-y平面为圆形, 而其他位置均不相等, 此时束团形状为椭圆形, 因此可以看出目标束团在准直空间形成的形状是一系列变化的椭圆. 选择一段包含准直平面在内的边界线可以作为准直器的孔径参数, 本文以z = 20 cm (前端)到z = 120 cm (后端)之间的部分作为准直器, 即准直器厚度为1 m.

      图  4  截断角为2 mrad、客体尺寸为5 cm时的目标束团边界线

      Figure 4.  Boundary lines of the target bunch when the angle-cut is 2 mrad and the object size is 5 cm.

      图5给出了客体尺寸为5 cm、截断角变化时的孔径值曲线, 可以看出孔径值与截断角大小成正比. 其中, 图5(a)是前端口孔径值的变化曲线, 可以得出xy方向的半轴不相等, 所以端口形状为椭圆形; 图5(b)是后端口孔径值的变化曲线, 可以得出端口形状为圆形, 此时准直器孔径的整体结构是椭圆台状的. 图6是客体半径为5 cm, 截断角为2 mrad时的准直器形状.

      图  5  端口处的孔径值随截断角的变化 (a) $z = 20\;{\rm{cm}}$; (b) $z = 120\;{\rm{cm}}$

      Figure 5.  Aperture size varies with the angle-cut at the ports: (a) $z = 20\;{\rm{cm}}$; (b) $z = 120\;{\rm{cm}}$.

      图  6  准直器孔径的形状 (a) x-y平面; (b) y-z平面

      Figure 6.  Shape of aperture of the collimator: (a) x-y plane; (b) y-z plane.

    • 利用Geant4[21,22]软件, 设计了1.6 GeV的质子成像系统. 该系统中探测器使用薄层介质板, 质子束流为单能质子. 模拟中分别使用理想型、拉伸型和利用本文模型设计的准直器称之为设计型, 并对比通过客体后的通量分布. 客体结构如图7所示, 图7(a)是厚度为1 mm、半径为5 cm的铜板. 图7(b)由一组同心球层组成, 中心区是空区, 空球半径为1 cm, 第二层和第三层分别是锂和钠, 球半径分别是4 cm和5 cm. 准直器分别以客体尺寸为5 cm, 截断角为2和3.5 mrad设计, 表2列出三类准直器的孔径参数.

      截断角/mrad准直器类型前端/cm 后端/cm厚度/m外半径/m材料
      xy xy
      2理想型 0.6110–73Al
      拉伸型0.61 0.611 & 0.43W
      设计型2.042.46 0.6113W
      3.5理想型 1.0710–73Al
      拉伸型1.07 1.071 & 0.43W
      设计型2.343.08 1.0713W

      表 2  准直器的孔径参数

      Table 2.  Aperture parameters of the collimator.

      图  7  客体示意图 (a)铜板; (b)同心球体

      Figure 7.  Diagram of the object: (a) The round copper plate; (b) the concentric spheres.

      使用不同的准直器时, 质子束流通过客体后的通量分布如图7图8所示. 图8(a)图8(b)分别是截断角为2和3.5 mrad时质子束流通过铜板的结果. 图9(a)图9(b)分别是截断角为2和3.5 mrad时质子束流通过同心球客体的结果. 在使用该模型设计的准直器时, xy方向的通量均与使用理想准直器时的结果相吻合, 使用拉伸型准直器且厚度为1 m时, 通量分布与使用理想准直器时的结果相差较大, 越靠近边缘, 相差越大; 当厚度为0.4 m时, 整体误差比1 m厚时的较小; 当截断角为3.5 mrad, 位置范围在(–10, 10) mm时, 0.4 m厚的拉伸型准直器造成的误差较小. 因此拉伸型准直器对于客体尺寸较小、截断角较大的情况下造成的通量误差较小, 反之误差较大. 准直器足够厚时, 利用该模型设计的准直器对客体的通量分布影响较小, 这有利于减小密度重建的误差, 提高材料诊断的准确性.

      图  8  通过铜板的通量分布 (a)截断角为2 mrad; (b)截断角为3.5 mrad

      Figure 8.  Flux distribution after passing the round copper plate: (a) Angle-cut of 2 mrad; (b) angle-cut of 3.5 mrad.

      图  9  通过同心球的通量分布 (a)截断角为2 mrad; (b)截断角为3.5 mrad

      Figure 9.  Flux distribution after passing the concentric spheres: (a) Angle-cut of 2 mrad; (b) angle-cut of 3.5 mrad.

    • 通过建立模型和理论分析, 给出了高能质子照相中角度准直器的设计方法. 利用该方法设计的角度准直器, 孔径结构是椭圆台状的, 它的厚度将不影响目标束团的通量. 基于Geant4程序模拟了1.6 GeV的质子成像系统, 系统中分别使用了理想型、拉伸型和该方法设计的准直器, 通过对比不同准直器下客体的通量分布, 得出在使用拉伸型时, 当客体尺寸小、截断角较大时, 通量差异较小, 反之误差较大. 在使用设计型准直器时, 与理想条件下的结果相吻合, 表明使用该方法设计的角度准直器时, 厚度将不影响目标束团的通量. 高能质子照相的能量一般在GeV量级以上, 需要更厚的准直器才能有效阻挡质子, 使用该方法设计的准直器有更大的优势. 这对于减小密度重建的误差和提高材料诊断的准确性有重要意义.

参考文献 (22)

目录

    /

    返回文章
    返回