降雨对斜坡稳定性的影响

1、 本科毕业论文（设计） 题目： 降雨对斜坡稳定性的影响 以南平市斜坡为例 姓 名： 刘莹 学号： 2114481 院 （系）： 环境学院 专业： 地下水科学与工程 指导老师： 湖沉 职称： 讲师 评 阅 人： 赵锐 职称： 2014 年 06 月本科生毕业论文（设计）原创性声明本人以信誉声明：所呈交的毕业论文（设计）是在导师指导下进行的研究工作及取得的研究成果，论文中引用他人的文献、数据、图件、资料均已明确标注出，论文中的结论和结果为本人独立完成，不包含他人成果及为获得中国地质大学或其他教育机构的学位或证书而使用过的材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示

2、了谢意。 毕业论文作者（签字）： 签字日期： 年 月 日摘 要 斜坡失稳破坏是多山地区常见的一种自然灾害，而降雨是斜坡失稳破坏最重要的触发因素。降雨通过以下方式影响斜坡的稳定性：降雨时，雨水入渗使原本的非饱和土变成饱和土，孔隙水压力增大也即基质吸力减小，使斜坡的抗剪强度降低；其次，雨水入渗会增加土壤的重度，加大下滑力；同时，长时间的雨水入渗会抬高地下水位，当地下水位达到或超过潜在滑动面时，地下水会软化滑面从而降低斜坡的抗剪强度。在此文中，利用Geostudio2007软件主要分析降雨条件下，孔隙水压力的变化（也即基质吸力的变化）对斜坡稳定性的影响。本文以福建省南平市西芹镇兴华街道后山比较有代表

3、性的区房后斜坡为例，首先对其进行斜坡变形的宏观分析及敏感性因素的分析并计算得到天然条件下的斜坡稳定性系数。后取区房后斜坡28-28号地质剖面，使用Geostudio2007软件建立了边坡饱和非饱和降雨入渗有限元计算模型，得到不同降雨条件下雨水入渗时孔隙水压力的变化情况，并计算得到不同降雨强度和不同降雨历时下的斜坡稳定性系数，进而分析降雨对斜坡稳定性的影响。经过分析得到如下结论：（1）区房后斜坡受降雨或者其他不利因素影响时，中上部岩土体仍存在崩塌或滑坡的可能。计算结果表明，在天然情况下，此斜坡能基本保持稳定，稳定性系数为1.294。（2）在相同降雨时长的条件下，降雨强度越大，斜坡的稳定性系数减小

4、得越快。但在降雨强度大于雨水渗透系数时，各种不同降雨强度对斜坡稳定性的影响趋于相近。（3）在相同降雨强度的条件下，降雨历时越长，斜坡的稳定性系数减小得越快。（4）通过模拟斜坡监测点中孔隙水压力的变化可知，随着降雨的持续，孔隙水压力会一直升高，而斜坡的稳定性系数下降。因此降雨引起的孔隙水压力增大是导致斜坡失稳的重要原因。关键词：降雨入渗 斜坡稳定性 Geostudio2007 孔隙水压力Abstract Slope instability and failure is a common natural disaster in mountainous areas, and rainfall is

5、the most important trigger factor. Rainfall affects the slope stability in the following ways. Firstly, rainwater infiltration will make the original unsaturated soil become saturated soil so that pore water pressure increases ( ie, matric suction decreases) which will then reduce the shear strength

6、 of the slope. Secondly, rainwater infiltration will increase the soils density and the pore water pressure and hence increase decline in force. Whats more, prolonged rainfall infiltration will raise the water table and when the water table reaches or exceeds the potential sliding surface, the groun

7、dwater will soften the slip surface which will reduce the shear strength. In this paper, Geostudio is used to analyze how the changes of pore water pressure affects the stability of slope. This paper takes the example of a representative district slope in Fujian Province, Nanping, Xiqin town, Xinghu

8、a Street. The macro-deformation of the slope and sensitivity factors are analyzed and the slope stability coefficient under the natural condition is calculated. By using Geostudio, a saturated-unsaturated finite element model of rainfall infiltration is built in taking the 28-28 geological section o

9、f the district slope. I get the results of changes in pore water pressure and slope stability coefficient under different rainfall intensity and duration. Then I analyze the rainfalls impact on slope stability. Following are the conclusions. 1. district slope is possible to collapse or landslide in

10、the upper and middle part of the slope after another rainfall or other unfavorable factors. The calculation results show that under natural condition, this slope can basically stay stable and the slope stability coefficient is 1.294. 2. Under the same rainfall duration, the higher the rainfall inten

11、sity, the faster that slope stability coefficient decreases. However, when the rainfall intensity is greater than rainwater permeability coefficient, the impact of various different rainfall intensity is the barely the same on the stability of the slope no matter what the rainfall intensity is. 3. u

