重复频率倍增光频梳时域互相关绝对测距*

梁旭 林嘉睿 吴腾飞 赵晖 邾继贵

(天津大学,精密测试技术及仪器国家重点实验室,天津 300072)

利用法布里-珀罗标准具对光纤光频梳的重复频率(重频)进行倍增,使光频梳重频从最初的250 MHz 提升至10 GHz,对应的脉冲间距从1200 mm 缩减至30 mm,极大地降低了脉冲互相关测距方法对参考臂扫描范围的需求.建立了重频倍增光频梳的时域互相关干涉信号数学模型,通过数值模拟分析了光源参数(重频、起始偏移频率)和法布里-珀罗标准具参数(色散、腔长、中心波长)对滤出光谱形状以及互相关信号的影响.在实验中,使用重频倍增后的光频梳进行脉冲互相关干涉绝对测距,与参考干涉仪对比,在210 mm 范围内获得优于4 µm 的测距精度.

1 引言

高精度绝对距离测量是科学研究和工业应用的先导和基础,直接影响着先进装备制造(大飞机、大型科学装置)和空天探测技术(卫星编队、引力波探测)等前沿科学领域的发展[1−3].非相干脉冲飞行时间测距法,受光电探测器响应速度限制,测距分辨力仅在毫米量级,无法满足上述测量场景下的高精度测量需求[4].基于单波长激光的零差、外差干涉方法,可在几十米范围内达到微米级测距精度,但受限于干涉条纹 2π 相位模糊问题,在大空间测量场景下存在测量光束连续引导、断光续接等应用问题,因此测量效率较低、灵活性较差[5].传统多波长干涉[6]、合成波长干涉[7]和光频扫描干涉[8]可通过构建合成波长的方式扩大非模糊范围,实现高精度绝对测距.然而,多波长方法需搭建多路稳频激光,使得光路系统整体结构复杂,测量实时性差且实施难度较大.光频扫描干涉测距可实现无非模糊范围高精度绝对测距,但该方法存在光频扫描非线性、振动误差放大效应,以及干涉信号数据处理量大等问题.

近二十年来,飞秒光频梳的出现为长度计量带来了全新的测量原理和解决思路[9].光频梳本质上可看作上百万个具有固定频率间隔单波长激光的相干叠加,在频域上表现为在宽光谱范围内,众多离散的、均匀分布的光学纵模,时域上对应为脉冲间距极其稳定的超短脉冲序列[10].光频梳本身作为可溯源频率源,当为其他测距光源提供高稳定度频率参考时,其测距结果可直接溯源至时间频率基准[11,12];光频梳本身亦可作为多波长光源,通过此单光源就可直接构建多级微波合成波长进行相位干涉测距[13,14];光频梳梳齿状离散宽光谱特性革新了传统白光光谱干涉测距方法,将原本白光光谱干涉的最大测距范围拓展定义为非模糊距离,在光频梳的相干长度范围内,实现大非模糊范围下的快速高精度绝对测距[15,16].此外,基于其脉冲相干特性,创新性地将脉冲飞行时间测距方法与干涉相干测距结合,以脉冲间距为刻尺(几十毫米至几米),依靠光学采样来降低对光电探测器带宽的需求,提高脉冲飞行时间的分辨力,理论上可达到优于半波长的测距精度.为实现参考脉冲扫描范围覆盖脉冲间距,依据实施采样方式的不同,衍生出机械扫描参考臂的脉冲互相关法[17]、腔长调节光学采样法[18,19],双光梳异步光学采样法[20,21].腔长调节光学采样法:在参考臂中引入长光纤作为脉冲存储器,在调节光频梳重复频率(重频)时,利用倍数放大效应实现测量脉冲对参考脉冲的大范围扫描采样,参考长光纤作为敏感器件,极易受到温度、压力、振动的影响而产生光程漂移,长光纤锁定技术虽能在一定程度上解决此问题,但锁定系统较为复杂(需额外的外差干涉仪、锁相环、可变光纤拉伸器等),锁定时间较短等问题使得该测量方法应用受限.双光梳异步光学采样法:由两台具有微小重频差的光频梳作为测距光源,系统无需任何移动部件便可自动完成极快的光学采样,实现快速、高精度任意测距,但双光梳作为多参数系统,重频、重频差、中心波长、光谱带宽的选择与优化,双光频梳的同步锁定与噪声抑制,使得构建稳定的双光梳系统难度较大且成本高昂.

