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

（1.河海大學(xué)地球科學(xué)與工程學(xué)院，江蘇南京211000；2.河南城建學(xué)院測(cè)繪學(xué)院，河南平頂山467036）

基于短基線集技術(shù)的城市地表沉降監(jiān)測(cè)研究

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

（1.河海大學(xué)地球科學(xué)與工程學(xué)院，江蘇南京211000；2.河南城建學(xué)院測(cè)繪學(xué)院，河南平頂山467036）

通過(guò)對(duì)合成孔徑雷達(dá)干涉測(cè)量的去相干分析，提出了基于短基線集（SBAS）技術(shù)進(jìn)行地表形變監(jiān)測(cè)的方法。該方法能夠克服DInSAR易受時(shí)間、基線去相干等因素的影響，通過(guò)對(duì)短時(shí)-空基線組合的影像對(duì)進(jìn)行干涉處理以提取高相干點(diǎn)，利用大氣延時(shí)相位與相位噪聲在頻域不同特性以達(dá)到二者分離的目的，最終獲取監(jiān)測(cè)區(qū)域長(zhǎng)時(shí)間緩慢地表形變的演變規(guī)律。本文利用2007—2011年16景南京地區(qū)ALOS數(shù)據(jù)進(jìn)行了短基線差分干涉試驗(yàn)，并通過(guò)實(shí)測(cè)水準(zhǔn)數(shù)據(jù)進(jìn)行了對(duì)比，結(jié)果表明該技術(shù)能夠較準(zhǔn)確地反演出地表的形變場(chǎng)及累積形變量。

短基線集（SBAS）；形變監(jiān)測(cè)；差分干涉；形變場(chǎng)

地表形變監(jiān)測(cè)是進(jìn)行地表形變分析研究并采取相應(yīng)防治對(duì)策的重要基礎(chǔ)。隨著科學(xué)技術(shù)的不斷提高，地表形變監(jiān)測(cè)從以往的以光學(xué)水準(zhǔn)測(cè)量為主，逐步向采用GPS技術(shù)，并進(jìn)而以合成孔徑雷達(dá)干涉測(cè)量（InSAR）為主要方式，從而使其在監(jiān)測(cè)范圍、監(jiān)測(cè)實(shí)效等方面都有了顯著的提高。

InSAR技術(shù)成功地融合了合成孔徑雷達(dá)成像原理和干涉測(cè)量技術(shù)，該技術(shù)能利用傳感器的系統(tǒng)參數(shù)、姿態(tài)參數(shù)和軌道之間的幾何關(guān)系等精確測(cè)量地表某一點(diǎn)的三維空間位置及微小變化。但常規(guī)差分干涉測(cè)量方法（DInSAR）受時(shí)間、空間去相干和大氣延時(shí)相位的影響，嚴(yán)重限制了DInSAR方法在地表形變監(jiān)測(cè)，尤其是長(zhǎng)時(shí)間緩慢地表形變監(jiān)測(cè)中的應(yīng)用［1-2］。短基線集（SBAS）技術(shù)充分利用短時(shí)空基線影像的相干信息，能有效地抑制相位噪聲對(duì)地形相位的影響，從而獲取長(zhǎng)時(shí)間緩慢地表形變的演變規(guī)律。本文以南京市為例對(duì)該技術(shù)在城市沉降監(jiān)測(cè)方面進(jìn)行深入探討和研究。

一、基于短基線集技術(shù)的沉降監(jiān)測(cè)方法

短基線集（SBAS）技術(shù)繼承了常規(guī)DInSAR方法的優(yōu)點(diǎn)，利用較短時(shí)空基線的影像對(duì)產(chǎn)生干涉圖提高相干性，下面對(duì)基于SBAS的地表沉降監(jiān)測(cè)方法進(jìn)行簡(jiǎn)要的論述［3-4］。

假設(shè)（t0，t1，…，tN）時(shí)刻獲取了覆蓋同一區(qū)域N+1幅SAR影像圖，并且假設(shè)所有SAR影像已配準(zhǔn)到同一坐標(biāo)系下，從而可以得到M對(duì)時(shí)空基線均小于某一閾值的多視差分干涉對(duì)，且M滿足以下條件

假設(shè)從tA、tB兩個(gè)時(shí)間獲得的SAR圖像產(chǎn)生第j幅干涉圖，并假設(shè)tB〉tA，去除地形相位后，可建立未考慮大氣相位、地形誤差及失相關(guān)噪聲等的簡(jiǎn)化模型，則干涉圖在像元x處的干涉相位可表示為

