中國VLBI網軟件相關處理機測地應用精度分析

劉 磊，鄭為民,2,3，張 娟，舒逢春,2,3，童鋒賢，童 力

1. 中國科學院上海天文臺, 上海 200030; 2. 中國科學院射電天文重點實驗室，江蘇 南京 210008； 3. 上海市導航定位重點實驗室，上海 200030

?

中國VLBI網軟件相關處理機測地應用精度分析

劉 磊1，鄭為民1,2,3，張 娟1，舒逢春1,2,3，童鋒賢1，童 力1

1. 中國科學院上海天文臺, 上海 200030; 2. 中國科學院射電天文重點實驗室，江蘇 南京 210008； 3. 上海市導航定位重點實驗室，上海 200030

為支持國際VLBI大地測量數據處理，中國VLBI網(CVN)軟件相關處理機完成了功能、性能升級，提高了信噪比，實現了以國際測地通用Mk4格式數據輸出結果，可以直接被VLBI通用測地后處理軟件用于時延數據解算。本文通過與國外軟件相關處理機DiFX的實測數據比對，系統(tǒng)地分析了CVN軟件處理機帶寬綜合殘余時延和時延率精度、帶寬綜合總時延和時延率精度、信噪比和VLBI站坐標解算值。數據顯示，CVN軟件相關處理機已經達到了測地數據處理的精度要求，可以用于IVS國際聯(lián)測數據處理。

VLBI; 軟件相關處理機; 大地測量; 信噪比; 時延; 時延率

甚長基線干涉測量(very long baseline interferom-etry，VLBI)技術誕生于20世紀60年代。得益于氫原子鐘、現代通信、電子計算機等技術的快速發(fā)展以及觀測方法的進步，VLBI測量能力進展迅速，其測量精確度已經達到亞毫角秒量級，是目前角分辨率最高的天文觀測技術[1]，在天體物理[2-3]、天體測量與大地測量(測地)[4]、深空探測[5-10]等多個領域發(fā)揮著重要作用。

VLBI系統(tǒng)通過分布在全球不同位置的VLBI射電望遠鏡(臺站)，在一段時間內連續(xù)地觀測一系列射電源，數據經相關處理后再進行條紋擬合，鐘差鐘速改正，大氣時延及電離層等傳播介質補償，最終得到的觀測量可用于天文、測地及深空探測研究[1]。VLBI系統(tǒng)已被用于建立天球參考架。與基于恒星位置的參考架相比，基于河外射電源的天球參考架因具有更高的穩(wěn)定性，更接近于理想慣性系，被國際天文學聯(lián)合會(IAU)采納為國際天球參考架(International Celestial Reference Frame, ICRF)[11-12]。此參考架確定了天球坐標系統(tǒng)的零點以及方向，并提供了一系列參考點，以確定其他天體在參考架中的精確位置。在測地方面，VLBI技術可用于研究地球運動狀態(tài)變化和確定精確的時間系統(tǒng)，包括進動、章動、極移、世界時等。通過數年至數十年的長時間監(jiān)測臺站坐標的變化，還可分析地球板塊運動[13-14]。

近年來國際VLBI大地測量與天體測量服務組織(IVS)[15-16]提出了滿足VLBI2010技術標準的全球觀測系統(tǒng)VGOS(VLBI2010 global observing system)，以適應未來天體測量與大地測量的高精度需求[17]。VGOS系統(tǒng)的預期目標是：2020年后，將在全球尺度上，實現臺站位置和地球定向參數的連續(xù)監(jiān)測，24 h快速數據處理，VLBI臺站位置精度達到1 mm。這些目標的實現需要對現有VLBI系統(tǒng)進行全面升級，主要包括：采用直徑13 m左右的新型小口徑望遠鏡，通過在單位時間內快速觀測更多的射電源(每個scan時間30 s左右)實現全天區(qū)更好的參考源覆蓋并獲得更高的天線參考點穩(wěn)定性；通過超寬帶接收機數據采集系統(tǒng)，以16～32 Gbps高數據采樣率增加小口徑天線觀測靈敏度。采用實時超寬帶軟件相關處理機、高速e-VLBI和自動化數據處理等關鍵技術以實現長期24 h不間斷觀測。

VLBI測地觀測與數據處理的基本流程是：①VLBI數據中心制定觀測綱要，確定觀測時間、臺站、頻率等；②臺站根據觀測綱要進行觀測，并將數據和日志文件傳到指定VLBI中心；③VLBI中心相關處理機完成相關處理，生成可見度數據，經過相位校正、條紋擬合后，得到VLBI帶寬綜合時延和時延率，并生成專用數據庫格式；④通過測地軟件解算后，獲得臺站坐標、地球定向參數等測量結果。

相關處理是整個VLBI測地數據處理的重要步驟，具有數據密集和計算密集的特點。中國科學院上海天文臺使用自主研發(fā)的CVN(Chinese VLBI network)測地軟件相關處理機，已成功處理了數十次測地數據[18-19]。但是，由于該處理機采用探月工程的專用輸出格式等原因，與國際測地VLBI處理機輸出數據格式不兼容。為支持國際VLBI大地測量數據處理，CVN軟件相關處理機進行了功能、性能升級，實現了輸出結果的國際測地通用Mk4數據格式輸出，提高了信噪比，使處理機結果可直接被通用VLBI測地后處理軟件用于時延數據解算。本文以DiFX(distributed FX correlator)處理機[20-21]作為參照，系統(tǒng)比較了CVN軟件相關處理機測地時延(率)。分析表明，CVN軟件相關處理機測地數據精度與DiFX一致，可以用于IVS國際聯(lián)測的數據處理以及臺站站坐標解算。

1 軟件相關處理機

1.1 VLBI相關處理機原理

FX型(一條基線涉及的兩臺站信號先從時域到頻域做傅里葉變換，再交叉相乘)VLBI相關處理機的輸入是基線兩端臺站獲得的時域信號，輸出是臺站間互相關譜(稱為可見度函數)。由圖1，FX型相關處理機數據處理流程所示，兩路信號進入相關處理機后首先按頻率通道進行解碼；然后根據預先計算的時延模型進行整數比特時延補償；經過條紋旋轉及快速傅里葉變換后，信號從時域變換至頻域，并完成分數比特時延補償；最終兩路信號進行共軛相乘及累加積分，得到可見度數據。為完成精確的測地數據處理，對部分失效數據引入的誤差必須進行精細的數據有效率控制，對在記錄和傳輸過程中缺失的數據也需專門處理，并進行可見度函數幅度校正。

由于光速有限，射電源發(fā)出的信號波前在不同時刻分別到達兩個望遠鏡，二者的時間差被稱為VLBI時延，表示為

(1)

圖1 FX型相關處理機數據處理流程Fig.1 Data processing flow of FX type correlator

1.2CVN軟件相關處理機

根據不同實現方式，VLBI相關處理機可分為軟件處理機和硬件處理機兩類。近年來，基于通用CPU集群的軟件相關處理機因其開發(fā)和調試方便、配置靈活、擴展性強等特點，已經成為VLBI相關處理機的主流形式。

CVN軟件相關處理機最早由中國科學院上海天文臺針對探月工程需求自主開發(fā)，已成功應用于嫦娥1號、2號、3號和探月三期CE-5T1試驗，在一系列的月球探測器測定軌任務中發(fā)揮了關鍵作用。該相關處理機用C語言開發(fā)，運行于Linux集群，采用OpenMP和MPI實現并行，并使用IntelIPP庫進行加速。處理機具備實時相關處理功能，并支持探測器信號的快速條紋搜索。在探月工程中，CVN處理機輸出的可見度數據采用專用CVN格式，包含一個多維矩陣，能夠輸出各臺站、基線的自相關和互相關函數及輔助信息。

為支持“九五”國家重大科學工程——中國地殼運動觀測網絡(陸態(tài)網)，在探月工程相關處理機基礎上研制的CVN測地軟件相關處理機系統(tǒng)[23]，已經成功支持了國內探月工程站坐標測量、陸態(tài)網絡大量VLBI觀測資料處理。2012年后，上海VLBI相關處理機被接納為IVS相關處理機，準備承擔國際觀測數據處理任務。上海VLBI中心從2015年開始采用國外DiFX軟件相關處理機承擔常規(guī)的IVS國際數據處理任務，同時對CVN軟件處理機進行功能與性能升級，使其能夠輸出國際通用數據格式，并提高了輸出結果信噪比。