直接使用高重频光频梳进行脉冲互相关测距,是最简单、直接且有效的时域互相关测量方法,电光调制频率梳、微腔光频梳的不断发展,使得应用高重频光频梳测距成为可能[22,23].但如今较为成熟的商用化钛宝石光频梳和光纤锁模光频梳,受激光器腔长的物理长度限制,重频一般较低,通常在几十MHz 至1 GHz 之间,对应脉冲间距为百毫米至几米之间.以德国公司Menlosystems 的超低相位噪声光频梳FC1500-250-ULN 为例,重频仅为250 MHz,对于脉冲互相关测距方法而言,这就要求机械扫描导轨范围长达600 mm,才能确保参考脉冲与测量脉冲在任意位置重合,实现无测量盲区测距.长距离的位移扫描将直接影响测量效率,扫描位移台的性能也将影响干涉条纹质量.针对此问题,科研人员提出借助外部设备来对低重频光源的重频进行扩增,将原本几百MHz 的重频提升至几十GHz,从而压缩脉冲间距,降低对参考臂扫描范围的需求.现已证实的重频倍增方法包括:法布里-珀罗(Fabry-Perot,FP)标准具法[24]、Mach-Zehnder 干涉仪法[25]、光纤谐振腔法[26].光纤谐振腔和Mach-Zehnder 干涉仪均需依靠多级串联的方式实现高倍数的重频倍增,这使得系统整体结构复杂,稳定性较差.FP 标准具作为一种光学滤波器,基于多光束干涉原理,可对光频梳的精密光谱进行选频,单腔即可完成高倍数重频倍增.但在使用FP标准具滤除光频梳梳齿过程中,光频梳光源与FP标准具的匹配程度,从根本上影响着滤除光谱的形状和强度,就脉冲互相关测距方法而言,光频梳光谱分布将直接决定互相关条纹的质量[27],并影响最终的测距精度.

本文分析了FP 标准具对光频梳进行重频倍增的原理,通过数值模拟详细讨论了光源参数(重频、起始偏移频率(偏频))和FP 标准具参数(腔长、中心波长、色散)对滤除光谱(形状、强度)的影响,并建立了较为全面的重频倍增光频梳时域干涉数学模型,分析了在不同透射光谱情况下的互相关干涉条纹.实验中,利用FP 标准具将250 MHz重频的光频梳倍增至10 GHz,将脉冲间距由原先的1200 mm 压缩至30 mm,并以10 GHz 的重频倍增光源进行脉冲互相关绝对测距实验.在210 mm范围内,与条纹计数干涉仪进行比对,最终测距结果与参考真值的一致性优于4 µm.

2 测量原理

2.1 脉冲互相关测距基本原理

如图1 所示,由光频梳发出的超短脉冲序列(红色脉冲)射入非等臂Michelson 干涉光路时,参考脉冲序列(绿色脉冲)Eref=∑ANexp[i(2πfNt)]射向扫描参考镜,测量脉冲序列(蓝色脉冲)Emea=∑ANexp[i(2πfN(t−τ))]射向目标镜,两束光在分束镜位置汇合并发生干涉,由光电探测器采集、示波器进行显示.为了简化分析,在此认为两臂测量脉冲光强相等,AN为第N根梳齿的光强幅值,fN为第N根梳齿的光频,可表示为fN=N·fr+f0.这里,N为正整数,对应光源光谱范围,通常取值区间为几万至几十万;fr为重频,表示相邻纵模的频率间隔;f0为偏频;τ=2nL/c为测量脉冲相对于参考脉冲的延迟时间,其中n为空气的折射率,L为测量脉冲相对参考脉冲的几何路径差,c为光在真空中的传播速度.在光电探测器的响应周期T内,由其测得的干涉信号可表示为

