基于中继测量的环月探测器定轨能力分析

嫦娥四号月球探测拟首次实现月球背面的软着陆,测控与数传依赖地月L2平动点的中继卫星,并有望获取四程测量与星间测量数据。对基于中继测量的环月探测器测定轨能力进行了仿真分析,结果表明,中继卫星可较好地实现环月探测器连续跟踪;在定轨能力方面,中继卫星自身轨道精度是制约环月探测器定轨精度的重要因素,当跟踪弧段达到5h以上时,定轨精度趋于稳定,但轨道精度较中继卫星的轨道精度相差1个量级;对于星间链路测量,除中继卫星自身的轨道精度外,星钟的稳定性是制约定轨精度的另一个重要因素,如果辅助以每天1h的地基跟踪亦可实现优于百m的定轨精度。     

月球重力场模型精度比较研究

从阶方差、谱域分析以及大地水准面自由空气重力异常3个方面分析了月球重力场模型LP150Q,LP165P与SGM150各自的特点及差异,研究结果表明,SGM150重力场在中高阶项上得到了明显改善,在远月面可以发现明显的重力异常正负交替的现象。另外,对"嫦娥二号"变轨试验期间的测轨数据进行了轨道解算和预报,比较不同重力场模型的精度影响,发现基于SGM150模型的解算结果优于基于LP150Q和LP165P的,在定轨弧段内有10m量级的精度提升,可用于高精度的月球探测器精密定轨。     

嫦娥探测器轨道测定中的科学与技术问题

嫦娥探测器轨道测定技术既是成功完成工程任务的基础,也对月球科学研究起着重要作用。嫦娥探测器的轨道测定与地球航天器的轨道测定有诸多不同之处,本文介绍了绕月球探测与绕地球探测在测定轨技术方面的差异,并介绍了基于"嫦娥一号"、"嫦娥二号"月球探测器建立月面地形以及月球重力场模型的研究进展,最后分析了"嫦娥二号"平动点飞行试验以及小行星探测试验的测定轨能力。这些技术不仅在月球探测工程中起到关键作用,对我国后续深空探测也有很好的借鉴意义。     

三向测量模式在嫦娥三号探测器中的应用

以嫦娥三号探测器的三向测量跟踪数据为依据,利用最小二乘方法分析站间钟差,在5阶拟合计算中实现了ns级的拟合精度;然后使用探测器精密星历对三向测轨数据进行标校,经过修正后的测距系统差达到10m量级,拟合噪声水平优于1m;最后将三向测量数据应用于探测器的定轨、着陆点定位和动力下降段的弹道计算,其中环月段100km×15km轨道计算结果与精密定轨相比较偏差百m级,120km×70km轨道计算结果与精密定轨相比较偏差10m量级,月面定位与双程测距加时延定位结果相比较偏差10m量级。     

影响环月飞行器定轨精度的误差源分析

以我国正在实施的探月计划“嫦娥1号”工程为背案,在现有测控网分布、观测弧段以及尽可能接近真实情况的误差源等前提下,利用仿真模拟的方法对影响环月飞行器定轨精度的误差源进行了初步探讨和分析。重点考察了月球重力场误差、观测量精度、初始时刻的先验轨道误差以及观测资料类型等对环月飞行器定轨精度的影响。     

利用DORIS测轨系统实现高精度定轨

仿真研究DORIS测轨系统,重点考察了大气密度模型误差、测量精度、测轨网分布对定轨精度的影响。仿真结果表明,除了测轨网的地理分布,动力学模型中的大气模型误差对中低轨卫星精密定轨结果影响也较大。对ENVISAT卫星的DORIS实测数据进行了定轨分析,结果表明实测数据的定轨精度比仿真精度大约低1个量级不到。综合仿真结果和实测数据进行精度分析,推断对800km高度的太阳卫星轨道,采用8个DORIS信标站布设,24h定轨,定轨三维位置精度可以达到29cm(1σ),其中径向误差为3.4cm(1σ)。若采用30个DORIS信标站布设,定轨精度可提高30%。     

USB与VLBI联合确定“探测一号”卫星轨道

