吳瑋瑩 王曉青 鄧 飛

?

基于高分衛(wèi)星遙感影像的地震應(yīng)急滑坡編目與分布特征探討——以2017年8月8日九寨溝7.0級地震為例1

吳瑋瑩 王曉青 鄧 飛

（中國地震局地震預(yù)測研究所，北京 100036）

地震應(yīng)急是減輕地震災(zāi)害的重要途徑之一。地震應(yīng)急工作具有時(shí)間緊迫、事關(guān)重大的特點(diǎn)。2017年8月8日四川九寨溝S7.0級地震發(fā)生后，為快速、準(zhǔn)確地提供地震引發(fā)的滑坡災(zāi)害分布，本研究基于震后第一天獲取到的高分辨率遙感影像（高分二號衛(wèi)星影像、北京二號衛(wèi)星影像），通過人工目視解譯的方法初步建立了四川九寨溝地震滑坡編目。結(jié)果表明，該地震至少觸發(fā)了622處同震滑坡，分布在沿使用影像邊界框定的面積為3919km2的區(qū)域內(nèi)。本研究還利用這個(gè)地震滑坡編目，統(tǒng)計(jì)了九寨溝地震滑坡數(shù)量和滑坡點(diǎn)密度（LND）與地形（坡度、坡向）、地震（地震烈度、震中距）等因素的關(guān)系。結(jié)果表明九寨溝地震滑坡多發(fā)生在坡度為20°—50°的區(qū)域內(nèi)，滑坡的易發(fā)性隨著坡度的增加而增加。受地震波傳播方向的影響，E、SE向是地震滑坡較易發(fā)生的坡向?；碌囊装l(fā)程度和地震烈度呈正相關(guān)，即隨著烈度的增大，滑坡易發(fā)性增大。滑坡易發(fā)性還隨著震中距增加而降低，這是由于地震波能量隨震中距的增加而衰減導(dǎo)致的。

四川九寨溝地震 地震應(yīng)急 同震滑坡分布

引言

地震是危害嚴(yán)重的自然災(zāi)害之一。提高地震監(jiān)測預(yù)報(bào)水平、加強(qiáng)震害防御能力和提高地震應(yīng)急救援水平是目前國內(nèi)外減輕地震災(zāi)害的3個(gè)主要途徑。在現(xiàn)有科技和經(jīng)濟(jì)條件下，地震應(yīng)急具有“較為現(xiàn)實(shí)、投入較少、見效較快、實(shí)效顯著”等優(yōu)點(diǎn)，受到世界各國的普遍重視。從理論上講，地震應(yīng)急過程實(shí)際上就是一系列在有限時(shí)間內(nèi)所作出的決策與實(shí)施過程的集合，具有時(shí)間緊迫、事關(guān)重大的特點(diǎn)。因此，在一定的時(shí)間約束條件下，地震應(yīng)急決策和行動是否科學(xué)合理和各種決策實(shí)施過程是否及時(shí)有效，是地震應(yīng)急工作成敗的關(guān)鍵。而決策的科學(xué)性和合理性以及決策實(shí)施的時(shí)效性和實(shí)際效用，根本上取決于決策和決策實(shí)施過程是否建立在全面、準(zhǔn)確、及時(shí)的地震災(zāi)情信息的基礎(chǔ)上（聶高眾等，2002，2012；蘇桂武等，2003）。能否快速提取和分析地震災(zāi)情信息并形成決策是地震應(yīng)急工作能否順利開展及發(fā)揮實(shí)際作用的核心。

