2. 中国科学院大学, 北京 100049;
3. 中国科学院空间目标与碎片观测重点实验室, 江苏 南京 210008
2. University of Chinese Academy of Sciences, Beijing 100049, China;
3. Key Laboratory of Space Object and Debris Observation, Chinese Academy of Sciences, Nanjing 210008, China
激光测距技术是一项综合技术,涵盖激光、电子、天文测量和卫星轨道计算等多个学科领域,在实际工作中是一个非常复杂的过程,受目标特性、望远镜指向误差、距离、接收孔径大小等多种因素的影响[1-2]。在众多影响因素中,空间目标的轨道预报准确性和稳定性直接影响测距望远镜发出的激光脉冲能否持续到达其表面,接收端的单光子探测器能否及时探测到回波光子对激光测距的成功率有很大影响[3]。
卫星激光测距系统广泛采用了距离门控技术,即时间滤波技术。该技术能有效防止噪声的干扰,是在预期回波到达时刻之前一瞬间给单光子探测器送开门脉冲,这样使受到噪声干扰的概率减到最小,使淹没在众多噪声中的单光子探测成为可能[4]。
非合作目标的激光测距预报,即载入跟踪望远镜的轨道预报基于双行根数(Two Line Elements, TLE),经过SGP4 (Simplified General Perturbations Version 4)轨道衍生器外推出来的①,在对目标进行跟踪时产生的较大偏差主要体现在望远镜视场中的位置偏差以及测站与目标之间的距离偏差。视位置偏差过大导致激光脉冲不能击中目标;测站和目标的预报距离偏差过大,使预期回波到达时刻与实际回波到达时刻出现大的偏差,导致单光子探测器的距离门开启期间探测不到回波光子。
① http://mobile.celestrak.com/NORAD/documentation/spacetrk.pdf
本文提出非合作目标激光测距预报实时修正方法,先利用双行根数原始预报对可见的非合作目标跟踪几分钟,获取其在望远镜跟踪视场中的脱靶量数据信息,根据这些信息解算时间根数偏差,并基于该时间根数偏差对轨道预报进行修正后,重新载入望远镜继续跟踪目标。激光测距卫星(例如Ajisai, Beaconc, Starlette, Stella等)具有双行根数和国际激光测距服务机构的官方预报星历(Consolidated Prediction Format, CPF)两种预报,文[5]将中低轨卫星的CPF星历与利用全球激光跟踪站的观测数据定轨出的精密轨道进行了比较,分析CPF星历精度在米级以内;文[6]分析了SGP4模型处理近地目标时,定轨精度在百米量级。为了排除望远镜指向误差等测距系统误差的影响,在理论上证明该方法的有效性,本文把中低轨道高度的激光测距卫星视为非合作目标进行研究分析,把某一时间段可见的激光卫星的CPF预报作为参考值,即模拟观测值。目标的脱靶量由同一时刻双行根数预报和CPF预报中的方位角和高度角偏差量代替,然后进行时间根数偏差的解算和轨道预报的修正。
1 激光测距预报的生成 1.1 时间系统的转换协调世界时(UTC)是一种修正过的原子时,空间目标的状态以及观测数据都使用UTC作为输入、输出的时间参数。格林尼治平恒星时(Greenwich sidereal time, GST) θG是坐标系转换的一个关键参数,要由协调世界时计算出θG,详见文[7]。
1.2 TEME坐标系与站心地平坐标系之间的转换利用预测模型外推出来的空间目标位置信息是在真赤道平春分点坐标系(True Equator Mean Equinox Frame, TEME)中表达的,而激光测距望远镜使用的轨道预报是在站心地平坐标系中表达的,所以要实现由真赤道平春分点坐标系向站心地平坐标系的转换,如图 1:
$ {\rm{图 1中, }}{R_{\rm{z}}}\left( {{{\bar \theta }_{\rm{G}}}} \right) = \left( {\begin{array}{*{20}{c}} {{\rm{cos}}{{\bar \theta }_{\rm{G}}}}&{{\rm{sin}}{{\bar \theta }_{\rm{G}}}}&0\\ {-{\rm{sin}}{{\bar \theta }_{\rm{G}}}}&{{\rm{cos}}{{\bar \theta }_{\rm{G}}}}&0\\ 0&0&1 \end{array}} \right); $ | (1) |
$ (MLT) = {R_{\rm{x}}}\left( {\frac{{\rm{ \mathsf{ π} }}}{2}-\phi } \right){R_{\rm{z}}}\left( {\frac{{\rm{ \mathsf{ π} }}}{2} + \lambda } \right), $ | (2) |
其中,
$ {R_{\rm{x}}}\left( {\frac{{\rm{ \mathsf{ π} }}}{2}-\phi } \right) = \left( {\begin{array}{*{20}{c}} 1&0&0\\ 0&{{\rm{sin}}\phi }&{{\rm{cos}}\phi }\\ 0&{-{\rm{cos}}\phi }&{{\rm{sin}}\phi } \end{array}} \right), $ | (3) |
$ {R_{\rm{z}}}\left( {\frac{{{\rm{ \mathsf{ π} }}}}{2} + \lambda } \right) = \left( {\begin{array}{*{20}{c}} {-{\rm{sin}}\lambda }&{{\rm{cos}}\lambda }&0\\ {-{\rm{cos}}\lambda }&{-{\rm{sin}}\lambda }&0\\ 0&0&1 \end{array}} \right), $ | (4) |
ϕ和λ分别为激光测距望远镜所在位置的纬度和经度;
当对中低轨空间目标进行激光测距时,测站和目标之间往往具有几千千米的距离,在激光脉冲的上行和下行期间,空间目标相对测站发生了较大变化,因此对于激光测距一般用发射时刻作为时标生成跟踪文件,即望远镜的实时指向有个前指量。针对本文的应用,由于下行光行差的原因,目标可见对应的预报应该滞后,即将激光接收方向作为模拟观测值。CPF预报中,前指量的算法如下:
对于某一时刻T,空间目标和测站在地固坐标系中的位置矢量为
$ {\rm{d}}A = A-{A_{\rm{B}}} = {\rm{ta}}{{\rm{n}}^{-1}}\left( {-\frac{{{\rho _y}}}{{{\rho _x}}}} \right) - {\rm{ta}}{{\rm{n}}^{ - 1}}\left( { - \frac{{{\rho _{{\rm{B}}y}}}}{{{\rho _{{\rm{B}}x}}}}} \right), $ | (5) |
$ {\rm{d}}E = E-{E_{\rm{B}}} = {\rm{ta}}{{\rm{n}}^{-1}}\frac{{{\rho _z}}}{{\sqrt {{\rho _x}^2 + {\rho _y}^2} }}-{\rm{ta}}{{\rm{n}}^{ - 1}}\frac{{{\rho _{{\rm{B}}z}}}}{{\sqrt {{\rho _{{\rm{B}}x}}^2 + {\rho _{{\rm{B}}y}}^2} }}. $ | (6) |
利用空间几何相关知识,下行光路的光行差引起的激光接收滞后量为dA,dE。
2 时间根数偏差的解算以及预报的修正 2.1 理论表述图 2中,阴影部分为地心天球坐标系O-XYZ的基本平面(XY平面),接近J2000.0平赤道面,X轴指向接近J2000.0平春分点,Z轴指向地球的北极。r为空间目标的轨道半径,
$ r = \frac{{a\left( {1-{e^2}} \right)}}{{1 + e{\rm{cos}}\left( {\theta-\omega } \right)}}, $ | (7) |
其中,a为半长径;e为偏心率;近地点幅角ω和极坐标变量θ都是从升交点方向起量的。当有摄动时,椭圆轨道随时间变化,升交点方向也在变化,ω从变化的升交点方向起量,而极坐标变量θ仍从一个定义的不变方向起量。
平近点角M定义为
$ M = n\left( {t-\tau } \right), $ | (8) |
从近地点开始起量,其中,n为平运动角速度。空间目标的平近点角是一个与时间相关的参数,能够确定对应时刻目标在椭圆轨道上的位置。当t=τ时,M=0,相应地,空间目标在近地点处,r达到极小值,τ是空间目标过近地点的时刻。
空间目标在轨运行时会受到多种摄动力的影响,包括地球的非球形引力、日月引力、大气阻力、太阳辐射压力等,对处于中低轨道的空间目标,受到大气阻力的影响比较大。双行根数考虑了地球非球形引力、日月引力的长期和周期摄动影响,以及大气阻力模型产生的引力共振和轨道衰退。双行根数是“平均”根数,它用特定的方法去掉了周期扰动项,预测模型SGP4用同样的方法重构了周期扰动项。在观测中,空间目标往往在沿迹方向与望远镜的指向中心有提前或者滞后的偏差。该偏差的主要原因是SGP4模型采用的大气模型简化为B*一个阻力参数,沿迹方向的偏差量可以简化归因于时间根数τ的偏差。
如果在跟踪目标时,适当调整τ的值,可以减小空间目标与望远镜指向中心的偏移量,修正空间目标与测站的距离量,进而提高跟踪的稳定性。
时间根数τ偏差的解算算法如下,记空间目标在站心地平坐标系中的方位角为A,高度角为E,离测站距离值为s,实测的序列记为(ti, Ai, Ei, si), i=1, 2, 3, …, n,在计算预报时适当调整τ的值,令τ + Δτ取代τ,Δτ是调整量,空间目标在真赤道平春分点坐标系中相应时刻点的位置
$ \overrightarrow r \left( t \right) = {R_{\rm{z}}}(\overline \theta {_{\rm{G}}}\left( t \right))\overrightarrow R {\rm{ }}(t + \Delta \tau ), \overrightarrow \rho \left( t \right) = \left( {MLT} \right)[\overrightarrow r \left( t \right)-{\overrightarrow r ^{{\rm{sta}}}}]; $ |
由此可得预报序列(ti, A′i, E′i, s′i; Δτ),i=1, 2, 3, …, n。记ΔAi=A′i-Ai,ΔEi=E′i-Ei,当Δτ的选值使
$ RM{S^*} = \sqrt {\frac{{{\sum _i}{{(\Delta {A_i}{\rm{cos}}{E_i})}^2} + {\sum _i}\Delta {E_i}^2}}{n}} $ | (9) |
为跟踪时间段上的极小值时,Δτ为预报的时间根数偏差。上述问题是一个复杂的非线性问题,直接求解比较困难,考虑到时间根数偏差精确到毫秒,因此一个简单直接的方法就是在一定范围内进行搜索。
当空间目标的指向信息得到修正时,距离信息也得到了相应的修正,进而能够提高单光子探测器的距离门控的精度,增大接收回波光子的概率。
2.2 实例分析Ajisai, Beaconc, Starlette, Stella等激光测距卫星位于近地轨道,与大部分的非合作目标的轨道高度相仿,空间环境类似,能够很好地替代仅有双行根数预报的非合作目标进行研究分析。
本文以Stella激光测距卫星为例进行分析。Stella卫星的轨道高度为792.3 × 803.1 km,处于非合作目标的密集区。当对空间目标进行激光测距时,测站和目标之间往往具有上千千米的距离,在激光脉冲的上行和下行期间,空间目标相对测站发生了较大变化,指向滞后量如图 3,横轴为日积秒,纵轴分别为方位角偏差和高度角偏差。从图 3可以看出,滞后量波动幅度比较大,最高时可达十几角秒,是解算时间根数偏差时不可忽视的因素。
图 4显示了在Stella的某次过境观测弧段,分别用前1 min,2 min,3 min以及全弧段跟踪数据进行时间根数搜索的结果,图中横轴为时间根数偏差Δτ,纵轴为对应的RMS*,时间根数搜寻范围为[-100, 200]。由图 4可见,无论用前1 min,2 min,3 min还是全弧段跟踪数据进行搜寻,RMS*都有一个极小值,对应时间根数偏差分别为26 ms,27 ms,29 ms,32 ms并无显著差异。对于其它的激光测距卫星,根据不同时长的跟踪弧段,搜索的时间根数偏差结果如表 1。
激光卫星 | 轨道高度/km | 过境时长 | 前一分钟/ms | 前两分钟/ms | 前三分钟/ms | 全弧段/ms |
Ajisai | 1 479 × 1 497 | 12 min 20 s | 35 | 35 | 35 | 34 |
Beaconc | 928.6 × 1 305 | 8 min 30 s | 30 | 28 | 26 | 23 |
Starlette | 805.5 × 1 107 | 6 min 9 s | 36 | 38 | 40 | 43 |
通过表 1可以看出,上述3个激光测距卫星在过境可见的时间段内,用前1 min,2 min,3 min以及整个弧段解算的时间根数偏差无显著差异。根据搜索时间根数偏差的思想可知,利用的跟踪数据段越长,解算的结果越精确,利用全弧段跟踪数据解算的结果在理论上最精确。考虑到中低轨目标过境时间短等问题,选择前1 min数据解算的时间根数偏差适用于整个弧段。
激光卫星Stella的轨道预报经过相应的修正后,视位置偏差明显减小,且波动幅度也很小。图 5 (a)中,实线表示Stella的双行根数原始预报中的方位角与CPF参考值的残差随观测时间的变化,下方的虚线分别表示双行根数预报经过26 ms,27 ms,29 ms,32 ms修正的方位角与CPF参考值的残差随观测时间的变化。类似地,图 5 (b)中,实线表示Stella的双行根数原始预报中的高度角与CPF参考值的残差随观测时间的变化,中间的虚线分别表示双行根数预报经过26 ms,27 ms,29 ms,32 ms修正的高度角与CPF参考值的残差随观测时间的变化。图 6是以ΔAcosE为横坐标,ΔE为纵坐标画出了整个跟踪弧段以0 ms (未修正)和26 ms(前一分钟)分别修正后的视位置偏差量的变化。综合以上可以看出,双行根数预报经过时间根数偏差改正后,视位置偏差量明显减小,且在某一固定值附近小范围波动,波动范围在几角秒之内。
表 2展示了处于近地轨道的4颗激光卫星的双行根数预报在未修正及基于前一分钟,整个弧段解算出的时间根数偏差分别修正后与CPF参考轨道在方位角、高度角的偏差均值以及相应的位置偏差RMS*。
激光星 | RMS0* | RMS1* | RMS2* | ||||||
Ajisai | -14.730 1 | 15.160 2 | -20.168 6 | 4.003 7 | -20.012 4 | 4.322 0 | 19.206 6 | 13.580 2 | 13.582 6 |
Beaconc | -13.451 4 | 4.064 5 | -9.397 2 | -8.887 3 | -10.350 1 | -5.871 8 | 12.174 0 | 8.769 7 | 4.160 2 |
Starlette | -69.325 8 | 15.458 1 | -62.306 6 | -0.363 6 | -60.936 5 | -3.439 9 | 89.690 7 | 84.207 6 | 87.890 5 |
Stella | 13.286 7 | 18.568 2 | 4.530 2 | 5.710 9 | 2.503 2 | 2.746 0 | 13.693 9 | 1.721 3 | 2.403 6 |
从图 7可以看出,激光卫星Stella的距离预报分别以26 ms,27 ms,29 ms,32 ms为时间根数偏差修正后,距离偏差由最大的二百多米减小到几十米,并且变化范围更稳定。
表 3展示了处于近地轨道的4颗激光卫星的轨道预报在未修正及基于前一分钟、整个弧段解算出的时间根数偏差分别修正后相较于CPF参考轨道在预报距离方面的偏差均值D0, D1, D2以及相应的RMS,即RMS0, RMS1, RMS2。可以看出,预报距离偏差均值由几百米减小到了几十米,甚至几米;均方根值RMS也减小了一个量级,说明偏差稳定性显著增强。上述结果说明基于时间根数偏差修正双行根数预报能够提高单光子探测器距离门控的精度。
激光星 | D0 | D1 | D2 | RMS0 | RMS1 | RMS2 |
Ajisai | -215.89 | -49.78 | -54.52 | 78.80 | 9.68 | 11.60 |
Beaconc | -154.19 | 8.72 | -29.21 | 46.16 | 9.11 | 4.35 |
Starlette | -214.18 | -4.56 | 36.19 | 56.31 | 14.84 | 7.67 |
Stella | -180.97 | -37.51 | -4.99 | 62.42 | 18.58 | 8.67 |
本文介绍了一种非合作目标激光测距预报实时修正的方法。在搜索时间根数偏差的过程中,以RMS*取得极小值时的Δτ值为时间根数偏差,用于空间目标的双行根数预报修正。用不同时长的跟踪数据搜寻出的时间根数偏差无显著差异,用前一分钟数据解算的时间根数偏差可以适用于整个弧段。修正后的双行根数预报更接近空间目标的实测运行轨迹且预报距离偏差值均由几百米减小到了几十米,甚至几米,能够提高单光子探测器距离门控的精度。
本文选取若干个空间目标的双行根数预报经过修正后,视位置值和预报距离值的偏差得到了明显的改善。
[1] | 翟东升, 汤儒峰, 李祝莲, 等. 激光脉冲特性对漫反射激光测距系统精度的影响研究[J]. 天文研究与技术, 2016, 13(3): 326–332 DOI: 10.3969/j.issn.1672-7673.2016.03.008 |
[2] | 唐美荣, 和丽娟, 翟东升, 等. 激光测距系统回波光子分布特性研究[J]. 天文研究与技术, 2017, 14(1): 32–38 |
[3] | 李语强, 李荣旺, 李祝莲, 等. 空间碎片激光测距应用研究[J]. 红外与激光工程, 2015, 44(11): 3324–3329 DOI: 10.3969/j.issn.1007-2276.2015.11.027 |
[4] | 李祝莲, 熊耀恒, 何妙婵, 等. 云南天文台人造卫星激光测距系统原理[J]. 天文研究与技术——国家天文台台刊, 2008, 5(3): 248–252 |
[5] | 陈国平, 何冰, 张志斌, 等. CPF星历精度分析[J]. 中国科学院上海天文台年刊, 2010(31): 35–44 |
[6] | 韦栋, 赵长印. SGP4/SDP4模型精度分析[J]. 天文学报, 2009, 50(3): 332–339 DOI: 10.3321/j.issn:0001-5245.2009.03.010 |
[7] | 黄珹, 刘林. 参考坐标系及航天应用[M]. 北京: 电子工业出版社, 2015. |
[8] | 刘林, 胡松杰. 航天器定轨理论与应用[M]. 北京: 电子工业出版社, 2015. |