斜坡InSAR變形監(jiān)測效果野外判定方法研究

梁京濤 趙聰 王軍

摘要：InSAR技術(shù)憑借其在地表形變監(jiān)測方面的優(yōu)勢，已逐漸廣泛應(yīng)用于地質(zhì)災(zāi)害領(lǐng)域。但由于氣候條件、茂密植被以及山區(qū)地形起伏過大等因素限制，地質(zhì)災(zāi)害InSAR變形監(jiān)測的應(yīng)用范圍及效果仍存在一定局限性。為有效評估InSAR變形監(jiān)測在某區(qū)域應(yīng)用效果優(yōu)劣，亟需建立較為系統(tǒng)的野外調(diào)查判定標(biāo)準(zhǔn)，以便于批量快捷開展InSAR變形監(jiān)測結(jié)果野外調(diào)查評估工作。從斜坡表部變形特征和孕災(zāi)地質(zhì)環(huán)境條件兩個(gè)方面總結(jié)了野外調(diào)查內(nèi)容，建立了一套斜坡InSAR變形監(jiān)測效果野外判定標(biāo)準(zhǔn)方法，并選取了滇西北地區(qū)典型實(shí)例，闡述了該套方法的具體應(yīng)用流程。研究結(jié)果表明：① 野外判定內(nèi)容包括斜坡表部變形特征和孕災(zāi)地質(zhì)環(huán)境條件兩個(gè)方面，其中斜坡表部變形特征是野外判定內(nèi)容的核心，而孕災(zāi)地質(zhì)環(huán)境條件則是地質(zhì)災(zāi)害發(fā)育的基礎(chǔ)條件。② 該套方法依據(jù)肉眼以及簡易測量等手段，實(shí)現(xiàn)了滇西北地區(qū)斜坡表部變形特征的快速識別，應(yīng)用效果較為理想，適宜于開展大批量斜坡InSAR變形監(jiān)測效果快速評估。

關(guān) 鍵 詞：斜坡地表形變; 地質(zhì)災(zāi)害; InSAR監(jiān)測; 野外調(diào)查; 滇西北地區(qū)

中圖法分類號： P642.22 ? 文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2022.01.019

0 引 言

InSAR是近年來逐步發(fā)展起來的一門對地觀測技術(shù)，最早應(yīng)用于獲取高精度的地形數(shù)據(jù)[1]。Gabriel等[2]于1989年首次發(fā)現(xiàn)利用InSAR技術(shù)可以獲取厘米級精度的形變，但是并沒有引起足夠的重視。直到Massonnet等[3]通過ERS-1數(shù)據(jù)獲取了1992年Landers地震的形變場，引起了國際地震界的震驚，此后InSAR獲取形變技術(shù)才得到了廣泛的關(guān)注。隨著InSAR技術(shù)的不斷發(fā)展，該技術(shù)已逐步引入于地面沉降、滑坡等諸多地質(zhì)災(zāi)害形變監(jiān)測中來[4-7]，現(xiàn)已成為地質(zhì)災(zāi)害早期識別領(lǐng)域主要技術(shù)手段之一。

目前，基于InSAR技術(shù)開展地質(zhì)災(zāi)害的應(yīng)用研究主要集中在不同類型雷達(dá)衛(wèi)星數(shù)據(jù)對比[8]、不同干涉處理方法對比[9-12]、典型區(qū)域應(yīng)用[13-14]等。相比而言，在滑坡InSAR變形監(jiān)測效果評價(jià)及野外判定方法方面的研究文獻(xiàn)較少。在野外調(diào)查中，由于肉眼對斜坡微變形（毫米級）觀測能力較弱甚至無法感知，此項(xiàng)工作一直很少涉及。雖然借助儀器可以現(xiàn)場測定，但是耗時(shí)耗力，效率低下且成本較高，不適用于大批量InSAR變形監(jiān)測點(diǎn)快速復(fù)核。因此，如何進(jìn)一步提高效率，方便快捷地開展InSAR變形監(jiān)測結(jié)果批量驗(yàn)證，是無法回避的技術(shù)問題，需要進(jìn)一步研究。

