三峡库区马家沟滑坡在库水位涨落影响下的稳定性机制与评价研究

三峡库区马家沟滑坡在库水位涨落影响下的稳定性机制与评价研究一、引言1.1研究背景与意义三峡工程作为世界上最大的水利枢纽工程之一，自建成蓄水以来，在防洪、发电、航运、水资源利用等方面发挥了巨大的综合效益，对中国的经济发展和社会稳定起到了至关重要的作用。然而，随着三峡水库水位的周期性涨落，库区周边地质环境发生了显著变化，滑坡等地质灾害的发生频率和规模也有所增加。三峡水库正常运行时，水位在145-175米之间周期性变化，年变幅达30米。这种大幅度的水位涨落，使得库岸岩土体长期处于饱水与风干的交替状态，力学性质不断劣化。同时，水位升降过程中产生的动水压力、渗透压力等，也对库岸稳定性产生了极大的影响。据不完全统计，三峡库区在175米库水位影响的范围内共有大小滑坡1190余个。2003年三峡水库蓄水后，诸多滑坡的部分区域已被水浸没，在库水位175-145米周期涨跌时，滑坡受到库水动力作用，岩土体力学强度降低，有效重力减轻，还会承受动水荷载，加速了滑坡的破坏进程。马家沟滑坡位于三峡库区，地处长江左岸卧牛山麓，长江支流吒溪河左岸，距长江河口2.1km。该滑坡面积约1.5平方公里，为黄土石积体滑坡，呈鱼鳞状，南北向展布，顶部海拔约725米，最大高度达150米。马家沟滑坡发育于第四纪黄土和更新世高洞段悬崖受侵蚀的基础上，向水库方向滑去。2003年长江三峡水库蓄水至135米后的3个月内，滑体后缘出现了拉张裂缝，虽其后拉裂变形趋于稳定，但表明该滑坡对水库蓄水反应敏感，在水位继续变化时，有复活的可能性。该滑坡体上有居民居住，一旦失稳滑动，将会对附近的居民点、公路、桥梁、通讯线路、架空线路、耕地、林木、河道等产生较大的危害，甚至可能危及附近其他滑坡和诱发河对岸的滑坡，严重威胁人民生命财产安全和基础设施的稳定运行。研究库水位涨落作用下马家沟滑坡的稳定性具有重大的现实意义。从保障人民生命财产安全角度看，准确掌握马家沟滑坡在不同库水位条件下的稳定性状况，能够为滑坡灾害的预警和防治提供科学依据，提前采取有效的防范措施，避免或减少人员伤亡和财产损失。在指导工程防治工作方面，通过对滑坡稳定性的研究，可以深入了解滑坡的变形机制和影响因素，从而优化防治工程的设计和施工方案，提高防治工程的效果和可靠性，降低治理成本。对马家沟滑坡稳定性的研究成果，还能为三峡库区其他类似滑坡的研究和防治提供参考和借鉴，丰富和完善库水位涨落作用下滑坡稳定性研究的理论和方法体系，推动地质灾害防治领域的科学技术进步。1.2国内外研究现状库水位涨落作用下滑坡稳定性的研究一直是国内外地质灾害领域的重点和热点问题。国外在该领域的研究起步较早，积累了丰富的理论和实践经验。20世纪60年代，意大利瓦依昂水库滑坡事件后，学者们开始重视库水位变化对滑坡稳定性的影响。众多学者通过物理模型试验、数值模拟等方法，对库水位涨落过程中滑坡体的渗流场、应力场变化规律以及稳定性演化机制进行了深入研究。有学者利用有限元软件对某水库滑坡进行模拟，分析了不同库水位升降速率下滑坡体的渗流特性和稳定性系数变化，结果表明水位下降速率越大，滑坡稳定性降低越明显。还有学者通过室内模型试验，研究了库水位周期性涨落对滑坡体岩土力学性质的影响，发现长期的水位变化会导致岩土体强度降低，进而影响滑坡稳定性。国内在库水位涨落与滑坡稳定性关系的研究方面也取得了丰硕成果。特别是三峡工程建设以来，大量科研人员针对三峡库区滑坡开展了全方位的研究。在理论研究上，完善了考虑库水作用的滑坡稳定性计算方法，如改进的极限平衡法、强度折减法等，并将其应用于实际工程分析中。在数值模拟方面，运用多种先进的数值计算软件，建立了复杂的滑坡模型，模拟库水位涨落、降雨、地震等多种因素耦合作用下的滑坡稳定性变化过程。在监测预警方面，建立了一套完善的滑坡监测体系，通过地面监测、卫星遥感、InSAR等技术手段，实时获取滑坡体的变形信息，为滑坡稳定性分析和灾害预警提供数据支持。针对马家沟滑坡稳定性的研究，国内已有一些相关成果。有研究通过现场勘查和室内试验，分析了马家沟滑坡的工程地质条件和岩土力学参数，运用剩余推力传递系数法计算了不同工况下的滑坡稳定性系数，得出了人类活动、降雨以及库水位升降等因素对滑坡稳定性的影响规律。还有学者利用有限元数值模拟方法，研究了库水位涨落过程中马家沟滑坡体内的渗流场和应力场变化，发现水位上升时，滑坡体前缘受到较大的浮托力和动水压力作用，稳定性降低；水位下降时，由于渗流作用，滑坡体后部产生较大的渗透力，也会对稳定性产生不利影响。尽管国内外在库水位涨落作用下滑坡稳定性研究以及马家沟滑坡稳定性研究方面取得了一定进展，但仍存在一些不足之处。现有研究多集中在单一因素对滑坡稳定性的影响，对于库水位涨落与降雨、地震等多种因素耦合作用下的滑坡稳定性研究还不够深入，耦合模型和计算方法有待进一步完善。在滑坡稳定性评价指标和判据方面，尚未形成统一的标准，不同方法计算结果存在一定差异，给实际工程应用带来困难。马家沟滑坡的研究在监测数据的长期连续性和全面性上还有所欠缺，对滑坡深部变形机制和破坏模式的认识还不够清晰，需要加强长期监测和深入研究，以提高对该滑坡稳定性的认识和预测能力。1.3研究内容与方法1.3.1研究内容本研究将深入剖析库水位涨落作用下马家沟滑坡的稳定性，具体研究内容如下：滑坡工程地质条件研究：对马家沟滑坡的地形地貌、地层岩性、地质构造、水文地质条件等进行详细调查与分析，通过现场勘查、地质测绘、钻探、物探等手段，获取滑坡体的几何形态、岩土体结构、地下水分布等基础数据，为后续的稳定性分析提供依据。通过现场勘查，准确绘制滑坡的边界、地形等高线，确定滑坡的规模和范围；利用钻探技术，获取不同深度的岩土体样本，进行室内物理力学性质试验，测定岩土体的密度、含水率、抗剪强度等参数；运用物探方法，如地质雷达、地震波速测试等，探测滑坡体内部的结构和地质缺陷。库水位涨落规律及对滑坡的作用机制研究：收集三峡水库历年的水位数据，分析库水位的涨落过程、变化幅度、升降速率等规律。研究库水位涨落过程中，滑坡体所受到的浮托力、动水压力、渗透压力等作用的变化规律，以及这些力对滑坡体应力应变状态的影响机制。建立库水位涨落与滑坡体渗流场、应力场之间的耦合关系模型，揭示库水位涨落作用下滑坡变形破坏的内在机制。滑坡稳定性计算与评价：基于上述研究，选取合适的滑坡稳定性计算方法，如极限平衡法、有限元强度折减法等，考虑库水位涨落、降雨、地震等多种因素，对马家沟滑坡在不同工况下的稳定性进行计算分析。确定滑坡的稳定状态、潜在滑动面位置、安全系数等关键指标，依据相关规范和标准，对滑坡的稳定性进行评价，判断其是否处于稳定状态以及失稳的可能性大小。滑坡稳定性影响因素敏感性分析：针对影响马家沟滑坡稳定性的主要因素，如库水位升降速率、岩土体抗剪强度参数、地下水位变化等，开展敏感性分析。通过改变各因素的取值，计算滑坡稳定性系数的变化情况，确定各因素对滑坡稳定性的影响程度和敏感性排序，找出对滑坡稳定性影响最为显著的因素，为滑坡防治措施的制定提供重点方向。滑坡稳定性预测与防治措施研究：结合滑坡的历史变形监测数据和稳定性分析结果，运用时间序列分析、灰色预测等方法，对马家沟滑坡在未来库水位涨落条件下的稳定性变化趋势进行预测。根据预测结果，提出针对性的滑坡防治措施，如排水工程、支挡工程、减重与反压等，并对防治措施的效果进行模拟分析和评估，确保防治措施的有效性和可靠性。1.3.2研究方法为实现上述研究目标，本研究拟采用以下多种研究方法：现场调查与监测：通过现场勘查，对马家沟滑坡的地形地貌、地质构造、岩土体特征等进行详细观察和记录，绘制地质草图和剖面图。在滑坡体上布置监测点，运用全站仪、GPS、水准仪等测量仪器，对滑坡体的地表位移、垂直位移进行定期监测；采用钻孔倾斜仪监测滑坡体深部的位移变化；利用孔隙水压力计监测孔隙水压力；通过地下水位监测井监测地下水位变化。通过长期的监测，获取滑坡体在库水位涨落过程中的变形数据，为稳定性分析和预测提供数据支持。室内试验：采集滑坡体的岩土体样本，在实验室进行物理力学性质试验，包括密度试验、含水率试验、颗粒分析试验、三轴压缩试验、直剪试验等，测定岩土体的基本物理参数和抗剪强度参数，为稳定性计算提供准确的参数依据。数值模拟：运用有限元软件，如ANSYS、ABAQUS等，建立马家沟滑坡的三维数值模型，模拟库水位涨落、降雨、地震等因素作用下，滑坡体的渗流场、应力场和位移场的变化规律，分析滑坡的稳定性变化过程。通过数值模拟，可以直观地展示滑坡体在不同工况下的力学响应，深入研究滑坡的变形破坏机制。理论分析：基于土力学、岩石力学、渗流力学等相关理论，对库水位涨落作用下滑坡的稳定性进行理论分析。运用极限平衡理论，计算滑坡在不同工况下的安全系数；利用渗流理论，分析库水位变化时滑坡体内的渗流特性；依据岩土力学本构模型，描述滑坡体的应力应变关系，为数值模拟和工程实践提供理论指导。综合分析与对比研究：将现场调查、监测数据、室内试验结果、数值模拟和理论分析成果进行综合分析，对比不同方法得到的结果，相互验证和补充，全面、准确地认识库水位涨落作用下马家沟滑坡的稳定性特征和变形破坏规律，为滑坡防治提供科学合理的建议。二、马家沟滑坡工程地质条件2.1区域地质背景马家沟滑坡所在的三峡库区，在大地构造位置上处于扬子准地台的上扬子台褶带，其地质构造复杂，经历了多期构造运动的叠加影响。区内褶皱构造发育，主要褶皱轴向呈北东-南西向，如秭归向斜等。这些褶皱构造控制了区域地层的展布和山体的形态，为马家沟滑坡的形成提供了宏观的地质背景。褶皱构造使得地层发生弯曲变形，在褶皱的翼部和转折端，岩石的完整性受到破坏，节理、裂隙发育，增加了岩土体的透水性和变形能力，从而降低了山体的稳定性，为滑坡的孕育创造了条件。在断裂构造方面，库区存在多条规模较大的断裂，如北北东向的水田坝断裂等。断裂的存在不仅破坏了岩石的连续性和完整性，还改变了区域的应力场分布。断裂附近的岩石破碎，形成了软弱结构面，当受到库水位涨落、降雨等外力作用时，这些软弱结构面容易成为滑坡的滑动面。此外，断裂活动还可能引发地震，进一步加剧滑坡的稳定性问题。例如，历史上的一些地震活动，使得断裂附近的山体出现松动，在后续的库水位变化过程中，更容易发生滑坡灾害。三峡库区的地层岩性较为复杂，出露的地层从元古界到新生界均有分布。马家沟滑坡所在区域主要出露的地层有第四系全新统滑坡堆积层（Q4del）、残坡积层（Q4el+dl），以及侏罗系上统沙溪庙组（J3s）地层。第四系全新统滑坡堆积层主要由块石、碎石、粉质黏土等组成，结构松散，透水性较强，在库水位涨落和降雨作用下，容易发生变形和失稳。残坡积层则由粉质黏土、碎石等组成，厚度变化较大，其力学性质相对较差，对滑坡的稳定性也有一定影响。侏罗系上统沙溪庙组地层主要为砂岩、泥岩互层，砂岩强度较高，但泥岩遇水后容易软化、崩解，导致岩体的整体强度降低。在库水位周期性涨落过程中，泥岩的软化作用会使岩体内部的应力分布发生改变，从而增加滑坡发生的风险。这种地层岩性的组合特征，使得马家沟滑坡所在区域的岩土体在力学性质上存在明显差异，在外部因素作用下，容易产生不均匀变形，进而引发滑坡。2.2马家沟滑坡基本特征马家沟滑坡位于湖北省秭归县归州镇澎家坡村8组，地处长江左岸卧牛山麓，长江支流吒溪河左岸，距长江河口2.1km。该滑坡地貌形态呈近东西向展布，坐落于向吒溪河倾斜的斜坡之上。坡体整体为顺坡向，坡面形态呈现为陡缓相间的折线型，大致呈上陡26°-35°，中缓10°-15°，下陡30°-45°的特征。滑坡前缘受河流侵蚀作用，形成了河流相冲积阶地。从规模上看，马家沟滑坡面积约1.5平方公里，属于规模较大的滑坡。其最大高度达到150米，顶部海拔约725米，顶部的最高处比滑坡头高出20余米，平均坡度10度。滑坡体呈鱼鳞状，南北向展布，这种独特的形态是在长期的地质演化过程中，受地形、岩土体性质以及外力作用等多种因素共同影响而形成的。在物质组成方面，马家沟滑坡为黄土石积体滑坡，主要由第四系全新统滑坡堆积层（Q4del）和残坡积层（Q4el+dl）构成。滑坡堆积层中包含块石、碎石、粉质黏土等，这些物质大小混杂，结构松散，透水性较强。块石和碎石的存在使得滑坡体具有一定的骨架结构，但粉质黏土的填充又增加了其孔隙性和软弱性。残坡积层则主要由粉质黏土和碎石组成，厚度变化较大，其力学性质相对较差，在外界因素作用下容易发生变形。从地质结构上分析，滑坡体下伏基岩为侏罗系上统沙溪庙组（J3s）地层，主要为砂岩、泥岩互层。砂岩具有较高的强度和抗风化能力，但泥岩遇水后容易软化、崩解，导致岩体的整体强度降低。这种软硬相间的地层结构，在库水位涨落和降雨等作用下，容易产生不均匀变形，进而引发滑坡。泥岩的软化会使岩体内部的应力分布发生改变，在砂岩与泥岩的接触面处，容易形成软弱结构面，成为滑坡潜在的滑动面。同时，滑坡体与下伏基岩之间的界面也可能因岩土体性质的差异和长期的地质作用，存在一定的软弱带，对滑坡的稳定性产生不利影响。2.3水文地质条件马家沟滑坡区域的地下水类型主要有松散岩类孔隙水和基岩裂隙水。松散岩类孔隙水赋存于滑坡堆积层和残坡积层中，由于这些地层结构松散，孔隙度较大，透水性较强，因此孔隙水含量较为丰富。孔隙水的补给来源主要为大气降水和地表水的入渗，在降雨过程中，雨水迅速渗入地下，补充孔隙水。同时，库水位的涨落也会对孔隙水产生影响，当库水位上升时，河水通过滑坡前缘渗入地下，使孔隙水水位升高；当库水位下降时，孔隙水则向河流方向排泄。基岩裂隙水主要赋存于下伏基岩的裂隙中，侏罗系上统沙溪庙组砂岩、泥岩互层的地层中，砂岩裂隙相对发育，是基岩裂隙水的主要储存空间。泥岩由于其隔水性能，对基岩裂隙水的分布和运移起到一定的阻隔作用，导致基岩裂隙水在空间上的分布不均匀。基岩裂隙水的补给来源同样包括大气降水和孔隙水的下渗，其排泄方式主要为向河流排泄和通过基岩裂隙向深部运移。三峡水库正常运行时，水位在145-175米之间周期性变化，年变幅达30米。这种大幅度的库水位涨落对马家沟滑坡的水文地质条件产生了显著影响。在库水位上升阶段，滑坡前缘被淹没，库水通过岩土体的孔隙和裂隙向滑坡体内渗透，导致滑坡体地下水位迅速上升。地下水位的上升使岩土体饱水程度增加，重度增大，有效应力减小，抗剪强度降低。同时，库水的渗透还会产生渗透压力，该压力方向指向滑坡体内部，对滑坡体产生一定的推力，不利于滑坡的稳定。当库水位下降时，滑坡体中的地下水向库水方向排泄，但由于岩土体的渗透性有限，地下水的排泄速度相对较慢，导致滑坡体内部形成水头差，产生渗透力。渗透力的方向与地下水的渗流方向一致，即从滑坡体内部指向库水方向，在滑坡体后部形成向上的渗透力，使岩土体的有效重量减轻，抗滑力降低；在滑坡体前缘形成向下的渗透力，增加了下滑力。这种由于库水位下降产生的渗透力，对滑坡的稳定性极为不利，容易引发滑坡的滑动。库水位的周期性涨落还会使滑坡体岩土体长期处于饱水与风干的交替状态，导致岩土体的物理力学性质发生劣化。长期饱水会使岩土体中的黏土矿物发生膨胀、软化，颗粒之间的胶结力减弱，抗剪强度降低；风干过程则会使岩土体产生干裂，增加其透水性和孔隙率，进一步降低岩土体的强度。这种反复的饱水与风干作用，使得滑坡体的稳定性逐渐降低，增加了滑坡发生的风险。三、库水位涨落作用下的滑坡稳定性影响因素3.1库水位涨落特征三峡水库自2003年开始蓄水运行，其水位涨落呈现出独特的规律。在正常运行情况下，每年的水位变化过程可分为消落期、蓄水期和平稳运行期。每年汛期（6-9月），为了预留防洪库容，应对长江上游可能发生的洪水，三峡水库会将水位降至145米左右，这一过程即为消落期。在消落期，库水位以一定的速率逐渐下降，平均下降速率约为0.3-0.5米/天。汛期过后，大约从9月底或10月初开始，水库进入蓄水期，通过控制下泄流量，使库水位逐步上升，直至11月份左右达到正常蓄水位175米。蓄水期的库水位上升速率相对较为稳定，平均上升速率约为0.3-0.4米/天。在11月至次年5月期间，三峡水库保持在较高水位运行，水位波动相对较小，处于平稳运行期。然而，即使在平稳运行期，库水位也并非完全静止不变，会受到上游来水、发电用水等因素的影响，在一定范围内波动，波动幅度一般在1-2米左右。三峡水库水位的年变幅高达30米，这种大幅度的水位变化在世界大型水库中较为罕见。如此大的水位变幅，使得库岸岩土体在一年中经历了从饱水到风干的剧烈变化过程。在水位上升阶段，岩土体逐渐被水淹没，长时间处于饱水状态；而在水位下降阶段，岩土体又逐渐露出水面，经历风干过程。这种周期性的饱水与风干交替作用，对岩土体的物理力学性质产生了深远影响。库水位的升降速率也是影响滑坡稳定性的重要因素之一。在消落期和蓄水期，库水位的升降速率虽然相对较为稳定，但不同年份之间仍存在一定差异。当库水位快速下降时，滑坡体内部的地下水来不及及时排出，形成较大的水头差，从而产生较大的渗透力，对滑坡的稳定性产生不利影响。相反，库水位上升时，虽然会使滑坡体受到浮托力作用，有效重量减轻，但上升速率相对较慢，对滑坡稳定性的影响相对较小。不过，如果库水位上升过程中伴随着强降雨等其他不利因素，也可能会加剧滑坡的变形和失稳。3.2岩土体物理力学性质马家沟滑坡的岩土体主要包括第四系全新统滑坡堆积层（Q4del）、残坡积层（Q4el+dl）以及下伏的侏罗系上统沙溪庙组（J3s）地层。这些岩土体在库水位涨落作用下，物理力学性质发生了显著变化。在抗剪强度方面，第四系全新统滑坡堆积层和残坡积层的抗剪强度受库水位影响较大。库水位上升时，岩土体饱水，孔隙水压力增大，有效应力减小，抗剪强度降低。相关研究表明，饱水状态下的粉质黏土抗剪强度指标内摩擦角可能降低10%-20%，黏聚力降低20%-30%。对于滑坡堆积层中的块石和碎石，其与粉质黏土之间的咬合作用也会因饱水而减弱，进一步降低了整个堆积层的抗剪强度。当库水位下降时，岩土体中的水分排出，土体产生收缩变形，内部结构发生改变，也会导致抗剪强度下降。长期的库水位涨落，使得岩土体反复经历饱水与风干过程，颗粒之间的胶结物逐渐被溶解或冲刷，抗剪强度持续劣化。侏罗系上统沙溪庙组地层中的砂岩抗剪强度相对较高，但泥岩的抗剪强度较低且对水的敏感性强。库水位变化时，泥岩遇水软化，其抗剪强度大幅降低，软化后的泥岩内摩擦角可降低30%-40%，黏聚力降低50%-60%。这种泥岩抗剪强度的降低，会使整个地层的稳定性受到影响，尤其是在砂岩与泥岩的接触面处，容易形成滑动面。岩土体的重度在库水位涨落过程中也发生变化。库水位上升，岩土体饱水，重度增大。对于第四系全新统滑坡堆积层和残坡积层，饱水后重度可增加10%-20%。重度的增加会使滑坡体的下滑力增大，对滑坡稳定性产生不利影响。而在库水位下降过程中，虽然岩土体中的水分排出，重度有所减小，但由于之前饱水状态下岩土体结构的改变，其力学性质已经发生劣化，对滑坡稳定性的改善作用有限。此外，库水位涨落还会影响岩土体的渗透性。长期的饱水与风干交替作用，使岩土体的孔隙结构发生变化。在饱水过程中，土体中的细小颗粒可能会被水携带迁移，堵塞部分孔隙，导致渗透性降低；而风干过程中，岩土体干裂，又会增加新的孔隙和裂隙，使渗透性增大。这种渗透性的复杂变化，进一步影响了滑坡体在库水位涨落过程中的渗流特性和稳定性。3.3渗流作用库水位涨落引发的渗流对马家沟滑坡稳定性有着至关重要的影响，其作用机制主要涉及动水压力和孔隙水压力等方面。在库水位变化过程中，动水压力是不可忽视的关键因素。当库水位上升时，水流速度相对较慢，动水压力主要表现为对滑坡体的浮托力。随着水位的升高，滑坡体前缘被水淹没的范围逐渐扩大，浮托力作用面积增大。以马家沟滑坡为例，在库水位从145米上升至175米的过程中，滑坡体前缘的浮托力逐渐增大，使得滑坡体的有效重量减轻。根据阿基米德原理，浮托力等于排开液体的重量，即F_{浮}=\rhogV，其中\rho为水的密度，g为重力加速度，V为排开液体的体积。在滑坡体被水淹没的部分，浮托力的增加导致滑坡体与下伏基岩之间的正压力减小，从而降低了抗滑力。相关研究表明，浮托力的作用可使滑坡体抗滑力降低10%-20%，对滑坡稳定性产生一定的不利影响。当库水位下降时，情况则更为复杂。由于水位下降速度相对较快，滑坡体内部的地下水来不及及时排出，形成了水头差，进而产生了动水压力。动水压力的方向与水流方向一致，即从滑坡体内部指向库水方向。在马家沟滑坡中，当库水位快速下降时，滑坡体后部的地下水向库水方向渗流，在滑坡体后部形成向上的动水压力，使岩土体的有效重量减轻，抗滑力降低；在滑坡体前缘形成向下的动水压力，增加了下滑力。有研究通过数值模拟分析得出，库水位快速下降时，动水压力可使滑坡体的下滑力增加15%-25%，对滑坡稳定性的影响十分显著。这种动水压力的作用，在一定程度上改变了滑坡体的应力分布状态，增加了滑坡失稳的风险。孔隙水压力在库水位涨落过程中也发生着明显变化。库水位上升时，水迅速渗入滑坡体，导致孔隙水压力快速升高。在马家沟滑坡的第四系全新统滑坡堆积层和残坡积层中，由于其结构松散，孔隙度较大，水的渗入使得孔隙水压力在短时间内大幅增加。孔隙水压力的升高会减小岩土体的有效应力，根据有效应力原理，\sigma'=\sigma-u，其中\sigma'为有效应力，\sigma为总应力，u为孔隙水压力。有效应力的减小直接导致岩土体抗剪强度降低，从而影响滑坡的稳定性。相关室内试验研究表明，库水位上升引起孔隙水压力增加，可使粉质黏土的抗剪强度降低20%-30%，对滑坡稳定性产生较大的负面影响。当库水位下降时，滑坡体中的孔隙水开始向外排出，但由于岩土体的渗透性有限，孔隙水压力并不会立即恢复到初始状态，而是存在一定的滞后现象。在马家沟滑坡中，这种滞后现象使得滑坡体内部在库水位下降后一段时间内仍保持较高的孔隙水压力。较高的孔隙水压力持续减小有效应力，降低岩土体抗剪强度，增加了滑坡滑动的可能性。研究发现，库水位下降后的滞后阶段，孔隙水压力对滑坡稳定性的影响可持续数天甚至数周，具体时长取决于岩土体的渗透性和水位下降速率等因素。动水压力和孔隙水压力的综合作用，使得马家沟滑坡在库水位涨落过程中面临着严峻的稳定性挑战。在库水位上升和下降的不同阶段，这两种力的相互作用导致滑坡体的应力应变状态不断变化，抗滑力和下滑力的平衡关系被打破，从而增加了滑坡失稳的风险。因此，深入研究渗流作用下动水压力和孔隙水压力的变化规律及其对滑坡稳定性的影响，对于准确评估马家沟滑坡的稳定性和制定有效的防治措施具有重要意义。3.4其他影响因素除了库水位涨落、岩土体物理力学性质和渗流作用等主要因素外，降雨、地震和人类活动等因素与库水位涨落共同作用，也对马家沟滑坡的稳定性产生着重要影响。降雨是影响马家沟滑坡稳定性的关键外部因素之一，它与库水位涨落的耦合作用较为复杂。降雨通过入渗增加滑坡体的含水量，使岩土体重度增大，下滑力增加。雨水入渗还会使孔隙水压力升高，有效应力减小，抗剪强度降低。当降雨与库水位上升同时发生时，滑坡体饱水程度加剧，地下水水位大幅上升，岩土体抗剪强度进一步降低，滑坡稳定性受到极大威胁。有研究表明，在库水位较高且持续降雨的情况下，滑坡发生的概率比单一因素作用时增加了30%-50%。若降雨发生在库水位下降阶段，由于库水位下降导致的渗透力与降雨入渗产生的孔隙水压力变化相互叠加，会进一步改变滑坡体的应力状态。降雨入渗使得滑坡体后部孔隙水压力升高，与库水位下降产生的向上渗透力共同作用，使后部岩土体有效重量减轻，抗滑力降低；而在滑坡体前缘，降雨入渗和库水位下降产生的向下渗透力叠加，增加了下滑力。这种情况下，滑坡更容易发生失稳滑动。相关案例分析显示，某类似滑坡在库水位下降过程中遭遇强降雨，最终发生了大规模滑动，造成了严重的灾害损失。地震作为一种突发性的强烈动力作用，与库水位涨落共同作用时，对马家沟滑坡稳定性的影响更为显著。地震产生的地震力会使滑坡体产生惯性力，增加下滑力。同时，地震还会使岩土体结构受到破坏，孔隙率增大，抗剪强度降低。当库水位处于较高位置时，地震作用下饱和岩土体的抗剪强度降低更为明显，滑坡更容易失稳。研究表明，在地震和库水位共同作用下，滑坡稳定性系数比静态条件下降低20%-40%。在地震作用下，库水位涨落引起的渗流场变化也会加剧。地震使岩土体的渗透性改变，导致地下水渗流更加复杂。库水位上升时，地震作用可能使岩土体中的孔隙和裂隙扩张，加速库水的渗入，进一步提高孔隙水压力；库水位下降时，地震造成的岩土体结构破坏会使地下水排泄不畅，增加滑坡体内部的水头差，增大渗透力。这种渗流场的变化与地震力的叠加，极大地增加了滑坡失稳的风险。历史上的一些地震引发的滑坡灾害中，库水位较高的区域滑坡发生的规模和危害程度明显更大。人类活动在马家沟滑坡区域也较为频繁，对滑坡稳定性产生了不可忽视的影响。在滑坡体上进行的工程建设，如修建房屋、道路等，改变了滑坡体的原始地形地貌和应力分布状态。建筑物的加载增加了滑坡体的重量，使下滑力增大；道路开挖可能切断了滑坡体的排水通道，导致地下水水位上升，影响滑坡稳定性。有研究指出，不合理的工程建设可使滑坡稳定性系数降低10%-20%。在滑坡体周边进行的农业灌溉活动，也会对滑坡稳定性产生影响。过量的灌溉水渗入地下，增加了滑坡体的含水量，改变了地下水的渗流状态。当灌溉与库水位涨落同时作用时，可能导致滑坡体地下水位长时间处于高位，岩土体长期处于饱水状态，抗剪强度持续降低。此外，植被破坏也是人类活动影响滑坡稳定性的一个方面。植被具有固土护坡、调节地表径流和涵养水源的作用，植被破坏后，滑坡体的抗侵蚀能力减弱，地表径流增加，雨水更容易渗入地下，对滑坡稳定性产生不利影响。四、库水位涨落作用下滑坡稳定性分析方法4.1极限平衡法极限平衡法是基于静力平衡原理来分析边坡在各种潜在破坏模式下的受力状态，通过考量边坡滑体上的抗滑力与下滑力之间的关系，以此来评价边坡的稳定性。该方法在工程领域应用广泛，其核心思想是假定岩土体中存在一个潜在的破坏面，将破坏面内的土体视为脱离体，通过计算作用于脱离体上的力系，使其达到静力平衡时所需的岩土抗力或抗剪强度，并与破坏面实际能够提供的抗力或抗剪强度进行对比，从而求得稳定性安全系数。在极限平衡法中，根据不同的假设条件和计算方式，衍生出了多种具体的分析方法，如费伦纽斯（Fellenius）法、毕肖普（Bishop）法、泰勒（Taylor）法、简布（Janbu）法、摩根斯顿-普赖斯（Morgenstern-Price）法、斯潘塞（Spencer）法、萨尔玛（Sarma）法、楔形体法、平面直线法、传递系数法，以及贝克-加伯（Baker-Garber）临界滑面法等。以费伦纽斯法为例，其假设滑动面为圆弧面，将滑体分成若干个垂直土条，不考虑土条间的相互作用力。对于每个土条，根据静力平衡条件，计算作用在土条上的下滑力和抗滑力。下滑力主要由土条的自重和可能存在的外部荷载产生，抗滑力则由土条滑动面上的黏聚力和摩擦力提供。整个滑体的稳定性系数即为抗滑力与下滑力对滑弧圆心的力矩之比。在实际应用中，通过试算不同的圆弧滑动面，找出最小稳定性系数对应的滑动面，该系数即为边坡的稳定性系数。对于马家沟滑坡，在运用极限平衡法进行稳定性分析时，首先需要通过现场勘查、钻探、室内试验等手段，获取滑坡体的详细几何参数，包括滑坡的边界、厚度、坡度等，以及岩土体的物理力学参数，如重度、黏聚力、内摩擦角等。这些参数的准确获取对于稳定性分析的结果至关重要。针对不同的工况，如天然状态、库水位上升、库水位下降、降雨、地震等，分别考虑相应的荷载作用。在库水位上升工况下，需要考虑滑坡体前缘受到的浮托力作用，浮托力会减小滑坡体的有效重量，从而降低下滑力，但同时也会减小抗滑力。在库水位下降工况下，除了考虑浮托力的变化外，还需考虑由于水位下降导致的渗透力作用，渗透力的方向和大小会影响滑坡体的应力分布，进而改变下滑力和抗滑力。降雨工况下，要考虑雨水入渗导致的岩土体重度增加、孔隙水压力升高以及抗剪强度降低等因素。地震工况则需考虑地震力对滑坡体的作用，地震力会增加滑坡体的下滑力。以某一具体工况为例，假设马家沟滑坡在库水位上升至175米时，运用传递系数法计算其稳定性系数。首先，根据滑坡的几何形状和岩土体参数，将滑坡体划分为若干个条块。对于每个条块，计算其自重、浮托力、下滑力和抗滑力。然后，根据传递系数法的计算公式，依次计算每个条块的剩余下滑力，并将其传递到下一个条块。最后，计算整个滑坡体的稳定性系数。经过计算，得到该工况下马家沟滑坡的稳定性系数为1.05。根据相关规范和经验，当稳定性系数大于1.05时，滑坡处于基本稳定状态；当稳定性系数在1.00-1.05之间时，滑坡处于欠稳定状态；当稳定性系数小于1.00时，滑坡处于不稳定状态。由此判断，在该工况下，马家沟滑坡处于基本稳定状态，但仍需密切关注其变形情况，因为稳定性系数接近临界值，一旦外部条件发生变化，如库水位继续上升、遭遇强降雨或地震等，滑坡的稳定性可能会降低，有失稳的风险。极限平衡法在马家沟滑坡稳定性分析中具有重要的应用价值，通过合理选择计算方法和准确获取参数，能够为滑坡的稳定性评价和防治措施的制定提供有力的依据。然而，该方法也存在一定的局限性，如对滑动面形状的假设较为理想化，未充分考虑岩土体的变形和应力-应变关系等。因此，在实际应用中，常结合其他方法，如有限元法等，进行综合分析，以提高分析结果的准确性和可靠性。4.2数值模拟方法4.2.1数值模拟软件介绍本研究选用Geo-Studio软件进行数值模拟分析。Geo-Studio是一套专业、高效且功能强大的适用于岩土工程和岩土环境模拟计算的仿真软件，在岩土工程领域应用广泛。该软件以Geo-SLOPE为主体，集成了多种专业分析软件，可对岩土工程中的渗流、稳定、应力变形、动力响应等多方面问题进行综合分析。其中，SEEP/W模块用于地下水渗流分析，能够全面而简洁地表述多孔渗水材料，如土体和岩石中的地下水渗流和超孔隙水压力消散问题，可分析从简单的饱和稳态问题到复杂的饱和-不饱和时变问题。在研究马家沟滑坡时，可利用该模块模拟库水位涨落过程中滑坡体内的渗流场变化，确定不同时刻孔隙水压力的分布状况。SLOPE/W模块是计算岩土边坡安全系数的主流软件产品，可同时用八种方法分析计算简单或复杂的边坡稳定问题。该模块使用极限平衡理论对不同土体类型、复杂地层和滑移面形状的边坡中的孔隙水压力分布状况进行建模分析，还提供多种土体模型，并支持确定性和随机输入参数方法，可进行随机稳定性分析。针对马家沟滑坡，通过SLOPE/W模块能计算不同工况下的滑坡稳定性系数，分析滑坡的稳定性变化。SIGMA/W模块用于岩土应力变形分析，可模拟岩土体在各种荷载作用下的应力应变状态。在马家沟滑坡的研究中，借助该模块能够分析库水位涨落、降雨等因素作用下，滑坡体内部的应力分布和变形情况。Geo-Studio软件具有多模块集成分析功能，所有模块可在同一界面下运行，用户只需建立一个几何模型，即可在所有分析中使用，便于综合分析岩土工程问题。其项目管理方式允许用户快速创建分析文档，方便进行多工况比选以及耦合分析，创建不同软件之间的耦合分析时，软件会自动调用上级目录的分析结果。软件界面友好，建模步骤和参数设置简单清晰，易于学习和使用，能节省分析时间。后处理功能丰富多样，可将计算结果绘制成多种图形并进行输出整理，还支持一键生成计算报告和动画输出瞬态计算结果等。4.2.2模型建立与参数选取根据马家沟滑坡的实际地形地貌、地层岩性和地质构造等条件，利用Geo-Studio软件建立三维数值模型。首先，通过现场勘查获取滑坡的边界、地形等高线等数据，使用专业的绘图软件绘制滑坡的二维剖面图，然后导入Geo-Studio软件中，利用软件的建模工具将二维剖面图转化为三维几何模型。在建模过程中，精确描绘滑坡体、下伏基岩以及周边地形的几何形状，确保模型能够准确反映实际情况。对于模型的参数选取，主要依据现场勘查、室内试验和工程地质类比法确定。岩土体的物理力学参数，如密度、重度、黏聚力、内摩擦角等，是数值模拟的关键参数。通过现场钻探获取滑坡体和下伏基岩的岩土体样本，在实验室进行物理力学性质试验。例如，采用三轴压缩试验测定岩土体的抗剪强度参数，通过直剪试验确定黏聚力和内摩擦角；利用密度试验和含水率试验测定岩土体的密度和含水率，进而计算出重度。对于一些难以通过试验直接获取的参数，参考工程地质类比法，结合三峡库区其他类似滑坡的研究成果和工程经验进行取值。以第四系全新统滑坡堆积层为例，根据室内试验结果，其密度取值为2.0-2.2g/cm^3，重度取值为20-22kN/m^3，黏聚力取值为15-20kPa，内摩擦角取值为18°-22°。侏罗系上统沙溪庙组砂岩的密度取值为2.5-2.7g/cm^3，重度取值为25-27kN/m^3，黏聚力取值为30-50kPa，内摩擦角取值为30°-35°；泥岩的密度取值为2.3-2.5g/cm^3，重度取值为23-25kN/m^3，黏聚力取值为10-15kPa，内摩擦角取值为15°-20°。在模拟库水位涨落时，根据三峡水库的实际运行数据，设定库水位的初始值、变化范围和升降速率。如初始库水位设定为145米，变化范围为145-175米，上升速率设定为0.3-0.4米/天，下降速率设定为0.3-0.5米/天。同时，考虑降雨因素时，根据马家沟滑坡所在地区的降雨统计资料，设定降雨强度和持续时间等参数。通过合理的模型建立和参数选取，为准确模拟库水位涨落作用下马家沟滑坡的稳定性变化奠定基础。4.2.3模拟结果分析通过Geo-Studio软件对库水位涨落作用下马家沟滑坡进行数值模拟，得到了滑坡体在不同工况下的应力、应变和位移变化规律。在应力变化方面，当库水位上升时，滑坡体前缘由于受到浮托力作用，有效应力减小。模拟结果显示，在库水位从145米上升至175米的过程中，滑坡体前缘部分区域的有效应力降低了10%-20%。而在滑坡体后部，由于库水的渗入，孔隙水压力升高，导致有效应力也有所减小。当库水位下降时，滑坡体内部产生渗流，形成渗透力，使得滑坡体后部的有效应力进一步减小，而前缘部分由于渗透力的作用，有效应力略有增加。在库水位快速下降阶段，滑坡体后部有效应力可降低20%-30%。这种应力的变化导致滑坡体内部的应力分布发生改变，在一些薄弱部位容易产生应力集中现象。在滑坡体与下伏基岩的接触面处，以及岩土体性质差异较大的部位，应力集中较为明显，这些部位的应力值可比周围区域高出15%-25%，增加了滑坡失稳的风险。从应变变化来看，库水位涨落过程中，滑坡体的应变主要集中在滑坡体的前缘和后缘。库水位上升时，前缘土体由于受到浮托力和库水压力的作用，产生向临空面的挤压变形，导致应变增大。模拟结果表明，库水位上升阶段，滑坡体前缘的水平应变可达到1.5×10^{-3}-2.5×10^{-3}。库水位下降时，后缘土体由于渗透力的作用，产生拉伸变形，应变也相应增大。在库水位快速下降时，滑坡体后缘的垂直应变可达到1.0×10^{-3}-2.0×10^{-3}。随着库水位的反复涨落，滑坡体的应变不断累积，当应变超过岩土体的极限应变时，岩土体就会发生破坏，从而引发滑坡。关于位移变化，库水位涨落引起滑坡体的位移主要表现为水平位移和垂直位移。库水位上升时，滑坡体前缘产生向库水方向的水平位移，垂直位移则表现为下沉。在库水位从145米上升至175米的过程中，滑坡体前缘的最大水平位移可达5-10厘米，垂直位移可达3-5厘米。库水位下降时，滑坡体后缘产生向坡外的水平位移，垂直位移表现为抬升。库水位快速下降时，滑坡体后缘的最大水平位移可达8-12厘米，垂直位移可达5-8厘米。位移的变化反映了滑坡体的变形趋势，当位移达到一定程度时，滑坡就可能发生滑动。通过对位移变化规律的分析，可以提前预测滑坡的发生，为滑坡防治提供依据。综合应力、应变和位移的模拟结果可知，库水位涨落对马家沟滑坡的稳定性产生了显著影响。在库水位上升和下降过程中，滑坡体的应力、应变和位移变化相互作用，导致滑坡体的稳定性逐渐降低。尤其是在库水位快速下降阶段，滑坡体内部的应力、应变和位移变化最为剧烈，滑坡失稳的风险最大。因此，在滑坡防治中，应重点关注库水位快速下降阶段，采取有效的防治措施，降低滑坡失稳的风险。4.3现场监测方法4.3.1监测内容与方法为全面掌握马家沟滑坡在库水位涨落作用下的稳定性变化情况，对其进行了系统的现场监测，监测内容涵盖多个关键方面，采用了多种先进且有效的监测方法。在位移监测方面，地表位移监测是重要环节。运用全站仪对滑坡体上预先设置的多个地表监测点进行定期观测，全站仪利用电磁波测距原理，通过测量仪器到监测点的距离以及水平角和垂直角，精确计算出监测点的三维坐标。在每次观测时，确保仪器的对中、整平精度，严格按照测量规范进行操作，以获取准确的坐标数据。通过对比不同时期的坐标数据，可计算出监测点在水平和垂直方向上的位移量，从而掌握滑坡体的地表位移变化情况。为提高监测精度，还结合了GPS监测技术，GPS接收机通过接收卫星信号，利用卫星定位原理确定监测点的位置。GPS具有全天候、高精度、实时性强等优点，能够提供更为准确的三维坐标信息。在滑坡体上均匀布置多个GPS监测点，与全站仪监测点相互补充，共同监测地表位移。定期对全站仪和GPS监测数据进行对比分析，确保数据的可靠性和一致性。对于滑坡体的深部位移监测，采用钻孔倾斜仪进行监测。在滑坡体上钻孔，将钻孔倾斜仪放入钻孔中，通过测量钻孔倾斜仪与铅垂线的夹角变化，来确定滑坡体深部不同深度处的水平位移情况。钻孔倾斜仪一般由测斜管、探头和读数仪组成，测斜管安装在钻孔中，探头在测斜管内移动进行测量。在安装测斜管时，严格控制其垂直度和安装深度，确保测量数据的准确性。定期对钻孔倾斜仪进行校准和维护，保证其正常工作。地下水位监测对于分析滑坡稳定性至关重要。在滑坡体上布置多个地下水位监测井，监测井采用PVC管制作，管壁设有滤水孔，以保证地下水能够顺利进入监测井。通过水位计测量监测井中的水位高度，水位计可采用压力式水位计或投入式水位计，压力式水位计通过测量地下水的压力来计算水位高度，投入式水位计则直接测量水位深度。定期记录水位数据，分析地下水位在库水位涨落过程中的变化规律。孔隙水压力监测也是重要的监测内容。在滑坡体不同位置埋设孔隙水压力计，孔隙水压力计一般采用振弦式或电阻应变式。振弦式孔隙水压力计通过测量钢弦的振动频率来确定孔隙水压力，电阻应变式孔隙水压力计则利用电阻应变片的变形来测量孔隙水压力。孔隙水压力计的埋设位置根据滑坡体的地质条件和渗流场特征进行合理选择，确保能够准确反映孔隙水压力的变化情况。定期读取孔隙水压力计的数据，分析孔隙水压力在库水位涨落过程中的变化趋势及其对滑坡稳定性的影响。在监测过程中，还运用了地面摄影测量技术。通过在滑坡体周边合适位置设置摄影站点，使用专业的测量相机对滑坡体进行定期拍摄。利用摄影测量软件对拍摄的照片进行处理，通过图像匹配和三维重建技术，获取滑坡体表面的三维模型和变形信息。地面摄影测量技术能够快速、全面地获取滑坡体的表面变形情况，且具有非接触、效率高的优点。将地面摄影测量结果与全站仪、GPS等监测数据进行对比分析，相互验证和补充，提高监测结果的准确性和可靠性。4.3.2监测数据处理与分析对现场监测获取的数据进行了严谨的处理和深入的分析，以验证数值模拟结果的准确性，全面掌握马家沟滑坡在库水位涨落作用下的稳定性变化规律。在数据处理阶段，首先对原始监测数据进行了质量检查，剔除明显异常的数据点。对于位移监测数据，若某个监测点的位移量突然出现大幅跳跃或与相邻监测点的变化趋势明显不符，且经检查仪器设备和测量过程无故障，可判断该数据为异常数据，予以剔除。对于地下水位和孔隙水压力监测数据，若出现超出合理范围的值，也进行仔细排查和甄别，确定其是否为异常数据。采用滤波算法对数据进行平滑处理，去除数据中的噪声干扰。常用的滤波算法有滑动平均滤波、中值滤波等。以位移监测数据为例，使用滑动平均滤波时，设定一个合适的窗口大小，如5个监测周期，对窗口内的数据进行平均计算，得到平滑后的位移值。通过滤波处理，使监测数据更加平稳，能够更清晰地反映滑坡体的真实变形趋势。对处理后的数据进行统计分析，计算位移、地下水位、孔隙水压力等参数的均值、方差、最大值、最小值等统计量。统计不同监测时段内位移的均值，可了解滑坡体在该时段内的平均变形程度；计算方差则能反映位移数据的离散程度，方差越大，说明位移变化越不稳定。通过这些统计量，对监测数据的整体特征有更直观的认识。将监测数据与数值模拟结果进行对比分析。在库水位上升阶段，对比监测得到的滑坡体前缘位移与数值模拟中该部位的位移结果。若监测位移值为5-8厘米，而数值模拟结果为6-9厘米，两者在数量级和变化趋势上基本一致，说明数值模拟能够较好地反映库水位上升时滑坡体前缘的位移变化情况。在库水位下降阶段，对比监测的地下水位和孔隙水压力与数值模拟结果。若监测显示地下水位在库水位下降后的2-3天内下降了1-2米，孔隙水压力相应降低，数值模拟也得到类似的变化趋势和数值范围，表明数值模拟在反映库水位下降时的渗流场变化方面具有较高的准确性。针对两者的差异进行深入分析。若在某些工况下，监测数据与数值模拟结果存在较大偏差，如在某次强降雨与库水位下降耦合作用的工况下，监测到的滑坡体后缘位移明显大于数值模拟结果。此时，从多个方面查找原因，检查数值模拟中降雨入渗参数的设置是否合理，是否充分考虑了强降雨对岩土体力学性质的影响；分析监测过程中是否存在外界干扰因素影响了监测数据的准确性；研究滑坡体实际地质条件与数值模拟模型中假设条件的差异，如是否存在未被发现的软弱结构面等。通过对差异的分析，不断改进数值模拟模型和参数设置，提高数值模拟的精度和可靠性。通过对监测数据的处理和分析，验证了数值模拟结果在一定程度上能够准确反映库水位涨落作用下马家沟滑坡的稳定性变化规律，但也存在需要改进和完善的地方。监测数据与数值模拟结果的相互验证和补充，为深入研究滑坡稳定性提供了有力支持，也为滑坡防治措施的制定提供了更可靠的依据。五、马家沟滑坡稳定性评价与预测5.1稳定性评价指标与标准马家沟滑坡稳定性评价采用稳定性系数和位移速率作为关键指标，各指标依据相关规范和研究成果确定评判标准。稳定性系数是衡量滑坡稳定性的重要量化指标，通常通过极限平衡法、有限元强度折减法等计算得出。在极限平衡法中，以抗滑力与下滑力的比值作为稳定性系数，当稳定性系数大于1时，表明抗滑力大于下滑力，滑坡处于稳定状态；等于1时，抗滑力与下滑力相等，滑坡处于极限平衡状态；小于1时，抗滑力小于下滑力，滑坡处于不稳定状态。有限元强度折减法通过不断折减岩土体的抗剪强度参数，直至滑坡体达到极限平衡状态，此时的折减系数即为稳定性系数。依据《建筑边坡工程技术规范》（GB50330-2013），对于永久性边坡，安全等级为一级时，稳定性系数不应小于1.35；安全等级为二级时，不应小于1.30；安全等级为三级时，不应小于1.25。马家沟滑坡由于其潜在的危害性较大，涉及居民点及重要基础设施，将其安全等级定为一级，因此稳定性系数应不小于1.35。若计算得到的稳定性系数在1.35及以上，表明滑坡处于稳定状态，发生滑动的可能性较小；在1.25-1.35之间，滑坡处于欠稳定状态，需密切关注，采取一定的监测和防范措施；小于1.25时，滑坡处于不稳定状态，随时可能发生滑动，必须立即采取有效的防治措施。位移速率反映了滑坡体的变形速度，是判断滑坡稳定性变化趋势的重要依据。通过现场监测获取滑坡体的位移数据，计算单位时间内的位移量，即得到位移速率。一般来说，位移速率越大，表明滑坡体的变形越剧烈，稳定性越差。根据相关研究和工程经验，对于马家沟滑坡，当位移速率小于5mm/月时，滑坡处于相对稳定状态。此时，滑坡体的变形较为缓慢，短期内发生滑动的可能性较低，但仍需定期进行监测，掌握其变形动态。位移速率在5-10mm/月之间时，滑坡处于基本稳定状态，但需加强监测频率。这表明滑坡体的变形有逐渐增大的趋势，可能受到库水位涨落、降雨等因素的影响，需要密切关注位移速率的变化情况，及时分析原因，采取相应的措施。若位移速率大于10mm/月，滑坡处于不稳定状态，应立即发出预警，组织人员撤离危险区域，并采取紧急防治措施。此时，滑坡体的变形急剧增加，可能即将发生滑动，必须高度重视，确保人民生命财产安全。在实际评价马家沟滑坡稳定性时，综合考虑稳定性系数和位移速率这两个指标。当稳定性系数满足安全要求且位移速率在相对稳定范围内时，可认为滑坡处于稳定状态；若其中一个指标不满足要求，就需要对滑坡的稳定性进行进一步分析和评估，结合其他监测数据和影响因素，判断滑坡的稳定性状况，并制定相应的防治策略。5.2不同工况下的稳定性评价5.2.1正常蓄水位工况在正常蓄水位175米工况下，运用极限平衡法中的传递系数法和有限元强度折减法，结合Geo-Studio软件进行数值模拟，对马家沟滑坡的稳定性展开深入分析。通过现场勘查、室内试验获取的岩土体物理力学参数，如第四系全新统滑坡堆积层的密度、重度、黏聚力、内摩擦角，以及侏罗系上统沙溪庙组砂岩和泥岩的相应参数，作为计算和模拟的基础数据。采用传递系数法计算时，将滑坡体划分为多个条块，依次计算每个条块的下滑力和抗滑力，并考虑条块间的作用力传递。计算结果表明，在正常蓄水位175米工况下，马家沟滑坡的稳定性系数为1.28。此稳定性系数略低于一级边坡稳定性系数不应小于1.35的标准，表明在该工况下，滑坡处于欠稳定状态。从计算过程可知，滑坡体前缘由于受到较大的浮托力作用，有效重量减轻，抗滑力降低，对稳定性产生了一定的不利影响。同时，滑坡体后部因库水的渗入，孔隙水压力升高，有效应力减小，也在一定程度上降低了滑坡的稳定性。运用有限元强度折减法进行分析时，通过不断折减岩土体的抗剪强度参数，直至滑坡体达到极限平衡状态。数值模拟结果显示，在正常蓄水位工况下，滑坡体的位移主要集中在滑坡体的前缘和后缘。前缘出现向库水方向的水平位移，最大位移量达到8厘米左右；后缘则产生向坡外的水平位移，最大位移量约为6厘米。滑坡体内部的应力分布也发生了明显变化，在滑坡体与下伏基岩的接触面处，以及岩土体性质差异较大的部位，出现了应力集中现象，应力集中区域的应力值比周围区域高出20%左右。综合位移和应力分析结果，该工况下的滑坡稳定性系数为1.26。这一结果与传递系数法计算结果相近，进一步验证了在正常蓄水位175米工况下，马家沟滑坡处于欠稳定状态。通过对正常蓄水位工况下的稳定性分析可知，虽然目前滑坡尚未处于不稳定状态，但稳定性系数接近临界值，且滑坡体已经出现了一定的位移和应力集中现象。若外部条件发生变化，如遭遇强降雨、地震等，滑坡的稳定性可能会进一步降低，有失稳的风险。因此，在该工况下，需密切关注滑坡的变形情况，加强监测和预警工作，提前制定应对措施，以保障滑坡体周边居民的生命财产安全和基础设施的稳定运行。5.2.2水位骤涨骤落工况在水位骤涨工况下，假设库水位在短时间内从145米快速上升至175米，上升速率达到0.5米/天。利用Geo-Studio软件中的SEEP/W模块和SLOPE/W模块进行数值模拟分析。在SEEP/W模块中，模拟库水位快速上升过程中滑坡体内的渗流场变化。结果显示，库水位快速上升时，滑坡体前缘的地下水迅速上升，在1-2天内，前缘部分区域的地下水位就可上升至与库水位接近的位置。由于水位上升速度快，滑坡体内部形成了较大的水头差，导致孔隙水压力在短时间内急剧升高。在滑坡体前缘，孔隙水压力升高幅度可达15-20kPa。通过SLOPE/W模块计算滑坡的稳定性系数。考虑到孔隙水压力升高导致有效应力减小，抗剪强度降低，以及浮托力对滑坡体有效重量的影响。计算结果表明，在水位骤涨工况下，马家沟滑坡的稳定性系数降至1.18。此稳定性系数明显低于一级边坡的安全标准，表明滑坡处于不稳定状态。从稳定性系数的降低幅度来看，与正常蓄水位工况相比，稳定性系数下降了约8%。这说明水位骤涨对滑坡稳定性的影响较为显著，在水位快速上升过程中，孔隙水压力和浮托力的共同作用，极大地降低了滑坡的稳定性。对于水位骤落工况，假设库水位在短时间内从175米快速下降至145米，下降速率为0.5米/天。同样利用Geo-Studio软件进行模拟。在SEEP/W模块中，模拟结果显示，库水位快速下降时，滑坡体内部的地下水排泄不畅，形成了较大的水头差，产生了较强的渗透力。在滑坡体后部，由于地下水向库水方向渗流，形成向上的渗透力，使岩土体的有效重量减轻。在滑坡体前缘，渗透力方向向下，增加了下滑力。在滑坡体后部，渗透力可使岩土体的有效重量减轻10%-15%；在滑坡体前缘，下滑力可增加15%-20%。通过SLOPE/W模块计算滑坡稳定性系数。考虑到渗透力对滑坡体应力分布和抗滑力、下滑力的影响。计算结果表明，在水位骤落工况下，马家沟滑坡的稳定性系数降至1.15。此稳定性系数进一步降低，说明水位骤落对滑坡稳定性的影响更为严重，滑坡处于不稳定状态。与水位骤涨工况相比，水位骤落时稳定性系数下降幅度更大，这是因为水位骤落产生的渗透力对滑坡体的应力状态改变更为显著，增加了滑坡失稳的风险。综上所述，水位骤涨骤落工况对马家沟滑坡的稳定性影响显著，在这两种工况下，滑坡均处于不稳定状态。在三峡水库运行过程中，应尽量避免库水位的骤涨骤落，合理控制水位变化速率，以减少对滑坡稳定性的不利影响。同时，针对水位骤涨骤落可能引发的滑坡失稳风险，需加强监测和预警，制定应急预案，确保在滑坡失稳时能够及时采取有效的应对措施，保障人民生命财产安全。5.2.3降雨工况在考虑降雨工况时，假设马家沟滑坡遭遇一场持续时间为3天，日降雨量达到50mm的强降雨。降雨通过入渗增加滑坡体的含水量，使岩土体重度增大，下滑力增加。雨水入渗还会使孔隙水压力升高，有效应力减小，抗剪强度降低。利用Geo-Studio软件进行数值模拟分析，首先在SEEP/W模块中，根据降雨强度和持续时间，设定入渗边界条件，模拟降雨过程中滑坡体内的渗流场变化。结果显示，降雨开始后，滑坡体上部的含水量迅速增加，在1-2天内，上部部分区域的饱和度可达到80%-90%。随着降雨的持续，水分逐渐向下渗透，地下水位上升。在降雨结束时，地下水位上升幅度可达2-3米。通过SLOPE/W模块计算滑坡的稳定性系数。考虑到降雨导致的岩土体重度增加、孔隙水压力升高以及抗剪强度降低等因素。计算结果表明，在降雨工况下，马家沟滑坡的稳定性系数降至1.20。此稳定性系数低于一级边坡的安全标准，表明滑坡处于不稳定状态。从稳定性系数的变化来看，与正常工况相比，稳定性系数下降了约6%。这说明降雨对滑坡稳定性有明显的不利影响，强降雨使得滑坡体的力学条件恶化，增加了滑坡失稳的风险。在实际情况中，降雨往往与库水位涨落相互作用，进一步影响滑坡的稳定性。若降雨发生在库水位上升阶段，滑坡体饱水程度加剧，地下水水位大幅上升，岩土体抗剪强度进一步降低，滑坡稳定性受到更大威胁。当降雨与库水位上升同时发生时，稳定性系数可比单纯降雨工况下再降低5%-8%。若降雨发生在库水位下降阶段，由于库水位下降导致的渗透力与降雨入渗产生的孔隙水压力变化相互叠加，会进一步改变滑坡体的应力状态，增加滑坡失稳的可能性。在库水位下降过程中遭遇降雨，稳定性系数可比单纯降雨工况下降低8%-10%。降雨工况对马家沟滑坡的稳定性产生了显著的不利影响，在降雨尤其是强降雨条件下，滑坡处于不稳定状态。在三峡库区的日常监测和管理中，应密切关注降雨情况，加强对滑坡体的监测频率和力度。在降雨来临前，提前做好预警工作，采取相应的防范措施，如加强排水系统的维护和疏通，确保降雨期间滑坡体的排水顺畅，降低滑坡失稳的风险。同时，在进行滑坡稳定性分析和防治措施制定时，应充分考虑降雨与库水位涨落的耦合作用，提高防治措施的针对性和有效性。5.2.4地震工况在考虑地震工况时，假设马家沟滑坡遭遇里氏5.0级地震。地震产生的地震力会使滑坡体产生惯性力，增加下滑力。同时，地震还会使岩土体结构受到破坏，孔隙率增大，抗剪强度降低。利用有限元软件ABAQUS建立马家沟滑坡的三维数值模型，考虑地震作用下的动力响应。在模型中，输入地震波加速度时程曲线，模拟地震过程中滑坡体的应力、应变和位移变化。模拟结果显示，在地震作用下，滑坡体的加速度响应在滑坡体的顶部和前缘较为明显，顶部的加速度峰值可达0.2g-0.3g，前缘的加速度峰值可达0.15g-0.25g。由于地震力的作用，滑坡体产生了较大的惯性力，导致下滑力增加。在滑坡体的后缘，惯性力使岩土体产生拉应力，当拉应力超过岩土体的抗拉强度时，会出现拉裂破坏。在滑坡体的前缘，惯性力使岩土体产生向临空面的挤压变形，导致应变增大。通过对地震作用下的应力、应变和位移分析，计算滑坡的稳定性系数。考虑到地震导致的下滑力增加、岩土体抗剪强度降低等因素。计算结果表明，在地震工况下，马家沟滑坡的稳定性系数降至1.10。此稳定性系数远低于一级边坡的安全标准，表明滑坡处于不稳定状态。与正常工况相比，稳定性系数下降了约14%。这说明地震对滑坡稳定性的影响非常显著，强烈的地震作用使得滑坡体的力学条件急剧恶化，极大地增加了滑坡失稳的风险。地震与库水位涨落共同作用时，对马家沟滑坡稳定性的影响更为复杂。当库水位处于较高位置时，地震作用下饱和岩土体的抗剪强度降低更为明显，滑坡更容易失稳。在库水位较高且遭遇地震时，稳定性系数可比单纯地震工况下再降低8%-10%。地震还会使库水位涨落引起的渗流场变化加剧，进一步影响滑坡的稳定性。地震使岩土体的渗透性改变，导致地下水渗流更加复杂。库水位上升时，地震作用可能使岩土体中的孔隙和裂隙扩张，加速库水的渗入，进一步提高孔隙水压力；库水位下降时，地震造成的岩土体结构破坏会使地下水排泄不畅，增加滑坡体内部的水头差，增大渗透力。这种渗流场的变化与地震力的叠加，极大地增加了滑坡失稳的风险。地震工况对马家沟滑坡的稳定性产生了严重的不利影响，在地震尤其是与库水位涨落共同作用时，滑坡处于高度不稳定状态。在三峡库区的地震监测和滑坡防治工作中，应加强地震监测和预警，提高对地震灾害的防范能力。对于马家沟滑坡，应制定针对地震工况的应急预案，在地震发生时，能够迅速采取有效的应对措施，如组织人员撤离、进行紧急抢险救援等，以减少地震引发的滑坡灾害损失。同时，在进行滑坡稳定性分析和防治工程设计时，应充分考虑地震与库水位涨落的耦合作用，提高防治工程的抗震能力和可靠性。5.3滑坡稳定性预测基于对马家沟滑坡历史数据的深入分析，以及不同工况下稳定性评价结果，运用时间序列分析和灰色预测等方法，对其在未来库水位涨落条件下的稳定性变化趋势进行预测。时间序列分析方法是通过对过去的稳定性系数、位移速率等数据进行建模，来预测未来的变化趋势。以稳定性系数为例，收集了过去10年马家沟滑坡在不同库水位工况下的稳定性系数数据，运用ARIMA（自回归积分滑动平均）模型进行分析。首先对稳定性系数数据进行平稳性检验，通过ADF检验发现原始数据不平稳，对其进行一阶差分后，数据达到平稳状态。然后利用自相关函数（ACF）和偏自相关函数（PACF）确定模型的参数，经过多次试验和比较，确定ARIMA(1,1,1)模型为最优模型。利用该模型对未来5年马家沟滑坡在正常库水位涨落工况下的稳定性系数进行预测。预测结果显示，在未来5年内，若库水位按照正常的涨落规律变化，马家沟滑坡的稳定性系数将呈现缓慢下降的趋势。在第1年，稳定性系数预计为1.26，到第5年，稳定性系数可能降至1.22左右。这表明随着时间的推移，在正常库水位涨落作用下，马家沟滑坡的稳定性逐渐降低，虽然下降幅度相对较小，但仍需持续关注。灰色预测方法则是基于灰色系统理论，对含有不确定因素的系统进行预测。对于马家沟滑坡稳定性预测，选取了滑坡体的位移数据作为预测指标。收集了过去5年的滑坡体月位移数据，建立GM(1,1)灰色预测模型。首先对原始位移数据进行累加生成处理，以弱化数据的随机性，得到累加生成序列。然后根据累加生成序列建立微分方程，通过最小二乘法估计方程中的参数，得到灰色预测模型。利用该模型对未来12个月的滑坡体位移进行预测。预测结果表明，未来12个月内，滑坡体的位移将呈逐渐增加的趋势。前3个月，位移增加较为缓慢，月位移增量约为3-5mm；随着时间的推移，从第4个月开始，位移增量逐渐增大，到第12个月，月位移增量预计达到8-10mm。这说明马家沟滑坡体的变形在未来一段时间内将持续发展，稳定性逐渐变差。综合时间序列分析和灰色预测结果，在未来库水位涨落条件下，马家沟滑坡的稳定性总体呈下降趋势。稳定性系数的降低和位移的增加，都表明滑坡体的稳定性逐渐变差，失稳的风险逐渐增大。尤其是在遭遇库水位骤涨骤落、强降雨、地震等不利工况时，滑坡的稳定性可能会急剧下降，随时有失稳滑动的危险。因此，根据预测结果，必须加强对马家沟滑坡的长期监测，建立完善的监测预警系统，及时掌握滑坡体的变形动态和稳定性变化情况。同时，提前制定有效的防治措施，如加强排水工程建设，降低地下水位，减小孔隙水压力；在滑坡体前缘设置支挡结构，增加抗滑力；对滑坡体后部进行减重处理，减小下滑力等。通过这些措施，降低滑坡失稳的风险，保障滑坡体周边居民的生命财产安全和基础设施的稳定运行。六、滑坡防治措施与建议6.1防治原则与目标马家沟滑坡防治遵循“以人为本、预防为主、综合治理、安全经济”的原则。“以人为本”是将保障人民生命财产安全置于首位，这是滑坡防治的核心宗旨。在制定和实施防治措施过程中，充分考虑滑坡体周边居民的安全，避免因滑坡失稳造成人员伤亡和重大财产损失。“预防为主”强调在灾害发生前采取积极有效的措施，降低滑坡发生的风险。通过加强监测预警，及时掌握滑坡体的变形动态和稳定性变化情况，提前发现潜在的滑坡隐患，为采取防治措施提供依据。同时，采取工程措施对滑坡体进行加固和防护，从源头上减少滑坡发生的可能性。“综合治理”要求综合运用多种防治手段，针对马家沟滑坡的特点和影响因素，全面考虑工程地质条件、水文地质条件、库水位涨落规律等，将工程措施与非工程措施相结合。工程措施包括排水工程、支挡工程、减重与反压等，通过这些措施改善滑坡体的力学条件，提高其稳定性。非工程措施如加强监测预警、制定应急预案、开展地质灾害科普宣传等，提高对滑坡灾害的应对能力和公众的防灾意识。“安全经济”原则是在确保防治措施有效性和安全性的前提下，合理控制防治成本，提高资金使用效率。对不同的防治方案进行技术经济比较，选择既能有效防治滑坡灾害，又具有良好经济效益的方案。在工程设计和施工过程中，优化设计方案，合理选用材料和施工工艺，避免不必要的浪费。马家沟滑坡防治的总体目标是确保滑坡体在各种工况下的稳定性，保障滑坡体周边居民的生命财产安全和基础设施的稳定运行。具体目标包括将滑坡体的稳定性系数提高到安全标准以上，根据《建筑边坡工程技术规范》（GB50330-2013），对于马家沟滑坡这种安全等级为一级的边坡，稳定性系数应不小于1.35。通过采取有效的防治措施，使滑坡体在正常库水位涨落、降雨、地震等工况下，稳定性系数均能满足安全要求。将滑坡体的位移控制在允许范围内，通过现场监测和数值模拟分析，确定滑坡体在不同工况下的位移变化规律，设定合理的位移控制指标。采取防治措施后，使滑坡体的地表位移和深部位移均能控制在允许范围内，确保滑坡体不会发生明显的变形和滑动。建立完善的监测预警系统也是防治目标之一，实时掌握滑坡体的变形动态和稳定性变化情况。通过布置多种监测设备，如全站仪、GPS、钻孔倾斜仪、孔隙水压力计等，对滑坡体的位移、地下水位、孔隙水压力等参数进行实时监测。利用监测数据，及时分析滑坡体的稳定性变化趋势，一旦发现异常情况，能够迅速发出预警信号，为人员撤离和采取应急措施提供充足的时间。制定科学合理的应急预案，提高应对滑坡灾害的能力。应急预案应包括灾害发生后的应急响应程序、人员疏散方案、抢险救援措施、物资调配等内容。定期对应急预案进行演练和修订，确保在灾害发生时能够迅速、有效地实施，最大限度地减少灾害损失。6.2防治措施针对马家沟滑坡的实际情况，采取工程措施与非工程措施相结合的方式进行防治。工程措施主要包括排水工程、支挡工程和减重与反压工程，旨在改善滑坡体的力学条件，提高其稳定性；非工程措施则主要是监测预警，通过实时监测滑坡体的变形动态，及时发出预警信号，为采取应急措施提供依据。排水工程对于降低马家沟滑坡地下水位、减小孔隙水压力和渗透力至关重要。在滑坡体表面，布置截水沟和排水沟，截水沟设置在滑坡体后缘及周边，拦截地表水，使其不流入滑坡体。排水沟则沿滑坡体坡面布置，将地表水迅速排出。截水沟和排水沟采用浆砌石或混凝土结构，沟底和沟壁应平整光滑，以保证排水顺畅。沟的尺寸根据汇水面积和降雨量计算确定，一般截水沟的断面尺寸为宽0.8-1.0米，深0.6-0.8米；排水沟的断面尺寸为宽0.5-0.8米，深0.4-0.6米。为防止沟壁坍塌和水流冲刷，沟壁设置伸缩缝和泄水孔，伸缩缝间距一般为10-15米，泄水孔采用直径5-10厘米的PVC管，间距1-2米，呈梅花形布置。对于滑坡体内部的地下水，采用仰斜排水孔进行疏干。在滑坡体中上部钻孔，钻孔方向与坡面成一定仰角，一般为10°-15°，使地下水能够顺利流入排水孔并排出。排水孔采用PVC管或钢管，管径一般为5-10厘米，长度根据滑坡体厚度和地下水位情况确定，一般为10-20米。为保证排水效果，排水孔内填充砾石或粗砂等滤水材料，防止孔壁坍塌和堵塞。在滑坡体中下部设置排水盲沟，盲沟采用片石或碎石填筑，外包土工布，以拦截和排除地下水。盲沟的断面尺寸一般为宽0.8-1.0米，深1.0-1.5米，长度根据实际情况确定。通过这些排水措施，可有效降低滑坡体的地下水位，减小孔隙水压力和渗透力，提高滑坡体的稳定性。支挡工程是提高马家沟滑坡稳定性的重要手段，主要通过在滑坡体前缘设置抗滑桩和挡土墙来实现。抗滑桩是一种常用的支挡结构，在马家沟滑坡前缘，根据滑坡体的规模和稳定性要求，合理布置抗滑桩。抗滑桩采用钢筋混凝土结构，桩径一般为1.2-1.5米，桩长根据滑坡体厚度和滑动面位置确定，一般为15-25米。桩间距根据滑坡推力和桩的承载能力计算确定，一般为3-5米。抗滑桩的施工采用人工挖孔或机械钻孔的方式，施工过程中要注意保证桩的垂直度和混凝土浇筑质量。在桩身达到设计强度后，将桩与滑坡体紧密连接，通过桩的抗滑力来阻止滑坡体的滑动。挡土墙设置在滑坡体前缘坡脚处，采用重力式挡土墙或悬臂式挡土墙。重力式挡土墙依靠自身重力来维持稳定，适用于小型滑坡或滑坡推力较小的情况。其墙身采用浆砌石或混凝土砌筑，墙高一般为3-5米，墙顶宽度不小于0.5米，