图1 脉冲互相关测距原理示意图(OFC,光频梳;BS,分束镜;Mr,参考镜;Mt,目标镜;M0,零点参考镜;PD,光电探测器;Scope,示波器)Fig.1.Schematic of the pulse cross-correlation ranging principle.OFC,optical frequency comb;BS,beam splitter;Mr,reference mirror;M0,zero position mirror;Mt,target mirror;PD,photodetector;Scope,oscilloscope.

由(1)可知,互相关干涉信号由直流分量IDC和交流分量IAC=两部分组成,其中N∈[Ni,Nf],Ni和Nf分别表示光谱范围内第一根梳齿和最后一根梳齿所对应的整数N.为了便于后续分析,将IAC表示为功率谱密度与余弦函数的形式

由(2)式可知,交流信号的强度主要由光源的功率谱密度S(fN)和两臂光程延迟τ所决定.根据Wiener-Khinchine 定理,自相关函数Γ(τ)和S(fN) 互为傅里叶变换关系,因此IAC可进一步表示为

其中,fc为光频梳的中心载波光频;k为测量脉冲相对参考脉冲的序数差;Δφceo为载波包络偏移相位,Δφceo=2πf0/fr;⊗和δ分别表示卷积函数和冲击函数;T为互相关干涉信号的重复周期,T=1/fr.对于任意的被测距离L可表示为

Lpp为脉冲间距,决定互相关测距的非模糊范围,Lpp=c/(fr·ng),其中ng为空气的群折射率;l表示测量脉冲相对参考脉冲的非整数脉冲间距部分,l=2L−kLpp.(3)式中Γ(τ)的定义域为 [−τW/2,+τW/2],其中τW表示飞秒脉冲的宽度.所以在目标移动过程中,只有测量脉冲与参考脉冲在空间重合位置的几十微米范围内才会出现高对比度的干涉条纹,其他位置几乎均为直流分量.为实现测量脉冲与参考脉冲在任意位置重合,这就要求参考位移台的扫描范围覆盖Lpp/2,提高光频梳的重频可间接缩短脉冲间距,降低对位移台扫描范围的需求.

2.2 FP 标准具滤梳齿仿真分析

FP 标准具可看作一种基于多光束干涉原理的光学滤波器,可对光频梳的光谱进行操控和选频.如图2 所示,FP 标准具通常由两片高反射率的镜片组成,一般情况下,FP 标准具的透过函数可表示为

图2 FP 标准具滤梳齿示意图Fig.2.Schematic of filtered frequency comb using a FP etalon.

式中R是反射镜的镜面反射率,f为输入光频的大小,Lc为FP 标准具的腔长,nc为FP 标准具的折射率.FP 标准具的自由光谱范围fFSR=c/(2ncLc),决定了透射光谱峰之间的频率间隔;精细度Fa==fFSR/fFWHM,决定了透射光谱峰的宽度,Fa主要由镜面反射率R所决定,fFWHM表示透射光谱峰值信号的半高全宽.当FP 标准具的精细度较高时,上述FP 标准具的透过函数T(f) 可看作是洛伦兹函数p(f) 与梳齿函数的卷积:

当光频梳入射至FP 标准具时,透射光谱可表示为

透射光谱EFP(f)表现为光谱形状为E(f)、受T(f)调制的梳状谱.接下来仿真分析光源参数(重频fr、偏频f0)、FP 标准具参数(腔长、色散、中心波长)对透射光谱的影响.仿真中,光频梳参数以Menlosystem 的超低相位噪声光频梳FC1500-250-ULN 为原型,光频梳的中心波长为1560 nm,光谱全宽约为100 nm,重频fr=250 MHz,为简化分析,偏频f0=0 MHz,光频梳的光谱形状设置为理想高斯型光谱,时域脉冲半高全宽为120 fs.FP 标准具的参数为:腔镜的反射率R=99%,对应精细度为312 左右,腔长Lc=15 mm,对应fFSR=10 GHz.