對于地質(zhì)災(zāi)害隱患的InSAR監(jiān)測而言，變形證據(jù)能否在野外有效獲取，直接關(guān)系到監(jiān)測效果的優(yōu)劣和該項(xiàng)技術(shù)的推廣應(yīng)用。目前，對地質(zhì)災(zāi)害隱患InSAR變形監(jiān)測結(jié)果的野外復(fù)核內(nèi)容主要集中在變形有無、變形程度大小、變形出現(xiàn)時(shí)間與InSAR監(jiān)測時(shí)間是否一致、斜坡表部變形現(xiàn)象之間關(guān)聯(lián)分析及斜坡穩(wěn)定性判定等方面。查閱國內(nèi)外文獻(xiàn)和現(xiàn)有InSAR行業(yè)規(guī)范可知[15]，該部分內(nèi)容和方法未得到明確。基于上述原因，筆者結(jié)合近年InSAR監(jiān)測效果野外工作實(shí)踐，將野外調(diào)查中存在的問題和解決方法進(jìn)行了初步總結(jié)，供同行參考。

1 野外判定內(nèi)容

InSAR變形監(jiān)測結(jié)果以矢量化圖形信息為載體，其干涉圖借助DSM可以實(shí)現(xiàn)3D表達(dá)，包含了空間信息和時(shí)間信息兩個(gè)方面。其中空間信息包括變形區(qū)位置、位移方向、變形尺度等;因數(shù)據(jù)處理方法不同，時(shí)間信息包含內(nèi)容略有差異，包括變形起止時(shí)間、持續(xù)時(shí)間、隨時(shí)間表現(xiàn)出的變形特征等。

1.1 野外判定內(nèi)容分類

根據(jù)監(jiān)測對象不同，InSAR變形監(jiān)測內(nèi)容主要包括兩部分。① 不考慮地質(zhì)災(zāi)害要素條件，僅針對斜坡變形特征與InSAR監(jiān)測結(jié)果吻合程度進(jìn)行復(fù)核，包括是否存在變形、變形程度大小和變形時(shí)間特征。② 考慮地質(zhì)災(zāi)害的要素條件，除了斜坡變形特征外，還包括斜坡臨空條件、物質(zhì)組成及威脅對象等孕災(zāi)地質(zhì)環(huán)境條件。后者是形成地質(zhì)災(zāi)害的必要條件，直接影響地質(zhì)災(zāi)害識別的有效性。

為考慮地質(zhì)災(zāi)害InSAR監(jiān)測識別的準(zhǔn)確率和有效性，野外判定內(nèi)容應(yīng)包括變形特征和孕災(zāi)地質(zhì)環(huán)境條件兩個(gè)方面。其中變形特征是野外判定內(nèi)容的核心，同時(shí)也是前提條件。一般來說，斜坡如果未見明顯變形特征，則可以不進(jìn)行孕災(zāi)地質(zhì)環(huán)境條件驗(yàn)證工作;如果只有變形特征卻不具備孕災(zāi)地質(zhì)環(huán)境條件的斜坡，則不屬于常規(guī)地質(zhì)災(zāi)害隱患范疇，變形特征沒有實(shí)際意義。

1.2 變形特征

斜坡變形特征包括空間變形信息和時(shí)間變形信息兩個(gè)方面?？臻g變形信息主要包含以下6個(gè)方面：① 斜坡表部裂縫的類型、數(shù)量、規(guī)模、組合、性質(zhì);② 斜坡表部固定建筑物如房屋墻面、院壩及地基開裂破壞情況;③ 斜坡表部線性工程如道路、水渠拉裂變形情況;④ 斜坡下坐程度、表部溜滑或局部滑動(dòng)變形等;⑤ 斜坡表部植被樹木破壞情況;⑥ 斜坡已有工程如擋墻、抗滑樁、攔石墻、防護(hù)網(wǎng)運(yùn)行及破壞情況等。時(shí)間變形信息包括斜坡變形出現(xiàn)時(shí)間、持續(xù)時(shí)間、結(jié)束時(shí)間及其與監(jiān)測周期的吻合性等。

1.3 孕災(zāi)地質(zhì)環(huán)境條件