式中，λ為雷達(dá)波長(zhǎng)；d（tB，x）、d（tA，x）分別為tA和tB時(shí)刻像元相對(duì)于初始時(shí)刻t0的LOS方向地表形變，即有d（t0，x）=0；假設(shè)相位δφ（x）為解纏后的相位，所有干涉圖經(jīng)過(guò)配準(zhǔn)，并選取相同的解纏起始點(diǎn)（穩(wěn)定點(diǎn)或形變已知點(diǎn)）。該方法對(duì)干涉圖進(jìn)行逐像元的時(shí)間序列分析，因此，以下討論均以某一像元為例來(lái)建立方程。

假設(shè)主影像時(shí)序集IE=｛IE1，IE2，…，IEM｝和從影像集IS=｛IS1，IS2，…，ISM｝，且滿足

則所有差分干涉圖相位可以組成如下觀測(cè)方程

對(duì)所有干涉圖，可將式（4）的線性模型表示為矩陣形式δφ=Aφ，其中A為M×N維矩陣。當(dāng)M≥N時(shí)，則該矩陣秩為N，可對(duì)式（4）進(jìn)行最小二乘法即可求解出φ的估計(jì)值 ^φ

通常為了減少基線去相干影響，會(huì)將干涉對(duì)進(jìn)行分組，這樣矩陣A的M常小于N時(shí)，相應(yīng)法方程系數(shù)陣秩虧，可采用奇異值分解法求解［5-7］。

采用較短基線（通常小于200 m）干涉紋圖集可以降低幾何去相干對(duì)它們的影響。此外，由于較大的高程模糊度也使得它們對(duì)DEM誤差的敏感性降低。通過(guò)累積差分干涉紋圖測(cè)量地殼形變，是該方法與永久散射體相干技術(shù)的共同點(diǎn)，但在去除大氣效應(yīng)時(shí)采用組合多景干涉紋圖平均去除仍然帶有一定的不確定性，因此增加了該方法的復(fù)雜度。該方法的優(yōu)點(diǎn)是可測(cè)量非線性形變，SBAS方法的數(shù)據(jù)處理流程如圖1所示。

圖1　短基線集技術(shù)流程

二、應(yīng)用分析

圖2為試驗(yàn)區(qū)對(duì)應(yīng)的SAR幅度影像，以南京市區(qū)為中心，北至八卦洲，南至將軍山，西至老山森林公園，東至仙林以西，試驗(yàn)范圍幾乎覆蓋整個(gè)南京市轄各區(qū)，如圖3所示。選擇日本ALOS衛(wèi)星的PALSAR 16景升軌數(shù)據(jù)，時(shí)間跨度從2007—2011年。監(jiān)測(cè)區(qū)域植被覆蓋大，地形復(fù)雜，包括山區(qū)、長(zhǎng)江漫灘平原及水體等，數(shù)據(jù)采集密度小，時(shí)間跨度大，易受大氣效應(yīng)影響。試驗(yàn)數(shù)據(jù)的詳細(xì)情況見(jiàn)表1。

圖2　試驗(yàn)研究區(qū)域的SAR幅度

圖3　試驗(yàn)研究區(qū)域的Google Earth影像

表1　ALOS SAR影像數(shù)據(jù)參數(shù)

經(jīng)過(guò)對(duì)各SAR影像進(jìn)行空間及時(shí)間基線的分析與比較，采用以2009年8月23日獲取的影像為主影像進(jìn)行配準(zhǔn)及干涉，從其余各數(shù)據(jù)以此景影像為參照得到的垂直基線數(shù)據(jù)可以看出，最短的垂直基線為328.96 m，最長(zhǎng)則達(dá)到2 982.11 m。由于基線較長(zhǎng)且時(shí)間跨度近5年，因此依據(jù)短基線集技術(shù)原理，在進(jìn)行差分干涉配對(duì)時(shí)進(jìn)行了空間與時(shí)間基線的約束，最終得到了94組干涉對(duì)，詳細(xì)的組合情況如圖4所示。為了充分考慮時(shí)序上相位變化信息，采用了三維解纏方法即加入時(shí)序維進(jìn)行空間時(shí)間三維解纏，SAR影像的Delaunay三維連接圖如圖5所示。

圖4　SAR影像時(shí)間空間分布圖

圖5　SAR影像Delaunay三維連接圖