首先分析光源参数对透过光谱的影响.在fr=250 MHz,f0=0 MHz,fFSR=10 GHz 情况下,光频梳的梳齿光频与FP 标准具透过函数满足如下关 系:fN=N·fr+f0=q·fFSR,40fr=fFSR,即光频梳第一根梳齿对应的光频等于整数倍的自由光谱范围,同时自由光谱范围等于40 倍的重频.如图3(a)和3(b)所示,在初始光谱的基础上每隔40 根梳齿滤出一根梳齿,滤出后单根梳齿的强度并未发生明显变化,光谱形状仍为严格的高斯型.图3(c)和图3(d)则表示fr=250.001 MHz 情况下的透射光谱,标准具的fFSR与fr的基础频率匹配偏差为 δf=(40fr−fFSR)=0.4 MHz .可以看出,随着光频逐渐的增大,基于游标倍增效应,频率匹配偏差q×δf变大,但仍处在 [−Δf/2,Δf/2] 之内,始终被涵盖在标准具的单透射峰内(Δf为FP 标准具透射峰的宽度),透射光谱整体仍可保持高斯型,但透射梳齿的强度整体变低.图3(e)和图3(f)表示fr=250.1 MHz 情况下的透射光谱.可以看出,在光频较小的情况下,频率匹配偏差q×δf∈[−Δf/2,Δf/2],但随着频率的增加,透射梳齿幅值越来越低,当,梳齿光强全部被抑制掉.随着光频进一步的增大,频率匹配偏差将落在下一个FP 标准具透射峰内,即q×δf∈[fFSR−Δf/2,fFSR+Δf/2],部分梳齿信号重新落在FP 标准具的透射峰内.因此当fr与fFSR的匹配偏差较大时,透射光谱将出现周期性、间断性的局部模式匹配,滤出后的光谱遭受严重调制.

图3 光频梳重频对透射光谱的影响(a)理想透射谱;(c) fr=250.001 MHz 时的透射谱;(e) fr=250.1 MHz 时的透射谱;(b),(d),(f)分别为(a),(c),(e)的局部放大图Fig.3.Influence of frequency comb repetition rate on transmission spectrum:(a) Ideal transmission spectrum;(c) the transmission spectrum at fr=250.001 MHz;(e) the transmission spectrum at fr=250.1 MHz;(b),(d),(f) partial enlarged views of (a),(c),(e).

接下来分析偏频变化对透射光谱的影响,图4(a)分别表示f0=0,5,15,25,35 MHz 情况下透射光谱整体的变化,光谱形状并未遭受较大的调制,仍保持高斯型,但偏频的变化使得整体梳齿信号相对FP 标准具透过函数进行平移,整体强度变低.图4(b)表示不同偏频下,局部透射光谱单梳齿强度的变化,随着偏频不断增大,透射梳齿强度越来越低.但总的来说,光频梳偏频变化对透射光谱的影响较小.

图4 光频梳偏频对透射光谱的影响,其中(b)为(a)的局部放大图Fig.4.Influence of frequency comb offset frequency on transmission spectrum.(b) is partial enlarged view of (a).

紧接着分析FP 标准具参数对滤出光谱的影响.由公式 δf=(40fr−fFSR) 可知,腔长变化引起的频率匹配偏差与重频变化引入的匹配偏差变化规律近似.图5(a)和图5(b)表示FP 标准具存在不同腔长误差(0—5 nm,以1 nm 为增量)下的透射光谱,透射光谱对腔长变化十分敏感,5 nm 腔长偏差下透射光谱强度已经降低很多.实际上,因光频梳包含众多频率成分,FP 标准具器件本身的色散将导致自由光谱范围fFSR=c/(2ncLc)不再是一固定的常数,其将会随着光频的变化而变化,这在一定程度上限制透射光谱的带宽.为方便理解,将(5)式中的干涉相位重新表示为

