唐歌实,曹建峰,韩松涛,胡松杰,任天鹏,陈略,孙靖,王美,李羿霏,李黎

(1.航天飞行动力学技术重点实验室,北京100094;2.北京航天飞行控制中心,北京100094)

基于CE-3的无线电测月研究

唐歌实1,2,曹建峰1,2,韩松涛1,2,胡松杰1,2,任天鹏1,2,陈略1,2,孙靖1,2,王美1,2,李羿霏1,2,李黎1,2

(1.航天飞行动力学技术重点实验室,北京100094;2.北京航天飞行控制中心,北京100094)

CE-3号探测器于2013年12月14日成功着陆于月表虹湾区,在预期1年的生存期内,着陆器能够相干转发下行X波段信号。基于这些信号,地面测控天线以及VLBI(very long baseline interferometry,甚长基线干涉测量)天线能够进行距离、速度、载波相位、VLBI时延以及时延率的测量。传统的激光测月技术只能提供距离测量信息,无线电测月技术能够同时提供更丰富的观测量类型,有可能获取优于激光测月技术的科学成果。为了分析LRM(lunar radio measurement,无线电测月)可能的贡献,MEKAS(Moon Earth kinematical analysis software)软件用于支撑无线电测月的数据分析处理,其能够模拟全部观测量类型,具备协方差分析以及CE-3着陆器定位、测站坐标解算等功能。模拟分析表明无线电测月技术在地月球科学方面具有广阔的应用前景。

嫦娥3号;玉兔号月球车;X波段转发器;无线电测月

0 引 言

激光测月(lunar laser ranging,LLR)是一种非常重要的空间大地测量技术,可确定月球星历、天平动、月球内部结构、地球自转参数以及月地动力学参数等重要参数[1-3]。LLR还是检验爱因斯坦广义相对论的有力工具[4-5]。但是,在所有的ILRS (International Laser Ranging Service,国际激光测距服务组织)观测站中,只有少数几个测站有能力开展激光测月观测,其中92%有效观测集中在美国的McDonald站和法国的Grasse站[6-8]。此外,LLR观测与ICRF(International Celestial Reference Frame,国际天球参考架)无关,很难把月球运动与准惯性天球参考架相连接。这些问题制约了LLR在地月科学研究中发挥更大作用,相比来说,利用CE-3号的全球联测及参考架连接优势可为地月科学研究提供新的机会。

CE-3探测器于2013年12月2日发射,经过12天的飞行,于14日成功着陆于月表虹湾区,设计寿命1年。着陆器安装了X波段转发器,能够以相干模式转发中国深空网(Chinese Deep Space Network,CDSN)上行测控信号。利用这些信号,中国深空网可以进行测速、测距,此外全球分布的IVS测站还可以开展VLBI观测。以上各类地月之间的测量信息统称为无线电测月(lunar radio measurement,LRM)。相对于激光测月而言,无线电测月具有如下优势:首先,大尺度的IVS天线分布提供了更好的PDOP因子;其次,VLBI是唯一的用于确定惯性天球参考架的技术,可以更好地连接月球轨道与天球参考架的关系。第三,可以同时获取到测距、多普勒测速、载波相位,时延、时延率数据。基于以上原因,无线电测月技术将有助于更好地了解月球行星科学,尤其是地月系统的动力学与运动学参考架。因而,对CE-3进行无线电测月技术研究具有重要科学价值。

1 系统组成

无线电测月技术体系包含两部分:一是地面跟踪测量系统,二是星上转发器。地面跟踪测量系统由CDSN、中国VLBI网(China VLB Net,CVN)、北京航天飞行控制中心(BACC)组成。深空网的佳木斯66 m天线、喀什35 m天线负责CE-3上行指令的发送以及遥测信号的接收,并兼顾测距、测速功能,此外深空站支持VSI、VSR格式的VLBI数据记录采集,干涉测量数据传输至北京中心相关处理中心,提取时延观测量。5个CVN网测站记录MARK5格式的VLBI数据,并传输至上海天文台相关处理中心用于提取干涉测量观测量。CDSN以及CVN网测站均配置有氢原子钟,GNSS接收机、水汽微波辐射计等辅助设备,用于时间、频率稳定性、传播介质时延等误差修正。

利用IVS全球布站的优势能够进一步提升位置精度因子。由于CE-3只在深空网可视跟踪弧段内具有下行信号,IVS分布在亚洲、欧洲、非洲以及大洋洲的测站具有更优的观测条件。此外由于CE-3下行信号相对于标校射电源功率更强,VGOS标准的12 m口径天线能够有效支持无线电测月观测。

