郭 睿，陳金平，朱陵鳳，唐桂芬，吳曉莉

1. 北京衛(wèi)星導航中心, 北京 100094； 2. 大地測量與地球動力學國家重點實驗室，湖北 武漢 430077

?

北斗衛(wèi)星超短弧運動學定軌方法優(yōu)化與試驗分析

郭 睿1,2，陳金平1，朱陵鳳1，唐桂芬1，吳曉莉1

1. 北京衛(wèi)星導航中心, 北京 100094； 2. 大地測量與地球動力學國家重點實驗室，湖北 武漢 430077

針對北斗衛(wèi)星姿軌控后的軌道快速確定難題，系統(tǒng)地研究了基于多項式擬合和基于星歷擬合兩種運動學定軌方法，推導建立了相應的運動學定軌模型。同時針對接收機系統(tǒng)差和頑固多徑問題，利用基于并置比對的接收機系統(tǒng)差解算方法和CNMC的多徑削弱方法，實現(xiàn)了超短弧跟蹤條件下接收機數(shù)據(jù)質(zhì)量的有效控制。利用北斗GEO/IGSO/MEO衛(wèi)星的實測偽距數(shù)據(jù)進行了試驗驗證，結果表明在10 min超短弧跟蹤條件下，GEO、IGSO和MEO衛(wèi)星的運動學定軌位置精度分別為3.27 m、8.19 m和5.90 m，實現(xiàn)了超短弧跟蹤條件下的北斗衛(wèi)星快速定軌，滿足了衛(wèi)星機動期間的北斗RDSS服務對軌道精度的需求，為北斗RDSS服務走向全球提供了技術支撐。

北斗衛(wèi)星導航系統(tǒng)；精密定軌；軌道機動；運動學法定軌；多徑誤差

北斗衛(wèi)星導航系統(tǒng)(BDS)作為全球四大衛(wèi)星導航系統(tǒng)之一，于2012年12月27日正式開始運行，主要為中國及周邊地區(qū)提供定位、授時和短報文服務[1]，首次同時實現(xiàn)了無線電定位業(yè)務(radio determination satellite service, RDSS)和無線電導航業(yè)務(radio navigation satellite service, RNSS)服務。RDSS服務是BDS系統(tǒng)的重要組成和鮮明特色，主要依托在軌5顆GEO(geostationary earth orbit)衛(wèi)星完成有源服務，具有不間斷和強實時性的特點[4]。

在BDS系統(tǒng)中，連續(xù)和穩(wěn)定的衛(wèi)星精密軌道測定與預報是衛(wèi)星導航系統(tǒng)服務性能及可用性的重要保證。受地球非球形引力等攝動因素的影響，衛(wèi)星軌道根數(shù)存在長期攝動影響，衛(wèi)星位置在經(jīng)度和緯度方向均存在漂移問題。但是受軌位和頻率資源的限制，北斗衛(wèi)星需要定期進行軌道機動。特別是對于軌位資源相對匱乏的GEO衛(wèi)星，其軌道機動的頻度更高，部分衛(wèi)星的機動周期達到25 d左右[7]。

實時、不間斷的現(xiàn)代導航服務需求對導航系統(tǒng)的運行控制能力提出了新的更高要求，傳統(tǒng)的精密定軌理論難以處理衛(wèi)星姿軌控期間和姿軌控后的精密軌道確定。

對于GPS系統(tǒng)，當衛(wèi)星機動后，通常將該衛(wèi)星從星座中剔除，機動前和機動期間的觀測數(shù)據(jù)將不再參與定軌，衛(wèi)星將重新積累數(shù)據(jù)，待衛(wèi)星定軌狀態(tài)穩(wěn)定之后，再重新進行星座組網(wǎng)處理。受地面跟蹤幾何條件、觀測數(shù)據(jù)量、測量誤差和上行注入等諸多因素的影響，在GPS系統(tǒng)地面運控端通常需要重新積累4～6 h的數(shù)據(jù)后再進行定軌處理，其定軌精度才能恢復到正常水平，因此在軌衛(wèi)星軌道機動后存在6 h左右的不可用狀態(tài)。

對于BDS衛(wèi)星機動期間定軌，文獻[10—11]中提出了機動力建模條件下的GEO衛(wèi)星機動期間定軌方法，但由于跟蹤弧段有限，同時衛(wèi)星在軌機動噴力數(shù)據(jù)存在誤差，基于機動力建模的動力學定軌難以實現(xiàn)米級的衛(wèi)星軌道測定，其精度大約為50 m。在衛(wèi)星機動后的恢復期間，文獻[9]提出了時間同步約束的快速動力學定軌方法，實現(xiàn)了機動后4 h的軌道快速恢復，其定軌位置精度優(yōu)于10 m，有效解決了衛(wèi)星頻繁軌道機動后的快速恢復問題。

但是從單顆衛(wèi)星可用度而言，北斗衛(wèi)星在機動期間(通常約2 h)與恢復4 h內(nèi)共6 h左右時間期間，衛(wèi)星廣播星歷參數(shù)無法達到正常服務水平，機動衛(wèi)星仍處于不可用狀態(tài)，北斗RNSS服務存在下降風險。對于北斗RDSS服務，雖然其對衛(wèi)星軌道精度的要求不高，但從北斗系統(tǒng)整體服務性能來看，仍然存在較大的提升空間。

