摘要:提出了一种基于动态光学频率梳啁啾脉冲干涉的绝对距离快速测量方法. 借助于重复频率的线性扫描, 可获取啁啾光谱干涉信号中最宽条纹的动态频率偏移, 从而完成被测距离的高精度测量. 动态重复频率能够延伸光谱仪探测极限, 相应地减小测量盲区, 并且基于合成波长也使得测量非模糊范围得到极大的拓展. 本文测距系统无需重复频率的锁定, 能够摆脱对锁相环等复杂应用模式的依赖, 在简化系统的同时借助铷钟提供精准时钟参考, 从而进一步提高测量精度. 此外, 电荷耦合器件成像帧速以及重复频率扫描速度的提升可实现干涉信号的快速采集, 弥补传统光谱干涉测量以及腔调谐方式在探测速度上的不足. 实验结果表明, 本文方法单点测距数据更新率为13.5 Hz, 相较于参考数值, 在20 m的测量范围内, 测量不确定度优于27 μm, 相对精度为1.35 × 10^–6.
关键词: 绝对距离测量/
动态光学频率梳/
啁啾脉冲干涉

English Abstract

全文HTML

--> --> -->

光学频率梳自诞生之初, 就在众多应用领域的发展中展现出划时代意义, 并体现重要的研究价值, 也因此受到国内外专家****的广泛关注^[1-3]. 其中光频梳在频率计量^[4,5]、光谱分析^[6,7]、太空梳^[8,9]和绝对距离测量^[10,11]等领域中, 依据其自身特性, 打破了传统激光应用的技术壁垒, 为各领域的方法革新提供了特有的解决方案. 光学频率梳在时域上表现为锁模脉冲序列, 在频域中具有频率间隔严格相等的频率纵模, 并保持极强的相干性. 光学频率梳其内在从射频到光频的频率纵模可溯源至相同的时钟基准, 这也意味着光频梳自身的稳定性和精度能够借助外部参考得以提升.
在光学频率梳众多应用领域中, 飞秒光学频率梳的绝对距离测量技术得到了迅猛发展, 并且借助于长度这个基本物理量所衍生出的相关计量技术也成为国内外课题组的热点研究方向. 目前基于光频梳的时频特性已提出多种较为成熟的绝对距离测量技术, 其中飞行时间测量基于脉冲互相关信号, 精确得出测量脉冲与参考脉冲的时间延迟, 将测量精度从mm级提升至nm级^[12]. 基于长延迟光纤的谐振腔扫描光学采样(OSCAT)技术在保证测距精度的同时, 可实现探测范围内的任意距离测量, 从而消除测量盲区^[13,14]. 此外, 光学频率梳的多波长相位探测相较于传统激光相位测量, 在保证亚波长精度的基础上极大地拓展了非模糊范围^[15], 另外基于光频梳模间拍频, 通过引入高次谐波可显著提高绝对距离测量精度^[16,17]. 对于依据光频及啁啾信息的色散干涉, 其在光纤色散测量和相干层析成像等领域的应用证明了该方法优秀的计量特性^[18-20]. 借助光学频率梳的频率溯源性, 色散干涉可提取解卷裹的光谱相位斜率^[21], 或探测啁啾脉冲的频率偏移完成绝对距离测量^[22]. 系统中无需运动机械结构即可获取包含距离信息的光谱干涉条纹, 这也使得色散干涉具备实现快速测量的理论可能. 基于异步光学采样的双光梳技术, 通过两台具有微小重频差异的光学频率梳, 获取双光梳干涉条纹以实现绝对距离测量. 双光梳测距系统在保证干涉测量精度的同时, 具备极高的数据更新率, 并且通过重复频率的调整也可得到极大的非模糊范围^[23-25].
应对激光雷达以及光学遥感等大尺寸应用的工程实践需要, 本文提出了基于动态啁啾脉冲干涉的绝对距离测量方法. 在重复频率动态扫描的过程中, 系统实时高速采集对应重频的啁啾脉冲干涉条纹, 从中依据中心频率的动态偏移量实现被测距离的解算. 基于动态重频的合成波长, 该方法极大地拓展了测量非模糊范围, 并且保证了方法的测距精度和测量速度. 此外, 本文方法在一定程度上放宽光谱仪分辨率的限制, 延伸啁啾脉冲干涉的可测范围, 依据动态频率梳参数, 当被测范围大于~250 m, 可实现任意绝对距离测量. 而且动态频率梳包含距离测量结果的动态信息, 其可消除无法直接分辨光谱仪处重合脉冲相对时延的局限性, 从而实现绝对距离的唯一确定. 测量结果以及长时间测距实验表明本文方法能够实现被测距离的高精度快速测量, 并且验证了方法的测量稳定性.

