基于多源遙感技術(shù)的紅層滑坡識(shí)別與監(jiān)測(cè)研究

王 慶 芳，鄭 志 軍，董 繼 紅，余 天 彬，劉 文，黃 細(xì) 超

(1.成都市地質(zhì)環(huán)境監(jiān)測(cè)站，四川 成都 610041； 2.自然資源部成都地質(zhì)災(zāi)害野外科學(xué)觀測(cè)研究站，四川 成都 610041； 3.四川省地質(zhì)調(diào)查院 稀有稀土戰(zhàn)略資源評(píng)價(jià)與利用四川省重點(diǎn)實(shí)驗(yàn)室，四川 成都 610081；4.四川省智慧地質(zhì)大數(shù)據(jù)有限公司，四川 成都 610081)

0 引 言

滑坡作為一種頻繁發(fā)生，破壞最為嚴(yán)重的自然災(zāi)害之一，具有隱蔽性強(qiáng)、破壞性大、突發(fā)性高等特點(diǎn)，嚴(yán)重危害著人民群眾的生命財(cái)產(chǎn)安全，同時(shí)也制約著區(qū)域經(jīng)濟(jì)的發(fā)展[1-3]。紅層是一種外觀以紅色為主色調(diào)的陸相碎屑巖沉積地層，龍泉山是典型的低山丘陵紅層區(qū)。地層主要是侏羅系、白堊系，巖性以紅色、紫紅色、磚紅色砂巖、粉砂巖、泥巖、礫巖為主[4-5]，該區(qū)域的植被覆蓋率較高，人類工程活動(dòng)強(qiáng)烈。地層具有強(qiáng)度低、親水性強(qiáng)、易崩解、軟化、膨脹等特性，易發(fā)生規(guī)模較大的順層巖性滑坡[6-8]。據(jù)調(diào)查顯示，自2008年以來(lái)，該區(qū)域發(fā)生多起大規(guī)模的滑坡災(zāi)害，如雷打石滑坡、五家墳滑坡、油榨房滑坡等。區(qū)內(nèi)滑坡具有分布廣、隱蔽性強(qiáng)、突發(fā)性高、識(shí)別難度大等特點(diǎn)。因此很有必要利用遙感技術(shù)對(duì)該區(qū)域的滑坡隱患開展大范圍的早期識(shí)別研究，對(duì)典型的滑坡開展時(shí)序監(jiān)測(cè)，分析變形規(guī)律及特征。

光學(xué)遙感技術(shù)因結(jié)果準(zhǔn)確率高、解譯精度高、獲取信息豐富、可以回溯歷史形變等特點(diǎn)廣泛用于滑坡識(shí)別、解譯分析等工作[9-10]，但是該技術(shù)自動(dòng)化程度低、工作量大。而合成孔徑雷達(dá)干涉測(cè)量(Interferometric Synthetic Aperture Radar，InSAR)技術(shù)因其非接觸、不受云霧天氣影響、空間覆蓋范圍廣、監(jiān)測(cè)精度高等優(yōu)勢(shì)，廣泛用于地表形變監(jiān)測(cè)[11-13]。張路等[14]以四川省丹巴縣為例，利用InSAR技術(shù)獲取的形變速率圖及差分干涉圖成功識(shí)別出17處滑坡隱患，同時(shí)總結(jié)了影響時(shí)序InSAR技術(shù)監(jiān)測(cè)效果的主要因素。戴可人等[15]通過SBAS技術(shù)獲取的形變速率圖進(jìn)行隱患的識(shí)別圈定，并討論了研究區(qū)域內(nèi)Sentinel-1數(shù)據(jù)可視性分布情況。朱賽楠等[16]利用光學(xué)遙感技術(shù)、InSAR技術(shù)結(jié)合物探及地表位移監(jiān)測(cè)等多個(gè)手段對(duì)金沙江流域色拉滑坡的特征、變形過程、形變機(jī)理進(jìn)行研究，并對(duì)其變形趨勢(shì)進(jìn)行了預(yù)判。李媛茜等[17]基于InSAR技術(shù)結(jié)合野外資料成功圈定了白龍江流域114處滑坡隱患，分析了滑坡發(fā)育特征，并通過對(duì)牙豁口滑坡結(jié)合無(wú)人機(jī)數(shù)據(jù)進(jìn)行分區(qū)研究，證實(shí)了InSAR技術(shù)在滑坡變形監(jiān)測(cè)的有效性和可靠性。同時(shí)一些學(xué)者也利用Stacking InSAR技術(shù)及SBAS InSAR技術(shù)進(jìn)行滑坡隱患的大范圍識(shí)別及監(jiān)測(cè)，取得了不錯(cuò)的結(jié)果[18-19]。

