Preliminary study on scatter quantification method for flash Multi-MeV radiography
-
摘要: 高能MeV闪光照相所针对的客体通常具有极高的面密度。当X射线穿过客体时,直穿X射线的强度将被极大衰减,到达成像面的直穿信号可能被散射“噪声”所淹没,若直接对图像进行反演将严重影响照相重建精度。从散射抑制角度出发,目前主要采用网柵相机即阵列型准直孔阻挡散射,但网栅相机的应用效果受光源位置稳定性影响较大,且网栅不易加工。提出了一种可实时定量散射强度的照相方案,该方案利用狭缝准直器对散射的抑制能力不随散射强度变化而改变这一特点,对现有照相布局进行小改进,利用已知客体实验结果标定狭缝准直器对散射的抑制能力,进一步自洽确定待测客体的散射量大小。基于蒙卡方法的仿真照相实验结果表明,当采用低面密度客体标定散射抑制系数时,高面密度客体散射强度的估计值与模拟真实值偏差可小于2%。Abstract: For multi-MeV X-ray flash radiography, the areal density of object can be obtained by the primary direct X-ray. Objects of flash radiography often have very high areal densities which greatly attenuate the intensity of direct X-rays emitted by the source. At this time, the direct penetration signal that can transmit the region of interest inside the object will be smaller than that of the scattered X-ray “noise”. If the captured image is reconstructed directly without scatter correction, it will affect the accuracy of reconstruction. The main method to reduce the scatter X-ray from a physical point of view is to use an anti-scatter grid, that is, an array-type collimation hole. However, the performance may be affected by the stability of the X-ray source spot, and the manufacture of such anti-scatter grid is very difficult. This paper proposes a new imaging method that does not rely on anti-scatter grid. This method only makes small improvements on the existing imaging layout, and can easily and self-consistently determine the amount of scattering for scattering correction. A Monte Carlo simulation is given to show the performance of scatter estimation, and the relative difference between the estimated intensity of scatter and the real value (provided by the simulation) for an object is less than 2% when a known object with less areal density is applied for calibration.
-
Key words:
- high-energy X-ray /
- flash photography /
- imaging plan /
- image inversion /
- scattering compensation
-
表 1 采用FTO-U30客体标定R用于估计其余客体散射值
Table 1. Calibration of R by FTO-U30 object to estimate the scattering value of other objects
object name variation before and after slit
collimation (measurement)simulated
scatter valuecalculated
scatter valuescatter after
collimation (simulation)R relative
deviation/%FTO-U30(for calibration) 0.0294 0.064 — 0.0059 0.0975 — FTO-U40 0.035 0.055 0.056 — 0.0975 +1.8 FTO-U45 0.035 0.055 0.056 — 0.0975 +1.8 表 2 采用FTO-U45客体标定R用于估计其余客体散射值
Table 2. Calibration of R by FTO-U45 object to estimate the scattering value of other objects
object name variation before and after slit
collimation (measurement)simulated
scatter valuecalculated
scatter valuescatter after
collimation (simulation)R relative
deviation/%FTO-U45(for calibration) 0.0354 0.055 — 0.0052 0.0945 — FTO-U30 0.0294 0.062 0.046 — 0.0945 −25.6 FTO-U35 0.034 0.060 0.053 — 0.0945 −11.7 -
[1] 施将君, 刘军, 刘进, 等. 闪光照相中源尺寸和散射对成像的影响[J]. 强激光与粒子束, 2006, 18(10):1727-1730Shi Jiangjun, Liu Jun, Liu Jin, et al. Effects of source spot size and scatter on imaging in radiography[J]. High Power Laser and Particle Beams, 2006, 18(10): 1727-1730 [2] 刘军, 张绚, 刘进, 等. 闪光照相中散射分布均匀性的影响因素[J]. 强激光与粒子束, 2009, 21(3):447-450Liu Jun, Zhang Xuan, Liu Jin, et al. Factors affecting scatter uniformity in flash X-ray radiography[J]. High Power Laser and Particle Beams, 2009, 21(3): 447-450 [3] Liu Jin, Liu Jun, Jing Yuefeng, et al. Decreasing the scatter effect in density reconstruction in high-energy X-ray radiography[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2013, 716: 86-89. [4] Cao Liangzhi, Wang Mengqi, Wu Hongchun, et al. Scattering correction for image reconstruction in flash radiography[J]. Nuclear Engineering and Technology, 2013, 45(4): 529-538. doi: 10.5516/NET.08.2012.074 [5] 刘军, 张绚, 刘进, 等. 高能闪光照相中网栅相机的性能分析[J]. 强激光与粒子束, 2016, 28:024003 doi: 10.11884/HPLPB201628.024003Liu Jun, Zhang Xuan, Liu Jin, et al. Performance analysis of anti-scatter grid camera in high energy flash radiography[J]. High Power Laser and Particle Beams, 2016, 28: 024003 doi: 10.11884/HPLPB201628.024003 [6] Deng Li, Li Gang, Zhang Baoyin, et al. JMCT Monte Carlo code with capability of integrating nuclear system feedback[C]//Proceedings of the 2018 2nd International Conference on Applied Mathematics, Modelling and Statistics Application (AMMSA 2018). 2018: 48-54. [7] Ekdahl C, Abeyta E O, Bartsch R, et al. DARHT-II long-pulse beam-dynamics experiments[C]//Proceedings of the 2005 Particle Accelerator Conference. Knoxville. 2005: 19-23. [8] Watson S A, Appleby M, Klinger J, et al. Design, fabrication and testing of a large anti-scatter grid for megavolt γ-ray imaging[C]//IEEE Nuclear Science Symposium Conference Record. 2005: 717-721.