1. 引言
在地下开挖过程中，由于岩体被结构面切割成大小不同、形状不一的各种岩块，在自然状态下，这些空间块体处于静力平衡状态 [1]。当进行隧道开挖时，使暴露在临空面上某些块体失去原始的静力平衡状态，因而造成某些块体首先沿着结构面滑移、失稳，进而产生连锁反应，造成整个隧道围岩的破坏 [2]。在岩体结构的研究中，各类结构面的分布情况、表面特征、连续性、密度、结构面之间的接触关系及其的空间组合形式为主要的研究内容，其中结构面之间的空间组合形式为最重要的研究内容。
由R. E. Goodman与G. H. Shif提出的块体理论 [3]，是近年来发展和完善起来的一种适合于工程岩体稳定性分析的有效方法。该理论提出后受到国内外学者的广泛重视，刘锦华和吕祖珩最早将块体理论介绍到中国，其对块体理论的基本原理和分析方法及在岩体工程中的应用进行了系统、全面的介绍 [4]。王思敬等应用矢量分析法对地下工程围岩块体进行分析 [5]。张子新等将赤平投影图解析化，提出了赤平投影解析法 [6] [7]。张奇华进行块体理论软件开发，并应用于工程 [8]。块体理论在发展过程中，在三峡等许多重大工程中都得到了推广应用。本文利用块体理论这一原理分析评价该地下洞室围岩的破坏机制，从而进行洞室围岩的稳定性评价。进行关键块体的搜索计算并判断其稳定性，进而为施工支护提供依据，并使洞库在开挖前能采取相应防治措施，达到提前预报的目的。
2. 工程实例
2.1. 工程概况及地质背景
某地下水封洞库设计库容300 × 104 m3，共设有8条断面，跨度为19 m (宽) × 24 m (高)，长934 m的储油洞室，每两个相邻洞室设为一个洞罐，共分为4个罐，每组洞罐的2条储油洞室之间通过连接巷道互相连通。截面形状为直墙圆拱形，洞室轴线走向为东西方向。
储油洞室围岩岩性主要为微、未风化花岗岩，岩性较为单一，局部见辉绿岩、细晶岩、角闪闪长玢岩等岩脉，场区无区域断裂和活动性断层分布，地壳稳定，局部发育有小断层和节理密集带，场区地下水主要为第四系松散岩类孔隙潜水和基岩裂隙水。
2.2. 结构面发育情况
根据洞室开挖地质素描可知，场区主要存在三组节理，产状分别为① 275˚∠30˚、② 140˚∠85˚、③240˚∠85˚，其中以前两组最为常见，节理张开度一般小于1 mm，节理面多平直、平滑。
① 275˚∠30˚，倾角较缓，在场区分布较为普遍，节理间距一般大于2 m，部分地段改组节理较密，节理面多平直、不规则或平滑，张开度一般小于1 mm，无充填，部分段改组节理张开度可达5~20 mm，粘土质成分或风化花岗岩碎块充填。场区大多数岩脉产状与该组节理相同。
② 140˚∠85˚，该组节理陡直展布，在场区分布最为普遍，节理间距一般大于1 m，部分地段改组节理较密，节理面多平直、平滑状，张开度一般小于1 mm，无充填，部分段改组节理张开度可达5~20 mm，可见绿泥石或风化花岗岩碎块充填。该组节理透水性相对较好，水幕巷道中多数出水点与该组节理有关，深层风化也多沿该组节理发育。场区大多数热液蚀变带和小型破碎带产状与该组节理相同。
③ 240˚∠70˚，该组节理在场区分布不如前两组普遍，节理间距一般大于1 m，局部地段改组节理较密，节理面多平直、平滑状，张开度一般小于1 mm，无充填，局部可见粘土或风化花岗岩碎块充填。场区只有少数岩脉、热液蚀变带和小型破碎带产状与该组节理相同。
2.3. 关键块体软件分析计算
2.3.1. Unwedge程序简介
Unwedge软件是E. HOEK等人基于石根华块体理论编制的地下洞室块体稳定性分析Unwedge程序 [9]。该程序是一种分析在坚硬岩体中开挖形成的块体稳定性的应用分析软件。它假定结构面相切的块体为四边形，即由三组结构面和开挖临空面组成，它仅考虑块体的重力及结构面的力学性质，并假定结构面为平面，岩体的变形仅为结构面的变形，结构体为刚体，结构面贯穿研究区域，且在保持产状不变的情况下任意移动，开挖断面沿轴向方向保持不变。
2.3.2. 关键块体软件分析计算
本文利用块体数值模拟软件Unwedge作为分析计算工具，对某地下储备库储油洞室的围岩块体进行数值模拟，对储油洞室的围岩稳定性系数进行计算根据施工勘察地质素描调查资料显示，该场区洞室围岩主要发育有三组优势结构面，其主要参数见表1，输入Unwedge软件结构面相关参数，得到可移动块体相关计算结果见表2及可移动块体分布结果见图1。
Table 1. Parameters of main structural planes
表1. 主要结构面参数表
Table 2. Calculation results of movable blocks
表2. 可移动块体计算结果
Figure 1. Distribution of movable blocks
图1. 可移动块体分布图
通过对计算结果进行分析，并结合可移动块分析图的相关内容可以看出，2号块体以及4号块体均为稳定块体，其中4号块体稳定性弱于2号块体，5号、7号以及8号块体的稳定性系数依次为1.876、1.457，0，由于地下储油洞室的施工过程中，围岩块体会受到工程爆破等人类工程活动的扰动，且扰动程度较为剧烈，加之地下储油洞室对于块体稳定性有较高的要求，因此为确保在服役周期内洞室的安全性，要求储油洞室围岩的稳定性系数需大于2.0，因此应该采取支护措施对稳定性系数小于2.0的储油洞室围岩块体进行加固。
2.3.3. 支护设计
受结构面切割形成的可移动块体，针对其发育的部位以及稳定性系数而采用不同强度的锚喷支护设计，相关支护参数见表3，支护效果图及支护前后稳定性系数变化分别见图2、表4。
Table 3. Spray anchor support parameter design
表3. 喷锚支持参数设计
Figure 2. Three-dimensional view of spray anchor support
图2. 喷锚支护立体图
Table 4. Changes of stability coefficient before and after support
表4. 支护前后稳定性系数变化表
通过比较支护前以及支护后围岩的稳定性系数可以看出，在采取支护措施后，洞室围岩块体的稳定性系数显著提高。其中，对仅采取锚杆支护作为支护措施的围岩的稳定性系数进行分析可以看出，5号、7号及8号块体在支护后，稳定性系数均满足工程要求，但其中5号、7号块体的稳定性系数仍低于8号块体；为进一步提高围岩稳定性，本文进一步计算了采取喷锚支护后围岩块体的稳定性系数，计算得出7号及8号块体的稳定性系数依次由2.390、3.687提高到4.268、22.203。通过数值模拟结果可以看出，采取喷锚支护可以显著提高洞室围岩的稳定性系数，考虑到喷锚支护具有工艺简单，施工快速的优点，因此在施工过推荐采用喷锚作为围岩支护手段，可以对稳定围岩起到明显而积极效果。
3. 结论
1) 本文根据勘察试验资料给出了不同产状的结构面抗剪强度参数指标，具有一定代表性，但考虑到试验结果的离散性较大，需要工程人员有足够的经验来结合具体地质条件确定相关参数，需要找到一种更好的方法来解决这个问题。
2) 在拟建地下水封洞库中，洞室顶部的关键块体所占比例最大，且优势明显；左侧边墙稳定性较差；右侧边墙稳定性较好。应重点防治洞室顶部和左侧边墙关键块体失稳。关键块体主要以结构面②、③面滑动，因此在洞库开挖施工过程中，针对这②、③结构面的出露处需加强监测。
3) 本文通过Unwedge软件直接生成某地下水封储库三维的关键块体图形，展示了受结构面切割形成的最大不稳定块体其大小和位置，并得出各关键块体的破坏模式及稳定性系数，通过采用针对性的喷锚支护设计使得关键块体稳定性显著提高。