基于动态啁啾脉冲干涉的实验原理如图1所示. 动态光学频率梳的实现基于环形腔长的线性扫描, 环形腔中的精密机械位移台通过高定位精度完成腔长的周期性调制. 此外, 整个光源系统处于隔温隔震的密闭环境, 保证了整体系统的稳定性. 测距干涉仪由图1中的分束镜(BS)分为两部分, 其中参考路光程固定并且通过一对光栅实现脉冲啁啾, 另外一路行进经过待测距离称为测量路. 两路光束通过回射系统反射, 回光在分束镜处合束, 随后借助光谱探测系统进行高速采集.
图 1 动态啁啾脉冲干涉实验原理图
Figure1. Schematic diagram of dynamic chirped pulse interferometry.

由于重复频率(f_rep)的快速扫描, 无论是传统的空间色散光谱仪还是傅里叶变换光谱仪均受限于其内部存在的机械组件, 因此无法在高分辨率模式下实现动态啁啾干涉信号的快速获取. 在本文方法中, 干涉信号通过光栅色散, 借助电荷耦合器件(CCD)实现高帧速采集, 弥补了传统光谱仪探测速度的不足^[26,27].
假设光频梳发出理想的高斯脉冲, 其光谱为理想高斯光谱, 则系统中引入啁啾的参考脉冲电场可表示为

$ {E_{{\text{ref}}}}(t) = \alpha {E_0}\exp \left( { - {a_1}{t^2}} \right)\exp \left[ { - {\text{i}}\left( {{\omega _{\text{c}}} + bt} \right)t} \right] , $

相对应的测量脉冲可表示为

$ {E_{{\text{mea}}}}(t) = \beta {E_0}\exp \big[ { - {a_2}{{\left( {t - \tau } \right)}^2}} \big]\exp \left[ { - {\text{i}}{\omega _{\text{c}}}\left( {t - \tau } \right)} \right] , $

其中α和β为功率因子, $ a_1=2?\ln(2/t_c^2) $, t_c为啁啾后的脉冲宽度, b为啁啾系数, $ b\approx \Delta \omega/(2t_c) $, Δω为脉冲的带宽, $ a_2=2\ln(2/t_0^2) $, t₀为未产生啁啾的脉冲宽度, 上述各项分别为表征系统光栅、光源自身性能的参数. 此外, ω_c为中心角频率, τ为测量脉冲与参考脉冲间的飞行时间延迟, 其指代系统所包含的被测距离信息. 则脉冲光谱干涉结果可表述为

$\begin{split}I(\omega ) =\;& {{\left| {\Im \left\{ {{E_{{\rm{ref}}}}(t)} \right\} + \Im \left\{ {{E_{{\rm{mea}}}}(t)} \right\}} \right|}^2}\\ =\;& {I_0}(\omega ) + {I_1}(\omega ){\rm{cos}}\left( {\varphi (\omega )} \right).\end{split}$

式中, $\Im $代表傅里叶变换; I₀为两脉冲光谱的叠加, 可看作直流项; 而另一项为交流项, 其中I₁为调制幅值. 则啁啾脉冲的光谱干涉信号相位?(ω)可表示为

