周東旭 袁運(yùn)斌 李子申 劉西鳳

(1)中國(guó)科學(xué)院測(cè)量與地球物理研究所動(dòng)力大地測(cè)量學(xué)重點(diǎn)實(shí)驗(yàn)室，武漢 430077 2)中國(guó)科學(xué)院研究生院，北京100049)

GPS接收機(jī)儀器偏差的長(zhǎng)期變化特性分析*

周東旭1，2)袁運(yùn)斌1)李子申1，2)劉西鳳1)

(1)中國(guó)科學(xué)院測(cè)量與地球物理研究所動(dòng)力大地測(cè)量學(xué)重點(diǎn)實(shí)驗(yàn)室，武漢 430077 2)中國(guó)科學(xué)院研究生院，北京100049)

基于歐洲定軌中心2001—2009年發(fā)布的部分IGS測(cè)站接收機(jī)儀器偏差估值，分析了不同類型接收機(jī)DCB的長(zhǎng)期變化特性。結(jié)果表明：各類接收機(jī)DCB的長(zhǎng)期時(shí)變特征均具備一定的周期性，其中以月和年周期最為明顯;接收機(jī)DCB的長(zhǎng)期變化中還存在一定的波動(dòng)，部分接收機(jī)DCB的穩(wěn)定性隨時(shí)間積累逐漸變差，甚至存在偏移。

接收機(jī)儀器偏差;長(zhǎng)期變化;全球?qū)Ш叫l(wèi)星系統(tǒng);歐洲定軌中心;周期性

1 引言

GPS觀測(cè)量中的儀器偏差(DCB)，是GPS衛(wèi)星和用戶接收機(jī)的硬件對(duì)信號(hào)產(chǎn)生延遲影響的統(tǒng)稱［1］。它是利用GPS觀測(cè)數(shù)據(jù)解算電離層總電子含量的重要誤差源，同時(shí)反映GPS接收機(jī)的硬件性能。分析接收機(jī)DCB隨時(shí)間的變化特性，有助于進(jìn)一步掌握接收機(jī)的工作狀態(tài)，以及TEC的精確確定及預(yù)報(bào)。國(guó)內(nèi)外學(xué)者基于GPS實(shí)測(cè)數(shù)據(jù)，采用不同方法對(duì)DCB進(jìn)行了探討［2-13］，其中，歐洲定軌中心CODE(Center for Orbit Determination in Europe)基于P1、P2碼觀測(cè)量對(duì)全球IGS基準(zhǔn)站的接收機(jī)DCB進(jìn)行了多年的研究和計(jì)算，其發(fā)布的接收機(jī)DCB估值具有較高的精度和可靠性，得到廣泛的認(rèn)可和應(yīng)用。但目前的研究大多集中于站星組合或衛(wèi)星DCB估計(jì)算法，鮮有專門針對(duì)接收機(jī)DCB變化特性的深入分析。

本文基于全球IGS基準(zhǔn)站2001—2009年的CODE接收機(jī)DCB估值，探討了接收機(jī)DCB的長(zhǎng)期變化特性，同時(shí)結(jié)合譜分析和濾波方法，研究了接收機(jī)DCB長(zhǎng)期變化的周期性特征。

2 GPS接收機(jī)DCB長(zhǎng)期變化特性分析

首先從CODE在2001—2009年期間發(fā)布的各IGS測(cè)站接收機(jī)DCB估值中，選取了不同區(qū)域IGS站的接收機(jī)DCB估值，以及同一區(qū)域內(nèi)部分IGS站的接收機(jī)DCB估值，探討接收機(jī)DCB連續(xù)9年的長(zhǎng)期變化特性;然后采用譜分析和濾波方法，具體分析和掌握接收機(jī)DCB的周期變化特征。

2.1 接收機(jī)DCB長(zhǎng)期變化分析