基于此，本文以四川省成都市龍泉山區(qū)域?yàn)閷?shí)驗(yàn)區(qū)域，收集了覆蓋該區(qū)域的ALOS PALSAR-2數(shù)據(jù)、Sentinel-1數(shù)據(jù)及高分光學(xué)數(shù)據(jù)，開展多源遙感技術(shù)在低山丘陵區(qū)紅層滑坡隱患的早期識(shí)別，并對(duì)比不同波段SAR數(shù)據(jù)的隱患識(shí)別情況；基于SBAS技術(shù)結(jié)合無(wú)人機(jī)數(shù)據(jù)及現(xiàn)場(chǎng)調(diào)查結(jié)果分析了典型滑坡的變形規(guī)律及特征，并給出了相應(yīng)建議。

1 研究區(qū)概況及數(shù)據(jù)來(lái)源

1.1 研究區(qū)概況

龍泉山位于四川省成都市東側(cè)，西接成都平原區(qū)，東部接川東丘陵區(qū)，地勢(shì)中部較高沿龍泉山山脊展布，東西兩側(cè)低，本次研究區(qū)主要分布在龍泉山中段(見圖1)，研究區(qū)地貌類型主要為構(gòu)造侵蝕低山地貌，屬于川東低山丘陵區(qū)。區(qū)內(nèi)屬中亞熱帶濕潤(rùn)氣候區(qū)，四季分明、濕潤(rùn)溫暖，降雨豐沛，雨季多集中于6～9月，降雨量約占全年的75%。調(diào)查分析發(fā)現(xiàn)降雨條件對(duì)區(qū)域的滑坡災(zāi)害具有較強(qiáng)的控制作用，滑坡發(fā)生的時(shí)間主要集中在強(qiáng)降雨期。

圖1 研究區(qū)位置Fig.1 Geographical location of the study area

1.2 數(shù)據(jù)收集與處理

本次研究采用多源遙感技術(shù)綜合識(shí)別，利用多源SAR數(shù)據(jù)和高分光學(xué)衛(wèi)星影像開展區(qū)域滑坡隱患早期識(shí)別。收集了覆蓋該區(qū)域2017～2020年的ALOS PALSAR-2降軌數(shù)據(jù)。同時(shí)也收集了可以免費(fèi)獲取的Sentinel-1數(shù)據(jù)進(jìn)行研究區(qū)內(nèi)典型滑坡隱患的形變規(guī)律監(jiān)測(cè)分析(見表1)。利用AW3D30 DSM數(shù)據(jù)來(lái)消除InSAR干涉處理中的地形相位及輔助SAR影像進(jìn)行地理編碼。Sentinel-1衛(wèi)星的POD精密軌道星歷(POD Precise Orbit Ephemerides)數(shù)據(jù)被用來(lái)輔助Sentinel-1數(shù)據(jù)的預(yù)處理和基線誤差改正。

表1 選用的遙感數(shù)據(jù)信息

同時(shí)在數(shù)據(jù)處理中對(duì)與高程不相干的隨機(jī)大氣擾動(dòng)誤差，利用李振洪教授團(tuán)隊(duì)研制的通用型衛(wèi)星雷達(dá)在線大氣改正系統(tǒng)(Generic Atmospheric Correction Online Service for InSAR，GACOS)來(lái)去除[20]。

光學(xué)衛(wèi)星數(shù)據(jù)主要包括國(guó)產(chǎn)高分二號(hào)和北京二號(hào)數(shù)據(jù)，用來(lái)進(jìn)行滑坡隱患的光學(xué)遙感解譯和識(shí)別驗(yàn)證。

2 低山丘陵區(qū)滑坡識(shí)別監(jiān)測(cè)方法

