Poroelastic heterogeneity in the presence of a fresh fault: experimental insights and numerical modelling

During natural and induced seismic activities, pore fluid pressure within fault zones and their surrounding rock may respond differently to stress variations, introducing additional complexities to seismic hazard assessment. While theoretical investigations have recognized the influence of such poroelastic heterogeneity on fault instability, incorporating phenomena like slip-induced dilation or compaction, the chosen poroelastic properties in these studies lack robust constraints from experimental measurements. Addressing this gap, our study focuses on quantifying the heterogeneity of poroelastic properties in the presence of a fresh fault, aiming to elucidate the coupling between poroelasticy and fault dilatancy during fault slip.

In our experimental investigation, we examined the evolving dynamics of pore pressure both on- and off-fault in initially intact Westerly granite samples. Applying confining stress of 100 MPa and a pore pressure of 60 MPa at two sample ends to replicate crustal settings, we induced a sliding fault plane through loading to failure under a constant strain rate. In the faulted samples, we measured the pore pressure response under sudden step loading in the direction of the maximum compression σ₁. Each loading step of around 5 MPa was imposed incrementally increasing the differential stress from 5 MPa to approximately 80 MPa (frictional resistance) after achieving pore pressure equilibrium. Detailed measurements, including displacement, bulk deformation, differential stress, local pore pressure and acoustic emissions were recorded throughout these tests. A spring-slider model coupled with 1-D fluid diffusion was used to try to simulate experimental observations.

Our results indicate that both the shear zone and the bulk exhibit a diminishing Δp/Δσ₁ with increasing differential stress. Measurements within the fault zone consistently yield positive values, surpassing those off the fault, with the discrepancy more pronounced at lower stress levels. In regions farther away from the shear zone, the off-fault response Δp/Δσ₁ presents a smaller value compared to locations proximal to the fault zone and may even exhibit slight negativity. During fault slip, on-fault measurements exhibit an instantaneous increase upon step loading followed by a gradual decrease, as a result of the interplay between poroelasticity and fault dilatancy. These observations were effectively reproduced by the numerical model integrating the poroelastic measurements and rate-and-state fault friction with slip-dependent dilatancy. The implications of this investigation extend to an enriched understanding of the heterogeneity in poroelastic responses between fault zones and host rocks, serving as valuable benchmarks for informing future numerical simulations, particularly in the context of naturally formed fresh faults.