剪胀对复杂裂隙岩体渗透性影响的离散元分析

裂隙岩体流固耦合问题是目前国内外研究热点之一,采用离散元软件UDEC对裂隙岩体发生节理剪胀的渗透性变化规律进行了模拟分析。基于现场调查的裂隙信息统计生成裂隙网络岩体模型。 通过固定垂直应力、不断增加应力比RS(RS=水平应力/垂直应力)使岩体出现剪胀,采用库伦滑移节理模式对岩体在剪胀过程中的渗透性变化情况进行模拟。结果发现:当应力比较小(RS3.1)时,节理水力隙宽、流速、渗透系数等参数都随着应力比的增加表现出明显的降低; 而当岩体出现剪胀现象之后(应力比大于3.1),发生剪切滑移和剪胀现象的节理控制着裂隙岩体的总体渗流行为,与不考虑节理剪胀的计算结果相比,岩体渗透能力出现了显著增长。这一结果表明,剪胀对裂隙岩体渗透性的影响是显著而不可忽视的。     

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.     

Stress-Dependent Fluid Flow and Permeability in Fractured Media: from Lab Experiments to Engineering Applications

Summary. Permeability is a physical property in rocks of extreme importance in energy engineering, civil and environmental engineering, and various areas of geology. Early on, fractures in fluid flow models were assumed to be rigid. However, experimental research and field data confirmed that stress-deformation behavior in fractures is a key factor governing their permeability tensor. Although extensive research was conducted in the past, the three-dimensional stress-permeability relationships, particularly in the inelastic deformation stage, still remain unclear. In this paper, laboratory experiments conducted on large concrete blocks with randomly distributed fractures and rock core samples are reported to investigate fluid flow and permeability variations under uniaxial, biaxial and triaxial complete stress-strain process. Experimental relationships among flowrate, permeability and fracture aperture in the fractured media are investigated. Results show that the flowrate and stress/aperture exhibit “cubic law” relationship for the randomly distributed fractures. A permeability-aperture relationship is proposed according to the experimental results. Based on this relationship, stress-dependent permeability in a set of fractures is derived in a three-dimensional domain by using a coupled stress and matrix-fracture interactive model. A double porosity finite element model is extended by incorporating such stress-dependent permeability effects. The proposed model is applied to examine permeability variations induced by stress redistributions for an inclined borehole excavated in a naturally fractured formation. The results indicate that permeability around underground openings depends strongly on stress changes and orientations of the natural fractures.     

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.     

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     

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 %.     

岩石裂隙的非饱和渗透特性及其演化规律研究

岩石裂隙的非饱和渗透特性是岩土、能源和环境等领域科学研究中的热点问题。采用三维激光扫描获取花岗岩裂隙的表面形貌特征,分析裂隙微观形貌特征对非饱和渗透特性的影响。研究在张拉、压缩、剪切等复杂荷载作用下裂隙开度分布的演化规律,建立复杂荷载作用下岩石裂隙非饱和毛细压力曲线演化模型。基于裂隙的微观形貌特征推导了岩石裂隙非饱和相对渗透系数模型,通过与试验数据对比,验证了模型的准确性和有效性,并在此基础上建立了复杂荷载作用下岩石裂隙非饱和相对渗透系数演化模型。研究成果对非饱和条件下裂隙岩体的水-力耦合机制研究具有一定指导意义。     

A simple failure criterion for rough joints and compound shear surfaces

The subject under investigation is the strength of a single shear plane which exhibits a regular, asymmetric roughness pattern. In the shear direction the asperities are so steeply inclined that the joint becomes mechanically non-effective with the result that the asperities are sheared off. Against the shear direction the asperities are only gently inclined. It is shown that this particular roughness pattern is of some general importance in geomechanics (examples: unconfined compression test; shear plane with secondary fractures).

