基于GIS的甘肃省东乡县锁南镇地质灾害易发性和危险性分区评价.pdf

传源南环境1SSN1672-9064环境大视野CN35-1272/TK基于GIS 的甘肃省东乡县锁南镇地质灾害易发性和危险性分区评价赵欣1陈宁2(1甘肃省有色工程勘察设计研究院2西北师范大学地理与环境科学学院甘肃兰州730000甘肃兰州730070)摘要随着GIS技术的发展，其强大的数据管理、空间分析和可视化能力被国内外学者广泛应用于地质灾害易发性和风险性评价中。研究基于CIS技术，利用野外调查和遥感解译数据，对研究区内发育的17 0 处地质灾害隐惠点进行统计分析，确定了地质灾害的相关评价指标并采用量化评分法进行了易发性评价，在易发性评价的基础上，结合斜坡失稳范围和条件概率确定了危险性分区。结果表明：地质灾害高、中、低易发区和不发育区分别占研究区面积的2 3.6 2%、55.50%、16.2 1%、4.6 6%，地质灾害高、中、低危险区分别占研究区面积的2 2.17%、58.8 7%、18.96%。通过后续排查结果验证，分区结果与实际情况较为吻合，同时建议锁南镇未来的防治重点应侧重于城镇交通枢纽两侧及重要基础设施区，未来可基于多源观测模型来提高风险区划的精度，最大程度的管理好地质灾害频发区域。关键词地质灾害CIS易发性危险性锁南镇中图分类号：X821;P694文献标识码：A文章编号：16 7 2-90 6 4(2 0 2 3)0 5-12 3-0 4地质环境是人类生存和发展的基础，由于环境恶化或人类不合理的开发利用导致的地质灾害严重破坏人类赖以生存的环境，并且造成巨大的生命财产损失1,尤其在我国西北地区发生频率较高2 ,地质灾害风险管理难度较大3。东乡族自治县位于甘肃省中部西南面，临夏回族自治州东北部，以东乡族发祥地和丝绸古道南路上的重要通道而闻名，是全国唯一的以东乡族为主体的少数民族自治县。然而，区内沟壑密集，地形坡度较大且坡形、坡向多变，植被覆盖率低，雨季集中，暴雨频发，水土流失严重，加之区内特殊的黄土一泥岩组合的岩土体发育，崩塌、滑坡等地质灾害发生频繁，严重影响着境内人民群众的生命财产安全4。因此，如何选择合适的灾害易发性评价方法，建立评价模型，划分研究区内危险区域和易发性等级，为防灾减灾提供科学的指导，是当前的研究热点之一。1研究区概况锁南镇位于东乡县中部黄河上级支流巴谢河上游，临夏盆地北缘，地势为黄土覆盖的残梁及沟谷呈放射状向四周延伸，区内出露地层主要为新近系临夏组泥岩、砂岩与第四系黄土。1.1气象研究区地处黄土高原和青藏高原过渡地带，属温带半干旱气候。区内多年平均降水量536 mm,24h最大降水量6 4.7 mm（19 7 0 年8 月18 日），1h最大降水量45mm，降水主要集中于6 一9 月，占全年降水总量的7 0.4%，多年平均蒸发量为1421.2mm,是年均降水量的2.6 倍。降水多以暴雨和连阴雨收稿日期：2 0 2 3-0 2-2 0基金项目：甘肃省国土资源厅临夏回族自治州东乡族自治县重要城镇地质灾害调查与风险区划项目（2 0 136 2 2 92 6 0 2)作者简介：赵欣(197 5一），男，研究生，硕士，高级工程师，研究方向为水工环地质。通讯作者：陈宁(1998 一),男，硕士在读，研究方向为生态环境遥感。123的形式出现，年均暴雨天数2.9d,最长连续降雨天数15d,标准冻土深度12 0 cm。1.2地形特征锁南镇地处东乡县中部，属黄土丘陵地貌区，高程2068.1,2460.3m,高差达150,2 0 0 m,基本形成以锁南镇为中心，黄土覆盖的残梁及沟谷呈放射状向四周延伸的地势特征。地貌类型决定了地形的此起彼伏，呈现出地形坡度、坡形、坡向的多变性。1.2.1坡度特征利用1:1万DEM数据对东乡县锁南镇地形进行坡度因子提取，东乡县锁南镇坡度0,5)的面积占9.8%;5,10)的面积占6.3%；10，15）的面积占17.9%；15，2 0）的面积占18.1%；2 0,2 5)的面积占12.7%；2 5,30）的面积占10.1%；30,35）的面积占9.1%；35，45）的面积占12.7%；45，55)的面积占3.1%；坡度55地段极少，多为陡崖、断壁，由于坡度过大，投影在平面上面积较小，仅占0.2%。1.2.2高差特征区内斜坡坡高以30,12 0 m居多，占全区域面积的86.64%,其中坡高 6 0 m的地形占全县总面积的47.46%,坡高60,120）m的地形占46.93%，坡高12 0 m的占5.6 1%。这一比例也侧面反映出本区沟谷发育，大部分地带斜坡体高差较大,具有临空面，为滑坡、崩塌的形成提供了非常有利的地形条件，是区内地质灾害多发的重要因素之一。2023.NO.5.(2)传源南环境环境大视野ISSN1672-9064CN35-1272/TK1.3地貌类型与特征研究区地处黄土丘陵地貌单元内，锁南镇及其周边呈沟梁相间的地形特征，梁顶起伏不大，总体呈缓坡状的起伏形态,高差一般在2 0,30 m,宽10 0,2 0 0 m。研究区四周沟谷较发育,其沟谷底部与梁顶相对高差150,2 0 0 】m,沟谷多呈“V型,两侧坡度多在2 0,45,受雨水冲刷作用影响，近沟脑地带侧蚀和向源侵蚀严重，梁顶完整性较差，大部分地段已成为弗状地形。原地形沟谷侵蚀切割强烈，梁顶完整性差，受人工改造作用大部分地段已成为阶状地形。1.4新构造运动区内新构造运动明显，以垂直升降运动为主，形成典型的黄土沟壑地形。区内各河谷形成多级阶地，阶地多呈基座式，且在各河谷两岸分布呈不对称性，表现为震荡性及差异性抬升。新构造运动为本区滑坡的形成奠定了基础。1.5地质灾害发育特征根据调查资料，研究区内共发育有滑坡（潜在滑坡）和崩塌(潜在崩塌)2 大类共17 0 处地质灾害隐患点,其中滑坡(潜在滑坡）94处，占比55%；崩塌（潜在崩塌）6 6 处，占比45%。滑坡(潜在滑坡)剖面形态多以“阶梯型”为主，按照物质组成可分为土质滑坡（主要为黄土滑坡和残坡积层滑坡）、岩土质滑坡和岩质滑坡（泥岩滑坡)3类；按照滑体厚度可分浅层滑坡(0,10）m)、中层滑坡(10,2 5）m)、深层滑坡(2 5,50）m)和超深层滑坡（50 m）。区内崩塌（潜在崩塌)剖面形态多以“直线型”为主，以人工斜坡为主，自然斜坡次之，人工斜坡坡宽数米序号影响因素1地层岩性0.3850单一黄土2斜坡结构0.2652土质斜坡类型3坡度/（）4坡高/m斜坡变形50.045 8破坏特征序号影响因素1坡度/（）地层岩性与20.20岩土体结构3地质构造4坡高/m5崩塌发生规模写发生频率0.20至数十米不等，坡高5,30 m,坡度一般50,8 5；自然斜坡一般坡宽数十米至数百米不等，坡高50,2 0 0 m不等。2地质灾害易发性评价2.1评价方法2.1.1滑坡易发性评价本文通过对研究区的已发生滑坡灾害的规模、大小、地形特征、坡度、岩性等信息进行统计分析并依据滑坡形成条件，选择地层岩性、斜坡结构类型、坡度、坡高和斜坡变形破坏特征等5项影响因素（表1)采用易滑分值法进行滑坡易发程度综合评判，滑坡易滑性分值计算见式(1)。(1)i=1式中：E潜表示滑坡易滑性分值;x;表示滑坡易发性影响因素；i表示x;影响因素的权重。根据E滑值将滑坡易发性分为4级：E滑2 5对应高易发区,E滑=2 0,2 5）对应中易发区,E滑=15,20)对应低易发区,E滑550.116 7100存在不稳定斜坡40变形体；存在正在坡存在活动的滑坡表2 崩塌易发程度量化评分表权重严重(A)0.3055裂隙发育的单一黄土：块状、厚层状的40坚硬岩体、岩体破碎或软硬相间陡峭斜坡走向与区域性断裂二致，几0.10组断裂交汇部位，40褶皱核部，褶皱轴向与坡面平行0.20100中型以上崩塌发生或存在体，常有小40规模崩塌落石2023.NO.5.得分40黄土一泥岩顺向层状结构斜40坡、崩滑残留体斜30斜向层状结构斜10坡；水平层状地层1坡坡4045,554040,100具复活条件的滑30已稳定的古老滑10无滑坡现象坡存在量级划分得分中等(B)4045,55软硬相间的层状30片状变质岩体结构碎裂岩体断裂密集分布，褶皱轴向与坡面斜30褶皱轴向垂直手10发育交坡向4040,100小型崩塌发生或存在小型以上崩塌的危岩体，危石30有落石发生存在124中等(B)得分30泥岩、砂岩互层3030,45 3030,40得分轻微(C)3030,45节理裂隙不发育，3030,40 无崩塌危岩体，时10无崩塌现象轻微(C)得分10单一岩性砂岩逆向层状结构斜斜坡100,30100,30得分一般(D)100,3010岩性单一厚层状1岩体，软弱岩体断裂褶皱构造不1100,30一般(D)得分1111得分111传源卤环境1SSN1672-9064CN35-1272/TK环境大视野表示x影响因素的权重。根据E崩值将崩塌易发性分为4级：相似的区域划分为一个评价单元,以研究区地理底图为基础，E前2 5对应高易发区，E册=2 0,2 5）对应中易发区,E前=充分考虑“区内相似、区际差异”原则和地质环境条件的差异15,19)对应低易发区,E崩 15对应不发育区。性，将研究区划分为不同的评价单元。对于地质环境条件相同2.1.3易发性分区的单元，可将单元合并；反之，地质环境条件变化明显的单元，根据易发性评价方法评判出每个评价单元的地质灾害易可将评价单元进一步细划，主要根据单元内致灾体的类型、特发程度,并将其划分为A、B、C、D 4类,分别表示地质灾害高、征、空间密度等，以使评价结果更为合理。中、低易发区和不易发区，若有2 种或2 种以上的地质灾害高按照上述划分原则，结合已确定的17 0 处地质灾害隐患易发区重叠时，则取值为4（表3）。点和区内的地质环境条件，以锁南镇规划区斜坡地带为评价重点，并兼顾现状地质灾害和潜在地质灾害发育特征，按照上表3地质灾害易发程度取值表述划分方法，将研究区划分成2 0 0 个评价单元。地质灾害易发区划分类型高易发区中易发区低易发区不易发区滑坡4崩42.2评价单元根据甘肃省重要城镇地质灾害调查与风险区划技术导则，将研究区的每个斜坡单元滑坡体和崩塌体视为一个完整对象，将灾害体以外的斜坡带、冲沟及其邻近的地质环境背景面积/易发性分区km2高易发区3.089主要分布于研究区毛毛沟的沟脑斜坡地带与沟脑右岸斜坡地带。地质环境脆弱，植被覆盖率较低，沟谷两侧斜坡低缓区开垦指数大，地质灾害的分布发育与人类工程活动密切相关主要包括县城北部锁南村毛毛村一带与前沟西侧黄土梁两侧地带。地形起伏变化较大，地质中易发区7.81低易发区2.44不易发区0.75地带3地质灾害危险性评价3.1评价方法基于目前各类危险性评价方法的优缺点，以及锁南镇地质灾害的发育特征，本文在地质灾害易发性评价的基础上，以斜坡单元作为评价重点，兼顾现状地质灾害和潜在地质灾害易发特征，通过对坡体失稳范围的计算，考虑条件概率和动态变化过程，得出危险性分区结果。3.1.1滑坡(潜在滑坡)威胁范围对于斜坡变形的危害主要为其前缘危害区和其本身，在评价时，确定不同规模、频率地质灾害危害区范围。滑坡直接或间接的最大影响范围包括滑动区、滑体堆积区、裂缝区和堵断河流的影响区。其堆积危害区范围确定可按公式估算法计算。(1)水平滑移距离。滑坡前缘滑移面的最大威胁距离，即为滑体影响范围的最远点，滑移的最大距离与坡高、滑速和滑移面上的摩擦系数有关，见式(3)。L(x)=nH/(0.5t a n g)式中：L(x）表示滑坡体最大水平滑移距离，m；H 表示滑坡体前、后缘高差，m;n表示滑坡体滑出条件系数，可参考表5取0.1,1.0;表示黄土的综合内摩擦角()。2.3评价结果3232环境脆弱，植被覆盖率较低，地质灾害的分布发育与人类工程活动密切相关主要分布于县城主城区已进行地质灾害综合治理的区域与研究区局部坡体较缓、地层较稳定的地带。地质灾害相对发育较少主要分布于县城主城区已进行地质灾害综合治理的区域与研究区局部坡体较缓、地层较稳定的表5滑坡体滑出条件系数表滑出条件滑坡堆积区地势狭窄，临空条件差，滑体前方有大型阻隔体剪出口位于坡脚附近，滑坡堆积区地表起伏，分布有中、小型阻隔体剪出口位于坡脚附近，滑坡堆积区地表稍有起伏，无大型阻隔体滑坡临空面开阔，滑坡前方地势平坦，无阻隔体剪出口临空高度高，滑坡前方地势平坦，无阻隔体2018年7 月18 日，甘肃临夏、定西等地发生短时强降雨，暴雨引发灾害，特别是临夏州东乡县受灾严重，经对灾后锁南镇隐患点排查的数据分析，预测值与实测值比较接近。(2)横向散开范围。滑坡横向散开的范围与地形有关，一般取前缘宽的1,2 倍。(3)3.1.2滑坡概率计算通过计算滑坡诱发因素的年超越概率来进行失稳概率分析，研究区滑坡的主要诱发因素为地震活动、降水条件和人类工程活动。根据建筑抗震设计规范（GB50011一2 0 10)（2 0 162023.NO.5.12511表4地质灾害易发性评价区划结果分布特征根据上述易发性评价方法，在定性分析和定量计算的基础上，斜坡区评价单元与区内灾害体评价单元采取“就高原则”，即当定性分析的和定量计算得出的评价结论不一致时，取二者中的较高等级。将研究区划分为地质灾害高易发区、中易发区、低易发区和不发育区。其中，地质灾害中易发区面积占比最大约7.8 1km，高易发区面积3.0 8 9km,低易发区面积2.44km,不发育区面积0.7 5km(表4)。面积百价单元分比/%3623.6210555.505716.2124.66系数(n)0.1,0.30.3,0.50.5,0.70.7,0.90.9,1.0涉及评传源南环境环境大视野ISSN1672-9064CN35-1272/TK版），东乡县设防烈度为VI度，地震动峰值加速度按50 a超越概率为10%计算，诱发滑坡的概率为0.0 0 2 1。暴雨强度重现期按50 a一遇计算，则诱发滑坡的概率为0.0 2,2 0 a一遇降雨，诱发滑坡的概率为0.0 5。3.1.3崩塌(潜在崩塌)与易发性评价相同,崩塌（潜在崩塌)的威胁范围和概率计算参照滑坡（潜在滑坡）评价方法，进而进行危险性评价。根据调查，区内发育的崩塌（潜在崩塌）平面形态多呈规则分布，剖面形态多为“直线型”。人工不稳定斜坡坡宽一般数米至数十米不等，坡高5,30 m，坡度50,8 5；自然不稳定斜坡坡宽一般数十米至数百米不等，坡高50,2 0 0 m,坡体总体上缓下陡,坡度2 5,40 。大多在雨季有小规模的溜土、塌现象发生，同时受地形及人类工程活动的影响，易失稳严重威胁斜坡影响范围内的人员及财产的安全。危险性分区面积/km地质灾害高危险区主要分布毛毛沟的沟脑斜坡地带与沟脑右岸斜坡地带，主要威胁坡面居民及坡高危险区2.907.70中危险区低危险区4结论(1)研究区目前共发育17 0 处地质灾害隐患点,其中滑坡（潜在滑坡）占55%，崩塌（潜在崩塌）占45%。(2)本文在传统地质灾害危险性评价方法的基础上,将GIS技术应用于锁南镇地质灾害易发性评价。根据锁南镇地质灾害发育特征选取评价因子进行易发性评价，将研究区划分为地质灾害高易发区、中易发区、低易发区和不发育区，其面积分别为3.0 8 9、7.8 1、2.44和0.7 5km，中易发区占比最大为55.5%高易发区主要分布于研究区毛毛沟的沟脑斜坡地带与沟脑右岸斜坡地带。(3)基于易发性评价，将研究区划分为地质灾害高危险区、中危险区和低危险区，其面积分别为2.90.7.7 0 和2.48 km，中危险区占比最大为58.8 7%，高危险区主要分布于毛毛沟的沟脑斜坡地带与沟脑右岸斜坡地带。(4)本文的评价结果与实际情况吻合，锁南镇未来的防治重点应侧重于城镇交通枢纽两侧、重要基础设施区的黄土质滑坡。针对相同类型的灾害区划研究，未来可结合多重数据分(上接第10 0 页）电力,2 0 19,52(3)：10-15.8王争荣，耿宣，汪洋，等.燃煤电厂湿烟羽消除设计方案对比分析J.华电技术,2 0 18,40(9):5-9.9孙金博.基于Excel实现烟气除湿脱白工艺参数计算J.当代化工研究,2 0 2 3(0 3):16 4-16 6.3.2评价单元基于易发性评价划分的2 0 0 个评价单元和斜坡影响范围，将每个斜坡单元在暴雨工况下失稳的覆盖范围作为1个评价单元、地震工况下失稳超出暴雨工况下的覆盖范围再作为1个评价单元，而其他评价范围不再单独划分危险性评价单元，故危险性评价共划分出6 0 0 个评价单元。3.3评价结果依据地质灾害易发性分区和斜坡影响范围的计算结果，对上述每个评价单元的危险性进行分析，对各评价单元危险性分级赋值,其中高、中、低危险区分别赋值为4、3、2。按照地质灾害危险性评价方法，结合每个评价单元危险性计算结果，以“区内相似、区间相异”的原则,将研究区划分为地质灾害高危险区、中危险区和低危险区。地质灾害高危险区面积为2.90 km，中危险区面积为7.7 0 km,低危险区面积为2.48 km（表6)。表6 地质灾害危险性评价区划结果分布特征脚工厂。主要包括115、118、12 0、137 等6 3个评价单元地质灾害中危险区主要分布在县城马蹄状黄土梁外侧斜坡地带，斜坡对象主要为斜坡坡顶及坡面上部居民。主要包括2 1、2 2、2 4、2 6 等2 7 3个评价单元地质灾害低危险区主要分布于县城主城区已进行地灾害综合治理的区域与研究区局部坡度较缓的2.48斜坡地带。地质灾害相对发育较少。主要包括1、2、3、4、5等2 6 6 个评价单元与环境,2 0 2 2(0 1):6 1-6 3,6 5.11周西华，梁茵,王小毛.饱和水蒸气分压经验公式的比较J.辽宁工程技术大学学报，2 0 0 7,2 6(0 3)：331-333.百分比/%22.1758.8718.96析、机器学习算法的多源观测模型来提高风险区划的精度，最大程度管理好地质灾害频发区域。参考文献1罗元华，张梁，张业成.地质灾害风险评估方法M.北京：地质出版社,1998,2 FAN X M,SCARINGI G,KORUP G,et al.Earthquake-induced chainsof geologic hazards:patterns,mecha-nisms and impactsJJ.Reviews ofGeophysics,2019,57(2):421-503.3唐亚明，张茂省,李政国，等.国内外地质灾害风险管理对比及评述J.西北地质,2 0 15,48(2):2 38-2 46.4赵欣，韵建武,汪玮玮，等.临夏回族自治州东乡族自治县锁南镇地质灾害调查与风险区划报告R.甘肃有色工程勘察设计研究院，2014:1-187.5李松，张川.甘肃省东乡县特大滑坡地质灾害形成机制与治理方法J.甘肃地质,2 0 13,2 2(1):6 5-7 0.6王宏睿.基于大数据平台的县级地质灾害监测一体化应急管理体系构建D.兰州：兰州理工大学,2 0 2 0.10郑鸿.烟气消白烟道流体动力学数值模拟流场优化设计J.能源2023.NO.5.126