線性插值在改善GPS 共視時間比對性能中的應用

馬慧穎 ,王威雄 ,武文俊,3

(1.中國科學院國家授時中心,西安 710600;2.中國科學院大學,北京 100049;3.中國科學院大學天文與空間科學學院,北京 100049)

0 引言

時間比對在導航定位、衛星發射、天文觀測、戰場調度等軍事、民生各領域都發揮著重要作用。當前的國際標準時間是協調世界時(UTC),它是將全球八十多個守時實驗室通過時間比對鏈路獲得的數據進行加權平均等一系列處理后得到的,因此時間比對也是計算協調世界時的重要一環[1]。

自二十世紀80 年代GPS 共視時間比對方法被提出后至今,GPS 共視時間比對方法以其覆蓋范圍廣,比對精密度高,設備便宜,操作簡單可靠的特點得到了國內外學者的廣泛研究[2],現已成為最常用的遠程時間比對技術之一。但使用共視時間比對技術的前提是兩個守時實驗室可以同時接收來自同一顆或多顆衛星的時間信號,因此該方法的比對精密度會受到基線長度的限制[3],而當基線過長時,甚至會出現無法進行共視時間比對的情況。

為了解決共視時間比對受基線限制的問題,國際計量局(BIPM)的江志恒等人提出了全視時間(AV)比對技術[4],但是此方法需要用到國際全球導航衛星系統服務組織(IGS)產品改善,計算方式較為復雜。雖然2019 年Verhasselt 等人提出通過多系統融合的方式提高共視時間比對精密度[5],但依然無法解決共視時間比對受到基線長度的限制。

本研究基于星地鐘差具有良好的線性特征,通過將鐘差估計值應用到共視時間比對的計算過程中以實現改善GPS 共視時間比對性能的目標。首先強調了改善長基線共視時間比對精密度的重要性,其次詳細介紹了基于線性插值的共視時間比對方法(Common-View time comparison based on Linear Interpolation,iCV)的原理,并給出了相應的試驗設計,然后對試驗結果進行了分析,最后對試驗結果的總結,并對下一步研究進行了展望。

1 基于線性插值的共視時間比對原理及試驗設計

1.1 基于線性插值的共視時間比對原理

眾所周知,搬運鐘時間比對方法是一種最簡單、直接的時間比對方法,此方法通常是閉環搬運,搬運方法如圖1 所示,首先在A 位置把搬運鐘與A本地鐘進行比對獲得鐘差,然后將搬運鐘搬運到B位置,獲得搬運鐘與B 本地鐘的鐘差,最后將搬運鐘從B 位置搬運回A 位置,再次獲得A 本地鐘與搬運鐘的鐘差。由于本地鐘和搬運鐘穩定性都較好,在整個往返過程中可以使用插值方法來估計搬運鐘的時差偏移量[6],進而間接得到兩個本地鐘之間的鐘差結果。

圖1 搬運鐘時間比對原理圖Fig.1 The principle of portable clock time comparison

在GPS 共視時間比對中,由于每個GPS 衛星在一個儒略日后通常會回到天空中同一位置,因此衛星鐘可以被視為搬運鐘,如圖2 所示,以GPS 系統PRN01、PRN15 衛星分別在KRISS、NTSC 站點上空時的星地鐘差為例,在電離層和對流層延時未被消除的情況下,GPS 系統的每個衛星與地面鐘的鐘差的行為是較難預測的,但是在消除電離層和對流層延時后,如圖3 所示,可以看出GPS 系統的每個衛星與地面鐘的鐘差都具有良好的線性行為,且KRISS 和NTSC 兩站點的時間系統都以穩定性良好的氫鐘作為主鐘,因此可以使用線性插值來估計GPS 衛星在非觀測周期內理想狀態下的星地鐘差。

