潘 超1,2? 劉慶會1? 鄭 鑫1,2 賀慶寶1,2 吳亞軍1,2

(1中國科學(xué)院上海天文臺上海200030) (2中國科學(xué)院大學(xué)北京100049)

基于差分相時延的月球車速率測定和動作分析?

潘 超1,2? 劉慶會1? 鄭 鑫1,2 賀慶寶1,2 吳亞軍1,2

(1中國科學(xué)院上海天文臺上海200030) (2中國科學(xué)院大學(xué)北京100049)

2013年12月14日嫦娥三號成功軟著陸,隨后開展兩器分離、兩器互拍等任務(wù).利用同波束VLBI(Very long baseline interferometry)技術(shù)同時觀測著陸器和玉兔月球車發(fā)射的信號,解算出的二者之間的差分相時延能夠反映月球車位置數(shù)厘米的微小變化.基于差分相時延監(jiān)測月球車動作的高靈敏度,利用月球車移動時的差分相時延數(shù)據(jù),求得月球車5次移動的速率,其平均值為0.056 m/s.利用差分相時延1階多項式擬合后殘差的頻譜,分析了月球車移動過程中的抖動情況及與月面地形的關(guān)系.

月球,技術(shù):干涉,技術(shù):其他諸多方面,方法:解析,方法:其他諸多方面

1 概述

嫦娥三號探測器于北京時間2013年12月14日21時11分在月球虹灣以東地區(qū)成功軟著陸,著陸點為西經(jīng)19.51?,北緯44.12?.著陸后開始執(zhí)行兩器分離和兩器互拍任務(wù).在兩器分離和兩器互拍過程中著陸器和月球車的距離很近,射電望遠鏡主波束可同時接收來自兩個探測器的信號,滿足同波束VLBI觀測條件.中國科學(xué)院VLBI網(wǎng)4架射電望遠鏡(上海天馬65 m(TM)、北京50 m(BJ)、昆明40 m(KM)和烏魯木齊25 m(UR))對整個過程采用同波束VLBI技術(shù)進行跟蹤觀測.利用觀測得到的著陸器和月球車的信號,提取積分時間0.983 04 s的相關(guān)相位,分別計算兩個探測器的殘余相時延,差分得到含整周模糊度的差分相時延.對兩個探測器的相時延進行差分,可以消除中性大氣、電離層和觀測裝置絕大部分的影響[1?5].由于著陸器在月面保持不動,差分相時延的變化主要反映月球車的動作及地球和月球的相對運動.

月球車在月面微小運動,差分相時延會隨之發(fā)生變化.月球車行走時,差分相時延不同基線相應(yīng)地快速變大或變小,月球車晃動時,差分相時延也來回晃動[6].基于差分相時延的高靈敏度,我們可以準(zhǔn)確地確定月球車動作開始和結(jié)束的時間點,進而根據(jù)兩個停泊點間的距離計算月球車移動的速率.本文利用月球車移動時的差分相時延數(shù)據(jù),求取了月球車玉兔在第1個月晝5次移動的速率,并利用差分相時延擬合殘差的頻譜,分析了月球車移動過程中的抖動情況及與月面地形的關(guān)系.本研究提供了一種基于地面測量進行月球車速率測定和動作分析的方法.

2 差分相時延監(jiān)視月球車動作靈敏度分析

圖1給出了VLBI差分相時延監(jiān)視月球車動作靈敏度的示意圖.如圖1(a)所示, cτ為VLBI幾何時延,cτ=B cosθ,因此,d(cτ)/dθ=?B sinθ.在此,θ約等于仰角.當(dāng)θ變化很小時,dθ≈Δθ≈?ΔSp/L,d(cτ)=Δ(cτ),因此,

式中B為基線長度,c為光速,τ為時延,L為地月間距離,ΔSp是在天球切面內(nèi)監(jiān)視月球車動作的靈敏度[7].

圖1 VLBI時延在天球切面內(nèi)(a)和月球表面(b)的監(jiān)視靈敏度Fig.1 The monitoring sensitivity of VLBI delay in the plane of the sky(a),and the lunar surface(b)

