scholarly journals Probabilistic earthquake locations of induced seismicity in the Groningen region, the Netherlands

2020 ◽  
Vol 222 (1) ◽  
pp. 507-516 ◽  
Author(s):  
Jonathan D Smith ◽  
Robert S White ◽  
Jean-Philippe Avouac ◽  
Stephen Bourne

SUMMARY The Groningen gas reservoir, situated in the northeast of the Netherlands, is western Europe’s largest producing gas field and has been in production since 1963. The gas production has induced both subsidence and seismicity. Seismicity is detected and located using the Koninklijk Nederlands Meteorologisch Instituut shallow-borehole array for the period 2015–2017, incorporating the back projection techniques of QuakeMigrate and the nonlinear location procedure to constrain earthquake locations and depths. The uncertainties on the estimated depths are estimated taking into account velocity model, changes in station array geometry and uncertainties in the measurement of arrival times of the P and S waves. We show that the depth distribution of seismicity is consistent with nucleation within the reservoir (28 per cent) or in the overburden (60 per cent) within ∼500 m from the top of the reservoir. Earthquakes with hypocentres in the overburden likely originate from overlying Zechstein anhydrite caprock. Based on their depth distribution, it seems like the earthquakes are primarily driven by the elastic strain in the reservoir and overburden, induced by the reservoir compaction. We estimate the probability of earthquakes nucleating beneath the reservoir in the underlying Carboniferous limestone and basement, to be no more than 12 per cent.

Author(s):  
Pauline P. Kruiver ◽  
Manos Pefkos ◽  
Erik Meijles ◽  
Gerard Aalbersberg ◽  
Xander Campman ◽  
...  

AbstractIn order to inform decision-making regarding measures to mitigate the impact of induced seismicity in the Groningen gas field in the Netherlands, a comprehensive seismic risk model has been developed. Starting with gas production scenarios and the consequent reservoir compaction, the model generates synthetic earthquake catalogues which are deployed in Monte Carlo analyses, predicting ground motions at a buried reference rock horizon that are combined with nonlinear amplification factors to estimate response spectral accelerations at the surface. These motions are combined with fragility functions defined for the exposed buildings throughout the region to estimate damage levels, which in turn are transformed to risk in terms of injury through consequence functions. Several older and potentially vulnerable buildings are located on dwelling mounds that were constructed from soils and organic material as a flood defence. These anthropogenic structures are not included in the soil profile models used to develop the amplification factors and hence their influence has not been included in the risk analyses to date. To address this gap in the model, concerted studies have been identified to characterize the dwelling mounds. These include new shear-wave velocity measurements that have enabled dynamic site response analyses to determine the modification of ground shaking due to the presence of the mound. A scheme has then been developed to incorporate the dwelling mounds into the risk calculations, which included an assessment of whether the soil-structure interaction effects for buildings founded on the mounds required modification of the seismic fragility functions.


Geology ◽  
2020 ◽  
Author(s):  
Berend A. Verberne ◽  
Suzanne J.T. Hangx ◽  
Ronald P.J. Pijnenburg ◽  
Maartje F. Hamers ◽  
Martyn R. Drury ◽  
...  

Europe’s largest gas field, the Groningen field (the Netherlands), is widely known for induced subsidence and seismicity caused by gas pressure depletion and associated compaction of the sandstone reservoir. Whether compaction is elastic or partly inelastic, as implied by recent experiments, is a key factor in forecasting system behavior and seismic hazard. We sought evidence for inelastic deformation through comparative microstructural analysis of unique drill core recovered from the seismogenic center of the field in 2015, 50 yr after gas production started, versus core recovered before production (1965). Quartz grain fracturing, crack healing, and stress-induced Dauphiné twinning are equally developed in the 2015 and 1965 cores, with the only measurable effect of gas production being enhanced microcracking of sparse K-feldspar grains in the 2015 core. Interpreting these grains as strain markers, we suggest that reservoir compaction involves elastic strain plus inelastic compression of weak clay films within grain contacts.


2017 ◽  
Vol 15 (9) ◽  
pp. 3555-3580 ◽  
Author(s):  
Pauline P. Kruiver ◽  
Ewoud van Dedem ◽  
Remco Romijn ◽  
Ger de Lange ◽  
Mandy Korff ◽  
...  

2017 ◽  
Vol 96 (5) ◽  
pp. s55-s69 ◽  
Author(s):  
Christopher J. Spiers ◽  
Suzanne J.T. Hangx ◽  
André R. Niemeijer

