王維，張志斌，KALIUZHNYI MyKolay,于亮，唐正宏，王廣利

微型VLBI系統(tǒng)的溫度效應(yīng)改正與精度評估

王維1，張志斌1,2，KALIUZHNYI MyKolay3,于亮4，唐正宏1，王廣利1

（1. 中國科學(xué)院 上海天文臺，上海 200030；2. 中國科學(xué)院大學(xué)，北京 100049；3.尼古拉耶夫天文臺，尼古拉耶夫 54030； 4. 西安測繪總站，西安 710054）

微型甚長基線干涉測量（VLBI）系統(tǒng)（MVLBI）因其出色的性價比，在地球靜止衛(wèi)星監(jiān)測領(lǐng)域展示出其巨大的應(yīng)用前景。通過在上海建造包括兩臺天線與3個終端的微型VLBI并置站和開展試觀測，發(fā)現(xiàn)了微型VLBI時延序列隨溫度變化而出現(xiàn)異常波動。提出一種線性分段擬合算法分析了時延的該特性，并采用溫控措施克服了溫度變化所引起的異常效應(yīng)且獲得了穩(wěn)定的時延序列?；诖?，通過在并置站開展如換線法等測試，標(biāo)校了系統(tǒng)誤差并評估了單站精度。結(jié)果表明時延序列可能會在溫度驟變時產(chǎn)生異常波動，這是由GNSS授時模塊的溫敏效應(yīng)引起的。此外，如不考慮系統(tǒng)差，微型VLBI單站的測量精度可達2.6 ns。

甚長基線干涉測量；時延；溫度；精度評估

0 引言

隨著地球靜止衛(wèi)星（geostationary earth orbit satellite，GEOS）數(shù)目的快速增長，地球靜止軌道資源已變得相當(dāng)稀缺。當(dāng)前可取的一種方法是開展GEOS的高精度監(jiān)測，從而使得地球靜止軌道帶上可以容納更多目標(biāo)，這便對GEOS的高精度全天時全天候監(jiān)測提出了需求。作為一項幾何相關(guān)測時技術(shù)，甚長基線干涉測量（very long baseline interferometry，VLBI）廣泛應(yīng)用于天體物理、大地測量以及深空探測中[1]。近年來，VLBI發(fā)展趨勢已趨于明朗[2]，即專鏡專用，一方面，大口徑天線因其靈敏度較高常被分配以觀測射電源或深空探測器等暗弱目標(biāo)[3]，另一方面，13 m口徑的寬頻速轉(zhuǎn)小天線被用以開展站址坐標(biāo)以及地球定向參數(shù)的監(jiān)測[4]。然而，還有一種VLBI天線系統(tǒng)尚未被歸類，該類天線有著小于2m的 口徑，即微型VLBI系統(tǒng)，亦稱為被動相關(guān)測距系統(tǒng)（passive correlation ranging，PaCoRa），該系統(tǒng)最初由歐洲空間局于2013年提出并開展了仿真及原型機設(shè)計[5]。作為一種非侵入式的相關(guān)測距方法，可用于接收多類下行鏈路信號。由于其無源特性，故無需向衛(wèi)星輸入跟蹤信號，從而將對航天器和地面運行的潛在影響降至最低。此外，因地面站不需要大型天線，也極大地降低了系統(tǒng)成本；同期，烏克蘭尼古拉耶夫天文臺基于相同原理，獨立開發(fā)出一套針對電視通訊GEOS的微型VLBI接收系統(tǒng)，所測時延可用以對Eutelsat-13B GEOS開展公里級精度的定位[6-7]；中國科學(xué)院上海天文臺也于2019年建成了微型VLBI觀測網(wǎng)，可對Apstar6C GEOS開展百米級精度的連續(xù)軌道監(jiān)測[8]。