當衛(wèi)星進行軌道機動后，筆者前期提出了基于多項式擬合[12]和基于星歷擬合[13-14]兩種運動學定軌方法，并利用北斗GEO衛(wèi)星和MEO衛(wèi)星的實測數(shù)據(jù)進行了驗證。試驗結果表明在10 min左右的超短弧條件下，GEO衛(wèi)星定軌精度約為20 m[13]，MEO衛(wèi)星定軌精度約為8 m[14]。受偽距多徑和接收機系統(tǒng)差等因素的制約，該定軌方法存在一定的優(yōu)化空間。同時需要對兩種運動學定軌方法在北斗GEO/IGSO/MEO混合星座的應用特性進行系統(tǒng)性分析，特別是超短弧條件下的軌道預報能力論證，有助于北斗衛(wèi)星精密定軌體制的確立與優(yōu)化。

本文立足北斗衛(wèi)星超短弧跟蹤條件下的快速軌道確定問題，系統(tǒng)地研究了基于多項式擬合和基于星歷擬合兩種超短弧運動學定軌方法，采用CNMC(code noise and multipath correction)算法的多徑削弱技術和基于并置比對的接收機系統(tǒng)差解算技術，有效實現(xiàn)了北斗機動衛(wèi)星測軌數(shù)據(jù)質(zhì)量控制，并利用北斗3種不同類型的衛(wèi)星實測數(shù)據(jù)進行了定軌試驗驗證，完成了北斗衛(wèi)星運動學定軌的系統(tǒng)論證，得到了一些有益結論，其中北斗GEO衛(wèi)星超短弧條件下的定軌精度得到大幅提升。

1 超短弧運動學定軌方法

超短弧運動學定軌方法主要有兩種，其一是基于多項式擬合的運動學定軌，其二是基于星歷擬合的運動學定軌。下面對兩種方法進行簡要介紹。

1.1 基于多項式擬合的運動學定軌模型

在短弧跟蹤條件下，采用單點定位只能得到一系列時間連續(xù)但不光滑的衛(wèi)星位置，而無法直接獲得衛(wèi)星速度。當跟蹤弧段很短時，衛(wèi)星運動軌道是平滑的，因此可以利用1個n階多項式來近似地描述衛(wèi)星的地固系位置，建立基于多項式擬合的短弧運動學法定軌模型，利用多項式與衛(wèi)星位置或星地距離之間的模型關系，對多項式系數(shù)進行估計，實現(xiàn)超短弧運動學定軌。當m個地面跟蹤站實現(xiàn)站間時間同步的前提下同時對衛(wèi)星進行跟蹤，其觀測模型如式(1)所示

(1)

式中，ρ(j,t)表示t時刻第j個跟蹤站觀測的偽距；(Xt,Yt,Zt)表示t時間衛(wèi)星的地固系位置；(Xsj,Ysj,Zsj)表示第j個跟蹤站的地固系位置；Ai(i=1,…,n)為衛(wèi)星X方向多項式系數(shù)；Bi(i=1,…,n)為衛(wèi)星Y方向多項式系數(shù)；Ci(i=1,…,n)為衛(wèi)星Z方向多項式系數(shù)；τ表示t時刻的衛(wèi)星鐘差。對觀測方程進行線性化處理，可以構建相應的誤差方程，依據(jù)最小二乘法解算待估參數(shù)[A0,A1,…,An,B0,B1,…,Bn,C0,C1,…,Cn,τ]；當衛(wèi)星鐘差通過其他手段(如衛(wèi)星雙向頻率傳遞)測定后，衛(wèi)星鐘差參數(shù)τ可以不需要估計，定軌解算的DOP(dilution of precision)值將得到較大改善。在具體的實現(xiàn)過程中，可以采用基于單點定位的多項式擬合法和基于偽距的多項式擬合法，詳細原理參考文獻[12]。為了提高定軌精度，本文主要采用基于偽距的多項式擬合法，即直接建立多項式與原始偽距觀測之間的數(shù)學關系。

1.2 基于星歷擬合的運動學定軌模型

在上述16參數(shù)的廣播星歷模型中，toe和6個開普勒軌道參數(shù)在描述衛(wèi)星軌道中是必需的，3個一階攝動參數(shù)能夠反映重要軌道參數(shù)的一階變化率，而6個調(diào)和攝動參數(shù)主要反映了各種類型的軌道攝動。由于調(diào)和改正項主要與衛(wèi)星軌道周期相關，在超短弧條件下這6個參數(shù)估計結果嚴重失真。

2 數(shù)據(jù)質(zhì)量控制

接收機觀測數(shù)據(jù)質(zhì)量是北斗衛(wèi)星運動學定軌精度的重要保證。受衛(wèi)星星座和觀測環(huán)境等因素的制約，北斗衛(wèi)星偽距觀測數(shù)據(jù)中存在較嚴重多徑問題(特別是GEO衛(wèi)星)[15-16]，同時接收機時延也存在一定的系統(tǒng)差偏差，并且隨著在軌運行存在長期漂移的問題。下面主要針對接收機系統(tǒng)差和偽距多徑問題，研究和設計接收機數(shù)據(jù)質(zhì)量控制方法。

2.1 基于并置比對的接收機系統(tǒng)差標校方法