孕災(zāi)地質(zhì)環(huán)境條件主要指形成地質(zhì)災(zāi)害的臨空條件、物質(zhì)組成及激發(fā)因素。臨空條件包括斜坡高度、坡度、走向、坡向、平面形態(tài)、剖面形態(tài)等微地貌特征;物質(zhì)組成是地質(zhì)災(zāi)害形成的載體，斜坡出露的地層巖性、產(chǎn)狀及斜坡結(jié)構(gòu)特征，是地質(zhì)災(zāi)害發(fā)生與否的內(nèi)在條件;激發(fā)因素主要指周邊人類工程及經(jīng)濟(jì)活動(dòng)類型，包括農(nóng)戶修建房屋、修建公路、水渠、土地整理、工程開礦、隧道開挖、開墾耕地等。

2 野外判定方法

2.1 野外調(diào)查驗(yàn)證方式

通過野外實(shí)地調(diào)查、目視判定、簡易設(shè)備測量和訪問等形式開展。調(diào)查通過攜帶GPS定位加圖紙位置識別的方式，基于目視觀測進(jìn)行，同時(shí)需要借助卷尺、羅盤等少量測量設(shè)備，以現(xiàn)場徒步穿插的形式沿路觀察、記錄描述，記錄方式包括相機(jī)拍攝、無人機(jī)航空攝影、現(xiàn)場文字記錄等。在尋找變形證據(jù)的同時(shí)，兼顧周邊地質(zhì)環(huán)境調(diào)查。在野外調(diào)查工作開展之前，應(yīng)盡量收集工作區(qū)已有地質(zhì)災(zāi)害資料，對于現(xiàn)場調(diào)查未見明顯變形特征的點(diǎn)，盡可能尋找當(dāng)?shù)厝罕娺M(jìn)行訪問復(fù)核，以保證調(diào)查精度。

2.2 斜坡變形程度量化及劃分類別

InSAR監(jiān)測變形區(qū)斜坡變形類型多樣、變形程度大小不一。為了評價(jià)InSAR監(jiān)測效果，便于理解和統(tǒng)計(jì)分析，本次研究在綜合考慮斜坡表部滑動(dòng)特征、裂縫發(fā)育特征、固定建筑物、線性工程變形程度、植被發(fā)育特征及已有工程破壞程度因素的基礎(chǔ)上，將斜坡變形程度分為明顯變形、輕微變形、未見明顯變形三類，并根據(jù)每一個(gè)因素具體發(fā)展程度，進(jìn)行指標(biāo)量化、綜合分類（見表1）。

2.3 變形時(shí)間獲取

斜坡變形起止時(shí)間與監(jiān)測周期吻合程度是滑坡早期識別的重要環(huán)節(jié)。在野外調(diào)查過程中，斜坡變形時(shí)間主要通過現(xiàn)場實(shí)地調(diào)查和訪問當(dāng)?shù)厝罕妰煞N方法獲取?，F(xiàn)場調(diào)查主要依據(jù)變形新鮮程度、發(fā)展趨勢、表部植被覆蓋情況等，來推斷斜坡變形時(shí)間。以拉張裂縫為例，可通過裂縫張開度、深度、植被和泥土填充情況進(jìn)行判斷，如裂縫張開度良好，具有一定深度且裂縫內(nèi)部無明顯填充物質(zhì)，可判定為新鮮裂縫，為近期發(fā)生;如裂縫明顯閉合、裂縫內(nèi)部被明顯填充，周邊植被恢復(fù)，僅留下坐臺坎，則該裂縫時(shí)間較長，為早期發(fā)生;如若進(jìn)一步判定裂縫發(fā)育的具體時(shí)間，則需要借助多期光學(xué)影像或進(jìn)一步走訪周邊當(dāng)?shù)卮迕?。一般來說，調(diào)查訪問周期要涵蓋InSAR干涉像對起止時(shí)間，同時(shí)還要適當(dāng)外延，進(jìn)一步確定變形發(fā)生時(shí)間與監(jiān)測周期吻合度。

2.4 InSAR變形監(jiān)測效果評價(jià)