AbstractThis paper describes a research programme recently initiated at Utrecht University that aims to contribute new, fundamental physical understanding and quantitative descriptions of rock and fault behaviour needed to advance understanding of reservoir compaction and fault behaviour in the context of induced seismicity and subsidence in the Groningen gas field. The NAM-funded programme involves experimental rock and fault mechanics work, microscale observational studies to determine the processes that control reservoir rock deformation and fault slip, modelling and experimental work aimed at establishing upscaling rules between laboratory and field scales, and geomechanical modelling of fault rupture and earthquake generation at the reservoir scale. Here, we focus on describing the programme and its intended contribution to understanding the response of the Groningen field to gas production. The key knowledge gaps that drive the programme are discussed and the approaches employed to address them are highlighted. Some of the first results emerging from the work in progress are also reported briefly and are providing important new insights.


2009 ◽  
Vol 1 (4) ◽  
Author(s):  
Sebastiano Imposa ◽  
Jean-Pierre Fourno ◽  
Rosario Raffaele ◽  
Antonio Scaltrito ◽  
Luciano Scarfi

AbstractA one-dimensional velocity model and station corrections for the Middle-Durance fault zone (south-eastern France) were computed by inverting P-wave arrival times recorded on a local seismic network of 8 stations. A total of 93 local events with a minimum of 6 P-phases, RMS 0.4 s and a maximum gap of 220° were selected. Comparison with previous earthquake locations shows an improvement for the relocated earthquakes. Tests were carried out to verify the robustness of inversion results in order to corroborate the conclusions drawn from our findings. The obtained minimum 1-D velocity model can be used to improve routine earthquake locations and represents a further step toward more detailed seismotectonic studies in this area of south-eastern France.


Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. KS35-KS47 ◽  
Author(s):  
Andreas Wuestefeld ◽  
Sonja M. Greve ◽  
Sven Peter Näsholm ◽  
Volker Oye

A priori assessment of the expected location accuracy of a sensor network is typically done through inversion of the traveltime spatial gradients. This approach assumes that the applied location algorithm successfully recovers the global minimum of the objective function. However, even for accurate and precise phase picks, complexity in the velocity model and limitations in the network layout may inhibit the finding of a global minimum. The location algorithms may end up in a local minimum instead. We compare the location performance of various objective functions and minimization approaches. Although most of the analyzed location approaches mostly lead to good location results, none of the analyzed approaches recovered the correct location for all event locations. This implies that (microseismic) event locations estimates include an additional inherent error, which is linked to the applied location approach. This aspect is often neglected when interpreting event locations. Our site comprises a sensor network of two vertical strings, in a 1D velocity model of the Groningen gas field (The Netherlands), complicated by two thin, very high-velocity layers. For a series of synthetic event locations, we calculate arrival times, add picking errors, and then feed these synthetic picks into a set of different location routines. We also determine a novel way to analyze this approach-dependent location error of the sensor network for the given velocity model and a set of event locations: We compare the distances between a set of assumed event locations to resulting locations in the target region of the subsurface. From the cumulative distribution function of that mislocation distance, we determine the [Formula: see text] and [Formula: see text] confidence distances for each method. This results in scalar values representing the location confidence distances for a given method. In turn, this can be used to easily compare the location capabilities of different sensor layouts.


2010 ◽  
Vol 89 (2) ◽  
pp. 91-107 ◽  
Author(s):  
R. Herber ◽  
J. de Jager

AbstractThe impact of oil and, in particular, gas fields discovered in the Dutch subsurface has been very significant. However, 50 years after the discovery of the giant Groningen gas field the Netherlands has become very mature for exploration of oil and gas, and the gas volume left to be discovered in conventional traps is insignificant compared to what has been found already. The total portfolio of conventional prospects held by the industry contains several 100s of billions of cubic metres (bcm), as reported by the Ministry of Economic Affairs, but many of these prospects are unattractive to drill because of their small size or other geologically unfavourable aspects. Hence, for planning purposes of future national gas production the risk should be taken into account that the size of the conventional portfolio is overestimated. The major E&P companies have reduced their exploration efforts and the number of wells drilled as well as the size and total volume of discovered gas reserves has seen a steady decline over the last 10 years. Some surprises may still be in store and can occasionally add a welcome addition of gas. But the follow-up potential of new play and trapping concepts has been disappointing for many years now, and it is concluded that this is unlikely to be different in the future. Remaining conventional discoveries will mainly be in small near-field targets that as a result of technological advances made in the last few decades can be drilled with high confidence, despite their small volumes.This leaves the so-called unconventional gas (UG) resources for a real and significant increase in the exploration potential of the Netherlands. UG resources occur outside conventional structural or stratigraphic traps in tight (low permeability) rocks and are of regional or sub-regional extent, without well-defined hydrocarbon-water contacts. The potential for Basin Centred Gas, Shale Gas and Coal Bed Methane is reviewed. As, according to present-day technology, development of UG requires very dense drilling at low costs with well spacing of a few 100s of metres, only the onshore potential can be commercial, even in the longer term.Recent geological uplift is a characteristic for all North American commercial UG developments. Uplift helps bringing the resources close to the surface and facilitates development of fractures, which are essential for achieving commercial flow rates. This significantly reduces the area where commercial UG resources may occur in the Netherlands. In addition, sweet spots, where commercial flow rates and ultimate recovery per well can be achieved, represent only a fraction of the total ‘play area’. The UG plays in the Dutch subsurface remain to be proven, and there is still a significant technical risk associated with these plays, on top of the commercial risk. Therefore, despite potentially enormous in-place gas volumes in these unconventional plays, recoverable volumes are much less. If UG resources can be proven and are commercially developable, their cumulative volume potential is estimated by us in the order of a few tens to one or two hundreds bcm of recoverable gas at best. Finally, as UG resources produce at very low rates and require large numbers of wells to develop, the environmental impact in a densely populated country like the Netherlands is enormous, and needs to be seriously considered, already in the exploration phase.In a mature area like the Netherlands, industry focus should be on technology development to reduce risk, increase recovery, reduce cost and minimize surface impact. Cooperation between Operators to build multi-well campaigns is therefore strongly recommended to reduce mobilisation cost. In addition, government incentives should be targeted at the development phase, in order to increase economic attractiveness for difficult reservoirs, both conventional and unconventional. In this way State and industry will both be able to maximize their returns on the remaining potential for gas and oil in the next two to three decades.


