鞏秀強(qiáng)，陳俊平，周善石，吳 斌

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

1 引 言

中國北斗衛(wèi)星導(dǎo)航系統(tǒng)(Beidou navigation satellite system,BDS)具備無線電星地雙向時間比對測量技術(shù)[1]。國內(nèi)外已有較多雙向時間比對技術(shù)的研究，星地雙向時間同步技術(shù)通過上下行偽距求差，消除了共有誤差影響，也減小與信號頻率有關(guān)的電離層延遲等誤差[2,3]。綜合考慮目前時間比對技術(shù)各誤差源的精度量級[4]，星地?zé)o線電雙向法的理論精度可達(dá)到約幾百皮秒的量級[5,6]。本文利用北斗雙向?qū)崪y數(shù)據(jù)[7]，分析了電離層誤差和衛(wèi)星相位中心誤差對雙向時間比對精度的影響及其改正方法，統(tǒng)計了雙向鐘差的隨機(jī)噪聲。

衛(wèi)星鐘差模型是廣播電文的重要組成部分，其精度直接影響導(dǎo)航系統(tǒng)的服務(wù)性能[8?10]。BDS 利用雙向鐘差點，采用短期、中期與長期相結(jié)合的多項式模型擬合生成偏差、鐘速和加速度共3 個參數(shù)，并通過廣播電文播發(fā)給用戶[11]。本文采用衛(wèi)星雙向時頻傳遞(two-way satellite-ground time transfer,TWSTT)鐘差測量作為參考，統(tǒng)計了BDS 廣播電文鐘差模型的精度，并分析了模型存在的問題。

2 模型與誤差

在星地雙向時間同步技術(shù)中，地面站與衛(wèi)星在同一鐘面時T接收對方發(fā)射的測距信號，并測量相應(yīng)的偽距，得到下行時延觀測值Psg和上行時延觀測值Pgs。下行和上行時延的觀測方程如下：

式中，?Tg和?Ts分別表示地面站和衛(wèi)星的鐘面時T相對于導(dǎo)航系統(tǒng)時間的鐘差；DTg和DRg分別為地面站發(fā)射通道時延和接收通道時延，DRs和DTs分別表示衛(wèi)星的發(fā)射通道時延和接收通道時延；ρgs,Igs,Tgs,?τgs,Rgs表示地面站到衛(wèi)星的幾何時延、電離層時延、對流層時延、地球自轉(zhuǎn)效應(yīng)和衛(wèi)星相對論效應(yīng)；ρsg,Isg,Tsg,?τsg,Rsg表示衛(wèi)星到地面站的幾何時延、電離層時延、對流層時延、地球自轉(zhuǎn)和衛(wèi)星相對論效應(yīng)；ETs和ERs表示衛(wèi)星發(fā)射和接收相位中心；εsg、εgs表示上行偽距和下行偽距的隨機(jī)誤差。

將式(2)減去式(1)再除以2,可得到衛(wèi)星鐘相對于地面時間系統(tǒng)北斗時(Beidou time,BDT)的鐘差。即：

由式(3)可知，星地雙向時間同步誤差由偽距觀測誤差、空間傳播誤差、測站與衛(wèi)星發(fā)射接收通道時延、衛(wèi)星發(fā)射與接收相位中心改正以及地球自轉(zhuǎn)效應(yīng)的不完全改正誤差決定。其中地球自轉(zhuǎn)效應(yīng)和相對論效應(yīng)可以采用模型精確計算，其誤差可以忽略[1]。下行與上行對流層延遲誤差Tsg和Tgs可以完全抵消[5]。我們采用了每小時更新的廣播星歷來求取幾何時延，由于上下行軌道誤差可以抵消，幾何時延不完全改正誤差也可以忽略。綜上所述，北斗系統(tǒng)雙向時間比精度取決于測量偽距隨機(jī)誤差、信號在空間雙向傳播的電離層誤差、衛(wèi)星相位中心修正誤差、測站和衛(wèi)星發(fā)射時延誤差以及測站和衛(wèi)星接收時延誤差。北斗時間同步站的發(fā)射與接收相位中心一致，可以相互抵消。

2.1 噪聲分析