由于我國地形條件復(fù)雜，地震往往導(dǎo)致滑坡、崩塌、泥石流等地質(zhì)災(zāi)害發(fā)生。地震地質(zhì)災(zāi)害會造成人員傷亡、房屋毀壞、道路中斷，極大地阻礙災(zāi)區(qū)快速開展地震應(yīng)急救援工作，因此快速獲得震區(qū)地震地質(zhì)災(zāi)害分布及其規(guī)律十分重要。在地震應(yīng)急工作中，傳統(tǒng)的方法在震后快速組織相關(guān)專家前往災(zāi)區(qū)調(diào)查存在耗時(shí)長、費(fèi)用高、無法保障調(diào)查人員安全等問題（熊德清，2009；趙紀(jì)生等，2008；唐川等，2013；劉剛等，2015）。隨著我國遙感技術(shù)的發(fā)展及其在地震應(yīng)急工作中的應(yīng)用，震后第一時(shí)間已能夠獲取到大范圍覆蓋震區(qū)的遙感影像，基于這些影像的人工目視解譯方法成為較為快速、準(zhǔn)確獲得震后地震地質(zhì)災(zāi)害及其分布信息的方法（胡國超，2009；曾濤等，2009；秦緒文，2012；劉睿等，2013）。

北京時(shí)間2017年8月8日21時(shí)19分，四川省阿壩藏族羌族自治州九寨溝縣發(fā)生7.0級地震。根據(jù)中國地震臺網(wǎng)測定，震中（33.20°N，103.82°E）位于九寨溝風(fēng)景區(qū)比芒村（距離九寨溝縣城39km），震源深度20km。此次地震震后第一天獲取到震中附近高分辨率的遙感影像，隨即基于遙感影像通過人工目視解譯的方法得到了地震應(yīng)急階段地震滑坡初步編目成果。本文利用該地震滑坡編目初步研究了九寨溝地震滑坡空間分布與地形（坡度、坡向）、地震（地震烈度、距震中的距離）等滑坡影響因素之間的關(guān)系，為地震應(yīng)急救援、震后重建，以及地質(zhì)災(zāi)害預(yù)防提供了可供參考的結(jié)果。

1 研究區(qū)概況

九寨溝位于四川省阿壩藏族羌族自治州九寨溝縣境內(nèi)，地處青藏高原東緣向四川盆地過渡地帶，是一條縱深50余千米的山溝谷地，總面積650.74km2。該區(qū)域內(nèi)地勢南高北低，山谷深切，高差懸殊，區(qū)域北緣九寨溝口海拔僅2000m，中部峰嶺均4000m以上，南緣達(dá)4500m以上。九寨溝地質(zhì)背景復(fù)雜，碳酸鹽巖分布廣泛，褶皺斷裂發(fā)育，新構(gòu)造運(yùn)動強(qiáng)烈，地殼抬升幅度大，造就了多種地貌（郭建強(qiáng)等，2000；楊更，2005；戴嵐欣等，2017）。四川九寨溝地震發(fā)生在南北地震帶中斷，青藏高原東緣巴顏喀拉塊體與東部華南地塊碰撞擠壓邊界帶上，震中及其鄰區(qū)發(fā)育有NWW向東昆侖斷裂東端塔藏?cái)嗔选NW向虎牙斷裂和近南北向岷江斷裂（圖1）（徐錫偉等，2017）。截至北京時(shí)間2017年8月18日09時(shí)00分，本次九寨溝地震共記錄到余震總數(shù)為5019個(gè)，其中4.0—4.9級地震3個(gè)，3.0—3.9級地震28個(gè)，最大余震為北京時(shí)間8月9日10時(shí)17分發(fā)生的九寨溝縣4.8級地震。

2 數(shù)據(jù)和方法

2.1 數(shù)據(jù)

本研究使用的高分辨率震后影像包括：①中國資源衛(wèi)星應(yīng)用中心提供的高分二號衛(wèi)星影像，共4景（圖2），包括分辨率為1m的全色影像與分辨率為4m的多光譜影像，影像的獲取時(shí)間均為2017年8月9日；②二十一世紀(jì)空間技術(shù)應(yīng)用股份有限公司提供的北京二號衛(wèi)星影像，共2景（圖2），包括分辨率為0.8m的全色影像與分辨率為3.2m的多光譜影像，影像獲取時(shí)間均為2017年8月9日。震后影像獲取時(shí)間為震后一天，能夠滿足地震應(yīng)急工作的時(shí)效性，保證地震滑坡編目工作能夠迅速開展。震前影像來源于Google Earth平臺，影像獲取時(shí)間為2011年12月10日，用于排除震前已經(jīng)存在的滑坡。將影像導(dǎo)入ArcGIS平臺中，并且在該平臺上完成地震滑坡編錄工作。地形數(shù)據(jù)利用30m分辨率的ASTER DEM（http://www. gscloud.cn/）以及從該DEM中提取的坡度、坡向。地震數(shù)據(jù)為中國地震局發(fā)布的四川九寨溝地震震中與地震烈度（http://www.cea.gov.cn/）。