2017 ◽  
Vol 96 (5) ◽  
pp. s117-s129 ◽  
Author(s):  
Rob M.H.E. van Eijs ◽  
Onno van der Wal

AbstractNot long after discovery of the Groningen field, gas-production-induced compaction and consequent land subsidence was recognised to be a potential threat to groundwater management in the province of Groningen, in addition to the fact that parts of the province lie below sea level. More recently, NAM's seismological model also pointed to a correlation between reservoir compaction and the observed induced seismicity above the field. In addition to the already existing requirement for accurate subsidence predictions, this demanded a more accurate description of the expected spatial and temporal development of compaction.Since the start of production in 1963, multiple levelling campaigns have gathered a unique set of deformation measurements used to calibrate geomechanical models. In this paper we present a methodology to model compaction and subsidence, combining results from rock mechanics experiments and surface deformation measurements. Besides the optical spirit-levelling data, InSAR data are also used for inversion to compaction and calibration of compaction models. Residual analysis, i.e. analysis of the difference between measurement and model output, provides confidence in the model results used for subsidence forecasting and as input to seismological models.


2017 ◽  
Vol 96 (5) ◽  
pp. s175-s182 ◽  
Author(s):  
Stephen J. Bourne ◽  
Stephen J. Oates

AbstractThis paper reviews the evolution of a sequence of seismological models developed and implemented as part of a workflow for Probabilistic Seismic Hazard and Risk Assessment of the seismicity induced by gas production from the Groningen gas field. These are semi-empirical statistical geomechanical models derived from observations of production-induced seismicity, reservoir compaction and structure of the field itself. Initial versions of the seismological model were based on a characterisation of the seismicity in terms of its moment budget. Subsequent versions of the model were formulated in terms of seismic event rates, this change being driven in part by the reduction in variability of the model forecasts in this domain. Our approach makes use of the Epidemic Type After Shock model (ETAS) to characterise spatial and temporal clustering of earthquakes and has been extended to also incorporate the concentration of moment release on pre-existing faults and other reservoir topographic structures.


Author(s):  
Ben Baker ◽  
Monique M. Holt ◽  
Kristine L. Pankow ◽  
Keith D. Koper ◽  
Jamie Farrell

Abstract Immediately following the 18 March 2020 Mww 5.7 Magna, Utah, earthquake, work began on installing a network of three-component, 5 Hz geophones throughout the Salt Lake Valley. After six days, 180 geophones had been sited within 35 km of the epicenter. Each geophone recorded 250 samples per second data onsite for ∼40 days. Here, we integrate the geophone data with data from the permanent regional seismic network operated by the University of Utah Seismograph Stations (UUSS). We use machine learning (ML) methods to create a new catalog of arrival time picks, earthquake locations, and P-wave polarities for 18 March 2020–30 April 2020. We train two deep-learning U-Net models to detect P waves and S waves, assigning arrival times to maximal posterior probabilities, followed by a two-step association process that combines deep learning with a grid-based interferometric approach. Our automated workflow results in 142,000 P picks, 188,000 S picks, and over 5000 earthquake locations. We recovered 95% of the events in the UUSS authoritative catalog and more than doubled the total number of events (5000 vs. 2300). The P and S arrival times generated by our ML models have near-zero biases and standard deviations of 0.05 s and 0.09 s, respectively, relative to corresponding analyst times picked at backbone stations. We also use a deep-learning architecture to automatically determine 70,000 P-wave first motions, which agree with 93% of 5876 hand-picked up or down first motions from both the backbone and nodal stations. Overall, the use of ML led to large increases in the number of arrival times, especially S times, that will be useful for future tomographic studies, as well as the discovery of thousands more earthquakes than exist in the UUSS catalog.


Sign in / Sign up

Export Citation Format

Share Document