基于SBAS-InSAR的成都平原地面沉降監(jiān)測(cè)

孫曉鵬， 魯小丫， 文學(xué)虎， 甄艷， 王蕾

(1.國(guó)家測(cè)繪地理信息局第六地形測(cè)量隊(duì)四川省地理國(guó)情監(jiān)測(cè)工程技術(shù)研究中心,成都　610500；2.西南民族大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院,成都　610041)

?

基于SBAS-InSAR的成都平原地面沉降監(jiān)測(cè)

孫曉鵬1， 魯小丫2， 文學(xué)虎1， 甄艷1， 王蕾1

(1.國(guó)家測(cè)繪地理信息局第六地形測(cè)量隊(duì)四川省地理國(guó)情監(jiān)測(cè)工程技術(shù)研究中心,成都610500；2.西南民族大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院,成都610041)

摘要：汶川地震后，余震活動(dòng)頻繁，加之成都平原內(nèi)城市發(fā)展迅速，容易誘發(fā)地面沉降； 對(duì)成都平原地面沉降進(jìn)行監(jiān)測(cè)，及時(shí)掌握沉降信息，可為相關(guān)決策提供科學(xué)依據(jù)?；贓NVISATASAR數(shù)據(jù)，采用小基線集(smallbaselinesubset,SBAS)-InSAR技術(shù)，對(duì)成都平原2008―2010年地面沉降進(jìn)行了監(jiān)測(cè)。結(jié)果表明，各主要城市在監(jiān)測(cè)時(shí)段內(nèi)的地表累積形變量在-8～14mm之間，總體形變量不大； 成都平原西部區(qū)域受地震影響呈抬升趨勢(shì)，沉降主要集中于成都市北側(cè)和德陽(yáng)市以南部分區(qū)域，最大沉降量為-22mm，沉降范圍隨時(shí)間推移呈擴(kuò)大趨勢(shì)。通過(guò)實(shí)測(cè)數(shù)據(jù)驗(yàn)證了監(jiān)測(cè)結(jié)果，精度達(dá)到2.9mm。成都平原不存在區(qū)域性沉降的構(gòu)造背景，且地下水資源豐富，沉降自然誘因不明顯，城市建設(shè)活動(dòng)可能為沉降的人為誘因。該成果可為今后成都平原主要城市更加精細(xì)的地面沉降監(jiān)測(cè)工作提供參考。

關(guān)鍵詞：SBAS-InSAR； 成都平原； 地面沉降監(jiān)測(cè)

0引言

地面沉降是在自然和人為因素作用下，因地殼表層土體壓縮而導(dǎo)致區(qū)域性地面標(biāo)高降低的一種環(huán)境地質(zhì)現(xiàn)象[1]，是一種累進(jìn)性的緩變地質(zhì)災(zāi)害，一旦形成便難以恢復(fù)[2]。引起地面沉降的因素除地表松散或半松散地層的自固壓密、地質(zhì)構(gòu)造作用、巖溶發(fā)育地區(qū)的巖溶塌陷等自然因素外，大量開(kāi)采地下資源(地下水、石油及天然氣等)、固體礦產(chǎn)資源和地?zé)豳Y源，大規(guī)模工程建設(shè)、軌道交通、城市地下空間開(kāi)發(fā)及高層建筑物對(duì)地基施加的動(dòng)態(tài)載荷是引起地面沉降的重要人為因素。《2011—2020年全國(guó)地面沉降防治規(guī)劃》(以下簡(jiǎn)稱《規(guī)劃》)指出，目前全國(guó)遭受地面沉降災(zāi)害的城市超過(guò)50個(gè)，分布于北京、天津、河北、山西及內(nèi)蒙古等20個(gè)省(市、自治區(qū))。全國(guó)累計(jì)地面沉降量超過(guò)200mm的地區(qū)達(dá)到7.9萬(wàn)km2，并有進(jìn)一步擴(kuò)大趨勢(shì)?！兑?guī)劃》提出了未來(lái)10a的地面沉降監(jiān)測(cè)總體目標(biāo)，即查明全國(guó)地面沉降災(zāi)害的現(xiàn)狀、發(fā)展趨勢(shì)、形成原因和分布規(guī)律。

