Whole-process simulation of NaI (TI) scintillator detector response in radiation field
-
摘要:
NaI(TI)探测器是典型的闪烁体辐射探测器，其探测过程涉及辐射能量沉积、可见光信号产生与输运、光电转换与信号处理等物理过程。首先利用蒙特卡罗方法、Birks公式及射线追迹程序，开展了射线粒子在晶体中转为可见光输出过程的计算分析，并结合光电倍增管和信号处理电路的指标参数进行模拟仿真，得出探测器最终输出的脉冲电压信号参数；然后，在137Cs源辐射场中采用Φ50 mm×50 mm NaI(TI)晶体耦合光电倍增管开展实验验证，实验测得探测器输出脉冲信号的上升/下降时间比为0.39，与模拟计算数值0.36相比，相差约7.69%，表明模拟计算模型的输出结果与实测数据基本符合，初步证明了论文的模拟计算模型及计算分析过程的正确性。论文提出的方法，对于深入理解辐射粒子激发的荧光可见光在晶体闪烁体中的传输规律和闪烁体辐射探测器系统的优化设计，具有一定参考意义。
Abstract:NaI (TI) detector is a typical scintillator radiation detector, and its detection process involves radiation energy deposition, generation and transport of visible light signal, photoelectric conversion, signal processing and other physical processes. Firstly, the process of ray particles converting into visible light signal output in crystals was simulated and analyzed by Monte Carlo method, Birks formula and ray tracing program. Combined with the index parameters of photomultiplier tube and signal processing circuit, the pulse voltage signal parameters for the final output of the detector were obtained. Then, the experiment was verified by using Φ50 mm×50 mm NaI (TI) crystal-coupled photomultiplier in 137Cs source radiation field. The experimental results show that the rise/fall time ratio of the output pulse signal of detector is 0.39, which is about 7.69% different than the simulated value of 0.36, indicating that the output results of the simulation calculation model are basically consistent with the measured data, which preliminarily proves the correctness of the simulation calculation model and calculation analysis process. The proposed method is of certain reference significance for deeply understanding transmission law of fluorescent visible light excited by radiation particles in crystal scintillator and optimization design of scintillator radiation detector system.
-
Keywords:
- radiation detection /
- NaI (TI) crystal /
- Monte Carlo /
- transport of visible light
-
引言
随着数字视频制作技术和数字投影机技术的快速发展，球幕数字投影越来越受到青睐，被广泛应用于天文、地理教学、展览展示等领域[1-3]。然而，市场上数字投影机的品牌和规格繁多，所需要配制的数字鱼眼镜头的规格也很多，检测为数字投影机配制的数字鱼眼镜头便成为了难题。目前，国内外还没有相关的检测设备对数字鱼眼镜头的性能、放映效果进行全面的测试，只能购买价格昂贵的与之匹配的数字投影机来进行性能、放映效果的测试，这样不仅增加了测试费用成本，而且通用性差。为此，研制一种通用型可检测各种类型数字投影鱼眼镜头的检测系统，变得十分迫切。本文研究讨论了一种鱼眼镜头检测系统，可实现芯片尺寸为0.55″~1.55″的1DLP、3DLP、3LCD、3LCOS的各类型的数字投影或放映鱼眼镜头的像方视场角、分辨率、彩色还原性能、色差等性能和放映效果进行测试，通用性强。
1 技术要求
设计一种检测系统，可以满足0.55″~1.55″的1DLP、3DLP、3LCD、3LCOS的各类型的数字投影或放映鱼眼镜头的性能参数和放映效果的测试，可实现的照明视场范围大于φ8.5 mm~φ35 mm，照明视场的照度均匀性应不低于85%，出射光通量应不小于5 000 lm。
2 设计思路
数字鱼眼镜头的设计是根据所配接的数字投影或放映机的关键参数和尺寸而定，如数字鱼眼镜头的像高由芯片尺寸而确定。为了更好地匹配系统整体像差，数字鱼眼镜头在设计时需要加入等效平行平板来模拟数字投影或放映机内部光学引擎中分合色棱镜[4-7]。
要满足不同规格的数字鱼眼镜头的性能(视场角、分辨率、彩色还原性能、色差)和放映效果的测试，检测系统需要具有不同种类和规格的检测标板，同时具有模拟数字投影或放映机内部分合色棱镜的功能，以及具有可以形成均匀且大小尺寸能改变的照明视场的功能。
3 设计过程
3.1 鉴别率标板的设计
数字投影或放映机的分辨率直接决定了所配接的数字鱼眼镜头的极限分辨率[4, 8]，由公式(1)得出：
$$ ~{{R}_{L}}=\frac{1}{2\times \text{ }\frac{{{L}_{P}}}{{{R}_{P}}}~}=\frac{{{R}_{P}}}{2\times {{L}_{P}}}~ $$ (1) 式中：RL为数字鱼眼镜头的极限分辨率；RL为数字投影或放映机芯片的横向分辨率；LP为数字投影或放映机芯片的横向尺寸。
表 1列出了主流数字投影或放映机对数字鱼眼镜头极限分辨率的要求。一般情况下，镜头设计时会将极限分辨率的设计值提高到理论值的1.3到1.5倍，所以，实际数字鱼眼镜头的最高极限分辨率可能会达到100 lp/mm，甚至会更高。考虑到检测系统的通用型，鉴别率标板的单个分辨率图案包含5组线对，分别为125 lp/mm、100 lp/mm、80 lp/mm、63 lp/mm、50 lp/mm。鉴别率标板的尺寸根据各种规格的数字鱼眼镜头的成像像高而确定，每种规格的鉴别率标板在0、0.3、0.5、0.7、0.85视场上均匀分布着单个分辨率图案，如图 1所示，给出了鉴别率标板的示意图。
表 1 不同数字投影机对应的镜头极限分辨率Table 1. Different types of digital projectors and corresponding ultimate resolutions of lenses芯片尺寸/in 芯片像素 极限分辨率/lp/mm 0.76 1 920×1 200 59 0.7 1024×768 36 0.67 1 920×1 200 66.5 0.65 1 920×1 080、1 280×800、1 024×768 66.8 /45.7 /38.8 1.2 2 048×1 080 38 1.38 4 196×2 160 67 3.2 照明系统的设计
目前投影显示系统中主要采用复眼透镜照明系统和方棒照明系统2种照明方式[9]。
复眼透镜照明系统采用双排复眼透镜，属于柯拉照明。照明系统由于光源的整个光束被分为多个细光束照明，而且对称位置上的细光束光斑互相叠加，所以形成了较为均匀的照明场。但是复眼透镜照明系统结构较为复杂，很难保证2组复眼透镜严格的同轴，加工较为困难。
而方棒照明系统采用方棒形式，光在光棒里面多次反射，每次反射都会形成虚拟光源像，多次反射形成二维的虚拟光源矩阵，从而使得光更加均匀。但是方棒照明系统存在能量利用率相对较低，而且系统结构较长。
为了克服现有投影照明系统的不足和缺点，研究讨论了采用柱体(六边形)蝇眼透镜照明系统来提供高均匀的照明场。相比双排复眼透镜结构，柱体蝇眼透镜结构不但避免了在加工工艺中产生的两组蝇眼镜不同轴的误差和前后组蝇眼之间的粘合剂对光的吸收、散射以及杂光的入侵问题，而且还减少了加工面，简化了机械固定结构；相比方棒结构，能量利用率更高，系统的结构更加紧凑，有利于小型化。
系统由抛物型的反光碗光源、柱体蝇眼透镜阵列、准直物镜(场镜和集光镜)组成。其工作原理图如图 2所示，抛物型的反光碗光源发出类平行光，类平行光进入柱体蝇眼透镜阵列被形成多点光源，多点光源经准直物镜的调制之后均匀叠加在系统片窗上，如图 3所示。
3.3 照明系统中关键元件基本参数的确定
准直物镜的焦距f准与照明场的大小y存在以下关系：
$$ {f_{准}} = \frac{{y/2}}{{{\rm{tan}}w}} $$ (2) 式中，w为准直透镜的视场角，如图 4所示。
准直物镜的焦距f准与柱体蝇眼透镜组的子午方向上下边缘2个柱体中心距H(如图 3所示)的关系：
$$ H = {f_{准}}\cdot{\rm{tan}}\theta $$ (3) 式中，θ是准直物镜光学系统的孔径角。
该投影照明系统中柱体蝇眼透镜组的焦平面与准直物镜的物平面重合，即单个柱体透镜的像方孔径角也是准直物镜的视场角，故有：
$$ {f_{柱}} = \frac{{h/2}}{{{\rm{tan}}w}}{\rm{ }} $$ (4) 式中：f柱是单个柱体透镜的焦距；h是单个柱体透镜上下边的距离。
柱体蝇眼透镜组的子午方向上下边缘2个柱体中心距H与单个柱体透镜上下边的距离h之间的关系如下：
$$ ~H=h(n-1) $$ (5) 式中n为柱体蝇眼透镜组子午方向柱体透镜的总个数。
照明系统中，柱体透镜的排列方式，如图 3所示，故，柱体蝇眼透镜组的小透镜的总个数为
$$ N=1+\sum\limits_{1}^{\left( n-1 \right)/2}{\text{ }6\cdot \frac{\left( n-1 \right)}{2}}\text{ }~ $$ (6) 考虑到缩短整个光路，以及有效地提高照度均匀度，取准直透镜的视场角w=6°，准直物镜光学系统的孔径角θ=4°，再由各参数之间的关系式(2)~(6)，最终得出关键元件的设计尺寸。
3.4 照明视场大小可变功能的实现
测试系统需要满足不同类型数字鱼眼镜头的检测，故所需的照明视场大小也不同，为了防止光能量的不必要损失，测试系统应该具有照明视场大小可调的功能。由(2)式知，当系统的准直透镜的视场角确定后，改变准直透镜的焦距，就可以改变照明场的大小。通过改变场镜和集光镜之间的距离便可以得到不同焦距值的准直透镜[10]，进而实现照明视场的大小变化。
4 设计结果
设计了一款通用型数字鱼眼镜头检测系统，如图 4所示，主要由照明系统、检测标板、等效平行平板系统和球幕组成。照明系统采用的柱体蝇眼透镜阵列是由91个六边形柱透镜组成，测试标板主要是鉴别率板，用于检测数字鱼眼镜头的分辨率；等效平行平板系统是由4种不同厚度的平行平板组成，不同的组合可以满足不同类型数字鱼眼镜头的测试。球幕内部刻有角度线，用于检测数字鱼眼镜头的视场角。
经过定量化性能测试，主要性能指标测试结果如表 2所示，满足了设计指标要求。
表 2 检测系统测试结果Table 2. Testing result of testing system主要性能 测试结果 参数指标 照明视场范围φ/mm 8.1~36.3 ＞8.5~35 照度均匀度/% 91 ≥85 出射光通量/lm ≥5 000，可调 ≥5 000 5 结论
本文设计的通用型数字鱼眼镜头检测系统采用的柱体蝇眼透镜阵列照明系统不但光能利用率高，可产生较高的均匀度，而且有较高的准直性和较好的光束控制。通过调整更换不同厚度的等效平行玻璃平板，用于模拟各种数字投影或放映机内部棱镜，同时在片门处调换不同类型、不同尺寸的测试标板，被检测的数字鱼眼镜头将测试标板成像在球幕上，可实现对芯片尺寸为0.55″~1.55″的1DLP、3DLP、3LCD、3LCOS的各种类型的数字投影或放映鱼眼镜头的像方视场角、分辨率、彩色还原性能、色差等性能和放映效果进行测试，通用性强，降低了鱼眼镜头产品的检测成本。
-
![]()
图 1 NaI(TI)晶体闪烁体探测器结构简图
1.铝外壳；2.磁屏蔽套；3.反射层；4.闪烁晶体；5.硅胶；6.光学玻璃（光导）；7.光电倍增管；8.信号处理电路
Figure 1. Brief diagram of NaI (TI) crystal scintillator detector structure
![]()
图 4 NaI(TI)晶体探头的MCNP计算模型
1.放射源 2.空气 3.铝外壳 4.MgO反射层 5.NaI(Tl)晶体6.硅胶 7.光学玻璃。
Figure 4. MCNP calculation model of NaI (TI) crystal probe
![]()
图 5 不同角度入射能量沉积分布
Figure 5. Distribution of incident energy deposition at different angles
![]()
图 6 部分可见光在闪烁体中的传播轨迹
Figure 6. Propagation trajectories of partial visible light in scintillator
![]()
图 7 不同角度入射时晶体输出可见光功率参数
Figure 7. Visible light power parameters for crystal output at different incidence angles
-
[1] 任国浩. 无机闪烁晶体在我国的发展史[J]. 人工晶体学报,2019,48(8):1373-1385. doi: 10.3969/j.issn.1000-985X.2019.08.001 REN Guohao. Development history of inorganic scintillation crystals in China[J]. Journal of Synthetic Crystals,2019,48(8):1373-1385. doi: 10.3969/j.issn.1000-985X.2019.08.001
[2] KOROLEVA T S, SHULGIN B V, PEDRINI Ch, et al. New scintillation materials and scintiblocs for neutron and γ-rays registration[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2005, 537(1/2): 415-423.
[3] 安少康. NPTool模拟NaI(Tl)和LaBr_3(Ce)闪烁体探测器试验研究[J]. 铀矿冶, 2019, 38(3): 198-203. AN Shaokang. Experimental study on NPTools Simulation of NaI(Tl) and LaBr3(Ce) scintillator detectors[J]. Uranium Minig and Metallurgy, 2019, 38(3): 198-203.
[4] 袁航. NaI(Tl)和CsI(Na)闪烁体探测器的探测性能研究[D]. 北京: 华北电力大学, 2019. YUAN Hang. Detection performance of NaI(Tl) and CsI(Na) scintillator detectors[D]. Beijing: North China Electric Power University, 2019.
[5] DEMIR N, KULUÖZTÜRK Z N. Determination of energy resolution for a NaI(Tl) detector modeled with FLUKA code[J]. Nuclear Engineering and Technology,2021,53(11):3759-3763. doi: 10.1016/j.net.2021.05.017
[6] MOUHTI I, ELANIQUE A, MESSOUS M Y, et al. Validation of a NaI(Tl) and LaBr3(Ce) detector’s models via measurements and Monte Carlo simulations[J]. Journal of Radiation Research and Applied Sciences,2018,11(4):335-339. doi: 10.1016/j.jrras.2018.06.003
[7] HIRANO Y, ZENIYA T, IIDA H. Monte Carlo simulation of scintillation photons for the design of a high-resolution SPECT detector dedicated to human brain[J]. Annals of Nuclear Medicine,2012,26(3):214-221. doi: 10.1007/s12149-011-0561-4
[8] MOUATASSIM S, COSTA G J, GUILLAUME G, et al. The light yield response of NE213 organic scintillators to charged particles resulting from neutron interactions[J]. Nuclear Instruments and Methods in Physics Research Section A:Accelerators, Spectrometers, Detectors and Associated Equipment,1995,359(3):530-536.
[9] BIRKS J B, FIRK F W K. The theory and practice of scintillation counting[J]. Physics Today,1965,18(8):60.
[10] BIZARRI G, MOSES W W, SINGH J, et al. An analytical model of nonproportional scintillator light yield in terms of recombination rates[J]. Journal of Applied Physics,2009,105(4):044507. doi: 10.1063/1.3081651
[11] PAYNE S A, CHEREPY N J, HULL G, et al. Nonproportionality of scintillator detectors: theory and experiment[J]. IEEE Transactions on Nuclear Science,2009,56(4):2506-2512. doi: 10.1109/TNS.2009.2023657
[12] GAO F, XIE Y, KERISIT S, et al. Yield, variance and spatial distribution of electron–hole pairs in CsI[J]. Nuclear Instruments and Methods in Physics Research Section A:Accelerators, Spectrometers, Detectors and Associated Equipment,2011,652(1):564-567. doi: 10.1016/j.nima.2010.08.063
[13] WIRTH S, METZGER W, PHAM-GIA K, et al. Impact of photon transport properties on the detection efficiency of scintillator arrays[C]//2006 IEEE Nuclear Science Symposium Conference Record. New York: IEEE, 2006: 2602-2603.
[14] MENG F, MA X, ZHAO H, et al. A study of polarized light propagation in turbid medium by monte carlo simulations and experiments[C]//2006 International Symposium on Biophotonics, Nanophotonics and Metamaterials. New York: IEEE, 2006: 132-135.
[15] WANG J G, WANG G Y, XU Z Z. Monte Carlo simulation of photon migration path in turbid media[J]. Chinese Optics Letters,2008,6(7):530-532. doi: 10.3788/COL20080607.0530
[16] 王建岗, 王桂英, 徐至展. 光在分层散射介质中传输行为的蒙特卡罗模拟研究[J]. 光学学报,2000,20(3):346-350. doi: 10.3321/j.issn:0253-2239.2000.03.011 WANG Jiangang, WANG Guiying, XU Zhizhan. Monte Carlo simulations for light propagationin striated scattering medium[J]. Acta Optica Sinica,2000,20(3):346-350. doi: 10.3321/j.issn:0253-2239.2000.03.011
[17] 郑兴荣, 宋小永, 郑燕飞, 等. 基于MATLAB的光传播特性的数值计算[J]. 陇东学院学报,2017,28(5):36-40. doi: 10.3969/j.issn.1674-1730.2017.05.008 ZHENG Xingrong, SONG Xiaoyong, ZHENG Yanfei, et al. Numerical calculation on the propagation characteristics of light based on Matlab[J]. Journal of Longdong University,2017,28(5):36-40. doi: 10.3969/j.issn.1674-1730.2017.05.008
[18] BRIESMEISTER J F. MCNPTM-A general Monte Carlo N-particle transport code, Version 4C[R]. USA: LA-13709-M, Los Alamos National Laboratory, 2000.
[19] 彭礼韬. MCNP和Geant4在伽马测井领域的应用对比研究[D]. 北京: 华北电力大学, 2021. PENG Litao. Comparative study on the application of MCNP and Geant4 in well logging[D]. Beijing: School of Nuclear Science and Engineering, 2021.
[20] 马继明, 朱宏权, 王奎禄, 等. 硅酸镥γ图像转换屏荧光弥散特性研究[J]. 应用光学,2010,31(5):752-756. doi: 10.3969/j.issn.1002-2082.2010.05.016 MA Jiming, ZHU Hongquan, WANG Kuilu, et al. Fluorescence dispersion of Lu2SiO5 crystal as a γ-image converter[J]. Journal of Applied Optics,2010,31(5):752-756. doi: 10.3969/j.issn.1002-2082.2010.05.016
[21] 田秀劳. 光波在左手材料中的菲涅尔公式和布儒斯特定律[J]. 光子学报,2006,35(7):1103-1106. TIAN Xiulao. Fresnel formulate and Brewster law of optical waves in the left-handed materials[J]. Acta Photonica Sinica,2006,35(7):1103-1106.
[22] BELL I. Introduction to circuit simulation with LTspice-Part 4[J]. Everday Partical Electronics,2019,48(1):48-51.
[23] KONSTANTINOU G, LECOQ P, BENLLOCH J M, et al. Metascintillators for ultrafast gamma detectors: a review of current state and future perspectives[J]. IEEE Transactions on Radiation and Plasma Medical Sciences,2022,6(1):5-15. doi: 10.1109/TRPMS.2021.3069624
-
期刊类型引用(1)
1. 路文文，陈善勇，翟德德，熊玉朋. 混合仿生鱼眼-复眼的广角高清成像系统. 应用光学. 2019(02): 311-315 . 
本站查看
其他类型引用(1)
下载:
计量
- 文章访问数: 344
- HTML全文浏览量: 60
- PDF下载量: 88
- 被引次数: 2
陕公网安备 61011302001501号 