章潔君，陳俊平,2，張益澤,3

(1. 中國科學(xué)院上海天文臺(tái)，上海 200030； 2. 中國科學(xué)院上海天文臺(tái)上海市空間導(dǎo)航與定位技術(shù)重點(diǎn) 實(shí)驗(yàn)室，上海 200030； 3. 同濟(jì)大學(xué)測繪與地理信息學(xué)院，上海 200092)

GNSS空間信號(hào)法時(shí)差監(jiān)測方法與結(jié)果分析

章潔君1，陳俊平1,2，張益澤1,3

(1. 中國科學(xué)院上海天文臺(tái)，上海 200030； 2. 中國科學(xué)院上海天文臺(tái)上海市空間導(dǎo)航與定位技術(shù)重點(diǎn) 實(shí)驗(yàn)室，上海 200030； 3. 同濟(jì)大學(xué)測繪與地理信息學(xué)院，上海 200092)

比較了3種基于GNSS空間信號(hào)進(jìn)行時(shí)差監(jiān)測的方法，其中，單站精密單點(diǎn)定位(PPP)時(shí)差監(jiān)測法依賴于衛(wèi)星軌道和鐘差的精度，高精度PPP一般采用事后處理的方法；網(wǎng)解時(shí)差監(jiān)測法受衛(wèi)星軌道誤差的影響較小，且不受衛(wèi)星鐘差誤差的影響，因此能夠用于高精度時(shí)差實(shí)時(shí)監(jiān)測。為分析不同方法時(shí)差監(jiān)測的差異，在距離不同的時(shí)間實(shí)驗(yàn)室之間進(jìn)行了時(shí)差試驗(yàn)，試驗(yàn)中將外接相同原子鐘的兩個(gè)測站的時(shí)差用于評(píng)定時(shí)差監(jiān)測的精度，而相距2000 km的測站用于驗(yàn)證高精度實(shí)時(shí)時(shí)差監(jiān)測的精度。結(jié)果表明：①實(shí)時(shí)網(wǎng)解與事后PPP時(shí)差監(jiān)測的精度相當(dāng)，達(dá)到了0.16 ns；②采用PPP方法，是否固定站坐標(biāo)時(shí)差結(jié)果差異的RMS為0.04 ns；③采用實(shí)時(shí)網(wǎng)解的方法，對(duì)于相隔數(shù)千千米的測站實(shí)時(shí)時(shí)差監(jiān)測，其與PPP后處理結(jié)果的差異約為0.14 ns。

時(shí)差監(jiān)測；空間信號(hào)法；實(shí)時(shí)網(wǎng)解；PPP

傳統(tǒng)的時(shí)差監(jiān)測方法主要有衛(wèi)星雙向時(shí)間頻率傳遞法(TWSTFT)和衛(wèi)星共視法等方法，其中，TWSTFT利用專用的時(shí)間同步鏈路，系統(tǒng)復(fù)雜且經(jīng)濟(jì)成本較高，衛(wèi)星共視方法是基于對(duì)GNSS衛(wèi)星的偽距觀測值，其精度較差[1-2]。結(jié)合GNSS數(shù)據(jù)處理技術(shù)的發(fā)展，近年來基于導(dǎo)航系統(tǒng)空間信號(hào)的時(shí)差監(jiān)測技術(shù)得到了較快發(fā)展，提高了時(shí)差監(jiān)測的實(shí)時(shí)性及監(jiān)測精度[3-5]。空間信號(hào)法時(shí)差監(jiān)測主要有兩種方法：PPP技術(shù)和多站網(wǎng)解技術(shù)[6]。單站精密單點(diǎn)定位技術(shù)在求取測站坐標(biāo)的同時(shí)求取測站鐘差參數(shù)；它要求固定衛(wèi)星軌道和鐘差，因此高精度的PPP時(shí)差測量一般采用后處理的方式[7-8]。多站網(wǎng)解技術(shù)固定了參考站的坐標(biāo)及鐘差參數(shù)，求解其他站的相對(duì)鐘差，同時(shí)衛(wèi)星鐘差也作為參數(shù)進(jìn)行求解；網(wǎng)解技術(shù)中衛(wèi)星軌道誤差對(duì)網(wǎng)解數(shù)據(jù)處理影響較小，并且其影響最大的徑向部分會(huì)被衛(wèi)星鐘差參數(shù)吸收，因此網(wǎng)解技術(shù)能夠?qū)崿F(xiàn)高精度的實(shí)時(shí)時(shí)差監(jiān)測[9]。