12、nder the same conditions of rainfall intensity, the longer that rainfall lasts, the faster stability of the slope coefficient decreases. 4. Through monitoring changes in simulated pore water pressure, we can know that with rainfall continues, pore water pressure increases and it leads to decrease of

13、 the slope stability coefficient. Therefore increase in pore water pressure caused by rainfall is an important factor that triggers the slope instability .Keywords: rainfall infiltration; slope stability; Geostudio2007; pore water pressure目 录第一章 前言- 1 -第一节 选题意义和依据- 1 -第二节 国内外研究现状- 1 -第三节 存在的问题和展望- 2

14、 -第四节 研究内容及技术路线- 3 -第二章 自然地理及地质概况- 5 -第一节 自然地理条件- 5 -第二节 斜坡基本特征- 6 -第三节 斜坡岩土体特征- 7 -第三章 斜坡稳定性分析- 9 -第一节 斜坡变形宏观分析- 9 -第二节 斜坡稳定性敏感因素分析- 9 -第三节 天然条件下的斜坡稳定性- 10 -第四章 降雨条件下的斜坡稳定性对比分析- 13 -第一节 基本原理- 13 -第二节 计算模型- 15 -第三节 数值模拟结果分析- 26 -第五章 结论和建议- 30 -第一节 结论- 30 -第二节 建议- 30 -致 谢- 31 -参考文献- 32 -中国地质大学(武汉)学士学

15、位论文第一章 前言第一节 选题意义和依据 不稳定的斜坡在内、外力的作用下，易发生斜坡的变形破坏。斜坡的变形破坏给工程建筑等带来的危害非常广泛，甚至造成生命财产的巨大损失，全世界各种典型实例不胜枚举。影响斜坡稳定性的因素复杂多样。其中内因主要包括斜坡岩土类型、岩土体结构、地质构造、水文地质条件。此外还有外因，例如风化作用、地表水和大气降水的作用、地震、人类工程活动等。其中降雨是影响斜坡稳定性和导致斜坡失稳的最主要和最普遍的环境因素。每年降雨导致的斜坡失稳破坏引起了数以千计的死亡和建筑物的严重破坏，因此在滑坡易发区建立模型并分析降雨对斜坡稳定性的影响是非常重要的。在许多案例中，陡峭土坡上的滑坡都是

16、由于降暴雨时雨水入渗到地下所引起的。增加的含水量导致了毛细吸力的下降以及土壤重量的上升，这两个过程均导致了斜坡的失稳。以下为降雨引起的斜坡失稳实例。2009年7月23日，云南凤庆发生山体滑坡，根据专家组的考察测定，此次滑坡是由于在不良地质条件下的连续降雨和库水位连续上升造成的，降雨降低了边坡稳定性，引起滑坡，造成灾害。2010年9月1日，云南保山市隆阳区发生山体滑坡，造成20户人家被掩埋。据了解，此次灾害的主要原因，首先是滑坡体本身的地形地质条件不良；其次是此地区之前大旱，后又出现降雨所致。2012年5月21日，一场强降雨造成重庆市多处发生山体滑坡事故，当地部分地区的交通因此而中断。2013年

17、5月贺州市出现长时间的强降水过程，导致某段公路出现多处滑坡。2014年4月12日，东乡县董岭乡出现山体滑坡现象，致使几名施工人员被埋。经调查分析，该区地形本就险峻，近期连续的降雨触发了此次滑坡。在我国，由于特殊的自然地理和地质条件所制约，斜坡地质灾害分布广泛，活动强烈，危害严重，是山区主要的工程动力地质作用。而降雨是引起滑坡的主要因素之一，因此有必要对降雨和斜坡稳定性的关系进行研究。 2010年6月12日至6月25日，福建省南平市出现长时间的持续性暴雨过程，过程降雨量达到252.6mm，西芹镇兴华街道后山山坡多处发生山体滑坡、崩塌等地质灾害，滑坡体导致兴华二路部分路段被埋，个别民房毁坏，严重威

18、胁当地居民生命财产安全，影响交通要道的畅通。本文以南平市西芹镇兴华街道区房后斜坡为例，用Geostudio2007软件的SLOPE/W模块和SEEP/W模块耦合计算各种不同降雨条件下斜坡的稳定性系数，通过对比稳定性系数的大小以及孔隙水压力的变化情况来分析降雨对斜坡稳定性的影响。第二节 国内外研究现状广义滑坡，即各类滑坡、崩塌、泥石流，可因多种外界因素的刺激而发生，如地震、火山、河流冲刷、融雪、降雨及不良人类活动等，其中以降雨诱发的滑坡（通常称为降雨型滑坡）在世界上分布最广，发生频率最高，给人类造成的危害最大。由于斜坡的变形破坏常威胁到人们的生命安全并会对经济活动带来巨大损失，而降雨又是引起斜坡