為支持測地數據的后處理，需要將CVN格式可見度數據結合時延模型、相位校正信號、鐘差鐘速等信息生成國際通用測地的Mk4格式，供測地數據后處理軟件HOPS(Haystackobservatorypostprocessingsystem)讀取并做條紋擬合。Mk4是VLBI處理機測地數據處理的輸出格式，由美國Haystack天文臺提出并作為HOPS的輸入格式。Mk4格式由一系列具有相同后綴的文件組成，有5種類型：

(1)root文件——對應scan的整體信息，包括源名稱，觀測臺站，觀測時間，頻率通道，等等。

(2)type-1文件——相關處理機輸出的可見度數據。每條基線對應一個文件。

(3)type-2文件——條紋擬合結果，包括不同頻段的帶寬綜合時延，單通道時延，時延率，信噪比等信息。每個type-1文件可以對應多次擬合，根據擬合序號被保留為多個type-2文件。

(4)type-3文件——臺站信息，包括相位校正信號、時延模型，等等。

(5)type-4文件——日志信息，可選。

圖2顯示了CVN處理機生成Mk4測地格式數據示意圖。處理機采用專用軟件讀取輸出的可見度數據數組、job文件和vex文件，以及可選的相位校正信號文件，生成root文件、type-1和type-3文件。HOPS軟件包中的fourfit程序讀取Mk4格式進行條紋擬合，將條紋擬合結果，包括殘余時延(率)、總時延(率)、信噪比等輸出至type-2文件。

圖2 CVN處理機輸出Mk4格式示意圖Fig.2 Demonstration of CVN correlator output to Mk4 format

1.3DiFX軟件處理機

為便于判斷CVN處理機結果正確性，采用國外DiFX處理機作為對比的參考對象。DiFX是目前世界上應用最廣泛的VLBI軟件相關處理機，由澳大利亞Swinburne大學、美國NRAO等聯(lián)合開發(fā)，支持天體物理、測地以及深空探測等多種應用。DiFX采用C++語言，高度模塊化，運行于Linux集群，使用IntelIPP進行加速。上海VLBI中心已于2015年建立了基于DiFX的天文與測地VLBI數據相關處理平臺，目前已成功處理近30次IVS觀測數據，得到了IVS組織認可[24]?；诂F有的測地觀測數據，通過將CVN和DiFX處理結果在不同后處理階段進行比較，可分析驗證CVN處理機的輸出精度，為CVN處理機后續(xù)進行IVS數據處理做好準備。

2 CVN軟件相關處理機與DiFX輸出精度比較

(2)

帶寬綜合時延是條紋擬合時，分別綜合S頻段與X頻段所有頻率通道數據擬合出的時延，較準確地反映了S頻段與X頻段的群時延，在參數解算時作為觀測量使用。本文使用HOPS中的alist程序讀取條紋擬合結果，并從中提取每個scan每條基線擬合得到的帶寬綜合殘余時延和時延率、帶寬綜合總時延、時延率和信噪比等信息。通過比較同一批觀測數據經過兩個相關處理機后的擬合結果的差異** 本節(jié)所討論的時延、時延率、信噪比差異含義為對應scan，波段和基線CVN結果減去DiFX結果的差值。，驗證了CVN軟件相關處理機的精度。在此基礎上，本文通過一次完整的測地觀測數據處理過程，比較了基于兩個處理機相關數據的VLBI臺站坐標解算結果，二者在誤差范圍內一致。

2.1 帶寬綜合殘余時延和殘余時延率精度比較

帶寬綜合殘余時延和殘余時延率由條紋擬合程序根據可見度數據擬合得到，其擬合結果僅依賴于處理機的相關處理結果，是衡量處理機數據質量最直接的指標。在輸入時延模型相同的情況下，兩個處理機的帶寬綜合殘余時延和殘余時延率可以直接進行比較。比較時采用了2014年IVS聯(lián)測數據，觀測代碼為k14349，共有挪威、中國(上海)、日本、德國4個臺站(Ny，Sh，Ts，Wz)，16個頻率通道。每個通道帶寬8MHz；觀測時長1h，共30個scan。其中Wz站記錄格式為Mark5a格式，其余3站為Mark5b格式。

