周飛

(江蘇蘇地仁合土地房地產(chǎn)資產(chǎn)評估測繪造價咨詢有限公司，江蘇 南京 210005)

1 引 言

截至2017年底，拉薩市區(qū)的建筑規(guī)模近 100 km2，為了覆蓋此區(qū)域，決定布設(shè)15個D級控制點及75個E級控制點，經(jīng)過法定程序，我們從西藏自治區(qū)測繪局獲取了拉薩附近的C級控制點資料，由于有些點位已經(jīng)被破壞，經(jīng)過仔細(xì)普查篩選，決定選取C001自治區(qū)園林局、C003德慶鎮(zhèn)政府、C090熱堆、C189殯儀館這4個點作為起算點。

2 控制網(wǎng)布設(shè)觀測及數(shù)據(jù)處理

2.1 控制網(wǎng)布設(shè)

本測區(qū)將布設(shè)的15個D級GPS點與4個C級點按同步圖形擴展式布網(wǎng)，同步環(huán)之間以邊連接方式進(jìn)行網(wǎng)的構(gòu)建，從而由這些同步圖形組成一個D級GPS觀測網(wǎng)。繼而將布設(shè)的75個E級GPS點與前期布設(shè)的15個D級GPS點按同步圖形擴展式布網(wǎng)，同步環(huán)之間以邊連接的方式進(jìn)行網(wǎng)的構(gòu)建，從而由這些同步圖形組成一個E級GPS觀測網(wǎng)。具體布設(shè)的網(wǎng)型如圖1所示。

圖1 拉薩市D、E級網(wǎng)布設(shè)示意圖

由于拉薩市屬于河谷地形，東西狹長，將控制點均勻分布各主要街道，點位距離合理，覆蓋全面。

2.2 控制網(wǎng)觀測

根據(jù)實際項目需要，使用了6臺高精度大地型雙頻GPS接收機，天線采用扼流圈天線，性能穩(wěn)定可靠。對投入使用的儀器設(shè)備均經(jīng)過專業(yè)計量檢定機構(gòu)的檢定。D、E級GPS采用同步環(huán)邊連接的方式，進(jìn)行靜態(tài)定位模式測量。在安置儀器時，嚴(yán)格對中、整平，用三腳架安置天線時，其對中誤差不應(yīng)大于 1 mm。精確量取天線高，使在天線互為120°方向量取的天線高互差不超過 3 mm，按要求填寫測量手簿。

GPS原始觀測數(shù)據(jù)為Topcon的tps格式，為方便基線處理和后續(xù)利用，將原始觀測數(shù)據(jù)轉(zhuǎn)換為標(biāo)準(zhǔn)的Rinex格式，并將點號、天線高、天線類型、觀測者等相關(guān)信息輸入到Rinex文件中，使外業(yè)觀測信息完整。由于原始觀測數(shù)據(jù)在轉(zhuǎn)換成Rinex格式時，須輸入4位點號，于是將點號的數(shù)字之前加“0”補位。

外業(yè)觀測均按觀測計劃要求完成，經(jīng)檢查，外業(yè)數(shù)據(jù)記錄均按照設(shè)計要求的技術(shù)指標(biāo)執(zhí)行，數(shù)據(jù)質(zhì)量滿足要求，外業(yè)觀測手簿記錄完整。

2.3 數(shù)據(jù)處理

先對D級網(wǎng)進(jìn)行數(shù)據(jù)處理，基線處理采用LGO軟件進(jìn)行，網(wǎng)平差采用COSA軟件進(jìn)行。GPS觀測經(jīng)基線解算、檢核、基線向量組網(wǎng)后，進(jìn)行GPS最小約束(無約束)平差，以檢驗空間基線向量網(wǎng)的內(nèi)符合精度，檢驗基線是否存在粗差和系統(tǒng)誤差。

經(jīng)計算檢核，D級網(wǎng)共有62個同步環(huán)、160個異步環(huán)與11條重復(fù)基線。E級網(wǎng)處理方法與D級網(wǎng)類似，E級網(wǎng)經(jīng)計算檢核，整網(wǎng)共有303個同步環(huán)、218個異步環(huán)與26條重復(fù)基線。由于市場中各儀器的基線長度精度標(biāo)稱精度存在虛高，我們?nèi)匀徊捎肎B/T 18314-2001規(guī)范中的相鄰點間基線長度精度σ計算公式，公式如下：

(1)

其中a為固定誤差，單位mm；b為比例誤差系數(shù)；d為相鄰點距離，單位mm。對于D級網(wǎng)a取10，b取10；對于E級網(wǎng)，a取10，b取20。

