Pre-failure strain localization in siliclastic rocks: a comparative study of laboratory and numerical approaches

Abstract

We combined novel laboratory techniques and numerical modeling to investigate (a)seismic preparatory processes associated with deformation localization during a triaxial failure test on a dry sample of Berea sandstone. Laboratory observations were quantified by measuring strain localization on the sample surface with a distributed strain sensing (DSS) array, utilizing optical fibers, in conjunction with both passive and active acoustic emission (AE) techniques. A physics-based computational model was subsequently employed to understand the underlying physics of these observations and to establish a spatio-temporal correlation between the laboratory and modeling results. These simulations revealed three distinct stages of preparatory processes: (i) highly dissipative fronts propagated towards the middle of the sample correlating with the observed acoustic emission locations; (ii) dissipative regions were individuated in the middle of the sample and could be linked to a discernible decrease of the P-wave velocities; (iii) a system of conjugate bands formed, coalesced into a single band that grew from the center towards the sample surface and was interpreted to be representative for the preparation of a weak plane. Dilatative lobes at the process zones of the weak plane extended outwards and grew to the surface, causing strain localization and an acceleration of the simulated deformation prior to failure. This was also observed during the experiment with the strain rate measurements and spatio-temporally correlated with an increase of the seismicity rate in a similar rock volume. The combined approach of such laboratory and numerical techniques provides an enriched view of (a)seismic preparatory processes preceding the mainshock.