大渡河某水电站左岸边坡开挖的稳定性分析.docxVIP

? ? 大渡河某水电站左岸边坡开挖的稳定性分析 ? ? 钱洪建 黄梦婷 李 彪 (1.国家能源集团金沙江旭龙水电开发有限公司，四川 成都 610041；2.西南石油大学地球科学与技术学院，四川 成都 610500) 0 引言 在“碳达峰、碳中和”背景下，水利能源作为一种清洁能源越来越被重视，水电站建设也受到越来越多人的关注。目前我国西南地区有一大批规划和在建的水电站，但是西南地区高山峡谷众多，地势险峻复杂，地貌十分特殊，其高陡岩石边坡的开挖存在许多不容忽视的问题[1]。 近年来，关于岩质边坡稳定性研究的方法大致有以下几种：数值模拟法、工程类比法和极限平衡法。其中极限平衡法的原理比较简单，使用也较为广泛，但是其缺点也同样明显，对边坡开挖引起的应力变化和位移变化的分析不够准确，存在局限性[2]。工程类比法即综合分析与实际工程相似的工程案例，但是不可能有完全相同的两个工程，各个工程的地质情况也各不相同，从而可能导致分析产生一定的误差，因此并未被广泛应用于工程中[3]。而数值模拟法能够对边坡分步开挖引起的应力和位移变化进行模拟分析，优越性突出，被研究者广泛应用于各类边坡工程稳定性研究工作当中。 随着学科的发展和计算理论的不断完善，数值模拟技术的使用越来越广泛。数值模拟法大致可以分为以下几类：有限元方法、离散元方法、有限差分方法、边界元方法和多场耦合法等。国内外研究者也基于各种方法对边坡稳定性的研究取得了许多成果。Griffiths等[4]采用有限元（FE）方法对一些边坡案例的稳定性进行研究分析；Dawson等[5]基于强度折减法对边坡稳定安全系数进行了分析与对比研究；孙超伟等[6]基于霍克-布朗有限元强度折减法，为评价边坡稳定性建立了岩质边坡的稳定性图表，并且通过算例来佐证了图表的合理性；张辉等[7]依托于含有软弱夹层的岩质边坡，采用强度折减法对岩质边坡稳定性进行了研究分析。 本文基于上述研究，依托于大渡河某水电站，采用数值模拟与经验分析结合的方法对坝轴线左岸坝肩稳定性分析以及结构面对高陡岩石边坡破坏失稳机制及稳定性研究，以期为支护及开挖施工提供理论参考，也为类似结构面环绕开挖部分的岩质边坡的稳定性分析研究提供一些思路。 1 工程地质概况 1.1 工程地质条件 大渡河某水电站位于金川县城以北约13km、大渡河与右岸支流新扎沟河汇合口以上长约1km的河段上，河谷平面上呈反“S”型。河谷呈“V”型，枯水期河水位高程约2153m，水面宽约50～90m，正常蓄水位2253m高程相应谷宽约245～310m。坝址区两岸多基岩裸露，谷坡陡峻，为横向谷。大渡河某水电站边坡坝轴线左岸剖面及开挖区域和层间结构面如图1所示。 图1 边坡坝轴线左岸剖面图 1.2 岩土体物理力学性质 左岸边坡高程2030～2450m，坝基岩性为T32 2(5)-T32 2(6)的细砂岩，为顺向坡岩层结构，其中断层、错动带、裂隙有许多，但考虑到其对开挖的稳定性影响，因此在模拟中仅考虑结构面g15（PD3）以及断层f20和f28。 左岸边坡第一主应力大小为6MPa左右，方向为NNW向，其倾角在浅部与地表近于平行，在深部则逐步转为垂直；第二主应力大小为3MPa左右，方向为N-NW向，其倾角在浅部与地表近于垂直，在深部则逐步转为水平；第三主应力大小为1MPa左右，方向为N-E向，其倾角在浅部与地表近于垂直，在深部则逐步转为水平，整体应力水平较低。 2 有限差分法数值模拟 2.1 FLAC3D基本原理 FLAC3D属于有限差分软件，其在求解偏微分方程时用有限差分公式对导数进行替换，将复杂的偏微分方程转化为更易求解的代数方程[8]。 模拟过程中，首先需要将求解区域的网格细分为有限差分网格，网格建立完成后需选择与工程对应的本构关系并根据材料性质进行赋值，同时还需定义模型的边界条件及初始受力情况，定义完成后需进行第一次初始力学平衡。第一次平衡完成后若结果能够接受，需清零位移及速度场，再进行后续诸如开挖等静力计算分析。 2.2 模型的建立 根据水电站左岸边坡轴线工程地质剖面图在FLAC3D中建立左岸边坡三维开挖模型，如图2所示。模型横河向取230m，底部高程为2050m，最高点高程为2450m。模型x轴为逆河向，y轴为横河指向河谷，z轴垂直向上。模型底面采用x、y、z三个方向固定约束，前后左右四个面采用法向约束，坡面为自由表面未添加约束。边坡模型厚度为35m，水平距离为230m，垂直距离约230m，共生成网格1190个，单元4344个。并且在模拟中选用了摩尔-库伦弹塑性本构模型，在不考虑支护的情况下对边坡开挖卸荷过程进行数值模拟，并进行稳定性分析。 图2 边坡数值模型 模型开挖区域每隔20m进行划分，分为6个阶段从上至下开挖。同样不考虑支护情况下对边坡开挖卸荷过程进行数值模拟，并进行边