高应力下花岗岩耦合节理在剪切过程中渗透率演化特性

节理滑移时其渗透率的演化特性是岩体水力特性重要的科学问题之一。为了深入了解花岗岩节理在高应力条件下剪切滑移时渗透率的变化规律,对花岗岩耦合节理剪切时的渗透率进行了试验研究。为了解决剪切时高压孔隙流体的密封问题及节理渗透率测量问题,设计了一种适用于高压孔隙流体作用下的剪切渗流试验系统。通过剪切渗流试验获得了渗透率随剪切位移的变化规律,以及剪切应力和法向变形与剪切位移的关系。试验结果表明:高应力条件下剪切峰后节理发生明显的剪胀,但其渗透率呈降低趋势;节理剪切过程中节理的剪胀、压缩和凸起的剪断磨损机制共同影响渗透率的变化;高应力一方面限制了节理的剪胀特性,另一方面也增强了凸起的剪断磨损机制,因此,节理剪切时节理渗透性被显著降低。     

节理岩体蠕变特性研究

在已知岩石和节理蠕变规律的前提下,推导了节理岩体蠕变模型的一般表达式。假定岩石体积变形和节理法向压缩变形为弹性变形,忽略节理的剪胀现象,认为只有岩石畸变和节理剪切滑移与时间有关,从而推导了含三组相交节理的岩体蠕变模型及其参数。根据反演出的岩石和节理蠕变模型,计算了含三组相交节理的岩体在单轴应力作用下一些节理参数对岩体单轴蠕变的影响。分析表明,节理间距、剪胀系数以及节理夹角都对岩体的单轴蠕变变形有明显影响。节理间距越大,剪胀系数越大,节理夹角越小,节理岩体的单轴蠕变柔量也就越小,岩体的蠕变变形也越小。     

岩桥直剪细观破坏特征与剪胀效应研究

节理岩体的剪切贯通机制影响着边坡的稳定性。为揭示锁固段型非贯通节理岩体在不同连通率和法向应力下的破坏特征,在室内直剪试验中结合高速摄影与AE特征参数分析其剪切全过程及剪胀效应。结果表明:节理岩体直剪试验中,法向应力的增大及节理连通率的下降会致使峰值剪切应力及峰值剪切位移增大;节理连通率与法向应力对其破坏特征具显著影响,表现为节理连通率较高且法向应力较小时呈直接剪断的特性,节理连通率降低后呈拉剪复合破坏,出现剪胀现象,而法向应力的增大使得剪胀效应呈波动现象;AE特征与岩桥贯通过程一致,事件数峰值随节理连通率的降低及法向应力的增大而增大且位于峰后。试验得到的岩桥细观破坏特征与剪胀效应对研究锁固段型岩质边坡的贯通破坏机制具指导意义。     

压剪作用下裂隙扩展过程渗流与应力耦合模型研究

裂隙的剪胀特性及扩展演化规律对岩体的渗流特性具有重要影响。为揭示裂隙剪胀及扩展演化对岩体渗流的影响,基于残余强度提出了一种能较好描述岩石硬化-软化特性的全剪切本构关系;结合剪切变形与裂隙开度的关系,利用最小势能原理和立方定律,建立了压剪作用下考虑裂隙剪胀特性的渗流应力耦合模型;假定压剪作用下裂隙发生Ⅰ型扩展,提出了伴有翼型裂纹的渗流模型,该模型不仅考虑了岩石的剪胀特性,更反映了裂隙扩展过程渗流的演化规律。对不同裂隙粗糙度的剪切应力-位移曲线进行分析,全剪切本构模型表现出较高的拟合精度。在剪切应力-位移关系基础上,通过剪切渗流试验数据对压剪作用下渗流模型进行验证,结果表明,该模型能较好地描述岩体剪胀阶段渗透性演化规律。利用等效裂隙简化裂隙网络,并通过试验数据进行验证,证实了裂隙扩展过程渗流-应力耦合模型的准确性与适用性。     

动载下节理岩体破坏过程的数值试验研究

