姚當(dāng),弓劍軍,3,馬浪明,劉佳,吳文雅,徐磊,呂林蔚,吳元偉,張秀忠,王廣利,舒逢春,楊旭海,3

(1.中國(guó)科學(xué)院 國(guó)家授時(shí)中心,西安 710600；2.中國(guó)科學(xué)院大學(xué),北京 100049；3.中國(guó)科學(xué)院 精密導(dǎo)航定位與定時(shí)技術(shù)重點(diǎn)實(shí)驗(yàn)室,西安 710600；4.中國(guó)科學(xué)院 上海天文臺(tái),上海 200030)

0 引言

甚長(zhǎng)基線干涉測(cè)量(very long baseline interferometry,VLBI)是實(shí)現(xiàn)天球參考架和地球參考架的關(guān)鍵技術(shù),主要由于它觀測(cè)的目標(biāo)為遙遠(yuǎn)的河外射電源,在空間上其位置是固定的。VLBI通過(guò)確定特定射電源的位置來(lái)實(shí)現(xiàn)國(guó)際天球參考架(ICRF),并通過(guò)標(biāo)定30個(gè)分布全球的射電望遠(yuǎn)鏡的位置和速度來(lái)實(shí)現(xiàn)國(guó)際地球參考架(ITRF)。VLBI通過(guò)測(cè)量確定地球定向參數(shù)(EOP),包括只能由其測(cè)定的UT1和章動(dòng),來(lái)連接ICRF和ITRF。

國(guó)際測(cè)地測(cè)天VLBI服務(wù)(IVS)組織協(xié)調(diào)全球天線進(jìn)行觀測(cè)研究,主要包括8～10個(gè)站的24 h全天觀測(cè),其世界時(shí)UT1測(cè)量精度在6～7 μs(相對(duì)于IERS C04 05系列),及2～4個(gè)站的每天1 h UT1加強(qiáng)觀測(cè),精度為15 μs[1]。為提高測(cè)量精度,IVS提出新一代VGOS(VLBI global observation system)系統(tǒng)：1 mm位置精度和0.1 mm/年的速度精度；連續(xù)測(cè)量站位置和地球定向參數(shù)及24 h內(nèi)得到初始測(cè)地結(jié)果[2]。通過(guò)仿真研究,給出了使用12 m口徑、轉(zhuǎn)速快的小天線,寬帶寬,靈敏度高,數(shù)據(jù)傳輸快等特點(diǎn)的天線系統(tǒng)來(lái)達(dá)到上述要求,同S/X雙頻觀測(cè)系統(tǒng)相比較,天線口徑明顯減小,S/X系統(tǒng)的絕大多數(shù)天線口徑都大于20 m；帶寬增大,數(shù)據(jù)率明顯提高,兩者具體比較可參見(jiàn)表1,其中Jy=10-26W·m-2·Hz-1。

表1 S/X系統(tǒng)與VGOS系統(tǒng)對(duì)比

如今在天線制造,數(shù)字電子設(shè)備及數(shù)據(jù)傳輸技術(shù)飛速發(fā)展形勢(shì)下,新一代射電天線成為可能。2010年10月,IVS VLBI計(jì)劃執(zhí)行組將部分可升級(jí)的S/X系統(tǒng)升級(jí)為VGOS系統(tǒng),計(jì)劃到2017年全球?qū)崿F(xiàn)16個(gè)VGOS寬帶系統(tǒng)的建設(shè),2019年將實(shí)現(xiàn)30個(gè),其中北半球23個(gè),南半球7個(gè)[3]。

如圖1所示,截至2017年初,全球陸續(xù)已有9臺(tái)符合VGOS系統(tǒng)要求的站址建成,分別為KOKEE,GGAO,WESTFORD,YEBES,WETTZELL和ISHIOKA以及國(guó)內(nèi)的吉林、喀什和三亞,還有11臺(tái)正處于建設(shè)改造中,且從圖1可發(fā)現(xiàn)絕大多數(shù)VGOS站集中于北半球,希望今后能有更多南半球VGOS站加入全球觀測(cè)網(wǎng)中進(jìn)行測(cè)地實(shí)驗(yàn)。

