基于红外多光谱相机分析长后焦距对无热化设计的影响*

徐睆垚 徐亮 沈先春 徐寒杨 孙永丰 刘文清 刘建国

1)(中国科学院合肥物质科学研究院安徽光学精密机械研究所,合肥 230031)

2)(中国科学技术大学,合肥 230026)

在不同环境温度下,热差对红外多光谱相机的成像质量造成一定的影响,基于此,建立了红外多光谱相机的无热化模型,此模型将红外多光谱相机等效为分离式双透镜光学系统.在定焦距的情况下,分析了后焦距变长对前后透镜光焦度的影响,从光焦度绝对值与正负值变化情况对材料的选择范围进行约束,实现光学材料的快速选择.采用该模型对波段为8—14 µm,焦距为50 mm,F 数为1.4的红外多光谱相机在–40—+60 ℃范围内进行无热化设计.通过仿真分析,各视场在奈奎斯特频率为30 lp/mm 处的值均达到0.39,接近衍射极限;弥散斑均方根半径均小于艾里斑半径19.17 µm;轴向像差均小于0.02 mm.采用通道为9.43—11.53 µm的红外多光谱相机对SF6 气体进行成像实验,实验结果表明,经过无热化的红外多光谱相机对SF6 气体的成像效果良好,设计方法正确可行.

1 引 言

随着社会的快速发展,突发性污染事故的发生日趋频繁,例如火灾、化工厂爆炸污染、毒气泄漏等事件使得有毒有害污染气体向大气排放.而传统的光学遥感和现场采样检测方式无法满足事故处置部门对污染现场快速定性和定量远程遥测的需求[1−4].近年来,红外多光谱相机以其快速识别、远距离、低成本等显著优势成为气体检测领域的研究热点之一[5].

由于红外多光谱相机需要在不同温度条件下工作,温度变化产生的热差会导致成像质量下降,这严重影响了仪器在不同温度条件下的使用效果[6].相比于可见光光学系统,红外光学系统的材料热系数高、波段范围大、热差与轴向色差导致的像面偏移更加明显[7].国内外学者针对红外光学系统的无热化设计进行了大量研究:1943 年,Perry[8]对均匀温度场中光学系统的热效应进行了研究,首次提出了无热化的概念;1994 年,Tamagawa 等[9]建立无热图挑选材料实现了密接多透镜组的无热化设计,并在1996 年拓展了这一理论,给出了分离多透镜组的无热化设计实例[10];2012 年,Schwertz等[11]在图解法的基础上,提出一种扩展的图解法来选择光学材料组合,适用于特定膨胀系数的外壳.随着加工技术的进步和衍射光学元件(diffractive optical element,DOE)特性的深入研究,折衍混合光学系统被越来越多地应用于高质量红外光学系统的无热化设计[12−15].但这些设计方法是普适的,针对长后焦距红外光学系统的无热化设计缺乏具体理论分析,毫无目的地选取材料会极大增加无热化设计的难度.鉴于此,本文基于红外多光谱相机分析了长后焦距对无热化设计的影响,通过对材料的选择范围进行约束,可以实现长后焦距红外光学系统材料的快速选择.同时,本文提供了一个长后焦距大相对孔径红外多光谱相机的设计实例,通过仿真分析以及实验结果验证了设计方法的可行性.

2 红外多光谱相机无热化模型

2.1 无热化原理

对于一个具有K块薄透镜的光学系统,分割出K－1 个空间,如图1 所示,Li为第i块透镜,Mi为第i块透镜的材料,hi为近轴光线在第i块透镜上的高度.

图1 具有K 块薄透镜的光学系统Fig.1.An optical system containing K thin lenses.

将热差方程通过后焦距fb随温度T变化的形式表现出来,可以得到

式中,di为第i块透镜与第i+1 块透镜之间的距离,nαi为空气的折射率,ααi为空气的线膨胀系数,ϕi为第i块透镜的光焦度,γi为Mi的热系数,ωi为Mi的色散系数.其中,γi和ωi的表达式为[9]

式中,Δλ是所选的波长范围,αi和ni分别是透镜材料的线膨胀系数以及中心波长处的折射率.(1)式中,大括号内的第1 项表示透镜光焦度随温度的变化,第2 项表示随着di的变化光学系统光程的变化.一般来说,透镜的光焦度变化远大于光学系统的光程变化,为了便于计算,大括号内的第2项可以忽略不计[10].

机械件材料的热胀冷缩同样会导致后焦距fb发生变化:

式中,αh为镜筒材料的线膨胀系数.所以消热差的本质就是通过材料的合理搭配使光学系统本身热效应引起的热差与机械件材料热胀冷缩引起的热差相抵消,再结合光学系统总光焦度φ方程和轴向色差方程可以获得无热化方程组:

2.2 无热化方程组求解

由于长波红外波段光线能量弱,像面照度公式为[16]

式中,N1为空气和材料的透射界面数,N2为光学材料中心厚度总和,ρ1为投射界面的反射损失系数,α为光学材料的吸收系数,G为物方光束亮度,sin2U′为光学系统的数值孔径.通过(6)式可以得出结论:增大相对孔径、减少透镜数量、减小透镜厚度有利于像面照度的提高.

