□ 王建營(yíng)張琦（.天津市測(cè)繪院，天津西青30038；.天津市武清區(qū)規(guī)劃建筑設(shè)計(jì)所，天津武清30700）

GPS在濱海新區(qū)地面沉降監(jiān)測(cè)中的應(yīng)用研究

□ 王建營(yíng)1張琦2
（1.天津市測(cè)繪院，天津西青300381；2.天津市武清區(qū)規(guī)劃建筑設(shè)計(jì)所，天津武清301700）

本文簡(jiǎn)要介紹了濱海新區(qū)沉降監(jiān)測(cè)中GPS數(shù)據(jù)處理方案，研究了高精度GPS數(shù)據(jù)處理參數(shù)設(shè)置、數(shù)據(jù)處理中涉及到的成果精度評(píng)定。結(jié)合GPS基準(zhǔn)站和一級(jí)GPS監(jiān)測(cè)點(diǎn)數(shù)據(jù)進(jìn)行區(qū)域地表形變監(jiān)測(cè)數(shù)據(jù)處理及數(shù)據(jù)分析，通過(guò)數(shù)據(jù)分析，簡(jiǎn)單論述了濱海新區(qū)地面沉降分布發(fā)育情況以及引起濱海新區(qū)地面沉降的影響因素，找出數(shù)據(jù)之間存在的規(guī)律性，為今后濱海新區(qū)地面沉降防治與研究提供數(shù)據(jù)支持。

GPS高程；濱海新區(qū)；地表沉降；形變監(jiān)測(cè)

0.引言

地面沉降是一種導(dǎo)致區(qū)域性地面高程降低的環(huán)境地質(zhì)災(zāi)害，常發(fā)生在人口密集、工業(yè)發(fā)達(dá)的城市，具有沉降速度相對(duì)較慢難以明顯察覺、危害面積廣、破壞嚴(yán)重、幾乎不可恢復(fù)性等特點(diǎn)，越來(lái)越引起人們的重視[1]。伴隨著各種測(cè)量?jī)x器、相關(guān)地面沉降監(jiān)測(cè)技術(shù)方法的不斷改進(jìn)，地面沉降的自動(dòng)化監(jiān)測(cè)工作取得了巨大的進(jìn)步，近年來(lái)興起的GPS（Global Position System）技術(shù)應(yīng)用于地面沉降監(jiān)測(cè)已經(jīng)成為一個(gè)熱點(diǎn)研究課題[2]。

隨著全球?qū)Ш叫l(wèi)星系統(tǒng) GNSS（Global Navigation Satellite System）技術(shù)的發(fā)展，很多國(guó)家和地區(qū)建立了GPS形變監(jiān)測(cè)網(wǎng)，用于研究區(qū)域性地面沉降[3]。利用GPS建立大區(qū)域的形變監(jiān)測(cè)網(wǎng)可有效地克服常規(guī)監(jiān)測(cè)網(wǎng)的弊端。它以地心坐標(biāo)系為參考系，很好地解決了各種形變監(jiān)測(cè)的基準(zhǔn)難題，而且GPS監(jiān)測(cè)成果不存在誤差累積[4]。

本文從理論上探討正高、正常高與大地高以及地面沉降量的關(guān)系，論述了由GPS和精密水準(zhǔn)測(cè)量所獲得的地面沉降量在理論上的一致性。

并將GPS測(cè)量技術(shù)應(yīng)用于濱海新區(qū)地面沉降監(jiān)測(cè)的實(shí)踐中，通過(guò)GPS測(cè)量結(jié)果與一等水準(zhǔn)測(cè)量結(jié)果的對(duì)比分析，從而論證了采用GPS測(cè)量方法可以用于大范圍的地面沉降監(jiān)測(cè)。

