It is of great importance to determine peak shear strength (PSS) of rock fractures, and data-driven criteria have showed advances in fitting capability in recent years. However, the generalization ability of existing data-driven criteria is limited by dataset size and fracture roughness characterization, which is negative to predictive power and robustness of models. Here we proposed a novel data-driven criterion to predict PSS of rock fractures, with high generalization ability on real experimental data. We first created large-scale low-fidelity dataset by discrete-element modeling, and small-scale high-fidelity dataset by laboratory direct shear tests. The numeric features include normal stress, mechanical properties (including PSS of intact and flat-fracture rock specimens), secondary properties (including internal friction angle, cohesion strength and basic friction angle), and the matrixed feature is topography data. We then established domain adaptation (DA) models for cross-domain knowledge transfer between the low- and high-fidelity datasets, and roughness features were automatically extracted by convolution kernels. The best DA-based model is weighting adversarial neural network, outranking other models by error indicator, and the average relative error on experimental data of new rock types is within 10.0 %. Finally, the sensitivity of input features is investigated, which further proves the promising potential of the developed data-driven PSS criterion of rock fractures in engineering practice.
Rock failure under external force is a process of energy conversion between the external environment and the rock system. This study aims to quantify rock damage and predict failure from an energy perspective. Infrared radiation (IR) and acoustic emission (AE) technologies were used to monitor the failure process of red sandstone during uniaxial loading experiments in real time. The energy evolution law during the rock failure process was analyzed. Based on the Stefan–Boltzmann law, a quantitative parameter, average cumulative radiation energy increment (), was proposed for IR indicators. A coupling mathematical model between elastic strain energy and was derived. The correlation between cumulative AE energy and dissipated strain energy was also analyzed. Results reveal that the rock failure process can be divided into four stages according to energy evolution: compaction, elastic, elastic–plastic, and failure stages. The proposed can serve as a basis for dividing these stages. A cubic polynomial relationship was found between and elastic strain energy. AE cumulative energy and dissipated strain energy showed similar variation trends. Furthermore, based on , AE cumulative energy, and energy evolution theory, a failure prediction indicator () was proposed. This indicator can effectively identify precursor points of rock failure. A quantitative indicator for rock damage evolution under combined IR and AE action was created using as the characterization parameter of the rock damage variable, demonstrating high reliability. This research provides strong support for estimating rock states and guiding the design of rock engineering structures.
Fractures control fluid flow, solute transport, and mechanical deformation in crystalline media. They can be modeled numerically either explicitly or implicitly via an equivalent continuum. The implicit framework implies lower computational cost and complexity. However, upscaling heterogeneous fracture properties for its implicit representation as an equivalent fracture layer remains an open question. In this study, we propose an approach, the Equivalent Fracture Layer (EFL), for the implicit representation of fractures surrounded by low-permeability rock matrix to accurately simulate hydromechanical coupled processes. The approach assimilates fractures as equivalent continua with a manageable scale (≫1 μm) that facilitates spatial discretization, even for large-scale models including multiple fractures. Simulation results demonstrate that a relatively thick equivalent continuum layer (in the order of cm) can represent a fracture (with aperture in the order of μm) and accurately reproduce the hydromechanical behavior (i.e., fluid flow and deformation/stress behavior). There is an upper bound restriction due to the Young's modulus because the equivalent fracture layer should have a lower Young's modulus than that of the surrounding matrix. To validate the approach, we model a hydraulic stimulation carried out at the Bedretto Underground Laboratory for Geosciences and Geoenergies in Switzerland by comparing numerical results against measured data. The method further improves the ability and simplicity of continuum methods to represent fractures in fractured media.
Fracability evaluation for unconventional reservoir is critical to the selection of candidate zones for post-frac productivity and plays a key role in fracturing design. Historically, the prevailing models for assessing fracability have been largely relied on brittleness indices. Brittleness indices focus mainly on rock fracture characteristics and offers limited assessment of fracture surface area and the complexity of fracture network, which are more relevant to the practical production. We explored a new fracability evaluation model for unconventional reservoirs from the perspective of fracturing performance, which comprehensively characterizes the rock's ability to generate larger fracture surface areas, more shear fractures and complex fracture networks. The new fracability index considers both the physical processes of rock failure and fracture propagation, and is directly associated with the dynamic production capacities of reservoir. According to the analysis of energy conservation during hydraulic fracturing, we quantify the rock fracture surface area using the KGD and the PKN models. The ability of rock formation to generate shear fractures is mainly influenced by Poisson's ratio and mode II fracture toughness. Brittle mineral content and mineral heterogeneity are two vital criteria that significantly affect the complexity of fracture networks. Based on the logging and production data, this fracability model was applied to two types of unconventional reservoirs. Preliminary results show that this fracability model has an improved correlation with the pay zones and actual production, which is beneficial for optimizing fracturing strategies and identifying production sweet spots.
In many engineering applications, understanding gas adsorption and its induced swelling in nanoporous materials is crucial. In this study, we propose a novel coarse-grained molecular dynamics (CGMD) model with gas-gas, solid-solid, and gas-solid interactions explicitly controlled to achieve the coupling between gas transport and solid deformation at the microscale. The CGMD model has the capability to recover solid and gas properties, including density, Young's modulus of the solid, and viscosity of the gas to generate a broad range of swelling ratios relevant to nanostructures by using the innovative bead-spring chain networks. A comparison is made between gas transport through deformable and non-deformable nanochannels of varying sizes (35.4–123.9 nm), which is also compared with the macroscopic Hagen-Poiseuille equation. The proposed model has been further tested in a simplified nanoporous medium composed of four randomly distributed spherical solids. The Kozeny-Carman equation can generally describe the relationship between permeability and porosity, but small deviations are observed in the case of swelling porous media. Our results justify the effect of swelling on reducing gas permeability and provide a new approach to modeling gas transport in swelling porous media at the microscale within the framework of CGMD, with potential applications spanning nanofluidics, energy storage technologies, and environmental nanotechnology.
Rockburst hazards exhibit different spatiotemporal characteristics in deep tunnel excavation. Failure characteristics and energy evolution process of delayed and instantaneous rockburst of basalt rock were investigated based on single-sided unloading experiments under true triaxial conditions. High-speed photography and acoustic emission (AE) monitoring were used, and computed tomography (CT) scanning, fractal theory, and crack classification were employed for failure analysis. A three-dimensional damage model considering variable stiffness of testing machine was established to calculate the energy evolution of rock-machine system during the entire process of rockbursts. Results show that delayed rockburst includes three stages of small particles ejection, rock slab buckling, and violent mixed ejection, while instantaneous rockburst is characterized by rock slab spalling accompanied with slight particles ejection. Delayed rockburst exhibits a progressive failure mode of large-scale expansion of tensile cracks (before failure) to small-scale penetration of shear cracks (upon failure), while instantaneous rockburst shows a large-scale shear failure and abrupt penetration of shear planes upon failure. Delayed rockburst consumes less energy, and most of dissipated energy is converted into kinetic energy of ejected rock fragments, causing a higher intensity level of rockburst; instantaneous rockburst consumes more energy, but almost all dissipated energy comes from internal friction energy of shear failure, causing a higher scale of rock damage. Before rockburst failure, elastic strain energy stored in rock remains basically unchanged, while the energy stored in testing machine continuously decreases, indicating that rockburst is triggered by energy release of loading system. Energy dissipation rate (EDR) can be used as a precursory index for rock failure induced by quasi-static loading such as delayed rockburst. High EDR means damage intensification, stress drop, active AE events, and acceleration of shear crack expansion inside the rock. The findings of this study can provide new perspectives for the mechanisms and early warning of rockbursts.
Direct measurement of the real contact area of rock joints under normal loading is crucial for comprehending the subsurface geological processes. However, measuring this phenomenon quantitatively at site-scale or laboratory-scale is challenging. Here, we investigate the evolution mechanism of the real contact area in rock joints by conducting closure tests on artificial and saw-cut sandstone joints under normal stresses up to 50 MPa. Geometrical shapes of contact patches are quantified by the pressure-sensitive film using the adaptive threshold method. An extensive range of contact stress within contact patches is innovatively measured by integrating the results from multi-type pressure-sensitive films. Experimental results demonstrate that the real contact area increases with the increasing normal stress hyperbolically. Such a nonlinear contact evolution behavior can be attributed to the coalescence of adjacent contact patches. The fractal dimension of composite surface governs the geometrical shapes of contact patches and the distribution of contact stress. The relationship between patch areas and bearing loads follows the Hertzian theory when the patches are small, while it gradually becomes linear with the increasing patch size. A power model with exponential cut-off is proposed to predict the size distribution of contact patches. This work can provide new insights for estimating the patch-dependent seismic nucleation length and slip stability of subsurface joints.
Microseismic source mechanisms in underground mines can provide information about the rock mass response to mining. Conventional approaches to such studies rely upon moment tensor solutions that are susceptible to modeling assumptions and need reliable information about source locations and high-resolution velocity models. We propose the application of unsupervised clustering to group microseismic events into different classes directly from the waveform data such that the events in a specific class have similar source mechanisms. Our method has three main steps, first using spectral decomposition to separate the source terms from the path-receiver contributions in the observed amplitude spectra of events occurring in spatially dense clusters. Second, reducing the number of features from the source spectra using independent component analysis (ICA). Third, applying a Gaussian mixture model (GMM) to the reduced feature matrix to obtain event clusters. To test our method, we generate synthetic waveform data using the receiver network and the recorded microseismic event locations in an underground potash mine in Saskatchewan. Results show the ability of our method to separate events into different classes corresponding to differences in source mechanisms. Application to field data recorded in the mine during February 2021 successfully discriminates between blasts and microseismic events. The data recorded between 1 March and 30 June 2021 that contain microseismic events only are divided into two dominant classes. Using known moment tensors (MT) of some of these events for labeling, we interpret one of the two classes as having dominant double-couple mechanisms as compared to the other which most likely corresponds to the linear dipole-tensile mechanisms. Our method, combined with some expert knowledge such as MT of some larger magnitude events, can offer an assessment of source types of large microseismic populations as often encountered in induced seismicity.
Axial studies on cable bolts can be conducted using various scale testing apparatuses. Large scale testing, while providing a powerful platform for testing, is expensive and time consuming. This study presents details of a small scale pull out testing campaign on cable bolts and investigates the results achieved. Six popular types of cable bolts were studied using an anti rotation apparatus while encapsulated in cementitious grout and resin. The resin samples were tested under both monotonic and cyclic loading patterns. The results showed that grouted bulbed cables require higher displacement to reach their maximum load capacity which is lost at failure, while plain cables tend to hold lower loads for a longer time. Resin samples provided strain softening behaviour with low capacities, particularly in absence of cable indentation or bulbs. Cyclic loading tended to adversely affect the post peak behaviour of the resin samples, especially in the bulbed cables. Failed samples inspected after the testing suggested a non-uniform damage profile along the cable with extensive damage at the exit point transitioning into almost no damage at the entry point.