P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Dynamic Model of Fracture Normal Behaviour and Application to Prediction of Stress Wave Attenuation Across Fractures

Summary.  The purpose of this paper is to establish a dynamic constitutive model of fracture normal behaviour, based on laboratory tests of artificial fractures cast by cement mortar. A series of tests are systematically carried out under quasi-static (10⁻¹ MPa/s) up to highly dynamic (10³ MPa/s) monotonic loading conditions. The normal stress-fracture closure response is measured at different loading rates. Based on the measured curves, a nonlinear (hyperbolic) dynamic model of fracture normal behaviour, termed as dynamic BB model, is proposed. The dynamic model is modified from the existing BB model of static normal behaviour of fractures by taking into account the loading-rate effect. Two important dynamic parameters of fractures, FSC _d (dynamic fracture stiffness constant, which describes the incremental ratio of dynamic initial stiffness) and FCC _d (dynamic fracture closure constant, which describes the decremental ratio of dynamic maximum allowable closure), are identified. They indicate the quantitative degree of loading-rate effect on fracture normal behaviour subjected to dynamic loads. For practical application, the new model is incorporated into the Universal Distinct Element Code (UDEC) and subsequently, UDEC modelling of normally incident P-wave transmission across single fractures with the dynamic BB model is conducted. Wave transmission coefficient is obtained for various combinations of fracture dynamic parameters, as well as different wave amplitudes and frequencies. The numerical results show that wave transmission coefficient for a fracture with the dynamic BB model is greater than that for a fracture with the static BB model. In addition, a fracture with higher values of FSC _d and FCC _d leads to higher transmission (lower attenuation). Author’s address: J. Zhao, Ecole Polytechnique Federale de Lausanne (EPFL), Rock Mechanics Laboratory, 1015 Lausanne, Switzerland     

Validating the ability of the discontinuous deformation analysis method to model normal P‐wave propagation across rock fractures

The effects of fractures on wave propagation problems are increasingly abstracting the attention of scholars and engineers in rock engineering field. This study aims to fully validate the ability of discontinuous deformation analysis (DDA) to model normal P‐wave propagation across rock fractures. The effects of a single fracture and multiple parallel fractures are all tested. The results indicate that DDA can accurately reflect the fracture effects, including the fractures stiffness, the fracture spacing and the fracture number, and the effects of incident wave frequency on one‐dimensional P‐wave propagation problems. Thus, DDA is able to deal well with normal incident P‐wave propagation problems. Copyright © 2017 John Wiley & Sons, Ltd.     

Displacement and stress distributions under a uniform inclined rectangular load on a cross‐anisotropic geomaterial

In practical engineering, an applied rectangular area load is not often horizontally or vertically distributed but is frequently inclined at a certain angle with respect to the horizontal and vertical axes. Thus, the solutions of displacements and stresses due to such a load are essential to the design of foundations. This article yields the analytical solutions of displacements and stresses subjected to a uniform rectangular load that inclines with respect to the horizontal and vertical axes, resting on the surface of a cross‐anisotropic geomaterial. The planes of cross‐anisotropy are assumed to be parallel to the horizontal ground surface. The procedures to derive the solutions can be integrated the modified point load solutions, which are represented by several displacement and stresses elementary functions. Then, upon integrations, the displacement and stress integral functions resulting from a uniform inclined rectangular load for (1) the displacements at any depth, (2) the surface displacements, (3) the average displacements in a given layer, (4) the stresses at any depth, and (5) the average stresses in a given layer are yielded. The proposed solutions are clear and concise, and they can be employed to construct a series of calculation charts. In addition, the present solutions clarify the load inclinations, the dimensions of a loaded rectangle, and the analyzed depths, and the type and degree of geomaterial anisotropy profoundly affect the displacements and stresses in a cross‐anisotropic medium. Parametric results show that the load inclination factor should be considered when an inclined rectangular load uniformly distributed on the cross‐anisotropic material. Copyright © 2008 John Wiley & Sons, Ltd.     

Transmission of Elastic P-waves across Single Fractures with a Nonlinear Normal Deformational Behavior