相比于主動式射電觀測[9-11]、光學(xué)觀測、傳統(tǒng)VLBI及連線干涉測量系統(tǒng)[12]等方法，微型VLBI系統(tǒng)具有低成本、全天時、全天候以及被動式觀測的特點，測站設(shè)計簡潔、布署便捷，且可實現(xiàn)長基線、長時段、無人值守的GEOS連續(xù)監(jiān)測。因其設(shè)計理念中遵循了最高性價比的建設(shè)原則，微型VLBI系統(tǒng)中省略了一些傳統(tǒng)VLBI所必須的一些部件或觀測方法，如相位校正或射電源校正等，這便對該系統(tǒng)的誤差校正方法提出了新需求，因此，在遠(yuǎn)端站布設(shè)前對微型VLBI并置站開展誤差標(biāo)校方法研究顯得尤為必要。為進一步精化微型VLBI觀測系統(tǒng)，需對其隨機性和系統(tǒng)性誤差開展深入研究。本文對微型VLBI原理、本地測試過程以及系統(tǒng)差評定方法作闡述，并提出一種用以分析溫度序列及相應(yīng)的時延序列的擬合方法，來解釋微型VLBI時延變化原因，最后給出系統(tǒng)精度評估方法及結(jié)論。

1 觀測原理及實驗

1.1 VLBI觀測GEO衛(wèi)星原理

如圖1所示，微型VLBI觀測GEOS原理為：GEOS轉(zhuǎn)發(fā)信號的波前依次到達測站1和測站2，通過對兩個測站記錄信號開展相關(guān)處理，可得到兩站間的觀測時延，進而通過構(gòu)建星站之間的測量模型，開展誤差方程求解，即可實現(xiàn)利用站間時延對GEOS開展定位和定軌。

圖1 微型VLBI觀測GEOS示意圖

1.2 微型VLBI單站信號采集及相關(guān)處理

2018年5月，項目組在上海佘山搭建了兩臺并置站天線，并研發(fā)了對應(yīng)的A和B兩套接收終端組成并置微型VLBI站。期間微型VLBI相關(guān)時延表現(xiàn)出異常變化，這種變化常出現(xiàn)在本地時間的中午時段。為進一步識別和鎖定問題，引入測溫儀來監(jiān)測時延序列和A、B終端不同模塊處所對應(yīng)溫度。2018年11月，項目組引入了第3套終端C，至此得到兩天線對應(yīng)3套接收終端的微型VLBI并置站，基于此開展了短基線與零基線的相關(guān)測試。

圖2 微型VLBI單站組成圖

觀測目標(biāo)選取GEOS Apstar6C（134°E；垂直極化；碼率：27 500 kbit/s；頻率：12.395 GHz），該衛(wèi)星為電視信號轉(zhuǎn)發(fā)衛(wèi)星。

2 時延與溫度的變化關(guān)系

2.1 時延序列與溫度序列間的相關(guān)性

圖3所示為時延序列與其對應(yīng)時刻的環(huán)境溫度。與轉(zhuǎn)發(fā)式測距設(shè)備中時延與溫度的基本呈線性變化[13]不同，由于并置站的微型VLBI時延觀測噪聲為幾十納秒，時延與溫度間的線性關(guān)系并不顯著，其表現(xiàn)為，當(dāng)溫度在短期內(nèi)劇烈上升時，時延序列出現(xiàn)波動。然而，這種相關(guān)并非完全相關(guān)。時延序列與溫度序列中變化的峰值位置是一一對應(yīng)的，但其他部分的變化卻是不一致的。通過計算皮爾遜積差相關(guān)系數(shù)可知兩者間相關(guān)度為0.33，這也表明兩種序列是中度相關(guān)。但從圖3來看，較大溫差不會引起更高的時延變化。

注：零刻度位于約化儒略日58 382日00:00:00

2.2 時延序列的擬合方法

為分析時延隨溫度的變化特點，提出一種非連續(xù)線性分段擬合方法（discontinuous piecewise linear fitting，DPWLF）來處理該1維時間序列。該算法中的關(guān)鍵是在序列中選取合適的分段位置?？紤]到計算效率，提出一種兩步法，如下所示：

