鲁甸“8.3”地震后红石岩堰塞湖危岩体：特征剖析与防治策略

鲁甸“8.3”地震后红石岩堰塞湖危岩体：特征剖析与防治策略一、绪论1.1研究背景与意义2014年8月3日16时30分，云南省昭通市鲁甸县发生6.5级地震，震源深度12km，震中位于昭通市鲁甸县龙头山镇、火德红镇和巧家县包谷垴乡附近，此次地震震级高、震源浅，破坏力大，地震最高烈度为Ⅸ度。鲁甸“8.3”地震震中地质条件复杂，山高谷深，余震不断，加之降雨影响，产生了许多地质灾害并在牛栏江干流等处形成堰塞湖。其中，最大的红石岩堰塞体位于红石岩水电站取水坝下游约600m处，堰塞体最大高度102m，顺河向底部长约910m，顶部长约262m，垂直河流方向平均顶宽301m，总体呈不对称的马鞍状，鞍部高程1222.00m，堰塞体总方量约1200万m³，堆积物系右岸基岩（白云岩和灰岩）崩滑形成，碎块石含量占85%以上，粒径大部分在0.5-5.0m之间。堰塞体形成后，湖水位17h上涨24.50m，涨幅1.42m/h，红石岩水电站取水坝完全被淹没。照此上涨速度，2d内可能漫顶，形成集水面积11832km²、库容2.6亿m³、回水长度25km的大型堰塞湖，直接影响上下游3万余人，下游天花板、黄角树等水电站安全，严重威胁周边4县区3万多群众的生命财产安全。灾情就是命令，经过解放军、武警和民兵预备役部队1100多名官兵连续奋战9昼夜，火速打通排水泄洪通道，险情基本排除，下游群众暂告安全。当年10月4号，红石岩堰塞湖后续处置任务完成，险情进一步缓解，但是当时的这个防洪渡汛，只能够满足两年一遇的渡汛标准，对下游的老百姓生命财产安全依然是一个巨大的威胁。由于堰塞体方量达1000万立方米，堆积体高约102米，加上周围交通不便，堰塞湖难以通过拆除方式恢复自然河流状态。在山体崩塌后，堰塞湖右岸形成600米高的负角度陡峭岩壁，危岩体不断坍塌和崩落，而且堰塞体构成复杂，还有大量的孤块石，架空现象严重。虽然受红石岩堰塞湖地质、地形的限制，不能按照常规方案来彻底消除隐患，但是堰塞湖已具备综合利用的条件。应急救援阶段通过开挖应急泄流槽，联动上下游水库电站进行腾库的工程措施，有效地缓解了堰塞湖持续上涨的水位带来的溃坝风险；后续处置阶段通过利用原红石岩电站保留的引水隧洞进行堰塞湖库水的排泄，也为后续的永久性整治创造了干预条件。在此背景下，对红石岩堰塞湖治理工程危岩体特征与防治措施的研究显得尤为重要。通过深入研究危岩体的特征，包括其分布规律、成因机制、失稳模式等，可以为后续的防治措施提供科学依据。而有效的防治措施不仅能够保障当地居民的生命财产安全，避免因危岩体崩塌引发的二次灾害，还能为堰塞湖的综合整治和后续开发利用创造条件，使得堰塞湖从一个潜在的灾害源转变为能够造福当地的水利资源，实现“除害兴利、变废为宝”的目标，对促进当地的可持续发展具有深远意义。1.2国内外研究现状1.2.1危岩体成因机制研究现状在危岩体成因机制研究方面，国内外学者已取得了一系列成果。国外学者早期多从地质构造运动角度出发，如Lajtai通过对岩石节理裂隙发育与地质构造关系的研究，指出构造应力的长期作用导致岩体内部结构破坏，是危岩体形成的重要基础。随着研究的深入，越来越多的因素被纳入考虑范围。例如，一些学者开始关注风化作用对危岩体形成的影响，发现长期的风化作用可降低岩体强度，促进危岩体的发育。国内对危岩体成因机制的研究起步相对较晚，但发展迅速。众多学者结合我国复杂的地质环境，从多方面开展研究。在西南地区，由于地形地貌复杂，高山峡谷众多，学者们发现地形地貌对危岩体的形成具有显著控制作用。在高陡边坡地区，岩体卸荷作用强烈，使得岩体完整性遭到破坏，进而形成危岩体。构造运动在危岩体形成过程中也扮演着重要角色，如在板块碰撞带附近，强烈的构造应力使岩体产生大量裂隙，为危岩体的发育创造了条件。地震作为一种突发的强烈动力作用，会瞬间破坏岩体结构，促使大量危岩体形成，2008年汶川地震后，震区周边出现了众多危岩体，便是典型例证。此外，水的作用不可忽视，降雨入渗会增加岩体重量，降低岩体抗剪强度，地下水的长期浸泡和侵蚀也会弱化岩体结构。尽管在危岩体成因机制研究上已取得一定进展，但仍存在不足。不同因素之间的耦合作用机制研究尚不够深入，尤其是在多种复杂因素共同作用的情况下，如何准确量化各因素对危岩体形成的贡献，还缺乏系统的研究方法。1.2.2危岩体失稳模式研究现状对于危岩体失稳模式，国内外学者进行了大量研究，总结出坠落式、滑移式、倾倒式、滚（滑）落式等常见模式。在坠落式失稳模式研究中，国外学者通过实验和数值模拟，分析了危岩体在自重和外部荷载作用下从母体脱离坠落的过程，建立了相应的力学模型，用于预测坠落的可能性和轨迹。国内学者则结合实际工程案例，对坠落式危岩体的稳定性进行分析，提出了基于结构面特征和岩体力学参数的稳定性评价方法。在滑移式失稳模式方面，国外学者从岩体结构面的力学特性出发，研究了危岩体沿结构面滑动的力学机制，提出了基于极限平衡理论的滑移稳定性分析方法。国内学者在此基础上，考虑了地震、降雨等动态荷载对滑移稳定性的影响，建立了更加符合实际情况的稳定性分析模型。对于倾倒式危岩体，学者们主要研究了其在倾覆力矩作用下的变形破坏过程，通过理论分析和数值模拟，探讨了倾倒的临界条件和影响因素。滚（滑）落式失稳模式研究相对较少，主要集中在危岩体滚（滑）落过程中的运动学特征和能量变化方面。然而，现有研究对于复合失稳模式的关注较少，在实际工程中，危岩体可能同时受到多种因素影响，出现多种失稳模式相互耦合的情况，这方面的研究还存在较大空白。此外，在不同地质条件和环境因素下，危岩体失稳模式的演变规律也有待进一步深入研究。1.2.3危岩体稳定性评价方法研究现状危岩体稳定性评价方法可分为地质分析法、定量评价方法和数值模拟方法。地质分析法主要通过现场地质调查，对危岩体的地质条件、结构特征、变形迹象等进行分析，从而定性评价其稳定性。这种方法直观、简便，但主观性较强，评价结果的准确性依赖于评价人员的经验和专业水平。定量评价方法中，刚体极限平衡法应用最为广泛，它基于力的平衡原理，通过计算危岩体在各种荷载组合下的稳定性系数来评价其稳定性。该方法计算简单，但未考虑岩体的变形特性，在一些复杂情况下计算结果与实际情况存在偏差。有限元法、离散元法等数值模拟方法的出现，为危岩体稳定性评价提供了更有效的手段。有限元法能够考虑岩体的非线性力学特性和复杂的边界条件，离散元法则适用于模拟非连续介质的大变形和破坏过程。通过数值模拟，可以更直观地了解危岩体在不同工况下的应力应变分布和变形破坏过程，为稳定性评价提供更全面的信息。虽然危岩体稳定性评价方法不断发展，但仍存在一些问题。各种评价方法都有其局限性，如何合理选择和综合运用多种方法，提高评价结果的准确性和可靠性，是当前研究的重点之一。此外，对于一些特殊地质条件下的危岩体，如强风化岩体、岩溶地区危岩体等，现有的评价方法还不能很好地适应，需要进一步研究和改进。1.2.4危岩体防治措施研究现状针对危岩体的防治，国内外采取了多种措施，主要包括主动防治和被动防治。主动防治措施旨在增强危岩体的稳定性，如锚固、支撑、灌浆等。锚固技术通过锚杆或锚索将危岩体与稳定岩体连接在一起，提高危岩体的抗滑和抗倾覆能力。支撑则是在危岩体下方设置支撑结构，承受危岩体的部分重量。灌浆可以填充岩体裂隙，提高岩体的整体性和强度。国外在锚固和支撑技术方面发展较为成熟，研发了多种高性能的锚固材料和先进的支撑结构。国内在借鉴国外经验的基础上，结合自身工程实践，对锚固和支撑技术进行了改进和创新，使其更适用于我国的地质条件。被动防治措施主要是在危岩体可能崩塌的路径上设置拦截结构，如拦石墙、拦石网等，以阻挡崩塌的危岩体。在这方面，国内外都有大量的工程应用案例，并且不断改进拦截结构的设计和材料性能，提高其拦截效果。此外，还有一些生态防治措施，如植被护坡等，通过植物根系的加固作用和对坡面径流的调节作用，增强边坡的稳定性，减少危岩体的形成和崩塌。然而，在危岩体防治措施研究中，仍存在一些不足之处。不同防治措施的适用性和有效性缺乏系统的对比研究，在实际工程中难以根据具体情况选择最优的防治方案。对于一些复杂地质条件下的危岩体，单一防治措施往往难以奏效，需要综合运用多种措施，但目前对于综合防治措施的协同作用机制研究还不够深入。1.3研究内容与方法1.3.1研究内容本研究以鲁甸“8.3”地震红石岩堰塞湖治理工程中的危岩体为研究对象，围绕其特征与防治措施展开多方面研究。具体内容如下：研究区地质环境条件分析：收集并分析研究区的自然地理概况、地形地貌、地层岩性、地质构造、水文地质条件以及地震及地震动参数等基础资料，全面掌握危岩体发育的地质环境背景。在自然地理方面，了解研究区的气候、降水、水系等特征，因为降水作为水作用的主要来源，对危岩体稳定性有重要影响。地形地貌上，关注区域内的山脉走向、河谷形态、坡度与高差等，高陡地形易形成危岩体。地层岩性研究包括岩石类型、岩石强度、岩体结构等，不同岩性和结构的岩体抗风化和抗变形能力不同。地质构造方面，分析褶皱、断层的分布与特征，它们控制着岩体的完整性和裂隙发育程度。水文地质条件重点研究地下水水位、水量、径流方向以及与地表水的水力联系，地下水的渗流和静水压力会改变危岩体的力学平衡。地震及地震动参数则为评估地震对危岩体稳定性的影响提供依据。危岩体发育分布规律与成因机制研究：采用常规地质调查、无人机低空摄影测量技术和三维激光扫描技术等多种手段，对危岩体进行详细调查。通过常规地质调查，直接观察危岩体的出露情况、规模大小、形态特征、结构面产状等，并记录相关数据。无人机低空摄影测量技术可快速获取研究区的高分辨率影像，通过影像解译，绘制危岩体分布图，分析其总体分布规律。三维激光扫描技术能够精确获取危岩体的三维空间信息，包括表面形态、体积等，为后续的稳定性分析和防治措施设计提供准确的数据支持。在此基础上，对研究区危岩体进行分类，分析不同类型危岩体在坝址区右岸、左岸以及珍珠岩边坡等区域的分布规律，并深入探讨危岩体发育的影响因素，包括地形地貌、地层岩性组合、地质构造、地震、水以及其他因素等，进而揭示危岩体的成因机制及失稳模式。危岩体稳定性评价：系统阐述危岩体稳定性评价的方法，包括地质分析法和定量评价方法。地质分析法通过对危岩体的地质条件、结构特征、变形迹象等进行现场观察和分析，定性判断其稳定性。定量评价方法则采用刚体极限平衡法等，计算危岩体在不同荷载组合下的稳定性系数。选取典型危岩体，如坠落式危岩体（以右岸WYT19为例）和滑移式危岩体（以右岸WYT13为例），进行稳定性评价，进而对研究区危岩体的整体稳定性做出评估。在评价过程中，考虑不同工况，如天然状态、地震状态等，分析危岩体在不同情况下的稳定性变化。典型危岩体失稳模式、运动特征数值模拟分析与危害性分析：运用非连续变形分析（DDA）法等数值模拟方法，对典型危岩体的失稳模式与威胁范围进行数值模拟分析。以WYT11危岩体为例，首先明确其地质特征，然后建立数值模型，选取合适的计算参数，模拟其在天然状态和地震状态下的稳定性。通过模拟结果，分析危岩体的失稳模式、运动轨迹和速度等运动特征，评估其对周边环境的危害性。此外，对溢洪洞进口边坡“10.18”崩塌的成因机理及失稳模式进行深入分析，对ZZWY1危岩体变形破坏特征进行数值模拟分析，并计算其失稳后可能引发的涌浪，评估涌浪的危害性。危岩体防治措施研究：根据危岩体的稳定性评价结果和失稳模式分析，结合研究区的实际情况，制定危岩体防治措施原则。针对右岸边坡、左岸边坡和珍珠岩边坡等不同区域的危岩体，分别提出具体的防治措施。对于右岸边坡危岩体，根据不同区域的震裂损伤程度、卸荷特征和危岩体稳定性分类，采用放坡开挖、锚喷支护、小扰动清危、主被动拦挡等不同措施进行治理。左岸边坡顶部裂缝发育部位采用放坡开挖+锚喷支护的措施处理，中下部考虑到危害对象为公路和下部堰塞坝，建议公路采用开凿隧道的方式通过，并在滑坡堆积体缓坡上设置落石槽和拦石墙，拦截可能威胁堰塞坝的危岩体。珍珠岩边坡危岩体也根据其具体特征，选择合适的防治措施，如锚固、支撑、灌浆等主动防治措施，以及拦石网、拦石墙等被动防治措施，确保危岩体得到有效治理，保障堰塞湖治理工程的安全和周边居民的生命财产安全。1.3.2研究方法文献资料分析法：广泛查阅国内外有关危岩体成因机制、失稳模式、稳定性评价方法和防治措施等方面的文献资料，了解该领域的研究现状和发展趋势，为本研究提供理论基础和研究思路。通过对文献的梳理，总结不同研究方法的优缺点，分析现有研究中存在的不足，从而确定本研究的重点和方向。例如，在危岩体稳定性评价方法的研究中，对比分析地质分析法、定量评价方法和数值模拟方法的适用条件和局限性，为后续研究中方法的选择提供参考。现场调查法：对红石岩堰塞湖治理工程区域进行详细的现场调查，包括危岩体的位置、规模、形态、结构面特征、风化程度等。在调查过程中，采用地质罗盘测量结构面产状，用皮尺、全站仪等测量危岩体的几何尺寸，观察危岩体的变形迹象，如裂缝宽度、深度和延伸方向等。同时，对研究区的地质环境条件进行实地考察，验证和补充收集到的资料。通过现场调查，获取第一手资料，为危岩体特征分析和防治措施研究提供真实可靠的数据。无人机低空摄影测量技术：利用无人机搭载高清相机，对研究区进行低空摄影测量。通过获取的高分辨率影像，进行影像拼接和三维建模，生成研究区的数字正射影像图（DOM）和数字高程模型（DEM）。利用这些成果，可直观地观察危岩体的分布范围和形态特征，绘制危岩体分布图，分析其分布规律。无人机低空摄影测量技术具有快速、高效、成本低等优点，能够获取大面积的地形和危岩体信息，弥补了传统地面调查的局限性。三维激光扫描技术：运用三维激光扫描仪对危岩体进行扫描，获取危岩体表面的三维坐标信息，生成高精度的三维点云模型。通过对三维点云模型的处理和分析，可以精确测量危岩体的体积、表面积、裂缝宽度等参数，了解危岩体的空间形态和结构特征。三维激光扫描技术能够实现非接触式测量，获取的数据精度高、信息全面，为危岩体的稳定性分析和防治措施设计提供准确的数据支持。数值模拟方法：采用非连续变形分析（DDA）法、有限元法（FEM）、离散元法（DEM）等数值模拟软件，对危岩体的稳定性、失稳模式和运动特征进行模拟分析。根据危岩体的地质条件和力学参数，建立合理的数值模型，模拟不同工况下危岩体的受力和变形情况。通过数值模拟，可以直观地展示危岩体在各种因素作用下的变化过程，预测危岩体的失稳可能性和破坏范围，为危岩体的防治提供科学依据。例如，利用DDA法模拟危岩体在地震作用下的失稳过程，分析其运动轨迹和速度，评估其对周边环境的影响。二、红石岩堰塞湖形成背景与地质环境2.1鲁甸“8.3”地震概况2014年8月3日16时30分，一场6.5级的强烈地震在云南省昭通市鲁甸县骤然爆发，震源深度12km，震中精准地落在昭通市鲁甸县龙头山镇、火德红镇和巧家县包谷垴乡附近。此次地震震级颇高，震源又极浅，这两大因素叠加，使得其破坏力惊人，地震最高烈度达到了Ⅸ度。此次地震的发生，在鲁甸县及周边地区引发了一系列极为严重的后果。截至2014年8月8日15时，地震无情地夺走了617人的生命，其中鲁甸县526人、巧家县78人、昭阳区1人、会泽县12人；112人不幸失踪，3143人受伤，受灾范围更是广泛，波及108.84万人，8.09万间房屋轰然倒塌。如此惨重的人员伤亡和巨大的财产损失，给当地民众带来了沉重的灾难，也让整个社会为之痛心。从地震的成因来看，它发生在北东向西鱼河－昭通断裂附近，震源机制结果清晰地显示为走滑型破裂。这种破裂方式使得地壳岩石在水平方向上发生相对错动，导致岩体结构遭受严重破坏。地震产生的强烈地震波在传播过程中，引发地面的剧烈震动，这种震动使得原本相对稳定的山体岩石内部产生了巨大的应力变化。在高陡地形区域，山体的稳定性本身就较为脆弱，地震产生的应力远远超出了岩石的承受极限，使得岩石内部产生大量裂隙，岩体的完整性被彻底打破，为后续的山体崩塌、滑坡等地质灾害埋下了祸根。鲁甸“8.3”地震对红石岩地区的地质结构造成了毁灭性的影响。在红石岩地区，地震引发了山体大规模的崩塌和滑坡。原本完整的山体在地震的强烈作用下，岩石破碎，土体松动，大量的岩土体沿着山坡滑落，堵塞了牛栏江，最终形成了红石岩堰塞湖。地震还导致该地区岩体的结构面进一步发育和扩展，结构面的抗剪强度大幅降低，使得岩体的稳定性急剧下降，为危岩体的形成创造了极为有利的条件。这些危岩体在后续的风化、降雨、地震等因素作用下，随时可能发生崩塌，对堰塞湖治理工程以及周边居民的生命财产安全构成了严重威胁。2.2红石岩堰塞湖形成过程与特征2014年8月3日，鲁甸“8.3”地震的巨大破坏力使得牛栏江两岸山体发生了大规模的垮塌。地震产生的强烈地震波在传播过程中，引发地面的剧烈震动，这种震动使得原本相对稳定的山体岩石内部产生了巨大的应力变化。在高陡地形区域，山体的稳定性本身就较为脆弱，地震产生的应力远远超出了岩石的承受极限，使得岩石内部产生大量裂隙，岩体的完整性被彻底打破。牛栏江两岸的山体在地震的强烈作用下，岩土体沿着山坡滑落，大量的岩土体在牛栏江干流红石岩水电站取水坝下游约600m处迅速堆积，短时间内便堵塞了牛栏江河道，从而形成了红石岩堰塞湖。从规模上看，红石岩堰塞湖堰塞体规模巨大，最大高度达到102m，这一高度使得堰塞体能够阻挡大量的河水，为堰塞湖的形成和蓄水提供了基础条件。顺河向底部长约910m，顶部长约262m，垂直河流方向平均顶宽301m，总体呈现出不对称的马鞍状，这种独特的形状也影响着堰塞湖的蓄水量和水流形态。堰塞体总方量约1200万m³，如此庞大的方量表明堰塞体具有较强的稳定性，但也增加了后续治理和处置的难度。堆积物系右岸基岩（白云岩和灰岩）崩滑形成，碎块石含量占85%以上，粒径大部分在0.5-5.0m之间，这种物质组成使得堰塞体的结构较为复杂，架空现象严重，进一步影响了堰塞体的稳定性和渗透性。在水位变化方面，堰塞体形成后，湖水位在17h内快速上涨24.50m，涨幅达到1.42m/h，这样的上涨速度极为惊人，反映出堰塞湖在形成初期水位变化的急剧性。红石岩水电站取水坝也因此完全被淹没，这不仅对水电站的设施造成了严重破坏，也改变了原有的水利格局。照此上涨速度推算，2d内湖水可能漫顶，一旦漫顶，将形成集水面积11832km²、库容2.6亿m³、回水长度25km的大型堰塞湖。在后续的发展过程中，通过对堰塞湖水位监测资料的分析，发现其水位变化呈现出涨落急剧的特征，水位变幅达45.06m，最低水位为1137.5m，最高水位达到1182.56m。堰塞湖持续高水位时间长，且峰谷交错，呈现出多峰型的水位变化特征。水位变化主要受上游暴雨洪水过程、泄流量的大小以及上游德泽水库的调蓄等多种因素的综合影响。在暴雨洪水季节，上游来水量大幅增加，导致堰塞湖水位迅速上升；而当泄流量增大或者德泽水库进行有效调蓄时，堰塞湖水位则会相应下降。2.3区域地质环境条件2.3.1地形地貌红石岩堰塞湖位于滇东北高原向四川盆地过渡的斜坡地带，区域内地形起伏剧烈，山高谷深，属于典型的高山峡谷地貌。研究区地势总体上呈现出西南高、东北低的态势，山脉走向多为北东-南西向，与区域地质构造方向基本一致。牛栏江蜿蜒贯穿其中，河流下切作用强烈，河谷呈“V”字形，谷坡陡峭，坡度多在35°-70°之间。在堰塞湖周边，地形地貌特征复杂多样。右岸山体高大陡峭，悬崖林立，岩体裸露，部分区域形成了负角度的陡峭岩壁，如地震后形成的600米高的负角度岩壁，这种特殊的地形地貌使得岩体卸荷作用明显，岩体完整性遭到严重破坏，为危岩体的形成提供了有利的地形条件。左岸地形相对较为平缓，但也存在一定规模的滑坡和崩塌体，这些不良地质现象改变了原有的地形地貌，增加了岩体的不稳定性。在河流阶地方面，研究区内可见多级河流阶地，阶地堆积物主要为第四系松散堆积物，包括砂、砾石、粉质黏土等。这些阶地堆积物在地震等外力作用下，容易发生变形和滑动，对周边的危岩体稳定性也会产生一定的影响。2.3.2地层岩性研究区出露的地层主要有震旦系、寒武系、奥陶系、二叠系和第四系。震旦系主要为一套浅变质的碎屑岩和火山岩，岩石致密坚硬，但经过长期的地质构造运动和风化作用，岩体内部发育有大量的节理裂隙，降低了岩体的完整性和强度。寒武系地层以碳酸盐岩和碎屑岩为主，其中碳酸盐岩岩溶发育，岩体溶蚀现象明显，进一步破坏了岩体的结构。奥陶系主要由页岩、砂岩和灰岩组成，页岩的抗风化能力较弱，在风化作用下容易破碎，砂岩和灰岩则相对较硬，但也存在不同程度的风化和裂隙发育。二叠系地层主要为峨眉山玄武岩，岩石致密，强度较高，但在构造应力作用下，岩体也产生了大量的破裂面。第四系地层分布广泛，主要包括残积层、坡积层、冲积层和崩滑堆积层等。残积层主要分布在山顶和山坡上部，由岩石风化残积而成，物质组成以粉质黏土和碎石为主，结构松散。坡积层多分布在山坡中下部，是由山坡上的岩石风化碎屑物在重力和坡面水流作用下堆积而成，其厚度和物质组成变化较大，稳定性较差。冲积层主要分布在河流两岸和河谷底部，由河流搬运和沉积作用形成，物质主要为砂、砾石和粉质黏土，分选性和磨圆度较好。崩滑堆积层则主要分布在滑坡和崩塌体区域，如红石岩堰塞湖堰塞体就是由地震引发的山体崩滑堆积形成，主要由白云岩和灰岩的碎块石组成，碎块石含量占85%以上，粒径大部分在0.5-5.0m之间，堆积体结构复杂，架空现象严重。2.3.3地质构造研究区位于扬子准地台西缘，处于小江断裂带和昭通-大关断裂带之间，地质构造复杂，新构造运动活跃。区域内主要发育有褶皱和断裂构造，褶皱形态以紧闭褶皱和倒转褶皱为主，轴向多为北东-南西向，褶皱的存在使得地层发生弯曲和变形，岩体内部产生了复杂的应力状态，增加了岩体的破碎程度。断裂构造对研究区的地质环境和危岩体发育影响显著。主要断裂有北东向的小江断裂和北西向的昭通-大关断裂，这些断裂在漫长的地质历史时期中多次活动，导致岩体破碎，节理裂隙密集发育，为危岩体的形成提供了良好的边界条件。在红石岩堰塞湖周边，受到断裂构造的影响，岩体完整性遭到严重破坏，形成了大量的构造破碎带，这些破碎带中的岩体抗风化能力弱，容易发生崩塌和滑坡。此外，断裂构造还控制着地下水的运移和富集，地下水在构造破碎带中流动，对岩体进行溶蚀和软化作用，进一步降低了岩体的稳定性。2.3.4水文地质条件研究区的水文地质条件较为复杂，地下水类型主要有松散岩类孔隙水、基岩裂隙水和岩溶水。松散岩类孔隙水主要赋存于第四系松散堆积层中，其富水性受堆积层的厚度、颗粒组成和地形地貌等因素控制。在河流阶地和沟谷地带，松散堆积层较厚，颗粒较粗，孔隙水相对较为丰富；而在山坡上部，堆积层较薄，颗粒较细，富水性较差。基岩裂隙水主要赋存于基岩的节理裂隙中，其富水性与岩体的裂隙发育程度、连通性以及地形地貌等密切相关。在构造破碎带和风化强烈的区域，基岩裂隙发育，连通性好，基岩裂隙水相对较丰富。岩溶水主要分布在寒武系和奥陶系的碳酸盐岩地层中，由于岩溶发育，岩溶水的分布极不均匀，存在着溶洞、溶蚀裂隙等岩溶通道，岩溶水在其中流动迅速，水位变化较大。在红石岩堰塞湖区域，地下水的补给主要来源于大气降水和地表水的入渗，大气降水通过地表径流和下渗的方式进入地下，补充地下水。地表水与地下水之间存在着密切的水力联系，在河流两岸和低洼地带，地表水与地下水相互补给。地下水的排泄方式主要有泉排泄、向河流排泄以及人工开采等。在红石岩堰塞湖形成后，湖水水位的变化对周边的水文地质条件产生了显著影响。湖水的淹没使得周边地下水位上升，改变了地下水的渗流场和水力梯度，导致岩体的饱水程度增加，有效应力降低，抗剪强度减小，从而影响危岩体的稳定性。此外，湖水的长期浸泡还会对岩体产生溶蚀和软化作用，进一步削弱岩体的强度。三、危岩体特征分析3.1危岩体调查方法与数据获取在对红石岩堰塞湖治理工程危岩体进行研究时，采用了多种先进且互补的调查方法，以全面、准确地获取危岩体的相关数据。这些方法的综合运用，为深入了解危岩体的特征奠定了坚实基础。常规地质调查是危岩体研究的基础手段，通过实地勘查，直接观察危岩体的出露情况、规模大小、形态特征、结构面产状等。在红石岩堰塞湖区域，地质调查人员利用地质罗盘对危岩体的结构面产状进行测量，获取其走向、倾向和倾角等关键信息，这些数据对于分析危岩体的稳定性和失稳模式至关重要。使用皮尺、全站仪等工具精确测量危岩体的几何尺寸，包括长度、宽度、高度等，详细记录危岩体的位置、周边地形地貌特征以及与其他地质构造的关系。通过仔细观察危岩体的表面特征，如风化程度、裂隙发育情况、是否存在变形迹象等，判断危岩体的稳定性状态。常规地质调查能够获取第一手资料，直观地了解危岩体的基本情况，但在面对高陡边坡、难以到达的区域时，存在一定的局限性。无人机低空摄影测量技术作为一种高效的调查手段，在危岩体调查中发挥了重要作用。该技术利用无人机搭载高清相机，对研究区进行低空飞行拍摄，获取高分辨率影像。在红石岩堰塞湖治理工程危岩体调查中，无人机能够快速覆盖大面积区域，获取研究区的整体影像资料。通过对影像的拼接和处理，生成研究区的数字正射影像图（DOM）和数字高程模型（DEM）。在DOM上，可以清晰地观察危岩体的分布范围、形态特征以及与周边地形的关系，准确绘制危岩体分布图，分析其总体分布规律。DEM则提供了危岩体区域的地形起伏信息，有助于进一步了解危岩体所处的地形条件对其稳定性的影响。无人机低空摄影测量技术具有快速、高效、成本低、机动性强等优点，能够弥补常规地质调查在大面积区域调查和高陡地形区域调查的不足，为危岩体调查提供了宏观的视角和全面的影像资料。三维激光扫描技术则为危岩体调查提供了高精度的三维空间信息。运用三维激光扫描仪对危岩体进行扫描，仪器发射激光束并接收反射信号，快速获取危岩体表面的三维坐标信息，生成高精度的三维点云模型。在对红石岩堰塞湖危岩体的扫描过程中，通过对三维点云模型的处理和分析，可以精确测量危岩体的体积、表面积、裂缝宽度、深度等参数，了解危岩体的空间形态和结构特征。还能通过点云数据提取危岩体的结构面信息，包括结构面的产状、间距、连通性等，为危岩体的稳定性分析提供准确的数据支持。三维激光扫描技术实现了非接触式测量，获取的数据精度高、信息全面，能够详细呈现危岩体的细微特征和复杂结构，解决了常规测量方法在危岩体复杂结构测量中的难题。在实际调查过程中，将这三种方法有机结合，相互补充。先利用无人机低空摄影测量技术进行大面积的初步调查，获取研究区危岩体的整体分布情况和宏观特征，为后续的详细调查提供方向和范围。然后，通过常规地质调查对重点危岩体区域进行实地勘查，进一步了解危岩体的具体特征和地质条件。最后，运用三维激光扫描技术对关键危岩体进行高精度测量，获取详细的三维空间信息，为危岩体的稳定性评价和防治措施设计提供精确的数据基础。通过多种调查方法的综合运用，全面、准确地获取了红石岩堰塞湖治理工程危岩体的相关数据，为后续的研究工作提供了有力保障。3.2危岩体发育分布规律通过对红石岩堰塞湖治理工程区域的详细调查和数据分析，发现该区域危岩体可依据其形成机制和形态特征分为坠落式危岩体、滑移式危岩体和倾倒式危岩体三大类。坠落式危岩体多由岩体受多组结构面切割，在重力作用下与母岩分离形成，常呈孤立块状，底部与母岩接触面积小，稳定性差。滑移式危岩体则是在重力和其他外力作用下，沿着特定的结构面发生滑动，其滑动面多为节理、裂隙或软弱夹层。倾倒式危岩体通常是由于岩体底部支撑不足，在重力和侧向力作用下发生倾倒，这类危岩体一般具有较高的高度和较小的底部支撑面积。从总体分布来看，研究区危岩体数量众多，共计335处，其中坝址区右岸分布最为密集，有239处，左岸有70处，珍珠岩边坡有26处。这种分布差异与各区域的地形地貌、地层岩性、地质构造等因素密切相关。坝址区右岸地形陡峭，山体高大，岩体在地震、风化等作用下更容易破碎和分离，形成大量危岩体。左岸地形相对较为平缓，危岩体数量相对较少，但仍存在一些因岩体结构破碎和风化作用而形成的危岩体。珍珠岩边坡由于其特殊的地层岩性和地质构造，危岩体数量相对较少，但也不容忽视。在坝址区右岸，危岩体主要集中在高程1200-1500m的区域，该区域地形陡峭，坡度多在50°-70°之间，岩体受地震影响震裂损伤严重，节理裂隙发育，为危岩体的形成提供了有利条件。根据震裂损伤程度和卸荷特征，右岸可进一步分为4个区域。B1-1区高程在1200-1300m之间，地震震裂损伤强烈，卸荷松弛明显，岩体完整性差，危岩体数量众多，且规模较大，以坠落式和滑移式危岩体为主。在该区域的一些悬崖峭壁处，可见大量块状危岩体悬挂在岩体表面，底部仅与母岩有少量接触，随时可能坠落。B1-2区高程在1300-1400m之间，震裂损伤和卸荷特征相对B1-1区较弱，但仍有较多危岩体分布，以坠落式危岩体居多。B1-3区高程在1400-1500m之间，震裂损伤和卸荷特征进一步减弱，但由于岩体风化作用较强，仍有部分危岩体存在，多为小型坠落式危岩体。B1-4区位于右岸顶部，高程大于1500m，虽然地震震裂损伤相对较弱，但由于长期的风化作用和岩体卸荷，也有一定数量的危岩体分布，以倾倒式危岩体较为常见。坝址区左岸危岩体主要分布在高程1150-1300m的区域，该区域地形相对右岸较为平缓，但仍存在一些高陡边坡。危岩体多分布在边坡的顶部和中部，以坠落式和滑移式危岩体为主。在左岸顶部裂缝发育部位，由于岩体受拉应力作用，形成了一些规模较大的滑移式危岩体。而在边坡中部，由于岩体风化和卸荷作用，部分岩体与母岩分离，形成坠落式危岩体。左岸危岩体的稳定性相对右岸较好，但仍对下方的公路和堰塞坝构成一定威胁。珍珠岩边坡危岩体主要分布在边坡的上部，高程在1350-1500m之间，以坠落式危岩体为主。该区域岩体岩性主要为白云岩和灰岩，由于岩石节理裂隙发育，在风化、降雨等作用下，岩体逐渐破碎，形成危岩体。部分危岩体底部与母岩之间的连接较为薄弱，在重力作用下容易发生坠落。珍珠岩边坡危岩体虽然数量相对较少，但由于其下方为重要的交通要道和居民区，一旦发生崩塌，可能会造成严重的人员伤亡和财产损失。3.3危岩体成因机制红石岩堰塞湖治理工程区域危岩体的形成是多种因素共同作用的结果，这些因素相互影响、相互制约，深刻改变了岩体的结构和稳定性，最终导致危岩体的发育。地形地貌是危岩体形成的重要基础条件。研究区处于高山峡谷地貌区，山高谷深，地形起伏剧烈，牛栏江两岸谷坡陡峭，坡度多在35°-70°之间。这种高陡的地形使得岩体在重力作用下，更容易产生变形和破坏。在陡峭的山坡上，岩体所受的重力沿坡面的分力增大，当超过岩体的抗滑力时，岩体就会发生滑动或崩塌，形成危岩体。右岸山体高大陡峭，悬崖林立，部分区域形成负角度陡峭岩壁，岩体卸荷作用明显，在地震等动力作用下，岩体完整性遭到严重破坏，大量岩体与母岩分离，形成坠落式和滑移式危岩体。高陡地形还使得雨水和地表径流对岩体的冲刷作用增强，进一步削弱了岩体的稳定性，促进了危岩体的形成。地层岩性组合对危岩体的发育有着重要影响。研究区出露的地层岩性多样，不同岩性的岩石力学性质差异较大。震旦系浅变质碎屑岩和火山岩、寒武系碳酸盐岩和碎屑岩、奥陶系页岩、砂岩和灰岩、二叠系峨眉山玄武岩等，在长期的地质作用下，形成了复杂的地层岩性组合。当软硬岩性互层时，软弱岩层在风化、水的作用下容易被侵蚀和软化，导致上部坚硬岩层失去支撑，从而形成危岩体。在寒武系地层中，页岩等软弱岩层在风化作用下破碎，上部的灰岩等坚硬岩层因失去支撑而悬空，形成危岩体。岩石的结构和完整性也影响危岩体的形成，节理裂隙发育的岩体，其完整性较差，抗风化和抗变形能力弱，容易在各种外力作用下形成危岩体。地质构造是危岩体形成的关键因素之一。研究区位于扬子准地台西缘，处于小江断裂带和昭通-大关断裂带之间，地质构造复杂，新构造运动活跃。褶皱和断裂构造使得岩体产生变形和破裂，形成大量的节理裂隙和破碎带。褶皱作用使地层发生弯曲和变形，在褶皱的轴部和翼部，岩体内部应力集中，容易产生裂隙和破碎。断裂构造的多次活动导致岩体破碎，节理裂隙密集发育，为危岩体的形成提供了良好的边界条件。在红石岩堰塞湖周边，受到断裂构造的影响，岩体完整性遭到严重破坏，形成了大量的构造破碎带，这些破碎带中的岩体抗风化能力弱，在风化、地震等作用下，极易发生崩塌和滑坡，形成危岩体。断裂构造还控制着地下水的运移和富集，地下水在构造破碎带中流动，对岩体进行溶蚀和软化作用，进一步降低了岩体的稳定性，促进危岩体的发育。地震是危岩体形成的重要触发因素。鲁甸“8.3”地震震级高、震源浅，产生的强烈地震波对研究区岩体造成了巨大的破坏。地震使得岩体内部产生大量的裂隙，岩体结构遭到严重破坏，完整性丧失，强度大幅降低。在地震作用下，岩体所受的地震力与重力、地应力等共同作用，当超过岩体的承载能力时，岩体就会发生崩塌、滑坡等地质灾害，形成危岩体。红石岩堰塞湖的形成就是地震引发山体崩塌和滑坡的结果，而在堰塞湖周边，地震造成的岩体破坏也为危岩体的发育创造了条件。地震还可能导致岩体内部结构的重新调整，使得原本处于稳定状态的岩体变得不稳定，进一步增加了危岩体的数量和规模。水的作用在危岩体形成过程中不可忽视。研究区降水丰富，地表水和地下水活动频繁。降雨入渗会增加岩体的重量，使岩体所受的重力增大，同时降低岩体的抗剪强度，容易引发岩体的滑动和崩塌。雨水还会沿着岩体的节理裂隙流动，对岩体进行冲刷和侵蚀，扩大裂隙的宽度和深度，促进危岩体的形成。在红石岩堰塞湖形成后，湖水水位的变化对周边岩体产生了显著影响。湖水的淹没使得周边地下水位上升，岩体饱水程度增加，有效应力降低，抗剪强度减小，导致危岩体的稳定性下降。湖水的长期浸泡还会对岩体产生溶蚀和软化作用，进一步削弱岩体的强度，加速危岩体的形成和发展。此外，风化作用、人类工程活动等其他因素也对危岩体的发育有一定影响。长期的风化作用使得岩体表面的岩石逐渐破碎，形成松散的碎屑物，降低了岩体的强度和稳定性，为危岩体的形成创造了条件。人类工程活动如开挖、爆破等，改变了岩体的原始应力状态，破坏了岩体的完整性，也可能诱发危岩体的形成。3.4危岩体失稳模式坠落式失稳模式在红石岩堰塞湖治理工程危岩体中较为常见。这类危岩体通常被多组结构面切割，底部与母岩的连接较为薄弱，在重力作用下，危岩体从母岩上脱离并坠落。在坝址区右岸的一些高陡边坡区域，部分危岩体呈孤立块状，底部仅与母岩有少量接触，随着风化、地震等作用的持续影响，底部连接逐渐破坏，最终导致危岩体坠落。从力学原理来看，当危岩体所受的重力对底部支撑点产生的力矩超过底部连接的抗倾覆力矩时，危岩体就会发生坠落。坠落式危岩体的运动轨迹通常较为简单，基本沿垂直方向下落，但在下落过程中可能会与坡面发生碰撞，改变运动方向，对下方的建筑物、道路、人员等造成严重威胁。滑移式失稳模式也是研究区危岩体常见的失稳方式之一。这类危岩体在重力和其他外力作用下，沿着特定的结构面发生滑动。在坝址区左岸，部分危岩体由于岩体内部存在节理、裂隙或软弱夹层等结构面，这些结构面的抗剪强度较低。当危岩体所受的下滑力超过结构面的抗剪强度时，危岩体就会沿着结构面发生滑动。地震、降雨等因素会增加危岩体的下滑力，同时降低结构面的抗剪强度，从而促进危岩体的滑动。在地震作用下，危岩体受到地震力的作用，下滑力增大；降雨入渗会使危岩体重量增加，同时软化结构面，降低其抗剪强度。滑移式危岩体的运动速度相对较慢，但由于其滑动距离可能较长，对下方的基础设施和居民安全同样构成较大威胁。倾倒式失稳模式在研究区的部分危岩体中也有出现。这类危岩体通常具有较高的高度和较小的底部支撑面积，在重力和侧向力作用下，危岩体底部的支撑力不足以维持其平衡，从而发生倾倒。在珍珠岩边坡的上部，一些危岩体由于长期的风化作用和岩体卸荷，底部支撑逐渐减弱，在重力和风力等侧向力的作用下，容易发生倾倒。从力学分析角度，当危岩体所受的重力和侧向力对底部支撑点产生的合力矩超过底部支撑的抗倾覆力矩时，危岩体就会发生倾倒。倾倒式危岩体在倾倒过程中，会对周围的岩体产生挤压和破坏作用，引发连锁反应，可能导致更大范围的岩体失稳。在红石岩堰塞湖治理工程区域，由于地质条件复杂，多种因素相互作用，危岩体的失稳模式可能并非单一存在，而是多种模式相互耦合。一些危岩体在坠落过程中，可能会与坡面发生碰撞，改变运动方向，产生滚动或滑动，形成坠落-滚（滑）落复合失稳模式；部分滑移式危岩体在滑动过程中，由于地形变化或遇到障碍物，可能会发生倾倒，形成滑移-倾倒复合失稳模式。这种复合失稳模式增加了危岩体失稳的复杂性和危害性，使得对危岩体的稳定性评价和防治措施设计更加困难。四、危岩体稳定性评价4.1稳定性评价方法概述在危岩体稳定性评价领域，地质分析法和定量评价方法是两种重要且应用广泛的手段，它们从不同角度对危岩体的稳定性进行评估，为后续的防治工作提供了关键依据。地质分析法作为一种定性的评价方法，具有直观、全面的特点。它主要基于现场地质调查，通过对危岩体所处的地质条件进行深入剖析，包括地层岩性、地质构造、地形地貌等方面，来综合判断危岩体的稳定性。在红石岩堰塞湖治理工程中，地质分析法发挥了重要作用。通过对研究区地层岩性的分析，了解到出露地层的岩石类型、岩石强度以及岩体结构等信息，为判断危岩体的稳定性提供了基础资料。震旦系浅变质碎屑岩和火山岩，虽岩石致密坚硬，但经过长期地质构造运动和风化作用，岩体内部节理裂隙发育，降低了岩体的完整性和强度，使得该地层中的危岩体稳定性相对较差。地质构造对危岩体稳定性的影响也不容忽视。在研究区，处于小江断裂带和昭通-大关断裂带之间，地质构造复杂，新构造运动活跃。褶皱和断裂构造使得岩体产生变形和破裂，形成大量的节理裂隙和破碎带，这些破碎带中的岩体抗风化能力弱，容易发生崩塌和滑坡，从而影响危岩体的稳定性。通过对地质构造的分析，能够明确危岩体与构造破碎带的关系，判断危岩体是否处于构造应力集中区域，进而评估其稳定性。在红石岩堰塞湖周边，受到断裂构造的影响，岩体完整性遭到严重破坏，形成了大量的构造破碎带，这些破碎带中的危岩体稳定性较差，容易发生失稳。地形地貌也是地质分析法中需要重点考虑的因素。红石岩堰塞湖位于高山峡谷地貌区，山高谷深，地形起伏剧烈，牛栏江两岸谷坡陡峭，坡度多在35°-70°之间。这种高陡的地形使得岩体在重力作用下，更容易产生变形和破坏，形成危岩体。通过对地形地貌的分析，能够确定危岩体所处的地形部位，判断其是否处于高陡边坡、悬崖峭壁等容易发生崩塌的区域，从而评估其稳定性。在右岸山体高大陡峭，悬崖林立，部分区域形成负角度陡峭岩壁，这些区域的危岩体稳定性较差，在地震、风化等作用下，容易发生坠落等失稳现象。定量评价方法则侧重于运用数学和力学原理，通过具体的计算来量化危岩体的稳定性。刚体极限平衡法是定量评价方法中应用最为广泛的一种。该方法基于力的平衡原理，假设危岩体为刚体，通过分析危岩体在自重、裂隙水压力、地震力等各种荷载组合作用下的受力情况，计算危岩体的稳定性系数，以此来评价其稳定性。在计算过程中，需要确定危岩体的几何形状、尺寸、结构面产状、岩体物理力学参数等。对于滑移式危岩体，要计算其沿滑动面的下滑力和抗滑力，根据两者的比值得到稳定性系数；对于倾倒式危岩体，则要考虑其倾覆力矩和抗倾覆力矩的关系。在红石岩堰塞湖治理工程危岩体稳定性评价中，刚体极限平衡法被广泛应用于典型危岩体的稳定性计算。通过现场调查和试验，获取危岩体的相关参数，如岩体的重度、结构面的粘聚力和内摩擦角等，然后代入刚体极限平衡法的计算公式中，计算不同工况下危岩体的稳定性系数。在计算坠落式危岩体稳定性时，要考虑危岩体的自重、地震力以及后缘裂隙的抗拉强度等因素，通过特定的公式计算稳定性系数。在天然状态下，某坠落式危岩体的稳定性系数为1.2，表明其处于基本稳定状态；但在地震工况下，稳定性系数降至0.9，说明此时危岩体处于不稳定状态，需要采取相应的防治措施。地质分析法和定量评价方法各有优缺点，在实际应用中，通常将两者结合起来，相互补充，以提高危岩体稳定性评价的准确性和可靠性。地质分析法能够提供危岩体的宏观地质背景和定性的稳定性判断，而定量评价方法则能通过具体的计算结果，为危岩体的稳定性提供量化的依据，两者的有机结合，有助于制定更加科学合理的危岩体防治措施。4.2典型危岩体稳定性计算以右岸WYT19坠落式危岩体为例，该危岩体位于坝址区右岸，高程约1350m，处于B1-2区域。它呈块状，底部与母岩的连接较为薄弱，后缘发育有陡倾裂隙。根据现场调查和测量，获取了危岩体的相关参数：危岩体体积约为50m³，重度取26kN/m³，因此危岩体自重W=50×26=1300kN。危岩体后缘裂隙深度h=3m，危岩体重心到潜在破坏面的水平距离a0=2m，重心到过潜在破坏面形心的铅垂距离b0=1.5m。岩体粘聚力标准值c=50kPa，内摩擦角标准值φ=30°，根据岩石抗拉强度标准值乘以0.2的折减系数确定危岩体抗拉强度标准值flk，假设岩石抗拉强度标准值为200kPa，则flk=200×0.2=40kPa。地震水平系数ζe取0.05，地震力Q=ζe×W=0.05×1300=65kN。对于后缘有陡倾裂隙的悬挑式坠落式危岩，按下列两个公式计算稳定性系数，取较小值：F_1=\frac{c（H-h）-Q\tan\varphi}{W}F_2=\frac{\zetaf_{lk}（H-h）}{W\timesa_0+Q\timesb_0}危岩体抗弯力矩计算系数ζ依据潜在破坏面形态取值，此处取1/12。代入数据计算：F_1=\frac{50×（5-3）-65×\tan30°}{1300}\approx0.03F_2=\frac{\frac{1}{12}×40×（5-3）}{1300×2+65×1.5}\approx0.002取较小值，F=0.002，远小于1.0，表明在当前工况下，右岸WYT19坠落式危岩体处于不稳定状态。再以右岸WYT13滑移式危岩体进行稳定性计算。该危岩体位于右岸高程1250m处，属于B1-1区域。它受两组结构面切割，形成了潜在的滑动面，滑动面倾向坡外，倾角α=35°。通过现场调查和试验，获取相关参数：危岩体体积为80m³，重度为26kN/m³，危岩体自重W=80×26=2080kN。后缘裂隙深度h=4m，裂隙充水高度hw取裂隙深度的1/3，即hw=4×1/3≈1.33m。水容重γw取10kN/m³，根据裂隙水压力公式V=1/2γwhw²，可得裂隙水压力V=1/2×10×1.33²≈8.85kN/m。地震水平系数ζe=0.05，地震力Q=ζe×W=0.05×2080=104kN。后缘裂隙粘聚力标准值c=30kPa，后缘裂隙内摩擦角标准值φ=25°。根据滑移式危岩体稳定性计算公式：F=\frac{(W\cos\alpha-Q\sin\alpha-V)\tan\varphi+c\timesl}{W\sin\alpha+Q\cos\alpha}假设结构面长度l=5m，代入数据计算：F=\frac{(2080×\cos35°-104×\sin35°-8.85)\tan25°+30×5}{2080×\sin35°+104×\cos35°}\approx0.78计算结果表明，右岸WYT13滑移式危岩体在当前工况下稳定性系数为0.78，小于1.0，处于不稳定状态。4.3研究区危岩体整体稳定性评价通过对研究区内典型危岩体稳定性的计算，能够进一步推断整个研究区危岩体的稳定性状况。在红石岩堰塞湖治理工程区域，危岩体数量众多，分布广泛，不同区域和类型的危岩体稳定性存在较大差异。根据对坝址区右岸、左岸以及珍珠岩边坡危岩体的调查和分析，结合典型危岩体稳定性计算结果，可将研究区危岩体整体稳定性划分为三个等级：不稳定、基本稳定和稳定。不稳定危岩体主要分布在坝址区右岸B1-1区，该区域地震震裂损伤强烈，卸荷松弛明显，岩体完整性差，如右岸WYT19坠落式危岩体和WYT13滑移式危岩体就处于该区域，计算结果表明它们在当前工况下稳定性系数远小于1.0，处于不稳定状态。这些危岩体在重力、地震、降雨等因素作用下，随时可能发生崩塌、坠落或滑动，对下方的建筑物、道路、人员等构成严重威胁。基本稳定危岩体在研究区内分布较为广泛，在坝址区右岸B1-2区、B1-3区，左岸以及珍珠岩边坡均有分布。这些区域的危岩体虽然受到地震、风化等作用影响，但程度相对较弱，岩体完整性相对较好。部分坠落式危岩体在天然状态下稳定性系数略大于1.0，处于基本稳定状态，但在地震等工况下，稳定性系数会降低，接近或小于1.0，存在失稳风险。一些滑移式危岩体在天然状态下稳定性系数处于1.0-1.3之间，处于基本稳定状态，但在暴雨、地震等不利工况下，稳定性系数会下降，可能转变为不稳定状态。稳定危岩体相对较少，主要分布在坝址区右岸B1-4区以及左岸和珍珠岩边坡的部分区域。这些区域的危岩体受到的地震震裂损伤和卸荷作用相对较弱，岩体完整性较好，结构面抗剪强度较高。在各种工况下，这些危岩体的稳定性系数均大于1.3，处于稳定状态，在自然条件下发生失稳的可能性较小。但仍需对其进行定期监测，以防因后期环境变化或人类工程活动影响其稳定性。研究区危岩体整体稳定性不容乐观，不稳定和基本稳定的危岩体占比较大，对堰塞湖治理工程和周边居民的生命财产安全构成了较大威胁。因此，需要根据危岩体的稳定性等级，采取相应的防治措施，对不稳定危岩体进行重点治理，对基本稳定危岩体加强监测和防范，以确保红石岩堰塞湖治理工程的安全运行和周边地区的稳定发展。五、危岩体失稳数值模拟与危害性分析5.1数值模拟方法与原理在对红石岩堰塞湖治理工程危岩体失稳进行研究时，非连续变形分析（DDA）法是一种极为重要的数值模拟方法。该方法由美籍华裔科学家石根华博士于1988年提出，专门用于模拟岩体的非连续变形行为。DDA法以离散的块体集合作为模拟对象，充分考虑了岩体变形的非连续和大变形这两个物理本质，这使得它在处理危岩体这类非连续介质问题时具有独特的优势。DDA法的基本原理基于最小势能原理建立总体平衡方程。在该方法中，将刚体位移和块体变形统一考虑，全部块体同步进行求解。具体来说，块体系统的大位移和大变形通过分步迭代计算的小位移和小变形累加来实现。由于每一步都是小位移，所以可以假设每一块体在每一步过程中具有常应力和常应变。块体任一点的位移（u,v）可用6个位移不变量来表示，即（u_0,v_0,\omega,\varepsilon_x,\varepsilon_y,\gamma_{xy}），其中，（u_0,v_0）是块体内特殊点（如块体的重心）的刚体位移；\omega是块体绕转动中心的转动角（以弧度为单位）；\varepsilon_x,\varepsilon_y,\gamma_{xy}是块体的法向应变和切向应变。考虑块体平移、转动（\omega）、正应变（\varepsilon_x,\varepsilon_y）和剪应变（\gamma_{xy}）变形的情况下，取块体系统的全一阶位移模式，块体内任一点的位移可写为：\begin{pmatrix}u\\v\end{pmatrix}=\begin{bmatrix}1&0&-(y-y_0)&1&0&0\\0&1&(x-x_0)&0&1&0\end{bmatrix}\begin{pmatrix}u_0\\v_0\\\omega\\\varepsilon_x\\\varepsilon_y\\\gamma_{xy}\end{pmatrix}其中，x,y为块体内任一点的坐标，x_0,y_0为块体特殊点的坐标。由此则可推导出块体的应变-位移关系矩阵。块体系统的总势能包括块体单元的应变能、初始应力的势能、点荷载和线荷载作用下的势能、体荷载势能、锚杆连接的势能、惯性力势能和粘性力势能等。由最小势能原理，在势能泛函取最小值时系统达到平衡。块体系统的总势能可写成一般形式：U=\frac{1}{2}[D]^T[K][D]-[D]^T[F]其中，[D]为块体的位移向量，[K]为刚度矩阵，[F]为荷载向量。非连续变形分析的平衡方程式由总势能最小化原理来建立，即由各种力和应力产生的总势能推导得到平衡方程式：\frac{\partialU}{\partial\{D\}}=0经推导可简化为[K]\{D\}=[F]。引入边界条件和块体系统的运动学条件，即可对上述方程求解，得到每一个块体的位移与变形状态。在危岩体模拟中，DDA法能够很好地模拟块体间的滑动、张开和闭合，这对于分析危岩体的失稳模式和运动特征至关重要。在模拟坠落式危岩体时，DDA法可以清晰地展示危岩体从母岩脱离后的坠落过程，包括坠落的轨迹、速度以及与坡面碰撞后的运动变化。对于滑移式危岩体，DDA法能够准确模拟其沿结构面的滑动过程，分析滑动过程中的受力情况和变形特征。通过DDA法的模拟，可以直观地了解危岩体在不同工况下的稳定性变化，为危岩体的防治措施设计提供科学依据。5.2典型危岩体失稳模拟以WYT11危岩体为例，其位于红石岩堰塞湖坝址区右岸，处于B1-1区域，高程约1280m。该危岩体呈块状，体积约为100m³，重度为26kN/m³，自重W=100×26=2600kN。它受两组结构面切割，底部与母岩的连接较为薄弱，后缘发育有陡倾裂隙，裂隙深度h=4m，岩体粘聚力标准值c=40kPa，内摩擦角标准值φ=28°。地震水平系数ζe取0.05，地震力Q=ζe×W=0.05×2600=130kN。利用非连续变形分析（DDA）法对WYT11危岩体进行模拟。在建立数值模型时，将危岩体离散为多个块体，根据危岩体的地质特征和现场测量数据，确定块体的几何形状、尺寸和初始位置。考虑到危岩体底部与母岩的连接条件以及后缘裂隙的情况，设置块体之间的接触条件和力学参数，如摩擦系数、粘聚力等。在天然状态下，模拟结果显示危岩体处于相对稳定状态，但后缘裂隙处存在一定的应力集中现象。随着时间的推移，后缘裂隙有逐渐扩展的趋势，但整体危岩体尚未发生失稳。当施加地震作用时，危岩体的稳定性发生明显变化。在地震力的作用下，危岩体开始产生位移，底部与母岩的连接部位受到的剪应力增大。随着地震力的持续作用，后缘裂隙迅速扩展，危岩体底部与母岩的连接逐渐破坏。最终，危岩体从母岩上脱离，发生坠落式失稳。从运动特征来看，危岩体在坠落过程中，速度逐渐增大，其运动轨迹基本沿垂直方向下落，但由于与坡面的碰撞，运动方向发生了一定的改变。通过模拟计算，得到危岩体坠落过程中的速度-时间曲线，在坠落初期，速度增长较慢，随着下落高度的增加，速度迅速增大，在即将落地时，速度达到最大值。对于溢洪洞进口边坡“10.18”崩塌，其成因主要是地震导致岩体结构破碎，在后续的风化、降雨等作用下，岩体的稳定性进一步降低。从失稳模式分析，该崩塌属于坠落-滚（滑）落复合失稳模式。在崩塌初期，岩体受多组结构面切割，在重力作用下从母岩上坠落，坠落过程中与坡面发生碰撞，部分岩体改变运动方向，沿着坡面滚动或滑动，最终形成了较大范围的崩塌堆积体。ZZWY1危岩体位于珍珠岩边坡，其变形破坏特征与岩体的结构和受力状态密切相关。通过数值模拟分析，发现该危岩体在重力和侧向力作用下，底部支撑逐渐减弱，有发生倾倒式失稳的趋势。在模拟过程中，考虑了岩体的自重、地震力以及侧向风力等因素，分析危岩体在不同工况下的稳定性。当危岩体发生失稳后，可能引发涌浪。根据相关理论和经验公式，计算ZZWY1危岩体失稳后可能引发的涌浪高度和传播距离。假设危岩体体积为V，落入水中的速度为v，根据动量守恒原理，可计算出涌浪的初始能量。再结合堰塞湖的水深、湖盆形状等因素，利用涌浪传播模型，计算涌浪在湖中的传播过程。计算结果表明，ZZWY1危岩体失稳后引发的涌浪高度可达3-5m，涌浪传播距离较远，可能对堰塞湖周边的建筑物和设施造成较大破坏。5.3危岩体危害性综合评估综合考虑危岩体失稳的可能性、运动特征和影响范围等多方面因素，对红石岩堰塞湖治理工程危岩体的危害性进行评估，能为后续制定科学合理的防治措施提供关键依据。从失稳可能性来看，研究区危岩体数量众多，共计335处，且分布广泛，坝址区右岸、左岸以及珍珠岩边坡均有分布。通过稳定性计算和数值模拟分析可知，部分危岩体处于不稳定或基本稳定状态，在地震、降雨、风化等因素作用下，有较高的失稳可能性。坝址区右岸B1-1区的危岩体，由于地震震裂损伤强烈，卸荷松弛明显，岩体完整性差，许多坠落式和滑移式危岩体的稳定性系数远小于1.0，处于不稳定状态，失稳可能性极大。在地震等动力作用下，这些危岩体随时可能发生崩塌、坠落或滑动。危岩体的运动特征对其危害性有着重要影响。不同失稳模式的危岩体运动特征各异，坠落式危岩体从母岩脱离后，通常沿垂直方向快速坠落，速度在坠落过程中不断增大，在即将落地时达到最大值。在与坡面碰撞后，运动方向会发生改变，可能会对碰撞点附近的区域造成破坏。滑移式危岩体沿结构面滑动，运动速度相对较慢，但滑动距离可能较长，对滑动路径上的建筑物、道路等基础设施会产生持续的破坏作用。倾倒式危岩体在倾倒过程中，会对周围的岩体产生挤压和破坏作用，引发连锁反应，导致更大范围的岩体失稳，其破坏范围不仅包括倾倒方向的区域，还可能对周边一定范围内的岩体稳定性产生影响。危岩体的影响范围也是评估其危害性的重要因素。红石岩堰塞湖周边地形复杂，人口和基础设施分布较为集中。危岩体一旦失稳，可能会对下方的村庄、公路、堰塞坝等造成严重破坏。坝址区右岸的危岩体如果发生崩塌，可能会直接冲击下方的村庄，造成人员伤亡和房屋损毁。左岸危岩体的滑动可能会堵塞公路，影响交通通行，对堰塞坝的安全也构成威胁，一旦堰塞坝受到破坏，可能引发洪水等次生灾害，影响范围将进一步扩大。珍珠岩边坡危岩体的崩塌可能会对下方的交通要道和居民区造成严重影响，导致交通中断，威胁居民的生命财产安全。综合以上因素，红石岩堰塞湖治理工程危岩体的危害性较大。对于处于不稳定状态的危岩体，如坝址区右岸B1-1区的部分危岩体，应作为重点防治对象，采取有效的防治措施，降低其失稳风险。对于基本稳定的危岩体，虽然在当前条件下相对稳定，但在不利工况下仍有失稳可能，也需要加强监测和防范，及时发现潜在的安全隐患，采取相应的措施进行处理。只有对危岩体的危害性进行全面、准确的评估，并采取针对性的防治措施，才能有效保障堰塞湖治理工程的安全和周边居民的生命财产安全。六、危岩体防治措施研究6.1防治措施原则与目标红石岩堰塞湖治理工程危岩体防治措施的制定，必须紧密围绕保障人民生命财产安全、确保工程设施稳定运行以及维护生态环境平衡等核心要点，遵循一系列科学合理的原则，以实现全面且有效的防治目标。安全第一原则是危岩体防治的首要准则，始终将人民生命财产安全置于首位。在制定和实施防治措施过程中，充分考虑危岩体失稳可能带来的严重后果，采取切实有效的措施，最大限度地降低危岩体对周边居民、建筑物以及基础设施的威胁。对于处于不稳定状态的危岩体，优先采取主动防治措施，增强其稳定性，防止崩塌、坠落等灾害的发生。在坝址区右岸B1-1区，该区域危岩体稳定性差，极易发生失稳，因此对该区域的危岩体进行重点治理，采用放坡开挖、锚喷支护等措施，消除安全隐患，保障居民的生命安全和建筑物的安全。技术可行原则要求防治措施在技术层面具备可操作性和可靠性。结合研究区危岩体的特征，包括分布规律、成因机制、失稳模式以及稳定性状况等，选择合适的防治技术和方法。对于坠落式危岩体，由于其底部与母岩连接薄弱，在重力作用下易坠落，可采用锚固、支撑等主动防治措施，增强其与母岩的连接，提高稳定性。对于滑移式危岩体，可根据其滑动面的特征，采用抗滑桩、锚杆等措施，增加抗滑力，阻止其滑动。在选择防治技术时，还需考虑研究区的地形地貌、地层岩性、地质构造等地质条件，确保技术的可行性。在地形陡峭、施工条件困难的区域，选择操作相对简便、对施工场地要求较低的防治技术。经济合理原则强调在保障防治效果的前提下，尽可能降低防治成本，提高经济效益。对不同的防治措施进行成本效益分析，选择成本低、效果好的方案。在防治措施实施过程中，合理安排施工顺序和施工工艺，避免不必要的浪费。对于一些规模较小、稳定性相对较好的危岩体，可采用简单的防护措施，如设置拦石网、拦石墙等被动防治措施，既能达到防治目的，又能降低成本。而对于规模较大、稳定性差的危岩体，则在保证安全的前提下，优化防治方案，合理选用材料和施工设备，降低工程造价。生态保护原则注重防治措施对生态环境的影响，尽量减少对自然环境的破坏。在施工过程中，采取生态友好的施工方法，减少对植被的破坏，保护生态平衡。对于一些位于生态敏感区域的危岩体，优先采用生态防治措施，如植被护坡等。通过种植根系发达的植物，利用植物根系的加固作用和对坡面径流的调节作用，增强边坡的稳定性，减少危岩体的形成和崩塌，同时还能美化环境，促进生态系统的恢复和保护。红石岩堰塞湖治理工程危岩体防治的目标是通过采取科学合理的防治措施，使危岩体的稳定性得到显著提高，将其对人民生命财产安全和工程设施的威胁降至最低程度。对于不稳定危岩体，通过治理使其达到基本稳定或稳定状态，消除安全隐患；对于基本稳定危岩体，加强监测和防范，确保其在各种工况下的稳定性，防止其转化为不稳定危岩体。通过有效的防治措施，保障堰塞湖治理工程的安全运行，包括溢洪洞、引水隧洞等工程设施的正常运行，避免因危岩体失稳导致工程设施受损，影响堰塞湖的防洪、发电等功能。维护周边地区的生态环境平衡，减少因危岩体防治工程对生态环境的破坏，实现地质灾害防治与生态环境保护的协调发展，为当地的可持续发展创造良好条件。6.2右岸边坡危岩体防治措施6.2.1B1-1区处理措施B1-1区高程在1200-1300m之间，地震震裂损伤强烈，卸荷松弛明显，岩体完整性差，危岩体数量众多且规模较大，以坠落式和滑移式危岩体为主。针对该区域危岩体的特点，考虑采用放坡开挖+锚喷支护的措施进行处理。放坡开挖是一种将基坑边坡进行开挖，形成一定坡度的斜坡，以增加基坑的稳定性和安全性的方法。在B1-1区，通过放坡开挖，可以减小危岩体的高度和坡度，降低其下滑力和倾倒力矩，从而提高危岩体的稳定性。在开挖过程中，根据地质条件和危岩体的实际情况，合理确定放坡坡度，一般可控制在1:1.5-1:2.0之间。采用机械开挖和人工修坡相结合的方式，确保开挖面的平整度和坡度符合设计要求。在开挖过程中，及时对坡面进行防护，防止雨水冲刷和风化作用对坡面的破坏。锚喷支护是通过喷射混凝土、钢筋网、钢绞线等材料对基坑坡面进行加固，以提高基坑的稳定性和安全性的方法。在B1-1区，放坡开挖后，立即进行锚喷支护。首先，在坡面上钻孔，插入锚杆，然后注入水泥浆液，使锚杆与岩体形成摩擦力，增强岩体的整体性和稳定性。锚杆的长度和间距根据危岩体的厚度和稳定性确定，一般长度为3-5m，间距为1.5-2.0m。在锚杆施工完成后，进行钢筋网片绑扎，钢筋网片的网格尺寸一般为200mm×200mm-250mm×250mm。最后，将混凝土喷射到锚杆和钢筋网片上，形成支护结构。喷射混凝土的强度等级一般为C20-C25，喷射厚度为100-150mm。在喷射混凝土过程中，严格控制喷射压力和喷射角度，确保混凝土与岩体紧密结合，提高支护效果。6.2.2B1-2区评价与处理措施B1-2区高程在1300-1400m之间，震裂损伤和卸荷特征相对B1-1区较弱，但仍有较多危岩体分布，以坠落式危岩体居多。针对该区域危岩体的特点，主要采用小扰动清危+锚喷支护+主被动拦挡等措施进行治理。小扰动清危是指在不破坏危岩体整体稳定性的前提下，对危岩体进行清理和加固。在B1-2区，对于一些体积较小、稳定性较差的坠落式危岩体，采用人工或小型机械进行清理，将其从母岩上拆除，消除安全隐患。在清理过程中，采用控制爆破、静态破碎等小扰动方法，避免对周边岩体造成过大的震动和破坏。对于一些体积较大、难以直接拆除的危岩体，采用锚杆、锚索等进行加固，增强其与母岩的连接，提高稳定性。在小扰动清危的基础上，对B1-2区的坡面进行锚喷支护，其施工方法和参数与B1-1区类似，但根据该区域危岩体的特点，可适当调整锚杆的长度和间距，以及喷射混凝土的厚度和强度等级。在锚喷支护完成后，为防止危岩体崩塌对下方造成危害，在坡脚设置主被动拦挡结构。被动拦挡结构可采用拦石墙，拦石墙采用浆砌片石或混凝土浇筑而成，墙高一般为2-3m，墙厚根据计算确定，一般为0.5-1.0m。拦石墙背设置填土缓冲层，以吸收危岩体崩塌时的冲击力。主动拦挡结构可采用钢丝绳网、铁丝格栅网等柔性防护网系统，将其固定在坡面上，抑制危岩体的崩塌和风化剥落，限制少量及局部落石的运动范围。6.2.3B1-3区评价与处理措施B1-3区高程在1400-1500m之间，震裂损伤和卸荷特征进一步减弱，但由于岩体风化作用较强，仍有部分危岩体存在，多为小型坠落式危岩体。针对该区域危岩体的特点，采用清理危石+坡面防护的措施进行治理。对于B1-3区的小型坠落式危岩体，采用人工清理的方式，将其从母岩上拆除，运至安全地带。在清理过程中，注意安全防护，防止危岩体坠落对人员造成伤害。清理完成后，对坡面进行防护，以防止风化作用进一步破坏岩体。坡面防护可采用抹面、捶面等方法，在坡面上涂抹水泥砂浆或水泥石灰砂浆，厚度一般为50-100mm。抹面、捶面可有效防止坡面岩体风化剥落，保护岩体的完整性。对于一些风化较为严重的区域，可在抹面、捶面的基础上，设置锚杆进行加固，增强坡面的稳定性。6.2.4B1-4区评价与处理措施B1-4区位于右岸顶部，高程大于1500m，虽然地震震裂损伤相对较弱，但由于长期的风化作用和岩体卸荷，也有一定数量的危岩体分布，以倾倒式危岩体较为常见。针对该区域危岩体的特点，采用支撑加固+坡面防护的措施进行治理。对于B1-4区的倾倒式危岩体，采用浆砌片石、混凝土支顶等方法进行支撑加固。在危岩体底部设置支撑结构，增加危岩体的底部支撑面积，提高其抗倾倒能力。支撑结构的尺寸和强度根据危岩体的重量和稳定性计算确定，确保能够承受危岩体的重力和倾倒力矩。在支撑加固的同时，对坡面进行防护，可采用护面墙的方式，护面墙采用浆砌片石或混凝土浇筑而成，墙高和墙厚根据坡面情况确定，一般墙高为1-2m，墙厚为0.3-0.5m。护面墙可防止坡面岩体风化剥落，保护危岩体的稳定性。为增强危岩体与支撑结构和坡面的连接，可设置普通长锚杆，将危岩体、支撑结构和坡面岩体连接在一起，提高整体稳定性。6.3左岸边坡危岩体防治措施左岸边坡危岩体的分布呈现出特定的规律，在顶部裂缝发育部位，由于岩体受拉应力作用，形成了一些规模较大的滑移式危岩体；而在边坡中部，由于岩体风化和卸荷作用，部分岩体与母岩分离，形成坠落式危岩体。这些危岩体虽然稳定性相对右岸较好，但下方为公路和堰塞坝，一旦失稳，将对交通和堰塞坝安全构成严重威胁。针对左岸边坡顶部裂缝发育部位的危岩体，采用放坡开挖+锚喷支护的措施进行处理。放坡开挖能够减小危岩体的高度和坡度，降低其下滑力和倾倒力矩，从而提高危岩体的稳定性。在开挖过程中，依据地质条件和危岩体的实际情况，合理确定放坡坡度，一般控制在1:1.5-1:2.0