考慮到空間及時(shí)間基線較大的數(shù)據(jù)特點(diǎn)，通過(guò)對(duì)數(shù)據(jù)集中最大空間基線的干涉對(duì)對(duì)應(yīng)的相干圖（如圖6所示）與最大時(shí)間基線的干涉對(duì)對(duì)應(yīng)的相干圖（如圖7所示）進(jìn)行相干分析，可以看出相干像元分布基本覆蓋南京市區(qū)，在兩種極限情況下城區(qū)依舊表現(xiàn)出較高的相干性，但是水體與有植被的山區(qū)相干性較差，均無(wú)法達(dá)到高相干點(diǎn)的提取要求。在數(shù)據(jù)處理過(guò)程中，為抑制相位噪聲與大氣相位帶來(lái)的影響，選擇相干閾值為0.35，小于該值的相干點(diǎn)認(rèn)為失相干嚴(yán)重，不予采用。這里采用的參考DEM為SRTM3，其精度優(yōu)于16 m，滿足去除地形相位的要求；解纏參考點(diǎn)依據(jù)以下原則進(jìn)行選取:①高相干且無(wú)形變點(diǎn)或已知形變點(diǎn)；②遠(yuǎn)離形變區(qū)，由于SAR數(shù)據(jù)覆蓋江北區(qū)域，為避免由于跨江造成的解纏誤差，此處分別選取兩處參考點(diǎn):紫金山附近地質(zhì)穩(wěn)定點(diǎn)、老山附近地質(zhì)穩(wěn)定點(diǎn)。由配準(zhǔn)后的干涉處理、去平地效應(yīng)及解纏后查驗(yàn)結(jié)果可以發(fā)現(xiàn)，存在嚴(yán)重的地形誤差及因時(shí)間失相干引起相位誤差，這些結(jié)果會(huì)影響線形形變信息的正確提取，因此需要從結(jié)果中予以去除。經(jīng)過(guò)篩查后，得到結(jié)果較為準(zhǔn)確的78組干涉對(duì)來(lái)進(jìn)行后續(xù)的大氣分離和形變反演。圖8為獲取的南京市及其周邊的形變速度場(chǎng)。

圖6　最長(zhǎng)空間基線干涉對(duì)相干圖

圖7　最長(zhǎng)時(shí)間基線干涉對(duì)相干圖

圖8　研究區(qū)域視線向形變速度場(chǎng)（mm/a）

從形變速度場(chǎng)圖中可以看出，存在4個(gè)較為明顯的沉降區(qū)域:C區(qū)為秦淮河以西的地區(qū)（簡(jiǎn)稱河西），D區(qū)為月牙湖附近，E為南京長(zhǎng)江二橋及南岸一帶，F(xiàn)為江北南京高新技術(shù)產(chǎn)業(yè)開(kāi)發(fā)區(qū)附近。為驗(yàn)證利用SBAS技術(shù)進(jìn)行研究區(qū)域地表沉降監(jiān)測(cè)的精度，本文重點(diǎn)結(jié)合河西區(qū)2007—2009年的地面監(jiān)測(cè)點(diǎn)的時(shí)序沉降數(shù)據(jù)進(jìn)行對(duì)比分析（如圖9所示），這里選擇集慶門(mén)西監(jiān)測(cè)點(diǎn)、濱江風(fēng)光帶監(jiān)測(cè)點(diǎn)以及鼓樓區(qū)匯賢居監(jiān)測(cè)點(diǎn)進(jìn)行累積沉降對(duì)比。從序列圖可以看出，3個(gè)監(jiān)測(cè)點(diǎn)的地面監(jiān)測(cè)結(jié)果與SAR累積沉降序列基本一致。由于監(jiān)測(cè)的時(shí)間與SAR數(shù)據(jù)獲取的觀測(cè)時(shí)間并不對(duì)應(yīng)，需要將地面水準(zhǔn)數(shù)據(jù)進(jìn)行插值處理，難免會(huì)出現(xiàn)插值誤差，且由于解纏參考點(diǎn)的選取及殘留大氣相位噪聲等誤差的存在，使SAR數(shù)據(jù)得到的累積沉降量無(wú)法與地面監(jiān)測(cè)數(shù)據(jù)完全吻合，但經(jīng)時(shí)序?qū)Ρ瓤梢钥闯觯袒€集技術(shù)能夠準(zhǔn)確地獲得試驗(yàn)區(qū)的沉降趨勢(shì)。

圖9　地面監(jiān)測(cè)結(jié)果與SAR累積沉降序列對(duì)比

