A Monte Carlo Method for Probabilistic Hazard Assessment of Induced Seismicity due to Conventional Natural Gas Production

Abstract A Monte Carlo approach to probabilistic seismic‐hazard analysis is developed for a case of induced seismicity associated with a compacting gas reservoir. The geomechanical foundation for the method is the work of Kostrov (1974) and McGarr (1976) linking total strain to summed seismic moment in an earthquake catalog. Our Monte Carlo method simulates future seismic hazard consistent with historical seismic and compaction datasets by sampling probability distributions for total seismic moment, event locations and magnitudes, and resulting ground motions. Ground motions are aggregated over an ensemble of simulated catalogs to give a probabilistic representation of the ground‐motion hazard. This approach is particularly well suited to the specific nature of the time‐dependent induced seismicity considered. We demonstrate the method by applying it to seismicity induced by reservoir compaction following gas production from the Groningen gas field. A new ground‐motion prediction equation (GMPE) tailored to the Groningen field has been derived by calibrating an existing GMPE with local strong‐motion data. For 2013–2023, we find a 2% chance of exceeding a peak ground acceleration of 0.57 g and a 2% chance of exceeding a peak ground velocity of 22  cm/s above the area of maximum compaction. Disaggregation shows that earthquakes of M w  4–5, at the shortest hypocentral distances of 3 km, and ground motions two standard deviations above the median make the largest contributions to this hazard. Uncertainty in the hazard is primarily due to uncertainty about the future fraction of induced strains that will be seismogenic and how ground motion and its variability will scale to larger magnitudes.

[1]  Julian J. Bommer,et al.  Sigma: Issues, Insights, and Challenges , 2009 .

[2]  Robin K. McGuire,et al.  Probabilistic seismic hazard analysis: Early history , 2008 .

[3]  Fabio Moia,et al.  A survey of the induced seismic responses to fluid injection in geothermal and CO2 reservoirs in Europe , 2012 .

[4]  Benjamin Edwards,et al.  Predicting ground motion from induced earthquakes in geothermal areas , 2013 .

[5]  R. Davies,et al.  Induced seismicity and hydraulic fracturing for the recovery of hydrocarbons , 2013 .

[6]  Stefan Wiemer,et al.  Building Robust Models to Forecast the Induced Seismicity Related to Geothermal Reservoir Enhancement , 2012 .

[7]  Julian J. Bommer,et al.  The Influence of Magnitude Range on Empirical Ground-Motion Prediction , 2007 .

[8]  Gail M. Atkinson,et al.  Observations on Regional Variability in Ground-Motion Amplitudes for Small-to-Moderate Earthquakes in North America , 2009 .

[9]  W. Ellsworth Injection-Induced Earthquakes , 2013, Science.

[10]  A. Mossop An Explanation For Anomalous Time Dependent Subsidence , 2012 .

[11]  C. Cornell,et al.  Disaggregation of seismic hazard , 1999 .

[12]  Julian J. Bommer,et al.  Induced seismicity associated with Enhanced Geothermal Systems , 2007 .

[13]  Mark D. Petersen,et al.  Time-independent and Time-dependent Seismic Hazard Assessment for the State of California: Uniform California Earthquake Rupture Forecast Model 1.0 , 2007 .

[14]  Norman A. Abrahamson,et al.  Ground-Motion Attenuation Model for Small-To-Moderate Shallow Crustal Earthquakes in California and Its Implications on Regionalization of Ground-Motion Prediction Models , 2010 .

[15]  John Douglas,et al.  Modeling the difference in ground-motion magnitude-scaling in small and large earthquakes , 2011 .

[16]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[17]  F. H. Goutbeek,et al.  Seismic hazard due to small-magnitude, shallow-source, induced earthquakes in The Netherlands , 2006 .

[18]  Matthew C. Gerstenberger,et al.  Induced seismicity and its implications for CO2 storage risk , 2011 .

[19]  Julian J. Bommer,et al.  The impact of epistemic uncertainty on an earthquake loss model , 2005 .

[20]  Julian J. Bommer,et al.  A risk-mitigation approach to the management of induced seismicity , 2015, Journal of Seismology.

[21]  J. Bommer,et al.  Adapting earthquake actions in Eurocode 8 for performance‐based seismic design , 2006 .

[22]  M. Hitzman,et al.  Induced Seismicity Potential of Energy Technologies , 2013 .

[23]  Charles S. Mueller,et al.  2014 Update of the United States National Seismic Hazard Maps , 2012 .

[24]  A. Zollo,et al.  Ground‐Motion Prediction Equations for The Geysers Geothermal Area based on Induced Seismicity Records , 2013 .

[25]  Luca Malagnini,et al.  Effect of time-dependence on probabilistic seismic hazard maps and deaggregation for the central apennines, Italy , 2009 .

[26]  A. Barka,et al.  Heightened odds of large earthquakes near istanbul: An interaction-based probability calculation , 2000, Science.