本文簡(jiǎn)要介紹了濱海新區(qū)GPS形變監(jiān)測(cè)網(wǎng)的布設(shè)和GPS數(shù)據(jù)處理方案，研究了高精度GPS數(shù)據(jù)處理參數(shù)設(shè)置、數(shù)據(jù)處理中涉及到的成果精度評(píng)定。結(jié)合濱海新區(qū)范圍內(nèi)的GPS基準(zhǔn)站和一級(jí)GPS監(jiān)測(cè)點(diǎn)數(shù)據(jù)進(jìn)行區(qū)域地表形變監(jiān)測(cè)數(shù)據(jù)處理及數(shù)據(jù)分析，通過(guò)數(shù)據(jù)分析，簡(jiǎn)單論述了全市地面沉降分布發(fā)育情況以及引起濱海新區(qū)地面沉降的影響因素，找出數(shù)據(jù)之間存在的規(guī)律性，為今后濱海新區(qū)地面沉降防治與研究提供數(shù)據(jù)支持。

1.數(shù)據(jù)采集

濱海新區(qū)地面沉降GPS監(jiān)測(cè)網(wǎng)絡(luò)由GPS基準(zhǔn)站網(wǎng)和一級(jí)GPS監(jiān)測(cè)網(wǎng)兩部分組成。濱海新區(qū)及周邊共有3座GPS基準(zhǔn)站和10個(gè)一級(jí)GPS監(jiān)測(cè)點(diǎn)，基本形成了GPS監(jiān)測(cè)的骨干網(wǎng)絡(luò)體系。3個(gè)GPS連續(xù)運(yùn)行基準(zhǔn)站10個(gè)一級(jí)GPS監(jiān)測(cè)點(diǎn)位置分布如圖1所示，圖中，CH02、KC01以及TJA1為連續(xù)運(yùn)行基準(zhǔn)站，其余為一級(jí)GPS監(jiān)測(cè)點(diǎn)。

圖1 濱海新區(qū)GPS監(jiān)測(cè)網(wǎng)1絡(luò)點(diǎn)位分布圖

3個(gè)GPS連續(xù)運(yùn)行參考站屬于市GPS連續(xù)運(yùn)行參考站網(wǎng)，實(shí)行全天候長(zhǎng)期監(jiān)測(cè)，數(shù)據(jù)自動(dòng)采集。10個(gè)一級(jí)GPS監(jiān)測(cè)點(diǎn)每年11月進(jìn)行外業(yè)測(cè)量1次。監(jiān)測(cè)網(wǎng)布設(shè)觀測(cè)采用參照《全球定位系統(tǒng)（GPS）測(cè)量規(guī)范》（GB/T 18314-2009）B級(jí)GPS網(wǎng)相關(guān)技術(shù)要求，基于GPS連續(xù)運(yùn)行基準(zhǔn)站網(wǎng)的點(diǎn)觀測(cè)模式。每個(gè)監(jiān)測(cè)點(diǎn)連續(xù)監(jiān)測(cè)3個(gè)時(shí)段，共計(jì)72小時(shí)，采樣率設(shè)定30秒[5]。各批次涉及點(diǎn)同步觀測(cè)，各站點(diǎn)平均距離約為50km。

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

數(shù)據(jù)處理包括監(jiān)測(cè)站基線解算以及平差處理兩部分。

2.1 基線處理

原始數(shù)據(jù)采集并整理完成以后，選擇適當(dāng)數(shù)量的IGS站數(shù)據(jù)，并使用精密星歷解算GPS數(shù)據(jù)，可求得高精度基線向量和測(cè)站坐標(biāo)。GPS基線解算是基于 Linux fedora 16平臺(tái)上運(yùn)行GAMIT軟件完成的。

表1 控制參數(shù)設(shè)置

2.2 平差處理

對(duì)GPS基線處理結(jié)果的平差方法是利用GLOBK軟件對(duì)GAMIT軟件計(jì)算的單天結(jié)果進(jìn)行平差處理，使用數(shù)據(jù)是GAMIT軟件解算產(chǎn)生的H文件，處理結(jié)果給出最終的GPS測(cè)站坐標(biāo)和精度評(píng)定及速率和精度評(píng)定。

此次GPS觀測(cè)結(jié)果的數(shù)據(jù)處理中，平差計(jì)算使用的軟件是GLOBK軟件5.18版本。計(jì)算方案為。

（1）對(duì)本期數(shù)據(jù)的GAMIT計(jì)算結(jié)果分別進(jìn)行單期的平差合并處理，計(jì)算出的GPS測(cè)站坐標(biāo)的最佳估計(jì)值和它的精度估計(jì)值。

