Geotechnical Engineering Online Library

It is well known that shear wave propagates slower across than parallel to a fracture, and as a result, a travelling shear wave splits into two directions when it encounters a fracture. Shear wave splitting and permeability of porous rock core samples having single fracture were experimentally investigated using a high-pressure triaxial cell, which can measure seismic shear wave velocities in two directions mutually perpendicular to the sample axis in addition to the longitudinal compressive wave velocity. A single fracture was created in the samples using a modified Brazilian split test device, where the cylindrical sample edges were loaded on two diametrically opposite lines by sharp guillotines along the sample length. Based on tilt tests and fracture surface profilometry, the method of artificially induced tensile fracture in the sample was found to create repeatable fracture surfaces and morphologies. Seismic velocities of the fractured samples were determined under different levels of stress confinement and fracture shear displacement or mismatch. The effective confining stress was varied from 0.5 MPa to 55 MPa, while the fractures were mismatched by 0 mm, 0.45 mm and 1 mm. The degree of matching of the fracture surfaces in the core samples was evaluated using the joint matching coefficient (JMC). Shear wave splitting, as measured by the difference in the magnitudes of shear wave velocities parallel (VS1) and perpendicular (VS2) to the fracture, is found to be insensitive to the degree of mismatching of the fracture joint surfaces at 2 MPa, and decreased and approached zero as the effective stress was increased. Simple models for the stress- and JMC-dependent shear wave splitting and fractured rock permeability were developed based on the experimental observations. The effects of the joint wall compressive strength (JCS), JMC and stress on the stress dependency of joint aperture were discussed in terms of hydro-mechanical response. Finally, a useful relationship between fractured rock permeability and shear wave splitting was found after normalization by using JMC.

Explicit solution techniques have been widely used in geotechnical engineering for simulating the coupled hydro-mechanical H-) interaction of fluid flow and deformation induced by structures built above and under saturated ground, i.e. circular footing and deep tunnel. However, the technique is only conditionally stable and requires small time steps, portending its inefficiency for simulating large-scale H-M problems. To improve its efficiency, the unconditionally stable alternating direction explicit (ADE) scheme could be used to solve the flow problem. The standard ADE scheme, however, is only moderately accurate and is restricted to uniform grids and plane strain flow conditions. This paper aims to remove these drawbacks by developing a novel high-order ADE scheme capable of solving flow problems in non-uniform grids and under axisymmetric conditions. The new scheme is derived by performing a fourth-order finite difference (FD) approximation to the spatial derivatives of the axisymmetric fluid–diffusion equation in a non-uniform grid configuration. The implicit Crank-Nicolson technique is then applied to the resulting approximation, and the subsequent equation is split into two alternating direction sweeps, giving rise to a new axisymmetric ADE scheme. The pore pressure solutions from the new scheme are then sequentially coupled with an existing geomechanical simulator in the computer code fast Lagrangian analysis of continua (FLAC). This coupling procedure is called the sequentially-explicit coupling technique based on the fourth-order axisymmetric ADE scheme or SEA-4-AXI. Application of SEA-4-AXI for solving axisymmetric consolidation of a circular footing and of advancing tunnel in deep saturated ground shows that SEA-4-AXI reduces computer runtime up to 42%–50% that of FLAC's basic scheme without numerical instability. In addition, it produces high numerical accuracy of the H-M solutions with average percentage difference of only 0.5%–1.8%.

Understanding microcracking near coalesced fracture generation is critically important for hydrocarbon and geothermal reservoir characterization as well as damage evaluation in civil engineering structures. Dense and sometimes random microcracking near coalesced fracture formation alters the mechanical properties of the nearby virgin material. Individual microcrack characterization is also significant in quantifying the material changes near the fracture faces (i.e. damage). Acoustic emission (AE) monitoring and analysis provide unique information regarding the microcracking process temporally, and information concerning the source characterization of individual microcracks can be extracted. In this context, laboratory hydraulic fracture tests were carried out while monitoring the AEs from several piezoelectric transducers. In-depth post-processing of the AE event data was performed for the purpose of understanding the individual source mechanisms. Several source characterization techniques including moment tensor inversion, event parametric analysis, and volumetric deformation analysis were adopted. Post-test fracture characterization through coring, slicing and micro-computed tomographic imaging was performed to determine the coalesced fracture location and structure. Distinct differences in fracture characteristics were found spatially in relation to the openhole injection interval. Individual microcrack AE analysis showed substantial energy reduction emanating spatially from the injection interval. It was quantitatively observed that the recorded AE signals provided sufficient information to generalize the damage radiating spatially away from the injection wellbore.

