地下水封洞库的原理早在上世纪30年代就提出,并在世界上其它国家截止目前已成功建造了60多处地下水封式储能洞库,但我国仍未完全掌握地下水封储气库的核心技术,究其原因,主要在于国外对于技术的保密。我国已建的汕头、宁波、珠海和黄岛四处地下水封储气库都是由瑞典、法国等国外公司主导建设,对一些核心技术不透露。因此,我国的地下水封储气库建设技术尚处于摸索和试验阶段,对水封条件下洞室群围岩变形破坏特点、规律的认识还不够,对围岩变形破坏的控制理论与方法还不成熟,因此,在我国当前能源地下储备建设热潮时期,开展地下水封储气库围岩变形破坏机制与锚喷支护研究具有重要的理论意义与工程价值。地下水封储气库的工作原理为：利用油水不相混合的特性,合理控制地下洞库围岩中的水压力,使其在洞壁各处始终大于油压或饱和蒸汽压。由此可知,该类工程具有独特的工程特性,例如水幕系统持续地注入高压力水、洞壁高压力水持续地往洞内入渗、液化石油气频繁抽放所产生的动态内压等。这些独特的工程特性决定了水封洞库的建设不能照搬已有的其它地下工程经验和技术,而必须在深入认识建库地质环境条件基础上,开展针对性、创新性且满足工程需求的基础科学问题研究。从工程岩体稳定性角度看,由于存在水幕系统的长期作用,因而稳定性分析中需要考虑水岩耦合作用。考虑到水封洞库都选择在岩性坚硬、结构面不太发育的较高等级质量岩体内,其结构面一般无充填、风化程度低,岩块与裂隙水之间的化学作用轻微,岩块和裂隙之间基本没有水量交换,即水岩耦合主要表现为裂隙渗流-岩石应力之间的力学耦合。再考虑到建库岩体结构面粗糙度一般较低,岩体的剪胀特性不明显,即对裂隙渗流特性影响明显的是其法向应力。地下水封洞库一般都处于浅埋状况,因为从水封原理和经济角度进行埋深设计时,只要地下水位面与洞库的高差大于储存介质饱和蒸汽压外加一定的安全超高即可,因而洞库围岩的应力强度比一般是较低的,即洞库围岩变形破坏模式主要受岩体结构所控制。由于上述地下水封洞库工程的地质特性,其锚喷支护方案也有别于其它地下工程,需要得到更具针对性和合理性的研究,从而为该类工程安全、经济地建设提供理论支撑和依据。本文以烟台万华地下水封储气库工程为基础,采用试验、监测和理论分析相结合,定性分析和定量分析互动并行的手段,进行地下水封储气库围岩变形破坏机制与锚喷支护研究。研究的核心是该类工程中地下水-岩体-支护结构之间的相互作用及由此产生的岩体变形、破坏规律。主要内容有：块状结构裂隙岩体的网络模型,并采用计算机技术对其进行表征：在现场水文试验、地应力测试及钻孔影像的基础上,通过理论分析,建立了洞库区结构面等效水力隙宽与法向应力之间的关系；基于现场监测与数值模拟,探讨水封条件下大型地下储气库围岩变形破坏模式、发生条件,揭示大型地下水封储气库围岩变形破坏机制;基于Q系统与数值模拟,优化围岩锚喷支护方案,并基于现场监测数据进行支护效果评价。论文完成的主要工作包括以下几个方面：(1)基于钻孔摄像原理与空间信息架构关系,对钻孔摄像技术应用于倾斜孔中进行结构面产状测量的解译方法进行了推导,得出了精确的解析计算公式,并编写了相应的计算软件。理论测试表明,推导得出的倾斜孔内结构面产状计算的解析公式及其程序是正确的,可以为相关工程应用提供支撑和参考。采用PVC管模拟钻孔进行了钻孔摄像技术的可靠性试验,结果表明：当PVC管的倾伏向270°、倾伏角25°时,结构面产状的可靠度为0.10；倾伏向176°、倾伏角60°时,结构面产状的可靠度为0.67；不同PVC管方位下的结构面产状可靠度差异主要为PVC管倾伏角的不同所导致；对一般岩体工程须在倾伏角不小于71°的钻孔内进行测量。(2)对迭代自组织数据分析技术算法进行改进,再将模糊隶属度和模糊加权指数引入其中,并在欧氏距离和标准偏差两个关键算式中增加指标权重,构建起可考虑多项分组指标及其重要程度的结构面自组织模糊聚类分组模型,推导了该分组模型下的结构面间距计算式,并编写了相应的计算软件。理论和实践表明该分组模型具有较好的普适性、自适应性和自组织性。(3)采用Excel生成伪随机数,AutoCAD进行三维结构面圆盘建模,成功地实现了结构面三维网络模拟,指出岩体结构均质区、基于工程地质条件与实测数据的网络模型校核的重要性,指出了采用实测数据对结构面三维网络模拟进行校核时,不能校核各个结构面的具体位置,而只能对比一些统计参数,因为结构面三维网络模拟只是统计意义上的仿真,并非真实场景的再现,后者错误认识是目前滥用结构面三维网络模拟的根源。(4)由于计算机硬件的限制,目前对于数量庞大的结构面模型一般无法运算或即使能运算但由于数值截短误差而难以得到可靠的结果,因此通过设置直径阈值、产状差异阈值与中心点距离阈值的方式剔除短小、密集的结构面,在离散元模型建立之后,通过设置体积阈值删除易影响计算速率、易畸形的小型块体,最后通过给同一结构面赋予不同参数的方式来表征虚拟与真实结构面,由此成功地实现了基于随机模拟的岩体结构模型离散元表征,此表征算法全部利用3DEC软件中的FISH语言编程实现。(5)提出了一种结构面等效水力隙宽反演方法。该方法以钻孔影像、钻孔压水试验和渗透张量理论计算为基础,首先基于试算法确定折减系数,初步确定压水试段各结构面的等效水力隙宽；再对不确定性较大的等效水力隙宽采用遗传算法进一步优化,以最终获取更为准确的数值。其中裂隙网络计算法以单个的结构面为单位进行计算,而无需对结构面进行分组,提高了计算结果的可靠性。该方法能同时反演大量的结构面等效水力隙宽。(6)统计分析表明,平均隙宽是初始等效水力隙宽70～100倍的压水试段数占总试段数的77.5%,表明结构面几何隙宽与等效水力隙宽差异悬殊,在岩体水力学及裂隙渗流研究中需要特别注意这种差别：场区结构面等效水力隙宽分布在10-100um之间的数量约占总数的99%,且分布最密集的区域是10～45um,约占总数的65%;当法向应力大于7MPa时,等效水力隙宽与法向应力的联系不明显,而小于7MPa时,等效水力隙宽随法向应力增大而变小的现象明显。(7)基于E. Hoek提出的GSI和Q值之间的关系式以及Hoek-Brown强度准则,结合室内岩石力学试验资料,提出不同Q值岩体的力学参数估算公式,并以Q值为4岩体为例,估算出其力学参数,再与施工巷道实际监测资料进行对比,验证了取值方法的合理与可靠性。基于室内结构面直剪试验,得出了相应的结构面力学参数。以含竖井、水幕巷道、支洞及主洞室的洞室群单元作为研究对象,考虑Q值为4的围岩状况,不考虑水岩耦合作用,按照实际开挖次序进行离散元数值模拟研究,结果表明,各地下结构之间的距离是合理的,相互之间基本无影响,围岩整体变形破坏机制主要受二次应力的影响,围岩位移与塑性区较小(8)以主洞室模型为研究对象,考虑裂隙流-应力耦合作用,进行离散元数值模拟研究,结果表明,主洞室左洞肩出现失稳块体的主要原因在于结构面与临空面之间的有利组合,失稳块体处的流场与其余地方明显不同,地下水流量与流速在该处急剧增大,这主要是隙宽在拉张应力作用下增大所导致的渗透系数增大的缘故,同时由于开挖面处水压力为0,组成失稳块体的结构面水力梯度较大,因而导致流速与流量变大。总体分析认为,出现失稳块体这种局部破坏的主要原因在于结构面的有利组合,地下水在其中起促进与加速失稳产生的辅助性作用。并着重指出,在处理失稳块体的同时及时采取合理的注浆堵水措施是十分必要的。(9)基于地下水封储气库工程条件,基于Q系统得出围岩锚喷支护参数,并分别从连续围岩与非连续围岩角度进行参数的校核,结果表明,连续围岩下Q系统确定的参数是安全的,非连续围岩下由于出现失稳块体,Q系统确定的参数是不安全的。综合认为,基于Q系统所得的锚喷支护参数不针对局部块体失稳问题,因而需结合离散元数值模拟等进行综合设计,并且在施工期要特别重视对可能发生的局部失稳问题进行提前甄别。(10)基于Q系统确定的锚喷支护参数,选定锚杆纵向间距、长度及喷层厚度作为优化指标,采用喷层位移、喷层受力、锚杆受力与节点破坏情况作为支护效果评价指标,基于正交设计得出试验方案,对每种方案进行数值模拟,基于进化神经网络建立起支护安全性评价指标与锚喷优化指标间的非线性关系,在确定支护安全性评价指标合理数值的基础上对锚喷优化指标进行优化,结合该工程实际进展情况,表明优化后的参数对于特定Q值区段是安全、可靠的,可为类似工程提供依据与参考。本文主要创新成果有：(1)目前钻孔摄像技术探测结构面产状时仅限于竖直孔,基于空间架构关系与钻孔摄像原理,提出了倾斜孔中结构面产状的解译算法,并试验研究了该技术的可靠性。在结构面组数划分中,解决了能考虑多项分组指标及其重要程度的分组模型问题。这极大地扩充了结构面调查与分析方法。(2)针对目前结构面三维网络模拟的应用仅停留在几何层面的问题,基于3DEC软件与FISH语言,提出了基于随机模拟的结构面三维网络模型的离散元表征算法,从而使三维网络模拟能与力学计算衔接,且算法无需多建虚拟结构面,提高了计算效率。(3)目前工程中对岩体结构面等效水力隙宽难以确定,提出了基于钻孔摄像技术与钻孔压水试验,采用试算法与遗传算法反演等效水力隙宽的方法,该方法可同时反演整个压水试段的结构面等效水力隙宽,可靠性较高。再通过场区地应力回归分析,总结出了结构面等效水力隙宽-法向应力之间的关系式,为裂隙流-应力耦合研究奠定了基础。(4)目前地下水封储气库围岩变形破坏机制问题在国内研究较少,结论不系统,主要原因在于此类工程很少,从洞室群单元模型、考虑裂隙流-应力耦合的主洞室模型两种角度出发,通过大量的离散元数值计算,结合现场监测资料,进行了围岩变形破坏机制分析,重点区分了二次应力、岩体结构与地下水在围岩变形破坏过程中的作用。(5)针对以往锚喷支护参数优化时易于陷入局部最优值的弊端,提出了基于Q系统初选锚喷支护参数,再通过正交设计、进化神经网络与数值计算进行优化的锚喷设计新方法。该种方法易于工程应用,在整个工程区内具有较好的普适性。