Summary This paper presents a theoretical study on normally incident elastic P-wave transmission across single dry fractures with a nonlinear normal deformational behavior. The effects of nonlinear fracture normal behavior on P-wave transmission are examined without the mixture of fracture shear behavior. The linear displacement discontinuity model for wave propagation across fractures is extended to a nonlinear model – the hyperbolic elastic model (BB model). Numeric solutions of magnitudes of transmission (|T _non|) and reflection (|R _non|) coefficients, for normally incident P-wave transmission across the nonlinear deformable fractures, are obtained and related to the closure behavior of fractures. Parametric studies are conducted to acquire an insight into the effects of the nonlinear fracture normal deformation on P-wave transmission, in terms of initial normal stiffness and the ratio of current maximum closure to maximum allowable closure of the fractures, as well as the incident wave amplitude and frequency. Comparisons between the linear and nonlinear models are presented. It is shown that, |T _lin| and |R _lin| for the linear model are special solutions of |T _non| and |R _non| for the nonlinear model, when the incident wave amplitude is so low that the current maximum closure of fracture incurred during the wave transmission is much smaller, relative to the maximum allowable closure. In addition, the nonlinear fracture behavior gives rise to a phenomenon of higher harmonics during the wave transmission across the fracture. The higher harmonics contribute to the increase of |T _non| from |T _lin|.     

Analysis of a deformable fracture in permeable material

The response of deformable fractures to changes in fluid pressure controls phenomena ranging from the flow of fluids near wells to the propagation of hydraulic fractures. We developed an analysis designed to simulate fluid flows in the vicinity of asperity‐supported fractures at rest, or fully open fractures that might be propagating. Transitions between at‐rest and propagating fractures can also be simulated. This is accomplished by defining contact aperture as the aperture when asperities on a closing fracture first make contact. Locations on a fracture where the aperture is less than the contact aperture are loaded by both fluid pressure and effective stress, whereas locations where the aperture exceeds the contact aperture are loaded only by fluid pressure. Fluid pressure and effective stress on the fracture are determined as functions of time by solving equations of continuity in the fracture and matrix, and by matching the global displacements of the fracture walls to the local deformation of asperities. The resulting analysis is implemented in a numerical code that can simulate well tests or hydraulic fracturing operations. Aperture changes during hydraulic well tests can be measured in the field, and the results predicted using this analysis are similar to field observations. The hydraulic fracturing process can be simulated from the inflation of a pre‐existing crack, to the propagation of a fracture, and the closure of the fracture to rest on asperities or proppant. Two‐dimensional, multi‐phase fluid flow in the matrix is included to provide details that are obscured by simplifications of the leakoff process (Carter‐type assumptions) used in many hydraulic fracture models. Execution times are relatively short, so it is practical to implement this code with parameter estimation algorithms to facilitate interpretation of field data. Copyright © 2006 John Wiley & Sons, Ltd.     

Manuel Rocha Medal Recipient Wave interaction with underground openings in fractured rock

Summary The objective of this work is to lead to improved models of seismic wave propagation around underground openings by studying the interaction of the waves with the fractured rock surrounding these openings. It demonstrates that seismic models can help in stability problems such as rockbursting in deep-level mining, or in the interpretation of micro-fracturing at waste storage sites. A significant emphasis is placed on comparing the models with observations from controlled experiments. These comparisons demonstrate that the wave propagation can be reliably and accurately modelled, and in so doing it motivates their application to the larger rock engineering problems. Seismic wave models are first applied to laboratory experiments on multiple fractures. Simulation through multiple displacement discontinuities yields strikingly similar waveforms to the experiments, while also identifying the need to build stress dependence into the fracture models, such as stress dependent fracture stiffness. The wave-fracture modelling is extended to in situ fractures in rock at the surface of a deep tunnel, using data collected during an acoustic emission experiment at the URL Mine-by tunnel. Waveforms from the velocity scans are compared against those from elastic models and various models of fracture, such as random assemblies of small open fractures (cracks) and larger fractures with fracture stiffness. Results indicate that it is possible to account for the wave-speeds and amplitudes using models with fractures. A generic method is then proposed for calculating the frequency variation of wave-speed and amplitude for any collection of cracks. The models of fracture are then applied to the rockburst problem, to investigate how the excavation affects the amplitude and the distribution of ground motion. The results provide important insights into the causes of the apparent amplification observed by researchers in this field. The thesis also covers the theory of the models used, including novel numerical work on dispersion and new grid schemes. The full detail of the work cannot be covered in this paper which instead seeks to summarize the main achievements.     

层状场地中局部不均体对平面P波的散射