月球表面在較小的范圍內(nèi)可近似用切平面來表示.由于天球切面內(nèi)的監(jiān)視靈敏度ΔSp是月球車實際移動距離ΔSs在天球切面上的分量,所以當(dāng)月球車在經(jīng)度λ、緯度?的月球表面移動時,監(jiān)視靈敏度可按下述方法計算.圖1(b)給出了緯度?對靈敏度的影響(經(jīng)度λ的影響類似).當(dāng)月球車在月球表面移動ΔSs的距離,其在天球切面內(nèi)的分量ΔSp為:ΔSp=ΔSscos?.因此,當(dāng)月球車處于經(jīng)度λ、緯度?的月球表面時,其靈敏度為:

以較短的北京-昆明基線(B=2160 km)為例,差分相時延隨機誤差0.5 ps相當(dāng)于0.15 mm[6],假設(shè)仰角θ為30?,地月間距離為380 000 km.嫦娥三號探測器著陸點位于經(jīng)度λ=?19.5124?=?0.34 rad,緯度?=44.1206?=0.77 rad,月球車動作的監(jiān)視靈敏度估算如下:

3 月球車移動時段的截取和速率的測定

根據(jù)差分相時延隨時間的變化趨勢,我們可以監(jiān)測月球車在月面移動的情況.由于差分相時延具有超高靈敏度,我們可以準(zhǔn)確地確定月球車移動的開始時間和結(jié)束時間,由此得出月球車的移動時間.然后由月球車在各停泊點的坐標(biāo)可以求出每次移動的距離,進而求得每次移動的平均速率.表1給出了月球車在各停泊點A、B、C、D、E、S1的相對于著陸器的坐標(biāo),其結(jié)果是利用2～3 h的差分相時延數(shù)據(jù)對月球車進行相對定位得出的,誤差為1 m[7].圖2給出了月球車在各停泊點之間移動時的北京-昆明基線的差分相時延(DPD).

圖2 月球車5次移動時的差分相時延Fig.2 The di ff erential phase delay during rover’s 5 movements

圖2(a)為12月15日月球車由A點移動到B點的差分相時延.由表1所示A點和B點的坐標(biāo)可以得到AB間距離為8.24 m.月球車移動了兩次,由A點移動到B點共用時142 s,平均移動速率為0.058 m/s.12月20日月球車有2次移動,分別為由B點移動到C點,和由C點移動到D點.其中B到C移動分2步完成,如圖2(b)所示.C到D移動1次完成,如圖2(c)所示.B到C和C到D的距離、移動時間和平均速率由表2給出.12月21日月球車由D點移動到E點,如圖2(d)所示,本次移動分2步完成.23日,月球車由E點移動至S1點,如圖2(e)所示.這3點間的距離、移動時間和平均速率由表2給出.這6點間移動的總平均速率為0.056 m/s.

表1 月球車在停泊點的坐標(biāo)(著陸器位置為坐標(biāo)原點)Table 1 The coordinates at the parking point of rover(the lander position as the origin of coordinates)

表2 6點間移動的距離、時間和速率Table 2 The moving distance,time,and speed between each two of six points

文獻[7]通過對比基于差分相時延測定的月球車上高增益天線與低增益天線之間的相對距離的結(jié)果與實際的設(shè)計距離的差異(約0.16 m)及其他分析,得出了基于差分相時延測定的月球車停泊點的定位誤差為1 m的結(jié)論.基于差分相時延的定位結(jié)果對應(yīng)著陸器和月球車天線相位中心間的相對位置,而用星載照相機視覺定位結(jié)果對應(yīng)于著陸器和月球車的機械中心參考點.因此兩種定位結(jié)果在參考系定義上存在1～2 m的差異[8].這個差異對停泊點間的各段移動速率會有一定影響,但對各段總的平均速率影響不大.根據(jù)文獻[9]視覺定位的結(jié)果,利用本文中求得的各段移動時間,計算得到的總的平均速率為0.053 m/s,與本文求得的平均速率0.056 m/s的差異為0.003 m/s.綜合上述討論,我們認(rèn)為,由于月球車每一次移動時間不少于100 s,各停泊點的定位精度為1 m,因此,基于差分相時延求得的總的平均速率的誤差小于0.01 m/s.

