摘要
为分析天琴计划激光测距台站卫星激光测距探测能力,设置单脉冲能量(平均功率)分别为0.15 mJ(0.015 W),0.4 mJ(0.04 W)和4 mJ(0.4 W)对同步轨道卫星qzs2、compass i3和compass i5分别于夜晚和白天进行激光测距实验。理论方面,分析了白天天光背景噪声强度,并结合雷达方程计算了不同平均功率条件下卫星激光测距的有效回波率,重点分析了平均功率对卫星激光测距探测能力的影响。实验方面,固定望远镜俯仰角(E=50°),旋转望远镜测量不同指向时的天光背景噪声强度。采用探测效率为60%(@1064 nm)的超导探测器阵列,并结合空间滤波、时间滤波和光谱滤波的方法抑制背景噪声,开展了小功率白天卫星激光测距。分析得到,白天对同步卫星激光测距的最小功率为0.04 W,夜晚对同步卫星激光测距的最小功率为0.015 W。天琴计划激光测距台站已具备白天和夜晚全时段常规卫星激光测距能力,将为今后天琴计划引力波探测卫星全时段激光测距奠定基础。
激光发射重复频率和单脉冲能量是卫星激光测距的重要指标,将直接决定激光器的平均功率。最早的卫星激光测距系统采用10~20 Hz,单脉冲能量30~150 mJ,即平均功率0.3 W的激光器进行卫星激光测距,其缺点是有效回波光子数少。近二十年来为提高激光测距回波率,国内外激光测距台站发展了高重复频率激光测距技术。白天激光测距方面,云南天文台的李祝莲等
卫星激光测距正向着高重频、低能量方向发展。本文以天琴计划激光测距台
卫星激光测距指对装有激光角反射器的卫星进行测量。其有效回波光子数可以由
, | (1) |
式中,λ为发射激光的波长,h和c分别为普朗克常数和光在真空中传播的速度,ET为发射激光单脉冲能量,为探测器探测效率,为发射系统效率,为接收系统效率,为大气透过率,为光束增益,为目标散射截面,R为目标与测站的距离,为望远镜有效接收面积。
将由能量表述的激光雷达方程转换为平均功率形式:
, | (2) |
, | (3) |
式中,为系统接收功率,为系统发射功率。
本文以天琴台站激光测距系统为例,其系统指标参数如
接收望远镜口径 | 1.2 m |
---|---|
发射系统光学效率 | 0.64 |
遮拦比 | 0.25 |
探测器探测效率 | 60%@ 1 064 nm |
接收系统光学效率 | 0.28 |
激光中心波长 | 1 064 nm |
激光发散角(全角) | 2″ |
大气透过率 | 0.5 |

图1 卫星激光测距有效回波光子数(激光重复频率100 Hz)
Fig. 1 The number of effective echo photons from satellite laser ranging (100 Hz)
激光雷达方程计算结果如
通常认为卫星激光测距有效回波光子数与激光单脉冲能量之间具有正相关性质,当单脉冲能量越大时有效回波光子数越多。然而,有效回波光子数本质上是由系统接收的平均功率所决定的,即系统接收到的平均功率越大,有效回波光子数也就越多。比较国内外高重复频率激光测距技术,提高重复频率的同时,发射激光单脉冲能量减小,但依旧保持较高的平均功率。因此,影响激光测距探测能力的是平均功率,而不仅仅是激光器单脉冲能量。
白天对激光测距有效回波率最大的影响因素之一就是背景噪声。因此,仿真分析背景噪声量级至关重要。望远镜可以接收的天光背景强度主要受到系统接收视场、天顶角与太阳夹角影响。天光背景辐射主要是由于太阳光线照射在大气分子、气溶胶分子及其他微粒上发生散射和吸收作用产生。当光束传输途径中的微粒半径大小小于太阳的辐射波长时,这种情况下的散射强度服从瑞利散射,如
, | (4) |
式中,表示观测方向与太阳光线入射方向的夹角。瑞利散射强度随观测方向与太阳夹角的增大而减小。