图5 FP 标准具参数对透射光谱的影响 (a)不同腔长偏差;(c)不同群延迟色散,(e)不同中心波长;(b),(d),(f)分别为(a),(c),(e)的局部放大图Fig.5.Influence of FP etalon parameter on the transmission spectrum:(a) Different cavity length deviations;(c) different group delay dispersions;(e) different center wavelengths;(b),(d),(f) partial enlarged views of (a),(c),(e).

公式(8)中的第一项,干涉相位与光频保持线性关系,第二项中的φ(w) 为因色散引入的非线性相位,将其在FP 标准具的中心角频率wd处以泰勒级数展开为

其中,常数相位项φ0所带来的影响相当于对整体标准具透射函数进行平移,可通过适当调节光频梳光源的偏频f0来修正;线性相位项φ1的影响相当于均匀地改变了标准具的fFSR,可通过适当调整光频梳的重频进行匹配;但后两项由二阶项群延迟色散(group delay dispersion,GDD)和高阶色散引入的额外相位则无法进行补偿,这将导致随着梳齿频率的增加,光频梳梳齿与FP 标准具fFSR的匹配程度越来越差,限制透射光谱的带宽.图5(c)和图5(d)表示FP 标准具中心波长为1560 nm,干涉相位中添加不同程度GDD 时的透射光谱.可以看出随着GDD 的不断增加,在远离腔镜中心波长的位置,频率匹配偏差越来越明显,透射梳齿信号的幅值不断变低,导致透射光谱宽度变窄.图5(e)和图5(f)表示GDD 为 60 fs2,中心波长分别在1550,1560,1570 nm 情况下的透射光谱.因FP 标准具中心波长与光频梳中心波长不同,透射光谱峰值点位置发生偏移,频率匹配偏差将以FP 标准具中心波长位置向两边递增,导致透射光谱的形状扭曲,呈现出非对称的形状.

图6 模拟了光源参数(重频、频偏)和FP 标准具参数(腔长、群延迟色散)对透射光谱强度的影响.图6(a)表示fr在250 MHz±30 kHz 范围内,以10 Hz 为步长递增时,透射光谱整体强度的变化情况.可以看出,重频在大范围变化过程中,光频梳重频与透射函数的匹配变化剧烈,透射光谱强度展现出振荡式的变化,当光频梳梳齿与fFSR匹配程度较差时,透射光强仅为入射光强的0.2%.因此,若要实现光频梳重频与FP 标准具自由光谱范围的严格匹配,应确保腔长匹配精度达到纳米级的水平.图6(b)表示f0在±35 MHz 范围内,以10 kHz为步长变化时透过光谱强度的变化情况.偏频变化相当于光频梳相对FP 标准具透射函数进行平移,强度变化趋势比较平缓和规律,当f0=35 MHz时,透射光强为入射光强的0.8%.图6(c)表示FP 标准具腔长Lc在15 mm±0.5 µm 范围内,以1 nm 为步长变化过程中,透射光谱强度的变化情况.由公式 δf=[40fr−c/(2Lc)]可知,频率匹配偏差由fr和Lc共同决定,因此腔长变化过程中透过光谱强度的变化规律与重频变化规律相同.图6(d)表示FP 标准具GDD 在±100 fs2范围内,以1 fs2为增量时的透射强度变化规律.随着GDD 的不断增加,透射光谱不断变窄,透射强度也随之降低.

图6 光源参数对透射光谱强度的影响 (a) 重频;(b) 偏频;(c) 腔长;(d) 群延迟色散Fig.6.Influence of light source parameter on the transmission intensity:(a) Repetition rate;(b) offset frequency;(c) cavity length;(d) group delay dispersion.

2.3 重频倍增光频梳时域互相关干涉

在确定光源参数、FP 标准具参数对透射光谱的影响后,仿真分析对应透射光谱下的时域互相关干涉信号.经过FP 标准具调制后的光频梳,时域互相关干涉信号交流部分可表示为