圖1(圖中圓圈表示該時(shí)段有接收機(jī)設(shè)備的更換，虛線方框標(biāo)注部分為DCB存在偏移的時(shí)段)顯示了全球不同區(qū)域內(nèi)8個(gè)IGS站9年的接收機(jī)DCB估值時(shí)間序列(IGS站具體信息見(jiàn)表1，其中時(shí)間序列空白部分為數(shù)據(jù)缺失)，其他IGS站的接收機(jī)DCB估值時(shí)間序列與之類似。為便于分析，X軸以年序號(hào)標(biāo)注。從圖1和表1可以看到，如不考慮接收機(jī)硬件設(shè)備更換引起的數(shù)據(jù)跳變，各站的接收機(jī)DCB估值時(shí)間序列除2001—2003年的波動(dòng)幅度相對(duì)較大外，其他時(shí)段整體變化比較穩(wěn)定，平均波動(dòng)幅度最大為1.88 ns。不同類型接收機(jī)的DCB估值和波動(dòng)幅度存在差異，平均波動(dòng)幅度的差異可達(dá)1 ns，具有同一類型接收機(jī)的KUNM、GUAM兩站，其DCB均值與標(biāo)準(zhǔn)差也存在差異，其中平均波動(dòng)幅度差異為0.3 ns;除NLIB站由于接收機(jī)硬件設(shè)備的更換造成時(shí)間序列中斷，難以辨析周期變化外，其他站都呈現(xiàn)出近似年變化周期，但各站的周期變化趨勢(shì)和變化幅度存在差異。

KERG站在2004—2009年的平均波動(dòng)幅度為1.77 ns，在沒(méi)有更換接收機(jī)硬件設(shè)備情況下(最近一次更換接收機(jī)硬件設(shè)備的時(shí)間為2004年1月13日)，該站的接收機(jī)DCB估值時(shí)間序列在2007—2009年出現(xiàn)了較大的波動(dòng)和明顯的偏移現(xiàn)象，偏移量達(dá)到2.79 ns。NLIB站在接收機(jī)硬件設(shè)備更換前的2004年也出現(xiàn)明顯偏移現(xiàn)象，偏移量為2.13 ns。可見(jiàn)，隨著工作時(shí)間的積累，GPS接收機(jī)硬件設(shè)備的老化，一些GPS接收機(jī)的DCB穩(wěn)定性會(huì)變差，甚至出現(xiàn)偏移現(xiàn)象。此外，MADR站由于接收機(jī)內(nèi)部硬件的更換，其時(shí)間序列在2003年也發(fā)生跳變。

圖1 不同地區(qū)GPS接收機(jī)DCB估值時(shí)間序列(2001—2009年)Fig.1 Time series of estimated value of receivers’DCB in different region(2001-2009)

表1 IGS站信息表Tab.1 Information of IGS stations

圖2為中國(guó)區(qū)域內(nèi)KUNM、BJFS、WUHN、TWTF 4站2001—2009年的接收機(jī)DCB估值時(shí)間序列(測(cè)站信息見(jiàn)表2)。圖2顯示，各站的接收機(jī)DCB估值除2001—2003年的波動(dòng)偏大外，其他時(shí)段整體變化平穩(wěn)，不同站的DCB估值和波動(dòng)幅度存在差異，其中平均波動(dòng)幅度差異最大為 1.04 ns;除TWTF站外，其余3站都有近似的年周期變化特征，并且KUNM、BJFS兩站時(shí)間序列的周期變化趨勢(shì)相反。

由圖1～2和表1～2的分析結(jié)果可知，接收機(jī)DCB在長(zhǎng)期變化過(guò)程中整體相對(duì)穩(wěn)定，但具有一定的波動(dòng)和周期性，不同IGS站之間，GPS接收機(jī)DCB的波動(dòng)和變化周期存在差異，并且在部分時(shí)段，個(gè)別接收機(jī)的DCB穩(wěn)定性隨時(shí)間積累變差，甚至出現(xiàn)偏移現(xiàn)象。

表2 KUNM、BJFS、WUHN、TWTF站信息Tab.2 Information of IGS stations(KUNM，BJFS，WUHN，TWTF)in the same region(China)

圖2 不同類型接收機(jī)DCB估值時(shí)間序列(2001—2009年)Fig.2 Time series of receivers’DCB of different type(2001-2009)

2.2 接收機(jī)DCB的周期特征分析

圖3顯示了KUNM站GPS接收機(jī)DCB估值的譜分析結(jié)果。其中，接收機(jī)DCB具有明顯的年變化周期，同時(shí)具有一定的月、季節(jié)、半年等變化周期。

圖3 KUNM站接收機(jī)儀器偏差估值譜分析結(jié)果(2001—2008年)Fig.3 Results of receiver’s DCB with spectrum analysis (2001-2008)