4 基于擬合殘差譜的月球車抖動分析

由于月球表面凹凸不平和月球車自身機械構(gòu)造等原因,月球車在月面移動或轉(zhuǎn)彎過程中,不可避免地會產(chǎn)生晃動.差分相時延具有極高的靈敏度,月球車數(shù)厘米微小的晃動也會反映到差分相時延中.月球車移動1次的時間為100多秒,而月球和地球相對運動引起的差分相時延變化的周期約為24 h.由圖2也可以看出,每次移動引起的差分相時延變化近似直線,故我們可對差分相時延進行1階多項式擬合,利用擬合殘差以及擬合殘差的頻譜來研究月球車在移動過程中的晃動情況.

圖3(a)為月球車從A點到B點第1次移動時的差分相時延及擬合直線,圖3(b)為擬合殘差.圖3(c)和(d)分別為月球車從A點到B點第2次移動時的差分相時延、擬合直線和擬合殘差.4幅圖均為BJ-KM基線數(shù)據(jù),其余基線具有類似的情況.從圖3(a)和(c)中可以看出,月球車從A點到B點的2次移動,第2次移動時的差分相時延比第1次的更加貼合擬合直線.從圖3(b)和(d)中可以看出,第1次移動時比第2次移動時的晃動頻率更高,幅度更大.

圖3 月球車由A點移動到B點的差分相時延、擬合直線及擬合殘差Fig.3 The di ff erential phase delay, fi tting lines,and fi tting residuals of the rover’s movement from point A to point B

圖4(a)和(b)分別給出了月球車由A點到B點第1次移動時的BJ-UR和BJ-TM基線的差分相時延擬合殘差的頻譜.由圖4(a)和(b)可以看出,在0.13 Hz附近存在相對幅值較大的譜線(頻譜圖中的低頻成分比較突出,這主要是由于月球車實際行走路線不完全是直線造成的,不反映月球車實際晃動情況),表明月球車存在頻率約為0.13 Hz的晃動.月球車晃動的原因可能是A點附近小石塊較多所致,如圖4(e)所示.圖4(c)和(d)給出了月球車由A點到B點第2次移動時的差分相時延擬合殘差的頻譜.由圖4(c)和(d)可以看出頻譜譜線分布相對均勻,沒有明顯突出的頻率成分.這說明月球車由A點到B點第2次移動過程中沒有明顯的周期性晃動,由圖4(e)也可看出后半段月面地貌相對平整.

圖5(a)和(b)分別給出了月球車由B點到C點第1次移動時的BJ-KM和BJ-UR基線的差分相時延擬合殘差的頻譜.由這兩幅頻譜圖可以看出在0.14 Hz附近的頻譜成分是相對突出的.圖5(c)和(d)給出的是月球車由B點到C點第2次移動時的差分相時延擬合殘差的頻譜.這兩幅頻譜圖中0.4 Hz附近的譜線幅值相對較高.月球車由B點到C點的這2次移動都存在不同頻率的晃動.圖5(e)給出了著陸器拍攝的月球車由B點到C點的移動路線月面地貌,可以看出這片區(qū)域中小石塊、小撞擊坑遍布,月面地形起伏較大.

圖6(a)和(b)給出了月球車由C點到D點移動時的BJ-KM和BJ-UR基線的差分相時延擬合殘差的頻譜.從兩幅頻譜圖中可以看出,兩條基線差分相時延擬合殘差頻譜除了低頻外,不存在其他突出的頻譜成分.月球車在由C點移動到D點的過程中沒有明顯的周期性晃動,顯示出C點到D點之間的月面地貌比較平整,如圖6(c)所示.

