光学学报, 2020, 40 (20): 2030001, 网络出版: 2020-10-16   

布里渊光谱仪的高精度波长标定方法研究 下载: 957次

Wavelength Calibration of Ultra-High Resolution Brillouin Spectrometer
作者单位
1 中电科仪器仪表有限公司, 山东 青岛 266555
2 中国电子科技集团有限公司第四十一研究所电子测试技术重点实验室, 山东 青岛 266555
摘要
针对最高分辨率可达0.1 pm的布里渊近红外光谱仪波长标定的难题,基于多次标定逐步降低波长误差的思路,提出了一种协同融合标定原理的标定流程:首先,基于布里渊光谱分析理论模型,结合泵浦信号波长、受激布里渊频移、增益谱线型函数等推导出超高分辨率光谱的理论波长;然后,采用FP标准具进行相对波长的标定;最后,基于气室完成绝对波长的标定。气室中充入HCN+ 12CO+ 13CO混合气体,以实现C+L波段的覆盖。针对可能存在的光谱鬼线,提出了一种基于受激布里渊增益谱与衰减谱共生特性的鬼线识别方法。实验结果表明:布里渊光谱仪的理论波长偏差超过了10 pm,相对标定可将波长不确定度降低到±2.5 pm,绝对标定可将波长不确定度进一步降低到±0.035 pm。
Abstract
To address the problem of wavelength calibration of a Brillouin NIR spectrometer with a maximum resolution of 0.1 pm, a calibration process using the principle of collaborative fusion calibration is proposed. This process is based on the idea of gradually reducing the wavelength error of multiple calibrations. The calibration process is as follows: first, the theoretical wavelength was calculated based on the Brillouin spectral analysis modal and relative parameters, such as the pump signal wavelength, Brillouin frequency shift, and gain in spectrum lineshape; then, the relative wavelength was calibrated using the Fabry-Perot etalon; finally, the absolute wavelength calibration based on the gas chamber was completed. The gas chamber was filled with a mixed gas of HCN, 12CO, and 13CO to achieve absolute wavelength calibration in the C and L bands. A ghost identification and correction algorithm was also presented based on the commensalism of gain and loss spectrum of stimulated Brillouin scattering. Experimental results show that the theoretical wavelength deviation is more than 10 pm, and wavelength uncertainty of ±2.5 pm can be achieved by relative wavelength calibration. Moreover, the wavelength uncertainty of the absolute wavelength calibration is ±0.035 pm.

1 引言

光谱分析是通信、传感、分子光谱学、微波生成等光学应用中的关键诊断步骤。随着光通信和传输网不断朝着超大容量、超高速率、超长传输距离方向发展,在波分复用等技术的基础上引入各种先进调制格式提升光网络的性能成为必然选择,而这会导致通道间隔越来越小。例如,光正交频分复用信号的通道间隔通常为20~50 MHz,这就要求光谱分析仪的光谱分辨率要达到pm量级甚至更高,以满足先进光传输系统的测试和监测,以及光放大器、光纤滤波器等有源和无源光器件的量化分析[1-2]

布里渊光谱仪基于光纤受激布里渊效应的选频进行放大,通过可调谐激光器的推扫实现感兴趣光谱区域的测量分析。得益于非线性光纤受激布里渊效应的独特优势,测量光谱数据可以实现高达40 dB的近峰动态范围,这就允许研究人员可以对靠近光学载波的信号成分进行精确测量。因此,布里渊光谱分析技术优于色散光栅、法布里-珀罗(FP)滤波器或外差相干等光谱分析技术,是一类非常有前景的光谱分析路线[3]

布里渊光谱仪具有超高的光谱分辨率,其测量波段范围可以覆盖C、L等光通信波段。因此,如何实现布里渊光谱仪的高精度波长标定是一个极具挑战性的难题。现有的波长标定方法,如激光光源法、卤素灯法、积分球光源法、FP标准具法、气体吸收法等[4-8]都存在一定的局限性,难以满足新一代光通信用超高分辨率光谱仪的波长标定需求。本文提出了一种协同融合标定原理的pm量级分辨率的波长标定流程:首先根据布里渊光谱分析模型得到理论波长,然后基于FP标准具进行相对波长的标定,最后基于气体吸收法完成绝对波长的标定。为提高标定精度,本文将宽工作波段分为若干个子区间,分别计算标定系数。实验结果表明,标定波长的不确定度优于±0.035 pm。此外,本文还识别校正了光纤受激布里渊效应引入的光谱鬼线。

2 布里渊光谱分析仪的原理

