降雨条件对吾隘滑坡稳定性的多维度解析与评估

降雨条件对吾隘滑坡稳定性的多维度解析与评估一、引言1.1研究背景与意义滑坡作为一种极具破坏力的地质灾害，给人类的生命财产安全和社会经济发展带来了巨大威胁。据统计，我国受滑坡灾害威胁的区域占全国陆地面积的1/5-1/4，西南地区发生的滑坡约占全国滑坡次数的一半以上。滑坡不仅会直接摧毁房屋、道路、桥梁等基础设施，导致交通中断、江河堵塞，还可能引发次生灾害，如泥石流、堰塞湖等，进一步加剧灾害的影响范围和破坏程度。在城镇地区，滑坡可能砸埋房屋，造成人员伤亡和财产损失，导致工厂、学校、机关单位等无法正常运转；在乡村，滑坡会毁坏农田、房舍，伤害人畜，破坏农业生产和生态环境。降雨是诱发滑坡的主要因素之一，约90%的滑坡由暴雨引发。河池市南丹县吾隘镇地处山区，地形地貌复杂，地质条件较为脆弱。吾隘滑坡的稳定性对当地居民的生命财产安全以及区域的可持续发展至关重要。降雨通过多种机制影响吾隘滑坡的稳定性。一方面，降雨入渗会使滑坡体的含水量增加，导致岩土体容重增大，从而增加下滑力。另一方面，雨水的入渗会使孔隙水压力升高，有效应力降低，削弱岩土体的抗剪强度，降低滑坡的稳定性。此外，长期的降雨还可能导致岩土体的软化、泥化，进一步破坏滑坡体的结构，增加滑坡发生的风险。研究不同降雨条件下吾隘滑坡的稳定性具有重要的现实意义。从防灾减灾角度来看，通过深入分析降雨与滑坡稳定性之间的关系，可以建立更加准确的滑坡预警模型，提前预测滑坡的发生，及时采取有效的防范措施，如人员撤离、工程加固等，从而最大限度地减少滑坡灾害造成的损失。从地质研究角度而言，该研究有助于深化对滑坡形成机理和演化过程的认识，丰富和完善滑坡稳定性理论，为类似地质条件下的滑坡研究提供参考和借鉴。1.2国内外研究现状在降雨入渗研究方面，国外起步较早，20世纪中叶，Richards建立了非饱和土的渗流理论，为降雨入渗研究奠定了理论基础。此后，众多学者在此基础上进行拓展和深化。Philip通过理论分析和实验研究，提出了Philip入渗公式，定量描述了降雨入渗过程中水分的运动规律。在数值模拟方面，有限差分法、有限元法等数值方法被广泛应用于降雨入渗过程的模拟。如Sutcliffe等利用有限差分法对非饱和土壤中的水分运动进行了数值模拟，取得了较好的效果。国内学者在降雨入渗研究方面也取得了丰硕成果。张有芷等通过室内试验，研究了不同土壤质地条件下降雨入渗的规律，分析了土壤孔隙结构、初始含水量等因素对降雨入渗的影响。近年来，随着计算机技术的发展，国内学者利用数值模拟软件，如GeoStudio、FLAC等，对复杂地质条件下的降雨入渗过程进行模拟，取得了许多有价值的成果。然而，目前降雨入渗研究仍存在一些不足，如对复杂地质条件下的降雨入渗机理认识不够深入，数值模拟中参数的选取缺乏足够的理论依据等。在滑坡稳定性分析方面，国外在理论研究和工程应用方面都取得了显著进展。极限平衡法是最早应用且目前仍广泛使用的方法，如Fellenius法、Bishop法等，这些方法通过对滑坡体进行力学分析，计算滑坡的稳定系数，评估滑坡的稳定性。随着计算机技术和数值计算方法的发展，数值分析法如有限元法、离散元法等逐渐应用于滑坡稳定性分析。如Zienkiewicz等将有限元法引入岩土工程领域，为滑坡稳定性分析提供了新的手段。国内学者在滑坡稳定性分析方面也进行了大量研究。孙广忠提出了岩体结构控制论，为滑坡稳定性分析提供了新的理论框架。许强等通过对滑坡变形破坏机制的研究，建立了基于变形监测的滑坡稳定性评价方法。尽管滑坡稳定性分析取得了很大进展，但仍存在一些问题，如对于复杂地质条件下的滑坡，各种分析方法的计算结果存在一定差异，缺乏统一的评价标准；对滑坡的长期稳定性研究相对较少等。在降雨入渗条件下滑坡稳定性研究方面，国外学者开展了大量的研究工作。Crozier通过对新西兰地区滑坡的研究，分析了降雨强度、持续时间与滑坡发生之间的关系，提出了降雨诱发滑坡的临界雨量阈值。在数值模拟方面，不少学者利用耦合模型对降雨入渗条件下的滑坡稳定性进行分析。如Wu等建立了考虑降雨入渗、地下水渗流和土体力学响应的耦合模型，对滑坡稳定性进行了数值模拟。国内学者在该领域也取得了重要成果。殷坤龙等通过对三峡库区滑坡的研究，分析了降雨入渗对滑坡稳定性的影响机制，建立了基于GIS的降雨型滑坡预警模型。叶俊以永泰县口乡后亭村滑坡体为例，采用GeoStudio软件建立滑坡体饱和－非饱和渗流数值模型，分析了降雨过程中滑坡体孔隙水压力等值线和稳定系数的变化。然而，目前该领域仍存在一些亟待解决的问题，如降雨入渗与滑坡稳定性之间的耦合关系还不够明确，缺乏考虑多因素耦合作用的精细化模型；对降雨型滑坡的预测精度有待进一步提高等。1.3研究内容与技术路线本研究以河池市南丹县吾隘滑坡为研究对象，深入探究不同降雨条件下其稳定性变化规律。具体研究内容如下：吾隘滑坡地质条件勘察：对吾隘滑坡进行详细的野外地质勘察，包括滑坡的地形地貌、地层岩性、地质构造、水文地质条件等。通过现场调查、测绘、钻探等手段，获取滑坡的基本地质信息，为后续的稳定性分析提供基础资料。例如，详细记录滑坡体的边界范围、坡度、高差等地形数据，分析地层的岩性组合、厚度及分布特征，确定地质构造对滑坡的影响，查明地下水的水位、流向、补给与排泄条件等。岩土体物理力学参数测定：采集滑坡体及滑带的岩土体样品，在实验室进行物理力学性质试验，测定岩土体的密度、含水量、孔隙比、抗剪强度等参数。这些参数是滑坡稳定性分析的关键数据，其准确性直接影响分析结果的可靠性。如通过直剪试验测定岩土体的内摩擦角和黏聚力，通过固结试验测定土体的压缩性等。降雨入渗规律研究：建立降雨入渗模型，利用数值模拟方法研究不同降雨强度、降雨持续时间条件下，雨水在滑坡体内的入渗过程，分析入渗深度、湿润锋推进速度、孔隙水压力分布等随时间的变化规律。考虑岩土体的非饱和特性，采用合适的渗流理论和数值算法，如基于Richards方程的有限元法，模拟降雨入渗过程，揭示降雨入渗对滑坡体渗流场的影响机制。滑坡稳定性分析：运用极限平衡法、数值分析法等多种方法，结合降雨入渗分析结果，计算不同降雨条件下吾隘滑坡的稳定系数，分析滑坡的稳定性变化趋势。极限平衡法采用常用的Fellenius法、Bishop法等，通过对滑坡体进行受力分析，计算滑坡的抗滑力和下滑力，进而得到稳定系数。数值分析法利用有限元软件如GeoStudio、FLAC等，建立滑坡的数值模型，考虑降雨入渗引起的孔隙水压力变化、土体强度降低等因素，模拟滑坡在不同降雨条件下的应力应变状态，评估滑坡的稳定性。滑坡稳定性影响因素敏感性分析：对影响吾隘滑坡稳定性的因素，如降雨强度、降雨持续时间、岩土体抗剪强度、孔隙水压力等进行敏感性分析，确定各因素对滑坡稳定性的影响程度，找出影响滑坡稳定性的关键因素。通过改变各因素的取值，计算滑坡的稳定系数，分析稳定系数随因素变化的敏感性，为滑坡的防治提供科学依据。例如，固定其他因素，逐步改变降雨强度，观察稳定系数的变化情况，确定降雨强度对滑坡稳定性的影响程度。本研究采用野外勘查、试验研究和数值模拟相结合的技术路线。首先，通过野外勘查，全面了解吾隘滑坡的地质条件和地形地貌特征，获取第一手资料。同时，采集岩土体样品，在实验室进行物理力学参数测定，为后续研究提供数据支持。其次，利用数值模拟软件建立降雨入渗模型和滑坡稳定性分析模型，对不同降雨条件下的滑坡稳定性进行模拟分析。在模拟过程中，充分考虑岩土体的物理力学性质、降雨特征以及渗流场与应力场的耦合作用。最后，对模拟结果进行分析和验证，结合实际情况，提出合理的滑坡防治建议。技术路线如图1-1所示：[此处插入技术路线图]二、吾隘滑坡区域自然地理与地质环境剖析2.1自然地理状况2.1.1地理位置与交通条件吾隘滑坡位于广西壮族自治区河池市南丹县吾隘镇境内，地处云贵高原南麓，地势西北高东南低。地理坐标约为东经[具体经度]，北纬[具体纬度]。该区域处于[具体山脉名称]山脉的延伸地带，周边山峦起伏，地形复杂，地势高差较大，为滑坡的形成提供了有利的地形条件。吾隘镇交通较为便利，[主要公路名称]公路穿镇而过，连接了周边的乡镇和县城，是当地居民出行和物资运输的主要通道。此外，吾隘镇距离[最近的铁路站点名称]铁路站点约[X]公里，通过公路与铁路站点相连，方便了人员和货物的长途运输。交通的便利性在一定程度上促进了当地的经济发展，但也增加了人类工程活动对地质环境的影响，如道路建设、切坡填方等工程活动，可能破坏山体的稳定性，引发滑坡等地质灾害。同时，交通设施在滑坡发生时也容易受到破坏，导致交通中断，影响救援和物资运输工作的开展。2.1.2气象与水文条件吾隘镇属于亚热带季风气候，夏季高温多雨，冬季温和少雨。多年平均气温为[X]℃，极端最高气温可达[X]℃，极端最低气温为[X]℃。年平均降水量为[X]毫米，降雨主要集中在5-9月，这期间的降水量占全年降水量的[X]%以上。降水的时空分布不均，暴雨频繁发生，短时间内的强降雨容易导致山体饱和，增加滑坡发生的风险。区内河流主要有[河流名称1]、[河流名称2]等，这些河流均属于[水系名称]水系。河流的流量受降水影响较大，在雨季时，河流水位迅速上涨，流量增大；在旱季时，河流水位下降，流量减小。河流对滑坡的影响主要体现在两个方面：一方面，河流的侵蚀作用可能削弱坡脚的抗滑力，使山体更容易发生滑坡；另一方面，河流水位的变化会导致地下水位的波动，进而影响滑坡体的稳定性。吾隘镇地下水类型主要有孔隙水、裂隙水和岩溶水。孔隙水主要赋存于第四系松散堆积层中，受大气降水和地表水的补给，水位和水量随季节变化明显。裂隙水主要分布在基岩的裂隙中，其富水性与岩石的裂隙发育程度、连通性以及地形地貌等因素有关。岩溶水主要存在于碳酸盐岩地层中，由于岩溶发育的不均匀性，岩溶水的分布和运动规律较为复杂。地下水的存在会增加岩土体的重量，降低其抗剪强度，同时，孔隙水压力的升高会使有效应力降低，进一步削弱滑坡体的稳定性。2.2地质环境特征2.2.1地形地貌吾隘滑坡所在区域总体地势起伏较大，属于构造侵蚀低山地貌。滑坡体位于山坡中下部，地形坡度较陡，一般在25°-40°之间，局部地段可达50°以上。山坡自然坡度较陡，使得滑坡体在重力作用下具有较大的下滑趋势。滑坡后壁高程约为[X]米，前缘高程约为[X]米，相对高差约为[X]米。滑坡体平面形态呈不规则状，长约[X]米，宽约[X]米，面积约为[X]平方米。该区域的地形地貌对滑坡的形成和发展具有重要作用。陡峭的地形坡度为滑坡的发生提供了势能条件，在降雨等外部因素的作用下，岩土体更容易沿着潜在滑动面下滑。地形起伏较大，使得地表水和地下水的径流速度加快，容易对山坡土体产生冲刷和侵蚀作用，削弱坡体的稳定性。此外，地形地貌还影响着降雨的分布和入渗条件，如山坡上部由于地形较高，降雨更容易形成地表径流，而山坡下部则更容易受到径流的冲刷和地下水的浸泡，从而增加了滑坡发生的可能性。2.2.2地层岩性滑坡区域内地层主要由第四系全新统坡残积层（Q4dl+el）和三叠系中统百逢组（T2bf）组成。第四系全新统坡残积层主要分布于滑坡体表层，岩性为粉质黏土夹碎石，厚度一般在2-5米之间，局部地段可达8米。粉质黏土呈可塑-硬塑状，碎石含量约为20%-40%，粒径一般在2-20厘米之间，成分主要为砂岩和页岩。该层土体结构松散，透水性较好，在降雨作用下容易饱和，导致土体容重增大，抗剪强度降低。三叠系中统百逢组为滑坡体的下伏基岩，岩性主要为砂岩、页岩互层。砂岩呈灰白色、灰黄色，中细粒结构，层状构造，岩石较坚硬，抗压强度较高。页岩呈黑色、灰黑色，页理发育，岩石较软弱，抗风化能力较弱。砂岩和页岩的互层结构使得岩体的力学性质存在明显差异，在地下水的作用下，页岩容易软化、泥化，形成相对软弱的滑动面，从而影响滑坡体的稳定性。地层岩性与滑坡稳定性密切相关。第四系坡残积层由于结构松散、透水性好，在降雨时容易吸收大量水分，增加土体重量，降低抗剪强度，成为滑坡发生的主要物质来源。而基岩中的页岩软弱夹层，在地下水的长期作用下，强度降低，容易形成滑动面，控制着滑坡的滑动方向和规模。例如，当降雨入渗导致坡残积层饱和后，土体的下滑力增大，而页岩夹层的抗滑力相对较弱，就容易引发滑坡的发生。2.2.3地质构造研究区位于[具体构造单元名称]构造单元内，区域内地质构造较为复杂，主要发育有褶皱和断裂构造。褶皱构造表现为一系列紧闭的背斜和向斜，轴向大致为[具体方向]。断裂构造主要有[断裂名称1]、[断裂名称2]等，断裂走向多为[具体走向]，倾向[具体倾向]，倾角一般在60°-80°之间。地质构造对滑坡稳定性产生重要影响。一方面，褶皱和断裂构造破坏了岩体的完整性，增加了岩体的透水性，使得降雨更容易入渗到地下，从而改变岩土体的物理力学性质，降低滑坡体的稳定性。例如，断裂带附近的岩体破碎，裂隙发育，降雨入渗后，地下水在裂隙中流动，对岩体产生动水压力和软化作用，削弱了岩体的抗剪强度。另一方面，地质构造控制着山体的形态和地形地貌，影响着滑坡的形成和发展。例如，背斜顶部由于岩层张应力作用，岩体破碎，容易遭受风化和侵蚀，形成陡崖和斜坡，为滑坡的发生提供了有利的地形条件。2.2.4新构造运动及地震研究区新构造运动较为活跃，主要表现为间歇性的地壳抬升和差异性升降运动。长期的地壳抬升使得山体不断隆升，地形高差增大，河流下切作用加强，从而导致山坡坡度变陡，岩土体稳定性降低，增加了滑坡发生的可能性。差异性升降运动使得区域内地形起伏变化较大，在地形突变部位容易形成应力集中，破坏岩土体的结构，引发滑坡等地质灾害。该区域地震活动相对较弱，历史上记录到的地震震级一般在4级以下。虽然地震活动不强，但地震产生的地震波会对岩土体产生震动作用，使岩土体的结构受到破坏，强度降低。在地震作用下，滑坡体的稳定性会受到显著影响，尤其是对于已经处于临界稳定状态的滑坡，地震可能成为诱发滑坡的直接因素。例如，地震产生的水平和垂直地震力会增加滑坡体的下滑力，同时破坏滑带土的结构，降低其抗剪强度，从而导致滑坡的发生。2.2.5水文地质条件吾隘滑坡区地下水类型主要为孔隙水和裂隙水。孔隙水主要赋存于第四系坡残积层中，受大气降水和地表水的补给，水位随季节变化明显。在雨季，大气降水充沛，孔隙水水位迅速上升；在旱季，随着蒸发和排泄作用，孔隙水水位逐渐下降。裂隙水主要存在于基岩的裂隙中，其富水性与裂隙的发育程度、连通性以及地形地貌等因素有关。一般来说，断裂带和褶皱轴部附近的裂隙较为发育，富水性较好。地下水的水位、流向和含水层分布对滑坡稳定性有着重要影响。当孔隙水水位上升时，会使滑坡体的含水量增加，容重增大，从而增加下滑力。同时，孔隙水压力的升高会使有效应力降低，削弱岩土体的抗剪强度。例如，在强降雨后，孔隙水水位迅速上升，滑坡体的稳定性会明显下降。地下水的流向也会影响滑坡的稳定性，如果地下水的流向与滑坡的滑动方向一致，会对滑坡体产生动水压力，进一步推动滑坡体下滑。含水层的分布则决定了地下水的赋存和运移条件，不同含水层之间的水力联系也会影响滑坡体的稳定性。如当存在多层含水层且水力联系密切时，地下水的运动更加复杂，对滑坡稳定性的影响也更为显著。三、吾隘滑坡基本特性探究3.1滑坡边界、形态及规模特点通过高精度的全站仪测量、航空遥感影像解译以及地面地质测绘等综合手段，对吾隘滑坡的边界进行了精确圈定。滑坡的后缘边界位于山体的较高部位，以明显的弧形拉张裂缝为标志，裂缝宽度在5-30厘米之间，延伸长度约为200米。后缘边界的地形相对较高，坡度较陡，一般在35°-45°之间。滑坡的两侧边界则依据地形地貌特征和岩土体的错动痕迹来确定，两侧边界的岩土体存在明显的错动和位移现象，错动距离在0.5-2米之间。滑坡的前缘边界较为平缓，以堆积的滑坡体和鼓张裂缝为标志，前缘边界与河流或道路相邻，对周边的基础设施构成一定威胁。在平面形态上，吾隘滑坡整体呈现出不规则的舌状，其长轴方向大致与山坡的走向一致，约为东北-西南向。滑坡的舌状形态在遥感影像上清晰可见，其前缘部分向下方突出，形似舌头，表明滑坡体在滑动过程中受到地形和岩土体性质的影响，呈现出一定的方向性。从剖面形态来看，滑坡体呈现出上陡下缓的特征，上部滑坡体的坡度一般在30°-40°之间，下部滑坡体的坡度逐渐变缓，在15°-25°之间。滑坡体的剖面形态与岩土体的结构和力学性质密切相关，上部岩土体由于受到重力和降雨的作用，更容易发生变形和滑动，因此坡度较陡；下部岩土体则受到上部岩土体的挤压和堆积作用，坡度相对较缓。吾隘滑坡的规模属于中型滑坡，其长度约为300米，宽度约为150米，平均厚度约为8米。根据体积计算公式V=L×W×H（其中V为体积，L为长度，W为宽度，H为厚度），可估算出滑坡体的体积约为36万立方米。近年来，由于持续的降雨和人类工程活动的影响，滑坡体的规模有逐渐扩大的趋势。例如，在2022年的强降雨过程中，滑坡体的后缘出现了新的拉张裂缝，裂缝长度增加了约50米，前缘部分也有一定程度的下滑，导致滑坡体的长度和宽度略有增加。人类工程活动如道路建设、切坡填方等，破坏了山体的原有稳定性，也使得滑坡体的规模不断扩大。3.2滑坡结构特征3.2.1滑体特征吾隘滑坡的滑体主要由第四系全新统坡残积层构成，岩性为粉质黏土夹碎石。粉质黏土呈可塑-硬塑状态，具有一定的黏性和可塑性，但在饱水状态下，其强度会显著降低。碎石含量在20%-40%之间，粒径范围为2-20厘米，成分主要为砂岩和页岩。这种物质组成使得滑体的结构较为松散，孔隙率较大，透水性相对较好。滑体的结构呈现出明显的不均一性，在垂直方向上，上部粉质黏土含量相对较高，结构较为细腻，而下部碎石含量逐渐增加，结构变得较为粗糙。在水平方向上，由于滑坡体的滑动和堆积作用，不同部位的物质组成和结构也存在差异。例如，滑坡体的前缘部分，由于受到滑动的挤压和摩擦作用，碎石含量相对较高，且颗粒排列较为紧密；而后缘部分则粉质黏土含量相对较多，结构较为松散。滑体厚度变化较大，在滑坡体的中部和后缘，厚度一般在5-10米之间，而在滑坡体的前缘，厚度相对较薄，一般在2-5米之间。滑体厚度的变化与地形地貌、岩土体性质以及滑坡的滑动过程密切相关。在地形坡度较陡的部位，滑体厚度相对较大，这是因为在重力作用下，岩土体更容易在这些部位堆积；而在地形较为平缓的部位，滑体厚度相对较薄。此外，岩土体的抗滑能力也会影响滑体的厚度，抗滑能力较弱的部位，滑体更容易发生滑动和堆积，从而导致厚度增加。滑体的物质组成、结构和厚度变化对降雨入渗和滑坡稳定性有着重要影响。由于滑体结构松散、透水性好，降雨时雨水能够迅速入渗到滑体内部，增加滑体的含水量和重量，从而降低滑体的抗剪强度。例如，当降雨强度较大时，大量雨水在短时间内入渗到滑体中，使滑体迅速饱和，容重增大，下滑力显著增加。同时，滑体厚度的变化也会影响滑坡的稳定性，厚度较大的部位，下滑力相对较大，更容易发生滑坡。而且，滑体结构的不均一性会导致雨水入渗的不均匀性，使得滑体内部不同部位的含水量和力学性质存在差异，进一步影响滑坡的稳定性。3.2.2滑床特征滑床为滑坡体滑动时所依附的下伏基岩，吾隘滑坡的滑床主要由三叠系中统百逢组的砂岩、页岩互层组成。砂岩呈灰白色、灰黄色，具有中细粒结构和层状构造，岩石较为坚硬，抗压强度较高，一般在30-50MPa之间。页岩呈黑色、灰黑色，页理发育，岩石相对软弱，抗风化能力较弱，抗压强度一般在10-20MPa之间。砂岩和页岩的互层结构使得滑床的力学性质存在明显的各向异性，在垂直于层面方向上的强度较高，而在平行于层面方向上的强度较低。滑床的产状对滑坡的滑动方向和规模具有重要控制作用。通过地质测绘和钻探资料分析，滑床的岩层倾向与山坡的坡向基本一致，倾角一般在20°-30°之间。这种产状使得滑坡体在重力作用下容易沿着滑床的层面发生滑动，滑动方向与岩层倾向一致。滑床的倾角大小影响着滑坡体的下滑速度和势能，倾角越大，滑坡体下滑的速度越快，势能也越大，对滑坡的稳定性影响也更为显著。滑床的力学性质对滑坡的滑动起到关键控制作用。由于页岩的强度较低，在地下水和滑坡体的长期作用下，页岩容易发生软化、泥化，形成相对软弱的滑动面。当滑坡体的下滑力超过滑床的抗滑力时，滑坡就会发生。例如，在强降雨或地震等因素的作用下，滑床的抗滑力会进一步降低，而滑坡体的下滑力会增加，从而导致滑坡的发生。此外，滑床的力学性质还会影响滑坡的滑动方式，当滑床的抗滑力在不同部位存在差异时，滑坡体可能会出现局部滑动或整体滑动的情况。3.2.3滑带特征滑带是滑坡体与滑床之间的滑动面，吾隘滑坡的滑带物质主要由粉质黏土、页岩碎屑和少量的砂岩颗粒组成。这些物质在滑坡滑动过程中受到强烈的剪切作用，结构被破坏，颗粒重新排列，形成了相对软弱的滑动带。滑带土的含水量较高，一般在25%-35%之间，处于饱和或近饱和状态。高含水量使得滑带土的抗剪强度显著降低，这是因为水分的存在会削弱颗粒之间的摩擦力和黏聚力。滑带的厚度相对较薄，一般在0.2-0.5米之间，但在局部地段可能会有所增厚。滑带厚度的变化与滑坡的滑动历史、滑体和滑床的物质组成以及地下水的活动密切相关。在滑坡滑动过程中，滑带土受到反复的剪切和挤压作用，其厚度可能会发生变化。例如，在滑坡的初始滑动阶段，滑带厚度可能较薄，随着滑动的持续进行，滑带土不断被剪切破坏，厚度可能会逐渐增加。此外，滑体和滑床物质的差异以及地下水的浸泡作用也会影响滑带的厚度。滑带的物理力学性质在降雨作用下会发生显著变化。当降雨入渗到滑带时，滑带土的含水量进一步增加，容重增大，抗剪强度降低。研究表明，滑带土的抗剪强度与含水量呈负相关关系，随着含水量的增加，内摩擦角和黏聚力都会减小。例如，当滑带土的含水量从25%增加到35%时，内摩擦角可能会从20°降低到15°，黏聚力可能会从10kPa降低到5kPa。这种变化会导致滑坡体的稳定性急剧下降，增加滑坡发生的风险。同时，降雨还可能导致滑带土的软化和泥化，进一步破坏滑带的结构，降低其抗滑能力。3.3滑坡变形特征通过现场勘查和监测数据可知，吾隘滑坡存在明显的变形迹象。滑坡后缘出现了多条拉张裂缝，这些裂缝呈弧形分布，长度从数米到数十米不等，宽度在5-30厘米之间，部分裂缝深度可达1-2米。拉张裂缝的出现表明滑坡体后缘受到了拉伸作用，土体发生了分离和位移，这是滑坡变形的重要标志之一。滑坡体表面还存在一些局部的塌陷和鼓胀现象，在滑坡体的中部和前缘较为明显。塌陷区域的土体下沉，形成凹坑，深度一般在0.5-1米之间；鼓胀区域的土体则向上隆起，形成小丘，高度在0.3-0.8米之间。这些塌陷和鼓胀现象是由于滑坡体内部的土体在滑动过程中受到挤压和拉伸，导致土体结构发生变化而形成的。吾隘滑坡的变形历史较为复杂，据当地居民回忆和相关资料记载，该滑坡在过去几十年间曾多次发生小规模的滑动。在20世纪80年代，由于连续的强降雨，滑坡体出现了轻微的下滑，导致后缘出现了一些细小的裂缝，但当时并未引起足够的重视。此后，随着时间的推移和人类工程活动的影响，滑坡体的变形逐渐加剧。在2010年前后，滑坡体的前缘开始出现局部的坍塌，一些房屋受到了不同程度的损坏。近年来，由于极端降雨事件的增多，吾隘滑坡的变形更加明显，2022年5月30日，在强降雨的作用下，滑坡发生滑动，滑动方量约70立方米，滑坡体堆积至房屋后墙，所幸未造成人员伤亡和经济损失。通过对监测数据的分析，可以发现吾隘滑坡的变形呈现出一定的趋势。在降雨期间，滑坡体的变形速率明显增大，尤其是在强降雨过后的几天内，变形速率达到峰值。例如，在2022年5月的强降雨过程中，滑坡体的水平位移速率在降雨后的3天内达到了每天10-15毫米，垂直位移速率也达到了每天5-8毫米。随着时间的推移，变形速率逐渐减小，但滑坡体仍处于持续变形的状态。从长期来看，滑坡体的变形有逐渐加剧的趋势，这可能与降雨的频繁发生、人类工程活动的不断增加以及岩土体性质的逐渐劣化等因素有关。如果不采取有效的防治措施，滑坡体可能会发生大规模的滑动，对周边居民的生命财产安全构成严重威胁。3.4滑坡渗透及强度特征3.4.1滑坡体原位双环渗透试验为准确获取吾隘滑坡体的渗透性能，采用原位双环渗透试验方法。该试验在滑坡体不同位置选取了[X]个具有代表性的试验点，分别位于滑坡体的前缘、中部和后缘，以确保试验结果能够全面反映滑坡体的渗透特性。试验点的选择充分考虑了地形地貌、岩土体物质组成等因素，避免了因局部特殊情况对试验结果产生偏差。试验过程严格按照相关规范进行。首先，在选定的试验点处，用人工挖掘一个直径约为[X]米、深度约为[X]米的圆形试坑，以保证试验环能够垂直嵌入且周围土体不受扰动。然后，将两个高度为[X]厘米、直径分别为[X]厘米和[X]厘米的铁环同心垂直嵌入试坑底部，两环之间的间距保持在[X]厘米左右，确保内环水的垂向入渗不受外环的干扰。试验时，使用高精度的量筒同时向内、外铁环内注水，通过调整注水速度，使内外环的水柱高度始终保持在同一高度，控制在[X]厘米左右。采用秒表精确记录注水时间和注水量，每隔[X]分钟记录一次数据，直至单位时间内的渗入水量基本稳定，即达到稳定入渗状态。在整个试验过程中，密切观察试坑周围土体的变化情况，确保试验的准确性和安全性。试验结果表明，吾隘滑坡体的渗透系数呈现出明显的空间差异性。在滑坡体前缘，由于岩土体结构相对松散，碎石含量较高，渗透系数较大，平均值约为[X]×10⁻⁴厘米/秒。中部区域的渗透系数相对较小，平均值约为[X]×10⁻⁵厘米/秒，这主要是因为中部粉质黏土含量相对较多，结构较为紧密。后缘的渗透系数介于前缘和中部之间，平均值约为[X]×10⁻⁴厘米/秒。通过对不同位置渗透系数的对比分析，可以看出滑坡体的渗透性能与岩土体的物质组成和结构密切相关。滑坡体的渗透性能对降雨入渗有着重要影响。较大的渗透系数意味着降雨能够更快地入渗到滑坡体内部，使滑坡体在短时间内含水量迅速增加。以滑坡体前缘为例，在强降雨条件下，由于其渗透系数较大，雨水能够快速入渗，导致前缘土体迅速饱和，容重增大，抗剪强度降低。同时，快速入渗的雨水还可能在滑坡体内部形成较高的孔隙水压力，进一步削弱土体的抗滑能力，增加滑坡发生的风险。而对于渗透系数较小的区域，如滑坡体中部，降雨入渗速度相对较慢，但随着降雨时间的延长，也会逐渐使土体饱和，对滑坡稳定性产生不利影响。因此，滑坡体的渗透性能是研究降雨入渗条件下滑坡稳定性的关键因素之一。3.4.2岩土体物理力学试验岩土体物理力学试验是获取滑坡稳定性分析所需参数的重要手段。本次试验主要包括密度试验、含水量试验、颗粒分析试验、直剪试验和压缩试验等。在滑坡体和滑带不同位置共采集了[X]组岩土体样品，确保样品具有代表性。密度试验采用环刀法，通过测量环刀内岩土体的质量和体积，计算出岩土体的天然密度。含水量试验则采用烘干法，将岩土体样品在105-110℃的烘箱中烘干至恒重，根据烘干前后的质量差计算含水量。颗粒分析试验采用筛分法和比重计法，对不同粒径的颗粒含量进行测定，以分析岩土体的颗粒级配。直剪试验采用快剪和固结快剪两种方法，测定岩土体的抗剪强度指标，包括内摩擦角和黏聚力。压缩试验则通过对岩土体样品施加不同的压力，测定其压缩性指标，如压缩系数和压缩模量。通过试验，获取了吾隘滑坡岩土体的物理力学参数。滑坡体粉质黏土的天然密度为[X]克/立方厘米，含水量为[X]%，内摩擦角为[X]°，黏聚力为[X]kPa，压缩系数为[X]MPa⁻¹，压缩模量为[X]MPa。滑带土的天然密度为[X]克/立方厘米，含水量为[X]%，内摩擦角为[X]°，黏聚力为[X]kPa，压缩系数为[X]MPa⁻¹，压缩模量为[X]MPa。这些参数为后续的滑坡稳定性分析提供了重要的数据支持。在使用极限平衡法计算滑坡稳定系数时，需要用到岩土体的容重、内摩擦角和黏聚力等参数，而这些参数正是通过物理力学试验获得的。在数值模拟分析中，也需要根据试验参数来设定岩土体的本构模型和力学参数，以保证模拟结果的准确性。因此，准确获取岩土体的物理力学参数对于深入研究吾隘滑坡的稳定性具有重要意义。四、吾隘滑坡变形影响因素与机理阐释4.1吾隘滑坡变形影响因素解析4.1.1内在因素分析岩土体性质对吾隘滑坡稳定性起着至关重要的控制作用。滑坡体主要由粉质黏土夹碎石组成，粉质黏土的可塑-硬塑状态使其在饱水条件下强度降低明显。碎石的存在虽在一定程度上增加了土体的骨架支撑，但由于其含量和分布的不均匀性，导致土体结构松散，孔隙率较大，透水性增强。这种物质组成使得滑坡体在降雨时容易饱和，容重增大，抗剪强度降低。如在实验室直剪试验中，当粉质黏土含水量从15%增加到25%时，内摩擦角从25°降低到20°，黏聚力从20kPa降低到15kPa。滑带土的特殊性质也是影响滑坡稳定性的关键因素，其高含水量和相对较低的抗剪强度，使得滑带成为滑坡滑动的薄弱环节。地质构造对吾隘滑坡的影响不容忽视。研究区处于复杂的地质构造区域，褶皱和断裂构造破坏了岩体的完整性。褶皱使得岩层发生弯曲变形，在褶皱轴部，岩体受拉伸和挤压作用，裂隙发育，岩石破碎，为降雨入渗提供了通道。断裂构造则进一步加剧了岩体的破碎程度，地下水沿着断裂带流动，对岩体产生软化和侵蚀作用，降低了岩体的强度。例如，在吾隘滑坡附近的[具体断裂名称]断裂带，岩体破碎带宽达数米，裂隙纵横交错，地下水活动频繁，导致该区域的岩土体稳定性极差。地质构造还控制着山体的地形地貌，影响着滑坡的形成和发展。地形地貌是吾隘滑坡形成的重要基础条件。滑坡所在区域地势起伏较大，山坡坡度较陡，一般在25°-40°之间，局部可达50°以上。陡峭的地形为滑坡提供了较大的势能，在重力作用下，岩土体具有向坡下滑动的趋势。地形的起伏变化还影响着地表水和地下水的流动，山坡上部的降雨容易形成地表径流，对坡面产生冲刷作用，破坏坡体的稳定性。而山坡下部由于汇水作用，地下水位较高，岩土体长期处于饱水状态，抗剪强度降低，增加了滑坡发生的风险。例如，在吾隘滑坡的前缘，由于地势较低，地表水和地下水汇聚，导致该区域的岩土体长期处于饱和状态，多次发生局部坍塌现象。4.1.2外在因素分析降雨是诱发吾隘滑坡的最主要外在因素。河池市降雨充沛，且降雨集中在5-9月，暴雨频繁。降雨通过多种机制影响滑坡稳定性。一方面，降雨入渗使滑坡体含水量增加，容重增大，下滑力增大。根据饱和-非饱和渗流理论，降雨入渗过程中，水分在土体孔隙中运移，当土体饱和度达到一定程度后，孔隙水压力升高，有效应力降低，岩土体抗剪强度随之降低。研究表明，当孔隙水压力增加10kPa时，岩土体的抗剪强度可能降低10%-20%。另一方面，长期降雨会使岩土体软化、泥化，尤其是滑带土，在雨水浸泡下，强度大幅下降，进一步降低了滑坡的稳定性。例如，在2022年5月的强降雨过程中，吾隘滑坡体的含水量急剧增加，孔隙水压力迅速上升，导致滑坡发生滑动。地震对吾隘滑坡稳定性也有一定影响。虽然该区域地震活动相对较弱，但地震产生的地震波会对岩土体产生震动作用，使岩土体结构受到破坏，强度降低。在地震作用下，滑坡体的稳定性会受到显著影响。地震力会增加滑坡体的下滑力，同时破坏滑带土的结构，降低其抗剪强度。根据地震动力学理论，地震加速度每增加0.1g，滑坡体的下滑力可能增加10%-15%。对于已经处于临界稳定状态的滑坡，地震可能成为诱发滑坡的直接因素。如在历史上的某次小震中，吾隘滑坡体出现了微小的位移和裂缝扩展现象。人类活动对吾隘滑坡稳定性的影响日益显著。随着当地经济的发展，人类工程活动不断增加，如道路建设、切坡填方、房屋修建等。道路建设过程中的切坡填方破坏了山体的原有稳定性，形成了新的临空面，增加了滑坡发生的可能性。切坡导致坡体应力重新分布，在坡脚处形成应力集中，容易引发滑坡。填方则增加了坡体的重量，使下滑力增大。房屋修建过程中，不合理的基础开挖和排水设施设置也会影响滑坡体的稳定性。此外，人类活动还可能改变地表水和地下水的径流路径，导致坡体含水量增加，进一步降低滑坡的稳定性。例如，吾隘镇在道路拓宽工程中，对山坡进行了大规模切坡，切坡后不久，在降雨作用下，切坡部位发生了小规模滑坡。4.2滑坡变形机理分析在降雨条件下，吾隘滑坡的变形破坏是一个复杂的过程，涉及多种物理力学机制。降雨初期，雨水首先在滑坡体表面形成径流，部分径流通过土体孔隙和裂隙入渗到滑坡体内部。由于滑坡体主要由粉质黏土夹碎石组成，结构松散，孔隙率较大，使得雨水能够较快地入渗。随着入渗的进行，滑坡体的含水量逐渐增加，容重增大，根据重力计算公式G=ρgV（其中G为重力，ρ为密度，g为重力加速度，V为体积），容重的增大导致滑坡体的重力增大，下滑力随之增大。同时，降雨入渗使孔隙水压力升高，有效应力降低。根据有效应力原理，有效应力等于总应力减去孔隙水压力。当孔隙水压力升高时，有效应力降低，岩土体颗粒间的摩擦力减小，抗剪强度降低。在非饱和土中，随着含水量的增加，基质吸力减小，也会导致抗剪强度降低。例如，当孔隙水压力从0增加到10kPa时，滑带土的抗剪强度可能降低10%-20%。这使得滑坡体在重力作用下更容易发生变形和滑动。随着降雨的持续，滑坡体内部的水分不断增加，岩土体逐渐软化、泥化，尤其是滑带土。滑带土在雨水的长期浸泡下，结构被破坏，颗粒间的连接力减弱，抗剪强度进一步降低。研究表明，滑带土的抗剪强度与含水量呈负相关关系，含水量每增加5%，抗剪强度可能降低10%-15%。当滑坡体的下滑力超过滑带土的抗滑力时，滑带土开始发生剪切破坏，形成滑动面。滑动面一旦形成，滑坡体就会沿着滑动面发生滑动，导致滑坡的发生。在滑坡变形过程中，岩土体的力学响应也十分复杂。当滑坡体受到外力作用时，会产生应力应变。在弹性阶段，岩土体的应力应变关系符合胡克定律，但随着外力的增加，岩土体进入塑性阶段，发生不可逆的变形。在降雨条件下，由于岩土体性质的变化，其力学响应也会发生改变。例如，含水量的增加会使岩土体的弹性模量降低，泊松比增大，导致岩土体更容易发生变形。此外，滑坡体内部的应力分布也会发生变化，在滑坡体的后缘和两侧，由于受到拉应力作用，容易出现拉张裂缝；在滑坡体的前缘，由于受到挤压应力作用，可能会出现鼓胀和坍塌现象。五、不同降雨条件下滑坡渗流场模拟与分析5.1降雨入渗基本理论5.1.1饱和—非饱和渗基本理论饱和—非饱和渗流理论是研究降雨入渗过程中水分在岩土体中运动规律的重要基础。在饱和状态下，岩土体孔隙完全被水充满，水的流动主要受重力和压力梯度的驱动，其渗流规律遵循达西定律，表达式为v=ki，其中v为渗流速度，k为渗透系数，i为水力梯度。达西定律表明，渗流速度与水力梯度成正比，渗透系数反映了岩土体的透水性能。然而，在实际的降雨入渗过程中，岩土体往往处于非饱和状态，孔隙中同时存在水和空气。非饱和渗流中，水的流动不仅受重力和压力梯度的影响，还受到基质吸力的作用。基质吸力是指非饱和土中孔隙水与空气交界面处由于表面张力而产生的一种吸力，它使得水分在孔隙中保持一定的势能，阻碍水分的流动。非饱和渗流的基本方程为Richards方程，其表达式为：C(h)\frac{\partialh}{\partialt}=\nabla\cdot(K(h)\nablah)+\nabla\cdot(K(h)\nablaz)-\frac{\partialK(h)}{\partialz}其中，h为压力水头，C(h)为比水容量，表示单位体积岩土体中，压力水头变化单位值时含水量的变化量，反映了岩土体对水分的储存能力；K(h)为非饱和渗透系数，是压力水头的函数，随压力水头的变化而变化，体现了非饱和状态下岩土体的透水性能；t为时间，z为垂直坐标。Richards方程综合考虑了重力、压力和基质吸力对水分运动的影响，能够更准确地描述非饱和渗流过程。非饱和渗透系数与基质吸力密切相关，一般来说，随着基质吸力的增大，非饱和渗透系数减小。这是因为基质吸力的增大使得孔隙水的弯液面曲率增大，水分在孔隙中的流动阻力增加。许多学者通过实验和理论分析，建立了非饱和渗透系数与基质吸力的关系模型，如VanGenuchten模型、Brooks-Corey模型等。以VanGenuchten模型为例，其表达式为：K(h)=K_sS_e^{\frac{1}{2}}\left[1-(1-S_e^{\frac{1}{m}})^m\right]^2其中，K_s为饱和渗透系数，S_e为有效饱和度，m为与土性有关的参数。该模型通过有效饱和度将非饱和渗透系数与基质吸力联系起来，能够较好地描述非饱和渗透系数随基质吸力的变化规律。5.1.2边坡降雨入渗过程降雨入渗边坡是一个复杂的动态过程，可分为以下几个阶段：溅渗阶段：降雨初期，雨滴以一定的速度撞击坡面，部分雨滴溅起，形成微小的水花，另一部分雨滴则直接渗入土体表面。在这个阶段，由于土体表面较为干燥，孔隙较大，雨水能够迅速入渗。但由于降雨强度可能大于土体的入渗能力，部分雨水会在坡面形成薄层水流，即坡面径流。例如，在小雨或短时间降雨的情况下，溅渗阶段可能较为短暂，大部分雨水能够迅速入渗到土体中。下渗阶段：随着降雨的持续，雨水不断向土体内部渗透，形成下渗水流。在非饱和土体中，水分主要通过孔隙和裂隙进行运移。由于土体颗粒表面存在吸附力，水分在入渗过程中会受到一定的阻力。随着入渗深度的增加，基质吸力逐渐增大，非饱和渗透系数逐渐减小，水分入渗速度逐渐减慢。在这个阶段，降雨入渗主要受土体的渗透性、初始含水量、降雨强度等因素的影响。例如，对于渗透性较好的砂土，雨水入渗速度较快；而对于渗透性较差的黏土，雨水入渗速度较慢。再分布阶段：当降雨停止后，土体中的水分开始重新分布。在重力和基质吸力的作用下，入渗到土体中的水分会继续向深层渗透，同时也会向周围扩散。在这个过程中，土体的含水量逐渐趋于均匀，但由于土体的非均质性，不同部位的含水量仍可能存在差异。例如，在坡体的低洼处，水分可能会相对集中，而在坡体的高处，水分可能会相对较少。蒸发和排泄阶段：在降雨停止后的一段时间内，土体中的水分还会通过蒸发和排泄的方式逐渐减少。蒸发主要发生在土体表面，受气温、湿度、风速等因素的影响。排泄则是指土体中的水分通过地下径流的方式排出坡体，受地下水位、水力梯度等因素的影响。随着蒸发和排泄的进行，土体的含水量逐渐恢复到初始状态。降雨入渗对滑坡渗流场产生显著影响。随着降雨入渗，滑坡体的含水量增加，地下水位上升，孔隙水压力增大。孔隙水压力的增大使得有效应力减小，根据有效应力原理\sigma'=\sigma-u（其中\sigma'为有效应力，\sigma为总应力，u为孔隙水压力），有效应力的减小会导致岩土体的抗剪强度降低。研究表明，当孔隙水压力增加10kPa时，岩土体的抗剪强度可能降低10%-20%。降雨入渗还会改变滑坡体的渗流路径和渗流速度，使得渗流场分布更加复杂。例如，在滑坡体的局部区域，由于岩土体的渗透性差异或地形的影响，可能会形成渗流集中区，导致该区域的孔隙水压力异常升高，进一步降低滑坡体的稳定性。5.2不同降雨条件下滑坡渗流场数值模拟5.2.1模拟方法简介本研究采用GeoStudio软件中的SEEP/W模块进行吾隘滑坡渗流场的数值模拟。GeoStudio是一款专门用于岩土工程分析的软件，具有强大的数值计算和图形显示功能，在滑坡渗流场研究领域得到了广泛应用。SEEP/W模块基于有限元方法，能够准确模拟饱和-非饱和渗流过程，考虑了重力、压力和基质吸力对水分运动的影响，能够较好地反映降雨入渗条件下边坡渗流场的变化规律。与其他数值模拟方法相比，GeoStudio软件的SEEP/W模块具有以下优势：一是其前处理功能强大，能够方便地导入CAD图形，快速建立复杂的地质模型。通过将吾隘滑坡的地质剖面图导入软件，能够准确地定义滑坡体、滑带和滑床的几何形状和位置，提高模型的准确性。二是软件提供了丰富的材料本构模型和参数设置选项，能够根据岩土体的实际物理力学性质进行合理设置。在吾隘滑坡的模拟中，可以根据室内试验获取的岩土体参数，如渗透系数、孔隙率、饱和度等，准确地定义材料属性，使模拟结果更加符合实际情况。三是SEEP/W模块能够考虑多种边界条件，如降雨入渗边界、定水头边界、零流量边界等，能够真实地模拟滑坡渗流场的边界条件。在吾隘滑坡的模拟中，根据其实际的水文地质条件，设置了坡体表面的降雨入渗边界、基岩面的零流量边界和滑坡后缘稳定地下水位处的定水头边界，使模拟结果更加可靠。四是该模块的计算结果可视化程度高，能够直观地展示渗流场的分布和变化情况。通过绘制孔隙水压力等值线图、地下水位线图和渗流速度矢量图等，可以清晰地观察到不同降雨条件下吾隘滑坡渗流场的变化特征，为分析滑坡稳定性提供有力支持。5.2.2模型的建立与方案设计在建立吾隘滑坡数值模型时，首先利用专业的测绘仪器对滑坡进行详细的地形测量，获取高精度的地形数据。同时，结合前期的地质勘察资料，包括钻孔数据、地质剖面图等，在GeoStudio软件中准确绘制滑坡的几何形状，包括滑坡体、滑带和滑床的位置和范围。模型的边界条件设置如下：坡体表面设置为降雨入渗边界，根据当地的降雨数据和研究需求，确定降雨强度和降雨历时；基岩面设置为隔水零流量边界，以模拟基岩的不透水特性；滑坡后缘稳定地下水位处设置为定水头边界，确保后缘地下水位的稳定。岩土体参数根据室内物理力学试验结果进行赋值，主要包括渗透系数、孔隙率、饱和度、比水容量等。对于渗透系数，采用原位双环渗透试验和室内渗透试验相结合的方法确定，考虑到岩土体的非均质性，在模型中对不同区域的渗透系数进行了合理的分区赋值。孔隙率和饱和度根据试验数据和经验取值，比水容量则根据相关理论公式计算得到。模型的网格划分采用自适应网格划分技术，在滑坡体和滑带等关键部位进行加密处理，以提高计算精度。经过多次调试和优化，最终建立了能够准确反映吾隘滑坡地质条件的数值模型。为了全面研究不同降雨条件下吾隘滑坡渗流场的变化规律，设计了以下模拟方案：一是不同降雨强度工况，设置小雨（降雨强度为10mm/d）、中雨（降雨强度为30mm/d）、大雨（降雨强度为50mm/d）、暴雨（降雨强度为100mm/d）等不同降雨强度，分别模拟在这些降雨强度下渗流场随时间的变化。二是不同降雨历时工况，设置降雨历时为1d、3d、5d、7d等，研究在不同降雨历时条件下渗流场的动态变化。三是相同总降雨量工况，通过不同的降雨强度和历时组合，使总降雨量相同，如降雨强度为50mm/d、历时2d和降雨强度为25mm/d、历时4d等，对比分析不同降雨组合下渗流场的差异。四是一次降雨历程工况，根据当地实际的降雨过程数据，模拟一次完整的降雨过程中吾隘滑坡渗流场的变化。通过这些模拟方案，可以系统地研究降雨强度、降雨历时以及降雨组合等因素对吾隘滑坡渗流场的影响。5.2.3不同降雨强度下吾隘滑坡渗流场分析通过对不同降雨强度下吾隘滑坡渗流场的模拟分析，得到了一系列关于地下水位、孔隙水压力等参数的变化规律。在小雨（降雨强度为10mm/d）条件下，降雨初期，雨水缓慢入渗，地下水位上升较为缓慢。随着降雨的持续，地下水位逐渐上升，但上升幅度较小。在降雨结束后，地下水位逐渐稳定，上升高度约为0.5-1米。孔隙水压力在滑坡体表层有所增加，但增加幅度较小，且随着深度的增加，孔隙水压力变化不明显。这是因为小雨条件下，降雨入渗量较小，对滑坡体内部的渗流场影响有限。当中雨（降雨强度为30mm/d）发生时，地下水位上升速度明显加快。在降雨初期，地下水位迅速上升，在降雨持续3-5天后，地下水位上升趋于稳定，上升高度约为1.5-2米。孔隙水压力在滑坡体表层和浅层区域显著增加，形成了一定范围的孔隙水压力升高区。在滑坡体前缘和后缘等部位，孔隙水压力增加较为明显，这是因为这些部位更容易受到降雨入渗的影响，且排水条件相对较差。大雨（降雨强度为50mm/d）时，地下水位上升速度更快，在短时间内地下水位迅速上升。降雨持续2-3天后，地下水位上升高度达到2.5-3米。孔隙水压力在滑坡体大部分区域都有显著增加，孔隙水压力等值线向深部扩展，形成了较大范围的高孔隙水压力区。此时，滑坡体的渗流场发生了较大变化，渗流速度明显增大，尤其是在滑坡体的浅层和坡脚部位，渗流速度增加更为显著。暴雨（降雨强度为100mm/d）情况下，地下水位急剧上升，在降雨1-2天后，地下水位上升高度可达4-5米。孔隙水压力在滑坡体内部全面大幅增加，整个滑坡体几乎都处于高孔隙水压力状态。滑坡体的渗流场发生了剧烈变化，渗流速度急剧增大，在滑坡体的前缘和后缘等部位，渗流速度甚至超过了正常情况下的数倍。这种高孔隙水压力和大渗流速度的状态，极大地削弱了滑坡体的稳定性，增加了滑坡发生的风险。综上所述，降雨强度对吾隘滑坡渗流场有显著影响。随着降雨强度的增大，地下水位上升速度加快，上升高度增加，孔隙水压力增大的范围和幅度也相应增加。降雨强度的增大使得滑坡体内部的渗流场更加复杂，渗流速度增大，从而对滑坡的稳定性产生更大的威胁。5.2.4不同降雨历时吾隘滑坡渗流场分析在不同降雨历时条件下，吾隘滑坡渗流场呈现出明显的动态变化特征。当降雨历时为1d时，降雨入渗主要集中在滑坡体表层，地下水位上升幅度较小，仅在表层0-1米范围内地下水位有所上升，上升高度约为0.3-0.5米。孔隙水压力在滑坡体表层有所增加，但增加范围和幅度有限，主要集中在表层0-0.5米范围内，孔隙水压力增量在5-10kPa之间。此时，滑坡体的渗流场基本保持原有状态，渗流速度变化不大。降雨历时为3d时，地下水位逐渐向深部扩展，在滑坡体表层和浅层区域地下水位上升明显，上升高度约为1-1.5米。孔隙水压力在滑坡体表层和浅层区域显著增加，形成了一定范围的孔隙水压力升高区，孔隙水压力增量在10-20kPa之间。渗流速度在滑坡体表层和浅层区域有所增大，尤其是在坡脚部位，渗流速度增加较为明显。降雨历时达到5d时，地下水位进一步上升，在滑坡体大部分区域地下水位都有不同程度的升高，上升高度约为1.5-2.5米。孔隙水压力在滑坡体内部广泛增加，孔隙水压力等值线向深部延伸，形成了较大范围的高孔隙水压力区，孔隙水压力增量在20-30kPa之间。渗流速度在滑坡体内部全面增大，渗流路径更加复杂，在滑坡体的后缘和两侧等部位，渗流速度相对较大。当降雨历时为7d时，地下水位上升趋于稳定，但仍维持在较高水平，上升高度约为2.5-3米。孔隙水压力在滑坡体内部达到较高值，整个滑坡体几乎都处于高孔隙水压力状态，孔隙水压力增量在30-40kPa之间。渗流速度在滑坡体内部保持较大值，渗流场处于较为稳定但复杂的状态。由此可见，降雨历时对吾隘滑坡渗流场的影响显著。随着降雨历时的增加，地下水位逐渐上升并向深部扩展，孔隙水压力不断增大，渗流速度逐渐增大，渗流场变得更加复杂。较长的降雨历时使得滑坡体有更多的时间吸收雨水，导致渗流场的变化更加明显，从而对滑坡的稳定性产生更为不利的影响。5.2.5相同总降雨量下吾隘滑坡渗流场分析在相同总降雨量不同降雨组合的模拟中，设置了两种典型的降雨组合：组合一为降雨强度50mm/d、历时2d，总降雨量为100mm；组合二为降雨强度25mm/d、历时4d，总降雨量同样为100mm。对于组合一，在降雨初期，由于降雨强度较大，雨水迅速入渗，地下水位快速上升。在第1天，地下水位上升高度达到1.5-2米，滑坡体表层和浅层区域孔隙水压力急剧增加，孔隙水压力增量在20-30kPa之间。渗流速度在滑坡体表层和浅层区域显著增大，尤其是在坡脚部位，渗流速度增加更为明显。随着降雨的持续，在第2天，地下水位继续上升，上升高度达到2.5-3米，孔隙水压力在滑坡体大部分区域都有进一步增加，形成了较大范围的高孔隙水压力区。组合二的降雨过程相对平缓，在第1天，由于降雨强度较小，地下水位上升较为缓慢，上升高度约为0.5-1米，孔隙水压力在滑坡体表层有所增加，增量在5-10kPa之间。随着降雨的持续，在第2天和第3天，地下水位逐渐上升，上升高度分别达到1-1.5米和1.5-2米，孔隙水压力在滑坡体表层和浅层区域逐渐增大。到第4天，地下水位上升高度达到2-2.5米，孔隙水压力在滑坡体内部广泛增加，但增加幅度相对组合一较小。对比两种降雨组合，虽然总降雨量相同，但由于降雨强度和历时的不同，渗流场存在明显差异。组合一由于降雨强度大，在短时间内导致地下水位快速上升，孔隙水压力急剧增加，渗流速度迅速增大，对滑坡体的稳定性产生较大冲击。而组合二降雨强度较小，历时较长，地下水位和孔隙水压力的增加相对较为平缓，渗流场的变化相对较为缓和。这表明，在相同总降雨量下，降雨强度和历时的组合对吾隘滑坡渗流场有重要影响，高强度短历时的降雨更容易引起渗流场的剧烈变化，增加滑坡的不稳定因素。5.2.6一次降雨历程吾隘滑坡渗流场分析根据吾隘镇当地的气象监测数据，选取了一次典型的降雨过程进行模拟。该次降雨过程历时5天，降雨强度呈现出先增大后减小的变化趋势。降雨初期，降雨强度较小，为10-20mm/d，随着时间的推移，降雨强度逐渐增大，在第3天达到最大值50mm/d，随后降雨强度逐渐减小，在第5天降雨结束。在降雨初期，由于降雨强度较小，雨水入渗缓慢，地下水位上升较为平缓。在第1天，地下水位上升高度约为0.3-0.5米，滑坡体表层孔隙水压力有所增加，增量在5-10kPa之间。随着降雨强度的增大，在第2天和第3天，地下水位迅速上升，上升高度分别达到1-1.5米和2-2.5米。孔隙水压力在滑坡体表层和浅层区域显著增加，形成了较大范围的高孔隙水压力区，孔隙水压力增量在15-30kPa之间。渗流速度在滑坡体表层和浅层区域明显增大，尤其是在坡脚和后缘等部位，渗流速度增加更为显著。在降雨强度达到最大值后，随着降雨强度的减小，地下水位上升速度逐渐减缓。在第4天，地下水位上升高度约为2.5-3米，孔隙水压力在滑坡体内部继续增加，但增加幅度有所减小。到第5天降雨结束时，地下水位基本稳定，上升高度达到3-3.5米，孔隙水压力在滑坡体内部维持在较高水平。在降雨结束后的一段时间内，地下水位逐渐下降，但下降速度较为缓慢。孔隙水压力也逐渐减小，但在滑坡体内部仍存在一定范围的高孔隙水压力区。渗流速度逐渐减小，但在滑坡体的某些局部区域，由于排水不畅，渗流速度仍然相对较大。通过对这次降雨历程的模拟分析，可以看出吾隘滑坡渗流场在实际降雨过程中的动态变化十分复杂。降雨强度的变化导致地下水位、孔隙水压力和渗流速度等参数不断发生变化，在降雨强度较大的时段，渗流场的变化更为剧烈。这种复杂的渗流场变化对滑坡的稳定性产生了显著影响，在降雨强度较大的时期，滑坡的稳定性明显降低，增加了滑坡发生的风险。5.3本章小结本章深入研究了不同降雨条件下吾隘滑坡渗流场的变化规律，通过理论分析与数值模拟，全面揭示了降雨对滑坡渗流场的影响机制。在理论层面，详细阐述了饱和—非饱和渗流基本理论，明确了达西定律在饱和渗流中的应用以及Richards方程对非饱和渗流的准确描述，深入剖析了非饱和渗透系数与基质吸力的紧密关系。系统阐述了边坡降雨入渗过程的溅渗、下渗、再分布和蒸发排泄四个阶段，以及降雨入渗对滑坡渗流场中地下水位、孔隙水压力和渗流速度等关键参数的影响。在数值模拟方面，选用GeoStudio软件的SEEP/W模块，基于有限元方法，充分发挥其前处理功能强大、材料本构模型丰富、边界条件设置灵活以及计算结果可视化程度高等优势，成功建立了吾隘滑坡的数值模型。通过精心设计不同降雨强度、降雨历时、相同总降雨量以及一次降雨历程等多种模拟方案，全面系统地研究了降雨因素对滑坡渗流场的影响。模拟结果表明，降雨强度越大，地下水位上升速度越快，上升高度越高，孔隙水压力增大的范围和幅度也越大，滑坡体内部的渗流场更为复杂，渗流速度显著增大，对滑坡稳定性的威胁更为严重。随着降雨历时的增加，地下水位逐渐上升并向深部扩展，孔隙水压力不断增大，渗流速度逐渐增大，渗流场变得更加复杂，对滑坡稳定性的不利影响更为显著。在相同总降雨量下，高强度短历时的降雨更容易引发渗流场的剧烈变化，增加滑坡的不稳定因素。一次降雨历程的模拟显示，降雨强度的动态变化导致地下水位、孔隙水压力和渗流速度等参数不断改变，在降雨强度较大的时段，渗流场的变化更为剧烈，滑坡的稳定性明显降低。本章的研究成果为深入理解降雨对吾隘滑坡稳定性的影响提供了重要的渗流场数据支持，为后续滑坡稳定性分析奠定了坚实基础，也为滑坡的防治和预警提供了关键的理论依据和技术支持。六、不同降雨条件下吾隘滑坡稳定性评估6.1不同降雨条件下滑坡稳定性极限平衡分析6.1.1计算原理极限平衡法是基于静力平衡原理来分析滑坡稳定性的经典方法，其核心在于对滑坡体进行力学剖析，通过求解抗滑力与下滑力的比值，即稳定性系数，来判断滑坡的稳定状态。在众多极限平衡法中，瑞典条分法和Bishop法应用广泛。瑞典条分法由Fellenius提出，是最早的条分法。该方法将滑坡体沿滑动面划分为若干土条，假定土条间的条间力相互抵消，只考虑土条的重力、滑动面上的抗滑力和滑动力。通过对每个土条进行力矩平衡分析，以整个滑动土体对滑动圆心的力矩平衡为基础，建立稳定性系数计算公式。假设第i个土条的重力为W_i，其在滑动面上的法向分力为N_i=W_i\cos\alpha_i，切向分力为T_i=W_i\sin\alpha_i，其中\alpha_i为土条底面与水平面的夹角。滑动面上的抗剪强度\tau_i=c_i+\sigma_{i}\tan\varphi_i，其中c_i为黏聚力，\sigma_{i}为有效应力，\varphi_i为内摩擦角。稳定性系数F_s的计算公式为：F_s=\frac{\sum_{i=1}^{n}(c_il_i+N_i\tan\varphi_i)}{\sum_{i=1}^{n}T_i}其中l_i为第i个土条滑动面的长度。瑞典条分法概念简单，计算方便，但由于忽略了土条间的相互作用力，计算结果相对保守，安全系数偏低。Bishop法是在瑞典条分法的基础上发展而来，它考虑了土条间的水平作用力，在不考虑条块间切向力的前提下，满足力多边形闭合条件。假设滑动面是以圆心为O，半径为R的滑弧，从中任取一土条i为分离体，其周边作用力包括土条重W_i引起的切向力T_i和法向反力N_i，分别作用于底面中心处；土条侧面作用法向力E_i、E_{i+1}。根据静力平衡条件和极限平衡状态时各土条力对滑动圆心的力矩之和为零等条件，可得Bishop法求土坡稳定系数的普遍公式。在实际应用中，常采用Bishop简化式，其稳定性系数F_s的计算公式为：F_s=\frac{\sum_{i=1}^{n}\frac{1}{m_{\alphai}}(c_ib_i+(W_i-u_ib_i)\tan\varphi_i)}{\sum_{i=1}^{n}W_i\sin\alpha_i}其中b_i为土条宽度，u_i为孔隙水压力，m_{\alphai}=\cos\alpha_i+\frac{\sin\alpha_i\tan\varphi_i}{F_s}。Bishop法由于考虑了条块间的水平力作用，计算结果相对更接近实际情况，安全系数比瑞典条分法略高。6.1.2参数取值岩土体物理力学参数的准确取值是滑坡稳定性分析的关键。本研究通过多种方法获取吾隘滑坡岩土体的物理力学参数。室内试验是获取参数的重要手段之一。在实验室对采集的滑坡体和滑带土样品进行了一系列物理力学试验。采用环刀法测定岩土体的天然密度，通过烘干法确定含水量，利用直剪试验获取抗剪强度指标内摩擦角和黏聚力，通过固结试验得到压缩性指标等。例如，对滑坡体粉质黏土样品进行直剪试验，在垂直压力分别为100kPa、200kPa、300kPa的条件下，测得内摩擦角分别为[X1]°、[X2]°、[X3]°，黏聚力分别为[Y1]kPa、[Y2]kPa、[Y3]kPa。原位测试能够获取更接近实际工况的岩土参数。采用标准贯入试验（SPT），通过击入土层测量土层的抗力，间接推算出土的密实度和强度参数。进行静载荷试验，使用重物对土层施加静态荷载，测量沉降和刚度，分析土体的弹性模量等参数。在吾隘滑坡现场进行标准贯入试验，在不同位置测得的标准贯入击数为[具体击数范围]，根据相关经验公式，推算出土体的密实度和强度参数。同时，参考已有文献和经验公式，结合吾隘滑坡的地质条件和岩土体类型，对参数进行合理的验证和调整。例如，根据当地类似地质条件下的工程经验，对室内试验和原位测试得到的参数进行对比分析，确保参数取值的合理性和可靠性。综合考虑以上因素，最终确定吾隘滑坡岩土体的物理力学参数如下：滑坡体粉质黏土的天然密度为[X]g/cm³，含水量为[X]%，内摩擦角为[X]°，黏聚力为[X]kPa，压缩系数为[X]MPa⁻¹，压缩模量为[X]MPa；滑带土的天然密度为[X]g/cm³，含水量为[X]%，内摩擦角为[X]°，黏聚力为[X]kPa，压缩系数为[X]MPa⁻¹，压缩模量为[X]MPa。这些参数为后续的滑坡稳定性分析提供了准确的数据基础。6.1.3计算剖面与计算方案为了准确分析吾隘滑坡在不同降雨条件下的稳定性，选取了具有代表性的计算剖面。通过对滑坡的地形地貌、地质构造、岩土体分布等因素的综合分析，确定了沿滑坡主滑方向的典型剖面作为计算剖面。该剖面能够较好地反映滑坡体的结构特征和地质条件，包括滑坡体的厚度变化、滑带的位置和形态、滑床的产状等。在计算剖面上，清晰地显示出滑坡体主要由粉质黏土夹碎石组成，滑带位于滑坡体与滑床之间，滑床为砂岩、页岩互层。针对不同降雨条件，设计了多种计算方案。考虑不同降雨强度工况，设置小雨（降雨强度为10mm/d）、中雨（降雨强度为30mm/d）、大雨（降雨强度为50mm/d）、暴雨（降雨强度为100mm/d）等不同降雨强度，分别计算在这些降雨强度下滑坡的稳定性系数。针对不同降雨历时工况，设置降雨历时为1d、3d、5d、7d等，研究在不同降雨历时条件下滑坡稳定性的变化。在相同总降雨量工况下，通过不同的降雨强度和历时组合，使总降雨量相同，如降雨强度为50mm/d、历时2d和降雨强度为25mm/d、历时4d等，对比分析不同降雨组合下滑坡的稳定性。考虑一次降雨历程工况，根据当地实际的降雨过程数据，模拟一次完整的降雨过程中滑坡稳定性的动态变化。通过这些计算方案，可以全面系统地研究降雨条件对吾隘滑坡稳定性的影响。6.1.4不同降雨强度吾隘滑坡稳定性分析通过极限平衡法计算不同降雨强度下吾隘滑坡的稳定性系数，结果表明，降雨强度对滑坡稳定性有显著影响。在小雨（降雨强度为10mm/d）条件下，滑坡的稳定性系数为[F1]，处于基本稳定状态。这是因为小雨条件下，降雨入渗量较小，对滑坡体的含水量和孔隙水压力影响有限，滑坡体的抗滑力仍能抵抗下滑力。当中雨（降雨强度为30mm/d）发生时，稳定性系数下降至[F2]，滑坡稳定性有所降低。中雨使滑坡体的含水量增加，容重增大，下滑力增大。同时，降雨入渗导致孔隙水压力升高，有效应力降低，抗剪强度减小。这些因素共同作用，使得滑坡的稳定性下降。大雨（降雨强度为50mm/d）时，稳定性系数进一步下降至[F3]，滑坡处于欠稳定状态。随着降雨强度的增大，滑坡体的含水量和孔隙水压力显著增加，抗滑力大幅降低，下滑力明显增大，滑坡的稳定性受到严重威胁。暴雨（降雨强度为100mm/d）情况下，稳定性系数降至[F4]，小于1，滑坡处于不稳定状态。暴雨导致滑坡体几乎完全饱和，孔隙水压力急剧升高，抗剪强度急剧降低，下滑力远远超过抗滑力，滑坡极易发生滑动。综上所述，随着降雨强度的增大，吾隘滑坡的稳定性系数逐渐减小，滑坡稳定性逐渐降低。降雨强度的增加使得滑坡体的力学性质发生显著变化，增加了滑坡发生的风险。6.1.5不同降雨历时下吾隘滑坡稳定性分析在不同降雨历时条件下，吾隘滑坡的稳定性呈现出明显的变化规律。当降雨历时为1d时，滑坡的稳定性系数为[F5]，处于稳定状态。由于降雨历时较短，降雨入渗对滑坡体的影响范围和程度有限，滑坡体的力学性质变化较小，抗滑力能够维持滑坡的稳定。降雨历时为3d时，稳定性系数下降至[F6]，滑坡稳定性有所降低。随着降雨历时的增加，降雨入渗逐渐深入滑坡体内部，滑坡体的含水量和孔隙水压力逐渐增加，下滑力增大，抗剪强度减小，导致滑坡的稳定性下降。降雨历时达到5d时，稳定性系数进一步下降至[F7]，滑坡处于欠稳定状态。较长的降雨历时使得滑坡体有更多的时间吸收雨水，含水量和孔隙水压力持续增加，抗滑力不断降低，下滑力不断增大，滑坡的稳定性受到较大影响。当降雨历时为7d时，稳定性系数降至[F8]，小于1，滑坡处于不稳定状态。长时间的降雨使滑坡体处于饱和状态，孔隙水压力达到较高水平，抗剪强度大幅降低，下滑力远远超过抗滑力，滑坡发生滑动的可能性极大。由此可见，随着降雨历时的增加，吾隘滑坡的稳定性系数逐渐减小，滑坡稳定性逐渐降低。降雨历时的延长使得滑坡体的力学性质逐渐恶化，增加了滑坡发生的风险。6.1.6相同总降雨量下吾隘滑坡稳定性分析在相同总降雨量不同降雨组合的计算中，对比了两种典型的降雨组合。组合一为降雨强度50mm/d、历时2d，总降雨量为100mm；组合二为降雨强度25mm/d、历时4d，总降雨量同样为100mm。对于组合一，计算得到的稳定性系数为[F9]。由于降雨强度较大，在短时间内导致滑坡体含水量迅速增加，孔隙水压力急剧升高，抗滑力大幅降低，下滑力明显增大，对滑坡稳定性产生较大冲击。组合二的稳定性系数为[F10]。降雨强度较小，历时较长，滑坡体的含水量和孔隙水压力增加相对较为平缓，抗滑力的降低和下滑力的增大也相对较缓，滑坡稳定性的变化相对较小。对比两种降雨组合，虽然总降雨量相同，但由于降雨强度和历时的不同，滑坡的稳定性存在明显差异。组合一的稳定性系数明显小于组合二，表明高强度短历时的降雨更容易引起滑坡稳定性的急剧下降，增加滑坡的不稳定因素。这说明在相同总降雨量下，降雨强度和历时的组合对吾隘滑坡稳定性有重要影响。6.1.7一次降雨历程吾隘滑坡稳定性分析结合一次降雨历程的渗流场分析结果，对吾隘滑坡的稳定性进行评估。根据当地实际的降雨过程数据，该次降雨历时5天，降雨强度呈现出先增大后减小的变化趋势。在降雨初期，由于降雨强度较小，滑坡的稳定性系数为[F11]，处于稳定状态。随着降雨强度的增大，在第3天达到最大值50mm/d时，稳定性系数下降至[F12]，滑坡处于欠稳定状态。此时，降雨入渗导致滑坡体的含水量和孔隙水压