[Miller] Modelling the Bohemian Earthquake Swarms: Application of a poro-thermo-elastic plastic model with weakening, damage and multiphase flow in GPU

German Title: Modelling the Bohemian Earthquake Swarms: Application of a poro-thermo-elastic plastic model with weakening, damage and multiphase flow in GPU

Abbreviation: 315

Current Status: completed

Main Applicant:Prof. Dr. Stephen Miller

Resources Recipient

Other Persons

Dr. Thomas Heinze

Conveyor Begin:
Conveyor End:
Conveyor Duration:
Year: 2013


We propose to apply a newly developed model to numerically investigate earthquake swarms in general and the Bohemian swarms in particular. The model simulates a pore-elastic plastic rheology coupled with a non-linear diffusion model for pressure propagation, where the nonlinearity arises through the effective stress-dependence of the permeability. A two-phase (water and gas) flow model is currently being implemented. Numerical cracks in the model nucleate and grow in response to; i) far-field stresses developed at the boundaries, ii) stress perturbations arising from internal crack growth, and iii) poro-elastic stressing. Recent developments that include hardening, softening, and a damage growth model reproduce the entire (load-failure-unload) stress-strain curve observed in laboratory rock experiments. In this project, we propose to: 1) Include a heat flow model and introduce thermo-elasticity; 2) extend the model to 3-dimensions, 3) implement the model on the Graphics Processing Unit (GPU) platform, and 4) apply the model to the Bohemian swarms. The GPU platform allows fast and very high-resolution simulations in comparison to the standard CPU architecture. The numerical model is ideally suited to study the Bohemian swarms because the dominant physics of the Bohemian earthquake swarms are simulated, namely the complex interactions and feedbacks of CO2 injection, seismic slip (with the concomitant permeability enhancement) fluid flow, and advective heat flow. The model will be used to address: a) triggering mechanisms that generate swarm-earthquake activity, b) characterize fault-related fluid transport processes and effects of the fluid/rock interaction during migration, and c) investigate the system’s dependence on pressure and temperature.

Related Publications

Heinze, Thomas, Galvan, Boris (2016). "Novel numerical strategy for solving strongly coupled elastoplastic damage models with explicit return algorithms: Application to geomaterials" International Journal of Solids and Structures 80 p64-72

Heinze, Thomas, Jansen, Gunnar, Galvan, Boris, Miller, Stephen A (2016). "Systematic study of the effects of mass and time scaling techniques applied in numerical rock mechanics simulations" Tectonophysics 684 p4-11

Heinze, Thomas, Galvan, Boris, Miller, Stephen Andrew (2015). "A new method to estimate location and slip of simulated rock failure events" Tectonophysics 651 p35-43

Heinze, Thomas (2015). "Development and application of coupled THM solvers to estimate rock failure events from laboratory to field scales" Dissertation

Heinze, Thomas, Galvan, Boris, Miller, Stephen Andrew (2015). "Modeling porous rock fracturing induced by fluid injection" International Journal of Rock Mechanics and Mining Sciences p133-141