圖1 四川九寨溝地震震中位置與主要活動斷層分布

2.2 方法

目前，建立同震滑坡編目的方法主要有5種（Keefer，2002；Guzzetti等，2012；Xu，2015）：①野外調(diào)查；②基于航片的目視解譯；③基于衛(wèi)星影像的目視解譯；④結(jié)合航片和衛(wèi)星影像的目視解譯；⑤基于遙感影像的滑坡自動提取。其中高分辨率衛(wèi)星影像目視解譯的方法具有能夠覆蓋較大的震區(qū)范圍、較強(qiáng)的中小規(guī)模滑坡的解譯能力，同時(shí)有獲取時(shí)間相對較短且價(jià)格低廉等優(yōu)點(diǎn)，成為當(dāng)下開展地震滑坡編錄的首選方法。野外調(diào)查可作為驗(yàn)證結(jié)果可靠性、提高滑坡編錄質(zhì)量的方法。本文首先建立同震滑坡的解譯標(biāo)志，采用目視解譯高分辨率遙感影像提取地震滑坡的方法建立地震應(yīng)急滑坡編目。同時(shí)，通過對比Google Earth平臺提供的震前遙感影像，排除了震前已存在的滑坡，從而建立了本次研究的九寨溝地震滑坡編目。在同震滑坡編目基礎(chǔ)上，以點(diǎn)要素表示單體滑坡，通過ArcGIS平臺顯示其空間展布，并研究其與地形因子（坡度、坡向）和地震因子（地震烈度、震中距）等同震滑坡空間分布影響因子之間的關(guān)系。選取滑坡數(shù)量（Landslide Number）和滑坡點(diǎn)密度（Landslide Number Density, LND）作為2個(gè)衡量滑坡相對發(fā)生程度（易發(fā)性）的參數(shù)，通過空間分析探索地震滑坡與地形因子和地震因子的關(guān)系。圖3給出了本研究處理的主要流程。

圖2 四川九寨溝地震等烈度線與用于地震滑坡解譯的震后遙感影像分布

3 結(jié)果

3.1 滑坡解譯標(biāo)志

解譯標(biāo)志是建立滑坡與遙感影像聯(lián)系的關(guān)鍵。在對高分辨率遙感影像進(jìn)行解譯前，首先需要建立同震滑坡的解譯標(biāo)志。建立解譯標(biāo)志并沒有統(tǒng)一的標(biāo)準(zhǔn)，解譯者基于經(jīng)驗(yàn)和對一系列能夠在影像上識別出來的特征進(jìn)行分析來發(fā)現(xiàn)和描繪滑坡的形態(tài)。特征主要包括形狀、大小、影像色彩、色調(diào)、陰影、紋理和物體的類型等。解譯者還可將遙感影像與數(shù)字高程模型結(jié)合，將影像以三維方式顯示，依據(jù)地勢、位置等信息對滑坡進(jìn)行識別和描繪（Guzzetti等，2012）。

圖3 本研究處理流程圖

