A quantitative and analytical approach is adopted to estimate two important parameters for coupled hydro-mechanical analysis at the scale of a fractured rock mass, namely the equivalent Biot effective stress coefficient α¯ and Skempton pore pressure coefficient B¯ . We derive formal expressions that estimate the two equivalent poroelastic coefficients from the properties of both the porous intact rock and the discrete fracture network, which includes fractures with different orientation, size, and mechanical properties. The coefficients are equivalent in the sense that they allow effectively predicting the volumetric deformation of the fluid-saturated fractured rock under an applied load in drained and undrained conditions. The formal expressions are validated against results from fully coupled hydro-mechanical simulations on systems with explicit representation of deformable fractures and rock blocks. We find that the coefficients are highly anisotropic as they largely vary with fracture orientations with respect to the applied stress tensor. For a given set of fracture and rock properties, B¯ increases with the ratio of normal to average stress undergone by the fractures, while the opposite occurs for α¯ . Additionally, both α¯ and B¯ increase with fracture density, which directly impacts the deformation caused by a load in undrained conditions. Because the effective stress variation is proportional to the applied load by (1 - α¯ B¯) , a factor that partly compensates for the decrease in equivalent rock stiffness caused by the fractures, a fully saturated fractured rock may deform less than an intact rock in undrained conditions, while the opposite occurs in dry conditions.

This work was funded by the Swedish nuclear fuel waste management company, Svensk Kärnbränslehantering AB (SKB), and the Nuclear Waste Management Organization (NWMO) in Toronto, Canada. The authors wish to thank Rima Ghazal for her help with the numerical models in 3DEC.

Peer reviewed