节理岩体是工程中最常见的一类岩体,其在地震、爆炸等动载下的力学响应及破坏过程对相关工程安全性的影响至关重要。采用基于有限元应力分析和统计损伤理论开发的动态版RFPA2D数值模拟软件,对动载下节理岩体的动态破坏过程进行了模拟,重点讨论了节理条数、节理贯通度、节理倾角及应力波峰值对岩体动态破坏过程的影响规律。计算结果表明,断续节理岩体动态破坏过程及破坏强度与节理构造形态、应力波峰值密切相关。相同动载下,随着节理条数的增加,岩体破坏程度以及应力波能量损失增强,但当节理条数数超过一定值后,岩体破坏程度及应力波能量损失逐渐趋于稳定;节理贯通度较小时,岩体破坏程度较低且破坏单元自上而下均匀分布。随着节理贯通度的增加,岩体破坏增强,且破坏主要出现于节理上部岩体;节理倾角较小时,节理上部岩体破坏严重,易形成次生贯通裂纹。随着节理倾角增加,破坏范围逐渐变大,不易形成次生贯通裂纹;倾角为45°～60°时,岩体破坏效果最佳;动载荷的峰值越大,试样的破坏越严重。当峰值达到一定值时,节理附近发育出多条裂隙并向上下方不断发展而导致岩体完全破坏。在不同节理贯通度工况下与岩石霍布金森压杆（SHPB）试验结果进行比较,结论吻合,证明该数值模拟的合可行性和结论的可靠性。     

裂隙岩体渗透性影响因素的研究

采用UDEC离散单元法中关于裂隙岩体开挖模拟及水力全耦合分析模型,分析裂隙岩体洞室开挖后,因围岩应力与水力耦合作用导致裂隙隙宽变化及渗流变化的过程。为了更直观地了解耦合作用对裂隙岩体渗透特性的影响,以隧洞开挖为例,用开挖后隧洞内总涌水量来表征岩体的渗透特性。利用数值试验的方法,研究了块体边界大小、初始应力比、裂隙隙宽和裂隙夹角对开挖后隧洞内涌水量变化的影响,进而可以看出它们对裂隙岩体渗透性的影响。并得出如下结论：随着块体尺寸和初始应力比的增大,隧洞内总涌水量减少;随初始隙宽的增大涌水量增加并当达到某一固定值时保持不变;隧洞涌水量在θ2/θ1=3.5,其中θ1=30°,即两组节理的夹角为75°处达到最大。     

采动裂隙岩体压剪渗透规律试验研究

根据采动过程中裂隙岩体的应力变化,用法向载荷、剪切载荷和渗透水压分别模拟开采过程中采场的水平应力、垂向应力和水头压力,应用JAM-600型剪切-渗流耦合试验系统对裂隙岩体进行压剪-渗流藕合试验,探讨在恒定法向荷载(CNL)和法向刚度(CNS)条件下,裂隙岩体的法向载荷、裂隙粗糙度与渗透水压对试样的位移、应力和渗流性的影响规律,分析剪切位移大小和岩体裂隙的剪胀特性对裂隙岩体的剪切应力、法向位移、节理水力开度及渗透率影响规律。研究表明:剪切应力和水头压力对裂隙水力开度起促进作用,水平地应力对水力开度变化起抑制作用。随着剪切位移变化,水力开度可分为变小或持平、增大、稳定3个阶段。裂隙表面粗糙度越大,裂隙岩体的刚度越小,则水力开度最终稳定值越大。由于裂隙岩体的剪胀作用渗透率先变小后增大,剪切位移增大,渗透率增大;法向荷载增大,试件的渗透率越小;裂隙表面越粗糙,其渗透率越大,其研究结果可为岩体透水通道形成时的孕育、萌生和爆发的导渗灾变演化过程提供理论基础。     

区域非稳定动力学环境下的岩体松动效应

区域非稳定动力学环境下,长期的构造变形、重力卸荷以及地震动力作用的共同影响,可以导致岩体发生大范围变形、松动。松动岩体内发育大量的软弱结构面,且表现出整体破碎、松弛严重、透水强烈,张性节理裂隙发育、地表裂缝较发育、岩体地震动力破坏信息反映明显等特征。岩体的变形松动可以分为卸荷变形松动、倾倒变形松动、顺层滑移松动、断层控制松动、节理裂隙控制松动等5种模式。岩体的物理振动试验结果表明,地震动力是造成岩体松动的主因。数值模拟结果表明,在单纯自重应力影响作用下,松动岩体不会出现大面积的失稳破坏现象,但在地震动力作用下,松动岩体会发生大面积屈服破坏。     

三维应力作用下岩体单个裂隙的渗流特性分析