結(jié)合工作區(qū)滑坡發(fā)生的特殊性質(zhì)，如分布廣泛、隱蔽性強(qiáng)、空間尺度較小、受降雨影響大以及突發(fā)性強(qiáng)等特點(diǎn)，在基于InSAR技術(shù)用于滑坡隱患識(shí)別的相關(guān)研究和前期數(shù)據(jù)處理結(jié)果分析基礎(chǔ)上，本文結(jié)合收集到的SAR數(shù)據(jù)及其他相關(guān)資料，針對(duì)低山紅層區(qū)突發(fā)性滑坡隱患識(shí)別、監(jiān)測(cè)研究，給出了本文的滑坡隱患識(shí)別與監(jiān)測(cè)處理流程(見圖2)。該流程的關(guān)鍵是對(duì)基于地表形變速率解譯的滑坡隱患從光學(xué)影像中尋找依據(jù)，通過SBAS技術(shù)判定其變形量及變形規(guī)律，通過野外調(diào)查、核查識(shí)別結(jié)果的準(zhǔn)確性，并對(duì)災(zāi)害庫(kù)進(jìn)行更新，同時(shí)積累該區(qū)域的解譯標(biāo)識(shí)庫(kù)。

圖2 滑坡隱患識(shí)別與監(jiān)測(cè)流程Fig.2 Flow chart of landslide hidden danger identification and monitoring

2.1 Stacking InSAR技術(shù)用于滑坡隱患編目

目前開展大范圍滑坡隱患早期識(shí)別工作，普遍采用Stacking InSAR技術(shù)，該技術(shù)具有使用影像數(shù)量少、結(jié)果獲取效率高等優(yōu)勢(shì)。因此本文在對(duì)比不同波段數(shù)據(jù)獲取結(jié)果解譯情況與已有滑坡隱患庫(kù)基礎(chǔ)上，建議在進(jìn)行Stacking計(jì)算中采用ALOS-2數(shù)據(jù)與Sentinel-1數(shù)據(jù)相結(jié)合的形式進(jìn)行。

其中Stacking技術(shù)[21]是由Sandwell等在1998年提出，該技術(shù)原理是對(duì)解纏相位進(jìn)行加權(quán)平均解算，以達(dá)到削弱空間上不相關(guān)的噪聲的影響，從而獲取研究區(qū)域的平均形變速率[22-23]。主要處理步驟如圖3所示。

圖3 STACKING-INSAR技術(shù)流程Fig.3 Stacking-InSAR technology process

2.2 基于光學(xué)遙感技術(shù)的滑坡隱患識(shí)別

結(jié)合工作區(qū)地質(zhì)災(zāi)害特征，對(duì)滑坡隱患點(diǎn)的識(shí)別，可采用對(duì)比分析方法和直接識(shí)別方法進(jìn)行判定。對(duì)比分析方法首先利用遙感影像對(duì)斜坡失穩(wěn)前地質(zhì)環(huán)境條件和變形跡象進(jìn)行解譯，然后對(duì)工作區(qū)內(nèi)地質(zhì)環(huán)境條件相似或者變形跡象明顯的地質(zhì)災(zāi)害進(jìn)行遙感解譯。直接識(shí)別方法的關(guān)鍵是建立遙感影像識(shí)別標(biāo)志，遙感解譯標(biāo)志包括了地質(zhì)災(zāi)害發(fā)育的地質(zhì)環(huán)境條件和斜坡變形特征解譯兩個(gè)方面，如表2所列。

表2 滑坡光學(xué)遙感解譯識(shí)別標(biāo)志

2.3 SBAS InSAR技術(shù)用于滑坡監(jiān)測(cè)分析

對(duì)于典型滑坡隱患采用SBAS InSAR技術(shù)進(jìn)行長(zhǎng)時(shí)間序列監(jiān)測(cè)，獲取其變形趨勢(shì)及變形量等信息，因?yàn)镾entinel-1數(shù)據(jù)時(shí)間分辨率高(時(shí)間間隔為6/12 d)、免費(fèi)獲取，所以本文利用Sentinel-1數(shù)據(jù)基于SBAS技術(shù)開展典型滑坡隱患的時(shí)間序列監(jiān)測(cè)分析。

SBAS技術(shù)是由Berardino等于2002年提出的，基本原理是假設(shè)在監(jiān)測(cè)時(shí)間t0～ts內(nèi)獲取了S+1幅SAR影像，根據(jù)干涉組合條件可形成M幅干涉組合，符合式(1)：