第二步的目的是在第一步所找出分段點的基礎(chǔ)上找出各段間最佳的分段位置。由第一步分段方法可知，每段固定步長的分段序列之間必有一點可使得該點兩側(cè)序列的擬后殘差值為最小。對所有固定步長序列中的點，都將計算其兩側(cè)序列的線性擬后殘差，并求出最小殘差對應(yīng)的點位，最終，這些點位即是DPLWF中的近似分段位置。

不同于連續(xù)線性分段擬合法，DPWLF可有效減小時間序列的擬后殘差，特別是對于某些分散的和快速變化的時間序列。圖4中給出一例，圖4（a）中是某次調(diào)試所得時延序列，圖4（b）和（c）分別給出了兩步法中的擬合線段。其中SF設(shè)為2.0，搜索步長為500點。通過DPWLF的兩步擬合，時延序列（單位為m）的擬后殘差由4.7 m降到3.6 m。兩步法中第二步的必要性可由圖4（d）說明，其最小殘差點位于第3段分段序列中379點處，即該點應(yīng)被選為序列2和3之間的分段點。

3 時延序列變化特性及解釋

3.1 時延序列的溫度特性

這里用小樣本來統(tǒng)計分析微型VLBI時延的溫度變化特性。結(jié)合DPWLF方法，對時延相關(guān)所得的微型VLBI時延序列做了擬合，其中，將穩(wěn)定時延序列定義為其時延率變化不超過特定值，如（7～8）m/min（時延單位已換算為m），考慮到基線與GEOS目標(biāo)方向近乎垂直的關(guān)系，將該值設(shè)為30 m/d。

如圖5（a）所示，利用DPWLF（SF = 1.5）方法，這里擬合出3段相對穩(wěn)定的時延序列。圖5（b）、（c）和（d）分別給出了這3段時延對應(yīng)的溫度變化。這些溫度序列均表現(xiàn)出典型的“U”形變化。當(dāng)溫度在較短時段內(nèi)（約2 h）且在超過5℃的范圍內(nèi)下降和上升后，此后的時延序列將出現(xiàn)隆起。因此可以推斷時延序列的不穩(wěn)定變化是由于短期內(nèi)溫度緩慢下降和急劇升高而引起的。此外，值得注意的是，時延序列的這種變化是不能用同一的溫度模型來改正的，因為不同的接收終端有著不同的溫敏特性。如圖6所示，與終端C有關(guān)的基線在第95～110 h出現(xiàn)了跳動，A、B基線卻保持平穩(wěn)。

3.2 時延序列變動原因

GEOS通常一周開展一次機動，一般不會連續(xù)3日在每日同一時段內(nèi)機動，因此首先排除GEOS目標(biāo)機動的原因，如僅從微型VLBI接收系統(tǒng)組件角度來看，其中3個組件可能會受溫度變化影響，即差分（或相對）電纜時延變化、GNSS授時精度隨溫度變化或者是示波器的模數(shù)轉(zhuǎn)換（ADC）模塊。

同軸電纜環(huán)境溫度變化會引起同軸電纜內(nèi)金屬綁線的伸縮，從而一定程度上引起同軸電纜電容的改變，進而改變電纜時延。文獻[14-15]分析了同類型號的電視電纜時延相對于溫度的變化，結(jié)果表明，當(dāng)溫差達到30℃時，幾米長的電纜所產(chǎn)生附加時延的量級在幾個皮秒。本試驗中，兩臺微型VLBI并置站間電纜長度差為14.0 m。我們引入了測溫設(shè)備，實時記錄了兩條電纜溫度，結(jié)果表明，針對兩條電纜所記錄的日最大溫差約15℃，由溫度變化所引起的相對電纜時延將不超過幾十皮秒（對應(yīng)厘米級），這顯然不是時延序列變化的主要因素。

為了進一步從GNSS授時模塊和計算機ADC中鎖定引起時延序列波動的部件，這里對所有天寶GNSS授時模塊做了隔熱處理，如圖7所示，并由此得到了平穩(wěn)的時延序列。據(jù)此確定了GNSS授時模塊的熱敏特性是引起時延序列變化的原因。