天然岩体中存在着大量的孔隙和裂隙,这些缺陷不仅改变了岩体的力学性质,也严重影响了岩体的渗流特性。在对现有裂隙岩体渗流特性研究成果进行分析的基础上,讨论了岩体单个裂隙的力学性质和渗流对单个裂隙岩体产生的力学作用,研究了岩体单裂隙渗透系数与岩体三维应力的关系,考虑了裂隙粗糙度对渗流的影响以及不同方向应力对渗透性影响的差异,分析了单个裂隙岩体在三向应力作用下的渗流特征,得出了裂隙所受三维应力与渗透系数关系式,认为垂直于裂隙面的应力对岩体渗透性起主导作用,岩体渗透系数随垂直裂隙面应力的增加而迅速减小。通过与渗流规律试验结果对比分析,证明了所得单个裂隙岩体渗透系数表达式的正确性。     

基于3D打印的含复杂节理岩石直剪特性及破坏机制研究

岩体中节理的几何形态及力学特性是影响其剪切力学特性及破坏模式的重要因素之一。基于3D打印技术,建立了不同节理粗糙系数(JRC)的节理模型、几何形态节理模型和复杂裂隙网络物理模型,通过开展室内直接剪切试验分析了各组试件的剪切强度及破坏模式。结果表明节理模型的抗剪强度随JRC波动性较大,波动幅值越高,峰值剪切位移越低;平面形节理模型的峰值抗剪强度最低,矩形节理模型的峰值抗剪强度最高,正弦形和三角形节理试件的抗剪能力相近;离散裂隙网络模型和粗糙裂隙网络模型的峰值抗剪强度显著低于实心试件,考虑了节理粗糙性的裂隙网络模型抗剪强度高于直线型节理模型;实心试件破坏模式为典型脆性剪切破坏,裂隙网络模型的破坏模式相对复杂,沿着剪切方向主剪切裂面波动萌生,破断面由多个节理面的交叉点破坏与沿节理面的滑移构成。研究成果可以为3D打印技术的推广和复杂节理岩体剪切力学特性的室内试验研究提供参考。     

Flow–stress coupled permeability tensor for fractured rock masses

In this paper a new analytical model is proposed to determine the permeability tensor for fractured rock masses based on the superposition principle of liquid dissipation energy. This model relies on the geometrical characteristics of rock fractures and the corresponding fracture network, and demonstrates the coupling effect between fluid flow and stress/deformation. This model empirically considers the effect of pre‐peak shear dilation and shear contraction on the hydraulic behavior of rock fractures and can be used to determine the applicability of the continuum approach to hydro‐mechanical coupling analysis. Results of numerical analysis presented in this paper show that the new model can effectively describe the permeability of fractured rock masses, and can be applied to the coupling analysis of seepage and stress fields. Copyright © 2008 John Wiley & Sons, Ltd.     

Polyaxial stress-dependent permeability of a three-dimensional fractured rock layer

A study about the influence of polyaxial (true-triaxial) stresses on the permeability of a three-dimensional (3D) fractured rock layer is presented. The 3D fracture system is constructed by extruding a two-dimensional (2D) outcrop pattern of a limestone bed that exhibits a ladder structure consisting of a “through-going” joint set abutted by later-stage short fractures. Geomechanical behaviour of the 3D fractured rock in response to in-situ stresses is modelled by the finite-discrete element method, which can capture the deformation of matrix blocks, variation of stress fields, reactivation of pre-existing rough fractures and propagation of new cracks. A series of numerical simulations is designed to load the fractured rock using various polyaxial in-situ stresses and the stress-dependent flow properties are further calculated. The fractured layer tends to exhibit stronger flow localisation and higher equivalent permeability as the far-field stress ratio is increased and the stress field is rotated such that fractures are preferentially oriented for shearing. The shear dilation of pre-existing fractures has dominant effects on flow localisation in the system, while the propagation of new fractures has minor impacts. The role of the overburden stress suggests that the conventional 2D analysis that neglects the effect of the out-of-plane stress (perpendicular to the bedding interface) may provide indicative approximations but not fully capture the polyaxial stress-dependent fracture network behaviour. The results of this study have important implications for understanding the heterogeneous flow of geological fluids (e.g. groundwater, petroleum) in subsurface and upscaling permeability for large-scale assessments.     

Permeability tensor and representative elementary volume of fractured rock masses