(1)

一般根據(jù)實(shí)際情況設(shè)置一定的垂直基線閾值和時(shí)間基線閾值對(duì)M的值域進(jìn)行限制，通過這種形式可以降低因垂直基線和時(shí)間基線過大引起的失相干現(xiàn)象[23-24]。利用差分干涉計(jì)算獲得差分干涉相位，用于第j幅干涉相位是根據(jù)tA和tB(tA

3 結(jié)果分析

3.1 基于Stacking InSAR技術(shù)的滑坡隱患早期識(shí)別與效果分析

分別利用ALOS-2數(shù)據(jù)和Sentinel-1數(shù)據(jù)進(jìn)行Stacking InSAR處理。對(duì)ALOS-2數(shù)據(jù)進(jìn)行距離向和方位向采用2×4的多視處理來(lái)提高信噪比，由于ALOS-2數(shù)據(jù)量較少，因此采用全組合的形式進(jìn)行差分干涉處理，并進(jìn)行Stacking InSAR解算獲得沿雷達(dá)衛(wèi)星視線向(Line of Sight，LOS)年平均形變速率結(jié)果(見圖4)。對(duì)于Sentinel-1數(shù)據(jù)，在處理中采用8×2進(jìn)行距離向和方位向多視處理，設(shè)置空間基線不大于200 m，時(shí)間基線不大于48 d進(jìn)行干涉對(duì)組合，進(jìn)而提高相干性和運(yùn)算效率，并通過Stacking InSAR技術(shù)獲取了形變速率結(jié)果。

圖4 利用ALOS-2數(shù)據(jù)獲取沿LOS向地表形變速率Fig.4 Surface deformation rate along the LOS direction obtained by ALOS-2 data

根據(jù)形變速率圖中的條紋形態(tài)判識(shí)疑似的地災(zāi)隱患[29]，結(jié)合地形條件及光學(xué)影像，對(duì)基于不同數(shù)據(jù)源獲取的地表形變監(jiān)測(cè)結(jié)果進(jìn)行滑坡隱患的識(shí)別解譯，其中ALOS-2數(shù)據(jù)共識(shí)別出9處隱患點(diǎn)(見圖4)，野外核查發(fā)現(xiàn)8處有明顯變形，1處無(wú)明顯變形。Sentinel-1升軌數(shù)據(jù)識(shí)別出8處，經(jīng)野外驗(yàn)證5處為滑坡隱患，Sentinel-1降軌數(shù)據(jù)識(shí)別出6處，經(jīng)野外驗(yàn)證2處為滑坡隱患，Sentinel-1數(shù)據(jù)經(jīng)驗(yàn)證正確的隱患點(diǎn)與ALOS-2數(shù)據(jù)識(shí)別的隱患點(diǎn)位置一致，具體滑坡隱患信息如表3所列。

表3 滑坡隱患信息

選擇了3處經(jīng)野外驗(yàn)證為滑坡隱患進(jìn)行多源遙感數(shù)據(jù)識(shí)別對(duì)比分析(見圖5)。其中獅子山滑坡和宋家灣滑坡為新增隱患點(diǎn)。通過基于ALOS-2數(shù)據(jù)獲取了這3處滑坡隱患的解譯信息，可以看到明顯的變形信息。相應(yīng)地，在Sentinel-1升軌數(shù)據(jù)中可以看到較為明顯的形變信息。其中獅子山滑坡和鄢家溝滑坡的變形信息與ALOS-2數(shù)據(jù)獲取的變形信息具有較好的一致性，只是Sentinel-1是升軌數(shù)據(jù)，ALOS-2是降軌數(shù)據(jù)，所以在速率圖中圖斑顏色正好相反。這3處滑坡隱患在Sentinel-1降軌數(shù)據(jù)中沒有變形信息，故沒有展示。

宋家灣滑坡隱患是基于ALOS-2數(shù)據(jù)獲取的新增隱患點(diǎn)，該滑坡在Sentinel-1升軌數(shù)據(jù)中變形信息并不明顯，為此依據(jù)圖2流程進(jìn)行驗(yàn)證分析。圖6為該滑坡隱患對(duì)應(yīng)的光學(xué)影像，從圖6中可以看到該滑坡隱患具備地質(zhì)災(zāi)害發(fā)育的地質(zhì)環(huán)境條件，但是斜坡變形特征并不明顯。為此，提取基于InSAR監(jiān)測(cè)變形強(qiáng)烈的部分，再基于SBAS InSAR技術(shù)提取其時(shí)間序列變形信息(見圖7)。從圖7可以看到：該滑坡一直處于變形狀態(tài)，變形具有明顯的周期性，每年7月前后處于加速變形狀態(tài)，變形時(shí)間為工作區(qū)雨季，同時(shí)2018年1月至2020年12月累積變形量達(dá)到了132 mm，佐證了該點(diǎn)為滑坡隱患點(diǎn)，后期通過野外調(diào)查驗(yàn)證，確定該點(diǎn)為一新增滑坡隱患點(diǎn)。

圖6 宋家灣滑坡光學(xué)衛(wèi)星影像解譯信息(北京2號(hào))Fig.6 Optical satellite image interpretation information of Songjiawan landslide(satellite of Beijing 2)

圖7 宋家灣滑坡時(shí)序點(diǎn)變形信息Fig.7 Deformation information of Songjiawan landslide sequence point

3.2 典型滑坡隱患監(jiān)測(cè)

選擇利用時(shí)序InSAR技術(shù)識(shí)別的一處新增隱患點(diǎn)——鄭家房子滑坡為例進(jìn)行隱患監(jiān)測(cè)分析。該滑坡地理位置為104°13′21.93301″E，30°23′0.09891″N；滑坡形態(tài)整體呈簸箕形，滑坡前緣高程685 m，后緣高程770 m，高差85 m，縱長(zhǎng)約260 m，橫寬約150 m，坡向約245°，坡度15°～20°，滑坡平面面積為6.25萬(wàn) m2，按10 m厚度估算，估算方量為6.25×105m3，為一大型緩傾滑坡。地層巖性為侏羅系中統(tǒng)沙溪廟組(J2s)淺紫紅、黃灰色厚-塊狀中細(xì)粒巖屑砂巖、巖屑長(zhǎng)石砂巖與紫紅色、灰紫色粉砂質(zhì)泥巖、粉砂巖互層，滑坡結(jié)構(gòu)為一順向坡，屬于典型的低山紅層區(qū)順層滑坡。

同時(shí)獲取了該區(qū)域的無(wú)人機(jī)正射影像圖(見圖8)，從光學(xué)影像圖中可以看到滑坡后緣呈圈椅狀，可見早期滑動(dòng)形成的洼地緩坡平臺(tái)；滑坡兩側(cè)以負(fù)地形凹地為界，可見溝道微地貌；滑坡前緣以緩坡平臺(tái)為界，堆積體外凸擠向臨空面，陡緩交界線明顯。

圖8 鄭家房子滑坡光學(xué)遙感解譯(據(jù)2021年無(wú)人航測(cè)數(shù)據(jù))Fig.8 Interpretation of Zhengjiafangzi landslide by optical remote sensing(according to unmanned aerial survey data in 2021)

通過對(duì)此滑坡進(jìn)行現(xiàn)場(chǎng)核查，發(fā)現(xiàn)滑坡變形區(qū)域與監(jiān)測(cè)識(shí)別吻合性較好。滑坡變形區(qū)主要集中在左側(cè)中部居民區(qū)處，房屋局部發(fā)生拉裂變形，坡體局部發(fā)生溜滑。

據(jù)訪問坡體上部右側(cè)邊界處的一戶居民，可知該區(qū)域斜坡在2018年至2021年底，每年降雨期坡表均有不同程度的溜滑，該居民房屋無(wú)變形，屋后堡坎局部發(fā)生溜滑。現(xiàn)場(chǎng)調(diào)查坡體中部發(fā)育裂縫，裂縫走向NE35°，裂縫延伸長(zhǎng)度10～25 m，寬10～20 cm，錯(cuò)落高度5～15 cm。該滑坡為一土質(zhì)滑坡，斜坡坡向與滑動(dòng)方向基本一致，屬于龍泉山典型的易滑斜坡區(qū)。

為了對(duì)鄭家房子滑坡的變形特征進(jìn)行分析，采用SBAS InSAR技術(shù)對(duì)Sentinel-1升軌數(shù)據(jù)進(jìn)行處理，獲得了該滑坡沿LOS向地表形變速率結(jié)果(見圖9)，正值表示滑坡位移靠近衛(wèi)星運(yùn)動(dòng)方向變形，負(fù)值表示遠(yuǎn)離衛(wèi)星飛行方向。

圖9 鄭家房子滑坡形變速率Fig.9 Deformation rate of Zhengjiafangzi landslide

將圖9與ALOS-2數(shù)據(jù)結(jié)果對(duì)比后發(fā)現(xiàn)，兩者在強(qiáng)變形區(qū)域具有較好的一致性，圖9中形變速率最大達(dá)到了-40 mm/a。在滑坡體上選擇特征點(diǎn)提取其在時(shí)間域的累積形變量結(jié)果(見圖10)。可以看到特征點(diǎn)的變形信息具有明顯的規(guī)律性，每年的5～7月之間發(fā)生變形，變形速率最大，其他時(shí)間段基本趨于穩(wěn)定，其變形時(shí)間與當(dāng)?shù)氐挠昙疚呛稀Ｆ渲蠦、C點(diǎn)形變速率最大，累積變形量也最大，其在2018年1月至2020年12月累積變形量達(dá)到134 mm。結(jié)合現(xiàn)場(chǎng)調(diào)查情況可知該滑坡變形主要集中在每年的汛期，在強(qiáng)降雨作用下，容易發(fā)生大規(guī)模滑坡，威脅坡體居民區(qū)、村道，因此需要在雨季加強(qiáng)持續(xù)性監(jiān)測(cè)。

