张沟水库上游渗透稳定性分析.docxVIP

张沟水库上游渗透稳定性分析 1 均质土石坝坝水岸工程 张沟水库位于陕西省罗南市城关市七岐口村南部的张沟。水库距离罗南市约5.0公里。张沟水库以灌溉和供水为主,兼顾防洪。主体工程由均质土石坝、放水洞、溢洪道组成,坝高15.0 m,坝顶高程为99.80 m,坝顶长91.5 m、宽2.0 m,正常蓄水位97.40 m,设计洪水位98.45 m,校核洪水位99.07 m,下游水位85.50 m,总库容40万m3,属小(2)型水利工程,主要建筑物等别为四级。水库投入运行以来,上游坝面塌陷,下游坝坡蠕变、拉裂、鼓凸明显,右坝肩渗水严重,威胁下游七岔口村安全。 2 计算方法 2.1 渗透系数及流量 大坝渗流计算采用有限元法进行。由于大坝运行多年,可认为土体完全固结,因此采用不可压缩渗流方程。符合达西定律的二维均质土石坝稳定渗流方程为 ??x∣∣kxx?h?x+kxy?h?y∣∣+??y∣∣kyx?h?x+kyy?h?y∣∣+Q=0(1)??x|kxx?h?x+kxy?h?y|+??y|kyx?h?x+kyy?h?y|+Q=0(1) 式中:h为水头;kxy、kxx、kyy、kyx为以x、y轴为主轴的渗透系数;Q为流量。 水头边界及流量边界分别为 h=h0(2)lxvx+lyvy=vn(3)h=h0(2)lxvx+lyvy=vn(3) 式中:h0为已知水头;vx、vy分别为沿x、y方向的流速;vn为沿合向量方向的流速;lx、ly分别为沿x、y方向的渗透坡降。 稳定渗流场的计算可以通过求解下式得到: kh=F(4)kh=F(4) 式中:F为已知常数项,由已知水头节点得出;k为渗透系数。 求解出各节点的h后,可进一步求得各节点的流速和渗透坡降,在渗流场的自由面上,压力p=0,且无渗流量通过,因此还必须满足h=y、vn=0。 2.2 滑裂面法在坝坡稳定性中的计算 根据《碾压式土石坝设计规范》(SL274—2001),采用Bishop法对大坝进行抗滑稳定计算。按照有效应力圆弧滑裂面法计算大坝边坡稳定性,浸润线数据采用大坝渗流稳定分析的结果。有效应力圆弧滑裂面法就是把孔隙压力作为滑裂面上的作用力计算坝坡稳定性,其公式为 K=∑[Cb+(W?U)tanφ]1mθ∑[Wsinθ+Dcos(α?θ)](5)mθ=cosθ+sinθtanφK(6)Κ=∑[Cb+(W-U)tanφ]1mθ∑[Wsinθ+Dcos(α-θ)](5)mθ=cosθ+sinθtanφΚ(6) 式中:K为安全系数;C、φ为土的抗剪强度有效应力指标黏聚力和内摩擦角;W为条块重力,浸润线以上取重度,浸润线以下取饱和重度;U为条块所受到的浮力;D为条块所受到的渗透力;α为条块渗透力与水平线的夹角;θ为条块重力线与通过此条块底面中点半径之间的夹角;b为土条宽度。 3 计算模型、场景和参数 3.1 坝体的设计深度 取竖直向上为Y轴正向,顺水向为X轴正向,原点位于0高程水平面上。根据参考文献,地基的有效深度取坝基长度的0.5～1.0倍。考虑张沟水库坝基的地质情况,计算时,沿坝基截取20 m深度计算,并分别沿坝体上游、下游河床截取20 m为计算区域。渗流计算有限元网格节点总数为786个,单元总数为1 650个。 3.2 坝体自然工况设计 渗流计算分别考虑正常蓄水位、设计洪水位、校核洪水位3种工况。根据规范和坝体实际情况,确定上游坡计算工况为正常蓄水位稳定渗流期、正常蓄水位+7级地震、校核洪水位骤降至正常蓄水位,下游坡计算工况为设计洪水位稳定渗流期、设计洪水位+7级地震、校核洪水位稳定渗流期。 3.3 计算参数 根据《张沟水库工程地质勘察报告》,坝体、地基的渗透系数和材料物理力学指标见表1。 4 结论分析 4.1 浸渍线和渗透坡降 根据以上计算原理、边界范围以及选定参数,采用二维有限元软件对张沟水库的3种工况进行渗流分析,得出坝体水力梯度、流速矢量图以及坝体的单宽渗流量和最大坡降,0+024(最大坝高断面)渗流计算结果见表2。 由计算结果可知,坝体单宽流量不是很大,最大坡降都在允许范围内,能满足抗渗要求,但张沟水库0+024(最大坝高断面)正常蓄水位、设计洪水位和校核洪水位3种工况下浸润线出溢点在下游坝体高程90 m附近,均比较高,若不能有效进行保护,则存在渗透破坏的可能。浸润线和渗透坡降较高的原因是,水库受到严重的人为破坏,内坡反滤层坝面全部被揭,下游坝外坡被耕种,坝脚排水棱体遭到破坏,失去排水能力;坝体填筑质量较差,压实度合格率为0,渗透系数大。若无良好的排水设施,则必然导致坝体内浸润线较高。 流速矢量图表明,中游、下游坡脚流速较快,出溢点高程90 m左右应重点进行防护。最大水力梯度位于坡脚处,最大值为0.42,说明下游坡坝体内很大区域是饱和的,对坝体的稳定极为不利,同时下游也有可能出现渗透破坏。