The Influence of Grain Crushing and Pore Collapse on the Formation of Faults

Abstract During an earthquake, slip occurs in a localized shear zone that features a heavily granulated fault core that can be characterized as a shear band. We study the formation of this fault core in a granular rock such as sandstone by developing a model of crushable granular media within the framework of Breakage Mechanics. This model accounts for the evolution of the grain size distribution, while also accounting for the co‐evolution of the solid fraction. An enrichment with the Cosserat continuum allows for the model to predict finite‐width shear bands. The model is then calibrated against experimental data taken from tests on Bentheim sandstone, and a parametric study of the mechanical parameters is conducted using linear stability analysis. We find that for deeply‐buried rocks the shear bands have a compactive component, and the initial value of the solid fraction does not play a strong role in the initial band thickness, but can influence the rate of delocalization of the band. Post‐localization behavior is studied with the finite element method, which shows the formation of zones of dilation outside the band in addition to the compaction within the band. Using a modified Kozeny–Carman permeability law, it is shown that within the band the permeability reduces by several orders of magnitude, but can increase outside the band. Our results highlight the importance of modeling grain size and solid fraction evolution as they exert a controlling influence on hydromechanical properties that play an important role in fault formation and seismic slip.