本文擬比較分析基于PPP及網(wǎng)解技術(shù)時(shí)差監(jiān)測的模型和數(shù)據(jù)處理方法，設(shè)計(jì)時(shí)間實(shí)驗(yàn)室的時(shí)差監(jiān)測試驗(yàn)，通過時(shí)差參數(shù)的精度及穩(wěn)定性情況評(píng)定這兩種方法時(shí)差監(jiān)測的精度，驗(yàn)證基于網(wǎng)解技術(shù)長距離實(shí)時(shí)時(shí)差監(jiān)測的可行性。

1 GNSS空間信號(hào)法時(shí)差監(jiān)測原理

GNSS空間信號(hào)法時(shí)差監(jiān)測方法一般采用非差無電離層組合模型，其偽距、相位觀測方程如下

(1)

時(shí)差解算中，衛(wèi)星位置一般可采用精密軌道進(jìn)行固定，測站坐標(biāo)可以作為參數(shù)估計(jì)或固定為精確已知值；接收機(jī)鐘差、模糊度等參數(shù)一般作為參數(shù)估計(jì)。采用PPP進(jìn)行時(shí)差監(jiān)測需要固定衛(wèi)星的鐘差，因此這種方法在實(shí)時(shí)處理中精度受到實(shí)時(shí)鐘差精度的約束，基于廣播星歷的鐘差精度能夠達(dá)到大約2～5 ns[11]。網(wǎng)解技術(shù)則將衛(wèi)星鐘差也作為待求參數(shù)，因此其精度不受輸入衛(wèi)星鐘差精度的影響，能夠?qū)崿F(xiàn)高精度的實(shí)時(shí)解算[12]。網(wǎng)解技術(shù)中測站和衛(wèi)星鐘差參數(shù)的同時(shí)解算存在秩虧，需要定義一個(gè)鐘差基準(zhǔn)，從而解算的鐘差參數(shù)都是相對(duì)該基準(zhǔn)的相對(duì)鐘差[13]。

以上的數(shù)據(jù)處理都是基于相位觀測，單站的相位處理由于模糊度無法固定，一般存在一定的收斂時(shí)間?；诰W(wǎng)解的相位處理可進(jìn)行雙差模糊度固定，因此其實(shí)時(shí)收斂的時(shí)間更短。

獲取的鐘差參數(shù)包含了接收機(jī)設(shè)備時(shí)延，對(duì)于定時(shí)接收機(jī)該時(shí)延值可在實(shí)驗(yàn)室進(jìn)行高精度標(biāo)定。對(duì)于外接時(shí)頻信號(hào)T(A)的GNSS接收機(jī)A，以上鐘差結(jié)果扣除設(shè)備時(shí)延后，得到的測站鐘差為外接時(shí)間信號(hào)T(A)與衛(wèi)星系統(tǒng)時(shí)間GNSST之差[14]

δA=GNSST-T(A)

(2)

不同接收機(jī)鐘差作差得到外接時(shí)頻信號(hào)的差異，即為時(shí)差測量值

ΔT=δA-δB

=GNSST-T(A)-GNSST-T(B)

=T(B)-T(A)

(3)

2 計(jì)算結(jié)果及分析

對(duì)GNSS信號(hào)法時(shí)差監(jiān)測方法進(jìn)行測試，設(shè)計(jì)了驗(yàn)證試驗(yàn)。試驗(yàn)中在成都某時(shí)頻實(shí)驗(yàn)室實(shí)驗(yàn)場架設(shè)5臺(tái)GNSS接收機(jī)，其中1臺(tái)接收機(jī)REF連接了該實(shí)驗(yàn)室基準(zhǔn)時(shí)頻信號(hào)，另外4臺(tái)(A、B、C、D)是為了此次測試在其附近臨時(shí)架設(shè)的，此外將距離測試場地約為2000 km的測站E也加入處理分析。以上測站中A、B站外接至同一臺(tái)原子鐘。試驗(yàn)測量了2015年第228和229年積日的數(shù)據(jù)，并分別使用事后PPP和實(shí)時(shí)網(wǎng)解技術(shù)獲取測站鐘差。實(shí)時(shí)網(wǎng)解中以實(shí)驗(yàn)室原子鐘作為基準(zhǔn)，求取其他測站相對(duì)于REF站的時(shí)差。

2.1 時(shí)差結(jié)果分析

由于PPP時(shí)差測量需要收斂，因此結(jié)果比較中不包含剛開始觀測的1小時(shí)數(shù)據(jù)。由于篇幅限制，僅圖1給出網(wǎng)解技術(shù)所得的測站鐘差(扣除了線性項(xiàng))的結(jié)果，圖2為各測站相對(duì)于A站的時(shí)差計(jì)算值，表1在圖1的基礎(chǔ)上統(tǒng)計(jì)了擬合殘差，作為時(shí)差測量的精度。

圖1 網(wǎng)解時(shí)差監(jiān)測測站鐘差序列