Experiments on rock joint behaviors have shown that joint surface roughness is mobilized under shearing, inducing dilation and resulting in nonlinear joint shear strength and shear stress vs. shear displacement behaviors. The Barton–Bandis (BB) joint model provides the most realistic prediction for the nonlinear shear behavior of rock joints. The BB model accounts for asperity roughness and strength through the joint roughness coefficient (JRC) and joint wall compressive strength (JCS) parameters. Nevertheless, many computer codes for rock engineering analysis still use the constant shear strength parameters from the linear Mohr–Coulomb (M−C) model, which is only appropriate for smooth and non-dilatant joints. This limitation prevents fractured rock models from capturing the nonlinearity of joint shear behavior. To bridge the BB and the M−C models, this paper aims to provide a linearized implementation of the BB model using a tangential technique to obtain the equivalent M−C parameters that can satisfy the nonlinear shear behavior of rock joints. These equivalent parameters, namely the equivalent peak cohesion, friction angle, and dilation angle, are then converted into their mobilized forms to account for the mobilization and degradation of JRC under shearing. The conversion is done by expressing JRC in the equivalent peak parameters as functions of joint shear displacement using proposed hyperbolic and logarithmic functions at the pre- and post-peak regions of shear displacement, respectively. Likewise, the pre- and post-peak joint shear stiffnesses are derived so that a complete shear stress-shear displacement relationship can be established. Verifications of the linearized implementation of the BB model show that the shear stress-shear displacement curves, the dilation behavior, and the shear strength envelopes of rock joints are consistent with available experimental and numerical results.

This paper presents an improved understanding of coupled hydro-thermo-mechanical HTM) hydraulic fracturing of quasi-brittle rok using the bonded particle model (BPM) within the discrete element method (DEM). BPM has been recently extended by the authors to account for coupled convective–conductive heat flow and transport, and to enable full hydro-thermal fluid–solid coupled modeling. The application of the work is on enhanced geothermal systems (EGSs), and hydraulic fracturing of hot dry rock HDR) is studied in terms of the impact of temperature difference between rock anda flowing fracturing fluid. Micro-mechanical investigation of temperature and fracturing fluid effects on hydraulic fracturing damage in rocks is presented. It was found that fracture is shorter with pronounced secondary microcracking along the main fracture for the case when the convective–conductive thermal heat exchange is considered. First, the convection heat exchange during low-viscosity fluid infiltration in permeable rock around the wellbore causes significant rock cooling, where a finger-like fluid infiltration was observed. Second, fluid infiltration inhibits pressure rise during pumping and delays fracture initiation and propagation. Additionally, thermal damage occurs in the whole area around the wellbore due to rock cooling and cold fluid infiltration. The size of a damaged area around the wellbore increases with decreasing fluid dynamic viscosity. Fluid and rock compressibility ratio was found to have significant effect on the fracture propagation velocity.

Fracture systems have strong influence on the overall mechanical behavior of fractured rock masses due to their relatively lower stiffness and shear strength than those of the rock matrix. Understanding the effects of fracture geometrical distribution, such as length, spacing, persistence and orientation, is important for quantifying the mechanical behavior of fractured rock masses. The relation between fracture geometry and the mechanical characteristics of the fractured rock mass is complicated due to the fact that the fracture geometry and mechanical behaviors of fractured rock mass are strongly dependent on the length scale. In this paper, a comprehensive study was conducted to determine the effects of fracture distribution on the equivalent continuum elastic compliance of fractured rock masses over a wide range of fracture lengths. To account for the stochastic nature of fracture distributions, three different simulation techniques involving Oda's elastic compliance tensor, Monte Carlo simulation (MCS), and suitable probability density functions (PDFs) were employed to represent the elastic compliance of fractured rock masses. To yield geologically realistic results, parameters for defining fracture distributions were obtained from different geological fields. The influence of the key fracture parameters and their relations to the overall elastic behavior of the fractured rock mass were studied and discussed. A detailed study was also carried out to investigate the validity of the use of a representative element volume (REV) in the equivalent continuum representation of fractured rock masses. A criterion was also proposed to determine the appropriate REV given the fracture distribution of the rock mass.