环月飞行器精密定轨的模拟仿真

以中国正在实施的探月计划“嫦娥1号”工程为背景，分析了在中国联合S波段(USB)测控网和甚长基线射电干涉(VLBI)跟踪网的现有空间分布、观测精度水平下的环月飞行器精密定轨．采用的方法是模拟仿真计算，即首先模拟观测数据，然后在计入各误差源的影响后进行求解，并对解算结果进行比较．模拟仿真的工具是美国宇航局哥达德飞行中心的空间数据分析软件系统GEODYN．环月飞行的主要误差源是月球重力场，为此首先讨论了目前精度最高的月球重力场模型JGL165P1的(形式)误差．在模拟了测距、测速以及VLBI的时延、时延率数据后，计入月球重力场的误差进行精密轨道确定．定轨时采用了减缩动力学(reduced dynamic)方法，即选用合适的经验加速度参数吸收重力场误差对定轨的影响．结果表明对于一个不将月球重力场作为主要科学目标的探月计划(如“嫦娥1号”)，减缩动力学方法是一个简单、有效地提高环月飞行器定轨精度的方法．     

Orbit determination of satellite “Tance 1” with VLBI data

The satellite “Tance 1” of the “Double-Star Program” is the first truly scientific experimentation satellite of China. Its orbit is the farthest so far launched in China, with a geocentric apogee reaching 78 thousand kilometers. The tracking of “Tance 1” and of more distant space targets, such as the lunar exploration craft, can be realized with the VLBI technique of radio astronomy. In order to test and verify the role which the VLBI technique plays in the lunar exploration program of China, Shanghai Astronomical Observatory organized the only 3 tracking stations in China (located at Shanghai, Urumqi and Kunming), to carry out test tracking of “Tance 1,” and used the time delay data obtained to determine the orbit of “Tance 1” over a two-day period, so providing a preliminary assessment of the possibility of VLBI orbit determination. The fitting error of the orbit so obtained is about 5.5 m in the time delay and about 2 cm/s in the delay rate (this for checking only), much better than is provided by the preliminary orbit (used merely for ensuring tracking) in which the corresponding figures are around 2 km and 15 cm/s. Further, if the orbit is determined by using both the time delay and time delay rate data (with weights according to their internal accuracies), then the residuals are 5.5 m in the time delay and 2 cm/s in the delay rate. For an appreciation of the true accuracy of the VLBI orbit determination, we used simulation data (of the observed two-day VLBI data) and found the results depended greatly on the error in the dynamic model of the satellite which, however, is difficult to assess, while the formal residuals are of the order of 1 kin in the delay and of cm/s in the delay rate. The simulation computation also indicates that a joint determination using both VLBI and USB data will have an improved accuracy.     

利用VLBI数据确定"探测一号"卫星的轨道

双星计划的“探测一号”是中国首颗真正严格意义上的科学实验卫星,其运行轨道为中国迄今所发射的卫星中距地球最远,远地点地心距达7．8万公里．采用射电天文的VLBI技术可以对“探测一号”以及更远的深空目标,如探月飞行器实现跟踪．为了验证VLBI技术在我国探月计划中的作用,上海天文台组织了国内目前仅有的上海、乌鲁木齐和昆明3个台站对“探测一号”进行试跟踪,利用对“探测一号”约两天的VLBI观测数据,确定“探测一号”卫星的轨道,对VLBI的定轨能力做初步的探讨．按照测控部门提供的初轨 (其精度仅保证跟踪)推算的轨道与VLBI时延的拟合误差平均约2 km,时延率的拟合误差平均约15 cm／s．而利用VLBI数据定轨后的拟合程度相对于初轨有了很大的改善,结果表明,单独利用VLBI时延定轨,时延的拟合精度约5．5 m,作为外部检核的VLBI时延率的拟合精度在2 cm／s左右．单独利用VLBI时延率定轨,时延率的拟合精度约为1．3 cm／s,作为外部检核的VLBI时延的拟合精度约为29 m．而若将时延和时延率数据联合定轨,采用其内符精度加权,VLBI时延和时延率的残差分别为5．5 m和 2 cm／s．为了合理地评估VLBI定轨的真实精度,利用模拟数据进行误差协方差分析,结果表明VLBI定轨精度受动力学模型误差的影响较大,由于"探测一号”卫星的动力学模型难以精确确定,所以利用两天弧段的VLBI数据确定“探测一号”卫星轨道的位置误差为km量级,而速度误差可达cm／s量级．模拟计算还表明, VLBI和USB数据联合定轨可以大大提高定轨精度．     