利用间接边界元方法,在频域内求解了层状场地中局部不均体对平面P波的散射。利用精确的土层动力刚度矩阵进行自由场反应分析,求得位移和应力响应。通过计算虚拟分布荷载的格林影响函数,求得相应位移和应力;根据边界条件确定虚拟分布荷载,最终得到问题的解答。研究了入射P波时,不均体宽度、埋深、厚度、入射角、入射频率度等参数对地表位移幅值的影响,并与相应自由场的结果进行了比较。不均体对P波散射有重要影响,在工程场地地震安全性评价中,应合理考虑局部不均体对场地设计地震动参数确定的影响。     

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.     

One‐dimensional analytical solution for mesoscopic flow induced damping in a double‐porosity dual‐permeability material

Heterogeneities, such as fractures and cracks, are ubiquitous in porous rocks. Mesoscopic heterogeneities, that is, heterogeneities on length scales much larger than typical pore size but much smaller than the wavelength, are increasingly believed to be responsible for significant wave energy loss in the seismic frequency band. When a compressional wave stresses a material containing mesoscopic heterogeneities, the more compliant parts of the material (e.g., fractures and cracks) respond with a greater fluid pressure than the stiffer portions (e.g., matrix pores). The induced fluid flow, resulting from the pressure gradients developed on such scale, is called mesoscopic flow. In the present study, the double‐porosity dual‐permeability model is adopted to incorporate mesoscopic heterogeneities into rock models to account for the attenuation of wave energy. Based on the model, the damping effect due to mesoscopic flow in a one‐dimensional porous structure is investigated. Analytical solutions for several boundary‐value problems are obtained in the frequency domain. The dynamic responses of infinite and finite porous layer are examined. Numerical calculations show that the damping effect of mesoscopic flow is significant on the pore pressure response and the resulting effective stress. For the displacement, the effect is seen only at the very low frequency range or near the resonance frequencies. Copyright © 2017 John Wiley & Sons, Ltd.     

Effects of fracture lineaments and in-situ rock stresses on groundwater flow in hard rocks: a case study from Sunnfjord, western Norway

Systematic field mapping of fracture lineaments observed on aerial photographs shows that almost all of these structures are positively correlated with zones of high macroscopic and mesoscopic fracture frequencies compared with the surroundings. The lineaments are subdivided into zones with different characteristics: (1) a central zone with fault rocks, high fracture frequency and connectivity but commonly with mineral sealed fractures, and (2) a damage zone divided into a proximal zone with a high fracture frequency of lineament parallel, non-mineralized and interconnected fractures, grading into a distal zone with lower fracture frequencies and which is transitional to the surrounding areas with general background fracturing. To examine the possible relations between lineament architecture and in-situ rock stress on groundwater flow, the geological fieldwork was followed up by in-situ stress measurements and test boreholes at selected sites. Geophysical well logging added valuable information about fracture distribution and fracture flow at depths. Based on the studies of in-situ stresses as well as the lineaments and associated fracture systems presented above, two working hypotheses for groundwater flow were formulated: (i) In areas with a general background fracturing and in the distal zone of lineaments, groundwater flow will mainly occur along fractures parallel with the largest in-situ rock stress, unless fractures are critically loaded or reactivated as shear fractures at angles around 30° to σ_H; (ii) In the influence area of lineaments, the largest potential for groundwater abstraction is in the proximal zone, where there is a high fracture frequency and connectivity with negligible fracture fillings. The testing of the two hypotheses does not give a clear and unequivocal answer in support of the two assumptions about groundwater flow in the study area. But most of the observed data are in agreement with the predictions from the models, and can be explained by the action of the present stress field on pre-existing fractures.     

Static deformation of an orthotropic multilayered elastic half-space by two-dimensional surface loads

The transfer matrix approach is used to solve the problem of static deformation of an orthotropic multilayered elastic half-space by two-dimensional surface loads. The general problem is decoupled into two independent problems. The antiplane strain problem and the plane strain problem are considered in detail. Integral expressions for displacements and stresses at any point of the medium due to a normal line load and a shear line load, acting parallel to a symmetry axis, are obtained. In the case of a uniform half-space, closed form analytic expressions for displacements and stresses are derived. The procedure developed is quite easy and convenient for numerical computations.     

On the mechanical interaction between a fluid-filled fracture and the earth's surface