圖2 電離層和對流層時延未被消除的情況下的星地鐘差圖Fig.2 The time difference between the local system and the satellite clock when the tropospheric refraction delay and ionospheric delay are not eliminated

圖3 電離層和對流層延時被消除的情況下的星地鐘差圖Fig.3 The time difference between the local system and the satellite clock when the tropospheric refraction delay and ionospheric delay are eliminated

GPS 系統中PRN01 衛星在A、B 兩個站點的觀測周期如圖4 所示,此處假設A、B 兩個站點無法在任一時刻同時接收到PRN01 衛星的信號。在第一個儒略日時間內,A 站點在t1所在觀測周期(一般情況下每個觀測周期時長為2～6 h,假設觀測周期時長均為τ)內接收到PRN01 的時間信號,并將本地鐘(TA)與衛星鐘(Tref)的鐘差記錄為數據集1:(TA-Tref)01(t1)。同理可獲得B 站點在t2所在觀測周期內的數據集2:(TB-Tref)01(t2),(TB代表B站點的本地鐘數據)。

圖4 PRN01 衛星在不同站點的觀測時間與觀測數據圖Fig.4 Observation time and observation data of the PRN01 satellite at different stations

到第二個儒略日時,PRN01 衛星再次出現在站點A 上空同一位置時,即在t3所在觀測周期內時,將本地鐘(TA)與衛星鐘(Tref)之間的鐘差記錄為數據集3:(TA-Tref)01(t3)。

基于星地鐘差具有良好的線性特征,可以通過數據集1 和數據集3 來估計在t2所在時間段內A站點本地鐘和PRN01 衛星鐘的鐘差。

A 站點和B 站點在t2所在觀測周期內的鐘差可以用式(1)表示

式中:I——估計值;R——觀測值。

通過B 站點在t2和t4所在觀測周期內的數據計算出在t3時間段內的估計值后,可以通過式(2)計算在t3所在時間段內A 站點和B 站點的鐘差。

重復此過程,可以生成PRN01 衛星在t1～t4內A 站點和B 站點的鐘差數據。

理論上通過此方法分別對剩余31 顆GPS 衛星與兩站點的本地鐘鐘差進行分析,可以得出在使用不同衛星鐘情況下A、B 兩站點的鐘差數據,然而在實際中,每個站點可以在同一時刻觀測到大約4～12個衛星,因此實際只需計算對應站點可以觀測到的衛星的數據即可。

將通過以上計算得到的PRN01 衛星在不同時間段內A、B 兩站點的數據進行綜合,可以得到四種類型的鐘差數據組合,如圖5 所示,其中上標I-R表示A 站點估計值和B 站點實際觀測值的差,上標R-I表示A 站點實際觀測值和B 站點估計值的差,I-I表示A 站點估計值和B 站點估計值的差,R-R表示A、B 兩站點實際觀測值的差,即共視時間比對數據,假設A、B 兩站點在此段時間內無法直接進行共視時間比對,因此此時R-R沒有數據。

圖5 四種進行共視時間比對計算的方法圖Fig.5 Four methods of comparing time for CV

將所有衛星數據分別處理后,需對每一時刻所有衛星的數據進行加權平均,以獲得每個時刻的共視時間比對結果,如式(3)所示

式中:IR-IR——使用線性插值方法后得到的最終共視時間比對結果;ti——不同時間;n——衛星的PRN 編號;wn——對應衛星的權重;IR0 -IR0——通過使用四種方法中的其中一種。

1.2 試驗設計及數據預處理

全球衛星導航時間傳遞標準文件(CGGTTS)中的數據以16 min 為一個周期,其中前2 min 用于捕獲衛星,中間13 min 采集數據,最后1 min 處理數據,該數據星地鐘差參數(REFSYS)是在獲得的偽距觀測值中有效扣除星地幾何距離及電離層、對流層等各項誤差所引起的測距時延后得到的,代表本地系統時間和衛星系統時間的差[7,8],因此使用CGGTTS 文件中的REFSYS 數據進行計算。

