曹炳強　劉智強　魯澤宇　簡程航　張雙成

關鍵詞：CORS系統(tǒng);地面沉降;大地高精度;監(jiān)測網(wǎng)型

近年來，國內(nèi)很多行業(yè)都建立了自己的衛(wèi)星定位連續(xù)運行網(wǎng)絡，建設完成、運行正常的GNSS基準站，可以滿足不同行業(yè)、不同用戶對定位精度的個性化要求，可以提供測繪、氣象、地震、交通、水利、應急保障和地質(zhì)災害監(jiān)測等多種服務（慕永燈，2019）。目前北京市有許多家企事業(yè)單位都進行了CORS系統(tǒng)的建設，實現(xiàn)了基準站網(wǎng)在北京市氣象局、北京市地震局、千尋網(wǎng)絡位置服務公司等不同行業(yè)間的數(shù)據(jù)共享，基準站在全市數(shù)量眾多且分布廣泛（北京市規(guī)劃和自然資源委員會，2019）。

地面沉降能夠造成交通和地下管線受損、河流防汛能力降低，同時也給城市建設的布局與規(guī)劃造成困難，已經(jīng)發(fā)展成為城市的主要地質(zhì)災害（楊艷等，2012;羅勇等，2013）。地面沉降成因復雜，在許多國家都有發(fā)生，其主要特點是持續(xù)周期長、影響范圍廣以及治理難度大。為了應對地面沉降問題，北京市很早就建立了以高精度水準測量、GNSS測量、InSAR測量和7個地面沉降監(jiān)測站構(gòu)成的地面沉降立體監(jiān)測體系，監(jiān)測范圍覆蓋主要沉降區(qū)（劉明坤等，2012;羅勇等，2014;北京市地質(zhì)礦產(chǎn)勘查局等，2008）。本文基于26個CORS站的監(jiān)測數(shù)據(jù)，探討了CORS系統(tǒng)在城市地面沉降監(jiān)測中的應用，也為下一步充分利用北京市已有的CORS站并將部分CORS站納入北京市地面沉降GNSS監(jiān)測系統(tǒng)，提供了參考數(shù)據(jù)和依據(jù)。

1北京市GNSS沉降監(jiān)測網(wǎng)現(xiàn)狀

目前北京市專門用于地面沉降監(jiān)測的GNSS監(jiān)測點共有105個，其中包括3個連續(xù)運行監(jiān)測站，分布在海淀區(qū)、通州區(qū)、昌平區(qū)、大興區(qū)、順義區(qū)、西城區(qū)、東城區(qū)、朝陽區(qū)、豐臺區(qū)、房山區(qū)等10個區(qū)，各點間距2～5 km不等，基本上覆蓋了北京市整個平原地區(qū)（北京市地質(zhì)礦產(chǎn)勘查開發(fā)局等，2008）。每個測量點均建為永久性混凝土觀測墩，觀測墩頂部有不銹鋼強制對中裝置，以消除天線對中誤差的影響。每年7—9月統(tǒng)一對100余個監(jiān)測點進行一次外業(yè)測量，接收機采用天寶雙頻接收機，接收機天線采用天寶扼流圈天線，以減弱多路徑效應。為了提高測量結(jié)果精度，削弱系統(tǒng)及周期性測量誤差，每期測量時每個監(jiān)測點對應的測量設備和測量時間都是固定的。根據(jù)多期的GNSS監(jiān)測成果顯示，GNSS監(jiān)測成果與一等精密水準測量、InSAR等其他監(jiān)測結(jié)果具有良好的符合性。近幾年隨著北京城市快速發(fā)展，導致少量監(jiān)測點遭到破壞，一部分監(jiān)測點觀測環(huán)境受到影響，降低了部分站點的監(jiān)測精度，也使得整個GNSS監(jiān)測網(wǎng)的精度受到影響。

2 CORS數(shù)據(jù)解算

2.1軟件介紹

GNSS基線解算采用GAMIT/GLOBK軟件包，其由美國麻省理工學院（MIT）與海洋研究所（SIO）研制。GAMIT軟件是開源軟件，高校及科研機構(gòu)可以通過官網(wǎng)申請，得到許可之后能夠免費使用，其自動化程度高，通過選取適當數(shù)量的IGS站，并使用精密星歷解算GNSS數(shù)據(jù)，可求得高精度基線向量和地面測站坐標，是目前國際、國內(nèi)科研及生產(chǎn)單位應用廣泛且解算精度較高的GNSS數(shù)據(jù)處理軟件之一（閔陽，2018;劉洋洋等，2018;呂旭陽，2015），最新版本已經(jīng)能夠解算包括我國北斗數(shù)據(jù)在內(nèi)的多系統(tǒng)數(shù)據(jù)。本文數(shù)據(jù)解算軟件為GAMIT/GLOBK 10.71，操作系統(tǒng)為Ubuntu 12.04。

