曹芬，楊旭海，李志剛，楊穎，孫保琪，馮初剛

?

TWSTFT站間鐘差殘差與測站環(huán)境溫度的相關(guān)性分析

曹芬1,2,3，楊旭海1,2，李志剛1,2，楊穎1,2,3，孫保琪1,2，馮初剛1,4

（1. 中國科學(xué)院國家授時中心，西安 710600；2.中國科學(xué)院精密導(dǎo)航定位與定時技術(shù)重點實驗室，西安 710600；3. 中國科學(xué)院大學(xué)，北京 100039；4. 中國科學(xué)院上海天文臺 上海 200030）

為了對衛(wèi)星雙向時間頻率傳遞（TWSTFT）站間鐘差殘差與測站環(huán)境溫度之間的關(guān)系進(jìn)行分析，利用2005年6月的C波段TWSTFT的觀測數(shù)據(jù)，對TWSTFT站間鐘差殘差與上海站、西安站的環(huán)境溫度的相關(guān)系數(shù)進(jìn)行了分析，并建立了聯(lián)結(jié)TWSTFT站間鐘差殘差與測站環(huán)境溫度的擬合函數(shù)。分析結(jié)果表明：測站環(huán)境溫度對TWSTFT站間鐘差殘差的影響是明顯的。對站間鐘差殘差進(jìn)行溫度擬合修正后，站間鐘差殘差RMS降低25%~50%。測站應(yīng)該保持觀測設(shè)備恒溫，這對于C波段TWSTFT的站間時間同步具有特殊意義。

環(huán)境溫度；衛(wèi)星雙向時間頻率傳遞（TWSTFT）；鐘差殘差；相關(guān)系數(shù)

0 引言

衛(wèi)星雙向時間頻率傳遞（TWSTFT）是現(xiàn)今世界上主要的遠(yuǎn)距離時間比對技術(shù)之一[1-2]。該方法利用傳播信號路徑的近似對稱性，原則上能抵消信號傳遞路徑上的大部分影響。

TWSTFT一般每周2~3次在固定時段內(nèi)觀測，每次觀測0.5 h，以5 min觀測為一組結(jié)果，每次觀測能給出2~3組結(jié)果，顯然這樣的觀測頻度無法發(fā)現(xiàn)觀測中的高頻部分誤差。隨著比對頻度的加大，在有些有條件的站，每天進(jìn)行24h觀測，并在TWSTFT比對結(jié)果中發(fā)現(xiàn)了嚴(yán)重影響比對精度的周日周期變化趨勢。因此，TWSTFT中的周日周期變化目前在國際上成為研究熱點[3]。

TWSTFT中的周日變化的成因是多因素的。單種誤差源不足以解釋TWSTFT中的周日變化。目前大多從實驗方法研究TWSTFT的誤差根源，以解釋部分周期性的變化，包括儀器系統(tǒng)誤差研究[4-6]和TWSTFT觀測原理誤差研究[7-9]等。對于TWSTFT，其主要的誤差來源包含設(shè)備發(fā)射與接收時延，衛(wèi)星轉(zhuǎn)發(fā)器時延，衛(wèi)星信號傳播路徑時延，時鐘時間信號的不穩(wěn)定，測時儀分辨率的影響等[9]。同時外界環(huán)境的變化也將影響時間比對精度，特別是溫度的影響不容忽視[3]。

有學(xué)者對測站環(huán)境溫度與TWSTFT精度的關(guān)系進(jìn)行了研究，并得出結(jié)論：在環(huán)境溫度超過12℃時，兩站的站間鐘差與環(huán)境溫度成正比關(guān)系；環(huán)境溫度低于12℃時，兩站的站間鐘差與環(huán)境溫度成反比關(guān)系[7]。美國國家標(biāo)準(zhǔn)技術(shù)研究所（NIST，National Institute of Standards and Technology）曾有實驗證明站間鐘差與環(huán)境溫度的擬合系數(shù)可達(dá)到（100±30）ps/k[4]。還有學(xué)者進(jìn)行研究而得出結(jié)論：環(huán)境溫度變化引起的站間鐘差的誤差最大[5]。

本文不討論測站環(huán)境溫度如何影響某種設(shè)備的具體機(jī)制，而利用2005年6月的TWSTFT觀測數(shù)據(jù)，分析站間鐘差殘差與測站環(huán)境溫度之間的相關(guān)系數(shù)，并建立聯(lián)結(jié)站間鐘差殘差與測站環(huán)境溫度的擬合函數(shù)，用以探求衛(wèi)星雙向時間頻率傳遞（TWSTFT）站間鐘差殘差與測站環(huán)境溫度之間的關(guān)系。

1 測站環(huán)境溫度與TWSTFT站間鐘差殘差的相關(guān)性分析

