快速检索        
  武汉大学学报·信息科学版  2016, Vol. 41 Issue (1): 93-99

文章信息

叶茂, 李斐, 鄢建国, 郝卫峰, 邵先远
YE Mao, LI Fei, YAN Jianguo, HAO Weifeng, SHAO Xianyuan
GRAIL月球重力场模型定轨性能分析
Orbit Determination Ability Analysis of the GRAIL Gravity Model
武汉大学学报·信息科学版, 2016, 41(1): 93-99
Geomatics and Information Science of Wuhan University, 2016, 41(1): 93-99
http://dx.doi.org/10.13203/j.whugis20150008

文章历史

收稿日期: 2015-01-13

GRAIL月球重力场模型定轨性能分析
叶茂1, 李斐1,2, 鄢建国1, 郝卫峰2, 邵先远1     
1. 武汉大学测绘遥感信息工程国家重点实验室, 湖北 武汉, 430079;
2. 武汉大学中国南极测绘研究中心, 湖北 武汉, 430079
摘要: 由最新GRAIL跟踪数据解算的月球重力场模型相较于之前的模型在分辨率和精度上有了巨大的飞跃,本文采用月球重力场模型GRGM660PRIM,通过与历史LP150Q和SGM150j模型对比,从轨道预报和精密定轨两个角度定量分析了GRGM660PRIM的定轨性能。计算结果表明,对于卫星轨道寿命预报,GRGM660PRIM并未体现出较大优势,与LP150Q和SGM150j相当;对于大倾角的Apollo 16子卫星的精密定轨,GRGM660PRIM体现出其高精度、高阶次的优势,双程Doppler测速残差RMS由1.36~1.48 mm• s降低至0.722 mm• s,三程Doppler测速残差RMS由1.67~1.69 mm• s降低至1.2 mm• s,定轨精度得到显著的提高。
关键词: 月球重力场     轨道预报     精密定轨     GRAIL     Apollo 16    
Orbit Determination Ability Analysis of the GRAIL Gravity Model
YE Mao1, LI Fei1,2, YAN Jianguo1, HAO Weifeng2, SHAO Xianyuan1     
1. State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan 430079, China;
2. Chinese Antarctic Center of Surveying and Mapping, Wuhan University, Wuhan 430079, China
First author: YE Mao, PhD candidate, specializes in planetary spacecraft precision orbit determination, gravity recovery and software development. E-mail: mye@whu.edu.cn.
Foundation support: The National Natural Science Foundation of China, Nos. 41174019, 41374024; Hubei Province Natural Science Foundation Innovation Group Project, No. 2015CFA011; Wuhan University PhD Short-Time Mobility Program.
Abstract: The new GRAIL gravity model GRGM660PRIM has advanced in the aspect of resolution and degree when compared with historical lunar gravity field models. In this paper, we present a quantitative analysis of the orbit determination ability of the GRGM660PRIM in orbit prediction and precise orbit determination (POD). In a comparison with the LP50Q and SGM150j models, results show that the GRGM660PRIM is close to the level of the LP150Q and SGM150 in satellite lifetime prediction; in the Apollo 16 subsatellite POD, GRGM660PRIM shows its advantages at high degrees and precision, the residuals RMS of two Doppler decreases from 1.36~1.48 mm• s to 0.722 mm• s, and the three way Doppler residuals RMS decrease from 1.67~1.69 mm• s to 1.2 mm• s;the precision of orbit determination improves significantly. These conclusions can provide a reference for the upcoming Chinese Chang'E-5 mission and further lunar missions.
Key words: lunar gravity field;orbit prediction;precise orbit determination;GRAIL;Apollo 16    

