Revisiting Stimulation Mechanism at Fenton Hill and an Investigation of the Influence of Fault Heterogeneity on the Gutenberg-Richter b-value for Rate-and-State Earthquake Simulations

In this work, we developed a numerical modeling framework for application to subsurface reservoir engineering problems related to fluid flow, heat transfer, and mechanical deformation in fractured porous media. An existing reservoir simulator was extended to include poroelastic and thermoelastic effects. We introduced a novel approach to couple poroelastic and thermoelastic stresses with fracture deformation calculations. In addition, we implemented a rate-and-state friction model capable of describing the earthquake rupture process on two-dimensional heterogeneous fault surfaces. A spatial random field model was used to generate heterogeneous, spatially correlated, fractal distributions of fault properties such as stress, friction, and permeability. We applied the numerical model to investigate reservoir stimulation efforts at the Fenton Hill enhanced geothermal system test site. A geologic conceptual model of the site was developed based on several independent datasets. Several scenarios were simulated to isolate the roles of fluid pressurization of the fractures, poroelastic stress, and thermal stress. Comparing the simulation results with the actual earthquake catalog data indicates that poroelastic and thermoelastic effects likely did play a role during the Fenton Hill stimulation treatments. Our analysis also suggests that shear stimulation did not contribute to enhanced permeability at the Fenton Hill reservoir. In another study, we hypothesized that by introducing heterogeneous and spatially correlated patterns of fault properties, it would be possible to simulate earthquake sequences that exhibit power law scaling of frequency-magnitude distributions using a rate-and-state friction model. We examined the influence of heterogeneous distributions of initial shear stress, initial normal stress, dynamic friction coefficient, and permeability. Our model was able to simulate earthquake sequences exhibiting Gutenberg-Richter-type frequencymagnitude distributions with b-values close to 1 when the fault had heterogeneous stress or friction. Permeability heterogeneity did not have an effect. This study has implications for understanding how the Gutenberg-Richter b-value is influenced by natural geologic conditions and injection well operations.