Hybrid implicit–explicit XFEM simulation of injection-induced seismicity: resolving multi-scale rupture nucleation and dynamics
- 1Technische Universität Berlin (Berlin, Germany)
- 2University of Bergen (Bergen, Norway)
Abstract
This study presents a novel hybrid implicit–explicit time-integration strategy for simulating injection-induced fault reactivation within a fully coupled poroelastodynamic extended finite element (XFEM) framework. Modeling these systems is inherently challenging due to the extreme separation of temporal scales – ranging from hours to days for fluid diffusion to milliseconds for dynamic rupture. Although the fully implicit formulation can provide enhanced numerical stability, stability alone does not guarantee sufficient temporal accuracy during nonlinear fault rupture. The results demonstrate close agreement between simulations with minimum time-step sizes of 10 ⁻⁵ and 10 ⁻⁶ s for the investigated slip and event-scale observables on the adopted spatial mesh, indicating temporal convergence for these quantities at 10 ⁻⁵ s. This time-step size is substantially smaller than the mesh-based explicit stability estimate. To reduce computational cost, the proposed strategy employs an implicit solver for quasi-static loading and nucleation and switches to a central-difference-based mechanical solver with implicit pressure integration once a prescribed slip-velocity threshold is exceeded. At a switching threshold of 0.01 m s ⁻¹ , the hybrid simulation requires approximately 10.0 h, compared with 26.4 h for the fully implicit reference, corresponding to a computational-time reduction of 62.12 %. The hybrid strategy reproduces the overall rupture sequence, stress evolution, final accumulated fault slip, and seismic moment. However, differences remain in the timing and magnitude of peak slip velocities, rupture-propagation speeds, and transient particle motion. These findings demonstrate the computational potential of the hybrid formulation for the investigated two-dimensional configuration, while highlighting the need for further assessment of transient dynamic accuracy and computational performance across different spatial discretizations and geological conditions.
Versions
Related
- Optimizing induced seismicity monitoring for CO2 injection at Trüllikon, Switzerland using risk-informed guidelines
- Development of a Combined Detection and Relocation Workflow (TM2DD) to Characterize Escalating Hydraulic Fracturing Induced Seismicity in the Karnes Trough Fault Zone of the Eagle Ford Shale Play, Texas
- Escalating soil CO ₂ degassing correlates with ground uplift and seismicity during the ongoing unrest of Campi Flegrei (2005–2025)
- Tectonic regionality and geophysical controls of the fractal dimension and b-value of seismicity in the Chinese mainland
- Temporal evolution and decline of seismicity associated with cyclic operations at the Serajeh underground gas storage, Iran
BibTeX
@article{sabah2026hybrid,
title = {{Hybrid implicit–explicit XFEM simulation of injection-induced seismicity: resolving multi-scale rupture nucleation and dynamics}},
author = {Mohammad Sabah and Mauro Cacace and Inga Berre and Iman R. Kivi and Hannes Hofmann},
journal = {Geoscientific Model Development},
year = {2026},
month = nov,
doi = {10.5194/gmd-19-10447-2026},
url = {https://gmd.copernicus.org/articles/19/10447/2026/},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Hybrid implicit–explicit XFEM simulation of injection-induced seismicity: resolving multi-scale rupture nucleation and dynamics
AU - Mohammad Sabah
AU - Mauro Cacace
AU - Inga Berre
AU - Iman R. Kivi
AU - Hannes Hofmann
PY - 2026
DA - 2026/11/04
JO - Geoscientific Model Development
DO - 10.5194/gmd-19-10447-2026
UR - https://gmd.copernicus.org/articles/19/10447/2026/
AB - This study presents a novel hybrid implicit–explicit time-integration strategy for simulating injection-induced fault reactivation within a fully coupled poroelastodynamic extended finite element (XFEM) framework. Modeling these systems is inherently challenging due to the extreme separation of temporal scales – ranging from hours to days for fluid diffusion to milliseconds for dynamic rupture. Although the fully implicit formulation can provide enhanced numerical stability, stability alone does not guarantee sufficient temporal accuracy during nonlinear fault rupture. The results demonstrate close agreement between simulations with minimum time-step sizes of 10 ⁻⁵ and 10 ⁻⁶ s for the investigated slip and event-scale observables on the adopted spatial mesh, indicating temporal convergence for these quantities at 10 ⁻⁵ s. This time-step size is substantially smaller than the mesh-based explicit stability estimate. To reduce computational cost, the proposed strategy employs an implicit solver for quasi-static loading and nucleation and switches to a central-difference-based mechanical solver with implicit pressure integration once a prescribed slip-velocity threshold is exceeded. At a switching threshold of 0.01 m s ⁻¹ , the hybrid simulation requires approximately 10.0 h, compared with 26.4 h for the fully implicit reference, corresponding to a computational-time reduction of 62.12 %. The hybrid strategy reproduces the overall rupture sequence, stress evolution, final accumulated fault slip, and seismic moment. However, differences remain in the timing and magnitude of peak slip velocities, rupture-propagation speeds, and transient particle motion. These findings demonstrate the computational potential of the hybrid formulation for the investigated two-dimensional configuration, while highlighting the need for further assessment of transient dynamic accuracy and computational performance across different spatial discretizations and geological conditions.
SN - 1991-9603
ER -