19、失稳破坏最重要的外界因素，因此关于降雨对斜坡稳定性影响的研究正日益受到重视，各国都有很多公开出版的相关研究文献。1999年吴宏伟和陈守义等分析了七个不同的因素对由于雨水入渗引起的斜坡暂态渗流场的影响，他们用有限元法对一个非饱和斜坡进行数值模拟，得到的结论是降雨强度、历时以及雨型对斜坡稳定性有较为明显的影响；土体的渗透系数影响显著；斜坡中相对隔水层的存在等因素对其也有影响。2003年张华、陈善雄通过对一维及二维非饱和土的降雨入渗分析及数值模拟，研究结果是Richards模型能很好地模拟二维降雨入渗过程以及地下水在土体中的运动；降雨入渗首先会使土体中压力升高，然后湿润锋面并向下推进，到达地下水面时

20、会使地下水位抬升；水在土体中运动分为入渗、水分重分布和排水三个过程。2003年兰恒星、周成虎等人以香港某地区降雨引起的浅层滑坡为例，在GIS的支持下，分析了瞬时降雨响应模型造成的瞬时孔隙水压力响应规律及滑坡危险性的时空分布特征，研究结果为短期集中降雨所引起的瞬时孔隙水压力变化是导致斜坡失稳的主要影响因素；降雨对浅层的低渗透性斜坡的稳定性影响很大，但对浅层的渗透性较好的斜坡稳定性影响较小。2007年在降雨对斜坡失稳的影响作用探讨一文中，高晓斐、陈植华等人讨论了降雨作用引发斜坡失稳的机理，并以施家梁子滑坡东变形体为实例进行分析，得出水的作用会增加斜坡体的下滑力同时会降低斜坡体的基质吸力、有效粘聚力

21、及有效内摩擦角, 从而降低斜坡体的抗剪强度。2010年程彬等采用了Geostudio软件中的边坡稳定性分析模块Slope/W、渗流分析模块Seep/W和应力应变分析模块Sigma/W对陕西一边坡为进行了应力场与渗流场的耦合分析及安全系数的计算。结论为斜坡稳定性在进行了应力场与渗流场的耦合时更符合实际。2013年程思、王运生在有限元法分析降雨条件下斜坡稳定性一文中采用工程地质条件分析和有限元分析法，进行了华蓥市马鞍坪滑坡的滑坡体在降雨条件下渗流与稳定性分析，结果表明降雨与斜坡稳定性的关系十分密切，降雨量是此区斜坡失稳下滑的主要触发因素，并且不同的降雨量对斜坡稳定性具有不同的影响。2013年陈芳和

22、田凯在降雨入渗作用下土质斜坡稳定性的数值分析一文中，运用有限元软件建立了一个数值模型，分析降雨入渗对斜坡稳定性影响。经过分析得到降雨入渗会使斜坡中的非饱和土变成饱和土，减小了土壤的基质吸力，降低其抗剪强度从而影响斜坡的稳定性。第三节 存在的问题和展望 （1）虽然关于降雨对斜坡稳定性影响的研究成果很多，但是仍有些机理没有完全弄清楚，很多文献都针对研究降雨入渗引起土壤基质吸力的变化导致的斜坡失稳，而各个不同研究方向的耦合分析是研究降雨型滑坡非常重要的一点，但是关于斜坡内具体的物质组成等一些岩土体的变化研究得不够深入。（2）二维饱和非饱和渗流控制方程、含水量变化方程等一些分析地下水渗流方面的研究方程

23、或者理论很多都是建立在许多假设条件下对实际情况进行的分析。（3）Geostudio2007有限元分析软件能耦合分析降雨和斜坡稳定性关系，由于计算机的飞速发展，以后完善的数值模拟技术和坚实的理论基础的结合将能更好地解决许多实际问题。第四节 研究内容及技术路线1.4.1 研究内容国内外许多学者对降雨与滑坡的关系进行了卓有成效的探索，他们在实验室中建模，使用数值模拟方法和在野外进行调查，研究说明了斜坡稳定性和降雨之间有一定的联系，并且针对不同的方面探讨了降雨对斜坡稳定性的影响。例如：许多非饱和土滑坡经常在降雨季节发生；降雨引起斜坡失稳破坏是由于雨水入渗到土壤中，改变了颗粒间的摩擦力和有效应力；降雨强

