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

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

1) (中国工程物理研究院研究生院, 北京 100088)

2) (北京应用物理与计算数学研究所, 北京 100094)

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

1 引 言

高能质子照相是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].由此可见, 角度准直器在质子成像中有着重要作用.

理想的角度准直器是没有厚度的, 但是在实际中要有效阻挡高能质子, 必须有一定的厚度.目前,实际的角度准直器一般以给定截断角时, 得到准直平面处的半径( R =m12θc, 其中 m12为传输矩阵元, θ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的质子成像系统进行了数值模拟, 得出利用该模型设计的准直器和使用理想准直器时的结果相吻合, 准直器的厚度将不影响截断角以内的质子通量.

图1 质子成像系统示意图Fig.1.Schematic diagram of proton radiography system.

2 角度准直器模型

通过数值模拟可以看出截断角以内质子束团的形状在放置角度准直器的漂移空间(准直空间)内是沿漂移位置变化的椭圆, 如图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 mFig.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.

根据质子成像系统结构, x和y方向传输矩阵形式一致, 因此将相空间 ( x,x′) 和 ( y,y′) 记为 ( u,u′) ,在相空间内质子初始坐标用 ( uo,) 表示, 多次库仑散射角以 θ 表示, 根据点对点成像要求坐标应满足：

其中 wu是u方向成像系统的匹配参数, 由传输矩阵决定.

束团通过传输矩阵R后的坐标记为 ( u,u′) ,则有

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

设初始束团尺寸为 rc, 截断角为 θc, 则有

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

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

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

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

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

图3 质子成像系统参数示意图Fig.3.Diagram of parameters of proton radiography sys−tem.

表1 1.6 GeV质子成像系统参数Table 1.Parameters of the proton radiography system of 1.6 GeV.

3 数值计算与模拟

3.1 数值计算

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

图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 =20cm ;(b)z=120cmFig.5.Aperture size varies with the angle−cut at the ports：(a) z =20cm ; (b) z =120cm .

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

图6 准直器孔径的形状 (a) x-y平面; (b) y-z平面Fig.6.Shape of aperture of the collimator：(a) x-y plane;(b) y-z plane.

3.2 蒙特卡罗模拟

利用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列出三类准直器的孔径参数.

图7 客体示意图 (a)铜板; (b)同心球体Fig.7.Diagram of the object：(a) The round copper plate; (b) the concentric spheres.

表2 准直器的孔径参数Table 2.Aperture parameters of the collimator.

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

图8 通过铜板的通量分布 (a)截断角为2 mrad; (b)截断角为3.5 mradFig.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 mradFig.9.Flux distribution after passing the concentric spheres：(a) Angle−cut of 2 mrad; (b) angle−cut of 3.5 mrad.

4 结 论

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