基于非饱和渗流理论的膨胀土岸坡稳定性分析.docxVIP

基于非饱和渗流理论的膨胀土岸坡稳定性分析 水位的快速变化对水库边缘的稳定性，尤其是对含膨胀土的通道边坡和岸壁的稳定性有显著影响。水位上升导致膨胀土的坍塌收缩，并改变坡体中的渗流场，这改变了岸壁的稳定性。特别是在一些差距发展的条件下，这是显著的。国内膨胀土岸坡、渠道因水位变化而导致失稳的实例较为常见。对于干湿循环条件下裂隙开展对膨胀土强度特性的影响,刘华强结合镇江南徐公路进行过研究;而对于水位升降时滑坡、边坡的稳定性,一般采用饱和稳定性渗流方法计算。实际上水位升降条件下,膨胀土岸坡体处于非饱和与饱和变化状态,采用非饱和渗流分析理论更接近实际情况。 笔者以含裂隙膨胀土岸坡为例,以非饱和渗流理论为基础,基于河流水位变化条件,运用有限元数值分析软件的SEEP/W模块进行水位升降条件下的瞬态渗流场模拟,研究了坡体内渗流场的变化规律,同时以渗流模拟得到的瞬态孔隙水压力为基础计算岸坡在水位升降和裂隙开展深度下的稳定系数,建立岸坡稳定性变化曲线,以期为库岸的治理和防治提供依据。 1 风化基岩带 以文献中临河岸坡外部条件为参考,建立图1所示的岸坡二维地质模型。岸坡上部为棕黄色粉质黏土(Qeol44eol),厚29.0～34.1 m,具有弱膨胀性,其中表层为强风化土层,发育宽0.1～0.3 cm、深0.5～6.0 m的网状裂隙;下部为粉沙质泥岩,产状近水平,上部3～6 m的为全风化层,其下岩体结构完整,极微透水,泥岩具有软化性,遇水易崩解,在水和空气交替作用下,风化剥蚀沿纵向延伸,上部土层极易沿风化面产生滑动破坏;河面静水位为11.42 m,岸坡地下水位埋深为14.2～21.0 m。 2 水位变化下岸壁梯度的渗透分析 2.1 土的渗流方程 水位的升降使岸坡体中地下水位发生变化,随着浸润线的变化,在岸坡土体中形成了非饱和区和饱和区。水在非饱和-饱和土中的渗流服从Darcy定律,在各向同性土中,渗流方程为 ??x(kw?hw?x)+??y(kw?hw?y)=mwγw?hw?t(1)??x(kw?hw?x)+??y(kw?hw?y)=mwγw?hw?t(1) 式中:kw为渗透系数,对于非饱和土,其为基质吸力或体积含水量的函数;hw为渗流水头;mw为土水特征曲线的斜率;γw为水的密度;x为水平方向;y为竖直方向;t为时间。 2.2 土水特征曲线 定义最左端节点为变水头边界,最右端为定水头边界,将不稳定状态作为岸坡稳定性的判别依据,选取汛期连续集中降雨引起河水位急剧上涨到21.42 m,水位在21.42 m运行2 d,其后降落到初始水位这一过程为研究对象,耦合分析中选取洪水位升降速度为5 m/d,同时不考虑动水条件。水位涨落曲线见图2。 该岸坡膨胀土与文献中所研究的弱膨胀土的性状相似,通过对文献中土水特征曲线参数进行调整,建立了土水特征曲线(图3),该曲线能够满足精度要求。渗透性函数采用FredlundXing模型进行估计,取饱和渗透系数为3.73 m/d;泥岩全风化层位于初始水位线以下,取饱和渗透系数为6.7 m/d;全风化层以下基岩完整,渗透系数取0.086 m/d。 2.3 最大水位下的各点孔隙水压力 为了研究水位升降过程中岸坡内外孔隙水压力的变化,进而研究渗流对岸坡内部土体的影响,在坡面选取1、2、3、4、5点,同时在岸坡内部选取a、b、c、d点(见图1),以5 m/d升降速度下的暂态渗流场模拟为基础,绘制孔隙水压力随时间变化曲线(见图4、图5)。 由图4、图5可以看出,在水位升降过程中,位于初始水位下的5点的孔隙水压力涨消和水位变化曲线一致,最终孔隙水压力恢复到初始状态;位于初始水位以上、最高水位以下的4点、d点孔隙水压力变化基本上和水位变化曲线一致,在水位恢复到初始水位时,集聚的超孔隙水压力最大,并且超孔隙水压力消散得最慢;位于最高水位以上的3点、c点孔隙水压力涨消与水位的涨落不同步,明显滞后于水位的变化,当水位达到最高水位21.42 m时,孔隙水压力才有明显增加,当水位恢复到初始水位时,也集聚一定的超孔隙水压力,不过在数值上没有4点、d点的大;远离最高水位的1点、2点、a点、b点孔隙水压力没有变化,没有产生超孔隙水压力。 库岸边坡既有一般山地边坡的共性,又有其特殊的一面,其活动与水位的升降有很大的关系,水位急剧升降使岸坡的渗流场发生改变,影响土体的稳定性。渗流场的变化对岸坡内部不同部位的土体影响是不均等的,有着明显的分带性。以岸坡内部a、b、c、d点孔隙水压力变化为参考,以超孔隙水压力为依据,对基岩全风化岩体以下不予考虑,近似将岸坡内部划分为Ⅰ、Ⅱ、Ⅲ区,见图6。远离河水位的岸坡后缘Ⅰ区土体基本不受水位的影响;在最高水位以上的Ⅱ区为受渗流影响明显区,孔隙水压力的增减和水位的升降不同步,在水位降落到初始水位时,产生的超孔隙水压力还没有完全消散