$ \phi (\omega ) = \frac{{b{{(\omega - {\omega _c})}^2}}}{{4({a_1}^2 + {b^2})}} - \tau \omega - \frac{1}{2}\arg ({a_1} + {\text{i}}b) . $

(4)式表明, 引入脉冲啁啾的相位函数与角频率为二次函数关系, 即光谱干涉条纹不再随着角频率的变化形成稳定振荡, 而是伴随着一个变化的调制频率. 当调制频率达到极小值时, 干涉条纹振荡最为缓慢且同时出现最宽振荡条纹. 则对(4)式求导可得

$ \frac{{\partial \phi (\omega )}}{{\partial \omega }} = \frac{{b(\omega - {\omega _{\text{c}}})}}{{2({a_1}^2 + {b^2})}} - \tau = 0 , $

其中飞行时延τ为被测距离在光速定义下的时域映射, 本文中定义c为真空中光的传播速度, n_g为脉冲在相应介质中的群折射率, l_d为动态啁啾脉冲干涉条纹所对应的被测距离. 因此基于对应关系τ = 2n_gl_d/c, 调制频率极小处对应的频率可表示为

$ \begin{split} {\omega _{\text{L}}} =\;& {\omega _{\text{c}}} + 2b\tau + {{2a_1 ^2\tau }}/{b}\\ \approx\;& {\omega _{\text{c}}} + {{4b{n_{\text{g}}}{l_{\text{d}}}}}/{c},\;{\text{ }}b \gg {a_1}. \end{split} $

依据本文的实验条件, 认为b(10²⁴)远远大于a₁(10¹⁹). 基于啁啾脉冲干涉特性, l_d可理解为距离测量的小数部分. 实际距离L为l_d与N倍脉冲间距l_pp的累加, 即$ L = (N\cdot l_{\rm pp})/2+l_{\rm d} $, 其中l_pp可表示为$ c/(n_{\rm g}\cdot f_{\rm rep}) $. 则实际被测距离L可表示为

$ L = \frac{N}{2} \cdot \frac{c}{{{n_{\text{g}}}{f_{{\text{drep}}}}}} + {l_{\text{d}}} , $

其中, f_drep为动态光频梳中某一时刻的重复频率. 当系统中引入动态光学频率梳时, 重复频率f_drep数值实现周期线性扫描. 由(7)式可知, 在被测距离L一定的前提下, 距离量l_d会随之产生动态变化, 并且光谱干涉条纹中最宽条纹的位置会发生往复偏移. 基于(6)式表征的啁啾频移关系, 距离量l_d可计算为

$ {l_{\text{d}}} = \frac{c}{{4b{n_{\text{g}}}}}({\omega _{{\text{dyn}}}} - {\omega _{\text{c}}}) = \frac{{c{\omega _{{\text{shift}}}}}}{{4{n_{\text{g}}}b}}{{,\;\; b}} \gg {a_1}, $

其中, ω_dyn为动态光频梳在某一重频时最宽条纹对应的频率, 并对应着一个频率偏移量ω_shift. 因此, 通过最宽条纹位置的变化即可解算被测距离L, 同时基于动态重复频率的合成波长能够拓展测量非模糊范围. 则最终实际被测距离L可表示为

$ L = \frac{N}{2} \cdot \frac{c}{{{n_{\text{g}}}{f_{{\text{drep}}}}}} + \frac{{c{\omega _{{\text{shift}}}}}}{{4{n_{\text{g}}}b}}{{, \;\; b}} \gg {a_1} . $

上述理论分析中, 假设光谱仪处测量臂脉冲较于参考臂脉冲存在一定时延, 然而传统色散干涉对于时延的正负存在模糊性, 而动态啁啾脉冲干涉基于条纹动态信息则可去除该时延模糊, 从而确定计算被测距离L时整数倍l_pp和l_d的加减对应关系, 使得绝对距离测量具有唯一性.