24、度的不同、降雨持续时间的不同等对不同渗透性系数的斜坡会产生不同的影响。这些研究成果为后人进行更深入的研究奠定了坚实的基础，本文主要研究了以下内容： (1) 通过对研究区前期资料的整理收集和成果汇总分析，掌握研究区地形地貌、水文、气象、构造、第四系土层结构等相关基础资料，分析外部环境和含水介质的特征。(2) 对南平市西芹镇兴华街道后山滑坡区房后斜坡的资料进行相关分析，完成斜坡变形宏观分析及敏感性因素的分析，并使用简化Bishop法计算得到天然条件下的斜坡稳定性系数。(3) 根据资料设计不同工况，运用Geostudio2007软件建立降雨入渗模型，得到不同降雨条件下孔隙水压力随时间的变化，计算不同

25、降雨条件下的斜坡稳定性系数，并对其进行对比分析，进而得出降雨会影响斜坡稳定性的结论。1.4.2 技术路线本人在前人研究的基础上，搜集了南平市西芹镇兴华街道后山滑坡的资料，进行了区房后斜坡变形的宏观分析和敏感性因素的分析，并选取某一典型斜坡地质剖面（本文为28-28 剖面），利用简易Bishop法计算了此剖面在天然条件下（降雨量为零时）的斜坡稳定性系数。之后通过使用Geostudio2007软件建立数值模型，耦合SEEP/W非饱和土渗流模块和SLOPE/W斜坡稳定性模块得到不同条件下雨水入渗时孔隙水压力的变化，并在SLOPE/W模块中同样使用简易Bishop法计算不同降雨强度和不同降雨持续时间下

26、的斜坡稳定系系数，再通过对比不同条件下孔隙水压力随时间的变化来分析降雨是如何对斜坡稳定性产生影响，并给出几条今后的研究建议。收集研究区有关地质、地貌、及自然地理条件等资料获取地下水水位线等相关要素的资料进行斜坡变形的宏观分析及敏感因素的分析计算天然条件下斜坡稳定性系数选取不同降雨条件的工况使用软件对所选滑坡区进行建模计算得到不同工况下的斜坡稳定性系数得到不同降雨条件下孔隙水压力随时间的变化图分析降雨对斜坡稳定性的影响图 1-1 技术路线图第二章 自然地理及地质概况第一节 自然地理条件2.1.1 地理位置西芹镇隶属于福建省南平市延平区，位于延城中心区10公里，人口7.1万。地理坐标为北纬2651

27、-2650，东经11750- 11820，境内自然资源丰富，交通方便， 205、316国道和京福衢南高速公路及鹰厦铁路纵横交汇。经济持续稳步发展，2010年，全镇实现农业总产值4.35亿元；工业总产值9.48亿元。图 2-1 西芹镇地理交通图2.1.2 气象水文该区域属中亚热带季风气候，全年温暖潮湿，四季分明。根据本地区多年资料，该区雨量充沛，属我省丰水区，年平均降雨量1656mm，最大年降水量2255mm，最小年降水量1131mm，过程最大降雨量为252mm。按降水量划分，23月为春雨季，46月为夏雨季（又称梅雨季），约占全年降雨量的52%，降雨量最多的是6月中旬，平均为141.2毫米。79

28、月为第二雨季，主要是台风雷雨，雨日30天左右，雨量340毫米左右，多年月平均降雨量134.08毫米，占年雨量8.3%。10月次年1月为旱季。多年平均气温7月最高，达28.4，最低位1月，为10.2，年平均气温为19.2。第二节 斜坡基本特征2.2.1 斜坡地形地貌 区房后斜坡位于路侧开挖斜坡往南约75m 4幢民房后方开挖坡体，斜坡总长约55m，开挖成台阶状，最顶部标高108.5m-110.8m，坡脚标高79.10m，高差约30m。坡顶为兴华街道山体后山，山顶坡面平缓延长，坡度约为20，坡面共分6级开挖而成，每级高2-7m不等，坡度40,坡脚4幢民房距离开挖坡面2-5m，3-5层，砖混结构，民房

29、前方为兴华街道。图 2-2 区房后斜坡2.2.2 斜坡空间形态区房后斜坡在2010年持续暴雨期间发生坡面崩塌，表层残坡积粘性土及坡面灌木一并下滑，崩塌体总宽度约60m，上下高差达40m，崩塌体厚1-5m，估算下滑土方量约4000m3。崩塌体直接危及到坡脚4幢民房内居民的安全。目前房后坡面分6级开挖成台阶状，每级高2-7m不等，总体坡度40。 2.2.3 斜坡物质组成及结构特征根据野外对区房后斜坡现有槽探开挖情况看，区房后斜坡崩塌体主要由表层残坡积粘性土组成，厚1-5m不等，坡面灌木随同崩塌体一并滑落。2.2.4 斜坡水文地质 斜坡水文地质条件较简单，以主山脊为分水岭，以各山沟为主要通道，坡面降

