張文峰，張旭東

（1.寧波市測(cè)繪設(shè)計(jì)研究院，浙江 寧波 315041）

寧波市似大地水準(zhǔn)面數(shù)字高程基準(zhǔn)建設(shè)及精度檢驗(yàn)

張文峰1，張旭東1

（1.寧波市測(cè)繪設(shè)計(jì)研究院，浙江 寧波 315041）

基于2000國(guó)家大地坐標(biāo)系和1985國(guó)家高程基準(zhǔn)，闡述了在NBCORS和寧波市基本高程控制網(wǎng)的三維框架下組織實(shí)施寧波市似大地水準(zhǔn)面數(shù)字高程基準(zhǔn)建設(shè)的方法，并進(jìn)行了似大地水準(zhǔn)面的正確性和實(shí)用性精度檢驗(yàn)以及精度區(qū)域劃分。

似大地水準(zhǔn)面；數(shù)字高程基準(zhǔn)；正確性檢驗(yàn)；實(shí)用性檢驗(yàn)

2009~2011年，寧波市測(cè)繪設(shè)計(jì)研究院組織實(shí)施了覆蓋全寧波市域及近海區(qū)域的“寧波市似大地水準(zhǔn)面數(shù)字高程基準(zhǔn)建設(shè)”項(xiàng)目。項(xiàng)目旨在探索“GNSS+似大地水準(zhǔn)面”測(cè)量新技術(shù)的應(yīng)用，以實(shí)現(xiàn)傳統(tǒng)高程基準(zhǔn)建立和維持模式的根本性轉(zhuǎn)變。

1 數(shù)字高程基準(zhǔn)建設(shè)

1.1 工藝流程

本項(xiàng)目中寧波市似大地水準(zhǔn)面數(shù)字高程基準(zhǔn)的建設(shè)工藝流程主要分為4個(gè)步驟：①高程異?？刂凭W(wǎng)建設(shè)，即一體化同步建立基于GNSS的基本控制網(wǎng)和全市基本高程控制網(wǎng)；②似大地水準(zhǔn)面計(jì)算，包括GNSS、水準(zhǔn)數(shù)據(jù)處理以及似大地水準(zhǔn)面模型的建立；③似大地水準(zhǔn)面精度檢驗(yàn)，包括采用正確性GNSS靜態(tài)檢驗(yàn)、實(shí)用性的網(wǎng)絡(luò)RTK檢驗(yàn)以及精度區(qū)域的劃分；④似大地水準(zhǔn)面的應(yīng)用，編制在線三維坐標(biāo)轉(zhuǎn)換系統(tǒng)，為用戶提供準(zhǔn)實(shí)時(shí)的GNSS高程轉(zhuǎn)換數(shù)字化服務(wù)。

從項(xiàng)目建設(shè)的工藝流程分析，本次精化主要解決4個(gè)問題：①統(tǒng)一性，即統(tǒng)一采用2000國(guó)家大地坐標(biāo)系和1985國(guó)家高程基準(zhǔn)作為平面和高程的起算基準(zhǔn)；②同步性，不僅保持水準(zhǔn)和GNSS測(cè)量的同步性，同時(shí)與全省CORS的觀測(cè)保持同步，以保障參考站點(diǎn)的CGCS2000大地坐標(biāo)的時(shí)效性；③取代傳統(tǒng)高程基準(zhǔn)，提出了使用似大地水準(zhǔn)面模型取代傳統(tǒng)水準(zhǔn)標(biāo)石并維持現(xiàn)代高程基準(zhǔn)的方法，研究和編制了采用GNSS高程測(cè)量代替各等級(jí)水準(zhǔn)測(cè)量的技術(shù)指標(biāo)，并經(jīng)過了理論推導(dǎo)和生產(chǎn)實(shí)踐應(yīng)用；④實(shí)用性，即項(xiàng)目利用似大地水準(zhǔn)面成果，結(jié)合網(wǎng)絡(luò)RTK技術(shù)，研制在線三維坐標(biāo)轉(zhuǎn)換系統(tǒng)，滿足城市建設(shè)定位需求。

