TBC和COSA_GPS GNSS數據處理軟件在延安引黃工程中的應用

劉 凱,彭方輝

(陜西省水利電力勘測設計研究院測繪分院,西安 710002)

文章編號：1006—2610(2015)03—0092—05

TBC和COSA_GPS GNSS數據處理軟件在延安引黃工程中的應用

劉 凱,彭方輝

(陜西省水利電力勘測設計研究院測繪分院,西安 710002)

筆者將TBC和COSA_GPS兩種數據處理軟件分別運用在工程前期的數據處理中，并通過充分的比較論證，得出兩者相結合的處理方法，即TBC用于基線解算、處理，COSA_GPS用于后期平差。討論了COSA_GPS在精度要求較高的隧洞工程中，軟件的一點一方位方法建立獨立坐標系的優越性。論述了這2種軟件的擬合GPS高程成果在水準檢核和小范圍測圖中的應用。 關鍵詞：TBC;COSA_GPS;獨立坐標系;高程擬合

1 概 述

1.1 工程任務

引黃工程是“十二五”期間延安市規劃建設的重大水源工程，該項目對于加快全市產業轉型和能源化工基地建設，促進經濟社會可持續發展，具有重大歷史意義。此引水工程始于黃河上游延川縣延水關鎮王家渠村，止于延安經濟技術開發區。取水量近期預計達到7 196萬m3/a，遠期有望將達到1億m3/a以上，為新建Ⅲ等中型供水工程，其由五大部分組成，分別為：取水樞紐工程(2處)、泵站工程(10座，其中加壓泵站8座，取水泵站2座)、輸水管(洞)線工程181.25 km、水庫工程(2處，調蓄水庫1處，事故應急水庫1處)和凈化水廠工程(2處)。

本次工程線路經由地段以山地為主，地形困難、起伏較大且極其破碎，土質都為濕陷性黃土，溝道較多，通視條件較差；加上水工建筑數量繁多，距離跨度較長，測量誤差大，若采用常規方法建網，不僅測量精度過低，工作量大，而且測量周期太長而影響施工進度。所以平面控制測量采用了高精度GPS網作為首級施工控制并加密的方法。具體分為骨架網和加密網，其中骨架網點按照每10 km布設1組對點；加密網點149個，以500 m為間距布設，并兩兩相互通視。骨架網除聯測已有國家C級點外，為了檢驗原有勘測設計階段GPS網成果點的穩定性，將原有部分成果點(基本都已破壞)納入骨架網。所有網點采用國家基準進行約束，對于隧洞工程建立了獨立坐標系，高程控制測量采用精密三等水準觀測，高差成果由高程擬合結果進行檢核。

1.2 數據處理軟件

Trimble Business Center(TBC)是Trimble的新一代后處理軟件，不僅能夠處理GNSS(包含GPS 和GLONASS)數據，還可以處理全站儀、水準儀、3D 掃描儀數據，集成了功能強大的可視工具和建模工具，利用多種視圖全面反映數據信息，全新的處理算法保證其處理速度，并提供了靈活的處理配置方案[1]。

COSA_GPS即“GPS工程測量控制網通用平差軟件包”，能讀取各種GPS接收機的基線向量解算文件，進行網的三維無約束平差和二維約束平差、在工程獨立坐標系中固定一點一方向的平差、高程擬合等功能[2]。

2 GPS基線解算

本次工程控制網由骨架網和加密網構成，骨架網全線共布設GPS點22個，其中現場澆注16個(C1…C16)，按照每10 km布設成對點；國家C級GPS已知點4個(F090、F091、F098、F012), 經查勘

樁志穩定，可以作為本次測量起算點；勘測設計階段布設的已知點1個(BD1，D級)，楊家山隧洞已知觀測墩1個(GP1，D級)，經查勘樁志穩定，可以作為本次測量檢核點。靜態觀測作業采用4臺Trimble雙頻(標稱精度為5×10-6mm)GPS接收機，以邊連式聯測各個樁志，全網共觀測12個時段，各項技術指標均按《水利水電工程測量規范》[3]要求執行，主要技術要求見表1。

表1 GNSS靜態測量的基本技術要求表

采用TBC數據處理軟件進行基線解算工作，在TBC中新建工程后導入Trimble R8的T01/T02文件后，按接收機和時段輸入對應的點號、接收機型號、天線類型以及量取方式和天線高等參數，確定無誤后進行整體基線處理，處理過程中隨時檢核解的類型是否固定，并以均方根RRMS值判斷觀測值的精度(RRMS越小代表觀測值的精度越高)。