本次處理使用的DiFX版本為2.4.1，其調用的IntelIPP庫版本為7.0。為精確比較殘余時延，兩個處理機需采用同樣的時延模型。為此，先將DiFX的時延模型替換為CVN的時延模型，然后再進行相關處理。DiFX的im文件中每個scan提供數個模型，彼此銜接，長度均為2min。CVN處理機每個scan提供一個模型，起始終止時刻與scan重合。替換時保持DiFX模型起始終止時刻不變，將CVN處理機的5次多項式系數變換到每個DiFX模型開始時刻。相關處理時兩處理機采用相同的設置，包括FFT點數、積分時間、鐘差鐘速等。在對兩處理機的輸出結果進行條紋擬合時，所有選項均采用fourfit默認設置，即采用相位校正信號自動調整和不剔除壞頻率通道。在進行比較時，剔除了所有帶有錯誤代碼以及數據質量(QC，qualitycode)低于5的數據，并對兩處理機的結果按照scan和基線進行匹配。此處的“數據質量”(QC)為fourfit程序對條紋擬合結果做的一個評估，分為0—9共10個等級(數字越高擬合質量越好)以及若干表示條紋擬合錯誤的字母。測地解算時通常采用數據質量不小于5的結果。

圖3是帶寬綜合殘余時延差異隨信噪比的變化。虛線和實線分別對應著X頻段和S頻段帶寬綜合殘余時延彌散的理論值，由如下公式給出[2]

(3)

式中，SSNR是信噪比；frms是對應頻段的各頻率通道帶寬均方根，根據k14349的頻率配置，X頻段和S頻段分別為298.3MHz和50.9MHz。由于S頻段相比于X頻段頻率通道帶寬均方根較小，相同信噪比條件下，其對應的帶寬綜合時延彌散較大。由圖中可以看出兩個處理機的帶寬綜合殘余時延差異隨著信噪比的增加而減小，并遠小于理論彌散值。圖中，X頻段和S頻段的帶寬綜合殘余時延的加權平均差異分別為0.02ps和0.74ps；彌散或加權均方根(weightedrootmeansquare，WRMS)分別為0.77ps和9.35ps。加權平均差異遠小于加權均方根，說明數據圍繞0值波動，無顯著系統(tǒng)性偏差。S頻段加權均方根差異超過9ps。此差異源于兩處理機實現細節(jié)的不同：Wz站數據采用Mark5a格式數據，DiFX輸出的最后一個積分周期結束時刻在scan結束時刻之后，而CVN軟件相關處理機所有積分周期均在scan范圍內，從而導致DiFX在處理某些scan數據時會比CVN處理機多輸出一個積分周期，條紋擬合時結果略有不同。

圖3 k14349帶寬綜合殘余時延差異隨信噪比變化Fig.3 Difference of multi-band residual delay for k14349 as a function of signal to noise ratio

圖4是殘余時延率隨信噪比的變化，同樣可以看出差異隨著信噪比的增加而減小，S頻段差異比X頻段要大，但是大部分在0.1ps/s之內。

圖4 k14349殘余時延率差異隨信噪比變化Fig.4 Difference of multi-band residual delay rate as a function of signal to noise ratio

圖5展示了k14349觀測CVN處理機相對于DiFX處理機結果信噪比的差異，以百分比表示。二者大部分數據信噪比相對差異在0.5 %以內。少量數據CVN結果相較于DiFX低約2 %。信噪比稍較低同樣是由于某些情況下DiFX會多輸出一個積分周期導致其對應條紋擬合結果信噪比較CVN結果稍高。

CVN結果相對于DiFX結果的差異圖5 k14349信噪比的相對差異Fig.5 Difference of signal to noise ratio for k14349

2.2 帶寬綜合總時延和總時延率精度比較

可見度數據經條紋擬合后，得到殘余時延(率)，通過將其與對應時刻模型時延(率)，以及鐘差(速)相加，可以得到總時延(率)。本節(jié)比較所采用的原始數據為CVN組織的烏魯木齊南山站(Ur)改造后的站坐標測定試驗，觀測代碼為cn1502。北京密云、上海佘山、云南昆明和烏魯木齊南山4個CVN臺站參加(代號Bj，Sh，Km，Ur)，采用S/X雙頻共16個頻率通道，每通道帶寬8MHz；觀測時長24h，共227個scan。觀測數據分別用DiFX和CVN軟件相關處理機處理，各自采用獨立的時延模型。DiFX版本為2.3，其調用的IntelIPP庫版本為7.0，模型計算使用Calc9.1。CVN模型計算程序由上海天文臺研發(fā)，并與國際上的主流模型進行了比較[25]。相關處理時兩處理機采用相同的設置，包括FFT點數、積分時間、鐘差鐘速等。處理機結果采用fourfit做條紋擬合，兩組數據剔除的相同的壞頻率通道。

