TEQC在鐵路工程測量中的應(yīng)用

湯偉堯

(中鐵第一勘察設(shè)計院集團有限公司，陜西西安 710043)

科學技術(shù)的發(fā)展日新月異，空間信息技術(shù)更是朝著多元化、實用化的方向不斷發(fā)展。關(guān)于海量空間信息的獲取手段，目前發(fā)展較為成熟的系統(tǒng)包括我國已基本建設(shè)完成的北斗全球定位系統(tǒng)(BDS)、美國GPS系統(tǒng)、俄羅斯GLONASS以及歐盟GALILEO系統(tǒng)。這些GNSS全球?qū)Ш蕉ㄎ幌到y(tǒng)都為用戶提供了海量、高精度、高分辨率的時空信息，并且已經(jīng)成功應(yīng)用于各個行業(yè)和領(lǐng)域。其中，在鐵路建設(shè)的勘察設(shè)計、施工運營等階段應(yīng)用較為廣泛。

關(guān)于多源海量空間位置信息，其全生命周期一般經(jīng)過數(shù)據(jù)采集、數(shù)據(jù)預(yù)處理和數(shù)據(jù)處理與分析三個重要過程。其中，數(shù)據(jù)預(yù)處理至關(guān)重要(主要包括數(shù)據(jù)編輯、加工和整理)。預(yù)處理工作是否完善，對后續(xù)數(shù)據(jù)平差計算以及平差結(jié)果的精度將產(chǎn)生重要的影響[1]。目前，鐵路工程測量中出現(xiàn)GNSS數(shù)據(jù)質(zhì)量較差的情況一般都返工重測，而對如何提高其利用率的探討研究較少。以下將研究利用預(yù)處理軟件TEQC的強大功能，探討并分析測量過程中出現(xiàn)短暫信號丟失情況時如何對數(shù)據(jù)進行后處理，從而避免返工重測，提高生產(chǎn)效率。

1 TEQC功能介紹

TEQC是UNAVCO Facility組織開發(fā)的用于GNSS觀測信息預(yù)處理的軟件[2]，相較于其他預(yù)處理軟件，其功能較多，使用命令簡單易懂，既可以檢核分析雙頻接收機的靜態(tài)觀測信息，也可以用于研究動態(tài)數(shù)據(jù)質(zhì)量。其原理主要是利用GNSS觀測信息的偽距和載波相位信息進行線性組合，從而實現(xiàn)對GNSS數(shù)據(jù)的誤差估計[3]。TEQC主要功能包括：格式轉(zhuǎn)換(Translate)、數(shù)據(jù)編輯(Edit)、質(zhì)量檢核(Quality Check)和單點定位(Coordinate)[4]。以下主要研究TEQC的數(shù)據(jù)編輯和質(zhì)量檢核功能。

1.1 質(zhì)量檢核

質(zhì)量檢核一般有l(wèi)ite模式和full模式(常用模式為full)。full模式：需要o文件和n文件，命令為teqc+qc -nav source. *.*n source.*.*o，執(zhí)行后文件夾里將產(chǎn)生9個文件。其中*.*S是分析結(jié)果匯總文件；*.ion是L1/L2電離層延遲；*.iod是電離層延遲的變化率；*.mp1是L1載波上的碼信息多路徑觀測誤差；*.mp2是L2載波上的p碼多路徑觀測誤差；*.sn1是L1觀測值的信噪比；*.sn2是L2觀測值的信噪比；*.ele是衛(wèi)星高度角；*.azi是衛(wèi)星方位角。

對數(shù)據(jù)質(zhì)量分析主要包括四點。

(1)數(shù)據(jù)完整性：數(shù)據(jù)完整性代表了實際觀測歷元和理論預(yù)計歷元的比值(possible obs/complete obs)，也就是預(yù)期歷元數(shù)和實際歷元數(shù)的比值。

(2)多路徑誤差分析：MP1和MP2主要用于分析GNSS接收機周圍環(huán)境是否造成多路徑效應(yīng)，進而評價觀測環(huán)境的質(zhì)量。其計算公式為

(1)

(2)

(3)

(4)

(5)

