多纵模激光器体制下直接多普勒测风技术的实现与仿真

上一篇

下一篇

高丰佳, 高飞, 赵婷婷, 汪丽, 李仕春, 闫庆, 华灯鑫. 多纵模激光器体制下直接多普勒测风技术的实现与仿真[J]. 物理学报, 2024, 73(20): 200701-1. doi: 10.7498/aps.73.20240949
引用本文: 高丰佳, 高飞, 赵婷婷, 汪丽, 李仕春, 闫庆, 华灯鑫. 多纵模激光器体制下直接多普勒测风技术的实现与仿真[J]. 物理学报, 2024, 73(20): 200701-1. doi: 10.7498/aps.73.20240949
Feng-Jia Gao, Fei Gao, Ting-Ting Zhao, Li Wang, Shi-Chun Li, Qing Yan, Deng-Xin Hua. Implementation and simulation of direct Doppler wind measurement technology under regime of multi-longitudinal mode laser[J]. Acta Physica Sinica, 2024, 73(20): 200701-1. doi: 10.7498/aps.73.20240949
Citation: Feng-Jia Gao, Fei Gao, Ting-Ting Zhao, Li Wang, Shi-Chun Li, Qing Yan, Deng-Xin Hua. Implementation and simulation of direct Doppler wind measurement technology under regime of multi-longitudinal mode laser[J]. Acta Physica Sinica, 2024, 73(20): 200701-1. doi: 10.7498/aps.73.20240949

多纵模激光器体制下直接多普勒测风技术的实现与仿真

    作者简介: 高丰佳. E-mail: 1093788469@qq.com .
    通讯作者: E-mail: gaofei@xaut.edu.cn.; 
  • 中图分类号: 07.05.Fb, 07.05.Tp, 42.79.Qx, 42.68.Wt