1.2 GNSS基本控制網(wǎng)測(cè)量

高程異?？刂凭W(wǎng)的GNSS網(wǎng)建立以ZJCORS為基礎(chǔ)，將NBCORS納入全省的CORS框架網(wǎng)中，與全省同步聯(lián)測(cè)，精確解算并確定NBCORS各參考站的CGCS2000坐標(biāo)。在此基礎(chǔ)上，再組織實(shí)施全市GNSS基本控制網(wǎng)。該網(wǎng)的特點(diǎn)如下：

①GNSS網(wǎng)點(diǎn)布設(shè)均勻，結(jié)構(gòu)合理。將GNSS點(diǎn)的布設(shè)與水準(zhǔn)路線選取統(tǒng)一考慮，減少了水準(zhǔn)聯(lián)測(cè)的工作量，使控制點(diǎn)上同時(shí)具有平面和高程成果。②采用8 h×2時(shí)段的C級(jí)GNSS同步環(huán)觀測(cè)技術(shù)，全網(wǎng)平均重復(fù)設(shè)站達(dá)4.32。③綜合了多種精密的數(shù)據(jù)處理方法，其GNSS網(wǎng)坐標(biāo)分量在水平方向上的平均精度為±0.003 5 m，大地高方向的平均精度為±0.006 7 m。

GNSS網(wǎng)的三維約束平差結(jié)果建立了該地區(qū)CGCS2000坐標(biāo)系下的平面測(cè)繪基準(zhǔn)，為高精度高分辨率似大地水準(zhǔn)面的建立提供了均勻、可靠和精確的地心坐標(biāo)數(shù)據(jù)。

1.3 基本高程控制網(wǎng)測(cè)量

高程異常控制網(wǎng)中的GNSS/水準(zhǔn)點(diǎn)的高程聯(lián)測(cè)也布設(shè)成水準(zhǔn)網(wǎng)，該網(wǎng)基于原有寧波市基本高程控制網(wǎng)，對(duì)其進(jìn)行改造，以國(guó)家一等水準(zhǔn)點(diǎn)為起算數(shù)據(jù)，全部布設(shè)成二等水準(zhǔn)網(wǎng)，旨在聯(lián)測(cè)并獲取GNSS/水準(zhǔn)點(diǎn)的高精度正常高。該網(wǎng)的特點(diǎn)如下：

① 水準(zhǔn)網(wǎng)優(yōu)化調(diào)整，全網(wǎng)全部改造并構(gòu)成28個(gè)閉合環(huán)的二等水準(zhǔn)網(wǎng)；②對(duì)于全部國(guó)家一等水準(zhǔn)點(diǎn)南北連貫采用一等水準(zhǔn)進(jìn)行檢測(cè)，確定了符合性起算點(diǎn)；③水準(zhǔn)觀測(cè)采用精密水準(zhǔn)測(cè)量技術(shù)，共完成了230 km一等水準(zhǔn)檢測(cè)、2 044 km二等水準(zhǔn)測(cè)量；④水準(zhǔn)數(shù)據(jù)處理時(shí)考慮了溫度改正，進(jìn)行了水準(zhǔn)標(biāo)尺長(zhǎng)度、正常水準(zhǔn)面不平行、重力異常以及環(huán)線閉合差改正。

全網(wǎng)平差后，二等水準(zhǔn)各種誤差均滿足限差要求，其每km水準(zhǔn)測(cè)量的全中誤差為±0.877 mm，符合設(shè)計(jì)和規(guī)范要求，為高精度高分辨率似大地水準(zhǔn)面的建立提供了精確的水準(zhǔn)高程數(shù)據(jù)。

1.4 數(shù)字高程基準(zhǔn)計(jì)算