接收機通道時延存在標定誤差，同時設備在軌運行存在時延長期漂移的問題。在超短弧跟蹤條件下，由于觀測數(shù)據(jù)少，測站和衛(wèi)星設備系統(tǒng)差難以建模，直接估計將導致定軌不收斂，因此接收機系統(tǒng)差的精確標校是制約超短弧運動學定軌的關鍵問題。

對于同一測站而言，為了保證地面系統(tǒng)運行的可靠性，通常安裝2—3臺接收機進行相互備份。接收機通常有外頻標和內(nèi)頻標驅(qū)動共兩種工作模式，在地面處理系統(tǒng)通常是采用外頻標驅(qū)動，所有接收機均采用相同的外部時頻系統(tǒng)進行驅(qū)動。接收機系統(tǒng)差主要包括不同頻率間相對硬件延遲和相對于外頻標1 pps的絕對硬件延遲兩部分。對于相對硬件延遲，可以通過選取夜間背景場和適當?shù)募s束條件進行解算，其精度可以優(yōu)于1 ns[17]。對于絕對硬件延遲，主要由接收機廠家在出廠前進行標定，但由于標定誤差和長期漂移，存在較大誤差，需要對其進行標校與解算。

在地面處理系統(tǒng)中，通常會在選擇一臺基準接收機作為參考基準，進行各類型業(yè)務處理，如精密定軌和差分完好性處理，并基于該基準接收機完成各測站接收機之間的時間同步。而基準接收機的時延主要通過高精度時延測量儀器測量得到，其長期變化量主要通過專門的監(jiān)測設備進行測量，同時配合系統(tǒng)差解算方法進行整網(wǎng)平差。在進行超短弧條件下定軌中需要高精度站鐘差的支持，需要將實際使用的定軌接收機與基準接收機進行時延標校，實現(xiàn)兩臺接收機之間的高精度時間同步。

在同一測站可以選取基準接收機作為參考基準，利用O-C并置比對的方法完成其他接收機的絕對時延標校。其主要原理如下

(2)

式中，R1、R2分別為接收機1和2的理論星地距離；ρ1和ρ2分別為接收機1和2的下行偽距觀測值；δts為衛(wèi)星鐘差；δtk為測站鐘差；ΔDtrop和ΔDion分別為對流層和電離層延遲誤差；ΔDrel為廣義相對論改正；ΔDtide為地球潮汐改正；ΔDSatTran為衛(wèi)星發(fā)射時延；τ21為接收機2的絕對硬件延遲；ε1和ε2分別為接收機1和2偶然誤差。通過計算兩個接收機O-C之差，可以完成接收機2絕對硬件延遲的解算與標校。

表1給出北京和海南兩個測站的接收機系統(tǒng)差多天解算結果，兩站均選取1號機作為參考基準，對2號和3號機進行系統(tǒng)差解算與標校。從結果中可以看出，與基準接收機相比較，不同接收機系統(tǒng)差大約為10～20 ns量級，單臺接收機多天的一致性較好，該誤差對運動學定軌精度的將達到數(shù)十米。不同接收機系統(tǒng)差平均RMS優(yōu)于1 ns，綜合多臺接收機系統(tǒng)差解算結果，平均標準差為0.6 ns，因此可以認為基于并置比對的接收機系統(tǒng)差解算精度優(yōu)于1 ns，可以有效削弱系統(tǒng)差對短弧運動學定軌的影響。

表1 接收機系統(tǒng)差標校結果

對于不同測站接收機系統(tǒng)偏差標校問題，如果異站接收機與基準接收機之間直接采用O-C并置比對法，難以消除公共的傳播路徑誤差。同時，由于兩個接收機不是同源時頻信號驅(qū)動，異站接收機時延與本地鐘差是耦合的，難以直接分離。因此，對于異站接收機系統(tǒng)差的標校，直接采用載波共視法或多星定軌法完成測站鐘差的估計。

2.2 基于CNMC的多徑削弱方法

WAAS(wide area augmentation system)系統(tǒng)提出一種采用CMC(code-minus-carrier)觀測量來改正GPS偽距觀測多路徑效應的方法[18-19],有效改善了WAAS偽距觀測量的精度。北斗導航系統(tǒng)采用GEO/IGSO/MEO混合星座設計，基于地面監(jiān)測站的實測數(shù)據(jù)分析表明，GEO和IGSO衛(wèi)星偽距觀測量中包含了明顯的多徑噪聲，且GEO的多徑噪聲尤為嚴重[20-21]，部分測站多徑接近1 m。

考慮到載波相位的多徑噪聲不超過3 cm，比偽距多徑噪聲小兩個數(shù)量級，文獻[22]根據(jù)偽距和雙頻載波相位的組合提取偽距觀測量中多徑噪聲，并提出偽距多徑噪聲改正的實時算法CNMC。該算法原理是利用雙頻相位組合計算每個歷元的電離層延遲，實現(xiàn)平滑單頻點偽距觀測的目的。該方法消除偽距多徑誤差的效果顯著,定位解算精度改善[23-24]。

下面選取國內(nèi)分布的7臺接收機對北斗GEO和IGSO衛(wèi)星實測數(shù)據(jù)，進行CNMC前后的多徑對比試驗，表2給出了多徑RMS統(tǒng)計對比結果。