成都平原地下水資源豐富，多年未見(jiàn)大范圍自然因素導(dǎo)致的地面沉降事件，但龍門山斷裂帶沿線地震活動(dòng)頻繁，對(duì)區(qū)域構(gòu)造造成重要影響； 人為因素導(dǎo)致的沉降甚至塌陷時(shí)有發(fā)生，如2014年7月9日成都暴雨致棕樹(shù)南街一處停車場(chǎng)發(fā)生塌陷事故，4輛轎車連同地面上的幾棵大樹(shù)一同墜入數(shù)米深的基坑內(nèi)。然而，目前針對(duì)成都平原形變監(jiān)測(cè)的相關(guān)研究還較鮮見(jiàn)，因此有必要對(duì)成都平原主要城市進(jìn)行地面沉降的監(jiān)測(cè)工作。

目前，使用合成孔徑雷達(dá)差分干涉測(cè)量(differentialinterferometrysyntheticapertureRadar,D-InSAR) 技術(shù)，能夠進(jìn)行大范圍、無(wú)接觸、面狀的mm級(jí)地表形變監(jiān)測(cè)。近年來(lái)，國(guó)內(nèi)研究者使用衛(wèi)星雷達(dá)干涉技術(shù)進(jìn)行了較多的沉降監(jiān)測(cè)研究，例如楊成生等[3]使用ASAR數(shù)據(jù)，采用小基線集(smallbaselinesubset,SBAS)監(jiān)測(cè)大同盆地的地面沉降，并進(jìn)行了地下水位與沉降之間關(guān)系的定量分析； 劉志敏等[4]使用SBAS技術(shù)進(jìn)行了長(zhǎng)治礦區(qū)2003—2010年地面沉降監(jiān)測(cè)。SBAS方法將SAR數(shù)據(jù)組成若干個(gè)子集，使用多景主圖像，能夠較好地克服時(shí)空失相干限制，獲得大范圍mm級(jí)精度的地表形變信息。因此，本文采用SBAS-InSAR技術(shù)監(jiān)測(cè)成都平原沉降的分布特征及時(shí)間演化規(guī)律，以期為成都平原主要城市更加精細(xì)的地面沉降監(jiān)測(cè)工作提供參考。

1數(shù)據(jù)處理與分析

1.1研究區(qū)及遙感數(shù)據(jù)

本研究的監(jiān)測(cè)范圍為成都平原部分地區(qū)。ENVISATASAR圖像覆蓋范圍如圖1所示。

圖1　ASAR圖像覆蓋范圍示意圖

用于監(jiān)測(cè)的ENVISARASAR圖像來(lái)源于歐空局(EuropeanSpaceAgency,ESA)，軌道號(hào)為Track18，時(shí)相分別為2008-05-28,2008-08-06,2008-09-10,2008-10-15,2009-01-28,2009-03-04,2009-04-08,2009-05-13, 2009-09-30,2010-01-13,2010-02-17,2010-03-24 ,2010-04-28 ,2010-06-02和2010-08-11。其時(shí)相與空間基線的關(guān)系見(jiàn)圖2。

圖2　ASAR時(shí)相與空間基線關(guān)系

1.2SBAS-InSAR技術(shù)原理

SBAS-InSAR是由Berardino和Lanari等研究人員提出的一種時(shí)間序列InSAR分析方法[5]。該方法將獲取的數(shù)據(jù)進(jìn)行適當(dāng)組合，得到一系列短空間基線差分干涉圖，能較好地克服空間去相關(guān)現(xiàn)象。在求解形變速率時(shí)，SBAS方法采用了奇異值分解(singularvaluedecomposition,SVD)法，可以將被較大空間基線分開(kāi)的孤立SAR數(shù)據(jù)集連接起來(lái)，從而提高觀測(cè)數(shù)據(jù)的采樣率。其主要步驟[6]如下。

1)在(t0,t1,…,tn)時(shí)間內(nèi)，對(duì)在同一地區(qū)獲取的(N+1)期SAR圖像，選取其中1景圖像作為主圖像，將其他SAR圖像與主圖像配準(zhǔn)。假設(shè)每期SAR圖像都有至少有1景圖像與之干涉，則(N+1)期圖像生成M景差分干涉圖，M滿足條件

(N+1)/2≤M≤N(N+1)/2 。

(1)

2)對(duì)于從圖像時(shí)刻tA和主圖像時(shí)刻tB(tB>tA)獲取的SAR圖像，生成的第j景差分干涉圖的干涉相位(忽略大氣延遲相位、殘余地形相位和噪聲相位)可以表示為