解譯標(biāo)志的建立還需考慮研究區(qū)的實(shí)際情況，得到該區(qū)域較為合適的解譯標(biāo)志。本研究基于四川九寨溝地區(qū)的實(shí)際情況，結(jié)合遙感影像特征，建立了以下兩種解譯標(biāo)志：①由于道路對于斜坡具有削切作用，因此當(dāng)受到地震波的影響，道路兩旁更易發(fā)生地震滑坡災(zāi)害，阻塞道路，對地震應(yīng)急救援工作造成阻礙。道路在遙感影像上呈現(xiàn)為清晰的線條狀，若道路兩旁的邊坡發(fā)生地震滑坡災(zāi)害，則滑坡的堆積物質(zhì)使得道路原本整齊的邊界遭到破壞，能夠清晰地從遙感影像中識別出阻塞點(diǎn)以及道路路基的坍塌。因此將道路作為此次地震的解譯標(biāo)志之一（如圖4（a）、4（b））。②地震滑坡災(zāi)害的發(fā)生會導(dǎo)致土地覆蓋類型的變化。由于地震發(fā)生正值夏季，四川九寨溝地區(qū)植被茂盛，滑坡災(zāi)害的發(fā)生會破壞山體表面的植被覆蓋，與植被覆蓋區(qū)域在遙感影像上反差較大。因此，將植被作為此次地震的滑坡解譯標(biāo)志之一（如圖4（c）、4（d））。

3.2 同震滑坡編目展示

依據(jù)上述研究思路，解譯了九寨溝地震震后高分遙感數(shù)據(jù)覆蓋范圍內(nèi)的同震滑坡，提取結(jié)果如圖3所示。結(jié)果表明九寨溝地震至少觸發(fā)了622處滑坡。根據(jù)所使用的影像范圍及解譯滑坡分布情況，沿所使用影像的外邊界框定出滑坡的分布范圍，面積為3919km2，去除云層覆蓋區(qū)域后，研究區(qū)面積為2506km2，滑坡點(diǎn)密度為622/2506km2＝0.24處/km2。該研究區(qū)長軸方向?yàn)楸蔽鳌蠔|方向，與中國地震局發(fā)布的地震烈度圈長軸方向一致，大多數(shù)滑坡發(fā)生在地震烈度為Ⅷ、Ⅸ區(qū)域內(nèi)（圖5）。

3.3 同震滑坡空間分布的影響因素

選擇坡度和坡向作為地形因素，統(tǒng)計(jì)他們與地震滑坡空間分布的關(guān)系（如圖6（a））。研究區(qū)坡度范圍為0°—80°。在GIS平臺中通過30m×30m分辨率的DEM提取得到坡度數(shù)據(jù)，并將其按照10°坡度間隔分為7個(gè)區(qū)間，統(tǒng)計(jì)每個(gè)區(qū)間內(nèi)滑坡的數(shù)量和區(qū)間面積，計(jì)算得到區(qū)間內(nèi)滑坡點(diǎn)密度（LND）。從圖6（a）中可以看出，隨著坡度的增大，滑坡的個(gè)數(shù)和LND值都呈現(xiàn)先增后減的趨勢，滑坡數(shù)量最大的區(qū)間為30°—40°，共有滑坡238處；LND值最大區(qū)間為60°—70°，可達(dá)1.15處/km2。總體上，滑坡多發(fā)生在坡度為20°～50°的區(qū)域內(nèi)，滑坡的易發(fā)性隨著坡度的增加而增加。在GIS平臺中從30m×30m分辨率的DEM中提取坡向，將坡向分為Flat、N、NE、E、SE、S、SW、W、NW共9類。圖6（b）顯示了每個(gè)分類區(qū)間內(nèi)滑坡的數(shù)量和LND的統(tǒng)計(jì)結(jié)果。結(jié)果表明N、E、SE向滑坡數(shù)量較多，分別為97、106、111處。滑坡點(diǎn)密度在E、SE向較大，分別為0.33處/km2和0.37處/km2，這說明E、SE向應(yīng)為地震滑坡較為易發(fā)的坡向。

（a）（b）為道路解譯標(biāo)志；（c）（d）為植被解譯標(biāo)志；（a）（c）為震前GF2號影像，分辨率為1m，拍攝時(shí)間為2017年7月1日；（b）（d）為震后GF2號影像，分辨率為1m，拍攝時(shí)間為2017年8月9日

