-
工程地质层组的划分是地质模型概化和岩土体质量评价的基础和首要环节。通过工程地质层组的划分，可以达到如下两个目的：(1)地质模型概化和参数研究，有助于设计师找出主要矛盾，更好地评价主要工程地质问题；(2)可以更为清晰的反应地质结构的变化规律[1]。数值仿真方法[2-6]可以采用施工前的地质调查资料、设计资料和现场监测资料进行模拟，了解工程设计的薄弱环节，达到改造设计的目的。欲使岩体内部响应(变形、位移和应变)逼近真解，其数值必须遵循几何、本构、力学和过程的仿真，这是变形体解答的唯一性定理所决定的[7]。数值仿真方法是一种近似算法，影响计算结果的准确性的因素较多，包括地质边界条件的准确性、地质概化模型的合理性、围岩力学参数的准确性、计算方法和本构模型的合理性、开挖和加固过程的仿真性等等，工程地质概化是进行洞室施工数值模拟的基础，因此，本文基于折减的围岩力学参数、设计支护参数和设计开挖方案，以某水电站大型调压井工程为研究背景，首先进行工程地质概化分析，然后建立地质概化模型和三维数值仿真模型进行数值模拟计算分析，对数值仿真的围岩应力、位移等力学性质进行分析，了解工程施工的围岩稳定性状况，获得工程区数值仿真合理的本构模型和初始地应力场，为洞室施工其它相关问题研究提供借鉴和指导。
HTML
-
某水电站尾水建筑物结构为三条尾水支洞连接一条尾水隧道，其交叉处连接一个尾水调压井共布置了三个尾水调压井。三个调压井按“一”字形布置，间距为102.0 m。计算考虑洞室有三个调压井、连通上室、九条尾水支洞、三条尾水隧洞。调压井开挖直径为29.3～34.3 m，净空高度92.0 m。调压井下部为五洞交叉口悬空结构，结构跨度大，对其洞室的支护提出了更高的要求。
-
当前岩体变形参数的取值常按风化划界或以某一大的软弱层带或岩性相差很大的界面为界，将不同区带岩体视为等值连续体而取一个固定不变的综合变形模量进行计算，这是目前计算中地质原型概化和参数取值常用的模式[8]，某水电站调压井工程区以表层均匀风化为主，由地表向深处风化程度逐渐减弱。将工程区风化程度划分为全风化、强风化、弱风化上部、弱风化下部和微风化～新鲜五级。全风化底界水平与垂直深度一般0.00～10.00 m，强风化底界约10.00～20.00 m，根据揭露的风化程度信息，全风化和强风化分界线不明显，因此将这两层风化岩体归类为一层—全强风化岩体，便于分析。弱风化上部(弱上岩体)底界水平深度一般10.00～80.00 m，垂直深度一般10.00～50.00 m；弱风化下部(弱下岩体)底界水平深度一般30.00～90.00 m，垂直深度一般30.00～50.00 m；弱下岩体以下为微风化岩体和新鲜岩体，简称为微新岩体。
根据岩性和风化程度，将调压井工程及围岩划分为四种介质。全风化和强风化岩层地质界线不明显，均呈疏松状态，为散体结构，作为一种介质，即全强风化岩体；弱风化带上部比较疏松，含少量完整结构，作为一种介质，即弱上岩体；弱风化带下部比较坚硬，含少量散体结构，作为一种介质，即弱下岩体；微风化和新鲜岩体岩体性质差别不大，仅沿裂隙有轻微的风化，作为一种介质，即微新岩体。调压井工程洞室布置于微新岩体中，围岩岩性较好。
-
将工程施工前就存在于岩体中的地应力，称之为初始应力或天然应力[9]。根据该工程区ZK448和ZK450钻孔地应力监测资料发现，该区域垂直向应力以自重为主，地应力转化至三维有限元计算模型坐标下的应力与埋深呈线性关系，如图1所示。
Figure 1. Stress and depth relationship of 3-D finite element simulation model
根据不同的计算工况来施加初始地应力场，分为只考虑自重应力场、自重应力和构造应力联合应力场。考虑自重应力和构造应力联合应力场时，垂直方向采用自重应力，构造应力根据实测应力回归线性方程模拟，即x向构造应力与埋深h的关系为：σx＝0.0377h＋4.1556，y向构造应力与埋深h的关系为：σy＝0.0305h＋3.4935。
-
通过尾水调压井地带的Ⅲ级断层有F20、F21、F22，属于Ⅲ级结构面的一般性断层、其余小断层、挤压面等Ⅳ级结构面发育平均间距约为20 m。调压井工程区主要断层详细情况见表1。
断层 产状 宽度/m 迹长/m 等级 F20 N20°～30°E，NW∠60°～64° 0.3～0.4 >200 Ⅲ F21 N10°～28°W，SW∠72°～82° 0.9 >300 Ⅲ F22 N0°～22°W，SW∠42°～58° 0.5～1.2 >660 Ⅲ Table 1.
Main faults of the surge tank engineering area details 经地质资料分析，该调压井工程区受F20、F21、F22断层影响较大，数值仿真时地质构造考虑这三条断层，断层厚度、岩性和方位根据实际情况进行概化。其他结构面因较细或距离开挖洞室较远，对洞室围岩稳定性影响较小，计算时不考虑。
-
许多岩土工程事故的分析表明，裂隙岩体渗流场和应力场耦合作用对工程安全的重要性不可忽视[10]。根据本工程的现场调查情况统计，微新岩体干燥或只有少量渗水，渗水量(钻孔压水试验求得岩体渗透系数，再由多布诺沃里斯基公式求得围岩的渗水量)为0.66～1.00 L·min-1·m-1；弱下岩体有渗水或滴水，渗水量为0.500～ 0.66 L·min-1·m-1；弱上岩体普遍有渗水滴水，渗水量为0.40～0.50 L·min-1·m-1；全强风化岩体普遍有渗水滴水，渗水量为0.40～0.50 L·min-1·m-1。
由于调压井洞室所处的微新岩体地下渗水较少，洞室施工亦未见围岩有明显地下水渗出，因此在数值计算时不考虑地下水对洞室施工的影响。
1.1 工程概况
1.2 围岩岩性概化
1.3 初始地应力场概化
1.4 结构面概化
1.5 地下水概化
-
根据上述地质概化分析，建立工程区的三维有限元计算模型。模型范围：x＝-349～-103 m，y＝0～450 m，z＝501～865 m，共划分单元184 336个，节点38 954个，采用FLAC3D软件进行模拟采用位移边界，实体单元模拟，钢筋和锚杆采用cable单元模拟。地下洞室布置及断层有限元计算模型见图2，围岩包括全强风化岩体、弱上岩体、弱下岩体、微新岩体，调压井地下洞室群处于微新岩体中，为Ⅱ级围岩，岩性较好。工程区三维地质模型如图3。洞室总共分38步开挖。初期支护参数的锚杆和锚索采用cable单元模拟，喷混凝土采用shell单元模拟，设计初期支护参数见表2，初始地应力采用侧压力系数法施加。
Figure 2. Underground cavern layout and faults finite element computation model
Figure 3. 3D Geological model of surge shaft project area
位置 部位 支护参数 锚杆 预应力锚索 喷混凝土 井身、顶拱及五岔口 顶拱 Φ28@2 m×2 m，L＝6 m 1 000 kN无粘结 C30，厚度0.2 m Φ36@2 m×2 m，L＝9 m 井身及五岔口边墙 Φ25@2 m×2 m，L＝4.5 m 1 000 kN无粘结 C20，厚度0.15 m Φ32@2 m×2 m，L＝9 m 锁口锚杆 预应力：125 KN — — Φ32@1 m×1 m，L＝9 m Φ32@2 m×2 m，L＝9 m 连通上室、尾水支洞、尾水隧洞 顶拱及边墙 Φ25@2 m×2 m，L＝4.5 m — C20，厚度0.15 m 锁口 预应力：125 kN — — Φ32@1.5 m×1 m，L＝9 m Table 2.
Design primary support parameters -
由围岩变形监测资料分析可知，施工过程中调压井围岩未出现流变特征，所以只考虑弹性本构和弹塑性本构模型对其进行模拟分析。
-
弹性本构模型为各向同性的线弹性本构，其主要计算参数为弹性模量E，泊松比μ和岩体密度。
-
弹塑性本构采用Mohr-Coulomb本构模型，所需要的参数包括弹性模量E、泊松比μ、岩体密度、抗拉强度、粘聚力c和内摩擦角φ。计算力学参数见表3。
围岩类别 微新岩体 弱下岩体 弱上岩体 全强风化岩体 断层 摩擦角φ/(°) 54.46 50.19 41.99 34.99 11.97 粘聚力c/MPa 1.70 1.00 0.30 0.20 0.03 弹性模量E/GPa 24.05 10.00 3.20 0.23 0.60 容重/(kN/m-3) 26.10 25.70 24.50 23.50 20.00 泊松比 0.22 0.24 0.26 0.28 0.33 抗拉强度/MPa 8.00 6.50 3.40 1.50 0.28 Table 3.
Mechanical parameters for finite element calculation
2.1 三维有限元模型建立
2.2 本构模型及其参数
2.2.1 弹性本构模型及其参数
2.2.2 弹塑性本构模型及其参数
-
根据分析目的，拟定四种计算工况进行数值模拟分析，计算工况见表4。
工况编号 本构模型 工况说明 GK1 线弹性 初始地应力只考虑自重应力场 GK2 线弹性 初始应力场考虑自重和构造应力场 GK3 弹塑性 初始地应力只考虑自重应力场 GK4 弹塑性 初始应力场考虑自重和构造应力场 Table 4.
Calculation conditions of numerical simulation
-
各个工况开挖完成后，围岩应力、位移计算结果分别见表5。为确定工程区洞室施工需采用的合理本构模型和初始地应力场，选取2号、5号、8号尾水支洞三个收敛断面(尾支0＋098.00、尾支0＋102.00、尾支0＋090.00)的表面收敛计BC测线位移进行对比分析，BC测线高程为573.25 m，表面收敛计BC测线布置如图4。各个工况下测线收敛位移与计算位移对比见图5～图7。
Figure 4. Measuring line BC of surface convergence gauge layout
Figure 5. Comparison of measuring line BC displacement of section 0＋098.00
Figure 6. Comparison of measuring line BC displacement of section 0＋102.00
Figure 7. Comparison of measuring line BC displacement of section 0＋090.00
工况编号 主应力范围/MPa 总位移/mm 应力位移计算结果分析 大主应力 小主应力 GK1 -25.13～-0.05 -2.99～1.26 0.00～8.74 (1)弹塑性模型的拉应力值较线弹性明显要小；初始地应力场对围岩的应力计算结果影响较大；相同本构模型下，构造应力场使得围岩拉应力增大；不同的本构模型下，围岩最大拉应力均较小，最大拉应力为2.70 MPa，远小于微新岩体的抗拉强度，围岩不会因为应力过大而发生破坏；调压井洞室拐角、断层穿越部位容易出现应力集中；(2)初始地应力场和本构模型对围岩变形影响均较大；对于相同的本构模型，施加构造应力场使得变形增大；弹性本构模型开挖造成的位移要比弹塑性小；在不同的本构模型下，围岩变形较小，最大位移17.20 mm，围岩较稳定；构造应力场对围岩位移分布影响很大；断层对位移有影响，断层附近位移相对较大。 GK2 -42.41～-2.46 -8.44～2.70 0.00～15.73 GK3 -26.52～-0.04 -4.45～0.83 0.00～13.21 GK4 -32.54～-0.08 -9.75～1.02 0.00～17.20 Table 5.
Statistical results of surrounding rock stress and displacement after calculation of each conditions 由测线位移对比可知，尾支0＋098.00断面处GK1、GK2、GK3、GK4计算值相对监测值变化值分别为15.18%、5.45%、13.14%、3.80%；尾支0＋102.00断面处GK1、GK2、GK3、GK4计算值相对监测值变化值分别为15.96%、5.25%、8.99%、2.23%；尾支0＋090.00断面处GK1、GK2、GK3、GK4计算值相对监测值变化值分别为21.30%、6.59%、9.86%、3.18%。
因此，GK4计算位移变化过程与监测的位移变化的相关性最好，该水电站大型调压井选择弹塑性本构模型并施加水平构造应力场进行施工期数值模拟计算，更加接近真实的工程区围岩特性。
-
结论如下：
1)洞室施工过程中的最大拉应力均小于调压井围岩抗拉强度，调压井围岩不会因为应力而发生破坏。
2)初始地应力场对围岩的应力和变形影响较大；调压井洞室拐角、断层穿越部位容易出现应力集中，结构面对围岩变形有较大影响；模型本构类型对应力位移均有较大影响。
3)围岩最大变形17.20 mm，说明该水电站大型调压井工程围岩较稳定，洞室开挖不会破坏围岩的稳定性。围岩与断层相交处以及断层附近围岩较松软，所以造成变形相对较大，需要加强监测和支护。
4)通过计算位移和监测位移对比分析，弹塑性本构模型并施加构造应力场更加接近真实的工程区围岩特性，对其它问题研究的数值计算具有重要的参考价值
DownLoad:
