库水位对心墙坝渗流稳定影响研究

2020-02-13陈章淼

水利技术监督 2020年1期

陈章淼，代彬

(贵州省水利水电勘测设计研究院，贵州贵阳 550002)

1 概述

坝体安全关系到水利工程的正常使用以及下游居民生产生活的安全，准确评价坝体安全是十分必要的[1]。影响坝体安全的因素众多，水是其中一个重要的因素，水可以降低坝体材料强度，同时水压力也是诱发坝体破坏的一个主要因素[2]。同时坝体自身结构设计也是影响其最终稳定性的一个重要因素[3]。目前，针对坝体渗流、稳定性研究的主要方法包括：数值模拟[4- 10]、理论计算[11- 13]；张寅寅等[14]使用有限元数值模拟分析不同排水条件下土坝渗流及稳定性情况；魏玉师[15]基于BiShop理论通过FORTRAN语言自编渗流分析软件，研究非饱和尺度土坡渗流及稳定性情况。

2 心墙坝计算模型

2.1 坝体结构

万营水库位于万营河距水城县城约80km，集水面积7.30km2。流域处于梯级状斜坡地带，万营水库工程为小(1)型水库。

水库大坝为粘土心墙防渗风化料坝壳坝，大坝坝顶设防浪墙，墙顶高程为1579.100m，高出坝顶1.0m，坝顶高程为1578.100m、长95.135m、宽6.0m，最大坝高41.1m。坝轴线方位角NE20°10′5″。

坝体材料分区、坝体填筑设计标准及大坝标准剖面如图1所示，填筑料见表1。

2.2 计算参数

大坝的采用粘土直心墙防渗，防渗料选用坝址附近的粘土料，上游及下游填筑区采用料场风化料、开挖洞渣料及开挖全风化料。材料参数见表2。

3 数值模拟模型及计算工况设计

3.1 数值模拟模型

根据心墙坝特征建立数值模拟模型(Geo-Studio)，模型地基自坝基至不透水层，建立水头变动时间参数模拟水头变化过程，数值模拟计算模型如图2所示。

3.2 数值模拟计算工况

为了分析不同水位变动工况(正常蓄水位到死水位，30年一遇洪水位到正常蓄水位，死水位到正常蓄水位，正常蓄水位到30年一遇洪水位)下心墙坝的响应特征，确定库水位变动速率为0.2、0.5、1.0m/d，通过渗流量、浸润线以及上、下游坝坡安全系数反映心墙坝的响应特征。水位变动坝体响应数值模拟计算工况见表3。

表3 计算工况

表1 坝体填筑标准

表2 数值模拟参数

图2 数值模拟计算模型

同时为了研究不同防渗墙深度对坝体防渗的影响，确定防渗墙深度分别为20、30、40m，研究库水位分别在1556、1565、1575、1576.93情况下的渗流情况。

3.3 数值模拟计算结果

3.3.1浸润线分布特征

不同库水位坝体浸润线如图3所示。从图3可知，在库水位恒定的情况下，不同库水位浸润线变化曲线分布规律较为一致，在心墙部位浸润线出现明显的下降，这是由于心墙渗透系数较低具有较好的防渗阻水效果。同时可以看出下游坝坡范围内浸润线分布较为密集，这表明上游坝坡浸润线高度与下游坝坡相比对库水位的变动响应更为敏感。

图3 20m防渗墙不同库水位浸润线

不同库水位上升速率坝体浸润线变化情况如图4所示。从图中可知，水位变化时浸润线变化规律与静态水位变化规律明显不同，浸润线突降部位位于心墙上游处，高水位与低水位浸润线在心墙下游边界处汇集，水位上升过程中，浸润线分布较为密集。

图4 不同水位上升速率浸润线分布特征

3.3.2防渗墙深度对渗流的影响

不同防渗墙深度渗流量计算结果如图5所示。从图中可知，心墙坝渗流量与库水位高程呈正相关关系，与防渗墙深度呈负相关关系；渗流量随着库水位的增高以及防渗墙深度的减小而增大。这表明随着防渗墙深度的增加，其防渗堵水效果越好。

图5 不同深度防渗墙渗流量

3.4 坝体稳定性计算结果

3.4.1不同库水位坝体稳定性分析

坝体稳定性是坝体安全性的直接反应，具有重要的研究意义和实际价值。上游坝坡、下游岸坡安全系数与库水位关系如图6所示。从图中可知库水位越高上游坝坡安全系数越低、下游坝坡安全系数越高，3种防渗墙深度，上游坝坡均处于安全情况，但下游坝坡在防渗墙深度为40m时处于安全情况。

图6 不同库水位上、下游坝坡安全系数

3.4.2水位变动坝体稳定性分析

图7 坝坡安全系数与时间关系

水位变动过程中，坝坡安全系数变化过程如图7所示。从图中可知，在初始状态下，坝坡安全系数均处于相对较高的水平，在水位变动过程中，坝坡安全系数处于一个动态变化的过程，随着库水位的上升，坝体安全系数逐渐减小，且水位变动速率越大，初期安全系数变化幅度越大，但水位达到最高值并稳定一段时间后，各个工况的安全系数趋于一致；之后随着水位降低，坝坡安全系数逐渐增大，且水位下降越慢到达同一水位高程时，坝坡的安全系数越高，这主要是因为，在库水位快速下降过程中，坝坡内部的水并未完全排出，坝坡材料处于饱水状态，强度较低，且库水快速下降降低了水压力对坝坡稳定产生的正向作用。