任政兆 朱照榮 張錫越 曾艷艷 張鳳錄 曹畢錚

(1. 北京市測(cè)繪設(shè)計(jì)研究院, 北京 100038;2. 城市空間信息工程北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100038)

0 引言

高程基準(zhǔn)是推算國家統(tǒng)一高程控制網(wǎng)中所有水準(zhǔn)高程的起算依據(jù),它包括一個(gè)水準(zhǔn)基準(zhǔn)面和一個(gè)永久性水準(zhǔn)原點(diǎn)。陳俊勇[1]指出為了保持我國高程控制網(wǎng)的現(xiàn)勢(shì)性,應(yīng)當(dāng)在15～20年復(fù)測(cè)一次。高程基準(zhǔn)雖經(jīng)過十多年的建設(shè)發(fā)展,但其基礎(chǔ)設(shè)施損壞嚴(yán)重,測(cè)量成果的現(xiàn)勢(shì)性難以滿足服務(wù)需求,屬性結(jié)果維持難度較大[2]。因此在新型基礎(chǔ)測(cè)繪背景下,利用連續(xù)運(yùn)行參考站(continuously operating reference stations,CORS)站數(shù)據(jù)解算的高差成果與傳統(tǒng)水準(zhǔn)測(cè)量的高差成果進(jìn)行分析將會(huì)有意義。

有研究表明地殼運(yùn)動(dòng)、地下流體資源開采、城市建設(shè)等人類活動(dòng)會(huì)引起城市局部地面沉降[3]。北京市存在多個(gè)沉降中心且不同程度的地面沉降造成了區(qū)域高程參考框架難以維持,同時(shí)在此框架下的測(cè)繪成果難以反映現(xiàn)勢(shì)情況。楊強(qiáng)等人[4]指出連續(xù)多年的全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system,GNSS)觀測(cè)數(shù)據(jù)在U方向精度有所提高。袁林果等人[5]提出利用GNSS連續(xù)運(yùn)行參考站進(jìn)行地殼形變監(jiān)測(cè)具有重要意義。張?jiān)娪竦热薣6]基于球諧函數(shù)定量分析了GNSS的大氣負(fù)荷效應(yīng),得到了大氣負(fù)荷是GNSS垂向分量季節(jié)性變化的因素之一的結(jié)論。有研究學(xué)者[7-9]分析了環(huán)境負(fù)載對(duì)GNSS測(cè)站序列的影響。即使GORS站解算坐標(biāo)速度更快、更簡(jiǎn)單便捷,但全球定位系統(tǒng)(global positioning system,GPS)觀測(cè)數(shù)據(jù)在高程方向受到的噪聲干擾遠(yuǎn)大于水平方向[10],所以在精密水準(zhǔn)測(cè)量、工程控制網(wǎng)中都要求使用傳統(tǒng)水準(zhǔn)測(cè)量傳遞高程信息,為解決因區(qū)域不均勻沉降造成的水準(zhǔn)點(diǎn)高程變化較大的問題,北京市測(cè)繪設(shè)計(jì)研究院每年以玉淵潭水準(zhǔn)原點(diǎn)作為基準(zhǔn)點(diǎn),進(jìn)行北京市沉降區(qū)高程復(fù)測(cè)及原點(diǎn)網(wǎng)監(jiān)測(cè)工程。

1 數(shù)據(jù)收集

收集2012年1月至2021年7月的北京市CORS站和中國周邊國際GNSS服務(wù)(international GNSS service, IGS)基準(zhǔn)站觀測(cè)數(shù)據(jù),包括:5個(gè)北京市CORS站點(diǎn)，即昌平站(CHPN)、東三旗站(DSQI)、西集站(XIJI)、牛欄山站(NLSH)、朝陽站(CHAO)；12個(gè)IGS站點(diǎn)，即aira,badg,bjfs,pimo,daej,pol2,shao,tcms,hyde,tixi,urum,yssk,數(shù)據(jù)采樣間隔為30 s。由于DSQI、CHAO 2012—2015年數(shù)據(jù)質(zhì)量存在問題(嚴(yán)重的周跳),只對(duì)該兩站2016—2021年的數(shù)據(jù)進(jìn)行處理。

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

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