取K=3,对无热化方程组进行求解:

图2 无热图图解 (a) L1的归一化光焦度值;(b)透镜的归一化光焦度正负值Fig.2.Graphic illustration of athermal chart:(a) Normalized power value of L1;(b) positive and negative values of normalized power of lens.

在无热化设计中,在任选材料的情况下,通过无热化方程组求出的光焦度正负值可能与初始结构的光焦度分配情况不一致,需要对材料的选择范围进行约束.假设选用库克三分离物镜作为红外多光谱相机的初始结构,库克三分离物镜将单凸透镜分为L1,L3,各安装在单凹透镜L2的前后一定距离处,形成大体对称式的结构[18].随机挑选材料Ge,ZnSe,ZnS 作为L1,L2,L3的材料,各透镜的权值为1,0.9,0.8,镜筒材料的线膨胀系数为23.6 ×10–6/℃,后焦距为40 mm,焦距为50 mm,表1列出了常用长波红外材料在8—14 µm 波段的热系数及色散系数,通过无热化方程组求出的为(−0.6624,+2.8833,−1.2208),与初始结构的光焦度分配情况相反.

表1 长波红外材料光学特性Table 1. Optical properties of long wave infrared materials.

可以看出,(ϕ1,ϕ2,ϕ3) 需 要找到一组对应的di才可以同时满足消热差且后焦距为fb.但是在实际光学设计中,受透镜厚度以及光学系统总长度的限制,di的变化是受限制的,在随机挑选材料的情况下,(9)式不一定能求解出一组合理的di,需要对材料的选择范围进一步约束.假设将Ge,ZnSe,ZnS替换为Ge,GaAs,AMTIR1,此时的为(+0.4724,−4.2006,+2.2738),由于各透镜的权值、焦距以及后焦距可通过初始结构确定,求出(ϕ1,ϕ2,ϕ3) 为(+0.0049,−0.0933,+0.0568) ,将ϕ1与ϕ2的值代入(9)式可以求出d1=204.08 mm,d2=−1.265mm.此时,L1 与L2的间隔远大于焦距,L2 与L3的间隔变为负值,不符合实际情况.

2.3 长后焦距对无热化设计的影响

对于长后焦距物镜,像方主面位于L3 附近,由L3 承担主要光焦度.为了简化分析,将L1,L2等效为前组透镜Lfront,将L3 作为后组透镜Lback,则原有的系统被等效为分离式双透镜光学系统[19].等效后的消色散方程可以简化为

通过消色散方程可以看出,承担主要光焦度的透镜需要选用较小色散系数的材料且前后组透镜光焦度比例越大,色散系数比例越小.分析不同后焦距下光焦度的变化情况,当后焦距fb增大时,像方主面H后移,得到如下关系式[16]:

以f′=50 mm 为例,结合(12)式与(13)式绘出不同后焦距下前组透镜归一化光焦度的变化情况,如图3 所示.

图3 不同后焦距下前组透镜的归一化光焦度变化情况Fig.3.Change of normalization power of front lens under different back focal length.

从图3 可以看出,通过增大或减小|d|可以抑制前组透镜归一化光焦度的变化趋势,但是需要对|d|进行大幅度的增减,当 Δfb取值较大时,|d|的取值不再合理.所以,当后焦距小于焦距时,随着后焦距的增大呈减小趋势,且后焦距越大趋势越明显,如图3 中的A,C 区域,此时变大,应当选择无热图中相隔较远的两种材料,如图4(a)所示.当后焦距大于焦距时,变化情况完全相反,如图3 中的B,D 区域,此时变小,应当选择无热图中相隔较近的两种材料,如图4(b)所示.

图4 无热图中长后焦距系统材料的选择 (a)后焦距小于焦距;(b)后焦距大于焦距Fig.4.Material selection of long back focal focus system on the athermal chart:(a) Back focal length is smaller than focal length;(b) back focal length is larger than focal length.

3 系统仿真与实验

3.1 设计指标

本文设计的红外多光谱相机选用响应波段为8—14 µm的非制冷型探测器,此探测器的像元尺寸为17 µm,像素为640 × 512,可在–40—+60 ℃环境下工作.机械材料选用铝合金.使用红外多光谱相机检测气体时,待测目标至少需要获得3 × 3个像素.假设需要在100 m远处检测到直径为0.1 m的气体云团,红外多光谱相机的设计指标如表2 所列.

表2 光学系统指标Table 2. Design specifications of optical system.

3.2 设计与仿真

根据2.2 节与2.3 节的分析,可以列出长后焦距红外光学系统的无热化设计流程,如图5 所示.选用库克三分离物镜作为红外多光谱相机的初始结构,在无热图中选择色散系数小的点作为L3的材料,作直线Z1连接,再以无热图中任意两点连接直线Z2,Z1与Z2的交点是等效透镜Lfront的材料.由于后焦距fb略小于焦距要离尽量远,结合2.2 节中的约束方法对材料选择范围进行约束,选择均大于的材料,或均小于的材料.