图2 太阳Rayleigh散射强度
Fig. 2 Sun Rayleigh scattering intensity
当光束传输途径中的微粒半径大小大于太阳的辐射波长时,这种情况下的散射强度服从Mie散射,如
. | (5) |

图3 太阳Mie散射强度
Fig. 3 Sun Mie scattering intensity
在大气粒子散射的基础上再考虑传输空间中粒子的吸收作用,Mie散射强度随观测方向与太阳的夹角增大而减小。
白天激光测距情况下,背景噪声的估计对于卫星激光测距系统回波率分析具有重要意义。以光子数描述天光背景噪声,采用激光雷达方程计算背景噪声,如
, | (6) |
式中,NB表示探测器接收到的背景噪声光子数,表示背景亮度,表示望远镜有效接收面积,表示接收视场角,表示窄带滤光片带宽,表示测距系统光学效率,表示探测器量子效率。在晴朗的白天时,式中,采用窄带滤光片,其带宽,天琴台站使用的激光测距望远镜有效接收面积,接收视场,测距系统光学效率,超导探测器的量子效率。故在晴朗白天条件下,激光测距系统接收到的背景噪声光子数。
天琴激光测距台站已经具备常规的卫星激光测距能力,使用小功率进行卫星激光测距主要目的在于分析激光测距系统的极限探测能力。本节采用4 mJ(0.4 W)、0.4 mJ(0.04 W)和0.15 mJ(0.015 W)的单脉冲能量(平均功率)对地球同步轨道卫星compass i3进行夜晚激光测距实验。实验时探测器工作温度2.30 K,4像元单光子探测器的偏置电流分别为12.0 A、12.5 A、12.5 A和12.5 A,各个像元的暗计数分别为138 cps、102 cps、147 cps和132 cps。

图4 0.15 mJ(0.015 W)单脉冲能量同步卫星激光测距
Fig. 4 0.15 mJ(0.015 W) monopulse energy synchronous satellite laser ranging

图5 0.4 mJ(0.04 W)单脉冲能量同步卫星激光测距
Fig. 5 0.4 mJ(0.04 W) monopulse energy synchronous satellite laser ranging

图6 4 mJ(0.4 W)单脉冲能量同步卫星激光测距
Fig. 6 4 mJ(0.4 W) monopulse energy synchronous satellite laser ranging
目标卫星 | 测距时长(s) | 目标仰角(°) | 激光单脉冲能量(mJ) | 重复频率(Hz) | 功率(W) | 回波率(photons/s) | 滤光片带宽(nm) |
---|---|---|---|---|---|---|---|
Compassi3 | 667 | 55 | 4 | 100 | 0.4 | 4.2 | / |
Compassi3 | 437 | 55 | 4 | 100 | 0.4 | 4.6 | / |
Compassi3 | 122 | 55 | 0.4 | 100 | 0.04 | 2.1 | / |
Compassi3 | 177 | 55 | 0.4 | 100 | 0.04 | 2.7 | / |
Compassi3 | 209 | 55 | 0.15 | 100 | 0.015 | 0.4 | / |
Compassi3 | 400 | 55 | 0.15 | 100 | 0.015 | 0.3 | / |
白天卫星激光测距的难点主要在于强天光背景噪声和望远镜的指向精度问题。针对强天光背景噪声,本文采取了空间滤波(小孔光阑)、时间滤波(距离门宽度为20 s)和光谱滤波(带宽为0.15 nm的窄带滤光片)作为噪声抑制手段;针对望远镜白天的指向误差问题,本文采取“边探测边调整”的策略,即由大逐步减小平均功率,并调整望远镜指向,使得接收有效回波率最高,最终确定指向。对同步轨道卫星qzs2和compassi5进行白天激光测距实验,实验结果如
PD = | (7) |
其中,回波光子数为s,噪声速率为vnoi,W= vnoi·(tgate-techo)·γ为回波分布范围之前的噪声光子数,tgate为门控时长,techo为回波展宽时长,γ为回波在门控中的位置,取值为0~1。n= vnoi·techo为回波分布范围内的噪声光子数。不同背景噪声强度对探测概率的影响如
目标卫星 | 测距时长(s) | 目标仰角(°) | 单脉冲能量(mJ) | 重复频率(Hz) | 功率(W) | 回波率(photons/s) | 滤光片带宽(nm) |
---|---|---|---|---|---|---|---|
Qz2 | 81 | 63 | 300 | 100 | 30 | 8.3 | 0.15 |
Qz2 | 55 | 63 | 300 | 100 | 30 | 7.2 | 0.15 |
Qz2 | 107 | 63 | 4 | 100 | 0.4 | 0.5 | 0.15 |
Qz2 | 134 | 63 | 4 | 100 | 0.4 | 0.4 | 0.15 |
Compassi5 | 123 | 58 | 300 | 100 | 30 | 7.8 | 0.15 |
Compassi5 | 713 | 58 | 4 | 100 | 0.4 | 0.2 | 0.15 |
Compassi5 | 675 | 58 | 0.4 | 100 | 0.04 | 0.02 | 0.15 |
Compassi5 | 438 | 58 | 0.4 | 100 | 0.04 | 0.01 | 0.15 |