月球交会对接VLBI差分时延研究

中国探月3期任务中,月球交会对接技术是任务成功的重要保障.利用嫦娥3号(CE03)绕月飞行的VLBI (Very Long Baseline Interferometry)时延数据,模拟仿真绕月交会对接过程中,同波束VLBI观测模式下,差分群时延的变化情况.仿真结果显示,在远程导引段,轨道器和上升组合体轨道距离保持100 km,持续半小时,差分群时延很好地反映了两者的轨道信息,可以用于定轨定位;自主控制段,上升组合体靠近轨道器,在轨道距离从5 km减小到20 m过程中,上升组合体加速追赶轨道器时,差分群时延快速趋近于0,上升组合体减速远离轨道器时,差分群时延绝对值快速变大.最后,利用嫦娥3号奔月段同时发射两个DOR (Differential One-Ranging)信号的VLBI时延数据,计算差分相时延,初步展示了月球交会对接过程中同波束VLBI差分相时延的误差情况.     

CVN硬盘系统和软件相关处理在e-VLBI试验中的应用

介绍了中国VLBI网（CVN）的e-VLBI技术研究进展．CVN包括上海佘山、乌鲁木齐南山2个固定观测站和云南昆明的流动站，以及上海天文台的2台站硬件相关处理机。2003年上海天文台自行研制了基于PC技术的VLBI数据记录、回放系统，命名为CVN硬盘系统，并成功将其安置于CVN观测站和处理机系统。硬件处理机经过改造后，已能处理来自硬盘和原有磁带系统的数据．从2003年至今，中国VLBI网采用该硬盘系统进行了多次VLBI观测和e-VLBI试验。在CVN硬盘系统基础上，软件相关处理技术的研究也得以开展。软件相关处理原型程序已经被用于台站条纹检测、卫星条纹搜索和数据处理中。该软件获得的计算结果被成功用于国内第一个3台站卫星VLBI的延迟和延迟率闭合试验，以及国内首次利用VLBI数据进行的卫星定轨试验。除此之外，该软件还用作硬件处理机的条纹引导器。为适应未来“嫦娥”月球探测工程，CVN将扩展成含有4个观测站和2个相关处理机（硬件、软件）的实时VLBI网。今后，e-VLBI将被应用于月球卫星导航以及测地和天体物理的VLBI观测。     

The Effect of GPS Ephemeris on the Accuracy of Precise Orbit Determination of LEO Satellite-borne GPS

The satellite-borne GPS receivers dedicated to precise orbit determination are now being carried by more and more low earth orbit (LEO) satellites and the satellite-borne GPS has become one of the main means for the precise orbit determination of low earth orbit satellites. The accuracy of satellite-borne GPS precise orbit determination depends on the accuracies of the GPS ephemeris and the clock error. Based on the orbit determination function of SHORDEIII zero-difference dynamics and using the observational data obtained by the GRACE satellites for the week from 2005 August 1 to 7 as an example, three versions of GPS ephemerides (igs, igr and igu) are used to carry out orbit determination under the same conditions and to estimate the effect of the GPS ephemeris accuracy on the accuracy of orbit determination of low earth orbit satellites. Our calculated results show that the two ephemerides, igs and igr, are equivalent to each other in orbit determination accuracy (about 9.5 cm), while igu is slightly less accurate, at about 10.5 cm. The effect produced by the data of the high frequency GPS satellite clock error on the accuracy of orbit determination is 1–6 cm.     

空间VLBI与天文地球动力学

介绍了空间VLBI的精密定轨及其在天文地球动力学应用研究中的最新进展．指出因为空间VLBI的时延和时延率观测量同时涉及到射电参考系、动力学参考系和地固参考系,所以特别适合于参考系的直接连接工作．对为评价参考系连接的精度而发展的协方差分析理论也作了介绍,还分析了将时延和时延率资料用于精密定轨时遇到的困难及其解决办法．     