布里渊光谱分析仪的原理模型可用数学公式描述为[9-10]

dIp(z,λ)dz=-αIp(z,λ)-gBAeffIs(z,λ)Ip(z,λ)-Iesp(z,λ),(1)dIs(z,λ)dz=αIs(z,λ)-gBAeffIs(z,λ)Ip(z,λ)-Iesp(z,λ),(2)dIB(z,λ)dz=gBλ-2πΩBAeffIp(z,λ)Isz,2πΩB,(3)

式中:z为介质光纤上的位置;λ为波长;IpIsIB分别为泵浦信号、待测信号和斯托克斯后向散射信号;Aeffα分别为介质光纤的有效面积和衰减系数;Iesp为由热波动引起的自发布里渊散射信号;ΩB为受激布里渊频移,ΩB为正时产生衰减谱,ΩB为负时产生增益谱;gB为布里渊增益系数。

待测信号Isz=L处注入待测光纤,泵浦信号由z=0处注入待测光纤,由(1)~(3)式所述受激布里渊过程可知,待测信号经受激布里渊放大后,在z=0处的输出信号为

Is(0,λ)=Is(0,λB)+Is(0,λλB)=Is(L,λB)exp[gBIp(0,λp)Leff]+Is(L,λλB)(4)

式中:λB为布里渊峰值波长;λp为泵浦光波长;Is(0,λB)表示待测信号中被受激布里渊放大的光谱组分;Is(0,λλB)表示待测信号中未被放大的带外光谱组分;Is(L,λB)表示z=L处经布里渊放大的光谱成分;Is(L,λλB)表示z=L处未被放大的带外光谱组分;Ip(0,λp)表示由z=0处注入光纤的泵浦信号功率;λp为泵浦波长,λp=λB+2πB;L表示介质光纤的长度;Leff表示介质光纤的有效长度。

当没有泵浦光注入时,在z=0处测量得到的输出信号为

Is~(0,λ)=Is(L,λB)+Is(L,λλB)(5)

布里渊光谱分析仪测量得到的分光光谱信号Sp(λ),可以通过移除待测信号中的未放大带外组分得到,即

Sp(λ)=Is(0,λ)-Is~(0,λ)=Is(L,λB){exp[gBIp(0,λp)L]-1}(6)

布里渊光谱分析仪的工作原理如图1所示。基于前面的数学描述,光谱测量流程如图2所示。利用光纤受激布里渊效应,通过选频放大得到待测信号在某一波长处的光谱成分所对应的分光光谱信号,并通过可调谐激光器连续推扫,获取待测信号在不同波长处的光谱成分所对应的全波段分光光谱,然后通过光谱重构算法得到待测信号的测量光谱数据Sm(λ),就可以实现感兴趣波长区域的光谱测量与分析。

图 1. 布里渊光谱仪的工作原理

Fig. 1. Principle of Brillouin spectrometer

下载图片 查看所有图片

图 2. 布里渊光谱仪的光谱测量流程。(a)待测信号;(b)泵浦信号;(c)受激布里渊效应放大后的输出信号;(d)基底移除后的分光光谱信号;(e)采用可调谐激光推扫得到的全波段分光光谱信号;(f)利用光谱重构算法得到的测量光谱

Fig. 2. Spectra measurement procedure of Brillouin spectrometer. (a) Signal under test (SUT); (b) pump signal; (c) output signal amplified by stimulated Brillouin scattering; (d) split spectra after baseline removing; (e) full band split spectra by tunable laser; (f) measured spectra obtained by spectrum reconstruction algorithm

下载图片 查看所有图片

3 波长标定流程

3.1 基于布里渊光谱分析模型计算理论波长

由布里渊光谱仪的工作原理可知,增益的谱型近似为高斯线型。测量光谱的理论波长可由泵浦信号波长和布里渊频移给出,即

λB=λp-2π/ΩB=λp-λp2nva,(7)

式中:n为光纤折射率;va为光纤中的声速。

需要指出的是,受泵浦信号波长漂移、仪器线型函数、器件缺陷等的影响,理论波长与真实波长之间存在一定偏差,需要进行校正。

当布里渊频差为负时会引入衰减谱,此时可能会在测量光谱中引入鬼线。因此,在相对和绝对波长标定前,需要根据受激布里渊效应增益谱与衰减谱的共生特性,对此类鬼线进行识别和校正。

3.2 基于FP标准具的相对波长标定