圖1 2017年全球VGOS站分布

1 VGOS系統(tǒng)

中國(guó)科學(xué)院國(guó)家授時(shí)中心和上海天文臺(tái)共同研制了“中科院VGOS系統(tǒng)”(如圖2所示),包括：三亞、吉林、喀什3個(gè)觀測(cè)站,和一個(gè)西安數(shù)據(jù)處理中心,其中喀什站是與新疆天文臺(tái)合作建成,吉林站是與中科院長(zhǎng)春人衛(wèi)站合作建成。各站配備標(biāo)準(zhǔn)的VGOS天線,采用方位俯仰式13 m拋物面天線,主副面板精度分別為0.3 mm和0.1 mm；方位角范圍為-270～+270°,轉(zhuǎn)速和加速度分別為12°/s和3°/s2；俯仰角范圍為5～90°,轉(zhuǎn)速與加速度分別為6°/s和3°/s2；配備寬帶接收機(jī)及寬帶記錄終端,并外接高穩(wěn)定的氫原子鐘,使用饋源制冷、寬帶相關(guān)等技術(shù),來(lái)提高測(cè)量精度。目前西安數(shù)據(jù)處理中心的觀測(cè)綱要編制軟件、三臺(tái)站寬帶VLBI硬件相關(guān)處理機(jī)、DIFX軟件相關(guān)處理機(jī)、測(cè)地與定軌解算軟件都已初步具備服務(wù)能力。

圖2 中科院VGOS系統(tǒng)組成

接收機(jī)是天線系統(tǒng)的重要組成部分,其性能的好壞直接影響到整個(gè)觀測(cè)系統(tǒng)。高靈敏度、高穩(wěn)定度、高可靠性的超寬帶接收機(jī)能夠大大提高系統(tǒng)觀測(cè)能力。饋源由四脊喇叭、耦合器、低噪聲放大器(LNA)、3 dB電橋及連接電纜組成,全部置于制冷杜瓦內(nèi),滿足臺(tái)站系統(tǒng)噪聲溫度小于50 K的要求,然而對(duì)于喀什站,由于電磁環(huán)境差,噪聲溫度在80 K左右。真空泵作為前期輔助工具使杜瓦內(nèi)的真空度低于0.001 pa,后通過(guò)氦氣壓縮機(jī)來(lái)維持杜瓦內(nèi)溫度處于12 K左右,以保障放大器的正常工作環(huán)境：低于70 K；接收機(jī)配備兩級(jí)放大器,左旋、右旋極化信號(hào)首先經(jīng)過(guò)一級(jí)低噪聲放大器,后經(jīng)過(guò)開(kāi)與關(guān)可控的二級(jí)放大器來(lái)切換衛(wèi)星和射電源觀測(cè)模式(衛(wèi)星信號(hào)太強(qiáng),容易引起系統(tǒng)飽和)之后經(jīng)過(guò)下變頻器轉(zhuǎn)換為中頻信號(hào)再由記錄設(shè)備采集數(shù)據(jù)并記錄；系統(tǒng)的10 MHz信號(hào)由氫鐘提供,相位校準(zhǔn)源從信號(hào)接收口注入相位信號(hào),以此拉平設(shè)備造成的相位跳動(dòng),具體可參見(jiàn)圖3。

圖3 VLBI 2010天線接收系統(tǒng)

每個(gè)測(cè)站天線的觀測(cè)頻段為1.2～9 GHz,每站配有4套相同的記錄終端,每一套記錄終端具有32個(gè)通道,每個(gè)通道32 MHz,每套記錄終端最大可記錄兩路16通道的512 MHz,2 bit量化或者一路32通道的1 024 MHz,1 bit量化的寬帶中頻信號(hào)采樣能力,記錄的數(shù)據(jù)首先以Mark5B格式存儲(chǔ)在本地磁盤陣列上,后通過(guò)專網(wǎng)傳回西安數(shù)據(jù)處理中心。

數(shù)據(jù)處理流程如圖4所示。數(shù)據(jù)傳回西安后,各臺(tái)站VLBI數(shù)據(jù)經(jīng)過(guò)10 GE SFP+接口將數(shù)據(jù)高速送入相關(guān)處理機(jī)的數(shù)據(jù)緩存服務(wù)器,經(jīng)分路處理后送入臺(tái)站時(shí)延處理單元的FPGA芯片進(jìn)行條紋旋轉(zhuǎn)、傅里葉變換和小數(shù)比特延遲補(bǔ)償運(yùn)算,然后再同步送入基線單元進(jìn)行交叉相乘得互相關(guān)函數(shù)的譜,并進(jìn)行短時(shí)標(biāo)的累加積分經(jīng)網(wǎng)口按規(guī)定的格式將結(jié)果送存儲(chǔ)服務(wù)器存儲(chǔ)。硬件相關(guān)處理機(jī)在基于觀測(cè)綱要生成的工作文件安排下運(yùn)行,后處理進(jìn)行帶寬綜合,提取觀測(cè)時(shí)延,數(shù)據(jù)分析軟件根據(jù)觀測(cè)時(shí)延數(shù)據(jù),利用最小二乘方法,解算并得到EOP及大氣等參數(shù)。

圖4 數(shù)據(jù)處理流程

2 實(shí)驗(yàn)與數(shù)據(jù)處理

2017年8月20日至8月22日對(duì)中星12號(hào)衛(wèi)星(定點(diǎn)東經(jīng)87.5°)進(jìn)行差分VLBI觀測(cè),射電源觀測(cè)頻率在3.5～4 GHz,衛(wèi)星為3.82 GHz,選取的參考射電源與目標(biāo)衛(wèi)星的角距小于15°[4]。因系統(tǒng)初建成,盡量選取C波段流量在1 Jy以上且在之前實(shí)驗(yàn)中系統(tǒng)能夠得到條紋的射電源,主要有1741-038,2008-159,2134+00,0003-066,0106+013,0420-014,0605-085,0742+103,1055+018,3C273B和1510-089等,具體參見(jiàn)表2。

表2 觀測(cè)射電源簡(jiǎn)介

因射電源的時(shí)延測(cè)量精度與SNR(signal noise ratio)和帶寬成正比,SNR與源流量、觀測(cè)時(shí)長(zhǎng)、帶寬成正比,所以在射電源觀測(cè)中,120 s為一個(gè)周期,射電源跟蹤記錄90 s,天線轉(zhuǎn)動(dòng)30 s,采用512 MHz帶寬2 bit量化,觀測(cè)頻率為3 522～4 034 MHz進(jìn)行實(shí)驗(yàn)觀測(cè)。

各站數(shù)據(jù)經(jīng)專網(wǎng)傳回西安數(shù)據(jù)處理中心后,通過(guò)硬相關(guān)機(jī)處理,得到各通道的頻譜值,后處理從相位信息中提取時(shí)延值,并經(jīng)過(guò)帶寬綜合處理,將3 522～4 034 MHz綜合到3 826 MHz頻率處,從中提取出射電源觀測(cè)群時(shí)延,再結(jié)合三站氣象站采集的溫度、濕度、氣壓數(shù)據(jù),生成分析軟件所需的NGS數(shù)據(jù)格式文件[5],如圖5所示。

圖5 數(shù)據(jù)處理流程

信號(hào)到達(dá)接收機(jī)的觀測(cè)時(shí)延主要包括幾何時(shí)延、對(duì)流層時(shí)延、電離層時(shí)延、天線形變引入的時(shí)延等。在觀測(cè)時(shí)延中扣除對(duì)流層時(shí)延、電離層時(shí)延、天線形變等非幾何時(shí)延后,同理論計(jì)算的幾何時(shí)延相減后進(jìn)行最小二乘求解所需參數(shù),具體流程如圖6所示。

圖6 解算流程圖