(2)

式中： φ為干涉相位； j為差分干涉圖的景號(hào)，j∈(1,…,M)； λ為信號(hào)的中心波長(zhǎng)； d(tA,x,r)和d(tB,x,r)分別為tA時(shí)刻和tB時(shí)刻相對(duì)于d(t0,x,r)=0的雷達(dá)視線方向累積形變量。

3)將式(2)中的相位表示為2個(gè)獲取時(shí)間之間的平均相位速度vj和時(shí)間的乘積，即

vj=(φj-φj-1)/(tj-tj-1) ,

(3)

則第j景干涉圖的相位值為

(4)

即各時(shí)段速度在主、從圖像時(shí)間間隔上的積分，寫(xiě)成矩陣形式為

Bv=δφ。

(5)

在求解系數(shù)矩陣B的過(guò)程中，由于SBAS的差分干涉圖采用了多主圖像策略，因此矩陣B可能會(huì)秩虧。采用SVD方法可以得到矩陣B的廣義逆矩陣，進(jìn)而得到速度矢量的最小范數(shù)解，最后通過(guò)各個(gè)時(shí)段速度的積分得到各個(gè)時(shí)間段的形變量。

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

SBAS-InSAR數(shù)據(jù)處理流程主要包括：

1)干涉對(duì)組合與干涉圖生成。根據(jù)設(shè)定的時(shí)間和空間基線閾值，將全部SAR圖像生成干涉數(shù)據(jù)對(duì)，并計(jì)算生成差分干涉圖，后續(xù)處理包括去平、濾波和解纏。

2)相干目標(biāo)選擇。以空間相干計(jì)算結(jié)果為參考，通過(guò)設(shè)定最小相干閾值，進(jìn)行目標(biāo)的分析與選擇。

3)差分干涉圖解纏。在相干性高且不存在干涉條紋(或者干涉條紋稀疏)的區(qū)域內(nèi)選擇GCP控制點(diǎn)，作為形變信息的對(duì)比參考點(diǎn)，使用Delaunay3D方法進(jìn)行解纏。

4)時(shí)間序列形變結(jié)果生成。利用SVD方法對(duì)解纏相位進(jìn)行解算，得到各個(gè)時(shí)段內(nèi)的形變速率，最后通過(guò)積分獲得各個(gè)時(shí)段的形變量(負(fù)值表示沉降量，正值表示抬升量)。

具體處理流程如圖3所示。

圖3　SBAS數(shù)據(jù)處理流程圖

1.4地表形變監(jiān)測(cè)結(jié)果

使用SABS-InSAR技術(shù)共獲得14期相對(duì)于2008-05-28獲取圖像的地表累積形變量圖，覆蓋范圍包括成都市主城區(qū)及新都區(qū)、郫縣、溫江區(qū)、青白江區(qū)、龍泉驛區(qū)、金堂縣、新津縣、雙流縣和彭州市，德陽(yáng)市旌陽(yáng)區(qū)、綿竹市、什邡市、廣漢市和羅江縣，綿陽(yáng)市涪城區(qū)、仙游區(qū)和江油市以及簡(jiǎn)陽(yáng)市。地表形變監(jiān)測(cè)結(jié)果如圖4所示。

(a) 2008-08-06(b) 2008-09-10(c) 2008-10-15(d) 2009-01-28 (e) 2009-03-04

(f) 2009-04-08(g) 2009-05-13(h) 2009-09-30(i) 2010-01-13(j) 2010-02-17

(k) 2010-03-24(l) 2010-04-28(m) 2010-06-02(n) 2010-08-11

圖4監(jiān)測(cè)范圍內(nèi)14期累積形變量

Fig.4Accumulateddeformationof14periodsinmonitoringarea

從圖4可以看出，監(jiān)測(cè)范圍內(nèi)地面沉降主要分布在成都市北側(cè)和德陽(yáng)市以南的區(qū)域，最大沉降量為-22mm，并且沉降范圍隨時(shí)間的推移逐漸擴(kuò)大。

以區(qū)縣、市為行政單元，統(tǒng)計(jì)了2008-05-28至2010-08-11每2a期間累積形變量的平均值(圖5)。

圖5　監(jiān)測(cè)范圍內(nèi)各個(gè)縣市平均累積形變量