[27]  Susan E. Hough,et al.  Shaking from Injection‐Induced Earthquakes in the Central and Eastern United States , 2014 .

[28]  J. D. Waal,et al.  On the rate type compaction behaviour of sandstone reservoir rock , 1986 .

[29]  Julian J. Bommer,et al.  The Influence of Strong-Motion Duration on the Seismic Response of Masonry Structures , 2004 .

[30]  BrianS-J. Chiou,et al.  An NGA Model for the Average Horizontal Component of Peak Ground Motion and Response Spectra , 2008 .

[31]  Julian J. Bommer,et al.  Modelling Seismic Hazard in Earthquake Loss Models with Spatially Distributed Exposure , 2006 .

[32]  B. Dost,et al.  Natural and induced seismicity , 2007 .

[33]  A. McGarr,et al.  Seismic moments and volume changes , 1976 .

[34]  G. Weatherill,et al.  OpenQuake Engine: An Open Hazard (and Risk) Software for the Global Earthquake Model , 2014 .

[35]  K. Assatourians,et al.  EqHaz: An Open‐Source Probabilistic Seismic‐Hazard Code Based on the Monte Carlo Simulation Approach , 2013 .

[36]  Jenny Suckale,et al.  Induced Seismicity in Hydrocarbon Fields , 2009 .

[37]  Bernard Dost,et al.  The August 16, 2012 earthquake near Huizinge (Groningen) , 2013 .

[38]  T. Van Eck,et al.  Scaling of peak ground acceleration and peak ground velocity recorded in the Netherlands , 2004 .

[39]  R. Mcguire Probabilistic seismic hazard analysis and design earthquakes: Closing the loop , 1995, Bulletin of the Seismological Society of America.

[40]  John B. Mander,et al.  Seismic performance of an unreinforced masonry building: An experimental investigation , 2009 .

[41]  N. Abrahamson,et al.  Alternative Approaches to Modeling Epistemic Uncertainty in Ground Motions in Probabilistic Seismic‐Hazard Analysis , 2014 .

[42]  S. Harmsen,et al.  Deaggregation of probabilistic ground motions in the central and eastern United States , 1999 .

[43]  S. Akkar,et al.  Empirical ground-motion models for point- and extended-source crustal earthquake scenarios in Europe and the Middle East , 2014, Bulletin of Earthquake Engineering.

[44]  Vincenzo Convertito,et al.  From Induced Seismicity to Direct Time-Dependent Seismic Hazard , 2011 .

[45]  Julian J. Bommer,et al.  THE PREDICTION AND USE OF PEAK GROUND VELOCITY , 2006 .

[46]  Christian D. Klose,et al.  Mechanical and statistical evidence of the causality of human-made mass shifts on the Earth’s upper crust and the occurrence of earthquakes , 2012, Journal of Seismology.

[47]  Julian J. Bommer,et al.  On the Selection of Ground-Motion Prediction Equations for Seismic Hazard Analysis , 2010 .

[48]  Roger M.W. Musson,et al.  The use of Monte Carlo simulations for seismic hazard assessment in the U.K. , 2000 .

[49]  V. V. Kostrov,et al.  AND SEISMIC FLOW OF ROCK , 1974 .

[50]  R. Puglia,et al.  Pan-European ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5 %-damped PSA at spectral periods up to 3.0 s using the RESORCE dataset , 2013, Bulletin of Earthquake Engineering.

[51]  Rodolfo Puglia,et al.  Erratum to: Pan-European ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5 %-damped PSA at spectral periods up to 3.0 s using the RESORCE dataset , 2014, Bulletin of Earthquake Engineering.

[52]  C. Cornell Engineering seismic risk analysis , 1968 .

[53]  Julian J. Bommer,et al.  The Variability of Ground-Motion Prediction Models and Its Components , 2010 .

[54]  Alan L. Kafka,et al.  A Monte Carlo approach to seismic hazard analysis , 1999 .

[55]  Julian J. Bommer,et al.  Consistent Source-to-Site Distance Metrics in Ground-Motion Prediction Equations and Seismic Source Models for PSHA , 2012 .

[56]  R. Musson DETERMINATION OF DESIGN EARTHQUAKES IN SEISMIC HAZARD ANALYSIS THROUGH MONTE CARLO SIMULATION , 1999 .

[57]  Brendon A. Bradley,et al.  Correlation of Significant Duration with Amplitude and Cumulative Intensity Measures and Its Use in Ground Motion Selection , 2011 .

[58]  Stewart w. Smith,et al.  SOURCE PARAMETERS OF EARTHQUAKES, AND DISCRIMINATION BETWEEN EARTHQUAKES AND NUCLEAR EXPLOSIONS. , 1968 .

[59]  S. Bourne,et al.  A seismological model for earthquakes induced by fluid extraction from a subsurface reservoir , 2014 .