30、雨顺地形汇集于各山沟中，由西向东径流，排泄于坡脚。一、场地地下水类型 根据地层岩性及地下水在含水介质中的赋存特征判断，地下水类型以孔隙水和基岩裂隙水为主。孔隙水分布于表层残坡积粘性土中，主要接受大气降水补给，由坡顶向坡脚渗透，含水量受季节、降雨量影响较大。基岩裂隙水赋存于下伏黑云母二长花岗岩节理裂隙和风化裂隙中，主要补给来源为大气降水、上部土层垂直补给与临近基岩裂隙水的侧向补给，多为潜水类型。路侧滑坡西南侧坡脚见有泉眼出露，流量约0.09L/s。1. 松散岩类孔隙含水岩组 场地松散岩类孔隙含水岩组主要为表层残坡积粘性土，地下水以孔隙水的形式赋存于孔隙中，富水性较差，主要接受大气降水补给，短途迳

31、流，由坡顶往坡脚排泄，为弱透水层。2. 基岩裂隙含水岩组场地基岩裂隙水主要赋存于花岗岩的节理裂隙和风化裂隙中，富水性差，主要补给来源为大气降水、上部含水层补给及临近基岩裂隙水的补给，短途迳流，向地势低洼处排泄，为弱透水层。二、地下水水位勘查期间在部分勘探孔内测得地下水位稳定水位埋深为0.50-31.20m，稳定水位标高为114.81-125.92m。水位年变幅约5.2m（没有区房后斜坡地下水水位埋深的数据）。表 2-1 场地地下水水位埋深孔号孔位水深孔号孔位水深孔号孔位水深ZK191.00ZK200.50ZK210.80ZK2313.50ZK2415.0ZK3012.40ZK333.90ZK3

32、716.70ZK3811.40ZK3915.20ZK4115.80ZK4317.80ZK4621.50ZK4731.20ZK4814.0ZK4918.30ZK5016.80 ZK5128.40ZK529.20第三节 斜坡岩土体特征根据钻孔揭露情况，场地岩土体自上而下可分为：残坡积粘性土、全风化花岗岩、散体状强风化花岗岩、碎块状强风化花岗岩、中风化花岗岩。第一层为残坡积粘性土，厚度大约在1.4018.40m左右，场地内均有分布。呈红褐色，成份主要为长石风化的粘粉粒，石英颗粒。粘性一般，韧性中等，干强度中等，无摇振反应，稍有光泽；第二层为全风化花岗岩，厚度大约在2.2013.40m左右，场地内均有

33、分布。呈黄褐色，矿物已基本风化为土状，有残余强度，属极软岩，极破碎岩体，岩石基本质量等级为类；第三层为散体状强风化花岗岩，厚度大约在2.036.90m左右，场地均有分布，民房后方及路侧滑坡坡脚局部有出露。整体呈 灰黄色，褐黄色，呈散体状，大部分矿物成份显著变化，花岗结构，具有泡水易软化，崩解，强度降低的特点，属极软岩，岩石基本质量等级为类；第四层为碎块状强风化花岗岩，厚度大约在0.411.40m左右，场地均有分布,局部民房后方有出露。整体呈 灰白色，矿物成份友长石、石英、云母组成，岩芯呈碎块状，岩石坚硬程度为较软岩，岩石基本质量等级为类；第五层为中风化花岗岩，场地均有分布。整体呈 灰白、浅灰色

34、，由长石、云母、石英、角闪石组成，中粗粒花岗结构，节理、裂隙较发育，岩芯呈短柱状，岩石坚硬程度属较软岩，岩体基本质量等级为级。第三章 斜坡稳定性分析斜坡稳定性评价方法可分为两大类，即定性评价和定量评价。定性评价方法包括成因历史分析法、工程地质类比法、赤平投影作图法等；定量评价方法包括极限平衡计算法、有限元分析法、破坏概率计算法等。本章选取区房后斜坡进行斜坡变形宏观分析及敏感性因素的分析，并使用简易Bishop法计算天然条件下斜坡的稳定性系数。第一节 斜坡变形宏观分析由于地质条件的复杂性和人们认识事物的局限性，工程地质定性评价在斜坡稳定性评价中仍然占有极其重要的地位。区房后斜坡现状：坡顶、坡脚未