FP标准具的光路图如图3(a)所示,其核心单元是两块固定的、具有很小楔角且涂有高反射涂层的平板玻璃。当波长连续变化的光入射后,光会在两个距离为h的平行平板G1、G2间发生多次反射和折射,出射光强度与波长的关系为

Iout(λ)=4Iin(λ)cos2(δ/2)=4Iin(λ)cos22πλκhcosi',(8)

式中:Iout(λ)为出射光强度;Iin(λ)为入射光强度;κ为平板玻璃G1、G2的折射率;i'为楔角;δ为光程差。对于不同波长的入射光,当满足δ=2π,4π,6π,…时,透射光强达到最大值;当满足δ=π,3π,5π,…时,透射光强为最小值。这样便可以锁定入射光的系列固定波长,如图3(b)所示,用作波长标定时的参考。图3(a)中的符号S、P为FP标准具的光纤耦合接口。

FP标准具在宽波段上拥有众多等间距的固定波长通道,但通道中心波长的准确度大概为数pm量级。因此,可采用FP标准具对理论波长进行相对波长的标定,以降低理论计算波长的误差。布里渊光谱仪的分光光谱信号经过FP标准具后可得到如图3(b)所示的光谱曲线,对曲线进行高斯拟合,由(7)式可得到各个拟合峰值点对应的理论波长λfit。理论波长λfit与FP标准具通道中心波长的标称值λFP之间的关系为

λFP=anλfitn+an-1λfitn-1++a0,(9)

式中:系数anan-1a0等为相对波长的标定系数,可由最小二乘拟合方法计算得到。然后,利用(9)式对理论波长进行标定,得到相对标定波长λrw

图 3. FP标准具的光路图和波长透过特性。(a)光路图;(b)波长透过特性

Fig. 3. Optical path and wavelength transmission property of FP etalon. (a) Optical path; (b) wavelength transmission property

下载图片 查看所有图片

3.3 基于气体吸收法的绝对波长标定

气体吸收法[11]是基于朗伯-比尔定律进行波长标定的。它的基本原理是当入射光经过气室时,光会因气体吸收而产生特征吸收谱线,该特征吸收谱线可用作波长标定的参考基准。经气体吸收后的特征谱线幅值为

I(λ)=I0(λ)exp[-κ(λ)cl],(10)

式中:I0(λ)与I(λ)分别为入射光和出射光的强度;κ(λ)为吸收系数;c为气体的体积分数;l为光与气体相互作用的长度。

由HITRAN光谱数据库[12]可知,HCN在1520~1565 nm之间有54条强吸收谱线,12CO在1560~1600 nm之间有41条强吸收谱线,13CO在1595~1630 nm之间有41条强吸收谱线,如图4所示。HCN+12CO+13CO混合气体可为C+L波段的绝对波长标定提供多个稳定的参考波长。

图 4. HCN+12CO+13CO混合气体的特征谱分布

Fig. 4. Characteristic spectrum of HCN+12CO+13CO mixed gas

下载图片 查看所有图片

进行绝对波长标定时,布里渊光谱仪的分光光谱信号经过气室后可以得到如图4所示的测量谱线;然后采用Voigt线型对测量谱线进行拟合,便可以得到由混合气体的吸收作用而产生的系列特征吸收谱线的波谷点。每个特征谱线波谷点的测量波长λm可采用(7)式给出理论值,并采用(9)式进行波长修正,以降低理论波长的偏差;同时,用作绝对定标的混合气体特征吸收谱线的波长(即混合气体同一特征吸收谱线的测量波长λm与标称波长的偏差,也就是布里渊光谱仪经相对波长标定后的波长偏差)是已知的。采用最小二乘拟合算法对混合气体系列特征吸收谱线的测量波长λm与标称波长进行多项式拟合,即可得到绝对波长标定时的波长标定系数。标称波长与测量波长的关系可描述为

λgas=bnλmn+bn-1λmn-1++b0,(11)

式中:λgas为气体特征吸收谱线的标称波长;系数bnbn-1b0等为绝对波长标定系数。然后,基于(12)式进行绝对标定,得到绝对标定波长λcal,此即为最终的测量光谱数据的标定波长。

λcal=bnλmn+bn-1λmn-1++b0(12)

4 分析与讨论

4.1 鬼线识别与校正

由受激布里渊效应可知,泵浦信号与待测信号间的布里渊频移为正时产生增益谱,为负时产生衰减谱[12],这会导致测量光谱中存在如图5(b)所示的鬼线。鬼线的存在不但会影响后续波长的标定,还会引入错误的虚假光谱信息。

图 5. 布里渊光谱仪产生鬼线的机制。(a)受激布里渊增益特性;(b)包含鬼线的测量数据

Fig. 5. Ghost lines generation mechanism in Brillouin spectrometer. (a) Gain property of stimulated Brillouin scattering; (b) measured data with ghost

下载图片 查看所有图片

鬼线识别算法利用的是受激布里渊效应增益谱与衰减谱共生的特性,而且增益谱与衰减谱之间的频差近似为2ΩB。鬼线识别流程如图6所示:1)通过AMPD算法[13]检出测量光谱中包含的波峰和波谷;2)对成对出现的波峰和波谷进行谱型拟合,如果拟合波峰与波谷间的频差近似为2ΩB,则判定波谷为鬼线;3)用鬼线邻近的基线采样点代替鬼线,完成鬼线的校正。

图 6. 鬼线识别校正流程。(a)波峰拟合;(b)波谷拟合;(c)鬼线校正后的测量光谱

Fig. 6. Ghost lines identification and correction procedure. (a) Peak fitting; (b) valley fitting; (c) spectra after ghost lines correction

下载图片 查看所有图片

4.2 基于FP标准具的相对波长标定实验

进行相对波长标定的实验装置如图7(a)所示,主要包括布里渊光谱仪实验装置和用作标定的FP标准具、偏振控制器、起偏器等。实验所用热稳定FP标准具的波长精度为±1 pm,自由光谱范围为50 GHz。对测量数据进行计算得到理论波长,移除鬼线后再进行插值拟合,以确定各个峰值的理论波长,然后将其与各个峰值的标称值(也就是对应的FP标准具通道的中心波长)代入(9)式计算相对波长标定系数,标定结果如图7(c)所示。完成相对波长标定后,标定波长的残差如图7(d)所示。考虑到FP标准具的波长精度,相对波长标定的不确定度在-2.5~2.5 pm以内。

图 7. 相对波长的标定。(a)实验装置;(b)测量数据;(c)波长标定的拟合结果;(d)标定波长的残差

Fig. 7. Calibration of relative wavelength. (a) Experimental setup; (b) measured data; (c) fitting result of wavelength calibration; (d) residual of calibration wavelength

下载图片 查看所有图片

4.3 基于多波段混合气体吸收谱的绝对波长标定实验

完成布里渊光谱仪相对波长的标定后,采用图8(a)所示的实验装置进行绝对波长的标定。实验装置主要包括布里渊光谱仪实验装置和用作标定的装有混合气体的气室、超连续光源等。气体特征吸收谱线具备非常高的波长精度和稳定性,是绝对波长标定的理想参考基准[14-15]。根据(10)式所示的朗伯-比尔定律得到的测量数据以及利用Voigt线型拟合后的结果如图8(b)所示。然后,基于特征谱线波峰的拟合值和标称值,利用(11)式计算标定系数并完成标定,结果如图8(c)所示。完成绝对波长定标后,将相对波长标定和绝对波长标定实验获得的标定系数,作为布里渊光谱仪测量光谱的波长标定系数,然后测量含有4种特征吸收峰的HCN标准气体的吸收光谱,比较测量光谱吸收峰值波长和基准值的偏差,结果如图8(d)所示。由图8(d)可见,最大的波长偏差约为0.025 pm。因此,可以认为本文所提方法的波长标定不确定度为±0.025 pm。考虑到气体特征谱线的波长不确定度约为±0.01 pm,可知布里渊光谱仪标定波长的不确定度约为±0.035 pm。

图 8. 绝对波长标定。(a)实验装置;(b)测量数据;(c)波长标定的拟合结果;(d)波长标定残差

Fig. 8. Calibration of absolute wavelength. (a) Experimental setup; (b) measured data; (c) fitting result of wavelength calibration; (d) residual of calibration wavelength

下载图片 查看所有图片

在完成布里渊光谱仪的波长标定后,分别对光纤光栅滤波器、正交调制电光调制器产生的调制信号进行测量,测得的光谱数据如图9所示。

图 9. 测量光谱。(a)光纤滤波器;(b)单边带调制信号

Fig. 9. Measured spectra. (a) Optical fiber filter; (b) single side-band modulation signal

下载图片 查看所有图片

5 结论

本文提出了一种基于FP标准具和气室的高精度波长标定方法,以满足新一代光通信测试所需的超高分辨率布里渊光谱分析仪的波长标定需求。本着多次标定逐步降低测量光谱波长误差的思路,首先基于布里渊光谱分析模型推导给出理论波长,然后基于FP标准具进行相对波长的标定,最后采用多波段混合气体吸收谱完成绝对波长的标定。实验结果表明,标定波长的不确定度可降低到±0.035 pm,实现了超高光谱分辨率、宽波段的布里渊光谱分析仪的波长标定。