圖10 鄭家房子滑坡體監(jiān)測(cè)點(diǎn)累計(jì)形變曲線Fig.10 Cumulative deformation curve of monitoring points of Zhengjiafangzi landslide mass

4 討 論

4.1 低山丘陵區(qū)雷達(dá)數(shù)據(jù)適用性及SAR數(shù)據(jù)選取分析

3.1節(jié)中提到Sentinel-1數(shù)據(jù)和ALOS-2數(shù)據(jù)識(shí)別滑坡隱患的數(shù)目及準(zhǔn)確率不一致，經(jīng)分析主要是受差分干涉相干性的影響。ALOS-2數(shù)據(jù)搭載L波段的傳感器，波長(zhǎng)較長(zhǎng)，對(duì)地表植被具有較好的穿透性。Sentinel-1數(shù)據(jù)搭載C波段，使得在植被茂密區(qū)域相干性差，尤其是夏季更甚。分別在5月與8月，每個(gè)月選取2景覆蓋研究區(qū)域時(shí)間間隔12 d的Sentinel-1數(shù)據(jù)進(jìn)行差分干涉處理，獲取了相干系數(shù)對(duì)比圖(見圖11)。從圖11中可以看到時(shí)間間隔12 d之內(nèi)，相關(guān)系數(shù)圖在山區(qū)相干性較差，大部分區(qū)域相關(guān)系數(shù)低于0.2，因此在進(jìn)行相位解纏的時(shí)候，絕大部分會(huì)被掩膜掉，同時(shí)滑坡易發(fā)區(qū)域位于山區(qū)，因此會(huì)造成監(jiān)測(cè)區(qū)域缺少有效信息。結(jié)合圖7可知該區(qū)域的滑坡隱患變形主要發(fā)生在5～7月，但是這段時(shí)間因?yàn)橹脖幻苁沟没赟entinel-1數(shù)據(jù)進(jìn)行差分干涉處理，獲取到的監(jiān)測(cè)信息效果較差。