基線解算中，主要是對于同步環(huán)閉合差、異步環(huán)閉合差和重復(fù)基線較差等技術(shù)指標(biāo)的處理。通過調(diào)整衛(wèi)星截止角，剔除觀測時間過短及周跳過多的衛(wèi)星數(shù)據(jù)，剔除受多路徑影響較大的觀測值，判別觀測值類型，觀察RATIO值限值等方法，使得各技術(shù)指標(biāo)在限差范圍內(nèi)。在處理數(shù)據(jù)過程中，我們明顯發(fā)現(xiàn)空間直角坐標(biāo)系中的Y值偏差明顯偏大，在低緯度地區(qū)意味著高程精度較差。為此我們在靜態(tài)觀測時延長了觀測時間，確保足夠的觀測數(shù)據(jù)。

滿足以上限差后，運用所得的合格基線數(shù)據(jù)進(jìn)行無約束平差，進(jìn)而加入4個C級控制點的坐標(biāo)成果予以固定進(jìn)行約束平差。平差結(jié)果的質(zhì)量評定我們主要以對單位權(quán)方差的檢查，對基線改正數(shù)的檢驗，對已知坐標(biāo)的檢驗。技術(shù)指標(biāo)均滿足規(guī)范要求后，出具各種報告。

3 中央子午線及投影面選擇計算

3.1 中央子午線選擇

拉薩市大地坐標(biāo)大致為東經(jīng)91°06′，北緯29°36′。最近的兩條標(biāo)準(zhǔn)中央子午線分別為東經(jīng)90°和東經(jīng)93°。那么拉薩城區(qū)某一點距離兩條中央子午線分別約為 80 km和 160 km，需考慮橢球面上的長度投影至高斯平面上產(chǎn)生的長度變形，公式如下：

(2)

其中，ym是兩端點間橫坐標(biāo)平均值，Rm是測區(qū)中點的平均曲率半徑，s0是歸算到橢球面上的長度。可以近似認(rèn)為Rm取值為地球平均半徑 6 371 km。由于測區(qū)東西方向大約為 20 km，若中央子午線穿過測區(qū)中心(即拉薩城市中心)，我們設(shè)定為91°08′，這里ym的最大取值為 10 km，那么根據(jù)式(2)，每千米最大長度變形僅為 1.23 mm，基本可以忽略不計。

3.2 投影面選擇

由于實地測量的邊長長度投影至橢球面會產(chǎn)生變形，可以根據(jù)最小二乘方法得到誤差方程式，按級數(shù)展開，取主項(其余可以忽略不計)，得公式如下：

(3)

其中，Hm是測區(qū)所在高程面對于橢球面的平均高程，Rc長度所在方向的橢球曲率半徑。這里Rc近似取值為地球平均半徑 6 371 km。由于拉薩海拔較高，如果還是按常規(guī)的大地高為0的投影面，那么根據(jù)式(3)，每千米的長度變形達(dá)到 56.5 cm，遠(yuǎn)遠(yuǎn)超過工程測量規(guī)范允許的1/40000，所以必須抬高投影面。選取大地高 3 600 m作為投影面，測區(qū)內(nèi)大地高最大為 3 700 m左右，每千米最大投影變形為 1.5 cm，符合規(guī)范要求。

綜合以上考慮，我們最終選擇了中央子午線91°08′，抵償高程投影面為大地高 3 600 m。

4 與西安80坐標(biāo)系轉(zhuǎn)換參數(shù)

為了方便CGCS2000坐標(biāo)系統(tǒng)與原西安80坐標(biāo)系統(tǒng)建立關(guān)聯(lián)，也為了方便當(dāng)?shù)夭粍赢a(chǎn)數(shù)據(jù)庫的銜接，故選取了遍布拉薩市區(qū)周邊與市中心區(qū)域的11個D級點作為轉(zhuǎn)換參數(shù)的起算點位，計算所得的七參數(shù)，由于其數(shù)據(jù)保密性質(zhì)，這里不做描述。其中，基礎(chǔ)點位轉(zhuǎn)換殘差最大值為 0.015 m，符合規(guī)范。所得西安80坐標(biāo)，方便施工單位進(jìn)行圖根點布設(shè)及接下來的測量需要。

5 結(jié) 論

拉薩市區(qū)地處西部高原，電離層活躍，觀測時間需要比對應(yīng)規(guī)范所規(guī)定等級網(wǎng)的時間適當(dāng)加長。由于地處城區(qū)，選取點位時應(yīng)避開高大建筑物，用雙拼扼流圈天線，避免多路徑效應(yīng)。在進(jìn)行GPS網(wǎng)解算時，三維坐標(biāo)Y值方向閉合差容易偏大，說明高原區(qū)域測取大地高需要加長觀測時間。在進(jìn)行嚴(yán)格的二維約束平差后，需用全站儀實地測量距離，與原坐標(biāo)進(jìn)行比對。

經(jīng)過嚴(yán)密的基線解算、GPS網(wǎng)平差，后又采用全站儀實地測距檢測，所有技術(shù)指標(biāo)符合規(guī)范要求。這套基于CGCS2000橢球的獨立坐標(biāo)系統(tǒng)將為拉薩市區(qū)的城市基礎(chǔ)建設(shè)，地籍管理測量需要等工作作出貢獻(xiàn)。