月球重力场是绕月探测器的重要摄动源,其精度对探测器的精密定轨及其着陆具有重要意义[1, 2]。美国于2011年9月发射了GRAIL卫星,分别解算了高达420阶次的月球重力场模型GL0420A、660阶次的模型GL0660B和GRGM660PRIM,以及900阶次的模型GL0900D和GRGM900C[3, 4, 5, 6, 7]。GRAIL之前的重力场模型主要为LP系列模型[8, 9]和SELENE系列模型[10, 11, 12],前者由于缺少远月面跟踪数据,在远月面的定轨性能不够;日本SELENE卫星采用高低跟踪模式,首次获得远月面重力场信号,提高了远月面重力场精度。但是由于LP和SELENE卫星轨道高度为100 km左右,无法进一步提取重力场的高阶信号,故其对绕月低轨卫星的定轨性能较差。GRAIL卫星采用了类似地球GRACE卫星跟踪卫星的模式,正常任务阶段轨道高度为50 km左右,扩展任务阶段为30 km左右,获取了全月球均匀分布的高精度星间Ka波段测量数据。与之前的重力场模型相比,GRAIL获取的重力场模型改正精度达3个量级,部分阶次(50~120阶)改进接近6个量级。中国于2007年成功发射嫦娥一号卫星,国内学者利用嫦娥一号的跟踪数据解算了50阶次的月球重力场CEGM-01[13],之后进一步融合LP、SELENE和其他历史跟踪数据解算得到100阶次的CEGM02[14],在重力场长波部分的精度有一定程度的提高。

高精度、高阶次的GRAIL月球重力场模型主要用于月球内部构造的研究,而其定轨性能目前并没有得以全面分析。因此,本文使用了660阶次的GRGM660PRIM模型,通过与历史LP150Q和SGM150j模型比较,从轨道预报和精密定轨两个角度[15, 16]分析了GRAIL重力场模型的定轨性能,着重分析了GRAIL对大倾角绕月卫星的定轨能力。

1 不同轨道高度、不同轨道倾角的近圆绕月探测器轨道仿真分析

为分析GRAIL重力场模型对绕月探测器轨道演化的影响,本文分别采用了20 km、50 km、100 km的轨道高度和10°、45°、89°、135°、170°的轨道倾角进行了轨道仿真计算。轨道计算使用GEODYNII[17]软件完成,采用的动力学模型为月球中心引力、月球非球形引力摄动、N体摄动、太阳光压摄动、固体潮摄动、相对论效应等。初始轨道根数时间为2017-10-01T00:00:00,月球平均半径取1 738 km,初始轨道偏心率为0.000 5,积分间隔为30 s。

表1~表5为采用LP150Q、SGM150j和GRGM660PRIM重力场模型在上述初始条件下的仿真计算结果,表1~表5中,Δa、Δe、Δrp、Δra表示运行终止时刻与起始时刻轨道半长轴、偏心率、近月点、远月点的差值;Δi表示运行终止时刻与起始时刻的平均变化率;Δt表示卫星在轨运行寿命。在仿真计算中,GRGM660PRIM重力场模型分别截断至500阶、420阶、150阶,发现截断后对轨道预报影响不大,但通过截断可以减少计算时间。限于篇幅所限,本文不列出全部计算结果。表1~表5为GRGM660PRIM截断至420阶次的结果。

表 1 倾角89°,不同重力场模型和不同高度下的卫星轨道根数变化 Tab. 1 Variation of Orbit Elements with Same Inclination 89° but Different Elevation for Various Gravity Field Models
100 km50 km20 km
LP150QSGM150jGRGMLP150QSGM150jGRGMLP150QSGM150jGRGM
Δa/m-907.855-911.048-807.057-780.975-811.129-1 768.72-802.117-734.218-958.376
Δe0.053 4800.053 4740.053 5300.027 1560.027 2260.026 5170.011 2740.011 3120.010 741
Δi/(°·d-1)-0.000 042-0.000 042-0.000 0520.010 6370.010 6350.010 817-0.019 43-0.019 40-0.021 114
Δrp/km-99.156-99.147-99.152-49.314-49.468-49.134-20.612-20.612-19.830
Δra/km97.34097.32597.53847.75247.84645.59719.00819.14417.914
Δt/d175.969 1175.969 1175.969 158.864 958.864 958.700 711.847 611.847 611.771 2

