鄢繼選，鄭志琴

（1. 甘肅農(nóng)業(yè)大學水利水電工程學院，甘肅蘭州730070；2. 甘肅農(nóng)業(yè)大學管理學院，甘肅 蘭州730070）

引言

我國是世界上地質(zhì)災害最嚴重、受威脅人口最多的國家之一，其中滑坡是地災中最嚴重的一種災害，西南地區(qū)是我國滑坡災害最嚴重的區(qū)域，滑坡災害能阻斷河流形成堰塞湖，間接影響人類安全，對地質(zhì)災害的監(jiān)測和防治工作，是社會治理的重要課題。2018年10月10日，習近平主持召開中央財經(jīng)委員會第三次會議，部署“九項重點工程”，其中有三項就是關于地質(zhì)災害的風險調(diào)查排查、綜合治理以及災害監(jiān)測預警信息化工程。由此可見地質(zhì)災害的防治工作，已經(jīng)上升到國家戰(zhàn)略層面。

對滑坡進行監(jiān)測和早期預報，是防止和減少滑坡災害的重要方法，傳統(tǒng)滑坡監(jiān)測和預警主要采用傾斜儀、應力計、水位計、位移計等儀器進行現(xiàn)場測量，以及運用TDR、GPS 技術等進行電子遠程監(jiān)測。但這些方法卻不適宜應用到人員難以進入的區(qū)域和廣面積區(qū)域監(jiān)測上，而且造價特別高。衛(wèi)星遙感技術的發(fā)展使大面積區(qū)域監(jiān)測成為可能，同時無需人員進入，在對大面積自然區(qū)域的滑坡監(jiān)測上有著傳統(tǒng)方法無法比擬的優(yōu)勢。但是，以前的遙感調(diào)查提供的解譯成果卻是宏觀的、半定量或定性的，缺少較精確的滑坡信息。合成孔徑雷達差分干涉測量（D-InSAR）技術是近些年來在雷達干涉基礎上發(fā)展起來的一種微波遙感技術，具有光學遙感不可比擬的優(yōu)點，而且對于地表微小形變探測具有厘米級甚至毫米級探測能力[1]，使其在對地表沉降、地震以及火山等大范圍地表形變的監(jiān)測中得到迅速而且廣泛的應用。

D-InSAR 技術應用于滑坡監(jiān)測最早的時間是1995年，法國學者Achache[2]等使用從ERS-1 衛(wèi)星獲取的SAR 數(shù)據(jù)，再利用D-InSAR 技術對法國南部saint-Etienne-de-Tinee 滑坡進行了研究，對比了D-InSAR 技術監(jiān)測滑坡的形變精度與地面監(jiān)測方法監(jiān)測的精度，最后得出兩種方法的監(jiān)測精度相一致。隨后許 多學者如Fruneau[3]，Vietmeier[4]，Squarzoni[5]，Rizzo and Tesauro[6]等也陸續(xù)開始對D-InSAR 技術應用于滑坡監(jiān)測進行了研究，此后，許多研究者也分別使用包括JERS-1、ENVISAT、RADARSAT-2 在內(nèi)的多種衛(wèi)星數(shù)據(jù)，利用D-InSAR 技術進行了有關滑坡監(jiān)測的研究，也取得了不錯的效果[7-11]。由此可見利用D-InSAR 技術進行滑坡災害研究具有非常重要的意義，目前已成為國際上的一個研究熱點[12]。

1 D-InSAR 測量原理

1.1 SAR 與InSAR 的基本原理