由于Bj站沒有注入相位校正信號，各臺站相位校正均采用手動模式，并對兩處理機數據對應通道采用相同的校正值。這樣做的原因是，由于兩處理機所采用不同的時延模型導致殘余條紋的差異，基于兩組數據各自進行相位校正會導致一個恒定的相位差異，進而導致條紋擬合時各基線10ns量級的恒定時延差異。該恒定時延差異完全由于相位校正時不同的初相造成，不影響解算結果。為比較時延的彌散情況，擬合采用了相同的相位校正值，以獲得同樣的初始相位。本次觀測4個臺站均采用Mark5b格式記錄，不會出現上節(jié)中DiFX處理機某些scan比CVN處理機多輸出一個積分周期的情況。由于開發(fā)環(huán)境、任務劃分策略等算法實現細以及所采用時延模型的不同，兩個處理機的擬合結果以及臺站坐標解算結果不完全相同，但分析顯示結果在精度范圍內一致。

圖6顯示了cn1502各基線精度差異。每幅圖從上到下依次是帶寬綜合總時延差異、總時延率差異、信噪比差異。由圖可知，WRMS與信噪比呈反相關關系，信噪比較低時WRMS較大。在本次觀測中，Ur站由于X頻段饋源兩個IF的極化均接反，導致Ur站與極化正確臺站的相關基線在X頻段上信噪均較低，普遍在20以下((f)、(j)圖)，但Sh站X頻段兩個IF的一個極化也接反，恰好與Ur站一致，導致對應Sh-Ur基線X頻段信噪比反而較高，總的效果介于全部接反和全部正確之間(圖l)。除Ur站以外的其余基線X頻段信噪比較高，普遍高于20，對應的WRMS均在3ps以內((b)、(d)、(h)圖)。數據表明，兩個處理機計算的總時延率符合很好，差異通常在0.1ps/s以內。

圖6 cn1502各基線精度差異比較Fig.6 Precision comparison of each baseline for cn1502

各分圖中從上到下依次為帶寬綜合總時延差異、總時延率差異和信噪比圖6 cn1502各基線精度差異比較(續(xù))Fig.6 Precision comparison of each baseline for cn1502(continude)

圖7展示了cn1502兩處理機結果數據信噪比的相對差異。由于參與相關的4個臺站均為Mark5b格式，兩處理機每個scan輸出有效積分周期個數基本相同，相較于k14349，信噪比的差異更小，基本集中在0.5%的范圍內。

圖7 cn1502信噪比的相對差異Fig.7 Difference of signal to noise ratio for cn1502

2.3 臺站坐標解算精度比較

在完成cn1502后處理之后，將其導入vgosDB數據庫進行Ur臺站坐標解算。本次解算使用了最新的測地VLBI后處理軟件nuSolve0.2.2版本，由NASA/GSFC開發(fā)[26]。nuSolve解算的算法基于廣泛使用的傳統(tǒng)測地軟件CALC/SOLVE。由于借鑒了現代軟件工程的開發(fā)方法，更加注重用戶體驗，nuSolve安裝和使用均比CALC/SOLVE便捷。針對兩個處理機的輸出結果，nuSolve附帶的vgosDbMake程序將條紋擬合得到的Mk4格式轉換為vgosDB版本1數據庫；vgosDbCalc程序調用Calc11生成vgosDB版本2數據庫；vgosDbProcLogs程序讀取4個臺站的log文件，生成vgosDB版本3數據庫。數據庫準備好之后，使用nuSolve將其分別導入，進行相同的數據操作，基本流程是：

(1) 選擇要導入的數據庫。

(2) 切換到單通道時延(SBdelay)，剔除壞數據點，逐基線調整跳秒。

(3) 切換到群時延(GRdelay)，解模糊度(ambig)。

