Simulation and design of in-line fiber optic etalon with Zemax
-
摘要:
从原理上分析了单层膜和法布里-珀罗谐振腔(Fabry-Parot cavity,FP腔),并使用Zemax软件的膜层设计工具对分光比膜和FP腔进行了建模,分析了不同折射率材料单层膜的反射率和透射率随波长的变化关系。对FP腔的模型进行简化分析,用3 mm间距的空气隙FP腔进行了反射率和透射率初步验证。设计了一款通讯常用的50 GHz的光学标准具,并介绍了一种基于自聚焦透镜准直器的在线光学标准具生产方案。通过设计,控制自聚焦透镜平面分光比膜的透射率和反射率,调整2个自聚焦透镜准直器之间的距离,达到控制FP腔的峰值、谷值透射率和自由光谱范围(free spectral range,FSR)等参数的目标。通过Zemax建模,发现模拟和样品的参数值基本一致。分析了样品间差异、光谱峰值透射率的差异和光谱相位差异。使用Zemax分析FP标准具方便直观,为FP的设计与分析提供了一种新方法。
Abstract:The single-layer coating and Fabry-Parot resonant cavity (FP cavity) were analyzed in principle, and the tap coating and FP cavity were modeled using the coatings tool of Zemax software. The reflectance and transmittance of single-layer coating with different refractive index materials were analyzed with respect to wavelength. The model of FP cavity was simplified and analyzed, and the reflectance and transmittance were preliminarily verified with an air-gap FP cavity with 3 mm gap. A fiber optical etalon of 50 GHz commonly used in communication was designed, and a production scheme of in-line optical etalon based on Glens collimator was introduced. Through the design, the transmittance and reflectance of the spectral ratio film of the Glens plane were controlled, and the distance between two Glens collimators was adjusted to achieve the control of the peak value, valley value transmittance and free spectral range (FSR) parameters of the FP cavity. Through the simulation, the parameters of the simulation and samples were obtained to be basically the same. The inter-sample differences, differences in spectral peak transmittance and spectral phase differences were also analyzed. Analysis of FP etalon by Zemax is convenient and intuitive, which provides a new method for FP design and analysis.
-
Keywords:
- Fabry-Parot etalon /
- Zemax /
- optical simulation
-
引言
飞机类目标作为现代战争的重要平台,广泛应用于侦察、攻击、运输等领域。基于此,其也成为红外探测系统、红外制导武器探测和打击的主要对象。准确获取来袭飞机的红外辐射特性,对于红外预警探测具有重要意义[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 Reflectance of different refractive indexes
n1 R/% 1.000 0.0 1.387 10.0 1.618 20.0 1.850 30.0 2.107 40.0 2.414 50.0 2.806 60.0 3.353 70.0 4.236 80.0 6.162 90.0 表 2 膜层FP的参数
Table 2 Parameters of coating FP
Material Thickness Is_absolute Loop_index A 0.25λ 0 0 Air 3000 μm1 0 A 0.25λ 0 0 表 3 样品和模拟参数
Table 3 Parameters of samples and simulation
Parameters Sample1 Sample2 Sample3 Sample4 Sample5 Sample6 Simulation fFSR/GHz 49.5 49.2 49.6 49.5 49.9 50.0 50.0 Peak loss/% 95.4 97.0 96.0 93.2 93.3 92.1 100.0 Bottom loss/% 11.3 11.3 11.4 11.2 9.5 9.2 11.1 -
[1] 余佳音, 樊静, 蓝旭辉, 等. 反射相移在Fabry-Perot标准具间距测量中的影响[J]. 激光与光电子学进展, 2020, 57(9): 091201. YU Jiayin, FAN Jing, LAN Xuhu, et al. Influence of reflection-induced retardance on the measurement of Fabry-Perot etalon interval[J]. Laser & Optoelectronics Progress, 2020, 57(9): 091201.
[2] 刘加庆, 刘磊, 刘雷, 等. 布里渊光谱仪的高精度波长标定方法研究[J]. 光学学报, 2020, 40(20): 2030001. doi: 10.3788/AOS202040.2030001 LIU Jiaqing, LIU Lei, LIU Lei, et al. Wavelength calibration of ultra-high resolution brillouin spectrometer[J]. Acta Optica Sinica, 2020, 40(20): 2030001. doi: 10.3788/AOS202040.2030001
[3] 李路, 庄鹏, 谢晨波, 等. 多普勒测风激光雷达的FP 标准具透过率曲线校准及风场观测[J]. 光子学报, 2020, 49(11): 1149016. LI Lu, ZHUANG Peng, XIE Chenbo, et al. Calibration of FP etalon transmittance curve and wind field observation based on Doppler wind measurement lidar[J]. Acta Photonica Sinica, 2020, 49(11): 1149016.
[4] 徐琼, 丁丽萍, 付号, 等. FP腔中驱动微悬臂的双稳效应的研究[J]. 光学与光电技术, 2012, 10(6): 49-52. XU Qiong, DING Liping, FU Hao, et al. Bistability of optical actuated micro-cantilever in a Fabry-Perot[J]. Optics & Optoelectronic Technology, 2012, 10(6): 49-52.
[5] 刘晓明, 张敏明, 刘德明. 基于法布里-珀罗滤波片的宽光谱可调谐滤波器[J]. 光学与光电技术, 2006, 4(1): 55-57. LIU Xiaoming, ZHANG Minming, LIU Deming. Widely tunable Febry-Perot etalon filter[J]. Optics & Optoelectronic Technology, 2006, 4(1): 55-57.
[6] 杨四刚, 潘晴, 胡必春, 等. 基于固体腔F-P 标准具的可调谐自适应光滤波器[J]. 半导体光电, 2003, 24(2): 107-109. YANG Sigang, PAN Qing, HU Bichun, et al. An adaptive tunable optical filter based on solid cavity F-P etalon[J]. Semiconductor Optoelectronics, 2003, 24(2): 107-109.
[7] 金文最, 赵洪霞, 安雪瑶, 等. 通信波段宽带吸收器设计研究[J]. 光电子·激光, 2020, 31(8): 800-805. JIN Wenzui, ZHAO Hongxia, AN Xueyao, et al. Research on the optimal design of a communication band broad-band absorber[J]. Journal of Optoelectronics·Laser, 2020, 31(8): 800-805.
[8] 江孝伟, 武华. 双通道窄带宽超材料吸收器[J]. 光学学报, 2021, 41(14): 1416002. doi: 10.3788/AOS202141.1416002 JIANG Xiaowei, WU Hua. Dual-channel narrow bandwidth metamaterial absorber[J]. Acta Optica Sinica, 2021, 41(14): 1416002. doi: 10.3788/AOS202141.1416002
[9] 胡文杰, 梅文婷, 黄峻峰, 等. 光栅外腔提高VCSEL微波调制效率实验研究[J]. 光电子·激光, 2021, 32(11): 1135-1139. HU Wenjie, MEI Wenting, HUANG Junfeng, et al. Experimental study on enhancement of microwave modulation efficiency in VCSEL laser using a grating external cavity[J]. Journal of Optoelectronics·Laser, 2021, 32(11): 1135-1139.
[10] 闫亚东, 何俊华, 仓玉萍, 等. 用ZEMAX模拟五棱镜误差对平行度检测的影响[J]. 应用光学, 2007, 28(5): 649-653. YAN Yadong, HE Junhua, CANG Yuping. Effect of penta prism error on parallelism detection with ZEMAX[J]. Journal of Applied Optics, 2007, 28(5): 649-653.
[11] 程义. 基于Zemax的C-lens光纤准直器点精度分析[J]. 机电工程技术, 2022, 51(6): 143-146. CHENG Yi. Analysis of point error for C-lens fiber collimator with Zemax[J]. Mechanical & Electrical Engineering Technology, 2022, 51(6): 143-146.
[12] 胡永宏. 用于成像的FP标准具研究[D]. 长春: 长春理工大学, 2010. HU Yonghong. Reserch on FP etalon for imaging[D]. Changchun: Changchun University of Science and Technology, 2010.
[13] 唐晋发, 顾培夫, 刘旭, 等. 现代光学薄膜技术[M]. 1版. 杭州: 浙江大学出版社, 2006: 61-64. TANG Jinfa, GU Peifu, LIU Xu, et al. Modern optical thin film technology[M]. 1st ed. Hangzhou: Zhejiang University Press, 2006: 61-64.
[14] 马科斯·玻恩, 埃米尔·沃尔夫. 光学原理: 光的传输、干涉和衍射的电磁理论[M]. 杨葭荪, 译. 7版. 北京: 电子工业出版社, 2009: 299-340. BORN M, WOLF E. Principles of optics: electromagnetic theory of propagation, interference and diffraction of light[M]. YANG J S, Transl. 7th ed. Beijing: Electronic Industry Press, 2009: 299-340.
[15] 程义. 一种基于FP标准具的光滤波器: CN212009020U [P]. 2020-11-24. CHENG Yi. An optical filter based on FP etalon: CN212009020U [P]. 2020-11-24.