本研究數據處理過程中主要分為3 個階段:

第一個階段為前期準備階段,即對觀測數據進行分析,剔除異常觀測數據,將處理過的觀測數據根據觀測周期計算平均值,為進行線性插值做準備。

如圖6 所示,以GPS 系統PRN01 衛星在KRISS站點不同高度角下的觀測數據為例,展示了不同高度角閾值下觀測數據的分布情況,其中星號、菱形、方形依次代表保留高度角在15°、30°、45°以上時的觀測數據,可以看出,當保留高度角在15°以上時的觀測數據時,仍存在較多異常值,而保留高度角在30°或45°以上時的觀測數據時均能有效剔除異常值,但保留高度角在30°以上時的觀測數據量明顯多于保留高度角在45°以上時的數據量,因此,選擇保留高度角在30°以上的數據進行分析。

第二階段是基于各站點本地時間與衛星系統時間線性、可預測的特點對觀測數據進行插值計算。由于PPP 數據的時間間隔是300 s,為了使計算結果能夠更好地與PPP 數據進行比對,因此以300 s 的時間間隔對觀測數據通過線性插值的方式進行采樣。以KRISS 上空的PRN01 衛星為例,給出了在觀測周期內線性插值的結果,如圖7 所示。

圖7 計算觀測周期內的估計值Fig.7 Calculate estimates over the observation period

圖7 中,圓形代表觀測數據,菱形代表線性抽樣估計值,其中估計值通過使用兩個相鄰最近的觀測數據進行線性插值計算獲得,具體如式(4)

式中:tu——通過式(4)獲得的兩個觀測點之間以300 s 為時間間隔的估計值;ti,ti+1——相隔為960 s(即16 min)的兩個觀測數據對應的時刻。

GPS 系統中某些衛星一個儒略日內僅在站點上空出現一次,如圖8 所示,以GPS 系統中每個儒略日在KRISS 上空出現一次的PRN01 衛星為例,圓形代表觀測數據,將每段觀測周期內的平均值作為觀測周期中點處的標準值,用菱形表示,然后在兩個標準值之間使用線性插值方法計算兩個相鄰觀測周期間的鐘差估計值,用圓點表示。

圖8 KRISS 上空出現的PRN01 衛星數據圖Fig.8 PRN01 satellite data appearing over KRISS station

另外還有一些GPS 衛星在同一個儒略日內出現在站點上空2 次或以上,如圖9 所示,以每天在KRISS 上空出現兩次的PRN29 衛星為例,圓形代表原始觀測數據,相鄰兩個觀測周期大約間隔0.5 天,每個觀測周期的平均值用菱形表示,圓點代表在兩個相鄰平均值之間使用線性插值方法而得到的估計值。此時由于每日周期性變化的影響,得到的插值結果相比實際結果誤差較大。因此,在進行插值計算時也以一個儒略日為周期,根據衛星在一個儒略日內出現的次序分別進行線性插值計算,最后將相同時間的插值結果進行平均作為最終結果。如圖10 所示,根據PRN29 將平均值分為兩個數據集H、L,對應數據標記為H1～H4,L1～L4,標記為菱形,兩個子集中每兩個相鄰數據(例如L2 和L3,H1 和H2)相差約為1 個儒略日。分別計算兩個子集的插值結果,然后將插值結果進行平均得到最終結果,用叉號表示。

圖9 KRISS 上空的PRN29 數據圖Fig.9 PRN29 satellite data appearing over KRISS station

圖10 PRN29 估計值計算圖Fig.10 The estimate of PRN29 is calculated

重復以上過程,獲取每顆GPS 衛星的鐘差數據估計值后,將相同時刻同一衛星的鐘差數據進行共視計算,最后將所有可用GPS 衛星的共視結果進行等權加權平均,得到最終的iCV 結果。