圖5 四川九寨溝地震滑坡分布圖

圖6 同震滑坡與坡度（a）坡向（b）烈度（c）震中距（d）分布關(guān)系

選擇地震烈度和震中距作為地震因素，統(tǒng)計(jì)他們與地震滑坡空間分布的關(guān)系。從圖2可以看出Ⅸ度區(qū)域面積較小。圖6（c）顯示研究區(qū)包含了Ⅵ、Ⅶ、Ⅷ、Ⅸ烈度區(qū)，滑坡數(shù)量隨著烈度的增大呈先增后減的趨勢，在Ⅷ度區(qū)達(dá)到最大，為383處?；旅芏入S著烈度的增大而增大，在Ⅸ度區(qū)達(dá)到最大，為0.60處/km2。從總體上看，滑坡的LND值與地震烈度呈現(xiàn)較好的相關(guān)性，即滑坡的LND值隨著烈度值的增大而增大，這也說明了烈度越大，越易發(fā)育滑坡。

在GIS平臺下以震中為圓心、10km為半徑間隔建立緩沖區(qū)文件，將整個(gè)研究區(qū)按照與震中的距離劃分為8個(gè)區(qū)間，然后開展滑坡分布與距離震中關(guān)系的統(tǒng)計(jì)（如圖6（d））。從中可以看出，滑坡數(shù)量先隨震中距的增加而增加，在距震中20km處達(dá)到最大，為243處。LND曲線隨著震中距的增加先減小，在震中距10km內(nèi)為最大值0.83處/km2，而后在距震中40km處出現(xiàn)局部峰值，為0.27處/km2，之后逐漸減少。從總體上看，隨著距離震中越遠(yuǎn)，同震滑坡易發(fā)性越低。

4 討論與結(jié)論

坡度對滑坡的發(fā)育具有重要的影響。斜坡的坡度與可能轉(zhuǎn)化為滑動面的坡體結(jié)構(gòu)面傾角決定了坡體臨空面能否成為滑坡發(fā)育的有效臨空面。同時(shí)，坡度應(yīng)力分布也決定著滑坡的穩(wěn)定性，地形坡度越陡則越容易引發(fā)滑坡（樊曉一，2013）。根據(jù)本文研究發(fā)現(xiàn)，四川九寨溝地震滑坡易發(fā)坡度范圍為20°—50°，滑坡的易發(fā)性隨著坡度的增加而增加。但在高坡度地區(qū)，滑坡數(shù)量與LND值減少，是因?yàn)樵撈露雀采w區(qū)域較少。這一現(xiàn)象在其他震例中也有類似的表現(xiàn)（Wang等，2007；Xu等，2015）。一些研究成果也表明相較于高坡度地區(qū)，滑坡更易發(fā)生在中等坡度的區(qū)域，這主要是由于在較大的斜坡梯度上沒有發(fā)生滑坡的物質(zhì)基礎(chǔ)（Masson等，2006；Reis等，2012）。

地震滑坡易發(fā)坡向受到區(qū)域內(nèi)坡向分布特點(diǎn)、地震波傳播方向以及區(qū)域其他地形地質(zhì)條件的影響。根據(jù)震后中國地震局地球物理研究所及北京大學(xué)張勇等（2017）對此次地震震源破裂過程的反演結(jié)果以及張旭等（2017）重新篩選波形資料并收集覆蓋震中區(qū)的InSAR資料對主震的震源破裂過程重新反演的分析結(jié)果表明，此次地震的破裂區(qū)可分為兩個(gè)主要凹凸體，一個(gè)較大且距起始破裂點(diǎn)較近，從起始破裂點(diǎn)北西較深的位置延展到起始破裂點(diǎn)南東較淺的位置；另一凹凸體較小且距起始破裂點(diǎn)較遠(yuǎn)，位于起始破裂點(diǎn)北西方向，正滑分量明顯。根據(jù)破裂過程可以得出，地震波是沿NW—SE向傳播，與地震波傳播方向一致的斜坡更易發(fā)生滑坡，即SE向?yàn)橐装l(fā)坡向。由于絕大多數(shù)滑坡的滑動和運(yùn)動方向基本與斜坡面傾向一致，研究區(qū)內(nèi)E向坡向分布范圍較廣，因此E向也為地震滑坡較為易發(fā)的坡向（許強(qiáng)等，2010）。

