基于气体探测器的X射线荧光计算机断层扫描成像扫描仪的模拟研究

《Frontiers in Detector Science and Technology》：Simulation of X-ray fluorescence computed tomography imaging scanner based on gaseous detectors

　　

在X射线荧光计算机断层扫描（X-ray Fluorescence Computed Tomography, XFCT）系统中，气体探测器为常用的半导体探测器提供了一种有吸引力的替代方案，主要归因于其在大有效面积和低成本方面的潜力。在本研究中，研究人员利用GATE软件（Geant4 Application for Tomographic Emission）通过蒙特卡罗（Monte Carlo）模拟研究了一种新型XFCT成像概念。该模拟系统集成了基于准直的光学系统以及一个充有氪气的位置和能量灵敏气体探测器，用于荧光光子探测。研究人员模拟了三种具有不同几何形状、尺寸、元素组成和衰减条件的模体（phantom），以评估系统性能。结果表明，该系统能够区分多种元素并重建截面和三维元素图像，从而证实了基于气体探测器的XFCT成像的可行性。这些发现凸显了该方法作为基于半导体的XFCT系统替代方案的潜力，支持进一步的优化和实验验证。

基于气体探测器的XFCT成像模拟研究解读

X射线荧光（X-ray Fluorescence, XRF）技术通过探测样品在电离辐射照射下发射的特征X射线来识别元素，具有多元素同时检测、无损和样品准备简单等优势。计算机断层扫描（Computed Tomography, CT）通过多角度投影重建三维结构。X射线荧光计算机断层扫描（X-ray Fluorescence Computed Tomography, XFCT）将两者结合，可无损地绘制样品内部三维元素分布，在生物医学成像和材料科学中具有应用价值。然而，多数XFCT系统依赖同步辐射单能X射线，限制其常规应用。实验室型系统虽有发展，例如Jones等人报道了用于定量金纳米粒子分布的锥束XFCT系统，Li等人提出了基于CZT线性光子计数阵列的全场扇形束系统，Moktan等人则通过Geant4模拟优化了检测设置，但这些系统多采用半导体探测器，存在复杂、成本高和供应短缺等问题。气体探测器因其成本低、可制成大面积、气体成分可调等优点，在XRF成像中显示出潜力，但尚未被用于XFCT。因此，研究人员提出探索基于气体探测器的XFCT成像概念，通过蒙特卡罗模拟评估其可行性。

研究人员利用GATE软件执行蒙特卡罗模拟，建立了集成了准直光学系统和位置及能量灵敏气体探测器的XFCT扫描仪模型。探测器以氪气为灵敏介质，前方配置铅多孔准直器。为全面评估系统性能，研究人员设计了三种不同复杂度、几何尺寸和衰减条件的模体，通过模拟获得了能谱、正弦图和重建截面图，并实现了三维元素分布重建。结果表明，该系统能够在不同条件下区分多种元素并重建二维和三维元素图像，证实了基于气体探测器的XFCT成像的可行性。该研究成果提供了半导体探测器系统的替代方案，具有显著的成本和规模化优势。论文发表在《Frontiers in Detector Science and Technology》。

在技术方法上，研究人员使用了GATE软件（Geant4 Application for Tomographic Emission）和蒙特卡罗模拟，采用Livermore低能物理列表模拟电磁相互作用。系统的核心部件包括一个4 × 4 × 1 cm³的氪气气体探测器和一个铅多孔准直器（孔直径1 mm或0.5 mm）。每种模体在旋转台上以1°/s的速度完成360°扫描。为降低计算开销，荧光X射线发射被近似为仅发生在物体表面，该近似基于光子透射计算得到验证。探测器能量分辨率设定为18%（@5.9 keV），符合已知实验值。投影数据采用滤波反投影（Filtered Back Projection, FBP）算法（Matlab iradon函数）重建截面图像，再堆叠成三维图像。空间分辨率通过提取径向和切向轮廓并拟合高斯函数得到全宽半最大值（Full Width at Half Maximum, FWHM）进行量化。

结果部分如下：

3.1 Phantom 1：该模体由空气中的铁和铜圆柱组成。能谱中清晰识别出氩（Ar）、铁（Fe）和铜（Cu）的特征峰，正弦图显示两条对应Cu和Fe的正弦轨迹。在旋转角接近90°时，Fe信号因Cu圆柱遮挡而消失。重建图像显示元素位置正确，但FWHM测量值大于理论直径（Fe径向1.76±0.03 mm，切向2.03±0.02 mm；Cu径向3.61±0.31 mm，切向4.43±0.20 mm），这主要归因于有限准直器孔径和重建模糊。两个圆柱间距离测得约10.2±0.1 mm，与设定的10 mm一致。三维重建成功组合了各元素的等值面，但存在由投影噪声引起的伪影。

3.2 Phantom 2：该模体为聚甲基丙烯酸甲酯（PMMA）圆柱内嵌入Fe、Cu、金（Au）三种元素的圆柱插件。能谱显示三个特征峰及Ar峰，正弦图可识别三个独立物体。重建截面中圆柱形状略显扁平，原因是某些旋转角度下元素插件互相阻挡，导致投影数据缺失。FWHM测量值均大于理论直径，且径向偏差大于切向，表明重建不确定性主要与系统空间分辨率和重建过程有关，而非元素类型。中心距离测量值略大于理论值（如Fe–Cu为7.68±0.10 mm，理论为6.00 mm），但元素相对排列仍然保持。三维图清晰显示了三种元素的分布。

3.3 Phantom 3：该模体包含一个水平Cu圆柱、一个垂直Au圆柱和两个Fe球体。能谱显示三个特征峰。由于正弦图沿垂直方向积分信号，两个Fe球在正弦图中表现为一条较宽的轨迹，二维截面图像无法区分它们，但在三维重建中，两个Fe球被清晰地分离并显示正确位置。重建结果存在一定噪声，在Fe球中更为明显，但整体几何和分布被准确恢复。

3.4 Detector point spread function characterisation：为单独评估探测器固有空间响应，研究人员在没有准直器和模体的情况下，使用单能笔形光束（每能量10⁵光子）垂直入射探测器中心，并在3–20 keV范围测量点扩散函数（Point Spread Function, PSF）。结果显示，探测器固有FWHM从3 keV的0.024 mm增加到15 keV的0.512 mm，在6.4、8.0和9.7 keV分别为0.12、0.20和0.31 mm。这一能量依赖展宽源于氪气中光子射程增加和电荷扩散效应。该结果证明探测器固有响应不是系统空间分辨率的限制因素，准直器孔径才是主导因素。

讨论部分指出，重建尺寸的系统性高估由三个因素共同造成：准直器有限孔径（1 mm或0.5 mm）是主要因素；FBP算法在物体边缘产生模糊；探测器固有空间响应贡献可忽略。PSF表征将各贡献分离并量化，确认探测器本身性能优异。尽管存在分辨率限制，所有模体的元素相对空间排列均被保持，表明系统能够恢复样品的整体几何。结论如下：本研究通过GATE模拟验证了基于气体探测器的XFCT成像的可行性，展示了多元素识别和三维元素分布重建能力，并强调了其作为半导体系统替代方案的潜力。未来工作将集中于实验验证（基于正在开发的气体探测器系统）、系统几何和重建过程的优化，以及对气体探测器和半导体光子计数探测器的定量比较。