岳　順，岳建平，邱志偉，汪學(xué)琴

(1. 河海大學(xué)地球科學(xué)與工程學(xué)院，江蘇 南京 210098；2. 水能資源利用關(guān)鍵技術(shù)湖南省重點實驗室，湖南 長沙 410014)

?

基于D-InSAR技術(shù)的地面沉降監(jiān)測方法

岳順1,2，岳建平1，邱志偉1，汪學(xué)琴1

(1. 河海大學(xué)地球科學(xué)與工程學(xué)院，江蘇 南京 210098；2. 水能資源利用關(guān)鍵技術(shù)湖南省重點實驗室，湖南 長沙 410014)

針對傳統(tǒng)地面沉降監(jiān)測方法的不足，D-InSAR技術(shù)因其全天候、空間分辨率高、精度高等優(yōu)點，在地面沉降監(jiān)測方面其表現(xiàn)得越來越重要。本文對香港地區(qū)的3幅ERS影像進行處理，提取了該地區(qū)高精度的三維沉降場，結(jié)果表明，該地區(qū)最大沉降值為-0.551 6 cm，且沉降趨勢較緩慢，分布平均，D-InSAR技術(shù)可以有效地提供高精度的地面沉降信息，在地面沉降監(jiān)測方面具有較高的實用性。

D-InSAR；地面沉降；ERS；高精度；三維沉降場

傳統(tǒng)的沉降監(jiān)測技術(shù)有水準(zhǔn)測量、全站儀、GPS測量等，但其具有勞動強度大、監(jiān)測成本大、空間分辨率低等缺陷，無法滿足現(xiàn)代地面沉降監(jiān)測的要求[1]。隨著對衛(wèi)星測量技術(shù)的不斷研究和應(yīng)用[2]，利用遙感測量技術(shù)來監(jiān)測地表沉降的方法也越來越完善。作為雷達遙感的重要分支，InSAR在過去的十幾年里取得了重大的進展[3-5]。該技術(shù)為大面積快速獲取高分辨率、高精度的DEM信息提供了一條新的路徑[6-8]。D-InSAR技術(shù)是由InSAR發(fā)展而來，通過若干幅影像干涉形成差分組合，消除地形影響，進而得到地表形變信息[7-10]。它能夠沿著雷達的視線向監(jiān)測地球表面厘米級或更加微小的形變[11]，以揭示地球許多現(xiàn)象，如地震形變、地面沉降、火山運動、冰川漂移以及山體滑坡等[12-14]。D-InSAR技術(shù)具有全天候、全天時、覆蓋面廣、高度自治化、觀測周期短、成本小和精度高等優(yōu)勢，現(xiàn)已成為地面沉降監(jiān)測的主要技術(shù)。利用D-InSAR技術(shù)對地面沉降進行監(jiān)測，有助于評估和預(yù)測地面沉降和殘余沉降災(zāi)害，進而減少或避免沉降災(zāi)害帶給人們和環(huán)境的潛在危險[15]，D-InSAR技術(shù)為地面沉降的監(jiān)測提供了全新、快速和有效的方法。

一、 D-InSAR進行地面沉降監(jiān)測方法研究

圖1為D-InSAR成像幾何示意圖，A1、A2分別表示衛(wèi)星在地表未發(fā)生形變前獲取的兩幅SAR圖像[16]，則由P點返回信號為

(1)

(2)

圖1　D-InSAR原理示意圖

在地表P點發(fā)生形變后，衛(wèi)星在A3位置獲取第3幅SAR圖像，此時P點返回的信號為

(3)

式中，ΔRd為視線向形變量。第1、3兩幅SAR圖像所引起的干涉紋圖的相位φ1既包含了區(qū)域的地形信息又包含了觀測期間地表的形變信息，由圖1可以看出，基線B、B′分別在雷達視線向上的投影分別為

B∥=Bsin(θ-α1)

(4)

(5)

將第1、3兩幅影像進行配準(zhǔn)后干涉共軛相乘得

(6)

得到干涉相位

(7)

如果要獲取地表的形變信息，則必須消除區(qū)域的地表信息，可以通過第1、2兩幅影像干涉相乘形成干涉紋圖，其干涉相位只包含地形信息

(8)

根據(jù)式(7)、式(8)所得最終的形變相位為

(9)

二、 實例數(shù)據(jù)分析

香港由九龍半島、香港島、新界內(nèi)陸地區(qū)及262個大小島嶼組成。隨著香港地區(qū)的快速發(fā)展，地面沉降現(xiàn)象越來越嚴(yán)重。常規(guī)的水準(zhǔn)和GPS監(jiān)測手段存在監(jiān)測時間長、成本高等缺點，而D-InSAR技術(shù)正好可以彌補這些缺點，正在逐漸應(yīng)用到地面沉降監(jiān)測。本文基于D-InSAR技術(shù)對香港地區(qū)3幅影像進行處理，所用數(shù)據(jù)是ERS衛(wèi)星數(shù)據(jù)，分別表示為影像24440、4767、6270。影像24440、4767為地面沉降之前的兩幅影像，作為地形像對，影像6270為地面沉降之后的影像，與影像24440作為形變像對，具體的主輔影像參數(shù)見表1。

表1　主輔影像參數(shù)

分別對地形像對和形變像對進行處理，得到兩像對的相干系數(shù)圖和干涉相位圖。