Based on a simulation of three-dimensional fracture networks and a superposition principle of liquid dissipation energy for fractured rock masses, a model of the fracture permeability tensor is proposed. An elastic constitutive model of rock fractures, considering fracture closure and dilation during shearing, is also proposed, based on the dilation angle of the fracture. Algorithms of flow-path searching and calculation of the effective flow coefficients for fracture networks are presented, together with a discussion on the influence of geometric parameters of the fractures (trace length, spacing, aperture, orientation and the number of fracture sets) on magnitude, anisotropy of hydraulic permeability and the size of a representative elementary volume (REV). The anisotropy of hydraulic permeability of fractured rock masses is mainly affected by orientation and the number of fracture sets, and the REV size is mainly influenced by trace length, spacing and the number of fracture sets. The results of studies on REV size and the influence of in-situ stress on hydraulic conductivity of the rock mass on the slope of Jinping-I hydropower station, China, are presented using the developed models and methods. The simulation results agreed well with the results obtained from field water-pressure measurements, with an error of less than 10 %.     

Hydraulic Characteristics of Rough Fractures in Linear Flow under Normal and Shear Load

Summary A hydro-mechanical testing system, which is capable of measuring both the flow rates and the normal and shear displacement of a rock fracture, was built to investigate the hydraulic behaviour of rough tension fractures. Laboratory hydraulic tests in linear flow were conducted on rough rock fractures, artificially created using a splitter under various normal and shear loading. Prior to the tests, aperture distributions were determined by measuring the topography of upper and lower fracture surfaces using a laser profilometer. Experimental variograms of the initial aperture distributions were classified into four groups of geostatistical model, though the overall experimental variograms could be well fitted to the exponential model. The permeability of the rough rock fractures decayed exponentially with respect to the normal stress increase up to 5 MPa. Hydraulic behaviours during monotonic shear loading were significantly affected by the dilation occurring until the shear stress reached the peak strength. With the further dilation, the permeability of the rough fracture specimens increased more. However, beyond shear displacement of about 7 to 8 mm, permeability gradually reached a maximum threshold value. The combined effects of both asperity degradation and gouge production, which prohibited the subsequent enlargement of mean fracture aperture, mainly caused this phenomenon. Permeability changes during cyclic shear loading showed somewhat irregular variations, especially after the first shear loading cycle, due to the complex interaction from asperity degradations and production of gouge materials. The relation between hydraulic and mechanical apertures was analyzed to investigate the valid range of mechanical apertures to be applied to the cubic law. Received June 12, 2001; accepted February 26, 2002 Published online September 2, 2002     

法向卸荷下贯通裂隙岩样压剪破坏特征研究

开展贯通裂隙岩样在法向卸荷下的压剪破坏研究,有助于更好地理解岩体裂隙在卸荷条件下的失稳特征。对单裂隙贯通岩样开展剪切应力峰前一系列法向卸荷破坏试验,分析了法向卸荷下试样裂隙的变形及能量演化规律。结果表明,法向卸荷下裂隙压剪失稳时的法向应力大于常规直剪试验的对应值,裂隙的抗剪能力减弱;剪切位移随法向应力减小而增大;卸荷变形比K随裂隙轮廓面积比Rs的增加而增大;卸荷过程中总变形能U0由减小转变为增加的现象可用于预测岩石裂隙卸荷失稳。研究结果对理解贯通裂隙岩体压剪卸荷失稳破坏具有一定参考意义。     

岩体裂隙网络等效渗透系数方向性的数值计算

采用渗流力学理论并结合Monte Carlo方法描述岩体裂隙的随机分布,研究渗流模型的尺寸效应并确定表征单元体积（REV）,得到了3种开口度分布形式的等效渗透系数椭圆曲线,建立了等效渗透系数方向性的判别标准。离散裂隙网络（DFN）模型假定流体只在岩体裂隙内部流动,而不通过岩体本身渗流。基于二维离散元程序UDEC并进行二次开发,建立DFN模型,通过改变流体的流动方向,得到不同流动方向下岩体裂隙网络的等效渗透系数,并分析不同的开口度分布形式对岩体裂隙网络等效渗透系数方向性的影响。计算结果表明,表征单元体积存在的条件是等效渗透系数保持稳定且渗透椭圆比较光滑。等效渗透系数的方向性受开口度分布形式的影响很大：当开口度-长度关联分布时,等效渗透系数各向异性;当开口度对数正态分布时,等效渗透系数各向同性;当开口度恒定分布时,等效渗透系数的特性介于二者之间。变化系数（CV）是否大于5%是判定岩体裂隙网络渗透系数是否具有方向性的判别标准。     