圖4 月球車由A點到B點移動時的差分相時延擬合殘差的頻譜以及著陸器拍攝的月球車移動路徑照片F(xiàn)ig.4 The spectrum of the fi tting residuals of the di ff erential phase delay when the rover moved from point A to point B,and the photo of the rover moving path taken by the lander

圖5 月球車由B點到C點移動時的差分相時延擬合殘差的頻譜以及著陸器拍攝的月球車移動路徑照片F(xiàn)ig.5 The spectrum of the fi tting residuals of the di ff erential phase delay when the rover moved from point B to point C,and the photo of the rover moving path taken by the lander

圖6 月球車由C點到D點移動時的差分相時延擬合殘差的頻譜以及著陸器拍攝的月球車移動路徑照片F(xiàn)ig.6 The spectrum of the fi tting residuals of the di ff erential phase delay when the rover moved from point C to point D,and the photo of the rover moving path taken by the lander

圖7(a)和(b)給出了月球車由D點到E點第1次移動時的BJ-UR和KM-UR基線的差分相時延擬合殘差的頻譜.圖7(c)和(d)給出了月球車由D點到E點第2次移動時的BJUR和KM-UR基線的差分相時延擬合殘差的頻譜.2次移動的頻譜圖中兩條基線的譜線分布都比較均勻,不存在明顯突出的頻譜成分.由D點到E點,月球車移動過程中無明顯周期性晃動.觀察圖7(e)可以看出,這片區(qū)域沒有明顯的裸露于月面地表的石塊,月面相對平整.

圖7 月球車由D點到E點移動時的差分相時延擬合殘差的頻譜以及著陸器拍攝的月球車移動路徑照片F(xiàn)ig.7 The spectrum of the fi tting residuals of the di ff erential phase delay when the rover moved from point D to point E,and the photo of the rover moving path taken by the lander

23日月球車由E點移動到S1點,這次移動分3步完成.第2次和第3次移動的時間較短,采樣點數(shù)太少,頻譜分辨率低,不作分析.由E點到S1點第1次移動的差分相時延擬合殘差的頻譜中無明顯突出的譜線,月球車沒有明顯的周期性抖動.根據(jù)前面的分析可以估計,這一段路線是比較平整的.

5 討論

差分相時延的快速變化反映了著陸器和月球車上的兩天線間的相對位置的變化.著陸器上的天線位置保持不變,而月球車上的天線會隨月球車本身的移動或抖動發(fā)生相應(yīng)的移動或抖動.當(dāng)月球車兩側(cè)車輪同時通過月球坑或小石塊,天線相位中心和車輪的位置變化幅度幾乎相同,約為月球坑或小石塊的尺寸;而當(dāng)月球車僅有單側(cè)車輪通過月球坑或小石塊時,天線相位中心較該側(cè)車輪的位置變化幅度會更大一些.月球車長度為150 cm,月球車前后車輪間間距約為75 cm.當(dāng)月球車以平均速率5.6 cm/s行駛,單側(cè)3個輪子依次越過一塊小石塊時,會引起3次震動,頻率為(5.6 cm/s)/(75 cm)=0.075 Hz.當(dāng)小石塊、月球坑分布密集且間距小于車輪間距75 cm時,震動頻率會增大(存在達到0.13 Hz、0.14 Hz和0.4 Hz的可能).14日月球車和著陸器分離時,有4次平移,平移距離分別為140 cm、10 cm、5 cm和5 cm,前面兩次大的平移時,差分相時延有明顯的變化,平移5 cm時,差分相時延反映不明顯,如圖8所示.這說明差分相時延監(jiān)視月球車動作的靈敏度在5～10 cm之間[6],與本文利用(2)式估算的結(jié)果基本符合.月球車有6個直徑30 cm、寬度15 cm的車輪,具備20?爬坡,20 cm越障能力,可以翻越使差分相時延產(chǎn)生明顯抖動的小石塊或者坑.在嫦娥三號中,對差分相時延擬合后求得殘差消除了大尺度月面高程變化的影響,擬合殘差主要反映10 cm左右的月面地貌的凹凸變化,如小的月球坑或石塊.擬合殘差的頻譜幅度越大,頻率成分越多,則說明月面地貌凹凸變化較多、較大,月面越不平整.文獻[9]給出了地形高程的相對變化曲線,僅從各點間曲線的小尺度變化分析可以看出,點N0102與點N0103之間曲線的前半部分明顯比后半部分變化劇烈;點N0103和點N0104之間曲線的凹凸變化較為劇烈;點N0104、點N0105和點N0106之間的兩段曲線凹凸變化的尺度較小,小于10 cm;點N0106和點N0107間的曲線過渡較為平滑,幾乎無小尺度的凹凸變化.點N0102、N0103、N0104、N0105、N0106、N0107分別對應(yīng)本文中點A、B、C、D、E、S1.上述結(jié)果與本文頻譜分析得到的結(jié)果有一定的相似性.

