胡光明，郭 睿，李曉杰，吳 杉

北京衛(wèi)星導(dǎo)航中心，北京， 100094

?

混合星座條件下的多星定軌及其精度

胡光明，郭 睿，李曉杰，吳 杉

北京衛(wèi)星導(dǎo)航中心，北京， 100094

為了實(shí)現(xiàn)與確保區(qū)域?qū)Ш椒?wù)能力，COMPASS設(shè)計(jì)采用GEO/IGSO/MEO混合星座，利用區(qū)域監(jiān)測(cè)網(wǎng)完成衛(wèi)星軌道測(cè)定與預(yù)報(bào)，并通過預(yù)報(bào)星歷確保導(dǎo)航服務(wù)。在區(qū)域布站條件下，難以實(shí)現(xiàn)對(duì)IGSO/MEO衛(wèi)星的全弧段跟蹤；混合星座條件下，不同衛(wèi)星的動(dòng)力學(xué)特性存在一定的差異，尤其是首次采用GEO衛(wèi)星，其靜地特性致使軌道與鐘差強(qiáng)相關(guān)，為動(dòng)力學(xué)模型選取帶來一定的困難。這些不利因素都給軌道確定帶來一定的挑戰(zhàn)。本文給出了區(qū)域監(jiān)測(cè)網(wǎng)條件下多星定軌的測(cè)量模型與動(dòng)力學(xué)模型，設(shè)計(jì)了軌道與鐘差參數(shù)的解算方案，系統(tǒng)地分析與論證了區(qū)域布站條件下的混合星座定軌精度。試驗(yàn)表明，衛(wèi)星三維位置精度優(yōu)于2m，基于激光數(shù)據(jù)檢驗(yàn)的GEO衛(wèi)星軌道預(yù)報(bào)2h視向精度優(yōu)于0.3m，IGSO衛(wèi)星軌道預(yù)報(bào)6h視向精度優(yōu)于0.4m。

COMPASS系統(tǒng)；多星定軌；衛(wèi)星激光測(cè)距數(shù)據(jù)；GEO；IGSO；MEO

1 引 言

COMPASS系統(tǒng)采用GEO/IGSO/MEO混合星座設(shè)計(jì)，以三個(gè)頻點(diǎn)碼分多址編碼方式，同時(shí)提供導(dǎo)航定位授時(shí)(Position-Navigation-Timing, PNT)等基本服務(wù)與差分服務(wù)[1]。此外，還通過星座中的GEO(Geostationary Earth Orbit)衛(wèi)星提供短報(bào)文通信和位置報(bào)告、 雙向定時(shí)等RDSS(Radio Determination Satellite Service)服務(wù)。

受限于地理分布有限的區(qū)域測(cè)軌監(jiān)測(cè)網(wǎng)，COMPASS系統(tǒng)無法實(shí)現(xiàn)對(duì)IGSO/MEO導(dǎo)航衛(wèi)星軌道的完全覆蓋。在區(qū)域監(jiān)測(cè)網(wǎng)內(nèi)，IGSO(Inclined GeoStationary Orbit)衛(wèi)星每天約有4個(gè)小時(shí)不可視，而MEO(Medium Earth Orbit)衛(wèi)星每天僅有30%的可視弧段，與全球布站的GPS(Global Positioning System)系統(tǒng)相比，有限的跟蹤弧段為高精度軌道測(cè)定與預(yù)報(bào)帶來諸多不利因素。

同時(shí)，COMPASS系統(tǒng)混合星座設(shè)計(jì)中，不同類型的衛(wèi)星高度、姿態(tài)控制方式、衛(wèi)星平臺(tái)和有效載荷均存在一定的差異，為動(dòng)力學(xué)模型的建模帶來一定的挑戰(zhàn)，特別是以太陽光壓模型為主的動(dòng)力學(xué)模型建模誤差，是制約混合星座定軌與預(yù)報(bào)精度的關(guān)鍵因素。

