劉 凱,彭方輝

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

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

TBC和COSA_GPS GNSS數(shù)據(jù)處理軟件在延安引黃工程中的應(yīng)用

劉 凱,彭方輝

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

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

1 概 述

1.1 工程任務(wù)

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

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

1.2 數(shù)據(jù)處理軟件

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

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

2 GPS基線解算

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

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

表1 GNSS靜態(tài)測量的基本技術(shù)要求表

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

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

圖1 TBC基線處理結(jié)果圖

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

3 GPS網(wǎng)平差

為了充分利用TBC和COSA_GPS各自的優(yōu)點(diǎn)，將同時運(yùn)用這2個軟件于同一組數(shù)據(jù)的處理過程中，從而獲取更為可靠、精確的GPS數(shù)據(jù)處理結(jié)果。

3.1 TBC軟件中GPS網(wǎng)平差

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

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

3.2 COSA_GPS中GPS網(wǎng)平差

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

(1) 新建工程

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

圖2 COSA_GPS新建工程圖

(2) 讀取同步基線數(shù)據(jù)并形成獨(dú)立基線文件

COSA_GPS以基線向量解算軟件按時段導(dǎo)出的基線向量文件作為觀測值輸入數(shù)據(jù)，由此形成的獨(dú)立基線向量作為平差時的觀測值，以輸入的已知點(diǎn)坐標(biāo)作為起算基準(zhǔn)，完成控制網(wǎng)平差計算。

(3) 按照規(guī)范檢查基線、閉合環(huán)是否超限

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

(1)

(2)

(3)

(4)

(5)

式中:n為閉合環(huán)中的邊數(shù)；W為環(huán)線全長閉合差，mm。

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

(4) 三維無約束平差

在WGS-84空間直角坐標(biāo)系中進(jìn)行三維向量網(wǎng)平差，從而驗(yàn)證基線向量解算結(jié)果的內(nèi)符合精度。在平差之前首先輸入一個三維已知點(diǎn)，平差之后得到各點(diǎn)的GPS大地高和二維基線向量，其中，GPS大地高可用于后續(xù)的高程擬合，二維基線向量作為二維約束平差的觀測值。

(5) 二維約束平差

設(shè)置中更改橢球類型，輸入和TBC中相同的國家C級網(wǎng)高斯坐標(biāo)成果進(jìn)行二維約束平差。結(jié)果主要的精度指標(biāo)為：

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

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

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

3.3 TBC和COSA_GPS平差結(jié)果分析

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

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

4 獨(dú)立坐標(biāo)系的建立及高程擬合

4.1 建立獨(dú)立坐標(biāo)系

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

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

圖3 COSA_GPS中獨(dú)立坐標(biāo)系的建立圖

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

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

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

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

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

5 結(jié) 語

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

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

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

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

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

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

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

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

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

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

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

[8] 劉明波，楊躍青，鹿恩鋒.南方GPS軟件中的獨(dú)立網(wǎng)平差方法探討[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”技術(shù)集成研究與開發(fā).

TP319

A

10.3969/j.issn.1006-2610.2015.03.026