雅安市名山區(qū)地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)的構(gòu)建

雅安市名山區(qū)地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)的構(gòu)建

【摘要】地質(zhì)災(zāi)害的發(fā)生受控于多種影響因子，同時(shí)影響因子中各屬性類別對(duì)地質(zhì)災(zāi)害的發(fā)育分別產(chǎn)生促進(jìn)或抑制作用，因此，構(gòu)建合適的指標(biāo)體系對(duì)滑坡危險(xiǎn)性評(píng)價(jià)具有重要意義。本文選取雅安市名山區(qū)作為研究區(qū)，參照前人研究，選取7個(gè)影響因子，并對(duì)各影響因子及其各屬性類進(jìn)行敏感性分析，建立評(píng)價(jià)指標(biāo)體系，得到名山區(qū)危險(xiǎn)性分區(qū)圖?！綤ey】地質(zhì)災(zāi)害敏感性分析危險(xiǎn)性評(píng)價(jià)指標(biāo)名山區(qū)本文選擇四川省雅安市名山區(qū)作為研究區(qū)，獲取名山區(qū)112個(gè)地質(zhì)災(zāi)害隱患點(diǎn)，并參照前人研究選取7個(gè)影響因子，根據(jù)唐川[1]在對(duì)德國波恩地區(qū)滑坡因子進(jìn)行分析時(shí)提出的敏感性分析方法，對(duì)各影響因子及其各屬性進(jìn)行敏感性分析，實(shí)現(xiàn)對(duì)名山區(qū)地質(zhì)災(zāi)害的危險(xiǎn)性評(píng)價(jià)。1研究數(shù)據(jù)與研究方法1.1研究數(shù)據(jù)本研究中地質(zhì)災(zāi)害隱患點(diǎn)數(shù)據(jù)來源于四川省雅安市人民政府?；诿絽^(qū)30m分辨率的DEM數(shù)據(jù)生成名山區(qū)坡度圖、坡向圖及河流水系圖;利用地理空間數(shù)據(jù)云上下載的2017年7月Landsat8影像（研究區(qū)內(nèi)均無云）計(jì)算NDVI;所用到的名山區(qū)的地質(zhì)年代圖來源于1：20萬的地質(zhì)圖矢量化;所涉及的路網(wǎng)數(shù)據(jù)是通過GoogleEarth下載得到。1.3敏感性指數(shù)IL分析法敏感性分析法是利用ArcGIS對(duì)屬性數(shù)據(jù)進(jìn)行預(yù)處理，繼而對(duì)地質(zhì)災(zāi)害發(fā)育的影響因子進(jìn)行篩選，得出關(guān)鍵影響因子，建立區(qū)域地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)指標(biāo)體系。式中：IL表示地質(zhì)災(zāi)害危險(xiǎn)性指數(shù);Pmn（L）表示影響因子m的n個(gè)屬性類別中存在地質(zhì)災(zāi)害隱患點(diǎn)的柵格個(gè)數(shù);Pm（L）表示影響因子m中的所有地質(zhì)災(zāi)害隱患點(diǎn)的柵格個(gè)數(shù);Pmn（G）表示影響因子m的n個(gè)屬性類別中不存在地質(zhì)災(zāi)害隱患點(diǎn)的柵格個(gè)數(shù);Pm（G）表示影響因子m所占的總柵格個(gè)數(shù)。2評(píng)價(jià)指標(biāo)體系的構(gòu)建2.1影響因子敏感性分析本文選取河流、NDVI、坡度、地質(zhì)年代、坡向、人類工程活動(dòng)以及高程7個(gè)影響因子[2]，并對(duì)這7個(gè)影響因子的各屬性類進(jìn)行相對(duì)頻率組合定量計(jì)算，最終確定其滑坡敏感性指數(shù)IL。具體計(jì)算過程是：①在ArcGIS中將選取的7個(gè)影響因子的屬性用30m×30m柵格單元表示;②將名山區(qū)內(nèi)112個(gè)地質(zhì)災(zāi)害隱患點(diǎn)圖層添加到ArcGIS中不同屬性圖層的柵格單元中，提取地質(zhì)災(zāi)害隱患點(diǎn)所處位置中各個(gè)影響因子的屬性數(shù)據(jù)，用柵格單元個(gè)數(shù)表示，得出影響因子m的n個(gè)類別中存在地質(zhì)災(zāi)害隱患點(diǎn)的柵格個(gè)數(shù)Pmn（L）、影響因子m中所有滑坡點(diǎn)的柵格個(gè)數(shù)Pm（L）、影響因子m的n個(gè)類別中不存在地質(zhì)災(zāi)害隱患點(diǎn)的柵格個(gè)數(shù)Pmn（G）以及影響因子m所占的總柵格個(gè)數(shù)Pm（G）;③根據(jù)公式（1）至公式（3）計(jì)算得出各個(gè)影響因子的各屬性類別的敏感性;④根據(jù)影響因子各屬性類別的敏感性可得出各影響因子的平均敏感性。根據(jù)敏感性計(jì)算方法及Gokceoglu[3]等的研究成果，說明某個(gè)因子中某一屬性類別的IL值大于1，則地質(zhì)災(zāi)害隱患點(diǎn)與該屬性類別成正相關(guān)，其對(duì)地質(zhì)災(zāi)害隱患點(diǎn)的發(fā)育具有較明顯的促進(jìn)作用，地質(zhì)災(zāi)害隱患點(diǎn)在該屬性類內(nèi)發(fā)生的概率較大;若IL值小于1，則反之。當(dāng)IL值接近于1時(shí)，說明該因子與地質(zhì)災(zāi)害隱患點(diǎn)相關(guān)性不大，接近于整個(gè)區(qū)域內(nèi)的平均水平。根據(jù)各影響因子各屬性類別的敏感性得出各影響因子的平均敏感性，若某個(gè)影響因子的平均敏感性值小于1，則不作為該地區(qū)地質(zhì)災(zāi)害隱患點(diǎn)危險(xiǎn)性評(píng)價(jià)的參評(píng)因子。3.2關(guān)鍵影響因子篩選利用各個(gè)影響因子的各屬性類別的IL值，對(duì)于選的7個(gè)影響因子進(jìn)行平均敏感性計(jì)算，得出：地質(zhì)年代（2.600）>坡向（2.020）>人類工程活動(dòng)（1.839）>高程（1.139）>河流（0.954）>NDVI（0.796）>坡度（0.619）。由于河流、NDVI和坡度的平均敏感度小于1，因此予以剔除。最終確定地質(zhì)年代、坡向、人類工程活動(dòng)以及高程4個(gè)因子作為名山區(qū)地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)的參評(píng)因子。4地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)根據(jù)敏感性分析篩選出4個(gè)參評(píng)因子，并對(duì)其進(jìn)行歸一化和數(shù)值化處理，每個(gè)因子按照4-5級(jí)危險(xiǎn)度劃分，由此可作出相應(yīng)因子的評(píng)價(jià)圖[4]。在ArcGIS中利用RasterCalculate的疊加工具進(jìn)行疊加，得到的數(shù)值范圍是0～8.5。在ArcGIS中將其重新分類，劃分為無危險(xiǎn)區(qū)（0-1.759）、輕度危險(xiǎn)區(qū)（1.759-2.246）、中度危險(xiǎn)區(qū)（2.426-4.089）以及高危險(xiǎn)區(qū)（4.089-8.414），最后得到名山區(qū)的地質(zhì)災(zāi)害危險(xiǎn)性區(qū)劃圖。5結(jié)論（1）本文利用地質(zhì)災(zāi)害的敏感性指數(shù)對(duì)影響名山區(qū)地質(zhì)災(zāi)害發(fā)育的影響因子進(jìn)行了詳細(xì)分析，最終得出，在坡度介于15-35°，距離高速公路1km以內(nèi)，在水系200m以內(nèi)，NDVI大于0.3，高程在544-650m、850-1470m，坡向介于90-135°、225-350°，地質(zhì)年代為三迭系、下三系、侏羅系以及白堊系，對(duì)研究區(qū)地質(zhì)災(zāi)害的發(fā)育起到促進(jìn)作用。（2）從得出的結(jié)論中可以看出，中、高危險(xiǎn)區(qū)主要集中在名山區(qū)的東南方向和西北方向，占總面積的18.19%，其中高危險(xiǎn)區(qū)僅占總面積的0.31%，故可得出名山區(qū)總體較為安全。【Reference】[1]唐川，JorgGrunert.滑坡災(zāi)害評(píng)價(jià)原理和方法研究[J].地理學(xué)報(bào)，1998（S1）：149-157.[2]齊識(shí).白龍江流域滑坡危險(xiǎn)度評(píng)價(jià)技術(shù)研究[D].蘭州大學(xué)，2014.[3]C.G?kceoglu，H.Aksoy.LandslidesusceptibilitymappingoftheslopesintheresidualsoilsoftheMengenregion（Turkey）bydeterministicstabilityanalys