Simple analytical considerations allow the formulation of a shear criterion, which is dependent on friction angle ø_m and cohesionc_m of the intact rock and on the inclination γ of the gently inclined parts of the asperities which are dipping against the shear direction. In the Mohr-diagram the criterion results in envelopes which converge at high normal stresses against the envelope of intact rock. Furthermore, the criterion expresses that both the slope of the envelopes and the dilation angle continuously decrease with increasing normal stress. Therefore the criterion adequately describes features which are regarded as most important when shearing rough joints or compound shear surfaces.     

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

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

粗糙裂隙水、气两相流相对渗透系数模型与数值分析

粗糙裂隙水、气两相流相对渗透系数是岩体工程多相渗流以及水力耦合分析的重要参数。从粗糙裂隙的细观结构出发，基于毛细吸持理论和立方定理，提出了粗糙裂隙水、气两相流相对渗透系数模型。通过与具有不同空间分布的粗糙裂隙水、气两相流试验数据对比分析，验证了模型的准确性。为进一步验证理论模型对不同粗糙程度裂隙的适用性，基于SRAM与Invasion Percolation模型，提出了粗糙裂隙的开度分布生成以及水、气两相流数值分析方法，计算结果表明理论模型与数值数据基本吻合一致，且优于X模型、V-C模型以及Corey模型。     

基于离散裂缝网络模型裂缝储层的随机模拟——以王窑区块为例

在燕山期、喜山期构造活动作用下,王窑区块发育北东向和南北向共轭剪切裂缝。针对该裂缝进行了统计分析,并利用离散裂缝网络模型对裂缝进行了随机模拟。模拟结果表明:研究区共轭裂缝北东向较为显著,且以高角度裂缝为主。裂缝对渗透率影响较为明显,将研究区的渗透率提高到20×10-3μm2,孔隙度增加约2%。     

Physical Modelling of Stress-dependent Permeability in Fractured Rocks

This paper presents the results of laboratory experiments conducted to study the impact of stress on fracture deformation and permeability of fractured rocks. The physical models (laboratory specimens) consisted of steel cubes simulating a rock mass containing three sets of orthogonal fractures. The laboratory specimens were subjected to two or three cycles of hydrostatic loading/unloading followed by the measurement of displacement and permeability. The results show a considerable difference in both deformation and permeability trends between the first loading and the subsequent loading/unloading cycles. However, the micrographs of the contact surfaces taken before and after the tests show that the standard deviation of asperity heights of measured surfaces are affected very little by the loadings. This implies that both deformation and permeability are rather controlled by the highest surface asperities which cannot be picked up by the conventional roughness characterization technique. We found that the dependence of flow rate on mechanical aperture follows a power law with the exponent n smaller or larger than three depending upon the loading stage. Initially, when the maximum height of the asperities is high, the exponent is slightly smaller than 3. The first loading, however, flattens these asperities. After that, the third loading and unloading yielded the exponent of around 4. Due to the roughness of contact surfaces, the flow route is no longer straight but tortuous resulting in flow length increase.     

A finite element porothermoelastic model for dual-porosity media

An existing dual-porosity finite element model has been extended to include thermo-hydro-mechanical coupling in both media. The model relies on overlapping distinct continua for the fluid and solid domains. In addition, conductive and convective heat transfers are incorporated using a single representative thermodynamics continuum. The model is applied to the problem of an inclined borehole drilled in a fractured formation subjected to a three-dimensional state of stress and, a temperature gradient between the drilling fluid and the formation. A sensitivity analysis has been carried out to study the impact of thermal loading, effect of heat transport by pore fluid flow and, the effect of parameters of the secondary medium used to represent the fractures. Copyright © 2004 John Wiley & Sons, Ltd.     

Discrete element modelling of material non‐coaxiality in simple shear flows