本文依据所提出方法的测量原理搭建相应的实验系统, 并在参考路中引入脉冲啁啾. 为验证原理, 并能够对性能参数做出准确分析, 本文进行了大尺寸的绝对距离测量实验. 测量过程中, 环境参数保持相对稳定, 其中温度为21.20 ℃, 气压为997.63 hpa, 相对湿度为45.1%. 依据Ciddor公式^[28], 在此环境条件下的群折射率为1.00026657.
本实验中, 飞秒脉冲激光器基于非线性偏振旋转锁模(NPR), 重复频率103.6 MHz, 动态重频可扫描范围 ± 300 kHz, 相应l_pp调谐范围为2.8846—2.9014 m, 该范围内扫描周期为1.16 s. 此外, 光源中心波长为1030 nm, 光谱宽度约为15 nm, 其光谱如图2所示.
图 2 脉冲激光光谱图
Figure2. Pulse laser spectrogram.

2

-->

3.1.动态啁啾脉冲测距系统标定

由于测量系统中, 频率偏移量、啁啾系数等参数难以做到准确获取, 并且实验条件的不确定性也将提高各项参数的获取难度. 在本文系统中, 基于动态重复频率就可使外部标定系统在无运动组件的情况下完成距离参数标定. 在标定系统中, 当L约为l_pp/2扫描范围内时, 随着重复频率的线性扫描, 最宽条纹位置会发生往复偏移. 基于重频扫描CCD标定原理, 搭建标定系统如图3所示.
图 3 CCD标定系统示意图
Figure3. Sketch of CCD calibration system.

其中商用光谱仪(Thorlabs, OSA201C)用来观察频率偏移量, 另一部分, 激光通过光栅色散, 经准直柱面镜(f = 200 mm)完成线激光准直, 借助CCD相机(PCO.pixelfly usb, 13.5 fps, 1392 × 1040 pixels)采集啁啾脉冲干涉条纹. 并且, 基于采集图像的灰度信息, 通过最宽条纹的插值拟合实现像素的光频校正, 从而获取重频变化对应的像素偏移量. 标定过程分别依次选取两个标定重频为103.631和103.710 MHz时的啁啾脉冲干涉条纹, 实验中使用光阑遮挡降低图像背景亮度, 从而提高条纹对比度. 相机标定结果如图4所示.
图 4 相机标定实验结果图　(a), (d) 两标定重复频率对应光谱图; (b), (e) 两标定重复频率对应相机采集图像; (c), (f) 图像灰度处理结果
Figure4. Experimental results of CCD calibration: (a), (d) Spectrograms for two calibrated repetition rates; (b), (e) captured camera images for two calibrated repetition rates; (c), (f) grayscale images processing results.

实验中受到振动及环境噪声影响, 无法直接获取最宽条纹处对应频率. 因此依据干涉条纹特点, 可获取最宽条纹两侧两个波峰或波谷的像素索引P_left和P_right, 从而确定两峰值处对应的中值位置, 完成最宽条纹处P_widest位置的精确测量. 该解算模式可表示为

$ {P_{{\text{widest}}}} = ({P_{{\text{left}}}} + {P_{{\text{right}}}})/2 . $

在此过程中, 可看作脉冲的空间位置实现了线性扫描, 则标定距离可计算为

$ {L_{{\rm{cal}}}} = \frac{c}{{{n_{\rm{g}}}}}\left( {\frac{{\Delta {f_{{\rm{rep}}}}}}{{{f_{{\rm{drep1}}}} \cdot {f_{{\rm{drep2}}}}}}} \right),$

其中Δf_rep为两个标定重频相应的重频差. 由于系统中啁啾量保持不变, 位移量与最宽条纹移动量保持线性关系, 因此对于单位像素的标定距离可以表示为

$ {L_{{\text{pPix}}}} = \dfrac{L_{{\text{cal}}}}{P_{{\text{widest1}}} - P_{{\text{widest2}}}} . $

基于本系统标定结果, 像素条纹偏移量的标定距离可计算为6.9 × 10^–6 m/pixel. 则某一重频下绝对距离L的计算可重写为

