基于GIS和证据权模型的山阳县地质灾害易发性评价.pdf
《基于GIS和证据权模型的山阳县地质灾害易发性评价.pdf》由会员分享,可在线阅读,更多相关《基于GIS和证据权模型的山阳县地质灾害易发性评价.pdf(8页珍藏版)》请在咨信网上搜索。
1、新疆地质XINJIANGGEOLOGY2023年6月Jun.2023第41卷 第2期Vol.41No.2中图分类号:P694文献标识码:A文章编号:1000-8845(2023)02-262-08项目资助:山阳县地质灾害风险调查评价项目(2249-194HDZB19203-02)资助收稿日期:2023-02-15;修订日期:2023-04-04;作者E-mail:第一作者简介:王伟中(1994-),男,甘肃渭源人,硕士,2021年毕业于长安大学地质工程专业,研究方向:地质灾害风险评价、地质灾害防治通讯作者:李树兴(1995-),男,甘肃永登人,硕士研究生,2021年毕业于长安大学地质工程专业,
2、从事地质灾害防治、地质灾害风险评价等工作;E-mail:基于GIS和证据权模型的山阳县地质灾害易发性评价王伟中1,李树兴2,杨成1,许涛1,宋飞3,曹小红1,李浩然3,王伟华4(1.新疆地质灾害防治重点实验室(新疆工程学院),新疆 乌鲁木齐 830000;2.甘肃省建筑设计研究院有限公司,甘肃 兰州730030;3.长安大学地质工程与测绘学院,陕西 西安 710054;4.兰州理工大学土木工程学院,甘肃 兰州 730050)摘要:地质灾害易发性涉及因素众多,评价结果取决于不同模型和参量权重赋值为当前研究热点。以陕西省山阳县为研究区域,在收集整理区域地质环境条件、地质灾害分布资料、相关性分析等基
3、础上,最终选取坡度、坡向、坡型、工程地质岩组、距断层距离、距河流距离、距道路距离等8个影响因子,采用证据权模型对8个影响因子进行分析,使用GIS空间分析功能对研究区开展地质灾害易发性评价。结果显示,山阳县地质灾害易发等级可划分为高易发、中易发、低易发和非易发4个等级,面积分别为262.2 km2、436.7 km2、1141.6 km2、1694.5 km2,占山阳县总面积的7.48%、12.35%、32.29%、47.94%。进一步采用ROC曲线方法检验评价等级结果的可靠性,得到AUC为0.824 2(精度达82.42%)。关键词:Arc GIS;地质灾害;易发性评价;证据权模型;ROC曲线
4、我国地域广阔,地质地理条件存在巨大差异,长期地球演化进程导致地质环境条件复杂多变1。位于陕南秦巴山区的山阳县,是陕西省地质灾害灾情严峻区域2,对该区开展地质灾害易发性评价,能为山阳县地质灾害防治提供参考依据。3S技术在地质灾害易发性、危险性和风险评价中得到广泛应用。周静静等基于ArcGIS平台,使用信息量模型对陕南秦巴山区进行地质灾害危险性区划2;冯显杰采用信息量模型、信息量-BP神经网络耦合模型等对怒江州高山峡谷地质灾害进行易发性评价3;白光顺等采用证据权法,结合GIS技术对昆明市五华区地质灾害进行易发性分区4;刘文涛等以阆中市滑坡灾害为例,采用层次分析法进行滑坡易发性评价5;毕结昂等基于G
5、IS平台,采用信息量-逻辑回归耦合模型对新疆玛纳斯河流域进行地质灾害易发性分区6;王雪冬等利用GIS技术,选取曲率、水力侵蚀等13个地形指标因子,基于Logistic回归与随机森林模型对和龙市地质灾害易发性进行预测7;李怡飞、易靖松等借助GIS平台对青藏高原高山峡谷区进行地质灾害危险性分区评价,获得该区域评价结果8-9。证据权法是基于贝叶斯准则,综合各种证据预测某事件发生概率大小的一种定量方法,应用到地质灾害评价中,取得较好效果10-13。证据权法综合各种影响因子,有效解决不同级别因子对地质灾害易发性影响14。本文以陕西省山阳县为研究区域,选取坡度、坡向、坡高、坡型、工程地质岩组、距断层距离、
6、距河流距离、距道路距离等8个评价因子,采用证据权模型,基于ArcGIS空间分析功能,对山阳县地质灾害进行易发性评价分析,为山阳县地质灾害防治提供参考。1 模型介绍证据权模型是贝叶斯理论模型之一,能综合各种证据预测预判某个事件发生的概率13。基本原理为:假设某研究区内已知有n个图层,区内已知灾害点的先验概率为P0。利用GIS将整个研究区划分为T个面积相等的栅格单元。有地质灾害发生栅格数为D,D=T-D表示无地质灾害发生的栅格数,B表示在某一证据因子下地质灾害发生的栅格数,在此证据下无地质灾害发生的栅格数为B=T-B。由上述假设得到P0=D/T,对任何一个证据因子,D的权重为:第41卷第2期王伟中
7、等:基于GIS和证据权模型的山阳县地质灾害易发性评价W+=lnP()B/DP()B/D-W=lnP()B-/DP()B-/D-(1)其中:W+、W分别表示证据因子存在区和不存在区的权重值,其大小表示与证据因子的密切程度。为说明各因子权重计算方法,选取第i个证据为例,设第i(i=1,2,n)个因子内发生灾害的单元数为Di,不发生灾害的单元数Di,则第i个证据因子对地质灾害发生的权重可用下式计算:W-i=lnDi/DDi/Di-W-=ln1-()Di/D1-()D-i/D-(2)为表征地质灾害发生与不发生之间的贡献值,可计算相对系数,公式为:Ci=W+-W-(3)对n个因子,研究区任一单元k为地质
8、灾害发生的可能性,即对数后验概率大小可表示如下14-15:lnR=W0+inWki()i=1,2,n(4)式中:W0=lnP0/()1-P0,其值为对数先验概率。Wki为第i个因子的权重值,其中若存在证据因子时为正,反之为负。最后计算后验概率为:P后检=R1+R(5)后验概率大小表示地质灾害易发高低,其值在01之间。2 研究区概况山阳县位于商洛市南部,地理位置E10932E11029,N3309N3342,东西长88.5 km,南北宽60.7 km,面积3 535 km2。研究区属亚热带向暖温带过渡的半湿润山地气候,四季特征明显,冬无极寒,夏无酷暑。区内呈南高北低地势特点,西南部郧岭、中部鹘岭
9、、北部流岭与中部河川区构成“三山夹两川”地势地貌格局。研究区地层出露较全,从老到新主要有青白口系、震旦系、寒武系、奥陶系、志留系、泥盆系、石炭系、二叠系、白垩系、新近系和第四系。震旦系、寒武系、奥陶系、志留系分布于东南部地区;泥盆系、石炭系大面积分布;二叠系分布于西南角,古近系局部出露,第四系广泛分布于河谷阶地、阶形斜坡缓坡地带等。研究区出露地层划分为5个工程地质岩组,即:坚硬块状坚硬侵入岩体(A)、中厚层状坚硬碳酸盐岩组(B)、薄层状较弱浅变质岩岩组(C)、中厚层状软弱碎屑岩岩组(D)和土体(E)。研究区位于祁秦地槽东秦岭褶皱系,境内褶皱、断裂构造发育。据史料记载,研究区共发生地震19次,多
10、为小于5级的有感地震或弱震,最大地震震级为4.5级。据野外调查,研究区内共发育地质灾害465处,主要为滑坡、崩塌、泥石流、地面塌陷4种,其中滑坡445处、崩塌8处、泥石流10处、地面塌陷2处,分别占灾害点总数的 95.69%、1.72%、2.15%、0.43%(图1)。3 基于GIS地质灾害易发性评价3.1 数据来源本文数据来源:12.5 m12.5 m的DEM栅格数据、山阳县1 5万地质图、野外调查资料(表1)。为方便后期评价需要,本文将研究区DEM重采样生成15 m15 m的栅格作为最小栅格单元分析评价,共划分成15 711 413个栅格单元。3.2 评价因子分析在对山阳县自然地理和地质灾
11、害发育特征分析图1 山阳县地质灾害隐患点分布图Fig.1 Map of potential geological hazards in Shanyang County1.地质灾害隐患点;2.山阳县263新疆地质2023年基础上,结合野外调查资料、前人研究成果,选取坡度、坡向、坡高、坡型、工程地质岩组、距断层距离、距河流距离、距道路距离等9个影响因子,进行分析。3.2.1 坡度坡度是控制斜坡稳定性、破坏及运移方式的重要因素7。坡度越陡斜坡越易失稳;坡度越缓,相对越稳定。通过GIS分析,将研究区坡度划分为5个区段。从图2-1(a)可见,地质灾害多发生在1545斜坡上。坡度在030时,地质灾害密度和
12、数量随坡度增大而增加;坡度处于3045时,地质灾害密度和数量随坡度增加而降低;坡度在4590时,地质灾害数量随坡度增加而降低,灾害密度随坡度增加而增加。3.2.2 坡向不同斜坡坡向太阳辐射强度不同,影响水分蒸发、植被覆盖、坡面侵蚀等,并影响斜坡中孔隙水压力分布和岩土物理力学参数,导致斜坡稳定性差异16。通过GIS将整个研究区坡向划分为 8 个区段。从图2-(b)可见,在坡向为E、SE、S、SW向斜坡中,地质灾害发育强烈;SE、SW 向斜坡中地质灾害点密度发育较高。3.2.3 坡高斜坡高度(坡高)一方面决定地质灾害破坏强度,另一方面影响灾害发生后运移距离。通过 GIS 将整个研究区坡高划分为5个
13、区段。从图2-(c)可见,地质灾害多集中发生于2080 m坡高范围内,且灾害点密度最高;随坡高增加,地质灾害密度和数量呈递增趋势,在5080 m范围内达最大。3.2.4 坡型斜坡坡型决定斜坡变形破坏方式及稳定性13。据野外实际调查情况,将研究区坡型分为凸型、直线型和凹型。在GIS中通过分析DEM斜坡地表曲率P基础数据DEM地质图地质环境数据数据来源中科图新科技有限公司收集野外调查数据数据类型12.5 m12.5 m栅格矢量矢量表1 研究所采用的数据来源和及其类型Table 1 The sources and types of data used in the study图2 评价因子地质灾害数
14、量、地质灾害点密度相关性统计图Fig.2 Correlation statistics of the number of geological hazards and thedensity of geological hazard points of the evaluation factors264第41卷第2期王伟中等:基于GIS和证据权模型的山阳县地质灾害易发性评价量化坡型。曲率P0时,为凸型坡;曲率P=0时,为直线型坡;曲率P0时,为凹型坡。本文将曲率介于-0.50.5的斜坡视作直线型坡,曲率大于0.5的斜坡视作凸型坡,曲率小于-0.5的斜坡视作凹形坡。通过GIS分析做相关统计得到图2
15、-(d),显示凸型和凹型斜坡更易产生滑坡灾害,地质灾害点密度相对较高。3.2.5 工程地质岩组不同类型岩石具不同软硬程度,岩体抗风化能力、应力不同,发生地质灾害的可能性与程度不同14。据区内出露地层工程地质特性,将研究区岩土体类型划为5类。从图2-(e)可见,薄层状较弱浅变质岩岩组(C)发育地质灾害数量最多,其地质灾害点密度最大;中厚层状坚硬碳酸盐岩组(B)发育地质灾害数量和地质灾害点密度位居其次。3.2.6 距断层距离断层造成岩土体破碎,断裂活动及对地下水控制作用直接影响斜坡稳定性,是地质灾害发育、发生的重要因素。在GIS中建立断层多环缓冲区,评价断裂构造对地质灾害发育的影响。从图2-(f)
16、看出:在距断层0100 m、5001 000 m范围内,地质灾害数量和地质灾害密度与距断层距离之间是负相关关系;在距断层100500 m、大于或等于1 000 m范围内,地质灾害数量和地质灾害密度与断层距离呈正相关关系,且在大于或等于1 000 m范围内发育的地质灾害数量最多,地质灾害密度最大。3.2.7 距河流距离河流对地质灾害的影响主要表现在侵蚀坡脚、软化岩土体、降低岩土体物理力学参数。在GIS中建立河流多环缓冲区,评价河流对地质灾害发育的影响。从图2-(g)看出,在距河流0200 m范围内,地质灾害数量和地质灾害密度与距河流距离之间是正相关关系;距河流2001000 m范围内,地质灾害数
17、量和地质灾害密度与河流距离存在负相关关系;大于或等于1 000 m的范围内发育地质灾害数量最多,地质灾害密度也最大。3.2.8 距道路距离道路对地质灾害的影响,主要表现为道路修建过程中开挖坡脚形成临空面,致使斜坡应力重分布,破坏斜坡原始稳定性,诱发地质灾害。在GIS中,建立主要道路多环缓冲区,用以评价道路对地质灾害发育的影响。从图2-(h)可见,在距道路0200 m范围内,地质灾害数量和密度与距道路距离之间。存在正相关关系;在距道路2001 000 m范围内,地质灾害数量和地质灾害密度与道路之间存在负相关关系;在大于或等于1 000 m范围内发育的地质灾害数量最多,密度却最低。3.2.9 地貌
18、类型地貌形态及变化是影响斜坡变形、失稳的重要因素。通过调查分析将研究区划分为河谷区、低山区、中山区3种地貌类型。通过GIS分析发现,地质灾害数量和密度与地貌类型呈负相关关系,海拔越高,分布的地质灾害数量越少,地质灾害密度越低;在低山区发育地质灾害数量最多,河谷区地质灾害密度最高。3.3 指标因子相关性分析为提高模型评价精度,本文采用相关系数量化各指标因子之间的相关关系。利用Arc GIS空间分析工具的多元分析模块对前文所述分级的指标因子进行波段值统计分析。通过计算其协方差与相关矩阵进行各指标因子间的相关性分析,并结合相关性划分表完成对各评价指标的独立性检验15。当各因子之间相关系数满足|R|0
19、.3时认为其弱相关或不相关。由表2,3,图3可见,地貌类型与距水系距离、距道路距离及坡高之间相关性较好,其余各评价指标指标因子工程地质岩组距河流距离坡向坡度距道路距离坡型坡高地貌类型距断层距离工程地质岩组1-0.098-0.004-0.184-0.1470.037-0.1710.2270.183距河流距离10.0110.0970.3120.0210.167-0.4440.151坡向10.0010.0090.012-0.0080.0190.009坡度10.1250.0140.118-0.134-0.002距道路距离10.0530.208-0.3510.066坡型10.156-0.163-0.05
20、5坡高1-0.7070.046地貌类型1-0.036距断层距离1表2 9个指标因子间的相关系数值Table 2 Value of correlation coefficient between 9 index factors265新疆地质2023年因子之间均满足相互独立要求,故将地貌类型指标因子剔除。3.4 证据层及后验概率计算据对上述各指标因子的分析统计,本文最终选取坡度、坡向等8个指标因子下共44个证据层,利用ArcGIS提取各指标因子不同状态分级下的栅格数据,计算各指标因子不同状态分级下的证据权重,进而计算各证据层的后验概率(表4)。4 地质灾害易发性评价及精度检验4.1 地质灾害易发性
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 GIS 证据 模型 山阳县 地质灾害 易发性 评价
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。