在監(jiān)測(cè)范圍內(nèi)，地面抬升區(qū)主要集中于成都平原西部，包括成都市新津縣，德陽(yáng)市所轄旌陽(yáng)區(qū)、羅江縣、綿竹市，綿陽(yáng)市所轄江油市、游仙區(qū)和涪城區(qū)，這是由于自汶川地震主震發(fā)生后至2010-09-08期間共發(fā)生4級(jí)以上的地震300余次，其中發(fā)震以逆沖型為主[7]，致使位于上盤一側(cè)的成都平原西部區(qū)域受到逆沖作用的影響而逐漸隆升。沉降區(qū)域分布于成都市所轄青白江區(qū)、溫江區(qū)、金牛區(qū)、成華區(qū)、新都區(qū)和雙流縣； 成都市青羊區(qū)、錦江區(qū)、武侯區(qū)、金堂縣、郫縣、什邡市、廣漢市和簡(jiǎn)陽(yáng)市的形變量都接近于零。盡管監(jiān)測(cè)范圍內(nèi)有不同程度的沉降和抬升，但2a期間累積形變量在-8～14mm之間，總體形變量不大。監(jiān)測(cè)范圍內(nèi)，平均形變速率在-10.3～11.1mm/a之間(圖6)。

圖6　平均形變速率

圖7示出成都、德陽(yáng)、綿陽(yáng)市各區(qū)縣的平均累積形變量。

(a) 成都市 (b) 德陽(yáng)市 (c) 綿陽(yáng)市

圖7不同成市各區(qū)縣平均累積形變量

Fig.7Averageaccumulateddeformationofvariouscountiesindifferentcities

從圖7可以看出，成都市青白江區(qū)、溫江區(qū)、金牛區(qū)、成華區(qū)、新都區(qū)和雙流縣的平均累積形變量在監(jiān)測(cè)時(shí)段內(nèi)為負(fù)，呈現(xiàn)沉降趨勢(shì)； 反之，新津縣呈現(xiàn)抬升趨勢(shì)； 其他行政區(qū)則在0值附近波動(dòng)。德陽(yáng)市所轄廣漢市為沉降區(qū)域，但是沉降值在-6～0mm之間，沉降量較小； 什邡市在2008年8月抬升后呈現(xiàn)沉降趨勢(shì)； 其他區(qū)縣為抬升區(qū)域。綿陽(yáng)市各行政區(qū)在監(jiān)測(cè)時(shí)段內(nèi)均呈抬升趨勢(shì)，但累積抬升量在2010年3月開(kāi)始有所下降。同時(shí)，以上3市在2008年下半年的平均累積形變量均存在較大波動(dòng)，這是因

為在2008-05-12至2008-10-31共發(fā)生4級(jí)以上余震240余次，地震活動(dòng)頻繁，地面受地震影響形變波動(dòng)較大。

2精度驗(yàn)證

選取覆蓋成都市的連續(xù)運(yùn)行衛(wèi)星定位服務(wù)綜合系統(tǒng)(continuousrunningsatelliteorientationserviceintegratedsystem,CROS)數(shù)據(jù)對(duì)地表形變監(jiān)測(cè)結(jié)果進(jìn)行驗(yàn)證，驗(yàn)證點(diǎn)的空間分布如圖8所示。

圖8　監(jiān)測(cè)驗(yàn)證點(diǎn)分布圖

本文收集到成都地區(qū)名稱為JITA,CDKC,LOQU和PUJI共4個(gè)CROS站點(diǎn)的數(shù)據(jù)。其中，PUJI站點(diǎn)位于監(jiān)測(cè)范圍之外；JITA和LOQU這2個(gè)站點(diǎn)位于山區(qū)，相干性較差，沒(méi)有形成有效的形變監(jiān)測(cè)驗(yàn)證； 故只能使用CDKC站數(shù)據(jù)進(jìn)行驗(yàn)證。另外，還收集到“汶川地震災(zāi)后重建——測(cè)繪專項(xiàng)1∶2 000測(cè)圖”項(xiàng)目提供名稱為QY10的水準(zhǔn)點(diǎn)作為驗(yàn)證點(diǎn)。驗(yàn)證數(shù)據(jù)點(diǎn)位的SBAS-InSAR形變監(jiān)測(cè)結(jié)果與實(shí)地測(cè)量形變結(jié)果的對(duì)比如表1所示。