從表中可以看出CNMC算法多徑改正效果明顯，其中GEO衛(wèi)星偽距多徑改善率達到60.17%，IGSO衛(wèi)星偽距多徑改善率達到39.54%，GEO衛(wèi)星的改善效果更加明顯?；贑NMC的多徑削弱方法，可以對接收機的偽距觀測質(zhì)量進行有效控制，有利于提升超短弧跟蹤條件下的運動學定軌精度。

表2 CNMC算法處理前后各星偽距觀測量中的多徑RMS

3 定軌試驗與分析

上文提出了兩種運動學定軌方法和接收機數(shù)據(jù)質(zhì)量控制方法，并詳細推導了具體實現(xiàn)過程。兩種運動學方法充分利用了高采樣率的測軌數(shù)據(jù)，減少了結果的噪聲，其優(yōu)點在于不需要累積長期測軌數(shù)據(jù)，可以實現(xiàn)近實時快速計算，克服了動力學法定軌發(fā)散和單點定位無法獲取速度信息的不足，實現(xiàn)了短弧跟蹤條件下的衛(wèi)星精密定軌。下面利用北斗系統(tǒng)GEO/IGSO/MEO衛(wèi)星的實測偽距觀測數(shù)據(jù)，對兩種定軌方法進行試驗驗證。

3.1 定軌方案與評估方法

為了驗證兩種超短弧運動學定軌方法的精度，采用北斗系統(tǒng)7個測站CNMC修正后的偽距數(shù)據(jù)進行定軌試驗，各測站之間通過所有在軌衛(wèi)星多星聯(lián)合定軌實現(xiàn)站間時間同步，各接收機的系統(tǒng)差基于2.1節(jié)的并置比對法完成了標校與修正，衛(wèi)星鐘差采用雙向時間頻率法測定的歷元鐘差值[25]-[27]。主要定軌處理策略如表3所示。

表3 超短弧運動學定軌處理策略

由于衛(wèi)星機動期間的定軌精度有限，為了評估多星多站條件下的機動定軌精度，采用非機動期間的數(shù)據(jù)進行定軌處理，利用長弧多星定軌的軌道產(chǎn)品作為評估標準，統(tǒng)計衛(wèi)星R、T、N 3個方向和位置誤差，表4給出了多星定軌精度的統(tǒng)計結果。從結果中可以看出，多星定軌的平均位置精度為1.28 m，能夠評估米級的超短弧運動學定軌精度。

表4 多星定軌精度統(tǒng)計(3 d定軌，重迭1 d)

Tab.4 Multi-satellite precise orbit determination accuracy (3 days tracking arc, 1 day overlap) m

3.2 定軌精度

下面采用2014年8月3日至4日C02(GEO)、C09(IGSO)和C11(MEO)衛(wèi)星的實測數(shù)據(jù)進行定軌試驗，每顆衛(wèi)星分別采用2～3個弧段數(shù)據(jù)進行定軌，每個弧段10 min，定軌方法分別為基于多項式擬合的運動學定軌和基于星歷擬合的運動學定軌。在基于多項式擬合的運動學定軌中，由于觀測數(shù)據(jù)只有10 min，同時考慮到待估參數(shù)的相關性和軌道預報精度的衰減，多項式的階數(shù)選擇為2階。

將定軌結果與多星定軌結果進行比較，統(tǒng)計兩個弧段R、T、N 3個方向和URE的誤差均值，同時為了進行結果比較，增加了單點定位結果的精度評估，試驗結果統(tǒng)計如表5所示。

表5 不同定軌方法的定軌精度-1(10 min弧長)

圖1 C02衛(wèi)星基于星歷擬合的定軌精度(弧段1)Fig.1 C02 satellite orbit determination accuracy based on ephemeris parameter fitting method (arc 1)

下面對上述運動學定軌精度進行分析：

(2) 從不同的定軌方法來看，綜合比較3類不同衛(wèi)星的定軌位置精度，基于星歷擬合的超短弧運動學定軌方法精度最高，其次是基于多項式擬合。究其原因，基于星歷擬合的定軌方法帶有衛(wèi)星軌道的部分動力學信息，能夠更好地描述衛(wèi)星軌道的物理特性，所以其精度更高。但從位置精度而言，基于多項式擬合的定軌方法同樣能夠?qū)崿F(xiàn)衛(wèi)星軌道的米級定軌精度，滿足衛(wèi)星軌道機動時的RDSS服務所需的衛(wèi)星星歷。

(3) 從不同類型衛(wèi)星的定軌位置精度來看，GEO衛(wèi)星雖然高度更高，但其位置精度最好，究其原因是定軌數(shù)據(jù)預處理中對衛(wèi)星系統(tǒng)差進行了有效標校，GEO衛(wèi)星的靜地特性使其相對于地球的運動速度小，衛(wèi)星和測站系統(tǒng)差標校精度高，而且隨時間變化緩慢。而IGSO和MEO衛(wèi)星相對于地球的運動速度更大，測站相對于衛(wèi)星的多普勒(偽距變化率)大，衛(wèi)星和測站系統(tǒng)差解算精度更低，隨時間變化較快。而衛(wèi)星和測站系統(tǒng)差通常是采用1d前的標定結果，因此不同類型衛(wèi)星表現(xiàn)不同的定軌精度。