Effects of the physical libration of the moon on lunar orbiters

Lunar physical libration, which is true oscillation of lunar equator in the space, alters the lunar gravitational field in the space coordinate system and affects the orbiting motion of lunar orbiters (hereafter called as lunar satellites) correspondingly. The effect is very similar to that of the precession and nutation on the earth satellites, and a similar treatment can be used. The variations in the gravitational force and in the orbit perturbation solution are clearly given in this paper together with numerical illustrations.     

Simulations of ultra-long wavelength interferometers in Earth orbit and on the lunar surface

We present simulations of interferometers in Earth orbit and on the lunar surface to guide the design and optimization of space-based ultra-long wavelength missions, such as those pioneered by China’s Chang’e Program. We choose parameters and present simulations using simulated data to identify inter-dependencies and constraints on science and engineering parameters. A regolith model is created for the lunar surface array simulation, and the results show that the lunar regolith will have an undesirable effect on the observations. We estimate data transmission requirements,calculate sensitivities for both cases, and discuss the trade-off between brightness temperature sensitivity and angular resolution for the Earth orbit array case.     

Lunar gravity derived from long-period satellite motion — A proposed method

A new method has been devised to determine the spherical harmonic coefficients of the lunar gravity field. This method consists of a two-step data reduction and estimation process. In the first step, a weighted least-squares empirical orbit determination scheme is applied to Doppler tracking data from lunar orbits to estimate longpperiod Kepler elements and rates. Each of the Kepler elements is represented by an independent function of time. The long-period perturbing effects of the Earth, Sun, and solar radiation are explicitly modeled in this scheme. Kepler element variations estimated by this empirical processor are then ascribed to the non-central lunar gravitation features. Doppler data are reduced in this manner for as many orbits as are available. In the second step, the Kepler element rates are used as input to a second least-squares processor that estimates lunar gravity coefficients using the long-period Lagrange perturbation equations.Pseudo Doppler data have been generated simulating two different lunar orbits. This analysis included the perturbing effects of the L1 lunar gravity field, the Earth, the Sun, and solar radiation pressure. Orbit determinations were performed on these data and long-period orbital elements obtained. The Kepler element rates from these solutions were used to recover L1 lunar gravity coefficients. Overall results of this controlled experiment show that lunar gravity coefficients can be accurately determined and that the method is dynamically consistent with long-period perturbation theory.     

日月轨道计算对TLE深空编目目标预报精度的影响