地震烈度是地震震動強(qiáng)度的客觀反映。在本研究對于地震烈度與地震滑坡數(shù)量的分析中，地震滑坡數(shù)量隨著烈度的增大先增加后減小。地震滑數(shù)量增加的原因是：根據(jù)Newmark累積位移量理論，坡體水平波動振蕩而產(chǎn)生的位移量與地震時(shí)地面的烈度呈相關(guān)關(guān)系。地震時(shí)地面烈度越高，地面水平加速度越大，坡體越接近于極限平衡狀態(tài)，則其它條件相同時(shí)，所產(chǎn)生的位移量也越大，越易發(fā)生滑坡（毛彥龍，1998）。地震滑坡數(shù)量在Ⅸ度區(qū)減少的原因可能是：①Ⅸ度區(qū)覆蓋面積較小；②Ⅸ度區(qū)在遙感影像中云覆蓋區(qū)域較多，影響解譯數(shù)量的準(zhǔn)確性。

通過探究此次地震滑坡在不同震中距下滑坡數(shù)量、LND值的變化趨勢，我們發(fā)現(xiàn)，從總體上看，隨著距離震中越遠(yuǎn)，同震滑坡易發(fā)性越低。這是由于地震波傳播過程中，能量隨震中距的增加而衰減，因而地面振動減弱，滑坡數(shù)量減少。在距震中40km處LND值出現(xiàn)小局部峰值與距震中80km處滑坡數(shù)量突然增多，可能是因?yàn)椋孩僬鹬卸ㄎ徊粔驕?zhǔn)確；②可能受到其它地形條件（坡度、坡向）等影響。

本研究利用震后第一天獲取到的高分辨率遙感影像（高分二號、北京二號）通過人工目視解譯的方法初步建立了四川九寨溝地震同震滑坡分布圖。結(jié)果表明，四川九寨溝地震至少觸發(fā)了622處滑坡，滑坡點(diǎn)密度為0.24處/km2。通過統(tǒng)計(jì)地形（坡度、坡向）、地震（烈度、距震中的距離）各分級區(qū)間內(nèi)同震滑坡數(shù)量和滑坡點(diǎn)密度（LAD），初步分析了各影響因子與同震滑坡空間分布之間的關(guān)系。結(jié)果表明，滑坡多發(fā)生在坡度為20°—50°的區(qū)域內(nèi)，滑坡的易發(fā)性隨著坡度的增加而增加。受地震波傳播方向的影響，E、SE向是地震滑坡較易發(fā)生的坡向?；碌囊装l(fā)程度和地震烈度呈正相關(guān)，即隨著烈度的增大，滑坡易發(fā)性增大?；乱装l(fā)性還隨著震中距增加而降低。本文研究存在影像未覆蓋整個(gè)地震震區(qū)、遙感影像存在云覆蓋區(qū)域等問題，但在地震應(yīng)急階段能夠滿足地震應(yīng)急的時(shí)效性要求，提供較為準(zhǔn)確的地震滑坡災(zāi)害分布，為震后救援、災(zāi)區(qū)重建、防治地震地質(zhì)災(zāi)害提供了具有價(jià)值的信息。

致謝：中國資源衛(wèi)星應(yīng)用中心和二十一世紀(jì)空間技術(shù)應(yīng)用股份有限公司提供了災(zāi)區(qū)衛(wèi)星遙感數(shù)據(jù)；地震預(yù)測研究所災(zāi)害信息研究中心科研團(tuán)隊(duì)成員參與了震后應(yīng)急滑坡提取工作。在此一并致謝！