图5 无热化设计流程图Fig.5.Flow chart of athermal design.

综上,将IRG204,ZnSe,IRG206 分别作为L1,L2,L3的材料,如图6(a)所示.由于透镜数量少,光学系统的高级像差难以校正,在L1 和L3 上各添加一个偶次非球面用来校正高级像差.利用Zemax 软件优化后,在–40—+60 ℃范围内,系统的热离焦均小于一倍焦深,校正了热差.但是无热图中的三角形面积较为扁平,光焦度分配不够均匀,系统依旧存在较大的轴向像差,如图7(a)所示.

图6 长波红外光学材料无热图 (a)折射系统;(b)折衍混合系统Fig.6.Athermal chart of long wave infrared optical materials:(a) Refraction system;(b) diffractive/refractive hybrid system.

考虑到DOE 具有特殊的热系数及色散系数[20]:

式中,αBOE为衍射面基底透镜材料的热膨胀系数,n0为环境介质的折射率,λd为所选波段的中心波长.在L3的偶次非球面基底上添加衍射面,根据(14)式和(15)式求出DOE的热系数和色散系数分别为γDOE=4.14×10−5,ωDOE=0.545 .DOE 和IRG206 等效成透镜Lj,等效材料Mj在无热图中的坐标轨迹在DOE 和IRG206的连线及其延长线上,如图6(b)所示,衍射元件为前组透镜承担一部分正色差,光焦度分配更加平滑,在校正色差的同时也减小了轴上点球差.添加衍射面后,对系统再次优化,系统的轴向像差从0.06689 mm 减小到0.01604 mm,如图7(b)所示.系统的热离焦由–12.46—+10.42 µm 减小到–2.52 —+3.86 µm,远小于一倍焦深,如图8 所示.

图7 轴向像差 (a)折射系统;(b)折衍混合系统Fig.7.Axial aberration:(a) Refraction system;(b) diffractive/refractive hybrid system.

图8 热离焦曲线Fig.8.Athermal defocus curve.

3.3 成像质量评价

经过无热化设计后的系统总长为86.53 mm,后焦距为40.53 mm,结构紧凑,如图9 所示.对其通过调制传递函数(modulation transfer function,MTF)和弥散斑均方根(root mean square,RMS)半径来评价成像质量[21].无热化后的红外多光谱相机在–40—+60 ℃范围内,各视场在奈奎斯特频率(30 lp/mm)处的MTF 值均达到0.39,接近衍射极限,如图10 所示.RMS 半径均小于艾里斑半径19.17 µm,如表3 所列.综上,该系统成像质量优异,基本不受环境温度的影响.

图9 光学系统结构图Fig.9.Structure chart of optical system.

图10 折衍射系统调制传递函数曲线图 (a) +20 ℃;(b)–40 ℃;(c) +60 ℃Fig.10.MTF chart of refractive/diffractive system:(a) +20 ℃;(b)–40 ℃;(c) +60 ℃.

表3 –40 —+60 ℃范围内弥散斑均方根半径Table 3. RMS radius in the range of–40—+60 ℃.

3.4 气体成像实验

在实验室中,采用通道为9.43—11.53 µm的红外多光谱相机对SF6 气体进行成像实验,选用SF6的原因是其不易燃、无毒、价格适中[22].气体成像实验于2020 年12 月31 日在中国安徽省合肥市的实验室中进行,实验室的环境温度经温度计测量是5 ℃,并于2021 年01 月27 日在同一地点对仪器进行温控实验,依旧采用9.43—11.53 µm的通道分别测试了仪器在20 和40 ℃情况下对SF6的成像结果.两次实验结果表明,在未采用任何图像处理的情况下,自研的红外多光谱相机可以在大温度范围内对SF6气体清晰成像,如图11 所示,光学系统性能优良,达到测试指标.

图11 红外多光谱相机对SF6 气体成像 (a)样机;(b) 2020 年12 月31 日实验背景;(c) 5 ℃环境下成像结果;(d) 2021 年01 月27 日实验背景;(e) 20 ℃环境下成像结果;(f) 40 ℃环境下成像结果Fig.11.SF6 gas imaging by infrared multispectral camera:(a) Prototype;(b) background of the experiment on December 31,2020;(c) gas imaging results in 5 ℃ environment;(d) background of the experiment on January 27,2021;(e) gas imaging results in 20 ℃environment;(f) gas imaging results in 40 ℃ environment.

4 结 论

本文建立了红外多光谱相机的无热化模型,此模型将红外多光谱相机等效为分离式双透镜光学系统.通过分析得出结论:由于前后组透镜的间隔受到约束,当后焦距小于焦距时,前后组透镜材料在无热图中的间隔随着后焦距的增大而增大;当后焦距小于焦距时,前后组透镜材料在无热图中的间隔随着后焦距的增大而减小.设计结果以及实验结果表明,该模型极大地简化了长后焦距红外光学系统的无热化过程,能够使红外多光谱相机在不同环境温度下工作.