從相干系數(shù)圖(如圖2、圖3所示)可以看出，亮度越強的地方相干系數(shù)越大，即相干性越好，另外可以看出地形像對的相干性比形變像對的相干性好。從圖4、圖5可以看出地形像對和形變像對的干涉圖發(fā)生了明顯的變化，說明在這期間地面發(fā)生了沉降。對得到的兩幅干涉圖進行相位解纏、差分、地理編碼，最終得到該地區(qū)的三維沉降圖(如圖6所示)。

圖2　地形像對相干系數(shù)圖

圖3　形變像對相干系數(shù)圖

圖4　地形像對干涉圖

圖5　形變像對干涉圖

從圖6中可以看出該地區(qū)的最大沉降量是-0.551 6cm，大部分在[-0.2,-0.4]區(qū)間內(nèi)，且沉降趨勢較緩慢，分布平均，說明該地區(qū)沉降量較小。

圖6　三維沉降圖

三、結(jié)束語

本文基于D-InSAR技術(shù)對香港地區(qū)3幅影像進行處理，得到該地區(qū)的三維沉降圖，實現(xiàn)了該地區(qū)地面沉降的動態(tài)監(jiān)測和分析，為地面災(zāi)害預(yù)測提供了新的監(jiān)測手段和方法。如果可以采集到沉降區(qū)域的地層條件、地下水等數(shù)據(jù)，可以利用D-InSAR對其地表移動參數(shù)進行反演，從地質(zhì)的角度解釋地面沉降的原理，完善D-InSAR技術(shù)監(jiān)測的結(jié)果分析。

[1]岳建平,方露.城市地面沉降監(jiān)控技術(shù)研究進展[J]. 測繪通報,2008(3):1-4.

[2]龍四春,張詩玉,馮濤,等. 公用主影像干涉圖加權(quán)疊加方法及其在地面沉降監(jiān)測中的應(yīng)用[J].測繪學(xué)報,2012，41(6):844-850

[3]侯建國,初禹. 合成孔徑雷達差分干涉測量技術(shù)在城市地面沉降監(jiān)測中的應(yīng)用[J]. 測繪工程,2014(8):40-44.

[4]張潔,胡光道. 基于合成孔徑雷達干涉測量技術(shù)的地面沉降研究綜述[J]. 地質(zhì)科技情報,2005(3):104-108.

[5]王桂杰,謝謨文,邱騁,等.D-INSAR技術(shù)在大范圍滑坡監(jiān)測中的應(yīng)用[J].巖土力學(xué),2010，31(4):1337-1344.

[6]舒寧. 雷達影像干涉測量原理[M]. 武漢: 武漢大學(xué)出版社, 2003:3-10.

[7]楊成生,瞿偉,季靈運,等.InSAR地理編碼處理中高程估算算法的研究[J]. 遙感信息,2007(6):15-17.

[8]陳蓓蓓, 李小娟, 張有全,等. 基于GPS和MERIS數(shù)據(jù)的InSAR大氣延遲改正[C]∥2008年測繪科學(xué)前沿技術(shù)論壇.南京：[s.n.]，2008.

[9]母景琴,姚國清. 干涉雷達時間序列分析方法在地面沉降監(jiān)測中的應(yīng)用[J]. 國土資源遙感,2007(3):28-30，46.

[10]張振生,孟昆, 谷延群.D-InSAR技術(shù)在礦山沉陷和地面沉降監(jiān)測中的應(yīng)用[J]. 華北地震科學(xué),2006，24(3):47-51.

[11]何亞倩,鄭宇,何晴倩.D-InSAR在礦區(qū)開采沉陷監(jiān)測中應(yīng)用的探析[J]. 現(xiàn)代礦業(yè),2009(9):76-78.

[12]ROGERSAE,INGALLSRP.Venus:MappingtheSurfaceReflectivitybyRadarInterferometry[J].Science, 1969,165(3895): 797-799.

[13]FARRTG,KOBRICKM.ShuttleRadarTopographyMissionProducesaWealthofData[J].Amer.Geophys.UnionEOS, 2000, 81(48): 583-585.

[14]NAKAGAWAH,MURAKAMIM,FUJIWARAS,etal.LandSubsidenceoftheNorthernKantoPlainsCausedbyGroundWaterExtractionDetectedbyJERS2InSARInterferometry[J].InternationalGeoscienceandRemoteSensingSymposium, 2000, 5: 2233-2235.

[15]周琦,趙文吉. 基于D-InSAR三軌法的北京平原區(qū)地表形變監(jiān)測[J]. 地理空間信息,2013，11(2):146-148，154.

[16]屈春燕,單新建,宋小剛,等.D-InSAR技術(shù)應(yīng)用于汶川地震地表位移場的空間分析[J]. 地震地質(zhì),2010，32(2):175-190.

The Methods of Ground Subsidence Monitoring Based on D-InSAR Technique

YUE Shun，YUE Jianping，QIU Zhiwei，WNAG Xueqin

10.13474/j.cnki.11-2246.2016.0160.

2015-06-25

國家自然科學(xué)基金(41174002)；湖南省重點實驗室開放研究基金(PHLHD201311)

岳順(1991—)，男，碩士，研究方向為3S集成與應(yīng)用。E-mail：yue_shun@163.com

P237

B

0494-0911(2016)05-0077-03

引文格式： 岳順，岳建平，邱志偉，等. 基于D-InSAR技術(shù)的地面沉降監(jiān)測方法[J].測繪通報，2016(5)：77-79.