Implementation and simulation of direct Doppler wind measurement technology under regime of multi-longitudinal mode laser

    Corresponding author: E-mail: gaofei@xaut.edu.cn.; 
  • MSC: 07.05.Fb, 07.05.Tp, 42.79.Qx, 42.68.Wt

  • 摘要: 考虑到单纵模激光器作为多普勒测风激光雷达的激励光源需要复杂的种子注入技术和高精度的稳频锁频技术, 本文提出了以自由运转的多纵模激光器作为激励光源的直接多普勒测风激光雷达, 以降低激励光源的实现难度和复杂性. 针对典型Nd:YAG脉冲激光器, 研究了不同激光辐射线宽、光学谐振腔长和中心波长条件下多纵模激光激励的大气弹性散射回波光谱的分布模式. 为了综合利用大气风场导致的多纵模激光中每个纵模所激励大气弹性散射回波光谱的多普勒频移, 利用具有周期性透过率曲线且四个输出通道相位各相差π/2的可调谐四通道马赫-曾德尔干涉仪, 实现对多纵模大气弹性散射回波光谱多普勒频移的高精度鉴别. 在此基础上, 构建了多纵模激光器体制下直接多普勒测风技术大气风场反演的数学模型. 仿真结果表明, 所提出的多纵模直接多普勒测风激光雷达能够实现对大气风场信息的高精度测量, 并且激光中心波长越大, 激光光学谐振腔长越短, 系统测风范围越大, 测风误差越小.
  • 加载中
  • 图 1  多纵模激光激励的大气弹性散射回波光谱示意图 (a) 激光辐射线宽为30 GHz, 激光光学谐振腔长为30 mm; (b) 激光辐射线宽为60 GHz, 激光光学谐振腔长为30 mm; (c) 激光辐射线宽为30 GHz, 激光光学谐振腔长为300 mm; (d) 激光辐射线宽为60 GHz, 激光光学谐振腔长为300 mm

    Figure 1.  Diagram of atmospheric elastic scattering echo spectra excited by multi-longitudinal mode laser: (a) Laser radiation linewidth is 30 GHz, laser optical resonator length is 30 mm; (b) laser radiation linewidth is 60 GHz, laser optical resonator length is 30 mm; (c) laser radiation linewidth is 30 GHz, laser optical resonator length is 300 mm; (d) laser radiation linewidth is 60 GHz, laser optical resonator length is 300 mm.

    图 2  不同激光中心波长时, 交叠系数随激光光学谐振腔长的变化曲线 (a) 多纵模气溶胶米散射回波光谱的交叠系数; (b) 多纵模大气分子瑞利散射回波光谱的交叠系数

    Figure 2.  Change curve of overlap coefficient with laser optical resonator length under different laser center wavelengths: (a) Overlap coefficient of multi-longitudinal mode aerosol Mie scattering echo spectra; (b) overlap coefficient of multi-longitudinal mode molecular Rayleigh scattering echo spectra

    图 3  直接单边缘检测技术鉴别多纵模大气弹性散射回波光谱多普勒频移的原理图

    Figure 3.  Principle diagram of direct single-edge detection technology for measuring the Doppler shift of atmospheric elastic scattering echo spetra excited by multi-longitudinal mode laser.

    图 4  基于可调谐四通道马赫-曾德尔干涉仪的多纵模直接多普勒测风激光雷达系统图. BS: 光束分离器; E: 准直扩束器; M: 反射镜; F: 光纤; CL准直透镜; QWP: 四分之一波片; DS: 纳米位移系统; GP: 格兰棱镜; L: 汇聚透镜; D: 光电探测器件

    Figure 4.  System diagram of multi-longitudinal mode direct Doppler wind lidar based on quadri-channel Mach-Zehnder interferometer. BS: beam splitter; E: expander; M: mirror; F: fiber; CL: collimating lens; QWP: quarter-wave plate; DS: displacement system; GP: Glan prism; L: convergent lens; D: photoelectric detector.

    图 5  四通道马赫-曾德尔干涉仪鉴别多纵模大气弹性散射回波光谱多普勒频移的原理图

    Figure 5.  Principle diagram of quadri-channel Mach-Zehnder interferometer for measuring Doppler shift of atmospheric elastic scattering echo spetra excited by multi-longitudinal mode laser.

    图 6  四通道马赫-曾德尔干涉仪鉴别多纵模大气弹性散射回波光谱多普勒频移的输出信号 (a) QMZI通道1输出信号; (b) QMZI通道2输出信号; (c) QMZI通道3输出信号; (d) QMZI通道4输出信号

    Figure 6.  Output signals of quadri-channel Mach-Zehnder interferometer for measuring Doppler shift of atmospheric elastic scattering echo spetra excited by multi-longitudinal mode laser: (a) Output signals of channel 1; (b) output signals of channel 2; (c) output signals of channel 3; (d) output signals of channel 4.

    图 7  激光光学谐振腔长为30 mm时, 不同激光中心波长下QMZI有效透过率随风速变化的仿真曲线 (a) 激光中心波长为1064 nm; (b) 激光中心波长为532 nm; (c) 激光中心波长为355 nm

    Figure 7.  Simulation curve of effective transmittance of QMZI with wind speed at different laser center wavelengths when the laser optical resonator length is 30 mm: (a) The laser center wavelength is 1064 nm; (b) the laser center wavelength is 532 nm; (c) the laser center wavelength is 355 nm.

    图 8  激光光学谐振腔长为300 mm时, 不同激光中心波长下QMZI有效透过率随风速变化的仿真曲线 (a) 激光中心波长为1064 nm; (b) 激光中心波长为532 nm; (c) 激光中心波长为355 nm

    Figure 8.  Simulation curve of effective transmittance of QMZI with wind speed at different laser center wavelengths when the laser optical resonator length is 300 mm: (a) The laser center wavelength is 1064 nm; (b) the laser center wavelength is 532 nm; (c) the laser center wavelength is 355 nm.

    图 9  激光光学谐振腔长为30 mm时, 不同激光中心波长下多纵模直接多普勒测风激光雷达系统的测风误差随风速变化的仿真曲线 (a) 激光中心波长为1064 nm; (b) 激光中心波长为532 nm; (c) 激光中心波长为355 nm

    Figure 9.  Simulation curve of wind measurement error of multi-longitudinal mode direct Doppler wind lidar with wind speed at different laser center wavelengths when the laser optical resonator length is 30 mm: (a) The laser center wavelength is 1064 nm; (b) the laser center wavelength is 532 nm; (c) the laser center wavelength is 355 nm.

    图 10  激光光学谐振腔长为300 mm时, 不同激光中心波长下多纵模直接多普勒测风激光雷达系统的测风误差随风速变化的仿真曲线 (a) 激光中心波长为1064 nm; (b) 激光中心波长为532 nm; (c) 激光中心波长为355 nm

    Figure 10.  Simulation curve of wind measurement error of multi-longitudinal mode direct Doppler wind lidar with wind speed at different laser center wavelengths when the laser optical resonator length is 300 mm: (a) The laser center wavelength is 1064 nm; (b) the laser center wavelength is 532 nm; (c) the laser center wavelength is 355 nm.

    表 1  多纵模激光器体制下直接多普勒测风激光雷达关键参数.

    Table 1.  Key parameters of multi-longitudinal mode direct Doppler wind lidar.

    参数 参数
    激光中心波长/nm 1064, 532, 355 理论风速/(m·s–1) –50—50 (间隔1 m·s–1)
    激光辐射线宽/GHz 30 大气温度/K 300
    激光光学谐振腔长/mm 30, 300 气溶胶后向散射比 2
    激光纵模个数 7, 61 QMZI自由光谱范围/GHz 5.0, 0.5
    激光纵模间隔/GHz 5.0, 0.5
    下载: 导出CSV
  • [1] Kumar D, Premachandran B 2019 Int. J. Therm. Sci. 138 263 doi: 10.1016/j.ijthermalsci.2018.12.010
    [2] Liu G Q, Perrie W 2013 Geophys. Res. Lett. 40 3150 doi: 10.1002/grl.50624
    [3] Gardiner B, Berry P, Moulia B 2016 Plant Sci. 245 94 doi: 10.1016/j.plantsci.2016.01.006
    [4] Yu L J, Zhong S Y, Bian X D, Heilman W E 2018 Int. J. Climatol. 39 1684
    [5] 马福民, 陈涌, 杨泽后, 周鼎富, 李晓锋, 陈春利, 冯力天, 余臣 2019 激光与光电子学进展 56 180003 doi: 10.3788/LOP56.180003 Ma F M, Chen Y, Yang Z H, Zhou D F, Li X F, Chen C L, Feng L T, Yu C 2019 Laser Optoelectron. Prog. 56 180003 doi: 10.3788/LOP56.180003
    [6] Reitebuch O, Lemmerz Ch, Nagel E, Paffrath U, Durand Y, Endemann M, Fabre F, Chaloupy M 2009 J. Atmos. Oceanic Technol. 26 2501 doi: 10.1175/2009JTECHA1309.1
    [7] Paffrath U, Lemmerz Ch, Reitebuch O, Oliver, Witschas B, Nikolaus I, Freudenthaler V 2009 J. Atmos. Oceanic Technol. 26 2516 doi: 10.1175/2009JTECHA1314.1
    [8] 储玉飞, 刘东, 王珍珠, 吴德成, 邓迁, 李路, 庄鹏, 王英俭 2020 量子电子学报 37 580 Chu Y F, Liu D, Wang Z Z, Wu D C, Deng Q, Li L, Zhuang P, Wang Y J 2020 Chin. J. Quantum Electron. 37 580
    [9] Jiang S, Sun D S, Han Y L, Han F, Zhou A R, Zheng J 2019 Curr. Opt. Photonics 3 466
    [10] Zhang Y P, Yuan J L, Wu Y B, Dong J J, Xia H Y 2023 Phys. Rev. Fluids 8 L022701 doi: 10.1103/PhysRevFluids.8.L022701
    [11] Zhang Y P, Wu Y B, Dong J J, Xia H Y 2022 IEEE Photonics J. 14 6047706
    [12] Liu Z L, Barlow J F, Chan P W, Fung J C H, Li Y G, Ren C, Mark H W L, Ng E 2019 Remote Sens. 11 2522 doi: 10.3390/rs11212522
    [13] 张亚飞, 冯玉涛, 傅頔, 畅晨光, 李娟, 白清兰, 胡炳樑 2022 物理学报 71 084201 doi: 10.7498/aps.71.20212086 Zhang Y F, Feng Y T, Fu D, Chang C G, Li J, Bai Q L, Hu B J 2022 Acta Phys. Sin. 71 084201 doi: 10.7498/aps.71.20212086
    [14] Vrancken P, Herbst J 2022 Remote Sens. 14 3356 doi: 10.3390/rs14143356
    [15] Kliebisch O, Uittenbosch H, Thurn J, Mahnke P 2022 Opt. Express 30 5540 doi: 10.1364/OE.445287
    [16] Wang L, Gao F, Wang J, Yan Q, Yan W X, Wang M, Hua D X 2019 Opt. Laser Eng. 121 61 doi: 10.1016/j.optlaseng.2019.03.020
    [17] Hill C 2018 Remote Sens. 10 466 doi: 10.3390/rs10030466
    [18] Shen F H, Wang B X, Shi W J, Zhuang P, Zhu C Y, Xie C B 2018 Opt. Commun. 412 7 doi: 10.1016/j.optcom.2017.11.073
    [19] 潘艺升, 闫召爱, 郭文杰, 徐轻尘, 胡雄 2016 激光技术 40 153 Pan Y S, Yan Z A, Guo W J, Xu Q C, Hu X 2016 Laser Technol. 40 153
    [20] Wu C T, Chen F, Dai T Y, Ju Y L 2015 J. Mod. Opt. 62 1535 doi: 10.1080/09500340.2015.1049232
    [21] 张明富, 杨天新, 葛春风 2022 红外与激光工程 51 20210435 doi: 10.3788/IRLA20210435 Zhang M F, Yang T X, Ge C F 2022 Infrared Laser Eng. 51 20210435 doi: 10.3788/IRLA20210435
    [22] 葛烨, 胡以华, 舒嵘, 洪光烈 2015 物理学报 64 020707 doi: 10.7498/aps.64.020702 Ge Y, Hu Y H, Shu R, Hong G L 2015 Acta Phys. Sin. 64 020707 doi: 10.7498/aps.64.020702
    [23] Bruneau D, Blouzon F, Spatazza J, Montmessin F, Pelon J, Faure B 2013 Appl. Opt. 52 4941 doi: 10.1364/AO.52.004941
    [24] 高飞, 南恒帅, 黄波, 汪丽, 李仕春, 王玉峰, 刘晶晶, 闫庆, 宋跃辉, 华灯鑫 2018 物理学报 67 030701 doi: 10.7498/aps.67.20172036 Gao F, Nan H S, Huang B, Wang L, Li S C, Wang Y F, Liu J J, Yan Q, Song Y H, Hua D X 2018 Acta Phys. Sin. 67 030701 doi: 10.7498/aps.67.20172036
    [25] Gao F, Nan H S, Zhang R, Zhu Q S, Chen T, Wang L, Chen H, Hua D X, Stanic S 2019 J. Quant. Spectrosc. Radiat. Transfer 234 10 doi: 10.1016/j.jqsrt.2019.05.021
    [26] 毛艳丽, 邱宏伟, 徐军, 邓佩珍, 干福熹 2001 光学学报 21 1264 Mao Y L, Qiu H W, Xu J, Deng P Z, Gan F X 2001 Acta Opt. Sin. 21 1264
    [27] Korb C L, Gentry B M, Weng C Y 1992 Appl. Opt. 31 4202 doi: 10.1364/AO.31.004202
    [28] Thompson B J, Wolf E 1957 J. Opt. Soc. Am. 47 895 doi: 10.1364/JOSA.47.000895
  • 加载中