表 2 倾角10°,不同重力场模型和不同高度下的卫星轨道根数变化 Tab. 2 Variation of Orbit Elements with Same Inclination 10° but Different Elevation for Various Gravity Field Models
100 km50 km20 km
LP150QSGM150jGRGMLP150QSGM150jGRGMLP150QSGM150jGRGM
Δa/m-287.566-115.455-100.151343.387308.615322.092579.092626.5981 346.589
Δe0.053 8350.053 9960.053 8530.027 8000.027 8360.027 7380.011 3080.011 2000.011 790
Δi/(°·d-1)-0.009 378-0.009 868-0.009 678-0.031 677-0.029 718-0.029 6690.190 6470.189 4270.187 569
Δrp/km-99.221-99.354-99.077-49.373-49.471-49.283-19.307-19.071-19.397
Δra/km98.64699.12398.87750.06050.08949.92820.46520.32322.090
Δt/d47.754 248.243 147.998 65.740 35.740 35.663 92.757 32.680 92.913 9

表 3 倾角170°,不同重力场模型和不同高度下的卫星轨道根数变化 Tab. 3 Variation of Orbit Elements with Same Inclination 170° but Different Elevation for Various Gravity Field Models
100 km50 km20 km
LP150QSGM150jGRGMLP150QSGM150jGRGMLP150QSGM150jGRGM
Δa/m-36.62437.982101.008133.271258.866-134.856-30.10240.137244.312
Δe0.053 8950.053 9230.053 9650.027 7140.027 8540.027 4300.011 1460.010 9780.010 982
Δi/(°·d-1)0.115 4240.115 1800.115 7260.123 6170.126 5330.117 022-0.225 713-0.238 139-0.234 347
Δrp/km-99.093-99.074-99.092-49.423-49.551-49.175-19.624-19.260-19.065
Δra/km99.02099.15099.29449.68950.06848.90619.56419.34019.553
Δt/d16.518 816.602 816.522 76.045 86.045 85.969 41.466 31.470 11.393 8

表 4 倾角45°,不同重力场模型和不同高度下的卫星轨道根数变化 Tab. 4 Variation of Orbit Elements with Same Inclination 45° but Different Elevation for Various Gravity Field Models
100 km50 km20 km
LP150QSGM150jGRGMLP150QSGM150jGRGMLP150QSGM150jGRGM
Δa/m3.424-592.993-552.711-371.038-237.673-167.188-187.599.798-980.307
Δe0.054 1050.053 8060.053 7510.027 3780.027 3470.027 3970.010 7630.010 8810.010 555
Δi/(°·d-1)-0.009 059-0.002 960-0.002 959-0.004 9700.010 3750.010 356-0.072 484-0.067 011-0.065 201
Δrp/km-99.442-99.455-99.317-49.313-49.127-49.148-19.107-19.119-19.525
Δra/km99.44998.27098.21248.57048.65148.81418.73219.13917.564
Δt/d72.427 876.189 976.189 925.070 426.170 426.170 48.960 18.116 08.112 4

表 5 倾角135°,不同重力场模型和不同高度下的卫星轨道根数变化 Tab. 5 Variation of Orbit Elements with Same Inclination 135° but Different Elevation for Various Gravity Field Models
100 km50 km20 km
LP150QSGM150jGRGMLP150QSGM150jGRGMLP150QSGM150jGRGM
Δa/m-462.106244.541-491.745-359.980-192.276-748.301-423.340-294.71041.306
Δe0.053 6730.054 1110.053 7270.027 5050.027 6520.027 1730.010 9100.010 7310.011 027
Δi/(°·d-1)-0.005 135-0.006 136-0.006 179-0.122 900-0.050 227-0.051 8440.031 7610.034 3770.032 140
Δrp/km-99.088-99.224-99.215-49.528-49.629-49.313-19.599-19.157-19.344
Δra/km98.16399.71398.23248.80849.24547.81618.75218.56819.427
Δt/d67.199 071.858 771.858 710.235 814.081 914.001 71.695 51.695 51.771 9

表1~表5可以看出,随着初始轨道高度降低,飞行器在轨时间急剧减少。极轨绕月卫星相对其他倾角卫星,在轨时间相对较长,对于大倾角、低轨的绕月卫星,例如170°倾角、20 km高度卫星,在轨时间极短,轨道调整难度比较大。整体上看,顺行轨道的卫星在轨时间比逆行轨道的在轨时间长,如10°、100 km轨道的卫星在轨时间为48 d左右,而170°、100 km的在轨时间仅为16 d。