由于初次實(shí)驗(yàn),并未深入研究設(shè)備時(shí)延、源結(jié)構(gòu)、天線熱形變等影響,所以在解算中這些誤差源都未扣除,解算時(shí)只扣除了電離層、對(duì)流層時(shí)延,其中電離層采用IGS(International GNSS Service Center)提供的電離層文件插值獲得,對(duì)流層則采用GMF(graphical definition model)模型。具體的時(shí)延模型公式為：

(1)

式(1)中,i表示觀測(cè)數(shù),j表示待解算參數(shù)個(gè)數(shù),Oi為觀測(cè)時(shí)延,vi為殘差,φi(x)為幾何時(shí)延,Ttrop,i為大氣時(shí)延,Tion,i為電離層時(shí)延,Tclock,i為鐘差,Tother,i為設(shè)備、射電源結(jié)構(gòu)等引入的時(shí)延。在不考慮設(shè)備、射電源結(jié)構(gòu)等誤差源的情況下,解算結(jié)果可表示如下：

Toc,i=Oi-φi(x0)-Tion,i-Ttrop,i-Tclock,i,

(2)

(3)

vi=Aδx-Toc,i。

(4)

式(2)～(4)中,Toc,i表示觀測(cè)值與理論值之差,A表示求解參數(shù)的偏導(dǎo)數(shù)矩陣。使用數(shù)據(jù)分析軟件VieVS(vienna VLBI and satellite software),采用單次解模式,射電源坐標(biāo)來(lái)自ICRF2[6],站坐標(biāo)由GPS測(cè)量提供,歲差章動(dòng)采用IAU 2006/2000A模型,EOP初值由IERS公布的EOP C04 08給出,每次觀測(cè)解算估計(jì)以下參數(shù)：各站濕大氣天頂距時(shí)延改正值ΔTzwd,站間鐘差改正值ΔTclock及一個(gè)UT1改正值ΔUT1[7],如式(5)和式(6)所示。

(5)

(6)

式(5)中,mw(ε)為濕大氣映射函數(shù)。由于此次觀測(cè)時(shí)間長(zhǎng),超過(guò)15 h,鐘差與大氣變化快,同國(guó)際上1 h加強(qiáng)觀測(cè)不同,其鐘差估計(jì)間隔為1 440 min,本文需對(duì)濕大氣和鐘差進(jìn)行分段估計(jì)。同時(shí)在解算過(guò)程中進(jìn)行高頻改正,采用拉格朗日方法將EOP初值序列插值到觀測(cè)歷元后,加入海洋潮汐影響和短周期UT1變化的高頻模型修正[8],便于與IERS C04比較。

3 結(jié)果與分析

經(jīng)過(guò)VieVS對(duì)3 d數(shù)據(jù)處理,解算結(jié)果參見(jiàn)表3,表中JL代表吉林站,SY代表三亞站,KS代表喀什站,ΔUT1表示解算結(jié)果與IERS C04 08數(shù)據(jù)文件中(UT1-UTC)參數(shù)值之差。從表3中可知3次的精度基本一致,為0.2 ms左右,且3次實(shí)驗(yàn)的解算結(jié)果在UT1約1σ情況下是吻合的。

表3 實(shí)測(cè)數(shù)據(jù)解算結(jié)果

為驗(yàn)證此次結(jié)果的合理性,采用蒙特卡羅仿真分析方法[9]對(duì)觀測(cè)進(jìn)行仿真測(cè)試,生成各個(gè)觀測(cè)的仿真時(shí)延,并得到仿真時(shí)延與理論時(shí)延之差,具體過(guò)程如圖7所示。仿真參數(shù)主要包括3部分：①鐘；②濕大氣；③白噪聲(表征系統(tǒng)誤差)。鐘的模擬采用Allan標(biāo)準(zhǔn)方差1×10-14@50 min,符合當(dāng)前站上氫鐘的頻率穩(wěn)定特性[10]；因系統(tǒng)不穩(wěn)定因素多,白噪聲設(shè)為30 ps；對(duì)流層模型根據(jù)湍流理論[11]及由Nilsson提出的策略[12]進(jìn)行仿真,三站對(duì)流層參數(shù)在仿真中設(shè)置一致。

圖7 仿真流程圖

