龙山水库大坝稳定性安全评价.docxVIP

龙山水库大坝稳定性安全评价 1 病险水库渗流 中国有8000多万座水库。目前，许多水库存在防洪模式低、技术破坏和维修不足的问题，导致了大量的保险水库。目前全国的病险水库约有30000万座,占水库总数的36%。这些病险水库严重威胁着下游人民生命财产的安全,不能充分发挥其兴利效益。因此,需要对病险水库进行安全鉴定,其中,对水库大坝进行渗流和稳定分析评价是非常重要的。本文对龙井水库大坝的渗流及坝体抗滑稳定进行了分析计算,以对其进行相应的安全评价。 2 坝体地质情况 龙井水库位于安徽省岳西县,枢纽工程等级为Ⅳ等,工程规模为小(1)型,永久建筑物级别为4级。水库设计防洪标准:50a一遇洪水设计,复核后设计洪水位564.09m(黄海高程系,下同),校核洪水位为564.71m,正常蓄水位为562.50m,最大下泄流量494.00m3·s-1。大坝为粘土心墙砂壳坝。 根据现场开挖探坑取样试验结果,坝体地质情况为:(1)坝壳料:以含细粒土砂为主,填土的干密度范围值1.37～1.55g·cm-3,现场渗透试验渗透系数平均值为8.2×10-3cm·s-1;(2)心墙料:以重粉质壤土、砂质粘土为主,软～可塑状,室内渗透试验渗透系数为1.31×10-5～6.94×10-7cm·s-1,填土的干密度范围值1.29～1.54g·cm-3;(3)碎裂岩,在坝体与坝基结合部位存在;下为花岗岩相对不透水层(图1)。 3 坝体稳定计算 根据Darcy定理及水流的连续性条件,当渗流系数主方向和坐标方向一致时,对于二维非均质各向异性土体渗流,考虑土的压缩性,非稳定渗流满足如下微分方程: ??x(kx?h?x)+??z(kz?h?z)=μs?h?t(1)??x(kx?h?x)+??z(kz?h?z)=μs?h?t(1) 式中,kx、kz、为x、z方向的渗透系数,h为水头函数,x、z为空间坐标,t是时间坐标,μs为贮水率。利用变分原理,可以将上述渗流问题转化为泛函在求解域内的极值问题,水库大坝非稳定渗流方程式的解应等于如下泛函所求的极值。 Ι(h)=∫∫Q{12kx(?h?x)2+12kz(?h?z)2+μs?h?th}dxdz+∫Γ2qhdΓ(2)q=μ?h*?tcosθ(3)I(h)=∫∫Q{12kx(?h?x)2+12kz(?h?z)2+μs?h?th}dxdz+∫Γ2qhdΓ(2)q=μ?h??tcosθ(3) 式中,I是h的函数,Q为求解域,Γ2为第二类边界条件,μ为给水度,h*表示自由面上各点的水头,θ为自由面向外的法线方向与竖直面的夹角,其余参数同上。经过求解得到支配方程: [Κ]{h}+[S]{?h?t}+[Ρ]{?h*?t}={f}(4)[K]{h}+[S]{?h?t}+[P]{?h??t}={f}(4) 式中,[K]、[S]、[P]为系数矩阵,{f}为已知常数列向量,其余参数同上。本文采用三节点三角形等参单元进行求解。 坝体稳定计算采用极限平衡法进行,极限平衡法是建立在摩尔-库仑强度准则基础上的,只考虑静力平衡条件和土体的摩尔-库仑破坏准则。极限平衡法采用条分法进行,这种方法假定边坡稳定问题是平面应变问题,以极限平衡理论为理论基础,设边坡的稳定安全系数为F,则当边坡体材料的抗剪强度(凝聚力c和摩擦因数tanφ)降低F倍后,边坡内某一最危险滑面上的滑体将濒于失稳的极限平衡状态。本文采用瑞典圆弧滑动法和简化Bishop法进行坝体稳定计算。 瑞典圆弧滑动法稳定性系数计算公式: Fs=∑[c′ili+(Wi-uibi)cosaitanφ′i]∑Wisinai(5)Fs=∑[c′ili+(Wi?uibi)cosaitanφ′i]∑Wisinai(5) 式中,Fs为安全系数,c′i、φ′i分别为第i土条底部滑动面处土的有效抗剪强度指标,bi、li分别为第i土条的宽度和底部滑弧长度,Wi为土条重量,浸润线上为天然容重,浸润线下为饱和容重,ai为第i土条底部圆弧切线与水平面夹角,ui为第i土条底部的孔隙水压力。 简化Bishop法稳定性系数计算公式: Fs=∑1mai[c′ibi+(Wi-uibi)tanφ′i]∑Wisinai+∑QieiR(6)mai=cosai+(tanφ′isinai)/Fs(7) 式中,Qi为i土条水平力;ei为第i土条的水平力至滑弧中心的力臂,其余参数同式(5)。 4 计算与分析 4.1 ．计算结果为 渗流计算的坝身各土层的计算指标见表1。 本文就以下工况进行计算:(1)上游正常蓄水位562.50m与下游相应的最低水位;(2)上游设计洪水位564.09m与下游相应的水位;(3)上游校核洪水位564.71m与下游相应的水位;(4)库水位由校核洪水位564.71m骤降至最不利水位542.78m。各种工况下的下游水位,由于无法知道其准确值,计算