此外,本文还利用受激布里渊增益谱和衰减谱的共生特性,对测量光谱中可能存在的鬼线进行了有效的识别校正。本文的波长定标方法对其他的超高分辨率光谱仪、可调谐光学滤波器、可调谐激光器等的波长标定,也具有一定的借鉴参考作用。

参考文献

[1] Al-Sammak M M A, Fyath R. Performance investigation of OFDM-based software-defined optical network[J]. International Journal of Networks and Communications, 2018, 8(1): 18-28.

[2] John C. Performance of coherent optical fiber transmission systems[J]. Frontiers of Optoelectronics, 2018, 11(2): 128-133.

[3] Budarnykh A, Vladimirskaya A, Lobach I, et al. High-resolution analysis of the wavelength stopping effect in a self-sweeping fiber laser[J]. Laser Physics Letters, 2019, 16(8): 085104.

[4] ManescauA, Araujo-HauckC, PasquiniL, et al. High resolution wavelength calibration: advancements with the laser frequency comb development[M]∥Astrophysics and Space Science Proceedings. Dordrecht: Springer Netherlands:411- 413.

[5] Rowe M A, Swann W C, Gilbert S L. Multiple-wavelength reference based on interleaved, sampled fiber Bragg gratings and molecular absorption[J]. Applied Optics, 2004, 43(17): 3530-3534.

[6] Duffin K. McGettrick A J, Johnstone W, et al. Tunable diode-laser spectroscopy with wavelength modulation: a calibration-free approach to the recovery of absolute gas absorption line shapes[J]. Journal of Lightwave Technology, 2007, 25(10): 3114-3125.

[7] Tateda M, Dong M H. Ferrule material dependence of axial force sensitivity of a tunable optical frequency filter made of fiber Fabry-Perot etalon[J]. Optical Review, 2011, 18(1): 39-41.

[8] Rivera E, Thomson D J. Accurate strain measurements with fiber Bragg sensors and wavelength references[J]. Smart Materials and Structures, 2006, 15(2): 325-330.

[9] Yeniay A, Delavaux J M, Toulouse J. Spontaneous and stimulated Brillouin scattering gain spectra in optical fibers[J]. Journal of Lightwave Technology, 2002, 20(8): 1425-1432.

[10] Fabrega J M, Sevillano P, Svaluto Moreolo M, et al. OFDM subcarrier monitoring using high resolution optical spectrum analysis[J]. Optics Communications, 2015, 342: 144-151.

[11] Swann W C, Gilbert S L. Pressure-induced shift and broadening of 1510-1540-nm acetylene wavelength calibration lines[J]. Journal of the Optical Society of America B, 2000, 17(7): 1263-1270.

[12] Gilbert S L, Swann W C, Dennis T. Wavelength standards for optical communications[J]. Proceedings of SPIE, 2001, 4269: 184-191.

[13] Rabin M W, Swann W C, Gilbert S L. Interleaved, sampled fiber Bragg gratings for use in hybrid wavelength references[J]. Applied Optics, 2002, 41(34): 7193-7196.

[14] 田鑫, 王蒙, 王泽锋. 基于倾斜光纤Bragg光栅的受激布里渊散射滤波器[J]. 光学学报, 2020, 40(10): 1006002.

    Tian X, Wang M, Wang Z F. Stimulated Brillouin scattering filters based on tilted fiber Bragg gratings[J]. Acta Optica Sinica, 2020, 40(10): 1006002.

[15] Gilbert SL, Etzel SM, Swann WC. Wavelength accuracy in WDM: techniques and standards for component characterization[C]∥Optical Fiber Communication Conference and Exhibit. 17-22 March 2002, Anaheim, CA, USA. New York: IEEE Press, 2002: 391- 393.

刘加庆, 刘磊, 刘雷, 李志增, 吴威, 胡垒军, 刘志明. 布里渊光谱仪的高精度波长标定方法研究[J]. 光学学报, 2020, 40(20): 2030001. Jiaqing Liu, Lei Liu, Lei Liu, Zhizeng Li, Wei Wu, Leijun Hu, Zhiming Liu. Wavelength Calibration of Ultra-High Resolution Brillouin Spectrometer[J]. Acta Optica Sinica, 2020, 40(20): 2030001.

本文已被 2 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!