后河水库坝肩稳定分析.docx

1、后河水库坝肩稳定分析后河水库坝肩稳定分析 - 农田水利1 主要工程问题及研究对策狮子滩水电站位于四川省长寿县境内，是龙溪河四级开发利用中最上游一级。电站装机容量48万kW，年发电量206亿kW&midd0t;h，水库总库容83亿m3，可对天然来流量作多年调节。主要水工建筑物由混凝土挡墙堆石坝、溢洪道、有压引水系统及电站厂房组成。溢洪道位于堆石坝左端条形山脊天然垭口，总宽1120 m，设计最大泄洪能力49l00 m3s，从溢流坝至挑水坎总长3580m，其中，首段为重力式混凝土溢流坝，共分5个坝段，坝段长220-230m， 闸墩宽40m，在闸墩中心线分沉陷缝，建基面高程从1号至5号坝段为32653

2、180m，溢流堰顶高程3400m，坝顶高程3500m，最大坝高32，0m，溢流坝上设置70 m180m(高宽)弧形闸门5孔。中段为泄水陡槽，长3225m，平均坡降495，末端为混凝土低鼻坎矩形差动式挑水坎。狮子滩水电站于l957年3月全部建成，至今已安全运行41年。在多年运行期间，水库最低蓄水位32509m，较死水位低314m，最高运行水位34758m，较正常高水位高085m，溢洪道实际下泄最大流量1750 m3s，仅达设计泄洪能力的356，溢洪道各建筑物还未经受设计最大泄流量考验，同时，经过多年蓄水运行，溢洪道相继出现一些影响安全运行的工程问题，主要包括：(1) 溢流坝上游铺盖被淘空。原设计

3、溢流坝上游铺盖为1：4坡度的粘土 铺盖，表面再铺块石保护，施工时作了修改，铺盖粘土只回填至3300m高程，3300m高程以上改用废渣回填，经过多次泄洪，下泄水流将回填废渣及3300m高程以下的粘土大量带走，铺盖大部分已失去保护，溢流坝前沿基础已经暴露在外，致使渗径缩短，增加了帐幕前坝基渗透压力，并可能拉裂帐幕，影响溢流坝体的整体安全。(2) 溢流坝稳定安全参数偏低。由于溢流坝基础大部分位于薄层砂岩上，砂岩厚度一般为3040m，以下为砂质粘土与砂岩互层，岩性变化大，节理裂隙发育，问有风化破碎带，在砂质粘土岩中分布有埋深较浅，产状近于水平且连通性较好的G4、G5软弱夹层，其中建基面以下80m的G5

4、夹层分布于整个溢流坝基。原设计根据竣工资料，对1号一5号坝段在校核洪水条件下的稳定性进行复核，1号一3号坝段抗滑稳定安全参数Kc为10108，4号坝段为091，5号坝段为094，由此可见4号、5号坝段的Kc值不能满足现行规范要求。为客观评价4号、5号坝段的稳定安全性，设计上建议在稳定分析中考虑尾部岩体的抗剪作用，这种计入尾部岩体抗力的分析方法是合理的，但对4号、5号坝段稳定安全参数量值上的变化设计上未作详细说明。(3) 溢流坝普遍出现贯穿性裂缝。据狮子滩溢洪道运行总结报告介绍，溢流坝体出现的贯穿性裂缝可分为三类：第一类是遗留的施工冷缝；第二类是由于混凝土浇筑过程中未采取严格温控措施和拆模过早等

5、因素形成的，先期(末蓄水发电前)裂缝；第三类则是在多年运行期间，由于坝面温度随外界气温变化所产生并逐步演变形成的。上述三类贯穿性裂缝的形成原因，以及对坝体安全性的影响程度，是工程上所关心的问题。针对狮子滩电站溢洪道所存在的上述三个主要 问题，结合溢洪道多年的运行与安全监测资料，以及 1990年狮子滩水电站大坝安全定检时所修定的最新地质水文成果，运用三维弹塑性有限方法，以上游铺盖、溢流坝及下游护坦和坝基岩体作为计算整体，对溢流坝体的变位与应力特征，以及上述(1)、(2)因素对坝体稳定性的影响进行定量分析评价，给出狮子滩电站溢流坝现今与进行工程处理后所具有的整体安全储备，同时，还对溢流坝出现贯穿性

