水位反弹下岸坡内渗流场分布及稳定性分析.docxVIP

水位反弹下岸坡内渗流场分布及稳定性分析 山坡是一种常见的自然灾害。中国每年因山坡原因造成的经济损失超过100亿元。统计表明,我国有新老滑坡约30万处,其中灾害性滑坡约1.5万处。据不完全统计,三峡库区在175m库水位影响的范围内共有大小滑坡1190余个,各类变形体更是广泛分布库水位骤降对岸坡稳定性会产生较大影响。由于国家经济和社会发展的需要,我国将在今后一段时间将水电站建设列为我国能源建设中的一个重要工作。地下水在岩体内流动,产生的渗透压力会改变岩体内的原始应力状态,使岩体的物理力学性质发生改变,从而影响到边坡的稳定性。 渗流对边坡稳定性具有显著的影响,因此在进行边坡稳定性分析时需考虑渗流的作用。毛昶熙等以三角形单元的渗透力取代条分法中垂直土条周边的水压力,提出了圆弧滑动有限元法;汪自力等以土体单元渗透力取代其周边孔隙水压力进行饱和—非饱和渗流作用下边坡稳定性分析;唐晓松等应用PLAXIS程序采用有限元强度折减法分析了水平自由面渗流作用下边坡的安全系数。 当库水位骤降时,由于土具有一定的持水性, 引起岸体中浸润线的瞬变,岸体内的自由水面或者浸润线高于外坡的库水位,因此具有较高的孔隙水压,同时受指向库内的渗透力作用,从而加速了岸坡的失稳破坏。本文结合有限元数值模拟软件GEO-SLOPE以及SEEP/W对均质岸坡在库水位骤降条件下进行岸坡稳定性分析,为库岸灾害防治提供参考依据以及指导。 1 渗流控制方程 在库岸工程实践中,经常发现部分潜在滑动面可能处于边坡地下水位线以上,即潜在滑动面可能通过具有负孔隙水压力的非饱和土;另一方面,由于库水位以上滑体处于非饱和状态,库水位以下处于饱和状态,随着库水位快速下降,岸坡土体内孔隙水压力也逐渐消散,土体由饱和变为非饱和状态,而非饱和区水的运动和饱和区水的运动是相互联系和影响的,将两者统一起来即所谓的饱和非饱和问题。 对于各向同性多孔介质,渗流控制方程为: ??x(kx?H?x)+??y(ky?H?y)+Q=mwγw?(H?y)?t(1)??x(kx?Η?x)+??y(ky?Η?y)+Q=mwγw?(Η-y)?t(1) 式中kx,ky,分别为x和y方向的导水率;Q为边界上水的流量;mw为比水容量;γw为水容重;该方程描述了渗流区单元内某时段流入和流出的总净流量等于体积含水率的变化量。 其边界条件为: (1) 水源限制 k?H?n|r1=h(x,y,t)(2)k?Η?n|r1=h(x,y,t)(2) (2) 初始条件的确定 k?H?n|r2=q(x,y,t)(3)k?Η?n|r2=q(x,y,t)(3) 初始条件通常为在开始时刻的渗流场水头分布。在进行非稳定渗流计算时,可先求得开始时刻稳定流场的水头分布,来作为已知的初始条件;也可取任一时刻的渗流状态作为初始条件。 求解方程,有限元程序采用的是伽辽金余量加权控制方程。利用二维渗流有限元分析程序GEO-SLOPE以及SEEP耦合可求得库水位骤降前后的岸坡的稳定性,并可对稳态和瞬态饱和非饱和渗流进行模拟。 2 分析均质岸边坡的渗透 2.1 研究计算模块 GeoStudio 是加拿大GeoStudio公司开发的一套功能强大、适用于岩土工程和岩土环境模拟计算的软件,在地质构造、土木工程、采矿工程、岩土工程、地下水分析中得到了广泛应用。SEEP/W二维渗流分析计算模块,可以通过渗流有限元模拟计算,对岸坡单宽剖面渗流路径、不同位置的流速、等水头线、通过指定流量截面的渗流量和各个节点的孔隙水压力值等问题进行分析。通过瞬态分析,可以得出不同时刻不同点的孔隙水压力状况,从而可对岸坡渗流进行分析。GeoStudio的各模块间可以结合运用,分析结果可以相互调用。同时结合GEO-SLOPE程序对岸坡的稳定性进行耦合计算。 2.2 织物初始浸渍线、岸坡系数变化规律 本文渗流分析使用有限元程序SEEP/W,利用土水特征曲线与渗透性函数模拟非饱和岸坡中的渗流情况,同时结合GEO-SLOPE耦合分析岸坡在库水位骤降过程中浸润线的变化以及岸坡安全系数的变化。数值模型如图1所示,岸坡体底宽55 m,顶宽20 m,岸高18 m,岸体边坡坡度为1∶1.82.岸体上游初始水位16 m,均质土石岸坡体内初始浸润线如图所示。岸坡面为变水头边界,通过介质运移模式来模拟库水位骤降的过程。 当库水位下降时,岸坡体内浸润面上的土体变为不饱和。随着基质吸力的进一步增加,水占土体空隙体积进一步减少。孔隙气压力假设为大气压,基质吸力值在数值上等于孔隙水压力,即可用孔隙水压力水头表示,将岸底选做基准面。 2.2 不同水位条件下的岸坡稳定性 对岸坡进行渗流分析,一般以稳态渗流分析的结果作为暂态渗流分析的初始条件,这与实际工程中的处理相似。库水位保持一定水位,在较长的一段时间里,岸体内渗流趋于稳定状态。稳态渗流