InSAR變形監(jiān)測效果的優(yōu)劣，需要通過斜坡有效變形率來進(jìn)行量化。所謂有效變形指經(jīng)過野外判定，斜坡在監(jiān)測周期內(nèi)發(fā)生變形，且具備發(fā)展成為滑坡的孕災(zāi)地質(zhì)背景條件。僅發(fā)生變形而不具備孕災(zāi)地質(zhì)環(huán)境條件或變形時(shí)間不在監(jiān)測周期之內(nèi)，均為無效InSAR變形。一般而言，人類工程活動(dòng)是造成無效變形的主要影響因素，比如開挖坡腳、修建公路、平整地基、土地整理、人工堆積棄土棄渣等。

斜坡有效變形率可以用該區(qū)域InSAR監(jiān)測的有效變形斜坡數(shù)量與InSAR解譯的斜坡總數(shù)量之比來確定，如式（1） 所示。若該區(qū)域InSAR監(jiān)測的有效變形斜坡越多，則斜坡有效變形率越高，即該區(qū)域InSAR變形監(jiān)測效果越理想。R=Ne（Ne+Nv）×100%

式中：R為有效變形率;Ne為有效變形斜坡的數(shù)量;Nv為無效變形靶區(qū)的數(shù)量。

3 典型實(shí)例驗(yàn)證

3.1 典型明顯變形斜坡野外識別

以車邑坪滑坡InSAR形變監(jiān)測點(diǎn)為例，該點(diǎn)位于云南省怒江傈僳族自治州蘭坪縣石登鄉(xiāng)車邑坪村，距離黃登水電站壩址沿江公路里程約39 km，地理坐標(biāo)為N 26°18′31″，E 99°07′27″。采用Stacking-InSAR技術(shù)的Sentinel-1數(shù)據(jù)（監(jiān)測周期為20170202～20191202）干涉結(jié)果提示該區(qū)域疑似存在變形（見圖1）。根據(jù)孕災(zāi)地質(zhì)背景條件分析，該區(qū)域?qū)儆诟呱胶庸鹊孛?，斜坡坡度?0°～30°，地表第四系（Q）覆蓋層分布較廣，下伏基巖為侏羅系中統(tǒng)花開左組（J2h），巖性為砂泥巖加粉砂巖;臨江附近有區(qū)域性次級斷裂通過，具備滑坡發(fā)育的地質(zhì)環(huán)境條件。

對該點(diǎn)室內(nèi)InSAR解譯變形區(qū)展開現(xiàn)場調(diào)查，該斜坡位于瀾滄江黃登水電站庫區(qū)左側(cè)，斜坡中部發(fā)育兩條較大規(guī)模沖溝（1號沖溝和2號沖溝），斜坡前緣坡腳位于黃登水電站庫區(qū)蓄水位線以下，兩條沖溝之間坡體為主變形區(qū)，斜坡體表部公路下部發(fā)育多條拉張裂縫，多為2019年6月水庫蓄水位變化過程中產(chǎn)生，裂縫主要位于斜坡中前部，以拉張裂縫為主，裂縫的產(chǎn)生多伴隨斜坡地面沉降、淺表層溜滑等，裂縫形態(tài)多呈弧形或圈椅狀展布。斜坡中部公路兩側(cè)為車邑坪村農(nóng)戶聚集區(qū)，多處房屋墻體及院壩可見拉裂破壞，裂縫寬度2～10 cm不等（見圖2）。

綜合分析可知，該斜坡表部多條裂縫成規(guī)模發(fā)育，坡體表部房屋多處拉裂變形，且斜坡發(fā)生形變時(shí)間正好處于InSAR監(jiān)測周期內(nèi)。根據(jù)表1斜坡變形程度量化及分類，該斜坡變形程度屬于明顯變形。該滑坡屬于水庫蓄水后形成的外動(dòng)力地質(zhì)作用引發(fā)的古滑坡復(fù)活，滑坡活動(dòng)性高，InSAR變形區(qū)解譯正確。該滑坡現(xiàn)處于基本穩(wěn)定狀態(tài)，在暴雨、水位驟降、地震等不利工況下發(fā)生進(jìn)一步變形失穩(wěn)可能性大。

