Optical design of ultraviolet continuous zoom lens for sulfur dioxide remote sensing monitoring
-
摘要:
为了满足大气二氧化硫遥感监测在不同距离下高清晰度成像的使用需求,设计了一种连续变焦紫外镜头。通过对变焦系统高斯光学解的理论分析,选取四组元正组补偿机械变焦初始结构,计算了二组元机械补偿变焦系统变倍组求解的3个区间,得到了变倍组与补偿组的光焦度分配;根据大气二氧化硫对紫外光的吸收特性,确定工作波段为250 nm~340 nm;利用Zemax光学设计软件设计了变倍比为10、F数为4、总长431.5 mm的紫外连续变焦光学系统,并使用Matlab软件绘制了变倍组与补偿组的运动曲线。模拟结果表明:系统的焦距范围为25 mm~250 mm,短焦视场范围0°~7.67°、中焦视场范围0°~1.45°、长焦视场范围0°~0.77°,变焦过程中传递函数在107 lp/mm下均大于0.3,畸变均小于1%,弥散斑半径小于5.5 μm。本系统适用于船舶尾气、工业生产烟囱尾气等多种场景的二氧化硫监测。
Abstract:In order to meet the needs of high-definition imaging of atmospheric sulfur dioxide remote sensing monitoring at different distances, a continuous zoom ultraviolet lens was designed. Based on the theoretical analysis of Gaussian optical solution of zoom system, four groups of positive compensation mechanical zoom initial structure were selected, three intervals of zoom group solution of two-group mechanical compensation zoom system were calculated, and the power distribution of zoom group and compensation group was obtained. According to the absorption characteristics of atmospheric sulfur dioxide to ultraviolet light, the working band was determined to be 250 nm~340 nm. An ultraviolet continuous zoom optical system with zoom ratio of 10, F number of 4 and total length of 431.5 mm was designed by using Zemax optical design software, and the motion curves of zoom group and compensation group were drawn by using Matlab software. The focal length range of the system is 25 mm~250 mm, the short focal field range is 0°~7.67°, the medium focal field range is 0°~1.45°, and the telephoto field range is 0°~0.77°. During zoom process, the transfer function is greater than 0.3 at 107 lp/mm, the distortion is less than 1%, and the radius of the speckle is less than 5.5 μm. This system is suitable for monitoring sulfur dioxide in various scenarios, such as ship exhaust, industrial production chimney exhaust and so on.
-
引言
飞机类目标作为现代战争的重要平台,广泛应用于侦察、攻击、运输等领域。基于此,其也成为红外探测系统、红外制导武器探测和打击的主要对象。准确获取来袭飞机的红外辐射特性,对于红外预警探测具有重要意义[1],因此,测量和研究飞机类目标的红外辐射特性是一项重要的基础工作,需长期进行。受保障条件限制,单独组织飞机类目标的红外辐射特性测量困难较大,一般会采取“搭车”方式测量。例如组织侦察探测实验时,利用已有的红外告警设备,结合实验开展目标红外辐射特性测量。由于红外告警设备并非专用的目标红外辐射特性测量设备,在具备原始图像采集的基础上,如何进行准确标定往往是测量前需要解决的关键问题[2]。
针对红外探测系统辐射标定的研究已有诸多报道,例如,吕原等提出了绝对辐射标定和基于多温度场的场地绝对辐射标定方法,实现了实时高精度绝对辐射标定[3]。蔡李靖等提出了一种无黑体式红外热像仪测温标定和温度补偿技术,并通过实验进行了验证[4]。凌玲等对傅里叶变换红外光谱仪开展了实验室标定,并分析了光谱仪测量过程中的不确定来源[5]。李宁等针对3 μm~5 μm红外焦平面阵列的标定拟合及评估,提出了利用剔除局外点的方法对拟合曲线进行优化处理,改善了辐射标定拟合的不确定度[6]。杨词银、李宪圣、李宁等针对大口径红外光电系统,分别提出了基于大面源黑体的辐射标定方法、用小口径照明光管标定大口径红外探测系统的方法和基于光栅单色仪和标准辐射计的相对光谱标定方法,并进行了实验验证及误差分析[7-10]。邹前进等采用不同的标定方法对焦距固定红外探测系统进行了标定实验,并对不同方法的适用情况进行了分析[11]。禄晓飞等基于红外标定参数受外场环境影响因素,提出了外场标定参数验证的方法流程,并进行了精度检验[12]。本文在此基础上,针对外场实际使用需求,分析验证了近距离直接扩展源法、远距离小源法对红外告警设备的标定误差,并给出了使用建议。
1 辐射标定及方法
标定一般采用红外探测系统测量标准红外辐射源的方式进行,目的是建立红外探测系统数字图像输出值(灰度)与其接收到的辐射亮度之间的对应关系,以便在目标特性测量过程中,准确得到目标的辐射亮度或强度。
标定过程首先将红外探测系统参数固定,并依次对不同温度的黑体进行探测,得到其在不同温度条件下的图像输出值(灰度)D,通过式(1)计算其辐射亮度响应度α和自身偏移值D0之间的对应关系,标定模型如下[13]:
$$ D = \alpha L + {D_0} $$ (1) 式中L为黑体在红外探测系统工作波段内的辐射亮度,不同温度对应的黑体辐射亮度通过式(2)计算得到[14]。
$$ L = \frac{\varepsilon }{\pi }\int\limits_{{\lambda _1}}^{{\lambda _2}} {\frac{{{{{c}}_1}}}{{{\lambda ^5}}} \cdot \frac{1}{{{{\mathrm{e}}^{{{{{{c}}_2}} / {\lambda T}}}} - 1}}{\rm{d}}\lambda } $$ (2) 式中:λ1~λ2为红外探测系统工作波段;ε为黑体发射率;c1、c2为辐射常数;T为黑体温度。
目前,针对红外探测系统的标定方法主要有近距离扩展源法、远距离小源法和平行光管法[15]。
近距离扩展源法一般采用大面源黑体作为标准红外辐射源,要求黑体辐射面覆盖红外探测系统光学镜头,并靠近红外探测系统光学镜头,以减小外界环境及背景的影响,该方法可以对红外探测系统探测器的每个像元进行标定。
当黑体辐射面不能覆盖红外探测系统光学镜头,或受现场条件限制,黑体无法靠近红外探测系统光学镜头时,可采用远距离小源法进行标定。远距离小源法一般将黑体置于距红外探测系统一定距离处,并保证红外探测系统能够对黑体清晰成像。
平行光管法需要使用大口径平行光管,标定过程首先将红外探测系统与平行光管对接,其次将黑体置于平行光管焦平面,形成无穷远辐射源。该方法由于标定系统较为复杂,一般不在外场使用。
本文主要探讨外场测量设备的标定,因此主要对近距离扩展源法和远距离小源法进行分析,不涉及平行光管法。
2 标定及结果
2.1 设备及参数
被标定设备为某长波红外告警设备,工作波段为8 μm~12 μm,像元尺寸为25 μm×25 μm,量化位数为16位。标准红外辐射源为以色列CI Systems公司生产的SR800扩展面黑体(见图1),参数见表1。
表 1 SR800扩展面黑体参数Table 1. Parameters of SR800 extended surface blackbody参 数 参数值 温度范围/℃ −40~170 温度精度/℃ ±0.02 辐射面积/mm2 300×300 发射率 0.97 标定时环境温度为28℃左右,为保证标定有效且标定温度分布在环境温度两侧,将标定温度设置在10℃~50℃之间。SR800扩展面黑体温度依次设置为10℃、20℃、30℃、40℃、50℃,利用长波红外告警设备对黑体进行测量,由式(1)、式(2)得到其辐射亮度响应度α和偏移值D0之间的对应关系。
2.2 近距离直接扩展源法标定结果
由于该长波红外告警设备采取旋转扫描方式对目标成像,标定过程中黑体无法紧贴其光学镜头,难以完全排除外界环境及背景的影响。近距离直接扩展源法的标定结果见图2,其中横坐标为黑体辐射亮度L,纵坐标为长波红外告警设备标定区域图像中心10×10个像元的输出灰度平均值。
对近距离直接扩展源法的标定结果进行线性拟合,得到长波红外告警设备的辐射亮度响应关系为D=225.07×L+
5973.7 。2.3 远距离小源法标定结果
标定时由于已将长波红外告警设备参数固定,距离太近无法对黑体清晰成像,距离太远大气及背景辐射对标定影响较大。因此,选取37 m和63 m两个能够清晰成像且距离较近的位置进行标定,37 m处黑体红外图像见图3。
标定结果如图4所示,其中横坐标为黑体辐射亮度L,纵坐标为长波红外告警设备标定区域图像中心10×10个像元的输出灰度平均值。
对远距离小源法的标定结果进行线性拟合,得到长波红外告警设备在37 m处的辐射亮度响应关系为D=198.63×L+
6031.9 ,在63 m处的辐射亮度响应关系为D=181.90×L+5701.2 。3 结果分析
从式(1)可以看出,在长波红外告警设备测量范围内,理论上2个不同的温度点就可以确定其辐射亮度响应度α和偏移值D0。表2~表4分别给出了由相邻温度点得到的近距离扩展源法、远距离小源法标定结果。
表 2 近距离直接扩展源法两点标定结果Table 2. Two points calibration results of short-range direct extended source method标定温度/℃ 10,20 20,30 30,40 40,50 响应度$\alpha $ 232.3 231.5 218.4 222.2 偏移值D0 5873.8 5885.0 6112.5 6034.1 表 3 远距离小源法(37 m)两点标定结果Table 3. Two points calibration results of long-distance small source method (37 m)标定温度/℃ 10,20 20,30 30,40 40,50 响应度$\alpha $ 193.0 206.3 194.4 199.2 偏移值D0 6100.5 5910.6 6117.7 6018.1 表 4 远距离小源法(63 m)两点标定结果Table 4. Two points calibration results of long-distance small source method (63 m)标定温度/℃ 10,20 20,30 30,40 40,50 响应度$\alpha $ 179.9 176.8 193.9 173.2 偏移值D0 5727.7 5772.1 5475.3 5907.5 表2中相邻温度点的辐射亮度响应度$\alpha $和偏移值D0与2.2节中拟合结果的相对误差分别在−2.96%~3.21%和−1.67%~2.32%之间。表3中相邻温度点的辐射亮度响应度$\alpha $和偏移值D0与2.3节中拟合结果的相对误差分别在−2.83%~3.86%和−2.01%~1.42%之间。表4中相邻温度点的辐射亮度响应度$\alpha $和偏移值D0与2.3节中拟合结果的相对误差分别在−4.78%~6.60%和−3.96%~3.62%之间。可以看出,近距离扩展源法的两点标定结果一致性较好,远距离小源法在37 m处的两点标定结果一致性也较好,但63 m处的两点标定结果一致性有所下降,说明随着距离的增加,外界环境及背景的影响有所增强。
在不考虑大气辐射传输的情况下,本次标定由设备带来的误差主要包括:SR800扩展面黑体光谱辐射出射度误差,长波红外告警设备光学系统透过率不确定度引起的误差,输出信号引起的误差,以及探测器响应不稳定引起的误差。SR800扩展面黑体温度精度为±0.02,结合文献[8,15-16]对黑体辐射出射度误差2%左右的分析,SR800扩展面黑体引起的测量误差${\sigma _1} $大约为2%。长波红外告警设备光学系统的光谱透过率实测误差${\sigma _2} $约为2%。根据当前长波红外告警设备技术状态,结合文献[15-16]的分析,其信号输出值D的不确定度${\sigma _3} $优于1%,探测器响应不稳定引起的误差${\sigma _4} $约为2%。
误差的合成标准不确定度可以表示为
$$ {\sigma }_{总}={\left({\sigma }_{1}^{2}+{\sigma }_{2}^{2}+{\sigma }_{3}^{2}+{\sigma }_{4}^{2}\right)}^{1/2} $$ (3) 经计算,由设备带来的总误差约为3.6%。
下面探讨不同标定方法带来的误差。为了验证2种标定结果的准确性,将SR800扩展面黑体作为目标布设于距长波红外告警设备140 m处进行辐射特性测量,利用式(4)对SR800扩展面黑体辐射亮度进行反演。
$$ {L_{\mathrm{F}}} = {{\left( {\frac{{D - {D_0}}}{\alpha } - {L_{{\mathrm{path}}}}} \right)} / {{\tau _\alpha }}} $$ (4) 式中:LF为反演的SR800扩展面黑体辐射亮度;τα为SR800扩展面黑体与长波红外告警设备之间的大气透过率;Lpath为SR800扩展面黑体与长波红外告警设备之间的大气程辐射,其结果通过大气辐射传输计算软件MODTRAN计算得到。测量期间,大气能见度为14.7 km,湿度为39%,温度为28℃,通过MODTRAN计算得到0 m、37 m、63 m距离的大气透过率分别为1、
0.9495 、0.9282 ;大气程辐射分别为0、0.5484 W∙m−2∙sr−1、0.7851 W∙m−2∙sr−1。表5~表7分别给出了采用近距离扩展源法、远距离小源法标定后的测量误差。表 5 近距离直接扩展源法测量误差Table 5. Measurement errors of short-range direct extended source method温度/℃ 黑体辐射亮度/
W∙m−2∙sr−1反演的辐射亮度/
W∙m−2∙sr−1相对
误差/%10 11.63 10.97 −5.67 20 14.30 13.00 −9.09 30 17.36 15.36 −11.52 40 20.81 18.08 −13.12 50 24.67 20.79 −15.73 表 6 远距离小源法(37 m)测量误差Table 6. Measurement errors of long-distance small source method (37 m)温度/℃ 黑体辐射亮度/
W∙m−2∙sr−1反演的辐射亮度/
W∙m−2∙sr−1相对
误差/%10 11.63 12.18 4.75 20 14.30 14.60 2.07 30 17.36 17.12 −1.38 40 20.81 20.19 −2.98 50 24.67 23.26 −5.72 表 7 远距离小源法(63 m)测量误差Table 7. Measurement errors of long-distance small source method (63 m)温度/℃ 黑体辐射亮度/
W∙m−2∙sr−1反演的辐射亮度/
W∙m−2∙sr−1相对
误差/%10 11.63 15.08 29.66 20 14.30 17.59 23.01 30 17.36 20.51 18.15 40 20.81 23.87 14.70 50 24.67 27.22 10.34 $$ \delta = \left( {\frac{{{L_{\mathrm{F}}} - {L_{\mathrm{H}}}}}{{{L_{\mathrm{H}}}}}} \right) \times 100{\text{%}} $$ (5) 式中LH为SR800扩展面黑体实际辐射亮度。
图5给出了采用近距离扩展源法、远距离小源法标定后的测量误差曲线。
1) 目标温度在10℃~50℃范围内,采用近距离直接扩展源法标定的测量误差在5.67%~15.73%之间。采用远距离小源法在37 m处标定的测量误差在1.38%~5.72%之间,在63 m处标定的测量误差在10.34%~29.66%之间。
2) 由于长波红外告警设备以周扫方式工作,无法对目标凝视成像,采用近距离直接扩展源法标定时,无法将黑体紧靠其光学系统,导致部分外界环境及背景辐射进入长波红外告警设备光学系统。另外,长波红外告警设备光学系统焦距较长,周扫过程中,由于角速度大对近距离放置的黑体探测稳定度较差。除3.6%的设备误差外,根据式(6),标定过程中环境及背景辐射、标定方法造成的误差可能在4.38%~15.31%之间。
$$ {\sigma }_{1}={\left({\sigma }_{总}^{2}\text-{\sigma }_{2}^{2}\right)}^{1/2} $$ (6) 3) 采用远距离小源法对长波红外告警设备进行标定,除3.6%的设备误差外,根据式(6),在37 m处,标定过程中环境及背景辐射、标定方法造成的误差优于4.45%,在63 m处,标定过程中环境及背景辐射、标定方法造成的误差可能在9.69%~29.44%之间。37 m处的标定误差明显优于63 m处的标定误差,由此可见,对于采取周扫方式工作的长波红外告警设备,远距离小源法是对其标定的有效手段。从验证结果可以看出,在保证长波红外告警设备能对黑体清晰成像的前提下,降低标定距离可有效减小环境及背景辐射对标定误差的影响。因此,在实际操作过程中,针对设备特点,找出其能对黑体清晰成像的最短距离,是远距离小源法标定的关键。
4 结论
本文利用长波红外告警设备、SR800扩展面黑体进行了外场标定实验,并验证了近距离直接扩展源法、远距离小源法的标定误差。结果表明,对于采取周扫方式工作的长波红外告警设备,近距离直接扩展源法误差较大,该方法适用性较差。远距离小源法是对周扫方式工作的长波红外告警设备进行标定的有效手段,但如何合理选择标定距离以提高其标定精度,是需要重点关注的问题。针对长波红外告警设备参数及特点,找出标定距离与标定误差的对应关系,总结规律并给出最优的标定方案,是未来研究的主要方向。
-
表 1 光学设计指标
Table 1 Optical design indexes
参数 指标 工作波段/nm 250~340 焦距/mm 25~250 F 数 4 CCD 像元尺寸/μm 4.65 有效靶面尺寸/pixels 1 024× 1024 MTF@107 lp/mm ≥0.3 畸变/% ≤1 总长/mm 431.5 变倍比 10 表 2 以衍射MTF子午为标准的蒙特卡罗公差分析结果
Table 2 Monte Carlo tolerance analysis results based on diffraction MTF meridian
蒙特卡罗
样本百分比%焦距/mm 25 133 250 $ \geqslant 90$ 0.273 3 0.301 0 0.148 9 $ \geqslant 80$ 0.288 8 0.362 1 0.187 6 $ \geqslant 50$ 0.324 8 0.452 1 0.275 7 $ \geqslant 20$ 0.359 0 0.544 2 0.367 9 $ \geqslant 10$ 0.372 4 0.594 2 0.422 8 表 3 以衍射MTF弧矢为标准的蒙特卡罗公差分析结果
Table 3 Monte Carlo tolerance analysis results based on diffraction MTF sagitta of arc
蒙特卡罗样
本百分比%焦距/mm 25 147.3 250 $ \geqslant 90$ 0.411 9 0.341 3 0.143 7 $ \geqslant 80$ 0.439 0 0.397 4 0.212 6 $ \geqslant 50$ 0.494 1 0.499 0 0.335 7 $ \geqslant 20$ 0.541 1 0.586 0 0.444 7 $ \geqslant 10$ 0.567 7 0.694 2 0.512 8 表 4 零件允许公差
Table 4 Allowable tolerances of parts
公差 范围 半径公差/光圈数 $ \pm $2 厚度公差/mm $ \pm $0.02 表面偏心公差/mm $ \pm $ 0.0075 表面倾斜公差/(°) $ \pm $ 0.0083 元件偏心公差/mm $ \pm $ 0.0075 元件倾斜公差/(°) $ \pm $ 0.0083 表面不规则度 $ \pm $0.2 折射率公差 $ \pm $ 0.0005 阿贝常数公差/% $ \pm $0.03 -
[1] 蒋斌. 环境空气中二氧化硫监测技术的发展[J]. 环境科学与管理, 2012, 37(1): 158-160. doi: 10.3969/j.issn.1673-1212.2012.01.039 JIANG Bin. Development of monitoring technology of sulfur dioxide in the air[J]. Environmental Science and Management, 2012, 37(1): 158-160. doi: 10.3969/j.issn.1673-1212.2012.01.039
[2] 张宇. 浅析大气污染现状及防治措施[J]. 清洗世界, 2021, 37(7): 90-91. ZHANG Yu. Analysis on the present situation of air pollution and its prevention measures[J]. Cleaning World, 2021, 37(7): 90-91.
[3] 张玉杰. 中国二氧化硫排放情况及其对环境的影响(待续)[J]. 硫磷设计与粉体工程, 2021(4): 30-35. ZHANG Yujie. Sulfur dioxide emission in China and its impact on the environment (to be continued)[J]. Sulfur and Phosphorus Design and Powder Engineering, 2021(4): 30-35.
[4] 张玉杰. 中国二氧化硫排放情况及其对环境的影响(续完)[J]. 硫磷设计与粉体工程, 2021(5): 43-48. ZHANG Yujie. Sulfur dioxide emission in China and its impact on the environment [J]. Sulfur and Phosphorus Design and Powder Engineering, 2021 (5): 43-48.
[5] 杜晓峰. 大气污染原因和环境监测治理技术研究[J]. 科技与创新, 2020(23): 135-136. DU Xiaofeng. Study on air pollution causes and environmental monitoring and control technology[J]. Science and Technology and Innovation, 2020(23): 135-136.
[6] 吴淑枞, 陈祖云. 工业废气中二氧化硫的治理[J]. 当代化工研究, 2022(13): 58-60. doi: 10.3969/j.issn.1672-8114.2022.13.021 WU Shucong, CHEN Zuyun. Treatment of sulfur dioxide in industrial waste gas[J]. Contemporary Chemical Research, 2022(13): 58-60. doi: 10.3969/j.issn.1672-8114.2022.13.021
[7] DALTON M P, WATSON I M, NADEAU P A, et al. Assessment of the UV camera sulfur dioxide retrieval for point source plumes[J]. Journal of Volcanology and Geothermal Research, 2009, 188(4): 358-366. doi: 10.1016/j.jvolgeores.2009.09.013
[8] 严韦, 任鑫, 刘建军, 等. 嫦娥三号极紫外相机几何定位方法及其应用[J]. 光子学报, 2016, 45(7): 119-124. YAN Wei, REN Xin, LIU Jianjun, et al. Geometric positioning method of Chang'e-3 extreme ultraviolet camera and its application[J]. Acta Photonica Sinica, 2016, 45(7): 119-124.
[9] 张英华, 李昂, 谢品华, 等. 污染气体浓度二维空间分布的紫外成像方法[J]. 光谱学与光谱分析, 2018, 38(5): 1476-1480. ZHANG Yinghua, LI Ang, XIE Pinhua, et al. Ultraviolet imaging method of two-dimensional spatial distribution of pollution gas concentration[J]. Spectroscopy and Spectral Analysis, 2018, 38(5): 1476-1480.
[10] VARNAM M, BURTON M, ESSE B, et al. Two independent light dilution corrections for the SO2 camera retrieve comparable emission rates at masaya volcano, nicaragua[J]. Remote Sensing, 2021, 13(5): 935. doi: 10.3390/rs13050935
[11] 段为民, 熊远辉, 陈振威, 等. 工业 SO2及碳黑颗粒紫外成像遥感监测技术[J]. 光子学报, 2020, 49(4): 153-161. DUAN Weimin, XIONG Yuanhui, CHEN Zhenwei, et al. Remote sensing monitoring technology of industrial SO2 and carbon black particles by ultraviolet imaging[J]. Acta Photonica Sinica, 2020, 49(4): 153-161.
[12] 熊远辉, 罗中杰, 陈振威, 等. SO2气体排放的紫外成像遥感监测[J]. 光谱学与光谱分析, 2020, 40(4): 1289-1296. XIONG Yuanhui, LUO Zhongjie, CHEN Zhenwei, et al. Monitoring SO2 emission by ultraviolet imaging remote sensing[J]. Spectroscopy and Spectral Analysis, 2020, 40(4): 1289-1296.
[13] 李世江. 8 mm~32 mm 双波段变焦安防镜头设计[D]. 西安: 西安工业大学, 2022. LI Shijiang. Design of 8 mm~32 mm dual-band zoom security lens[D]. Xi'an: Xi'an University of Technology, 2022.
[14] 高子英. 大口径变焦投影系统结构设计[D]. 长春: 长春理工大学, 2017. GAO Ziying. Structural design of large-aperture zoom projection system [D]. Changchun: Changchun Univer-sity of Science and Technology, 2017.
[15] 宋鹏飞. 10 倍变焦距镜头设计[D]. 长春: 长春理工大学, 2022. SONG Pengfei. Design of 10x zoom lens[D]. Changchun: Changchun University of Science and Technology, 2022.
[16] 曾振煌, 林峰. 30 mm~110 mm大孔径红外变焦热像镜头设计[J]. 应用光学, 2016, 37(5): 752-758. ZENG Zhenhuang, LIN Feng. Design of 30 mm~110 mm large aperture infrared zoom thermal image lens[J]. Journal of Applied Optics, 2016, 37(5): 752-758.
[17] 罗锐, 梁秀玲. 大孔径宽光谱变焦镜头设计[J]. 应用光学, 2022, 43(5): 839-845. doi: 10.5768/JAO202243.0501002 LUO Rui, LIANG Xiuling. Design of large aperture and wide spectrum zoom lens[J]. Journal of Applied Optics, 2022, 43(5): 839-845. doi: 10.5768/JAO202243.0501002