戴嵐欣，許強(qiáng)，范宣梅等，2017．2017年8月8日四川九寨溝地震誘發(fā)地質(zhì)災(zāi)害空間分布規(guī)律及易發(fā)性評價(jià)初步研究．工程地質(zhì)學(xué)報(bào)，25（4）：1151—1164．

樊曉一，2013．西南地區(qū)地震滑坡坡度分布特征研究．四川地質(zhì)學(xué)報(bào)，33（3）：328—331．

郭建強(qiáng)，彭東，曹俊等，2000．四川九寨溝地貌與第四紀(jì)地質(zhì)．四川地質(zhì)學(xué)報(bào)，20（3）：183—192．

胡國超，2009．遙感技術(shù)在“5·12”地震重災(zāi)區(qū)汶川縣地質(zhì)災(zāi)害調(diào)查中的應(yīng)用．成都：成都理工大學(xué)．

劉剛，董樹文，侯春堂等，2015．對四川蘆山2013年S7.0地震的初步認(rèn)識——中國地質(zhì)科學(xué)院地震應(yīng)急響應(yīng)與地質(zhì)災(zāi)害排查．地質(zhì)論評，61（S1）：947—949．

劉睿，蘇偉，張曉東等，2013．蘆山地震震后次生滑坡災(zāi)害風(fēng)險(xiǎn)評價(jià)研究．震災(zāi)防御技術(shù)，8（4）：423—433．

毛彥龍，胡廣韜，趙法鎖等，1998．地震動觸發(fā)滑坡體滑動的機(jī)理．西安工程學(xué)院學(xué)報(bào)，20（4）：45—48．

聶高眾，陳建英，李志強(qiáng)等，2002．地震應(yīng)急基礎(chǔ)數(shù)據(jù)庫建設(shè)．地震，22（3）：105—112．

聶高眾，安基文，鄧硯，2012．地震應(yīng)急災(zāi)情服務(wù)進(jìn)展．地震地質(zhì)，34（4）：782—791．

秦緒文，2012．汶川地震災(zāi)區(qū)次生地質(zhì)災(zāi)害航天遙感應(yīng)急調(diào)查與綜合研究．中國科技成果，（18）：64—65．

蘇桂武，聶高眾，高建國，2003．地震應(yīng)急信息的特征、分類與作用．地震，23（3）：27—35．

唐川，常鳴，2013．2013．4.20蘆山S7.0級地震誘發(fā)次生地質(zhì)災(zāi)害應(yīng)急調(diào)查與分析．見：第十四屆海峽兩岸地貌學(xué)研討會論文集．青島：南京大學(xué)，臺灣大學(xué)．

熊德清，2009．四川省理縣5·12汶川地震誘發(fā)地質(zhì)災(zāi)害應(yīng)急排查與危險(xiǎn)性評價(jià)．成都：成都理工大學(xué)．

許強(qiáng)，李為樂，2010．汶川地震誘發(fā)滑坡方向效應(yīng)研究．四川大學(xué)學(xué)報(bào)（工程科學(xué)版），42（S1）：7—14．

徐錫偉，陳桂華，王啟欣等，2017．九寨溝地震發(fā)震斷層屬性及青藏高原東南緣現(xiàn)今應(yīng)變狀態(tài)討論．地球物理學(xué)報(bào)，60（10）：4018—4026．

楊更，2005．九寨溝景觀地質(zhì)背景及成因研究．成都：成都理工大學(xué)．

曾濤，楊武年，簡季，2009．無人機(jī)低空遙感影像處理在汶川地震地質(zhì)災(zāi)害信息快速勘測中的應(yīng)用．測繪科學(xué)，34（S2）：64—65，55．

張旭，馮萬鵬，許力生等，2017．2017年九寨溝S7.0級地震震源過程反演與烈度估計(jì)．地球物理學(xué)報(bào)，60（10）：4105—4116．

張勇，許力生，陳運(yùn)泰，2017．http://www.csi.ac.cn/manage/eqDown/05LargeEQ/201708082119M7.0/zonghe. html#_圖16、地震破裂過程圖

