库水位升降作用下边坡失稳机理研究.pdfVIP

库水位升降作用下边坡失稳机理研究

摘要：介绍了饱和-非饱和渗流模型的基本理论，并用GRO-SLOPE分析软

件进行边坡失稳机理研究。

关键词：饱和-非饱和渗流模型GRO-SLOPE失稳机理研究

1引言

边坡在水的作用下可能会产生滑坡，“十个滑坡九个水”这句话充分反映

水是产生边坡失稳的重要条件之一。崩塌、滑坡、泥石流等边坡失稳现象的发生

和发展多受水等因素的控制。本文采用边坡稳定性分析软件GEO-SLOPE，计算

了在库水位下降8m期间，滑坡体稳定性随不同库水位下降速度的变化情况。

2基本理论

2.1饱和-非饱和渗流模型

水位涨落下边坡中的渗流场是饱和－非饱和渗流场，非饱和土区水的流

动服从达西定律。达西定律是多孔隙介质中流体流动的运动方程，质量守恒是物

质运动的普遍规律，将质量守恒定律运用于多孔介质中的流体运动即为渗流控制

方程[5]：

(1)

式中：、分别为x、y方向的饱和渗透系数；为介质相对渗透率，在饱和

区为1，非饱和区介于0～1之间；h为总水头，h=hp+y，hp为压力水头，y为

位置水头；Q为源汇项；为容水度，在饱和区为0，非饱和区容水度为,为体积

含水量；β为选项系数，在饱和区为1，非饱和区为0；Ss为饱和土体单位贮水

率，在饱和区为一常数，非饱和区为0；t为时间变量。

2.2饱和－非饱和渗流有限元算法

将渗流计算区域进行单元离散，单元内水头在空间和时间域上的变化近

似地表达为：(2)

式中：m为单元节点数；hm为单元节点m的测压管水头；Nm为单元节

点m的形函数。

由Galerkin加权余量法及格林公式可得求解渗流场的有限元法矩阵方程

为：

(3)

式中：；

；

三.计算模型与参数

3.1计算模型

计算软件选用加拿大SLOPE公司开发的边坡工程分析系列软件包中的

SLOPE/W和SEEP/W软件模块。其中SLOPE/W能够利用极限平衡原理来计算

含水的岩土边坡的安全系数。SEEP/W是一种既能分析饱和渗流又能分析非饱和

渗流的渗流分析通用程序。采用SEEP/W建立三角形单元和四边形单元相结合

的渗流计算网格,分别计算水位上升和下降的情况下，某库水位条件下的库岸坡

体内的渗流场分布，然后将渗流计算结果导入SLOPE/W中计算该水位条件下的

安全系数。

水位涨落差设定为h=18m，边坡土体和基础均为饱和-非饱和材料，渗透

系数为随含水量变化的函数，饱和土体渗透系数为1e-2m/day。基础与边坡土体

材料参数如表1所示，计算模型如图1所示。

3.2运用GRO-SLOPE运算求解

边界与水位涨落速度设置为：边界1、边界2为总水头随时间变化函数

并设置为可渗流边界，水位降落速度为2m/day。荷载考虑土体自重以及库岸水

压力的作用。路堤边坡失稳时剪入口位于边坡坡顶，剪出口位于边坡坡脚，滑动

体的土条数为30，迭代次数为2000，安全系数采用Morgenstern-Price求解，精

度为0.001。

四.求解和分析

4.1库水位下降

水位下降时安全系数随时间变化柱状图如图2所示。从图2可以看出，

安全系数随着水位的下降而逐渐降低，但是降低到一定程度以后又有所回升。安

全系数从第1天的水位高度18m开始下降，到第3天末水位高度12m为最低，

第四天开始安全系数逐渐缓慢增大。其原因主要是水位正常下降时，边坡内部渗

流场将随着边界条件的变化而进行不断调整，这一过程伴随着坡体内部的地下水

位下降。另外，所建模型渗透系数相对较大，当库水位缓慢下降时，坡体内部的

地下水位近似与库水位同步下降。与此同时，库水位的下降将对坡面施加水压力

作用，其大小为hwγw（hw为库水位深度，γw为水容重），水位下降时，水流

在土体内部额外增加了一个向外的渗透力，容易导致路基的失稳破坏。一方面坡

体的浮托力减小，另一方面由于坡体内孔隙水来不及排出而使坡体水位高于库水

位，由此产生的渗流作用是潜在滑动面的抗滑能力降低，