谢晔,荀晓慧,魏玉峰,聂德新
(成都理工大学,成都 610059)
受地质运动和构造应力作用,天然岩体中存在着大量可见或微裂隙[1]。对岩体结构和其物理力学特征进行研究是岩体力学的基本问题。经过数百万年的长期地质演化和复杂构造,岩石受到累积应力的作用,存在着各种不同尺度的不连续面,从中、微裂缝到大型断层和褶皱。这些软弱带由于其强度、刚度和抗变形能力较低,对岩石的完整性乃至整体力学性能有很大的负面影响。尤其针对碎裂岩体此类力学性质较差的岩体,对其工程特性进行研究,是评价岩体是否适用于工程实践的关键[2]。
在过去的几十年里,国内外许多学者们对于碎裂岩体的工程特性有着广泛关注。聂德新等[3-5]首先于1994年提出了岩体强度再生理论,认为作用于软弱结构面上的正应力可以一定程度上改善软弱夹层的力学性质,生成与新的地质环境相适应的结构,该结构特征下岩体的力学性质将会恢复至新的强度;郭启良[6]认为,尽管岩体的结构与完整程度各不相同,但在原地状态下大都具有一定的承载力;Heping Xie等[7]和Pei J L等[8]对天然裂隙大理岩进行了单轴压缩试验。认为天然裂隙岩石的破坏特征受裂隙分布和裂隙形态的控制,是一个局部渐进的破坏过程,最终是突然破坏。尽管天然岩石比类岩石材料更难达到理想的效果,但天然岩石的试验研究能更好地揭示实际岩石的工程特性。
在岩体的等效力学特征研究上,Singh[9-10]提出了节理岩体的连续介质表征方法,并给出了等效连续介质各向异性岩体弹性模量的估算公式。Zienkiewicz等人[11]使用了等效连续体方法(称为多层模型)来模拟不连续岩体;Gerrard[12]还使用了等效连续体方法,将单元的柔度表示为完整岩石的柔度和单个节理组的柔度之和;Fossum[13]提出了估算随机节理岩体有效弹性性质的关系式;Chen[14]基于弹性岩石基质和具有非线性法向和剪切响应的节理组之间应变分配的连续介质假设,提出了节理介质应力-应变响应的速率方程;Duncanh和Goodman[15]将节理岩体描述为等效的各向异性连续体;Amadei和Goodman[16]在对节理岩体建模的尝试中,在为裂隙岩体开发的本构关系中隐含地考虑了节理的影响;Wei和Hudson[17]开发了一种等效材料方法,根据不连续面的几何和力学性质估算变形模量。
在关于碎裂岩体的力学试验研究方面,目前最常用的方法就是采用相似材料进行配比来代替岩石材料进行试验。如廖美春[18]采用石膏模型来代替岩体对深埋隧洞(厂房)开挖时产生的分区破裂现象进行物理模拟;也有学者[19-20]根据相似原理设置物理模拟材料,通过插入铝皮设置节理缝隙,探究了经锚杆加固后碎裂岩体的力学特征;印度的Tiwari[21]通过石灰与细砂经过一定配比制作岩石试样进行三轴试验。
对于非原岩材料制样的模拟实验来说,试验结果并不能真实地反映岩体的真实参数,怎样才能得到岩体的真实参数是一个需要解决的问题。对于大多数岩石工程,现场地质条件总是复杂多变的,现场原位试验总是难以控制[22]。此外,由于不同位置工程岩体的地质特征不同,现场试验结果往往反映局部场地,即试验结果仅对试验区域有效,而不适用于整个工程区域。因此,在有关天然状态下碎裂岩体的工程特性方面的研究,目前少有考虑岩体的赋存环境。本文通过分析碎裂岩体的成岩过程,结合国内外沉积物压密试验成果,以公伯峡水电站碎裂岩体为研究对象,开展室内等效压密实验,从而对碎裂岩体的等效工程特性进行探讨。
图1 向家坝左岸挤压带碎裂岩体
根据谷德振先生在《岩体工程地质力学基础》[23]提出的岩体结构理论认为,受构造或风化作用影响、由岩屑和碎屑组成的结构面间距不超过50 cm的岩体为碎裂岩体。而《水力发电工程地质勘察规范》(GB50287-2006)[24]将碎裂结构岩体又进一步细分为两大类:结构面间距在10~30 cm之间的为块裂结构岩体和结构面间距小于10 cm的为碎裂结构。上述两种碎裂岩体结构划分依据均为岩体结构面间距,而在构造或风化作用强烈地区的碎裂岩体,其结构特征往往受大量发育、不规则的微小结构面控制,这些微小的裂纹尚未贯通扩展为裂隙面,但在开挖、爆破等冲击震动作用下极易相互贯通,形成形状各异的短小裂隙,裂隙之间经过相互交切,使岩体显得十分破碎。针对此类破碎较强的岩体,其结构面间距测量工作开展往往相对困难,为了进一步获取碎裂岩体的结构尺度特征,本文对向家坝左岸挤压带(图1)、溪洛渡玄武岩层间错动带(图2)和白鹤滩水电站柱状节理碎裂岩体(图3)进行取样分析,通过室内筛分实验,获取不同电站碎裂岩体级配特征如图4~6。
图2 溪洛渡玄武岩错动带
图3 白鹤滩柱状节理碎裂玄武岩
图4 白鹤滩水电站碎裂岩体级配曲线
图5 溪洛渡水电站碎裂岩体级配曲线
图6 向家坝水电站碎裂岩体级配曲线
从不同水电站坝址区碎裂岩体级配曲线中可以看出,试样中块度在10 cm以下的土重含量均大于80%,且基本不含黏粒,该试验结果与上述关于碎裂岩体块度的规定有着很好的对应性。因此,结合现场调查结果,对于碎裂岩体结构也可以通过岩体块度指标进行定义为:由构造作用形成的、块度在10 cm以下的岩块或碎屑;由风化卸荷形成的块度在10~30 cm之间的岩块或碎屑。
碎裂岩体由于其特殊的结构性易受扰动而解体,因此在取样过程中稍有碰撞或冲击就会破坏其原状性质。为了尽量获取岩体真正原位状态下的物理指标,避免钻探过程中冲洗液浸泡导致岩体含水率增大,利用减震法将错动带上层完整岩体削去,采用团队新开发的原位取样技术,获取溪洛渡水电站和白鹤滩水电站坝址区层内错动带以及向家坝水电站左坝肩破碎带碎裂岩体未松弛状态下的原状试样(图7)。通过上述取样方法获取的试样是在打开岩盖后即刻进行的,且采用了特质的材料进行包裹,能够有效避免岩体长时间暴露后吸水或者松弛,通过室内比重试验、及烘干实验,获取不同坝址区碎裂岩体原状试样密度及含水率如表1~3所示。
图7 碎裂岩体原状样取样过程
由上述试验结果可以看出,原位状态下的碎裂岩体具有较高的密度和较低的含水率,天然密度在2.3~2.5 g/cm3之间,干密度在2.3~2.4 g/cm3之间含,水率基本在3%~6%之间。且从溪洛渡水电站收集到的试样室内实验结果(图7)可以看出,随着距离洞壁深度的增加,岩体的含水率逐渐降低,密度逐渐增大。这是由于位于深部的岩体未受扰动,避免了表部风华剥蚀,能够基本保持早期构造或重力压密的状态。由此可以看出,岩体的物理指标很大程度受岩体所处的地质环境影响,深埋状态下的碎裂岩体应具有相对较高的密度和较低的含水率特性。
表1 溪洛渡碎裂岩体物理性质
表2 白鹤滩碎裂岩体物理性质
表3 向家坝碎裂岩体物理性质
构造运动过程中岩体在高地应力条件下会产生褶皱和断层,断层错动继而引起岩体发生变形和破碎,在剪胀作用下碎裂的岩块之间出现空隙,新出现的空隙使原有的岩体体积变大,随着时间推移,当后期的构造运动不再剧烈或构造应力达不到形成断层时的应力量值时,此时岩体受上覆新的沉积物自重应力作用而被逐渐填充、挤压,从而体积减小、密度增大。通过总结国内外有关自重压密的研究成果,获取不同沉积物盖层厚度-干密度关系曲线图、孔隙比-压力曲线图如图8。
在不同盖层厚度下孔隙比和盖层压力关系曲线如图8所示,方程对应的干密度和孔隙比指标见表4。
表4 国外不同沉积物孔隙比、压力关系汇总表
图8 国内外不同盖层厚度物理指标关系曲线
从收集的资料可以看出,当沉积物干密度在2.3~2.4 g/cm3之间时,所承受的压力在50~60 MPa左右,因此可以认为,现场获取的试样在历史上曾受到的压力达50~60 MPa。
由于组成碎裂岩体的岩块或碎屑本身的抗变形能力较大,通常认为其是不可压缩的刚体,因此碎裂岩体压缩的实质可以被看作是岩块或岩屑之间空隙被填充或挤压的过程,岩体的孔隙比计算公式可以表示为:孔隙比=空隙体积/块体体积。
通过对公伯峡水电站右岸坝址区揭露的岩层进行现场调查发现,临近古风化花岗岩中随处可见嵌入的第三系红色砾砂岩(图9),对二者进行物理指标参数试验发现其密度基本一致。从两种不同的岩性成岩过程来看,红色砾砂岩是在临夏组地层沉积时填充至断层或褶皱中,随古风化花岗岩一同经历重力压密作用,即使花岗岩无论从结构还是构造上密度都要比砾砂岩要低,但经历了相同历史时间的压密和重分布,最终导致两者在密度和物理力学性质上保持高度的相似性。由此可以看出,虽然天然压密条件下存在一定的时间效应,原位碎裂岩体的强度和压缩量更大,但同时可以作为室内等效压密实验结果的安全冗余度来考虑,即天然碎裂岩体的室内等效压密与天然地质条件下的压密结果应该是具有等效性的。
图9 公伯峡水电站右岸揭露的岩层
根据实际工程调查结果可以看出,原岩状态下的碎裂岩体即使被大量发育的随机节理相互切割,但在长时间重力压密条件下,如果未发生错动,其结构应当仍具有一定的完整性和较高的强度及模量。通过室内高精度压密试验,理论上可以对相同碎裂程度的岩块进行等效压密恢复至原岩状态下岩体的物理状态,即给予岩体历史上经受的天然压力进行压密,应该能够获取与野外原位状态下相近的密度、孔隙比等物理指标。根据胡克定律:
σ=E·Es
(1)
式中,σ为应力;E为弹性模量;Es为压缩模量。
(2)
式中,P为压力(MPa);S为不同压力下对应的压缩量(应变)。
根据前述等效压密原理,在室内高精度加压条件下,是否能获得与天然地质条件下碎裂岩体相同或接近的压缩模量,是检验室内等效试验结果的可靠依据。
试验采用是团队所在单位国家重点实验室引进的美国MTS岩石试验机,仪器主要由伺服机操作系统、加压系统组成,其中,试验伺服控制系统可以记录试样加载全过程压力、变形曲线。试验加压装置采用特质的具有极高强度、抗变形能力达1 000 MPa的模具及加压活塞,模具外径250 mm,内径与加压活塞尺寸相同,为50 mm。加压过程如图10所示。由图10(e)、(f)可以看出,加载后的碎裂岩体试样在高压作用下展现出与原状岩体极为相似的形貌。
图10 试验过程
试验过程中的砾砂岩与花岗岩各项物理力学指标变化曲线如图11所示。从图11中的过程曲线可以看出,无论是孔隙比、干密度、压缩模量和变形模量,随压力的变化曲线规律性均较好,拟合出的曲线方程相关度均较高。通过查阅相关地质资料,公伯峡水电站盖层分布厚度在1 900~2 000 m左右,根据前述物理指标试验结果,若取干密度在2.4~2.5 g/cm3之间,则压力可在48~50 MPa之间,考虑该地区曾受喜马拉雅构造运动作用,以10 MPa构造应力为冗余度,可认为公伯峡水电站第三系地层及古全风化花岗岩曾受50~60 MPa压力。将室内试验结果、国内外研究资料以及现场研究资料中50~60 MPa压力下各项指标进行对比,如表5所示。
表5 等效实验、现场测量与国内外研究成果对比
从表5中可以看出,在50~60 MPa压力下,等效试验中古全风化花岗与砾砂岩各项物理指标基本接近,证明二者在相同量级的压力下可以最终达到相同的物理力学性质;室内试验结果与现场测试获取的结果具有较好的对应性,无论是干密度与变形模量基本能够保持一致。其中等效试验结果中无论是花岗岩还是砾砂岩,均与国外资料中的下线值有着较高的匹配度,从而证明了开展室内等效实验压密实验的准确性。
本文通过不同水电站碎裂岩体的物理指标和等效压密试验,对碎裂岩体的等效工程特性进行研究,主要获得了以下研究成果:
(1) 获取了不同电站的碎裂岩体粒度成分、密度、含水率、孔隙比等物理指标,从块度角度重新定义了碎裂岩体结构的划分标准。
(2) 收集了国内外岩体自重压密试验成果,获取了不同沉积物厚度与密度、压力、孔隙比等指标的关系曲线。
图11 等效压密试验力学指标过程曲线
(3) 通过对公伯峡水电站古风化花岗岩和第三系砾砂岩进行室内等效压密试验,验证了室内压密试验获取的各项物理指标与原位状态下碎裂岩体的物理指标具有一定的等效性。