留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

百kHz重复频率卫星激光测距

程少宇 龙明亮 张海峰 吴志波 秦思 张忠萍

程少宇, 龙明亮, 张海峰, 吴志波, 秦思, 张忠萍. 百kHz重复频率卫星激光测距[J]. 红外与激光工程, 2022, 51(11): 20220121. doi: 10.3788/IRLA20220121
引用本文: 程少宇, 龙明亮, 张海峰, 吴志波, 秦思, 张忠萍. 百kHz重复频率卫星激光测距[J]. 红外与激光工程, 2022, 51(11): 20220121. doi: 10.3788/IRLA20220121
Cheng Shaoyu, Long Mingliang, Zhang Haifeng, Wu Zhibo, Qin Si, Zhang Zhongping. Research on satellite laser ranging at pulse repetition frequency of 100 kHz[J]. Infrared and Laser Engineering, 2022, 51(11): 20220121. doi: 10.3788/IRLA20220121
Citation: Cheng Shaoyu, Long Mingliang, Zhang Haifeng, Wu Zhibo, Qin Si, Zhang Zhongping. Research on satellite laser ranging at pulse repetition frequency of 100 kHz[J]. Infrared and Laser Engineering, 2022, 51(11): 20220121. doi: 10.3788/IRLA20220121

百kHz重复频率卫星激光测距

doi: 10.3788/IRLA20220121
基金项目: 国家自然科学基金(12003056);上海市自然科学基金面上项目(20ZR1467500);上海市仪器共享配套操作与应用技术研究项目(21142201900)
详细信息
    作者简介:

    程少宇,男,硕士生,主要从事卫星激光测距数据处理方面的研究

  • 中图分类号: P228.5

Research on satellite laser ranging at pulse repetition frequency of 100 kHz

  • 摘要: 对卫星激光测距(Satellite Laser Ranging, SLR)回波数与重复频率、脉冲能量及功率关系进行分析,表明单位时间内相同激光回波数,重复频率越高所需激光脉冲能量和平均功率越低;同时对SLR单次测量精度及标准点数据精度进行分析,表明标准点时长内测距点数越多,SLR标准点精度越高。提出点火脉冲群与门控脉冲群收发交替的工作模式,解决超高重复频率后向散射光噪声对激光回波干扰问题。开发多缓冲区存储模式,使测量软件数据实时处理与储存效率提升4~6倍。基于中国科学院上海天文台60 cm口径SLR系统,以快速事件计时器、脉冲群生成器、低噪声单光子探测器等,采用脉冲间隔5 μs、单脉冲能量80 μJ的皮秒激光,收发交替脉冲群模式下实现100 kHz重复频率低轨至高轨卫星的SLR测量,近地星Hy2b标准点精度达到28.55 μm,远地星Galileo218标准点精度达到136.51 μm,为发展更高重频和高精度空间目标激光测距提供了有效方法。
  • 图  1  同样测量回波点数下(D=100)的关系:(a) 单脉冲能量与重复频率;(b) 激光平均功率与重复频率

    Figure  1.  Relationship on the same laser echo (D=100): (a) Between single pulse energy and repetition rate; (b) Between laser average power and repetition rate

    图  2  激光脉冲收发交替工作模式时序模式

    Figure  2.  Time sequence of alternated working mode of laser emission and reception

    图  3  激光与单光子探测器的脉冲群外触发信号时序

    Figure  3.  External trigger pulse trains signal timing for laser and single photon detector

    图  4  固定低频脉冲群周期下的激光损失:(a) 固定低频周期小于实际周期;(b) 固定低频周期大于实际周期

    Figure  4.  Laser loss under the fixed pulse trains: (a) The fixed period is shorter than the actual period; (b) The fixed period is longer than the actual period

    图  5  卫星激光测距系统数据储存模式: (a) 无缓冲区同步模式;(b) 单缓冲区同步模式;(c) 多缓冲区异步模式

    Figure  5.  Data storage mode of SLR system: (a) No buffer in synchronous mode; (b) Single buffer in synchronous mode; (c) Multi-buffer storage in asynchronous mode

    图  6  100 kHz重复频率SLR系统框图

    Figure  6.  Block diagram of SLR system with 100 kHz PRF

    图  7  脉冲群周期信号

    Figure  7.  Pulse trains periodic signal

    图  8  Ajisai卫星的100 kHz SLR激光回波探测实时显示(a)与输出处理(b)

    Figure  8.  Real-time display (a) and data processing (b) for 100 kHz SLR laser echoes from the Ajisai satellite

    图  9  Glonass105卫星100 kHz重复频激光数据处理结果

    Figure  9.  Data processing result for 100 kHz SLR for Glonass105 satellite

    表  1  100 kHz重复率激光测距系统观测数据结果

    Table  1.   Data result from 100 kHz satellite laser ranging system

    Sate nameNumber of laser echoesRMS/cmNumber of normal pointNormal point RMS/μm
    Starlette3588120.97432.39
    Ajisai2034651.59249.85
    Beaconc3937091.85672.22
    Lageos2640212176.75
    Stella2254661.1123.17
    Glonass10597362.031205.73
    Galileo20384011.271138.56
    Galileo21047471.11159.66
    Galileo21851540.981136.51
    Glonass13367062.011245.45
    Hy2b3630510.86428.55
    下载: 导出CSV
  • [1] Lucchesi D, Anselmo L, Bassan M, et al. General relarivity measurements in the field of earth with laser-ranged state of the art and perspectives [J]. Universe, 2019, 5(6): 141. doi:  10.3390/universe5060141
    [2] Zhang Zhongping, Cheng Zhien, Zhang Haifeng, et al. Global laser ranging observation of beidou satellites and data application [J]. Chinese Journal of Lasers, 2017, 44(4): 0404004. (in Chinese) doi:  10.3788/CJL201744.0404004
    [3] Shao Kai, Gu Defeng, Ju Bing, et al. Analysis of Tiangong2 orbit determination and prediction using onboard dualfrequency GNSS data [J]. GPS Solut, 2020, 24(1): 11.
    [4] Deng Huarong, Zhang Haifeng, Long Mingliang, et al. Research on 4 kHz repetition rate satellite laser ranging system and its application [J]. Acta Optica Sinica, 2019, 39(3): 0314002. (in Chinese) doi:  10.3788/AOS201939.0314002
    [5] Sung Ki-Pyoung, Choi Eun-Jung, Lim Hyung-Chul, et al. Development of operation software for high repetition rate satellite laser ranging [J]. Journal of the Korean Society for Aeronautical and Space Sciences, 2016, 44(12): 1103-1111. doi:  10.5139/JKSAS.2016.44.12.1103
    [6] Wu Zhibo, Zhang Zhongping, Zhang Haifeng, et al. The preliminary results of SLR with 10 kHz laser system at Shanghai Station [C]//International Workshop on Laser Ranging, 2013.
    [7] Gen Rengfang, Wu Zhibo, Tang Kai, et al. Precise control of laser emission epoch in laser time transfer with high repetition rate [J]. Infrared and Laser Engineering, 2021, 50(10): 20200473. (in Chinese) doi:  10.3788/IRLA20200473
    [8] Long Mingliang, Zhang Haifeng, Meng Linlin, et al. Satellite laser ranging at high-repetition 10 kHz in all day [J]. Journal of Infrared and Millimeter Waves, 2020, 39(6): 778-785. (in Chinese)
    [9] Bai Zhenxu, Chen Hui, Gao Xiaoqiang, et al. Highly compact nanosecond laser for space debris tracking [J]. Optical Materials, 2019, 98: 109470. doi:  10.1016/j.optmat.2019.109470
    [10] Deng Huarong, Long Mingliang, Zhang Haifeng, et al. Experiment of satellite laser ranging in daytime based on 1064 nm wavelength [J]. Infrared and Laser Engineering, 2020, 49(10): 20200021. (in Chinese) doi:  10.3788/IRLA20200021
    [11] Li Zhulian, Zhai Dongsheng, Zhang Haitao, et al. Superconductivity detector applied to daytime satellite laser ranging experiment and research [J]. Infrared and Laser Engineering, 2020, 49(8): 20190536. (in Chinese) doi:  10.3788/IRLA20190536
    [12] Zhang Haitao, Li Zhulian, Tang Rufeng, et al. Application of array detection technology in laser ranging [J]. Infrared and Laser Engineering, 2020, 49(10): 20200006. (in Chinese)
    [13] Hampf D, Schafer E, Sproll F, et al. Satellite laser ranging at 100 kHz pulse repetition rate [J]. Ceas Space Journal, 2019, 11(4): 363-370. doi:  10.1007/s12567-019-00247-x
    [14] Hampf D, Wagner P, Schafer E, et al. Concept for a new minimal SLR system [C]//21st International Laser Ranging Workshop, 2018.
    [15] Dequal D, Agnesi C, Sarrocco D, et al. 100 kHz satellite laser ranging demonstration at matera laser ranging observatory [C]//ILRS Technical Workshop, 2019.
    [16] Courde C, Mariey H, Chabé J, et al. High repetition rate SLR at GRSM [C]//ILRS Technical Workshop, 2019.
    [17] Wang P, Franz K, Georg K, et al. Contributions to sub-MHz SLR in Graz [C]//ILRS Technical Workshop, 2019.
    [18] Wang P, Steindorfer M A, Koidl F, et al. Megahertz repetition rate satellite laser ranging demonstration at Graz observatory [J]. Optics Letters, 2021, 46(5): 937-940. doi:  10.1364/OL.418135
    [19] Degnan J J. Satellite laser ranging: Current status and future prospects [J]. IEEE Trans Geosci Remote Sens, 1985, 23: 398-413.
    [20] Zhang Zhongping, Zhang Haifeng, Long Mingliang, et al. High precision space debris laser ranging with 4.2 W double-pulse picosecond laser at 1 kHz in 532 nm [J]. Optik, 2019, 179: 691-699. doi:  10.1016/j.ijleo.2018.10.219
    [21] Tyahla L J. Normal point data [EB/OL]. (2018-04-02) [2021-10-28]. https://ilrs.gsfc.nasa.gov/data_and_products/data/ npt/index.html.
    [22] Wu Chenyi, Wang Linli, Shi Haotian, et al. Single-photon ranging with hundred-micron accuracy [J]. Acta Physica Sinica, 2021, 70(17): 174201. (in Chinese)
    [23] Wu Zhibo, Deng Huarong, Zhang Haifeng, et al. Analysis and improvement on stability of satellite laser ranging system [J]. J Infrared Millim Waves, 2019, 38(4): 479-484, 492. (in Chinese)
  • [1] 林海声, 吴志波, 郑敏, 龙明亮, 耿仁方, 喻荣宗, 张忠萍.  卫星激光测距系统皮秒准确度时延标定研究及应用 . 红外与激光工程, 2023, 52(10): 20230070-1-20230070-12. doi: 10.3788/IRLA20230070
    [2] 安宁, 关博文, 张旖伦, 高健, 温冠宇, 董雪, 马磊, 范存波.  卫星激光测距数据处理方法研究进展 . 红外与激光工程, 2021, 50(8): 20200408-1-20200408-9. doi: 10.3788/IRLA20200408
    [3] 耿仁方, 吴志波, 汤凯, 李岩, 孟文东, 张忠萍.  高重复率激光时间比对激光发射时序精确控制 . 红外与激光工程, 2021, 50(10): 20200473-1-20200473-9. doi: 10.3788/IRLA20200473
    [4] 郑立, 汪会波, 田文龙, 张大成, 韩海年, 朱江峰, 魏志义.  LD泵浦的高重复频率全固态飞秒激光器(特邀) . 红外与激光工程, 2020, 49(12): 20201069-1-20201069-13. doi: 10.3788/IRLA20201069
    [5] 邓华荣, 龙明亮, 张海峰, 吴志波, 汤凯, 张忠萍.  1 064 nm波长白天卫星激光测距试验 . 红外与激光工程, 2020, 49(10): 20200021-1-20200021-6. doi: 10.3788/IRLA20200021
    [6] 李翔, 汪宏, 乔忠良, 张宇, 牛智川, 佟存柱, 刘重阳.  2 μm GaSb基被动锁模激光器重复频率变化的研究(特邀) . 红外与激光工程, 2020, 49(12): 20201054-1-20201054-5. doi: 10.3788/IRLA20201054
    [7] 孙敬华, 孙克雄, 林志芳, 孙继芬, 晋路, 徐永钊.  高功率高重复频率飞秒掺镱光纤激光频率梳的研究(特邀) . 红外与激光工程, 2019, 48(1): 103001-0103001(9). doi: 10.3788/IRLA201948.0103001
    [8] 吕华昌, 陈念江, 钟声远, 李楠楠, 李长桢, 郭丽娜, 吴健, 耿园园.  试验卫星激光反射器的设计和试验 . 红外与激光工程, 2018, 47(8): 806005-0806005(9). doi: 10.3788/IRLA201847.0806005
    [9] 温冠宇, 王爽, 安宁, 董雪, 韩兴伟.  光行差对高轨卫星激光测距的影响分析 . 红外与激光工程, 2018, 47(9): 906001-0906001(5). doi: 10.3788/IRLA201847.0906001
    [10] 邓华荣, 吴志波, 李朴, 汤凯, 张海峰.  偏振调能技术对SLR数据偏差的改善研究 . 红外与激光工程, 2017, 46(9): 917005-0917005(4). doi: 10.3788/IRLA201746.0917005
    [11] 吴志波, 邓华荣, 张海峰, 汤凯, 张忠萍.  高重复率卫星激光测距中后向散射干扰及规避 . 红外与激光工程, 2017, 46(2): 206002-0206002(6). doi: 10.3788/IRLA201746.0206002
    [12] 黄珂, 黄超, 赵柳, 马连英, 唐影, 朱峰, 易爱平, 钱航, 李高鹏, 冯国斌, 赵军.  流场均匀性对重频HF激光能量稳定输出的影响 . 红外与激光工程, 2016, 45(1): 106007-0106007(6). doi: 10.3788/IRLA201645.0106007
    [13] 吴志波, 邓华荣, 张海峰, 汤凯, 张忠萍.  卫星激光测距中光束亮度的偏振影响及应用 . 红外与激光工程, 2016, 45(3): 306005-0306005(5). doi: 10.3788/IRLA201645.0306005
    [14] 武伟, 陈桂明, 赵娜, 樊博璇.  激光在高速钢表面加工沟槽表面织构的实验研究 . 红外与激光工程, 2016, 45(2): 206008-0206008(8). doi: 10.3788/IRLA201645.0206008
    [15] 李语强, 李荣旺, 李祝莲, 翟东升, 伏红林, 熊耀恒.  空间碎片激光测距应用研究 . 红外与激光工程, 2015, 44(11): 3324-3329.
    [16] 何洋, 陆君, 陈飞, 李殿军, 郑长彬, 谢冀江, 郭劲.  声光调Q CO2激光器的实验研究 . 红外与激光工程, 2015, 44(8): 2280-2285.
    [17] 代福, 熊胜明.  高重复频率DPL 激光对光学薄膜元件损伤实验研究 . 红外与激光工程, 2014, 43(7): 2074-2080.
    [18] 郑义军, 谭荣清, 王东雷, 张阔海, 黄文武, 刘世明, 李能文, 孙科, 卢远添, 刁伟伦.  新型高重复频率脉冲CO2 激光器 . 红外与激光工程, 2013, 42(10): 2697-2701.
    [19] 陈豪, 石磊, 马丽华, 彭志明, 徐志燕, 乔忆炜.  脉冲数目和重复频率对纳秒激光推进的影响 . 红外与激光工程, 2013, 42(9): 2403-2408.
    [20] 汤伟, 郭劲, 邵俊峰, 郭汝海, 王挺峰.  激光重频对脉冲非稳腔TEA CO2激光远场传输特性的影响分析 . 红外与激光工程, 2013, 42(9): 2380-2385.
  • 加载中