采用2016年9月1日北斗衛(wèi)星星地雙向?qū)崪y鐘差進(jìn)行分析，對雙向鐘差序列作1 min窗口的滑動線性階擬合，分析擬合的殘差。為了增加衛(wèi)星可視弧段，除主控站外，北斗區(qū)域衛(wèi)星導(dǎo)航系統(tǒng)布設(shè)多個外場站，并通過C 波段雙向時頻傳遞技術(shù)保持與主控站的時間同步。通過外場站星地雙向時頻傳遞，我們可測得衛(wèi)星鐘差與外場站相對于主控站的相對鐘差之和，并得到衛(wèi)星相對于主控站BDT 的鐘差。如圖1所示，紅點表示主控站直接測得衛(wèi)星鐘差噪聲，綠點表示外場站經(jīng)站間時間同步歸算后的鐘差噪聲。統(tǒng)計得出主控站測得衛(wèi)星鐘差精度RMS為0.08 ns。外場站測得衛(wèi)星鐘差的誤差明顯放大，RMS達(dá)到0.15 ns；這說明外場站與主控站間的站間歸算處理增大了觀測誤差。圖2給出了擬合殘差的統(tǒng)計直方圖，顯示出雙向鐘差噪聲呈正態(tài)分布的特征，具有隨機(jī)誤差特性。

圖1 衛(wèi)星鐘差的誤差序列

圖2 鐘差誤差的直方圖

2.2 電離層延遲改正分析

在雙向測量中，與頻率有關(guān)的電離層延遲誤差雖然減小很多，其值仍可達(dá)4 ns[7]，是影響雙向鐘差精度的主要誤差源。根據(jù)式(3)，電離層延遲對雙向鐘差的影響?tiono表示為：

式中，fdn和fup表示上行和下行的頻率值，TEC為地面站至衛(wèi)星路徑上電離層電子總含量。北斗衛(wèi)星導(dǎo)航系統(tǒng)時間同步站布設(shè)有并置監(jiān)測接收機(jī)[8]，可利用接收機(jī)雙頻實測數(shù)據(jù)求取時間同步站至衛(wèi)星的電離層延遲量[9]。根據(jù)式(1)，將B1 與B2 頻點的偽距觀測值P1和P2做差，可求得同步站至衛(wèi)星路徑上電離層電子總含量TEC，為：

式中f1和f2表示B1 和B2 兩頻點的頻率，tgd12、ifb12分別表示衛(wèi)星與接收機(jī)B1 和B2 頻點的硬件時延之差。相位平滑偽距可以有效地降低偽距噪聲，消除偽距多路徑效應(yīng)[10]。BDS實時采用偽距噪聲和多路徑改正(code noise and multipath correction,CNMC)方法計算接收機(jī)的相位平滑偽距，并供電離層延遲解算使用[11,12]。

相位數(shù)據(jù)觀測精度在毫米級，我們利用雙頻相位觀測數(shù)據(jù)計算電離層延遲，其結(jié)果作為參考，來評估利用雙頻相位平滑偽距計算電離層延遲的精度。采用2016年9月1日北斗監(jiān)測接收機(jī)雙頻實測數(shù)據(jù)，我們分別利用相位平滑偽距與相位觀測數(shù)據(jù)計算電離層延遲對雙向鐘差的影響，將兩者之差畫出。從圖3a)可以看出，對于不同GEO (Geostationary Earth Orbit)衛(wèi)星(不同顏色的線條)，雙頻相位平滑偽距法求得電離層延遲對雙向鐘差影響在0.05 ns 以內(nèi)。當(dāng)相位平滑重新收斂時，電離層延遲的解算精度會受影響，誤差達(dá)到0.1 ns。從圖3b)可以看出，IGSO (Inclined Geosynchronous Satellite Orbit)衛(wèi)星在剛?cè)刖硶r，由于觀測高度角低、偽距質(zhì)量差[13]，電離層延遲的解算精度較低，因此雙向鐘差的誤差可以達(dá)0.5 ns以上。入境一段時間后，誤差逐漸縮小到0.05 ns 以內(nèi)。當(dāng)移動衛(wèi)星剛?cè)刖硶r，雙向鐘差精度較低，在利用雙向鐘差進(jìn)行廣播鐘差參數(shù)擬合時，要將此段數(shù)據(jù)剔除。

圖3 雙頻相位平滑偽距計算的電離層延遲對雙向鐘差的影響

另外一種電離層延遲改正方法為采用電離層模型進(jìn)行計算，其不存在以上的收斂問題。進(jìn)一步分析采用實時電離層模型進(jìn)行電離層誤差的改正。采用的模型為歐洲定軌中心(Center Orbit Determination Europe,CODE)提供的全球電離層模型(Global Ionosphere Map,GIM)，其在中國范圍內(nèi)精度可以達(dá)到10 TECU (TECU=1016m?2)以內(nèi)[9]。以雙頻相位觀測數(shù)據(jù)計算電離層延遲參考來評估模型計算電離層延遲的精度。從圖4a)可以看出，不同GEO 衛(wèi)星(不同顏色線條)的電離層誤差趨勢基本一致，說明GIM 模型在1 h 內(nèi)精度波動明顯，無法表達(dá)電離層變化的高頻部分。對雙向鐘差影響波動幅度達(dá)到0.6 ns，不能滿足雙向鐘差的精度要求。從圖4b)可以看出，IGSO 衛(wèi)星出入境時，由于高度角低，GIM 模型投影函數(shù)誤差較大，這也會造成電離層改正精度大幅降低。

