均质土坝渗流场计算分析.docxVIP

均质土坝渗流场计算分析 人类在河流洪水的斗争中创造了文明，促进、消除劣势的历史过程形成了河流防洪工程体系的今天。由于渗流引起的大小险情难以胜数, 这表明堤坝渗流稳定是制约堤防防洪能力的主要因素, 由此引发的堤身稳定问题是堤防工程面临的重大课题。当前具有代表性的工程渗流问题主要表现在以下几个方面:①在水利工程中, 闸坝渗漏损失引起下游浸没, 土坝中渗流作用引起的冲蚀和滑坡等破坏。②在建筑工程中, 渗流引起的基坑边坡失稳, 基坑支护结构的破坏;渗透变形和破坏以及引起周围土体移动, 引发的地面建筑和地下设施破坏。③在采矿工程中, 矿坑突水造成矿井淹没, 矿山排水造成的塌陷等。④在农业工程中, 灌渠渗漏造成的工程效益下降等。⑤与渗流有关的环境问题, 地震引起的砂土液化, 水库蓄水诱发的地震, 大量抽取地下水或开采液体矿引起的地面沉降等。⑥库水位变化时, 引起上游边坡发生滑坡等问题。 本文以某均质土坝为研究目标, 建立堤坝渗流计算的数学模型, 设计多种库水位升降工况, 利用非饱和稳定渗流理论, 计算分析了其对渗流场的影响, 并在此基础上考虑了渗流场变化对上下游坝坡稳定性的影响。为水库调水、蓄水和泄水管理提供参考。 1 计算模型和分析方案的设计 1.1 大坝的模型计算 某均质土坝, 坝高50 m, 顶宽10 m, 上游坝坡为1∶2.5, 下游坝坡为1∶2, 均质坝填筑沙壤土料的渗透系数0.293 76 m/d, 孔隙率为0.464, 大坝的模型如图1所示。 1.2 水库泄水速度对大坝稳定性的影响 考虑到蓄洪时库水位上升速度不同对坝体渗流场和坝坡稳定性的影响不同;同时水库泄水速度的不同同样会对大坝的稳定性带来不利影响。为定量分析上述水库使用过程中可能出现的水位升降对坝体稳定性的影响, 为水库安全运行及科学管理制度的制定提供参考依据, 特设计了如下水位升降计算方案, 具体见表1。 2 结论分析 2.1 自由水域的变化 2.1.1 降速速度对边界条件的影响 计算边界条件:上游库水位从45 m分别以速度4 m/d、2 m/d和1 m/d下降到5 m水位, 下游水位保持在5 m的水位, 大坝底部按不透水边界处理, 初始渗流场为大坝上游水位保持在45 m时的饱和稳定渗流计算结果, 库水位下降工况下不同降速计算结果整理成图2, 时间单位为d。 从图2可以看出, 下降速度的不同对上游半边坝体渗流场影响较大, 图中表现为自由面上半段随着降速的减小下弯程度减弱;下游半边坝体渗流场则变化不大, 这跟筑坝材料的渗透系数有关, 渗透系数越大, 降速影响范围越大。而且20d之内自由面变化由密变疏, 说明降速较小时, 坝体内水分可以慢慢疏干, 降速大初期会形成超孔隙水压力, 坝体内水分来不及排出, 不利坝体稳定。不同水位下降速度, 自由面在100d以后变化逐渐变慢, 直到趋于稳定, 这是因为经过长时间的渗流, 上游水位早已降到某一确定的水位, 坝体内的孔隙水压力已逐渐消散。孔隙水压力消散的时间与坝体的渗透系数有很大的关系, 渗透系数小, 需要较长时间孔隙水压力才能逐渐消散, 自由面下降的速度会减慢很多, 自由面在初期变化很慢, 但水向上游“倒流”会严重得多。总的来说, 水位下降速度快, 在上游面附近会形成比较严重的“倒流”, 这样会严重危及上游坝坡的稳定。 2.1.2 不同水位上升速度计算结果 上游水位从0 m分别以速度2.0 m/d和1.0 m/d上升到45.0 m水位。下游水位以1.0 m/d速度上升到5.0 m的水位。库水水位上升工况下不同水位上升速度计算结果整理成图3。 从图3可以看出, 不同水位上升速度,60d之内的自由面靠近上游面的半边整体变化差别很大, 下游面的半边则变化不大。升速越大初期自由面的上弯越大, 这是水位上升坝体内渗流场来不及调整造成的, 筑坝材料的渗透系数越大, 升速影响范围越大, 但是影响程度越小, 即图中60d以内的自由面上弯程度减小。随着时间的延长, 坝体内渗流场调整到位,80d以后自由面开始稳定变化, 图中表现为80d以后自由面位置基本同幅演化。 2.2 坝坡稳定性分析 在上述非稳定渗流计算结果基础上, 也分析了各计算工况条件下的坝坡稳定性。整个计算过程在Geo-studio下完成, 其中渗流计算用seep/w模块实现, 然后把不同时段的渗流计算结果调入slope/w模块计算瞬态坝坡稳定性。 2.2.1 滑弧安全滑体及基层滑动的稳定性 由于库水位下降时上游面坝体内水分不能及时排出, 存在超孔压, 对坝坡稳定不利;下游面在该计算工况下安全系数只会越变越大, 因为下游坝坡水分补给来源消除, 残留水分还一直往坡外排出, 所以坡内水位降低, 坝坡逐步变为干坡。在此只分析最不利工况, 上游面坝坡稳定性计算结果整理成图4。 图4(a) 表明,