图(9) / 表(1)
计量
  • 文章访问数:  193
  • HTML全文浏览量:  20
  • PDF下载量:  59
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-02-24
  • 修回日期:  2022-03-23
  • 录用日期:  2022-05-16
  • 刊出日期:  2022-11-30

百kHz重复频率卫星激光测距

doi: 10.3788/IRLA20220121
    作者简介:

    程少宇,男,硕士生,主要从事卫星激光测距数据处理方面的研究

基金项目:  国家自然科学基金(12003056);上海市自然科学基金面上项目(20ZR1467500);上海市仪器共享配套操作与应用技术研究项目(21142201900)
  • 中图分类号: P228.5

摘要: 对卫星激光测距(Satellite Laser Ranging, SLR)回波数与重复频率、脉冲能量及功率关系进行分析,表明单位时间内相同激光回波数,重复频率越高所需激光脉冲能量和平均功率越低;同时对SLR单次测量精度及标准点数据精度进行分析,表明标准点时长内测距点数越多,SLR标准点精度越高。提出点火脉冲群与门控脉冲群收发交替的工作模式,解决超高重复频率后向散射光噪声对激光回波干扰问题。开发多缓冲区存储模式,使测量软件数据实时处理与储存效率提升4~6倍。基于中国科学院上海天文台60 cm口径SLR系统,以快速事件计时器、脉冲群生成器、低噪声单光子探测器等,采用脉冲间隔5 μs、单脉冲能量80 μJ的皮秒激光,收发交替脉冲群模式下实现100 kHz重复频率低轨至高轨卫星的SLR测量,近地星Hy2b标准点精度达到28.55 μm,远地星Galileo218标准点精度达到136.51 μm,为发展更高重频和高精度空间目标激光测距提供了有效方法。