35、见变形迹象，但受雨水冲刷影响，坡面可见多条冲沟，且开挖坡体总高约30m，坡面岩土体裸露，仅在坡体下部用简易干砌挡墙进行支挡，受降雨等不利因素影响，中上部岩土体仍存在崩塌或滑坡的可能。如图所示，区斜坡失稳破坏后的斜坡坡度仍然较大，并且坡面岩土体裸露，在降雨尤其是暴雨时，雨水入渗冲刷引起再次的滑坡崩塌可能性很大。 图3-1 区斜坡失稳破坏后第二节 斜坡稳定性敏感因素分析3.2.1 地形地貌兴华街道后山山体属于侵蚀剥蚀丘陵地貌，地势西南高，东北低，高差起伏较大。尤其是原滑坡-泥石流处坡面，大都三面环山、一面出口，两侧及后缘坡度陡，中间坡度较缓，空间上呈向下倾斜的簸箕状、漏斗状，形成有利的汇水条件。降

36、雨时，坡面雨水能迅速汇集到下部沟口，且沟口又直对民房。区房后斜坡开挖呈台阶状，总体高差约30m，但上部台阶开挖坡度57-65，仍无法满足放坡稳定坡度要求。3.2.2 岩土体特征场地覆盖层表面主要由残坡积粘性土构成，土层厚度较厚，多在1.40-18.40m，区最厚6.20m，呈可塑硬塑状态。当遭遇降雨时，雨水渗入坡体，软化各岩土层，受雨水软化后，岩土层物理力学性质降低，尤其表现为C值粘聚力的下降，不利于斜坡的稳定。3.2.3 降雨降雨是影响斜坡稳定的重要外界因素，南平市延平区年平均降雨量约1664mm，最大过程降雨量可超过150mm。在2010年6月中旬遭遇百年一遇的短时性特大暴雨，13天内降雨

37、量达684.8mm。其中典型特大暴雨有6月18日延平大横12点至13点的1小时降雨量达到61.5mm。降雨一方面渗入岩土体，增加了土体的重量，加大下滑力，另一方面受雨水浸泡，岩土体的力学强度降低，从而降低了斜坡的稳定性，也就是在这一轮强降雨后，斜坡发生了滑塌地质灾害。可以说强降雨是此次边坡产生滑塌地质灾害的主要外在诱因。3.2.4 人类工程活动在斜坡区域的勘查中发现对斜坡有产生影响的常见人类工程活动为居民日常耕作和对原始山体坡脚进行大面积的开挖。区域内山体坡面普遍被居民来当日常耕作用地，种植蔬菜、毛竹、灌木等浅层根茎农作物，一方面破坏了山坡的植被覆盖，另一方面日常劳作中需要对山坡表层土体进行松

38、动，致使山坡表层土体孔隙增大，有利于雨水的渗入；对斜坡稳定影响最大的人类活动是对原始山体坡脚进行大面积的开挖，使得山体坡脚形成高陡临空面，破坏了原始山体的应力平衡，不利于边坡的稳定。第三节 天然条件下的斜坡稳定性 本节通过简易Bishop法来计算天然条件下（降雨量为零）斜坡的稳定性并对其进行简单分析。瑞典条分法忽略了土条侧面的作用力，算出的安全系数可能偏低10%20%。这种假定方法不是很精确，它是将不平衡的问题按极限平衡的方法来考虑并且未能考虑有效应力下的强度问题。在同样的边坡条件下，简化Bishop法的安全系数一般比瑞典条分法高67。在此先介绍简易Bishop法的基本原理：假定滑动面是以圆心

39、为O，半径为R的圆弧，从中任取一土条i为分离体，其分离体的周边作用力为：土条的重力Wi，切向力Ti和法向反力Ni；土条侧面作用法向力Pi、Pi+1和切向力Hi、Hi+1。根据静力平衡条件和极限平衡状态时各土条力对滑动圆心的力矩之和为零，整理得到Bishop法的公式为： （3-1）其中： （3-2）图 3-2 毕肖普法条块作用力分析简化Bishop法中土坡稳定系数是指整个滑动面上土的抗剪强度f与实际产生的剪应力之比，考虑了各土条侧面间存在的作用力，但忽略了条间切向力。即： （3-3）将式3-3带入式3-1，于是,式3-1进一步简化为式3-4： （3-4）选取区房后斜坡28-28剖面。此剖面实际长

40、度为280m，以5m为单位剖分土条，每隔5m搜索新圆心，以5m为半径画圆弧滑面。场地岩土体物理力学指标值见下表，再以上述简易Bishop公式来计算，得到此剖面的斜坡稳定性系数为F=1.294，因此区房后斜坡在天然条件下（采用天然参数）处于基本稳定状态。表 3-1 场地岩土体物理力学指标值名称天然重度kN/m3粘聚力C (kPa)内摩擦角坡残积粘土层191620全风化花岗岩层201525散体状强风化花岗岩层213024碎块状强风化花岗岩层213228中风化花岗岩层图 3-3 区房后斜坡28-28地质剖面图第四章 降雨条件下的斜坡稳定性对比分析Geostudio2007软件是加拿大GEO-SLOP