图7 不同天区天空背景噪声亮度变化
Fig. 7 Variation of sky background noise brightness in different sky regions

图8 不同噪声强度对探测概率的影响
Fig. 8 Influence of different noise intensities on detection probability
卫星激光测距实验中选择qzs2和compassi5卫星做为激光测距实验目标。相对于测站,qzs2和太阳位置如

(a)

(b)
图9 测卫星和太阳相对位置 (a)卫星qzs2位置,(b)卫星compassi5位置
Fig. 9 he position of satellite and the Sun (a)position of qzs2, (b)position of compassi5

(a)

(b)
图10 空背景亮度 (a) 北京时间14:00(2020.12.27)天空亮度,(b)北京时间15:00(2021.02.03)天空亮度
Fig. 10 Sky background brightness (a)Beijing time 14:00 (2020.12.27) sky brightness, (b)Beijing time 15:00 (2021.02.03) sky brightness
实验中激光器发射能量衰减为0.4 mJ(0.04 W),其能量测量和激光测距回波残差散点图如


图11 0.4 mJ白天同步卫星激光测距
Fig. 11 Satellite laser ranging daytime with 0.4 mJ


图12 4 mJ白天同步卫星激光测距
Fig. 12 Satellite laser ranging daytime with 4 mJ
上文分别在不同功率条件下对同步轨道卫星于白天和夜晚进行了激光测距实验。通过降低系统发射平均功率的方法,分析天琴激光测距台站卫星激光测距探测能力极限。白天以同步轨道卫星为测距目标,成功实现测距的最小平均功率为0.04 W,夜晚以同步轨道卫星为测距目标,成功实现测距的最小平均功率为0.015 W。与国内外其他激光测距台站相比,本文所分析的系统发射单脉冲能量并非最小,其他测站使用的单脉冲能量量级在J~mJ之间,但是国外激光测距系统的重复频率为kHz,而本系统的激光重复频率仅为100 Hz,因此本文系统的平均功率最低。
国外kHz激光测距系统的单脉冲能量已经达到J量级,但正是由于高重复频率的原因,系统发射激光的平均功率并不低,达到百毫瓦量级,故可以实现远距离激光测距,相关数据统计如
单脉冲 能量 | 激光重复频率 | 平均功率 | 测量距离 | 测站/系统 |
---|---|---|---|---|
100 J | 2 kHz | 200 mW | 22,000 km |
SLR 2000 |
400 J | 2 kHz | 80 mW | 22.000 km |
Graz |
0.6 mJ | 10 kHz | 6 W | 36,000 km |
Shanghai |
400 mJ | 100 Hz | 40 W | 36,000 km | Yunnan |
150 J | 100 Hz | 15 mW | 36,000 km | TianQin |
单脉冲能量 | 激光重复频率 | 平均功率 | 测量距离 | 测站/系统 | 实验时段 |
---|---|---|---|---|---|
400 mJ | 100 Hz | 40 W | 36 000 km |
云南天文 | 白天 |
190 J | 10 kHz | 1.9 W | 36 000 km |
上海天文 | 白天 |
3 mJ | 1 kHz | 3 W | 20 000 km | 长春人卫站 | 白天 |
80 mJ | 200 Hz | 16 W | 1 500 km空间碎片 |
Gra | 白天 |
400 J | 2 kHz | 0.8 W | 6 000 km Lageos |
Gra | 白天 |
1.5 mJ | 10 kHz | 15 W | 36 000 km |
TROS100 | 白天 |
7.8 J | 1 MHz | 7.8 W | 36 000 km |
Gra | 白天 |
0.4 mJ | 100 Hz | 40 mW | 36 000 km | TianQin | 白天 |
为分析研究天琴计划激光测距台站卫星激光测距探测能力极限,在白天和夜晚分别采用小功率进行实验,成功实现对同步轨道卫星(轨道高度约36,000 km)的测量。通过降低单脉冲能量的策略,分析研究了平均功率对激光测距探测能力的影响,并总结出白天和夜晚的极限激光发射功率。天琴激光测距台站已经具备白天和夜晚全时段卫星激光测距能力,今后可更好地服务于深空目标探测任务。
References
LI Zhu-Lian, Zhai Dong-Sheng, Zhang Hai-Tao, et al. Superconductivity detector applied to daytime satellite laser ranging experiment and research [J]. Infrared and Laser Engineering, 2020, v.49;No.310(08):127-132. (李祝莲,翟东升,张海涛,等. 基于超导探测器的白天卫星激光测距试验与研究[J]. 红外与激光工程), 2020, v.49;No.310(08):127-132. [百度学术]
LI Ming, XUE Li, HUANG Chen, et al. Estimation of detection range for space debris laser ranging system based on efficient echo probability [J]. Optics and Precision Engineering, 2016, 24(2): 260-267. [百度学术]
李明, 薛莉, 黄晨, 等. 基于有效回波概率估计空间碎片激光测距系统作用距离[J]. 光学精密工程, 2016, 24(2): 260-267. [百度学术]
DENG Hua-Rong, LONG Ming-Liang, ZHANG Hai-Feng,et al. Experiment of satellite laser ranging in daytime based on 1064nm wavelength [J]. Infrared and Laser Engineering, 2020, v.49;No.312(10):134-139. 邓华荣,龙明亮,张海峰,等. 1 064 nm波长白天卫星激光测距试[J]. 红外与激光工程, 2020, v.49;No.312(10):134-139. [百度学术]
LONG Ming-Liang, ZHANG Hai-Feng, MEN Lin-Lin,et al. Satellite laser ranging at 10 kHz repetition rate in all day [J]. J. Infrared Millim.WavesJ. Infrared Millim.Waves , 2020,39(06):778-785. [百度学术]
龙明亮,张海峰,门琳琳,等. 10 kHz重复率全天时卫星激光测距[J].红外与毫米波学报, 2020,39(06):778-785. 10.11972/j.issn.1001-9014.2020.06.016 [百度学术]
LI Xin, WANG Pei-Yuan, ZOU Tong, et al, Experiment on kHz laser ranging at Wuhan satellite laser ranging station [J]. High Power Laser and Particle Beams, 2011,23(02):367-370. [百度学术]
李欣,王培源,邹彤,等. kHz激光器在武汉卫星观测站的测距实验[J]. 强激光与粒子束, 2011,23(02):367-370. [百度学术]
Kirchner G, Koidl F. Laser Ranging to Space Debris from Graz Laser Station [J]. VGI, 2015, 103(2-3):151-155. [百度学术]
Kucharski D, Otsubo T, Kirchner G, et al, Spin axis orientation of Ajisai determined from Graz 2 kHz SLR data [J]. Advances in Space Research, 2010, 46(3):251-256. [百度学术]
Kucharski D, Kirchner G , Koidl F . Spin parameters of nanosatellite BLITS determined from Graz 2 kHz SLR data [J]. Advances in Space Research, 2011, 48(2):343-348. [百度学术]
DENG Hua-Rong, ZHANG Hai-Feng, LONG Ming-Liang, et al. 4 kHz Repetition rate satellite laser ranigng system and its application [J]. Acta Optica Sinica, 2019. 39(003), 225-231. (邓华荣, 张海峰, 龙明亮,等. 4 kHz重复频率卫星激光测距系统及其应用[J]. 光学学报), 2019. 39(003), 225-231. [百度学术]
Steindorfer M A, Kirchner G, Koidl F , et al. Daylight space debris laser ranging [J]. Nature Communications, 2020, 11(1), 1-6. [百度学术]
Degnan J J, Mcgarry J F. SLR2000: Eyesafe and autonomous single photoelectron satellite laser ranging at kilohertz rates [J]. Proceedings of Spie the International Society for Optical Engineering, 1997, 3218:63-77. [百度学术]
Isyanova Y, Wall K F, Flint J H, et al. High-power, short-pulse, compact SLR2000 laser transmitter [J]. 2004, p. MB14. [百度学术]
Degnan J J. SLR2000 technical overview, status, and schedules [J]. 2001, June 1. [百度学术]
Lim H C, Sung K P, Yu S Y, et al. Satellite laser ranging system at Geochang station[J]. Journal of Astronomy and Space Science, 2018, 35(4):253-261. [百度学术]
Lim H, Choi M, Park E, et al. Technical Aspects and Progress of Korean SLR Systems [C]// 18th International Workshop on Laser Ranging. 2013. [百度学术]
Choi E J, Bang S C, Sung K P, et al. Design and Development of High-Repetition-Rate Satellite Laser Ranging System [J]. Journal of Astronomy and Space Science, 2015, 32(3):209-219. [百度学术]
Jianwei Mei, Yan-Zheng Bai, Jiahui Bao, et al. The TianQin project: current progress on science and technology [J]. Progress of Theoretical and Experimental Physics, 2021, (5), 05A107. [百度学术]
Tang Mei-Rong, Li Zhu-Lian, Li Yu-Qiang, et al. Influence of Atmospheric Turbulence Intermittency on Return Photon Number in Laser Ranging [J]. Acta Optica Sinica, 2017. 037(002), 1-7. (唐美荣, 李祝莲, 李语强, 等. 大气湍流间歇性对激光测距回波光子数的影响[J]. 光学学报), 2017. 037(002), 1-7. [百度学术]
Han Guang-Yu, Ju Feng, Guo Jin, et al. Analysis and realization of day-time ranging of SLR [J]. Chinese Journal of Scientific Instrument, 2012. (04), 885-890.韩光宇, 瞿锋, 郭劲, 等. 卫星激光测距中白天测距的分析与实现[J]. 仪器仪表学报, 2012. (04), 885-890. [百度学术]
J. J. Degnan, “SLR 2000” in Satellite Laser Ranging in the 1990’s: Report of the 1994 Belmont Workshop [C]. NASA Conference Publication, 1994, 3283, pp. 101-106. [百度学术]
Kucharski D, Kirchner G, Lim H Cet al. Spin parameters of High Earth Orbiting satellites Etalon-1 and Etalon-2 determined from kHz Satellite Laser Ranging data [J]. Advances in Space Research, 2014, 54(11), 2309–2317. [百度学术]
Kirchner G, Koidl F. Graz kHz SLR system: design, experiences and results[C]//Proc. 14th Int. Workshop Laser Ranging. 2004: 501-505. [百度学术]
Guo T, Wang P, Li X, et al. Progress of the satellite laser ranging system TROS1000[J]. Geodesy and Geodynamics, 2015, 6(1): 67-72. [百度学术]
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. [百度学术]