圖11 不同時(shí)間段相關(guān)系數(shù)(間隔12 d)Fig.11 Correlation coefficient diagram for different time periods(interval of 12 days)

為了提高該區(qū)域的隱患識(shí)別能力和監(jiān)測(cè)效果，經(jīng)過上面的對(duì)比和結(jié)合文獻(xiàn)[30]，認(rèn)為在該區(qū)域進(jìn)行滑坡隱患的早期識(shí)別，應(yīng)選擇分辨率更高或者波長(zhǎng)更長(zhǎng)的SAR數(shù)據(jù)，比如日本的ALOS-2衛(wèi)星數(shù)據(jù)，阿根廷的SAOCOM-1A衛(wèi)星數(shù)據(jù)，以及即將發(fā)射的國(guó)產(chǎn)SAR衛(wèi)星。而Sentinel-1數(shù)據(jù)因高時(shí)間分辨率、免費(fèi)獲取可以用來(lái)對(duì)典型滑坡隱患開展時(shí)序監(jiān)測(cè)分析。

4.2 Stacking InSAR技術(shù)適用性分析

在本文中用Stacking 技術(shù)進(jìn)行隱患識(shí)別的前提是假設(shè)識(shí)別對(duì)象變形為線性形變，進(jìn)而獲取了地表形變結(jié)果。但是對(duì)該區(qū)域的滑坡通過SBAS InSAR技術(shù)進(jìn)行時(shí)序監(jiān)測(cè)，發(fā)現(xiàn)其變形具有一定的周期性，這與假設(shè)相矛盾，仍然可以識(shí)別出滑坡隱患的原因有：① 收集到的ALOS-2數(shù)據(jù)集，在夏季有影像數(shù)據(jù)，使得可以很好獲取降雨期間的變形信息。② Sentinel-1數(shù)據(jù)在5～8月雖然相干性較差，但是時(shí)間分辨率高，差分干涉組合多，因此在加權(quán)平均過程中可以獲取到有效變形信息。