似大地水準(zhǔn)面數(shù)字高程計(jì)算融合了規(guī)劃區(qū)已有建設(shè)成果，共使用了3 673個(gè)點(diǎn)重力數(shù)據(jù)和174個(gè)GNSS/水準(zhǔn)點(diǎn)資料，以EIGEN03C地球重力場(chǎng)模型作為參考重力場(chǎng)，由第二類Helmert凝集法完成了似大地水準(zhǔn)面的計(jì)算。其特點(diǎn)如下：

1）在利用原有規(guī)劃區(qū)精化成果時(shí)，首先將規(guī)劃區(qū)與全市高程異?？刂凭W(wǎng)平面和高程采用統(tǒng)一起算點(diǎn)約束，計(jì)算其CGCS2000大地坐標(biāo)和水準(zhǔn)高程，經(jīng)過重復(fù)點(diǎn)的分析和選擇，再進(jìn)一步進(jìn)行高程異常計(jì)算。

2）格網(wǎng)重力異常的內(nèi)插和推估計(jì)算，在航天飛機(jī)雷達(dá)地形測(cè)繪使命（SRTM）高分辨率7.5"×7.5"數(shù)字地面模型基礎(chǔ)上，利用地形均衡歸算，通過移去-還原原理計(jì)算。其中地形均衡歸算采用Airy-Haiskanen均衡模型，地形改正和均衡改正采用了考慮地球曲率的嚴(yán)密球面積分公式，積分半徑為300 km。

3）似大地水準(zhǔn)面的計(jì)算采用了第二類Helmert凝集法，公式是考慮了地球曲率影響的嚴(yán)密球面積分公式。

通過計(jì)算，174個(gè)GNSS/水準(zhǔn)資料與重力似大地水準(zhǔn)面獨(dú)立比較精度為±0.018 m。利用球冠協(xié)調(diào)和分析方法將GNSS水準(zhǔn)與重力似大地水準(zhǔn)面聯(lián)合求解得出的2'×2'格網(wǎng)似大地水準(zhǔn)面精度達(dá)到±0.010 m。

2 似大地水準(zhǔn)面精度檢驗(yàn)

在似大地水準(zhǔn)面精化之后，為了檢測(cè)似大地水準(zhǔn)面成果的正確性和實(shí)用性，寧波市測(cè)繪設(shè)計(jì)研究院先后組織了不同等級(jí)、不同區(qū)域的GNSS精度測(cè)試，對(duì)寧波市高精度高分辨率似大地水準(zhǔn)面成果進(jìn)行了質(zhì)量檢驗(yàn)。

2.1 正確性檢驗(yàn)

正確性檢測(cè)采用外部檢核方式，總共選取了26個(gè)符合要求的點(diǎn)位作為檢測(cè)點(diǎn)，主要均勻分布在寧波市的慈溪、余姚、寧海、象山、奉化等5個(gè)代表性的特征區(qū)域。每個(gè)區(qū)域均布設(shè)8個(gè)檢測(cè)點(diǎn)，其中包含2～3個(gè)已知點(diǎn)。檢測(cè)點(diǎn)布設(shè)均考慮了利用具有二等水準(zhǔn)成果的點(diǎn)位或布設(shè)于已知二等水準(zhǔn)點(diǎn)附近，既符合GNSS觀測(cè)要求，又便于水準(zhǔn)聯(lián)測(cè)。每個(gè)檢測(cè)點(diǎn)的GNSS、水準(zhǔn)測(cè)量等級(jí)與似大地水準(zhǔn)面精化時(shí)高程異常控制點(diǎn)的測(cè)量要求一致，數(shù)據(jù)處理也采用相同的軟件和處理方法，檢測(cè)結(jié)果見表1。

表1 正確性靜態(tài)檢測(cè)結(jié)果統(tǒng)計(jì)表/m