(4) 解電離層改正值(IonoC)。

(5) 解大氣時延(zenithdelay)和臺站坐標(station)。兩組數據均選擇Sh站為參考站。

(6) 解基線鐘差(baselineclocks)。兩組數據均選擇Bj站為參考站。

(7) 解時間變化率(dUT1rate)、章動(nutation)。

(8) 做權重調整(reweighting)。

(9) 做數據修整(outliner)，首先剔除(elimination)模式，再恢復(restoration)模式，然后重復，直到沒有數據被剔除。

(10) 解算完成，生成解算報告。

對兩個處理機的對應數據根據上述步驟進行了解算、評估(見表1)。解算過程中所涉及的觀測量均為時延，其中有效觀測量是指通過fourfit進行條紋擬合得到的數據質量在5—9之間的數據。表1顯示，兩組數據的觀測量數目及擬合質量一致。由前文可知，兩處理機條紋擬合結果差異的均方差約10ps，小于解算結果的均方差(40ps)，即CVN處理機輸出結果的條紋擬合精度可以滿足解算的需要。本次解算主要關心Ur站坐標，兩組數據的解算結果X、Y、Z、U、E、N分量結果分別羅列于表2中。解算結果表明，在精度范圍內，兩組數據對于Ur站坐標的解算結果無顯著差異。

表1 cn1502解算結果評估參數

表2 cn1502 解算Ur站坐標結果

3 結 論

為支持VLBI國際測地數據處理，CVN軟件相關處理機進行了升級改造。通過對CVN和DiFX軟件相關處理機的測地觀測結果的條紋擬合及站坐標解算，分析了CVN軟件相關處理機測地應用精度。結果表明，兩處理機輸出的帶寬綜合殘余時延的加權平均差異無系統(tǒng)性偏差，S頻段加權均方差在10 ps之內，X頻段加權均方差在1 ps以內。S頻段帶寬綜合總時延差異的加權均方差在9 ps以內，X頻段在信噪比較高時(大于20)加權均方差在3 ps以內。殘余時延率和總時延率差異通常在0.1 ps/s以內。信噪比差異通常在0.5%以內。兩個處理機相關結果分別經參數解算得到的烏魯木齊南山站坐標一致，差異在誤差范圍內。綜合對比試驗數據可知，CVN軟件相關處理機已經達到了測地數據處理要求，具備處理IVS國際聯(lián)測數據的能力。

[1] THOMPSON A R, MORAN J M, SWENSON G W JR. Interferometry and Synthesis in Radio Astronomy[M]. New York: Wiley-Interscience Publication, 2001.

[2] FISH V L, JOHNSON M D, DOELEMAN S S, et al. Persistent Asymmetric Structure of Sagittarius A on Event Horizon Scales[J]. The Astrophysical Journal, 2016, 820(2): 90.

[3] WALTER F, BRINKS E, DE BLOK W J G, et al. THINGS: The H I Nearby Galaxy Survey[J]. The Astrophysical Journal, 2008, 136(6): 2563.

