云南省施甸县旧城乡集镇泥石流风险评估与防治策略研究

云南省施甸县旧城乡集镇泥石流风险评估与防治策略研究一、引言1.1研究背景与意义泥石流作为一种极具破坏力的地质灾害，常常在短时间内释放出巨大的能量，给人类生命财产和基础设施带来严重威胁。它是指在山区或者其他沟谷深壑、地形险峻的地区，因为暴雨、暴雪或其他自然灾害引发山体滑坡，并携带有大量泥沙及石块的特殊洪流，具有突发性、流速快、流量大、物质容量大和破坏力强等特点。据统计，我国是世界上泥石流灾害最为严重的国家之一，近十年来，平均每年造成的直接或间接经济损失达100多亿元，死亡近千人，严重制约社会经济发展。云南省施甸县旧城乡地处山区，特殊的地理位置与复杂的地质条件，使其深受泥石流灾害的困扰。旧城乡属于勐波罗河大断裂带上与怒江断裂带以北，区内分布地质灾害隐患点14个，其中大型泥石流冲沟4条、中型3条，小型滑坡7个，威胁人口865人，威胁资产约3720万元。每年雨季，强降雨天气频繁，加上地形地貌复杂，山体坡度较大，岩土体稳定性差，为泥石流的发生创造了极为有利的条件，致使该区域泥石流灾害频发。这些灾害不仅冲毁了大量的农田、房屋、道路等基础设施，给当地居民的生命财产安全带来了巨大损失，也对当地的生态环境造成了严重破坏，阻碍了区域经济的可持续发展。对旧城乡集镇进行泥石流风险评价具有重大的现实意义。精确的风险评估能够为当地政府及相关部门提供科学的决策依据，助力合理规划土地利用，从源头上降低泥石流灾害的潜在威胁。通过明确不同区域的风险等级，可指导在低风险区域进行合理开发建设，而在高风险区域避免不必要的投资，实现资源的优化配置。同时，风险评估结果有助于制定科学有效的防灾减灾措施，如加强监测预警系统建设，及时发布泥石流灾害预警信息，使居民能够提前做好防范准备，最大程度减少人员伤亡；优化工程治理方案，对泥石流沟道进行加固、疏导等工程处理，增强其抵御灾害的能力；开展防灾减灾宣传教育，提升居民的自我保护意识和应对灾害的能力，从而保障社会的稳定和可持续发展，促进旧城乡的经济繁荣与人民生活的安宁。1.2国内外研究现状在泥石流风险评估领域，国外起步相对较早。20世纪60年代起，欧美、日本等国家和地区就开始关注泥石流灾害，并开展了一系列研究。早期主要侧重于泥石流的形成机制与运动规律探索，如日本高桥堡在1991年编著的《泥石流》，系统阐述了泥石流的基本特征、运动机理、发生发展及堆积过程，为后续风险评估研究奠定了理论基础。随着研究的深入，逐渐从定性描述向定量评估转变。奥地利、瑞士等欧洲国家较早采用类似交通信号灯红、黄、绿三色含义，划分泥石流危险区、潜在危险区和无危险区，初步实现了风险区域的划分。在风险评估方法上，国外不断创新。基于概率统计的方法被广泛应用，通过对历史泥石流事件的发生频率、规模等数据统计分析，预测未来灾害发生的可能性和危害程度。数值模拟技术也得到了快速发展，如FLO-2D、RAMMS等模型，能够模拟泥石流的运动过程，包括流速、流量、泥深、堆积范围等参数，为风险评估提供了更精确的数据支持。近年来，机器学习、深度学习等人工智能技术逐渐应用于泥石流风险评估，通过对大量数据的学习和分析，建立更准确的风险评估模型，提高评估的精度和效率。国内对泥石流风险评估的研究始于20世纪80年代。1986年，谭炳炎提出了泥石流沟严重程度的数量化综合评判方法，开启了我国泥石流风险定量评估的先河。此后，刘希林、唐川等学者在泥石流危险性分析评价方面做了大量工作。1988年，刘希林首次提出泥石流危险度的判定方法；1993年，唐川、刘希林等应用泥石流二维非恒定流理论建立危险度评价的数学模型，并初步应用于实际。随着地理信息系统（GIS）、全球定位系统（GPS）、遥感（RS）等技术在我国的普及，泥石流风险评估得到了极大的推动。这些技术能够快速获取地形地貌、地质构造、土地利用等多源数据，并进行空间分析和处理，直观展示泥石流风险的空间分布特征，为风险评估提供了强大的数据支持和分析手段。尽管国内外在泥石流风险评估方面取得了丰硕成果，但仍存在一些不足。现有评估模型在复杂地形和地质条件下的适应性有待提高，部分模型对数据的要求较高，而实际应用中数据获取往往存在困难，导致模型的推广和应用受到限制。不同评估方法之间的对比和验证研究相对较少，缺乏统一的评估标准和规范，使得评估结果的可比性和可靠性受到影响。对泥石流灾害的不确定性研究还不够深入，难以准确评估灾害发生的概率和损失程度，无法满足实际防灾减灾工作的需求。针对以上不足，本研究以云南省施甸县旧城乡集镇为研究区域，综合考虑当地复杂的地形地貌、地质条件、气象因素以及人类活动等影响，运用多源数据和多种评估方法，构建适用于该区域的泥石流风险评估模型，旨在提高风险评估的准确性和可靠性，为当地泥石流灾害的防治提供科学依据。1.3研究内容与方法本研究针对云南省施甸县旧城乡集镇的泥石流风险评估，主要围绕以下几个方面展开研究内容。1.3.1研究内容地质条件分析：全面调查旧城乡集镇区域的地质构造特征，包括地层岩性、断裂构造、褶皱等情况，深入分析这些地质因素对泥石流形成的控制作用。研究发现，旧城乡属于勐波罗河大断裂带上与怒江断裂带以北，区内分布地质灾害隐患点14个，其中大型泥石流冲沟4条、中型3条，小型滑坡7个，地质构造复杂，为泥石流的发生提供了潜在的地质条件。对区域内岩土体的物理力学性质进行详细测试，如岩土体的抗剪强度、渗透性、容重等参数，评估岩土体在强降雨等条件下的稳定性，明确易引发泥石流的岩土体类型和分布区域。地形地貌特征研究：利用高分辨率的遥感影像和地形数据，精确提取旧城乡集镇的地形地貌信息，如坡度、坡向、沟谷密度、相对高差等。通过分析这些地形地貌因素与泥石流发生的相关性，确定有利于泥石流形成和运动的地形条件。研究表明，坡度较大、沟谷纵横的区域更容易发生泥石流灾害。绘制详细的地形地貌图，直观展示区域内地形地貌的空间分布特征，为后续的泥石流风险评估提供基础数据。气象因素分析：收集旧城乡集镇多年的气象数据，包括降雨量、降雨强度、降雨历时、气温、风速等，分析气象因素的时空变化规律，重点研究强降雨事件的发生频率、强度和持续时间，以及其与泥石流发生的时间关系。通过统计分析发现，每年雨季（6-9月）是泥石流的高发期，且降雨量和降雨强度与泥石流的发生密切相关。结合气象数据和历史泥石流灾害记录，建立气象因素与泥石流发生的统计模型，预测在不同气象条件下泥石流发生的可能性。人类活动影响评估：调查旧城乡集镇的人类工程活动，如道路建设、房屋建筑、采矿活动、土地开垦等，分析这些活动对地表植被、地形地貌和岩土体稳定性的破坏程度，评估人类活动对泥石流形成和发展的影响。研究发现，不合理的工程建设和土地开垦活动，如在山坡上随意开挖、破坏植被等，会增加泥石流发生的风险。评估人类活动对泥石流风险的影响程度，提出合理的人类活动规划和管理建议，以减少因人类活动引发的泥石流灾害。泥石流灾害历史分析：系统收集旧城乡集镇历史上发生的泥石流灾害资料，包括灾害发生的时间、地点、规模、造成的损失等信息，整理分析泥石流灾害的发生规律和发展趋势，为风险评估提供历史依据。通过对历史灾害的分析，可以总结出泥石流灾害的发生频率、危害程度等特征，为制定防灾减灾措施提供参考。利用历史灾害数据，对风险评估模型进行验证和校准，提高模型的准确性和可靠性。泥石流风险评估模型构建：综合考虑地质、地形地貌、气象和人类活动等因素，选取合适的评估指标，构建适用于旧城乡集镇的泥石流风险评估模型。模型的构建将采用层次分析法（AHP）、模糊综合评价法等方法，确定各评估指标的权重，实现对泥石流风险的定量评估。利用构建的风险评估模型，对旧城乡集镇不同区域的泥石流风险进行计算和评价，划分风险等级，绘制泥石流风险分布图，直观展示风险的空间分布特征。风险评估结果分析与应用：对泥石流风险评估结果进行深入分析，明确不同风险等级区域的分布范围和特征，分析高风险区域的形成原因和潜在危害。根据风险评估结果，为旧城乡集镇的土地利用规划、防灾减灾措施制定、应急救援预案编制等提供科学依据。例如，在高风险区域，应避免进行大规模的建设活动，加强监测预警和工程治理措施；在中低风险区域，也应加强风险管理，制定相应的防范措施。对风险评估结果的应用效果进行跟踪和评估，不断完善风险评估模型和防灾减灾措施，提高应对泥石流灾害的能力。1.3.2研究方法资料收集与整理：广泛收集旧城乡集镇的地质、地形地貌、气象、水文、社会经济等相关资料，包括地质勘察报告、地形图、气象数据、遥感影像、土地利用规划图等，对收集到的资料进行系统整理和分析，筛选出与泥石流风险评估相关的信息，建立基础数据库。野外调查与勘查：组织专业人员对旧城乡集镇进行野外实地调查，详细了解区域内的地质构造、地形地貌、岩土体特征、植被覆盖情况、人类工程活动以及泥石流灾害遗迹等。通过实地勘查，获取第一手资料，补充和验证资料收集阶段的数据，为后续的分析和评估提供可靠依据。在调查过程中，采用GPS定位技术、地质罗盘测量、现场采样等方法，对重要的地质和地形地貌特征进行准确记录和分析。数据分析与统计：运用统计学方法对收集到的气象、地质、地形地貌等数据进行分析，计算相关参数的平均值、标准差、变异系数等，分析数据的分布特征和变化规律。通过相关性分析、回归分析等方法，研究不同因素之间的相互关系，确定影响泥石流发生的关键因素。利用地理信息系统（GIS）技术对空间数据进行处理和分析，如坡度分析、坡向分析、地形起伏度分析等，提取与泥石流风险相关的地形地貌信息，直观展示数据的空间分布特征。数值模拟与模型构建：选用合适的泥石流数值模拟软件，如FLO-2D、RAMMS等，对旧城乡集镇的泥石流运动过程进行数值模拟。通过输入地形地貌、地质条件、气象数据等参数，模拟泥石流在不同工况下的流速、流量、泥深、堆积范围等参数，预测泥石流可能造成的危害范围和程度。基于层次分析法（AHP）、模糊综合评价法等方法，构建泥石流风险评估模型。通过专家打分、问卷调查等方式确定各评估指标的权重，利用模糊数学理论对风险进行综合评价，实现对泥石流风险的定量评估。对比分析与验证：将不同方法得到的泥石流风险评估结果进行对比分析，评估不同方法的优缺点和适用性，选择最适合旧城乡集镇的风险评估方法和结果。利用历史泥石流灾害数据对风险评估模型进行验证，通过比较模型预测结果与实际灾害情况，检验模型的准确性和可靠性。对模型进行调整和优化，提高模型的精度和适应性。综合评估与制图：综合考虑地质、地形地貌、气象、人类活动等因素，以及数值模拟和模型计算的结果，对旧城乡集镇的泥石流风险进行全面评估，划分风险等级，确定高、中、低风险区域的范围和边界。利用GIS技术绘制泥石流风险分布图，图中直观展示不同风险等级区域的空间分布情况，以及与地形地貌、土地利用、人口分布等因素的关系，为风险管理和决策提供可视化依据。1.4技术路线本研究的技术路线旨在系统、全面地对云南省施甸县旧城乡集镇的泥石流风险进行评估，为防灾减灾提供科学依据。具体流程如下：数据收集与整理：通过多种渠道广泛收集旧城乡集镇的地质、地形地貌、气象、水文、社会经济等相关资料，如地质勘察报告、高分辨率遥感影像、地形图、多年气象数据、土地利用规划图以及历史泥石流灾害记录等。对收集到的资料进行系统整理，筛选出与泥石流风险评估相关的数据，建立基础数据库，为后续分析提供数据支持。野外调查与勘查：组织专业的地质、地理等领域的技术人员，对旧城乡集镇进行全面细致的野外实地调查。实地观察区域内的地质构造特征，包括地层岩性的分布、断裂构造的走向和规模、褶皱的形态等；详细记录地形地貌特征，如山坡的坡度、坡向、沟谷的密度和形态、相对高差等；实地勘查岩土体的类型、结构和物理力学性质；查看植被的覆盖类型、覆盖程度及其分布状况；调查人类工程活动的类型、规模和分布，如道路建设、房屋建筑、采矿活动、土地开垦等；仔细观察泥石流灾害遗迹，包括泥石流的堆积物特征、冲蚀痕迹、灾害影响范围等。在调查过程中，运用GPS定位技术准确记录调查点的地理位置，使用地质罗盘测量地质构造要素，采集岩土体样本用于实验室测试分析，确保获取的数据真实可靠，为后续研究提供第一手资料。因素分析与指标选取：对收集到的数据和野外调查结果进行深入分析，研究地质条件、地形地貌、气象因素、人类活动等对泥石流形成和发展的影响机制。基于分析结果，选取能够准确反映泥石流风险的评估指标，如地层岩性、断裂密度、坡度、地形起伏度、年平均降雨量、降雨强度、植被覆盖率、工程活动强度等，为构建风险评估模型奠定基础。风险评估模型构建与计算：采用层次分析法（AHP）、模糊综合评价法等方法，确定各评估指标的权重，构建适用于旧城乡集镇的泥石流风险评估模型。利用构建的模型，输入相关指标数据，对旧城乡集镇不同区域的泥石流风险进行计算，得到每个区域的风险值。数值模拟与结果验证：选用FLO-2D、RAMMS等泥石流数值模拟软件，输入地形地貌、地质条件、气象数据等参数，对旧城乡集镇的泥石流运动过程进行数值模拟，得到泥石流在不同工况下的流速、流量、泥深、堆积范围等参数。将数值模拟结果与风险评估模型计算结果进行对比分析，相互验证，提高评估结果的准确性和可靠性。风险等级划分与制图：根据风险评估模型计算结果和数值模拟结果，采用自然断点法、等间距法等方法，对旧城乡集镇的泥石流风险进行等级划分，划分为高风险、中风险、低风险等不同等级。利用地理信息系统（GIS）技术，将不同风险等级区域在地图上进行可视化表达，绘制泥石流风险分布图，直观展示风险的空间分布特征。结果分析与防治策略提出：对泥石流风险评估结果进行深入分析，研究不同风险等级区域的分布规律、形成原因和潜在危害。针对不同风险等级区域，结合当地实际情况，提出科学合理的泥石流防治策略，包括工程治理措施，如修建挡土墙、拦砂坝、排导槽等；生物治理措施，如植树造林、恢复植被等；监测预警措施，如建立泥石流监测站点、完善预警系统等；应急管理措施，如制定应急预案、开展应急演练等，为旧城乡集镇的泥石流灾害防治提供科学指导。具体技术路线流程如图1.1所示：@startumlstart:收集地质、地形地貌、气象、水文、社会经济等资料，建立基础数据库;:进行野外调查与勘查，获取第一手资料;:分析地质、地形地貌、气象、人类活动等因素对泥石流的影响，选取评估指标;:采用AHP、模糊综合评价法等构建风险评估模型，计算风险值;:选用FLO-2D、RAMMS等软件进行数值模拟，验证评估结果;:采用自然断点法、等间距法等划分风险等级，绘制风险分布图;:分析风险评估结果，提出泥石流防治策略;stop@endumlstart:收集地质、地形地貌、气象、水文、社会经济等资料，建立基础数据库;:进行野外调查与勘查，获取第一手资料;:分析地质、地形地貌、气象、人类活动等因素对泥石流的影响，选取评估指标;:采用AHP、模糊综合评价法等构建风险评估模型，计算风险值;:选用FLO-2D、RAMMS等软件进行数值模拟，验证评估结果;:采用自然断点法、等间距法等划分风险等级，绘制风险分布图;:分析风险评估结果，提出泥石流防治策略;stop@enduml:收集地质、地形地貌、气象、水文、社会经济等资料，建立基础数据库;:进行野外调查与勘查，获取第一手资料;:分析地质、地形地貌、气象、人类活动等因素对泥石流的影响，选取评估指标;:采用AHP、模糊综合评价法等构建风险评估模型，计算风险值;:选用FLO-2D、RAMMS等软件进行数值模拟，验证评估结果;:采用自然断点法、等间距法等划分风险等级，绘制风险分布图;:分析风险评估结果，提出泥石流防治策略;stop@enduml:进行野外调查与勘查，获取第一手资料;:分析地质、地形地貌、气象、人类活动等因素对泥石流的影响，选取评估指标;:采用AHP、模糊综合评价法等构建风险评估模型，计算风险值;:选用FLO-2D、RAMMS等软件进行数值模拟，验证评估结果;:采用自然断点法、等间距法等划分风险等级，绘制风险分布图;:分析风险评估结果，提出泥石流防治策略;stop@enduml:分析地质、地形地貌、气象、人类活动等因素对泥石流的影响，选取评估指标;:采用AHP、模糊综合评价法等构建风险评估模型，计算风险值;:选用FLO-2D、RAMMS等软件进行数值模拟，验证评估结果;:采用自然断点法、等间距法等划分风险等级，绘制风险分布图;:分析风险评估结果，提出泥石流防治策略;stop@enduml:采用AHP、模糊综合评价法等构建风险评估模型，计算风险值;:选用FLO-2D、RAMMS等软件进行数值模拟，验证评估结果;:采用自然断点法、等间距法等划分风险等级，绘制风险分布图;:分析风险评估结果，提出泥石流防治策略;stop@enduml:选用FLO-2D、RAMMS等软件进行数值模拟，验证评估结果;:采用自然断点法、等间距法等划分风险等级，绘制风险分布图;:分析风险评估结果，提出泥石流防治策略;stop@enduml:采用自然断点法、等间距法等划分风险等级，绘制风险分布图;:分析风险评估结果，提出泥石流防治策略;stop@enduml:分析风险评估结果，提出泥石流防治策略;stop@endumlstop@enduml@enduml图1.1技术路线流程图二、旧城乡集镇概况及泥石流灾害现状2.1自然地理概况2.1.1地理位置旧城乡位于施甸县南部，地处东经99°05′05″—99°19′34″，北纬24°10′24″—24°30′33″，距县城60公里。其地理位置独特，呈“S”形走向，与多个地区接壤，东北与昌宁县和姚关镇相邻，东南与永德县隔勐波罗河相望，南与龙陵县隔怒江对峙，西与酒房乡相连，形成了与两市、四县相交的特殊区位格局。这种特殊的地理位置，使得旧城乡成为区域交通和经济交流的重要节点。勐波罗河和怒江两条重要水系流经境内，不仅为当地的农业灌溉、居民生活用水提供了保障，也在一定程度上影响了区域的交通布局和经济发展模式。然而，特殊的地理位置也使旧城乡面临着诸多挑战。由于地处山区，地形复杂，周边山脉纵横，河流深切，地质条件不稳定，为泥石流等地质灾害的发生创造了条件。同时，受季风气候影响，降水集中且强度大，在雨季时，山区容易遭受暴雨侵袭，进一步增加了泥石流灾害的发生风险。此外，与多个地区接壤的地理位置，使得人员流动和经济活动频繁，人类工程活动对地表植被和地形地貌的破坏也较为严重，加剧了地质灾害的潜在威胁。2.1.2地形地貌旧城乡地形地貌复杂多样，以山地和河谷为主。辖区内最高海拔达2595.2米，最低海拔仅583.88米，相对高差较大，地势起伏显著。这种高差大的地形条件，使得地表径流速度快，水流侵蚀能力强，为泥石流的形成提供了强大的动力条件。当强降雨发生时，短时间内大量的降水迅速汇聚成地表径流，在陡峭的山坡上急速下泄，容易携带大量的泥沙、石块等固体物质，从而引发泥石流灾害。山地面积广阔，坡度较陡，大部分地区坡度在25°以上，部分区域甚至超过45°。陡峭的山坡使得岩土体在重力作用下处于不稳定状态，容易发生滑坡和崩塌等地质灾害，这些灾害产生的松散物质又为泥石流的形成提供了丰富的物源。此外，山区的地形起伏还导致了沟谷发育，沟谷密度较大，且多呈“V”字形，这种沟谷形态有利于水流的汇聚和快速流动，进一步促进了泥石流的形成和发展。河谷地带地势相对较低平，是勐波罗河和怒江的流经区域。河谷两岸阶地发育，阶地沉积物厚度较大，主要由砂、砾石和黏土等组成。这些沉积物在长期的水流作用下，结构较为松散，抗冲刷能力较弱。在雨季，河水水位上涨，水流速度加快，对河谷两岸的阶地产生强烈的侵蚀作用，导致阶地沉积物被冲刷、搬运，为泥石流提供了物质来源。同时，河谷地区人口密集，经济活动频繁，人类工程活动如道路建设、房屋建筑等，破坏了原有的地形地貌和地表植被，进一步加剧了泥石流灾害的发生风险。2.1.3气候特征旧城乡属于典型的干热河谷气候，具有炎热少雨、干湿季分明、雨热同期、日照时间长、年积温高的特点。年平均气温为21.2℃，最高气温可达42℃，全年有效积温8358℃，热量资源丰富，有利于农作物的生长，农作物可一年三熟。然而，这种气候条件也对泥石流的发生有着重要影响。干湿季分明是旧城乡气候的显著特征之一。每年5-10月为雨季，降水集中，降雨量占全年的80%以上，且多以暴雨形式出现。据统计，多年平均降雨量为1084.6毫米左右，其中日最大降雨量可达150毫米以上。强降雨是引发泥石流的主要触发因素之一，短时间内大量的降雨使得地表径流迅速增加，土壤含水量饱和，岩土体的抗剪强度降低，容易导致山体滑坡和泥石流的发生。而11月至次年4月为干季，降水稀少，气候干燥，蒸发量大，导致地表植被生长受到限制，植被覆盖率较低，岩土体缺乏植被的保护，稳定性较差，也为雨季泥石流的发生埋下了隐患。日照时间长，年日照时数达1992小时，充足的光照使得河谷地区气温较高，加剧了水分的蒸发，进一步加重了干旱程度。在干热的气候条件下，岩土体风化作用强烈，岩石破碎，形成大量的松散碎屑物质，这些物质在强降雨作用下，容易被水流携带，形成泥石流。此外，干热河谷气候还导致了当地生态环境脆弱，植被恢复困难，一旦发生泥石流灾害，对生态环境的破坏将更为严重，恢复也更加困难。2.1.4地质条件旧城乡地质构造复杂，处于勐波罗河大断裂带上与怒江断裂带以北，区内断裂构造发育，主要有勐波罗河断裂、怒江断裂及其次级断裂等。这些断裂带的存在使得区域内地层破碎，岩石完整性遭到破坏，岩土体的力学性质变差，抗剪强度降低，为泥石流的形成提供了有利的地质条件。断裂活动还可能导致山体岩石发生错动、变形，产生大量的裂隙和破碎带，增加了山体的不稳定性，在强降雨等因素的作用下，容易引发山体滑坡和泥石流灾害。地层岩性主要有砂岩、页岩、泥岩、灰岩等。其中，砂岩和页岩的抗风化能力较弱，在长期的风化作用下，容易破碎形成松散的碎屑物质；泥岩遇水后容易软化、崩解，降低了岩土体的强度和稳定性；灰岩虽然硬度较大，但岩溶发育，溶洞、溶蚀裂隙等岩溶现象广泛存在，使得岩石的完整性受到破坏，也为泥石流的形成提供了物源。此外，不同岩性地层的组合和接触关系，也影响着岩土体的稳定性。例如，砂岩与页岩互层的地层，由于两者的力学性质差异较大，在地下水和外力作用下，容易发生层间滑动，引发滑坡和泥石流。区内岩土体类型多样，包括残积土、坡积土、洪积土等。残积土主要分布在山坡顶部和缓坡地带，厚度较薄，结构松散，颗粒大小不均，抗剪强度较低；坡积土多分布在山坡中下部，是由山坡上的岩土体在重力和水流作用下搬运堆积而成，其物质组成和结构较为复杂，稳定性较差；洪积土主要分布在河谷和沟口地带，是由洪水携带的泥沙、石块等物质堆积而成，颗粒较大，分选性差，透水性较强，但在饱水状态下，其强度会显著降低。这些岩土体在强降雨、地震等因素的作用下，容易发生变形、滑动和崩塌，为泥石流的形成提供了丰富的物质来源。2.2社会经济概况旧城乡全乡辖8个村（社区）委员会，81个村民小组，86个自然村，总人口4548户，16408人，有彝族、傣族、布朗族、傈僳族、白族、佤族等少数民族。人口分布呈现出明显的不均衡特征，河谷地带地势平坦，水源充足，交通便利，经济相对发达，吸引了大量人口聚居，人口密度较大；而山区地形复杂，交通不便，经济发展相对滞后，人口分布较为稀疏。在产业结构方面，旧城乡以农业为主要产业，充分利用当地干热河谷气候优势，发展特色农业。农作物总播种面积达100171亩，粮食作物以水稻、玉米为主，经济作物则以蔬菜、甘蔗、瓜果等作物为主。蔬菜种植面积12179亩，总产量191821吨，冬季蔬菜凭借独特的气候优势和上乘的品质，畅销全国各地。甘蔗种植也颇具规模，是当地的重要经济作物之一。此外，以牛、羊、猪、马、鸡、鸭、鹅饲养为主的畜牧业已形成规模，全乡肉牛养殖存栏6883头，出栏3065头，生猪存栏7.8679万头，出栏6.12万头，畜牧业在当地经济中占据重要地位。近年来，随着经济的发展和产业结构的调整，旧城乡的工业和服务业也逐渐兴起。但总体来说，工业基础较为薄弱，主要以农产品加工、小型制造业等为主，规模较小，技术水平相对较低，对劳动力的吸纳能力有限。服务业发展相对滞后，主要集中在餐饮、零售、交通运输等传统领域，新兴服务业态如电子商务、文化旅游、现代物流等发展缓慢，尚未形成规模效应。社会经济因素与泥石流灾害之间存在着复杂的相互作用。一方面，不合理的人类经济活动会加剧泥石流灾害的发生风险。随着人口的增长和经济的发展，人类对土地资源的需求不断增加，山区的过度开垦、滥砍滥伐现象较为严重，导致地表植被遭到破坏，水土流失加剧，为泥石流的发生提供了丰富的物源。在山区进行的道路建设、房屋建筑、采矿等工程活动，往往缺乏科学规划和有效的防护措施，破坏了山体的稳定性，容易引发山体滑坡和泥石流灾害。另一方面，泥石流灾害的发生又会对当地的社会经济发展造成严重的破坏。泥石流冲毁农田、房屋、道路、桥梁等基础设施，导致农业减产、居民财产受损、交通中断，给当地的经济发展带来巨大损失。灾害还会威胁到居民的生命安全，导致人员伤亡，影响社会的稳定和和谐。因此，在旧城乡的社会经济发展过程中，必须充分考虑泥石流灾害的影响，采取科学合理的措施，加强灾害防治，实现社会经济的可持续发展。2.3泥石流灾害现状2.3.1历史泥石流灾害事件回顾历史上，旧城乡发生过多起泥石流灾害，给当地带来了严重的破坏。2005年8月，旧城乡遭遇持续强降雨，降雨量在短时间内超过200毫米，引发了大规模的泥石流灾害。此次泥石流主要发生在勐波罗河沿岸的多个村庄，如旧城村、芒别村等。泥石流流速极快，携带大量泥沙、石块和树木，形成一股巨大的洪流，以排山倒海之势冲向下游。据统计，此次灾害导致20余间房屋被冲毁，大量农田被淹没，农作物受灾面积达500余亩，直接经济损失超过100万元。更令人痛心的是，造成了3人死亡，5人失踪，给当地居民的生命财产安全带来了巨大损失。2012年7月，旧城乡再次遭受泥石流灾害的侵袭。由于前期降雨较多，土壤含水量饱和，在一场暴雨的触发下，山区多条沟谷同时爆发泥石流。其中，大山村附近的泥石流规模较大，冲毁了村内的部分道路和桥梁，导致交通中断，村民出行受阻。此次灾害还对当地的水利设施造成了严重破坏，部分灌溉渠道被堵塞，影响了农田的灌溉，导致农作物减产。据不完全统计，此次灾害造成直接经济损失约80万元，受灾人口达200余人。2018年9月，受台风“山竹”外围云系影响，旧城乡出现强降雨天气，引发了多处泥石流灾害。在一些坡度较陡、植被覆盖率较低的区域，泥石流灾害尤为严重。如大田坝村的一条泥石流沟，在强降雨作用下，泥石流迅速形成，冲毁了沟口的几户民房，导致居民财产受损。灾害还造成了部分电力和通信线路中断，影响了当地居民的正常生活。此次泥石流灾害共造成直接经济损失约60万元，受灾人口150余人。这些历史泥石流灾害事件表明，旧城乡的泥石流灾害具有突发性强、破坏力大、损失严重等特点。每次灾害的发生，不仅给当地居民的生命财产安全带来了巨大威胁，也对当地的基础设施、农业生产和生态环境造成了严重破坏，制约了当地社会经济的发展。通过对这些历史灾害事件的回顾和分析，可以总结经验教训，为今后的泥石流灾害防治工作提供参考依据。2.3.2现有泥石流隐患点分布根据施甸县自然资源局的调查资料，旧城乡境内目前分布有14个地质灾害隐患点，其中包括大型泥石流冲沟4条、中型3条，小型滑坡7个，威胁人口865人，威胁资产约3720万元。这些隐患点主要分布在山区和河谷地带，具体分布情况如下：在山区，泥石流隐患点多集中在坡度较陡、沟谷发育的区域。如大山村、大田坝村等地，由于地形起伏大，岩土体稳定性差，在强降雨条件下，极易发生泥石流灾害。其中，位于大山村的一条大型泥石流冲沟，沟谷长度超过2公里，沟床比降大，两侧山坡坡度多在35°以上，沟内松散固体物质丰富，一旦遭遇强降雨，泥石流爆发的可能性极大。该冲沟下游有多个村民小组，居住人口约200人，房屋100余间，还有大量的农田和果园，一旦发生泥石流灾害，将对当地居民的生命财产安全造成严重威胁。河谷地带的泥石流隐患点主要分布在勐波罗河和怒江沿岸。这些区域地势较低，河流流速快，在洪水期，河水容易漫溢，携带大量泥沙和石块，形成泥石流。例如，在旧城村靠近勐波罗河的区域，存在一处中型泥石流隐患点。该隐患点周边有部分居民房屋和农田，由于河流的冲刷作用，河岸土体松散，稳定性差。在雨季，当河水上涨时，容易引发河岸滑坡和泥石流灾害。据估算，该隐患点威胁人口约150人，资产约500万元。此外，在一些人类工程活动频繁的区域，如道路建设、采矿活动区等，也存在一定数量的泥石流隐患点。这些区域由于工程建设破坏了原有的地形地貌和岩土体稳定性，增加了泥石流发生的风险。例如，在旧城乡的一处采矿活动区，由于长期的采矿作业，山体被开挖，形成了大量的松散堆积物，这些堆积物在降雨作用下，容易引发泥石流灾害。该隐患点威胁到附近采矿工人和周边居民的安全，威胁资产约300万元。这些泥石流隐患点的存在，犹如一颗颗“定时炸弹”，时刻威胁着旧城乡居民的生命财产安全和社会经济的稳定发展。因此，加强对这些隐患点的监测和治理，是旧城乡泥石流灾害防治工作的当务之急。通过定期的监测，及时掌握隐患点的动态变化情况，提前发布预警信息，能够为居民的撤离和防范提供宝贵时间；采取有效的工程治理措施，如修建拦砂坝、排导槽等，能够降低泥石流的危害程度，保障人民群众的生命财产安全。三、泥石流风险评估指标体系构建3.1评估指标选取原则在构建云南省施甸县旧城乡集镇泥石流风险评估指标体系时，需遵循一系列科学合理的原则，以确保评估结果的准确性和可靠性。这些原则对于全面、客观地反映泥石流风险状况，以及为后续的风险评估和防灾减灾决策提供有力支持具有至关重要的意义。科学性是首要原则，要求所选取的评估指标必须基于泥石流形成的科学原理和相关理论知识。地质条件中的地层岩性、断裂构造等指标，它们与泥石流的物源和动力条件密切相关。地层岩性决定了岩土体的抗风化能力和稳定性，易风化破碎的岩石如砂岩、页岩等，在长期的风化作用下，会形成大量松散碎屑物质，为泥石流提供丰富的物源。断裂构造的存在破坏了岩石的完整性，降低了岩土体的力学强度，增加了山体的不稳定性，在强降雨等触发因素作用下，容易引发泥石流灾害。地形地貌指标如坡度、地形起伏度等，直接影响地表径流的速度和汇流情况，进而影响泥石流的形成和运动。只有基于科学原理选取这些指标，才能准确反映泥石流风险与各影响因素之间的内在联系，使评估结果具有科学依据和可信度。全面性原则强调评估指标应涵盖影响泥石流风险的各个方面，避免遗漏重要因素。除了地质和地形地貌因素外，气象因素如降雨量、降雨强度等也是引发泥石流的关键因素。强降雨会使地表径流迅速增加，土壤含水量饱和，导致岩土体抗剪强度降低，从而触发泥石流。人类活动因素同样不可忽视，山区的不合理工程建设、过度开垦、滥砍滥伐等活动，会破坏地表植被和地形地貌，增加松散物质来源，改变地表径流路径，进而加大泥石流发生的风险。因此，只有综合考虑地质、地形地貌、气象和人类活动等多方面因素，选取相应的评估指标，才能全面反映泥石流风险的形成机制和影响因素，为风险评估提供全面的数据支持。可操作性原则要求所选取的评估指标应易于获取和量化。在实际研究中，数据的获取往往受到各种条件的限制，如果选取的指标难以获取或无法准确量化，将导致评估工作无法顺利进行。对于地形地貌指标，可以利用高分辨率的遥感影像和数字高程模型（DEM），通过地理信息系统（GIS）软件进行快速提取和分析，获取坡度、坡向、沟谷密度等参数。气象数据可以从当地气象部门获取多年的降雨、气温等观测数据，经过整理和统计分析，得到降雨量、降雨强度等指标。地质数据可以通过地质勘察报告、野外调查和实验室测试等方式获取。人类活动指标可以通过实地调查、问卷调查等方式进行统计和量化。只有保证评估指标的可操作性，才能确保风险评估工作的顺利开展，提高评估效率和准确性。独立性原则要求各评估指标之间应相互独立，避免指标之间存在重复或高度相关的情况。如果指标之间存在重复信息，会导致评估结果的偏差和信息冗余，影响评估的准确性和科学性。在选取地质指标时，地层岩性和断裂构造是两个相互独立的因素，分别从岩土体的物质组成和地质结构方面影响泥石流风险。地形地貌指标中，坡度和地形起伏度虽然都与地形有关，但它们从不同角度反映地形特征，坡度主要影响地表径流的流速，地形起伏度则反映地形的复杂程度，两者相互独立。通过确保指标的独立性，可以提高评估指标体系的质量，使各指标能够独立地反映泥石流风险的不同方面，从而提高评估结果的准确性和可靠性。灵敏性原则要求评估指标对泥石流风险的变化具有较高的敏感性，能够及时准确地反映风险的变化情况。一些指标在泥石流风险变化时，其数值或状态会发生明显改变，这些指标就具有较高的灵敏性。降雨量和降雨强度是对泥石流风险变化非常敏感的指标，当降雨量和降雨强度超过一定阈值时，泥石流发生的可能性会急剧增加。坡度的变化也会对泥石流风险产生显著影响，坡度增大，岩土体的稳定性降低，泥石流发生的风险相应增加。通过选取灵敏性高的指标，可以及时捕捉到泥石流风险的变化趋势，为提前采取防灾减灾措施提供预警信息，降低灾害损失。3.2具体评估指标分析3.2.1地形地貌指标地形地貌在泥石流的形成和发展过程中扮演着关键角色，其相关指标对泥石流风险评估具有重要意义。坡度是影响泥石流发生的重要地形因素之一。当坡度较大时，岩土体在重力作用下的稳定性降低，更容易发生滑动和崩塌，为泥石流提供丰富的固体物质来源。同时，较大的坡度使得地表径流速度加快，水流的侵蚀和搬运能力增强，能够更有效地将松散物质带入沟谷，激发泥石流的形成。研究表明，坡度在25°-45°之间的区域，泥石流发生的概率相对较高。在旧城乡集镇，山区部分区域坡度超过35°，这些区域在强降雨等条件下，泥石流发生的风险显著增加。坡向主要通过影响光照、降水和植被生长等因素，间接对泥石流产生作用。阳坡接受的太阳辐射较多，气温较高，蒸发量大，植被生长相对较差，岩土体风化作用强烈，松散物质较多，在一定程度上增加了泥石流发生的可能性。阴坡则相对湿润，植被生长较好，对岩土体有一定的保护作用，泥石流发生的风险相对较低。此外，当坡向与盛行风向一致时，降水更容易在坡面汇聚，增加了泥石流发生的风险。在旧城乡，一些阳坡区域由于植被覆盖率低，岩土体稳定性差，在雨季时泥石流发生的风险较高。高差反映了地形的起伏程度，高差越大，地表径流的势能越大，水流速度越快，对岩土体的侵蚀和搬运能力越强，越有利于泥石流的形成和发展。在高差较大的山区，山体岩石在长期的风化、侵蚀作用下，容易破碎形成大量的松散物质，这些物质在强降雨等触发因素作用下，很容易被水流携带形成泥石流。旧城乡境内最高海拔达2595.2米，最低海拔仅583.88米，相对高差较大，山区地形起伏明显，这为泥石流的发生提供了有利的地形条件。沟谷密度和沟谷形态也与泥石流密切相关。沟谷密度大，意味着地表径流有更多的汇聚通道，容易形成集中的水流，增加了泥石流发生的可能性。“V”字形沟谷由于其狭窄陡深的特点，水流在沟谷中流速快，对沟谷两侧和底部的侵蚀作用强，容易导致沟谷两侧的岩土体崩塌、滑落，为泥石流提供物质来源。而“U”字形沟谷相对较为平缓，水流速度较慢，泥石流发生的风险相对较低。在旧城乡集镇，山区沟谷密度较大，且部分沟谷呈“V”字形，这些沟谷在雨季时是泥石流的高发区域。地形起伏度综合反映了地形的复杂程度，地形起伏度越大，地形越复杂，岩土体的稳定性越差，泥石流发生的风险越高。复杂的地形使得降水在地表的分布不均匀，容易形成局部的强降雨区域，同时也增加了地表径流的汇聚和流动路径的复杂性，为泥石流的形成创造了条件。通过对旧城乡集镇地形起伏度的分析，可以确定地形复杂、泥石流风险较高的区域，为风险评估和防治提供依据。3.2.2地质条件指标地质条件是泥石流形成的内在基础，对泥石流的发生起着关键的控制作用，其中岩土体类型和地质构造是两个重要的指标。岩土体类型决定了泥石流固体物质的来源和性质。不同类型的岩土体具有不同的物理力学性质，其抗风化、抗侵蚀能力以及稳定性存在显著差异。在旧城乡集镇，砂岩、页岩和泥岩分布较为广泛。砂岩硬度相对较低，抗风化能力较弱，在长期的风化作用下，容易破碎形成大量的砂质碎屑，这些碎屑颗粒较小，在水流作用下容易被搬运，为泥石流提供了丰富的细颗粒物质。页岩具有明显的页理构造，遇水后容易软化、崩解，强度大幅降低，其结构的不稳定性使得页岩地区在强降雨等条件下容易发生滑坡和泥石流灾害。泥岩的黏性较大，透水性差，遇水后会形成泥浆状物质，增加了泥石流的黏性和破坏力。此外，残积土、坡积土和洪积土等松散堆积物在旧城乡也有一定分布，这些堆积物结构松散，颗粒之间的黏聚力较小，在重力和水流作用下，极易发生滑动和崩塌，为泥石流的形成提供了直接的物质来源。地质构造对泥石流的发生和发展有着重要影响。断裂构造是岩石中的破裂面，它破坏了岩石的完整性，降低了岩土体的力学强度。在断裂带附近，岩石破碎，裂隙发育，地下水活动频繁，使得岩土体的稳定性受到严重影响。当断裂带处于活动状态时，会引发地震等地质活动，进一步加剧山体的不稳定，增加泥石流发生的风险。旧城乡位于勐波罗河大断裂带上与怒江断裂带以北，区内断裂构造发育，这些断裂带的存在使得区域内地层破碎，为泥石流的形成提供了有利的地质条件。褶皱构造改变了地层的原始产状，形成了复杂的地形地貌。在褶皱的轴部和翼部，岩石受力集中，容易发生破裂和变形，导致岩土体的稳定性降低。此外，褶皱构造还会影响地下水的流动和分布，使得局部地区地下水水位升高，岩土体饱水，抗剪强度降低，从而增加了泥石流发生的可能性。地层的产状和厚度也与泥石流风险密切相关。倾斜的地层在重力作用下，岩土体有向下滑动的趋势，当受到外部因素如强降雨、地震等作用时，更容易发生滑动和崩塌。地层厚度较大时，意味着有更多的岩土体参与泥石流的形成，可能会导致泥石流的规模更大，破坏力更强。在旧城乡集镇，部分山区地层倾斜角度较大，且厚度变化明显，这些区域在地质条件和外部触发因素的共同作用下，泥石流发生的风险较高。3.2.3气象水文指标气象水文条件是泥石流发生的重要触发因素，降雨量、降雨强度、河流分布等指标对泥石流的形成和发展具有显著影响。降雨量和降雨强度是引发泥石流的关键气象因素。大量的降雨是泥石流形成的必要条件之一，充足的水源能够使岩土体饱水，降低其抗剪强度，同时形成强大的地表径流，为泥石流的发生提供动力。当降雨量达到一定程度时，土壤孔隙被水充满，水分开始在地表积聚，形成坡面径流。随着降雨量的增加，坡面径流的流量和流速不断增大，对地表的侵蚀作用增强，能够将地表的松散物质带入沟谷。降雨强度对泥石流的影响更为直接，短时间内的高强度降雨会使地表径流迅速增加，来不及下渗和蒸发，从而形成强大的洪流。这种高强度的水流具有巨大的能量，能够迅速冲刷和搬运大量的泥沙、石块等固体物质，激发泥石流的发生。研究表明，当短时间内降雨量超过一定阈值，如日降雨量超过50毫米，或小时降雨量超过20毫米时，泥石流发生的概率会显著增加。在旧城乡，每年雨季（5-10月）降雨集中，且多暴雨天气，日最大降雨量可达150毫米以上，这些强降雨事件是引发泥石流的主要原因。降雨历时也不容忽视，长时间的降雨会使岩土体持续饱水，逐渐软化，稳定性不断降低。即使降雨强度不大，但如果降雨历时较长，也可能导致地下水位上升，岩土体处于饱和状态，在重力作用下容易发生滑动和崩塌，进而引发泥石流。例如，连续降雨3天以上，且累计降雨量达到100毫米以上时，旧城乡山区部分区域就可能出现岩土体失稳的情况，增加泥石流发生的风险。河流分布对泥石流的形成和发展有着重要影响。河流是地表径流的主要通道，其流量和流速的变化直接影响着泥石流的发生。在河流的上游和支流区域，地形坡度较大，河流流速快，对河岸和河床的侵蚀作用强烈，容易导致河岸崩塌和河床物质的搬运，为泥石流提供物质来源。当河流遇到强降雨时，河水水位迅速上涨，流量和流速急剧增大，可能会漫溢到河岸两侧，携带大量泥沙和石块，形成泥石流。此外，河流的弯曲程度和河谷形态也会影响泥石流的发生。在河流弯曲处，水流速度不均匀，外侧河岸受到的冲刷力较大，容易发生崩塌和滑坡，为泥石流提供物质。河谷狭窄且两岸坡度较陡的区域，河流在洪水期的能量集中，对河谷两侧的侵蚀作用更强，更容易引发泥石流。旧城乡境内有勐波罗河和怒江等主要河流，这些河流的存在使得河谷地带成为泥石流的高发区域。河网密度反映了区域内河流的密集程度，河网密度大的区域，地表径流更容易汇聚，增加了泥石流发生的可能性。在河网密集的地区，降雨后地表径流能够迅速通过众多的河流通道汇聚，形成强大的水流，对地表的侵蚀和搬运能力增强。同时，河网的存在也使得地下水与地表水的交换更加频繁，可能导致地下水位升高，岩土体饱水，稳定性降低。在旧城乡集镇，山区部分区域河网密度较大，这些区域在雨季时泥石流发生的风险相对较高。3.2.4人类活动指标人类活动对泥石流的诱发作用日益显著，工程建设、农业开发等活动改变了地表的自然状态，增加了泥石流发生的风险。工程建设活动如道路修建、房屋建筑、采矿等对地形地貌和岩土体稳定性产生了深刻影响。在山区进行道路修建时，往往需要开挖山体、填方筑路，这会破坏原有的地形地貌和植被，形成大量的松散堆积物。这些堆积物如果没有得到妥善处理，在降雨等条件下，容易发生滑坡和泥石流灾害。道路的修建还可能改变地表径流的路径和汇流条件，使得水流集中，增加了对道路两侧山体的冲刷力。例如，在旧城乡的一些山区道路建设过程中，由于开挖山体产生的弃土随意堆放在山坡上，没有采取有效的防护措施，在雨季时这些弃土被雨水冲刷，引发了小规模的泥石流灾害，对道路和周边居民造成了一定的影响。房屋建筑活动也会对泥石流产生影响。在山坡上不合理地选址建房，会破坏山体的稳定性，增加滑坡和泥石流的风险。为了平整地基，可能会对山坡进行开挖和填方，改变了山体的自然坡度和岩土体结构。此外，房屋建设过程中产生的建筑垃圾如果随意堆放，也会为泥石流提供物质来源。在旧城乡的一些村庄，部分居民在山坡上自行建房，没有进行科学的规划和地质勘察，导致房屋周边山体在降雨后出现滑坡迹象，威胁到居民的生命财产安全。采矿活动对地质环境的破坏更为严重。采矿过程中，大量的岩土体被开挖和扰动，形成了众多的采空区和废石堆。采空区的存在使得山体内部结构变得不稳定，容易引发地面塌陷和山体滑坡。废石堆通常堆积在矿区周边，这些废石结构松散，缺乏植被覆盖，在降雨作用下，极易被冲刷形成泥石流。例如，旧城乡的一些小型采矿场，由于缺乏规范的管理和有效的环保措施，废石随意堆放，在暴雨后曾多次引发泥石流灾害，不仅对矿区造成了破坏，还对下游的村庄和农田构成了威胁。农业开发活动如开垦耕地、砍伐森林等也会加剧泥石流的发生风险。开垦耕地会破坏原有的植被，使得土壤失去植被的保护，抗侵蚀能力降低。在山区进行陡坡开垦，容易导致水土流失，大量的泥沙被雨水冲刷进入沟谷，为泥石流提供物质来源。砍伐森林会削弱森林的水土保持功能，森林植被具有截留雨水、涵养水源、固持土壤等作用，森林被砍伐后，这些功能丧失，地表径流增加，土壤侵蚀加剧，泥石流发生的可能性增大。在旧城乡，部分山区存在过度开垦和滥砍滥伐的现象，导致植被覆盖率下降，水土流失严重，在雨季时泥石流发生的频率明显增加。此外，人类活动还可能改变区域的水文条件。如修建水库、灌溉渠道等水利设施，可能会改变地下水的水位和流向，影响岩土体的稳定性。不合理的灌溉方式可能导致土壤水分过多，岩土体饱水，增加泥石流发生的风险。在旧城乡，一些农田灌溉渠道由于设计不合理或维护不善，在雨季时出现漏水现象，导致周边土壤含水量过高，引发了小规模的滑坡和泥石流灾害。3.2.5承灾体指标承灾体是泥石流灾害的承受对象，人口密度、建筑物分布等承灾体因素与泥石流灾害损失密切相关，直接影响着泥石流风险的评估结果。人口密度是衡量承灾体易损性的重要指标之一。在人口密度大的区域，一旦发生泥石流灾害，受到威胁的人口数量众多，人员伤亡和社会影响的风险也相应增加。旧城乡集镇的河谷地带人口密度较大，这些区域如果发生泥石流灾害，可能会造成大量人员伤亡和财产损失。相比之下，山区人口相对稀疏，虽然泥石流发生的风险较高，但灾害造成的人员伤亡和社会影响相对较小。此外，人口的分布特征也会影响灾害损失。如果人口集中分布在泥石流沟道下游或地势低洼等危险区域，那么在灾害发生时，人员疏散和救援难度较大，灾害损失将更为严重。建筑物分布是影响泥石流灾害损失的关键因素。建筑物的数量、结构类型和抗震性能等都会对灾害损失产生重要影响。在泥石流灾害发生时，建筑物可能会受到泥石流的冲击、掩埋和倒塌等破坏。数量众多的建筑物意味着更大的财产损失风险。不同结构类型的建筑物对泥石流的抵抗能力不同，砖混结构和框架结构的建筑物相对较为坚固，在一定程度上能够抵御泥石流的冲击；而土木结构和简易结构的建筑物则较为脆弱，容易在泥石流灾害中受损或倒塌。建筑物的抗震性能也与灾害损失密切相关，抗震性能差的建筑物在泥石流引发的地面震动和冲击作用下，更容易发生倒塌，造成人员伤亡和财产损失。在旧城乡集镇，部分老旧房屋多为土木结构，且建筑年代久远，抗震性能较差，这些房屋在泥石流灾害中面临着较大的风险。基础设施如道路、桥梁、电力、通信等的分布和状况也对泥石流灾害损失有着重要影响。道路和桥梁是交通的重要枢纽，在泥石流灾害中，它们可能会被冲毁或掩埋，导致交通中断，影响救援物资的运输和人员的疏散。电力和通信设施的损坏会导致停电和通信中断，影响应急救援工作的开展和居民的正常生活。在旧城乡，一些山区道路和桥梁在历史上的泥石流灾害中曾多次受损，导致交通中断，给当地居民的生产生活带来了极大不便。此外，承灾体的经济价值也是评估泥石流风险的重要因素。经济价值高的区域，如商业区、工业园区等，一旦遭受泥石流灾害，经济损失将更为巨大。旧城乡虽然工业和商业相对不发达，但一些重要的农业生产区域和居民聚居区具有较高的经济价值，这些区域如果发生泥石流灾害，将对当地的经济发展造成严重影响。四、泥石流风险评估方法选择与应用4.1常用风险评估方法概述泥石流风险评估方法种类繁多，每种方法都有其独特的原理、适用范围和优缺点，在实际应用中需根据具体情况进行合理选择。经验法是一种基于专家经验和历史数据的评估方法。其原理主要是依靠专家对泥石流形成条件、运动特征以及历史灾害案例的了解和认识，通过实地考察、资料分析等方式，对泥石流的风险进行定性或半定量的评估。在评估某一区域的泥石流风险时，专家会综合考虑该区域的地形地貌、地质条件、气象因素以及历史上泥石流的发生情况等，凭借自身经验判断该区域泥石流发生的可能性和危害程度。这种方法的优点是简单易行，不需要复杂的数学模型和大量的数据支持，能够快速得出评估结果。在一些数据匮乏的地区，经验法可以在短时间内对泥石流风险做出初步判断，为后续的防灾减灾工作提供参考。然而，经验法也存在明显的局限性。其评估结果主观性较强，不同专家的经验和判断标准可能存在差异，导致评估结果缺乏一致性和可比性。该方法难以对泥石流风险进行精确的量化评估，对于一些复杂的泥石流灾害情况，仅凭经验可能无法准确评估其风险。半定量法是在经验法的基础上，引入一些量化指标和模型，对泥石流风险进行相对定量的评估。层次分析法（AHP），它通过构建层次结构模型，将复杂的泥石流风险评估问题分解为多个层次和因素，然后通过两两比较的方式确定各因素的相对重要性权重，最后综合计算得出泥石流风险的相对等级。模糊综合评价法，它利用模糊数学的理论，将模糊的风险因素进行量化处理，通过建立模糊关系矩阵和综合评价模型，对泥石流风险进行综合评价，得出风险等级。半定量法的优点是在一定程度上克服了经验法的主观性，通过量化指标和模型的运用，使评估结果更加客观、准确。它能够综合考虑多个因素对泥石流风险的影响，比经验法更全面、系统。半定量法仍然存在一些不足。其评估结果的准确性依赖于指标的选取和权重的确定，若指标选取不合理或权重确定不准确，会影响评估结果的可靠性。该方法对于一些难以量化的因素，如人类活动对泥石流风险的潜在影响等，处理起来存在一定困难。定量法是运用数学模型和大量的数据，对泥石流的发生概率、规模、危害程度等进行精确量化评估的方法。基于概率统计的方法，通过对历史泥石流事件的发生频率、规模等数据进行统计分析，建立概率模型，预测未来泥石流发生的概率和可能造成的损失。数值模拟方法，如FLO-2D、RAMMS等模型，利用计算机模拟技术，根据研究区域的地形地貌、地质条件、气象数据等，模拟泥石流的运动过程，包括流速、流量、泥深、堆积范围等参数，从而精确评估泥石流的风险。定量法的优点是能够对泥石流风险进行精确量化，评估结果具有较高的科学性和可靠性。它可以为防灾减灾决策提供具体的数据支持，有助于制定更加科学合理的防治措施。然而，定量法对数据的要求较高，需要大量准确的历史数据和实时监测数据，数据获取难度较大且成本较高。模型的建立和求解过程较为复杂，需要专业的技术人员和先进的计算设备，在实际应用中受到一定限制。4.2评估方法选择依据旧城乡集镇的泥石流风险评估选用层次分析法（AHP）和模糊综合评价法相结合的方式，主要基于以下考虑。旧城乡地形地貌复杂，山区面积广大，坡度陡峻，沟谷纵横交错，且地质构造处于勐波罗河大断裂带上与怒江断裂带以北，地层破碎，岩性多样，这些因素相互交织，使得泥石流风险影响因素众多且关系复杂。同时，该区域气象条件多变，干湿季分明，雨季降水集中且多暴雨，人类活动如工程建设、农业开发等对地表扰动较大，进一步增加了泥石流风险评估的复杂性。在数据获取方面，虽然可以通过地质勘察、地形测量、气象观测等手段获取部分数据，但仍存在一些难以精确量化的数据，如人类活动对泥石流风险的潜在影响等。层次分析法（AHP）具有高度的逻辑性、系统性、简洁性与实用性，适用于目标结构复杂且缺乏必要数据的情况。它能够将复杂的泥石流风险评估问题分解为多个层次和因素，通过专家经验和判断，对各因素进行两两比较，确定其相对重要性权重。在旧城乡的泥石流风险评估中，利用AHP可以将地质条件、地形地貌、气象因素、人类活动等多个影响因素进行层次划分，构建层次结构模型，从而有效地处理各因素之间的复杂关系。对于地层岩性、断裂构造、坡度、降雨量等因素，可以通过专家打分的方式，确定它们在泥石流风险评估中的相对重要性，为后续的综合评价提供权重依据。模糊综合评价法能够很好地处理模糊性和不确定性问题。在旧城乡的泥石流风险评估中，许多因素具有模糊性，如泥石流发生的可能性、危害程度等，难以用精确的数值来描述。模糊综合评价法利用模糊数学的理论，将这些模糊因素进行量化处理，通过建立模糊关系矩阵和综合评价模型，对泥石流风险进行综合评价。对于泥石流发生可能性的评估，可以将其划分为“高”“中”“低”等模糊等级，然后通过模糊关系矩阵，综合考虑各种影响因素，得出最终的风险评价结果。将层次分析法（AHP）和模糊综合评价法相结合，既能充分发挥AHP在处理复杂因素权重确定方面的优势，又能利用模糊综合评价法处理模糊性和不确定性问题的能力，从而更全面、准确地评估旧城乡集镇的泥石流风险。这种方法组合在类似复杂地质和地理条件的区域泥石流风险评估中也具有较高的适用性和参考价值。4.3基于层次分析法和模糊综合评价法的旧城乡集镇泥石流风险评估过程4.3.1数据收集与整理数据收集与整理是泥石流风险评估的基础环节，其准确性和完整性直接影响后续评估结果的可靠性。在本次对旧城乡集镇泥石流风险评估中，通过多渠道广泛收集各类数据。地质数据主要来源于施甸县地质矿产局的地质勘察报告、旧城乡的地质图以及相关的地质研究文献。这些资料详细记录了区域内地层岩性的分布情况，如砂岩、页岩、泥岩、灰岩等各类岩石的分布范围和厚度；断裂构造的位置、走向和规模，明确了勐波罗河断裂、怒江断裂及其次级断裂在旧城乡的具体分布；岩土体的物理力学性质参数，包括岩土体的抗剪强度、抗压强度、容重、渗透性等，为分析地质条件对泥石流的影响提供了关键数据。地形地貌数据借助高分辨率的遥感影像和数字高程模型（DEM）获取。利用遥感影像能够直观地识别地形地貌特征，如山脉、沟谷、河流的分布情况。通过DEM数据，在地理信息系统（GIS）软件中进行处理和分析，精确提取坡度、坡向、高差、沟谷密度、地形起伏度等地形地貌指标。在提取坡度时，通过GIS的空间分析功能，对DEM数据进行计算，得到每个栅格单元的坡度值，从而生成坡度图；在提取沟谷密度时，利用水文分析工具，根据DEM数据生成水流方向、汇流累积量等数据，进而识别出沟谷网络，计算出沟谷密度。气象水文数据从施甸县气象局和水文站收集。获取多年的降雨量、降雨强度、降雨历时、气温、风速、河流流量、河网密度等数据，并对其进行整理和统计分析。将多年的降雨量数据按年份、月份进行统计，分析降雨量的年际变化和年内分配规律；对降雨强度数据进行统计，确定不同量级降雨强度的发生频率；通过对河流流量数据的分析，了解河流的水文特征及其与泥石流发生的关系。人类活动数据通过实地调查、问卷调查和相关部门的统计资料收集。实地调查主要观察区域内工程建设、农业开发、采矿等活动的情况，记录工程建设的规模、位置，农业开垦的面积、方式，采矿活动的范围、开采方式等信息。问卷调查面向当地居民和相关企业，了解人类活动对地表植被、地形地貌的破坏程度，以及对泥石流灾害的认知和防范措施。相关部门的统计资料，如土地利用规划图、矿产资源开发规划等，提供了人类活动的宏观信息，为评估人类活动对泥石流风险的影响提供了全面的数据支持。社会经济数据从旧城乡政府的统计报表、经济发展规划等资料中获取。收集人口密度、建筑物分布、基础设施分布、经济发展水平等数据，了解区域内承灾体的情况。通过人口普查数据，获取不同区域的人口密度信息；根据建筑物普查资料，掌握建筑物的数量、结构类型和分布情况；从基础设施建设规划中，了解道路、桥梁、电力、通信等基础设施的分布和建设情况。对收集到的各类数据进行整理和预处理。对数据进行清洗，去除重复、错误和缺失的数据；对不同格式的数据进行统一格式转换，使其便于后续的分析和处理；将各类数据按照评估指标体系的要求进行分类存储，建立数据库，为后续的风险评估提供数据支持。4.3.2指标量化与标准化指标量化与标准化是将不同性质、不同量纲的评估指标转化为具有可比性的数值，以便于后续的分析和计算。对于地形地貌指标，坡度可直接从DEM数据中提取得到，取值范围一般为0°-90°。为了便于分析和计算，将其进行标准化处理，采用极差标准化方法，公式为：x_{i}^{*}=\frac{x_{i}-x_{min}}{x_{max}-x_{min}}其中，x_{i}^{*}为标准化后的指标值，x_{i}为原始指标值，x_{min}和x_{max}分别为该指标的最小值和最大值。在旧城乡集镇，坡度最小值为0°，最大值为60°，若某区域原始坡度值为30°，则标准化后的坡度值为：x_{i}^{*}=\frac{30-0}{60-0}=0.5坡向取值范围为0°-360°，同样采用极差标准化方法进行处理。将坡向划分为不同的区间，如0°-90°为阳坡，90°-180°为半阳坡，180°-270°为阴坡，270°-360°为半阴坡。对不同区间赋予不同的分值，阳坡赋值为1，半阳坡赋值为0.75，阴坡赋值为0.5，半阴坡赋值为0.25，然后进行标准化处理。高差通过计算区域内最高海拔与最低海拔的差值得到。假设旧城乡集镇最高海拔为2595.2米，最低海拔为583.88米，某区域高差为1000米，则标准化后的高差指标值为：x_{i}^{*}=\frac{1000-0}{2595.2-583.88}\approx0.497沟谷密度通过计算单位面积内沟谷的长度得到。若某区域面积为10平方公里，沟谷总长度为5公里，则沟谷密度为0.5公里/平方公里。将沟谷密度进行标准化处理，使其取值范围在0-1之间。地形起伏度通过计算一定窗口范围内的高程标准差得到。在GIS软件中，设置合适的窗口大小，如500米×500米的窗口，计算每个窗口内的地形起伏度。然后采用极差标准化方法对地形起伏度进行处理。对于地质条件指标，岩土体类型可根据其稳定性和抗风化能力进行赋值。砂岩、页岩等抗风化能力较弱、稳定性较差的岩土体赋值为0.8，泥岩、灰岩等赋值为0.6，花岗岩等稳定性较好的岩土体赋值为0.4。然后进行标准化处理。地质构造根据断裂构造的规模、活动性等因素进行赋值。规模较大、活动性较强的断裂构造附近区域赋值为1，规模较小、活动性较弱的断裂构造附近区域赋值为0.6，远离断裂构造的区域赋值为0.2。再进行标准化处理。对于气象水文指标，降雨量和降雨强度采用极差标准化方法进行处理。假设旧城乡集镇多年平均降雨量最小值为800毫米，最大值为1200毫米，某区域年降雨量为1000毫米，则标准化后的降雨量指标值为：x_{i}^{*}=\frac{1000-800}{1200-800}=0.5降雨历时可根据不同的时长范围进行赋值。降雨历时小于1小时赋值为0.2，1-3小时赋值为0.5，3-6小时赋值为0.8，大于6小时赋值为1。然后进行标准化处理。河流流量和河网密度同样采用极差标准化方法进行处理。对于人类活动指标，工程建设、农业开发、采矿等活动的强度可根据活动的规模、范围等因素进行赋值。大规模的工程建设、过度开垦和采矿活动区域赋值为1，中等规模的活动区域赋值为0.6，小规模的活动区域赋值为0.2。然后进行标准化处理。对于承灾体指标，人口密度根据不同的人口数量范围进行赋值。人口密度大于500人/平方公里赋值为1，200-500人/平方公里赋值为0.6，小于200人/平方公里赋值为0.2。然后进行标准化处理。建筑物结构类型根据其抗震性能和抗泥石流冲击能力进行赋值。砖混结构和框架结构赋值为0.4，土木结构和简易结构赋值为0.8。然后进行标准化处理。通过以上指标量化与标准化处理，使各评估指标具有了可比性，为后续的风险评估提供了数据基础。4.3.3权重确定权重确定是泥石流风险评估中的关键步骤，它反映了各评估指标对泥石流风险的相对重要程度。本研究运用层次分析法（AHP）来确定各指标的权重。首先，构建层次结构模型。将泥石流风险评估问题分解为目标层、准则层和指标层。目标层为旧城乡集镇泥石流风险评估；准则层包括地形地貌、地质条件、气象水文、人类活动和承灾体五个方面；指标层则包含坡度、坡向、高差、沟谷密度、地形起伏度、岩土体类型、地质构造、降雨量、降雨强度、降雨历时、河流流量、河网密度、工程建设强度、农业开发强度、采矿强度、人口密度、建筑物结构类型等具体指标。然后，通过专家问卷调查的方式，邀请地质、地理、气象、灾害防治等领域的专家，对准则层和指标层各因素进行两两比较，判断其相对重要程度。采用1-9标度法进行量化，1表示两个因素同等重要，3表示一个因素比另一个因素稍微重要，5表示一个因素比另一个因素明显重要，7表示一个因素比另一个因素强烈重要，9表示一个因素比另一个因素极端重要，2、4、6、8为上述相邻判断的中间值。例如，在判断地形地貌和地质条件对泥石流风险的相对重要程度时，若专家认为地形地貌比地质条件稍微重要，则赋值为3。根据专家打分结果，构建判断矩阵。以准则层为例，假设地形地貌、地质条件、气象水文、人类活动和承灾体分别用B1、B2、B3、B4、B5表示，则判断矩阵A为：A=\begin{pmatrix}1&3&2&4&3\\\frac{1}{3}&1&\frac{1}{2}&2&1\\\frac{1}{2}&2&1&3&2\\\frac{1}{4}&\frac{1}{2}&\frac{1}{3}&1&\frac{1}{2}\\\frac{1}{3}&1&\frac{1}{2}&2&1\end{pmatrix}计算判断矩阵的特征向量和最大特征根。利用方根法或和积法等方法进行计算。以方根法为例，计算步骤如下：计算判断矩阵A中每行元素的乘积M_{i}：M_{1}=1\times3\times2\times4\times3=72M_{2}=\frac{1}{3}\times1\times\frac{1}{2}\times2\times1=\frac{1}{3}M_{3}=\frac{1}{2}\times2\times1\times3\times2=6M_{4}=\frac{1}{4}\times\frac{1}{2}\times\frac{1}{3}\times1\times\frac{1}{2}=\frac{1}{48}M_{5}=\frac{1}{3}\times1\times\frac{1}{2}\times2\times1=\frac{1}{3}计算M_{i}的n次方根\overline{W}_{i}（n为判断矩阵的阶数，此处n=5）：\overline{W}_{1}=\sqrt[5]{72}\approx2.27\overline{W}_{2}=\sqrt[5]{\frac{1}{3}}\approx0.80\overline{W}_{3}=\sqrt[5]{6}\approx1.43\overline{W}_{4}=\sqrt[5]{\frac{1}{48}}\approx0.53\overline{W}_{5}=\sqrt[5]{\frac{1}{3}}\approx0.80对\overline{W}_{i}进行归一化处理，得到特征向量W_{i}：W_{1}=\frac{\overline{W}_{1}}{\sum_{i=1}^{n}\overline{W}_{i}}=\frac{2.27}{2.27+0.80+1.43+0.53+0.80}\approx0.39W_{2}=\frac{\overline{W}_{2}}{\sum_{i=1}^{n}\overline{W}_{i}}=\frac{0.80}{2.27+0.80+1.43+0.53+0.80}\approx0.14W_{3}=\frac{\overline{W}_{3}}{\sum_{i=1}^{n}\overline{W}_{i}}=\frac{1.43}{2.27+0.80+1.43+0.53+0.80}\approx0.25W_{4}=\frac{\overline{W}_{4}}{\sum_{i=1}^{n}\overline{W}_{i}}=\frac{0.53}{2.27+0.80+1.43+0.53+0.80}\approx0.09W_{5}=\frac{\overline{W}_{5}}{\sum_{i=1}^{n}\overline{W}_{i}}=\frac{0.80}{2.27+0.80+1.43+0.53+0.80}\approx0.14计算最大特征根\lambda_{max}：\lambda_{max}=\sum_{i=1}^{n}\frac{(AW)_{i}}{nW_{i}}其中，(AW)_{i}为判断矩阵A与特征向量W的乘积的第i个元素。经计算，\lambda_{max}\approx5.03。进行一致性检验。计算一致性指标CI：CI=\frac{\lambda_{max}-n}{n-1}=\frac{5.03-5}{5-1}=0.0075查找平均随机一致性指标RI，当n=5时，RI=1.12。计算一致性比例CR：CR=\frac{CI}{RI}=\frac{0.00