A microstructural rock model based on the distinct element method employing the Subspring Network contact model with rigid, Breakable, Voronoi-shaped grains (SNBV model) is proposed. The model consists of a mesh (3D Voronoi tessellation) of rigid, breakable, Voronoi blocks. The SNBV model is a microstructural rock model because it is a discrete model that can mimic rock microstructure at the grain scale. SNBV material mimics the microstructure of angular, interlocked, breakable grains with interfaces that may have an initial gap and can sustain partial damage. The model embodies the microstructural features and damage mechanisms that occur at the grain scale: initial microcrack fabric; heterogeneity-induced local tension; and intergranular and transgranular damage. The heterogeneity-induced local tension can be introduced in a controlled fashion that is not tied directly to the shape and packing of the grains and the interface stiffnesses. The synthetic material exhibits behavior during direct-tension and triaxial compression tests that matches the behavior of compact rock. The material can be calibrated to match the standard material properties and characteristic stresses of pink Lac du Bonnet granite. The material properties consist of Young's modulus and Poisson's ratio corresponding with uniaxial compression and Young's modulus corresponding with direct tension, as well as tensile strength, crack-closure stress, crack-initiation stress, secondary crack-initiation stress to mark the onset of grain breakage, crack-damage stress, and compressive strengths up to 4 MPa confinement. The model is suitable for studying the grain-scale micromechanics of brittle rock fracture.

Rock fracture toughness is a critical parameter for optimizing reservoir stimulation during deep resource extraction. This index characterizes the in situ resistance of rocks to fracture and is affected by high temperature, in situ stress, thermal shock, and chemical corrosion, etc. This review comprehensively examines research on rock fracture properties in situ environments over the past 20 years, analyses the influences of various environmental factors on rock fracture, and draws the following conclusions: (i) Environmental factors can significantly affect rock fracture toughness through changing the internal microstructure and grain composition of rocks; (ii) While high temperature is believed to reduce the rock strength, several studies have observed an increase in rock fracture toughness with increasing temperature, particularly in the range between room temperature and 200 degrees C; (iii) In addition to a synergistic increase in fracture toughness induced by both high temperature and high in situ stress, there is still a competing effect between the increase induced by high in situ stress and the decrease induced by high temperature; (iv) Thermal shock from liquid nitrogen cooling, producing high temperature gradients, can surprisingly increase the fracture toughness of some rocks, especially at initial temperatures between room temperature and 200 degrees C; and (v) Deterioration of rock fracture toughness occurs more rapidly in acidic environments than that in alkaline environments. In addition, this review identified current research trends and suggested some potential directions to provide suggestions for deep subsurface resource extraction. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/

The geometric properties of fracture surfaces significantly influence shear-seepage in rock fractures, introducing complexities to fracture modelling. The present study focuses on the hydro-mechanical behaviours of rough rock fractures during shear-seepage processes to reveal how dilatancy and fracture asperities affect these phenomena. To achieve this, an improved shear-flow model (SFM) is proposed with the incorporation of dilatancy effect and asperities. In particular, shear dilatancy is accounted for in both the elastic and plastic stages, in contrast to some existing models that only consider it in the elastic stage. Depending on the computation approaches for the peak dilatancy angle, three different versions of the SFM are derived based on Mohr-Coulomb, joint roughness coefficient-joint compressive strength (JRC-JCS), and Grasselli's theories. Notably, this is a new attempt that utilizes Grasselli's model in shear- seepage analysis. An advanced parameter optimization method is introduced to accurately determine model parameters, addressing the issue of local optima inherent in some conventional methods. Then, model performance is evaluated against existing experimental results. The findings demonstrate that the SFM effectively reproduces the shear-seepage characteristics of rock fracture across a wide range of stress levels. Further sensitivity analysis reveals how dilatancy and asperity affect hydraulic properties. The relation between hydro-mechanical properties (dilatancy displacement and hydraulic conductivity) and asperity parameters is analysed. Several profound understandings of the shear-seepage process are obtained by exploring the phenomenon under various conditions. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/ licenses/by-nc-nd/4.0/).