基線處理軟件使用美國麻省理工學(xué)院和Scripps研究所共同研制的流體力學(xué)模擬基礎(chǔ)軟件(GAMIT)軟件,該軟件是世界上最優(yōu)秀的GNSS數(shù)據(jù)處理軟件之一。利用GAMIT 10.70版本和GLOBK 5.12版本軟件,附加合理的約束對(duì)北京市CORS站網(wǎng)進(jìn)行數(shù)據(jù)處理。

2.2 參數(shù)設(shè)置

數(shù)據(jù)處理參數(shù)設(shè)置具體情況如表1所示。

表1 GAMIT參數(shù)設(shè)置

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

首先對(duì)北京市CORS系統(tǒng)數(shù)據(jù)進(jìn)行預(yù)處理,主要批量修改2018年之前的基準(zhǔn)站數(shù)據(jù)文件頭接收機(jī)和天線信息格式。以GAMIT/GLOBK(一個(gè)研究高精度全球定位系統(tǒng)測(cè)量數(shù)據(jù)處理等領(lǐng)域的綜合分析軟件包)軟件為主進(jìn)行GNSS數(shù)據(jù)處理。通過對(duì)CORS站和IGS站每天的GNSS觀測(cè)量數(shù)據(jù)(觀測(cè)站位置、接收機(jī)鐘差與衛(wèi)星鐘差等多種參數(shù))解算,利用GLOBK軟件獲得CORS站的單日區(qū)域松弛解(h-file),經(jīng)過坐標(biāo)約束后,得到基于國際地球參考框架(international terrestrial reference frame,ITRF)2014下所有基準(zhǔn)站位置的高精度坐標(biāo)和時(shí)間序列。

2.4 精度檢核

由于GAMIT軟件采用的是網(wǎng)解(即全組合解),同步環(huán)閉合差在基線解算的同時(shí)已經(jīng)進(jìn)行了分配。為了檢核GAMIT軟件基線解的同步環(huán),可以把時(shí)間序列的國際GNSS服務(wù)(international GNSS service,IGS)作為同步環(huán)質(zhì)量好壞的一個(gè)指標(biāo),一般要求NRMS值不能大于0.5且不等于0。2012—2021年基線解算NRMS值大小如下圖1所示。

圖1 基線解算NRMS值

通過圖1分析,除2016年外其他年份基線解算的NRMS值均小于0.5,滿足精度要求。經(jīng)查閱,2016年前221天DSQI站由于接收機(jī)和天線設(shè)備故障,數(shù)據(jù)觀測(cè)質(zhì)量較差,存在嚴(yán)重的周跳變化,導(dǎo)致基線解算的NRMS值存在超過限值的情況。從2016年第222天開始,DSQI站更換天線、接收機(jī)設(shè)備后基線解算的NRMS值均小于0.5且趨于穩(wěn)定狀態(tài)。

3 高程基準(zhǔn)變化分析

3.1 時(shí)間序列分析

CHAO站和DSQI站位于在北京市嚴(yán)重沉降檢測(cè)區(qū)內(nèi),圖2是基準(zhǔn)站的原始時(shí)間序列,在垂直方向明顯發(fā)現(xiàn),CHAO站在2015—2018年之間是不斷沉降的,2019—2021年沉降速度有了明顯減緩。CHPN站相對(duì)比較穩(wěn)定,年周期變化相對(duì)明顯。NLSH站2012—2017年下沉較大,但是2018年開始沉降逐漸減緩,甚至2020—2021年該站有抬升的趨勢(shì)。DSQI站在2016年時(shí)間序列出現(xiàn)階躍,相對(duì)位置存在一定偏移,經(jīng)核實(shí)是該站更換了天線,該站同樣存在一定的沉降。XIJI站在2014—2021年階段中處于均勻沉降。

圖2 基準(zhǔn)站原始時(shí)間序列

3.2 基站速度分析