图( 10) 表( 1)
计量
  • 文章访问数:  246
  • HTML全文浏览数:  246
  • PDF下载数:  2
  • 施引文献:  0
出版历程
  • 收稿日期:  2024-07-09
  • 刊出日期:  2024-10-20

多纵模激光器体制下直接多普勒测风技术的实现与仿真

    通讯作者: E-mail: gaofei@xaut.edu.cn.; 
    作者简介: 高丰佳. E-mail: 1093788469@qq.com
  • 1. 西安理工大学机械与精密仪器工程学院, 西安 710048
  • 2. 陕西省现代装备绿色制造协同创新中心, 西安 710048
  • 3. 国家市场监管重点实验室(计量光学及应用), 西安 710048

摘要: 考虑到单纵模激光器作为多普勒测风激光雷达的激励光源需要复杂的种子注入技术和高精度的稳频锁频技术, 本文提出了以自由运转的多纵模激光器作为激励光源的直接多普勒测风激光雷达, 以降低激励光源的实现难度和复杂性. 针对典型Nd:YAG脉冲激光器, 研究了不同激光辐射线宽、光学谐振腔长和中心波长条件下多纵模激光激励的大气弹性散射回波光谱的分布模式. 为了综合利用大气风场导致的多纵模激光中每个纵模所激励大气弹性散射回波光谱的多普勒频移, 利用具有周期性透过率曲线且四个输出通道相位各相差π/2的可调谐四通道马赫-曾德尔干涉仪, 实现对多纵模大气弹性散射回波光谱多普勒频移的高精度鉴别. 在此基础上, 构建了多纵模激光器体制下直接多普勒测风技术大气风场反演的数学模型. 仿真结果表明, 所提出的多纵模直接多普勒测风激光雷达能够实现对大气风场信息的高精度测量, 并且激光中心波长越大, 激光光学谐振腔长越短, 系统测风范围越大, 测风误差越小.