3.2 典型輕微變形斜坡野外識別

以獨(dú)家村三組滑坡InSAR形變監(jiān)測點(diǎn)為例，該點(diǎn)位于云南省怒江傈僳族自治州蘭坪縣中排鄉(xiāng)獨(dú)家村三組，地理坐標(biāo)為N 26° 51′26″，E 99°10′25″。采用Stacking-InSAR技術(shù)的Sentinel-1數(shù)據(jù)（監(jiān)測周期：20170202～20191202）干涉結(jié)果提示該區(qū)域疑似存在變形（見圖3）。該區(qū)域孕災(zāi)地質(zhì)背景條件與車邑坪滑坡相似，屬于高山河谷地貌，斜坡平均坡度25°，地表第四系（Q）覆蓋層分布較廣，下伏基巖為侏羅系中統(tǒng)花開左組（J2h），巖性為砂泥巖加粉砂巖，同樣具備滑坡發(fā)育的地質(zhì)環(huán)境條件。

進(jìn)一步對該點(diǎn)進(jìn)行了現(xiàn)場調(diào)查。該斜坡位于瀾滄江黃登水電站庫區(qū)右側(cè)，整體平面形態(tài)呈圈椅狀，斜坡坡體后緣略微下凹，中前部外凸擠壓明顯，剖面形態(tài)呈階梯型，具備較為典型的滑坡地貌特征。該斜坡坡腳至坡體上部盤山公路水泥路面普遍開裂，最大拉伸寬度1～2 cm，長約20～100 m不等，局部路段路面橫向拉斷（水泥路面于2020年4月重新補(bǔ)修，原裂縫于2020年9月現(xiàn)場調(diào)查期間未見進(jìn)一步開裂現(xiàn)象）。此外，斜坡體中后部房屋普遍存在小規(guī)模輕微開裂跡象，裂縫寬度多小于0.5 cm，長度5～20 m不等（見圖4）。訪問當(dāng)?shù)剞r(nóng)戶可知，坡體表部公路及房屋裂縫集中出現(xiàn)于2018年黃登水電站庫區(qū)正式蓄水后，推測該處斜坡的變形與庫區(qū)水位變化緊密相關(guān)。

綜合分析可知，該斜坡表部公路路面出現(xiàn)多條小規(guī)模裂縫，居民房屋輕微開裂，且斜坡發(fā)生形變時(shí)間正好處于InSAR監(jiān)測周期內(nèi)。根據(jù)表1斜坡變形程度量化及分類，該斜坡變形程度屬于輕微變形。該處斜坡屬于庫區(qū)蓄水位變化引發(fā)的滑坡，滑坡活動(dòng)性中等，現(xiàn)處于基本穩(wěn)定狀態(tài)，威脅黃登水電站庫區(qū)以及滑坡體表部的農(nóng)戶、鄉(xiāng)村公路，潛在風(fēng)險(xiǎn)較高。

3.3 典型未見明顯變形斜坡野外識別

以福坪村斜坡InSAR形變監(jiān)測點(diǎn)為例，該點(diǎn)位于云南省怒江傈僳族自治州蘭坪縣金頂鎮(zhèn)福坪村，地理坐標(biāo)為N 26°27′16″，E 99°26′58″。采用Stacking-InSAR技術(shù)的Sentinel-1數(shù)據(jù)（監(jiān)測周期：20170202～20191202）干涉結(jié)果提示該區(qū)域疑似存在變形（見圖5）。根據(jù)孕災(zāi)地質(zhì)背景條件分析，該區(qū)域?qū)儆谥猩降孛?，斜坡平均坡度約15°，斜坡兩側(cè)發(fā)育沖溝，坡體表部廣泛分布第四系（Q）覆蓋層，下伏基巖為三疊系上統(tǒng)三合洞組（T3s），巖性為灰?guī)r夾泥灰?guī)r、砂巖，斜坡前緣處有區(qū)域性次級斷裂通過，斷裂走向近南北向，該點(diǎn)具備滑坡發(fā)育的地質(zhì)環(huán)境條件。