驗(yàn)證結(jié)果表明，監(jiān)測(cè)結(jié)果與實(shí)測(cè)結(jié)果有較好的一致性，精度達(dá)到2.971mm。

表1　SBAS-InSAR形變監(jiān)測(cè)結(jié)果與實(shí)地測(cè)量形變差異比較

3問(wèn)題與討論

1)是否存在地面沉降的自然誘因。就沉降原因而言，根據(jù)對(duì)成都平原形成及演化的相關(guān)研究[8-11]可知： ①喜山運(yùn)動(dòng)以來(lái)，四川盆地(成都平原)處于擠壓構(gòu)造背景下的非拉張環(huán)境； ②剛性塊體被周緣構(gòu)造帶所圍限，而成都平原并非前陸坳陷區(qū)，因此不具備區(qū)域沉降的構(gòu)造條件； ③成都平原基底正斷裂不發(fā)育或不活動(dòng)，不存在區(qū)域性沉降的誘因； ④成都平原中、新生代經(jīng)歷海盆—湖盆—陸盆的長(zhǎng)期演化，陸源碎屑巨厚沉積，沉積環(huán)境穩(wěn)定，更新世以來(lái)未經(jīng)歷區(qū)域沉降歷史； ⑤汶川地震及余震以逆沖型為主，致使位于上盤一側(cè)的成都平原有逐漸抬升的趨勢(shì)。因此，監(jiān)測(cè)時(shí)段內(nèi)成都平原不存在地面沉降的自然誘因。

2)區(qū)域降水對(duì)地面沉降的貢獻(xiàn)。從中國(guó)氣象數(shù)據(jù)網(wǎng)收集到的監(jiān)測(cè)范圍內(nèi)編號(hào)為56196和56188的氣象臺(tái)站在2008年1月至2010年11月期間監(jiān)測(cè)的月降水?dāng)?shù)據(jù)(圖9)來(lái)看，2個(gè)臺(tái)站在監(jiān)測(cè)時(shí)段內(nèi)監(jiān)測(cè)到的平均年降水量分別為832mm和1 036mm，降水較為豐富； 且成都市水資源承載力與年降雨量呈顯著的正相關(guān)關(guān)系，從總量上來(lái)說(shuō)，基本上能滿足成都市需水總量的要求[12]。另外，成都平原多年平均地下水資源模數(shù)量為72.02萬(wàn)m3/(km2·a)，地下水系統(tǒng)很穩(wěn)定[13]，因而不存在地下水超采導(dǎo)致地面沉降的問(wèn)題。

圖9　2個(gè)氣象臺(tái)站監(jiān)測(cè)的月降水量

3)是否存在地面沉降的人為誘因?？紤]到成都市城市建設(shè)發(fā)展迅速，近年來(lái)地鐵和城市二環(huán)立交等施工活動(dòng)頻繁； 因此，城市建設(shè)活動(dòng)有可能成為地面沉降的人為誘因，今后需要對(duì)其密切關(guān)注。

4結(jié)論

利用SBAS-InSAR方法，得到了成都平原地區(qū)2008―2010年地面沉降信息。監(jiān)測(cè)結(jié)果表明：

1)成都平原累積形變量平均值在-8～14mm之間，汶川地震后總體形變量不大。

2)受汶川地震及其余震影響，成都平原西部區(qū)域受到逆沖作用的影響而逐漸隆升； 沉降則主要集中于成都市北側(cè)和德陽(yáng)市以南部分地區(qū)，監(jiān)測(cè)期間區(qū)內(nèi)尚沒(méi)有出現(xiàn)沉降漏斗，沉降量較小，但沉降范圍呈逐漸擴(kuò)大的趨勢(shì)。

需要說(shuō)明的是，因未能收集到研究區(qū)域內(nèi)較為詳盡的與地質(zhì)、地下水相關(guān)的數(shù)據(jù)，故未能針對(duì)沉降成因進(jìn)行深入分析。然而，與地震活動(dòng)相關(guān)的資料較好地吻合了研究范圍內(nèi)的抬升區(qū)域，而沉降區(qū)則位于城市發(fā)展迅速的成都市周邊； 同時(shí)，監(jiān)測(cè)結(jié)果與實(shí)地測(cè)量結(jié)果有較好的一致性，其結(jié)果可為今后成都平原主要城市更加精細(xì)的地表形變監(jiān)測(cè)工作提供參考。

