Numerical modeling of fracturing in permeable rocks via a micromechanical continuum model

Summary The paper presents a micromechanical approach to describe the failure of low‐permeability brittle rocks as a multiscale fracturing process based on a poroelastic microcrack‐damage model. Failure is formulated deep down at the fine pore scale as a material degradation phenomenon driven by microcrack growth that also impacts upon hydromechanical properties. A set of damage tensors describes the effect of dual‐scale porosities (nanopores and microcracks) on both the hydraulic and poroelastic rock properties. Essentially, the multiscale model reconstructs the coupling effect of hydromechanical forces at the continuum level from the ground up through the upscaling of multiphase interactions at the fundamental structural level of the material. As a result, many macroscopic characteristics emerge naturally such as friction angle, fracture properties, and most importantly, Biot's coefficient taking on a tensorial form that is generally anisotropic. The model is validated within the framework of finite elements to illustrate various baseline constitutive features such as the effect of microcrack growth on the nonlinear stress‐strain response and the induced anisotropy in the context of lab experimental tests and boundary value problems. Heterogeneities of the rock samples were incorporated by choosing material properties to be stochastic following Weibull and lognormal distributions. Numerical results appropriately replicated typical experimental observations where fracture localization and propagation are shown to be a multiscale phenomenon emerging from microcrack growth and coalescence at the microscale, with concomitant enhancement in fluid conductivity.