Induced Seismicity in the Groningen Field a Statistical Assessment of Tremors Along Faults in a Compacting Reservoir
Total Page:16
File Type:pdf, Size:1020Kb
Induced seismicity in the Groningen field A statistical assessment of tremors along faults in a compacting reservoir Rick Wentinck Datum November 2015 Editors Jan van Elk & Dirk Doornhof General Introduction In this report, a statistical model is presented to better understand and predict the tremor rates and magnitudes from trends in the past at a local scale following from oil/gas producing reservoirs. The local scale can be a small gas field or a region of several square kilometres of a large field, such as the Groningen field. The state-driven model supposes that the tremors are generated along the faults in the region of interest. These faults become critically stressed during the compaction of the reservoir in which the reservoir pressure drops. Herewith, the state which triggers the tremors is the mean stress-state along the faults in this region. The model is executed as a Monte Carlo simulation and uses stochastic variables which are related to probabilities that a tremor occurs and that it has a certain magnitude. The document explores further investigations into the seismological model analyzing the historical record of earthquakes. This should be read along-side the reports on the activity rate model: Papers: 1. Bourne, S. J., S. J. Oates, J. van Elk, and D. Doornhof (2014), A seismological model for earthquakes induced by fluid extraction from a subsurface reservoir, J. Geophys. Res. Solid Earth, 119, 8991–9015, doi:10.1002/2014JB011663.) 2. Bourne S.J., S.J. Oates, J.J. Bommer, B. Dost, J. van Elk, D. Doornhof, A Monte Carlo method for probabilistic hazard assessment of induced seismicity due to conventional gas production, Bulletin of the Seismological Society of America, V.105, no. 3, June 2015 in press. Reports: 3. Stephen Bourne and Steve Oates, An activity rate model of induced seismicity within the Groningen Field, (Part 1), February 2015. 4. An activity rate model of seismicity induced by reservoir compaction and fault reactivation in the Groningen gas field (Part 2), S.J.Bourne, S.J. Oates, June, 2015 5. Review of “An Activity rate model of seismicity induced by reservoir compaction and fault reactivation in the Groningen gas fields”, Ian Main, September 2015 Title Induced seismicity in the Groningen field - A statistical Date November 2015 assessment of tremors along faults in a compacting Initiator NAM reservoir Author(s) Rick Wentinck Editors Jan van Elk Dirk Doornhof Organisation Shell Organisation NAM Place in the Study Theme: Seismological Model Study and Data Comment: Acquisition Plan Understanding of earthquake historical record and prediction of earthquakes. Directliy linked (1) Geomechanical modelling of ruptures, research (2) Statistics of earthquakes (3) Ground Motion Prediction (incl. duration and spectral acceleration)m. Used data Historical record of earthquakes based on measurements from geophone and accelerometer network. Associated organisation Assurance Internal. Restricted SR.15.11335 Induced seismicity in the Groningen field - statistical assessment of tremors along faults in a compacting reservoir by H.M. Wentinck This document is classified as Restricted. Access is allowed to Shell personnel, designated Associate Companies and Contractors working on Shell projects who have signed a confidentiality agreement with a Shell Group Company. ‘Shell Personnel’ includes all staff with a personal contract with a Shell Group Company. Issuance of this document is restrictedto staff employed by a Shell Group Company. Neither the whole nor any part of this document may be disclosed to Non-Shell Personnel without the prior written consent of the copyright owners. Copyright Shell Global Solutions International, B.V., 2015. Shell Global Solutions International B.V., Rijswijk Further electronic copies can be obtained from the Global Information Centre. SR.15.11335 - I - Restricted Executive Summary In this report, we present a statistical model to understand and predict the tremor rates and magnitudes from trends in the past at a local scale following from oil/gas producing reservoirs. The local scale can be a small gas field or a region of several square kilometers of a large field, such as the Groningen field. The state-driven model supposes that the tremors are generated along the faults in the region of interest. These faults become critically stressed during the compaction of the reservoir in which the reservoir pressure drops. Herewith, the state which triggers the tremors is the mean stress-state along the faults in this region. The model is executed as a Monte Carlo simulation and uses stochastic variables which are related to probabilities that a tremor occurs and that it has a certain magnitude. The relative likelihood of a seismic fault failure depends on the Weibull distribution func- tion. The tremor rates are normalised to the observed tremor rates using the Poisson probability distribution function. The seismic moments of the tremors follow from a Pareto distribution function, as is commonly done when analysing of natural and man- made induced seismicity. The model includes that the tremors reduce the mean stress on the faults in proportion to the seismic moment released. The model does not differentiate between a-seismic and seismic fault slip based on field or fault attributes. This could be seen as a shortcoming but we have no clue what this could be, so far. The model has been applied to 6 relatively small regions in the Groningen field and in the Annerveen and Eleveld fields south of the Groningen field. The model fit parameters have been derived from these regions by comparing the modelled and observed tremors. A sensitivity study has been done and the predictive potential of the model has been explored. The fit parameters of the model, τfail and cM0 have values which can be understood and constrained from geomechanical calculations. τfail is a typical mean shear stress along the faults for which the relative likelihood for seismic fault failure substantially increases. cM0 relates the mean stress reduction along the faults by the tremors generated by uniaxial compaction. Remarkable is the high value of the Weibull shape parameter kW in relation to variations in parameters for rock properties and fault geometry which determine the onset of shear- type fault failure. The Weibull shape parameter determines the period of the transition from no tremors to a constant tremor rate in the case of a constant compaction rate. It’s high value may indicate that the tremors are predominantly determined by only a few of these parameters or by another mechanism that differs from a-seismic slip. The present model cannot reproduce the drastic reduction in the tremor rate in the Eleveld field by a mean field stress relaxation process. This could point to a shortcoming SR.15.11335 - II - Restricted in the submodel for the stress relaxation or to a shortcoming in the submodel for the generation of the tremors. The latter submodel is driven by the Poisson process and the Weibull distribution probability distribution function for the relative likelihood of a tremor. Presumably, the Poisson process expression which normalises the number of tremors should include a term related to the rate of compaction. With reservations, we illustrate that the model has some predictive power about tremor rates in a certain region, using the observed tremors in this region in the previous period. In essence, the state-driven model consistently reproduces the gradual onset of the tremor rate followed by a more or less constant tremor rate. Taking a constant and uniform pres- sure drop over the Groningen field, it must be noted that the model does not capture a small but significant increase in the tremor rate over the last three to four years in the regions around Ten Boer and Loppersum. We found no correlations between the regional variation of the observed tremor den- sity and the fault dip angle, fault throw and fault azimuth angle in the Groningen field. A remarkable observation is that most tremors with a minimum magnitude of M = 2.5 can be associated with faults with high dip angles of about 70◦ and throws of about 100 m, see Appendix C. These faults could potentially release most of the stored energy from compaction and the imbalance between the field stresses, see Appendix B, B.2. § We recommend to follow the effects of local gas production reductions in the Gronin- gen field to validate and improve this model, especially in relation to the possible effect of stress relaxation by tremors. The present model does not differentiate between seismic and a-seismic fault slip because we do not know yet which field or fault attributes are determining here. If this knowledge becomes available, it could be used to constrain the model parameters further. Table of Contents 1 Introduction 6 2 Model 9 2.1 Assumptionsandequations . 9 2.1.1 Assumptions............................... 9 2.1.2 Stressstate ............................... 10 2.1.3 Therelativelikelihoodofatremor . 11 2.1.4 Frequency-magnitude-relation-of-tremors . ........ 13 2.2 Energy release by tremors and mean stress reduction . ........ 13 2.2.1 Liberatedenergy,slipandstress . 13 2.2.2 Reduction in the mean shear stress along the faults in the region of interestbytremors .. ... .. .. .. .. ... .. .. .. .. ... 17 2.2.3 Relationwithobservations . 18 2.3 Relation between frequency-magnitude relationship and the asperity dis- tribution..................................... 19 2.3.1 Otherb-values ............................. 20 2.4 Pareto, Weibull and Poisson probability distributions