English Abstract

    • 卫星激光测距(Satellite Laser Ranging,SLR)通过精确测定激光脉冲从地面观测点到卫星反射器往返时间间隔,从而获取卫星与站点之间的高精确距离[1]。北斗导航卫星激光测距定轨残差均方根值(Root Mean Square, RMS)好于1 cm[2],天宫二号飞行器激光数据轨道检核精度达到1.7 cm[3]。SLR测距频率由早期使用低重复频率(小于20 Hz)已发展到目前普遍使用的千赫兹[4-9],由于高重频激光测距具有目标搜索快、激光回波数据量大的特点,能够增大测量标准点时长内的回波数,获得更高精度激光标准点数据,在国外得到快发展,并向更高的重复频率激光测距发展。2018年,中国科学院上海天文台激光测距站在4 kHz重复频率激光测距下分析测距数据量与标准点精度,表明4 kHz比1 kHz重复频率的测量数据量和标准点数据精度分别提高了约2.62倍和1.62倍,4 kHz的SLR对低轨卫星标准点精度平均为135 μm,对高轨卫星标准点精度平均为405 μm [4]。次年,通过升级激光及距离门控,使用激光点火信号移动的方式进行后向散射激光规避,实现全天时从低轨到同步轨道的10 kHz的激光测距测量[6-9]。然而,测距数据中依旧存在明显的后向散射激光对激光回波数据影响情况[8]。随着激光技术及激光探测技术的不断发展,低噪声的红外1064 nm得到应用,2020年,中国科学院上海天文台通过红外相机监视光束指向实现了白天1064 nm的SLR[10];中国科学院云南天文台使用超导纳米线单光子探测器实现了白天SLR[11],并于次年通过阵列探测技术提高了弱回波信号的探测效率[12]

      目前,国际上超高重频(100 kHz及以上)激光测距主要采用脉冲调频burst脉冲群及自由模式运转的单光子探测器来实现。2018年,德国Daniel Hampf等人采用体积小的光纤纳秒激光器实现了激光发射脉冲数100 k对20000 km高轨SLR,测得标准点精度达5~15 mm[13],接着完成了小尺寸超高重复频率100 kHz激光测距测系统,其尺寸仅为1.8 m ×1.2 m×1.6 m,质量200 kg[14]。2019年,在德国斯图加特(Stuttgart)国际激光测距网组织的全球技术大会上,意大利D. Dequal等人分析了100 kHz激光测距标准点精度有望突破亚毫米量级[15];法国C. Courde等人报告了100 kHz 激光测距测量[16];奥地利Graz的SLR站报告了500 kHz的SLR[17],并于2020年基于收发交替频率实现了1 MHz重复频率白天SLR[18]。然而,对于超高重复频率激光测距,激光脉冲数目多,测量数据需要高速采集和储存,且激光脉冲之间脉冲间距短,激光后向散射光极易对卫星激光回波信号干扰,增大了激光回波信号识别难度。

      文中从激光测距雷达方程对平均回波数进行分析,给出了超高重复频率技术优势,提出了单光子探测器门控模式下收发交替的脉冲群方式,分析了收发交替周期随卫星距离变化的关系,解决了大气后向散射噪声光对激光回波干扰。研究了不同数据储存模式的时效性,以保证数据采集处理的实时性,最后以中国科学院上海天文台SLR系统为平台,采用工作在Geiger模式的低噪声外触发单光子探测器,实现了100 kHz重复频率532 nm波长SLR。

    • 激光测距雷达方程是评估激光测距系统测量能力的基础,其公式如下[19]

      $$ {n_0} = \frac{{16}}{{{\pi ^2}}} \cdot \frac{{\lambda {\eta _q}}}{{hc}} \cdot \frac{{{E_t}{A_r}{A_s}}}{{{R^4}\theta _t^2\theta _s^2}} \cdot {T^2} \cdot {K_t} \cdot {K_r} \cdot \alpha $$ (1)

      式中:$ {n_0} $为系统接收平均激光回波光电子数;$ \lambda $为激光波长;$ h $为普朗克常量;$ c $为光速;$ {n_q} $为探测器探测效率;$ {E_t} $为激光器单脉冲能量;$ {A_r} $为系统有效接收面积;$ {A_s} $为目标反射器反射面积;$ {T_2} $为双程大气透过率;$ {K_t} $为发射光学系统效率;$ {K_r} $为接收光学系统效率;$ {\theta _t} $为目标反射器发散角;$ {\theta _s} $为激光束发散角;$ R $为目标卫星到测距站点的距离;$ \alpha $为大气衰减因子。激光测距一般采用单光子探测器,其光电转换产生的光电子数服从泊松分布,由此产生一个以上光电子的概率为:

      $$ P(0,N) = \sum\limits_1^\infty {p(m,N)} = 1 - P(0,N) = 1 - {{\rm{e}}^{ - N}} $$ (2)

      激光测距过程中不考虑后散信号的光子数,当单光子探测器探测到为背景噪声与暗噪声光电子时,卫星激光回波的光电子将得不到探测,故单个激光脉冲的卫星激光回波探测概率为:

      $$ P(0,{n_0}) = {{\rm{e}}^{ - ({n_1} + {n_2})}}(1 - {{\rm{e}}^{ - {n_0}}}) $$ (3)

      式中:$ {n_1} $为背景光噪声产生的光电子数;$ {n_2} $为单光子探测器自身暗噪声光电子数。由上式可知,为了提高卫星激光回波的探测概率,应降低$ {n_1} $$ {n_2} $。在重复频率$ f $下,结合公式(1),单位时间内的回波平均点数$ D $为:

      $$ D{\text{ = }}f{{\rm{e}}^{ - ({n_1} + {n_2})}}\left(1 - {{\rm{e}}^{ - \tfrac{{16}}{{{\pi ^2}}} \cdot \tfrac{{\lambda {\eta _q}}}{{hc}} \cdot \tfrac{{{E_t}{A_r}{A_s}}}{{{R^4}\theta _t^2\theta _s^2}} \cdot {T^2} \cdot {K_t} \cdot {K_r} \cdot \alpha }}\right) $$ (4)

      由公式(4)得,为了增大单位时间内的回波平均点数D,可以增加激光的重复频率,降低背景光噪声$ {n_1} $与暗噪声的光电子数$ {n_2} $、系统有效接收面积$ {A_r} $、激光器单脉冲能量$ {E_t} $。由此,在$ D $恒定的情况下,$ f $可以降低激光器单脉冲能量及系统有效接收面积,即可用口径小的接收望远镜以及高重复频率、低脉冲能量的激光器系统,降低系统成本,利于系统小型化发展。在中国科学院上海天文台千赫兹激光测距系统测量中[20]λ=532 nm,h=6.6260693×10−34 J·s,c= 3×108 m/s,nq=20%,T=0.6,Ar=0.251 m2Kt=0.6,Kr=0.4,α=0.2,θt=30 μrad,As=0.1 m2R=20000 km,θs=30 μrad,在不考虑噪声的情况下,D=100点内,Etf的关系为:

      $$ {E_t}{{ = }}{{- 0.000\;207}}{\rm{ln}}\left(1 - \frac{{{\text{100}}}}{f}\right) $$ (5)

      根据公式(5)得到的关系曲线如图1(a)所示,可以看出,重复频率越高,在单位时间内获得相同回波数时所需要单脉冲能量越低。由公式(5)可得激光输出平均功率随重复频率的变化为:

      $$ P = f{E_t}{{ = - 0}}{{.000\;207}}f{\rm{ln}}\left(1 - \frac{{{\text{100}}}}{f}\right) $$ (6)

      根据公式(6)得到的关系曲线如图1(b)所示,可以看出,重复频率越高,在单位时间内获得相同回波数时所需要激光输出的平均功率较低。

      图  1  同样测量回波点数下(D=100)的关系:(a) 单脉冲能量与重复频率;(b) 激光平均功率与重复频率

      Figure 1.  Relationship on the same laser echo (D=100): (a) Between single pulse energy and repetition rate; (b) Between laser average power and repetition rate

    • 由参考文献[13],对SLR系统的单次测量精度$ \sigma _{{\text{single}}}^{} $进行估计:

      $$ {\text{ }} {\sigma _{{\text{single}}}} = \sqrt {\sigma _L^2 + \sigma _{D1}^2 + \sigma _{D2}^2 + \sigma _{ET}^2 + \sigma _S^2} $$ (7)

      式中:$ \sigma _L^{} $为由激光脉冲引起的时间偏差(与脉冲持续时间有关);$\sigma _{{{D}}1}^{}$为探测器开门时间的偏差,$\sigma _{{{D}}2}^{}$为探测器关门时间的偏差;$\sigma _{{{ET}}}^{}$为事件计时器响应的偏差;$\sigma _{{S}}^{}$为测距目标脉冲响应偏差。目前国际激光测距规范采用标准点反映数据精度[21],标准点精度$ \sigma $计算公式如下[4,13,22]

      $$ \sigma = \frac{{{\sigma _{{\text{single}}}}}}{{\sqrt N }} $$ (8)

      式中:N为激光标准点时长内数据点数。由公式(7)可知,激光测量的工作频率基本不影响$ \sigma _{{\text{single}}}^{} $,由此,增加激光测距的工作频率可增加卫星激光标准点时长内的回波数(低轨卫星标准点时长为30 s,高轨卫星标准点时长为300 s)[21],从而可获得更高精度的标准点精度。

    • 超高重复频率下,激光脉冲之间的间距短,激光在大气中的后向散射光比较强,此时若探测器打开,一方面很容易导致探测器饱和而损坏,另一方面后向散射激光会干扰卫星激光回波,增大卫星激光回波的识别难度,以收发交替脉冲群模式进行激光后向散射规避。相比常规1 kHz重复频率SLR测距数据,超高重复频率(100 kHz)的SLR测距数据量将提高两个数量级,常规数据储存模式将造成计算机操作系统产生大量阻塞而崩溃,采用多缓冲区异步数据储存模式有效避免阻塞,提高大数据量的处理效率。

    • 采用收发交替脉冲群模式进行后向散射光规避,其工作模式如图2所示,激光器输出激光脉冲,发送一定脉冲数后停止发送。由于激光已经停止输出,待其在大气散射光至望远镜的光子返回后,开启单光子探测器收集远距离卫星激光回波,避免后向散射光的干扰;单光子探测器接收对应的回波后,停止工作,此时激光发射,以此循环往复。

      图  2  激光脉冲收发交替工作模式时序模式

      Figure 2.  Time sequence of alternated working mode of laser emission and reception

      为了确定脉冲群周期,并考虑激光脉冲输出与单光子探测器工作模式,设计的触发时序如图3所示。

      图  3  激光与单光子探测器的脉冲群外触发信号时序

      Figure 3.  External trigger pulse trains signal timing for laser and single photon detector

      图3中,激光脉冲群外触发信号最大持续工作时间为:

      $$ \Delta {t_1} = \frac{{2L}}{c} $$ (9)

      式中:c为光的传播速度;L为卫星距SLR站的距离,测距过程中L随测量卫星的方位、俯仰而变。假定T为激光器发射脉冲间距,故在$ \Delta {t_1} $内激光发射的脉冲数最大为:

      $$ {{{N}}_{{L}}} = \frac{{2L}}{{c{{T}}}} $$ (10)

      此时$ \Delta {t_2} $内无外触发信号激光脉冲未发射,而单光子探测器由激光脉冲群对应的门控脉冲群信号开启,接收激光回波信号,则对应单光子探测器开启最大持续工作时间$ \Delta {t_2} = \Delta {t_1} $,此时激光输出脉冲群周期TL为:

      $$ {{{T}}_{{L}}} = \Delta {t_1} + \Delta {t_2} = \frac{{4L}}{c} $$ (11)

      $ \Delta {t_1} $内激光发射脉冲数的最大占空比DR为:

      $$ DR = \frac{{\Delta {t_1}}}{{{T_L}}} = 50 {\text{%}} $$ (12)

      即发射激光的最大允许脉冲数为激光工作重复频率的一半。对不同轨道距离的卫星,由公式(10)、(11)可得对应的脉冲群周期,例如对某一时刻距SLR站距离为3750 km近地卫星,脉冲群周期TL为50 ms,对应脉冲群收发交替频率为20 Hz,相应脉冲群内的脉冲发射最大时间为25 ms。而对于37500 km的同步卫星,脉冲群周期为TL为0.5 s,对应收发交替频率为2 Hz,对应脉冲发射最大时间为TL/2,即0.25 s。实际测距过程中,卫星距离随着卫星轨道与SLR站的夹角变化,TL也对应变化。对不同卫星轨道预报距离,可以根据公式(10)、(11)计算脉冲群的重复频率及脉冲群内的脉冲数目,文中结合已有SLR系统中的固定低频门控与脉冲群生成器的方式实现脉冲群的收发交替模式,根据目标预报的特定距离确定低频脉冲群周期,实现超高重复频率100 kHz测量。使用低频门控触发脉冲群生成器获得脉冲群,部分激光脉冲将得不到在对应有效的距离门控下的单光子探测器的外触发工作,因此得不到探测而带来损失。图4展示了脉冲群周期的激光脉冲数损失,TR表示固定的脉冲群周期,ΔT表示实际能接收到回波的距离门控下开启单光子探测器的有效接收窗口大小。

      图  4  固定低频脉冲群周期下的激光损失:(a) 固定低频周期小于实际周期;(b) 固定低频周期大于实际周期

      Figure 4.  Laser loss under the fixed pulse trains: (a) The fixed period is shorter than the actual period; (b) The fixed period is longer than the actual period

      图4可知,脉冲群周期固定时,回波损失率Qloss为:

      $$ Q_{\text {loss }}=\left\{\begin{array}{l}1-\dfrac{2 {T}_{{R}}-{T}_{{L}}}{{T}_{{L}}}, \;\;\;\;\; {T}_{{R}}<{T}_{{L}} \\ 1-\dfrac{2 {T}_{{R}}-{T}_{{L}}}{{T}_{{L}}},\;\;\;\;\; {T}_{{R}} \geqslant {T}_{{L}}\end{array}\right.$$ (13)

      以脉冲群周期为50 ms的近地卫星为例,若实际距离为3500 km,由公式(11)可知对应的脉冲群周期TL应为47 ms,由公式(13)可知发射的脉冲群内将有6.38%的脉冲因得不到探测而损失。

    • 常规千赫兹重复率卫星激光测距系统中使用了无缓冲区或单缓冲区同步模式实时采集、储存测距数据,其工作原理如图5(a)和图5(b)所示。无缓冲区模式每次采集测距数据后直接写入磁盘,单缓冲区模式将测距数据暂存在缓冲区,缓冲区满后持久化到磁盘,等待写入完毕后继续采集。无缓冲区或单缓冲区都在写入磁盘时造成采集过程阻塞,降低了数据采集效率。使用中国科学院上海天文台激光测距系统测试[23],缓冲区与单缓冲区采集模式对应最大采集频率,无缓冲区模式在激光发射频率为20 kHz时软件出现阻塞,单缓冲区模式在激光发射频率为80 kHz时软件出现阻塞,不能满足高重复频率100 kHz的SLR,必须采用新的采集、储存方法。

      图  5  卫星激光测距系统数据储存模式: (a) 无缓冲区同步模式;(b) 单缓冲区同步模式;(c) 多缓冲区异步模式

      Figure 5.  Data storage mode of SLR system: (a) No buffer in synchronous mode; (b) Single buffer in synchronous mode; (c) Multi-buffer storage in asynchronous mode

      文中提出采用多缓冲区异步数据储存模式进行数据实时储存,其工作原理如图5(c)所示,具有一个前台缓冲区和多个后台缓冲区,每个缓冲区在不同线程中工作。数据储存过程中,采集线程向前台缓冲区存入测距数据,同时储存线程将后台缓冲区的数据写入到磁盘,若后台缓冲区均为空,则储存线程将休眠。前台缓冲区接收数据达到设定的阈值时,将空的后台缓冲区与前台缓冲区交换,此时原本的前台缓冲区变为后台缓冲区,而原本空的后台缓冲区变为新的前台缓冲区,接收采集线程传输的数据。由于数据写入磁盘在后台储存线程中进行,采集线程不会受限于磁盘IO带来时间损失。使用计算机读取实际测距数据至内存中,进行仿真比较三种模式储存1~200 k的等量数据所消耗时间,测试使用的计算机配置为:Intel i5-9300H八核CPU,8 GB RAM,测得多缓冲区异步数据储存模式储存200 k测距数据消耗时间分别是无缓冲区和单缓冲区的0.1452与0.2238倍,经测试,多缓冲区异步数据储存模式可满足100 kHz重复率激光测量要求。

    • 以中国科学院上海天文台激光测距系统为平台,根据参考文献[23]中激光测距基本框架,100 kHz重复频率激光测距系统框图如图6所示。

      图  6  100 kHz重复频率SLR系统框图

      Figure 6.  Block diagram of SLR system with 100 kHz PRF

      图6中,相对于参考文献[23]中常规千赫兹重复率激光测距系统,脉冲群发生器基于低频SLR的卫星轨道预报值确定给出对应预报的低频SLR激光与单光子探测器开启触发信号,低频SLR激光开启触发信号触发脉冲群生成器1 (Pulse group generator 1,型号为DG645)产生脉冲间距5 μs (200 kHz)的点火脉冲群信号,占空比50%,激光器经点火脉冲群触发,发射脉冲群激光脉冲;同样,低频SLR单光子探测器开启触发信号给脉冲群生成器2 (Pulse group generator 2,型号为DG645)产生脉冲群信号触发,再触发单光子探测器。由此,对应单位时间内总脉冲数为100 k, 此时频率为100 kHz。采用华东师范大学提供的低暗噪声APD单光子探测器,同时对接收到的卫星激光回波进行光谱滤波及空间滤波,降低背景光噪声,提高激光回波信号探测成功率。激光脉冲发射时,单光子探测器关闭;单光子探测器开启时,激光脉冲停止发射,由此实现脉冲群收发交替模式,避免单光子探测器在开启时,激光在大气中后向散射光被单光子探测器探测接收,实现后向散射激光的规避。事件计时器 (产地:拉脱维亚,型号USB,A033)记录激光脉冲的主波与单光子探测器的回波时刻,并传输到计算终端,计算终端进行信号识别后,使用多缓冲区数据储存模式存储数据。

      经过低频SLR触发产生的脉冲群触发信号如图7所示,其中波形1为低频SLR产生的低频SLR周期信号,波形2为触发脉冲群生成器1生成的脉冲激光点火脉冲群信号,触发激光器输出激光脉冲,波形3为触发脉冲群生成器2生成的单光子探测器触发信号,触发单光子探测器工作。波形2与波形3中的脉冲群内脉冲之间的间距为5 μs。采用重复频率200 kHz (脉冲间距5 μs)的激光器,激光发射、回波探测交替工作模式 (100 kHz)下,激光平均输出功率为8 W,单脉冲能量80 μJ,各脉冲之间幅值相等。为了降低回波损失率,根据卫星预报距离的平均值,对近地星(LEO)采用40 Hz低频,群内脉冲数为2500个,远地星(HEO)采用4 Hz,群内脉冲数为25000个,即对应单位时间的脉冲数为100 k。

      图  7  脉冲群周期信号

      Figure 7.  Pulse trains periodic signal

    • 基于中国科学院上海天文台60 cm口径SLR系统,由图6建立100 kHz测距系统,对近地轨道、远地轨道合作目标开展了100 kHz重复率激光测量,数据测量及处理结果如表1所示。

      表 1  100 kHz重复率激光测距系统观测数据结果

      Table 1.  Data result from 100 kHz satellite laser ranging system

      Sate nameNumber of laser echoesRMS/cmNumber of normal pointNormal point RMS/μm
      Starlette3588120.97432.39
      Ajisai2034651.59249.85
      Beaconc3937091.85672.22
      Lageos2640212176.75
      Stella2254661.1123.17
      Glonass10597362.031205.73
      Galileo20384011.271138.56
      Galileo21047471.11159.66
      Galileo21851540.981136.51
      Glonass13367062.011245.45
      Hy2b3630510.86428.55

      图8为Ajisai卫星的脉冲群超高重复频率激光测距激光回波探测实时软件测量显示(图8(a))与测量数据处理(图8(b)),图9为远地星Glonass105卫星的测量数据处理,图8图9的数据显示及处理呈现周期性状态,表明单光子探测器在工作时是启停交替方式,体现脉冲群的收发交替模式。

      图  8  Ajisai卫星的100 kHz SLR激光回波探测实时显示(a)与输出处理(b)

      Figure 8.  Real-time display (a) and data processing (b) for 100 kHz SLR laser echoes from the Ajisai satellite

      图  9  Glonass105卫星100 kHz重复频激光数据处理结果

      Figure 9.  Data processing result for 100 kHz SLR for Glonass105 satellite

      根据参考文献[21]中标准点时长以及表1中测距点数可以推算出每圈观测数据能够生成的标准点个数及每个标准点内测距数据点数,由公式(8)可以计算出每圈观测数据的标准点精度。对近地星的测量,精度最好的结果为Hy2b卫星的28.55 μm;对远地星的测量,Galileo218卫星的精度达到136.51 μm,相比德国Daniel Hampf等人采用光纤纳秒激光器测量标准点精度[13],实现了数量级的提高,体现了皮秒激光在超高重复百千赫兹卫星激光测距上的优势。

    • 从激光测距雷达方程出发,得出单位时间内的回波平均点数与重复频率及回波探测概率关系,表明测距重复频率越高,可降低对激光输出功率、激光脉冲能量、接收望远镜的尺寸的需求。通过重复频率100 kHz的SLR发射时间段与卫星激光回波接收时间段分隔开,解决后向散射光对卫星激光回波的干扰;提出多缓冲区异步数据储存模式,满足100 kHz高重复频率数据采集要求。基于中国科学院上海天文台卫星激光测距系统,设置激光收发交替控制,实现100 kHz重复频率激光测距,近地星标准点测距精度最优结果可达到28.55 μm,远地星最优标准点测距精度达到136.51 μm,验证了脉冲群收发交替对避免后向散射光的有效性。后期进一步通过高速的距离门控编程实现脉冲群可变,避免脉冲损失。通过文中的研究与试验,脉冲群交替模式将成为超高重复频率激光测距实现的有效方法,有力促进了激光测距技术的发展与应用。

      致谢 该文得到大恒新纪元科技股份公司在超高重复频率皮秒激光器与华东师范大学在单光子探测器上的大力支持,在此表示感谢。

参考文献 (23)

目录

    /

    返回文章
    返回