We investigate the quasi‐static simple shear flow of a two‐dimensional assembly of cohesionless particles using discrete element method (DEM) simulations. We focus on the unsteady flow regime where the solid would experience significant evolution of stresses, mobilised shear strength and dilation. We construct the DEM model using a discretised‐wall confined granular cell where the apparent boundary is allowed to dilate or contract synchronously with the confined solid. A rather uniform simple shear field is achieved across the whole assembly, which benefits rheological studies in generalising constitutive laws for continuum methods. We examine two aspects of the simple shear behaviour: macroscopic stress and strain rate evolution, particularly the non‐coaxiality between the principal directions of the two; and micromechanics such as evolution of fabric. For an initially anisotropic specimen sheared under constant normal pressure, the direction of principal stress rotates towards that of the principal strain rate, gradually reducing the degree of non‐coaxiality from about 45° to fluctuating around 0°. The rate in approaching coaxiality is slower in samples with larger initial porosity, stress ratio and mean stress. Generally, a faster rate in approaching coaxiality in simple shear is observed in a more dilatant sample, which often shows a larger degree of mobilised fabric anisotropy, suggesting the possible important role of instantaneous internal friction angle. The evolution of principal fabric direction resembles that of the principal stress direction. © 2013 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons, Ltd.     

Faulting, fracturing and in situ stress prediction in the Ahnet Basin, Algeria — a finite element approach

Many low-efficiency hydrocarbon reservoirs are productive largely because effective reservoir permeability is controlled by faults and natural fractures. Accurate and low-cost information on basic fault and fracture properties, orientation in particular, is critical in reducing well costs and increasing well recoveries. This paper describes how we used an advanced numerical modelling technique, the finite element method (FEM), to compute site-specific in situ stresses and rock deformation and to predict fracture attributes as a function of material properties, structural position and tectonic stress. Presented are the numerical results of two-dimensional, plane-strain end-member FEM models of a hydrocarbon-bearing fault-propagation-fold structure. Interpretation of the modelling results remains qualitative because of the intrinsic limitations of numerical modelling; however, it still allows comparisons with (the little available) geological and geophysical data.

In all models, the weak mechanical strength and flow properties of a thick shale layer (the main seal) leads to a decoupling of the structural deformation of the shallower sediments from the underlying sediments and basement, and results in flexural slip across the shale layer. All models predict rock fracturing to initiate at the surface and to expand with depth under increasing horizontal tectonic compression. The stress regime for the formation of new fractures changes from compressional to shear with depth. If pre-existing fractures exist, only (sub)horizontal fractures are predicted to open, thus defining the principal orientation of effective reservoir permeability. In models that do not include a blind thrust fault in the basement, flexural amplification of the initial fold structure generates additional fracturing in the crest of the anticline controlled by the material properties of the rocks. The folding-induced fracturing expands laterally along the stratigraphic boundaries under enhanced tectonic loading. Models incorporating a blind thrust fault correctly predict the formation of secondary syn- and anti-thetic mesoscale faults in the basement and sediments of the hanging wall. Some of these faults cut reservoir and/or seal layers, and thus may influence effective reservoir permeability and affect seal integrity. The predicted faults divide the sediments across the anticline in several compartments with different stress levels and different rock failure (and proximity to failure). These numerical model outcomes can assist classic interpretation of seismic and well bore data in search of fractured and overpressured hydrocarbon reservoirs.     

A multiscale DEM-LBM analysis on permeability evolutions inside a dilatant shear band

This paper presents a multiscale analysis of a dilatant shear band using a three-dimensional discrete element method and a lattice Boltzmann/finite element hybrid scheme. In particular, three-dimensional simple shear tests are conducted via the discrete element method. A spatial homogenization is performed to recover the macroscopic stress from the micro-mechanical force chains. The pore geometries of the shear band and host matrix are quantitatively evaluated through morphology analyses and lattice Boltzmann/finite element flow simulations. Results from the discrete element simulations imply that grain sliding and rotation occur predominately with the shear band. These granular motions lead to dilation of pore space inside the shear band and increases in local permeability. While considerable anisotropy in the contact fabric is observed with the shear band, anisotropy of the permeability is, at most, modest in the assemblies composed of spherical grains.     

Permeability evolution in fractured coal: The roles of fracture geometry and water-content