双行根数(Two Line Element, TLE)作为一类广泛使用的空间物体编目数据, 其预报精度和误差特性是TLE编目 在空间碎片研究中所要关注的问题之一. TLE编目需要配合SGP4/SDP4 (Simplified General Perturbations 4/Simplified Deep Space 4)模型进行轨道预报, 对深空物体来说, 主要考虑带谐项$J_2$、$J_3$、$J_4$摄动、 第三体日月摄动和特殊轨道共振问题修正等. 其中, SGP4/SDP4模型第三体摄动计算时, 对日月轨道近似采用了长期进动根数和 简单平运动的方式, 在外推10d时存在约2$^\circ$--3${^\circ     

The Study on the Group Delay Jumping of CE-4 VLBI System

The phase delay software is a key component of VLBI (Very Long Baseline Interferometry) orbit determination system in CE-4 mission, which is critical to deep space exploration. When processing X-band DOR (Differential One-way Ranging) signals of probes in the past, the group delay jumping occasionally occurred in the observed scan. This paper consists of two parts. First, according to the VLBI station signal reception and VLBI center data processing steps, we analyzed the frequency, phase, power, and self-correlation spectrum of the signal, and located the influencing factor at the abnormal analog signal of Kunming (KM) station. Second, the group delay jumping was revised by the post-correction method, which was based on the analysis of the abnormal unwrap phase. We used the corrected phase to get the new group delay, and used orbit determination software to verify the effectiveness of the post-correction method, which upgraded the phase delay processing software.     

嫦娥四号VLBI群时延跳变研究

相时延处理软件是嫦娥四号(CE-4)任务甚长基线干涉测量(Very Long Baseline Interferometry,VLBI)测轨分系统软件配置项之一,利用其处理着陆巡视器X波段DOR(Differential One-way Ranging)信号时,观测弧段(scan)内偶尔存在残余群时延跳变问题.论文主要分为两部分:一是根据VLBI信号接收和数据处理流程,从相关相位、频率、幅度和功率方面进行分析,最终将影响因素定位在昆明测站模拟信号异常;二是以scan内时延跳变影响因素分析为基础,通过修正异常基线相关相位,研究残余群时延跳变事后改正方法,并利用定轨软件验证其有效性,升级了相时延处理软件.     

A Proposed lunar orbiting gravity gradiometer experiment

Analysis of the gravity gradiometer developed by R. L. Forward and C. C. Bell at the Hughes Research Laboratories suggest than an accuracy, in the range 0.1 to 0.5 EU can be expected in a lunar orbiter application. This accuracy will allow gradient anomalies associated with mascons to be mapped with 1% accuracy and should reveal a great deal of new information about the lunar gravity field.The proposed experiment calls for putting such a gradiometer into a closely circular polar orbit at an average height of about 30 km above the lunar surface. This orbit allows the entire lunar surface to be covered in fourteen days, the gradiometer to be checked twice per revolution and results in successive passes above the lunar surface being spaced at about the resolution limit of about 30 km set both by the satellite altitude and instrumental integration time. Doppler tracking will be employed and the spacecraft will carry an electromagnetic altimeter. Gradient and altitude data from the far side of the Moon can be stored for replay when communication is re-established.     

Future Satellite Gravimetry for Geodesy

After GRACE and GOCE there will still be need and room for improvement of the knowledge (1) of the static gravity field at spatial scales between 40 km and 100 km, and (2) of the time varying gravity field at scales smaller than 500 km. This is shown based on the analysis of spectral signal power of various gravity field components and on the comparison with current knowledge and expected performance of GRACE and GOCE. Both, accuracy and resolution can be improved by future dedicated gravity satellite missions. For applications in geodesy, the spectral omission error due to the limited spatial resolution of a gravity satellite mission is a limiting factor. The recommended strategy is to extend as far as possible the spatial resolution of future missions, and to improve at the same time the modelling of the very small scale components using terrestrial gravity information and topographic models.We discuss the geodetic needs in improved gravity models in the areas of precise height systems, GNSS levelling, inertial navigation and precise orbit determination. Today global height systems with a 1 cm accuracy are required for sea level and ocean circulation studies. This can be achieved by a future satellite mission with higher spatial resolution in combination with improved local and regional gravity field modelling. A similar strategy could improve the very economic method of determination of physical heights by GNSS levelling from the decimeter to the centimeter level. In inertial vehicle navigation, in particular in sub-marine, aircraft and missile guidance, any improvement of global gravity field models would help to improve reliability and the radius of operation.     

Gravity field and rotation state of Mercury from the BepiColombo Radio Science Experiments

The ESA mission BepiColombo will include a Mercury Planetary Orbiter equipped with a full complement of instruments to perform Radio Science Experiments. Very precise range and range-rate tracking from Earth, on-board accelerometry, altimetry and accurate angular measurements with optical instruments will provide large data sets. From these it will be possible to study (1) the global gravity field of Mercury and its temporal variations due to tides, (2) the medium to short scale (down do 300400 km) gravity anomalies, (3) the rotation state of the planet, in particular the obliquity and the libration with respect to the 3/2 spin orbit resonance and (4) the orbit of the center of mass of the planet.With the global gravity field and the rotation state it is possible to tightly constrain the internal structure of the planet, in particular to determine whether the solid surface of the planet is decoupled from the inner core by some liquid layer, as postulated by dynamo theories of Mercury's magnetic field. With the gravity anomalies and altimetry it is possible to study the geophysics of the planet's crust, mantle and impact basins. With the orbit of the planet closest to the Sun it is possible to constrain relativistic theories of gravitation.The possibility of achieving these scientific goals has been tested with a full cycle numerical simulation of the Radio Science Experiments. It includes the generation of simulated tracking and accelerometer data, and the determination, by least squares fit, of a long list of variables including the initial conditions for each observed arc, calibration parameters, gravity field harmonic coefficients, and corrections to the orbit of Mercury. An error budget has been deduced both from the formal covariance matrices and from the actual difference between the nominal values used in the data simulation and the solution. Thus the most complete error budget contains the effect of systematic measurement errors and is by far more reliable than a formal one. For the rotation experiment an error budget has been computed on the basis of dedicated studies on each separate error source.The results of the full cycle simulation are positive, that is the experiments are feasible at the required level of accuracy. However, the extraction of the full accuracy results from the data will be by no means trivial, and there are a number of open problems, both in the data processing (e.g., the selection of the orbital arc length) and in the mission scheduling (e.g., the selection of the target areas for the rotation experiment).     

The Experiment and Analysis on VLBI Observations of Space Passive Hydrogen Maser

China plans to establish a lunar orbital VLBI (Very Long Baseline Interferometer) station around 2025, which will carry a space passive hydrogen maser as the time and frequency reference. Since it is the first time to use a space passive hydrogen maser for VLBI observation, its feasibility needs to be studied and verified. Therefore, we carried out VLBI observations using the space passive hydrogen maser as the frequency reference. In the experiment, the active hydrogen atomic clock and space passive hydrogen maser were used as the frequency standard, and the alternate VLBI observations of China’s Mars probe TW1 (Tianwen 1) were carried out using the 25 m radio telescope at Sheshan, Shanghai, and other VLBI stations. The results of data processing and analysis show that the standard deviation of VLBI residual group delay based on both active hydrogen atomic clock and space passive hydrogen maser are within 0.5 ns, which indicates that the space passive hydrogen maser can meet the accuracy requirements of VLBI measurement for deep space exploration, and verify its feasibility as the frequency standard of lunar orbital VLBI stations.     

利用国内VLBI网跟踪大椭圆轨道卫星

2004年7月,昆明VLBI站经过改造,由上海、乌鲁木齐和昆明站组成的中国VLBI网(CVN)采用统一的MARK4格式编制器和CVN硬盘记录系统,对大椭圆轨道卫星“探测1号”的2圈轨道的共同可视弧段进行了跟踪观测．软件相关处理程序已成功地用于检测卫星遥测信号的干涉条纹和数据相关处理．采用基于条纹幅度的加权最小二乘条纹拟合方法,获得了卫星VLBI观测量及其精度估计,完成了卫星VLBI观测量的3基线闭合误差检验．应用河外射电源校准方法和多频点相位校正信号提取方法,进行了台站钟差和仪器延迟等系统误差改正．经系统差改正后的卫星VLBI观测量序列已用于“探测1号”卫星的轨道确定．     

基于双频星载GPS数据的LEO卫星运动学定轨研究

运动学定轨是星载GPS特有的定轨方法,该方法不依赖于任何力学模型(地球重力场、大气阻力及太阳辐射压等),尤其适用于受大气阻力影响严重的低轨卫星定轨.基于双频星载GPS数据,研究了运动学定轨原理,讨论了数据预处理方法,建立了一套非差运动学定轨算法.并以GRACE (Gravity Recovery And Climate Experiment)-A、B卫星2008年2月实测数据作为试算验证了本研究方法的有效性和可靠性.GRACE 卫星实测数据计算结果表明:运动学定轨能达到5 cm精度(相对于SLR (Satellite Laser Ranging)),与动力学和简化动力学定轨精度相当.     

星载氢钟VLBI观测实验及分析

中国计划于2025年左右建立月球轨道VLBI (Very Long Baseline Interferometer)测站,将会搭载被动型星载氢钟作为时间频率标准.由于是首次在VLBI观测中使用星载氢钟,需要研究和验证其可行性.因此,利用星载氢钟作为频率基准开展了VLBI观测.实验时,分别使用主动型地面氢钟和被动型星载氢钟作为频率基准,利用上海天文台佘山25 m射电望远镜和其他测站对我国火星探测器天问一号进行了交替VLBI观测.数据处理分析结果表明,基于地面氢钟与星载氢钟的VLBI残余群时延标准差均在0.5 ns以内,表明星载氢钟可满足深空探测VLBI测定轨的精度要求,验证了其作为月球VLBI测站频率基准的可行性.