心墙坝渗透稳定分析.docxVIP

心墙坝渗透稳定分析 1 广东飞来峡水利工程 在我国的中小型水利水电工程中，心墙粘土水库和心墙土石坝是常见的水库类型。对水库进行水位分析和坡面稳定的评估是安全评估的重要内容之一。它的目标是确定不同区域水库和不同土层之间的渗透、稳定和渗透流量。当坝基是软土地基时,由于软土的低强度和大变形亦可能成为控制大坝稳定的重要因素,所以,对于这一类大坝的安全评估必须将坝体和坝基视为整体,以全面分析各种工况下大坝在自重和渗流的耦合作用下的渗透稳定性。 广东飞来峡水利枢纽位于北江干流中游,主要由泄水和挡水坝段、船闸、主土坝和副坝等组成。其中右岸主土坝全长1 778 m,坝顶高程34.8 m,最大坝高18 m,坝顶宽为8 m。该主坝跨越不同的地质结构,大约有400 m长的坝段基于软黏土和淤泥质黏土,施工时分别采用开挖换土和砂井排水固结的方法进行了地基处理,坝体填筑时采用控制填土速率、加设反压平台的方法保证其稳定。 由于该坝段在施工过程中曾由于软基沉降产生过坝体裂缝、在运行过程中也发现坝体有不同程度的沉陷趋势,为对大坝的安全性作出总体评价,本文根据大坝实际运行的上下游水位和坝体应力状态,选定了多种工况作为计算内容对该主土坝典型断面进行了渗透稳定分析。 2 分析模式 2.1 渗流微分方程解析解的方程法法传统方程,在较 稳定渗流微分方程为 式中:kx、kz分别为土体水平、铅直方向的渗透系数;h为土体中任一点的水头高度。结合式(2)所示的边界条件,即可求出式(1)的定解。 式中:h*为自由面边界上的水头,为保证式(1)的解具有唯一性,应满足: 式(1)所示的渗流微分方程对于复杂的边界条件和非均质介质较难求得解析解,故在实际运用中往往采用数值解。根据变分原理和泛函求极值的方法可得到上述渗流微分方程的有限元计算公式为 式中:[K]为与渗透系数kx、kz有关的矩阵;{f}为已知常数项,由已知水头结点得出。 2.2 土体稳定验算 渗流作用下,土坝的破坏原因可分为两种:一种是因出逸坡降导致的土体渗透变形(管涌或流土);另一种是因渗流场存在的孔隙水压力导致的土体的滑动,所导致的后果有局部失稳定和整体失稳。 如果定义Ms为抗滑力矩,Mτ为滑动力矩,则根据塑性极限平衡原理,可得土体稳定安全系数η(包括局部稳定和整体稳定)的定义式为 依据圆弧滑动和条分法思路,可推导出考虑渗流作用力的坝坡稳定有限元计算公式为 式中:c′、?′为土体有效应力指标;l为土条底部弧长;?为所划分的单元面积(具体计算时采用了三角形单元);γ为土体重度,浸润线以下取有效重度;Jz和Jx分别为单元在z、x方向的渗透坡降,其定义式如式(7)所示,也可用其结点的坐标和水头值表示。 其余几何符号的意义见图1。 采用式(6)进行坝坡稳定分析实际是根据渗流分析结果,将作用于各结点的压力水头值换算成各单元渗透力(图1中的T或Fx和Fz,T=γwJ),分析中不再计算边坡外所受的静水压力,也不计算各单元体接触边界上的孔隙水压力。 2.3 土体渗透系数随孔隙率的变化规律 式(1)所示的渗流方程中渗透系数被假定为常数,不随固结过程发生变化,这对于软土、特别是压缩特性比较明显的软土来说是有较大误差的。如果考虑压缩性对渗流场的影响,可将式(1)改写为 式中:kx(σij)、kz(σij)分别为土体水平、铅直方向的渗透系数,与压缩应力σ有关。 根据现有研究,可采用下述经验关系描述渗透系数随孔隙率变化的规律: 式中:k0为初始渗透系数;α为试验常数,对于软黏土,可取α=4.0~6.0;e为孔隙比;n为孔隙率。 假设土体的体积应变全部产生于孔隙体积应变(饱和土体),有: 则式(9)可表示为 由于εv与应力有唯一的关系,所以,式(11)所示的渗透系数最终可以表示为应力的函数,即k=k(σij)。将式(11)代入式(4),即可进行考虑软土层渗透特性变化的心墙坝渗流分析。 3 计算当前和参数 3.1 计算噪声 根据大坝的设计标准和实际运行情况选定了13种计算工况,各工况对应的上下游水位见表1。 3.2 原状土层土性指标 软基上心墙坝的典型断面见图2(主土坝B2+330),根据不同的土层将断面划分为18个区域,各区域土性指标通过原状土的室内试验获得,参见表2。 4 计算与成就分析 4.1 心墙下游逸出位势随水位的变化 假设土体的水平渗透系数kx与表2中kz相等,计算各种工况下心墙坝典型断面的渗流特性(图3),可以发现心墙的防渗效果与上下游水位有关: (1)当上游水位低于心墙高程时,在工况1所示的低水位时(图3(a)),心墙下游逸出位势接近0,表明心墙几乎削减了全部水头;随着水位的增加,心墙下游逸出位势逐渐增大,在工况2所示的水位时(图3(b))削减水头约为70%; (2)当上游水位高于心墙高程时,心墙下游逸出位势一般在30%~60%,同等