（2）坐標(biāo)參考框架為ITRF2005全球坐標(biāo)系統(tǒng)框架，使用ITRF2005系統(tǒng)的核心觀測(cè)站的坐標(biāo)和坐標(biāo)運(yùn)動(dòng)速率為坐標(biāo)系統(tǒng)平差控制參數(shù)。對(duì)核心站的坐標(biāo)給予強(qiáng)約束，東西向和南北向?yàn)?mm，豎向?yàn)?0mm，對(duì)坐標(biāo)運(yùn)動(dòng)速率給予東西向和南北向約束為2mm/a，豎向?yàn)?mm/a[6]。本次數(shù)據(jù)處理中ITRF2005系統(tǒng)的核心觀測(cè)站有6個(gè)，它們分別是：SHAO、BJFS、USUD、IRKT、YSSK、POL2。對(duì)其余的GPS觀測(cè)站的坐標(biāo)和坐標(biāo)運(yùn)動(dòng)速率給予松約束，坐標(biāo)東西向和南北向約束為10m，豎向?yàn)?0m，速度約束為東西向和南北向?yàn)?0mm/a，豎向?yàn)?00mm/a[7]。

3.沉降提取

水準(zhǔn)測(cè)量測(cè)得的是正常高的變化，而GPS測(cè)得的是大地高的變化，同一個(gè)點(diǎn)上的正常高與大地高是不同的，兩者不能直接比較。同一個(gè)點(diǎn)上正常高的變化Δh與大地高的變化ΔU兩者存在如下關(guān)系。

其中cosa是該點(diǎn)上垂線偏差的余弦。由于天津地區(qū)a很?。ㄒ话阍趲资氲牧考?jí)），因此可以直接用ΔU代替Δh表示地面沉降的變化，即直接將GPS測(cè)得的大地高的變化代替水準(zhǔn)測(cè)量正常高的變化。于是，可用周期監(jiān)測(cè)的大地高變化量來(lái)表示當(dāng)年該點(diǎn)的沉降。

2013年沉降量值由2013年與2012年GPS站點(diǎn)大地高相減所得。根據(jù)2012～2013年GPS解算成果年度沉降值表，做出2012～2013年沉降速率等值線圖，如圖2所示。

圖2 2012-2013年濱海新區(qū)地面沉降速率等值線圖

從圖2可以得知，2012～2013年下沉量最大的GPS點(diǎn)為TJA1，沉降值達(dá)到了-85.4mm，下沉值最小的為TJ10，沉降值為-6.2mm，濱海新區(qū)所有GPS點(diǎn)沉降量平均值為-35.83mm。

4.成果比較

由上節(jié)可知，兩期GPS測(cè)量所獲得大地高變化量ΔU表示地面沉降量，與應(yīng)用兩期精密水準(zhǔn)測(cè)量方法所獲得的正常高變化量Δh表示地

面沉降量，在變化趨勢(shì)上基本一致。為了進(jìn)一步比較2012～2013年GPS測(cè)量成果與水準(zhǔn)測(cè)量成果的一致性，特從相關(guān)部門收集到2012年與2013年部分GPS監(jiān)測(cè)點(diǎn)水準(zhǔn)測(cè)量成果，通過(guò)對(duì)比分析，二者較差見表2。

表2 GPS測(cè)量與水準(zhǔn)測(cè)量成果較差表

從表2可知，GPS測(cè)量成果與水準(zhǔn)測(cè)量成果二者較差小于±5mm的站點(diǎn)有5個(gè)，占比為50%；較差介于5-10mm之間的站點(diǎn)有3個(gè)，占比30%；較差介于10-20mm之間的站點(diǎn)有2個(gè)，占比20%;較差大于20mm的站點(diǎn)有2個(gè)。統(tǒng)計(jì)以上數(shù)據(jù)可知，80%的站點(diǎn)較差小于10mm，充分說(shuō)明數(shù)據(jù)解算的高精度和可靠性。兩種測(cè)量方法成果較差分段統(tǒng)計(jì)圖如圖3所示。

圖3 GPS測(cè)量成果與水準(zhǔn)測(cè)量成果較差占比圖