雷達遙感影像是有源系統(tǒng)，不需要外部照明，而且是全氣候，穿透云、雨、干沙和部分植被可進行高精度形變測量。SAR（Synthetic Aperture Radar）合成孔徑雷達的基本原理是把雷達放在運動平臺上，獲得目標距離—方位二維高分辨微波圖像。SAR 系統(tǒng)分類有星載SAR、機載SAR、地基SAR，其中星載SAR 和地基SAR 適用于形變測量和變化檢測。 InSAR（Interferometry Synthetic Aperture Radar）合成孔徑雷達干涉測量技術是在SAR 技術的基礎上發(fā)展而來的，它是通過兩幅天線同時觀測或兩次近平行的觀測而獲取地面同一區(qū)域兩幅影像的方法，通過SAR 以及InSAR 技術能進行廣域監(jiān)測，發(fā)現(xiàn)新的地質(zhì)災害區(qū)域，能提供目標形變的歷史變化序列，以及災害點治理前后的變化數(shù)據(jù)與目標變化的趨勢性信息，SAR 形變觀測具備時空大數(shù)據(jù)的特征。InSAR 處理流程包括：基線估算、干涉圖生成、干涉圖去平、干涉圖自適應濾波和相干生成、相位解纏、軌道重定義、高程/形變轉(zhuǎn)換[13，14]。

1.2 D-InSAR 基本原理

差分干涉合成孔徑雷達測量（D-InSAR）技術是傳統(tǒng)InSAR 技術的擴展，該技術利用兩幅干涉圖或者一幅干涉圖結(jié)合高精度數(shù)字高程模型（DEM）來測量地表的微小形變信息[15]。D-InSAR 技術的流程包括：基線估算、干涉圖生成、自適應濾波及相干性生成、相位解纏、軌道精煉和重去平、形變圖生成等。

2 滑坡災害監(jiān)測案例

2018年10月11日，金沙江白格發(fā)生山體滑坡，金沙江的干流河道堵塞，形成了堰塞湖。危及附近地區(qū)人民群眾的生命與財產(chǎn)安全，同時又對多個水電站造成威脅。2018年11月3日，白格村原山體滑坡點發(fā)生二次滑坡。災害發(fā)生后，“地質(zhì)災害InSAR 技術研究中心”第一時間獲取高分辨率衛(wèi)星影像進行了災情遙感影像解譯分析，并利用合成孔徑雷達干涉測量技術展開了滑前形變探測等工作。

2.1 基于災前Sentinel-1 衛(wèi)星分析

利用4 期災前（2018年9月4日～10月10日，每間隔12 天一期）Sentinel-1 衛(wèi)星SAR 數(shù)據(jù)，開展D-InSAR 技術處理，發(fā)現(xiàn)該滑坡一直在持續(xù)發(fā)生形變。9月4日～9月16日，該山體斜坡中上部形變特征明顯（圖1a）；至9月28日，斜坡上部持續(xù)保持變形（圖1b）；至10月10日，整個斜坡體形變特征更為明顯，形變量級及范圍均有顯著增加，對岸亦有形變特征出現(xiàn)（圖1c），直至11日凌晨發(fā)生山體滑坡。

其中圖 a 為 Sentinel-1 第一幅干涉圖（20180904-20180916）， 圖 b 為 第 二 幅 干 涉 圖（20180916-20180928）， 圖 c 為 第 三 幅 干 涉 圖（20180928-20181010）；紅線代表InSAR 觀測得到的山體斜坡變形邊界。

圖1 金沙江白格山體滑坡災前Sentinel-1 衛(wèi)星形變監(jiān)測結(jié)果

2.2 基于災前ALOS-2 衛(wèi)星分析

利用4 期災前（2017年7月24日～2018年7月23日）ALOS-2 衛(wèi)星SAR 數(shù)據(jù)，開展D-InSAR 技術處理，發(fā)現(xiàn)該滑坡在監(jiān)測時段均具有不同程度的形變特征。其中，2017年7月至11月歷時4個月，斜坡上部形變特征明顯，中部失相干嚴重，可能是該時段內(nèi)形變量過大所致（圖2a）；2017年11月至2018年5月歷時6個月，斜坡形變特征依然明顯，形變量大（圖2b）；2018年5月至7月歷時近2個月，整個斜坡出現(xiàn)嚴重失相干，可能是斜坡體持續(xù)形變所致（圖2c）。這與Sentinel-1 衛(wèi)星得到的形變結(jié)論一致，且較好地反應了該斜坡在更早期就在發(fā)生持續(xù)形變。