为验证小口径天线对CE-3的观测性能,北京中心航天飞行动力学技术重点实验室改造一个12 m口径天线,其采用主动相位补偿技术向下变频等模块提供高稳定的频率基准。

CE-3着陆器、巡视器作为月球表面的信标机资源,具有如下特点:1)只在月昼期间工作;2)多个下行点频信号(DOR测音以及测距音);3)下行信号由北京中心遥控指令控制;4)巡视器遥测信号全程开启。

2 观测试验及精度评估

在CE-3任务期间,获取了多种测量模式的观测量:1)深空网测站获取的测距、测速观测量;2)由相关中心获取的不同基线时延、时延率;3)由相关中心获取的SBI双差分时延、时延率;4)由基带设备以及原始采集数据处理分析获取的载波相位。

为计算时延、时延率,测站原始数据需要传递至相关处理中心,目前软相关处理器以其灵活可扩展的优势被广泛应用[9]。北京中心研发了BSCS (BACC software correlator system)软相关处理系统,用于干涉测量的数据处理分析。软件采用C语言开发,运行于高性能集群系统。

图1 Delta-DOR测量精度分析(时延0.5 ns,时延率0.33 ps/s) Fig.1 Analysis of Delta-Dor measurement

下面几幅图显示了CE-3多种观测模式的精度分析情况。图1统计分析表明,Delta-DOR时延、时延率精度分别为0.5 ns以及0.33 ps/s;图2所示SBI测量结果精确反映了巡视器相对着陆器的运动状态,可分辨巡视器毫米级的位置变化。

图3统计分析表明,深空网载波相位测量精度为0.1 rad,12 m天线测量精度0.5 rad,在X波段等价于1～2 mm量级的距离测量精度。

图2 SBI测量结果分析Fig.2 Analysis of SBI measure result

图3 多普勒以及载波相位测量结果Fig.3 Measure result of Doppler and DOR

3 科学目标

通过CE-3着陆器的无线电测量技术,我们可以获取到高精度和更多采样点的地—月测量数据。基于该技术,通过CE-3项目有望进一步增强对地球与月球行星科学的认识。

为分析月球无线电测量数据,北京航天飞行控制中心开发了用于地月运动分析的软件系统(MEKAS)。该软件包括一系列基本模块,使用FORTRAN语言编码。MEKAS的4个主要基本功能包括:模拟观测数据、观测偏导数计算、待估参数的解算与协方差分析。软件可处理的观测数据类型包括双程/三程测距、多普勒测速、VLBI时延与时延率、激光测月数据等。可解算分析的参数包括着陆器的位置,站址坐标,地球自转参数,观测数据系统差,月球勒夫数等。参数的选取与设置非常灵活,MEKAS的输入设置为一个名为“mmekas.dat”的输入文件。在输入文件中,可以设置以下项:1)月球着陆器或是月球激光反射阵的位置与不确定度; 2)定位计算所使用的观测数据与权重设置;3)定位计算中需解算的参数设置;4)一个或多个观测数据文件。

利用MEKAS,我们进行了协方差分析。利用中国深空网与中国VLBI网对“嫦娥3号”着陆器测量数据进行定位分析,表1给出了可能的定位精度。

表1 着陆器定位精度Table 1 The positioning accuracyof the lunar lander

在定位计算中,测距数据的系统性偏差是影响制约定位精度的主要误差源,而且该偏差不易解算,因而在分析中考虑了高程方向约10 m的先验不确定度。如表1所示,如果不考虑站址误差,定位精度将随着测量数据的增加而提高,当跟踪数据达到5天后,结果表明可以实现米级的定位精度。如果考虑15 cm的站址误差,精度将会下降至10 m量级。

如果南半球的VLBI测站能够参与对“嫦娥3号”着陆器的跟踪测量,使用CDSN、CVN、IVS的测量数据进行了仿真分析,表2给出了IVS站参与跟踪的前提下可能的定位精度。

表2 着陆器定位精度(IVS站参与跟踪测量)Table 2 The positioning accuracy of the lunar lander(IVS involved in the measurement)

LLR激光反射阵在主轴坐标系下的坐标在行星历表建立的过程中同时解算。比较DE421与DE430,反射阵在x方向的坐标分量约有1 m的偏移。这表明着陆器的精度有可能达到与激光反射阵同等的精度水平,因而着陆器可以作为一个月球的基本控制点。

通常地球自转参数在IERS组织同组合估值的方式以时间序列的形式给出,并在公告中延迟30天发布。其发布的地球自转参数中极移的精度大约为0.1 mas,UT1精度为0.02 ms。精确确定了着陆器的位置,则能获取与IERS发布的相同精度的EOP参数。表3给出了考虑1 m不确定度的着陆器的位置误差解算的EOP参数精度。