由表1可知，應(yīng)用似大地水準(zhǔn)面成果計(jì)算的GNSS高程與直接水準(zhǔn)高程較差中誤差為±0.015 m，達(dá)到了較高精度水平，滿足了精化預(yù)期的目標(biāo)。

2.2 實(shí)用性檢驗(yàn)

由于本項(xiàng)目實(shí)施的似大地水準(zhǔn)面精化范圍較大，實(shí)用性檢測(cè)主要基于NBCORS系統(tǒng)，在全市區(qū)按照8 km×8 km的密度，均勻布設(shè)網(wǎng)絡(luò)RTK高程檢測(cè)點(diǎn)，總共布設(shè)了132個(gè)網(wǎng)絡(luò)RTK檢測(cè)點(diǎn)。布設(shè)的實(shí)用性檢測(cè)點(diǎn)位均具有近期觀測(cè)的四等及其以上等級(jí)的已知水準(zhǔn)成果，這些點(diǎn)位穩(wěn)定可靠，沒有進(jìn)行重新測(cè)量。

RTK觀測(cè)執(zhí)行《衛(wèi)星定位城市測(cè)量技術(shù)規(guī)范》中的一級(jí)GNSS控制點(diǎn)觀測(cè)要求，觀測(cè)4個(gè)測(cè)回，測(cè)回間的垂直坐標(biāo)分量均不大于3 cm，取各測(cè)回結(jié)果的平均值作為最終觀測(cè)成果。最后再利用LGO軟件，在軟件中加載似大地水準(zhǔn)面精化模型，導(dǎo)入觀測(cè)數(shù)據(jù)，可直接獲取網(wǎng)絡(luò)RTK正常高，其檢測(cè)結(jié)果見表2。

表2 實(shí)用性網(wǎng)絡(luò)RTK動(dòng)態(tài)檢測(cè)結(jié)果統(tǒng)計(jì)表/cm

從表2可以看出，網(wǎng)絡(luò)RTK動(dòng)態(tài)檢測(cè)的模型計(jì)算高程與其水準(zhǔn)高程比較平均較差為0.011 m，中誤差為±0.022 m，均沒有超出相關(guān)規(guī)范規(guī)定的±6 cm限差，可以滿足現(xiàn)有國(guó)家規(guī)范中代替等外水準(zhǔn)、圖等水準(zhǔn)的技術(shù)指標(biāo)要求。

2.3 精度區(qū)域劃分

考慮到測(cè)繪生產(chǎn)活動(dòng)的行政區(qū)域性，為了指導(dǎo)寧波市各個(gè)行政區(qū)域內(nèi)利用似大地水準(zhǔn)面精化成果進(jìn)行高程測(cè)量，本項(xiàng)目將精化模型精度區(qū)域劃分到行政區(qū)，共6個(gè)區(qū)域，分布情況見圖1。

圖1 全市似大地水準(zhǔn)面精度區(qū)域劃分圖

1）內(nèi)符合精度區(qū)域劃分：該精化模型在寧波市各個(gè)行政區(qū)域內(nèi)的內(nèi)符合精度在±9.2 mm～±12.4 mm，山區(qū)精化模型內(nèi)符合精度為±12.0 mm，平地為±9.7 mm，全市范圍內(nèi)的內(nèi)符合精度為±10 mm。依據(jù)相關(guān)規(guī)范，在寧波市域范圍內(nèi)，精度均滿足利用四等GNSS高程測(cè)量代替四等水準(zhǔn)測(cè)量的技術(shù)指標(biāo)要求。

2）外符合精度區(qū)域劃分：全市實(shí)用性靜態(tài)GNSS高程精度為±18.0 mm，實(shí)用性網(wǎng)絡(luò)RTK高程精度為±22.0 mm，可滿足寧波城市建設(shè)各領(lǐng)域的高程測(cè)量應(yīng)用需求。

3 結(jié) 語