對于部分平面或者垂直精度太差的基線選擇單基線處理，可通過在對應基線的時段編輯器中禁用不健康的衛星、屏蔽質量較差的衛星觀測值、調整觀測時間等方法處理之后保存結果完成基線處理，其結果如圖1所示，隨后在工具欄中查看基線處理以及閉合差報告。

圖1 TBC基線處理結果圖

通過報告可見該網整體精度較高，完全滿足規范要求。此外，在TBC基線解算過程中發現，基線向量的解算成果與測站環境因素、觀測時間、時段及基線長度都有著密切的關系，如9:00左右測區出現星況不良導致返測、基線過短相對誤差過大而超限等情況，通過調整觀測時間以及延長觀測時間的方法使得上述情況明顯改善。

3 GPS網平差

為了充分利用TBC和COSA_GPS各自的優點，將同時運用這2個軟件于同一組數據的處理過程中，從而獲取更為可靠、精確的GPS數據處理結果。

3.1 TBC軟件中GPS網平差

在所有基線都合格以后就可以三維無約束平差，在菜單欄中依次進行“測量”-“平差網”即可執行自由網平差操作，平差結束后在平差網對話框的“結果”選項卡中查看本次平差結果。如果顯示“卡方檢驗(95%)”失敗則通過將上次平差的參考因子加權后再次平差，若檢驗通過則自由網平差完成。

在約束網平差之前需要建立對應的坐標系統，本次工程采用國家80坐標系，中央子午線為111°，在菜單欄中選擇“工具”-“坐標系統管理器”建立相應的坐標系統后即可將此坐標系加載到本次工程。在已知的4個公共點(BD1與GP1用于檢核)中添加當地坐標后即可進行二維約束平差，方法和上述一致,平差完后查看平差報告是否滿足限差要求，并將此結果與后續COSA_GPS成果進行比較。

3.2 COSA_GPS中GPS網平差

在COSA_GPS中處理由TBC導出的基線數據,骨架網共觀測了12個時段，從TBC中共導出12個后綴為*.asc格式的基線向量觀測文件，以此作為COSA_GPS觀測值輸入文件。軟件解算過程如下。

(1) 新建工程

新建工程界面見圖2，這里需要注意的是獨立基線的條數，省缺值為-1，即認為選定的基線全部為獨立基線，則平差后的精度指標比實際值偏高，在此輸入獨立基線的實際條數，軟件將對平差后的精度指標進行修正，從而與獨立基線平差結果的精度指標基本一致[2]。其他文本框按照坐標系統實際情況輸入即可。

圖2 COSA_GPS新建工程圖

(2) 讀取同步基線數據并形成獨立基線文件

COSA_GPS以基線向量解算軟件按時段導出的基線向量文件作為觀測值輸入數據，由此形成的獨立基線向量作為平差時的觀測值，以輸入的已知點坐標作為起算基準，完成控制網平差計算。

(3) 按照規范檢查基線、閉合環是否超限

在工具中分別生成同步環閉合差、異步環閉合差和重復基線差文件，檢查文件是否滿足限差規定[3]，其中，異步環限差如下所示。

(1)

(2)

(3)

(4)

(5)

式中:n為閉合環中的邊數；W為環線全長閉合差，mm。

同步環各坐標分量閉合差的限差值為異步環閉合差限差的1/2。從生成的異步環文件可見，異步環閉合差的坐標分量平均值保持在1 cm左右，對于某些長度較短的基線，盡管相對閉合差略大但仍符合要求。

(4) 三維無約束平差

在WGS-84空間直角坐標系中進行三維向量網平差，從而驗證基線向量解算結果的內符合精度。在平差之前首先輸入一個三維已知點，平差之后得到各點的GPS大地高和二維基線向量，其中，GPS大地高可用于后續的高程擬合，二維基線向量作為二維約束平差的觀測值。

(5) 二維約束平差

設置中更改橢球類型，輸入和TBC中相同的國家C級網高斯坐標成果進行二維約束平差。結果主要的精度指標為：

1) 邊長精度。邊長中誤差0.9～3.5 mm，平均值2.4 mm，相對中誤差3.26×10-6～3.12×10-7，平均相對誤差0.40×10-6，最弱邊邊長中誤差2 mm。