圖8 月球車14日4次平移差分相時延(KM-UR基線)Fig.8 The di ff erential phase delay when the rover moved 4 times on December 14(KM-UR baseline)

6 總結(jié)

本文基于差分相時延的高靈敏度,對月球車的移動動作進行了分析.通過分析差分相時延的變化,可以準(zhǔn)確確定月球車移動動作的時間,進而得到月球車的移動速率,其平均值為0.056 m/s.差分相時延擬合殘差反映的是月球車移動過程中的抖動,當(dāng)擬合殘差的頻譜中存在比較突出的頻譜成分,則說明月球車的抖動有周期性.分析結(jié)果表明,月球車由B點移動至C點時,抖動較大,對應(yīng)的月面地形起伏也較大.

[1]劉慶會,陳明,熊蔚明,等.中國科學(xué)G輯,2010,40:253

[2]Liu Q H,Kikuchi F,Tsuruta S,et al.ITAP,2007,55:1466

[3]Liu Q H,Kikuchi F.AdSpR,2007,40:51

[4]Liu Q H,Kikuchi F,Goossens S,et al.JGSJ,2009,55:243

[5]Chen M,Liu Q H.SCSMP,2011,54:2284

[6]鄭鑫,劉慶會,吳亞軍,等.中國科學(xué)G輯,2014,44:872

[7]Liu Q H,Zheng X,Huang Y,et al.RaSc,2014,49:1080

[8]黃勇,昌勝騏,李培佳,等.科學(xué)通報,2014,59:2268

[9]Liu Z Q,Di K C,Peng M,et al.SCSMP,2015,58:1

Speed Measurement and Motion Analysis of Chang’E-3 Rover Based on Di ff erential Phase Delay

PAN Chao1,2LIU Qing-hui1ZHENG Xin1,2HE Qing-bao1,2WU Ya-jun1,2
(1 Shanghai Astronomical Observatory,Chinese Academy of Sciences,Shanghai 200030) (2 University of Chinese Academy of Sciences,Beijing 100049)

On 2013 December 14,the Chang’E-3 made a successful soft landing on the lunar surface,and then carried out the tasks of separating the lander and the rover, and taking the photos of each other.With the same beam VLBI(Very long baseline interferometry)technique to observe the signals transmitted by the lander and the rover simultaneously,the di ff erential phase delay between them is calculated,which can re fl ect a minor change of the rover’s position on a scale of a few centimeters.Based on the high sensitivity of di ff erential phase delay,the rover’s speeds during 5 movements are obtained with an average of 0.056 m/s.The relationship between the rover’s shake in moving process,and lunar terrain is analyzed by using the spectrum of the residual of the di ff erential phase delay after the fi rst-order polynomial fi tting.

Moon,techniques:interferometric,techniques:miscellaneous,methods: analytical,methods:miscellaneous

P164;

:A

2014-12-22收到原稿,2015-02-04收到修改稿

?國 家 自 然科 學(xué) 基 金 項 目(11273049,11473059)、 上 海 市 導(dǎo) 航 和 定 位重 點 實 驗室(3912DZ22733001)資助

?panchao@shao.ac.cn

?liuqinghui@shao.ac.cn

10.15940/j.cnki.0001-5245.2015.04.005