上式中，P1和P2分別代表L1和L2載波上的測距碼；φ1和φ2顯示載波相位觀測值；M1和M2是雙頻偽距的多路徑效應(yīng)；m1和m2是雙頻載波相位的多路徑效應(yīng)；n1和n2是整周模糊度；f1和f2是頻率；λ1和λ2是波長；假如沒有產(chǎn)生周跳，B1和B2是常量，m1、m2遠小于M1和M2；α是L1、L2波段頻率之比的平方。

(3)信噪比(SNR)：信噪比主要用于分析接收信息與高度角變化的關(guān)系，還可以分析評價接收機對高度角較低衛(wèi)星信號的追蹤接收能力。

(4)觀測歷元數(shù)與周跳數(shù)的比值(o/slps)：經(jīng)常用CRS=1 000(o/slps)表示。

(5)電離層延遲微分周跳(IOD)：一般用這個指標來衡量相位觀測值中是否存在周跳，如果IOD變化率大于400 cm/min，應(yīng)對其進行重點關(guān)注，有

(6)

根據(jù)IGS跟蹤站的數(shù)據(jù)統(tǒng)計，一般要求MP1<0.5，MP2<0.75，o/slps>100[5]。

1.2 數(shù)據(jù)編輯

常見的合并功能應(yīng)用于將多日數(shù)據(jù)分割成單日文件，也可以將單時段數(shù)據(jù)進行合并，形成一個長時間序列的觀測文件。通過給定一個時間范圍，可以對標準格式數(shù)據(jù)進行該時間窗口的切割。例如：teqc-dm 20 source.**o > result.**o表示截取觀測文件最后20 min數(shù)據(jù)。如果合并兩個時間段的RINEX數(shù)據(jù)，可以使用teqc urum.135？.10o>urum.1350.10o。

2 TEQC工程應(yīng)用實例

以2019年8月3日某鐵路精密工程測量實驗為例。平面坐標系采用2000國家大地坐標系，長半軸為6 378 137，扁率為298.257 222 101。按照高斯窄帶投影的方法建立工程獨立坐標系，滿足規(guī)范中投影變形不大于10 mm/km的要求[6-8]。本實驗研究范圍中央子午線為103°，投影面大地高為2 610 m。截止衛(wèi)星高度角設(shè)置為15°，數(shù)據(jù)采樣率為15 s，PDOP小于等于6，觀測時間長度為120 min。其中， CPI57-1點在一個時段中陸續(xù)停電3次，測試具體概況如表1和表2所示。

表2 CPI57-1四個時段概況

表1 CPI57-1四個時段數(shù)據(jù)質(zhì)量檢核

執(zhí)行teqc+qc -nav CPI02150.19o CPI02150.19n，得到每一時段的匯總報告文件。其中包含的參數(shù)信息有每顆衛(wèi)星的接收數(shù)據(jù)情況、觀測時間、首歷元、末歷元、MP1(L1載波的多路徑影響)、MP2(L2載波的多路徑影響)、SN1(L1載波的信噪比)、SN2(L2載波的信噪比)、IOD(電離層延遲變化率)、ION(電離層延遲誤差)、AZI(方位角)、ELE(衛(wèi)星高度角)。由表1可以看出，四個時段的衛(wèi)星數(shù)都多于10顆，觀測條件良好；根據(jù)一般性原則，重點關(guān)注MP1和MP2，如果其值分別小于0.35和0.45，說明站點周圍的多路徑影響效應(yīng)較弱[9-12]，本研究中的MP1和MP2遠小于限差值，數(shù)據(jù)有效率都在90%以上，說明多路徑、遮擋以及電磁干擾等影響非常小。

由表2可以看出， CPI57-1點第一時段和第二時段相差只有11 s，由于靜態(tài)測量中設(shè)置的采樣率為15 s，所以斷電小于15 s的情況下觀測數(shù)據(jù)的損失很小。第二時段和第三時段間斷42 s，丟失3個歷元的數(shù)據(jù)，其影響也相對較小。第三個時段和第四個時段間斷時間相對較長(共有3 min57 s的間斷時間)，相對觀測數(shù)據(jù)損失較多。鑒于觀測過程中間斷時間短暫，采樣率、儀器高等參數(shù)不變，故利用TEQC的文件合并功能，執(zhí)行teqc file1 file2>myfile，將以上四個時段數(shù)據(jù)合并為一個文件。

