基于透镜阵列多元聚焦的连续晶体三维位置灵敏探测器设计和仿真.pdf

第58 卷第1期2024年1月原子能科学技术Atomic Energy Science and TechnologyVol.58,No.1Jan.2024基于透镜阵列多元聚焦的连续晶体三维位置灵敏探测器设计和仿真李婷1,2.3，庄凯1-2.3，李道武1-3，梁秀佐1-3,刘彦韬1-3,张译文1-3，孔令钦1-2 3，章志明1-2.34，帅磊1.34，魏 1.2.34（1.中国科学院高能物理研究所，北京市射线成像技术与装备工程技术研究中心，北京10 0 0 49；2.中国科学院大学核科学与技术学院，北京10 0 0 49；3.国家原子能机构核技术（核探测与核成像）研发中心，北京10 0 0 49；4.济南中科核技术研究院，山东济南2 50 131)摘要：连续晶体闪烁体探测器可以实现射线作用点的三维位置定位。本文提出一种新的探测器结构，引入透镜阵列实现作用点定位。探测器采集多元聚焦图像，从而改变对连续晶体闪烁光分布探测的需求。文中建立连续晶体、透镜阵列和光子计数器阵列的模型进行模拟研究，建立多元聚焦的光路反演重建算法定位射线作用点的三维位置。通过重建结果分析评估基于透镜阵列的连续晶体探测器性能，这种结构具有较好的定位能力和位置分辨。对于尺寸为48 mmX48mmX45mm的硅酸钇(LYSO)闪烁晶体，实现了cy平面位置分辨优于1.54mm，之方向位置分辨优于3.13mm。基于连续晶体和透镜阵列结构的探测器还可对多作用点的情形进行重建，分析三组选定位置的重建结果，该结构对双作用点定位具有可行性。关键词：连续晶体探测器；透镜阵列；作用点定位；位置灵敏探测器中图分类号：TL99;0572.212doi:10.7538/yzk.2023.youxian.0042Design and Simulation of Three-dimensional Position SensitiveDetector with Monolithic Scintillation Crystal Based on Lens ArrayLI Tingl,ZHUANG Kail-23,LI Daowu,LIANG Xiuzuol*,LIU Yantaol*,ZHANG Yiwen,KONG Lingqinl-2-3,ZHANG Zhiming.-2.3.,(1.Beijing Engineering Research Center of Radiographic Techniques and Equipment,Institute of High Energy Physics,Chinese Academy of Sciences,Beijing 100049,China;2.School of Nuclear Science and Technology,University of Chinese Academy of Sciences,Beijing 100049,China;3.CAEA Center of Ercellence on Nuclear TechnologyApplications for Nuclear Detection and Imaging,Beijing 100049,China;4.Jinan Laboratory of Applied Nuclear Science,Jinan 250131,China)Abstract:The monolithic scintillation crystal detector has the ability to position the文献标志码：Afor Multiple Focus PointsSHUAI Leil.3.*,WEI Long.文章编号：10 0 0-6 9 31(2 0 2 4)0 1-0 2 39-0 9收稿日期：2 0 2 3-0 2-0 6；修回日期：2 0 2 3-0 2-2 4基金项目：国家自然科学基金（117 0 52 2 4，118 0 52 13）*通信作者：帅磊240three-dimension interaction location of the gamma ray in the detector depending on accu-rate acquisition of scintillation light distribution.Depth estimation of light field can besolved by a lens array in light field technology.A new detector model was designed inthe paper which introduced a lens array for positioning the interaction in the monolithicdetector.The scintillation light is focused by each lens unit in the lens array.Photons insome discrete pixels of the photon counter array where the scintillation light is focus onare detected to form an image of multiple focus points instead of the scintillation lightdistribution detection based on continuous crystal.The position of the scintillation lightsource in the crystal can be reconstructed from the focus image through optical pathinversion.The detector model consisting of a monolithic crystal,a lens array and amulti-pixel photon counter was established in the paper.The monolithic scintillator wasa lutetium-yttrium silicate（L Y SO)s c in t illa t o r w it h d im e n s io n s 48 m m X48 m m X45 mm.The lens was hemispherical with a radius of 15 mm supposed to be uniformlydistributed in 3X3 array and the material was as same as the crystal.The crystal andthe lens array were considered as a whole in the simulation.The multi-pixel photoncounter was set at 44.4 mm from the crystal surface.The pixel size of the photon counterarray was 3 mm and the quantum efficiency was 30%.The axis of symmetry of thethree parts was same.The reconstruction algorithm was based on optical path inversionfrom multiple focus points for locating the 3D position of gamma-ray interaction.Theintersection points of the inverted optical paths and the set plane with different depthsin the crystal was calculated in the reconstruction process.When the intersection pointsconverged towards the same point in the zo plane,the depth of interaction was zo and theintersection point in the zo plane was the interaction point.The reconstruction resultswere analyzed to evaluate the performance of the monolithic scintillation crystal detectorbased on the lens array.The x/y resolution is better than 1.54 mm and the z resolutionis better than 3.13 mm.The two interaction points when coincidence occurs are recon-struced in the condition that the energy is the same.The detector is feasible for the posi-tioning of two interaction points by the reconstruction results of three groups of selectedpositions.The monolithic scintillation detector based on a lens array has both good posi-tioning ability and position resolution.The detector is suitable for the position recon-struction of multiple interaction points.Key words:monolithic crystal detector;lens array;interaction position;position sensi-tive detector射线成像技术广泛应用于生物医学、天文学、辐射环境监测、国土安全以及公共安全等领域1-5。位置灵敏探测器作为射线探测的关键设备，是决定射线成像系统性能指标的重要因素。基于闪烁体的位置灵敏探测器阻止本领高，位置分辨、能量分辨和时间分辨优异，还具有操作简易、稳定性好等优点，使得闪烁体探测器在射线成像应用中起着重要的作用。闪烁体探测器利用闪烁晶体将射线能量转化为闪烁原子能科学技术第58 卷光，再由光子计数器转化为电信号，电信号内包含射线和晶体作用的位置信息和能量信息。目前大多数高分辨的探测器采用晶体阵列耦合光子计数器阵列来实现。当以缩小晶体单元尺寸的方式进一步提高晶体阵列探测器的空间分辨时，受到了机械加工的限制，难以继续提升分辨性能，并且晶体间隔的填充降低了探测器的探测效率，这些是晶体阵列探测器存在的问题。第1期李婷等：基于透镜阵列多元聚焦的连续晶体三维位置灵敏探测器设计和仿真基于连续晶体的探测器可以避免晶体阵列带来的问题，具有设计结构简单、能量分辨率高、探测器效率高以及低造价等优点。随着新的闪烁晶体、硅光电倍增管（SiPM)和数字电子学技术的发展，连续晶体探测器有了新的可行性方案。连续晶体探测器通过所探测到每个事例的闪烁光分布来计算射线在晶体内的作用位置。因此当采用连续晶体结构的探测器时，高分辨的射线位置定位算法是不可缺少的。目前已有多种高分辨的定位算法被研究和应用开发，如神经网络定位算法6-7 以及K-NN定位算法8-10，这些算法都依赖于精确获取闪烁光输出分布。光学成像技术中，在普通成像系统主透镜的一次像面处插入透镜阵列，每个透镜单元及其后对应的传感器区域记录光线在场景中相同241用中，通过透镜阵列中的透镜单元在不同视场对同一个物体成像，相当于将三维物体的深度信息转化为二维透镜阵列所成像的角度信息14-15，每个透镜单元所成的像都包含着深度信息。经过透镜中心的光线传输方向不变，以此简化光线的传播，可将透镜阵列简化为小孔阵列，光场成像深度估计原理如图1所示。不同深度的光源A、B的光线在经过透镜阵列后，在像平面的单元图像中位置不同，反映出A、B光源在物体空间中的位置差异。采集到的光线物体A3部分在不同视角下所成像的集合，因此采用二维透镜阵列能得到同时包含位置和传播方向在内的四维光场数据，解决光场深度估计和光源定位问题。本文引人透镜阵列耦合连续晶体光输出面，提出一种新型的连续晶体探测器结构，实现对射线作用点的三维位置探测，并有望实现多个作用点的定位重建。在闪烁光传输过程中，透镜阵列中的各个透镜单元对其接收到的闪烁光聚焦，在光子计数器阵列平面上形成含有多个聚焦点的图像，每个聚焦点包含有闪烁光源位置和传播方向的信息，由聚焦点图像通过光路反演即可重建闪烁光源在晶体内的位置。这种结构的探测器不需要对连续晶体闪烁光分布进行精确探测，通过探测聚焦位置即可实现作用点定位重建。本文通过建立模型进行蒙特卡罗模拟研究，提出一种连续晶体对射线作用点的三维定位的解决方法，并通过重建算法结果分析来评估这一技术的可行性，并对晶体内两个作用点的情形进行研究。1重建算法原理透镜阵列应用在光场成像领域中，通常单元透镜直径非常小，因而又可称为微透镜阵列。微透镜阵列作为一种基础的阵列光学元器件，对光信息有很好的聚焦、准直、交换、多重成像和综合成像的能力，已经应用于光场成像3D技术、深度估计等研究中1-13。在深度估计应1单元图像位于不同深度的物点小孔阵列图1透镜阵列用于光场成像深度估计原理示意图Fig.1Schematic diagram of depth estimationin light field imaging based on lens array将透镜阵列应用于连续晶体探测器，射线在晶体内作用转化的闪烁光传输到晶体光输出面，再经由透镜单元传输到透镜与空气的界面上发生折射而改变闪烁光传输方向，闪烁光在空气中向靠近光轴的方向继续传输。经过透镜阵列时闪烁光的传播示意图如图2 所示，蓝色光线表示晶体沉积射线能量后转化的闪烁光，图2 经过透镜阵列的闪烁光传播示意图Fig.2Schematic diagram of propagationof scintillation light through lens array242红色光线表示从透镜阵列出射的闪烁光，出射闪烁光继续传输直到被光电器件探测，透镜对闪烁光传输方向的调制作用使得其接收到的闪烁光到达光子计数器阵列的位置更加集中。光子计数器阵列的光电器件采集图像并输出电信号，被后续电子学转化为数字信号形成数字图像。当调节光子计数器阵列与透镜阵列的间距，使得闪烁光被探测的位置集中形成的光斑小于光电器件像素大小时，采集图像可视为一组聚焦点的像素化图像，一次事例形成的图像中含有多个成像亮点。这些光斑或亮点仍包含有射线作用点的位置和沉积能量信息，成像亮点在图像中的位置与作用点的位置有关，成像亮点的光子数与射线的沉积能量、透镜对射线作用点的立体角有关。采集到多个聚焦亮点图像进行射线作用点位置重建时，主要基于闪烁光路的光线对应关系和光的可逆性原理。射线能量不大时射线作用点可视为一个点光源，因而在闪烁光的传播路径上，闪烁光源、任一透镜中心和该透镜对应的图像亮点在一条光线上。将光路反演，图像亮点和透镜中心的连线反推到晶体内，反推的光线必定经过闪烁光源。算法中通过预设晶体内不同深度的平面，示意图如图3所示，计算反推光线与平面的交点S；，当各个透镜反推光线的交点在晶体的某一深度平面内汇聚为一点时，该点即为射线与晶体的作用点，此时的平面深度即为作用点的值，再计算当前深度平面中反推光线与平面的交点即为作用点的、y值。由此，基于透镜阵列的连续晶体探测器通图3多点聚焦的光路反演重建算法示意图Fig.3 Schematic diagram of reconstructionalgorithm based on optical path inversionfrommultiplefocus points原子能科学技术第58 卷过多点聚焦和光路反演，实现了晶体探测器对射线作用点的位置重建。2探测器模拟研究及算法处理流程2.1探测器蒙特卡罗模型及模拟参数设置使用FLUKA蒙特卡罗模拟软件，建立连续晶体及透镜阵列模型以计算光子计数器阵列的采集图像。闪烁体选用尺寸为48 mmX48mmX45mm的连续晶体，材料为LYSO,折射率为1.8 2。晶体表面除出光面耦合透镜阵列端面外，均涂有黑色吸收层以防止多次反射造成多重影像。透镜采用半球型，半径为15mm，选取3X3阵列，球心均匀分布，材料与晶体相同，在模型中可将透镜阵列与连续晶体视为一个整体，对闪烁光没有介质界面。光子计数器阵列采用 SiPM阵列,像素大小为3 mm，量子效率为30%，设置在距离晶体出光面44.4mm的位置。连续晶体、透镜阵列、SiPM阵列的对称轴在一条直线上，垂直于晶体出光面。采集图像中每个像素的数值为当前像素探测到的闪烁光子计数。设置射线能量为511keV，在晶体内仅发生一次作用。探测器的作用过程为：当射线与闪烁晶体发生作用时，沉积部分或全部能量产生闪烁光。闪烁光在连续晶体和透镜阵列的一体化结构中传输，当传输到透镜与空气的界面时发生折射，并沿折射方向继续传播至光子计数器阵列被探测成像，成像结果如图4所示。为研究晶体内某选定作用点位置的重建结果统计特性，将事例转换为质子事例，质子的出射位置即为事例的作用点位置。以晶体入射面中心点（0，0，一45）（单位mm，下同）开始设置作用点位置，在E一45，一2 0 、E0,8、y E,8（单位mm，下同)晶体空间内，间隔1mm设置作用点。在E12,2 4，y E,2 4 晶体空间内，间隔4mm设置作用点，方向间隔仍为1mm。本文将设置在晶体边沿1的作用点、y 坐标取为2 3mm,而非2 4mm。2.2多元聚焦的光路反演重建算法流程算法处理流程分为以下3个步骤。1）确定亮点像素位置采集图像以亮点为中心形成9 个分立的区域，区域内除亮点外仍有少量的光子。因采用第1期李婷等：基于透镜阵列多元聚焦的连续晶体三维位置灵敏探测器设计和仿真45Fa4035F30F25F20H15F10F55101520 25300.35.4045方向像素a中心对称轴；b偏离中心对称轴图4作用点在晶体中心对称轴和偏离中心对称轴时成像结果示意图Fig.4 Schematic diagram of multiple focus point image when interactionposition is on or departure from axis of crystal symmetry的半球型透镜对“点”物的成像存在像差，部分闪烁光子分布在焦点区域外，亮度低于焦点。首先对采集图像进行区域划分，区域划分采用K-Means聚类算法。划分区域后,在每个区域内采用峰值法确定亮点像素位置，记录为（，y,），i=19。2）计算重建位置在2.1节所述的模型下，闪烁光源、透镜中心、图像亮点所在的光线为一条直线，在光路反演时可直接采用几何比例关系，计算反演光路光线与晶体内深度为zo的平面交点S，的位置为(i,y)z o，i=l9。计算所有交点之间的距离，以式(1)最小距离平方和来确定作用点的深度Smin值，如图5所示Smin可估计为一40 mm。y值由深度为smin的平面上各交点的平均值(,)zsmin确定。600500400300F2001000-45图5重建时计算最小距离平方和参数结果示意图Fig.5 Schematic diagram of minimum distancesquare sum parameter in reconstruction24345b12401035F830F256520415F210F5F0151015202530354045方向像素f(z)=Z Il S;(z)-S,(z)I 23）判选重建结果固定一个作用点位置进行多次事例重建了后,对重建位置结果的、之3个方向分别做拟合分析，、y方向重建结果为高斯分布，超过3c的重建结果认定为异常值，在后续分析中剔除。在算法程序中，可用式（1)的f(zsmin）作为参数进行判选。在单个作用点时，f（z s m i n）的大小可以表征成像亮点位置确定的准确程度，从而可反映重建结果偏离真实值的程度。文中数据分析时设定f（z s m i n）值超过40 时，重建结果偏离真实值，不进行统计分析。3重建结果分析3.1单个作用点时探测器定位性能以2.2 节重建流程计算得到2.1节模型中作用点位置设置下10 0 0 次事例的重建结果，重建结果绘制为散点图，如图6 所示，以作用点位置为（4，4，一2 0）时的重建结果示意，画图时将取绝对值（下同）。分析重建结果的各方向统计特性，如图7 所示。以下对选取位置的重建结果分析时发现方向的结果不是标准的高斯分布，在下文的计算中仍以高斯分布进行计算，对不同选取位置的重建-40-3530-25-202平面深度/mm20151050结果进行对比分析。选取表1中6 个作用点位置，分别计算得出重建结果的、3个方向平均值和标准差来评估当前结构探测器的定位能力。结果表明，(1)181？244原子能科学技术第58 卷25a20151010y方向重建位置/mm50-5-5x方向重建位置/mm21b2019OO800105y方向重建位置/mm422x方向重建位置/mm64a一一未进行判选；b一一经过判选图6 未进行判选和经过判选的重建结果散点图Fig.6 Scatter diagram of reconstruction result before or after data filtering120a1008060402002.53.03.54.04.55.05.5x方向像素a图7 重建结果在、y、之方向分布情况Fig.7 Distribution of reconstruction result in,y and z directions透镜阵列的引入，能解决连续晶体的作用点定位问题，定位精度较好。表1固定在不同作用点位置的重建结果Table 1Reconstruction resultof interaction at different positions设置值/mm重建值/mmy重建值/mm(0,0,-25)0.010.35(0,0,-40)0.040.29(0,4,-25)0.050.34(0,4,-40)0.060.31(4,4,25)4.110.29(4,4,-40)4.290.26对中心透镜所对的E一45，一2 0、Eo,8 、y E,8晶体空间区域的作用点重建结果的、y、3个方向误差进行统计分析，结果列于表2。结果表明，基于连续晶体和透镜阵列结构的探测器在、方向的定位准确度较高，在方向的定位准确度较弱。140b12010080事604020H2.53.03.54.04.55.05.5y方向像素方向;by方向;c方向表2中心区域重建结果定位误差情况Table 2Position error of reconstruction resultof interaction in central region重建类别误差/mm最小值-0.80中位值0.07重建值/mm最大值-0.010.3225.220.640.040.31-38.060.504.070.31-24.660.574.070.28-39.410.724.070.3124.410.554.260.26-40.970.94100r80604020F04171.00选取y=0 mm平面,以=0、2、4、6、8 m m时重建结果的、方向平均值做分布图，如图8 所示。重建结果的位置是不均匀、非线性60(0,0)(0,2)(0,4)(0,6)(0,8)uu/z40204一1图：y=0mm平面固定作用点重建结果Fig.8Reconstruction result of interactionin y=0 mm planeC182方向像素重建误差/mm0.840.071.0213x/mm19520重建误差/mm-2.060.612.857921第1期李婷等：基于透镜阵列多元聚焦的连续晶体三维位置灵敏探测器设计和仿真的，在虚线框选的位置分别出现了之方向的重建位置重叠和y平面内的重建位置重叠。这种现象的出现是由于采集图像的像素化和重建算法中对图像像素点的定位网格化带来的。3.2单个作用点时探测器位置分辨性能表3列出在不同区域的重建结果，分别计Table 3Position error and resolution of reconstruction result of interactions in different regions文半高宽J半高宽半高宽半高宽半高宽半高宽作用点误差区域范围/mm中心透镜-0.801.000-0.841.02-2.062.850.241.120边线透镜-0.730.299-0.831.04-0.043.21边角透镜0.600.600-0.680.59245算作用点定位在、3个方向的误差及半高宽，3个区域分别对应中心、边线、边角位置上的透镜。其中，中心透镜对应坐标值范围为E0,8 、E,8 ;边线透镜对应坐标值范围为E0,8 、y E12,2 4；边角透镜对应坐标值范围为E12,24、y E,24。表3不同区域的作用点定位误差及位置分辨误差误差范围/mm范围/mm0.531.390.401.540.813.130.763.700.551.200.571.220.972.68可重建范围/范围/mmmm0.221.12范围/中位值/中位值/中位值/mmmm0.542.960.680.780.77事例数mmmm0.691.460.711.740.822.012957951042910183中心透镜对应的区域内误差及分辨性能优于边线透镜，边角透镜对应的区域内作用点重建结果误差及分辨最差。可重建事例数为数据判选后的事例数，数值高可表征发生在该作用点的事例容易被重建，可见晶体内可重建范围集中在晶体人射面、靠近晶体中心轴。以各方向的半高宽最大值计，y平面位置分辨优于1.54mm，z 方向位置分辨优于3.13mm。以各方向的半高宽中位值来看，均匀分布在晶体横向截面上的作用点，一半数量的作用点位置分辨优于0.8 2 mm；均匀分布在方向上的作用点，一半数量的作用点位置分辨优于2.0 1mm。探测器的位置分辨性能良好。3.3两个作用点时探测器定位性能晶体内发生康普顿散射事例或偶然符合事例时，在两个以上作用点位置产生闪烁光，这些闪烁光在一次采集时间窗内不能被区分，因此光子计数器阵列采集到的图像是闪烁光的叠加图像。连续晶体通过光分布重建作用点定位的重建结果仅能给出一个作用点位置，且可能不能定位到任何一个作用点的真实位置上。透镜阵列的引入，将光子计数器阵列对光分布的精确连续探测转换成亮点的离散探测，不仅可以实现单个作用点的位置重建，还有望实现多作用点的位置重建。本文研究发生偶然符合事例时双作用点的情况，射线能量为单一能量。2.1节模型中固定作用点的图像为一次光电事例的图像，当晶体内发生偶然符合的两次光电事例，成像结果数据可由单个作用点的图像进行加和来获得。选取3组事例分析双作用点重建的可行性，位置组合分别为（0,1，一40）十（0,1，一30）（0,1，一40)+(0,4,30)、(8,8，一40)+(8,8,一30),重复10 0 次事例。表4列出3组事例重建结果的统计情况。(0,1,一40)+(0,1,一3 0)、（8,一8，一40)+（8,8，一30）两组事例的重建结果平均值及标准差与单个作用点的情况一致，（0，1，一40）十(0,4，一30)这组事例的位置重建结果较差。通过分析y重建值的分布发现，按照2.2 节重建流程的判选条件仍有偏离较大的结果未被剔除，这部分偏离结果是由于亮点位置定位不准确带来的。重复10 0 0 次事例后分析y重建值分布，超过3c的重建结果认定为偏离结果，剔除这些偏离结果后这组事例的重建结果平均值及标准差与单个作用点的情况一致，但可重建事例率仅为7%，即当前组合位置的准确重建程度较低。图 9为 3组事例的重建结果散点图，其中图9 b中已剔除偏差较大的结果。以上重建是基于两个作用点沉积的能量接近，当能量差异较大时重建准确度会下降。考虑应用到康普顿散射时，如设定射线人射方向为垂直晶体，（0，1，一40）十（0，4，一30）这组事例的散射角约为16.7（一8，一8，一40)十246位置组合(0,1,-40)+(0,1,-30)(0,1,-40)+(0,4,-30)(8,-8,40)+(8,8,-30)原子能科学技术第58 卷表4连续晶体内两个作用点重建结果统计情况Table 4Reconstruction result of two interactions in monolithic crystala1重建值/mmy1重建值/mm21重建值/mm2重建值/mmy2重建值/mm2重建值/mm0.020.301.310.120.270.331.260.320.270.331.070.28-7.630.25-7.600.23-39.670.3639.430.5939.240.69-38.190.290.010.220.010.350.050.247.630.161.010.052.880.863.520.137.640.16-29.900.1628.411.0929.190.18-28.890.29403530252.0y方向重建位置/mm1.51.00.3-1x方向重建位置/mma8800088040b353025y方向重建位置/mm4240C35080%30256y方向重建位置/mm00-1x方向重建位置/mm100-10-10 x方向重建位置/mma(0,1,40)+(0,1,30);b(0,1,40)+(0,4,30);c(8,8,40)+(8,8,30)图9连续晶体内两个作用点重建结果散点图Fig.9 Scatter diagram of reconstruction result of two interactions in monolithic crystal（8,8，一30)这组事例的散射角约为6 6.2。康普顿散射角较大的事例重建准确度较高，且康普顿散射角较大的事例散射点能量和吸收点能量相近，利于双作用点的定位。本文重点研究探测器的定位性能，后续研究如能进一步解决能量探测问题，透镜阵列结构的连续晶体探测器便可为康普顿射线成像技术提供一种新型的探测器构型。4结论透镜阵列应用于连续晶体探测器时，射线沉积能量转换的闪烁光传输到透镜与空气的界面上发生折射，改变其传输方向，在透镜单元接收范围内的闪烁光出射到空气后均向靠近该透镜单元光轴的方向继续传输，直至被光子计数器阵列采集，形成光斑阵列的离散图像。光子计数器阵列采集位置在透镜焦距附近时，采集的光斑图像会形成聚焦点样式的图像。透镜阵列的引入，改变了连续晶体探测器中对出射闪烁光连续分布的精确探测，转换为对离散的闪烁光聚焦点阵列的探测。通过多元聚焦的光路反演重建能解决连续晶体探测器对射线作用点的定位问题，实现射线作用点的三维位置探测。文中模拟采用像素大小为3mm的光子计数器阵列采集聚焦点图像，y平面重建位置误差在1.04mm内，方向重建位置误差在3.7 0 mm内,y平面位置分辨优于1.54mm，方向位置分辨优于3.13mm，探测器性能良好。本文采用数据为连续晶体和透镜阵列一体化结构的较理想情况下的结果，实际探测时透镜阵列和连续晶体之间存在的界面将影响聚焦点图像的亮暗情况，从而影响“亮点”的定位，后续研究中需要改进“亮点”像素位置算法。透镜阵列的引人还能实现对连续晶体内双作用点的同时定位，有望应用于康普顿散射探测，但同样需要改进“亮点”像素位置算法并解决能量探测的问题。参考文献：1 JONES T,TOWNSEND D.History and futuretechnical innovation in positron emission tomo-graphyJJ.Journal of Medical Imaging,2017,4(1):011013.2MCGREGOR D S.Materials for gamma-rayspectrometers:Inorganic scintillatorsJ.Annu-al Review of Materials Research,2018,48(1):245-277.3NIKL M,YOSHIKAWA A.Recent R&Dtrends in inorganic single-crystal scintillatormaterials for radiation detectionJ.Advanced第1期李婷等：基于透镜阵列多元聚焦的连续晶体三维位置灵敏探测器设计和仿真Optical Materials,2015,3(2):463-481.4 SSIVELS C B,MCINTYRE J I,BOWYER T W,et al.A review of the developments of radiox-enon detectors for nuclear explosion monitoringJ.Journal of Radioanalytical and NuclearChemistry,2017,314(2):829-841.5YANAGIDA T.Inorganic scintillating materialsand scintillation detectorsJ.Proceedings of theJapan Academy Series B-Physical and BiologicalSciences,2018,94(2):75-97.6周宇，王永纲基于连续晶体PET探测器模块电子学设计JI核电子学与探测技术，2 0 14，34(3):307-311.ZHOU Yu,WANG Yonggang.The readoutelectronics for the continuous crystal PET detec-tor moduleJ.Nuclear Electronics&.DetectionTechnology,2014,34(3):307-311(in Chinese).7DU J W,WANG Y G,ZHANG L J,et al.Per-formance study of neural network position esti-mators for the monolithic scintillator PET detec-tor modulesCI/IEEE Nuclear Science Sympo-suim&Medical Imaging Conference.S.1.:s.n.,2010.8van DAM H T,SEIFERT S,VINKE R,et al.A practical method for depth of interaction deter-mination in monolithic scintillator PET detectorsJ.Physics in Medicine and Biology,201l,56(13):4 135-4 145