6、裂缝的原因及其影响作了初步分析。2 有限元数学模型及计算方法21 有限元计算范围及结构离散选取最大坝高，抗滑稳定安全参数最低的4号坝段作为分析对象，并沿上下游及坝基方向延伸一定范围构成三维计算域，其中：垂直水流方向(Y向)：溢流坝1+813l溢流坝l+10331m；平行水流方向(X向)：0-3000-0+8000m：垂直方向(Z向)：建基面以下320m至坝顶，绝对高程290.0 3500m，相对高差600m。该区域包括4号坝段左右闸墩，溢流坝体、上游前坦300m与下游护坦500m，以及坝基以下薄层砂岩、砂质粘土岩互层与软弱破碎带G5。根据溢流坝受力特点及坝基地层分布界面，对计算域进行剖分，其中

7、坝体和护坦混凝土结构，以及坝基岩体采用空间8节点等参元模拟，G5软弱夹层以8节点夹层单元模拟，为反映坝体与坝基岩体接触面相对剪切变形，该接触面用接触界面单元模拟。22 材料物理力学参数参考原设计采用材料物理力学参数，并结合1990年狮子滩水电站溢洪道坝基岩石力学试验成果，本次有限元分析计算选用材料参数见表l。23 荷载及其组合本文有限元分析考虑两种荷载组合：基本荷载组合：岩体自重+坝体自重+正常水载=泥沙压力+扬压力特殊荷载组合：岩体自重+坝体自重+校核水载+泥沙压力+扬压力其中水载及泥沙压力基本参数为：工作桥每孔自重1010 t，交通桥每孔自重1421t，闸门自重505t，均以均布力作用于左

8、右闸墩上，工作桥与交通桥面活荷重作为安全储备，略去不计。岩体自重在地质历史时期已形成，是一种环境力场，坝体自重在坝体完建后业已形成，是一种初始力场，有限元分析中只计入两者形成的自重应力。考虑到溢流坝上游铺盖已淘至坝基混凝土齿墙，坝基灌浆帐幕前渗透压力按10倍H计，帷幕线及排水孔线位置渗透压力水头折减系数按设计规范取为05与03，下游坝趾处渗透水头参考扬压力实测资料取为025倍H。24 岩体强度及本构模型在外荷载组合作用下，溢流坝及坝基岩体应力集中或材料强度较低的部位首先出现开裂或塑性破坏，由此将引起这些部位承载能力下降，超过其承载能力的部分荷载(超余应力)将转移至附近坝体或岩体单元，进而可能引

9、起附近单元破坏。因此，溢流坝及坝基岩体的变位、应力及破坏是不断调整最终趋于稳定(收敛)的渐进过程，非线性弹塑性有限元法能较好地模拟分析这一客观过程。溢流坝混凝土及坝基岩体材料开裂条件用宏观强度描述：iiRt(i1，2，3) (21)式中 ii吨代表应力特征点的3个主应力，以拉为正；Rt表示混凝土或岩体材料抗拉强度，分析中可能呈现单向，双向及三向开裂情况，由程序自行校核并进行刚度修正。坝体混凝土或岩体是否进入塑性状态，按DruckerParger准则判别：FI1+J2-k (22)式中 I1和J2分别为应力张量的第一不变量和应力偏张量的第二不变量；、k是与材料摩擦系数tg和凝聚力C有关的常数。据

10、关联流动法则导出材料的弹塑性本构阵：坝体建基面及坝基G5软弱错动带，按层面法向不抗拉材料分析，其剪切滑移按mohrCoulomb条件复核。25 坝基抗滑稳定安全系数计算方法采用整体强度储备安全系数表征坝基的抗滑稳定性，基本方法是假定基本或特殊组合荷载不变，通过不断下浮坝基与岩体接触面，G5软弱夹层及尾部抗力体(砂岩)抗剪强度指标，分析强度参数下浮倍数与潜存滑移通道中破坏面积之间的演变关系，寻求坝基滑移面整体移滑时，相应的Kc值即为整体抗滑稳定安全系数。3 坝基抗滑稳定安全度评价针对溢流坝上游铺盖被淘空条件，首先研究正常蓄水与校核洪水两种工况，坝体沿建基面与G5软弱破碎带两种滑动面的抗滑稳定性(

11、文中G5滑动面包含坝基下近于水平分布的G5破碎带及坝趾尾部砂岩抗剪断段)，其成果见表2中方案1-方案8，在此基础上，分析坝体前坦深槽回填处理后，以及坝基扬压力两种分布形式相应的坝体抗滑稳定性，其成果见表2方案9方案l0。表2中沿坝基面或G5软弱破碎带滑移通道的破裂面积比是指滑移通道中产生破坏的面积与总面积之比，该参数可以定量评价坝体抗滑稳定安全系数的可靠程度(保证率)。分析表2可以看出： 在不计滑移通道尾部岩体抗剪作用及正常蓄水位条件下，坝体抗滑安全系数12时，沿G5滑移通道的破裂面积比达948，已接近整体滑移极限(破裂面积比l00)。因此，在正常蓄水位及不计尾部抗剪作用条件下，坝体的最大抗滑