為了研究測站環(huán)境溫度與站間鐘差殘差的相關(guān)性，首先，根據(jù)TWSTFT原理得到站間鐘差，經(jīng)線性擬合后得到站間鐘差殘差?？紤]到站間鐘差殘差的周日周期變化應(yīng)當(dāng)與兩個站的環(huán)境溫度有關(guān)，我們對TWSTFT站間鐘差殘差與兩個測站的環(huán)境溫度進(jìn)行相關(guān)分析，并建立站間鐘差殘差與兩站環(huán)境溫度之間的關(guān)系，其流程如圖1所示。

圖1 站間鐘差殘差與測站環(huán)境溫度的相關(guān)性分析流程圖

1.1 站間鐘差的計算與擬合

1.1.1 TWSTFT比對原理

1.1.2 站間鐘差的線性擬合

不考慮頻率漂移，站間鐘差可以用線性函數(shù)擬合：

1.2 鐘差殘差與測站環(huán)境溫度的相關(guān)系數(shù)分析

1.3 建立聯(lián)結(jié)鐘差殘差與測站環(huán)境溫度的擬合函數(shù)

2 結(jié)果與分析

我們使用C波段TWSTFT的觀測數(shù)據(jù)分析測站環(huán)境溫度對TWSTFT站間鐘差殘差的影響。分析計算中選取2005年6月7日、6月9日、6月13日、6月15日和6月17日西安站和上海站的TWSTFT連續(xù)觀測資料和氣象數(shù)據(jù)。

2.1 TWSTFT鐘差計算

利用TWSTFT得到上海與西安的站間鐘差，并進(jìn)行線性擬合，用3倍均方差進(jìn)行粗差剔除，得到鐘差線性擬合后的殘差序列。表1給出了站間鐘差線性擬合系數(shù)及相應(yīng)的1值。

表1 上海站與西安站的站間鐘差線性擬合系數(shù)及RMS1值

注：，即式（2）中的兩系數(shù)。

從表1中可以看出，上海站與西安站的站間鐘差擬合的一次項系數(shù)較為平穩(wěn)，約為-0.000 560 ns/s。除6月7日的站間鐘差殘差均方根稍大外，其他天的鐘差殘差均方根均小于1 ns。

具有代表性的6月15日和17日的上海站與西安站的站間鐘差曲線及其線性擬合曲線，如圖2和3所示。

圖3 2005年6月17日上海與西安的站間鐘差曲線及其線性擬合曲線

2.2 站間鐘差殘差與測站環(huán)境溫度之間的相關(guān)系數(shù)

對TWSTFT站間鐘差殘差與上海站，西安站的環(huán)境溫度之間的相關(guān)系數(shù)的計算結(jié)果示于表2。

表2 TWSTFT站間鐘差殘差與測站環(huán)境溫度的相關(guān)系數(shù)

由表2可見，鐘差殘差與西安站環(huán)境溫度之間的相關(guān)系數(shù)約為0.607 5，與上海站環(huán)境溫度之間的相關(guān)系數(shù)約為0.104 4。因此TWSTFT站間鐘差殘差與TWSTFT比對鏈路的測站環(huán)境溫度存在一定的相關(guān)性。

2.3 建立聯(lián)結(jié)TWSTFT站間鐘差殘差與測站環(huán)境溫度的擬合函數(shù)及經(jīng)溫度擬合修正后的鐘差殘差

圖4和5給出了TWSTFT站間鐘差殘差及上海站，西安站的環(huán)境溫度變化趨勢。雙向時間比對的測站環(huán)境溫度在一天之內(nèi)變化大約10℃，并呈現(xiàn)周期性的變化。鐘差殘差也呈現(xiàn)一定的周期性變化趨勢，與測站環(huán)境溫度呈現(xiàn)相似的變化趨勢。

圖4 2005年6月15日上海與西安的站間鐘差殘差及上海站、西安站的環(huán)境溫度

由于TWSTFT站間鐘差殘差與測站環(huán)境溫度存在相關(guān)性，建立聯(lián)結(jié)站間鐘差殘差與測站環(huán)境溫度的擬合函數(shù)（如式（6）所示）。表3中給出了擬合函數(shù)（式（6））的各項系數(shù)及溫度擬合修正后的鐘差殘差2。

表3 站間鐘差殘差與測站環(huán)境溫度之間的擬合函數(shù)的系數(shù),2等

圖6和圖7給出了鐘差殘差，溫度擬合后的鐘差殘差及經(jīng)過溫度擬合修正后的鐘差殘差（）。

由圖6和7可以看出，利用測站環(huán)境溫度與鐘差殘差之間的擬合函數(shù)，我們能夠較好地模擬出與測站環(huán)境溫度相關(guān)的鐘差殘差曲線。由表3可見，經(jīng)過溫度擬合修正后，鐘差殘差的均方根降低了約25%~50%。這體現(xiàn)了鐘差殘差與測站環(huán)境溫度存在相關(guān)性。

3 結(jié)論

本文對C波段TWSTFT的觀測數(shù)據(jù)進(jìn)行處理，分析了站間鐘差殘差與測站環(huán)境溫度之間的相關(guān)系數(shù)，并建立了聯(lián)結(jié)鐘差殘差與測站環(huán)境溫度的擬合函數(shù)。分析結(jié)果表明：TWSTFT站間鐘差殘差與測站環(huán)境溫度存在相關(guān)性。經(jīng)過溫度擬合修正后，站間鐘差殘差RMS降低了約25%~50%。

