-
激光测距雷达方程是评估激光测距系统测量能力的基础,其公式如下[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 m2,Kt=0.6,Kr=0.4,α=0.2,θt=30 μrad,As=0.1 m2,R=20000 km,θs=30 μrad,在不考虑噪声的情况下,D=100点内,Et与f的关系为:$$ {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)所示,可以看出,重复频率越高,在单位时间内获得相同回波数时所需要激光输出的平均功率较低。
-
由参考文献[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中,相对于参考文献[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。
-
基于中国科学院上海天文台60 cm口径SLR系统,由图6建立100 kHz测距系统,对近地轨道、远地轨道合作目标开展了100 kHz重复率激光测量,数据测量及处理结果如表1所示。
表 1 100 kHz重复率激光测距系统观测数据结果
Table 1. Data result from 100 kHz satellite laser ranging system
Sate name Number of laser echoes RMS/cm Number of normal point Normal point RMS/μm Starlette 358812 0.97 4 32.39 Ajisai 203465 1.59 2 49.85 Beaconc 393709 1.85 6 72.22 Lageos2 6402 1 2 176.75 Stella 225466 1.1 1 23.17 Glonass105 9736 2.03 1 205.73 Galileo203 8401 1.27 1 138.56 Galileo210 4747 1.1 1 159.66 Galileo218 5154 0.98 1 136.51 Glonass133 6706 2.01 1 245.45 Hy2b 363051 0.86 4 28.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],实现了数量级的提高,体现了皮秒激光在超高重复百千赫兹卫星激光测距上的优势。
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,为发展更高重频和高精度空间目标激光测距提供了有效方法。Abstract: The relationship between the number of satellite laser ranging (SLR) echoes and the repetition rate, laser pulse energy and power are analyzed, it shows that under the same laser echoes, the higher the repetition rate, the lower the laser pulse energy and laser average power are required. Meanwhile, single-shot accuracy and normal point (NP) accuracy of SLR are analyzed, it shows that the more measuring points within the NP time are, the higher the NP accuracy would be. The working mode of firing pulse trains and receiving pulse trains triggered by a fixed period range gate is proposed to solve the interference problem of laser echo caused by back-scattered laser noise with ultra-high pulse repetition frequency (PRF). The multi-buffer storage mode is developed to improve the real-time processing and storage efficiency of data in measurement software by 4-6 times. Based on the telescope aperture of 60 cm SLR in Shanghai Astronomical Observatory, CAS, the PRF of 100 kHz SLR is realized by using fast event timing, ultra-high pulse trains generator, low noise single photon detector and picosecond laser with pulse spacing of 5 μs, single pulse energy of 80 μJ. Measured data of satellites by the PRF of 100 kHz SLR is from low earth orbit (LEO) and high earth orbit (HEO). The NP accuracy for Hy2b satellite in LEO is 28.55 μm, and that of Galileo218 satellite in HEO is 136.51 μm. This study provides an effective method for higher frequency and high-accuracy space target laser ranging.
-
表 1 100 kHz重复率激光测距系统观测数据结果
Table 1. Data result from 100 kHz satellite laser ranging system
Sate name Number of laser echoes RMS/cm Number of normal point Normal point RMS/μm Starlette 358812 0.97 4 32.39 Ajisai 203465 1.59 2 49.85 Beaconc 393709 1.85 6 72.22 Lageos2 6402 1 2 176.75 Stella 225466 1.1 1 23.17 Glonass105 9736 2.03 1 205.73 Galileo203 8401 1.27 1 138.56 Galileo210 4747 1.1 1 159.66 Galileo218 5154 0.98 1 136.51 Glonass133 6706 2.01 1 245.45 Hy2b 363051 0.86 4 28.55 -
[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)