GEO衛(wèi)星的采用為COMPASS混合星座定軌帶來新的難題。影響GEO衛(wèi)星定軌精度的因素主要有兩方面：一是高軌衛(wèi)星的幾何跟蹤條件受到局部跟蹤網(wǎng)的限制，同時(shí)GEO衛(wèi)星的靜地特性使得跟蹤幾何幾乎不變，導(dǎo)致地面站對(duì)衛(wèi)星的動(dòng)力學(xué)約束強(qiáng)度非常弱；二是目前廣泛采用的GEO測(cè)軌技術(shù)具有較顯著的系統(tǒng)誤差，對(duì)定軌誤差影響較大，并且預(yù)報(bào)精度隨時(shí)間衰減較為明顯[2]。

常用的定軌模式包括單星定軌和多星定軌。其中，單星定軌主要利用地面站對(duì)單顆衛(wèi)星測(cè)軌數(shù)據(jù)獨(dú)立進(jìn)行軌道處理，多星定軌可以聯(lián)合地面站對(duì)所有在軌衛(wèi)星的測(cè)軌數(shù)據(jù)進(jìn)行綜合處理。前者由于測(cè)軌數(shù)據(jù)有限，同時(shí)，單星條件下不利于各類誤差的分離，因此定軌精度稍低，主要用于衛(wèi)星入網(wǎng)或機(jī)動(dòng)后快速恢復(fù)；后者能夠充分利用整網(wǎng)觀測(cè)數(shù)據(jù)，對(duì)衛(wèi)星鐘差和測(cè)站鐘差進(jìn)行建?；蚣s化處理，同時(shí)，對(duì)單個(gè)測(cè)站或單顆衛(wèi)星的公共誤差進(jìn)行建模，有利于各類誤差的分離，從而提高定軌精度，是目前最常用的定軌模式。

因此，在區(qū)域布站和混合星座的特殊條件下，如何實(shí)現(xiàn)COMPASS混合星座的多星整網(wǎng)軌道測(cè)定，是一個(gè)需要深入研究的問題。本文給出了多星定軌的數(shù)學(xué)模型，設(shè)計(jì)了多星整網(wǎng)定軌解算策略，利用實(shí)測(cè)數(shù)據(jù)對(duì)區(qū)域網(wǎng)絡(luò)跟蹤條件下的定軌精度進(jìn)行了分析論證，得到了一些有益的結(jié)論。

2 多星定軌原理

為了實(shí)現(xiàn)全星座條件下的多星定軌，設(shè)計(jì)采用中國(guó)區(qū)域內(nèi)地面站的偽距、相位無電離層組合觀測(cè)量進(jìn)行綜合處理，在軌道確定的同時(shí)，能夠?qū)崿F(xiàn)對(duì)衛(wèi)星鐘差和測(cè)站鐘差的有效估計(jì)。

2.1 測(cè)量模型及動(dòng)力學(xué)模型

測(cè)量型接收機(jī)不僅能夠進(jìn)行偽距觀測(cè)，而且可以獲得更高精度的相位數(shù)據(jù)。因此，多星定軌可以綜合采用偽距與相位數(shù)據(jù)進(jìn)行聯(lián)合定軌，在軌道確定的同時(shí)，實(shí)現(xiàn)對(duì)衛(wèi)星鐘差與測(cè)站鐘差的估計(jì)。多星定軌中一般采用批處理最小二乘方法進(jìn)行精密軌道確定。

多星定軌采用的數(shù)據(jù)為偽距、相位雙頻無電離層組合，其觀測(cè)方程可以表示為[3]：

(1)

(2)