The mechanical interaction between a fluid-filled fracture (e.g., hydraulic fracture joint, or igneous dike) and the earth's surface is analyzed using a two-dimensional elastic solution for a slit of arbitrary inclination buried beneath a horizontal free surface and subjected to an arbitrary pressure distribution. The solution is obtained by iteratively superimposing two fundamental sets of analytical solutions.For uniform internal pressure the slit behaves essentially as if it were in an infinite region if the depth-to-center is three times greater than the half-length. For shallower slits interaction with the free surface is pronounced: stresses and displacements near the slit differ by more than 10% from values for the deeply buried slit. The following changes are noted as the depth-to-center decreases:

1. (1) the mode I stress intensity factor increases for both ends of the slit, but more rapidly at the upper end;

2. (2) the mode II stress-intensity factor is significantly different from zero (except for vertical slits) suggesting propagation out of the original plane of the slit;

3. (3) displacements of the slit wall are asymmetric such that the slit gaps open more widely near the upper end. Similar changes are noted if fluid density creates a linear pressure gradient that is smaller than the lithostatic gradient. Under such conditions natural fractures should propagate preferentially upward toward the earth's surface requiring less pressure as they grow in length.

If deformation near the surface is of interest, the model should account explicitly for the free surface. Stresses and displacements at the free surface are not approximated very well by values calculated along a line in an infinite region, even when the slit is far from the line. As depth-to-center of a shallow pressurized slit decreases, the following changes are noted:

1. (1) displacements of the free surface increase to the same order of magnitude as the displacements of the slit walls,

2. (2) tensile stresses of magnitude greater than the pressure in the slit are concentrated along the free surface. The relative surface displacements over a shallow vertical slit are downward over the slit and upward to both sides of this area. The tensile stress acting parallel to the free surface over a shallow vertical slit is concentrated in two maxima adjacent to a point of very low stress immediately over the slit.

The solution is used to estimate the length-to-depth ratio at which igneous sills have gained sufficient leverage on overlying strata to bend these strata upward and form a laccolith. The pronounced mode II stress intensity associated with shallow horizontal slits explains the tendency for some sills to climb to higher stratigraphie horizons as they grow in length. The bimodal tensile stress concentration over shallow vertical slits correlates qualitatively with the distribution of cracks and normal faults which flank fissure eruptions on volcanoes. The solution may be used to analyze surface displacements and tilts over massive hydraulic fractures in oil fields and to understand the behavior of hydraulic fractures in granite quarries.     

Origin of sheet structure, 1. Morphology and boundary conditions

Sheet structure, or large-scale exfoliation, is the division of a rock mass into lenses, plates or “sheets” approximately parallel to the earth's surface. Sheet fractures, which separate the plates, are characterized by surface markings resembling those formed during the brittle fracture of metals, glass and ceramics and those on joints in rock. Lineations, or hackle marks, parallel to the direction of fracture propagation, are common. Branching occurs during propagation so that the fractures appear as many echelon segments in profile. The opposing surfaces of sheet fractures observed in quarry walls and natural exposures are typically in contact. Damage, if any, to rock adjacent to sheet fractures is generally limited to a zone less than one centimeter wide. Sheets tend to parallel preferred orientations of microcracks in many rock masses. Most sheets that are exposed today are of prehistoric origin, but in some rock masses, such as Chelmsford granite in Massachusetts, new sheets continue to form.

Sheet structure forms in environments of high differential stress, dominated by large-scale compression parallel to an exposed rock surface. In parts of the Chelmsford granite the magnitude of this compression is greater than 30 MPa in the zone of sheeting. High differential stresses near the ground surface can result from several natural agents, including contemporary tectonic forces, vertical unloading of a rock mass that formed at depth under high triaxial compression, and suppression of expansion that would otherwise result from temperature increases or chemical alteration of the rock. Exfoliation of boulders during forest fires and spheroidal weathering of boulders appears to occur under similar states of stress, although the agents responsible for the stresses differ.     

黏弹性节理岩体的位移不连续模型研究

选用三参数标准线性固体作为岩石本构,提出了一种考虑岩体黏弹性的位移不连续模型;根据一维黏弹性波的特征线法,推导了节理处质点速度、应力和应变递推公式。首先,基于分离式霍普金森压杆（SHPB）对砂层进行试验,得到其应力-应变关系,并换算出砂层节理的法向刚度;接着,通过一维强间断黏弹性波的波速公式、高频波衰减系数以及任一频率下的衰减系数,确定出数值算法中的三参数。最后,基于自制的摆锤装置,探讨了一维应力波在节理岩体中的传播规律,试验中以两根长1 000 mm、直径为68.50 m的岩杆作为入射和透射杆,以3 mm砂层模拟节理。试验和数值结果吻合度良好,进一步验证了该方法的可靠性。     

