CASRock—a code to simulate THMC processes in fractured geological media

<p>This work presents the advancement of a self-developed Cellular Automata Software for engineering Rockmass fracturing processes (CASRock, http://en.casrock.cn/) in the applications of thermo-hydro-mechanical-chemical (THMC) processes in fractured geological media. It contains a series of previous developed numerical systems, namely EPCA for simulation of heterogeneous rock failure process, VEPCA for visco elastoplastic analysis, D-EPCA for rock dynamic response simulation, THMC-EPCA for coupled THMC processes in geological media and RDCA for simulation of rock cracking process from continuity to discontinuity. In CASRock, the non-isothermal, unsaturated fluid flow, mechanical process and chemical reaction are sequentially coupled by updating all the state variables using cellular automaton technique and finite difference method on spatial and temporal scale, respectively. The Lagrangian method is used to simulate the particle transport. The control equations, coupling scheme and numerical implementation are briefly introduced. Several applications, including  in the background of high level nuclear waste disposal are provided to show the abilities of CASRock in the simulation of coupling processes between physical fields. These applications include, (1) stability analysis of engineering rockmass under mechanical loading, (2) numerical study on coupled TM processes in hard rock pillar, (3) study on coupled THM processes in engineering barrier, (4) simulation of the THMC process in fractured rock.</p>

A Numerical Investigation on the Hydraulic Fracturing Effect of Water Inrush during Tunnel Excavation

When a high-pressure water source is located near a tunnel under excavation, water inrush is commonly associated with a hydraulic fracturing effect. To study the hydraulic fracturing effect of water inrush (HFEWI), flow-rock failure process analysis (F-RFPA2D) was adopted to simulate the water inrush process. The simulated results indicated that a stress disturbance area formed in front of the excavation face and that a hydraulic fracture zone formed in front of the karst cavity. Similarly, stress concentrations formed in front of the excavation face and the karst cavity. The hydraulic fracturing effect was characterized by stress concentration, and the local hydraulic crack propagation was the result of stress concentration. In addition, a pore pressure gradient formed in the crack-free area of the surrounding rock, and the occurrence of hydraulic cracking was the root cause of the significant change in water flow. When the hydraulic cracks initially formed and expanded, the zone of crack activity was large. As the cracks continued to expand, the range of activity decreased and finally concentrated directly in front of the excavation face. Additionally, the shapes of the water inrush channel obtained by the experimentation and numerical simulation were basically the same: semielliptical. During the evolution of hydraulic crack initiation, expansion, and penetration, the bottom of the excavated borehole was initially dry and then experienced seepage and water inrush. Finally, the minimum safe thickness of the rock wall was calculated to provide a safety guideline for this type of water inrush.

A Crack Propagation Control Study of Directional Hydraulic Fracturing Based on Hydraulic Slotting and a Nonuniform Pore Pressure Field

Hydraulic fracturing techniques for developing deeply buried coal reservoirs face routine problems related to high initial pressures and limited control over the fracture propagation direction. A novel method of directional hydraulic fracturing (DHF) based on hydraulic slotting in a nonuniform pore pressure field is proposed. A mechanical model is used to address crack initiation and propagation in a nonuniform pore pressure field, where cracks tend to rupture and propagate towards zones of high pore pressure for reducing the effective rock stress more. The crack initiation pressure and propagation morphology are analyzed by rock failure process analysis software. The numerical results show that the directional propagation of hydraulic fracturing cracks is possible when the horizontal stress difference coefficient is less than or equal to 0.5 or the slotting deviation angle is less than or equal to 30°. These findings are in good agreement with experimental results, which support the accuracy and reliability of the proposed method and theory.

Numerical Simulation of Rock Failure Process with a 3D Grain-Based Rock Model

A grain-based rock model was developed and applied to study mechanical characteristics and failure micromechanics in thick-walled cylinder and wellbore stability tests. The rock is represented as an assembly of tetrahedral blocks with bonded contacts. Material heterogeneity is modeled by varying the tensile strength at the block contacts. This grain-based rock model differs from previous disk/sphere particle-based rock models in its ability to represent a zero (or very low) initial porosity condition, as well as highly interlocked irregular block shapes that provide resistance to movement even after contact breakage. As a result, this model can reach higher uniaxial compressive strength to tensile strength ratios and larger friction coefficients than the disk/sphere particle-based rock model. The model captured the rock fragmentation process near the wellbore due to buckling and spalling. Thin fragments of rock similar to onion skins were produced, as observed in laboratory breakout experiments. The results suggest that this approach may be well suited to study the rock disaggregation process and other geomechanical problems in the rock excavation.

Relationship between Movement Laws of the Overlaying Strata and Time Space of the Mined-Out Volume

The study and accurate prediction of the movement of overburden rock mass and surface subsidence are crucial for a safe production in metal mines. This study investigates the relationship between the movement laws of overlaying strata and the time space of a mined-out volume using Rock Failure Process Analysis (RFPA) System. Furthermore, the movement, deformation, and failure laws of overlaying strata are examined in different positions when a goaf volume is certain and the failure behavior of the overlaying strata. This study analyzes the similarities and differences of the overlaying strata comparatively. Results show that, regardless of the movement range or subsidence value of the overlying rock mass, a power function relationship is observed between them and working face advancement. Setting the equation shows that the scope of the overlying rock mass is significant when the ratio of a certain position distance roof to the working face distance is small. The results provide a reference for controlling the displacement of the overlying rock mass and treating goaf.