衛(wèi)星精密定軌所采用的動(dòng)力學(xué)模型是影響定軌精度的重要因素。衛(wèi)星動(dòng)力學(xué)模型也稱攝動(dòng)力模型，用于描述衛(wèi)星受到的攝動(dòng)力對(duì)衛(wèi)星運(yùn)動(dòng)產(chǎn)生的加速度。衛(wèi)星受到的攝動(dòng)力可分為兩大類：保守力攝動(dòng)和非保守力攝動(dòng)。其中，保守力攝動(dòng)又分為N體攝動(dòng)、與地球引力位相關(guān)的攝動(dòng)、相對(duì)論效應(yīng)攝動(dòng)三種，與地球引力位相關(guān)的攝動(dòng)又可分為地球形狀攝動(dòng)、固體潮攝動(dòng)、海潮攝動(dòng)、大氣潮攝動(dòng)、地球自轉(zhuǎn)形變攝動(dòng)五類；非保守力攝動(dòng)包括太陽輻射壓攝動(dòng)和周期性RTN攝動(dòng)[4]。對(duì)于GEO與IGSO衛(wèi)星來說，其動(dòng)力學(xué)模型的最大誤差來自于太陽輻射壓模型的誤差。本文在進(jìn)行精密定軌研究時(shí)，采用的主要?jiǎng)恿W(xué)模型見表1。

表1 動(dòng)力學(xué)模型

攝動(dòng)力模型N體攝動(dòng)DE403太陽、行星攝動(dòng)重力場(chǎng)EIGEN-GRACE02S10×10太陽輻射壓GPST20地球固體潮IERS-Conventions1996經(jīng)驗(yàn)力軌道T/N方向調(diào)和函數(shù)

2.2 定軌解算策略

多星定軌試驗(yàn)采用國(guó)內(nèi)7個(gè)地面站的觀測(cè)數(shù)據(jù)，定軌弧長(zhǎng)為3天，數(shù)據(jù)采樣率為60秒。估計(jì)參數(shù)包括初始時(shí)刻衛(wèi)星位置、速度，太陽輻射壓攝動(dòng)系數(shù)，輻射壓y-bias，T/N方向經(jīng)驗(yàn)力參數(shù)，各監(jiān)測(cè)站大氣天頂延遲，相位模糊度以及所有歷元的鐘差參數(shù)。

IGS(International GNSS Service)計(jì)算表明，采用經(jīng)驗(yàn)太陽輻射壓模型可以模制輻射壓對(duì)衛(wèi)星軌道的主要攝動(dòng)影響，但仍存在一定殘余誤差，因此，精密定軌還需要估計(jì)經(jīng)驗(yàn)力參數(shù)對(duì)殘余誤差進(jìn)行模制。從軌道誤差看，受地面站的約束，軌道徑向誤差可以得到較好的控制，但軌道的沿跡方向和軌道面法向的誤差較難控制，因此，需要在軌道跡向和法向增加經(jīng)驗(yàn)力參數(shù)加以模制，吸收光壓模型的殘余誤差。

采用多星定軌還可以估計(jì)觀測(cè)模型參數(shù)，如接收機(jī)的大氣天頂延遲。由于常用的大氣折射誤差模型SAASTAMOINEN可以較好地模制干大氣折射誤差，而濕大氣部分受局部水汽壓影響較難精確模制，因此，多星定軌一般采用分段估計(jì)大氣天頂延遲方式來降低大氣模型誤差對(duì)精密定軌的影響。多星定軌估計(jì)的大氣參數(shù)為大氣模型在常溫常壓下計(jì)算修正值的修正因子。假設(shè)在3天的定軌弧段內(nèi)一個(gè)接收機(jī)估計(jì)M個(gè)天頂延遲參數(shù)，假定某時(shí)刻ti滿足ttrop(k)

(3)

式中，dt為歸一化的時(shí)差參數(shù)，trop_ori為采用SAASTAMOINEN模型在常溫常壓下計(jì)算的斜路徑大氣折射改正。

xtrop(k)(k=1,…M)為該接收機(jī)估計(jì)的第k個(gè)大氣參數(shù)。

ttrop(k)(k=1,…M)為第k個(gè)大氣參數(shù)對(duì)應(yīng)的時(shí)刻。

對(duì)大氣參數(shù)進(jìn)行線性化，第ti時(shí)刻對(duì)第k和k+1個(gè)大氣參數(shù)偏導(dǎo)數(shù)為：