Analytical solution of the problem for a set of shear fractures with Coulomb friction

A simplified solution of a two-dimensional problem on a set of interaction shear fractures with Coulomb friction along their wings is proposed. In the classical statement of this problem, the equality between shear stresses and friction stresses must necessarily hold in a new equilibrium state at each point along the fracture plane. The solution of the problem is reduced to the solution of a set of singular integral equations with respect to unknown functions of a shear jump on fractures. In reality, our knowledge of fracture conditions is rather approximate. In addition, in the problems of seismology and tectonophysics it is sufficient to approximately estimate dynamic (taking into account the inertia forces) and static perturbations from a fracture in the form of the first terms of a true solution series. With this aim in view, point models of an earthquake source or continual representation of discontinuous deformations are used. Within the framework of the requirements of such approaches, it is assumed that the condition of equality between the sum of shear stresses and friction stresses on the shear fractures in a new equilibrium state is met. In addition, the complex potential function is calculated for each fracture based on the function of the jump of its wings obtained in the problem for a solitary fracture. Such statement of the problem of a set of neighboring and even intersecting fractures makes it possible to reduce its solution to a set of linear algebraic equations with respect to shear stresses relieved at each fracture and the average along its length, they being the unknown values.     

Displacement efficiency of faults and fractures

Faults and fractures evolve kinematically as Somigliana dislocations. As a result, displacement gradients are always present, with displacements reducing to zero at fault/fracture tips. We present a graphical method for the comparison of displacement gradients from fractures of various magnitudes and types. Displacement efficiency of a fault/fracture is defined as its ability to maximise displacement across its surface.     

Multi-Scattering of Plane P Waves by Several Fluid-Filled Cavities

The conformal mapping method of complex functions and expansion method of wave functions are adopted, the stresses and displacements at the boundary between the cavities and inner fluids are considered as continuous, and then the theoretical solutions are obtained about the scattering problems of incident P waves by several fluid-filled cavities. Two circular cavities are taken as examples, the influences of frequencies of incident P waves, whether there are inner fluids or not and the spacing distances between the two cavities on the scattering amplitudes are studied, and the results reveal that: the scattering amplitudes of two cavities decrease when the distance between the two cavities increase and they reach constants when the spacing distance exceeds 7.5 times of cavities radius, and the inner fluids would absorb part of energy of the scattered waves and the scattering amplitude of the cavities filled with fluids would slightly smaller than those without fluids.     

Theory of chocolate tablet boudinage

Two-dimensional boudinage in a flattening type of bulk deformation, with equal layer-parallel extensions in all directions, leads to the development of roundish or polygonal outlines of boudins in plan-view. As combined experimental and theoretical studies show, chocolate tablet boudinage with two sets of mutually perpendicular boudin axes may form in different ways. (1) Unequal layer-parallel extension in the matrix results in one set of extension fractures forming perpendicular to the greatest principal stress in the matrix. Once these long narrow boudins are formed, the greatest principal stress in the brittle layer becomes approximately parallel to the long axis of the boudin. As a result a second set of fractures forms normal to the first set. (2) In lineated rocks the anisotropy of tensile strength leads to the sequential formation of two sets of extension fractures, parallel and perpendicular to the lineation. Depending on the orientation of the lineation the boudin axes may or may not be parallel to the principal stresses in the matrix. (3) Boudins with rectangular plan-view may also form when two successive events of unidirectional boudinage are superposed on one another. Irrespective of the direction of principal extensional strain rate in the matrix, the second generation extension fractures are likely to form approximately perpendicular to the first generation boudin axes.     

基于qp 波及qs 波衰减系数检测裂缝的方法

裂缝储层预测是当前国内外油气勘探中的重点和难点,而裂缝储层预测最关键的问题是裂缝发育带的检测。基于qp 波、qsh 波、qsv 波的衰减系数对裂缝发育带的检测,利用Hudson 理论模型经数值计算,分别研究了裂缝密度、入射角及频率对qp 波、qsh 波、qsv 波衰减系数的影响。理论模型试算表明,频率和裂缝密度是影响衰减系数的主要参数,而裂缝中充填物不同时,衰减系数也表现出很大的差异,这为利用qp 波、qsh 波、qsv 波衰减属性检测裂缝发育带提供了理论依据。