由(10)式可知,经FP 标准具后的互相关信号周期由光频梳的重频fr以及自由光谱范围fFSR共同决定.当fFSR=m·fr,且fN=N·fr+f0=q·fFSR时,互相关信号周期可缩减m倍.接下来,在绝对距离为12 m 位置附近分析光频梳重频、偏频、FP标准具色散、精细度,对脉冲互相关干涉信号的影响.仿真中在12 m±50 µm 范围内,以30 nm 为增量扫描参考臂,图7(a)和图7(b)为fr=250 MHz+20 Hz,250 MHz+60 Hz,以及250 MHz+100 Hz情况下的互相关条纹,重频微小的变化便使得互相关条纹幅值迅速降低,条纹对比度变差;图7(c)和图7(d)为f0=5,20,35 MHz 情况下的互相关条纹,因偏频对透射光谱形状和强度的影响较小,较大偏频的变化对互相关干涉条纹的影响也较弱.

图7 光源参数对互相关信号的影响 (a) 不同重频偏差下的干涉条纹;(c) 不同偏频偏差下的干涉条纹;(b),(d) 分别为(a),(c)的局部放大图Fig.7.Influence of light source parameter on cross-correlation signal:(a) Interference fringe under different repetition rate;(c) interference fringe under different offset frequency;(b),(d) partial enlarged views of (a),(c).

由图6(a)和图6(c)的仿真可知,标准具腔长变化与重频变化对透射光谱影响近似,在此省略了标准具腔长变化对互相关信号的影响.图8(a)和图8(b)表示标准具GDD 分别为 20,60,100 fs2情况下的互相关条纹,结合图5(c)中的仿真结果容易理解,色散导致滤出光谱宽度变窄,基于光谱宽度与脉冲宽度的反比关系,对应时域信号的脉冲宽度变宽,互相关条纹也随之展宽;图8(c)和图8(d)表示标准具精细度Fa=210,520,1050 情况下的互相关条纹,条纹质量和对比度均未发生明显变化.基于上述四种因素的分析,单纯重频变化对互相关信号的影响最为明显.然而在实际实验中,基于FP 标准具的重频倍增脉冲互相关系统,上述多种因素将同时作用于脉冲互相关信号,实际情况要比上述分析复杂得多.

图8 FP 标准具参数对互相关信号的影响 (a) 不同群延迟色散下的干涉条纹;(c) 不同精细度下的干涉条纹;(b),(d) 分别为(a),(c)的局部放大图Fig.8.Influence of FP etalon parameter on cross-correlation signal:(a) Interference fringe under different group delay dispersion;(c) interference fringe under different finesse;(b),(d) partial enlarged views of (a),(c).

3 实验及分析

测距系统如图9 所示,共由两部分组成,分别是重频倍增光频梳光源和互相关绝对测距光路.重频倍增光源部分:由德国MenloSystem 公司的超低相位噪声光频梳FC1500-250-ULN 作为种子光源,fr=250 MHz (可调范围±2.5 MHz),f0=35 MHz(可调范围±5 MHz),两参数均被锁定至原子钟频率基准(Symmetricom,8040 C),由其发出的脉冲序列经准直器CL1 (Thorlabs,RC02APC-P01)准直后射入重频倍增透镜组,经过模式匹配透镜1(f=100 mm)、FP 标准具(Thorlabs,SA210-B)、模式匹配透镜2 (f=100 mm)扩展为fr=10 GHz的光频梳光源.在FP 标准具的滤梳齿作用下,透射后的光强遭到极大地缩减,在此将透射光经铒镱共掺光纤放大器(TOP Photonics,EYDFA)放大至16 dBm.绝对测距部分共由三组Michelson 光路组成,首先经倍增后的光频梳光源射入到第一组非等臂Michelson 干涉光路中,测量光路中的半反半透薄膜SRM (Thorlabs,BP145B3)作为测距的绝对零位参考,将扫描参考镜CR1 放置到高精度扫描位移台(PI,M-521.DD.1)上,由目标镜CR2、绝对零位参考镜反射的光与扫描参考镜反射的光在BS1 位置合束并发生干涉,由光电探测器(Thorlabs,PDB470C)接收.搭建的第二路干涉仪主要用来确定第一路干涉仪中两干涉峰之间的距离,所用光源为RIO ORION Laser Module,波长为1550.12 nm、输出功率约为20 mW,所使用的光电探测器为Thorlabs 的PDA10CS-EC,为尽可能地减小阿贝误差,目标镜CR4和扫描参考镜CR1 背靠背放置在扫描位移台上.第三路干涉仪为雷尼绍干涉仪XL-80,其目标镜CR6与互相关干涉光路中的目标镜CR1 背靠背摆放,XL-80 的输出值作为测量真值来验证系统的测距精度.雷尼绍的环境参数传感器XC-80 来记录实验过程中的温度、湿度、气压,用以计算空气折射率并修正空气中的波长[28].最终两光电探测器采集到的干涉信号由示波器(Tektronix,MDO4050C)进行显示.