我国绕月探测工程“嫦娥一号”卫星将以统一S波段(USB)为主,辅以甚长基线干涉仪(VLB I)测轨分系统来完成测控任务。由于“探测一号”卫星轨道与“嫦娥一号”调相轨道段相似,有关单位于2005年3月17日—20日进行了USB和VLB I联合跟踪“探测一号”试验。通过对联合测轨数据的处理,研究了USB—VLB I联合定轨方法,分析了联合定轨和预报精度,得出了一些结论。     

环月快速交会调相策略设计与任务分析

针对未来载人登月任务中登月前环月轨道交会对接与组装问题,基于我国现有近地两天交会对接飞行方案,设计了环月轨道一天快速降轨交会任务的调相变轨方案,采用四脉冲修正特殊点变轨算法进行求解。分析了快速交会的调相终端控制精度、最优初始相位角范围等任务特性参数,给出了满足调相段终端控制精度所需要的定轨精度,分析了环月轨道倾角、调相段终端瞄准相对状态和调相时间对最优初始相位角范围和变轨总速度增量的影响规律。仿真结果表明,实施环月快速交会任务,要求追踪器与目标器的定轨精度均较高,但追踪器的最优初始相位角范围较大。     

月球重力场及低轨环月卫星轨道保持仿真分析

为了在满足精度要求的前提下节省月球重力场模型的计算时间,依据Kaula准则分析了目前国际上公认的最精确的两个重力场模型GLGM-2和LP165P,提出了在一定阶次截断重力场模型的问题.通过仿真不同阶次重力场模型作用下轨道高度为50 km的圆形极轨道环月卫星轨道特征的变化,验证了50 km以上高度卫星非球形摄动分析时可以将重力场模型截断至一定阶次的结论,并利用截断至70阶次的重力场模型仿真分析200 km圆轨道卫星一年内轨道下降程度.最后在仿真结果的基础上得到了200 km高度环月卫星需要每50天进行一次轨道保持控制的结论并完成一次轨道保持控制仿真.研究结论可以为我国低轨环月卫星轨道保持提供参考.     

金星探测器测定轨系统设计与试验验证

针对我国未来自主金星探测活动中的测量与导航实施问题,重点研究了地基测定轨系统的设计与实现,阐述了测距、测速和甚长基线干涉测量等能力,提出了基于我国金星探测测定轨系统设计方案,建立了高精度测量模型和定轨策略,并完成了软件实现。利用与欧空局联合开展的金星快车探测器跟踪试验,对该系统进行了初步验证,实测数据定轨结果与欧空局事后精密轨道的位置偏差优于485m,速度偏差优于5.5cm/s。试验结果验证了该系统的数据处理和轨道计算能力达到欧空局金星快车探测器同等水平,初步检验了该系统测量数据处理和轨道确定部分的正确性、有效性,可为后续我国自主金星探测提供测定轨技术支持。     

Phase Residual Estimations for PCVs of Spaceborne GPS Receiver Antenna and Their Impacts on Precise Orbit Determination of GRACE Satellites

In-flight phase center systematic errors of global positioning system(GPS) receiver antenna are the main restriction for improving the precision of precise orbit determination using dual-frequency GPS.Residual approach is one of the valid methods for in-flight calibration of GPS receiver antenna phase center variations(PCVs) from ground calibration.In this paper,followed by the correction model of spaceborne GPS receiver antenna phase center,ionosphere-free PCVs can be directly estimated by ionosphere-free carrier phase post-fit residuals of reduced dynamic orbit determination.By the data processing of gravity recovery and climate experiment(GRACE) satellites,the following conclusions are drawn.Firstly,the distributions of ionosphere-free carrier phase post-fit residuals from different periods have the similar systematic characteristics.Secondly,simulations show that the influence of phase residual estimations for ionosphere-free PCVs on orbit determination can reach the centimeter level.Finally,it is shown by in-flight data processing that phase residual estimations of current period could not only be used for the calibration for GPS receiver antenna phase center of foretime and current period,but also be used for the forecast of ionosphere-free PCVs in future period,and the accuracy of orbit determination can be well improved.     

国内外深空探测器精密定轨软件研究综述及WUDOGS简介