[4] 喬書波, 李金嶺, 柴洪洲, 等. 天球參考架的穩(wěn)定源選取及其最新實現[J]. 測繪學報, 2010, 39(2): 120-128. QIAO Shubo, LI Jinling, CHAI Hongzhou, et al. Newly Realization of Celestial Reference Frame by Stable Sources Selection[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(2): 120-128.

[5] 喬書波, 李金嶺, 孫付平. VLBI在探月衛(wèi)星定位中的應用分析[J]. 測繪學報, 2007, 36(3): 262-268. DOI: 10.3321/j.issn:1001-1595.2007.03.004. QIAO Shubo, LI Jinling, SUN Fuping. Application Analysis of Lunar Exploration Satellite Positioning by VLBI Technique[J]. Acta Geodaetica et Cartographica Sinica, 2007, 36(3): 262-268. DOI: 10.3321/j.issn:1001-1595.2007.03.004.

[6] 錢志瀚, 李金嶺. 甚長基線干涉測量技術在深空探測中的應用[M]. 北京: 中國科學技術出版社, 2012. QIAN Zhihan, LI Jinling. The Application of Very Long Baseline Interferometry in Deep Space Navigation[M]. Beijing: China Science and Technology Press, 2012.

[7] 陳世杰, 杜蘭, 旦增曲英, 等. 航天器姿軌控始末時刻的動力學標定方法[J]. 測繪學報, 2013, 42(6): 824-830. CHEN Shijie, DU Lan, DANZENG Quying, et al. Dynamic Detection Algorithm for Spacecraft Attitudinal Control and Orbital Control Epoch[J]. Acta Geodaetica et Cartographica Sinica, 2013, 42(6): 824-830.

[8] 童鋒賢, 鄭為民, 舒逢春. VLBI相位參考成像方法用于玉兔巡視器精確定位[J]. 科學通報, 2014, 59(34): 3362-3369. TONG Fengxian, ZHENG Weimin, SHU Fengchun. Accurate Relative Positioning of Yutu Lunar Rover Using VLBI Phase-referencing Mapping Technology[J]. Chinese Science Bulletin, 2014, 59(34): 3362-3369.

[9] 周歡, 童鋒賢, 李海濤, 等. 深空探測器同波束相位參考成圖相對定位方法[J]. 測繪學報, 2015, 44(6): 634-640. DOI: 10.11947/j.AGCS.2015.20140240. ZHOU Huan, TONG Fengxian, LI Haitao, et al. Relative Position Determination between Deep-space Probes Based on Same Beam Phase-referencing Imaging Technique[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(6): 634-640. DOI: 10.11947/j.AGCS.2015.20140240.

[10] 朱柯, 鄭為民, 童鋒賢, 等. 深空探測器VLBI成圖技術研究[J]. 中國科學: 物理學 力學 天文學, 2016, 46(8): 089513. ZHU Ke, ZHENG Weimin, TONG Fengxian, et al. A Study of VLBI Imaging Techniques for Deep Space Probes[J]. Scientia Sinica: Physica, Mechanica & Astronomica, 2016, 46(8): 089513.

[11] MA C, ARIAS E F, EUBANKS T M, et al. The International Celestial Reference Frame as Realized by Very Long Baseline Interferometry[J]. The Astronomical Journal, 1998, 116(1): 516.

[12] MA C, ARIAS E F, BIANCO G, et al. The Second Realization of the International Celestial Reference Frame by Very Long Baseline Interferometry[R]. IERS Technical Note No. 35. Frankfurt am Main: Verlag des Bundesamts für Kartographie und Geod?sie, 2009.

[13] 李金嶺, 王廣利. 上海VLBI站相對于歐亞板塊穩(wěn)定區(qū)域的運動[J]. 中國科學院上海天文臺年刊, 2001(22): 57-62. LI Jinling, WANG Guangli. The Relative Motion of Shanghai VlbI Station to the Stable Part of Eurasia Plate[J]. Annals of Shanghai Observatory Academia Sinica, 2001(22): 57-62.

[14] FEIGL K L, AGNEW D C, BOCK Y, et al. Space Geodetic Measurement of Crustal Deformation in Central and Southern California, 1984-1992[J]. Journal of Geophysical Research, 1993, 98(B12): 21677-21712.

[15] BEHREND D. Data Handling within the International VLBI Service[J]. Data Science Journal, 2013, 12: WDS81-WDS84.

[16] SCHUH H, BEHREND D. VLBI: A Fascinating Technique for Geodesy and Astrometry[J]. Journal of Geodynamics, 2012(61): 68-80.

[17] PETRACHENKO W, BEHREND D, HASE H, et al. The VLBI2010 Global Observing System (VGOS)[C]∥General Assembly European Geosciences Union. Vienna, Austria: EGU, 2013.

[18] SHU F, ZHENG W, ZHANG X, et al. Shanghai Correlation System Upgrade for Geodetic Application[C]∥Proceedings of the 19th European VLBI for Geodesy and Astrometry Working Meeting. Bordeaux, Frane: VLBI, 2009.

[19] WANG Guangli, XU Minghui, ZHANG Zhibin, et al. The CVN Geodetic Observation and Its Result[C]∥Proceedings of the 22nd European VLBI Group for Geodesy and Astronomy Working Meeting. Azores, Portugal: EVGA, 2015.

[20] DELLER A T, TINGAY S J, BAILES M, et al. DiFX: A Software Correlator for Very Long Baseline Interferometry Using Multiprocessor Computing Environments[J]. Publications of the Astronomical Society of the Pacific, 2007, 119(853): 318-336.

[21] DELLER A T, BRISKEN W F, PHILLIPS C J, et al. DiFX-2: A More Flexible, Efficient, Robust, and Powerful Software Correlator[J]. Publications of the Astronomical Society of the Pacific, 2011, 123(901): 275-287.

[22] TAKAHSASHI F, KONDO F, TAKAHASHI Y, et al. Very Long Baseline Interferometer (Wave Summit Course)[M]. Tokoyo: Ohmsha Ltd., 2000.

[23] ZHENG Weimin, QUAN Ying, SHU Fengchun, et al. The Software Correlator of the Chinese VLBI Network[C]∥Proceedings of the Sixth International VLBI Service for Geodesy and Astronomy. Hobart, Tasmania, Australia: VLBI, 2010: 157-161.

[24] SHU Fengchun, ZHENG Weimin, JIANG Wu, et al. Shanghai VLBI Correlator[C]∥International VLBI Service for Geodesy and Astrometry 2014 Annual Report. 2015: 171.

[25] KLOPOTEK G, ARTZ T, BELLANGER A, et al. Results from the VLBI Analysis Software Comparison Campaign 2015[C]∥International VLBI Service for Geodesy and Astrometry 2016 General Meeting. Johannesburg, South Africa: VLBI, 2016.

[26] BOLOTIN S, BAVER K, GIPSON J, et al. Current Status of nuSolve[C]∥Proceedings of the 21st Meeting of the European VLBI Group for Geodesy and Astronomy. Espoo, Finland: VLBI, 2013: 77-80.

(責任編輯：叢樹平)

Precision Analysis of Chinese VLBI Network Software Correlator for Geodetic Applications

LIU Lei1,ZHENG Weimin1,2,3,ZHANG Juan1,SHU Fengchun1,2,3,TONG Fengxian1,TONG Li1

1. Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China; 2. Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, Nanjing 210008, China; 3. Shanghai Key Laboratory of Space Navigation and Positioning Techniques, Shanghai 200030, China

To support international VLBI geodetic data processing, Chinese VLBI Network (CVN) software correlator has been updated for better performance and more features, including geodetic Mk4 data format output and improvement of signal to noise ratios. After updating, the visibility output of the software correlator can be sent to the standard VLBI geodetic postprocessing software to derive observables. To verify the correlation result, we have carried out correlation with both CVN software correlator and the popular DiFX correlator by using the same data set. We have compared the derived multi-band residual delay/rate, total delay/rate, signal to noise ratio, and the solved VLBI site positions. The comparison results demonstrate that the CVN software correlator has achieved the required precision, and therefore it can be used for IVS geodetic data processing.

VLBI; software correlator; geodesy; signal to noise ratio; delay; delay rate

The National Natural Science Foundation of China (Nos. 11373061; 11303077; 11573057); The Shanghai Program of Shanghai Academic Research Leader (No. 14XD1404300); The Shanghai Outstanding Academic Leaders Plan; CAS Key Technology Talent Program; Ten Thousand Talent Program

LIU Lei (1984—), male, PhD, assistant researcher, majors in VLBI geodetic data processing, development of VLBI software correlator.

ZHENG Weimin

劉磊，鄭為民，張娟，等.中國VLBI網軟件相關處理機測地應用精度分析[J].測繪學報，2017,46(7)：805-814.

10.11947/j.AGCS.2017.20160600. LIU Lei,ZHENG Weimin,ZHANG Juan,et al.Precision Analysis of Chinese VLBI Network Software Correlator for Geodetic Applications[J]. Acta Geodaetica et Cartographica Sinica,2017,46(7):805-814. DOI:10.11947/j.AGCS.2017.20160600.

P228

A

1001-1595(2017)07-0805-10

國家自然科學基金(11373061; 11303077; 11573057); 上海市優(yōu)秀學術帶頭人項目(14XD1404300); 上海市領軍人才項目; 中國科學院關鍵技術人才項目; 國家萬人計劃

2016-11-22

劉磊(1984—)，男，博士，助理研究員，研究方向為VLBI測地數據處理與VLBI軟件相關處理機開發(fā)。

E-mail： liulei@shao.ac.cn

鄭為民

E-mail： zhwm@shao.ac.cn

修回日期： 2017-03-21