圖4為KUNM站接收機(jī)DCB估值的各主要變化周期的濾波結(jié)果，其他IGS基準(zhǔn)站的結(jié)果與之類似，在此不一一列舉。由圖4可見(jiàn)，月周期與年周期整體變化顯著，振幅較大，季節(jié)與半年周期項(xiàng)振幅相對(duì)較小;其中月周期和季節(jié)周期項(xiàng)振幅隨時(shí)間變化較劇烈，半年與年周期項(xiàng)振幅變化相對(duì)平緩。

為進(jìn)一步分析接收機(jī)DCB估值的月、季節(jié)、半年、年等變化周期的大小與變化情況，選取KUNM等IGS站的接收機(jī)DCB各變化周期9年的最大振幅和平均振幅作為統(tǒng)計(jì)指標(biāo)(表3)進(jìn)行分析。由表3可見(jiàn)，月周期變化中，各IGS站的最大振幅均大于1.1 ns，而平均振幅最大為0.69 ns，二者的平均差異達(dá)到1.66 ns，表明月周期變化隨時(shí)間變化比較劇烈;季節(jié)周期和半年周期變化的最大振幅、平均振幅以及二者的平均差異都相對(duì)較小，即季節(jié)周期和半年周期變化在GPS接收機(jī)DCB長(zhǎng)期變化中不顯著;年周期變化中，各站的最大振幅和平均振幅相對(duì)其他周期變化整體偏大，平均振幅的均值為0.86 ns，最大振幅和平均振幅的平均差異為0.45 ns，表明年周期變化在GPS接收機(jī)DCB周期變化中最為顯著，并且隨時(shí)間變化相對(duì)平緩。

3 結(jié)論

1)GPS接收機(jī)DCB的變化具有一定的周期性，主要包括月、季節(jié)、半年、年等變化周期，其中月周期變化隨時(shí)間變化劇烈，年周期變化在整體變化過(guò)程中最為明顯，平均振幅可達(dá)ns級(jí)。并且不同類型GPS接收機(jī)DCB的周期變化特性存在差異。

2)GPS接收機(jī)DCB變化整體呈比較穩(wěn)定的特性，但也存在一定的波動(dòng)，并且一些接收機(jī)DCB的穩(wěn)定性隨著工作時(shí)間的積累會(huì)逐漸變差，甚至產(chǎn)生偏移現(xiàn)象。

圖4 KUNM站接收機(jī)DCB估值濾波結(jié)果(2001—2008年)Fig.4 Filtered results of receiver’s DCB(2001-2008)

表3 周期變化統(tǒng)計(jì)表(單位：ns)Tab.3 Statistics of periodic change terms(unit：ns)

3)目前，常借助電離層模型進(jìn)行接收機(jī)DCB的精確確定，但這是一種后處理方法，并且當(dāng)觀測(cè)數(shù)據(jù)量不足時(shí)，接收機(jī)DCB的解算精度會(huì)嚴(yán)重降低，此時(shí)，利用GPS接收機(jī)DCB長(zhǎng)期變化特征監(jiān)測(cè)得到的預(yù)報(bào)值進(jìn)行實(shí)時(shí)導(dǎo)航、定位、TEC的提取具有重要意義，并且有助于掌握接收機(jī)的工作狀態(tài)。

由于CODE接收機(jī)DCB估值是與電離層TEC值同時(shí)確定的，本文所分析的周期變化特征是否受電離層變化周期的影響有待進(jìn)一步分析。

致謝 感謝張寶成、蘇曉慶博士研究生的幫助，感謝CODE提供儀器偏差值。

1 袁運(yùn)斌.基于GPS的電離層監(jiān)測(cè)及延遲改正理論與方法的研究［D］.中國(guó)科學(xué)院測(cè)量與地球物理研究所，2002.(Yuan Yunbin.Study on theories and methods of correcting ionospheric delay and monitoring ionosphere based on GPS［D］.Institute of Geodesy and Geophysics of Chinese Academy of Sciences，2002)

2 Brunini C and Azpilicueta F.Accuracy assessment of the GPS-based slant total electron content［J］.Journal of Geodesy，2009，83：773-785.

3 Brunini C and Azpilicueta F.GPS slant total electron content accuracy using the single layer model under different geomagnetic regions and ionospheric conditions［J］.Journal of Geodesy，2010，84：293-304.