深空探测器精密定轨软件系统的研制在深空探测活动中是一个非常重要的环节,一直受到各大航天机构的重视。针对国内外深空探测器精密定轨软件平台的研究现状,重点介绍了具有代表性的美国JPL(Jet Propulsion Laboratory,喷气推进实验室)的DPTRAJ/ODP(Double Precision TRAJectory program/Orbit Determination Program,双精度轨道程序/定轨程序)和MONTE(Mission analysis,Operations,and Navigation Toolkit Environment,任务分析、操作和导航工具箱环境),GSFC(Goddard Space Flight Center,戈达德航天飞行中心)的GEODYN-II以及法国CNES(Centre National dEtudes Spatiales,国家空间研究中心)的GINS(Géodésie par Intégrations Numériques Simultanées,同步数值积分大地测量)软件系统,对这些软件的结构与功能进行了总结。之后对武汉大学自主研制的深空探测器精密定轨软件系统WUDOGS(Wuhan University Deep space Orbit determination and Gravity recovery System,武汉大学深空探测器精密定轨与重力场解算软件系统)的主要模块与功能进行了介绍,通过与GEODYN-II的交叉对比验证,表明:对于探测器的轨道预报,WUDOGS与GEODYN-II的1个月位置差异小于0.3mm,2d位置差值小于5×10~(-3) mm;双程测距、双程测速的理论计算值和GEODYN-II的差值RMS(Root Mean Square,均方根)分别在0.06mm,0.002mm/s的水平;WUDOGS目前已初步具备了月球和火星探测器精密定轨能力。最后对WUDOGS的下一步发展方向进行了展望。     

Mars Global Surveyor Observations of Solar Wind Magnetic Field Draping Around Mars

We examine the magnetic field in the martian magnetosheath due to solar wind draping. Mars Global Surveyor provided 3-D vector magnetic field measurements at a large range of altitudes, local times, and solar zenith angles as the spacecraft orbit evolved. We choose orbits with very clean signatures of draping to establish the nominal morphology of the magnetic field lines at local times of near-subsolar and near-terminator. Next, using a compilation of data from Mars Global Surveyor, we determine the average magnetic field morphology in the martian magnetosheath due to the solar wind interaction. The topology of the field is as expected from previous observations and predictions. The magnetic field magnitude peaks at low altitude and noon magnetic local time and decreases away from that point. The magnetic field has an inclination from the local horizontal of 5.6° on average in the dayside magnetosheath and 12.5° on the nightside. The inclination angle is closest to zero at noon magnetic local time and low altitude. It increases both upward and to later local times. The magnetic field in the induced magnetotail flares out from the Mars—Sun direction by 21°. Finally, we compare the observations to gasdynamic model predictions and find that the shocked solar wind flow in the martian magnetosheath can be treated as a gasdynamic flow with the magnetic pileup boundary as the inner boundary to the flow.     

Collocation of two GEO satellites and one inclined GSO satellite

Three collocation strategies are planned and analyzed for the cluster of two geostationary orbit (GEO) satellites and one inclined geosynchronous orbit (GSO) satellite in the same longitude control band of 116°E±0.05° . The longitudinal control bands are allocated for the two GEO satellites and one inclined GSO satellite with seven-day East/West station-keeping maneuver cycle. The latitudinal control bands are allocated for the two GEO satellites with fourteen-day North/South station-keeping maneuver cycle. One inclined GSO satellite is allowed for natural inclination drift. The coordinated eccentricity vector and inclination vector separation method is applied for the collocation, and the maneuver schedule is planned to minimize the operational load by avoiding simultaneous maneuvers. A total of six months of station-keeping maneuver simulations are performed for the three different strategies.     

The Lunar Gravity Ranging System for the Gravity Recovery and Interior Laboratory (GRAIL) Mission

The Lunar Gravity Ranging System (LGRS) flying on NASA’s Gravity Recovery and Interior Laboratory (GRAIL) mission measures fluctuations in the separation between the two GRAIL orbiters with sensitivity below 0.6 microns/Hz^1/2. GRAIL adapts the mission design and instrumentation from the Gravity Recovery and Climate Experiment (GRACE) to a make a precise gravitational map of Earth’s Moon. Phase measurements of Ka-band carrier signals transmitted between spacecraft with line-of-sight separations between 50 km to 225 km provide the primary observable. Measurements of time offsets between the orbiters, frequency calibrations, and precise orbit determination provided by the Global Positioning System on GRACE are replaced by an S-band time-transfer cross link and Deep Space Network Doppler tracking of an X-band radioscience beacon and the spacecraft telecommunications link. Lack of an atmosphere at the Moon allows use of a single-frequency link and elimination of the accelerometer compared to the GRACE instrumentation. This paper describes the implementation, testing and performance of the instrument complement flown on the two GRAIL orbiters.     