41、E公司开发的一款应用于岩土、环境、交通等领域的仿真软件，其中使用得最多的是SLOPE/W、SEEP/W、SIGMA/W模块。本章内容中，SEEP/W模块用来模拟不同的降雨条件，并得到孔隙水压力随时间变化的曲线图；SLOPE/W模块能够运用多种不同方法计算斜坡的稳定性，并具有引入SEEP/W模块的功能。因此将SLOPE/W模块与SEEP/W模块耦合，计算得到不同降雨条件下斜坡的稳定性系数，并对其计算结果进行对比，同时根据孔隙水压力随时间变化的曲线图分析降雨对斜坡稳定性的影响，并分析是如何影响斜坡的稳定性的。分析降雨是通过影响斜坡的哪些方面如何影响斜坡稳定性分析在不同的工况下，降雨条件和斜坡稳定性

42、系数变化的关系与SLOPE/W模块耦合，计算出不同工况下斜坡的稳定性系数SEEP/W模块模拟降雨功能，同时得到孔隙水压力随时间的变化图 4-1 本章流程图第一节 基本原理4.1.1 SEEP/W的基本原理SEEP/W是一款分析地下水渗流的有限元软件。它不仅可以分析简单的饱和土渗流也可以分析与时间相关的饱和-非饱和土渗流问题，通过分析在时间变化下的孔隙水压力从而可以研究斜坡稳定性和时间的关系。 水在饱和土及非饱和土中的渗流均符合达西定律。与在饱和土中的渗流不同，水在非饱和土中渗流的渗透系数K不再是一个常数，而是一个随着土壤含水率或者土壤基质吸力变化的函数。下式为二维饱和一非饱和渗流控制方程，由水

43、流的连续性条件和达西定律相结合得出。 (4-1)式中 H水头； x,z相对水平距离和高程； t时间； Q流量边界条件； Kx(h)、Kz(h)水平方向和垂直方向的渗透系数，在渗流计算中假设渗透系数各向相等，即Kx(h)=Kz(h)=K(h)； h为用水头表示的基质吸力，当Hz时，坡体处于非饱和阶段，如果h不超过设定的最大基质吸力hmax，h=z-H，否则h=hma；当Hz时，坡体处于饱和阶段，h=0； C(h)比水容重，非饱和条件下为土水特征曲线（SWCC）中单位基质吸力的变化引起的含水量变化： (4-2)式中：C(h)代表非饱和土的储水能力，为体积含水量。4.1.2 SLOPE/W的基本原理

44、SLOPE/W是一款计算斜坡稳定性系数，分析斜坡稳定性的软件。该软件能准确、快速地求解各种不同危险滑动面和最小安全系数，最符合的滑动面则再根据经验和实际情况选取。SLOPE/W采用的是极限平衡法（条分法）。条分法的基本原理为将滑动土体竖直分成若干土条，把土条当成刚体，分别求作用于各土条上的力对圆心的滑动力矩和抗滑力矩，土坡的稳定安全系数Fs为抗滑力矩除以滑动力矩。SLOPE/W中常使用的方法有Ordinary法，简易Bishop法，Janbu法，Morgenstem-Price法，Spence法和GLE分析方法。各个方法的静力平衡条件和土条间力的性质如下表所示。本章使用SLOPE/W计算斜坡稳

45、定性的时候采用的是简易Bishop法。表4-1静力平衡条件方法是否有法向力（E）是否有切向力（X）X-E性质和关系Ordinary法否否没有条间力简化的Bishop法是否水平简化的Janbu法是否水平Spencer法是是固定值Morgenstem-Price法是是变量 用户定义表4-2土条间力的性质和关系方法是否满足力矩平衡是否满足力平衡滑动面形状Ordinary法是否圆弧简化的Bishop法是否圆弧简化的Janbu法否是任意Spencer法是是任意Morgenstem-Price法是是任意 第二节 计算模型 本节首先使用Geostudio2007软件的SEEP/W模块，通过导入CAD中的斜坡

46、剖面图建立计算模型，对不同岩层进行分区并进行参数赋值（主要参数为饱和渗透系数、饱和体积含水量和残余含水量），并添加边界条件和初始条件。最重要的步骤是选取不同降雨条件的工况。通过考虑，选择了如下几种工况：降雨时长均为5天，降雨强度分别为50 mm/d、100 mm/d、200 mm/d；降雨时长均为10天，降雨强度分别为50 mm/d、100 mm/d、200 mm/d。之后再在SEEP/W模块下建立SLOPE/W模块，在SLOPE/W模块中对不同的岩层进行参数赋值（主要参数为天然重度、粘聚力和内摩擦角），最后通过耦合两个模块计算出不同条件下的斜坡稳定性系数。4.2.1 基本模型 选取区房后斜坡