其中圖 a 為 ALOS-2 第一幅干涉圖（20170724-20171127）， 圖 b 是 第 二 幅 干 涉 圖（20171127-20180528）， 圖 c 是 第 三 幅 干 涉 圖（20180528-20180723），分別表示災前不同時段內(nèi)斜坡變形特征；紅線代表InSAR 觀測得到的山體斜坡變形邊界。

圖2 金沙江白格山體滑坡災前ALOS-2 衛(wèi)星形變監(jiān)測結(jié)果

2.3 基于災中、災后COSMO-SkyMed 衛(wèi)星監(jiān)測分析

利用2018年10月14日至11月27日的COSMO-SkyMed 衛(wèi)星影像（圖3）對滑坡進行遙感解譯，發(fā)現(xiàn)滑坡體總面積達1.62 km2，其中滑源區(qū)面積0.77 km2，堰塞壩面積0.85 km2?；磪^(qū)長約1350 m，平均寬約600 m，頂部高程約3680 m，底部高程約2980m，高差達700 m（高程測自ALOS-1 DEM）；堰塞體順河長約1600 m，平均寬約490 m，最大沖高達150 m，堰塞體平均厚度估算約40 m，體積約3400萬m3（圖4）。

圖3 COSMO-SkyMed 衛(wèi)星影像數(shù)據(jù)

圖4 COSMO-SkyMed 衛(wèi)星形變監(jiān)測結(jié)果

2.4 COSMO-SkyMed 衛(wèi)星監(jiān)測與實測數(shù)據(jù)交叉驗證

從首次滑坡至12月份，持續(xù)對滑坡區(qū)域進行了持續(xù)的數(shù)據(jù)獲取、處理和分析，對滑坡體的分布、動態(tài)及變化規(guī)律等進行了動態(tài)監(jiān)測如圖5，為現(xiàn)場應急處理等提供了參考。同時對衛(wèi)星監(jiān)測值與實測數(shù)據(jù)進行比較如圖6，得出兩種方法形變趨勢基本一致，監(jiān)測到的形變量級不同，SAR 監(jiān)測形變量級要比實測數(shù)據(jù)大。SAR 監(jiān)測面狀形變信息，可監(jiān)測到人員不易到達的危險區(qū)域。地面實測得到點狀形變信息，一般布設在人員可達的滑坡體邊緣處，形變可能相對較小。

圖5 不同時期衛(wèi)星形變監(jiān)測

圖6 不同時期的衛(wèi)星監(jiān)測值與實測值的比較

3 結(jié)論

利用D-InSAR 技術對滑坡區(qū)域進行動態(tài)監(jiān)測為現(xiàn)場應急處理提供了參考，同時D-InSAR 監(jiān)測值與實測數(shù)據(jù)比較得出兩種方法形變趨勢基本一致，只是兩種方法監(jiān)測到的形變量級不同，地面實測可得到點狀形變信息，這些點一般布設在人員可達的滑坡體邊緣處，形變可能相對較小。D-InSAR 監(jiān)測形變量級要比實測數(shù)據(jù)大而且D-InSAR 可監(jiān)測面狀形變信息，能夠監(jiān)測到人員不易到達的危險區(qū)域，這為無人區(qū)滑坡檢測提供了很好的監(jiān)測路徑與方法。

4 結(jié)語

高分辨率星載SAR 監(jiān)測，可對滑坡體形變空間分布、發(fā)展動態(tài)及運動規(guī)律進行監(jiān)測，并為地面測量點位布設提供指導，為驗證滑坡體動力模型提供數(shù)據(jù)參考。金沙江白格滑坡區(qū)斜坡在發(fā)生滑坡災害前，已存在明顯的、可被遙感技術感知和發(fā)現(xiàn)的形變信息。如果能在此滑坡發(fā)生前，利用高精度遙感和D-InSAR 技術對相關區(qū)域進行地質(zhì)災害隱患“普查”，提前發(fā)現(xiàn)滑坡的變形跡象，識別出滑坡潛在隱患，對該滑坡實施監(jiān)測預警和風險評估（包括對滑坡堵河的可能性及危害性進行預測評估），便可掌握防災減災的主動權，變事后應急救災為事前主動防災。由此可見提高滑坡災害動態(tài)監(jiān)測技術的水平對于防災減災具有重要的意義。