The FedSat Microsatellite Mission

An Australian research microsatellite, FedSat with a complement of four payloads was launched from Tanegashima, Japan on 14 December 2002 into a near-circular sun synchronous 10:30 LT polar orbit at an inclination of 98.7° and altitude 800 km. Scientific experiments include a triaxial fluxgate magnetometer with a frequency response up to 100 Hz and a GPS receiver to monitor total electron content (TEC) and provide a precise orbit determination. Communications experiments include a Ka-band transponder and a UHF packet data service. A high performance computer payload will test reconfigurable computing technology. This revised version was published online in August 2006 with corrections to the Cover Date.     

A geopause satellite system concept

The forthcoming 10 cm range tracking accuracy capability holds much promise in connection with a number of Earth and ocean dynamics investigations. These include a set of earthquake-related studies of fault motions and the Earth's tidal, polar and rotational motions, as well as studies of the gravity field and the sea surface topography which should furnish basic information about mass and heat flow in the oceans. The state of the orbit analysis art is presently at about the 10 m level, or about two orders of magnitude away from the 10 cm range accuracy capability expected in the next couple of years or so. The realization of a 10 cm orbit analysis capability awaits the solution of four kinds of problems, namely, those involving orbit determination and the lack of sufficient knowledge of tracking system biases, the gravity field, and tracking station locations. The Geopause satellite system concept offers promising approaches in connection with all of these areas. A typical Geopause satellite orbit has a 14 hour period, a mean height of about 4.6 Earth radii, and is nearly circular, polar, and normal to the ecliptic. At this height only a relatively few gravity terms have uncertainties corresponding to orbital perturbations above the decimeter level. The orbit s, in this sense, at the geopotential boundary, i.e., the geopause. The few remaining environmental quantities which may be significant can be determined by means of orbit analyses and accelerometers. The Geopause satellite system also provides the tracking geometery and coverage needed for determining the orbit, the tracking system biases and the station locations. Studies indicate that the Geopause satellite, tracked with a 2 cm ranging system from nine NASA affiliated sites, can yield decimeter station location accuracies. Five or more fundamental stations well distributed in longitude can view Geopause over the North Pole. This means not only that redundant data are available for determining tracking system biases, but also that both components of the polar motion can be observed frequently. When tracking Geopause, the NASA sites become a two-hemisphere configuration which is ideal for a number of Earth physics applications such as the observation of the polar motion with a time resolution of a fraction of a day. Geopause also provides the basic capability for satellite-to-satellite tracking of drag-free satellites for mapping the gravity field and altimeter satellites for surveying the sea surface topography. Geopause tracking a coplanar, drag-free satellite for two months to 0.03 mm per second accuracy can yield the geoid over the entire Earth to decimeter accuracy with 2.5° spatial resolution. Two Geopause satellites tracking a coplanar altimeter satellite can then yield ocean surface heights above the geoid with 7° spatial resolution every two weeks. These data will furnish basic boundary condition information about mass and heat flows in the oceans which are important in shaping weather and climate.     

自主机动条件下的实时定轨方法

中继卫星在跟踪自主机动用户目标时,由于机动轨道未知,需要利用中继卫星下传的星载GNSS(Global Navigations Satellite System,全球导航卫星系统)数据进行实时轨道确定与预报,为中继卫星跟踪提供实时的引导信息,以方便中继卫星快速捕获目标和连续稳定跟踪。针对该类用户目标的任务需求,讨论了基于星载GNSS数据自主机动条件下的实时定轨方法,建立了连续推力机动力学模型。以某一型号卫星的实测数据进行分析验证,并对轨道机动进行辨识,计算的机动加速度和机动时间与试验单位提供的结果一致。针对卫星不同机动情况,5min的观测数据定轨预报10min的弧段,最大位置误差小于8km,可以为中继卫星快速捕获提供高精度的引导信息。