47、28-28剖面图，基于下图所示水文地质概化模型采用二维渗流问题有限元分析程序进行暂态的饱和一非饱和渗流计算。图1为入渗模型的地质概化剖面斜坡，图2为渗流计算采用的有限元网格模型，此层状非均质斜坡计算模型从上到下一共分层残坡积粘性土层、全风化花岗岩层、散体状强风化花岗岩层、碎块状强风化花岗岩层、中风化花岗岩层。网格剖分时大小约为1m，剖分得到有限元网格模型共有节点2218个，单元2178个。 图4-2 区房后斜坡模拟地质剖面图图4-3 边坡饱和非饱和降雨入渗有限元计算模型4.2.2 计算参数一、饱和渗透系数、饱和体积含水量及残余含水量坡残积粘土层中的渗透系数通过现场渗坑实验确定，由于只在区和区做

48、过渗坑实验，区的渗透系数只能够根据区和区的数据来估算。而全风化花岗岩层、强风化花岗岩层及中风化花岗岩层的渗透系数没有确定，因此只能依靠经验值估算。表4-3 场地渗透系数估算值名称渗透系数 m/days坡残积粘土层0.05全风化花岗岩层0.1散体状强风化花岗岩层0.02碎块状强风化花岗岩层0.005中风化花岗岩层0.003残坡积粘性土和全风化花岗岩层、强风化花岗岩层以及中风化层的饱和体积含水量及残余含水量由实验室的土工试验结果得到。表4-4 场地饱和体积含水量及残余含水量名称饱和体积含水量m 3/m3残余含水量m3/m3坡残积粘土层0.60.28全风化花岗岩层0.60.15散体状强风化花岗岩层0

49、.40.16碎块状强风化花岗岩层0.30.11中风化花岗岩层0.30.12二、场地岩土体物理力学指标值 建立SLOPE/W模块时需要给每层赋值，最重要的三个指标是天然重度、粘聚力和内摩擦角。各值见第三章第三节表3-1场地岩土体物理力学指标值。三、各土层的土水特征曲线土水特征曲线的研究非常重要，因此很多学者提出了大量有关土水特征曲线的数学模型，其中以Fredlund&Xing模型及Van Genuchten模型效果较好。1、 Fredlund&Xing模型原理 通过对土体孔径分布曲线的研究，用统计分析理论推导出非饱和土体积含水率与基质吸力之间的对数函数型式的土水特征曲线，其归一化表达式如下： （

50、4-3）式中：为基质吸力；为与基质吸力相对应的体积含水率，其取值范围为r ,s，r 和s分别为残余含水率和饱和含水率；n为归一化的体积含水率(有效饱和度)；a、n、m为拟合参数（a是与进气值有关的吸力值，对应于曲线的拐点；n是与孔径分布有关的参数，对应于曲线的斜率；m是与曲线拐点处的不对称性有关的参数)；e为自然对数的底。2、Van Genuchten模型原理 通过对土水特征曲线的研究，得出非饱和土体积含水率与基质吸力之间的幂函数形式的土水特征曲线，其归一化表达式如下： （4-4）通过阅读文献，Fredlund&Xing模型可较好地表征非饱和粘性土土壤含水量和基质吸力的相关关系,文献中从拟合曲

51、线可知Fredlund&Xing模型和Van Genuchten模型在大部分含水量范围内与试验数据吻合较好，但在接近于完全干燥状况时，Fredlund&Xing模型计算的含水量接近于0，更符合理论情况，而此时Van Genuchten模型的误差较大。因此本文在使用Geostudio2007软件(有Fredlund&Xing模型和Van Genuchten模型这两种选项)设置参数计算时，选择的是Fredlund&Xing模型。以下为各个土层的含水率与基质吸力的关系图（也即土水特征曲线）。 图 4-4 第一层岩土层土水特征曲线 图 4-5 第二层岩土层土水特征曲线 图 4-6 第三层岩土层土水特征曲线 图 4-7 第四层岩土层土水特征曲线图 4-8 第五层岩土层土水特征曲线4.2.3 计算模型中的边界条件和初始条件根据第二章的地下水水位埋深表大概画出相应的地下水水位线。模型的左侧边界为分水岭，将其设置为零流量边界；模型底部为不透水层边界，也设为零流量边界；模型的右侧设置潜在渗流面（potential seepage pace）；斜坡表层设置降雨条件，见图4-3。降雨条件分别按如下几种工况设置：降雨时长均为5天，降雨强度分别为50 mm