黑水县俄瓜沟泥石流：特征剖析与风险性精准评价

黑水县俄瓜沟泥石流：特征剖析与风险性精准评价一、引言1.1研究背景与意义泥石流作为山区常见且极具破坏力的地质灾害，严重威胁着人类生命财产安全与生态环境稳定。中国是世界上泥石流灾害最为严重的国家之一，众多山区深受其扰，泥石流灾害频发，每年都造成巨大的人员伤亡和经济损失。据统计，2005-2020年间，全国共发生泥石流1.3万起，因灾死亡失踪超4000人，尤其在西南山区，突发、群发性泥石流频率高、范围广，极易引发群死群伤事件，严重阻碍当地社会经济发展。黑水县位于四川省西北部的阿坝藏族羌族自治州中部，地处青藏高原东南缘横断山脉中段北端的岷江上游。境内山高坡陡，地形地貌复杂，地质构造活跃，降水集中且多暴雨，这些自然条件使得黑水县成为泥石流灾害的高发区域。俄瓜沟作为黑水县内的一条典型泥石流沟，其泥石流活动频繁，给当地带来了严重的灾害威胁。2012年7月9日，黑水县双溜索乡双溜索村俄瓜热十多沟流域发生了较大规模的泥石流灾害。此次灾害由突降暴雨引发，暴雨冲刷坡面，洪水揭底冲刷、掏蚀沟岸，携带大量松散固体物源从沟道内倾泻而出。泥石流堆积于沟口并进入黑水河，形成平均高度达2.66m的堰塞坝，一次冲出固体物质约10.28×10⁴m³。沟内地形陡峻，相对高差达2155m，主沟平均纵坡降为268.6‰，纵向上沟谷形态呈折线型，发育12道跌坎，具备爆发泥石流的多元松散物源和势能动力条件，并且管道增大效应明显。该泥石流形成堰塞湖长约1000m，淹没上游部分安置区房屋，冲毁对岸及下游部分房屋，造成一人失踪，沿线乡村公路约1000m被冲毁，省道S302线公路两侧民房不同程度被淹。一旦泥石流堰塞坝溃坝，还将直接威胁黑水河下游两岸居民和双溜索乡镇的安全。2024年6月20日下午17时50分左右，黑水县芦花镇片区发生短时强降水，最大降雨量达59.3mm，降雨造成县城新区出现内涝、多处区域发生泥石流，部分车辆受损，部分区域断电，尽管暂无人员伤亡，但也再次敲响了泥石流灾害防范的警钟。对黑水县俄瓜沟泥石流的基本特征及风险性进行深入研究，具有极其重要的意义。从防灾减灾角度来看，准确掌握俄瓜沟泥石流的形成条件、运动特征、堆积规律等基本特征，能够为泥石流灾害的监测、预警提供科学依据。通过建立科学合理的风险性评价体系，评估泥石流灾害发生的可能性和危害程度，可以提前制定有效的防灾减灾措施，如合理规划居民点和基础设施布局，避免在高风险区域建设；加强对泥石流沟的治理，采取工程措施和生物措施相结合的方式，减少泥石流的发生频率和规模；制定应急预案，提高应对泥石流灾害的能力，从而最大限度地减少人员伤亡和财产损失。从区域发展规划角度而言，研究俄瓜沟泥石流有助于为黑水县的土地利用规划、交通规划、水利设施建设等提供决策支持。在进行区域发展规划时，充分考虑泥石流灾害的风险，可以避免因规划不合理而导致的灾害损失，保障区域经济社会的可持续发展。例如，在交通规划中，合理选择公路、铁路等交通线路的走向，避开泥石流高风险区域，减少交通设施遭受泥石流破坏的可能性；在水利设施建设中，考虑泥石流对水库、大坝等的影响，采取相应的防护措施，确保水利设施的安全运行。1.2国内外研究现状泥石流作为一种复杂的地质灾害，其基本特征与风险性评价一直是国内外学者研究的重点领域。在国外，泥石流研究起步较早，发展较为成熟。早期研究主要集中在泥石流的形成机制方面，如日本学者足立胜治等在1977年开展了泥石流发生危险度的判定研究，从地貌条件、泥石流形态和降雨三方面判定泥石流发生率，为后续研究开辟了道路。随着研究的深入，逐渐拓展到泥石流的运动特征和堆积规律。在运动特征研究中，借助先进的监测技术和数值模拟方法，对泥石流的流速、流量、冲击力等参数进行精确测定和模拟分析。例如，利用多普勒雷达、高速摄像机等设备，实时监测泥石流的运动过程，获取其运动参数；通过数值模拟软件，如FLUENT、DEM等，对泥石流在不同地形条件下的运动路径和速度变化进行模拟，为泥石流灾害的防治提供理论支持。在堆积规律研究中，注重分析泥石流堆积物的粒度分布、结构特征以及堆积形态与地形、流量等因素的关系。通过对大量泥石流堆积现场的实地调查和样品分析，建立了堆积物特征与泥石流形成条件之间的关联模型，为泥石流灾害的评估和预测提供了重要依据。在泥石流风险性评价方面，国外学者建立了多种评价模型和方法。其中，基于统计分析的评价方法，通过对历史泥石流灾害数据的收集和整理，运用概率统计理论，分析泥石流发生的概率和危害程度。如美国学者Holling在1981年将洪水危险度的研究模式引入泥石流危险度研究中，通过对大量历史数据的统计分析，建立了泥石流危险度与地形、降雨等因素之间的统计关系模型。基于物理过程的评价方法，则从泥石流的形成、运动和堆积的物理机制出发，建立数学模型来模拟泥石流的发生和发展过程，进而评估其风险性。例如，欧洲一些国家的学者利用数值模拟技术，建立了考虑地形、地质、水文等多因素的泥石流物理过程模型，对泥石流的运动路径、堆积范围和危害程度进行精确预测。地理信息系统（GIS）和遥感（RS）技术在国外泥石流研究中也得到了广泛应用。通过RS技术获取大面积的地形、地貌、植被等信息，利用GIS强大的空间分析功能，对这些信息进行整合和分析，实现对泥石流灾害的快速识别、监测和评估。例如，利用卫星遥感影像解译泥石流沟的分布范围和形态特征，结合地形数据和降雨数据，在GIS平台上进行空间分析，绘制泥石流危险性分区图，为区域规划和灾害防治提供科学依据。国内的泥石流研究始于20世纪中叶，经过多年的发展，取得了丰硕的成果。在泥石流基本特征研究方面，国内学者结合我国山区的地质、地形和气候条件，对泥石流的形成条件进行了深入分析。研究发现，我国泥石流的形成主要受地形地貌、地层岩性、地质构造、水文气象等因素的控制，其中，地形陡峻、岩石破碎、降水集中是导致泥石流发生的关键因素。在运动特征研究中，通过野外实地观测和室内模拟实验，对泥石流的运动速度、流量、冲击力等参数进行了大量的研究。例如，在云南蒋家沟、四川映秀等泥石流频发地区，建立了长期的观测站，对泥石流的运动过程进行实时监测，获取了大量宝贵的数据；通过室内水槽实验，模拟泥石流在不同坡度、流量和物源条件下的运动特征，分析其运动规律。在堆积规律研究中，对泥石流堆积物的特征和分布规律进行了详细研究，发现堆积物的粒度组成、分选性和磨圆度等特征与泥石流的运动过程和形成条件密切相关，堆积物的分布范围和厚度也受到地形和泥石流流量的影响。在风险性评价方面，国内学者针对我国泥石流灾害的特点，提出了一系列适合我国国情的评价方法和指标体系。早期的评价方法主要是定性分析，通过对泥石流沟的地形地貌、物源条件、历史灾害记录等因素的综合分析，对泥石流的危险性进行初步评估。随着科学技术的发展，逐渐向定量评价转变。例如，采用层次分析法（AHP）、模糊综合评价法等数学方法，对影响泥石流危险性的多个因素进行量化分析，确定各因素的权重，从而建立泥石流危险性评价模型。同时，结合我国的实际情况，建立了一系列具有针对性的评价指标体系，如考虑地形地貌、地层岩性、地质构造、水文气象、人类工程活动等多因素的泥石流危险性评价指标体系，以及考虑承灾体类型、数量、价值等因素的泥石流易损性评价指标体系。此外，国内学者还将3S技术（GIS、RS、GPS）与泥石流研究相结合，实现了对泥石流灾害的动态监测和实时预警。通过建立泥石流灾害监测预警系统，利用卫星遥感和地面监测设备，实时获取泥石流的相关信息，如降雨、水位、地声等，通过数据分析和模型计算，及时发布泥石流灾害预警信息，为防灾减灾工作提供了有力支持。尽管国内外在泥石流基本特征及风险性评价方面取得了众多成果，但仍存在一些不足之处。在基本特征研究中，对于一些特殊类型的泥石流，如冰川泥石流、黄土泥石流等，其形成机制和运动规律尚未完全明确，需要进一步深入研究。不同地区泥石流特征的对比研究相对较少，缺乏对全球范围内泥石流特征的系统性总结和归纳。在风险性评价方面，评价模型和指标体系的通用性和适应性有待提高，不同模型和方法之间的对比和验证工作还不够充分。对于一些不确定性因素，如气候变化、人类活动对泥石流风险性的影响，还难以准确评估和预测。此外，在泥石流灾害的监测和预警技术方面，虽然取得了一定的进展，但仍存在监测精度不够高、预警时效性不强等问题，需要进一步改进和完善。1.3研究内容与方法1.3.1研究内容本研究聚焦黑水县俄瓜沟泥石流，围绕其基本特征与风险性评价展开多维度探索。在基本特征研究板块，从地形地貌、地质构造、水文气象等方面剖析泥石流的形成条件。通过高精度地形图和实地勘测，详细分析俄瓜沟流域的地形起伏、坡度变化、沟谷形态等，明确其地形地貌对泥石流形成的控制作用；深入研究地层岩性、断裂构造等地质构造特征，探究岩石的抗风化能力、节理裂隙发育程度以及地质构造对松散物质的控制和影响；结合长期的气象数据和水文监测资料，分析降水特征、洪水过程以及冰雪融化等水文气象因素在泥石流形成中的诱发作用。针对泥石流的运动特征，运用现场监测和数值模拟相结合的方法，对泥石流的流速、流量、冲击力等关键参数进行研究。在泥石流沟内设置监测点，安装先进的监测设备，如雷达流速仪、超声波流量计等，实时获取泥石流的运动参数；利用专业的数值模拟软件，如FLO-2D、DAN3D等，构建泥石流运动模型，模拟不同工况下泥石流的运动路径、速度变化和堆积范围，深入揭示其运动规律。在堆积特征研究中，对泥石流堆积物的粒度分布、结构特征以及堆积形态进行详细分析。通过现场采样和室内实验，运用激光粒度分析仪、扫描电子显微镜等设备，研究堆积物的粒度组成、分选性和磨圆度等特征；结合地形数据和遥感影像，分析堆积物的堆积形态和分布规律，探讨堆积特征与泥石流运动过程和形成条件之间的内在联系。风险性评价方面，基于泥石流的基本特征和历史灾害数据，构建科学合理的风险性评价指标体系。选取地形地貌、地质构造、水文气象、物源条件、人类工程活动等作为危险性评价指标，综合考虑承灾体的类型、数量、价值以及脆弱性等因素，建立易损性评价指标体系。采用层次分析法（AHP）、模糊综合评价法、灰色关联分析法等多种数学方法，确定各评价指标的权重，构建俄瓜沟泥石流风险性评价模型。运用该模型对俄瓜沟泥石流的风险性进行定量评价，划分风险等级，绘制风险区划图，直观展示不同区域的风险程度。同时，对风险性评价结果进行不确定性分析，评估模型的可靠性和稳定性，为泥石流灾害的防治提供科学、可靠的决策依据。1.3.2研究方法为实现上述研究目标，本研究采用多种研究方法。实地调查法是获取一手资料的重要途径，深入俄瓜沟流域，对泥石流沟的地形地貌、地质构造、物源条件、水文气象等进行详细的现场勘查。利用全站仪、GPS等测量仪器，精确测量沟谷的坡度、高差、沟道长度和宽度等地形参数；观察地层岩性、断裂构造的出露情况，采集岩石样本进行室内分析；调查松散固体物质的分布范围、储量和稳定性，记录降雨、水位等水文气象数据；走访当地居民，了解泥石流的历史发生情况和危害程度。资料收集与整理也是必不可少的环节，广泛收集黑水县的地质、地形、气象、水文等相关资料，包括区域地质调查报告、地形图、气象站数据、水文监测资料等。同时，查阅俄瓜沟泥石流的历史文献、研究报告和灾害记录，对收集到的资料进行系统整理和分析，为后续研究提供数据支持和参考依据。实验分析法则用于对采集的样本进行深入研究，在室内对岩石样本进行物理力学性质测试，如密度、抗压强度、抗剪强度等，分析岩石的工程地质特性；对泥石流堆积物样本进行粒度分析、矿物成分分析和化学成分分析，了解堆积物的组成和结构特征；通过模拟实验，如泥石流水槽实验、降雨入渗实验等，研究泥石流的形成机制和运动规律。在现代技术应用方面，充分利用地理信息系统（GIS）强大的空间分析功能，对地形、地质、气象、水文等多源数据进行整合和分析。通过构建数字高程模型（DEM），提取地形地貌特征参数，如坡度、坡向、地形起伏度等；利用遥感（RS）技术获取大面积的地表信息，解译泥石流沟的分布范围、形态特征和物源分布情况；运用全球定位系统（GPS）进行精确的定位和测量，确保实地调查数据的准确性。通过3S技术的集成应用，实现对俄瓜沟泥石流的动态监测和可视化表达。数值模拟方法在研究泥石流的运动特征和风险性评价中发挥着重要作用，运用FLO-2D、DAN3D等专业数值模拟软件，建立泥石流运动模型和风险性评价模型。通过输入地形、地质、水文气象等参数，模拟泥石流的运动过程，预测其流速、流量、冲击力和堆积范围；利用风险性评价模型，对泥石流的危险性和易损性进行定量评估，分析不同因素对风险性的影响程度，为灾害防治提供科学依据。统计学方法用于对大量的数据进行分析和处理，通过对历史泥石流灾害数据和相关影响因素数据的统计分析，运用概率统计理论，研究泥石流发生的概率和频率，分析各因素与泥石流发生之间的相关性。采用主成分分析、因子分析等多元统计方法，对影响泥石流风险性的多个因素进行降维处理，提取主要影响因子，简化评价指标体系，提高评价模型的精度和可靠性。二、黑水县俄瓜沟泥石流概述2.1地理位置与地质背景黑水县俄瓜沟位于四川省阿坝藏族羌族自治州黑水县双溜索乡，地处青藏高原东南缘横断山脉中段北端的岷江上游，地理坐标介于东经102°35′—103°30′，北纬31°35′—32°38′之间。该区域处于我国地势第一级阶梯向第二级阶梯的过渡地带，地势西北高、东南低，平均海拔3544米，相对高差较大，地形起伏剧烈。俄瓜沟流域内山峦起伏，沟谷纵横，属于典型的高山峡谷地貌。俄瓜沟所在区域的地质构造复杂，处于松潘-甘孜褶皱带与龙门山断裂带的交汇部位，受到强烈的构造运动影响，地层褶皱、断裂发育。区域内主要发育有近南北向和北西向的断裂构造，这些断裂构造不仅控制了山体的隆升和地形地貌的演化，还使得岩石破碎，节理裂隙发育，为泥石流的形成提供了丰富的松散固体物质来源。由于断裂构造的活动，山体稳定性降低，在降雨、地震等因素的触发下，极易发生山体滑坡、崩塌等地质灾害，这些灾害产生的松散物质进一步增加了泥石流的物源储备。地层岩性对泥石流的形成也具有重要影响。俄瓜沟流域出露的地层主要有三叠系、二叠系和侏罗系等。三叠系地层主要为砂板岩、千枚岩等，岩石质地较软，抗风化能力弱，在长期的风化作用下，易形成大量的松散碎屑物质。二叠系地层以石灰岩、白云岩为主，岩石坚硬，但岩溶发育，岩石破碎带和溶蚀裂隙较多，也为松散物质的产生提供了条件。侏罗系地层则主要由砂岩、泥岩组成，泥岩遇水易软化，砂岩抗风化能力相对较强，但在差异风化作用下，也会导致岩石破碎，增加物源量。不同岩性地层的组合和分布，使得流域内的物源条件复杂多样，为泥石流的形成提供了丰富的物质基础。在地形地貌方面，俄瓜沟流域地势陡峭，沟谷狭窄，相对高差大，主沟平均纵坡降为268.6‰。这种陡峭的地形使得水流速度快，侵蚀能力强，能够迅速汇集地表径流，为泥石流的形成提供强大的动力条件。流域内沟谷形态复杂，呈折线型，发育有12道跌坎。跌坎的存在增加了水流的能量损耗和紊动程度，使得水流对沟床和沟岸的冲刷作用增强，促进了松散物质的启动和搬运，有利于泥石流的形成和发展。沟谷两岸山坡坡度多在45°以上，局部地段甚至达到70°，山坡稳定性差，在降雨、地震等因素的作用下，极易发生滑坡、崩塌等地质灾害，为泥石流提供大量的固体物质。俄瓜沟流域内的地形地貌还具有明显的垂直分带性。从沟口到沟源，海拔逐渐升高，地形地貌特征也发生明显变化。沟口附近地势相对较缓，堆积作用较为明显，形成了一定规模的堆积扇；沟道中游地段，地形陡峻，谷深壁峭，以侵蚀作用为主，是泥石流的主要流通区；沟源地区则多为高山峡谷，地势起伏大，积雪冰川发育，为泥石流的形成提供了丰富的水源和物源。这种垂直分带性的地形地貌特征，使得俄瓜沟泥石流的形成和运动过程更加复杂多样。黑水县俄瓜沟特殊的地理位置决定了其复杂的地质背景，地质构造、地层岩性和地形地貌等因素相互作用、相互影响，为泥石流的形成创造了有利条件。深入研究这些地质背景因素，对于理解俄瓜沟泥石流的形成机制、运动特征和堆积规律具有重要意义，也为泥石流灾害的防治提供了关键的地质依据。2.2气候与降水特征黑水县属于季风高原型气候区，其气候具有旱、雨季分明，日照充足，年温差较小，日温差大的显著特点。这种独特的气候特征，对区域内的自然环境和人类活动产生了深远的影响。从年平均气温来看，黑水县年平均气温为9.5℃，最低极端气温可达-14.4℃。较低的年平均气温使得区域内的热量条件相对有限，对农业生产和植被生长的热量需求满足程度有一定限制。例如，在农作物种植方面，一些对热量要求较高的作物难以在此生长，农业种植结构主要以适应冷凉气候的作物为主。然而，较小的年温差使得气温变化相对平稳，减少了因气温骤变对生态系统和人类生活造成的不利影响。而较大的日温差则为水果、蔬菜等农作物的生长提供了有利条件。白天充足的日照和较高的温度有利于作物进行光合作用，积累有机物质；夜晚较低的温度则减少了作物的呼吸作用消耗，使得果实更加香甜，蔬菜口感更佳，提高了农产品的品质。在降水方面，黑水县降雨分布极为不均，夏季降水集中，秋季阴雨连绵。这一降水特点与泥石流的发生密切相关。据统计，黑水县年平均降雨量为620.2毫米，且雨量的日变化明显，夜间降雨占全年总降雨量的70.2%。夏季集中降水，使得大量雨水迅速汇集，形成强大的地表径流。地表径流在山区陡峭的地形条件下，流速极快，对坡面和沟道的冲刷能力极强。它能够迅速携带大量的泥沙、石块等松散固体物质，为泥石流的形成提供了充足的动力和物质条件。以2012年7月9日黑水县双溜索乡俄瓜热十多沟流域发生的泥石流灾害为例，此次灾害正是由突降暴雨引发。短时间内的强降雨使得坡面被暴雨冲刷，洪水揭底冲刷、掏蚀沟岸，大量松散固体物源被洪水携带，从沟道内倾泻而出，最终导致了大规模泥石流的发生。秋季的阴雨连绵也不容忽视。长时间的降雨使得土壤含水量持续饱和，土体的抗剪强度降低，稳定性变差。尤其是在山区，原本就较为脆弱的山体在长时间雨水的浸泡下，更容易发生滑坡、崩塌等地质灾害。这些地质灾害产生的松散物质，进一步丰富了泥石流的物源。同时，持续的降雨还会使得沟道内的水位不断上升，水流能量不断积聚，一旦超过一定的阈值，就可能引发泥石流灾害。从降水的时空分布来看，空间上，受地形地貌和大气环流的影响，不同区域的降水量存在明显差异。山区迎风坡由于地形的抬升作用，降水相对较多；而背风坡则降水较少。这种降水的空间差异，导致了泥石流灾害发生的空间分布也不均匀。在降水较多的区域，泥石流发生的频率相对较高，灾害风险也更大。时间上，降水主要集中在5-10月，这一时期也是泥石流灾害的高发期。随着全球气候变化的影响，极端降水事件的发生频率和强度都有增加的趋势，这进一步加大了黑水县泥石流灾害的风险。例如，近年来，黑水县出现了多次短时强降雨事件，这些强降雨事件往往在短时间内产生大量降水，引发了多起泥石流灾害，对当地的生命财产安全和生态环境造成了严重的破坏。2.3历史泥石流灾害事件回顾俄瓜沟历史上曾多次遭受泥石流灾害的侵袭，这些灾害给当地的生态环境、基础设施和居民生活带来了严重的影响。通过对历史泥石流灾害事件的回顾和分析，能够更深入地了解俄瓜沟泥石流的活动规律和危害程度，为灾害的防治提供宝贵的经验和教训。2012年7月9日，黑水县双溜索乡双溜索村俄瓜热十多沟流域发生了一次大规模的泥石流灾害。此次灾害由突降暴雨引发，短时间内的强降雨使得坡面被暴雨冲刷，洪水揭底冲刷、掏蚀沟岸，携带大量松散固体物源从沟道内倾泻而出。泥石流堆积于沟口并进入黑水河，形成平均高度达2.66m的堰塞坝，一次冲出固体物质约10.28×10⁴m³。沟内地形陡峻，相对高差达2155m，主沟平均纵坡降为268.6‰，纵向上沟谷形态呈折线型，发育12道跌坎，具备爆发泥石流的多元松散物源和势能动力条件，并且管道增大效应明显。该泥石流形成堰塞湖长约1000m，淹没上游部分安置区房屋，冲毁对岸及下游部分房屋，造成一人失踪，沿线乡村公路约1000m被冲毁，省道S302线公路两侧民房不同程度被淹。一旦泥石流堰塞坝溃坝，还将直接威胁黑水河下游两岸居民和双溜索乡镇的安全。此次灾害发生后，当地政府迅速组织力量进行抢险救灾工作，疏散转移受灾群众，抢修受损的基础设施，努力降低灾害造成的损失。同时，相关部门也对此次泥石流灾害进行了深入的调查和分析，总结经验教训，为今后的灾害防治工作提供参考。除了2012年的这次灾害，俄瓜沟在历史上还发生过多次小规模的泥石流灾害。这些灾害虽然规模相对较小，但也对当地的生态环境和居民生活造成了一定的影响。例如，一些小规模的泥石流灾害导致沟道内的树木被冲倒，植被遭到破坏，影响了生态系统的平衡；部分泥石流灾害还冲毁了一些农田和灌溉设施，影响了当地的农业生产。这些小规模的泥石流灾害也提醒我们，即使是规模较小的灾害，也不能掉以轻心，需要加强监测和防范，及时采取有效的措施进行应对。对俄瓜沟历史泥石流灾害事件的回顾，让我们清晰地认识到泥石流灾害的突发性、破坏性和多发性。这些灾害不仅给当地带来了直接的人员伤亡和财产损失，还对生态环境和社会经济发展产生了长期的负面影响。为了减少泥石流灾害的发生和降低其危害程度，我们需要加强对泥石流灾害的研究和监测，提高灾害预警能力，制定科学合理的防治措施，加强宣传教育，提高公众的防灾减灾意识，共同应对泥石流灾害的威胁。三、黑水县俄瓜沟泥石流基本特征3.1泥石流形成条件3.1.1物源条件俄瓜沟泥石流的物源丰富多样，主要来源于山体岩石的风化破碎、崩塌滑坡堆积物以及沟道内的松散沉积物。该区域的地层岩性复杂，主要出露的地层有三叠系、二叠系和侏罗系等。三叠系地层中的砂板岩、千枚岩等岩石质地较软，抗风化能力弱。在长期的风化作用下，岩石表面逐渐破碎，形成大量的碎屑物质。这些碎屑物质在重力和水流的作用下，不断向沟道内汇聚，为泥石流的形成提供了重要的物源基础。例如，在俄瓜沟的上游地区，三叠系地层广泛分布，山坡上随处可见风化破碎的岩石，这些岩石在降雨的冲刷下，很容易被带入沟道，增加了泥石流的物源量。二叠系地层以石灰岩、白云岩为主，虽然岩石本身较为坚硬，但岩溶发育，岩石破碎带和溶蚀裂隙较多。这些破碎带和裂隙使得岩石的整体性遭到破坏，在外界因素的作用下，容易产生崩塌和剥落，从而为泥石流提供松散物质。在俄瓜沟的中游地段，部分沟段两侧的山体为二叠系石灰岩，由于岩溶作用，山体表面布满了溶蚀裂隙和洞穴，岩石破碎严重。在暴雨等强降雨条件下，这些破碎的岩石极易被雨水冲刷下来，进入沟道，成为泥石流的物源。侏罗系地层的砂岩、泥岩组合也对物源的形成有重要影响。泥岩遇水易软化，强度降低，在水流的冲刷和浸泡下，容易发生剥落和坍塌。砂岩虽然抗风化能力相对较强，但在差异风化作用下，与泥岩的接触部位容易形成破碎带，导致岩石破碎。在俄瓜沟的一些区域，侏罗系地层出露，泥岩的软化和砂岩的破碎使得该区域的山坡稳定性较差，为泥石流提供了丰富的物源。除了地层岩性因素外，山体的崩塌和滑坡也是俄瓜沟泥石流物源的重要来源。该区域地质构造活跃，断裂、褶皱发育，地震活动频繁，这些因素导致山体岩体破碎，稳定性降低。在降雨、地震等因素的触发下，极易发生崩塌和滑坡。例如，2012年7月9日俄瓜沟发生的泥石流灾害，就是由突降暴雨引发的山体滑坡和崩塌，大量的滑坡和崩塌堆积物进入沟道，与水流混合，形成了大规模的泥石流。据现场调查，此次灾害中，滑坡和崩塌的堆积物总量达到了数万立方米，极大地丰富了泥石流的物源。沟道内的松散沉积物也是泥石流物源的一部分。在长期的水流侵蚀和搬运作用下，沟道内逐渐堆积了大量的泥沙、砾石等松散物质。这些物质在平时处于相对稳定的状态，但在遇到强降雨等激发条件时，很容易被水流重新启动，参与到泥石流的形成过程中。在俄瓜沟的沟道内，尤其是在沟道的弯曲段和缓坡段，松散沉积物的堆积较为明显，这些区域成为了泥石流物源的重要储备区。此外，人类工程活动也对俄瓜沟泥石流的物源产生了一定的影响。随着区域经济的发展，在俄瓜沟流域内进行了一些公路建设、采矿等工程活动。这些活动破坏了山体的自然植被和岩土结构，产生了大量的弃土、弃渣。例如，在公路建设过程中，开挖山体产生的土石方如果没有得到妥善的处理，就会随意堆积在沟道两侧或沟道内，成为泥石流的潜在物源。在俄瓜沟的部分沟段，由于公路建设和采矿活动的影响，沟道两侧的山坡上布满了弃土和弃渣，这些松散物质在降雨的作用下，很容易被冲刷进入沟道，增加了泥石流发生的风险。3.1.2地形条件俄瓜沟流域的地形条件为泥石流的形成和运动提供了有利的基础。沟谷形态对泥石流的形成和发展具有重要影响。俄瓜沟主沟长8.3km，流域面积达29.10km²，相对高差为2176m。沟谷整体呈折线型，在纵向上发育12道跌坎。这种折线型的沟谷形态使得水流在沟道内的流动过程中不断改变方向，增加了水流的紊动程度和能量损耗。跌坎的存在进一步加剧了水流的能量变化，当水流经过跌坎时，会产生强烈的下蚀和侧蚀作用，对沟床和沟岸进行冲刷，使得沟道内的松散物质更容易被启动和搬运，从而促进了泥石流的形成。在沟道的上游形成区，地形较为开阔，呈三面环山、一面出口的瓢状或漏斗状，这种地形有利于水和碎屑物质的集中汇聚。大量的降水和松散物质在这里聚集，为泥石流的形成提供了充足的物质和能量条件。例如，在暴雨天气下，山坡上的水流迅速向沟道内汇集，携带大量的泥沙、石块等物质，在沟道的开阔区域形成强大的洪流，一旦洪流的能量超过了沟道的承载能力，就会引发泥石流。中游流通区的地形则多为狭窄陡深的峡谷，谷床纵坡降大，主沟平均纵坡降为268.6‰。这种狭窄陡深的峡谷地形使得泥石流在流通区能够迅猛直泻，流速加快。较大的纵坡降为泥石流的运动提供了强大的势能动力，使得泥石流能够携带大量的固体物质快速向下游流动。峡谷两侧的陡壁在泥石流的冲击下，也容易发生崩塌和剥落，进一步增加了泥石流的物源。在2012年的泥石流灾害中，中游流通区的狭窄峡谷地形使得泥石流的流速急剧加快，冲击力增强，对沟道两岸的山体和建筑物造成了严重的破坏。下游堆积区地形为开阔平坦的山前平原或河谷阶地，为泥石流携带的碎屑物提供了堆积场所。当泥石流到达下游堆积区时，由于地形变缓，流速降低，泥石流的搬运能力减弱，携带的固体物质逐渐沉积下来，形成堆积扇或堆积体。这些堆积物在一定程度上改变了下游地区的地形地貌，也可能对下游的居民点、农田和交通设施等造成威胁。在俄瓜沟的下游地区，已经形成了较大规模的泥石流堆积扇，堆积扇上的堆积物厚度较大，颗粒大小混杂，对当地的土地利用和生态环境产生了一定的影响。坡度和高差也是影响泥石流形成和运动的重要地形因素。俄瓜沟流域内山坡坡度多在45°以上，局部地段甚至达到70°，这种陡峭的山坡使得土体和岩石在重力作用下处于不稳定状态，容易发生滑坡和崩塌等地质灾害。大量的滑坡和崩塌堆积物为泥石流提供了丰富的物源。同时，陡峭的山坡也使得水流在山坡上的流速加快，侵蚀能力增强，能够迅速将山坡上的松散物质冲刷到沟道内，促进泥石流的形成。较大的高差使得水流具有较大的势能，在水流从高处向低处流动的过程中，势能转化为动能，为泥石流的运动提供了强大的动力。俄瓜沟相对高差达2155m，这种巨大的高差使得泥石流在运动过程中能够获得较高的速度和能量，具有更强的破坏力。3.1.3水源条件水源是泥石流形成的关键因素之一，它既是泥石流的重要组成部分，又是泥石流的重要激发条件和搬运介质的动力来源。俄瓜沟泥石流的水源主要来自降水、冰雪融水等。降水是俄瓜沟泥石流最主要的水源。黑水县降雨分布不均，夏季降水集中，秋季阴雨连绵。夏季集中的降水，尤其是暴雨，能够在短时间内产生大量的地表径流。这些地表径流在山区陡峭的地形条件下，流速极快，对坡面和沟道的冲刷能力极强。它能够迅速携带大量的泥沙、石块等松散固体物质，为泥石流的形成提供了充足的动力和物质条件。以2012年7月9日黑水县双溜索乡俄瓜热十多沟流域发生的泥石流灾害为例，此次灾害正是由突降暴雨引发。短时间内的强降雨使得坡面被暴雨冲刷，洪水揭底冲刷、掏蚀沟岸，大量松散固体物源被洪水携带，从沟道内倾泻而出，最终导致了大规模泥石流的发生。据气象资料统计，当日该地区的降雨量达到了50mm以上，且降雨强度大，在短时间内形成了强大的地表径流，成为泥石流爆发的直接诱因。秋季的阴雨连绵也不容忽视。长时间的降雨使得土壤含水量持续饱和，土体的抗剪强度降低，稳定性变差。尤其是在山区，原本就较为脆弱的山体在长时间雨水的浸泡下，更容易发生滑坡、崩塌等地质灾害。这些地质灾害产生的松散物质，进一步丰富了泥石流的物源。同时，持续的降雨还会使得沟道内的水位不断上升，水流能量不断积聚，一旦超过一定的阈值，就可能引发泥石流灾害。在一些年份的秋季，由于长时间的阴雨天气，俄瓜沟流域内发生了多起小规模的泥石流灾害，虽然规模相对较小，但也对当地的生态环境和居民生活造成了一定的影响。除了降水，冰雪融水也是俄瓜沟泥石流的重要水源之一。俄瓜沟流域内海拔较高，部分区域存在积雪和冰川。在气温升高的季节，尤其是春季和夏季，积雪和冰川开始融化，形成大量的融水。这些融水在向山下流动的过程中，会与山坡上的松散物质混合，形成泥石流。冰雪融水的产生具有一定的季节性和周期性，其融化量和融化速度受到气温、日照等因素的影响。在气温较高、日照时间长的年份，冰雪融水的量会相对较大，泥石流发生的风险也会相应增加。在一些年份的春季，随着气温的快速回升，俄瓜沟上游地区的积雪大量融化，融水携带大量的泥沙和石块，形成了泥石流，对下游地区的交通和农田造成了一定的破坏。此外，沟道内的地下水在一定条件下也可能参与到泥石流的形成过程中。当地下水位升高，地下水溢出地表时，会增加土体的含水量，降低土体的抗剪强度，使得土体更容易被水流冲刷和搬运。在一些沟段，由于地下水的作用，沟岸的土体变得松软，在降雨或地表径流的作用下，容易发生崩塌和滑坡，为泥石流提供物源。同时，地下水的流动也可能对沟道内的松散物质产生侵蚀和搬运作用，促进泥石流的形成。3.2泥石流运动特征3.2.1流速与流量为准确掌握俄瓜沟泥石流的流速和流量变化规律，本研究综合运用实地观测与数据模拟两种方法。在实地观测方面，于2024年6月20日黑水县芦花镇片区发生短时强降水引发泥石流灾害时，迅速在俄瓜沟泥石流沟内设置多个监测点，安装先进的雷达流速仪和超声波流量计。通过雷达流速仪发射电磁波，利用电磁波与泥石流表面颗粒相互作用产生的多普勒频移，精确测量泥石流的流速。超声波流量计则通过测量超声波在泥石流中的传播时间差，计算出泥石流的流量。在此次观测中，记录到泥石流在形成初期，由于松散固体物质的逐渐加入和水流能量的不断积聚，流速和流量呈现快速上升的趋势。在某监测点，泥石流流速在短时间内从1m/s迅速攀升至5m/s，流量也从初始的5m³/s增加到20m³/s。随着泥石流向下游运动，遇到沟道的转弯和狭窄地段时，流速和流量发生明显变化。在一处沟道转弯处，由于水流受到侧向挤压，流速瞬间增大到8m/s，流量也相应增加到30m³/s，强大的冲击力对沟岸造成了严重的冲刷和破坏。除了实地观测，还运用专业的数值模拟软件FLO-2D进行数据模拟。FLO-2D是一款基于有限差分法的二维泥石流模拟软件，能够考虑地形、物源、水流等多种因素对泥石流运动的影响。在模拟过程中，将俄瓜沟的高精度地形数据、物源分布数据以及不同降雨条件下的水文数据输入到软件中，构建准确的泥石流运动模型。通过多次模拟不同工况下的泥石流运动过程，得到了丰富的流速和流量数据。模拟结果显示，在中等降雨强度下，泥石流在沟道上游的平均流速约为3-4m/s，流量为10-15m³/s；进入中游流通区后，由于地形陡峻，纵坡降大，流速迅速增大到6-8m/s，流量也增加到20-30m³/s；当泥石流到达下游堆积区时，地形变缓，流速逐渐降低到2-3m/s，流量也减小到5-10m³/s。不同区域的流速和流量变化与实地观测结果基本相符，进一步验证了模拟结果的可靠性。通过对实地观测和数据模拟结果的综合分析，发现俄瓜沟泥石流的流速和流量变化受到多种因素的影响。地形条件是影响流速和流量的重要因素之一，沟道的坡度、纵坡降以及沟道的宽窄变化都会导致流速和流量的改变。在坡度较陡、纵坡降大的地段，泥石流能够获得更大的势能，从而加速运动，流速和流量相应增大；而在地形平缓、沟道变宽的区域，泥石流的流速和流量则会减小。物源条件也对流速和流量有显著影响，松散固体物质的含量和颗粒大小会改变泥石流的粘性和流动性。当物源中粗颗粒物质较多时，泥石流的粘性减小，流动性增强，流速和流量可能会增大；反之，当细颗粒物质含量较高时，泥石流的粘性增大，流速和流量则可能会减小。降水强度和历时也是影响流速和流量的关键因素，强降雨能够在短时间内提供大量的水源，增加泥石流的流量，同时也会增强水流的动力，提高流速。3.2.2运动路径与影响范围借助地理信息系统（GIS）强大的空间分析功能，确定俄瓜沟泥石流的运动路径和可能影响范围。首先，收集俄瓜沟流域的高精度数字高程模型（DEM）数据，该数据分辨率达到10m×10m，能够精确反映地形的起伏变化。利用ArcGIS软件对DEM数据进行处理，提取沟道的中心线和边界，生成沟道的矢量数据。通过对沟道矢量数据的分析，结合实地调查结果，确定了泥石流在正常情况下的主要运动路径。在确定运动路径的基础上，运用GIS的水文分析工具，对流域的水流方向、汇流累积量等进行计算。通过水流方向分析，明确了水流在地形作用下的流动方向，从而推断出泥石流可能的运动轨迹。汇流累积量分析则帮助确定了不同区域的水流汇聚情况，为判断泥石流的流量和冲击力提供了依据。结合历史泥石流灾害的现场照片和相关记录，在GIS平台上对泥石流的运动路径进行了可视化表达，直观地展示了泥石流在不同时期的运动轨迹。为了预测泥石流在不同工况下的可能影响范围，采用FLO-2D数值模拟软件与GIS相结合的方法。在FLO-2D中模拟不同降雨强度、物源条件下的泥石流运动过程，得到泥石流在不同时刻的位置和范围数据。将这些模拟结果导入到GIS中，与地形数据、土地利用数据等进行叠加分析。通过叠加分析，可以清晰地看到泥石流可能影响到的区域，包括居民点、农田、公路、铁路等重要设施。在模拟一次强降雨引发的泥石流灾害时，结果显示泥石流可能沿着主沟道向下游运动，影响范围主要集中在沟道两侧200m以内的区域。在该区域内，有多个居民点和部分农田受到威胁。一些位于沟道下游的公路路段也可能被泥石流冲毁，导致交通中断。对于可能受到影响的居民点，根据模拟结果制定了详细的疏散方案，明确了疏散路线和安置地点，以确保居民的生命安全。对于农田和基础设施，提出了相应的防护措施，如修建挡土墙、加固路基等，以减少泥石流造成的损失。3.3泥石流堆积特征3.3.1堆积物组成与结构通过对俄瓜沟泥石流堆积物的现场采样和室内实验分析，发现其堆积物组成复杂多样，结构特征明显。在颗粒组成方面，运用激光粒度分析仪对堆积物样本进行分析，结果显示堆积物的粒径范围广泛，从小于0.002mm的粘粒到大于200mm的巨砾均有分布。其中，粒径小于2mm的细颗粒物质（粉粒和粘粒）含量约占30%-40%，粒径在2-60mm之间的砂粒和砾石含量约占40%-50%，粒径大于60mm的粗颗粒物质（块石和巨砾）含量约占10%-20%。这种宽粒径分布特征表明泥石流在运动过程中具有很强的搬运能力，能够携带各种大小的颗粒物质。在2012年泥石流灾害的堆积物中，发现了大量直径超过1m的巨砾，这些巨砾在泥石流的强大动力作用下被搬运到下游堆积区，充分体现了泥石流的搬运能力。分选性是衡量堆积物颗粒均匀程度的重要指标。俄瓜沟泥石流堆积物的分选性极差，分选系数大于2.0。这是因为泥石流是一种快速流动的特殊洪流，在运动过程中，不同粒径的颗粒物质被强烈搅拌混合，没有足够的时间进行分选。分选性差的堆积物结构较为松散，孔隙率较大，这也使得堆积物的稳定性相对较差。在堆积区，由于堆积物的分选性差，容易受到水流和风力的侵蚀作用，导致堆积物的再次搬运和流失。磨圆度反映了颗粒在搬运过程中受到磨损的程度。俄瓜沟泥石流堆积物的颗粒磨圆度普遍较差，多呈棱角状和次棱角状。这是因为泥石流的搬运距离相对较短，颗粒在搬运过程中与其他物体的摩擦和碰撞时间较短，磨损程度较小。与河流沉积物相比，河流中的颗粒经过长时间的搬运和磨蚀，磨圆度较好，多呈圆形或椭圆形。而俄瓜沟泥石流堆积物的棱角状和次棱角状颗粒表明其形成过程较为迅速，搬运过程中的磨蚀作用不明显。堆积物的结构特征还包括颗粒的排列方式和孔隙特征。通过扫描电子显微镜观察发现，堆积物中的颗粒排列杂乱无章，没有明显的定向性。这是由于泥石流在堆积过程中，颗粒受到复杂的水流和动力作用，随机堆积在一起。堆积物的孔隙大小不一，分布不均匀，大孔隙主要由粗颗粒之间的空隙形成，小孔隙则存在于细颗粒之间。孔隙率的大小对堆积物的透水性和稳定性有重要影响，较大的孔隙率使得堆积物的透水性较强，但也降低了堆积物的抗剪强度，增加了滑坡和坍塌的风险。堆积物的组成与结构特征反映了泥石流的形成过程。丰富的物源条件使得泥石流能够携带各种粒径的颗粒物质，地形条件决定了泥石流的搬运距离和动力大小，从而影响颗粒的分选性和磨圆度。分选性差、磨圆度低的堆积物特征表明泥石流形成过程迅速，搬运过程中颗粒没有充分的时间进行分选和磨蚀。堆积物的结构特征也与泥石流的运动和堆积方式密切相关，杂乱的颗粒排列和不均匀的孔隙分布是泥石流快速堆积的结果。3.3.2堆积形态与分布规律俄瓜沟泥石流在沟口和下游地区形成了独特的堆积形态，其分布规律受到多种因素的影响。在堆积形态方面，通过现场调查和遥感影像分析，发现泥石流在沟口主要形成了扇形堆积体，即泥石流堆积扇。堆积扇的平面形态呈扇形，从沟口向四周展开，半径可达数百米。堆积扇的纵剖面呈上凸下凹的弧形，顶部较陡，坡度可达10°-15°，向扇缘逐渐变缓，坡度一般在2°-5°。堆积扇的横剖面则呈现出中间厚、两侧薄的特点，中心部位的堆积物厚度可达数米甚至数十米，向两侧逐渐变薄。在2012年泥石流灾害后，形成的堆积扇面积较大，覆盖了沟口附近的大片区域，对当地的土地利用和交通造成了严重影响。堆积扇的形成与泥石流的运动特征密切相关。当泥石流冲出沟口后，由于地形突然开阔，流速迅速降低，搬运能力减弱，携带的大量固体物质逐渐沉积下来。在堆积过程中，粗颗粒物质首先沉积在扇顶，形成堆积扇的骨架，细颗粒物质则随着水流向扇缘搬运，在扇缘逐渐沉积，从而形成了上粗下细的堆积结构。堆积扇的坡度和形态也受到泥石流流量、流速以及地形条件的影响。流量大、流速快的泥石流形成的堆积扇坡度较陡，面积较小；而流量小、流速慢的泥石流形成的堆积扇坡度较缓，面积较大。地形条件也会对堆积扇的形态产生影响，在沟口附近地形较为平坦的区域，堆积扇能够较为完整地向四周展开；而在地形有起伏或受到障碍物阻挡的区域，堆积扇的形态会发生变形，可能出现偏态堆积或局部堆积加厚的情况。在下游地区，泥石流的堆积物沿着河谷呈条带状分布。这是因为下游河谷地形相对狭窄，泥石流在河谷中流动时受到河谷两岸的约束，只能沿着河谷方向延伸。堆积物的厚度和宽度在不同地段有所变化，一般在河谷的转弯处和河谷变宽的地段，堆积物厚度较大，宽度也较宽。这是因为在这些地段，泥石流的流速降低，搬运能力减弱，堆积作用增强。在河谷的直段和河谷变窄的地段，堆积物厚度相对较薄，宽度也较窄。泥石流堆积物的分布规律还受到地形地貌、水流条件和物源供应的影响。地形地貌决定了泥石流的流动路径和堆积场所，水流条件影响着泥石流的流速和搬运能力，物源供应则决定了堆积物的数量和规模。在俄瓜沟流域，地势较高的区域，泥石流的堆积物相对较少，主要是因为这些区域物源相对较少，且水流速度较快，不利于堆积。而在地势较低的沟口和下游河谷地区，物源丰富，水流速度降低，有利于堆积物的沉积，因此堆积物分布较为广泛。水流条件的变化也会导致堆积物分布的差异，当水流流量较大时，泥石流能够搬运更远的距离，堆积物的分布范围也会相应扩大；当水流流量较小时，堆积物则主要集中在沟口附近。物源供应的变化同样会影响堆积物的分布，物源丰富时，堆积物的数量和规模都会增加，分布范围也会更广；物源减少时，堆积物的数量和规模也会相应减小，分布范围也会缩小。四、黑水县俄瓜沟泥石流风险性评价指标体系构建4.1风险性评价的概念与内涵泥石流风险性评价是对泥石流灾害可能造成的损失进行综合评估的过程，旨在量化泥石流灾害发生的可能性及其对人类社会和自然环境造成的潜在危害程度。其核心内涵在于将泥石流灾害的危险性、易损性以及承灾体的特征相结合，全面分析泥石流灾害的风险状况，为灾害防治和风险管理提供科学依据。泥石流灾害的危险性是风险性评价的重要组成部分，它反映了泥石流发生的可能性及其潜在的破坏能力。危险性的大小取决于泥石流的形成条件和运动特征，包括物源条件、地形条件、水源条件、流速、流量、运动路径等因素。丰富的物源、陡峭的地形和充足的水源会增加泥石流发生的可能性和规模，从而提高其危险性；流速快、流量大的泥石流具有更强的破坏力，危险性也更高。在黑水县俄瓜沟，复杂的地层岩性、活跃的地质构造以及特殊的地形地貌，使得该区域具备了丰富的物源和有利于泥石流形成的地形条件，加上降水集中的气候特征，导致泥石流发生的危险性较高。2012年7月9日俄瓜沟发生的大规模泥石流灾害，就是在强降雨的触发下，丰富的物源和陡峭的地形条件相互作用的结果，充分体现了该区域泥石流的高危险性。易损性则关注承灾体在面对泥石流灾害时的脆弱程度，即承灾体受到泥石流破坏的可能性和损失程度。易损性受到承灾体自身的特性以及其所处环境的影响。不同类型的承灾体，如居民点、农田、基础设施等，由于其结构、功能和价值的差异，对泥石流灾害的抵抗能力各不相同。居民点中的建筑物结构坚固程度、人口密度，农田的土壤质量、农作物种类，基础设施的材质、设计标准等，都会影响其易损性。承灾体所处的地理位置、周边地形以及防护措施的完善程度等环境因素，也会对易损性产生重要影响。位于泥石流沟下游、地势较低且缺乏有效防护措施的居民点，在面对泥石流灾害时，易损性较高，遭受破坏的可能性和损失程度也更大。承灾体是指可能受到泥石流灾害影响的人类社会和自然环境的各个要素，包括人口、建筑物、基础设施、农田、生态系统等。人口的数量、分布和年龄结构等因素，会影响到人员在泥石流灾害中的伤亡风险；建筑物的类型、结构和建设标准，决定了其在泥石流冲击下的稳定性；基础设施如公路、桥梁、电力设施等的损坏，会对区域的交通、能源供应等造成严重影响，进而影响社会经济的正常运转；农田的受损会影响农业生产，威胁粮食安全；生态系统的破坏则会对生物多样性、水土保持等产生长期的负面影响。在黑水县俄瓜沟流域，居民点、农田、省道S302线公路等都是重要的承灾体。2012年的泥石流灾害冲毁了部分居民点和农田，导致居民财产损失和农业生产受阻，省道S302线公路两侧民房不同程度被淹，交通中断，对当地的社会经济造成了严重影响。4.2评价指标选取原则在构建黑水县俄瓜沟泥石流风险性评价指标体系时，评价指标的选取遵循了一系列科学合理的原则，以确保评价结果的准确性和可靠性。科学性原则是评价指标选取的首要原则。要求所选取的指标能够客观、准确地反映泥石流风险性的本质特征，具有明确的科学内涵和理论依据。指标的定义、计算方法和数据来源都必须科学合理，经得起理论和实践的检验。在选取地形条件相关指标时，坡度、高差、沟谷纵坡降等指标的确定都是基于地貌学和工程地质学的理论基础，这些指标能够准确地反映地形对泥石流形成和运动的影响。通过对坡度的测量和分析，可以判断山坡的稳定性，坡度越大，山坡越不稳定，越容易发生滑坡和崩塌，从而为泥石流提供物源；高差和沟谷纵坡降则直接影响水流的势能和流速，进而影响泥石流的运动能力。代表性原则强调选取的指标应具有广泛的代表性，能够全面反映泥石流风险性的各个方面。从泥石流的形成条件、运动特征到堆积特征，以及承灾体的易损性等多个角度选取指标，确保评价指标体系的完整性。在物源条件方面，选取了岩石风化程度、崩塌滑坡堆积物储量等指标，这些指标能够代表不同类型的物源，全面反映物源对泥石流风险性的影响。岩石风化程度反映了岩石在长期风化作用下的破碎程度，风化程度越高，岩石越破碎，为泥石流提供的物源就越多；崩塌滑坡堆积物储量则直接体现了因崩塌和滑坡产生的松散物质的数量，储量越大，泥石流发生时可利用的物源就越丰富。可操作性原则是评价指标选取的重要原则之一。要求选取的指标数据易于获取、测量和计算，在实际应用中具有可行性。优先选择那些可以通过实地调查、监测数据或现有资料获取的指标，避免选取过于复杂或难以获取数据的指标。在气象条件指标选取中，选择年平均降雨量、最大日降雨量等易于从气象站获取的指标，这些指标能够较好地反映降水对泥石流的影响，同时数据获取相对容易。通过气象站的长期监测记录，可以直接获取这些指标的数据，为泥石流风险性评价提供可靠的数据支持。独立性原则要求各个评价指标之间应相互独立，避免指标之间存在过多的信息重叠。这样可以提高评价指标体系的效率和准确性，避免因指标相关性过高而导致评价结果的偏差。在选取指标时，通过相关性分析等方法，对指标之间的相关性进行检验，确保所选指标之间的相关性较低。在地质构造指标选取中，选择断裂密度和褶皱强度作为不同的指标，它们分别从不同角度反映地质构造对泥石流的影响，断裂密度主要反映断裂构造的发育程度，而褶皱强度则体现了地层褶皱的强烈程度，两者之间独立性较强，能够更全面地反映地质构造对泥石流风险性的影响。动态性原则考虑到泥石流风险性会随着时间和环境的变化而变化，评价指标也应具有一定的动态性。及时更新和调整指标，以适应不同时期和不同环境条件下的泥石流风险性评价。随着全球气候变化和人类工程活动的加剧，泥石流的形成条件和风险性也在不断变化。在评价指标体系中，应适时增加一些反映气候变化和人类活动影响的指标，如气温变化趋势、土地利用变化等。通过对这些动态指标的监测和分析，可以及时了解泥石流风险性的变化情况，为灾害防治提供更及时、准确的决策依据。4.3具体评价指标确定4.3.1危险性指标在构建黑水县俄瓜沟泥石流风险性评价指标体系时，危险性指标的选取至关重要，它直接反映了泥石流发生的可能性及其潜在的破坏能力。泥石流规模是一个关键的危险性指标，它主要通过一次泥石流冲出的固体物质总量来衡量。规模越大，泥石流携带的固体物质越多，其冲击力和破坏力也就越强。在2012年7月9日俄瓜沟发生的泥石流灾害中，一次冲出固体物质约10.28×10⁴m³，如此巨大的规模使得泥石流在运动过程中能够摧毁沿途的一切障碍物，对沟道两岸的山体、建筑物和基础设施造成了严重的破坏。大量的固体物质堆积在沟口和下游地区，形成了大面积的堆积扇，改变了当地的地形地貌，对土地利用和交通产生了长期的影响。发生频率也是衡量泥石流危险性的重要指标之一。俄瓜沟泥石流的发生频率相对较高，这与该区域的地质、地形和气候条件密切相关。复杂的地质构造、陡峭的地形以及降水集中的气候特征，使得该区域具备了泥石流频繁发生的条件。较高的发生频率意味着当地居民和基础设施面临着更大的威胁，增加了灾害发生的不确定性和风险。频繁发生的泥石流会对当地的生态环境造成持续的破坏，影响植被的生长和恢复，导致水土流失加剧，进一步恶化地质条件，形成恶性循环。流速和流量同样是重要的危险性指标。流速反映了泥石流运动的速度，流量则表示单位时间内通过某一断面的泥石流体积。流速和流量越大，泥石流的动能和势能就越大，其破坏力也就越强。在2024年6月20日黑水县芦花镇片区发生的泥石流灾害中，通过实地观测和数值模拟，记录到泥石流在形成初期，流速在短时间内从1m/s迅速攀升至5m/s，流量也从初始的5m³/s增加到20m³/s。在沟道的转弯和狭窄地段，流速瞬间增大到8m/s，流量也相应增加到30m³/s，强大的冲击力对沟岸造成了严重的冲刷和破坏。高流速和大流量的泥石流能够迅速冲毁房屋、桥梁等建筑物，掩埋农田，阻断交通，对当地的社会经济造成巨大的损失。地形坡度和高差也对泥石流的危险性有重要影响。俄瓜沟流域内山坡坡度多在45°以上，局部地段甚至达到70°，相对高差达2155m。陡峭的山坡使得土体和岩石在重力作用下处于不稳定状态，容易发生滑坡和崩塌等地质灾害，为泥石流提供丰富的物源。较大的高差则使得水流具有较大的势能，在水流从高处向低处流动的过程中，势能转化为动能，为泥石流的运动提供强大的动力。在这样的地形条件下，泥石流一旦发生，就具有很强的突发性和破坏力，难以预测和防范。松散固体物质储量也是危险性评价的重要指标。俄瓜沟泥石流的松散固体物质储量丰富，主要来源于山体岩石的风化破碎、崩塌滑坡堆积物以及沟道内的松散沉积物。丰富的松散固体物质为泥石流的形成提供了充足的物质基础，增加了泥石流发生的可能性和规模。在2012年的泥石流灾害中，山体滑坡和崩塌产生的大量堆积物进入沟道，与水流混合，形成了大规模的泥石流。这些松散固体物质在泥石流的运动过程中，不断被搬运和堆积，进一步扩大了泥石流的影响范围和破坏程度。4.3.2易损性指标易损性指标的确定对于准确评估黑水县俄瓜沟泥石流的风险性至关重要，它主要反映了承灾体在面对泥石流灾害时的脆弱程度。受威胁人口数量是一个关键的易损性指标。在俄瓜沟流域内，有多个居民点分布在泥石流沟道附近，受威胁人口数量较多。这些居民点的房屋大多为普通砖混结构，抗震、抗泥石流冲击能力较弱。一旦发生泥石流灾害，受威胁人口的生命安全将受到严重威胁。2012年的泥石流灾害导致部分居民点房屋被冲毁，居民财产遭受重大损失，甚至造成一人失踪。大量受威胁人口的存在，不仅增加了灾害发生时人员伤亡的风险，也加大了灾害救援和人员疏散的难度。在灾害发生时，如何快速、安全地疏散大量受威胁人口，是防灾减灾工作面临的一个重要挑战。建筑物密度也是衡量易损性的重要指标之一。俄瓜沟流域内一些区域的建筑物密度较大，尤其是在居民集中居住的区域。建筑物密度大意味着在泥石流灾害发生时，更多的建筑物将受到冲击和破坏，造成更大的财产损失。密集的建筑物还会阻碍泥石流的流动，导致泥石流在局部区域堆积，进一步加剧灾害的危害程度。在一些建筑物密度较大的地段，泥石流发生时，建筑物之间的空间狭小，使得泥石流的冲击力无法有效分散，导致建筑物更容易被冲垮，造成严重的人员伤亡和财产损失。基础设施价值是易损性评价的重要方面。俄瓜沟流域内有省道S302线公路等重要基础设施，这些基础设施对于当地的经济发展和居民生活至关重要。省道S302线公路是连接黑水县与外界的重要交通通道，一旦该公路因泥石流灾害受损，将导致交通中断，影响物资运输和人员往来，对当地的经济发展造成严重影响。公路两侧的桥梁、涵洞等设施在泥石流的冲击下也容易受损，修复这些基础设施需要耗费大量的人力、物力和财力。一些电力、通信等基础设施在泥石流灾害中也可能受到破坏，导致区域停电、通信中断，给居民的生活带来极大不便。土地利用类型对易损性也有显著影响。该区域内的土地利用类型主要包括农田、林地和建设用地等。农田是当地农业生产的重要基础，一旦遭受泥石流破坏，将导致农作物受损，影响农业产量，威胁当地的粮食安全。2012年的泥石流灾害冲毁了部分农田，使得当年的农作物减产，给农民带来了经济损失。林地在保持水土、涵养水源方面发挥着重要作用，但在泥石流灾害中，林地也可能遭受破坏，导致植被覆盖率下降，水土流失加剧，进一步恶化生态环境。建设用地的易损性则主要体现在建筑物和基础设施的破坏上，不同类型的建设用地，如商业用地、工业用地和居住用地等，其易损程度也有所不同。商业用地和工业用地通常集中了大量的商业设施和工业设备，一旦遭受泥石流灾害，不仅会造成建筑物的破坏，还会导致商业活动和工业生产的中断，带来巨大的经济损失。4.3.3承灾体指标承灾体指标在黑水县俄瓜沟泥石流风险性评价中具有关键作用，它全面反映了可能受到泥石流灾害影响的人类社会和自然环境的各个要素的特征。人口密度是一个重要的承灾体指标。在俄瓜沟流域，部分居民点的人口密度较大，这使得在泥石流灾害发生时，单位面积内受到威胁的人口数量增多。一旦泥石流来袭，人员疏散难度加大，伤亡风险显著提高。例如，在一些人口密集的居民点，狭窄的街道和有限的疏散空间，使得居民在灾害发生时难以迅速撤离，增加了被泥石流掩埋或冲击的危险。较大的人口密度还会导致救援工作的难度增加，救援队伍和物资难以快速抵达受灾区域，影响救援效率。经济发展水平也对承灾体的应对能力产生重要影响。俄瓜沟所在的黑水县整体经济发展水平相对较低，这限制了当地在泥石流灾害防治方面的投入能力。在灾害预防阶段，经济实力不足使得难以建设完善的监测预警系统和有效的防护工程。由于缺乏资金，无法安装先进的泥石流监测设备，不能及时准确地掌握泥石流的发生迹象，导致预警信息的发布存在延迟，影响居民的提前防范。在灾害发生后的救援和恢复阶段，经济发展水平低意味着可用于救灾和恢复重建的资金有限，这将延缓受灾地区的恢复速度，影响居民的生活和生产秩序的尽快恢复。有限的资金使得在修复受损的基础设施、重建房屋等方面面临困难，受灾群众可能需要长时间才能恢复正常生活。建筑物结构类型是影响承灾体脆弱性的关键因素。俄瓜沟流域内的建筑物结构类型多样，包括砖混结构、木结构和简易结构等。其中，木结构和简易结构的建筑物抗泥石流冲击能力较弱，在泥石流的强大冲击力下，很容易倒塌或被冲毁，对居民的生命安全构成严重威胁。在2012年的泥石流灾害中，一些木结构和简易结构的房屋被泥石流瞬间冲垮，导致屋内居民伤亡。相比之下，砖混结构的建筑物在一定程度上具有较好的抗冲击能力，但如果结构设计不合理或施工质量不达标，也难以抵御大规模泥石流的冲击。建筑物的抗震性能、墙体的强度和稳定性等因素，都会影响其在泥石流灾害中的受损程度。生态系统的稳定性同样是重要的承灾体指标。俄瓜沟流域的生态系统较为脆弱，森林覆盖率相对较低，植被的水土保持能力较弱。一旦发生泥石流灾害，生态系统的稳定性将遭到破坏，导致水土流失加剧，进一步恶化地质条件，形成恶性循环。泥石流会冲毁大量的植被，破坏土壤结构，使得土壤的保水保肥能力下降，影响植物的生长和恢复。生态系统的破坏还会对生物多样性造成威胁，许多动植物失去了栖息地，导致物种数量减少，生态平衡被打破。生态系统的不稳定也会增加泥石流再次发生的风险，因为植被的减少使得山坡更容易受到雨水的冲刷，为泥石流的形成提供了更多的物源和动力条件。五、黑水县俄瓜沟泥石流风险性评价方法与模型5.1常用风险性评价方法概述在泥石流风险性评价领域，多种方法被广泛应用，每种方法都有其独特的原理、优势和适用场景。层次分析法（AHP）是一种定性与定量相结合的多准则决策分析方法，由美国运筹学家萨蒂（T.L.Saaty）教授于20世纪70年代初期提出。该方法的核心在于将复杂问题分解为若干层次和因素，通过各因素之间的两两比较和计算，得出不同方案的权重，为决策者提供定量化的决策依据。在构建层次结构模型时，需明确决策问题，并将其分解为目标层、准则层和方案层。目标层是决策的最终目的，准则层包含影响决策的各种因素，方案层则是可供选择的方案或策略。在确定黑水县俄瓜沟泥石流风险性评价指标权重时，可将泥石流风险性评价作为目标层，将危险性指标、易损性指标和承灾体指标作为准则层，各具体指标作为方案层。在准则层，决策者需要根据各准则之间的相对重要性进行两两比较，并采用1-9标度法给出定量的判断值，这些值构成一个判断矩阵。若认为危险性指标比易损性指标略为重要，则在相应矩阵位置填入3，反之则填入1/3。为保证判断矩阵的一致性，需要计算一致性指标CI和一致性比率CR，CR=CI/RI，其中RI是随机一致性指数，与矩阵阶数有关，通常CR值小于0.1，则认为判断矩阵的一致性是可以接受的。通过求解判断矩阵的最大特征值所对应的特征向量，经过归一化处理后，得到各准则的权重向量。层次分析法能够将人的主观判断用数量形式表达和处理，有效减少个人主观臆断所带来的弊端，使评价结果更具可信度，适用于具有多个层次的综合评价中。但该方法也存在一定局限性，判断矩阵的构建依赖专家经验，主观性较强，当因素较多时，判断矩阵的一致性检验难度较大。模糊综合评价法（FCEM）于1965年由美国自动控制专家查德（L．A．Zadeh）教授提出，是一种基于模糊数学的综合评价方法。该方法根据模糊数学的隶属度理论把定性评价转化为定量评价，对于受到多种因素制约的事物或对象做出一个总体的评价。在对俄瓜沟泥石流风险性进行评价时，首先要确定评价因素集，将影响泥石流风险性的各种因素构成一个集合U，如U={泥石流规模，发生频率，流速，流量，地形坡度，高差，松散固体物质储量，受威胁人口数量，建筑物密度，基础设施价值，土地利用类型，人口密度，经济发展水平，建筑物结构类型，生态系统的稳定性}。确定评价等级集，将评价结果划分成不同的等级，构成集合V，如V={高风险，较高风险，中等风险，较低风险，低风险}。确定模糊关系矩阵R，对每个评价因素ui，评估其在各个评价等级vj上的隶属度，形成模糊关系矩阵R，矩阵R的元素rij表示第i个评价因素属于第j个评价等级的隶属度。通过专家打分或利用隶属函数计算等方式确定各评价因素的权重A，满足ai相加为1。通过模糊合成运算B=AR，得到综合评价结果B，B是一个1m的行向量，其元素bj表示评价对象属于第j个评价等级的综合隶属度。模糊综合评价法能较好地处理模糊性和不确定性问题，使评价结果更接近实际情况，结果清晰，系统性强，适用于解决各种非确定性问题。然而，该方法计算相对复杂，对指标权重矢量的确定主观性较强，当指标集U较大时，在权矢量和为1的条件约束下，相对隶属度权系数往往偏小，权矢量与模糊矩阵R不匹配，可能会出现超模糊现象，分辨率很差，无法区分谁的隶属度更高，甚至造成评判失败。灰色关联分析法是一种多因素统计分析方法，它是以各因素的样本数据为依据用灰色关联度来描述因素间关系的强弱、大小和次序。其基本思想是根据序列曲线几何形状的相似程度来判断其联系是否紧密，曲线越接近，相应序列之间的关联度就越大，反之就越小。在黑水县俄瓜沟泥石流风险性评价中，将泥石流风险性作为参考序列，将各评价指标数据作为比较序列。先对数据进行无量纲化处理，消除量纲和数量级的影响。计算各比较序列与参考序列的关联系数，关联系数反映了某一时刻比较序列与参考序列的关联程度。通过对关联系数进行加权平均，得到各评价指标与泥石流风险性的关联度。关联度越大，说明该指标对泥石流风险性的影响越大。灰色关联分析法对数据要求较低，可在数据量少、信息不完全的情况下进行分析，能有效处理复杂系统中各因素之间的非线性关系。但该方法在确定各因素的权重时，通常采用主观赋值法，存在一定主观性，且关联度的计算结果受数据变换方法的影响较大。5.2本研究采用的评价方法与模型选择综合考虑黑水县俄瓜沟泥石流的特点以及数据的可获取性，本研究选用层次分析法（AHP）确定指标权重，运用模糊综合评价法进行风险性评价。层次分析法能够将复杂的决策问题分解为有序层次结构，通过数学方法进行定量分析，有效处理多准则决策问题。在确定俄瓜沟泥石流风险性评价指标权重时，它可以将泥石流风险性评价作为目标层，危险性指标、易损性指标和承灾体指标作为准则层，各具体指标作为方案层，通过各准则之间的两两比较和计算，得出各指标的权重。这种方法将人的主观判断用数量形式表达和处理，减少了个人主观臆断的影响，使评价结果更具可信度，适用于具有多个层次的综合评价，与本研究构建的多层次评价指标体系相契合。模糊综合评价法基于模糊数学的隶属度理论，把定性评价转化为定量评价，对于受到多种因素制约的事物或对象做出一个总体的评价。在对俄瓜沟泥石流风险性进行评价时，该方法能充分考虑到影响泥石流风险性的各种因素的模糊性和不确定性，使评价结果更接近实际情况。通过确定评价因素集、评价等级集、模糊关系矩阵以及权重集，进行模糊合成运算，最终得到综合评价结果。其结果清晰，系统性强，能较好地解决模糊的、难以量化的问题，适合本研究中对泥石流风险性这种具有不确定性问题的评价。将层次分析法与模糊综合评价法相结合，既能通过层次分析法合理确定各评价指标的权重，又能利用模糊综合评价法处理评价过程中的模糊性和不确定性，从而构建出科学、合理的黑水县俄瓜沟泥石流风险性评价模型，为泥石流灾害的防治提供可靠的决策依据。5.3评价模型的构建与参数确定5.3.1层次分析法确定权重运用层次分析法（AHP）确定黑水县俄瓜沟泥石流风险性评价指标的权重，具体步骤如下：构建层次结构模型：将黑水县俄瓜沟泥石流风险性评价这一复杂问题，分解为清晰的目标层、准则层和指标层。目标层明确为黑水县俄瓜沟泥石流风险性评价，它代表了整个评价的最终指向。准则层包含危险性指标、易损性指标和承灾体指标，这些准则从不同维度影响着泥石流风险性。指标层则由泥石流规模、发生频率、流速、流量、地形坡度、高差、松散固体物质储量、受威胁人口数量、建筑物密度、基础设施价值、土地利用类型、人口密度、经济发展水平、建筑物结构类型、生态系统的稳定性等具体指标构成，它们是对准则层各因素的进一步细化和量化。通过这样的层次结构构建，将复杂的评价问题转化为有序的层次体系，便于后续的分析和计算。建立判断矩阵：在准则层和指标层，邀请地质、气象、地理信息系统等领域的10位资深专家，依据各准则及指标之间的相对重要性进行两两比较。采用1-9标度法来赋予定量的判断值，进而构成判断矩阵。1-9标度法中，1表示两个因素具有同等重要性；3表示一个因素比另一个因素略重要；5表示一个因素比另一个因素明显重要；7表示一个因素比另一个因素强烈重要；9表示一个因素比另一个因素极端重要；2、4、6、8则是上述相邻判断的中间值。若认为危险性指标比易损性指标明显重要，那么在判断矩阵中对应的位置就填入5，反之则填入1/5。通过这种方式，将专家的主观判断转化为具体的数值，为后续的权重计算提供数据基础。计算权重向量：利用方根法计算判断矩阵的特征向量和最大特征根。以准则层判断矩阵为例，先计算判断矩阵每列元素的乘积，再对乘积开n次方（n为判断矩阵的阶数），得到一个列向量。将该列向量归一化处理，即每个元素除以列向量元素之和，得到的结果就是各准则的权重向量。计算出的危险性指标权重为0.52，易损性指标权重为0.30，承灾体指标权重为0.18，这表明在泥石流风险性评价中，危险性指标的影响最为显著，易损性指标次之，承灾体指标相对较小。在指标层，同样按照方根法计算各指标的权重，得到泥石流规模权重为0.18，发生频率权重为0.15，流速权重为0.12等。一致性检验：为确保判断矩阵的一致性，计算一致性指标CI和一致性比率CR。CI=（λmax-n）/（n-1），其中λmax是判断矩阵的最大特征根，n为判断矩阵的阶数。CR=CI/RI，RI是随机一致性指数，与矩阵阶数有关。当CR值小于0.1时，认为判断矩阵的一致性是可以接受的。对于准则层判断矩阵，计算得到的CR值为0.05，小于0.1，说明专家的判断具有较好的一致性，权重计算结果可靠。通过一致性检验，可以保证层次分析法确定的权重具有合理性和科学性，为后续的模糊综合评价提供准确的权重数据。5.3.2模糊综合评价法确定隶属度采用模糊综合评价法对黑水县俄瓜沟泥石流风险性进行评价，具体步骤如下：确定评价因素集：将影响黑水县俄瓜沟泥石流风险性的各种因素构成评价因素集U，即U={泥石流规模，发生频率，流速，流量，地形坡度，高差，松散固体物质储量，受威胁人口数量，建筑物密度，基础设施价值，土地利用类型，人口密度，经济发展水平，建筑物结构类型，生态系统的稳定性}。这些因素全面涵盖了泥石流形成条件、运动特征、承灾体特征等方面，是进行风险性评价的重要依据。确定评价等级集：将评价结果划分为不同的等级，构成评价等级集V，本研究将其分为高风险、较高风险、中等风险、较低风险、低风险五个等级，即V={高风险，较高风险，中等风险，较低风险，低风险}。这样的等级划分便于直观地展示泥石流风险性的程度差异，为风险管理和决策提供明确的参考。确定模糊关系矩阵R：对于每个评价因素ui，通过专家打分和隶属函