新集水库均质土坝渗流稳定性分析.docxVIP

新集水库均质土坝渗流稳定性分析 水库是水库中心的主要建筑之一，主要用于水库的蓄水。由于地质条件要求低、成本低，因此已被广泛使用。在土坝研究中, 土坝的安全一直是人们关注的重点, 而渗流是威胁土坝安全的主要因素之一, 因此土坝的渗流问题越来越引起人们的关注。水库在运行期间, 当遇到当地雨季来临或者暴雨相对集中的时候, 为了水库的安全运行, 就要根据库容量实时泄水, 保证水库的安全, 此时库水位就会以不同的速度下降。水库水位下降会改变坝体正常的渗流状态, 不同的水位下降速度下, 渗流场的变化也不同, 坝体渗流场的改变, 会影响上游迎水坝坡的稳定性, 尤其当水库长期蓄水之后水位的下降, 这种影响更加明显。如1968年湖南省流光岭水库均质土石坝由于水位降落, 坝体内渗流场变化导致上游坝坡严重滑坡;1970年10月, 江西省七一水库粘土心墙坝, 由于在坝体右岸开挖引水隧洞的影响, 致使库水位开始快速的下降, 进而改变了坝体内的渗流场, 最终导致迎水面坝坡发生严重的破坏;1979年湖北省狮子岩水库土坝上游坝坡受库水位骤降影响产生滑坡;1972年—1981年间福建省发生的一些土石坝滑坡, 其中有12座是由库水位降落引起的, 在美国、日本30%~40%的坝体滑坡也是由于库水位的下降而引起的, 因此对于土坝有必要进行库水位下降工况下坝体内渗流场的分析, 进而为坝体稳定性分析以及安全评价提供参考。 本文以新集水库均质土坝为研究对象, 考虑了土坝粉质粘土的渗透系数、体积含水率与基质吸力的函数关系, 利用有限元数值分析软件Geo-slope中的Seep/w模块分析在不同的库水位下降情况下, 土坝非稳定渗流场的变化规律, 得到了坝体内不同时刻浸润线的变化趋势, 然后将计算得出的不同时刻坝体渗流孔隙水压力值应用于Slope/w模块, 分析了渗流影响下不同水位下降速度迎水面坝坡的稳定性, 计算了渗流场变化时的坝坡稳定安全系数。 1 非饱和土壤侵蚀分析原理 1.1 非饱和土体渗流理论 水库水位以不同的速度下降时, 由于坝体材料透水性的不同, 就会引起坝体内浸润线有不同的变化趋势, 孔隙水压力值也随之变化, 考虑《碾压式土石坝设计规范》(SL274-2001) 中库岸边坡内自由水面的曲线方程, 为了更加准确的分析坝体渗流及稳定性, 采用有限元数值分析法进行渗流场和稳定性分析。 非饱和土体渗流理论公式基于饱和与非饱和土体渗流的达西定律, 表述如下: 式中:q为单位体积的流量;k为渗透系数;i为总水头梯度。 非饱和渗流基本微分方程: 若假设总应力与孔隙气压力对时间的导数为零, 则可得: 式中:为与基质吸力变化有关的水的体积变化系数, 可将上式变为: 非饱和土体的渗流方程从形式上看与饱和土体完全类似, 但不同的是非饱和渗透系数不是常数, 而是随着点的位置的变化而变化, 为所在点的体积含水率或饱和度或者基质吸力的函数, 将上式进一步展开, 可得如下非线性方程: 1.2 边界水头h和边界边界面的计算 土坝非稳定渗流的定解条件包括初始条件和边界条件, 初始条件为: 根据渗流场情况边界条件可以分为三类: (1) 场变量 (水头H) 的边界条件: (2) 水流梯度 (流量Q) 的边界条件: (3) 混合边界条件: 式中:α为正常数, α和β都是边界条件中各点的已知量, 在计算的过程中运用迭代的方法依此去满足边界水头h和的已知关系。 1.3 算例2:总体渗流矩阵式,代次加乘子法 进行饱和-非饱和渗流的数值求解时, 将渗流场划分为小的单元体, 依此进行迭代计算, 此时考虑渗流计算的迭代方程为: 式中:D为总体渗流矩阵式, 是由De依次叠加组成, De矩阵中的各个元素可以表示为, E是由各单元矩阵Ee依次叠加组成, Ee中的各个元素可以表示为 2 计算有限值和渗透分析 利用Geo-slope有限元分析软件中的Seep/W模块, 模拟在不同的库水位下降的情况下坝体内浸润线的变化规律。 2.1 计算模型及网格剖分 新集水库枢纽工程位于甘肃省灵台县百里乡政府上游的达溪河上, 水库正常蓄水位35.0 m, 工程设计总库容为3.21×107m3, 回水长度约6.4 km, 库周长25.7 km, 库水面积3.2 km2, 最大坝高41.5m, 大坝坝顶高程为1 068.50 m, 坝顶宽8.0 m, 坝体上游边坡为1∶3和1∶3.5, 在高程为1 059.00 m处设宽度为2.0 m的马道, 下游坝坡比也为1∶3和1∶3.5, 在高程为1 049.00 m处设宽度为3.0 m的马道, 上游坡采用现浇C20F250混凝土板护坡, 下游坡采用C20混凝土网格草皮护坡, 设置贴坡排水, 坝型主要为均质土坝, 工程等别为Ⅲ等, 工程规模为中型, 根据Geo-slope提供的网格剖分功能, 将坝体划