GRGM660PRIM模型对20 km、50 km、100 km高度卫星的轨道寿命预报并未表现出优越性,和LP150Q、SGM150j的预报结果相当。值得一提的是,LP150Q模型对轨道倾角为45°、135°,轨道高度为100 km、50 km的卫星的轨道寿命预报和SGM150j、GRGM660PRIM的计算结果相差比较大,其原因可能为LP150Q重力场的解算数据中缺少月球背面的轨道跟踪数据,在月球背面误差比较大。图1为LP150Q重力异常图和135°倾角、50 km高度的卫星星下点轨迹图,可以看出LP150Q在月球背面存在较大的条带误差,当卫星经过月球背面时,这些误差将直接导致预报的偏差。

图 1 135°倾角,50 km高度的卫星星下点轨迹图(南北纬45°之间的红色部分) Fig. 1 Sub-satellite Points (Between 45°N and 45°S) for 135° Inclination and 50 km Height Satellites
2 Apollo 16子卫星的精密定轨分析

GRAIL卫星采用了卫卫跟踪模式,其星间Ka波段跟踪数据精度为0.1 μm/s,相比于传统S波段大约1 mm• s的精度,提高了4个数量级,基于这些全球均匀分布的高精度Ka波段星间测距数据,得到了高阶次高精度的GRAIL重力场模型,其解算目前并未使用历史跟踪数据。作为一个检核,文献[18]利用GRAIL重力场模型对LP、SELENE和LRO探测器进行了精密定轨,得到了精度更高的重建轨道,这些精密轨道将有利于其他科学数据的精化,如激光测高数据。但是这些卫星均为倾角约90°的极轨卫星,本文采用轨道倾角为170°的Apollo 16子卫星进行了精密定轨。该子卫星是Apollo 16登月任务完成后,返回地球前由命令舱在轨释放,其初始轨道高度为100 km,在轨运行34 d后撞月,期间无任何轨道机动[9, 19]。通过对Apollo 16子卫星精密定轨分析,可进一步反映出GRAIL重力场模型对大倾角低轨卫星的定轨性能。

在对Apollo 16子卫星的精密定轨中,估计参数包括6个轨道根数、光压系数以及测量系统偏差,测量数据包括双程和三程Doppler测速。图2给出了利用不同模型得到的测速残差,表6为各残差的统计信息,图3为Apollo 16子卫星星下点覆盖图。

图 2 Apollo 16子卫星双程和三程Doppler测速残差图 Fig. 2 Two Way and Three Way Doppler Residuals of Apollo 16 Subsatellite

表 6 不同模型精密定轨残差统计信息/(mm• 5s-1) Tab. 6 Statistical Results of POD Residuals for Different Models/(mm• 5s-1)
LP150QSGM150jGRGM660PRIM
双程Doppler测速残差
均值 -7.5×10-4 1.37×10-3 3.14×10-4
均方根 1.36 1.48 0.722
三程Doppler测速残差
均值 1.37×10-3 6.8×10-4 -1.58×10-4
均方根 1.67 1.69 1.2

图 3 Apollo 16子卫星星下点覆盖图 Fig. 3 Sub-satellite Points of Apollo 16 Subsatellite

图2可以看出,利用GRGM660PRIM模型解算的残差水平明显低于LP150Q和SGM150j,双程Doppler测速残差RMS由1.36~1.48 mm• s降低至0.722 mm• s,三程Doppler测速残差RMS由1.67~1.69 mm• s降低至1.2 mm• s。LP150Q和SGM150j对Apollo 16子卫星具有相近的定轨水平,SGM150j在月球背面精度的提高并未对Apollo 16子卫星的精密定轨产生显著影响,其原因为其背面精度的提高主要来自极轨卫星,而对于大倾角卫星的精密定轨不会产生显著的影响。同时,这两个模型的残差在5月10日至5月12日期间对应大的跳变,而利用GRGM660PRIM模型解算时没有出现这一跳变现象,主要是因为在此期间Apollo 16子卫星的近月点比较低,接近10 km(见图4),GRGM660PRIM模型的阶次较高,更有利于较低轨道处的精密定轨。