4 袁運(yùn)斌，歐吉坤.GPS觀測(cè)數(shù)據(jù)中的儀器偏差對(duì)確定電離層延遲的影響及處理方法［J］.測(cè)繪學(xué)報(bào)，1999，28(2)：110-114.(Yuan Yunbin and Ou Jikun.The effects of instrumental bias in GPS observations on determining ionospheric delays and the methods of its calibration［J］.Acta Geodaetica et Cartographica Sinica，1999，28(2)：110-114)

5 Sardon E and Zarraoa N.Estimation of the transmitter and receiver differential biases and the ionospheric total electron content from global positioning system observations［J］.Radio Science，1994，19(3)：577-586.

6 Ma G and Maruyama T.Derivation of TEC and estimation of instrumental biases from GEONET in Japan［J］.Annales Geophysicae，2003，21：2 083-2 093.

7 Gao Y，et al.Modeling and estimation of C1-P1 bias in GPS receivers［J］.Journal of Geodesy，2001，74：621-626.

8 Ciraolo L，et al.Calibration errors on experimental slant total electron content(TEC)determined with GPS［J］.Journal of Geodesy，2007，81：111-120.

9 蔡昌盛，李征航，張小紅.GPS系統(tǒng)硬件延遲修正方法的探討［J］.測(cè)繪通報(bào)，2002，4：15-16.(Cai Changsheng，Li Zhenghang and Zhang Xiaohong.A study of calibration method of GPS satellite and receiver instrumental biases［J］.Bulletin of Surveying and Mapping，2002，4：15-16)

10 常青，等.GPS系統(tǒng)硬件延遲修正方法［J］.科學(xué)通報(bào)，2000，45(15)：1 676-1 680.(Chang Qing，et al.Correct way of GPS hardware delay［J］.Chinese Science Bulletin，2000，45(15)：1 676-1 680)

11 章紅平，施闖，唐衛(wèi)明.地基GPS區(qū)域電離層多項(xiàng)式模型與硬件延遲統(tǒng)一解算分析［J］.武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)，2008，33(8)：805-809.(Zhang Hongping，Shi Chuang and Tang Weiming.United solution to polynomial VTEC modeling and DCB analysis using ground based GPS observations［J］.Geomatics and Information Science of Wuhan University，2008，33(8)：805-809)

12 宋小勇，等.GPS接收機(jī)碼間偏差(DCB)的確定［J］.大地測(cè)量與地球動(dòng)力學(xué)，2009，(1)：127-131.(Song Xiaoyong，et al.Determination of GPS receiver’s DCB［J］.Journal of Geodesy and Geodynamics，2009，(1)：127-131)

13 安家春，等.基于單站的硬件延遲求解方法［J］.大地測(cè)量與地球動(dòng)力學(xué)，2009，(2)：86-90.(An Jiachun，et al.Methods for estimation of differential code biases by using single station GPS measurements［J］.Journal of Geodesy and Geodynamics，2009，(2)：86-90)

ANALYSIS OF LONG-TERM VARIATIONS OF GPS RECEIVERS’DIFFERENTIAL CODE BIAS

Zhou Dongxu1，2)，Yuan Yunbin1)，Li Zishen1，2)and Liu Xifeng1)

(1)Institute of Geodesy and Geophysics，Chinese Academy of Sciences，Key Laboratory of Dynamic Geodesy，Wuhan 430077 2)Graduate University of Chinese Academy of Sciences，Beijing100049)

The long-term variations of a variety of GPS receivers’DCB were analyzed based on the DCB published by the Center for Orbit Determination in Europe during 2001—2009.The results show that there are periodic behaviors in the time-series of the change receivers’DCB，among then the monthly and yearly periods are the most predominant.In addition，certain fluctuations underlying the long-term behaviors in GPS receivers’DCB can also be found，and some of the GPS receivers may suffer the deterioration of DCB stability，or even somewhat deviation.

DCB(Differential Code Bias);long-term variations;GPS;CODE(Center for Orbit Determination in Europe);periodicity

周東旭，男，1982年生，碩士研究生，主要研究方向?yàn)榛贕NSS的電離層研究.E-mail：zhoudongxu12@163.com

1671-5942(2011)05-0114-05

2011-03-10

國(guó)家自然科學(xué)基金(40890160，41021003，40625013);中國(guó)科學(xué)院研究生科技創(chuàng)新實(shí)踐項(xiàng)目(Y005381038);導(dǎo)航年會(huì)青優(yōu)資助課題(CSNC2010-QY-013);國(guó)家863計(jì)劃項(xiàng)目(2008AA12Z303)

P227;P207

A