因系統(tǒng)誤差影響是隨機(jī)的,所以對(duì)2017年8月20號(hào)的數(shù)據(jù)進(jìn)行30次仿真(3 d觀測(cè)綱要基本一致),得到一個(gè)可信的統(tǒng)計(jì)結(jié)果來(lái)分析比較實(shí)際解算結(jié)果。在理想情況下(無(wú)站坐標(biāo)、電離層誤差,只含鐘及大氣誤差)統(tǒng)計(jì)結(jié)果如圖8所示,圖8(a)表示VieVS解算得到的ΔUT1(UT1-UTC)值與IERS C04 08中事后精密值比較,實(shí)心點(diǎn)表示VieVS值,中間橫線表示IERS值；圖8(b)表示偏差值ΔUT1VieVS-ΔUT1IERS,其偏離IERS(UT1-UTC)參數(shù)平均為0.0 761 ms,方差為0.1 399 ms。

圖8 理想情況下仿真結(jié)果(2017-08-20)

目前三站站坐標(biāo)采用GPS測(cè)量獲得,同國(guó)際上測(cè)地觀測(cè)方式得到的精確站坐標(biāo)存在十幾厘米的誤差,這對(duì)解算結(jié)果中的偏差有很大影響[13],為更好分析其對(duì)結(jié)果影響,在喀什站站坐標(biāo)加入Δx=10 cm,Δy=10 cm,Δz=10 cm偏差,仿真結(jié)果如圖9所示,平均偏差為0.188 7 ms,方差為0.139 9 ms,可見(jiàn)站坐標(biāo)誤差主要影響偏差。

圖9 站坐標(biāo)偏差影響(2017-08-20)

3 d觀測(cè)都是在C波段單頻觀測(cè),不同于國(guó)際上S/X雙頻觀測(cè),它能有效地扣除電離層時(shí)延,本文在解算過(guò)程中則采用IGS提供的電離層網(wǎng)格插值獲得,這勢(shì)必將存在幾個(gè)TECU的誤差,仿真中加入0.1倍電離層時(shí)延到觀測(cè)時(shí)延中,結(jié)果如圖10所示,平均偏差為0.048 3 ms,方差為0.278 9 ms,可見(jiàn)電離層對(duì)方差影響大。

綜上分析,在站坐標(biāo)不精確,無(wú)法有效扣除電離層情況下,經(jīng)仿真可知,該模式下的解算精度在0.2 ms左右,同實(shí)測(cè)結(jié)果在一個(gè)量級(jí)上,可證實(shí)在觀測(cè)數(shù)據(jù)少的情況下,實(shí)測(cè)解算結(jié)果是合理的,同時(shí)也為后續(xù)開(kāi)展相關(guān)工作提供了經(jīng)驗(yàn)。

圖10 電離層影響(2017-08-20)

4 結(jié)語(yǔ)

當(dāng)前中科院首先在國(guó)內(nèi)建立3臺(tái)VGOS天線,并有相應(yīng)的數(shù)據(jù)處理中心；近期利用3臺(tái)天線開(kāi)展了射電源與衛(wèi)星的觀測(cè),通過(guò)上述實(shí)測(cè)數(shù)據(jù)結(jié)果及分析可知,本次解算結(jié)果與IERS事后精密值在≈1σ下基本吻合,這表明三站能獨(dú)立自主地開(kāi)展UT1測(cè)量,這對(duì)授時(shí)單位來(lái)講是非常重要的,同時(shí)在衛(wèi)星定軌精度上也有改進(jìn)。后續(xù)將在C波段開(kāi)展專門的UT1觀測(cè)模式,測(cè)試系統(tǒng)的測(cè)量精度。

致謝:在此感謝上海天文臺(tái)的錢志翰老師,張秀中老師,舒逢春老師,王廣利老師等在觀測(cè)、硬相關(guān)及數(shù)據(jù)處理上的指導(dǎo),并對(duì)合作單位新疆天文臺(tái)與長(zhǎng)春人衛(wèi)站表示感謝,正是在他們的運(yùn)營(yíng)維護(hù)下,基礎(chǔ)設(shè)施才有了保障。