2) 點位精度。點位中誤差0～5.4 mm，平均值4.3 mm，其中X坐標中誤差的平均值為3.2 mm，Y坐標中誤差的平均值為2.9 mm，最弱點的點位中誤差僅為5.4 mm。

此平差精度較好，完全可以滿足施工控制的需求。

3.3 TBC和COSA_GPS平差結果分析

由基于相同基準的二維約束平差得出，TBC和COSA_GPS兩種軟件的處理結果較為一致，無顯著差異，在縱坐標(x)平均相差0.7 mm，最大坐標差為1.5 mm，在橫坐標(y)平均相差0.9 mm，最大坐標差為2.9 mm，此差異小于水利水電施工測量規范要求的點位誤差(如特高精度水電站施工控制網規定最弱點精度為5 mm，平均點位精度3 mm)，可見這2個軟件的二維約束平差結果沒有本質的區別且外部符合精度較高。

與TBC的平差結果文件相比，COSA_GPS的輸出成果內容全，且同步環閉合差、異步環閉合差和重復基線差結果文件能夠兼容Excel的CSV格式，方便用戶進行閉合差和重復基線差的統計分析[4,5]；加之平差結果文件中包括了基線向量殘差、可靠性，平差后方位角、邊長及精度，最弱邊以及相鄰點坐標分量及其相對中誤差等結果，可以作為最終的成果資料提供給甲方，故在本次應用中采用TBC處理基線并導出合格的基線數據，利用COSA_GPS進行平差得到最終成果并提供資料。

4 獨立坐標系的建立及高程擬合

4.1 建立獨立坐標系

引黃工程沿線有姚家山、新舍古、盧草梁、文安驛、薛家溝和曹家園則等隧洞，隧洞總長度達25.017 km，平均長度都在1 km以上，根據《水利水電工程測量規范》[3]要求，對于水利樞紐地區及重要工程建筑物地區，當測區內投影長度變形值不大于5 cm/km，均應采用國家坐標系統，按統一的高斯投影3°分帶；對于投影長度變形值大于5 cm/km時，可采用一個國家大地點的坐標和該點至另一大地點的方位角作為起始數據的獨立坐標系統。采用國家80國家坐標系時，姚家山隧洞投影變形為5.84 cm/km，新舍古隧洞投影變形為5.9 cm/km。如果按照此坐標成果進行貫通，則會造成橫向貫通誤差偏大。故工程采用2套坐標系統：針對隧洞工程必須使用獨立坐標系統，其他施工項目使用國家坐標系統。

這里主要采用COSA_GPS中的一點一方位方法建立獨立坐標系，該過程首先將三維GPS基線向量投影至工程橢球(在設置中將橢球類型改為80)，此時工程橢球球心和扁率與80參考橢球相同，但長半軸要考慮引水工程的平均高程面的海拔高度以及大地水準面與80橢球的差距[6]；然后根據設置的中央子午線投影到高斯平面上再轉化至工程獨立坐標系。如圖3的對話框所示，固定點信息有點名、高斯平面坐標和正常高、大地坐標、投影面正常高。對于獨立網來說，已知點、已知方位最好是選在網的重心附近，這里已知點選擇處于測區中心C8點的國家高斯平面坐標，使得平差成果掛靠在國家坐標系下，固定方位角為該點至另一已知點的方向。投影至高斯平面的中央子午線取110°01′09″(固定點附近)，投影面正常高取900(測區平均高程面)。實踐證明，通過這種方法建立的獨立坐標系的各項綜合變形都遠遠小于上限，滿足了隧洞工程對基準的要求。

圖3 COSA_GPS中獨立坐標系的建立圖

4.2 TBC與COSA_GPS在高程擬合中的差異

本次工程沿線的二等水準路線共計146 km，且已知水準點間相距較遠，為保證觀測質量和數據的可靠性，防止外業成果出現錯誤，采取擬合GPS高程的方法對聯測對點之間的高差進行檢核；同時，為了滿足施工過程中線路整改和工程棄渣場等測圖的需求(比例尺：1∶1000)，對水準難以到達的地段考慮采用擬合高程成果作為點校正的起始數據。

TBC中采用EGM96大地水準面模型；在COSA_GPS中的擬合方法選擇平面擬合，擬合模型為：N=f(B,L)=f1+f2·B+f3·L。在此只試驗選擇整網中的一段(C5-C8、D33-D52)，得到最后結果對比如圖4所示。