參考文獻(xiàn)(References)：

[1]鄭銑鑫,武強(qiáng),侯艷聲,等.關(guān)于城市地面沉降研究的幾個(gè)前沿問(wèn)題[J].地球?qū)W報(bào),2002,23(3):279-282.

ZhengXX,WuQ,HouYS,etal.Somefrontierproblemsonlandsubsidenceresearch[J].ActaGeoscientiaSinica,2002,23(3):279-282.

[2]殷躍平,張作辰,張開(kāi)軍.我國(guó)地面沉降現(xiàn)狀及防治對(duì)策研究[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào),2005,16(2):1-8.

YinYP,ZhangZC,ZhangKJ.LandsubsidenceandcountermeasuresforitspreventioninChina[J].TheChineseJournalofGeologicalHazardandControl,2005,16(2):1-8.

[3]楊成生,劉媛媛,敖萌.基于SBAS時(shí)序分析的大同地面沉降與地下水活動(dòng)研究[J].國(guó)土資源遙感,2015,27(1):127-132.doi:10.6046/gtzyyg.2015.01.20.

YangCS,LiuYY,AoM.StudyoflandsubsidenceandgroundwateractivityusingSBAStime-seriesanalysis[J].RemoteSensingforLandandResources,2015,27(1):127-132.doi:10.6046/gtzyyg.2015.01.20.

[4]劉志敏,李永生,張景發(fā),等.基于SBAS-InSAR的長(zhǎng)治礦區(qū)地表形變監(jiān)測(cè)[J].國(guó)土資源遙感,2014,26(3):37-42.doi:10.6046/gtzyyg.2014.03.06.

LiuZM,LiYS,ZhangJF,etal.AnanalysisofsurfacedeformationintheChangzhiminingareausingsmallbaselineInSAR[J].RemoteSensingforLandandResources,2014,26(3):37-42.doi:10.6046/gtzyyg.2014.03.06.

[5]尹宏杰,朱建軍,李志偉,等.基于SBAS的礦區(qū)形變監(jiān)測(cè)研究[J].測(cè)繪學(xué)報(bào),2011,40(1):52-58.

YinHJ,ZhuJJ,LiZW,etal.GroundsubsidencemonitoringinminingareausingDInSARSBASalgorithm[J].ActaGeodaeticaetCartographicaSinica,2011,40(1):52-58.

[6]BerardinoP,FornaroG,LanariR,etal.AnewalgorithmforsurfacedeformationmonitoringbasedonsmallbaselinedifferentialSARinterferograms[J].IEEETransactionsonGeoscienceandRemoteSensing,2002,40(11):2375-2382.

[7]廖明生,王騰.時(shí)間序列InSAR技術(shù)與應(yīng)用[M].北京:科學(xué)出版社,2014:84-86.

LiaoMS,WangT.TimeSeriesInSARTechnologyandApplication[M].Beijing:SciencePress,2014:84-86.

[8]林向東.汶川地震以來(lái)龍門山斷裂帶地震矩張量時(shí)空演化特征[D].北京:中國(guó)地震局地球物理研究所,2014:84-86.

LinXD.LongmenShanFaultSeismicMomentTensorCharacteristicTemporalEvolutionSincetheEarthquake[D].Beijing:InstituteofGeophysics,ChinaEatrhquakeAdministration,2014:84-86.

[9]錢洪,唐榮昌.成都平原的形成與演化[J].四川地震,1997(3):1-7.

QianH,TangRC.OntheformationandevolutionoftheChengduPlain[J].EarthquakeResearchinSichuan,1997(3):1-7.

[10]何銀武.試論成都盆地(平原)的形成[J].中國(guó)區(qū)域地質(zhì),1987(2):169-175.

HeYW.OntheformationoftheChengduBasin(Plain)[J].RegionalGeologyofChina,1987(2):169-175.

[11]李永昭,郭兵.成都平原的晚新生代構(gòu)造[J].成都理工大學(xué)學(xué)報(bào):自然科學(xué)版,2008,35(4):371-376.

LiYZ,GuoB.CenozonictectonicsofChengduPlain,Sichuan,China[J].JournalofChengduUniversityofTechnology:Science&TechnologyEdition,2008,35(4):371-376.

[12]歐洋銘,白若男,朱國(guó)宇.成都市水資源承載力動(dòng)態(tài)變化及成因[J].南水北調(diào)與水利科技,2012,10(6):110-114.