We report laboratory experiments that investigate the permeability evolution of an anthracite coal as a function of applied stress and pore pressure at room temperature as an analog to other coal types. Experiments are conducted on 2.5 cm diameter, 2.5-5 cm long cylindrical samples at confining stresses of 6 to 12 MPa. Permeability and sorption characteristics are measured by pulse transient methods, together with axial and volumetric strains for both inert (helium (He)) and strongly adsorbing (methane (CH₄) and carbon dioxide (CO₂)) gases. To explore the interaction of swelling and fracture geometry we measure the evolution of mechanical and transport characteristics for three separate geometries — sample A containing multiple small embedded fractures, sample B containing a single longitudinal through-going fracture and sample C containing a single radial through-going fracture. Experiments are conducted at constant total stress and with varied pore pressure — increases in pore pressure represent concomitant (but not necessarily equivalent) decreases in effective stress. For the samples with embedded fractures (A and C) the permeability first decreases with an increase in pressure (due to swelling and fracture constraint) and then increases near-linearly (due to the over-riding influence of effective stresses). Conversely, this turnaround in permeability from decreasing to increasing with increasing pore pressure is absent in the discretely fractured sample (B) — the influence of the constraint of the connecting fracture bridges in limiting fracture deformation is importantly absent as supported by theoretical considerations. Under water saturated conditions, the initial permeabilities to all gases are nearly two orders of magnitude lower than for dry coal and permeabilities increase with increasing pore pressure for all samples and at all gas pressures. We also find that the sorption capacities and swelling strains are significantly reduced for water saturated samples — maybe identifying the lack of swelling as the primary reason for the lack of permeability decrease. Finally, we report the weakening effects of gas sorption on the strength of coal samples by loading the cores to failure. Results surprisingly show that the strength of the intact coal (sample A) is smaller than that of the axially fractured coal (sample B) due to the extended duration of exposure to CH₄ and CO₂. Average post-failure particle size for the weakest intact sample (A) is found to be three times larger than that of the sample B, based on the sieve analyses from the samples after failure. We observe that fracture network geometry and saturation state exert important influences on the permeability evolution and strength of coal under in situ conditions.     

Modelling of steady‐state fluid flow in 3D fractured isotropic porous media: application to effective permeability calculation

In this paper, 3D steady‐state fluid flow in a porous medium with a large number of intersecting fractures is derived numerically by using collocation method. Fluid flow in the matrix and fractures is described by Darcy's law and Poiseuille's law, respectively. The recent theoretical development presented a general potential solution to model the steady‐state flow in fractured porous media under a far‐field condition. This solution is a hypersingular integral equation with pressure field in the fracture surfaces as the main unknown. The numerical procedure can resolve the problem for any form of fractures and also takes into account the interactions and the intersection between fractures. Once the pressure field and then the flux field in fractures have been determined, the pressure field in the porous matrix is computed completely. The basic problem of a single fracture is investigated, and a semi‐analytical solution is presented. Using the solution obtained for a single fracture, Mori‐Tanaka and self‐consistent schemes are employed for upscaling the effective permeability of 3D fractured porous media. Copyright © 2012 John Wiley & Sons, Ltd.     

Changes in permeability of sand during triaxial loading: effect of fine particles production

Various mechanisms can affect the permeability of dense unconsolidated sands: Volumetric dilation can lead to permeability increase, whereas strain localization in shear bands may increase or decrease the permeability depending on the state of compaction and on the level of grains breakage inside the band. To investigate these various mechanisms, an experimental study has been performed to explore the effect of different factors such as grain size and grain shape, confining pressure, level of shear, stress path, and formation of one or several shear bands on the permeability of dense sands under triaxial loading. The experimental results show a reduction of permeability during the consolidation phase and during the volumetric contraction phase of shear loading, which can be related to the decrease of porosity. The experimental results also show that, depending on the confining pressure, the permeability remains stable or decreases during the volumetric dilation phase despite the increase of total porosity. This permeability reduction is attributed to the presence of fine particles, which result from grains attrition during pre-localization and grains breakage inside the shear band during the post-localization phase.