從圖4可以看出，利用COSA_GPS和TBC軟件的不同模型在此次工程中整體效果差異較小，最大相差3.6 mm。2種結果的最大誤差都位于D36和D50處，其原因主要是這2個點的高程沒有和它附近的已知高程點保持同樣的變化趨勢，例如，D36的高程值遠大于它左右的已知高程點。同時，運用高程擬合所得結果和已知水準高程在10 km范圍內最大只有2.4 cm(COSA_GPS)的差異；而COSA_GPS和TBC擬合結果平均誤差分別為4.2 mm與4.8 mm，標準差分別為8.9 mm與8.1 mm，所以在小范圍內高程擬合精度可以滿足規范中測圖比例尺對高程的要求。

圖4 COSA_GPS和TBC高程擬合對比圖

5 結 語

在水利水電施工控制網的布設中，上述過程都是經常會碰到的，本文所討論的軟件使用方法以及得出的結論具有較高的探索和實際參考價值，主要包括以下幾個方面。

(1) COSA_GPS和TBC在國家坐標系下的二維約束平差結果沒有本質上的區別，故成果都可作為最終資料交付施工單位使用，但實際運用中，由于COSA_GPS的成果的可讀性強，便于理解，故在工程運用中我們采取以TBC軟件用于基線解算和處理，COSA_GPS用于平差計算。

(2) 在隧洞、橋梁及大壩等對相對精度較高的控制網，應該使控制點成果的投影長度變形值滿足一定上限，同時，固定工程網中的某一特定方向的方位角，例如隧洞控制網的隧洞軸線方向、大壩控制網的壩軸線方向等。為此，可使用一點一方位的方法建立獨立坐標系。

(3) 本工程地勢復雜，測區內多為山區地帶，水準測量進度緩慢而已知水準點之間距離較遠，實驗中COSA_GPS和TBC不同的高程擬合模型所得結果差異不大，故兩者的高程擬合成果都可對每天水準測段高差成果進行檢核，同時驗證了擬合成果也滿足了范圍內的測圖需求。

[1] 吳偉.基于TBC的高精度GPS數據處理若干問題探討[D].桂林：桂林理工大學測繪地理信息學院,2012.

[2] 郭際明,羅年學.GPS工程測量網通用平差軟件包使用說明書[R].武漢大學測繪學院,2010.

[3] SL 197-97,水利水電工程測量規范[S].北京：中國水利水電出版社,2013.

[4] 郭際明,羅年學.精密工程測量數據處理綜合系統講座.第一講.科傻(COSA)系統構成及其在工程測量中的應用[J].測繪信息與工程,2010,35(1):53-54.

[5] 郭際明,羅年學.精密工程測量數據處理綜合系統講座.第三講.COSA_GPS及在大規模高精度GPS控制網數據處理中的應用[J].測繪信息與工程,2010,35(3):19-20.

[6] 張項鐸,張正祿.淺談特長隧道GPS網的布測及其精度-以秦嶺隧道工程為例[J].勘察科學技術.2000(5).46-50.

[7] 李祖鋒,賀麗娟,張成增,周濤.高角度信號遮擋區域GPS測量控制網精度控制[J].西北水電,2014,(1): 19-23.

[8] 劉明波，楊躍青，鹿恩鋒.南方GPS軟件中的獨立網平差方法探討[J].西北水電， 2012，(3)：13-15.

Application of Data Processing Software of TBC and COSA_GPS GNSS in Yan'an Yellow River Diversion Project

LIU Kai, PENG Fang-hui

(Shaanxi Province Institute of Water Resources and Electric Power Investigation and Design,Survey & Mapping Branch, Xi'an 710002,China)

Both data Processing Software of TBC and COSA_GPS GNSS are applied in the data process of the early engineering stage. Furthermore, the processing method by combination of two software is derived from sufficient comparison and demonstration. Namely, TBC is applied for baseline calculation and processing, COSA_GPS for adjustment at late stage. In the paper, superiority of independent coordinate system established by one-point one-azimuth method of the software of COSA_GPS applied in tunneling with high precision requirement are studied. Application of the fitted GPS elevation by the two software in level verification and small-scope mapping are demonstrated.

TBC; COSA_GPS; independent coordinate system; elevation fitting

2014-09-19

劉凱(1988- ),男,陜西省榆林市人,助理工程師,從事“3S”技術集成研究與開發.

TP319

A

10.3969/j.issn.1006-2610.2015.03.026