OuYM,BaiRN,ZhuGY.DynamicvariationofwaterresourcescarryingcapacityanditscausesinChengdu[J].South-to-NorthWaterTransfersandWaterScienceandTechnology,2012,10(6):110-114.

[13]鄭義加.成都平原地下水資源分布的初步分析[J].四川水利,2005,26(4):28-31.

ZhengYJ.TentativeanalysisofgroundwaterresourcesdistributionintheChengduPlain[J].SichuanWaterResources,2005,26(4):28-31.

(責(zé)任編輯： 劉心季)

Monitoring of ground subsidence in Chengdu Plain using SBAS-InSAR

SUN Xiaopeng1， LU Xiaoya2， WEN Xuehu1， ZHEN Yan1， WANG Lei1

(1. Geographic National Condition Monitoring Engineering Research Center of Sichuan Province， No.6 Topographic Survey Team，National Mapping Geographic Information Bureau， Chengdu 610500， China； 2. School of Computer Science and Technology， Southwest University for Nationalities， Chengdu 610041， China)

Abstract:AfterWenchuanearthquake,theaftershockshappenedfrequently.Coupledwiththerapiddevelopmentofcities,theyarelikelytocausegroundsubsidence.SomonitoringthesurfaceofChengduPlainandobtainingthedeformationinformationcouldprovidescientificbasisfortherelevantdecisions.Inthispaper,basedonENVISATASARdata,theauthorsmonitoredthegroundsubsidenceofChengduPlainfrom2008to2010usingSBAS-InSARtechnology.Theresultsshowthattheaveragesurfacedeformationswerebetween-8to14mminmajorcitiesintheChengduPlainduringthemonitoringperiod,thedeformationisnotprominent,andthewesternpartoftheplainshowedauplifttrendcausedbyearthquake.ThesubsidenceareainthenorthoftheChengduCityandsouthtotheDeyangCitywasupto-22mmwiththeexpansionofthesubsidencearea.Themonitoringresultswerevalidatedbymeasureddataandtheaccuracyis2.9mm.TheChengduPlainhasnoregionaltectonicsettingofsubsidenceandhasabundantgroundwaterresources,sothenaturalcauseofsubsidenceisnotobvious;thecitybuildingactivitymightbethemajorcauseofsubsidence.

Keywords:SBAS-InSAR；ChengduPlain；groundsubsidencemonitoring

doi:10.6046/gtzyyg.2016.03.20

收稿日期：2015-04-08；

修訂日期：2015-05-30

基金項(xiàng)目：四川省地理國(guó)情監(jiān)測(cè)工程技術(shù)研究中心開(kāi)放基金項(xiàng)目“基于時(shí)間序列InSAR技術(shù)的紫平鋪水庫(kù)壩體及消落帶形變監(jiān)測(cè)方法研究”(編號(hào)：GC201508)、2014年中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)基金項(xiàng)目“基于D-InSAR技術(shù)的汶川地震重點(diǎn)區(qū)域地表形變監(jiān)測(cè)關(guān)鍵技術(shù)的研究”(編號(hào)： 2014NZYQN28)和數(shù)字制圖與國(guó)土信息應(yīng)用重點(diǎn)實(shí)驗(yàn)室開(kāi)放課題項(xiàng)目“基于開(kāi)源環(huán)境的地理國(guó)情監(jiān)測(cè)統(tǒng)計(jì)成果數(shù)據(jù)可視化展示方法研究與技術(shù)實(shí)現(xiàn)”(編號(hào)：DM2013SC09)共同資助。

中圖法分類號(hào)：TP79

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1001-070X(2016)03-0123-07

第一作者簡(jiǎn)介：孫曉鵬(1986-),男,碩士研究生,助理工程師,主要從事地理國(guó)情監(jiān)測(cè)和相關(guān)科研工作。Email：sunxiaopeng05@163.com。

引用格式： 孫曉鵬,魯小丫,文學(xué)虎,等.基于SBAS-InSAR的成都平原地面沉降監(jiān)測(cè)[J].國(guó)土資源遙感,2016,28(3)：123-129.(SunXP,LuXY,WenXH,etal.MonitoringofgroundsubsidenceinChengduPlainusingSBAS-InSAR[J].RemoteSensingforLandandResources,2016,28(3)：123-129.)