分析圖1和圖2可知，測站鐘差解算兩天的吻合度較好。C站未連接氫原子鐘，其本身噪聲較大，其他3個(gè)測站的鐘差值均表現(xiàn)出較小的噪聲。從表1可以看到，對(duì)于外接了原子鐘的測站，3種方法的測站鐘差精度都在0.2 ns以內(nèi)，平均為0.16 ns。對(duì)于A與B站的時(shí)差，其頻率差異為-1.28×10-15，由于A、B外接了相同的時(shí)頻信號(hào)，因此該頻率的差異主要反映了接收機(jī)的噪聲及性能。而C站3種技術(shù)求得的測站鐘差RMS均大于0.3 ns，由于未接時(shí)頻信號(hào)，其時(shí)差噪聲較大。

表1 3種時(shí)差監(jiān)測方法測站鐘差結(jié)果RMS值比較 ns

圖2 網(wǎng)解測站時(shí)差序列

2.2 實(shí)時(shí)網(wǎng)解與PPP事后處理時(shí)差結(jié)果的比較

將網(wǎng)解方法與事后PPP方法獲取的時(shí)差測量值作比較，如圖3所示。從圖3中可看到，網(wǎng)解結(jié)果與PPP結(jié)果吻合較好，其差異的RMS為0.04 ns。PPP處理中對(duì)于坐標(biāo)已知的測站可以進(jìn)行坐標(biāo)的固定，將其與常規(guī)PPP的結(jié)果進(jìn)行比較，圖4列出了不同測站這兩種方法的差異，其結(jié)果基本吻合，差值的RMS約為0.04 ns。因此，是否固定測站坐標(biāo)對(duì)解算結(jié)果影響不大。

網(wǎng)解方法把衛(wèi)星鐘差與測站鐘差一起作為參數(shù)估計(jì)，可實(shí)現(xiàn)長距離實(shí)時(shí)的時(shí)差觀測。圖5比較了距離達(dá)2000 km的E站與A站的實(shí)時(shí)時(shí)差測量值與PPP事后處理的結(jié)果。兩者結(jié)果基本吻合，且兩天的結(jié)果一致，差值的RMS均值為0.14 ns。因此可以認(rèn)為，對(duì)于長距離的測站，實(shí)時(shí)網(wǎng)解可與事后PPP解算精度一致。

2.3 實(shí)時(shí)網(wǎng)解與PPP事后處理時(shí)差頻率穩(wěn)定度的比較

在時(shí)間領(lǐng)域中，使用Allan方差來評(píng)估利用GNSS方法得到的產(chǎn)品的時(shí)域穩(wěn)定度，在實(shí)際測量中，Allan方差是描述隨機(jī)噪聲時(shí)域統(tǒng)計(jì)特性的通用標(biāo)準(zhǔn)[15]。其公式為

圖3 網(wǎng)解時(shí)差監(jiān)測與PPP時(shí)差監(jiān)測的比較

圖4 PPP是否固定站坐標(biāo)對(duì)時(shí)差監(jiān)測的影響

(4)

式中，xi為歷元i的觀測值；N為觀測序列的長度；τ為采樣間隔，τ=nτ0，τ0為原始觀測序列采樣間隔；n為正整數(shù)(1≤n≤N)。本文中，xi是歷元i的鐘差解，N為鐘差解序列樣本總數(shù)(即歷元數(shù))，τ0為觀測值采樣間隔30 s。

圖5 長距離實(shí)時(shí)網(wǎng)解與PPP事后處理時(shí)差結(jié)果的比較

基于以上數(shù)據(jù)處理的結(jié)果，分析了不同采樣率下的Allan方差，如圖6所示。

圖6 Allan方差時(shí)差序列

結(jié)合圖6可知，3種方法計(jì)算得到的Allan方差統(tǒng)計(jì)結(jié)果相當(dāng)，說明這3種時(shí)差解算方法的頻率穩(wěn)定度一致。其中10 000 s間隔的頻率穩(wěn)定度在1×10-13～1.5×10-13左右。

3 結(jié) 論

針對(duì)以上關(guān)于3種時(shí)差監(jiān)測方法的研究和驗(yàn)證，得出以下結(jié)論：

(1) PPP、PPP固定和實(shí)時(shí)網(wǎng)解時(shí)差監(jiān)測的結(jié)果均能夠達(dá)到優(yōu)于0.2 ns的精度，頻率萬秒的穩(wěn)定度為1×10-13～1.5×10-13。

(2) 實(shí)時(shí)網(wǎng)解與PPP事后處理時(shí)差監(jiān)測的精度相當(dāng)，都優(yōu)于0.2 ns。距離達(dá)2000 km的實(shí)時(shí)網(wǎng)解時(shí)差測量精度與PPP差異為0.14 ns。

(3) 采用PPP方法，是否固定站坐標(biāo)造成時(shí)差監(jiān)測結(jié)果差異的RMS為0.04 ns。

