剩余推力法在边坡稳定性评价中的应用.docxVIP

剩余推力法在边坡稳定性评价中的应用 1 剩余推力法的计算困境 剩余动力法是最广泛的工程应用的不确定性的边坡稳定评价方法。经过90多年的国内和国外专家和科学家的努力，剩余动力法得到了不断改进和完善。但由于工程地质条件的复杂性和剩余推力法本身理论上的欠缺,其计算有时并不十分合理,尤其是在处理一些动态因素的影响时,这些局限性益发显露出来。地下水位动态变化,就是一个典型的影响斜坡稳定性的动态因素,但在剩余推力法中还没能很好的考虑。尤其对于水库斜坡来说,在人类活动的影响下,水库水位的变化可能是影响地下水位的动态变化的主导因素。本文基于上述思路,先对库水位变化引起地下水位的变化进行研究,将地下水位在水库蓄水和下降时的浸润线方程引入到原剩余推力法的计算公式和程序中,使之能考虑水位波动的影响,并结合一个具体的工程实例说明其应用。 2 渗流压力变化 以往一般都以静态的观点考虑水压力,即根据地下水位的分布,计算斜坡的稳定性系数。而渗透压力实际上反映了地下水渗流在渗透过程中,总水头损失的那一部分孔隙水压力转化为作用在水流方向上的有效压力。因此,有必要对地下水位对库水位波动的响应进行研究,并找出岸坡可能失稳的最不利时刻。 2.1 库水位上升时地下水位浸渍线的计算 水库蓄水后,低于蓄水位的斜坡体内的地下水位将随水库蓄水位的高低、某种水位运行时间的长短等因素而产生相应的变化。根据《重庆市三峡库区滑坡危岩勘查规定(试行)》建议采用下列的公式计算库水位上升时滑体中的地下水位线(浸润线)(见图1,a)。假定斜坡的隔水底板为水平,当库水位上升时,地下水位浸润线可按下式计算。 y=y2p+s′(h2?h2p)/s??????????????√(1)y=yp2+s′(h2-hp2)/s(1) 式中,yp为水位上升后的库水位(从隔水层顶面起算)(m);hp为水位上升前的库水位(从隔水层顶面起算)(m);s、s′为水位上升前、上升后计算点至库岸距离(m);y为水位上升后计算点的地下水位(从隔水层顶起算)(m);h为水位上升前计算点的地下水位(从隔水层顶面起算)(m)。 滑体内的地下水位的响应与物质的渗透性有很大关系,比如渗透系数K10-3cm/ s时,地下水位的上升基本上与库水位的上升同步,并且与库水位持平。如果滑体物质的渗透性较小,滑体内的地下水位将随时间而变,但只要某一库水位维持的时间足够长,滑体内的地下水位也将与这个库水位持平。从这点来看,库水位不再上升的那一时刻渗透力最大,也是对斜坡稳定性最不利时刻。 2.2 渗流自由面位置最高的确定 库水位下降,有两种情况:①库水位骤降;②库水位缓降。所谓水位骤降,一般是指水位降落很块,以致在降落过程中没有显著的水量从土体中排出,土体仍全部饱和,斜坡内自由面或浸润线基本保持不变。所谓缓降,是指在水位的降落过程中,自由面有显著下降。 一般来说,库水位骤降,对斜坡来说是最危险的。但是最危险水位的确定仍然是存在的一个难题。Neuman教授(1974)采用Galerkin加权残值法对二维饱和-非饱和渗流问题进行了大量地计算分析和深入的研究,取得了较满意的成果。彭良泉等(2003)采用了Neuman教授的分析方法研究三峡库区水位降落时的水库斜坡的渗流情况提出了渗流自由面位置最高点确定的经验计算公式。 h0H=1?(0.28?0.33)tT(kμv)14(2)h0Η=1-(0.28-0.33)tΤ(kμv)14(2) 式中,k,μ为土的渗透系数和给水度;v为水位的降速;h0为自由面最高点水头;H为库水头的最大降距;t为库水位降低的任意时刻;T为库水位降落所需要的时间。式(2)表明,若已知斜坡组合参数k/(μv),则h0/H与t/T就是直线关系, 如此, 就能把库水位降落停止在任意所需要的位置上。 渗流自由面下降过程中的最高点位置h0确定以后,就可以根据巴甫洛夫斯基-达赫勒分段法计算渗出点E的高度He和自由面的位置(见图1,b),从最高点到渗出点间坡身段的渗流量为: qk=(h1+h0)2?h2e2(L?m1he)(3)qk=(h1+h0)2-he22(L-m1he)(3) 由渗出点E到坡脚C的坡面渗流量为: qk=(he?h1)m1(1+Lnhehe?h1)(4)qk=(he-h1)m1(1+Lnhehe-h1)(4) 式中,h1为计算下降水面时的水深;m1为坡率;L为计算的自由面最高点到坡脚的水平距离。使上两式相等, 求出he, 然后算出q/k, 再利用式(3)可以导出离开自由面最高点的水平距离x处的自由面高度h: h=(h1+h0)2?2x(q/k)????????????????√(5)h=(h1+h0)2-2x(q/k)(5) 由式(2)可知,当t=T时,库水位降落所引起的浸润面最陡,也即所动水力坡度最大,所以此时应是对斜坡稳定