在DE430的构建中,勒夫数h2是固定的,而k2设置为GRAIL确定的值。比较DE421与DE430的h2值,其偏差大约为0.01。使用长期的着陆器跟踪数据,可以进行月球勒夫数的解算。表4为基于协方差分析给出的勒夫数估计精度。当跟踪数据达到30天,h2解算精度可达0.005 3,l2可达0.002 3。

表3 利用着陆器测量解算EOP参数的精度Table 3 The uncertainty of EOP estimation with lunar lander

表4 月球勒夫数解算精度Table 4 The uncertainty of lunar love numbers estimation

4 结束语

CE-3着陆器以及巡视器定点于月球表面并发射下行无线电信号,信标定位精度提升至亚米级可以有效支持地月科学研究。MEKAS软件分析结果表明,在观测弧段足够长的情况下,CE-3信标定位精度可达到1 m。利用IVS测站天线以及12 m口径天线能够接收CE-3下行信号,相对LLR技术而言具有极高的性价比,在地月科学领域具有广阔的应用前景。

VLBI测量技术对射电源进行时延与时延率测量,可以用于高精度天球参考架的建立,EOP参数的解算。通过无线电测月技术不仅可以获取时延、时延率数据,还可获取测距、测速数据,如果结合激光测月计算,有望在地球与月球行星科学方面取得更多进展。包括月球历表,月球物理天平台,内部结构,不同参考架以及地球定向参数等等。

[1]Dickey J O.Lunar laser ranging:a continuing legacy of the Apollo program[J].Science,1994(265):482-490.

[2]Murphy T W.Laser ranging to the lost Lunokhod 1 reflector [J].Icarus,2011(2):1103-1108.

[3]Dehant V.Geodesy instrument package on the Moon for improving our knowledge of the Moon and the realization of reference frames[J].Planetary and Space Science,2012 (68):94-104.

[4]Adelberger E G.Eötvös experiments,lunar ranging and the strong equivalence principle[J].Nature,1990(347):261-263.

[5]Rambaux N,Williams J G.The Moon's physical librations and determination of their free modes[J].Celestial Mechanics and Dynamical Astronomy,2011(109):85-100.

[6]Bender P L.The Lunar laser ranging experiment[J]. Science,1973(182):229-238.

[7]Murphy T W.APOLLO:millimeter lunar laser ranging[R]. [S.l.]:Classical and Quantum Gravity(IOP Publishing),2012.

[8]Martini M.Moon LIGHT:a USA-Italy lunar laser ranging retroreflector array for the 21st century[J].Planetary and Space Science,2012(74):276-282.

[9]Deller A T,Tingay S J,Bailes M.DiFX:a software correlator for very long baseline interferometry using multiprocessor computing environments[J].Publications of the Astronomical Society of the Pacific,2007(119):318 336.

通信地址:北京市北清路26号院(100094)

电话:(010)66365934

E-mail:tanggeshi@bacc.org.cn

[责任编辑:宋宏]

Research on Lunar Radio Measurement on Chang'e-3

TAGN Geshi1,2,CAO Jianfeng1,2,HAN Songtao1,2,HU Songjie1,2,REN Tianpeng1,2, CHEN Lue1,2,SUN Jing1,2,WANG Mei1,2,LI Yifei1,2,LI Li1,2
(1.National Key Laboratory of Science and Technology on Aerospace Flight Dynamics,Beijing 100094,China; 2.Beijing Aerospace Control Center,Beijing 100094,China)

After the successful launch on December 2,2013,Chang'e-3 performed soft landing on lunar surface on December 14.Which is designed to stand in place for more than 1 year,and it can transmit X-band signal when the transponder is switched on.With their signals,ranging,Doppler,carrier phase,VLBI delay and delay rate can be acquired from ground TT&C antennas and VLBI antennas.These techniques can be called Lunar Radio Measurement(LRM)collectively.By LRM technique,more rich observations can be acquired nearly at the same time,so it could be more effective in contributing to space geodesy than LLR(Lunar Laser Ranging)which only provides ranging observation.To analyze the LRM technique,MEKAS(Moon Earth Kinematical Analysis Software)is developed,which can simulate all kinds of observations mentioned above,carry out covariance analysis and determine the parameters including CE-3 position on the lunar surface,ground site coordinates UT1,EOP and love number.Simulation results show that LRM technique has wide prospects in earth and lunar science.

Chang'e-3;Yutu Rover;X-band transponder;lunar radio measurement

V556.6

:A

:2095-7777(2014)03-0236-05

10.15982/j.issn.2095-7777.2014.03.013

唐歌实(1969—),男,研究员,主要研究方向:航天测控技术。

2014-07-26

2014-08-08

国家自然科学基金(61401014)