12、稳定安全系数为120，考虑一定安全储备，若以沿滑移通道破裂面积比90控制(下同)，坝基实际具有抗滑稳定安全系数110-115。由方案5和方案6可看出，在校核洪水下，若不计入尾部岩体抗剪作用，抗滑安全系数分别为100 和110时，坝体沿G5滑移通道的破裂面积比均已 超过90，表明坝体实际具有安全系数小于10。 方案3和方案4表明，计入滑移通道尾部岩体抗剪作用，正常蓄水时抗滑安全系数120时，沿G5滑移通道的破裂面积比仅867，表明坝体的整体抗滑安全系数已达到120。方案7和方案8表明，计入滑移通道尾部抗剪作用，校核洪水当抗治安全系数110时，坝体沿G5滑移通道的破裂面积比超过90，而安全系数为1

13、0时，破裂面积比仅达675，因此，校核洪水时计入尾部抗剪作用的坝基抗滑安全系数为105左右。 考虑坝体前坦深槽回填处理，方案9给出了不计水平防渗处理效果(坝基扬压力较回填前不变)，坝基抗滑安全系数为110时，沿G5滑移通道的破裂面积比较方案8下降了近10(即安全度保证率提高l0)，方案10表明计入水平防渗效果(坝在扬压力减小)，沿G5滑移通道的破裂面积比较方案9又进一步降低，即安全度进一步提高。综合方案9、方案10来看，对坝体前坦深槽采用回填并作粘土水平铺盖处理，校核洪水时计入尾部抗剪作用的坝基抗滑安全系数可达110。 从方案l至方案10中可看出，正常蓄水位与校核洪水两种工况，以及是否计入尾部

14、抗剪作用多种组合条件下，沿G5滑移通道的破裂面积比均大于坝体建基面，表明埋深较浅产状近于水平的G5软弱破碎带是坝基抗滑稳定的关键。4 坝体裂缝形成原因浅析前面已提及狮子滩水库溢流坝存在的贯穿性裂缝可以分为三类，第一类是坝体施工过程中遗留的冷缝，第二类是由于混凝土浇筑过程未采取严格温控措施等因素形成的先期裂缝；第三类是坝面温度随外界气温变化演变形成的裂缝。对于第一、二类裂缝可采用灌环氧胶等工程措施一次性处理，而第三类裂缝在水库长期运行过程中，可能重复出现，尤应引起重视。对第三类裂缝产生的条件，分别用三维有限元法和重力坝设计规范建议方法进行计算分析。41 坝体运行期温度应力三维有限元计算结果狮子滩

15、电站溢流坝已运行多年，坝体内部混凝土稳定温度接近建坝区多年平均气温，而上下游坝面随外界气温呈年际变化，坝面与坝体内部的温差将引起坝体变形，这种坝体变形由于受到坝基岩体的约束，在坝体内将形成一定的温度应力。参考坝区气温、水温和地温资料，拟定以下参数计算夏季和冬季两种情况：(1) 坝体内部稳定温度18，下游坝面气温37，上游水温在3400347.0m高程内取25，3300-3400m高程内取20，地温取l0。(12) 坝体内部混凝土稳定温度18，下游坝面气温-2，上游水温在3400-3470 m取2，3300-3400m高程取6。经计算在溢流坝段(35剖面)堰顶及下游坝面出现主拉应力最大，其方向近于平行坝轴线，量值达08-12MPa。42 按重力坝设计规范分析(1) 外界温度年变化引起的应力。最大应力一般出现在坝表面，值按下式估算：若假定坝面混凝土允许抗拉(裂)强度t10MPa，由上式计算，坝面允许的温度变幅为15 ，换言之，若坝面温度