有研究文章指出，通過(guò)十年GPS沉降監(jiān)測(cè)數(shù)據(jù)試驗(yàn)，全天津市GPS測(cè)量大地高變化與水準(zhǔn)測(cè)量大沽高變化之差的均值約為10mm。因此，表2里80%的對(duì)比結(jié)果正常。對(duì)比結(jié)果大于15mm只有TJA1。仔細(xì)分析TJA1站GPS測(cè)量與水準(zhǔn)測(cè)量之差差異大的原因，有以下三條：一是GPS觀測(cè)數(shù)據(jù)質(zhì)量不好，GPS信號(hào)產(chǎn)生多路徑效應(yīng)，影響了數(shù)據(jù)觀測(cè)質(zhì)量；二是地質(zhì)環(huán)境改變，TJA1周圍2013年5月有大型建筑工地開工，沉降加??；三是水準(zhǔn)觀測(cè)與GPS觀測(cè)不完全同步。由此造成水準(zhǔn)觀測(cè)與GPS測(cè)量之間差異。

通過(guò)不同監(jiān)測(cè)方法比較，GPS監(jiān)測(cè)精度滿足要求，可以應(yīng)用于濱海新區(qū)地面沉降監(jiān)測(cè)，一是GPS監(jiān)測(cè)能實(shí)現(xiàn)自動(dòng)化，大大節(jié)省了人力物力財(cái)力；二是它以地心坐標(biāo)系為參考系，很好地解決了形變監(jiān)測(cè)的基準(zhǔn)難題；三是GPS監(jiān)測(cè)不像水準(zhǔn)測(cè)量存在誤差累積的現(xiàn)象，使得監(jiān)測(cè)成果更接近于沉降真值。

5.結(jié)論

本文通過(guò)同期水準(zhǔn)監(jiān)測(cè)數(shù)據(jù)和GPS監(jiān)測(cè)數(shù)據(jù)對(duì)比，經(jīng)過(guò)統(tǒng)計(jì)分析，驗(yàn)證了GPS監(jiān)測(cè)技術(shù)應(yīng)用于地面沉降監(jiān)測(cè)的可行性。

從2012～2013年濱海新區(qū)沉降等值線圖，得出濱海新區(qū)沉降概況。基本情況是以TJA1為中心，沉降逐漸擴(kuò)散，中間沉降大，周邊沉降逐步減小。

通過(guò)GPS連續(xù)運(yùn)行基準(zhǔn)站觀測(cè)結(jié)果可知，基準(zhǔn)站所在地區(qū)地面沉降呈線性下沉，無(wú)法恢復(fù)。

[1]許才軍，張朝玉等．地殼形變測(cè)量與數(shù)據(jù)處理[M]．武漢：武漢大學(xué)出版社，2009（10）．

[2]許曙光，劉樹臣等．國(guó)內(nèi)外地面沉降研究動(dòng)向[A]．第六屆地面沉降論文集，2001．

[3]崔振東，唐益群．國(guó)內(nèi)外地面沉降現(xiàn)狀與研究[J]．西北地質(zhì)學(xué)報(bào)，2009，29（03）：275-278．

[4]GB/T 18314-2009．全球定位系統(tǒng)（GPS）測(cè)量規(guī)范[S]．北京：中國(guó)標(biāo)準(zhǔn)出版社，2009．

[5]董克剛，易長(zhǎng)榮．利用GPS監(jiān)測(cè)天津市地面沉降的可行性研究[J]．大地測(cè)量與地球動(dòng)力學(xué)，2008，28（04）：68-71．

[6]黃立人，武艷強(qiáng)等．GPS連續(xù)站在地面沉降監(jiān)測(cè)中的應(yīng)用[J].工程地質(zhì)學(xué)報(bào)，2005，S1(13)：94-97．

[7]任立生，黃立人等．天津GPS沉降監(jiān)測(cè)的十年試驗(yàn)結(jié)果[J]．測(cè)繪科學(xué)，2006，31(04)：17-19．

P228.3

B

2095-7319（2015）03-0062-05

王建營(yíng)（1981—），男，雙學(xué)士，工程師，天津市測(cè)繪院，現(xiàn)在主要從事GNSS應(yīng)用研究。