(4) 從衛(wèi)星的速度確定精度來看，IGSO和MEO衛(wèi)星只能獲得6cm/s左右的速度精度，GEO衛(wèi)星可以獲得1cm/s左右的速度精度，與IGSO/MEO衛(wèi)星相比提高5倍，究其原因，GEO衛(wèi)星相對于地球幾乎不動，運動速度小，更有利于速度信息的建模。

(5) 從運動學定軌在北斗系統(tǒng)的應用來看，與傳統(tǒng)的機動期間單點定位方法相比，采用運動學定軌方法后軌道精度得到明顯改善，提高了1倍左右。特別是GEO衛(wèi)星，可以實現(xiàn)3.27m左右的定軌精度，由于GEO衛(wèi)星機動期間需要提供RDSS服務所需的星歷，運動學定軌方法能夠滿足其對GEO衛(wèi)星軌道精度的需求，從而確保連續(xù)的RDSS服務。

(6) 與文獻[12—13]中定軌精度相比較，MEO衛(wèi)星定軌精度相當，而GEO衛(wèi)星定軌精度得到較大改善，這主要得益于有效的偽距數(shù)據(jù)質(zhì)量控制，對接收機系統(tǒng)差和多徑誤差進行了有效削弱。

為了進一步說明運動學定軌方法的普適性，需要利用更多的衛(wèi)星進行短弧運動學定軌試驗。由于GEO衛(wèi)星相對于地面跟蹤幾何條件幾乎不變，下面僅利用2014年8月3日C06、C08和C12衛(wèi)星進行分析，試驗方法和條件與前面相同。表6給出了不同方法的統(tǒng)計結果。

表6 不同定軌方法的定軌精度-2(10 min弧長)

綜合表5和表6的結果可以看出，表6中的整體定軌精度優(yōu)于表5，原因是不同跟蹤弧段中衛(wèi)星相對地面測站的高度角不同。其中表5中各弧段的衛(wèi)星高度角大約為10°～40°，而表6中衛(wèi)星高度角在60°左右，因此定軌精度略高。從整體而言，運動學定軌方法對于提高衛(wèi)星姿軌控期間的定軌精度具有較好的應用價值。

3.3 預報精度

導航用戶采用的是衛(wèi)星預報軌道，因此軌道預報精度非常關鍵。為了分析兩種超短弧運動學定軌方法的預報性能，采用2014年8月3日C02(GEO)、C09(IGSO)和C11(MEO) 3種類型衛(wèi)星的實測數(shù)據(jù)進行定軌與預報，定軌弧段為10 min，預報弧段為若干分鐘。由于C06、C08和C12衛(wèi)星幾何跟蹤條件更好，定軌精度也略高，在此不再進行預報精度分析與統(tǒng)計。下面將定軌結果與多星定軌結果進行比較，統(tǒng)計R、T、N 3個方向和URE的誤差，試驗結果統(tǒng)計如下。

下面對上述基于多星定軌的運動學軌道預報精度評估結果進行分析：

(1) 對于C02衛(wèi)星，在10 min左右的超短弧跟蹤條件下，綜合考慮多星定軌誤差和誤差傳播規(guī)律，運動學定軌方法預報10 min的位置精度和URE精度分別優(yōu)于為7.2 m和1 m，從位置和URE精度來看，機動衛(wèi)星與非機動衛(wèi)星的定軌精度相當。因此，從運動學定軌的預報精度來看，能夠同時滿足北斗RNSS和RDSS服務性能的需求，從而提高單顆GEO衛(wèi)星的可用度，實現(xiàn)北斗系統(tǒng)可用度的整體提升。

表7 不同定軌方法的預報精度(10 min弧長)

(2) 對于C09和C11衛(wèi)星，在10 min左右的超短弧跟蹤條件下，兩種運動學定軌方法預報10 min的軌道位置精度表現(xiàn)了很大差異，其中基于多項式擬合方法達到百米量級，而星歷擬合方法為30～50 m精度水平，究其原因，星歷擬合方法包含了衛(wèi)星軌道的部分動力學信息，更加有利于提升軌道的預報能力。

(3) 對比不同類型衛(wèi)星的軌道預報精度，為了實現(xiàn)優(yōu)于3 m的URE精度，在10 min的跟蹤弧段條件下，GEO衛(wèi)星能夠達到10 min，甚至更長，而IGSO和MEO衛(wèi)星只能預報3～6 min，其中基于星歷擬合法優(yōu)勢更加明顯。

(4) 兩種運動學定軌方法克服了短弧跟蹤條件下動力學法定軌發(fā)散問題，同時解決了單點定位不光滑和無法獲得衛(wèi)星速度信息等問題，實現(xiàn)了短弧跟蹤條件下的北斗衛(wèi)星運動學定軌與預報。從預報性能來看，基于星歷擬合的運動學定軌方法更加有效，有利于提高軌道的預報能力。

(5) 從北斗系統(tǒng)RDSS服務來看，目前GEO衛(wèi)星機動期間由于定軌精度不高，軌道預報精度衰減非?？欤诰€系統(tǒng)通常是1 min左右進行一次幾何法定軌和擬合預報，而采用本文提出的運動學定軌方法，不僅提高了衛(wèi)星定軌和預報精度，而且可以降低系統(tǒng)的處理頻度和星歷注入頻度，降低地面系統(tǒng)的處理強度和負荷，以提升超短弧跟蹤條件下RDSS的服務性能。