(4)

將大氣參數(shù)與其他參數(shù)一起進(jìn)行最小二乘平差，可得到大氣估計(jì)值。

考慮到衛(wèi)星鐘差和測(cè)站鐘差可能存在不規(guī)則變化，多星定軌解算需對(duì)衛(wèi)星鐘差和測(cè)站鐘差進(jìn)行單歷元建模，即每個(gè)歷元設(shè)置一個(gè)衛(wèi)星或測(cè)站鐘差參數(shù)。對(duì)同一歷元，多星定軌只需要估計(jì)M+N(M為衛(wèi)星鐘差個(gè)數(shù)，N為測(cè)站鐘差個(gè)數(shù))個(gè)鐘差參數(shù) (實(shí)際由于秩虧只能獲得M+N-1個(gè)獨(dú)立的鐘差參數(shù)估計(jì))。而對(duì)整個(gè)定軌弧段，多星定軌需要估計(jì)Nepoch(M+N-1)個(gè)鐘差參數(shù), 其中Nepoch為數(shù)據(jù)歷元數(shù)。以60秒數(shù)據(jù)采樣計(jì)算，1天多星定軌需要估計(jì)多達(dá)1440×(M+N-1)個(gè)鐘差參數(shù)(未考慮不可視弧段)。

大量鐘差參數(shù)的估算需要耗費(fèi)大量的計(jì)算資源，為提高處理效率，多星定軌需采用鐘差約化算法，逐歷元約化鐘差參數(shù)。其基本思想是將待估參數(shù)分為全局參數(shù)和鐘差參數(shù)兩類，首先，通過逐歷元對(duì)鐘差參數(shù)進(jìn)行約化，消去衛(wèi)星鐘差和測(cè)站鐘差，僅解算所有動(dòng)力學(xué)參數(shù)、大氣延遲尺度參數(shù)和相位模糊度等全局參數(shù)；然后，將全局參數(shù)回代解算所有測(cè)站鐘差和衛(wèi)星鐘差，從而減小法方程大小，降低存儲(chǔ)空間，加快求解速度。解算鐘差時(shí)需設(shè)置基準(zhǔn)，原則上可以選取任一定軌測(cè)站的站鐘作為時(shí)間基準(zhǔn)。相關(guān)原理詳見文獻(xiàn)[5]。

3 多星定軌試驗(yàn)及精度分析

采用多星定軌方法對(duì)COMPASS系統(tǒng)混合星座進(jìn)行整網(wǎng)處理。采用2013年3月3日～6日共3天弧長(zhǎng)的數(shù)據(jù)進(jìn)行定軌試驗(yàn)，在軌衛(wèi)星包括5顆GEO、5顆IGSO以及4顆MEO，共14顆衛(wèi)星，地面站共7個(gè)站的偽距相位數(shù)據(jù)參與測(cè)定軌。

本文采用內(nèi)符合與外符合評(píng)估相結(jié)合的方法對(duì)多星定軌精度進(jìn)行綜合評(píng)估。其中，內(nèi)符合精度評(píng)估方法包括定軌殘差和重疊弧段法，外符合精度評(píng)估指利用衛(wèi)星激光測(cè)距數(shù)據(jù)(Satellite Laser Ranging, SLR)評(píng)估軌道視向精度的方法。

3.1 定軌殘差精度分析

定軌得到每顆衛(wèi)星與每個(gè)測(cè)站的偽距、相位殘差。例如GEO-1衛(wèi)星，其定軌殘差為7個(gè)測(cè)站對(duì)該星的定軌殘差統(tǒng)計(jì)值，參考站的定軌殘差為14顆衛(wèi)星對(duì)該站的定軌殘差統(tǒng)計(jì)值。各測(cè)站與各衛(wèi)星的偽距、相位無電離層組合殘差如圖1所示。3天弧段區(qū)域網(wǎng)定軌偽距殘差RMS(Root Mean Square)優(yōu)于50cm，相位殘差RMS優(yōu)于1cm，與區(qū)域GPS網(wǎng)定軌的殘差水平相當(dāng)。