$ L = \frac{1}{2}\left( {N \cdot {l_{{\text{pp}}}} + {P_{{\text{shift}}}} \cdot {L_{{\text{pPix}}}}} \right) , $

式中P_shift为该重复频率条件下, 对应实际被测距离L的最宽条纹像素偏移量. 上述标定方法简化了光谱干涉信号相机标定的方式, 在有效地提高标定速度的同时, 摆脱了传统方式中对位移台定位精度的依赖.
2 -->

3.2.实验结果

如图5所示, 基于Michelson干涉仪结构, 本文将动态啁啾脉冲干涉应用于大尺寸绝对距离测量当中, 并以增量式干涉仪(Agilent 5519B)的测距结果作为参考数值, 其相应的稳定度在20 m处, 10 min的测量时间内优于0.3 μm. 实验中仔细调试了相机及光栅位置, 以确保获取图像的完整和清晰. 此外, 脉冲动态重频的精准获取, 以及对应时刻相机信号的同步采集是决定探测精度的关键步骤. 因此在信号探测过程中, 铷钟(Microsemi 8040)被用作频率计数器(Keysight 53230A)的外部时钟参考, 并且借助信号源(Tektronix AFG3102C)提供极为精准的同步触发信号, 为频率计数器门控和相机的同步采集提供统一的时间基准.
图 5 实验装置图
Figure5. Diagram of experimental setups.

为了能够获取包含距离信息的光谱干涉信号, 实验中, 测量参考点以l_pp为间距均匀分布在光学长导轨上. 其中当动态重复频率为103.374 MHz, 在2.9 和17.4 m处的光谱干涉条纹如图6所示. 随着重复频率的动态扫描, 可直接获得最宽条纹的左右偏移, 从而完成正负时延的确定, 得到唯一确定的被测距离. 由图6可知, 长距离实验易受到环境振动及空气扰动等影响, 且光路的准直等高也会随着被测距离的增大而变化, 对干涉条纹的调制深度产生影响.
图 6 啁啾脉冲干涉实验结果图　(a) 2.9 m处光谱图; (b) 2.9 m处相机采集条纹; (c) 17.4 m处光谱图; (d) 17.4 m处相机采集条纹
Figure6. Experimental results of chirped pulse interferometry: (a) Spectrogram at 2.9 m; (b) captured camera fringes at 2.9 m; (c) spectrogram at 17.4 m; (d) captured camera fringes at 17.4 m.

在光路调试过程中, 测量光路与He-Ne激光干涉仪光路严格对准, 使之平行于光学导轨. 通过步进导轨不断地调试光路, 找出光路的最优状态, 从而抑制系统中阿贝误差和余弦误差对测量结果的影响. 在20 m的测量范围内, 每个测量位置完成8次重复性测量, 其中选取2个重复频率为一组, 共获取4组测量数据. 相较于参考干涉仪, 实验结果的残差和标准差如图7所示. 实验数据表明, 在20 m范围内测量不确定度优于27 μm, 相对精度为1.35 × 10^–6. 相较于双光梳或模间拍频等测距技术, 本文方法摆脱了系统对锁相环的依赖, 但频率稳定性及光谱分辨率则影响着最终的测距精度. 然而, 本文方法机理上与波长信息相连接, 保证了方法的测量精度, 并且借助最宽条纹的中值解算模式, 可有效地抑制环境噪声对测距结果的影响. 此外, 该方法基于动态啁啾信息可实现非模糊范围的快速拓展.
图 7 距离测量结果
Figure7. Results of distance measurements.

为验证本文方法长时间测量的稳定性, 实验在17.4 m处进行了长时间测量, 实验结果残差散点图及阿伦方差如图8所示. 为便于观察, 散点图为1/10的系统抽样结果, 数据标准差为8.4 μm, 采集数据更新率为13.5 Hz, 在1和100 s的平均时间内阿伦方差分别可达到2.3 μm及265 nm.
图 8 长时间距离测量实验结果　(a) 数据散点图; (b) 测距结果阿伦方差
Figure8. Long-term experimental results of the distance measurement: (a) Data scatter diagram; (b) Allan deviation of ranging results.