表2給出了基準(zhǔn)站VN,VE,VU的三個(gè)方向近十年的年均運(yùn)動(dòng)速率,基準(zhǔn)站在南北方向每年10 mm左右,東西方向每年的運(yùn)動(dòng)量大約在30 mm,垂直方向是根據(jù)不同地形運(yùn)動(dòng)量各不相同,但各站的運(yùn)動(dòng)速度均為負(fù)值,表示基準(zhǔn)站均呈現(xiàn)沉降趨勢(shì),其中在嚴(yán)重沉降區(qū)內(nèi)的DSQI站沉降速率最大,每年可以達(dá)到50.54 mm,CHAO站明顯低于DSQI站,每年以約22 mm的速度下沉,NLSH次之,年沉降速率約為17 mm,CHPN站屬于北京市沉降區(qū)外的山區(qū),該點(diǎn)的高程比較穩(wěn)定,大約為每年1.79 mm的速度下降。

表2 基準(zhǔn)站年均運(yùn)動(dòng)速率 單位：mm/a

3.3 結(jié)合水準(zhǔn)復(fù)測(cè)分析

為了驗(yàn)證CORS站解算高程基準(zhǔn)變化的準(zhǔn)確性,收集了近6年北京市沉降區(qū)高程水準(zhǔn)復(fù)測(cè)的高差結(jié)果,二者對(duì)比如圖3。由于2013—2015年北京市沒有對(duì)CORS基準(zhǔn)站進(jìn)行聯(lián)測(cè),選取2016年的基準(zhǔn)站高程作為基準(zhǔn)做差值。

雖然部分測(cè)站的GNSS結(jié)果與水準(zhǔn)成果具有一定的差距,比如CHAO站2020年高差二者相差4.8 cm,DSQI站兩個(gè)高差結(jié)果平均差11.3 cm,但不可否認(rèn)的是GNSS高差與水準(zhǔn)高差的趨勢(shì)項(xiàng)相似度較高,其中NLSH站和XIJI站二者相差毫米級(jí),除CHPN站比較穩(wěn)定,其他四個(gè)基準(zhǔn)站均在2016年下沉比較嚴(yán)重,2019—2021年期間基準(zhǔn)站沉降趨勢(shì)明顯減弱,其中CHAO站在2021年水準(zhǔn)測(cè)量成果中抬升4 mm,NLSH站從2018年開始一直以8 mm/a的速度做抬升運(yùn)動(dòng),XIJI站處于北京市沉降區(qū)內(nèi),保持下沉趨勢(shì),但運(yùn)動(dòng)速度明顯減緩。圖3也從側(cè)面反映出北京市“海綿城市”綠地建設(shè)的有效成果。堅(jiān)持生態(tài)為本、自然循環(huán),增強(qiáng)城市綠地海綿體功能,采取嚴(yán)格限制地下水開采、南水北調(diào)的工程等措施從一定程度上緩解了北京市區(qū)域沉降的趨勢(shì),取得了良好的整治成果。

圖3 GNSS高差與水準(zhǔn)高差對(duì)比

4 結(jié)束語

本文對(duì)北京CORS基準(zhǔn)站2012—2021年的數(shù)據(jù)進(jìn)行了處理,獲取了5個(gè)基準(zhǔn)站高程基準(zhǔn)變化的時(shí)間序列,結(jié)合了近6年來北京市水準(zhǔn)測(cè)量情況,得出以下結(jié)論：

(1)CORS站垂直方向運(yùn)動(dòng)普遍存在明顯的周期變化。

(2)將CORS站數(shù)據(jù)解算的高差結(jié)果與傳統(tǒng)水準(zhǔn)測(cè)量的結(jié)果相比較,雖然部分測(cè)站的GNSS結(jié)果與水準(zhǔn)成果具有一定的差距,但GNSS高差與水準(zhǔn)高差的趨勢(shì)項(xiàng)相似度較高,驗(yàn)證了利用CORS站獲取垂直方向的準(zhǔn)確性。

(3)地面沉降對(duì)北京城市建設(shè)和基礎(chǔ)設(shè)施造成了隱患,北京市采取嚴(yán)格限制地下水開采、北京市“海綿城市”綠地建設(shè)、南水北調(diào)的工程等措施從一定程度上緩解了北京市沉降的趨勢(shì),取得了良好的整治成果,以CHAO站、NLSH站尤為明顯,甚至出現(xiàn)了抬升的趨勢(shì)。