為進(jìn)一步分析偽距相位殘差，表2給出了各測(cè)站所有14顆衛(wèi)星在3天定軌弧段內(nèi)偽距、相位無電離層組合的殘差RMS統(tǒng)計(jì)表。

圖1 多星定軌殘差統(tǒng)計(jì)結(jié)果

表2 衛(wèi)星多星定軌殘差統(tǒng)計(jì)表(單位:m)

衛(wèi)星類型衛(wèi)星號(hào)偽距殘差RMS殘差均值衛(wèi)星號(hào)相位殘差RMS殘差均值GEOGEO-10．2920GEO-20．2340GEO-30．3555GEO-40．3239GEO-50．32570．3062GEO-10．0088GEO-20．0077GEO-30．0083GEO-40．0082GEO-50．00950．0085IGSOIGSO-10．3118IGSO-20．2779IGSO-30．3216IGSO-40．2885IGSO-50．27330．2946IGSO-10．0096IGSO-20．0097IGSO-30．0110IGSO-40．0104IGSO-50．00910．0099MEOMEO-30．3960MEO-40．4664MEO-50．3907MEO-60．50690．4400MEO-30．0101MEO-40．0095MEO-50．0095MEO-60．00990．0098

圖2 1021站GOE-4、IGSO-1、MEO-3殘差圖

從表2可以看出，所有三類衛(wèi)星的相位殘差基本處于一個(gè)量級(jí)，高軌衛(wèi)星(GEO/IGSO)的偽距殘差也大體相近，而MEO的偽距殘差稍大，其主要原因在于，MEO觀測(cè)弧段較短，衛(wèi)星出入境前后數(shù)據(jù)質(zhì)量較差。圖2為1021站GEO-4、IGSO-1及MEO-3的殘差圖，從圖中可以看出，由于區(qū)域網(wǎng)無法覆蓋IGSO/MEO衛(wèi)星的全部軌道，所以該站對(duì)IGSO-1與MEO-3衛(wèi)星的可視弧段均存在空隙。按10°高度角截止，我國(guó)境內(nèi)站每24h對(duì)IGSO約有4h的不可視弧段，對(duì)MEO衛(wèi)星約有16h的不可視弧段。比較GEO與IGSO/MEO衛(wèi)星的殘差還可發(fā)現(xiàn)，在IGSO/MEO衛(wèi)星的出入境弧段，偽距的測(cè)量噪聲明顯增大。

3.2 重疊弧段精度分析

重疊弧段精度評(píng)估標(biāo)準(zhǔn)反映了兩次定軌的內(nèi)符度；同時(shí)，通過多天重疊弧段的比較，能夠反映定軌的穩(wěn)定性。評(píng)估方法為3天弧段重疊兩天，由重疊弧段在徑向、切向、法向上的誤差(ΔRi,ΔTi,ΔNi)，計(jì)算得到位置誤差：

(5)

以及URE(User Range Error)值(式6適用于GEO/IGSO衛(wèi)星，式7適用于MEO衛(wèi)星)：

(6)

(7)

統(tǒng)計(jì)得到各離散點(diǎn)處軌道之差的均方根誤差：

(8)

作為重疊弧段結(jié)果(其中Δσ可以代表R、T、N、Pos、URE中的任一變量含義)。

表3給出了將所有在軌衛(wèi)星于2013年3月2日～8日共7天的定軌結(jié)果進(jìn)行重疊弧段比較得到的R、T、N、Pos、URE精度。

表3 衛(wèi)星多星定軌重疊弧段統(tǒng)計(jì)表(單位:m)