進(jìn)一步對該點(diǎn)進(jìn)行現(xiàn)場調(diào)查，該斜坡位于蘭坪縣城東側(cè)，距離約2 300 m，整體平面形態(tài)呈舌形，斜坡坡體后緣略微下凹，坡體中部為福坪村農(nóng)戶聚集區(qū)，前部及后部可見大量光伏發(fā)電板分布。根據(jù)現(xiàn)場調(diào)查，該斜坡區(qū)域內(nèi)水泥公路、房屋、空心磚墻體以及田間小路均未見明顯裂縫以及局部溜滑、垮塌等變形跡象發(fā)育，訪問當(dāng)?shù)卮迕竦弥?，? a來斜坡區(qū)域內(nèi)均未見明顯變形現(xiàn)象（見圖6）。該斜坡現(xiàn)整體處于穩(wěn)定狀態(tài)，由于該電站于2019年12月建成完工，大規(guī)模人類工程活動(dòng)時(shí)間正好處于InSAR監(jiān)測周期內(nèi)，因此，InSAR形變解譯異常推測與修建光伏發(fā)電站的人類工程活動(dòng)相關(guān)。

綜合分析可知，該斜坡表部未見局部滑動(dòng)跡象，建筑物、公路等均完好，未見明顯裂縫發(fā)育，根據(jù)表1斜坡變形程度量化及分類，該斜坡變形程度屬于未見明顯變形。

4 討 論

（1） InSAR技術(shù)已逐漸廣泛應(yīng)用于地質(zhì)災(zāi)害隱患識別調(diào)查領(lǐng)域，國家、省部級等相關(guān)部門也大力推行各種不同方式類型、不同精度的地質(zhì)災(zāi)害隱患InSAR監(jiān)測識別項(xiàng)目，取得了較為理想的應(yīng)用效果。同時(shí)也認(rèn)識到由于受氣候條件、茂密植被以及山區(qū)地形起伏過大等因素限制，地質(zhì)災(zāi)害InSAR變形監(jiān)測應(yīng)用范圍及效果仍存在一定局限性[16]。評估某種技術(shù)手段或處理方法在地質(zhì)災(zāi)害隱患識別調(diào)查中應(yīng)用效果的優(yōu)劣，需要結(jié)合野外調(diào)查工作來進(jìn)行驗(yàn)證。此外，由于InSAR處理過程中不可避免地會存在一些誤差，InSAR監(jiān)測變形數(shù)據(jù)不能等同于斜坡變形的直接證據(jù)，同樣需要野外調(diào)查工作進(jìn)行輔助評估。因此，對于地質(zhì)災(zāi)害隱患識別而言，基于“變形靶區(qū)”的微變形InSAR監(jiān)測和斜坡宏觀變形特征野外調(diào)查獲取相結(jié)合，或許能得到更好的識別效果。

（2） 本文所論述的斜坡InSAR變形監(jiān)測效果野外驗(yàn)證方法，主要是基于簡易測量手段以及肉眼對斜坡變形特征可以識別判定的前提下開展的。在野外調(diào)查驗(yàn)證工作中，必然也存在著少量簡易測量手段和肉眼無法識別的斜坡微變形，一定程度上影響了驗(yàn)證精度和監(jiān)測效果評估。

5 結(jié) 論

（1） InSAR變形監(jiān)測信息包含了空間信息和時(shí)間信息兩個(gè)方面，對室內(nèi)監(jiān)測結(jié)果的野外復(fù)核應(yīng)圍繞這兩個(gè)方面信息展開?？紤]到地質(zhì)災(zāi)害InSAR監(jiān)測的準(zhǔn)確率和早期識別效果，野外判定內(nèi)容應(yīng)包括斜坡表部變形特征和孕災(zāi)地質(zhì)環(huán)境條件兩個(gè)方面，其中斜坡表部變形特征是野外判定內(nèi)容的核心，孕災(zāi)地質(zhì)環(huán)境條件是地質(zhì)災(zāi)害早期識別的必要條件。