English Abstract

    • 风是太阳辐射引起的一种空气流动现象, 其与地球上的大气对流和海陆循环等自然规律息息相 关[1,2]. 高精度探测大气风场信息, 不论是在气候学和气象学等理论研究方面, 还是在环境监测、航空安全、风能利用、农业生产等国民经济建设方面都具有非常重要的意义[3,4]. 近年来, 多普勒测风激光雷达因其探测精度高、探测距离大、可通过布控实现大范围风场监测而被广泛应用于大气风场信息探测[57]. 多普勒测风激光雷达是基于光的多普勒效应, 利用随大气风场运动的气溶胶和大气分子对激光信号的弹性散射作用获得多普勒频移, 进而实现对大气风场信息的测量[8].

      多普勒测风激光雷达根据激励光源的不同可以分为连续波多普勒激光雷达[9]和脉冲多普勒激光雷达[10]. 连续波多普勒激光雷达的距离信息通过几何方法获得, 其最大探测距离通常小于1 km[11]. 而脉冲多普勒激光雷达的距离信息由激光脉冲的飞行时间来确定, 其最大探测距离较远, 通常可以达到几公里到几十公里[12]. 多普勒测风激光雷达根据鉴频技术的不同可以分为外差多普勒测风激光雷达[13]和直接多普勒测风激光雷达[14]. 外差和直接探测技术都是通过测量激光发射信号和大气回波信号之间的多普勒频移来实现风场信息反演的. 不同的是, 外差探测技术是利用激光的相干性质和电子信号处理技术解调激光发射信号和大气回波信号之间的差频信号来获得多普勒频移[15], 而直接探测技术是利用光谱鉴别器鉴别气溶胶米散射和大气分子瑞利散射回波信号的多普勒频移所引起的能量变化, 再反演获得多普勒频移[16]. 相比之下, 外差探测技术虽然可以获得更高的探测精度, 但同时也对光学准直性有更高的要求, 并且其只能实现对窄带气溶胶米散射回波光谱多普勒频移的鉴别[17]. 而直接探测的边缘检测技术对于鉴别窄带气溶胶米散射回波光谱多普勒频移和宽带大气分子瑞利散射回波光谱多普勒频移都比较有效[18]. 无论是基于哪种探测技术的多普勒测风激光雷达, 都需要通过复杂的种子注入技术[19,20]和高精度的稳频、锁频技术[21,22]来提供高稳定度的单纵模激光作为激励光源, 这极大地增加了多普勒测风激光雷达系统的复杂性和技术难度.

      考虑到单纵模激光作为激励光源的实现难度, 2013年Bruneau等[23]提出了直接采用自由运转的多纵模激光作为多普勒测风激光雷达系统的激励光源. 多纵模激光光谱可以看作是多个单纵模激光光谱在频域上的等距排列分布, 因此多纵模激光与单纵模激光激励的大气弹性散射回波光谱分布存在较大差异[24,25]. 为了高精度测量多纵模激光 激励的大气弹性散射回波光谱的多普勒频移, 本文探讨了不同激光辐射线宽、光学谐振腔长、中心波长条件下的多纵模激光激励的大气弹性散射回波光谱的分布特性. 在此基础上, 提出了一种基于可调谐四通道马赫-曾德尔干涉仪(quadri-Mach-Zender interferometer, QMZI)的多纵模直接多普勒测风激光雷达系统, 分析了QMZI鉴别多纵模激光激励的大气弹性散射回波光谱多普勒频移的技术原理, 构建了多纵模光谱模式下大气风场的数据反演模型, 并对所提出的多纵模直接多普勒测风激光雷达系统进行了仿真验证.

    • Nd:YAG激光器增益介质中钕浓度较高[26], 其单位体积的工作物质能提供比较高的激光功率, 且其通过倍频技术可以提供1064 nm(基频)、532 nm(二倍频)和355 nm(三倍频)三种中心波长的激光, 因此被广泛地应用于激光雷达系统中. 自由运转的Nd: YAG激光器输出的激光具有多纵模光谱模式, 其可以看作是多个单纵模激光光谱在频域上的等距排列分布[24,25]. 激光与大气颗粒发生相互作用会产生大气散射回波信号, 通常利用大气弹性散射(气溶胶米散射和大气分子瑞利散射)回波信号来分析反演大气风场信息. 多纵模激光激励的大气弹性散射回波光谱可以表示为[24,25]

      式中, ν0代表激光中心频率; βa, βm分别代表气溶胶米散射和大气分子瑞利散射的后向散射系数; Ra代表气溶胶后向散射比. MSa(νν0)和MSm(νν0)分别代表多纵模激光激励的气溶胶米散射和大气分子瑞利散射回波光谱, 其受到激光辐射线宽、光学谐振腔长、中心波长和大气温度的影响[24,25]. 而MS(νν0)代表由气溶胶米散射和大气分子瑞利散射回波光谱叠加获得的多纵模大气弹性散射回波光谱, 因此其与激光辐射线宽、光学谐振腔长、中心波长和大气温度、气溶胶后向散射比均相关. 为了探讨激光器自身参数与其光谱的关系, 分析了激光辐射线宽Δνe为30, 60 GHz, 光学谐振腔长L为30, 300 mm, 中心波长λ0为1064, 532, 355 nm, 大气温度T为300 K, 气溶胶后向散射比Ra为2时, 多纵模激光激励的大气弹性散射回波光谱如图1所示.

      图1的四幅图中均可以看出, 多纵模激光激励的大气弹性散射回波光谱可以看作是激光辐射线宽内所有激光纵模激励的气溶胶米散射和大气分子瑞利散射回波光谱有序叠加的结果. 激光中心波长越小, 则单纵模大气分子瑞利散射回波光谱的谱宽越宽. 对比图1(a)图1(b)图1(c)图1(d)可以看出, 激光光学谐振腔长相同的情况下, 多纵模激光激励的大气弹性散射回波光谱的纵模间隔相同, 而此时纵模个数与辐射线宽相关, 辐射线宽越大, 则纵模个数越多. 对比图1(a)图1(c)图1(b)图1(d)可以看出, 激光辐射线宽相同的情况下, 多纵模激光激励的大气弹性散射回波光谱的纵模个数和纵模间隔均受激光光学谐振腔长影响, 激光光学谐振腔长越长, 则纵模个数越多, 纵模间隔越小. 从图1的四幅图中还可以看出, 当多纵模激光激励的大气弹性散射回波光谱的纵模个数较多、纵模间隔较小时, 相邻纵模之间会相互交叠, 而相邻纵模之间的相互交叠会影响多普勒激光雷达鉴别多纵模大气弹性散射回波光谱多普勒频移的性能. 本文采用交叠系数来描述多纵模激光激励的大气弹性散射回波光谱相邻纵模之间的交叠程度, 其可以用下式定义:

      式中, OaOm代表多纵模气溶胶米散射和大气分子瑞利散射回波光谱的交叠系数; Δνa和Δνm代表单纵模气溶胶米散射和大气分子瑞利散射回波光谱的半高全宽; Δνq代表纵模间隔; kB代表玻尔兹曼常数; T代表大气温度; M代表平均分子质量; c代表光速; L代表光学谐振腔长; λ0代表激光中心波长. 图2为不同激光中心波长时, 交叠系数随激光光学谐振腔长的变化曲线. 若交叠系数≤0.5, 则说明多纵模激光激励的大气弹性散射回波光谱相邻纵模之间的交叠较弱, 其对多普勒激光雷达鉴频性能的影响几乎可以忽略; 若交叠系数>0.5, 则说明多纵模激光激励的大气弹性散射回波光谱相邻纵模之间的交叠较明显, 其会导致多普勒激光雷达鉴频性能明显下降.

      图2可以看出, 多纵模气溶胶米散射回波光谱的交叠系数仅与激光光学谐振腔长相关, 其随激光光学谐振腔长的增大而增大, 但直到激光光学谐振腔长为750 mm(对应的纵模个数为151)时, 交叠系数才等于0.5, 即此时多纵模气溶胶米散射回波光谱的相邻纵模之间的交叠才开始影响多普勒激光雷达的鉴频性能. 而多纵模大气分子瑞利散射回波光谱的交叠系数与激光光学谐振腔长和中心波长均相关, 其随激光光学谐振腔长的增大而增大, 随激光中心波长的增大而减小. 激光中心波长分别为1064, 532, 355 nm时, 在激光光学谐振腔长分别为57, 29, 19 mm (对应的纵模个数为11, 7, 5)的位置, 交叠系数等于0.5, 即此时多纵模大气分子瑞利散射回波光谱的相邻纵模之间的交叠开始影响多普勒激光雷达的鉴频性能. 因此, 当激光中心波长分别为1064, 532, 355 nm时, 纵模个数较少时(N < 11, 7, 5时), 多纵模气溶胶米散射和大气分子瑞利散射回波光谱的多普勒频移均可被鉴别, 进而反演大气风场信息. 当纵模个数较多时(N > 11, 7, 5时), 仅多纵模气溶胶米散射回波光谱的多普勒频移可以被鉴别去反演大气风场信息.

    • 考虑到大气风场会给多纵模激光中每个纵模激励的大气弹性散射回波光谱都引入多普勒频移, 因此最优的光谱鉴频器需要选择具有周期性透过率曲线的光学干涉仪, 并且参数需要满足多纵模激光的纵模间隔(Δνq)与光谱鉴频器透过率曲线的自由光谱范围(FSR)相匹配的条件. 满足上述条件的基础上, 借鉴单纵模光谱模式下直接单边缘检测技术[27], 可以给出直接单边缘检测技术鉴别多纵模激光激励的大气弹性散射回波光谱多普勒频移的原理图, 如图3所示, 图中以激光辐射线宽Δνe为30 GHz, 光学谐振腔长L为30 mm, 中心波长λ0为532 nm, 大气温度T为300 K, 气溶胶后向散射比Ra为2, 纵模个数为7时的多纵模激光激励的大气弹性散射回波光谱为例说明了多纵模模式下回波光谱多普勒频移的鉴别原理, 其他多纵模模式下鉴别原理是相同的.

      图3中可以看出, 多纵模激光的中心频率被锁定在光谱鉴频器透过率曲线上升沿的最陡峭边缘处[27], 由于多纵模激光的纵模间隔和光谱鉴频器透过率曲线的自由光谱范围是相匹配的, 因此多纵模激光激励的大气弹性散射回波光谱中每个纵模的峰值频率都被锁定在了光谱鉴频器透过率曲线上升沿的最陡峭边缘处. 当多纵模激光激励的大气弹性散射回波光谱无多普勒频移时, 大气回波光谱的中心频率与发射激光中心频率一致, 光谱鉴频器输出端的能量不发生变化; 而当多纵模激光激励的大气弹性散射回波光谱存在多普勒频移时, 大气回波光谱的中心频率将会发生偏移, 光谱鉴频器输出端的能量发生变化, 因此通过此能量变化可以计算得到多普勒频移, 进而反演出大气风场信息.

      法布里-珀罗干涉仪和迈克耳孙干涉仪都可以提供单通道的周期性透过率曲线, 以实现对多纵模激光激励的大气弹性散射回波光谱多普勒频移的单边缘检测. 马赫-曾德尔干涉仪可以提供具有互补输出的双通道周期性透过率曲线, 实现对多纵模回波光谱多普勒频移的双边缘检测. 而四通道马赫-曾德尔干涉仪可以提供四个周期相同、相位各相差π/2的周期性透过率曲线, 能够以更高的探测灵敏度实现对多纵模激光激励的大气弹性散射回波光谱多普勒频移的鉴别.

    • 为了实现对多纵模激光激励的大气散射回波光谱多普勒频移的高精度鉴别, 设计了如图4所示的基于可调谐四通道马赫-曾德尔干涉仪的多纵模直接多普勒测风激光雷达(multi-longitudinal mode direct Doppler wind lidar, MLM DDWL)系统, 其利用自由运转的大功率Nd: YAG脉冲激光器向大气中传输多纵模激光. Nd: YAG脉冲激光器发出的多纵模激光被光束分离器(beam separator, BS)3分为两路, 其中一路通过光纤(fiber, F)传输至望远镜作为参考信号; 另一路经准直扩束器(expander, E)后入射到大气中与气溶胶和大气分子发生相互作用并产生多纵模气溶胶米散射和大气分子瑞利散射回波信号, 其中后向散射信号被垂直放置的大口径望远镜接收. 然后, 回波信号通过光纤传输至具有四通道周期性透过率曲线的QMZI来进行多普勒频移的高精度鉴别.

      图4中红色实线框标出的部分即为可调谐光程差的QMZI的结构示意图, 其四通道输出是由光束分离器、格兰棱镜(Glan prism, GP)、四分之一波片(quarter-wave plate, QWP)和纳米位移系统(displacement system, DS)等核心器件来实现的. BS可以将信号以1:1的比例分为反射信号和透射信号, 并且能够为介质涂层方向入射信号的反射信号引入π的相移. GP可以将信号以1:1的比例分为寻常光信号和非常光信号, 而QWP可以为寻常光信号和非常光信号引入π/2 (或π/2奇数倍)的相位差. 纳米位移系统可以实现对QMZI光程差在纳米级别上的调整.

      当大气弹性后向散射回波信号被望远镜接收并通过光纤传输至QMZI时, 其首先被准直透镜(collimating lens, CL)准直扩束, 然后被BS1分束. BS1的反射信号通过光路①, 其经过反射镜(mirror, M)M1反射、QWP透射、M2反射后到达BS2. BS1的透射信号通过光路②, 其经过经安装在DS上的M3与M4反射后到达BS2. 两路信号在BS2处汇聚并发生干涉现象, 进而产生两路干涉信号. 而这两路干涉信号又分别到达GP1和GP2, 并被其分为寻常光信号和非常光信号, 因此得到四路相位各相差π/2的干涉信号输出. 四路干涉信号经过四个汇聚透镜(convergent lens, L)分别聚焦、四个光电探测器件(detector, D)分别接收和光电转换, 最后传输至示波器进行数据接收与处理.

      为了分析QMZI的透过率函数, 假设入射至QMZI的信号为单频线偏振光, 电矢量的振幅为A, 辐射频率为ν, 初始相位为$\varphi $, 则其电场强度分布E0(ν, t)可以表示为

      式中, ν为光谱频率. 则入射信号的干涉光强度分布I0(ν)可以表示为

      从(4)式中可以看出, 干涉光强度分布I0(ν)中包含了原频率项和倍频项的分量. 通过对光电倍增管型号的选择与输出电路的设计可以使得光电倍增管不响应干涉光强度分布I0(ν)中倍频项的分量, 因此实际有I0(ν) = A2/2.

      假设Eij(ν, t)代表QMZI入射信号以不同光路到达不同通道时的电场强度, 其下标i = 1, 2 表示通过的光路, j = 1, 2, 3, 4表示到达的通道. 通过光路①到达D1的非常光信号, 经过BS1反 射(光线由介质涂层侧入射,反射信号相位延迟π)、M1反射(整体半波损失)、QWP透射(非常光信号无影响, 寻常光信号相位延迟π/2)、M2反射(整体半波损失)、BS2透射(光线由介质涂层侧入射,透射信号无相移)、GP1透射,相位损失为3π, 则此部分信号电场强度分布E11(ν, t)为

      RBS1, RM1, TQWP, RM2, TBS2, TGP1分别代表BS1的反射率、M1的反射率、QWP的透射率、M2的反射率、BS2的透射率、GP1的透射率, d1代表 光路①的光学长度. 按照光路②到达D1的非常光信号, 经过BS1透射 (光线由介质涂层侧入射, 透射信号无相移)、M3反射 (整体半波损失)、M4反射 (整体半波损失)、BS2反射 (光线非介质涂层侧入射,透射与反射时均无相移)、GP1透射,其相位损失为2π,则此部分信号电场强度分布E21(ν, t)为

      式中, TBS1, RM3, RM4, RBS2分别代表BS1的透射率、M3的反射率、M4的反射率、BS2的反射率; d2代表光路②的光学长度. 按照不同光路到达D1的两路光线相位差为π, 则D1接收信号的干涉光强度分布I1(ν)为

      因此, D1通道的透过率函数T1(ν)为

      忽略光路损耗和光学元件瑕疵时, D1通道的理想透过率函数为

      同理可以得到其他通道的理想透过率函数, 因此QMZI四通道的理想透过率函数Ti(ν) (i = 1, 2, 3, 4)为

      式中, αi (i = 1, 2, 3, 4)表示QMZI透过率函数中三角项的附加相位因子.

      图5图6分别为四通道马赫-曾德尔干涉仪鉴别多纵模激光激励的大气弹性散射回波光谱多普勒频移的原理图与输出信号. 图中提供了无频移和具有高频方向多普勒频移的大气弹性散射回波光谱.

      图5可以看出, 多纵模激光每个纵模的峰 值频率都被锁定在了QMZI通道2透过率曲线的上升沿的最陡峭边缘频率处. 由于QMZI四个 通道的透过率函数相位各相差π/2, 因此多纵模激光每个纵模的峰值频率同样与QMZI通道1透过率曲线的谷值频率、通道3透过率曲线的峰值 频率、通道4透过率曲线下降沿的最陡峭边缘频率相对应.

      图6可以看出, 无频移的大气弹性散射回波光谱通过QMZI后, 通道1和通道3的输出信号分别为四个通道中的最小值和最大值, 通道2和通道4的输出信号是互补的且其强度相等. 相比于无频移的大气弹性散射回波光谱, 具有高频方向多普勒频移的大气弹性散射回波光谱通过QMZI后, 四个通道的输出信号分别增大、增大、减小、减小. 因此QMZI四个通道输出信号的变化情况正好反映了多普勒频移的变化, 以及与之相对应的大气风场信息.

    • 当入射至QMZI的信号为多纵模激光激励的大气弹性散射回波信号时, QMZI对该光谱的有效透过率可以由QMZI透过率函数和多纵模激光激励的大气弹性散射回波光谱卷积获得

      式中, ETiν)(i = 1, 2, 3, 4)代表QMZI对带有多普勒频移的多纵模激光激励的大气弹性散射回波信号的有效透过率; Δν代表多普勒频移; MS(ν–(v0ν))代表带有多普勒频移的多纵模激光激励的大气弹性散射回波光谱分布. 由于(11)式中的卷积计算极其复杂, 而多纵模激光及其激励大气弹性散射回波光谱均属于准单色光[24,25], 因此, 根据准单色光干涉的部分相干理论与光波偏振效应[28], 多纵模大气弹性散射回波信号入射至QMZI后获得的四路干涉信号可以表示为

      式中, Io1Io2代表GP1和GP2反射的寻常光干涉光强分布; Ie1Ie2代表GP1和GP2透过的非常光干涉光强分布; ItoIro代表BS1透过和反射的寻常光干涉光强分布; IteIre代表BS1透过和反射的非常光干涉光强分布; |γM|与αM分别表示多纵模大气弹性散射回波信号的复相干度的模数因子与相位因子, 可以表示为

      式中, |γMa|和|γMm|分别表示多纵模气溶胶米散射和多纵模大气分子瑞利散射回波信号的复相干度的模数因子, 其可由下式求得

      由于GP1与GP2的分光比均为1∶1, 且四个光电探测器件接收到的干涉光强I1ν), I2ν), I3ν), I4ν)分别对应Ie1, Io1, Ie2, Io2, 因此QMZI四个通道对带有多普勒频移的多纵模大气弹性散射回波信号的有效透过率可以求出:

      借鉴双边缘检测技术的鉴频原理, 可以将QMZI鉴频的响应函数定义为

      式中, Qν)代表QMZI的响应函数. 对其进行反演可以得到多纵模大气弹性散射回波信号的多普勒频移为

      进而通过多普勒频移和多普勒测风技术原理[8]可以反演得到大气风场信息:

      式中, V代表激光束方向的风速.

    • 为了验证多纵模激光器体制下基于QMZI的直接多普勒测风激光雷达的可行性与探测性能, 本文采用可以对光电探测系统进行分析的Zemax光学仿真软件进行测风仿真实验. 在Zemax非序列模式下搭建如图4中红色实线框出所示的QMZI仿真模型, 自定义Zemax高斯光源中的SPCD光谱文件以构建直接多普勒测风激光雷达中多纵模激光激励的大气弹性散射回波光谱分布. 本文按照表1的关键参数配置了6组不同风速下的多纵模大气弹性散射回波信号的SPCD光谱文件. 需要注意的是, 由于Zemax光学仿真软件的限制, 用于配置光源信号的SPCD光谱文件中离散点不能超过200个. 此外, 风速值为正时代表风速方向为靠近光源方向, 风速值为负时代表风速方向为远离光源方向.

      将配置的不同激光中心波长(λ = 1064 nm, λ = 532 nm, λ = 355 nm)、不同激光光学谐振腔长(L = 30 mm, L = 300 mm)的SPCD光谱文件按照风速的大小顺序分别依次输入至QMZI仿真模型中, 然后在Zemax中进行基于蒙特卡罗原理的光线追踪, 最终获得6组QMZI四个通道的输出信号. 将这6组输出信号分别代入(17)式即可获的QMZI对带有多普勒频移的多纵模大气弹性散射回波信号的有效透过率(如图7图8所示), 进而通过(18)式可以获得仿真测量的风速值. 将仿真测量风速与理论风速作差, 即可求得多纵模激光器模式下基于QMZI的直接多普勒测风雷达的测风误差, 如图9图10所示.

      图7图8可以看出, 不同激光中心波长、不同激光光学谐振腔长条件下, QMZI四个通道的有效透过率随风速变化呈现出相似的变化趋势. 当风速为正时, QMZI四个通道的有效透过率随风速绝对值增大分别呈现出增大、增大、减小、减小的趋势, 而风速为负时, QMZI四个通道的有效透过率随风速绝对值增大分别呈现出增大、减小、减小、增大的趋势. 此外, 激光中心波长越大, QMZI四个通道有效透过率的变化程度越大. 激光光学谐振腔长越大, 纵模个数越多, QMZI四个通道有效透过率的变化程度越小.

      图8的三个子图可以发现, 激光中心波长较小时, 系统的测风范围也较小. 由(16)式可以看出QMZI响应函数具有周期性且其周期受到多普勒频移影响, 即QMZI只能鉴别一定范围内的多普勒频移. 而由(18)式可以看出, 激光中心频率越小(即激光中心波长越大)的条件下, 相同多普勒频移对应的测风范围越大. 但是在图7的三个子图中, 三种激光中心波长条件下, 系统的测风范围未显示出差异. 这是由于多纵模激光器模式下的直接多普勒测风技术需要多纵模激光纵模间隔与QMZI自由光谱范围相匹配, 而QMZI自由光谱范围随其自身光程差增大而减小. 由(16)式可知QMZI响应函数周期受到光程差的影响, 即QMZI鉴频范围随其光程差增大而减小, 随其自由光谱范围增大而增大. 图7中激光光学谐振腔长较小, 即纵模个数较小, 其对应的纵模间隔和相匹配的QMZI自由光谱范围均较大, 因此三种激光中心波长条件下系统测风范围均超出了±50 m/s.

      图9图10可以看出, 不同激光中心波长、不同激光光学谐振腔长条件下, 系统测风误差基本都呈现出随风速绝对值增大而波动增大的趋势. 激光光学谐振腔长为30 mm(对应纵模个数为7), 激光中心波长分别为1064, 532, 355 nm, 风速在–50—50 m/s之间变化时, 系统测风误差分别可以控制在1.0, 2.0, 4.0 m/s. 激光光学谐振腔长为300 mm(对应纵模个数为61), 激光中心波长分别为1064, 532, 355 nm, 风速在±50, ±30, ±20 m/s之间变化时, 系统测风误差均可以控制在4.0 m/s. 激光光学谐振腔长为300 mm(对应纵模个数为61), 风速在±20 m/s之间变化时, 不同激光中心波长时系统测风误差分别可以控制在2.5, 3.0, 4.0 m/s. 由此对比可以发现, 当激光中心波长越大时, 系统测风误差越小; 当激光光学谐振腔长越大时, 纵模个数越多, 系统测风误差越大.

    • 多纵模激光器体制下直接多普勒测风激光雷达利用QMZI四通道透过率曲线周期相同、相位各相差π/2的特性, 通过探测QMZI四通道输出信号的能量变化来实现对多纵模大气弹性散射回波光谱多普勒频移的高精度鉴别. 大气风场会给多纵模激光中每个纵模激励的大气弹性散射回波光谱都产生多普勒频移, 为了综合利用所有纵模多普勒频移带来的能量变化, MLM DDWL系统参数需要满足“多纵模激光纵模间隔与QMZI透过率曲线自由光谱范围相匹配”的条件. 此外, 多纵模大气弹性散射回波光谱纵模间隔较小时, 相邻纵模之间会相互交叠, 因此导致多普勒频移带来的能量变化变小, 从而导致系统对多普勒频移的鉴别能力下降. 因此MLM DDWL系统参数还需要满足“多纵模激光纵模间隔大于单纵模气溶胶米散射或大气分子瑞利散射回波光谱谱宽”的条件. 多纵模激光的纵模间隔与激光光学谐振腔长、中心波长、辐射线宽、大气温度均相关, 单纵模气溶胶米散射回波光谱的谱宽的激光单纵模光谱谱宽相关, 单纵模大气分子瑞利散射回波光谱的谱宽与激光中心波长、大气温度相关, 而QMZI透过率曲线的自由光谱范围受QMZI光程差的影响. 因此, 要实现MLM DDWL对多纵模大气弹性散射回波光谱多普勒频移的高精度鉴别, 则有必要从气溶胶米散射和大气分子瑞利散射两种条件下去分别考虑实现方案和参数设计. 此外, 多纵模激光器的光谱分布会受到激光器内在因素(如谐振腔的几何形状、材料、尺寸、电磁场分布等因素)和外部环境因素(如环境温度、空气湍流、机械振动等因素)的影响, 这不仅会导致多纵模激光的中心频率发生漂移, 而且会导致多纵模激光不同纵模之间的纵模间隔并非严格一致. 而QMZI的透过率曲线会也受到其自身光学元件瑕疵和外部环境因素的影响, 进而导致多纵模激光纵模间隔与QMZI透过率曲线自由光谱范围的匹配出现误差, 因此需要对MLM DDWL开展误差分析与不确定度评估以解决此问题. 今后将从这些方面继续开展多纵模激光器体制下直接多普勒测风技术的研究工作.

    • 本文在探讨自由运转的Nd: YAG脉冲激光器输出的多纵模激光激励的大气弹性散射回波光谱分布模式的基础上, 借鉴单纵模光谱模式下的直接边缘检测技术, 提出了多纵模激光器体制下直接多普勒测风激光雷达系统, 其利用透过率函数周期相同、相位各相差π/2的可调谐QMZI作为光谱鉴频器. 基于准单色光干涉的部分相干理论与光波偏振效应, 建立了QMZI对多纵模大气弹性散射回波信号有效透过率的数学模型, 并提出了QMZI鉴别多纵模大气弹性散射回波信号多普勒频移的数据反演算法.

      为了验证多纵模激光模式下直接多普勒测风激光雷达技术的可行性及其探测精度, 本文基于Zemax光学仿真软件搭建了QMZI仿真模型, 并在不同纵模个数和不同激光中心波长条件下进行了测风仿真实验. 仿真实验结果表明, 在不同激光中心波长、不同激光光学谐振腔长条件下, 系统的测风误差均随所测风速绝对值的增加呈波动增大的趋势. 激光光学谐振腔长为30 mm(对应纵模个数为7), 激光中心波长分别为1064, 532和355 nm, 风速在–50—50 m/s范围内变化时, 系统的测风误差可分别控制在1.0, 2.0和4.0 m/s以内. 激光光学谐振腔长为300 mm(对应纵模个数为61), 激光中心波长分别为1064, 532和355 nm, 风速在–20—20 m/s范围内变化时, 系统的测风误差分别控制在2.5, 3.0和4.0 m/s以内. 即激光中心波长越大, 激光光学谐振腔长越短, MLM DDWL系统测风精度越高, 测风范围越大. 因此, 多纵模激光模式下直接多普勒测风激光雷达可以实现对大气风场的高精度探测.

    参考文献 (28)

目录

/

返回文章
返回