衛(wèi)星類型衛(wèi)星號(hào)RTNPosUREGEOGEO-10．2671．2081．5031．9460．366GEO-20．1930．7621．0501．3110．257GEO-30．2191．3201．5142．0200．341GEO-40．6062．5481．8613．2130．733GEO-50．2970．8681．0911．4250．348IGSOIGSO-10．2621．5050．4911．6050．334IGSO-20．4342．0370．6412．1790．516IGSO-30．4881．1001．3251．7900．537IGSO-40．2401．1970．7141．4150．301IGSO-50．3901．8430．5861．9730．464MEOMEO-30．2401．1650．6461．3540．297MEO-40．2421．2060．6771．4040．302MEO-50．2690．9790．7311．2510．313MEO-60．2511．1360．6261．4190．324均值—0．3141．3480．9611．7360．388

可以發(fā)現(xiàn)，多星定軌算法穩(wěn)定，所有衛(wèi)星徑向重疊弧段精度(均值)優(yōu)于0.4m，位置精度優(yōu)于2m，URE精度優(yōu)于0.4m。還可以看出，各衛(wèi)星T、N方向的誤差明顯大于R方向誤差，其原因在于，我們所采用的是中國(guó)境內(nèi)站數(shù)據(jù)，所有跟蹤站均偏于一側(cè)所致。

3.3 衛(wèi)星激光測(cè)距檢驗(yàn)軌道視向精度分析

通過衛(wèi)星激光測(cè)距手段，獲得由地面站至衛(wèi)星間的幾何距離常被用作外符合手段對(duì)軌道視向精度進(jìn)行評(píng)估。該測(cè)量方法通過測(cè)量地面發(fā)射與接收激光信號(hào)之間的時(shí)差進(jìn)行距離測(cè)量，其測(cè)量手段決定了激光數(shù)據(jù)對(duì)衛(wèi)星與接收機(jī)鐘差不敏感。同時(shí)，由于激光測(cè)距的信號(hào)頻率較高，傳播路徑上的介質(zhì)時(shí)延量級(jí)比較小且可以較精確模制(厘米級(jí))，因此可以采用SLR進(jìn)行軌道精度評(píng)估。具體步驟是：對(duì)激光測(cè)距觀測(cè)值進(jìn)行測(cè)站和衛(wèi)星偏心改正、大氣折射誤差修正得到誤差修正后的觀測(cè)值O；利用多星定軌解算衛(wèi)星軌道和激光站坐標(biāo)可以計(jì)算給定時(shí)刻激光雙程距離理論值C，將計(jì)算O-C得到激光數(shù)據(jù)的殘差作為軌道的視向精度。

圖3給出了2013年3月4日利用北京的激光數(shù)據(jù)評(píng)估多星定軌的軌道視向精度殘差圖，其中3幅子圖分別表示GEO-1、IGSO-3、IGSO-5共3顆衛(wèi)星的軌道視向精度殘差圖。由于廣播星歷的更新周期是1h，考慮到IGSO衛(wèi)星每天存在約4h的出境弧長(zhǎng)，故本文采用激光數(shù)據(jù)評(píng)估GEO衛(wèi)星預(yù)報(bào)2h、IGSO衛(wèi)星預(yù)報(bào)6h的軌道視向精度。

圖3 SLR評(píng)估多星定軌不同衛(wèi)星視向精度殘差圖

可以看出，GEO衛(wèi)星預(yù)報(bào)2h、IGSO衛(wèi)星預(yù)報(bào)6h的SLR軌道視向精度均優(yōu)于0.4m，不同衛(wèi)星的激光測(cè)距殘差有所差別，各衛(wèi)星定軌弧段及預(yù)報(bào)弧段的SLR評(píng)估結(jié)果RMS均值詳見表4。

表4SLR評(píng)估的多星定軌視向精度(單位:m)

衛(wèi)星精度GEO-1IGSO-3IGSO-5定軌精度0．2270．1990．159預(yù)報(bào)精度0．2750．2070．320

4 小 結(jié)