圖7 GNSS授時模塊隔熱示意圖

4 系統(tǒng)誤差校正與系統(tǒng)精度評估

在對各接收終端的GNSS授時模塊采取溫控措施后，可利用平穩(wěn)時延來校正設(shè)備間的系統(tǒng)誤差。不同于傳統(tǒng)VLBI方法利用射電源校正或相位校正方法來觀測GEOS目標(biāo)，單一功能的微型VLBI站因暫不能改變天線指向，只能觀測一定“視場”內(nèi)的GEOS目標(biāo)，這樣，在將微型VLBI測站部署到遠(yuǎn)端站點之前，必須嚴(yán)格標(biāo)定測站間的相對系統(tǒng)差，這包括部件間相對誤差標(biāo)定以及系統(tǒng)開關(guān)機測試。這里通過換線法來標(biāo)記測站間的系統(tǒng)差，具體如下一節(jié)所述。

4.1 設(shè)備時延校正

式（5）中：中包括兩個星站距離和間的差，即幾何時延（或）；兩條電視電纜和間的線纜相對延遲（或）；表示線纜與終端相對終端間的時延；表示線纜與終端相對終端間的時延；和統(tǒng)稱為設(shè)備延遲。線纜和設(shè)備名稱具體如圖9所示。

表1 利用圖8（f）（g）（h）所測時延解算的并置站間幾何時延、線纜相對時延和設(shè)備間相對時延

4.2 系統(tǒng)精度評估

表2 3條基線時延與其閉合時延平均偏差與RMS統(tǒng)計

5 結(jié)語

系統(tǒng)差的識別與改正是微型VLBI系統(tǒng)在部署到遠(yuǎn)端站點前最為關(guān)鍵的步驟。試驗結(jié)果表明，GNSS授時模塊的溫敏特性會引起時延序列的異常波動。通過提出連續(xù)分段線性擬合法，了解到系統(tǒng)環(huán)境溫度的緩慢下降與迅速上升會引起該異常波動，且這種異常隨不同設(shè)備而異，因此無法構(gòu)建統(tǒng)一模型來改正該異常，只能通過溫控的方法來獲得平穩(wěn)時延?；诖嗽O(shè)計了換線法和閉合時延分別評估了系統(tǒng)精度和單站接收精度，并驗證了評估結(jié)果，這些結(jié)果均與估計值有著較好符合。本研究為后續(xù)微型VLBI長基線布站奠定了基礎(chǔ)。

致謝：感謝中國科學(xué)院上海天文臺李巖、毛銀盾等老師對系統(tǒng)組裝與觀測試驗所提供的大力支持。

[1] Schuh H, Behrend D. VLBI: a fascinating technique for geodesy and astrometry[J]. Journal of Geodynamics, 2012(61)： 68-80.

[2] 方成. 關(guān)于我國未來天文大設(shè)備的一點戰(zhàn)略思考[J]. 中國科學(xué)(物理學(xué) 力學(xué) 天文學(xué)), 2012, 42(12): 1265-1273.

[3] 張波.FAST的VLBI觀測與科學(xué)[J]. 中國科學(xué)(物理學(xué) 力學(xué) 天文學(xué)), 2017, 476(6): 7-16.

[4] NIELL A, WHITNEY A, PETRACHENKO B, et al. VLBI2010: current and future requirements for geodetic VLBI systems[J]. International VLBI Service for Geodesy and Astrometry, 2005: 13-40.

[5] Rodriguez L, Krier G, Thill M, et al. Passive ranging for geostationary satellites: on a novel system and operational benefits[C]//13th International Conference on Space Operations, AIAA: Pasadena, 2014: 1857-1863.

[6] Bushuev F, Kalyuzhny M A, Slivinsky A P, et al. Ranging geostationary telecommunication satellites with satellite TV signals[J]. Telecommunications & Radio Engineering, 2016, 53(5): 5-16.