測站環(huán)境溫度在一天中變化較大，尤其具有明顯大陸性氣候的西安，日溫差約為10℃，與TWSTFT站間鐘差殘差有較明顯的相關(guān)性。因此工作站應(yīng)該保持觀測設(shè)備恒溫，這對于C波段TWSTFT的站間時間同步具有特殊意義。

[1] Rec. ITU-R TF.1153-2. The operational use of two-way satellite time and frequency transfer employing PN codes[S]. 2003.

[2] GUREVICH E L. Two-way time transfer experiment via the Horizon satellite[J]. Mrtrolrigiu, 1995, 32(1): 51-53.

[3] PARKER T E, ZHANG V S, MCKINLEY A, et al. Investigation of instabilities in two-way time transfer[C]//34th annual precise time and time interval (PTTI) meeting. 2002: 381-388.

[4] ASCARRUNZ F G, JEFFERTS S R, PARKER T E, etc. Earth station errors in two-way time and frequency transfer [J]. IEEE transactions on instrumentation and measurement, 1997, 46(2): 205-208.

[5] 張虹, 李煥信, 李志剛. 雙向衛(wèi)星時間與頻率傳遞溫度效應(yīng)[J]. 宇航計測技術(shù), 2006, 26 (1): 24-28.

[6] ASCARRUNZ F G, JEFFERTS S R, PARKER T E. A delay calibration system for two-way satellite time and frequency transfer[C]//Proc. of IEEE International Frequency Control Symp. Pasadena: [s.n.], 1998: 250-253.

[7] 孫宏偉, 李志剛, 李煥信, 等. 衛(wèi)星雙向時間比對原理及比對誤差估算[J]. 宇航計測技術(shù), 2001, 21(2): 55-58.

[8] 武文俊, 李志剛, 楊旭海, 等. 衛(wèi)星雙向時間頻率傳遞中的Sagnac效應(yīng)[J]. 宇航學(xué)報, 2012, 33(7): 936-941.

[9] 葉玲玲, 石明華, 樓楊. 基于衛(wèi)星雙向法的時間比對系統(tǒng)構(gòu)建與分析[J]. 現(xiàn)代電子技術(shù), 2011, 34(24): 141-143.

[10] 李志剛, 程宗頤, 馮初剛, 等. 電離層預(yù)報模型研究[J]. 地球物理學(xué)報, 2007, 50(2): 327-337.

[11] 武文俊, 衛(wèi)星雙向時間頻率傳遞的誤差研究[D]. 西安: 中國科學(xué)院國家授時中心, 2012.

[12] 孔垚, 楊旭海, 孫保琪, 等. 基于“北斗”觀測數(shù)據(jù)的C波段雙向衛(wèi)星時間頻率傳遞中的電離層改正[J]. 時間頻率學(xué)報, 2012, 35(3): 148-155.

[13] 數(shù)學(xué)手冊編寫組. 數(shù)學(xué)手冊[M]. 北京: 人民教育出版社, 1979: 789-790.

Analysis of correlation between clock offset and environmental temperatures of stations for TWSTFT

CAO Fen1,2,3, YANG Xu-hai1,2, LI Zhi-gang1,2, YANG Ying1,2,3, SUN Bao-qi1,2, FENG Chu-gang1,4

(1. National Time Service Centre, Chinese Academy of Sciences, Xi′an 710600, China;2. Key Laboratory of Precision Navigation and Timing Technology, National Time Service Center,Chinese Academy of Sciences, Xi′an 710600, China;3. Graduate University of Chinese Academy of Sciences, Beijing 100039, China;4. Shanghai Astronomical Observatory of Chinese Academy of Sciences, Shanghai 200030, China)

For analyzing the correlation between the clock offset and the environmental temperatures of stations for TWSTFT(two-way satellite time frequency transfer), the correlation coefficients between the clock offset and the environmental temperatures of stations for TWSTFT are calculated and analyzed, and a fitted function that relates the residual of clock offset with the environmental temperatures of stations is constructed by using the observation data of TWSTFT (in C-band, for Shanghai-Xi′an link, during the June of 2005). The results show that the environmental temperatures obviously affect the residuals of clock offset. The residuals of clock offset are revised with the fitted function and the RMS(root-mean-square) values of the revised residuals of clock offset reduce by 25%~50%. The observation devices should be kept in an environment with constant temperature, which is of special significance to the time synchronization between the stations for TWSTFT in C-band.

environmental temperature; TWSTFT; residual of clock offset; correlation coefficient

TM935

A

1674-0637(2013)03-0148-08

2012-09-03

國家自然科學(xué)基金重點資助項目（11033004）；西部之光聯(lián)合學(xué)者資助項目（2007LH01）

曹芬，女，博士研究生，主要從事衛(wèi)星定軌與時間同步研究。