Published July 2025 | Published
Journal Article Open

Flow2Quake, an integrated multiphase flow, geomechanical and seismicity model for efficient forecasting of injection and extraction induced earthquakes

  • 1. ROR icon California Institute of Technology
  • 2. ROR icon Virginia Tech
  • 3. ROR icon École Polytechnique

Abstract

Efforts to secure and decarbonize the energy sector are driving various subsurface reservoir operations. These operations carry a risk of inducing surface deformation and earthquakes. To assess these risks, modeling tools integrating fluid flow, geomechanical and seismicity modeling are needed. Here, we demonstrate the use of an efficient Vertical Flow Equilibrium (VFE) multiphase fluid flow model in an integrated framework for deformation and seismicity modeling both under fluid extraction or injection configurations. The VFE-computed spatio-temporal pressure evolution is fed to a geomechanical module to compute surface deformation and stress changes in and around the reservoir. Stress changes feed a seismicity module to calculate earthquake probabilities. First, we apply the benchmarked model to gas extraction from Groningen. There, we can reduce the variance of pressure measurements by 38% with respect to a pre-existing single phase flow model while remaining computationally efficient. The surface deformation and seismicity simulations show remarkable agreement with observed data. Second, we study induced seismicity due to CO sequestration in the Decatur phase 1 project. We find that, for the Decatur phase 1 project, poroelastic stress changes can account for most of the non-clustered observed seismicity within modeling uncertainties. Finally we simulate scenarios for CO sequestration using the Quest field. The sloping reservoir topography significantly impacts the predicted position of the CO plume but the effects on geomechanical deformation (and seismicity) are minimal. Incorporating VFE models with geomechanical and seismicity forecasts with real-world case applications can allow real-time hazard assessment and mitigation procedures.

Copyright and License

© 2025 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

Acknowledgement

This study was supported by the NSF/IUCRC Geomechanics and Mitigation of Geohazards (National Science Foundation, USA award # 1822214). T.L. acknowledges École Polytechnique, FR for its financing support through the Chaire Energies Durables. M.A. Acknowledges funding from the Swiss National Science Foundation through grant P2ELP2195127 and from Caltech’s Resnick sustainability institute, USA. We gratefully acknowledge data and support from Shell Global Solutions, USA . We gratefully acknowledge discussions with Guanli Wang and Adrian Moure. Authors declare that they have no competing interests.

Data Availability

All data and codes can be found in: doi: 10.17605/OSF.IO/WSH32.

Supplemental Material

MMC S1. The appendix contains information on (A) Fluid properties (B) The multiphase flow model derivation (C) Benchmarking of the MP-Flow model (PDF)

Files

1-s2.0-S1750583625000866-main.pdf
Files (11.8 MB)
Name Size Download all
md5:c08725dedf9b143f398da1599fca932c
9.5 MB Preview Download
md5:c1b42b39455480048ab4a4eaa218f820
2.3 MB Preview Download

Additional details

Created:
June 6, 2025
Modified:
June 6, 2025