2 -->

3.3.测量不确定度分析

为对本文方法测量不确定度做出评价, 研究系统各项参数对测距结果的精度影响, (13)式可重写为

$ L = \frac{1}{2} \cdot \frac{c}{{{n_{\rm{g}}}}} \cdot \left( {\frac{N}{{{f_{{\rm{drep}}}}}} + \frac{{{\rm{\pi }} \cdot {f_{{\rm{shift}}}}}}{b}} \right),$

则测量结果的合成不确定度可表示为

$ \begin{split}u_{\rm{L}}^{\rm{2}} =\;& {\left( {\frac{L}{{{n_{\rm{g}}}}} \cdot {u_{{n_{\rm{g}}}}}} \right)^2} + {\left( {\frac{{c \cdot N}}{{2{n_{\rm{g}}} \cdot f_{{\rm{drep}}}^2}} \cdot {u_{{f_{{\rm{drep}}}}}}} \right)^2} \\&+ {\left( {\frac{{c \cdot {\rm{\pi }}}}{{2{n_{\rm{g}}} \cdot b}} \cdot {u_{{f_{{\rm{shift}}}}}}} \right)^2}.\end{split}$

本文方法测量不确定度主要包含群折射率n_g、动态重复频率f_drep、频率偏移量f_shift等影响因素. 首先, 对于群折射率的确定需考虑传感器固有不确定度、环境稳定性以及在测量范围内的环境不均匀性. 基于上述分析, 实验中通过环境传感器进行环境监测和参数采集, 借助Ciddor公式对空气折射率进行实时补偿, 其中温度、大气压以及相对湿度等环境参数波动分别为42.7 mK, 12.9 Pa, 和 3.1%, 对应的测量不确定度分别为4 × 10^–8, 3.4 × 10^–8和2.9 × 10^–8, 则合成不确定度为6 × 10^–8·L. 第二项中, 动态重复频率的采集依靠频率计数器, 并且与CCD的采集触发信号严格锁定至同一时间基准上, 其中动态重频不确定度为5.38 Hz(标准差), 则第二项对应测量不确定度为9 × 10^–7 m. 第三项与最宽条纹频率偏移量相关, 在本文系统中也可对应为像素偏移量的测量不确定度. 通常在实验中, 像素分辨率、像素尺寸、干涉条纹调制深度以及最宽条纹定位精度均会对测量不确定度造成一定的贡献. 基于相机标定结果, 光频与像素存在空间映射关系, 依据重复性测量可得频率偏移量不确定度为0.063 THz, 则该项测量不确定度为26.9 μm. 基于上述分析, 本文方法合成测量不确定度为[(27 μm)² + (6 × 10^–8·L)²]^1/2. 此外, 在啁啾系数一定的情况下, 系统光路经过调试和优化始终保持出色的调制深度, 并且所采用的最宽条纹定位提取算法也能实现频率偏移量的准确测量.

本文提出一种快速绝对距离测量方法, 基于重复频率的线性扫描, 实现动态啁啾脉冲干涉. 本文方法突破了光学频率梳重复频率锁定的传统应用模式, 借助于重频的动态变化极大地延伸了测量非模糊范围, 并且能有效拓展光谱仪探测分辨极限. 在简化系统的同时, 光谱干涉信号也能够实现快速采集, 满足了几何量计量的工程实践需求. 基于实验结果, 数据更新率为13.5 Hz的长时间测量验证了本文方法的测量稳定性, 并且在20 m范围内, 测量不确定度优于27 μm, 相对精度为1.35 × 10^–6. 本文系统在参考臂引入脉冲啁啾, 而在今后的研究中也可通过大气或水体色散在测量臂实现脉冲啁啾, 从而将影响因素转变为测量信息, 以提高方法对实际环境的适用性. 并且, 本文工作可为大尺寸绝对距离测量、面型三维重构、以及水体等强色散介质中的目标探测提供相应的技术储备.