[7] Bushuev F, Kaliuzhnyi M, Sybiryakova Y, et al. Results of the ongoing monitoring of the position of a geostationary telecommunication satellite by the method of spatially separated basis receiving of digital satellite television signals[J]. Latvian Journal of Physics & Technical Sciences, 2016, 53(5): 5-16.

[8] 張志斌, 王維, 楊鵬, 等. 一種用以同步微型監(jiān)測的微型VLBI網(wǎng)[EB/OL]. (2020-04-23)[2020-08-25]. http://gb.oversea.cnki.net/KCMS/detail/11.1859.V.20200422.1839.008.html.

[9] Peltzer K E. Apollo unified S-band system[EB/OL]. (1996-07-29)[2020-08-05]. https://ntrs.nasa.gov/citations/19660018739.

[10] Harles G, Siebert P, Eberlein E, et al. A novel ranging method using DVB-S transport stream packets[C] // Proceeding of the Seventh International Workshop on Digital Signal Processing Techniques for Space Communications, Sesimbra: [s.n.], 2001.

[11] 雷輝, 李志剛, 楊旭海, 等. 基于轉(zhuǎn)發(fā)式的北斗衛(wèi)星導(dǎo)航系統(tǒng)地球靜止軌道衛(wèi)星精密定軌試驗[J]. 測繪學(xué)報, 2011, 40(S1): 31-33.

[12] 杜蘭, 李曉杰, 王若璞. 基于同波束CEI的GEO共位衛(wèi)星相對軌道監(jiān)視[J]. 大地測量與地球動力學(xué), 2012, 32(3): 54-58.

[13] 武文俊, 李志剛, 李孝輝, 等. 轉(zhuǎn)發(fā)式測距設(shè)備時延中的溫度效應(yīng)[J]. 時間頻率學(xué)報, 2013, 36(2): 113-119.

[14] Czuba K, Sikora D. Temperature stability of coaxial cables[J]. Acta Physica Polonica A, 2011, 119(4): 553-557.

[15] Parashare C R, Richard F B. 75?Transmission System[EB/OL]. (2011-09-26)[2020-08-05]. http://eor.berkeley.edu/wpcontent/uploads/2011/09/p011.cparashare.pdf.

Thermal effect correction and accuracy evaluation on MVLBI

WANG Wei1, ZHANG Zhi-bin1,2, KALIUZHNYI MyKolay3, YU Liang4,TANG Zheng-hong1, WANG Guang-li1

(1. Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China;2. University of Chinese Academy of Sciences, Beijing 100049, China;3. Mykolayiv Astronomical Observatory, Research Institute, Mykolayiv 54030, Ukraine;4. Xi’an Division of Surveying and Mapping, Xi’an 710054, China)

With the advantages of being cost effective, the micro VLBI system (MVLBI) shows a great application prospect in geostationary satellite monitoring. We built a co-location station with two antennas and three terminals in Shanghai, and found that the time delay series fluctuates abnormally with temperature changes. A linear piecewise fitting algorithm is developed to analyse this characteristic. Besides, temperature controls are taken to overcome the abnormal effect caused by temperature changing, then the stable delay series is obtained. Thus, by performing tests such as the cable changing method at the co-location station, systematic errors are calibrated and the accuracy of single station is evaluated. It shows the delay series may vary abnormally with the sharp temperature change, due to the thermosensitive in GNSS timing module. In addition, the tested accuracy of the single MVLBI system can reach up to 2.6 ns regardless of systematic errors.

VLBI; time delay; temperature; accuracy evaluation

10.13875/j.issn.1674-0637.2021-02-0142-11

王維, 張志斌,KALIUZHNYI M, 等. 微型VLBI系統(tǒng)的溫度效應(yīng)改正與精度評估[J]. 時間頻率學(xué)報, 2021, 44(2): 142-152.

2020-11-16；

2020-12-05

國家自然科學(xué)基金資助項目（1703067；11873077）；探月工程三期測控系統(tǒng)關(guān)鍵技術(shù)研究資助項目