4 結 論

本文針對北斗衛(wèi)星姿軌控后超短弧跟蹤條件下的軌道快速確定難題，系統(tǒng)地研究了基于多項式擬合和基于星歷擬合兩種運動學定軌方法，給出了定軌解算模型與實現(xiàn)過程，設計了北斗衛(wèi)星數(shù)據(jù)質(zhì)量控制方法，利用北斗GEO/IGSO/MEO衛(wèi)星的實測偽距數(shù)據(jù)進行了定軌試驗，主要結論如下：

(1) 兩種數(shù)據(jù)質(zhì)量控制方法有效削弱了系統(tǒng)差和多徑誤差對定軌精度的影響，其中基于并置比對的接收機系統(tǒng)差標定精度優(yōu)于1 ns，GEO和IGSO衛(wèi)星偽距多徑改善率分別達到60.17%和39.54%。

(2) 從定軌位置精度來看，基于多項式擬合和基于星歷擬合兩種運動學定軌方法是有效的，在10 min跟蹤條件下運動學定軌位置精度，GEO、IGSO和MEO衛(wèi)星定軌精度分別為3.27 m、8.19 m和5.90 m，經(jīng)過數(shù)據(jù)質(zhì)量控制，GEO衛(wèi)星的運動學定軌精度得到較大改善。

(3) 從定軌速度精度來看，基于星歷擬合的定軌方法帶有軌道的動力學信息，更好地描述衛(wèi)星軌道物理特性，其速度精度更高。GEO衛(wèi)星的靜地特性有利于速度建模，IGSO和MEO衛(wèi)星只能獲得6 cm/s左右的速度精度，GEO衛(wèi)星可以獲得1 cm/s左右的速度精度，與IGSO/MEO衛(wèi)星相比提高5倍。

(4) 從定軌預報精度來看，GEO衛(wèi)星的超短弧運動學定軌預報10 min能夠?qū)崿F(xiàn)優(yōu)于7.2 m的軌道位置精度，URE精度優(yōu)于1 m，其精度水平與非機動衛(wèi)星相當。在相同條件下為實現(xiàn)優(yōu)于3 m的URE精度，IGSO和MEO衛(wèi)星只能預報3～6 min。

(5) 與傳統(tǒng)的機動期間定軌方法相比，采用運動學定軌方法后軌道精度得到明顯改善，提高了1倍左右，能夠有效滿足GEO衛(wèi)星機動期間的RDSS服務對軌道精度的需求，確保北斗RDSS服務連續(xù)和穩(wěn)定。

總之，兩種運動學法定軌方法克服了短弧跟蹤條件下傳統(tǒng)動力學法定軌和單點定位中諸多問題，實現(xiàn)了超短弧條件下的北斗衛(wèi)星軌道快速確定，能夠適用于不同軌道類型的北斗衛(wèi)星。本文的方法和試驗成果能夠有效支撐北斗RDSS向全球系統(tǒng)擴展，為北斗RDSS服務走向全球提供技術支撐。

