![点击分享此内容可以赚币 分享](/master/images/share_but.png)
基于信息量、加权信息量与逻辑回归耦合模型的云南罗平县崩滑灾害易发性评价对比分析.pdf
《基于信息量、加权信息量与逻辑回归耦合模型的云南罗平县崩滑灾害易发性评价对比分析.pdf》由会员分享,可在线阅读,更多相关《基于信息量、加权信息量与逻辑回归耦合模型的云南罗平县崩滑灾害易发性评价对比分析.pdf(11页珍藏版)》请在咨信网上搜索。
1、DOI:10.16031/ki.issn.1003-8035.202208030杨得虎,朱杰勇,刘帅,等.基于信息量、加权信息量与逻辑回归耦合模型的云南罗平县崩滑灾害易发性评价对比分析J.中国地质灾害与防治学报,2023,34(5):43-53.YANG Dehu,ZHU Jieyong,LIU Shuai,et al.Comparative analyses of susceptibility assessment for landslide disasters based oninformation value,weighted information value and logistic
2、 regression coupled model in Luoping County,Yunnan ProvinceJ.TheChineseJournalofGeologicalHazardandControl,2023,34(5):43-53.基于信息量、加权信息量与逻辑回归耦合模型的云南罗平县崩滑灾害易发性评价对比分析杨得虎1,朱杰勇1,2,刘帅1,马博1,代旭升2(1.昆明理工大学国土资源工程学院,云南昆明650093;2.自然资源部高原山地地质灾害预报预警与生态保护修复重点实验室,云南昆明650093)摘要:以罗平县崩滑地质灾害为研究对象,选取工程岩组、坡度、坡向、高程、起伏度、曲率
3、、地貌类型、距河流距离、距断裂距离 9 个评价因子,基于共线性诊断和相关性分析对其进行独立性检验。然后采用信息量法计算各评价因子分类分级的信息量值,采用层次分析法和逻辑回归法对各评价因子进行权重的定量计算,从而构建信息量、加权信息量和信息量-逻辑回归耦合易发性评价模型并进行对比分析。基于 GIS 的自然断点法将评价结果划分为非、低、中和高 4 个等级,并采用 ROC 曲线对其精度进行检验。结果表明:3 种评价模型的 AUC 值分别为 0.757、0.723 和 0.852,信息量-逻辑回归耦合模型的精度最高,模型结果分区与崩滑地质灾害点的分布较吻合,其非、低、中和高的面积(分级比)分别为 77
4、1.1km2(25.55%)、836.6km2(27.73%)、864.36km2(28.64%)和 545.94km2(18.08%)。关键词:崩塌;滑坡;地质灾害;信息量模型;加权信息量模型;信息量-逻辑回归耦合模型;易发性评价中图分类号:P694文献标志码:A文章编号:1003-8035(2023)05-0043-11Comparative analyses of susceptibility assessment for landslidedisasters based on information value,weighted information valueand logisti
5、c regression coupled model in Luoping County,Yunnan ProvinceYANGDehu1,ZHUJieyong1,2,LIUShuai1,MABo1,DAIXusheng2(1.Faculty of Land Resource Engineering,Kunming University of Science and Technology,Kunming,Yunnan650093,China;2.Key Laboratory of Geohazard Forecast and Geoecological Restoration in Plate
6、au Mountainous Area,MNR,Kunming,Yunnan650093,China)Abstract:ThisstudyfocusesonlandslidesusceptibilityassessmentsinLuopingCounty,where9evaluationfactors,includingengineeringrockgroup,slope,slopeaspect,elevation,undulation,curvature,landformtype,distancefromrivers,anddistancefrom fault,were selected a
7、s the research variables.After conducting collinearity diagnosis and correlation analysis,the收稿日期:2022-08-22;修订日期:2022-10-14投稿网址:https:/ model;weighted information model;information-logisticregressioncoupledmodel;susceptibilityevaluation0 引言云南以红河深大断裂为界,滇西为高山纵谷区,滇东为喀斯特高原1,滇东北为地质灾害高易发区,滇西地质灾害易发性南相对较弱2
8、。罗平县地处云南滇东喀斯特高原地区,是典型的岩溶山区。目前,崩滑易发性制图主要分为定性、定量和机器学习 3 种方法34。定性分析是通过对成因机制的全面认识,基于专家经验和知识确定评价因子权重,定量分析方法通过数学或数值算法估计滑坡易发性5。定性分析方法主要有层次分析法68,定量分析方法主要有频率比法910、信息量法1112,机器学习方法有逻辑回归法、随机森林法、K 近邻、支持向量机和神经网络等1315。国内外学者进行研究分析,并采用多种研究方法进行对比,吉日伍呷等15通过逻辑回归、K 近邻、朴素贝叶斯和随机森林算法对鲁甸进行地震滑坡易发性评价,并得出统计建模地更多的是寻找变量之间的可解释关系;
9、樊芷吟等16通过信息量法+Logistic 回归模型对汶川县进行易发性评价;张晓东17通过定量信息量法和确定性系数法分别与 Logistic 回归耦合对宁夏盐池县进行易发性评价,得出耦合模型结果均优于单一模型评价结果。根据以上研究,易发性评价模型首先选取崩滑评价影响因子,再通过分类分级进行计算,根据不同方法得出的结果基于 ArcGIS 进行等级划分。信息量法只考虑了评价因子分类分级状态下的权重,其优点在于原理简单,易于实现;缺点在于无法体现所选取的评价因子的权重。因此选取层次分析法和逻辑回归法分别赋予评价因子的权重,其中层次分析法依据崩滑灾害成因分析,通过专家法构建矩阵计算得出评价因子的权重,
10、逻辑回归法是依据样本数据和连接方法,通过两种方法得出评价因子的权重分别与信息量法进行耦合,其耦合模型将评价因子的权重和分类分级的权重叠加得到综合权重,降低了单一评价模型人为主观性因素的影响,论证不同方法评价结果的准确性。本文以罗平县作为研究区域,基于野外详实的地质灾害调查成果,综合分析孕灾地质条件和崩滑点分布规律,选取岩土体(工程岩组)、地形地貌(坡度、坡向、高程、起伏度、曲率、地貌类型)、地质构造(距断裂距离)、气象水文(距河流距离)等评价因子。采用信息量法、加权信息量法、信息量-逻辑回归耦合法构建易发性评价模型进行对比分析,并对评价结果进行精度检验分析,选取精度最高模型易发性分布图,可为罗
11、平县今后地质灾害治理提供参考依据,对今后城市的发展和防灾减灾有重要意义,也可为岩溶地区地质灾害易发性评价提供参考。1 研究区概况及数据来源1.1研究区概况研究区(罗平县)位于云南省曲靖市东部。东西最大横距 75km,南北最大纵距 99km。相对高差为1705m,全县面积 3018km2,山区面积占 78%,坝区面积占 22%(图 1)。研究区西部和北部属于岩溶盆地地貌和岩溶低中山地貌,中部属岩溶断陷湖形盆地,东部和南部受九龙河和南盘江流域侵蚀切割,形成峰林洼地和岩溶中山地貌。区内地层出露主要有古生界泥盆系(D)浅灰、深灰色中厚层状灰岩、泥灰岩、泥质白云岩;石炭系(C)深灰、灰黑色块状灰岩、白云
12、质灰岩、泥质灰岩;古生界二叠系(P)灰、深灰色厚层块状、生物碎屑灰岩,结晶灰岩夹虎斑状灰岩及白云岩;中生界三叠系(T)上统为黄褐色粉砂岩、泥质粉砂岩及细砂岩、中统为深灰色灰岩夹泥质灰岩、中上部为黄色白云岩、下统为紫红色含长石粉细砂夹泥灰岩页岩及含铜44中国地质灾害与防治学报第5期页岩;新生界古近系(E)+新近系(N)褐黄紫红色砾岩、细砂岩及粉砂质泥岩、底部砾岩;新生界第四系(Q)细砂、砂砾石及砂质黏土。主要构造体系和构造型式有北东向构造、新华夏系构造、网状构造等。北东向构造为区内主导构造,是研究区内最重要的构造成分之一,主要断裂有:金鸡山断裂、长家湾断裂和腊庄断裂等。其次为新华夏系构造,多发育
13、在褶皱边缘、密集成束、规模大、延伸远、呈舒缓波状,主要分布在西部及南盘江两岸。主要断裂有:洒土革断裂、大水塘断裂、罗格断裂等。010 km5崩塌滑坡河流罗平县边界高程/m断裂2 420715N图 1 研究区概况Fig.1 Overview of the study area1.2数据来源本研究数据主要包括:(1)12.5m 分辨率数字高程模型(DEM)收集自 ASF,用于提取坡度、坡向、起伏度、曲率等评价因子;(2)120 万地质图收集自全国地质资料馆,用于提取岩性、断裂等因子;(3)15 万地理数据库提取水系;(4)历史崩滑数据:主要来自地矿眉山工程勘察院 15 万全区调查结果,共 154
14、个崩滑灾害点的数据。2 研究方法2.1信息量模型信息量模型是对崩滑历史数据进行统计分析,将影响崩滑的各因子的实测值转化为信息量值,来衡量崩滑的易发性18。首先计算各评价因子的信息量值,再对各因子信息量值进行总和,作为崩滑易发性的综合指标19。单因子信息量计算公式为:I=lnNj/NSj/S(1)式中:I评价因子 j 下的信息量;Nj评价因子 j 内发生的崩滑数;N研究区崩滑总数;Sj评价因子 j 下所占栅格数;S研究区总栅格数。将每个评价单元各分类分级进行叠加计算,其地质灾害发生的总信息量计算公式为:Ij=ni=1lnNj/NSj/S(2)式中:Ij总信息量,为地质灾害易发性指数;Ij值越大且
15、为正值则表示该单元内有利于崩滑发生。2.2层次分析模型层次分析模型是一种将决策者定性判断和定量计算有效结合起来的分析方法。通过比较相邻影响因子的重要性19,根据专家法构建判断矩阵20:A=(aij)=a11a1n.an1ann(3)式中:A要素判断矩阵;aij因子 i 和因子 j 重要性比较的结果,有以下性质:aij=1aji,aij,0,aii=1,i,j=1,2,n(4)为保证求得的权重的正确性及合理性,还需要进行一致性检验。CI=maxnn1(5)CR=CI/RI(6)式中:CI一致性指标;n判断矩阵的阶数;max判断矩阵的最大特征值;CR随机一致性比;当其0.1 时一致性检验通过;RI
16、随机一致性指标。2023 年杨得虎,等:基于信息量、加权信息量与逻辑回归耦合模型的云南罗平县崩滑灾害易发性评价对比分析452.3逻辑回归模型逻辑回归模型是一种研究二分类因变量常用的统计方法16,21。通过研究崩滑易发性与评价因子之间的关系,预测崩滑发生的概率。其中自变量为评价因子指标值(x1,x2,xn),是否发生地质灾害作为因变量(分别用 1 和 0 代表崩滑点和非崩滑点)。逻辑回归函数如下:Z=+1x1+2x2+nxn(7)P(y=1)=11+ez(8)式中:常数项;x1,x2,xn自变量;1,2,n回归系数;Z崩滑发生的可能性与各评价因子之间的关系;P崩滑灾害发生的概率,范围 01。2.
17、4加权信息量模型根据层次分析法得出各评价因子的权重值,结合信息量法各评价因子分类分级的信息量值,两者相乘得出加权信息量值,其计算公式可表示为:Ij=ni=1iIi=ni=1ilnNi/NSi/S(9)式中:Ij加权信息量;i每个评价因子的权重;Ii评价因子 i 的信息量值。2.5信息量-逻辑回归耦合模型将信息量模型与逻辑回归模型进行耦合,通过逻辑回归确定评价因子的权重,可降低信息量模型评价因子分级的主观性影响。其原理将信息量模型中评价因子分类分级的信息量值作为逻辑回归模型中的自变量,建立回归方程进行逻辑回归运算,得出各评价因子的回归系数,以此为依据建立信息量-逻辑回归耦合模型。3 评价因子选取
18、分级3.1评价因子选取分级本文在罗平县资料收集和野外地质调查的基础上,选取岩土体(工程岩组)、地形地貌(坡度、坡向、高程、起伏度、曲率、地貌类型)、地质构造(距断裂距离)、气象水文(距河流距离)等评价因子进行分析。根据 12.5m12.5m 栅格单元作为易发性评价的制图单元,通过对研究区评价因子与崩滑点数据进行归纳分析,得出各评价因子的分类分级处理(图 2)。3.2评价因子共线性诊断进行逻辑回归时,需确保所选评价因子之间的相互独立,相关性高会出现多重共线性2223。采用容忍度(tolerance,TOL)和方差膨胀因子(varianceinflationfactor,VIF)对自变量进行多重共
19、线性诊断:VIF=11R2(10)式中:R2以 xi为因变量时对其他自变量回归的复测 定系数;TOL 是 VIF 的倒数,当 TOL 大于 0.1 且 VIF 小于10 时,则不存在多重共线性。根据 308 个独立属性样本,提取每个样本的各类级信息量值,在 SPSS 软件中进行多重共线性诊断。结果显示对所选 9 个评价因子其 VIF 值在 11.5(表 1)。其VIF5,表明各因子之间相互独立,不存在共线性。3.3评价因子相关性分析崩滑的易发性与评价因子之间存在一定的相关性。为了保证各评价因子间的相互独立性和结果的可靠性,进行因子相关性检验24。结果显示各评价因子之间的相关系数均0.3(表 2
20、),评价因子之间的相关性较小,所以 9 个评价因子均可以进入模型。4 易发性评价结果4.1信息量模型评价结果信息量模型中,崩滑的易发性与因子信息量值有关,信息量值越大且为正值则表示单元内崩滑越容易发生2528。根据已有 154 个地质灾害进行重分类统计,根据公式(1)计算各评价因子分类分级的信息量值(表 3)。4.2加权信息量模型评价结果根据加权信息量法构建模型,对研究区地质灾害及其背景因素和影响因素的相对重要性进行分析,依据所选取的评价因子,按照专家法对选取的评价因子根据式(3)构建判断矩阵计算出各评价因子的权重值,根据式(5)求出每个判断矩阵的一致性指标 CI,并通过式(6)进行一致性检验
21、(表 4),各评价因子的权重值与各评价因子分类分级的信息量值根据式(9)得出加权信息量值(表 3)。4.3信息量-逻辑回归耦合模型评价结果根据研究区已有的 154 个崩滑点,并随机选取等量的非崩滑点,共计有 308 个独立属性样本。将全部样本点依次赋予相应评价因子的信息量值,导入 SPSS25 软件进行二项逻辑回归分析(表 5),各评价因子的信息量46中国地质灾害与防治学报第5期值作为自变量,是否发生地质灾害作为因变量(1 和 0代表崩滑点和非崩滑点)。在逻辑回归分析结果中,sig 值越小,代表评价因子的显著性越高,表(3)中 sig 小于 0.05,说明 9 个因子均有统计意义。基于模型分析
22、结果中的各因子系数值根据公式(7)得逻辑回归公式如下:0 5高程/m7158608601 2001 2001 3501 3501 5001 5001 6501 6501 8001 8001 9501 9502 4200 5地貌类型侵蚀谷地岩溶中山岩溶低中山地岩溶盆地峰丛洼地峰林谷地断块上升岩溶断坳盆地构造侵蚀剥蚀构造侵蚀岩溶石丘(垅岗)0 5起伏度/m04488151523233030383850502200 5曲率0000 5距断裂距离/m06006001 2001 2001 8001 8002 4002 4003 0003 0000 5距河流距离/m06006001 2001 2001 8
23、001 8002 4002 4003 0003 000010 km5工程岩组软硬相间碳酸盐岩夹碎屑岩岩组块状结构坚硬玄武岩岩组坚硬层状碳酸盐岩岩组第四系冲洪积松散岩组010 km5坡度()066121218182424303036366060900 5坡向平面北东北东东南南西南西西北(a)工程岩组(b)坡度(c)坡向(d)高程(e)地貌类型(f)起伏度(g)曲率(h)距断裂距离(i)距河流距离10 km10 km10 km10 km10 km10 km10 kmNNNNNNNNN图 2 崩滑评价因子分类分级图Fig.2 Classification map of landslide susce
24、ptibility evaluation factors2023 年杨得虎,等:基于信息量、加权信息量与逻辑回归耦合模型的云南罗平县崩滑灾害易发性评价对比分析47Z=0.165+0.698x1+0.031x2+0.761x3+0.309x4+0.171x5+0.641x6+1.523x7+0.528x8+0.458x9P=11+ez(11)式中:x1x9分别为地层岩组、坡度、坡向、高程、地貌类型、起伏度、曲率、距断裂距离、距河流距离的信息量值;运用 ArcGIS 的栅格计算器功能将 z 值代入式(8)得到崩滑灾害发生的概率 p。信息量模型根据信息量法求出各评价因子分类分级的信息量,然后进行叠加
25、分析,加权信息量模型根据表 1 评价因子VIF计算结果表Table 1 Calculation results of VIF for evaluation factors评价因子TOLVIF工程岩组0.8181.222坡度0.6561.524坡向0.9541.048高程0.9041.107地貌类型0.7131.402起伏度0.6691.495曲率0.9701.031距断裂距离0.9451.058距河流距离0.7171.396表 2 评价因子之间的相关系数矩阵Table 2 Correlation coefficient matrix of evaluation factors评价因子工程岩组坡
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 信息量 加权 逻辑 回归 耦合 模型 云南 罗平县 灾害 易发性 评价 对比 分析
![提示](https://www.zixin.com.cn/images/bang_tan.gif)
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。