Non-Sequential Simulation Analysis of Narcissus Effect for Cooled Infrared Optical System
-
摘要: 在制冷红外光学系统中，制冷探测器通过前面的光学表面反射，使得探测器探测到自身的像，形成边缘亮而中心暗的黑斑现象被称为“冷反射”现象。冷反射严重影响制冷红外光学系统的像面非均匀性，在进行光学设计时需对其严格控制。本文以某中波制冷红外光学系统为例，提出了能有效解决或降低冷反射效应的技术途径，对原设计改进优化，利用ZEMAX光学设计软件建立非序列数学模型，进行仿真分析，经分析改进后光学系统冷反射由40%下降到13%，并开展了工程样机成像实验，实验现象与仿真结果基本一致，表明该分析方法能准确地反映红外光学系统中的冷反射实际情况，可作为制冷红外光学系统研制生产前的判定依据。Abstract: In a cooled infrared optical system, the cooled detector can be reflected through the front optical surface, such that the detector can detect its own reflecting image. The phenomenon of a dark spot with a bright edge and a dark center is called "Narcissus, " which has a significant influence on image nonuniformity of the cooled infrared optical system. Narcissus should be strictly controlled in the process of optical design. In this study, a medium wavelength infrared optical system is taken as an example and a technical approach is proposed to effectively solve or reduce narcissus. A non-sequential mathematical model was established based on ZEMAX optical design software and simulation was conducted for analysis. The narcissus of the optical system was reduced from 40% to 13% after optimization and an engineering prototype imaging experiment was conducted. The analysis results of the simulation were consistent with the experimental phenomenon. This analysis method can accurately reflect the actual situation of Narcissus and can be used as a basis for judgment before the development and production of a cooled infrared optical system.
-
Keywords:
- infrared optical system /
- narcissus /
- MWIR /
- non-sequential
-
0. 引言
冷反射是制冷红外系统中一种普遍的杂光效应[1]，也是评价和计算红外系统性能的重要指标之一[2]。冷反射是制冷型红外系统中探测器能够“看到”“自身反射”的“冷像”，它是从系统中各透镜表面反射回来进入光敏面的辐射[3]。因此，在制冷红外光学系统设计中，需合理选择光学结构形式，约束透镜形状、镜面像高、光线入射角、膜层透过率等指标，但上述指标与光学系统成像质量及边界尺寸要求相互制约，在含头罩窗口的光学成像制导系统中影响尤其明显。因此，在制冷红外光学系统研制生产前，通过建立数学仿真模型，进行光线追迹。对冷反射影响强度进行合理评价非常重要。
在进行制冷红外光学系统设计时，通常利用YNI（指近轴边缘光线在该面上的投射高度Y、入射角I和折射率N的乘积）与I/I-BAR（I是近轴边缘光线的入射角，I-BAR是近轴主光线的入射角，I/I-BAR反映的是该面产生的冷像噪声在系统视场范围内的变化情况）两个参数作为光学系统中各表面产生的冷反射强度的评价标准[4]。近年来，已经有相关文献讨论了相关的设计方法，如文献[5-7]推导并总结了冷反射引入温差（Narcissus induced equivalent temperature difference, NITD）和成像系统参数及环境参数的关系式，并进行了定量分析。文献[8]提出了差分阈值制冷红外系统冷反射校正法，文献[9]分析了一种基于光学优化设计的冷反射抑制的设计方法，文献[10]开展了制冷型红外热像仪冷反射效应的数理分析，文献[11-13]讨论了中波红外系统中冷反射的抑制等相关工作。
本文针对某中波红外制冷光学成像制导系统的需求，在满足系统指标要求的前提下合理布局，对系统的冷反射强度进行仿真分析和评估，利用工程样机进行冷反射对比实验，经实测仿真结果与实验现象基本一致，仿真方法准确有效。
1. 冷反射分析方法及流程
在制冷红外光学系统中，可以通过光学设计软件进行光线追迹仿真分析，根据分析结果进一步优化改进设计或者通过后期图像处理算法对图像进行标定校正进而对冷反射进行抑制。本文提出了一种通过ZEMAX光学设计软件，建立非序列数学仿真模型，通过逆向光线追迹方法，对由于镜片表面反射所形成的系统冷反射进行相对强度仿真，并给出仿真结果。冷反射示意图如图 1。
制冷红外光学系统冷反射分析流程如下所述：
1）由于冷反射现象为一种特殊的鬼像，所以在进行冷反射分析前，应先利用ZEMAX鬼像分析功能Ghost focus generator对系统进行鬼像分析，根据鬼像距离探测器像面的距离、成像大小、放大倍率等，找出鬼像比较严重的表面，对表面曲率半径进行修正优化，通过鬼像的筛选，给仿真分析提出方向性建议，减少后期ZEMAX的数据分析量；
2）在ZEMAX中建立非序列仿真模型，光学系统中透镜的透过率依据镀膜经验数据曲线设置，材料自身的透过率按照ZEMAX软件自定义进行设置。系统中所涉及的结构件如镜筒、隔圈、压圈建模导入系统中，并按照零件表面反射率进行赋值；
3）在探测器像面处设置黑体辐射光源及接收探测器，进行真实光线追迹，光源出射能量经透镜表面反射，在探测器光敏面上形成能量分布，即由于系统冷反射所引起的像面非均匀性。
通过上述方法建立数学仿真模型，模型示意图如图 2所示。
2. 消除或降低系统冷反射的技术手段
通过上述分析可知，消除或降低冷反射的主要技术途径有以下3种：
① 在进行光学设计时，严格控制透镜形状和光线偏折角度，提高冷光线在各个透镜表面上的反射角度，使其经过多次震荡返回光线无法反射到探测器像面处，不能形成明显光斑，并通过鬼像分析快速确定造成较严重冷反射的问题面，进行优化设计，在优化过程中自动减小或消除其影响。
② 对结构件尤其是镜筒、隔圈、压圈等进行表面氧化发黑或涂覆消光漆处理，降低结构件表面的反射率，加强对杂散光的散射，降低系统的冷反射。
③ 根据空间尺寸布局，合理选择结构形式，减少镜片数量，提升系统透过率，同时提高光学镜片的膜层透过率，降低单面反射率。从而降低杂散辐射光源经透镜表面反射进入系统形成冷反射。
3. 分析实例
3.1 光学系统设计指标及改进措施
以某中波红外制冷光学系统为例，原设计冷反射现象明显，在满足系统指标要求的情况下，对系统改进设计，进行冷反射仿真分析，同时开展冷反射实验进行对比，表 1为该光学系统的技术指标要求。图 3为改进前后光学系统图。
表 1 光学系统技术指标要求Table 1. The technical requirements of optical systemProject Numerical value Focal length 27 mm FOV 20°×16° F# 2 Format 640×512 Pixel pitch 15 μm×15 μm Operating temperature -40℃~+70℃ 针对光学系统冷反射问题主要采取如下改进措施：
1）优化光学结构，增加光学非球面，减少镜片数量，提高增透膜透过率要求，单个镜片透过率由98%提升至99%，降低由于透镜镀膜剩余反射所产生的冷反射现象，光学系统透过率由88.6%提高至96.1%（不含整流罩）。
2）控制透镜曲率半径，增大各光学表面的YNI数值，最小值由改进前的0.3585增加至1.2798，降低系统的冷反射。
3）利用ZEMAX软件建立非序列仿真模型，反向追迹，对系统各表面进行冷反射追迹，确定影响严重表面，在满足成像质量要求前提下，修改曲率半径，降低系统的冷反射。
改进前后光学系统对比表如表 2所示。
表 2 改进前后光学系统对比Table 2. The comparison of optical system before and after improvementProject Optical system (exclude gairing) Before improvement After improvement Lens number 6 4 YNImin 0.3585 (optical lens 5) 1.2798 (optical lens 1) Optical transmittance 88.6% 96.1% Aspheric surface 0 1(optical lens 3) 分析改进后光学系统各个透镜表面的冷反射影响，经分析透镜1（optical lens 1）的前表面影响较大，在满足成像质量要求前提下，修改曲率半径，进行迭代分析优化，降低系统的冷反射，透镜1优化前后分析结果如图 4。
![]()
图 4 透镜1改进前后分析对比Figure 4. The comparative analysis of optical lens 1 before and after improvement3.2 冷反射分析
来自于光学系统外部非目标产生的杂散辐射及工作过程中内部机械件由于升温产生的杂散辐射是影响像面非均匀性的主要因素。为了降低由于镜筒内壁反射及自身辐射所产生的杂散光影响，在光线经过的镜筒内壁增加消光螺纹和消光漆。假设镜筒内壁为镜面反射，反射率为3%，经过5次反射到达探测器像面，根据公式(1)可以计算到达像面能量。
$$ {R}_{系统}={R}^{n} $$ (1) 式中：R系统为系统杂光总反射率；R为镜筒内壁反射率，R＝3%；n为反射次数[14]，n＝5；经计算，R系统＝2.43×10－8。
通过ZEMAX光学设计软件对改进前后的光学系统进行对比仿真分析（环境参数赋值相同），建立非序列仿真，该模型主要包括整流罩、光学系统、机械件及探测器光敏面4部分。其中光学系统、探测器将模型直接导入，整流罩膜层透过率为0.85，其他光学镜片膜层透过率单面为0.98，结构支撑件反射率设为5%，光源的辐射能量设为1 W，进行实际光线追迹，非序列仿真模型如图 5所示，分析结果如图 6、图 7所示。
![]()
图 6 改进前后像面非均匀性灰度图分析对比Figure 6. The comparison of non-uniformity grayscale image before and after improvement![]()
图 7 改进前后像面非均匀性曲线分析对比Figure 7. The comparison of non-uniformity curve of the image before and after improvement从分析结果可以看出，改进前系统像面非均匀性灰度图有明显凸台边界，改进后系统的像面非均匀性灰度图弥散性好，无明显边界突变；改进前系统像面非相干照度非均匀性为40%，冷反射较明显，改进后系统像面非相干照度非均匀性曲线平滑，中心与边缘差异较小，像面均匀性较好，非均匀性为13%，冷反射现象基本消除。
4. 冷反射成像实验
完成精密装校的工程样机如图 8所示，成像效果如图 9所示。同时开展了工程样机冷反射成像实验，改进前后冷反射成像效果如图 10所示，从实验结果可以看出，改进后系统成像质量优良，冷反射成像效果与仿真结果基本一致，进一步证明仿真方法的有效性。
5. 结论
本文利用ZEMAX光学设计软件建立非序列数学仿真模型，提出了进行冷反射分析的分析流程及有效解决和降低冷反射的效应的技术途径。并以某中波红外光学系统为例，对光学系统进行改进优化设计，经分析改进后光学系统冷反射由40%下降到13%。并开展了工程样机研制和冷反射实验，实验与仿真结果基本一致。表明该分析方法能准确地反映红外光学系统中的冷反射实际情况，可作为制冷红外光学系统研制生产前的判定依据。
-
![]()
图 4 透镜1改进前后分析对比
Figure 4. The comparative analysis of optical lens 1 before and after improvement
![]()
图 6 改进前后像面非均匀性灰度图分析对比
Figure 6. The comparison of non-uniformity grayscale image before and after improvement
![]()
图 7 改进前后像面非均匀性曲线分析对比
Figure 7. The comparison of non-uniformity curve of the image before and after improvement
表 1 光学系统技术指标要求
Table 1 The technical requirements of optical system
Project Numerical value Focal length 27 mm FOV 20°×16° F# 2 Format 640×512 Pixel pitch 15 μm×15 μm Operating temperature -40℃~+70℃ 
下载: 导出CSV
表 2 改进前后光学系统对比
Table 2 The comparison of optical system before and after improvement
Project Optical system (exclude gairing) Before improvement After improvement Lens number 6 4 YNImin 0.3585 (optical lens 5) 1.2798 (optical lens 1) Optical transmittance 88.6% 96.1% Aspheric surface 0 1(optical lens 3)
下载: 导出CSV
-
[1] J M Lloyd. Thermal Imaging Systems[M]. New York: Plenum Press, 1975.
[2] 任国栋, 张良, 兰卫华, 等. 红外成像系统冷反射的定量分析[J]. 红外技术, 2016, 38(4): 290-295. http://hwjs.nvir.cn/article/id/hwjs201604003 REN Guodong, ZHANG Liang, LAN Weihua, et al. Quantitative analysis of the narcissus of infrared imaging system[J]. Infrared Technology, 2016, 38(4): 290-295. http://hwjs.nvir.cn/article/id/hwjs201604003
[3] 刘欣, 潘枝峰. 红外光学系统冷反射分析和定量计算方法[J]. 红外与激光工程, 2012, 41(7): 1684-1688. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201207003.htm LIU Xin, PAN Zhifeng. Analysis and quantitative calculating methods for Narcissus of infrared optical system[J]. Infrared and Laser Engineering, 2012, 41(7): 1684-1688. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201207003.htm
[4] 张鹏, 罗长江, 熊钟秀, 等. 制冷型红外光学系统冷反射的逆光路分析[J]. 电光与控制, 2013, 20(s6): 21-27. https://www.cnki.com.cn/Article/CJFDTOTAL-DGKQ201306017.htm ZHANG Peng, LUO Changjiang, XIONG Zhongxiu, et al. Analysis of inverse path tracing rays of Narcissus for cooled infrared optical system[J]. Electronics Optics & Control, 2013, 20(s6): 21-27. https://www.cnki.com.cn/Article/CJFDTOTAL-DGKQ201306017.htm
[5] 刘洋, 安晓强. 制冷型红外焦平面系统冷反射效应的分析与控制[J]. 光学学报, 2012, 32(2): 0222007-1. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB201202043.htm LIU Yang, AN Xiaoqiang. Analysis and control of Narcissus effect of cooling IR focal plane system[J]. Acta Optical Sinica, 2012, 32(2): 0222007-1. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB201202043.htm
[6] 刘洋, 安晓强, 王茜. 基于光学离焦量的致冷型长波变焦红外成像系统冷反射效应的分析与控制[J]. 光学学报, 2012, 32(4): 0422003-1. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB201204037.htm LIU Yang, AN Xiaoqiang, WANG Qian. Analysis and control of narcissus effect of long-wavelength cooling infrared imaging system based on optical defocusing[J]. Acta Optical Sinica, 2012, 32(2): 0422003-1. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB201204037.htm
[7] 刘志祥, 马冬梅, 胡明鹏, 等. 凝视型红外成像系统中冷像的仿真分析[J]. 红外与激光工程, 2008, 37(4): 702-705. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ200804036.htm LIU Zhixiang, MA Dongmei, HU Mingpeng, et al. Simulation analysis of the narcissus in the staring infrared imaging system[J]. Infrared and Laser Engineering, 2008, 37(4): 702-705. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ200804036.htm
[8] 梅朝阳, 崔庆丰, 胡洋, 等. 差分阈值致冷红外系统冷反射校正法[J]. 光子学报, 2021, 50(11): 1111003-1. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB202111020.htm MEI Chaoyang, CUI Qingfeng, HU Yang, et al. Correction method of Narcissus effect of cooled infrared system based on differential threshold[J]. Acta Photonica Sinica, 2021, 50(11): 1111003-1. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB202111020.htm
[9] 王珊珊, 梁宇宏, 向阳. 一种基于光学优化设计的冷反射抑制[J]. 长春理工大学学报, 2012, 35(3): 17-18. https://www.cnki.com.cn/Article/CJFDTOTAL-CGJM201203006.htm WANG Shanshan, LIANG Yuhong, XIANG Yang. A inhibition of narcissus based on optical optimization design[J]. Journal of Changchun University of Science and Technology, 2012, 35(3): 17-18. https://www.cnki.com.cn/Article/CJFDTOTAL-CGJM201203006.htm
[10] 李忠. 制冷型红外热像仪冷反射效应的数理分析[J]. 舰船科学技术, 2006, 28(6): 88-89. https://www.cnki.com.cn/Article/CJFDTOTAL-JCKX200606026.htm LI Zhong. Mathematical analysis of cooled infrared imaging's reflection cool effect[J]. Ship Science and Technology, 2006, 28(6): 88-89. https://www.cnki.com.cn/Article/CJFDTOTAL-JCKX200606026.htm
[11] 顿雄, 陶玉, 孟军合. 双视场红外扫描成像系统冷反射抑制[J]. 红外与激光工程, 2010, 39(4): 732-735. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201004037.htm DUN Xiong, TAO Yu, MENG Junhe. Narcissus suppression in dual-FOV IR scanned imaging systems[J]. Infrared and Laser Engineering, 2010, 39(4): 732-735. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201004037.htm
[12] 谢洪波, 孟庆斌, 杨磊, 等. 中波红外光学系统无热化设计和冷反射抑制[J]. 应用光学, 2017, 38(3): 353-356. https://www.cnki.com.cn/Article/CJFDTOTAL-YYGX201703003.htm XIE Hongbo, MENG Qingbin, YANG Lei, et al. Athermalization and suppression of narcissus for medium-wave infrared optical system[J]. Journal of Applied Optics, 2017, 38(3): 353-356. https://www.cnki.com.cn/Article/CJFDTOTAL-YYGX201703003.htm
[13] 马永龙. 红外热成像系统中冷反射的分析和控制[J]. 光学与光电技术, 2010, 8(6): 68-70. https://www.cnki.com.cn/Article/CJFDTOTAL-GXGD201006027.htm MA Yonglong. Analysis and restrain of narcissus effect in scanning infrared system[J]. Optics & Optoelectronic Technology, 2010, 8(6): 68-70. https://www.cnki.com.cn/Article/CJFDTOTAL-GXGD201006027.htm
[14] 贺磊, 张建隆, 杨振, 等. 一种小型化滚-仰式长波红外光学系统设计[J]. 红外技术, 2018, 40(12): 1142-1148. http://hwjs.nvir.cn/article/id/hwjs201812005 HE Lei, ZHANG Jianlong, YANG Zhen, et al. Design of a small rolling-pitching long-wave infrared optical system[J]. Infrared Technology, 2018, 40(12): 1142-1148. http://hwjs.nvir.cn/article/id/hwjs201812005
-
期刊类型引用(1)
1. 郭来功，刘安坤，黄晨瑞，穆朝民. 高速红外测温系统. 红外技术. 2025(01): 36-43 . 
本站查看
其他类型引用(1)
计量
- 文章访问数: 217
- HTML全文浏览量: 79
- PDF下载量: 64
- 被引次数: 2