進(jìn)行監(jiān)測(cè)點(diǎn)對(duì)比驗(yàn)證后，通過(guò)對(duì)所有地面監(jiān)測(cè)數(shù)據(jù)進(jìn)行空間插值得到河西地區(qū)地表沉降趨勢(shì)分布圖，與SAR獲取的地表沉降速度場(chǎng)進(jìn)行對(duì)比驗(yàn)證，結(jié)果如圖10所示。從圖中可以看出，河西地區(qū)的沉降存在自北向南較明顯的條帶分布特點(diǎn)，從龍江經(jīng)莫愁湖到集慶門(mén)都有較為顯著的沉降趨勢(shì)，因此，地面監(jiān)測(cè)結(jié)果與SAR獲取的沉降趨勢(shì)基本一致，但由于SAR獲取的信息豐富，能夠得到更多的沉降細(xì)節(jié)，相對(duì)而言，地面監(jiān)測(cè)點(diǎn)由于點(diǎn)位分布稀少，通過(guò)插值獲取的分布圖會(huì)有較多的孤島區(qū)域，與實(shí)際情況不符。

圖10　地面監(jiān)測(cè)沉降趨勢(shì)與SAR地表形變場(chǎng)對(duì)比

結(jié)合本地區(qū)的地質(zhì)結(jié)構(gòu)，可以得到地面沉降有以下幾個(gè)主要原因:

1）C、E地區(qū)主要是長(zhǎng)江漫灘，軟土層較厚，固結(jié)沉降仍然在進(jìn)行中，具有較小的地基承載力。C地區(qū)主要是由于城市不斷發(fā)展，地上建筑迅速增加，地面荷載對(duì)其影響明顯，且伴隨著地下空間的利用，致使該地區(qū)產(chǎn)生較大的沉降變形；E地區(qū)是南京市新生圩港，常有萬(wàn)噸級(jí)貨輪在此進(jìn)行裝卸活動(dòng)，并且臨近長(zhǎng)江，易受江水水位變化影響。

2）D、F地區(qū)主要是由于城市建設(shè)較快，建筑荷載和工程施工中的抽取地下水可能引起地基不均勻的變形。需特別注意的是F區(qū)臨近山區(qū)，地質(zhì)結(jié)構(gòu)理應(yīng)穩(wěn)定，但依然出現(xiàn)沉降，可能與當(dāng)?shù)厣襟w開(kāi)挖活動(dòng)有關(guān)。

3）從沉降速率圖中可以發(fā)現(xiàn)，盡管南京一橋與二橋兩座橋體上的高相干點(diǎn)較少，但不難發(fā)現(xiàn)均出現(xiàn)較明顯的沉降，有關(guān)單位應(yīng)引起重視加強(qiáng)監(jiān)測(cè)，以免出現(xiàn)安全隱患。

三、結(jié)束語(yǔ)

本文利用短基線集技術(shù)對(duì)南京地區(qū)進(jìn)行了地表沉降分析，通過(guò)地面監(jiān)測(cè)數(shù)據(jù)進(jìn)行對(duì)比驗(yàn)證，得到了較為準(zhǔn)確的該地區(qū)沉降速率分布圖和時(shí)序沉降累積量。試驗(yàn)結(jié)果證明該技術(shù)能夠大面積、高密度地獲取地表形場(chǎng)，克服SAR數(shù)據(jù)基線及時(shí)間去相干的影響，提取出有效的形變信息。但由于技術(shù)理論方面的限制，監(jiān)測(cè)精度還有待進(jìn)一步提高。

［1］ 何秀鳳，何敏.InSAR對(duì)地觀測(cè)數(shù)據(jù)處理方法［M］.北京:科學(xué)出版社，2012:20-30.

［2］ 岳建平，方露，黎昵.變形監(jiān)測(cè)理論與技術(shù)研究進(jìn)展［J］.測(cè)繪通報(bào)，2007（7）:26-31.

［3］ 盧麗君.基于時(shí)序SAR影像的地表形變檢測(cè)方法及其應(yīng)用［D］.武漢:武漢大學(xué)，2008.

［4］ STROZZI T，F(xiàn)ARINA P，CORSINI A，et al.Survey and Monitoring of Landslide Displacements by Means of L-

Research on Urban Surface Subsidence Monitoring Based on SBAS

QIU Zhiwei，YUE Jianping，WANG Xueqin，YUE Shun，LIU Bin

10.13474/j.cnki.11-2246.2016.0217.

P237

B

0494-0911（2016）07-0025-05

2015-07-17；

2015-12-11

河南城建學(xué)院科研基金（2013JZD008）；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金（2014B38314）；水能資源利用關(guān)鍵技術(shù)湖南省重點(diǎn)實(shí)驗(yàn)室開(kāi)放研究基金（PKLHD201311）

邱志偉（1985—），男，博士生，講師，主要研究方向?yàn)槔走_(dá)干涉測(cè)量等。E-mail:qiuzhiwei-2008@163.com

引文格式：邱志偉，岳建平，汪學(xué)琴，等.基于短基線集技術(shù)的城市地表沉降監(jiān)測(cè)研究［J］.測(cè)繪通報(bào)，2016（7）:25-29.