Closure of Fracture Due to Cover Stress Re-establishment After Coal Mining

In situ measurements of deformations, stresses, and closure of fractures, affecting water inflow following coal mining, are challenging due to the inaccessibility of fractured rock. In this paper, the authors studied the closure process of the fractured rock mass with the cover stress re-establishment based on a theoretical analysis and a scale model testing. A quantitative analysis is used to study the fracture distribution in the fractured zone. A function to describe a fracture aperture distribution in the fractured zone is proposed, which takes into account the curvature and thickness of the fractured rock. The theoretical analysis and a scale model testing both indicate that the cover stress re-establishment with mining distance increasing and the relationship between the fracture closure and cover stress re-establishment both satisfy a logarithmic function. The scale model test also shows the following features: (1) the fracture ratio (which is the fracture area divided by the total area of fracture and intact rock with a unit width in the vertical or horizontal direction) in the lower part of the fractured rock mass is greater than that in the upper part; (2) the initially fast decreased of fracture ratios is then followed by a slower decrease during the cover stress re-establishment process; (3) in the upper part of the rock mass, the vertical directional fractures with small apertures are being closed with cover stress re-establishment, which indicates an increase in the water resistance reducing the seepage from these parts of the fractured zone. This study improves the general understanding of the fracture closure process and cover stress re-establishment in the fractured rock mass after coal mining ceased, and provides a theoretical basis for water resource protection in case of underground coal mining.     

A shear‐dilation‐based model for evaluation of hydraulically stimulated naturally fractured reservoirs

The role of shear dilation as a mechanism of enhancing fluid flow permeability in naturally fractured reservoirs was mainly recognized in the context of hot dry rock (HDR) geothermal reservoir stimulation. Simplified models based on shear slippage only were developed and their applications to evaluate HDR geothermal reservoir stimulation were reported. Research attention is recently focused to adjust this stimulation mechanism for naturally fractured oil and gas reservoirs which reserve vast resources worldwide. This paper develops the overall framework and basic formulations of this stimulation model for oil and gas reservoirs. Major computational modules include: natural fracture simulation, response analysis of stimulated fractures, average permeability estimation for the stimulated reservoir and prediction of an average flow direction. Natural fractures are simulated stochastically by implementing ‘fractal dimension’ concept. Natural fracture propagation and shear displacements are formulated by following computationally efficient approximate approaches interrelating in situ stresses, natural fracture parameters and stimulation pressure developed by fluid injection inside fractures. The average permeability of the stimulated reservoir is formulated as a function of discretized gridblock permeabilities by applying cubic law of fluid flow. The average reservoir elongation, or the flow direction, is expressed as a function of reservoir aspect ratio induced by directional permeability contributions. The natural fracture simulation module is verified by comparing its results with observed microseismic clouds in actual naturally fractured reservoirs. Permeability enhancement and reservoir growth are characterized with respect to stimulation pressure, in situ stresses and natural fracture density applying the model to two example reservoirs. Copyright © 2002 John Wiley & Sons, Ltd.     

The role of hydromechanical coupling in fractured rock engineering

This paper provides a review of hydromechanical (HM) couplings in fractured rock, with special emphasis on HM interactions as a result of, or directly connected with human activities. In the early 1960s, the coupling between hydraulic and mechanical processes in fractured rock started to receive wide attention. A series of events including dam failures, landslides, and injection-induced earthquakes were believed to result from HM interaction. Moreover, the advent of the computer technology in the 1970s made possible the integration of nonlinear processes such as stress–permeability coupling and rock mass failure into coupled HM analysis. Coupled HM analysis is currently being applied to many geological engineering practices. One key parameter in such analyses is a good estimate of the relationship between stress and permeability. Based on available laboratory and field data, it was found that the permeability of fractured rock masses tends to be most sensitive to stress changes at shallow depth (low stress) and in areas of low in-situ permeability. In highly permeable, fractured rock sections, fluid flow may take place in clusters of connected fractures which are locked open as a result of previous shear dislocation or partial cementation of hard mineral filling. Such locked-open fractures tend to be relatively insensitive to stress and may therefore be conductive at great depths. Because of the great variability of HM properties in fractured rock, and the difficulties in using laboratory data for deriving in-situ material properties, the HM properties of fractured rock masses are best characterized in situ. Electronic Publication