图 4 Apollo 16子卫星在轨运行阶段轨道高度图 Fig. 4 Orbit Height of Apollo 16 Subsatellite

为进一步反映GRGM660PRIM的定轨性能,本文选取了观测值相对密集的5月10日至5月12日三天进行了重复弧段分析,重复弧段时长为6 h,表7给出了3个重复弧段差异在径向、沿迹方向和法向3个方向上的残差RMS值。选取上述3 d的弧段进行重复轨道比较的原因是相比其他时段,这3 d的轨道跟踪数据量较为密集。由表7的结果可以看出,GRGM660PRIM在径向方向的变化幅度最小,达到亚米级,相比于LP150Q和SGM150j解算结果,残差RMS值显著减少。5月10日至5月12日期间,Apollo 16子卫星最低轨道高度达10 km,LP150Q和SGM150j模型的阶次不足以吸收该高度轨道处包含的重力场信息,定轨精度相对较差,而GRGM660PRIM的阶次比较高,可以取得较好的定轨精度。此外,利用GRAIL模型处理小倾角卫星,可以发现仍然具有很好的结果,这与文献[20]中的倾角卫星精密定轨仿真结果有所出入,也从侧面证实了GRAIL星间跟踪数据对月球重力场模型精度的显著改进。

表 7 不同模型重复弧段差异RMS/m Tab. 7 Statistics of Orbit Overlap Difference for Different Models/m
弧段径向沿迹法向
LP150QSGM150jGRGMLP150QSGM150jGRGMLP150QSGM150jGRGM
5/100.4178.1690.1781.65317.4270.8824.42132.7506.326
5/112.7409.9810.5755.29521.3601.79226.46040.7322.666
5/125.0944.6980.38510.89510.8511.3227.38263.69518.730

综合上述分析可知,GRGM660PRIM对绕月低轨卫星具有更好的定轨性能,这点对我国即将发射的嫦娥5号探测器具有重要的参考意义。

3 结 语

本文从轨道预报和精密定轨两个角度分析了最新高精度月球重力场模型GRGM660PRIM的定轨性能,结果表明:对于轨道预报而言,GRGM660PRIM模型对20 km、50 km、100 km高度卫星的轨道寿命预报并未表现出巨大的优越性;对于大倾角低轨卫星的精密定轨,可以显著地提高其定轨精度。我国即将发射的嫦娥5号探测器需要在月球附近完成交会对接,采样返回地球,对定轨提出了更高要求,本文结论可以为我国嫦娥5号及之后的探月任务提供一定程度的参考。

参考文献
[1] Li Fei, Yan Jianguo. Principle and Method of Lunar Gravity Field Determination and Project on Self-determinational Lunar Grdavity Field[J]. Geomatics and Information Science of Wuhan University, 2007, 32(1):6-10(李斐, 鄢建国. 月球重力场的确定及构建我国自主月球重力场模型的方案研究[J]. 武汉大学学报·信息科学版, 2007,32(1):6-10)
[2] Yan Jianguo. Lunar Gravity Field Research and Lunar Satellite Precise Orbit Determination[D]. Wuhan:Wuhan University, 2007(鄢建国. 月球重力场研究及绕月卫星精密定轨[D]. 武汉:武汉大学, 2007)
[3] Zuber M T, Smith D E, Watkins M M, et al. Gravity Field of the Moon from the Gravity Recovery and Interior Laboratory (GRAIL) mission[J]. Science, 2013, 339(6 120):668-671
[4] Lemoine F G, Goossens S, Sabaka T J, et al. High-Degree Gravity Models from GRAIL Primary Mission Data[J]. Journal of Geophysical Research:Planets, 2013, 118(8):1 676-1 698
[5] Lemoine F G, Goossens S, Sabaka T J, et al. GRGM900C:A Degree 900 Lunar Gravity Model from GRAIL Primary and Extended Mission Data[J]. Geophysical Research Letters, 2014, 41(10):3 382-3 389
[6] Konopliv A S, Park R S, Yuan D N, et al. The JPL Lunar Gravity Field to Spherical Harmonic Degree 660 from the GRAIL Primary Mission[J]. Journal of Geophysical Research:Planets, 2013, 118(7):1 415-1 434
[7] Konopliv A S, Park R S, Yuan D N, et al. High-resolution Lunar Gravity Fields from the GRAIL Primary and Extended Missions[J]. Geophysical Research Letters, 2014, 41(5):1 452-1 458
[8] Konopliv A S. LP150Q.A01 (JGL150Q1.SHA) Lunar Prospector Derived Spherical Harmonic Model for the Moon Submitted to the Planetary Data System[EB/OL]. http://pds-eosciences.wustl.edu/lunar/lp-l-rss-5-gravity-v1/lp_1001/sha/jgl150 q1.lbl, 2000
[9] Konopliv A S, Asmar S W, Carranza E, et al. Recent Gravity Models as a Result of the Lunar Prospector Mission[J]. Icarus, 2001, 150(1):1-18
[10] Matsumoto K,Goossens S, Ishihara Y, et al. An Improved Lunar Gravity Field Model from SELENE and Historical Tracking Data:Revealing the Farside Gravity Features[J]. Journal of Geophysical Research:Planets, 2010, 115(E6):1 991-2 012
[11] Goossens S J, Matsumoto K, Kikuchi F, et al. Improved High-Resolution Lunar Gravity Field Model from SELENE and Historical Tracking Data[C]. AGU Fall Meeting, San Francisco, 2011
[12] Yan J G,Zhong Z, Li F, et al. Comparison Analyses on the 150×150 Lunar Gravity Field Models by Gravity/Topography Admittance, Correlation and Precision Orbit Determination[J]. Advances in Space Research, 2013, 52(3):512-520
[13] Yan Jianguo, Li Fei, Ping Jingsong, et al. Lunar Gravity Field Model CEGM-01 Based on Tracking Data of Chang'E-1[J]. Chinese Journal of Geophysics, 2011, 53(12):2843-2851(鄢建国, 李斐, 平劲松,等. 基于"嫦娥一号"跟踪数据的月球重力场模型CEGM-01[J]. 地球物理学报,2011,53(12):2 843-2 851)
[14] Yan J G,Goossens S, Matsumoto K, et al. CEGM02:An Improved Lunar Gravity Model using Chang'E-1 Orbital Tracking Data[J]. Planetary and Space Science, 2012, 62(1):1-9
[15] Wang Wei, Yan Jianguo, Shi Xian, et al. Numerical Analysis of Lunar Satellite Orbit Revolutionary[J]. Geomatics and Information Science of Wuhan University, 2007, 32(1):19-23(王威, 鄢建国, 史弦, 等.绕月飞行器近圆形轨道演化的数值分析[J].武汉大学学报·信息科学版, 2007, 32(1):19-23)
[16] Zhong Zhen, Li Fei, Yan Jianguo, et al. Comparison and Analysis on Main and Newly Lunar Gravity Field Models[J]. Geomatics and Information Science of Wuhan University, 2013, 38(4):390-393(钟振, 李斐, 鄢建国, 等.新近月球重力场模型的比较与分析[J].武汉大学学报·信息科学版, 2013, 38(4):390-393)
[17] Pavlis D E, Wimert J, McCarthy J. J. GEODYN II system description[R]. SGT Inc., Greenbelt, MD, 2013
[18] Mazarico E, Lemoine F G, Goossens S J, et al. Improved Precision Orbit Determination of Lunar Orbiters from the GRAIL-Derived Gravity Models[C]. AAS/AIAA Space Flight Mechanics Conference, Kauai, Hawaii, 2013
[19] Sjogren W L, Wimberly R N, Wollenhaupt W R. Lunar Gravity via the Apollo 15 and 16 Subsatellites[J]. Earth, Moon, and Planets, 1974, 9(1):115-128
[20] Li Fei, Yan Jianguo, Ping Jingsong, et al. Lunar Gravity Field Simulation Based on Big Inclination Orbiter[J]. Chinese Journal of Geophysics, 2011, 54(3):666-672(李斐, 鄢建国, 平劲松, 等.基于大倾角卫星轨道跟踪数据的月球重力场模型仿真解算[J].地球物理学报, 2011, 54(3):666-672)