[1] 江志恒. GPS全視法時(shí)間傳遞回顧與展望[J]. 宇航計(jì)測技術(shù),2007(S1):53-71.

[2] SENIOR K,KOPPANG P,RAY J.Developing an IGSTime Scale[J].IEEE Transactions on UltrasonicFerroelectric Frequency Control,2003,50(6):585-593.

[3] KOUBA J,SPRINGER T.New IGS Station and Satellite Clock Combination[J].GPS Solutions,2001，4(4):31-36.

[4] SENIOR K.New IGS Timescaleamp;New UT1/LODSeries[C]∥Analysis Center Workshop.Florida：[s.n.]，2008.

[5] 周鈺，熊耀恒.利用衛(wèi)星進(jìn)行雙向時(shí)間傳遞[J]. 天文研究與技術(shù).國家天文臺(tái)臺(tái)刊,2006,3(1):28-34.

[6] ORGIAZZI D,TAVELLA P,LAHAYE F.Experimental Assessment of the Time Transfer Capability of Precise Point Positioning(PPP)[C]∥IEEE International Frequency Control Symposiumand amp; Exposition.Vancouver,BC,Canada：IEEE,2005.

[7] 楊文輝，黃新寶，羅銳.衛(wèi)星導(dǎo)航時(shí)間傳遞性能的分析與比較[J].現(xiàn)代導(dǎo)航,2013(1):14-18.

[8] HEERING T.Impact of Absolute Phase Center Models on GPS Reference Frames[C]∥European Geosciences Union General Assembly,Massachusetts.Cambridge:[s.n.],2007.

[9] 李浩軍,王解先,陳俊平，等. 基于GNSS網(wǎng)絡(luò)的實(shí)時(shí)精密單點(diǎn)定位及精度分析[J]. 地球物理學(xué)報(bào),2010,53(6)：1302-1307.

[10] 劉經(jīng)南,葉世榕. GPS非差相位精密單點(diǎn)定位技術(shù)探討[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2002,27(3):234-240.

[11] 張小紅,程世來,李星星,等.單站GPS載波平滑偽距精密授時(shí)研究[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2009,34(4):463-465.

[12] 劉曉剛,吳曉平,張傳定. 衛(wèi)星雙向共視法時(shí)間比對(duì)計(jì)算模型及其精度評(píng)估[J]. 測繪學(xué)報(bào),2009,38(5):415-421.

[13] 陳俊平,張益澤,謝益炳,等. 超大觀測網(wǎng)絡(luò)及多GNSS系統(tǒng)的快速數(shù)據(jù)處理[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2014,39(3):253-257.

[14] 辜聲峰. 多頻GNSS非差非組合精密數(shù)據(jù)處理理論及其應(yīng)用[D].武漢：武漢大學(xué),2013.

[15] 張小紅,蔡詩響,李星星,等.利用GPS精密單點(diǎn)定位進(jìn)行時(shí)間傳遞精度分析[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2010,35(3):274-278.

ResearchandResultsAnalysisofMethodsofGNSSSpatialSignalTimeMonitoring

ZHANG Jiejun1,CHEN Junping1,2,ZHANG Yize1,3

(1. Shanghai Astronomical Observatory,Chinese Academy of Sciences, Shanghai 200030, China; 2. Shanghai Key Laboratory of Space Navigation and Positioning Technology,ShanghaiAstronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China; 3. College of Surveying and Geo-informatics, Tongji University, Shanghai 200092, China)

Methods using GNSS signal in space to monitor time offset among timing stations are dicussed. Three methods namely the traditional PPP, PPP with fixed station coordinates and the real-time network solutions are compared. To validate and assess the performance of SIS time offset monitoring, data of 5 GNSS stations with different configurations are analyzed.Results show that: (1)Precision of real time network solutions is 0.16 ns, which is at the same level of post-processed PPP.(2)Station coordinates estimation has impact on the determination of station clocks，which is about 0.04 ns. (3) Precision of time transfer for stations depart 2000 km could reach 0.14 ns in the real-time network solution.

time monitoring； spatial signal; real-time network solution;PPP

P224

A

0494-0911(2017)01-0026-04

章潔君，陳俊平，張益澤.GNSS空間信號(hào)法時(shí)差監(jiān)測方法與結(jié)果分析[J].測繪通報(bào)，2017(1)：26-29.

10.13474/j.cnki.11-2246.2017.0006.

2016-05-05；

2016-05-27

中科院百人計(jì)劃；國家863計(jì)劃(2013AA122402；2014AA123102)；國家自然科學(xué)基金(11273046)

章潔君(1991—)，女，碩士，研究方向?yàn)镚NSS精密數(shù)據(jù)處理。E-mail：jjzhang@shao.ac.cn