（2） 斜坡InSAR變形程度可以分為明顯變形、輕微變形、未見明顯變形三類，劃分標(biāo)準(zhǔn)可在參照斜坡表部滑動(dòng)特征、裂縫發(fā)育特征、固定建筑物、線性工程變形特征、表部植被發(fā)育特征及已有治理工程破壞程度因素的基礎(chǔ)上，量化指標(biāo)，定性分類。

（3） InSAR變形監(jiān)測效果的優(yōu)劣，可通過斜坡有效變形率來進(jìn)行量化。斜坡在監(jiān)測周期內(nèi)發(fā)生變形，且具備滑坡的孕災(zāi)地質(zhì)背景條件的為有效變形;僅發(fā)生變形而不具備孕災(zāi)地質(zhì)環(huán)境條件或變形時(shí)間不在監(jiān)測周期之內(nèi)，均為無效InSAR變形。

參考文獻(xiàn)：

[1] GRAHAM L C.Synthetic interferometer radar for topographic mapping[J].Proceedings of the IEEE，1974，62（6）：763-768.

[2] ANDREW K，RICHARD M，HOWARD A.Mapping small elevation changes over large areas：differential radar interferometry[J].Journal of Geophysical Research.Part B：Solid Earth，1989，94：9183-9191.

[3] MASSONNET D，ROSSI M，CARMONA C，et al.The displacement field of the Landers earthquake mapped by radar interferometry[J].Nature，1993，364（6433）：138-142.

[4] 張潔，胡光道，羅寧波.InSAR技術(shù)在滑坡監(jiān)測中的應(yīng)用研究[J].工程地球物理學(xué)報(bào)，2004，1（2）：148-153.

[5] 葛大慶，王艷，范景輝，等.地表形變D-InSAR監(jiān)測方法及關(guān)鍵問題分析[J].國土資源遙感，2007，19（4）：14-22.

[6] 葛大慶，戴可人，郭兆成，等.重大地質(zhì)災(zāi)害隱患早期識別中綜合遙感應(yīng)用的思考與建議[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2019，44（7）：949-956.

[7] 廖明生，唐婧，王騰，等.高分辨率SAR數(shù)據(jù)在三峽庫區(qū)滑坡監(jiān)測中的應(yīng)用[J].中國科學(xué)（地球科學(xué)），2012，42（2）：217-229.

[8] 康亞，趙超英，張勤，等.InSAR滑坡探測技術(shù)研究：以金沙江烏東德水電站段為例[J].大地測量與地球動(dòng)力學(xué)，2018，13（10）：1053-1057.

[9] BATUHAN O，F(xiàn)ILIZ S，SHIMON W，et al.Time series analysis of InSAR data：methods and trends[J].ISPRS Journal of Photogrammetry and Remote Sensing，2016（115）：90-102.

[10] 程滔，葛春青，陶舒，等.小基線集合成孔徑雷達(dá)干涉測量算法及其應(yīng)用[J].測繪科學(xué)，2015，40（11）：96-99.

[11] 林昊.基于D-InSAR和Offset Tracking技術(shù)的滑坡形變場提取研究[D].北京：中國地質(zhì)大學(xué)（北京），2014.

[12] 王之棟，文學(xué)虎，唐偉，等.聯(lián)合多InSAR技術(shù)的龍門山：大渡河區(qū)域地災(zāi)隱患早期探測[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2020，45（3）：451-459.

[13] 陸會燕，李為樂，許強(qiáng)，等.光學(xué)遙感和InSAR結(jié)合的金沙江白格滑坡上下游滑坡隱患早期識別[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2019，44（9）：1342-1354.

[14] 黃潔慧，謝謨文，王立偉.基于SBAS-InSAR技術(shù)的白格滑坡形變監(jiān)測研究[J].人民長江，2019，50（12）：101-105.

[15] 中國地質(zhì)災(zāi)害防治工程行業(yè)協(xié)會.地質(zhì)災(zāi)害InSAR監(jiān)測技術(shù)指南（試行）：T/CAGHP 013-2018[S].武漢：中國地質(zhì)大學(xué)出版社，2018.

[16] 許強(qiáng).對地質(zhì)災(zāi)害隱患早期識別相關(guān)問題的認(rèn)識與思考[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2020，45（11）：1651-1659.

（編輯：劉 媛）