2.2數(shù)據(jù)來源

本文搜集了26個北京市連續(xù)運行監(jiān)測站的數(shù)據(jù)，包含3個專門的北京市地面沉降監(jiān)測連續(xù)運行站，23個測繪及氣象部門的連續(xù)運行站（圖1）。由于監(jiān)測站在高程方向上的最大周期性變化能達到4 cm（王利，2015），為了削弱監(jiān)測站在高程方向上的周期性誤差，所采用的數(shù)據(jù)采集日期均為2019年8月1日—8月7日的數(shù)據(jù)。

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

基線解算是將3個基巖基準點（BJFS、BJSH、JIXN）與26個連續(xù)運行站的數(shù)據(jù)聯(lián)網(wǎng)計算。在解算時對基巖點和連續(xù)運行站分別進行約束：基巖點為0.05 m的水平方向約束，0.1 m的垂直方向約束;連續(xù)運行站水平方向和垂直方向的約束分別為0.1 m和0.2 m。此外，對于海洋潮汐、地球固體潮、極潮、大氣負載潮等引起的站臺漂移，以及太陽光壓、相對論效應引起的測量誤差都進行了適當?shù)奶幚恚ǜ吆?014;黨亞民等，2007;李征航等，2016;劉大杰等，1996）。具體的基線解算參數(shù)見表1。

3解算結(jié)果分析

對26個連續(xù)運行站2020年8月1日—8月7日的數(shù)據(jù)進行了基線解算與平差處理，分別得到了監(jiān)測點的空間直角坐標和大地坐標，然后對其精度分別進行了統(tǒng)計分析。一般情況下，單時段解算得到的標準化均方根殘差（NRMS）值越小，基線解算精度越高;其值越大，則基線解算精度越低，一般要求NRMS值小于0.3 mm（曹炳強等，2015;王樹東等，2018）。最后對解算結(jié)果文件（Q文件）中的NRMS值進行了統(tǒng)計。NRMS值統(tǒng)計結(jié)果：最大值為0.1813，最小值為0.1612，平均值為0.1686，解算得到的NRMS值為0.16～0.19 mm，滿足基線解算要求。圖2為CORS站三維坐標方向的精度，圖3為大地坐標系下CORS站在BLH方向上的精度。

從圖2可以看到，CORS站解算結(jié)果在三維方向上的坐標精度都在2mm之內(nèi)，其中X、Y、Z方向的誤差平均值分別為0.47 mm、0.81 mm、0.71 mm，解算結(jié)果在三維方向上的坐標精度都在1mm之內(nèi)的CORS站占84.6%。

北京市地面沉降GNSS監(jiān)測站點沉降變化信息主要通過比較各期監(jiān)測站點的大地高H得到，所以GNSS監(jiān)測站點大地高H精度的高低直接影響到北京市地面沉降GPS監(jiān)測的可靠性。圖3可以看到大地坐標系下，CORS站解算結(jié)果大地高方向的誤差都在2 mm之內(nèi)，其中有23個站的高程誤差小于1 mm，占總數(shù)的88.5%。具體統(tǒng)計結(jié)果見表2。

綜合以上分析得出，采用CORS站數(shù)據(jù)解算，可以高精度獲取包括大地高在內(nèi)的坐標信息，可以高精度獲取地面沉降監(jiān)測成果。

4結(jié)論與建議

通過對26個北京市連續(xù)運行監(jiān)測站數(shù)據(jù)進行處理，可以得到基于CORS站解算的結(jié)果大地高方向誤差都在2 mm之內(nèi)，而且有23個站的高程誤差小于1 mm，能夠得到高精度的解算成果。

目前各城市都建有大量的CORS站，分布廣泛，聯(lián)合多部門實現(xiàn)CORS系統(tǒng)在不同行業(yè)間的數(shù)據(jù)共享，將不同行業(yè)的CORS站納入城市地面沉降監(jiān)測系統(tǒng)中，這樣可以借助CORS系統(tǒng)增加地面沉降監(jiān)測面積和提高監(jiān)測點的覆蓋密度，增強現(xiàn)有GNSS監(jiān)測網(wǎng)的精度與可靠性，提高解算成果的精度，也更容易實現(xiàn)城市地面沉降監(jiān)測的自動化。