需要指出來(lái)的是Stacking技術(shù)在該區(qū)域進(jìn)行滑坡隱患識(shí)別時(shí)，可能因?yàn)槠骄?jì)算使得獲取的地表形變形變速率結(jié)果值偏小，部分變形速率較小的滑坡隱患存在漏判現(xiàn)象，如文中的宋家灣滑坡，在基于Sentinel-1數(shù)據(jù)的Stacking結(jié)果中該滑坡隱患變形效果并不明顯，但是在累積形變量中均有較強(qiáng)烈的表現(xiàn)。

5 結(jié) 論

本文以龍泉山區(qū)域?yàn)槔捎脮r(shí)序InSAR技術(shù)利用多源SAR數(shù)據(jù)，結(jié)合光學(xué)遙感技術(shù)利用高分光學(xué)數(shù)據(jù)對(duì)該區(qū)域滑坡隱患開展識(shí)別研究，基于SBAS技術(shù)對(duì)典型滑坡隱患開展時(shí)序監(jiān)測(cè)分析，并對(duì)識(shí)別的隱患進(jìn)行野外驗(yàn)證。主要有以下結(jié)論：

(1) 利用Stacking技術(shù)、結(jié)合光學(xué)遙感技術(shù)成功識(shí)別出8處滑坡隱患，并通過野外驗(yàn)證驗(yàn)證了其可靠性，發(fā)現(xiàn)ALOS-2數(shù)據(jù)的識(shí)別效果優(yōu)于Sentinel-1數(shù)據(jù)，也證明了時(shí)序InSAR技術(shù)用于該區(qū)域的滑坡隱患識(shí)別具有可行性。

(2) 對(duì)基于ALOS-2數(shù)據(jù)識(shí)別的鄭家房子滑坡，通過時(shí)序監(jiān)測(cè)發(fā)現(xiàn)其最大形變速率為-40 mm/a，最大累積形變量為134 mm，并且其變形規(guī)律與汛期相吻合。通過野外調(diào)查發(fā)現(xiàn)，該滑坡在強(qiáng)降雨作用下，容易發(fā)生大規(guī)模失移，威脅坡體居民區(qū)、村道，因此需要在雨季加強(qiáng)持續(xù)性監(jiān)測(cè)，同時(shí)結(jié)果表明多種時(shí)序InSAR技術(shù)相結(jié)合用于滑坡隱患識(shí)別監(jiān)測(cè)具有較強(qiáng)可行性和可靠性。

(3) 通過對(duì)比不同波段SAR數(shù)據(jù)的識(shí)別結(jié)果及影響因素分析，發(fā)現(xiàn)影響該區(qū)域識(shí)別效果主要是干涉相干性，因此建議該區(qū)域采用分辨率更高或者波長(zhǎng)較長(zhǎng)的SAR數(shù)據(jù)開展滑坡隱患識(shí)別工作。

(4) 通過討論Stacking InSAR技術(shù)的適用性，分析了部分滑坡隱患漏判的原因，建議采用多種時(shí)序InSAR技術(shù)相結(jié)合形式開展滑坡隱患識(shí)別與監(jiān)測(cè)工作。