趙紀(jì)生，魏景芝，吳景發(fā)等，2008．汶川8.0級地震滑坡、崩塌機(jī)制．震災(zāi)防御技術(shù)，3（4）：379—383．

Guzzetti F., Mondini A. C., Cardinali M., et al., 2012. Landslide inventory maps: new tools for an old problem. Earth-Science Reviews, 112(1—2): 42—66.

Keefer D. K., 2002. Investigating landslides caused by earthquakes – a historical review. Surveys in Geophysics, 23(6): 473—510.

Masson D. G., Harbitz C. B., Wynn R. B., et al., 2006. Submarine landslides: processes, triggers and hazard prediction. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 364(1845): 2009—2039.

Reis S., Yalcin A., Atasoy M., et al., 2012. Remote sensing and GIS—based landslide susceptibility mapping using frequency ratio and analytical hierarchy methods in Rize Province (NE Turkey). Environmental Earth Sciences, 66(7): 2063—2073.

Wang H. B., Sassa K., Xu W. Y., 2007. Analysis of a spatial distribution of landslides triggered by the 2004 Chuetsu earthquakes of Niigata Prefecture, Japan. Natural Hazards, 41(1): 43—60.

Xu C., 2015. Preparation of earthquake—triggered landslide inventory maps using remote sensing and GIS technologies: principles and case studies. Geoscience Frontiers, 6(6): 825—836.

Xu C., Xu X. W., Shyu J. B. H., 2015. Database and spatial distribution of landslides triggered by the Lushan, Chinaw6.6 earthquake of 20 April 2013. Geomorphology, 248: 77—92.

Compilation and Spatial Analysis of Co-seismic LandslideInventory by Using High-resolution Remote SensingImages in Earthquake Emergency Response—An Exampleof the JiuzhaigouS7.0 Earthquake on August 8, 2017

Wu Weiying, Wang Xiaoqing and Deng Fei

(Institute of Earthquake Science, China Earthquake Administration, Beijing 100036, China)

Earthquake emergency response is one of the significant ways to mitigate earthquake disasters. After aS7.0 earthquake occurred in Jiuzhaigou, Sichuan Province on August 8, 2017, we used high-resolution remote sensing images (Gaofen-2 satellite image, Beijing-2 satellite image) obtained at the second day after earthquake, primarily established the inventory of landslide triggered by the JiuzhaigouS7.0 earthquake through artificial visual interpretation in order to obtain distribution caused by earthquake rapidly and accurately. The results show that the earthquake triggered at least 622 co-seismic landslides along the boundary of the image covered by 3,919km2. We also statistically analyzed the relations of the number of landslides with topography factors (gradient, aspect) and seismic factors (seismic intensity, distance from epicenter). It is found that landslides mostly occurred at the slope of 20°—50°, and the susceptibility of occurrence increases with the increase of slope gradient. The aspects of E and SE are prone to landslide occurrence because the direction of gradient may be related to seismic wave propagation. It is well correlated between the susceptibility of occurrence and seismic intensity, which suggests that as the intensity increases, the landslide increases. The susceptibility of landslide decreases with the increase of distance from the epicenter. The results provide valuable information for earthquake relief, and are useful in reconstruction and prevention of earthquake-induced geological disasters.

Jiuzhaigou earthquake; Earthquake emergency response; Distribution of co-seismic landslide

10.11899/zzfy20170410

科技部重點(diǎn)研發(fā)課題（編號：2017YFB0504104）

2017-10-02

吳瑋瑩，女，生于1993年。碩士研究生。主要從事地震滑坡危險(xiǎn)性研究。E-mail：wuweiying512@126.com

吳瑋瑩，王曉青，鄧飛，2017．基于高分衛(wèi)星遙感影像的地震應(yīng)急滑坡編目與分布特征探討——以2017年8月8日九寨溝7.0級地震為例．震災(zāi)防御技術(shù)，12（4）：815—825．