图9 重频倍增光频梳时域互相关绝对测距原理示意图(Frequency comb,光频梳;CL1—3,准直器;Lens1 和Lens2,模式匹配透镜;Cavity mirror,FP 标准具的腔镜;EYDFA,铒镱共掺光纤放大器;M1,平面反射镜;BS1—3,光学分束器;CR1—6,角锥棱镜;PD1 和PD2,光电探测器;SRM,半反半透薄膜;RIO,单频激光器;XL-80,雷尼绍干涉仪;Scope,示波器;黄色线,单模光纤;深蓝色线,电学线缆;红色线,光频梳出射的激光光束;浅蓝色线,RIO 单频激光器出射的激光光束;绿色线,雷尼绍干涉仪出射的激光光束)Fig.9.Schematic of cross-correlation ranging based on repetition rate multiplying optical frequency comb.OFC,optical frequency comb;CL1—3,collimator;Lens1 and Lens2,mode matching lens;EYDFA,erbium ytterbium doped fiber amplifier;M1,plane mirror;BS1—3,beam splitter;CR1—6,retroreflector;PD1 and PD2,photodetector;SRM,semi-reflective film;RIO,RIO single wavelength laser;XL-80,Renishaw commercial interferometer;Scope,oscilloscope;yellow line,single-mode fiber;dark bule line,the electrical cable;red line,the laser beam emitted by the optical frequency comb;light blue line,the laser beam emitted by the RIO single wavelength laser;green line,the laser beam emitted by the Renishaw interferometer.

图10(a)和图10(b)为光频梳经FP 标准具前后,由光谱仪(Yokogawa,AQ6370D)采集的光谱.起始状态下,光频梳光谱全带宽约为100 nm,经FP 标准具滤除后的光谱明显变窄,主要原因在于光频梳经FP 标准具后光强遭到了极大的缩减.在使用EYDFA 进行功率放大时,实验中所使用的EYDFA 的放大波段为1530—1570 nm,因此透射光谱变窄许多.图10(c)和图10(d)为光频梳经FP 标准具前后,由30 GHz 带宽高速光电探测器(Alphalas,UPD-15-IR2-FC)搭配26.5 GHz 带宽信号分析仪(Agilent,N9020A)获取的射频谱.可以看出,滤出光频梳射频谱的边模抑制比优于20 dB.

图10 光频梳经过FP 标准具前后的光谱和射频谱 (a) 未经FP 标准具的光谱;(b) 经过FP 标准具后的光谱;(c) 未经FP 标准具的射频谱;(d) 经过FP 标准具后的射频谱.Fig.10.Optical spectrum and the radio frequency spectrum of the optical frequency comb before and after passing through the FP cavity:(a) Optical spectrum before FP cavity;(b) optical spectrum after FP cavity;(c) radio frequency spectrum before FP cavity;(d) radio frequency spectrum after FP cavity.

