基于GIS的岷江上游汶川-叠溪河段滑坡灾害危险性精准评价与防控策略

基于GIS的岷江上游汶川—叠溪河段滑坡灾害危险性精准评价与防控策略一、引言1.1研究背景与意义1.1.1研究背景岷江上游汶川—叠溪河段地处青藏高原东缘，位于松潘—甘孜褶皱带、西秦岭造山带以及龙门山构造带的结合部位，地质构造复杂，新构造运动强烈。该区域山高谷深，地形起伏大，海拔高度变化显著，地势落差悬殊，河流深切，河谷狭窄。同时，岷江作为长江上游水量最丰富的支流之一，其上游水流湍急，河水的侵蚀与搬运作用强烈，进一步塑造了该地区复杂的地形地貌。此外，该区域属于亚热带季风气候区，干湿季分明，夏季降水集中且多暴雨，为滑坡灾害的发生提供了丰富的水源条件。在多种不利因素的综合影响下，岷江上游汶川—叠溪河段成为滑坡灾害的频发区。如1933年叠溪发生7.5级强烈地震，引发大规模山体滑坡，滑坡体堵塞岷江形成叠溪海子，导致河道断流，淹没了叠溪古城及周边村庄，造成了巨大的人员伤亡和财产损失。2017年6月24日，茂县叠溪镇新磨村新村组富贵山山体突发高位垮塌，约1800万立方米滑坡体将村庄掩埋，岷江支流松坪沟2公里河道堵塞，再次给当地带来了沉重的灾难。频繁发生的滑坡灾害不仅严重威胁着当地居民的生命财产安全，还对交通、水利等基础设施造成了极大的破坏，制约了区域的经济发展和社会稳定。面对如此严峻的滑坡灾害形势，开展滑坡灾害危险性评价显得尤为必要。通过对该区域滑坡灾害危险性进行科学、准确的评价，可以深入了解滑坡灾害的形成机制、分布规律和发展趋势，为制定有效的防灾减灾措施提供科学依据，从而降低滑坡灾害带来的损失。1.1.2研究意义从防灾减灾角度来看，准确的滑坡灾害危险性评价能够帮助相关部门提前识别出高危险区域，从而有针对性地制定防灾减灾规划。比如在高危险区域设置警示标识，加强监测预警，提前做好人员疏散和物资转移等准备工作，当滑坡灾害发生时，能够最大限度地减少人员伤亡和财产损失。保障人民生命财产安全是研究的重要目标。岷江上游汶川—叠溪河段居住着众多居民，滑坡灾害的威胁时刻存在。通过危险性评价，为居民提供安全的居住环境建议，引导居民合理选址建房，避免在滑坡易发区域建设，切实保障人民群众的生命财产安全。对于区域可持续发展而言，滑坡灾害危险性评价为区域的土地利用规划、基础设施建设等提供科学指导。在进行土地开发和建设时，充分考虑滑坡灾害的危险性，避开高危险区域，合理布局产业和基础设施，促进区域经济的可持续发展，实现人与自然的和谐共生。1.2国内外研究现状国外在基于GIS的滑坡灾害危险性评价方面开展了大量研究。早期，研究主要侧重于利用GIS的基本功能对滑坡相关数据进行管理和简单分析，如存储和展示地形、地质等数据。随着技术的不断发展，学者们开始运用各种数学模型与GIS相结合进行危险性评价。例如，利用概率模型分析滑坡发生的概率，通过将地形、地质、水文等多源数据输入模型，计算不同区域的滑坡发生概率，从而划分危险性等级。在欧洲阿尔卑斯山区，研究人员通过收集长时间序列的地形、气象等数据，运用统计分析方法与GIS技术相结合，对该区域的滑坡灾害进行了系统的危险性评价，为山区的防灾减灾和土地利用规划提供了重要参考。国内的研究也取得了丰硕成果。在数据收集与处理方面，充分利用高分辨率遥感影像、地理信息数据等，获取更准确的滑坡信息和影响因子数据。在评价模型方面，除了借鉴国外的模型外，还结合国内实际情况进行了改进和创新。例如，在三峡库区，研究人员运用层次分析法确定各评价因子的权重，再结合GIS的空间分析功能，对库区的滑坡灾害危险性进行评价，为三峡工程的安全运行和库区的生态保护提供了有力支持。同时，国内还注重多学科交叉融合，将地质学、气象学、地理学等学科知识与GIS技术相结合，提高评价的科学性和准确性。然而，当前研究仍存在一些不足与空白。在评价因子选取方面，部分研究对一些新出现的影响因素考虑不够全面，如人类工程活动对滑坡的影响日益显著，但在某些评价中未能充分体现。在模型构建方面，不同模型之间的比较和验证还不够深入，缺乏统一的标准来评估模型的优劣，导致在实际应用中难以选择最合适的模型。此外，对于滑坡灾害危险性的动态变化研究相对较少，不能及时反映由于环境变化、人类活动等因素导致的危险性变化情况。1.3研究内容与方法1.3.1研究内容本研究首先进行多源数据收集，涵盖地形数据，如数字高程模型（DEM），以获取该区域的地形起伏、坡度、坡向等信息；地质数据，包括地层岩性、地质构造等，用于分析地质条件对滑坡的影响；水文数据，如河流水系分布、降雨量等，明确水文因素与滑坡的关系；以及遥感影像数据，通过解译获取滑坡的分布范围和形态等信息，同时分析植被覆盖等情况。评价因子选取是关键环节，依据研究区域特点和已有研究成果，确定地形地貌因子（如坡度、坡向、高程、地形曲率等）、地质因子（地层岩性、断层距离等）、水文因子（河流水系距离、降雨量等）和植被因子（植被覆盖度）等作为评价因子，并对各因子进行量化和分级处理。构建并验证评价模型，选择层次分析法确定各评价因子的权重，结合信息量模型计算滑坡危险性指数。通过将研究区域划分为一定数量的评价单元，利用GIS的空间分析功能，计算每个评价单元的危险性指数，完成模型构建。然后，运用已知的滑坡样本数据对模型进行验证，评估模型的准确性和可靠性。进行危险性分区，根据计算得到的滑坡危险性指数，采用自然断点法等方法将研究区域划分为不同的危险性等级，如低危险性区、较低危险性区、中等危险性区、较高危险性区和高危险性区，并绘制滑坡灾害危险性分区图。对评价结果进行深入分析，研究不同危险性等级区域的分布规律，探讨各评价因子对滑坡危险性的影响程度，结合区域实际情况，分析滑坡灾害的形成机制和潜在风险，为防灾减灾提供针对性建议。1.3.2研究方法运用GIS空间分析方法，利用其强大的空间数据处理和分析能力，对收集到的多源数据进行处理和分析。通过叠加分析，将不同类型的图层（如地形图层、地质图层、水文图层等）进行叠加，获取综合信息，分析各因子之间的相互关系；利用缓冲区分析，确定河流、断层等因素的影响范围；通过坡度分析、坡向分析等获取地形地貌特征信息，为滑坡灾害危险性评价提供基础数据支持。采用层次分析法确定评价因子权重，该方法将与决策总是有关的元素分解成目标、准则、方案等层次，在此基础上进行定性和定量分析。首先构建层次结构模型，将滑坡灾害危险性评价目标作为最高层，评价因子作为中间层，评价单元作为最底层；然后通过专家打分等方式构造判断矩阵，计算各评价因子相对于目标层的相对权重，从而确定各因子对滑坡灾害危险性的影响程度。借助信息量模型进行滑坡危险性计算，该模型基于信息论原理，通过计算各评价因子在滑坡发生区域和非发生区域的信息量，来衡量因子与滑坡之间的相关性。信息量越大，表明该因子对滑坡发生的影响越大。通过计算每个评价单元内各评价因子的信息量，并进行累加，得到该评价单元的滑坡危险性指数，以此评估滑坡灾害危险性。二、研究区概况与数据来源2.1研究区概况2.1.1地理位置与地质构造岷江上游汶川—叠溪河段位于四川省西北部，地理坐标介于东经102°50′-103°40′，北纬31°30′-32°10′之间。该区域处于青藏高原东缘向四川盆地的过渡地带，是松潘—甘孜褶皱带、西秦岭造山带以及龙门山构造带的交汇部位，地质构造极为复杂。区域内断裂构造发育，主要有岷江断裂、虎牙断裂、龙门山断裂等。岷江断裂呈北西-南东向贯穿研究区，是控制该区域地质构造格局的重要断裂之一，其新构造运动强烈，历史上曾多次引发地震活动。虎牙断裂位于研究区东北部，同样具有较强的活动性，对区域内地层的变形和隆升起到了重要作用。龙门山断裂带则在研究区东南部，是一条规模巨大的逆冲推覆断裂带，2008年汶川8.0级特大地震就发生在龙门山断裂带的映秀-北川断裂上。这些断裂的活动导致研究区内地层褶皱、错动，岩石破碎，为滑坡灾害的发生提供了地质条件。断裂活动使得山体岩石的完整性遭到破坏，节理、裂隙发育，降低了岩石的抗剪强度，在重力、降雨、地震等因素的作用下，容易引发滑坡。如1933年叠溪7.5级地震就是由于岷江断裂的活动引发的，地震导致山体大规模滑坡，堵塞岷江形成叠溪海子。2.1.2地形地貌特征研究区属于高山峡谷地貌，地势西北高、东南低，海拔高度在1200-4500米之间，地势落差悬殊。山脉纵横交错，峡谷深切，河谷狭窄，坡度陡峭，大部分区域坡度在30°-60°之间。高山峡谷地貌的形成与区域地质构造运动和河流侵蚀作用密切相关。新构造运动导致山体隆升，而岷江及其支流的强烈下切侵蚀作用进一步塑造了峡谷地貌。这种地形地貌条件使得该区域成为滑坡灾害的高发区。陡峭的山坡在重力作用下，岩土体本身就处于不稳定状态，一旦受到外界因素的干扰，如地震、降雨等，就容易发生滑坡。峡谷地形使得滑坡体在滑动过程中能够获得较大的势能，加速滑动速度，增加滑坡的破坏力。2.1.3气象水文条件研究区属于亚热带季风气候区，受地形和季风的影响，气候垂直变化明显。年平均气温在10℃-15℃之间，气温随海拔升高而降低。降水主要集中在5-10月，约占全年降水量的80%以上，且多以暴雨形式出现，日最大降雨量可达100毫米以上。岷江是长江上游水量最丰富的支流之一，在研究区内蜿蜒而过，其支流众多，水系发达。河流水位受降水影响较大，雨季时河水暴涨，水流湍急，对河岸的侵蚀作用增强。河流的侵蚀作用会削弱河岸岩土体的稳定性，掏空坡脚，为滑坡的发生创造条件。降雨入渗会增加岩土体的重量，降低其抗剪强度，当降雨强度和持续时间达到一定程度时，就容易引发滑坡。2.1.4历史滑坡灾害情况历史上，岷江上游汶川—叠溪河段发生过多次严重的滑坡灾害。1933年8月25日，叠溪发生7.5级强烈地震，引发大规模山体滑坡，滑坡体堵塞岷江形成叠溪海子，淹没了叠溪古城及周边村庄，造成6800多人死亡，是该地区历史上最为严重的滑坡灾害之一。2008年5月12日，汶川发生8.0级特大地震，研究区受到强烈影响，引发了大量的滑坡、崩塌等次生地质灾害，对当地的基础设施和居民生命财产造成了极大的破坏。2017年6月24日，茂县叠溪镇新磨村新村组富贵山山体突发高位垮塌，约1800万立方米滑坡体将村庄掩埋，造成重大人员伤亡和财产损失。通过对历史滑坡灾害事件的分析可以发现，该区域滑坡灾害具有以下规律和特点：一是多与地震活动密切相关，强烈地震往往会引发大规模的滑坡灾害；二是集中发生在雨季，降雨是滑坡灾害的重要触发因素之一；三是滑坡规模较大，破坏力强，常常造成严重的人员伤亡和财产损失。2.2数据来源与预处理2.2.1数据来源地形数据主要来源于美国地质调查局（USGS）提供的30米分辨率的数字高程模型（DEM）数据，该数据能够准确反映研究区的地形起伏、坡度、坡向等信息，为地形地貌分析和滑坡灾害危险性评价提供了基础数据。地质数据包括地层岩性、地质构造等，主要通过收集四川省地质矿产勘查开发局编制的1:20万地质图件以及相关的地质调查报告获取。这些资料详细记录了研究区内地层的分布、岩性特征以及断裂构造的位置和性质等信息，对于分析地质条件对滑坡的影响至关重要。气象数据从中国气象数据网获取，包括研究区及周边气象站点多年的降水、气温、风速等数据。降水数据用于分析降雨对滑坡的触发作用，气温数据可辅助判断积雪融化等对滑坡的影响，风速数据在一定程度上也会影响滑坡的发生和发展。滑坡灾害数据通过实地调查、遥感影像解译以及相关部门的灾害统计资料收集。实地调查主要对研究区内已发生的滑坡灾害进行现场勘查，记录滑坡的位置、规模、形态、滑动方向等信息；遥感影像解译利用高分辨率卫星遥感影像，识别滑坡的分布范围和边界；相关部门的灾害统计资料则提供了历史滑坡灾害的时间、地点、造成的损失等信息，为滑坡灾害危险性评价提供了实际案例支持。2.2.2数据预处理对地形数据进行预处理，利用ArcGIS软件的空间分析工具对DEM数据进行拼接、裁剪，使其范围与研究区一致。通过填充洼地、计算坡度、坡向、地形曲率等操作，提取地形地貌特征信息，为后续的评价因子计算做准备。地质数据的预处理主要是对地质图件进行数字化处理，将纸质图件转化为矢量数据，方便在GIS软件中进行编辑和分析。对地层岩性数据进行分类和编码，统一数据格式，使其能够与其他数据进行叠加分析。对于地质构造数据，准确绘制断裂构造的位置和走向，并建立拓扑关系，以便分析其对滑坡的影响范围。气象数据的预处理包括数据的质量控制和格式转换。对气象数据进行异常值检查和修正，确保数据的准确性和可靠性。将气象数据从原始格式转换为GIS软件能够识别的格式，如将降水数据按照空间位置进行插值处理，生成栅格数据，以便与其他数据进行空间分析。滑坡灾害数据的预处理是对实地调查和遥感影像解译获取的数据进行整理和核实。对滑坡的位置信息进行精度验证，确保其准确性；对滑坡的规模、形态等数据进行标准化处理，使其具有可比性。将滑坡灾害数据与其他数据进行关联，如将滑坡位置与地形、地质数据进行叠加，分析滑坡与各影响因子之间的关系。三、基于GIS的滑坡灾害危险性评价指标体系构建3.1评价因子选取3.1.1地形因子地形因子是影响滑坡发生的重要因素之一，主要包括坡度、坡向、高程、地形起伏度等。坡度对滑坡的影响显著。当坡度较小时，岩土体在重力作用下相对稳定；随着坡度增大，重力沿坡面的分力逐渐增大，岩土体的下滑力增强，抗滑力相对减小，滑坡发生的可能性增加。研究表明，岷江上游汶川—叠溪河段大部分滑坡发生在坡度25°-60°的区域。通过对该区域历史滑坡数据与坡度数据的叠加分析发现，在坡度35°-45°的范围内，滑坡发生的频率最高。这是因为在这个坡度区间，岩土体既要承受自身重力，又要抵御外界因素的干扰，稳定性较差，一旦受到降雨、地震等触发因素的作用，就容易发生滑坡。坡向影响着太阳辐射、降水和风化作用的差异。不同坡向的岩土体干湿程度、温度变化和风化程度不同，进而影响其稳定性。在岷江上游地区，阳坡接受的太阳辐射较多，温度较高，岩土体风化作用较强，岩石破碎程度较大；阴坡相对湿度较大，岩土体含水量较高，抗剪强度降低。因此，阳坡和阴坡都有可能成为滑坡的易发区域。例如，在一些向阳的山坡，由于风化作用导致岩土体结构松散，在降雨条件下容易发生滑坡；而在背阴山坡，长期的高湿度环境使岩土体软化，也增加了滑坡的风险。高程与滑坡的关系较为复杂。随着高程的增加，地形起伏度增大，地质构造活动往往更加剧烈，岩石破碎程度增加，为滑坡的发生提供了物质基础。同时，高程的变化还会导致气候条件的差异，如气温降低、降水增多等，这些因素都可能影响滑坡的发生。在岷江上游汶川—叠溪河段，高海拔地区的滑坡发生率相对较高。这是因为高海拔地区的岩土体受到低温、冻融等作用的影响，结构更加不稳定，而且降水较多，容易引发滑坡。此外，高海拔地区的地形陡峭，一旦发生滑坡，滑坡体的势能较大，破坏力更强。地形起伏度反映了区域内地形的变化程度。地形起伏度大的区域，地势高差悬殊，岩土体在重力作用下的稳定性较差，容易发生滑坡。在该研究区域，高山峡谷地貌导致地形起伏度较大，是滑坡灾害频发的重要原因之一。如一些峡谷地段，两侧山体高差可达数百米，山坡陡峭，在地震、降雨等因素的作用下，极易发生大规模的滑坡。通过计算地形起伏度，并与滑坡分布数据进行对比分析，可以发现地形起伏度较大的区域，滑坡的发生概率明显增加。3.1.2地质因子地质因子是滑坡灾害发生的内在基础，主要包括岩土类型、地质构造、地震烈度等。岩土类型是滑坡发生的物质基础，不同类型的岩土体具有不同的物理力学性质，其抗剪强度、抗风化能力等差异显著，对滑坡的发生和发展有着重要影响。在岷江上游汶川—叠溪河段，广泛分布着砂泥岩互层、花岗岩、变质岩等岩土类型。其中，砂泥岩互层由于砂岩和泥岩的力学性质差异较大，泥岩遇水易软化，导致岩体整体抗剪强度降低，在外界因素作用下容易发生滑坡。花岗岩虽然质地坚硬，但经过长期的风化作用，表层形成了风化壳，风化壳中的松散岩土体在降雨等条件下也可能引发滑坡。变质岩由于其复杂的岩石结构和矿物成分，在构造应力和风化作用下，岩石破碎，节理裂隙发育，同样为滑坡的发生创造了条件。地质构造对滑坡的控制作用十分明显。断裂、褶皱等地质构造破坏了岩土体的完整性，使岩石节理、裂隙发育，降低了岩土体的抗剪强度。同时，地质构造还影响着地下水的运移和分布，进一步改变了岩土体的稳定性。研究区内的岷江断裂、虎牙断裂等活动断裂，是滑坡灾害的重要控制因素。断裂活动导致山体岩石破碎，形成了大量的破碎带和软弱结构面，这些区域在地震、降雨等触发因素作用下，极易发生滑坡。如1933年叠溪地震引发的大规模滑坡，就是由于岷江断裂的强烈活动，使得山体岩石破碎，在地震波的作用下，岩土体失稳滑动。地震烈度是衡量地震对地面影响和破坏程度的指标，与滑坡的发生密切相关。强烈地震产生的地震波会使山体岩土体受到强烈的震动和冲击，导致岩土体结构破坏，抗剪强度降低，从而引发滑坡。地震烈度越高，滑坡发生的可能性和规模越大。在岷江上游地区，历史上多次发生的强烈地震都引发了大量的滑坡灾害。2008年汶川8.0级特大地震，研究区地震烈度高达Ⅺ度，地震引发了大量的滑坡、崩塌等次生地质灾害，滑坡分布范围广，规模巨大，对当地的生态环境和基础设施造成了极大的破坏。通过对历史地震和滑坡数据的分析发现，地震烈度大于Ⅷ度的区域，滑坡发生的频率和规模明显增加，而且地震引发的滑坡往往具有突发性和群发性的特点，危害更为严重。3.1.3水文因子水文因子在滑坡灾害的发生发展过程中起着关键作用，主要包括河流分布、地下水位等。河流作为区域水文系统的重要组成部分，对滑坡的影响主要体现在两个方面。一方面，河流的侵蚀作用会削弱坡脚的稳定性。河水的长期冲刷使坡脚岩土体被掏空，导致上部岩土体失去支撑，在重力作用下容易发生滑坡。在岷江上游，河流蜿蜒曲折，水流湍急，对河岸的侵蚀作用强烈。一些河岸地段，由于长期受到河水的冲刷，坡脚逐渐变陡，岩土体松动，成为滑坡的高发区域。另一方面，河流的水位变化也会影响滑坡的发生。在雨季，河流水位迅速上升，浸泡坡脚岩土体，使其饱和软化，抗剪强度降低；当河流水位下降时，坡体内外产生水头差，形成动水压力，进一步加剧了坡体的不稳定。例如，在一些靠近河流的山坡，当河流水位上涨并长时间浸泡坡脚后，一旦遇到强降雨等触发因素，就容易引发滑坡。地下水位的变化对滑坡的影响不容忽视。地下水位上升时，岩土体含水量增加，重度增大，导致下滑力增大；同时，地下水对岩土体的软化作用也会降低其抗剪强度。此外，地下水位的变化还会产生孔隙水压力和动水压力，改变岩土体的应力状态，从而影响滑坡的稳定性。在岷江上游汶川—叠溪河段，由于地形起伏大，地下水的分布和运移较为复杂。在一些低洼地段或含水层丰富的区域，地下水位较高，岩土体长期处于饱水状态，稳定性较差。当遇到持续降雨或其他因素导致地下水位进一步上升时，就容易引发滑坡。通过对该区域地下水位监测数据与滑坡发生情况的对比分析发现，地下水位上升幅度较大的时期，滑坡发生的概率明显增加。3.1.4气象因子气象因子是滑坡灾害的重要触发因素，主要包括降雨、气温等。降雨是诱发滑坡的最主要气象因素之一。降雨对滑坡的触发机制主要体现在以下几个方面：一是雨水的入渗使岩土体含水量增加，重度增大，从而增加了下滑力；二是降雨入渗导致孔隙水压力升高，有效应力降低，抗剪强度减小；三是持续降雨或暴雨可能使岩土体达到饱和状态，形成超孔隙水压力，导致岩土体液化，进一步降低其稳定性。在岷江上游地区，降水主要集中在5-10月，且多以暴雨形式出现。大量的降雨数据和滑坡灾害记录表明，该地区大部分滑坡都发生在降雨过程中或雨后短时间内。如在一场暴雨过后，山坡上的岩土体迅速饱和，孔隙水压力急剧上升，抗剪强度大幅降低，导致滑坡的发生。通过对降雨强度、持续时间与滑坡发生关系的研究发现，当降雨强度超过一定阈值，且持续时间达到一定时长时，滑坡发生的概率会显著增加。气温对滑坡的影响主要通过积雪融化和风化作用来实现。在高海拔地区，冬季积雪较多，春季气温升高时，积雪融化产生大量的水流。这些水流渗入岩土体中，增加了岩土体的含水量和重量，同时也可能导致地下水位上升，从而引发滑坡。此外，气温的变化还会影响岩石的风化作用。昼夜温差和季节温差的变化使岩石发生热胀冷缩，加速岩石的风化破碎，降低岩土体的强度，为滑坡的发生创造条件。在岷江上游的高山区域，由于气温的垂直变化明显，高海拔地区的积雪融化和风化作用对滑坡的影响较为突出。3.1.5人为因子随着人类工程活动的日益频繁，人为因子对滑坡灾害的影响也越来越大，主要包括工程建设、植被破坏等。工程建设活动如道路修建、房屋建设、露天采矿等，往往会改变山体的原有地形地貌和岩土体结构，破坏山体的稳定性。在道路修建过程中，开挖坡脚、填方等工程活动会使坡体的应力状态发生改变，增加坡体的下滑力。如果在施工过程中没有采取有效的支护和排水措施，就容易引发滑坡。例如，一些山区公路在修建时，为了节省成本，随意开挖坡脚，导致坡体失衡，在降雨等因素的作用下，经常发生滑坡灾害，影响公路的正常通行。露天采矿活动会导致山体岩石被大量开采，形成采空区，使山体结构遭到严重破坏，岩土体稳定性降低，容易引发滑坡和崩塌等地质灾害。植被对山体具有重要的保护作用，它可以通过根系固土、截留雨水、减少地表径流等方式增强山体的稳定性。然而，人类的过度砍伐、开垦等活动导致植被破坏严重，削弱了植被对山体的保护功能。植被破坏后，地表裸露，水土流失加剧，岩土体失去植被根系的锚固作用，抗剪强度降低，在降雨等因素的作用下，容易发生滑坡。在岷江上游地区，由于部分地区的过度开垦和森林砍伐，植被覆盖率下降，山体的稳定性受到影响，滑坡灾害的发生率有所增加。通过对植被覆盖度与滑坡发生关系的研究发现，植被覆盖度较低的区域，滑坡发生的概率相对较高。3.2评价因子权重确定3.2.1层次分析法原理与步骤层次分析法（AnalyticHierarchyProcess，AHP）是一种将定性分析与定量分析相结合的多准则决策分析方法，由美国运筹学家萨蒂（T.L.Saaty）于20世纪70年代提出。其基本原理是将复杂的决策问题分解为若干个层次，每个层次包含若干个因素，通过对同一层次内各因素相对于上一层次中某一准则的重要性进行两两比较，构造判断矩阵，计算各因素的相对权重，从而确定各因素对目标的影响程度。运用层次分析法确定因子权重的步骤如下：建立层次结构模型：将滑坡灾害危险性评价问题分解为目标层、准则层和指标层。目标层为滑坡灾害危险性评价；准则层包括地形因子、地质因子、水文因子、气象因子和人为因子等；指标层则是各准则层下具体的评价因子，如坡度、岩土类型、河流分布等。构造判断矩阵：针对准则层中的每个准则，对其下一层的指标进行两两比较，判断它们对于该准则的相对重要性。采用1-9标度法来量化这种相对重要性，其中1表示两个因素同等重要，3表示一个因素比另一个因素稍微重要，5表示一个因素比另一个因素明显重要，7表示一个因素比另一个因素强烈重要，9表示一个因素比另一个因素极端重要，2、4、6、8则表示上述相邻判断的中间值。例如，对于地形因子准则下的坡度和坡向两个指标，若认为坡度比坡向稍微重要，则在判断矩阵中对应位置赋值为3。计算相对权重：通过计算判断矩阵的特征向量，得到各指标相对于上一层次中某一准则的相对权重。常用的方法有特征根法、和积法等。以特征根法为例，计算判断矩阵A的最大特征值λmax及其对应的特征向量W，将特征向量W归一化后得到各指标的相对权重。一致性检验：由于判断矩阵是基于专家主观判断构造的，可能存在不一致性。因此，需要进行一致性检验，以确保判断结果的合理性。一致性指标CI=(λmax-n)/(n-1)，其中n为判断矩阵的阶数。随机一致性指标RI可通过查表得到，一致性比例CR=CI/RI。当CR<0.1时，认为判断矩阵具有满意的一致性，否则需要对判断矩阵进行调整。3.2.2权重计算与一致性检验邀请地质、地理、灾害防治等领域的专家，根据层次分析法的原理和步骤，对各评价因子进行两两比较，构造判断矩阵。以地形因子准则下的坡度、坡向、高程、地形起伏度四个指标为例，构造的判断矩阵如下：A=\begin{pmatrix}1&3&5&7\\1/3&1&3&5\\1/5&1/3&1&3\\1/7&1/5&1/3&1\end{pmatrix}利用特征根法计算该判断矩阵的最大特征值λmax和特征向量W。通过计算得到λmax=4.12，特征向量W=(0.52,0.27,0.14,0.07)T。将特征向量W归一化后，得到各指标的相对权重为：坡度0.52、坡向0.27、高程0.14、地形起伏度0.07。计算一致性指标CI=(4.12-4)/(4-1)=0.04，随机一致性指标RI（4阶矩阵）=0.90，一致性比例CR=0.04/0.90=0.04<0.1，说明该判断矩阵具有满意的一致性。按照同样的方法，分别计算地质因子、水文因子、气象因子和人为因子准则下各指标的权重，并进行一致性检验。最终得到各评价因子的权重，如下表所示：准则层指标层权重地形因子坡度0.52坡向0.27高程0.14地形起伏度0.07地质因子岩土类型0.45地质构造0.35地震烈度0.20水文因子河流分布0.60地下水位0.40气象因子降雨0.70气温0.30人为因子工程建设0.65植被破坏0.35通过层次分析法确定的各评价因子权重，反映了不同因子对滑坡灾害危险性的影响程度。在后续的滑坡灾害危险性评价中，将根据这些权重对各评价因子进行综合分析，从而更准确地评估滑坡灾害的危险性。四、基于GIS的滑坡灾害危险性评价模型构建与应用4.1信息量模型原理与方法4.1.1信息量模型基本原理信息量模型是基于信息论发展起来的一种用于评估滑坡灾害危险性的有效方法，其核心在于通过量化各影响因子与滑坡之间的相关性，来确定每个因子对滑坡发生所提供的信息量大小。在信息论中，信息量是衡量事件不确定性的一个重要指标，对于滑坡灾害而言，某一影响因子在滑坡发生区域和非发生区域的分布差异越大，其所包含的关于滑坡发生的信息量就越大。假设研究区内存在n个影响滑坡发生的因子，分别为X_1,X_2,\cdots,X_n，对于每个因子X_i，可以将其划分为m个状态X_{i1},X_{i2},\cdots,X_{im}。设滑坡发生区域为A，非发生区域为\overline{A}，则因子X_{ij}对滑坡发生的信息量I_{ij}可以通过以下公式计算：[I_{ij}=\ln\frac[I_{ij}=\ln\fracI_{ij}=\ln\frac五、岷江上游汶川—叠溪河段滑坡灾害危险性评价结果分析5.1评价结果可视化5.1.1危险性等级划分通过层次分析法确定各评价因子权重，结合信息量模型计算得到研究区每个评价单元的滑坡危险性指数。为了直观地展示研究区滑坡灾害危险性的分布情况，采用自然断点法对滑坡危险性指数进行分级，将研究区滑坡灾害危险性划分为五个等级，分别为低危险性、较低危险性、中等危险性、较高危险性和高危险性。具体分级标准及对应的危险性指数范围如下表所示：危险性等级危险性指数范围低危险性[0,0.2)较低危险性[0.2,0.4)中等危险性[0.4,0.6)较高危险性[0.6,0.8)高危险性[0.8,1.0]低危险性区域表示该区域发生滑坡灾害的可能性较小，地质条件相对稳定，地形地貌较为平缓，岩土体的抗滑能力较强，受到外部触发因素影响的敏感度较低。较低危险性区域发生滑坡的可能性相对较低，但在特定条件下，如强降雨、地震等，仍有可能发生滑坡灾害，该区域的地形、地质条件存在一定的不稳定性因素。中等危险性区域处于中间状态，滑坡发生的可能性适中，地形、地质、水文等因素相互作用，使得该区域具有一定的滑坡潜在风险。较高危险性区域发生滑坡的可能性较大，地形坡度较陡，岩土体较为破碎，地质构造复杂，容易受到降雨、地震等因素的影响而引发滑坡。高危险性区域则是滑坡灾害极易发生的区域，地形地貌条件恶劣，地质构造活动强烈，岩土体稳定性极差，一旦遇到合适的触发条件，就可能发生大规模的滑坡灾害。5.1.2绘制危险性分区图利用GIS强大的制图功能，将研究区按照上述划分的五个危险性等级进行可视化表达，绘制出岷江上游汶川—叠溪河段滑坡灾害危险性分区图。在制图过程中，为了使不同危险性等级区域更加清晰直观，采用不同的颜色和符号来表示各个等级。例如，低危险性区域用浅蓝色表示，较低危险性区域用绿色表示，中等危险性区域用黄色表示，较高危险性区域用橙色表示，高危险性区域用红色表示。同时，在图上添加图例，对不同颜色和符号所代表的危险性等级进行说明，以便读者能够快速准确地理解图中信息。此外，为了增强地图的可读性和实用性，还在图上标注了主要的地理要素，如河流、山脉、城镇、道路等，这些地理要素与滑坡灾害危险性分区相结合，可以更全面地反映研究区的地理环境特征以及滑坡灾害危险性与人类活动、基础设施之间的关系。通过滑坡灾害危险性分区图，可以直观地看出研究区内不同区域滑坡灾害危险性的高低分布情况，为后续的结果分析和防灾减灾决策提供了直观的依据。5.2结果分析与验证5.2.1各等级区域分布特征低危险性区域主要分布在研究区的东南部和西北部的部分地区。在东南部，地势相对较为平缓，地形起伏度较小，坡度大多在20°以下，且岩土类型多为较为稳定的花岗岩和厚层砂岩，这些岩石抗风化能力强，结构致密，为岩土体提供了较好的稳定性。此外，该区域距离主要断裂构造较远，地质构造活动相对较弱，地震烈度较低，对滑坡的触发作用较小。西北部的低危险性区域主要位于一些山间盆地和河谷平原地区，这些地方地势平坦，地下水水位较深，水文条件对滑坡的影响较小，同时植被覆盖度较高，植被根系对岩土体起到了很好的锚固作用，进一步增强了山体的稳定性。较低危险性区域在研究区内呈斑块状分布，主要集中在低危险性区域的周边以及一些坡度在20°-30°之间的山坡地带。这些区域的地形条件相对低危险性区域稍差，但仍处于相对稳定的状态。岩土类型多为砂泥岩互层，但砂岩和泥岩的比例相对较为均衡，泥岩的软化作用对整体稳定性的影响相对较小。此外，该区域受到河流侵蚀和地质构造活动的影响也相对较弱，虽然存在一定的滑坡潜在风险，但在正常情况下发生滑坡的概率较低。中等危险性区域主要分布在研究区的中部和东北部，这些区域的地形坡度一般在30°-45°之间，地形起伏度较大，岩土体在重力作用下的稳定性相对较差。岩土类型以砂泥岩互层为主，且泥岩含量相对较高，遇水容易软化，导致岩土体抗剪强度降低。同时，该区域距离主要断裂构造较近，地质构造活动较为频繁，岩石节理、裂隙发育，增加了滑坡发生的可能性。此外，河流对该区域的侵蚀作用也较为明显，部分河岸地段坡脚被掏空，进一步削弱了山体的稳定性。较高危险性区域集中分布在岷江及其主要支流的两岸以及一些地势陡峭的山区。岷江及其支流的长期侵蚀作用使得河岸坡脚岩土体稳定性降低，河水的涨落还会导致岩土体含水量和应力状态的变化，增加了滑坡的发生风险。在地势陡峭的山区，坡度大多在45°-60°之间，岩土体在重力作用下极易失稳滑动。这些区域的岩土类型复杂，岩石破碎程度高，地质构造活动强烈，地震烈度较高，对滑坡的触发作用显著。此外，人类工程活动如道路修建、采矿等在这些区域也较为频繁，进一步破坏了山体的稳定性，加剧了滑坡灾害的危险性。高危险性区域主要位于历史上发生过强烈地震的区域以及一些大型断裂构造附近。如1933年叠溪地震和2008年汶川地震的震中及周边地区，这些区域在地震的强烈作用下，山体岩石破碎，形成了大量的松散堆积物，为滑坡的发生提供了丰富的物质基础。同时，地震还导致山体内部结构破坏，节理、裂隙大量发育，降低了岩土体的抗剪强度。大型断裂构造附近的区域，由于断裂活动频繁，岩石破碎，应力集中，也是滑坡灾害的高发区。此外，高危险性区域的地形条件极为恶劣，坡度陡峭，地形起伏度大，降雨容易形成地表径流，快速汇聚的水流对山体的冲刷作用强烈，进一步增加了滑坡发生的可能性。5.2.2与历史滑坡灾害的对比验证将滑坡灾害危险性评价结果与研究区历史滑坡灾害数据进行对比分析，以验证评价结果的准确性和可靠性。历史滑坡灾害数据包括1933年叠溪地震、2008年汶川地震以及其他年份发生的滑坡灾害的位置、规模、发生时间等信息。通过对比发现，历史上发生的滑坡灾害大多集中在危险性评价结果中的较高危险性和高危险性区域。在1933年叠溪地震引发的滑坡灾害中，滑坡主要分布在岷江断裂附近以及叠溪海子周边地区，这些区域在本次危险性评价中均被划分为高危险性区域。2008年汶川地震引发的滑坡灾害同样主要发生在龙门山断裂带附近以及地形陡峭、岩石破碎的区域，与评价结果中的较高危险性和高危险性区域高度吻合。在其他年份发生的滑坡灾害中，也有超过80%的滑坡点位于较高危险性和高危险性区域。这表明本次基于GIS的滑坡灾害危险性评价结果能够较好地反映研究区实际的滑坡灾害分布情况，评价结果具有较高的准确性和可靠性。同时，对于一些历史上发生滑坡灾害但危险性评价等级相对较低的区域，进一步分析发现，这些区域可能在评价时某些影响因子的数据存在一定误差，或者在滑坡发生时受到了特殊的触发因素影响，如局部暴雨强度过大、人类工程活动的突然扰动等。总体而言，通过与历史滑坡灾害数据的对比验证，证明了本次评价结果能够为研究区的滑坡灾害防治提供科学有效的依据。5.2.3敏感性分析为了确定各评价因子对滑坡灾害危险性评价结果的影响程度，对选取的地形、地质、水文、气象和人为等五大类评价因子进行敏感性分析。敏感性分析采用单因素变化法，即固定其他评价因子不变，仅改变某一个评价因子的值，观察滑坡灾害危险性指数的变化情况。在地形因子中，坡度对滑坡灾害危险性的影响最为敏感。当坡度从20°增加到40°时，滑坡危险性指数显著增大，尤其是在坡度30°-40°之间，危险性指数增长幅度最大。这表明坡度的变化对滑坡发生的可能性影响很大，随着坡度的增大，岩土体在重力作用下的稳定性急剧下降，滑坡灾害危险性迅速增加。坡向的敏感性相对较低，但不同坡向对滑坡危险性仍有一定影响。阳坡和阴坡由于光照、湿度等条件的差异，岩土体的物理力学性质有所不同，从而影响滑坡的发生。高程和地形起伏度的敏感性相对适中，随着高程的增加和地形起伏度的增大，滑坡灾害危险性也呈现出逐渐增加的趋势，但增长幅度不如坡度明显。地质因子中，岩土类型和地质构造对滑坡灾害危险性的影响较为敏感。不同岩土类型的抗剪强度和抗风化能力差异显著，如砂泥岩互层的岩土体在遇水后容易软化，导致滑坡危险性增加。地质构造活动强烈的区域，岩石破碎，节理、裂隙发育，滑坡危险性明显增大。地震烈度的敏感性相对较低，但在高地震烈度区域，地震对滑坡的触发作用仍然不可忽视，随着地震烈度的增加，滑坡灾害危险性也会相应增加。水文因子中，河流分布对滑坡灾害危险性的影响较为敏感。靠近河流的区域，由于河水的侵蚀作用和水位变化的影响，滑坡危险性较高。当河流距离评价单元从500米减小到100米时，滑坡危险性指数明显增大。地下水位的敏感性相对较低，但地下水位的上升会增加岩土体的含水量和重度，降低其抗剪强度，从而对滑坡灾害危险性产生一定影响。气象因子中，降雨对滑坡灾害危险性的影响最为敏感。随着降雨量的增加和降雨强度的增大，滑坡危险性指数迅速增大。当降雨强度从10毫米/小时增加到50毫米/小时时，滑坡危险性指数增加了数倍。这表明降雨是触发滑坡灾害的关键因素之一，强降雨极易引发滑坡。气温的敏感性相对较低，主要通过积雪融化和风化作用间接影响滑坡灾害危险性。人为因子中，工程建设对滑坡灾害危险性的影响较为敏感。不合理的工程建设活动，如开挖坡脚、填方等，会改变山体的应力状态，增加滑坡危险性。当工程建设活动强度增加时，滑坡危险性指数明显增大。植被破坏的敏感性相对较低，但植被覆盖度的降低会削弱植被对山体的保护作用，增加滑坡发生的风险。通过敏感性分析可知，在岷江上游汶川—叠溪河段滑坡灾害危险性评价中，坡度、岩土类型、地质构造、河流分布和降雨等因子对评价结果的影响较为显著，是影响滑坡灾害危险性的关键因子。在今后的滑坡灾害防治工作中，应重点关注这些关键因子的变化，加强对高敏感性区域的监测和防范。六、滑坡灾害防治对策与建议6.1工程防治措施6.1.1滑坡治理工程抗滑桩是滑坡治理中常用的工程措施之一，其工作原理是通过桩身将上部承受的坡体推力传给桩下部的侧向土体或岩体，依靠桩下部的侧向阻力来承担边坡的下推力，使边坡保持平衡或稳定。抗滑桩类型多样，按材料可分为木桩、钢板桩、钢筋混凝土桩等，其中钢筋混凝土桩因其强度高、耐久性好等优点应用最为广泛；按截面形式可分为圆形、方形、矩形等。在岷江上游汶川—叠溪河段，对于浅层滑坡且滑坡推力较小的情况，可选用木桩或小型钢筋混凝土桩；而对于深层滑坡和滑坡推力较大的区域，则需采用大直径、高强度的钢筋混凝土抗滑桩。如在某滑坡治理工程中，根据滑坡的规模和推力计算，选用了直径1.5米、桩长20米的钢筋混凝土抗滑桩，桩间距为3米，呈梅花形布置。施工时，先进行桩孔开挖，然后放置钢筋笼，最后浇筑混凝土，施工完成后，有效阻止了滑坡的进一步滑动，保障了周边区域的安全。挡土墙也是一种有效的滑坡治理手段，它通过抵抗土体的侧向压力，防止土体滑坡或坍塌。常见的挡土墙类型有重力式挡土墙、悬臂式挡土墙、扶壁式挡土墙等。重力式挡土墙依靠自身重力来维持稳定，结构简单，施工方便，但对地基承载力要求较高，适用于高度较低、地基条件较好的滑坡地段；悬臂式挡土墙利用悬臂结构来抵抗土体压力，对地基承载力要求相对较低，适用于土质较差的区域；扶壁式挡土墙则在悬臂式挡土墙的基础上增加了扶壁，提高了墙体的抗弯性能，适用于高度较大的滑坡地段。在该研究区域，对于一些小型滑坡且地基承载力较好的地段，可采用重力式挡土墙，用块石或混凝土砌筑而成，墙高一般在3-5米，墙体厚度根据滑坡推力计算确定。对于中型滑坡且地基条件一般的区域，可选用悬臂式挡土墙，墙高可达8米左右，通过合理设计悬臂长度和墙体配筋来满足抗滑要求。而对于大型滑坡和高陡边坡，常采用扶壁式挡土墙，墙高可超过10米，扶壁间距一般为2-3米，以增强挡土墙的稳定性。排水工程对于滑坡治理至关重要，它能有效降低孔隙水压力和动水压力，防止岩土体的软化及溶蚀分解，消除或减小水的冲刷和浪击作用。排水工程包括地表排水和地下排水。地表排水主要通过修建截水沟、排水沟等设施，拦截和排除地表水，防止其渗入滑坡体。截水沟应设置在滑坡体边界以外一定距离处，以拦截周边区域的地表水，沟底和沟壁应进行防渗处理，防止渗漏。排水沟则应根据地形和水流方向合理布置，将地表水引至安全地带。地下排水可采用水平钻孔疏干、垂直孔排水、竖井抽水、隧洞疏干、支撑盲沟等方法。在岷江上游地区，由于降雨量大且集中，地表排水尤为重要。对于山坡较陡的滑坡体，可在滑坡体后缘和两侧设置截水沟，沟深一般为0.5-1米，沟宽0.4-0.6米，采用浆砌片石或混凝土浇筑而成。在滑坡体表面，每隔一定距离设置横向排水沟，将地表水引入截水沟。对于地下水位较高的区域，可采用垂直孔排水的方法，在滑坡体内钻孔，孔深根据地下水位和滑坡体厚度确定，一般为5-10米，孔内安装排水管，将地下水引出滑坡体。6.1.2监测预警系统建设利用现代技术建立滑坡灾害监测预警系统，对于及时发现滑坡灾害迹象、提前发出预警、减少灾害损失具有重要意义。监测系统主要通过部署各种监测设备，实时获取滑坡体的变形、位移、地下水位、降雨量等数据。位移监测可采用全球定位系统（GPS）、全站仪等设备，能够高精度地测量滑坡体的水平和垂直位移。在岷江上游汶川—叠溪河段的一些重要滑坡隐患点，设置了多个GPS监测点，通过卫星信号实时传输位移数据，监测精度可达毫米级。当滑坡体位移超过设定的预警阈值时，系统会及时发出警报。地下水位监测通常使用地下水观测井和水位传感器，能够准确掌握地下水位的变化情况。在滑坡体周边和内部设置地下水观测井，井深根据地层情况确定，一般为10-20米，通过水位传感器将地下水位数据实时传输到监测中心。如某滑坡隐患点，通过地下水位监测发现，在雨季时地下水位迅速上升，当水位超过警戒水位时，滑坡发生的风险显著增加。降雨量监测利用雨量传感器来实现，能够实时监测降雨强度和降雨量。在研究区内分布多个雨量监测站点，雨量传感器将采集到的降雨数据通过无线传输方式发送到监测中心。通过对降雨量数据的分析，结合历史滑坡灾害与降雨的关系，确定降雨触发滑坡的临界值，当降雨量达到或超过该临界值时，及时发出预警。预警系统则基于监测数据，通过建立科学的预警模型，对滑坡灾害的发生进行预测和预警。预警模型可采用经验模型、数值模型或机器学习模型等。经验模型是根据历史滑坡灾害数据和相关影响因素建立的简单预测模型，如基于降雨强度和持续时间的滑坡预警经验公式。数值模型则通过模拟滑坡体的力学行为和变形过程，预测滑坡的发生和发展。机器学习模型利用大量的历史数据进行训练，自动学习滑坡灾害的特征和规律，从而实现对滑坡灾害的准确预测。当监测数据达到预警阈值时，预警系统通过短信、广播、警报器等多种方式及时向当地居民和相关部门发出预警信息，以便采取相应的应急措施。如在某滑坡预警案例中，监测系统实时监测到滑坡体位移和地下水位异常变化，预警系统通过分析判断，提前2小时向周边居民发送短信预警，并启动广播和警报器进行通知。当地政府接到预警信息后，迅速组织居民疏散，避免了人员伤亡和财产损失。通过建立完善的滑坡灾害监测预警系统，能够实现对滑坡灾害的实时监测、准确预警和有效应对，为保障人民生命财产安全和区域可持续发展提供有力支持。6.2非工程防治措施6.2.1土地利用规划调整根据滑坡灾害危险性评价结果，对研究区的土地利用规划进行合理调整，是预防滑坡灾害的重要非工程措施。在高危险性区域，应严格限制人类工程活动，避免进行大规模的建设项目。对于已有的建筑物和基础设施，应进行安全评估，如存在安全隐患，应采取搬迁或加固等措施。例如，在岷江上游汶川—叠溪河段的高危险性区域，有一些居民点和小型工厂，这些区域的滑坡灾害风险极高。通过危险性评价后，当地政府制定了搬迁计划，将居民点搬迁至安全区域，并对工厂进行关闭或搬迁，有效降低了滑坡灾害对人员和财产的威胁。对于中等危险性区域，应加强土地利用的管控，合理规划建设项目。在进行建设时，应充分考虑滑坡灾害的影响，采取相应的工程措施进行防护。如在某中等危险性区域进行道路建设时，设计单位根据滑坡灾害危险性评价结果，对道路线路进行了优化，避开了滑坡风险较高的地段。同时，在道路建设过程中，采用了抗滑挡土墙、排水等工程措施，确保道路的安全稳定。在低危险性和较低危险性区域，虽然滑坡灾害发生的可能性相对较小，但仍需进行合理的土地利用规划。鼓励发展生态农业、林业等对土地扰动较小的产业，保护植被，提高山体的稳定性。如在低危险性区域，推广种植经济林和生态林，既增加了农民的收入，又改善了生态环境，增强了山体的抗滑能力。通过科学合理的土地利用规划调整，能够从源头上减少滑坡灾害的发生风险，实现土地资源的合理利用和区域的可持续发展。6.2.2宣传教育与应急管理加强滑坡灾害宣传教育，提高公众的防灾减灾意识和自救互救能力，是减少滑坡