1）項(xiàng)目綜合運(yùn)用了各種現(xiàn)代大地測(cè)量技術(shù)，確定了覆蓋寧波全市域、2'×2'高分辨率、±1.0 cm高精度似大地水準(zhǔn)面模型；建立了統(tǒng)一基準(zhǔn)、統(tǒng)一地勢(shì)、統(tǒng)一精度、無縫接邊、幾何與物理一體化的寧波市似大地水準(zhǔn)面數(shù)字高程基準(zhǔn)；實(shí)現(xiàn)了傳統(tǒng)高程基準(zhǔn)建立和維持模式的根本性轉(zhuǎn)變，完善了寧波市現(xiàn)代測(cè)繪基準(zhǔn)體系的建設(shè)成果。

2）似大地水準(zhǔn)面同精度正確性檢測(cè)結(jié)果表明，寧波市似大地水準(zhǔn)面的外符合精度為±0.015 m，精化成果比較理想，達(dá)到了預(yù)期精化的目標(biāo)。實(shí)用性動(dòng)態(tài)檢測(cè)結(jié)果表明，網(wǎng)絡(luò)RTK高程測(cè)量中誤差為±0.022 m，可以滿足相應(yīng)技術(shù)規(guī)范規(guī)定代替等外水準(zhǔn)、圖等水準(zhǔn)的技術(shù)指標(biāo)要求，這為寧波市解決各類測(cè)繪工作的高程控制提供了基本技術(shù)保障，為下一步基于NBCORS的網(wǎng)絡(luò)RTK高程精度研究打下堅(jiān)實(shí)的基礎(chǔ)。

3）精度區(qū)域劃分表明，寧波市似大地水準(zhǔn)面高程異常模型在不同區(qū)域、不同地勢(shì)的精度基本無差異，但采用靜態(tài)和動(dòng)態(tài)測(cè)量GNSS高程精度與相應(yīng)的GNSS等級(jí)精度有關(guān)。

[1] 席青驥. 城市似大地水準(zhǔn)面精化成果的檢核及其應(yīng)用[J].地理空間信息，2007，5(4):90-92

[2] 肖建華，王厚之，李江衛(wèi)，等. 武漢市區(qū)域似大地水準(zhǔn)面精化成果應(yīng)用分析[J].地理空間信息，2008，6(2):1-3

[3] 章傳銀，黨亞民，柯寶貴，等. 高精度海岸帶重力似大地水準(zhǔn)面的若干問題討論[J].武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2012(10):1 150-1 154

[4] 符華年，張旭東，史秀保. 基于NBCORS的似大地水準(zhǔn)面精化成果檢測(cè)與分析[J].城市勘測(cè)，2012(2):93-95

[5] 曹學(xué)禮，張旭東，史秀保. 似大地水準(zhǔn)面精化成果的檢測(cè)與探討[J].測(cè)繪標(biāo)準(zhǔn)化，2008(4):22-24

[6] 陳達(dá)玉，譙勇，劉俊，等. 似大地水準(zhǔn)面精化方法研究與精度分析[J].鐵道勘察，2011(3):21-23

[7] 盧榮. 廣州市似大地水準(zhǔn)面精度檢核研究[J]. 礦山測(cè)量，2011(4):56-57

Construction and Testing Accuracy of Ningbo Quasi-geoid Digital Vertical Datum

byZHANG Wenfeng

Based on the 2000 national geodetic coordinate system and the 1985 National Vertical Datum, this paper described the method which organized the construction of Ningbo quasi-geoid digital vertical datum in NBCORS and Ningbo basic elevation control network 3D framework. At the same time, it conducted the accuracy test of correctness and practicality and the accuracy region's partition.

quasi-geoid, digital vertical datum, correctness testing,practicality testing

P223.0

B

1672-4623（2013）02-0157-03

10.11709/j.issn.1672-4623.2013.02.052

2012-12-05。

張文峰，工程師，主要從事大地測(cè)量技術(shù)的研究和應(yīng)用。