最后一個階段為驗證階段,將得到的結果與傳統共視時間比對結果(CV)、PPP 比對結果進行比較,驗證方法的可行性與準確性。

2 試驗結果及分析

本試驗中使用MJD 60035 至60064 連續30 天的數據,分別計算NTSC-PTB、KRISS-PTB 兩條長基線鏈路的時間比對數據,為對試驗結果進行驗證,給出了相同鏈路的傳統CV 與PPP 時間比對結果,如圖11 所示,PPP 時間比對結果標記為加號,iCV 結果標記為圓形,CV 標記為星號。由于PPP 時間比對采用載波相位數據其準確度較高,因此在本試驗中使用PPP 時間比對為參考來評估線性插值比對結果。

圖11 不同鏈路的iCV 結果圖Fig.11 The results of the iCV on different baselines

從圖11 可知,三種時間比對方法結果一致性較好,CV 結果的短期噪聲較大,iCV 結果的噪聲水平明顯小于傳統CV 結果。

兩條鏈路iCV 結果減去PPP 時間比對結果獲得的雙差結果(iCV-PPP)如圖12 所示。NTSC-PTB 鏈路iCV-PPP 的平均值為0.005 4 ns,RMS 為0.53 ns;KRISS-PTB 鏈路iCV-PPP 的平均值為0.035 ns,RMS 為0.8 ns;兩條鏈路iCV 與PPP 結果差值的均值均接近0,說明iCV 結果與PPP 結果一致性較好。

圖12 不同鏈路iCV 結果與PPP 時間比對結果差圖Fig.12 The results of the iCV minus the PPP time alignment result on different links

分別統計了兩條鏈路iCV-PPP 與CV-PPP 的RMS 結果,如表1 所示。

表1 不同基線長度下CV、iCV 的RMS表Tab.1 The RMS of CV and iCV at different baseline lengths

從表1 數據顯示,NTSC-PTB 和KRISS-PTB 這兩條鏈路通過傳統CV 得到的RMS 分別為1.03 ns和1.74 ns,而通過iCV 方法得到的RMS 分別為0.53 ns 和0.8 ns,與傳統CV 相比,RMS 分別提升了48%和54%。這一結果表明線性插值能夠提高時間比對結果的準確性。

為對穩定度提升進行評估,給出了iCV 結果與PPP、傳統CV 時間比對的穩定度對比結果,如圖13所示。

圖13 不同鏈路使用不同時間比對方法時的時間穩定度結果圖Fig.13 Time stability results for different links using different time comparison methods

從圖13 中可以看出:平均時間小于3 天時iCV得到的時間穩定性優于傳統CV 結果,PPP 結果的時間穩定性最好。當平均時間大于3 天時,iCV 結果的時間穩定性和傳統CV、PPP 結果接近。iCV 結果的時間穩定性一直優于傳統CV 的時間穩定性,說明線性插值能夠提高共視時間比對結果的穩定性。綜上所述,長基線情況下,iCV 與傳統CV 相比,在穩定性和精密度兩方面有所提升。

3 結束語

為提高共視時間比對的精密度,實現了一種基于線性插值的共視時間比對方法(iCV)。該方法使用標準CGGTTS 文件,計算得到NTSC-PTB、KRISSPTB 兩條長基線的iCV 結果,并與GPS 精密單點定位時間比對結果(PPP)對比,其互差結果的RMS 分別為0.53 ns、0.8 ns,與傳統共視時間比對(CV)結果相比RMS 結果提升了48%、54%,同時對兩條鏈路時間比對結果的時間穩定性進行了分析,結果表明,iCV 結果的時間穩定性一直優于CV 結果的時間穩定性,說明該方法可以從穩定度和精密度兩方面改善長基線共視時間比對結果的性能。下一步將繼續完善基于線性插值的共視時間比對方法,并將該方法應用于格洛納斯、北斗、伽利略多系統融合的時間比對。