图11(a)为经FP 标准具后的透射光谱结合(10)式仿真获得的脉冲互相关干涉信号.仿真中,位移步长设置为50 nm,扫描范围为±100 µm.图11(b)为实验中在等臂位置附近,扫描位移台获取的脉冲互相关干涉信号.从结果易知,实验获得干涉图样与仿真结果具有较好的一致性.图11(c)则为在等臂位置附近60 mm 范围内连续扫描参考臂获得的四组脉冲互相关干涉条纹.

图12(a)为测距实验获得的互相关干涉图样,其中峰值较低的互相关信号由零位镜产生,而峰值较高的互相关信号由目标镜所产生,紫色干涉条纹由单波长激光所产生.图12(b)为目标镜互相关信号峰值位置附近的放大图样,因经过FP 标准具,光源光谱变窄、对应时域信号脉冲宽度变宽,峰值位置附近的多个载波信号幅值接近.如图12(c)先获取载波信号多个峰值点位置,再通过高斯曲线拟合的方式来确定互相关信号包络峰值位置.在测量臂与参考臂等臂位置附近的210 mm 范围内,以30 mm 为步长,连续对8 个位置进行测量,单点快速采集5 组数据.图13 为长度比对实验结果,与XL-80 的 最大残差优于4 µm,单点重复性优于1.5 µm.

图12 重频倍增光频梳脉冲互相关测距图样 (a)互相关干涉信号;(b) 互相关干涉信号峰值附近位置放大图样;(c) 包络峰值信号提取Fig.12.Cross correlation pattern based on repetition rate multiplying frequency comb:(a) Cross correlation signal;(b) magnified view of the horizontal axis near the envelope peak for (a);(c) extracted envelope and the Gaussian fitting results.

图13 长度比对实验结果(红色圆圈表示五次测量结果的平均值;误差棒表示测量结果的标准差)Fig.13.Length comparison experiment result. The red circle represents the average result of five repeated measurements.Error bar indicates the standard deviation of the measurement result.

对影响测距精度的因素进行分析,误差主要来源于以下四个方面:首先,因FP 标准具腔长并未锁定至光频梳的重频上,FP 标准具的自由光谱范围与光频梳匹配较差且不够稳定,温度变化、微小机械振动、光机件的蠕变,都将影响互相关干涉条纹的质量,使得在利用高斯拟合算法提取干涉包络峰值位置时,存在一定的峰值拟合误差,影响脉冲重合位置的判读精度;其次,在FP 标准具和EYDFA 的共同作用下,滤出光谱宽度变窄、光谱形状不规则,对应时域脉冲宽度展宽,在进行峰值提取时,多个载波条纹峰值强度近似,峰值提取精度很难达到半波长;此外,由辅助干涉仪确定零位镜和目标镜所产生两干涉峰之间的距离时,辅助干涉仪条纹计数也会存在半波长的截取误差;最后,因放置目标镜的滑轨存在一定的漂移和蠕动,在单点数据采集过程中,会因与参考干涉仪数据采集不同步而引入一定的误差.针对上述提及的问题,后续将对光路进行优化,将FP 标准具的腔长严格锁定至光频梳重频上,以尽可能地保证系统的稳定性和干涉条纹质量,将测距精度提升至半波长以内,之后也将在长距离下来验证该方法的测距精度.

4 结论

就飞秒脉冲互相关绝对测距方法而言,光频梳的重频直接决定非模糊距离的大小,针对现有商用光纤光频梳重频较小、脉冲间距过大的问题.本文应用FP 标准具对光梳重频实施倍增,对标准具滤梳齿的原理进行了详细分析,建立了较为全面的重频倍增光频梳脉冲互相关干涉模型,分析了光源参数(重频、偏频)、标准具参数(腔长、色散、中心波长)对滤出光谱(形状、强度)和脉冲互相关信号的影响.在实验中,使用FP 标准具将光纤光频梳的重频从250 MHz 倍增至10 GHz,对应脉冲间距从1200 mm 压缩至30 mm,并利用倍增后的光源进行脉冲互相关绝对测距实验,在210 mm 测量范围内与条纹计数干涉仪进行精度比对,测量结果与参考真值的一致性优于4 µm.