為了研究合并后數(shù)據(jù)的可用性，將CPI57-1點與其他點組成基線，并且同步觀測6次，以進行對比分析。觀測數(shù)據(jù)的基線解算采用徠卡隨機商用軟件LGO8.4。重復(fù)基線、閉合環(huán)檢驗、整網(wǎng)平差計算利用武漢大學研制的COSAGPS軟件。圖1為CPI57-1與其他點位的位置分布。

圖1 CPI57-1與其他點位位置關(guān)系

2.1 異步環(huán)統(tǒng)計

根據(jù)工程測量規(guī)范及要求，外業(yè)觀測結(jié)束后，首先對基線進行處理和質(zhì)量分析，然后對重復(fù)基線、環(huán)閉合差進行檢驗，獨立觀測過程中的異步環(huán)坐標分量也要符合相關(guān)規(guī)范中的技術(shù)指標[8]，如表3所示。

表3 基線向量環(huán)閉合差限差

對6次獨立平差的基線向量網(wǎng)獨立環(huán)閉合差進行計算，其獨立環(huán)最大閉合差統(tǒng)計結(jié)果如表4所示。

表4 閉合環(huán)最大閉合差 mm

由表4可以看出，與CPI57-1組成的基線向量網(wǎng)中，各獨立環(huán)的最大閉合差在X、Y、Z以及基線長度分量上都能滿足限差要求，說明合并后的數(shù)據(jù)質(zhì)量相對較好。

2.2 重復(fù)基線較差統(tǒng)計

表5 重復(fù)基線最大較差 mm

2.3 三維自由網(wǎng)平差

圖2 同一組基線不同時段三維基線向量殘差對比

由圖2可以看出， CPI57-1參與構(gòu)成的基線向量6個時段的基線向量殘差互差都比較小。其中，CPI57-1—CPI57-2規(guī)范要求限差為1.51 cm，CPI57-1—CPI57-3限差為11.51 cm，CPI57-1—CPI059限差為5.6 cm，CPI57-1—CPI060限差為5.79 cm。與CPI57-1有關(guān)的三維基線向量殘差都小于限差。

2.4 二維約束平差

選擇相同的起算點，得到6個時段CPI57-1的二維約束平差值，如表6和圖3所示。

表6 不同時段二維約束平差坐標 m

圖3 不同時段二維約束成果對比

由表6和圖3可知，將6個時段聯(lián)合平差的結(jié)果作為參考值，將每個時段的二維約束平差值與其進行對比，可以看出第一時段和第四時段解得的CPI57-1坐標差相對較大。由于第一時段的數(shù)據(jù)為多個斷開時段合并的數(shù)據(jù)，多余觀測較少，且第一時段觀測過程中該點位周邊車輛和人員來往較密集，致使觀測點位穩(wěn)定性較差，故第一時段CPI57-1點位精度較差。第四時段可能是由于數(shù)據(jù)采集過程中多路徑影響或者衛(wèi)星星況較差造成的：在該區(qū)域內(nèi)，18:00～19:00時間段內(nèi)PDOP值較大，鎖星情況不好。此外，在聯(lián)合平差時,若將6個時段數(shù)據(jù)同時平差，可能會將觀測時段不好的數(shù)據(jù)引入，造成第一和第四時段相對精度較低。由于每個時段解算的數(shù)據(jù)都是單時段，所以多余觀測較少，精度較弱，但是整體來看，TEQC合并后的數(shù)據(jù)延長了同步觀測時間，優(yōu)化了觀測數(shù)據(jù)質(zhì)量，解算的點位坐標精度能夠滿足工程測量的要求。

3 結(jié)論

基于TEQC軟件強大的質(zhì)量檢核和數(shù)據(jù)編輯功能，對GNSS接收機在現(xiàn)場測量中由于各種復(fù)雜原因出現(xiàn)的短暫數(shù)據(jù)丟失情況進行后處理。對比分析表明，經(jīng)TEQC處理后的數(shù)據(jù)能夠延長同步觀測時間，滿足工程測量基線解算要求，解算的基線和平差后的點位精度符合鐵路工程測量的要求，不僅可以避免不必要的返工測量，提高工程測量的效率，節(jié)約成本，而且對于優(yōu)化基線解算和網(wǎng)平差也有著重要的意義。