圖4 GIM 求得電離層延遲對雙向鐘差的影響

2.3 衛(wèi)星發(fā)射與接收相位中心改正分析

由于北斗衛(wèi)星下行發(fā)射信號與上行接收信號的相位中心不一致，在進(jìn)行雙向時間比對時，衛(wèi)星相位中心引起的誤差不能抵消，因此需要對衛(wèi)星發(fā)射和接收相位中心進(jìn)行改正。從式(3)中可以看到，衛(wèi)星相位中心改正對雙向鐘差的影響?tecc等于衛(wèi)星下行發(fā)射相位中心改正值減去衛(wèi)星接收相位中心改正值，表示如下：

其中，ETs和ERs表示衛(wèi)星發(fā)射和接收相位中心。衛(wèi)星相位中心改正取決于兩個方面：衛(wèi)星瞬時姿態(tài)以及衛(wèi)星相位中心在星固系中的坐標(biāo)。GPS 衛(wèi)星只有動偏模式，而北斗區(qū)域?qū)Ш较到y(tǒng)衛(wèi)星采用了動偏與零偏兩種姿態(tài)模式[14,15]。若在計算北斗衛(wèi)星相位中心時不考慮此差異，會影響雙向時間比對精度。

采用2016年9月1日數(shù)據(jù)，在動偏與零偏兩種姿態(tài)模式下，我們把衛(wèi)星相位中心改正對雙向鐘差的影響進(jìn)行了比較。由圖5可知，在不同姿態(tài)控制模式下，北斗衛(wèi)星相位中心改正對鐘差影響可達(dá)0.9 ns，綠線代表IGSO 衛(wèi)星，藍(lán)線代表MEO 衛(wèi)星。

圖5 不同姿態(tài)下相位中心誤差對鐘差的影響

此外，相位中心的標(biāo)定值也是影響時間同步結(jié)果的一個因素。目前，北斗衛(wèi)星天線相位中心有三種來源，分別是衛(wèi)星出廠標(biāo)定值、國際多模GNSS (global navigation satellite system)實驗工程(Multi-GNSS Experiment and Pilot Project,MGEX)[16]、ESA (European Space Agency)解算值[17]。出廠標(biāo)定值與MGEX 給出的值基本一致。ESA 解算的衛(wèi)星天線相位中心與出廠標(biāo)定值存在較大差異，X方向誤差可達(dá)0.5 m，Z方向達(dá)到2 m。

針對以上情況，我們對相位中心X和Z方向存在的誤差進(jìn)行了仿真，計算了不同量級的相位中心誤差對雙向鐘差的影響。由圖6可知，當(dāng)Z方向存在1 m 誤差時，對鐘差影響在0.02 ns 以內(nèi)；存在2 m 誤差時，對鐘差影響在0.04 ns 以內(nèi)。因此Z方向誤差對鐘差的影響較小。當(dāng)X方向存在0.1 m 的誤差時，對鐘差影響為0.02 ns；當(dāng)X方向存在0.5 m 誤差時X方向誤差可以達(dá)到0.08 ns，說明相位中心X方向的誤差對雙向鐘差精度影響較顯著。

圖6 不同來源衛(wèi)星相位中心對雙向鐘差影響

3 廣播鐘差參數(shù)精度分析

BDS 直接利用雙向鐘差序列建模生成廣播鐘差參數(shù)，然后上注到衛(wèi)星。我們在鐘差參數(shù)的擬合中，直接利用了短時間雙向鐘差點進(jìn)行a0 和a1 兩個鐘差參數(shù)的線性擬合，把a(bǔ)2 項設(shè)為0。衛(wèi)星位于境內(nèi)時，廣播鐘差參數(shù)每小時更新一次；出境后廣播鐘差參數(shù)不更新；重新入境后，由于預(yù)報時間長，預(yù)報誤差會比較大。我們采用雙向?qū)崟r觀測的鐘差對廣播鐘差參數(shù)的預(yù)報精度進(jìn)行評估。表1統(tǒng)計了2016年1月至6月廣播鐘差參數(shù)預(yù)報精度(RMS，包括境內(nèi)與重新入境時段)。從表中可以看出，除2 號星和11 號星外，其余衛(wèi)星的鐘差預(yù)報精度均在2 ns 以內(nèi)。其中GEO 衛(wèi)星1,3,4 號星的精度在1 ns 以內(nèi)。GEO2 由于星鐘信號減弱造成衛(wèi)星鐘性能較低，因此參數(shù)精度較低。整體來看，GEO 衛(wèi)星廣播鐘差預(yù)報精度高于IGSO 和MEO 衛(wèi)星。

