基于DDA某水电站消力池高边坡稳定性分析.docVIP

基于DDA某水电站消力池高边坡稳定性分析 　　摘要：针对某水电站消力池高边坡，利用FLAC3D软件非线性有限差分法初步计算，得出坡体潜在滑动面的位置，建立了基于块体理论的非连续变形分析（DDA）块体系统分析模型。利用该模型分析了现状边坡和加固后边坡的位移与强度参数之间的关系，得出了边坡在加固前后的极限内摩擦角及相应的稳定安全系数。结果表明，该边坡在加固前安全系数较小，存在失稳的可能性，采用推荐方案预应力锚索加固后，边坡的稳定安全系数有显著提高，满足正常运行的要求。 　　关键词：岩体边坡； DDA；极限摩擦角；安全系数；稳定 　　中图分类号：TV223文献标识码：A文章编号2013 　　某水电站消力池高边坡后缘高程为3 240 m，前缘位于改线公路2 860 m高程附近，高差约380 m，坡角约42°，坡体中断层、裂隙及结构面发育，岩体破碎，为边坡变形破坏提供了潜在的??界条件。由于前缘处改线公路边坡开挖形成了高20～30 m的临空面，边坡高程2 860～3 010 m范围一带曾发生较大规模塌滑破坏。至今遇大雨或暴雨时，该边坡仍存在不同规模的塌滑，上部坡体的拉裂变形仍在持续。由于该边坡潜在变形体规模大，分布高程高，前缘塌滑后，临空条件变好，更有利于坡体变形的扩展，一旦遇连续暴雨可能再发生规模较大的破坏，影响电站消力池的正常运行。 　　本文采用基于块体理论的非连续变形分析（DDA）[14]方法，通过分析降雨控制工况下该边坡位移和强度参数之间的关系，得到边坡在加固前后的极限内摩察角及相应的稳定安全系数，为优化除险加固设计提供依据。研究中选择了有代表性的消力池高边坡22剖面进行计算分析，工程地质剖面见图1。 　　1岩坡块体系统模型 　　1.1DDA模型 　　根据该边坡实际情况及工程地质资料，选择典型剖面22进行计算分析。在分析该边坡岩体断层、裂隙及结构面产状等信息的基础上，依据所讨论问题的性质和要求、DDA法的建模原则，以及该边坡断层、裂隙的产状以及三组主要结构面的分布、连通情况，建立了块体系统的非连续变形分析概化离散模型。模型水平方向长5500 m，竖直方向高5400 m，左侧、右侧及底面均为固定边界，边坡表面为自由边 　　1.2岩体物理力学参数 　　边坡区出露的地层为上三叠统侏倭组（T3zh）变质岩和第四系（Q4）松散堆积物。边坡表部主要为崩坡积碎块石土，其下为风化较强的碎裂结构岩体，岩体拉裂、架空现象明显，完整性较差。除少量小断层及裂隙外，连坡区还发育有4组结构面：J1组，岩层层面，产状56°～74°∠53°～62°，反倾坡内，为主要结构面；J2组，产状160°～173°∠75°～82°，为顺坡向陡倾结构面；J3组，产状208°～239°∠25°～37°，为顺坡向缓倾结构面；J4组，产状273°～292°∠48°～58°，局部发育。 　　坡体具有明显的分层特征。结合现场调查、钻孔及洞探、物探揭露的坡体结构、变形程度，以及地表拉裂缝分布特征，将该坡体分为深层变形区Ⅰ和浅层变形区Ⅱ（图2）。Ⅰ区分布于Ⅱ区下部，位于边坡上梁后缘。Ⅱ区位于整个斜坡中下部，坡面发育多条延伸较长的拉张裂缝，变形破坏较Ⅰ区明显，主要表现为坡体中下部的浅表部塌滑和中后部的拉裂变形。 　　对于由裂隙、断层及结构面切割而成的非连续边坡来说，在无混凝土护坡及排水孔等措施的情况下，降雨入渗是影响其稳定的重要因素之一[56]。降雨使边坡表面以下一定深度岩体含水量升高，局部出现暂态饱和区，岩体系统的强度参数将大幅度降低，潜在滑动面抗滑力降低，岩体饱和引起自重的变化会导致滑动力增加，边坡的稳定安全系数降低[7]。 　　应用饱和-非饱和渗流理论建立边坡渗流分析模型，对降雨入渗作用下边坡渗流场进行模拟分析。可知，在降雨控制工况下，边坡降雨入渗的整体影响范围和边坡浅层变形体Ⅱ区形态较为接近。为简化DDA模型划分，保证一定的计算精度，将边坡浅层变形体Ⅱ区设定为降雨入渗暂态饱和区，该部分岩体采用饱和状态的强度参数。 　　根据提供的资料和参考工程经验拟定边坡岩土体、裂隙、断层及结构面的物理力学参数见表1。 　　2.1滑动面 　　在分析边坡稳定时DDA方法与其他数值计算方法有较大差别，有限单元法在建模时只需要对边坡岩层、断层构造等材料属性进行模拟，不必事先假设潜在滑动面，计算程序会自动根据节点受力等情况输出最后的单元变形及位移数据。DDA方法不能模拟单个块体在受力后的二次分割，只能模拟变形和位移，因此在边坡稳定分析前，首先需要确定抗滑稳定安全系数最小的潜在滑动面[8]。常规方法是根据边坡工程地质条件、节理面的产状及强度参数，参照已有滑坡体的滑动历史记录描述来确定边坡的最危险滑动面，但是这种方法受经验、勘测精度等因素影