四川旺苍县杜家营滑坡工程地质特征剖析与治理方案的科学构建

四川旺苍县杜家营滑坡工程地质特征剖析与治理方案的科学构建一、引言1.1研究背景与意义滑坡作为一种常见且危害巨大的地质灾害，常常给人类的生命财产安全和生态环境带来严重的威胁。四川旺苍县杜家营滑坡位于广元市旺苍县东部五权镇北端约200m处，其前缘有场镇公路经过，公路东侧紧邻南江河一级支流五郎河。该滑坡的变形历史可追溯至20世纪70年代中期，此后每逢雨季，便会出现不同程度的地面开裂等现象。杜家营滑坡对当地居民的生命财产安全构成了直接威胁。滑坡前缘居住着43户共179人，一旦滑坡发生，这些居民将首当其冲，极有可能遭受房屋损毁、人员伤亡等严重后果。滑坡前缘与五郎河距离极近，这使得滑坡具有引发河道阻塞的潜在风险。一旦滑坡体冲入河道，堵塞河流，将会导致河水水位迅速上涨，进而引发洪水灾害，间接威胁到五权场镇2500余居民的生命财产安全。从经济角度来看，该滑坡潜在经济损失可达5000万元以上，这不仅包括居民房屋、财产的损失，还涵盖了基础设施损坏、农业生产受影响以及后续救援、重建等方面的费用。滑坡灾害的发生，还会对当地的生态环境造成极大的破坏。它可能导致植被被掩埋、土壤侵蚀加剧、生物栖息地遭到破坏，进而影响区域生态平衡。此外，滑坡还可能引发一系列的次生地质灾害，如泥石流等，进一步加剧对生态环境和人类社会的危害。研究杜家营滑坡的工程地质特征，能够深入了解滑坡的形成机制、影响因素以及演化过程。通过掌握这些信息，我们可以更为准确地评估滑坡的稳定性，预测其发展趋势，为制定科学合理的治理方案提供坚实的数据支持和理论依据。而制定有效的治理方案，对于保障当地居民的生命财产安全、保护生态环境以及促进区域社会经济的可持续发展具有至关重要的意义。它可以从根本上消除或减轻滑坡灾害的威胁，避免人员伤亡和财产损失，维护社会的稳定与和谐。因此，对四川旺苍县杜家营滑坡工程地质特征及治理方案的研究迫在眉睫，具有极高的现实价值和重要的社会意义。1.2国内外研究现状国外对滑坡的研究起步较早，从20世纪初开始，许多发达国家就对滑坡问题进行了深入研究，在滑坡监测预警、风险评估和防治措施等方面取得了显著成果。在滑坡监测预警方面，不少研究者利用地球物理学、地质学和计算机科学等技术手段，建立了许多有效的监测预警系统。如美国采用分布式光纤传感技术，对科罗拉多州一处大型滑坡进行了实时监测，成功预测了滑坡发生的时间和位置。在风险评估方面，意大利和法国等国家制定了详细的地质灾害风险评估指南和规范，明确了评估流程和方法。一些研究者还提出了基于GIS技术的滑坡危险性评估模型，如“CLUE-S”模型等，这些模型综合考虑了地形、地质、气象等因素，可对滑坡危险性进行快速评估。在防治措施方面，国外的滑坡防治主要包括工程措施和非工程措施两大类。工程措施主要包括加固不稳定岩土体、修筑挡土墙和排水设施等；非工程措施则包括植被护坡、土地利用规划和管理、公众教育和紧急应对计划等，这些措施均在不同程度上取得了成功，为当地减缓滑坡灾害的发生提供了有力支持。我国是一个地质灾害多发的国家，滑坡灾害的防治工作一直受到高度重视。20世纪80年代以来，我国在滑坡防治方面取得了显著进展。在监测预警方面，我国研究者利用遥感技术、自动监测设备和数值模拟等方法，对多处大型滑坡进行了有效监测。例如，2008年汶川地震后，中国地质调查局在四川绵竹市开展了基于遥感的滑坡监测预警示范区建设，为灾后重建提供了科学依据。在风险评估方面，我国研究者结合遥感技术、GIS和数值模拟等方法，针对不同地区和不同类型的滑坡灾害，建立了相应的风险评估模型和规范。例如，中国科学院成都山地所提出了基于GIS的西南山区滑坡危险性评估方法，为该地区的防灾减灾工作提供了有力支持。在防治措施方面，我国研究者针对不同类型和规模的滑坡灾害，提出了多种有效的防治措施。对于小型滑坡，可采用排水、削坡和加固等工程措施；对于大型滑坡，可能需要采用综合性防治措施，如植被护坡、土地利用规划和管理以及公众教育和紧急应对计划等。尽管国内外学者在滑坡防治方面已取得显著进展，但仍存在以下不足之处：现有研究多针对特定地区或特定类型的滑坡灾害，缺乏普适性的滑坡防治方案，不同地区的地质条件、气候因素、地形地貌等存在差异，使得某一地区成功的防治方案难以直接应用到其他地区；现有研究多从单一角度出发，如监测预警或风险评估，缺乏对滑坡全过程的综合研究，滑坡的形成、发展到最终失稳是一个复杂的过程，涉及多个因素和环节，单一角度的研究难以全面揭示其本质和规律；现有研究多侧重于技术手段的开发和应用，缺乏对滑坡防治的社会经济和政策方面的研究，滑坡防治不仅是技术问题，还涉及社会经济发展、土地利用规划、政策法规等多个方面；现有研究多集中于纯粹的理论分析和数值模拟，缺乏足够的实地试验和验证，理论和模拟结果在实际应用中可能会受到各种因素的影响，需要通过实地试验来检验其有效性和可靠性。四川旺苍县杜家营滑坡具有独特的地质条件和环境背景，以往的研究成果难以直接应用于该滑坡的治理。因此，有必要针对杜家营滑坡的具体情况，深入研究其工程地质特征，综合考虑各种因素，提出适合该滑坡的治理方案，以填补在该特定区域滑坡研究方面的空白，为当地的滑坡防治工作提供科学依据和实践指导。1.3研究内容与方法本研究将深入剖析四川旺苍县杜家营滑坡的工程地质特征，全面分析其稳定性，并制定科学合理的治理方案，具体研究内容如下：滑坡工程地质特征研究：通过详细的野外实地勘查，精准掌握滑坡的平面形态、坡面特征、地形坡度等基本情况。同时，结合勘探资料、降雨数据、室内试验和野外调查数据，深入探究滑体、滑带以及滑床的物质组成，从而深入分析该滑坡的滑坡机理，并经过物理力学试验得出各岩土体的物理力学性质相关数据，为后续研究提供坚实基础。滑坡稳定性分析：运用工程地质法对滑坡稳定性进行定性分析，从地质构造、岩土体特性、地形地貌等方面综合判断滑坡的稳定状态。利用极限平衡法进行定量分析，通过计算滑坡体的下滑力和抗滑力，确定滑坡的稳定系数，精确评估滑坡的稳定性程度，为治理方案的制定提供科学依据。滑坡治理方案研究：根据滑坡的工程地质特征和稳定性分析结果，提出多种可行的治理方案。对每种方案进行详细的技术经济分析，综合考虑工程的安全性、可靠性、经济性以及对环境的影响等因素，筛选出最优的治理方案。对最优方案进行详细的设计计算，确定抗滑桩、挡土墙以及截水沟等工程结构的具体参数，并对其力学稳定性进行严格计算和验证，确保治理方案的可行性和有效性。为实现上述研究目标，本研究将采用以下研究方法：实地勘查法：对杜家营滑坡进行全面、细致的野外实地勘查，包括地质测绘、地形测量、水文地质调查等，详细记录滑坡的各种特征和相关信息，获取第一手资料。通过实地勘查，直观了解滑坡的现状和周边环境，为后续的分析和研究提供基础数据。试验分析法：采集滑体、滑带和滑床的岩土样本，在室内进行物理力学性质试验，如密度试验、含水率试验、剪切强度试验等，准确测定岩土体的各项物理力学参数。分析降雨数据，研究降雨与滑坡变形之间的关系，为滑坡机理分析和稳定性评价提供科学依据。数值模拟法：利用专业的岩土工程数值模拟软件，如PLAXIS2D等，建立杜家营滑坡的数值模型。通过模拟不同工况下滑坡的变形和稳定性情况，分析治理方案的效果，对治理方案进行优化和验证。数值模拟可以直观展示滑坡的演化过程和治理效果，为决策提供可视化的参考依据。理论分析法：运用工程地质、土力学、岩石力学等相关理论，对滑坡的形成机制、稳定性评价和治理方案进行深入分析。结合极限平衡理论、岩土体本构模型等，进行滑坡稳定性的定量计算和分析，确保研究结果的科学性和可靠性。二、杜家营滑坡区域地质背景2.1地形地貌杜家营滑坡处于四川盆地北缘的低山丘陵地带，该区域整体地势呈现出北高南低的态势。区域内山体走势较为复杂，山峦起伏，山脉纵横交错，山体的走向多为东北-西南向。从海拔高度来看，杜家营滑坡所在区域的海拔范围大致在500-800米之间，相对高差可达300米左右。这种相对较大的高差使得地形起伏较为明显，为滑坡的形成提供了一定的地形条件。在滑坡区域，地形坡度变化较大，滑坡中后部地形坡度相对较缓，约为20°，呈阶梯状分布，阶坎高度大约在1-4米之间。而滑坡中前部地形坡度则较为陡峭，坡度在20°-30°之间。这种上缓下陡的坡面形态，使得滑坡体上部的岩土体在重力作用下更容易产生向下的滑动趋势，增加了滑坡发生的可能性。滑坡平面形态呈现出典型的“圈椅”状，这是滑坡地貌的一种常见形态特征。“圈椅”状的平面形态表明滑坡体在滑动过程中，后部受到的拉力较大，导致岩土体向后退缩，形成了类似圈椅靠背的滑坡后壁；而两侧则由于滑动过程中的侧向挤压和剪切作用，形成了类似圈椅扶手的滑坡侧壁；滑坡前缘则向前突出，形如圈椅的前沿，这种形态特征在一定程度上反映了滑坡的运动方式和受力状态。坡面总体呈凸坡面，这种凸形坡面使得雨水在坡面上的径流速度较快，容易对坡面岩土体产生较强的冲刷作用，进一步破坏岩土体的结构，降低其抗滑能力。同时，凸形坡面在重力作用下，岩土体的稳定性相对较差，更容易发生变形和滑动。滑坡区域内沟谷发育，这些沟谷多为季节性流水通道，在雨季时水流较大，对沟谷两侧的坡脚产生强烈的冲刷和侵蚀作用，削弱了坡脚的支撑力，从而引发滑坡。此外，沟谷的存在还使得地下水的径流条件发生变化，地下水在沟谷附近容易汇聚，进一步软化岩土体，增加了滑坡发生的风险。该区域的地形地貌条件，如相对高差、坡度、坡面形态以及沟谷发育等因素，相互作用，共同影响着杜家营滑坡的形成和发展，为后续对滑坡工程地质特征和稳定性的分析提供了重要的地形基础。2.2地层岩性杜家营滑坡区域出露的地层岩性主要为侏罗系中统上沙溪庙组（J2s）砂岩、泥岩互层，以及第四系全新统滑坡堆积层（Q4del）、残坡积层（Q4el+dl）。侏罗系中统上沙溪庙组（J2s）广泛分布于滑坡区域，是滑坡体的主要组成部分。该地层岩性主要为一套陆相碎屑沉积岩，以砂岩和泥岩互层为特征。砂岩呈灰白色、浅灰色，主要由石英、长石等矿物组成，颗粒结构，中-细粒为主，分选性和磨圆度中等，钙质胶结，岩石的抗压强度较高，一般在30-50MPa之间，但由于受风化作用和构造运动的影响，岩石的完整性较差，节理裂隙较为发育。泥岩呈紫红色、灰绿色，主要由黏土矿物组成，具页理构造，岩石的抗压强度较低，一般在5-15MPa之间，遇水易软化、崩解，抗风化能力弱。砂岩和泥岩互层的结构特点使得岩体的力学性质差异较大，在风化、降雨等因素的作用下，泥岩易被侵蚀、软化，形成相对软弱的结构面，为滑坡的形成提供了潜在的滑动面。第四系全新统滑坡堆积层（Q4del）主要分布于滑坡体上，是滑坡发生后形成的堆积物。其物质组成主要为砂岩、泥岩的碎块石以及粉质黏土，碎块石含量约占40%-60%，粒径大小不一，一般在0.2-2m之间，大者可达5m以上，呈棱角状-次棱角状，分选性差，无明显层理，结构松散。粉质黏土充填于碎块石之间，其含水率较高，一般在20%-30%之间，塑性指数为15-20，具有一定的黏聚力和可塑性，但强度较低。滑坡堆积层的厚度变化较大，在滑坡后缘较薄，一般为1-3m，向前缘逐渐增厚，最厚处可达10m左右。该堆积层的工程性质较差，在自重和外部荷载作用下，容易产生变形和滑动。第四系全新统残坡积层（Q4el+dl）主要分布于滑坡体表层及周边斜坡地带，是基岩经风化、剥蚀后，在重力和坡面水流作用下堆积形成的。其物质组成主要为粉质黏土夹少量碎石，碎石含量约占10%-30%，粒径一般小于0.5m，呈次棱角状-次圆状，分选性和磨圆度较好，结构较松散。粉质黏土的含水率一般在15%-25%之间，塑性指数为12-18，黏聚力和内摩擦角相对较小。残坡积层的厚度一般为0.5-2m，局部地段可达3m左右。该层与下伏基岩呈不整合接触，由于其结构松散、强度较低，在降雨、地震等因素的作用下，容易发生滑动，进而引发滑坡。不同地层岩性的物理力学性质差异显著，对滑坡的形成和发展产生了重要影响。侏罗系中统上沙溪庙组砂岩、泥岩互层的结构特点，形成了潜在的滑动面；第四系全新统滑坡堆积层和残坡积层的工程性质较差，强度低、结构松散，在外部因素的作用下，容易产生变形和滑动，为滑坡的发生提供了物质基础。这些地层岩性特征是分析杜家营滑坡工程地质特征和稳定性的重要依据，也为后续治理方案的制定提供了关键的参考信息。2.3地质构造杜家营滑坡区域位于扬子准地台西缘龙门山褶断带与米仓山台穹的过渡地带，区域地质构造较为复杂，经历了多期构造运动，主要包括加里东运动、海西运动、印支运动和燕山运动，这些构造运动对区域内地层的沉积、褶皱、断裂以及岩体的力学性质产生了深远影响。区内褶皱构造较为发育，主要为一系列轴向近南北向的紧闭褶皱，褶皱形态以箱状褶皱和尖棱褶皱为主。这些褶皱使得地层发生强烈的弯曲变形，岩层产状变化较大。在滑坡区域附近，岩层倾向多为东-南东向，倾角一般在30°-50°之间。褶皱的核部和转折端部位，由于岩石受到强烈的挤压和拉伸作用，节理裂隙极为发育，岩石破碎，完整性遭到严重破坏，力学强度大幅降低，为滑坡的形成提供了有利的地质结构条件。断层构造在该区域也较为常见，主要有北东向和北西向两组断层。北东向断层规模较大，延伸长度可达数千米，断层走向一般为北东30°-60°，倾向南东，倾角在60°-80°之间。这些断层多为逆断层或逆冲断层，具有较强的活动性，在地质历史时期曾多次发生构造活动，导致断层两侧的岩体破碎，形成了宽度不等的断层破碎带。北西向断层规模相对较小，延伸长度一般在数百米至千米左右，断层走向为北西300°-330°，倾向北东，倾角在50°-70°之间。这些断层与北东向断层相互切割、错动，进一步破坏了岩体的完整性，使得岩体的结构更加复杂。构造运动对杜家营滑坡的控制作用主要体现在以下几个方面：褶皱和断层构造破坏了岩体的完整性，形成了大量的节理裂隙，为地下水的运移和储存提供了良好的通道和空间。地下水在这些裂隙中流动，不断软化和侵蚀岩土体，降低了岩土体的抗剪强度，增加了滑坡发生的可能性。褶皱和断层造成的地层岩性变化和岩体结构差异，使得斜坡岩土体的力学性质不均一，在重力作用下，容易产生应力集中现象，当应力超过岩土体的强度极限时，就会引发滑坡。构造运动产生的地应力场改变了斜坡的应力状态，使斜坡岩土体处于复杂的应力环境中，增加了斜坡的不稳定性。在区域构造应力的作用下，斜坡岩土体可能发生变形和破裂，从而为滑坡的发生创造了条件。杜家营滑坡区域的地质构造特征，如褶皱、断层的分布与特征，对滑坡的形成和发展具有重要的控制作用。这些构造因素与地形地貌、地层岩性等因素相互作用，共同影响着滑坡的稳定性，是研究杜家营滑坡工程地质特征和制定治理方案时不可忽视的重要因素。2.4气象与水文条件2.4.1气象条件杜家营滑坡所在的四川旺苍县属于亚热带湿润季风气候区，气候温暖湿润，四季分明。该区域多年平均气温为16.2℃，其中1月平均气温最低，约为5.1℃；7月平均气温最高，可达26.6℃。年温差相对较大，这使得岩土体在温度的反复变化作用下，容易产生胀缩效应，从而导致岩石的结构逐渐破坏，节理裂隙进一步发育，为滑坡的发生创造了一定的条件。在降水方面，旺苍县年平均降水量丰富，约为1200mm，降水主要集中在5-9月，这期间的降水量约占全年降水量的75%。降水集中且强度较大，容易形成地表径流，对坡面产生强烈的冲刷作用。大量的雨水渗入地下，使地下水位迅速上升，增加了岩土体的重量，同时还会软化岩土体，降低其抗剪强度，极大地增加了滑坡发生的可能性。例如，2018年7月，旺苍县遭遇了连续强降雨天气，降雨量在短时间内超过了300mm，导致杜家营滑坡区域出现了多处地面开裂和小规模的滑坡现象。强降雨等极端气象条件与滑坡发生具有显著的相关性。当降雨量达到一定阈值时，滑坡发生的概率会急剧增加。通过对历史降雨数据和滑坡事件的统计分析发现，当连续5日降雨量超过150mm，且日最大降雨量超过50mm时，杜家营滑坡区域发生滑坡的可能性明显增大。强降雨引发滑坡的机制主要在于，大量雨水迅速渗入地下，在滑带处形成超孔隙水压力，有效应力降低，使得滑带土的抗剪强度大幅下降。强降雨还会导致坡面径流增大，对坡体产生冲刷和侵蚀作用，削弱坡体的稳定性。年平均蒸发量约为1000mm，主要集中在气温较高的夏季。蒸发作用会使得岩土体中的水分散失，导致岩土体收缩干裂，增加了雨水的入渗通道，间接影响滑坡的稳定性。相对湿度较高，年平均相对湿度在75%左右，这使得空气较为潮湿，岩土体长期处于湿润状态，进一步降低了其强度。杜家营滑坡区域的气象条件，如气温变化、降水分布以及蒸发和相对湿度等因素，相互作用，对滑坡的形成和发展产生了重要影响。尤其是强降雨等极端气象条件，是触发滑坡的关键因素之一，在研究滑坡工程地质特征和制定治理方案时，必须充分考虑气象条件的影响。2.4.2水文条件杜家营滑坡区域的地表水主要来源于大气降水形成的坡面径流以及五郎河。五郎河作为南江河的一级支流，从滑坡前缘东侧流过，河水流量随季节变化明显。在雨季，由于降水丰富，五郎河水位迅速上涨，流量增大，对滑坡前缘坡脚产生强烈的冲刷和侵蚀作用。长期的冲刷会导致坡脚土体被掏空，削弱坡体的支撑力，使坡体上部处于临空状态，从而引发滑坡。在2015年雨季，五郎河水位暴涨，对杜家营滑坡前缘坡脚的冲刷作用加剧，导致坡脚局部坍塌，滑坡体出现了明显的变形迹象。滑坡区域内地表水系发育，冲沟众多，这些冲沟呈树枝状分布，将坡面分割成多个小块。在降雨过程中，坡面径流会沿着冲沟迅速汇集，形成较大的水流，对冲沟两侧的坡体产生强烈的冲刷作用，破坏坡体的稳定性。冲沟的存在还使得地表水更容易渗入地下，增加了地下水的补给量，进一步影响滑坡的稳定性。地下水主要赋存于第四系松散堆积层和侏罗系中统上沙溪庙组砂岩、泥岩互层中。第四系松散堆积层透水性较好，地下水在其中的径流速度相对较快，水位变化受降水影响明显。在雨季，随着降水量的增加，地下水水位迅速上升；而在旱季，地下水水位则会逐渐下降。侏罗系中统上沙溪庙组砂岩、泥岩互层中，砂岩透水性相对较好，为含水层；泥岩透水性差，为相对隔水层。这种地层结构导致地下水在砂岩中富集，并在泥岩顶面形成相对隔水边界，使得地下水在泥岩顶面附近汇聚，软化泥岩，降低其抗剪强度，形成潜在的滑动面。地下水的水位变化对滑坡稳定性有着重要影响。当地下水水位上升时，会增加岩土体的重量，使下滑力增大；同时，地下水对滑带土的软化作用会降低滑带土的抗剪强度，减小抗滑力，从而降低滑坡的稳定性。据监测数据显示，在2019年雨季，杜家营滑坡区域地下水水位上升了2-3m，滑坡体的位移量明显增加，稳定性系数降低了0.1-0.2。地下水的径流条件也会影响滑坡的稳定性。当地下水径流不畅时，会在坡体中积聚，增加孔隙水压力，导致有效应力减小，削弱岩土体的强度。滑坡区域内存在一些相对隔水层和局部地形低洼处，容易造成地下水的汇聚和径流不畅，增加了滑坡发生的风险。杜家营滑坡区域的水文条件，包括地表水和地下水的分布、水位变化以及径流条件等，对滑坡的稳定性产生了多方面的影响。地表水的冲刷和侵蚀作用以及地下水的软化、增重和改变应力状态等作用，是导致滑坡发生和发展的重要因素，在研究和治理滑坡时，必须充分考虑水文条件的作用，采取有效的排水等防治措施，以提高滑坡的稳定性。三、杜家营滑坡工程地质特征3.1滑坡基本特征3.1.1滑坡形态杜家营滑坡平面形态呈现出典型的“圈椅”状，这种形态是滑坡在滑动过程中，由于后部岩土体受到拉力，两侧受到侧向挤压和剪切，前缘向前突出而形成的。从整体上看，滑坡的后缘较为圆滑，宛如圈椅的靠背，两侧边界清晰，似圈椅的扶手，前缘则较为开阔，如同圈椅的前沿。这种形态特征不仅反映了滑坡的运动方式，还在一定程度上揭示了滑坡的受力状态。滑坡的长度约为350m，宽度约为200m，由此可计算出滑坡的面积大约为70000m²。滑坡体的厚度变化较为明显，在滑坡后缘厚度相对较薄，约为2-3m，这是因为后缘岩土体在滑动过程中受到的拉伸作用相对较小，堆积相对较少。而向前缘方向，滑坡体厚度逐渐增加，在滑坡前缘最厚处可达12m左右，这是由于前缘岩土体在滑动过程中受到的挤压和堆积作用较强。坡面总体呈现出凸坡面的形态，这种凸形坡面使得雨水在坡面上的径流速度较快。当降雨发生时，大量雨水迅速汇聚成地表径流，对坡面岩土体产生强烈的冲刷作用。在长期的冲刷作用下，坡面岩土体的结构逐渐被破坏，颗粒之间的连接变得松散，抗滑能力大幅降低。坡面还具有上缓下陡的特征，滑坡中后部地形坡度相对较缓，约为20°，且呈阶梯状分布，阶坎高度大约在1-4m之间。这种阶梯状地形是由于滑坡体在滑动过程中，不同部位的运动速度和位移量存在差异所导致的。中前部地形坡度则较为陡峭，坡度在20°-30°之间，这种陡峭的坡度使得岩土体在重力作用下更容易产生向下的滑动趋势，增加了滑坡发生的风险。3.1.2滑坡边界杜家营滑坡的周界范围较为清晰，通过实地勘查和地形测绘等手段，可以准确地确定滑坡体与周围稳定岩土体的分界线。滑坡周界的确定对于评估滑坡的影响范围和制定治理方案具有重要意义。滑坡后壁位于滑坡的后缘，是滑坡体与后部稳定岩土体脱离开的分界面，在平面上多呈圈椅状。滑坡后壁的高度一般在3-5m之间，壁面上可以观察到明显的顺坡向擦痕，这些擦痕是滑坡体在滑动过程中与后壁岩土体相互摩擦产生的，它们记录了滑坡滑动的方向和过程。滑坡台阶是由于各段滑体运动速度的差异而在滑体上部形成的滑坡错台。在杜家营滑坡中，滑坡台阶较为明显，分布在滑坡体的中后部。这些台阶的宽度和高度各不相同，宽度一般在5-10m之间，高度在1-3m之间。滑坡台阶的存在反映了滑坡体在滑动过程中的阶段性和不均匀性。滑坡舌又称滑坡前缘或滑坡头，位于滑坡的前部，形如舌状伸入沟谷或河流，甚至有可能越过河对岸。杜家营滑坡的滑坡舌较为突出，向前延伸至五郎河附近，对河道的稳定性构成了潜在威胁。滑坡舌部位的岩土体较为松散，在河水的冲刷和浸泡作用下，容易发生坍塌和滑动，进一步加剧滑坡的发展。滑坡体两侧多以沟谷或裂缝为界。在实地勘查中可以发现，滑坡体两侧发育有明显的冲沟，这些冲沟是地表水长期冲刷形成的。冲沟的存在不仅分割了滑坡体与周围稳定岩土体，还为地表水的汇聚和下渗提供了通道，进一步影响了滑坡体的稳定性。滑坡体两侧还存在一些裂缝，这些裂缝是由于滑坡体在滑动过程中受到侧向拉力和剪切力作用而产生的。裂缝的宽度和深度不一，宽度一般在0.1-0.5m之间，深度可达数米。裂缝的存在使得滑坡体的整体性遭到破坏，增加了滑坡发生的可能性。3.2滑体、滑带及滑床特征3.2.1滑体物质组成滑体物质的组成对滑坡的性质和稳定性有着重要影响。通过对杜家营滑坡的勘探和试验分析，滑体主要由第四系全新统滑坡堆积层（Q4del）和残坡积层（Q4el+dl）组成。第四系全新统滑坡堆积层（Q4del）的物质成分较为复杂，主要包括砂岩、泥岩的碎块石以及粉质黏土。碎块石含量约占40%-60%，其粒径大小不一，呈现出较大的离散性。小粒径的碎块石一般在0.2-0.5m之间，大粒径的碎块石可达2m以上，甚至有部分超过5m。这些碎块石的形状多为棱角状-次棱角状，分选性差，表明其在搬运和堆积过程中缺乏有效的分选作用，是在较为复杂的地质作用下快速堆积形成的。这种形状和分选特征使得碎块石之间的排列较为松散，相互之间的嵌锁作用较弱，从而影响了滑体的整体结构稳定性。粉质黏土充填于碎块石之间，其含水率较高，一般在20%-30%之间。较高的含水率使得粉质黏土具有一定的可塑性，但也降低了其强度。通过室内试验测定，粉质黏土的塑性指数为15-20，具有一定的黏聚力，这在一定程度上能够约束碎块石的相对运动，增强滑体的整体性。然而，由于粉质黏土本身强度较低，在受到较大外力作用或地下水的长期浸泡软化时，其黏聚力和抗剪强度会显著降低，导致滑体的稳定性下降。第四系全新统残坡积层（Q4el+dl）主要由粉质黏土夹少量碎石组成。碎石含量约占10%-30%，粒径一般小于0.5m，相对较小。这些碎石呈次棱角状-次圆状，分选性和磨圆度较好，说明其在搬运过程中受到了一定程度的磨蚀和分选作用，经历了相对较长距离的搬运和堆积过程。这种特征使得残坡积层的结构相对较为紧密，但由于碎石含量较少，其整体强度仍然相对较低。粉质黏土的含水率一般在15%-25%之间，塑性指数为12-18，黏聚力和内摩擦角相对较小。在降雨等因素的作用下，含水率的增加会进一步降低其抗剪强度，使得残坡积层在坡面上容易发生滑动。通过物理力学试验，得到滑体的一些关键物理力学性质数据。滑体的天然密度一般在1.8-2.0g/cm³之间，这反映了滑体物质的紧密程度和质量分布情况。天然重度在18-20kN/m³之间，重度的大小直接影响着滑体在重力作用下的下滑力大小。内摩擦角在18°-22°之间，内摩擦角是衡量滑体抗剪强度的重要指标之一，其大小决定了滑体在受到剪切力时抵抗滑动的能力。黏聚力在10-15kPa之间，黏聚力则体现了滑体内部颗粒之间的连接强度，对滑体的整体稳定性也有着重要影响。这些物理力学性质数据为后续分析滑坡的稳定性和制定治理方案提供了重要的基础数据。它们反映了滑体物质的组成、结构和力学特性，对于理解滑坡的变形机制和预测滑坡的发展趋势具有关键作用。例如，内摩擦角和黏聚力的大小直接影响着滑坡的抗滑力，而滑体的密度和重度则决定了下滑力的大小，通过对这些参数的分析，可以评估滑坡在不同工况下的稳定性状态，为采取有效的防治措施提供科学依据。3.2.2滑带特征滑带是滑坡体与滑床之间的一个特殊部位，其特征对于滑坡的发生和发展起着至关重要的作用。杜家营滑坡的滑带位于滑体与下伏基岩之间，是滑坡滑动的主要控制面。滑带的位置相对稳定，通过勘探和地质测绘确定，其深度在滑坡后缘较浅，一般为2-3m，向前缘逐渐加深，在滑坡前缘最深处可达10-12m。这种深度变化与滑坡体的厚度变化趋势一致，反映了滑带的形成与滑坡体的变形过程密切相关。在滑坡体的滑动过程中，滑带受到了强烈的剪切作用，其物质组成和结构特征发生了显著变化。滑带的厚度一般在0.5-1.5m之间，相对较薄，但却是滑坡稳定性的关键控制因素。滑带物质主要由粉质黏土和少量碎块石组成，其中粉质黏土含量较高，约占70%-80%。这些粉质黏土颗粒细小，在长期的剪切作用下，颗粒定向排列明显，形成了独特的结构特征。碎块石含量较少，一般占20%-30%，粒径相对较小，多在0.1-0.3m之间，它们分散在粉质黏土中，进一步影响了滑带的力学性质。滑带土的结构特征表现为明显的定向排列和揉皱现象。在显微镜下观察，可见粉质黏土颗粒沿着滑动方向呈定向排列，形成了类似层理的结构。这种定向排列使得滑带土在平行于滑动方向上的强度较低，容易发生剪切破坏。滑带土中还存在大量的揉皱和挤压痕迹，这是由于滑带在滑动过程中受到强烈的剪切和挤压作用所导致的。这些揉皱和挤压痕迹进一步破坏了滑带土的结构完整性，降低了其抗剪强度。滑带土的抗剪强度是衡量滑坡稳定性的关键指标之一。通过室内直剪试验和三轴剪切试验测定，滑带土的内摩擦角在10°-15°之间，黏聚力在5-10kPa之间。与滑体相比，滑带土的内摩擦角和黏聚力明显较低，这是因为滑带在滑动过程中经历了强烈的剪切变形，土体结构被破坏，颗粒之间的连接强度降低。较低的抗剪强度使得滑带成为滑坡体中最薄弱的环节，在外部因素（如降雨、地震等）的作用下，滑带土的抗剪强度进一步降低，容易引发滑坡的滑动。滑带土的含水率一般在25%-35%之间，相对较高。高含水率使得滑带土处于饱水状态，颗粒之间的摩擦力减小，抗剪强度降低。地下水在滑带中的存在还会产生孔隙水压力，进一步削弱滑带土的有效应力，增加滑坡的不稳定性。当孔隙水压力达到一定程度时，滑带土可能会发生液化现象，导致滑坡迅速滑动。滑带的特征，包括位置、厚度、物质组成、结构特征以及抗剪强度和含水率等，对杜家营滑坡的稳定性产生了决定性影响。深入研究滑带特征，对于理解滑坡的形成机制和制定有效的治理方案具有重要意义。在治理方案的设计中，应充分考虑滑带的薄弱性，采取相应的措施提高滑带的抗剪强度，如加固滑带土、降低地下水位等，以增强滑坡的稳定性。3.2.3滑床特征滑床作为滑坡体滑动的依托面，其特征对于滑坡的稳定性和发展具有重要影响。杜家营滑坡的滑床主要由侏罗系中统上沙溪庙组（J2s）砂岩、泥岩互层组成。滑床的岩石类型主要为砂岩和泥岩，砂岩呈灰白色、浅灰色，主要矿物成分为石英、长石，中-细粒结构，颗粒分选性和磨圆度中等，钙质胶结。砂岩具有较高的抗压强度，一般在30-50MPa之间，这使得滑床在一定程度上能够承受滑坡体的压力。然而，由于受到风化作用和构造运动的影响，砂岩的完整性较差，节理裂隙较为发育。这些节理裂隙的存在降低了砂岩的强度和整体性，为地下水的运移提供了通道，同时也增加了滑坡体与滑床之间的摩擦力，影响了滑坡的滑动过程。泥岩呈紫红色、灰绿色，主要由黏土矿物组成，具页理构造。泥岩的抗压强度较低，一般在5-15MPa之间，遇水易软化、崩解，抗风化能力弱。在滑坡区域，泥岩往往作为相对软弱的夹层存在于砂岩之间，这种软硬相间的地层结构使得岩体的力学性质差异较大。在风化、降雨等因素的作用下，泥岩容易被侵蚀、软化，形成相对软弱的结构面，成为滑坡滑动的潜在控制面。滑床的岩体结构较为复杂，由于受到多期构造运动的影响，岩层产状变化较大。在滑坡区域，岩层倾向多为东-南东向，倾角一般在30°-50°之间。这种倾斜的岩层结构使得滑坡体在重力作用下更容易沿着滑床发生滑动。在褶皱和断层的影响下，滑床岩体的完整性遭到严重破坏，节理裂隙密集发育，形成了各种不同规模和方向的结构体。这些结构体之间的相互作用和错动，进一步影响了滑床的稳定性。滑床的稳定性对滑坡的发展具有重要制约作用。由于滑床的岩石强度和岩体结构存在差异，其抵抗滑坡体滑动的能力也不同。在滑坡体的推动下，滑床岩体可能发生变形、破裂，甚至产生新的滑动面。当滑床的稳定性降低到一定程度时，滑坡体将加速滑动，导致滑坡灾害的发生。滑床的岩石类型、岩体结构以及稳定性等特征，与滑体和滑带相互作用，共同影响着杜家营滑坡的稳定性和发展趋势。在研究滑坡的工程地质特征和制定治理方案时，必须充分考虑滑床的这些特征，采取相应的措施增强滑床的稳定性，如对滑床岩体进行加固、改善地下水的排泄条件等，以防止滑坡的进一步发展。3.3滑坡变形破坏特征3.3.1变形迹象通过对杜家营滑坡的长期监测和实地勘查，发现滑坡体上存在多种明显的变形迹象，这些迹象反映了滑坡的变形过程和发展趋势。地面开裂是滑坡体上最为显著的变形迹象之一。在滑坡后缘，出现了多条平行于滑坡后壁的拉张裂缝，裂缝宽度在0.1-0.5m之间，深度可达数米。这些拉张裂缝的产生是由于滑坡体后部在滑动过程中受到拉力作用，岩土体被拉伸断裂所致。随着滑坡的发展，拉张裂缝有逐渐加宽和加深的趋势。在2020年雨季，滑坡后缘的一条拉张裂缝宽度从原来的0.2m扩展到了0.35m，深度也增加了1-2m。在滑坡体中部两侧，还分布着剪切裂缝，这些裂缝与滑坡的滑动方向呈一定夹角，是由于滑坡体在滑动过程中受到剪切力作用而产生的。剪切裂缝的宽度相对较小，一般在0.05-0.1m之间，但它们的存在进一步破坏了滑坡体的整体性，增加了滑坡的不稳定性。滑坡体上的建筑物也出现了明显的变形。位于滑坡前缘的一些房屋墙体出现了倾斜和开裂现象，裂缝宽度在0.05-0.2m之间，部分房屋的门窗变形严重，无法正常开关。这些建筑物的变形是由于滑坡体的滑动导致地基不均匀沉降，从而使建筑物受到不均匀的外力作用而发生损坏。根据对建筑物变形的监测数据，在2018-2020年期间，滑坡前缘某房屋墙体的倾斜角度从2°增加到了5°，裂缝宽度也不断增大。树木歪斜也是滑坡变形的一个重要表现。在滑坡体上，可以观察到许多树木呈现出倾斜生长的状态，部分树木甚至出现了弯曲变形，形成了“马刀树”的现象。这是因为滑坡体的缓慢滑动使得树木的根系受到不均匀的拉力作用，导致树木生长方向发生改变。通过对树木歪斜程度的测量和分析，可以大致推断滑坡体的滑动方向和滑动速率。在滑坡体中前部，一些树木的倾斜角度达到了15°-30°，表明该区域的滑坡变形较为明显。这些变形迹象并非孤立出现，而是相互关联、相互影响的。地面开裂为雨水的渗入提供了通道，进一步软化岩土体，加剧滑坡的变形；建筑物的变形反映了滑坡对人类工程设施的破坏，同时也会影响滑坡体的应力分布；树木歪斜则直观地展示了滑坡体的缓慢滑动过程。通过对这些变形迹象的综合分析，可以发现它们具有一定的发展规律。在雨季，由于降雨的作用，变形迹象会明显加剧，地面开裂加宽加深，建筑物变形加速，树木歪斜程度增大；而在旱季，变形迹象的发展相对缓慢，但仍会持续进行。这些变形迹象的发展过程和规律，为我们判断滑坡的稳定性和预测滑坡的发展趋势提供了重要依据。3.3.2破坏模式根据对杜家营滑坡的变形破坏过程的详细分析，确定其破坏模式为降雨触发型牵引式滑坡。这种破坏模式的形成机制与该地区的地质条件、气象条件以及岩土体特性密切相关。杜家营滑坡所在区域的地层岩性主要为侏罗系中统上沙溪庙组砂岩、泥岩互层，这种软硬相间的地层结构使得岩体的力学性质差异较大。泥岩作为相对软弱的夹层，遇水易软化、崩解，抗风化能力弱。在长期的风化、降雨等作用下，泥岩逐渐被侵蚀、软化，形成了相对软弱的结构面，成为滑坡滑动的潜在控制面。该地区属于亚热带湿润季风气候区，年平均降水量丰富，降水主要集中在5-9月，且降雨强度较大。大量的雨水迅速渗入地下，使得地下水位急剧上升。地下水位的上升增加了岩土体的重量，使下滑力增大；同时，地下水对滑带土的软化作用显著，降低了滑带土的抗剪强度，减小了抗滑力。当降雨强度和持续时间达到一定程度时，滑坡体下部的岩土体首先达到极限平衡状态，开始发生滑动。滑坡体下部的滑动是整个破坏过程的起始点。由于下部岩土体的滑动，使得上部岩土体失去了支撑，在重力作用下，上部岩土体也逐渐开始变形滑动。这种由下部向上部发展的滑动过程，呈现出牵引式的特点，使得滑坡范围逐渐扩大。在滑动过程中，滑坡体沿着预先形成的软弱结构面（滑带）滑动，滑带土在强烈的剪切作用下，结构被破坏，颗粒定向排列，抗剪强度进一步降低，加速了滑坡的滑动。降雨触发型牵引式滑坡的特点主要包括以下几个方面：其发生与降雨密切相关，降雨是主要的触发因素，强降雨事件往往会导致滑坡的发生或加剧滑坡的发展；滑动过程具有明显的牵引特征，从下部开始滑动，逐渐向上部扩展，形成上小下大的塔式外貌；滑动速度相对较慢，一般多具横向张性裂隙发育，表面多呈阶梯状或陡坎状，这是由于各段滑体运动速度的差异所导致的。这种破坏模式对滑坡的稳定性和发展趋势产生了重要影响。由于滑坡是由下部向上部逐渐牵引滑动，一旦下部岩土体开始滑动，滑坡就很难自行停止，容易导致滑坡规模的不断扩大。降雨的不确定性使得滑坡的发生时间和发展程度难以准确预测，增加了防治工作的难度。在制定治理方案时，需要充分考虑降雨触发型牵引式滑坡的特点，采取有效的排水措施，降低地下水位，增强滑带土的抗剪强度，以阻止滑坡的进一步发展。四、杜家营滑坡稳定性分析4.1分析方法选择为了全面、准确地评估杜家营滑坡的稳定性，本研究选用工程地质法和极限平衡法相结合的方式进行分析。工程地质法主要基于地质学原理，通过对滑坡区域的地形地貌、地层岩性、地质构造、气象水文条件以及滑坡的变形破坏迹象等方面进行综合分析，定性地判断滑坡的稳定性状态。这种方法能够充分考虑地质环境因素对滑坡稳定性的影响，从宏观角度把握滑坡的整体特征和发展趋势。例如，通过观察滑坡后壁的形态、滑体的裂缝分布以及滑带的物质组成和结构特征等，可以初步判断滑坡的稳定性程度和可能的滑动模式。极限平衡法是依据静力平衡原理，分析边坡在各种破坏模式下的受力状态，通过计算滑坡体上的抗滑力和下滑力之间的关系，来定量评价边坡的稳定性，其计算结果以稳定性系数的形式呈现。在实际工程中，常用的极限平衡法包括费伦纽斯（Fellenius）法、毕肖普（Bishop）法、简布（Janbu）法等。本研究采用传递系数法，该方法是一种较为常用的极限平衡分析方法，它假设滑坡体沿某一滑动面滑动，将滑坡体划分为若干个垂直条块，通过依次传递各条块的剩余下滑力，计算整个滑坡体的稳定性系数。这种方法考虑了条块间的作用力，计算过程相对简单，能够较好地适用于杜家营滑坡这种具有折线形滑动面的情况。工程地质法能够提供滑坡的地质背景和定性信息，而极限平衡法可以给出具体的稳定性系数，实现定量分析。两者相互补充，能更全面、准确地评估杜家营滑坡的稳定性，为后续治理方案的制定提供可靠依据。4.2工程地质法分析从地形地貌来看，杜家营滑坡处于低山丘陵地带，地势起伏较大，相对高差可达300米左右。滑坡区域地形坡度变化明显，中后部地形坡度约为20°，呈阶梯状分布，阶坎高度在1-4米之间；中前部地形坡度较陡，在20°-30°之间。这种上缓下陡的坡面形态，使得滑坡体上部岩土体在重力作用下，更容易产生向下的滑动趋势。坡面总体呈凸坡面，雨水在坡面上的径流速度较快，对坡面岩土体的冲刷作用较强，容易破坏岩土体结构，降低其抗滑能力。滑坡平面形态呈“圈椅”状，这种形态表明滑坡体在滑动过程中，后部受到较大拉力，两侧受到侧向挤压和剪切，前缘向前突出，反映了滑坡体的受力状态和运动方式，也暗示了滑坡的不稳定性。地层岩性对滑坡稳定性有着重要影响。杜家营滑坡区域出露的地层岩性主要为侏罗系中统上沙溪庙组砂岩、泥岩互层，以及第四系全新统滑坡堆积层、残坡积层。侏罗系中统上沙溪庙组砂岩、泥岩互层中，砂岩抗压强度较高，但完整性较差，节理裂隙发育；泥岩抗压强度低，遇水易软化、崩解，抗风化能力弱。这种软硬相间的地层结构，使得岩体力学性质差异较大，在风化、降雨等因素作用下，泥岩易被侵蚀、软化，形成相对软弱的结构面，为滑坡的形成提供了潜在的滑动面。第四系全新统滑坡堆积层和残坡积层结构松散，工程性质较差。滑坡堆积层物质组成主要为砂岩、泥岩碎块石及粉质黏土，碎块石含量高，粒径大小不一，分选性差，结构松散；粉质黏土含水率高，强度较低。残坡积层主要由粉质黏土夹少量碎石组成，结构较松散，强度也较低。这些松散的堆积层在自重和外部荷载作用下，容易产生变形和滑动，增加了滑坡的不稳定性。地质构造对滑坡稳定性的控制作用也不容忽视。杜家营滑坡区域位于扬子准地台西缘龙门山褶断带与米仓山台穹的过渡地带，地质构造复杂，经历了多期构造运动。区内褶皱构造发育，主要为轴向近南北向的紧闭褶皱，褶皱核部和转折端部位节理裂隙极为发育，岩石破碎，完整性遭到严重破坏，力学强度大幅降低。断层构造常见，主要有北东向和北西向两组断层，这些断层相互切割、错动，进一步破坏了岩体的完整性，使得岩体结构更加复杂。构造运动破坏了岩体的完整性，形成大量节理裂隙，为地下水的运移和储存提供了通道，地下水的软化和侵蚀作用降低了岩土体的抗剪强度；褶皱和断层造成的地层岩性变化和岩体结构差异，使得斜坡岩土体力学性质不均一，容易产生应力集中现象；构造运动产生的地应力场改变了斜坡的应力状态，增加了斜坡的不稳定性。气象与水文条件也是影响滑坡稳定性的重要因素。杜家营滑坡所在地区属于亚热带湿润季风气候区，年平均降水量丰富，降水主要集中在5-9月，且降雨强度较大。强降雨等极端气象条件与滑坡发生具有显著的相关性，大量雨水渗入地下，使地下水位迅速上升，增加了岩土体的重量，同时软化岩土体，降低其抗剪强度，极大地增加了滑坡发生的可能性。地表水主要来源于大气降水形成的坡面径流以及五郎河，五郎河对滑坡前缘坡脚的冲刷作用，以及坡面径流对冲沟两侧坡体的冲刷作用，都破坏了坡体的稳定性。地下水主要赋存于第四系松散堆积层和侏罗系中统上沙溪庙组砂岩、泥岩互层中，地下水水位的上升会增加岩土体重量，降低抗滑力，地下水径流不畅会积聚孔隙水压力，削弱岩土体强度，从而降低滑坡的稳定性。通过工程地质法分析可知，杜家营滑坡在地形地貌、地层岩性、地质构造以及气象与水文条件等多种因素的综合作用下，处于不稳定状态，具有较高的滑坡风险，且在未来有进一步发展的趋势，需要采取有效的治理措施来提高其稳定性。4.3极限平衡法分析4.3.1计算模型建立根据杜家营滑坡的实际几何形状、地层岩性分布以及勘探得到的滑带位置和形状等信息，建立极限平衡分析的计算模型。将滑坡体视为由多个垂直条块组成，各条块之间相互作用。为准确模拟滑坡的实际情况，模型的计算边界确定如下：左右边界选取滑坡体两侧相对稳定的岩土体边界，以确保边界条件的合理性；上边界为滑坡体的地表，直接与大气接触；下边界为滑带，作为滑坡体滑动的控制面。在确定计算条件时，充分考虑了多种工况。天然工况下，仅考虑滑坡体自身重力作用，不考虑外部荷载和地下水的动态变化，以此作为基础工况来评估滑坡在自然状态下的稳定性。暴雨工况下，增加了由于降雨导致的地下水水位上升的影响。地下水水位上升会使岩土体重度增加，同时降低滑带土的抗剪强度，通过合理设定地下水水位上升的高度和相关参数，模拟暴雨对滑坡稳定性的影响。地震工况下，考虑了地震作用产生的惯性力。根据该地区的地震设防烈度和相关规范，确定地震加速度的取值，将地震惯性力以水平和竖向力的形式施加到滑坡体上，模拟地震对滑坡稳定性的影响。通过对不同工况的模拟分析，可以全面评估滑坡在各种可能情况下的稳定性，为治理方案的制定提供更全面、准确的依据。4.3.2参数选取为确保极限平衡法分析的准确性，通过多种手段获取岩土体的物理力学参数。在室内试验方面，对采集的滑体、滑带和滑床的岩土样本进行了系统的物理力学性质试验。对于滑体样本，进行了密度试验，采用环刀法测定其天然密度，得到滑体的天然密度在1.8-2.0g/cm³之间；通过烘干法测定含水率，确定其含水率范围；进行直剪试验，测定滑体的抗剪强度参数，得到内摩擦角在18°-22°之间，黏聚力在10-15kPa之间。对于滑带样本，同样进行了密度、含水率试验，滑带土的含水率较高，一般在25%-35%之间。采用三轴剪切试验测定滑带土的抗剪强度，得到内摩擦角在10°-15°之间，黏聚力在5-10kPa之间。对于滑床的砂岩和泥岩样本，进行了抗压强度试验，砂岩的抗压强度一般在30-50MPa之间，泥岩的抗压强度在5-15MPa之间。现场测试也是获取参数的重要手段。利用原位测试技术，如标准贯入试验、动力触探试验等，对滑体和滑带的力学性质进行现场测定。通过标准贯入试验，得到滑体和滑带土的贯入击数，进而推算其力学参数，与室内试验结果相互验证和补充，提高参数的准确性。参考类似地质条件下的工程案例，进行工程地质类比。查阅相关文献资料，收集与杜家营滑坡地质条件相似地区的滑坡工程案例，对比分析其岩土体参数，对本次研究获取的参数进行合理性验证和调整，确保参数符合实际地质情况。综合室内试验、现场测试和工程地质类比的结果，确定了用于极限平衡法分析的岩土体物理力学参数。这些参数准确反映了杜家营滑坡岩土体的特性，为后续的稳定性计算提供了可靠的数据支持。4.3.3稳定性计算运用传递系数法进行稳定性计算，该方法是极限平衡法中的一种常用方法，适用于具有折线形滑动面的滑坡稳定性分析。其基本原理是将滑坡体沿滑动面划分为若干个垂直条块，依次计算每个条块的下滑力和抗滑力，并考虑条块间的作用力，通过传递系数将前一个条块的剩余下滑力传递到下一个条块，最终计算出整个滑坡体的稳定性系数。计算公式如下：F_s=\frac{\sum_{i=1}^{n}\left(c_il_i+\left(W_i\cos\alpha_i-U_i-e_i\sin\alpha_i\right)\tan\varphi_i\right)\prod_{j=1}^{i-1}\psi_j}{\sum_{i=1}^{n}\left(W_i\sin\alpha_i+e_i\cos\alpha_i\right)}其中，F_s为滑坡稳定性系数；c_i为第i条块滑带土的黏聚力（kPa）；l_i为第i条块滑带长度（m）；W_i为第i条块的自重（kN）；\alpha_i为第i条块滑带的倾角（°）；U_i为第i条块滑带的孔隙水压力（kN）；e_i为第i条块的水平地震力（kN）；\varphi_i为第i条块滑带土的内摩擦角（°）；\psi_j为第j条块的剩余下滑力传递至第j+1块时的传递系数。根据建立的计算模型和选取的参数，分别计算杜家营滑坡在天然工况、暴雨工况和地震工况下的稳定性系数。在天然工况下，不考虑孔隙水压力和地震力，仅计算滑坡体自重作用下的稳定性系数；在暴雨工况下，考虑由于降雨导致的地下水位上升，增加孔隙水压力的影响，重新计算稳定性系数；在地震工况下，根据该地区的地震设防烈度，确定水平地震力的大小，加入地震力的作用，计算稳定性系数。计算结果表明，在天然工况下，滑坡的稳定性系数为1.05，处于基本稳定状态，但安全储备相对较低；在暴雨工况下，由于地下水位上升，岩土体重度增加，滑带土抗剪强度降低，稳定性系数降至0.98，滑坡处于欠稳定状态，有发生滑动的风险；在地震工况下，考虑地震力的作用后，稳定性系数进一步降低至0.92，滑坡的稳定性明显下降，处于不稳定状态，在地震作用下极有可能发生滑动。通过对不同工况下稳定性系数的计算和分析，可以清晰地了解滑坡在各种情况下的稳定性状态，为制定针对性的治理方案提供了关键的定量依据。五、杜家营滑坡治理方案研究5.1治理原则与目标杜家营滑坡治理遵循安全、经济、合理、可行的原则，力求全面消除滑坡隐患，保障居民生命财产安全。安全原则是治理工作的首要准则，通过全面、科学地分析滑坡的工程地质特征、稳定性以及变形破坏规律，精准识别潜在的安全风险，采取可靠的工程措施，确保滑坡治理工程在施工过程中和运营期的绝对安全，杜绝因治理工程不当引发的二次灾害，为居民创造一个安全稳定的居住环境。经济原则要求在保证治理效果的前提下，合理控制工程成本。充分考虑当地的经济发展水平和资金投入能力，对不同的治理方案进行详细的技术经济分析，对比各方案的工程投资、施工难度、后期维护成本等因素，优化设计，选择性价比高的方案，避免不必要的浪费，提高资金的使用效率。合理原则强调治理方案应与滑坡的实际情况紧密结合，充分考虑地形地貌、地层岩性、地质构造、气象水文等因素的影响。依据滑坡的变形破坏模式和稳定性分析结果，选择合适的治理措施，确保治理方案的科学性和合理性，使治理工程能够有效地针对滑坡的关键问题，达到预期的治理效果。可行原则注重治理方案在技术和施工方面的可操作性。所选用的技术和工艺应成熟可靠，施工设备和材料易于获取，施工过程安全可控。同时，考虑施工过程中可能遇到的各种问题，制定合理的施工组织设计和应急预案，确保治理工程能够顺利实施。治理目标明确为彻底消除滑坡隐患，使滑坡体达到长期稳定状态。通过采取有效的治理措施，提高滑坡的稳定性系数，使其满足相关规范和标准的要求。在天然工况、暴雨工况和地震工况等各种可能情况下，滑坡稳定性系数均应达到安全阈值以上，确保滑坡不会发生滑动，保障滑坡前缘43户共179名居民以及五权场镇2500余居民的生命财产安全。保障滑坡前缘场镇公路的正常通行也是重要目标之一。通过治理工程，防止滑坡对公路造成破坏，确保公路的结构安全和交通流畅，维持当地的正常交通运输秩序，促进区域经济的稳定发展。保护滑坡前缘的五郎河河道畅通同样关键。避免滑坡体冲入河道导致河道阻塞，引发洪水等次生灾害，维护河道的自然生态环境和水利功能，保障河流上下游居民的生命财产安全和生态环境的稳定。在治理过程中，注重环境保护，减少对周边生态环境的破坏。采取生态友好的治理措施，尽量降低工程施工对植被、土壤等生态要素的影响，保护生物多样性，实现滑坡治理与生态环境保护的协调发展。5.2治理方案设计5.2.1方案一：三级抗滑桩+截水沟在方案一中，三级抗滑桩的布置是关键环节。第一级抗滑桩布置在滑坡的中前部，此处是滑坡体下滑力较大的区域，通过设置抗滑桩可以有效阻挡滑坡体的下滑趋势。桩间距设计为4m，这样的间距既能保证抗滑桩之间的相互协同作用，又能合理控制工程成本。桩长根据滑坡体的厚度和滑带的深度确定，一般为15-18m，确保抗滑桩能够深入到滑带以下稳定的岩土体中，提供足够的锚固力。桩截面尺寸为2m×2.5m，这种较大的截面尺寸可以增强抗滑桩的抗弯和抗剪能力，更好地抵抗滑坡体的推力。第二级抗滑桩位于滑坡中部，该区域的下滑力相对第一级抗滑桩处稍小，但仍然对滑坡的稳定性有重要影响。桩间距设置为5m，桩长为12-15m，截面尺寸为1.5m×2m。这样的参数设计既能满足该区域的抗滑需求，又能在保证工程效果的前提下，适当降低工程成本。第三级抗滑桩布置在滑坡后缘，主要作用是阻止滑坡体后缘的进一步变形和滑动。桩间距为6m，桩长8-12m，截面尺寸为1m×1.5m。后缘的滑坡体厚度相对较薄，下滑力也较小，因此抗滑桩的参数相对较小，但同样能够有效地发挥其抗滑作用。抗滑桩采用钢筋混凝土结构，以确保其具有足够的强度和耐久性。混凝土强度等级为C30，这种强度等级的混凝土能够满足抗滑桩在长期使用过程中承受各种荷载的要求。钢筋配置根据抗滑桩的受力计算确定，主筋采用HRB400级钢筋，直径为25mm，沿桩身均匀布置，以提供足够的抗拉强度。箍筋采用HPB300级钢筋，直径为10mm，间距为200mm，用于增强桩身的抗剪能力和约束混凝土的横向变形。截水沟的设计对于减少地表水对滑坡体的渗入至关重要。在滑坡后缘，沿滑坡边界开挖截水沟，其走向与滑坡后壁大致平行，这样可以最大限度地拦截从滑坡后部来的地表水，防止其流入滑坡体。截水沟的断面尺寸为梯形，上底宽1.2m，下底宽0.8m，深度为1m。这种断面尺寸能够保证在降雨量大时，截水沟有足够的过水能力，及时排除地表水。截水沟采用M7.5浆砌片石砌筑，片石应质地坚硬、无风化，以保证截水沟的耐久性。沟壁厚度为0.3m，沟底厚度为0.35m，确保截水沟能够承受水流的冲刷和压力。在截水沟的内侧和底部，采用1:2水泥砂浆抹面，厚度为2cm，以防止水流渗漏，提高截水沟的防渗性能。在截水沟的适当位置设置沉砂池和出水口，沉砂池用于沉淀水流中的泥沙，防止泥沙淤积堵塞截水沟，出水口则将截水沟中的水引至安全地点排放，避免对滑坡体和周边环境造成影响。5.2.2方案二：两级抗滑桩+截水沟+挡土墙在方案二中，两级抗滑桩和挡土墙相互配合，共同提高滑坡的稳定性。第一级抗滑桩布置在滑坡前缘，这里是滑坡体下滑的前沿阵地，承受着巨大的下滑力。桩间距设计为3.5m，桩长根据滑坡前缘的地质条件和滑带深度确定，一般为18-20m，以确保抗滑桩能够深入稳定地层，提供强大的锚固力。桩截面尺寸为2.5m×3m，这种较大的尺寸能够增强抗滑桩的承载能力和抗滑性能，有效阻挡滑坡体的下滑。抗滑桩采用钢筋混凝土结构，混凝土强度等级为C35，主筋采用HRB400级钢筋，直径为32mm，沿桩身均匀布置，以承受巨大的拉力；箍筋采用HPB300级钢筋，直径为12mm，间距为150mm，增强桩身的抗剪能力。第二级抗滑桩位于滑坡中部，桩间距为4.5m，桩长15-18m，截面尺寸为2m×2.5m。混凝土强度等级为C30，主筋采用HRB400级钢筋，直径为28mm，箍筋采用HPB300级钢筋，直径为10mm，间距为200mm。中部抗滑桩主要分担滑坡体中部的下滑力，与前缘抗滑桩协同作用，共同稳定滑坡体。截水沟的设计与方案一类似，但在一些细节上有所优化。同样在滑坡后缘沿边界开挖截水沟，走向与滑坡后壁平行。断面尺寸为梯形，上底宽1.5m，下底宽1m，深度为1.2m，过水能力更强。采用M10浆砌片石砌筑，沟壁厚度为0.35m，沟底厚度为0.4m，比方案一的截水沟更加坚固耐用。内侧和底部采用1:2.5水泥砂浆抹面，厚度为2.5cm，防渗性能更好。沉砂池和出水口的设置更加合理，沉砂池的容积适当增大，以更好地沉淀泥沙，出水口的排水能力也相应提高，确保截水沟排水畅通。挡土墙设置在滑坡前缘抗滑桩之间，起到辅助抗滑和防止坡脚被冲刷的作用。挡土墙采用重力式挡土墙，墙高根据滑坡前缘地形和抗滑需求确定，一般为4-6m。墙顶宽度为1m，墙底宽度根据墙体高度和地基承载力计算确定，以保证挡土墙的稳定性。墙体采用C25混凝土浇筑，在墙体内设置泄水孔，间距为2m，呈梅花形布置，以排除墙后积水，减小墙后土压力。在挡土墙的基础底部铺设一层厚0.3m的碎石垫层，以提高地基的承载能力和稳定性。两级抗滑桩、截水沟和挡土墙在该方案中协同工作。抗滑桩主要承受滑坡体的下滑力，将其传递到稳定地层；截水沟拦截地表水，减少其对滑坡体的渗入，降低滑坡体的重量和孔隙水压力，从而提高滑坡的稳定性；挡土墙则辅助抗滑桩，增强坡脚的稳定性，防止坡脚被地表水冲刷和侵蚀，进一步保障滑坡的整体稳定性。5.3方案比选从工程投资角度来看，方案一布置了三级抗滑桩，桩的数量较多，且桩长和截面尺寸在不同位置虽有差异，但总体工程量较大，导致抗滑桩的混凝土用量、钢筋用量以及施工费用较高。截水沟的建设费用相对固定，但综合抗滑桩和截水沟的费用，方案一的总工程投资相对较高。而方案二布置两级抗滑桩，桩的数量相对较少，尤其是在滑坡后缘未设置抗滑桩，减少了一部分工程费用。挡土墙的建设费用相对抗滑桩较低，虽然截水沟在尺寸和砌筑材料上有所优化，费用略有增加，但总体而言，方案二的工程投资比方案一低，具有一定的经济优势。在施工难度方面，方案一的三级抗滑桩布置，使得施工过程中需要在不同位置、不同深度进行桩的施工，施工工序相对复杂。抗滑桩的长度和截面尺寸变化较多，对施工技术和设备的要求较高，施工过程中需要频繁调整施工参数，增加了施工的难度和不确定性。方案二的两级抗滑桩布置相对简单，桩的参数变化较少，施工过程相对容易控制。挡土墙的施工工艺相对成熟，技术难度较低。截水沟的优化设计虽然对施工精度有一定要求，但整体施工难度相比方案一有所降低。从治理效果分析，方案一通过三级抗滑桩的布置，能够在滑坡的不同部位提供抗滑力，对滑坡体的稳定性有较好的控制作用。截水沟能够有效拦截地表水，减少其对滑坡体的渗入，降低滑坡体的重量和孔隙水压力，从而提高滑坡的稳定性。然而，由于抗滑桩布置相对分散，在抵抗较大规模滑坡推力时，可能存在局部抗滑能力不足的问题。方案二的两级抗滑桩和挡土墙协同工作，能够在滑坡前缘和中部形成较强的抗滑体系，有效阻挡滑坡体的下滑。截水沟的优化设计提高了排水能力，能更好地减少地表水对滑坡体的影响。挡土墙还能增强坡脚的稳定性，防止坡脚被地表水冲刷和侵蚀，相比方案一，方案二的治理效果更加全面和可靠。后期维护方面，方案一的三级抗滑桩数量多，维护检查的工作量较大，且不同位置的抗滑桩参数不同，维护难度较高。截水沟的维护相对简单，但总体而言，方案一的后期维护成本较高。方案二的两级抗滑桩和挡土墙数量相对较少，维护检查工作相对容易开展，后期维护成本较低。综合比较，方案二在工程投资上具有优势，施工难度相对较低，治理效果更加全面可靠，后期维护成本也较低。因此，方案二更适合作为杜家营滑坡的治理方案，能够在满足治理要求的前提下，实现经济效益和社会效益的最大化。六、治理方案设计计算与数值模拟6.1治理方案二设计计算6.1.1抗滑桩设计计算抗滑桩作为方案二中重要的抗滑结构，其设计计算需综合考虑滑坡推力、地质条件等多方面因素。首先，通过传递系数法计算滑坡推力，具体计算公式为：P_i=P_{i-1}\psi_{i-1}+KW_i\sin\alpha_i-W_i\cos\alpha_i\tan\varphi_i-c_il_i其中，P_i为第i条块的剩余下滑力（kN）；P_{i-1}为第i-1条块传递至第i条块的剩余下滑力（kN）；\psi_{i-1}为第i-1条块剩余下滑力传递至第i条块时的传递系数；K为滑坡推力安全系数，取1.3；W_i为第i条块的自重（kN）；\alpha_i为第i条块滑带的倾角（°）；\varphi_i为第i条块滑带土的内摩擦角（°）；c_i为第i条块滑带土的黏聚力（kPa）；l_i为第i条块滑带长度（m）。根据地质勘查资料，将滑坡体划分为多个条块，依次计算各条块的剩余下滑力，得到滑坡在不同位置的推力分布。经计算，在滑坡前缘，最大滑坡推力约为800kN/m；在滑坡中部，最大滑坡推力约为500kN/m。依据滑坡推力，拟定抗滑桩的桩长。对于前缘抗滑桩，考虑到此处滑坡推力较大，滑带较深，桩长确定为18-20m，其中锚固段长度不小于6m，以确保抗滑桩能够深入稳定地层，提供足够的锚固力。中部抗滑桩桩长为15-18m，锚固段长度不小于5m。桩径的确定需考虑桩身的承载能力和抗弯、抗剪性能。前缘抗滑桩采用矩形截面，尺寸为2.5m×3m，这种较大的截面尺寸能够有效增强桩身的承载能力，抵抗巨大的滑坡推力。中部抗滑桩截面尺寸为2m×2.5m，既能满足抗滑需求，又能合理控制工程成本。抗滑桩的配筋计算至关重要。根据《混凝土结构设计规范》（GB50010-2010），采用极限状态设计方法。对于前缘抗滑桩，主筋采用HRB400级钢筋，直径为32mm，沿桩身均匀布置，每侧布置8根，以承受巨大的拉力。箍筋采用HPB300级钢筋，直径为12mm，间距为150mm，增强桩身的抗剪能力。中部抗滑桩主筋采用HRB400级钢筋，直径为28mm，每侧布置6根；箍筋采用HPB300级钢筋，直径为10mm，间距为200mm。通过上述设计计算，确保抗滑桩在满足承载能力和稳定性要求的同时，能够有效抵抗滑坡推力，保障滑坡治理工程的安全。6.1.2挡土墙设计计算挡土墙主要承受土压力和滑坡推力，其设计计算需准确分析这些荷载。土压力计算采用朗肯土压力理论，对于墙后填土为粉质黏土的情况，主动土压力系数K_a=\tan^2(45°-\frac{\varphi}{2})，其中\varphi为粉质黏土的内摩擦角，取20°，则K_a=\tan^2(45°-\frac{20°}{2})\approx0.49。主动土压力强度计算公式为：p_a=\gammazK_a-2c\sqrt{K_a}，其中\gamma为填土重度，取18kN/m³；z为计算点深度；c为填土黏聚力，取10kPa。当z=0时，p_{a0}=-2c\sqrt{K_a}=-2×10×\sqrt{0.49}\approx-14kPa，出现负值，说明墙顶处存在拉应力，但实际工程中，墙顶处填土可能会产生裂缝，拉应力可忽略不计。当z=6m（挡土墙高度）时，p_{a6}=18×6×0.49-2×10×\sqrt{0.49}=52.92-14=38.92kPa。单位长度挡土墙上的主动土压力合力E_a=\frac{1}{2}×38.92×6=116.76kN/m，作用点位于墙底以上\frac{6}{3}=2m处。考虑滑坡推力作用时，假设滑坡推力呈三角形分布，在墙顶处为0，墙底处达到最大值，根据滑坡稳定性分析结果，墙底处滑坡推力约为200kN/m。墙体尺寸设计方面，采用重力式挡土墙，墙高根据滑坡前缘地形和抗滑需求确定为6m。墙顶宽度取1m，以保证施工和使用的安全；墙底宽度根据墙体高度、土压力和地基承载力计算确定，根据公式b\geq\frac{KE_ah}{f}，其中K为抗滑稳定安全系数，取1.3；E_a为主动土压力合力；h为合力作用点距墙底的距离；f为地基承载力特征值，取150kPa。代入数据计算可得b\geq\frac{1.3×116.76×2}{150}\approx2.03m，考虑一定的安全储备，墙底宽度取2.5m。基础埋深根据地质条件和抗倾覆稳定性要求确定，由于滑坡前缘地下水位较高，基础埋深取1.5m，以保证基础位于地下水位以上，且处于稳定的地基土层中。稳定性系数计算包括抗滑稳定性和抗倾覆稳定性。抗滑稳定性系数K_s=\frac{(G+E_{az})\tan\varphi+cA}{E_{ax}}，其中G为挡土墙自重；E_{az}为主动土压力垂直分力；E_{ax}为主动土压力水平分力；\varphi为地基土的内摩擦角，取25°；c为地基土的黏聚力，取15kPa；A为基础底面面积。经计算，K_s\approx1.5>1.3，满足抗滑稳定性要求。抗倾覆稳定性系数K_t=\frac{Gx_0+E_{az}x_f}{E_{ax}z_f}，其中x_0为挡土墙重心距墙趾的水平距离；x_f为主动土压力垂直分力作用点距墙趾的水平距离；z_f为主动土压力水平分力作用点距墙趾的垂直距离。经计算，K_t\approx2.0>1.6，满足抗倾覆稳定性要求。通过上述设计计算，确保挡土墙能够稳定地承受土压力和滑坡推力，发挥其抗滑和保护坡脚的作用。6.1.3截水沟设计计算截水沟的设计流量是关键参数，采用《滑坡防治工程设计与施工技术规范》（DZ/T0219-2006）中的小汇水面积设计流量公式计算，公式为：Q=\varphiSF其中，Q为设计频率地表水汇量（m³/s）；\varphi为当地径流系数，根据杜家营滑坡区域的地形、植被等条件，取0.6；S为设计降雨强度，根据该地区的气象资料，十年一遇1h降雨强度取70mm/h，换算为\frac{70}{3600}m/s；F为汇水面积，通过地形测绘和分析，确定滑坡后缘截水沟的汇水面积为0.3km²，即300000m²。将参数代入公式计算可得：Q=0.6×\frac{70}{3600}×300000\approx35m³/s截水沟采用梯形断面，根据流量计算过水断面相关参数。湿周长X=b+2\sqrt{h^2+(mh)^2}，其中b为沟底宽；h为水深；m为边坡系数，取1.5。水力半径R=\frac{bh+mh^2}{X}。沟床糙率n根据采用的M10浆砌片石材料，取0.017。纵向坡降i根据地形条件，取0.005。谢才系数C=\frac{R^{\frac{1}{6}}}{n}，流速v=C\sqrt{Ri}，流量Q=(bh+mh^2)v。通过试算，假设沟底宽b=1m，当水深h=1.2m时，湿周长X=1+2\sqrt{1.2^2+(1.5×1.2)^2}\approx6.2m，水力半径R=\frac{1×1.2+1.5×1.2^2}{6.2}\approx0.46m，谢才系数C=\frac{0.46^{\frac{1}{6}}}{0.017}\approx52.3，流速v=52.3\sqrt{0.46×0.005}\approx2.5m/s，流量Q=(1×1.2+1.5×1.2^2)×2.5=7.2m³/s，不满足设计流量要求。调整沟底宽b=1.5m，水深h=1.5m时，湿周长X=1.5+2\sqrt{1.5^2+(1.5×1.5)^2}\approx8.1m，水力半径R=\frac{1.5×1.5+1.5×1.5^2}{8.1}\approx0.63m，谢才系数C=\frac{0.63^{\frac{1}{6}}}{0.017}\approx57.4，流速v=57.4\sqrt{0.63×0.005}\approx3.2m/s，流量Q=(1.5×1.5+1.5×1.5^2)×3.2=21.6m³/s，仍不满足要求。继续调整沟底宽b=2m，水深h=1.8m时，湿周长X=2+2