表1 北斗衛(wèi)星鐘差預(yù)報精度(RMS)統(tǒng)計表 ns

為了更加細(xì)致分析GEO 衛(wèi)星的廣播鐘差預(yù)報精度，將1 號星的廣播鐘差預(yù)報誤差畫出。圖7a)表示GEO 1 號星鐘差的預(yù)報誤差?？梢钥闯鲈趨?shù)更新的1 h 內(nèi)，預(yù)報誤差逐漸增大；當(dāng)參數(shù)更新時，預(yù)報誤差回到0 附近。預(yù)報誤差成“鋸齒狀”，說明廣播鐘差參數(shù)a1 項存在偏差。同時發(fā)現(xiàn)一天中存在兩處預(yù)報誤差明顯偏大的地方。為了進(jìn)一步發(fā)現(xiàn)這兩處預(yù)報誤差較大的原因，對該天的雙向鐘差數(shù)據(jù)進(jìn)行二次擬合，殘差如圖7b)所示。鐘差存在類似軌道周期的波動[10,18]，預(yù)報誤差較大的地方正是鐘差斜率發(fā)生變化的拐點處。

圖7 廣播鐘差預(yù)報誤差和GEO 單天擬合殘差

北斗IGSO/MEO 存在出入境的情況，IGSO 出境時間約6 h，MEO 出境時間約12 h，因此還需采用長期數(shù)據(jù)二次擬合參數(shù)，給出a0,a1,a2。由于每顆衛(wèi)星星鐘性能不同，擬合采用的數(shù)據(jù)長度也會不同。圖8給出IGSO 和MEO 衛(wèi)星的預(yù)報誤差，其中，紅色表示主控站觀測的衛(wèi)星鐘差，綠色和藍(lán)色分別表示另外兩個外場站觀測的衛(wèi)星鐘差。衛(wèi)星剛?cè)刖硶r，誤差明顯變大。IGSO 的預(yù)報誤差達(dá)到10 ns，MEO 預(yù)報誤差達(dá)到8 ns。衛(wèi)星出境時，由于站間歸算的原因，鐘差噪聲會放大，預(yù)報誤差也會變大。

圖8 廣播鐘差預(yù)報誤差

4 結(jié) 論

中國北斗衛(wèi)星導(dǎo)航采用星地雙向時頻傳遞的方法直接測量衛(wèi)星鐘相對于地面主鐘BDT的差，并用于導(dǎo)航電文鐘差參數(shù)的預(yù)報建模。與定軌法不同，雙向法能夠得到頻率為1 Hz 高頻實測鐘差點，且鐘差與軌道沒有相關(guān)性。本文得到的主要結(jié)論如下。

(1)電離層延遲誤差是影響北斗雙向時間同步精度的主要誤差源。雙頻相位平滑偽距法求得電離層延遲對雙向鐘差影響在0.1 ns 以內(nèi)，當(dāng)衛(wèi)星剛?cè)刖硶r雙向鐘差精度較低，在利用雙向鐘差進(jìn)行廣播鐘差參數(shù)擬合時，要將此段數(shù)據(jù)剔除。

(2)電離層模型無法表達(dá)電離層變化的高頻部分，其精度無法滿足雙向鐘差的精度要求。

(3)北斗衛(wèi)星存在動偏與零偏兩種姿態(tài)模式，不同姿態(tài)下相位中心對造成鐘差誤差達(dá)到0.9 ns。不同來源的相位中心值，對雙向鐘差誤差較小。

(4)實測數(shù)據(jù)分析結(jié)果表明，星地雙向衛(wèi)星鐘差內(nèi)符合精度(RMS)優(yōu)于0.15 ns。

(5)北斗廣播鐘差參數(shù)預(yù)報精度在2 ns 以內(nèi)，GEO 衛(wèi)星廣播鐘差預(yù)報精度高于IGSO 與MEO 衛(wèi)星。IGSO/MEO 衛(wèi)星出入境時，預(yù)報誤差較大，可達(dá)10 ns。