本文深入研究了COMPASS系統(tǒng)導(dǎo)航衛(wèi)星多星軌道確定的特點(diǎn)與困難，詳細(xì)給出了混合星座條件下導(dǎo)航衛(wèi)星多星軌道確定的測(cè)量模型與動(dòng)力學(xué)模型，設(shè)計(jì)了軌道與鐘差參數(shù)的解算策略，并分別使用內(nèi)符合評(píng)估法(定軌殘差法、重迭弧段法)及外符合評(píng)估法(衛(wèi)星激光測(cè)距檢驗(yàn)軌道視向精度評(píng)估法)對(duì)多星定軌精度進(jìn)行了分析與評(píng)估。試驗(yàn)及分析結(jié)果表明：區(qū)域網(wǎng)混合星座條件下,多星定軌偽距殘差優(yōu)于0.5m，相位殘差優(yōu)于1cm；重疊弧段精度在徑向優(yōu)于0.4m，位置精度優(yōu)于2m；激光檢驗(yàn)的軌道視向精度優(yōu)于0.4m。

[1]北斗衛(wèi)星導(dǎo)航系統(tǒng)官網(wǎng).北斗衛(wèi)星導(dǎo)航系統(tǒng)簡(jiǎn)介[EB/OL]. [2010-01-15].http://www.beidou.gov.cn/xtjs.html.

[2]GuoRui,HuXiaogong,TangBo,etal.Preciseorbitdeterminationforthegeostationarysatellitewithmultipletrackingtechiques[J].Chin.Sci.Bull., 2010,55(6):428-434.

[3]周善石.基于區(qū)域監(jiān)測(cè)網(wǎng)的衛(wèi)星導(dǎo)航系統(tǒng)精密定軌方法研究[D]．上海:中國(guó)科學(xué)院上海天文臺(tái)，2011.

[4]DepartmentofDefense.GlobalPositioningSystemStandardPositioningServicePerformanceStandard[EB/OL].[ 2008-09-15].http://www.gps.gov/technical/ps/2008-SPS-performance-standard.pdf.

[5]周善石.區(qū)域衛(wèi)星導(dǎo)航系統(tǒng)GEO/IGSO混合星座的精密定軌和時(shí)間同步[J].中國(guó)科學(xué)G輯: 物理學(xué) 力學(xué) 天文學(xué)，2011，54(6): 1089-1097.

The Multi-satellite Precise Orbit Determination and It’s Accuracy for COMPASS Hybrid Constellation

Hu Guangming, Guo Rui, Li Xiaojie, Wu Shan

Beijing Satellite Navigation Center, Beijing 100094，China

In order to realize and ensure COMPASS regional navigation system services, COMPASS system adopts a hybrid constellation composed of GEO, IGSO and MEO, completes orbit determination and prediction by a regional tracking network and ensure navigation services through forecasting ephemeris. It is difficult for the regional tracking stations to make whole-arc tracking for IGSO/MEO. Meanwhile, for the hybrid constellation condition different type of satellites have differences in dynamic characteristics. Especially GEO satellites are first adopted for navigation purpose, whose stationary characteristics leads to the orbit and clock strongly correlated, which brings certain difficulty for selecting dynamic model parameters. This paper proposes a measurement model and a dynamic model for orbit determination in regional tracking network, designs a multi-satellite precise orbit determination (MPOD) for solving orbit and clock parameters simultaneously and analyzes the radial accuracy of hybrid constellation in regional stations. The experiments show that the 3-dimensional position accuracy is better than 2m, the radial accuracy of 2-hour orbital prediction for GEO satellites based on laser data validations is better than 0.3m and the accuracy of 6-hour orbital prediction is better than 0.4m.

COMPASS system; multi-satellite precise orbit determination(POD)； satellite laser ranging (SLR) data；geostationary earth orbit (GEO)；inclined geostationary orbit (IGSO)；medium earth orbit (MEO)

2015-03-16。

國(guó)家自然科學(xué)基金資助項(xiàng)目(41204022)， 上海市空間導(dǎo)航與定位技術(shù)重點(diǎn)實(shí)驗(yàn)室基金資助項(xiàng)目(12DZ2273300)。

胡光明(1974—)，男，工程師，主要從事衛(wèi)星導(dǎo)航數(shù)據(jù)處理和精密定軌方面的研究。

P228

A