[1] 楊元喜. 北斗衛(wèi)星導航系統(tǒng)的進展、貢獻與挑戰(zhàn)[J]. 測繪學報, 2010, 39(1): 1-6. YANG Yuanxi. Progress, Contribution and Challenges of Compass/BeiDou Satellite Navigation System[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(1): 1-6.

[2] YANG Yuanxi, LI Jinlong, XU Junyi, et al. Contribution of the Compass Satellite Navigation System to Global PNT Users[J]. Chinese Science Bulletin, 2011, 56(26): 2813-2819. DOI: 10.1007/s11434-011-4627-4.

[3] YANG Yuanxi, LI Jinlong, WANG Aibing, et al. Preliminary Assessment of the Navigation and Positioning Performance of BeiDou Regional Navigation Satellite System[J]. Science China Earth Sciences, 2014, 57(1): 144-152. DOI: 10.1007/s11430-013-4769-0.

[4] 譚述森. 廣義衛(wèi)星無線電定位報告原理及其應用價值[J]. 測繪學報, 2009, 38(1): 1-5.DOI:10.3321/j.issn:1001-1595.2009.01.001. TAN Shusen. Theory and Application of Comprehensive RDSS Position and Report[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(1): 1-5.DOI:10.3321/j.issn:1001-1595.2009.01.001.

[5] 譚述森. 衛(wèi)星導航定位工程[M]. 北京: 國防工業(yè)出版社, 2007. TAN Shusen. The Engineering of Satellite Navigation and Positioning[M]. Beijing: The National Defense Industry Press, 2007.

[6] 劉穎, 劉根友, 陳曉峰, 等. 北斗二代衛(wèi)星導航系統(tǒng)IGSO星增加RDSS功能問題探討[J]. 大地測量與地球動力學, 2012, 32(3): 72-75. LIU Ying, LIU Genyou, CHEN Xiaofeng, et al. Discussion on Adding RDSS to IGSO Satellite in 2ND Generation BeiDou Navigation Satellite System[J]. Journal of Geodesy and Geodynamics, 2012, 32(3): 72-75.

[7] 黃勇, 胡小工, 黃珹, 等. 利用CAPS測距數(shù)據(jù)確定GEO衛(wèi)星變軌期間的軌道[J]. 中國科學 G輯: 物理學 力學 天文學, 2008, 38(12): 1750-1758. HUANG Yong, HU Xiaogong, HUANG Cheng, et al. Precise Orbit Determination of a Maneuvered GEO Satellite Using CAPS Ranging Data[J]. Science in China Series G: Physics, Mechanics and Astronomy, 2009, 52(3): 346-352.

[8] 楊旭海, 李志剛, 馮初剛, 等. GEO衛(wèi)星機動后的星歷快速恢復方法[J]. 中國科學 G輯: 物理學 力學 天文學, 2008, 38(12): 1759-1765. YANG Xuhai, LI Zhigang, FENG Chugang, et al. Methods of Rapid Orbit Forecasting after Maneuvers for Geostationary Satellites[J]. Science in China Series G: Physics, Mechanics and Astronomy, 2009, 52(3): 333-338.

[9] GUO Rui, ZHOU Jianhua, HU Xiaogong, et al. Precise Orbit Determination and Rapid Orbit Recovery Supported by Time Synchronization[J]. Advances in Space Research, 2015, 55(12): 2889-2898. DOI: 10.1016/j.asr.2015.03.001.

[10] 郭睿, 李曉杰, 周建華, 等. 機動力建模條件下的GEO衛(wèi)星機動期間定軌[J]. 測繪科學技術學報, 2013, 30(5): 465-470. GUO Rui, LI Xiaojie, ZHOU Jianhua, et al. Precise Orbit Determination for the GEO Satellite Maneuver Based on the Thrust Force Model[J]. Journal of Geomatics Science and Technology, 2013, 30(5): 465-470.

[11] 李曉杰, 黃勇, 郭睿, 等. 基于機動力模型的GEO衛(wèi)星恢復期間定軌[J]. 空間科學學報, 2015, 35(6): 730-737. LI Xiaojie, HUANG Yong, GUO Rui, et al. Orbit Determination Method Based on Thrust Force Model for GEO Satellite in the Period of Orbit Recovery[J]. Chinese Journal of Space Science, 2015, 35(6): 730-737.

[12] 郭睿, 胡小工, 黃勇, 等. 導航衛(wèi)星的短弧運動學定軌[J]. 飛行器測控學報, 2010, 29(1): 60-67. GUO Rui, HU Xiaogong, HUANG Yong, et al. Navigation Satellite Short-arc Kinematic Orbit Determination[J]. Journal of Spacecraft TT&C Technology, 2010, 29(1): 60-67.

[13] 郭睿, 周建華, 胡小工, 等. 一種地球靜止軌道衛(wèi)星的快速恢復定軌方法[J]. 測繪學報, 2011, 40(S): 19-25, 30. GUO Rui, ZHOU Jianhua, HU Xiaogong, et al. A Strategy of Rapid Orbit Recovery for the Geostaionary Satellite[J]. Acta Geodaetica et Cartographica Sinica, 2011, 40(S): 19-25, 30.

[14] 郭睿, 胡小工, 黃勇, 等. 基于星歷擬合的短弧運動學定軌[J]. 宇航學報, 2010, 31(2): 416-422. GUO Rui, HU Xiaogong, HUANG Yong, et al. Short-arc Kinematic Orbit Determination Based on Ephemeris Fitting[J]. Journal of Astronautics, 2010, 31(2): 416-422.

[15] 馮曉超, 金國平, 范建軍, 等. GNSS接收機偽距測量中的多徑效應試驗分析[J]. 現(xiàn)代電子技術, 2013, 36(5): 77-81. FENG Xiaochao, JIN Guoping, FAN Jianjun, et al. Experimentation and Analysis of Multipath Effect in Pseudo-range Measurement of GNSS Receiver[J]. Modern Electronics Technique, 2013, 36(5): 77-81.

[16] FENG Xiaochao, WU Xiaoli, ZHANG Zhixue, et al. Multipath Mitigation Technique Based on Modifications to GNSS Monitor Station Antennas Field[C]∥Proceedings of China Satellite Navigation Conference (CSNC) 2012. Berlin Heidelberg: Springer, 2012: 655-666.

[17] 吳曉莉, 平勁松, 劉利, 等. 區(qū)域衛(wèi)星導航系統(tǒng)硬件延遲解算[J]. 武漢大學學報(信息科學版), 2011, 36(10): 1218-1221. WU Xiaoli, PING Jinsong, LIU Li, et al. Hardware Delay Solution of Regional Satellite Navigation System[J]. Geomatics and Information Science of Wuhan University, 2011, 36(10): 1218-1221.

[18] SHALLBERG K, SHENG Fang. WAAS Measurement Processing: Current Design and Potential Improvements[C]∥Proceedings of 2008 IEEE/ION Position, Location and Navigation Symposium. Monterey, CA: IEEE, 2008: 253-262.

[19] SCHEMPP T, BURKE J, RUBIN A. WAAS Benefits of GEO Ranging[C]∥Proceedings of the 21st International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS 2008). Savannah, GA: ION, 2008: 12-15.

[20] WANNINGER L, BEER S. BeiDou Satellite-induced Code Pseudorange Variations: Diagnosis and Therapy[J]. GPS Solutions, 2015, 19(4): 639-648. DOI: 10.1007/s10291-014-0423-3.

[21] HAUSCHILD A, MONTENBRUCK O. A Study on the Dependency of GNSS Pseudorange Biases on Correlator Spacing[J]. GPS Solutions, 2016, 20(2): 159-171. DOI: 10.1007/s10291-014-0426-0.

[22] WU Xiaoli, ZHOU Jianhua, WANG Gang, et al. Multipath Error Detection and Correction for GEO/IGSO Satellites[J]. Science China Physics, Mechanics and Astronomy, 2012, 55(7): 1297-1306.

[23] CAO Yueling, HU Xiaogong, WU Bing, et al. The Wide-area Difference System for the Regional Satellite Navigation System of Compass[J]. Science China Physics, Mechanics and Astronomy, 2012, 55(7): 1307-1315.

[24] 常志巧, 胡小工, 郭睿, 等. CNMC與Hatch濾波方法比較及其在北斗相對定位中的精度分析[J]. 中國科學: 物理學 力學 天文學, 2015, 45(7): 079508. CHANG Zhiqiao, HU Xiaogong, GUO Rui, et al. Comparison between CNMC and Hatch Filter & Its Precision Analysis for BDS Precise Relative Positioning[J]. Scientia Sinica: Physica, Mechanica & Astronomica, 2015, 45(7): 079508.

[25] 劉利, 韓春好. 衛(wèi)星雙向時間比對及其誤差分析[J]. 天文學進展, 2004, 22(3): 219-226. LIU Li, HAN Chunhao. Two Way Satellite Time Transfer and Its Error Analysis[J]. Progress in Astronomy, 2004, 22(3): 219-226.

[26] 劉利, 韓春好. 地心非旋轉(zhuǎn)坐標系中的TWSTT計算模型[J]. 測繪學院學報, 2004, 21(2): 96-98. LIU Li, HAN Chunhao. Model of TWSTT in the Geocentric Non-rotating Coordinate System[J]. Journal of Institute of Surveying and Mapping, 2004, 21(2): 96-98.

[27] 劉利, 朱陵鳳, 韓春好, 等. 星地無線電雙向時間比對模型及試驗分析[J]. 天文學報, 2009, 50(2): 189-196. LIU Li, ZHU Lingfeng, HAN Chunhao, et al. The Model of Two-way Radio Time Transfer between the Earth and Satellites and Analysis of Its Experiment[J]. Acta Astronomica Sinica, 2009, 50(2): 189-196.

(責任編輯：叢樹平)

Kinematic Orbit Determination Method Optimization and Test Analysis for BDS Satellites with Short-arc Tracking Data

GUO Rui1,2,CHEN Jinping1,ZHU Lingfeng1,TANG Guifen1,WU Xiaoli1

1. Beijing Satellite Navigation Center, Beijing 100094, China; 2. State Key Laboratory of Geodesy and Earth’s Dynamics, Wuhan 430077, China

Rapid orbit recovery is a puzzle for the BDS satellites after orbit maneuvers. Two kinematic orbit determination methods are studied, with two orbit determination models being established. The receiver system error and serious multipath error exist in the BDS system. The co-location method is proposed to estimate and calibrate the receiver system errors. A CNMC (code noise and multipath correction) method is introduced to weaken the multipath error. Therefore the data quality is controlled efficiently for the receivers in the short tracking arc. The GEO/IGSO/MEO real data is emploied to carry out tests and validation. Using 10 min short tracking arc, the kinematic precise orbit determination accuracy is about 3.27 m for the GEOs, and 8.19 m for the IGSOs, and 5.9 m for the MEOs. Rapid orbit determination is achieved, which satisfying the orbit requirements from the BDS RDSS services. The kinematic precise orbit determination method also supports the RDSS service walking up to the global world.

BDS; precise orbit determination; orbit maneuver; kinematic precise orbit determination; multipath error

The National Natural Science Foundation of China (Nos. 41204022; 41374038); The State Key Laboratory of Geodesy and Earth’s Dynamics Foundation of China (No. SKLGED2017-3-3-E)

GUO Rui(1982—), male, PhD, senior engineer, majors in satellite navigation, precise orbit determination and time synchronization.

CHEN Jinping

郭睿，陳金平，朱陵鳳,等.北斗衛(wèi)星超短弧運動學定軌方法優(yōu)化與試驗分析[J].測繪學報，2017,46(4)：411-420.

10.11947/j.AGCS.2017.20160361. GUO Rui, CHEN Jinping, ZHU Lingfeng,et al.Kinematic Orbit Determination Method Optimization and Test Analysis for BDS Satellites with Short-arc Tracking Data[J]. Acta Geodaetica et Cartographica Sinica,2017,46(4):411-420. DOI:10.11947/j.AGCS.2017.20160361.

P228

A

1001-1595(2017)04-0411-10

國家自然科學基金(41204022; 41374038)；大地測量與地球動力學國家重點實驗室開放研究基金(SKLGED2017-3-3-E)

2016-07-15

郭睿(1982—),男,博士,高級工程師,研究方向為衛(wèi)星導航、精密定軌與時間同步。

E-mail： shimbarsalon@163.com

陳金平

E-mail： jinping_chen@sina.com

修回日期： 2017-02-11