EGU22-5398, updated on 27 Mar 2022
https://doi.org/10.5194/egusphere-egu22-5398
EGU General Assembly 2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.

Characterizing induced seismic events in the Groningen gas field using an efficient Hamiltonian Monte Carlo sampler: a case study

La Ode Marzujriban Masfara1, Cornelis Weemstra1,2, and Thomas Cullison3
La Ode Marzujriban Masfara et al.
  • 1Delft University of Technology, Delft, The Netherlands (l.o.m.masfara@tudelft.nl)
  • 2Royal Netherlands Meteorological Institute, Utrecht, The Netherlands (c.weemstra@tudelft.nl)
  • 3Utrecht University, Utrecht, The Netherlands (t.a.cullison@uu.nl)

In May 2019, an earthquake with a magnitude of 3.4 (local magnitude) hit the area of the Westerwijtwerd village in the province of Groningen, the Netherlands. The event is the result of the gas extraction in the Groningen gas field and is one of the largest events to date. To better understand the source characteristics of the event, we apply a probabilistic full-waveform inversion technique that we recently developed to the event's recordings. Specifically, we use a variant of the Hamiltonian Monte Carlo (HMC) algorithm. When sampling high-dimensional model spaces, HMC is proven to be more efficient than the generic Metropolis-Hasting algorithm. Compared to probabilistic inversions of tectonic events, two main challenges arise while applying the algorithm. First, the prior information of the event is usually incomplete and inaccurate. That is, the only available information is (an estimate of) the hypocenter and origin time. Second, the frequency content of the induced event's seismograms is higher than that of typical tectonic events. This implies a higher non-linearity, which in turn complicates the ability of a probabilistic inversion algorithm to sample the model spaces, particularly when considering the first challenge. Consequently, to address both challenges, first, we develop a procedure to estimate the necessary prior information and use it as input to the HMC variant. Second, we run our HMC algorithm iteratively to mitigate the non-linearity. Using the relatively detailed velocity model of the Groningen gas field, we eventually estimate ten posteriors of the source parameters. The latter being the hypocenter (three parameters), the moment tensor (six independent parameters), and the origin time.  

How to cite: Masfara, L. O. M., Weemstra, C., and Cullison, T.: Characterizing induced seismic events in the Groningen gas field using an efficient Hamiltonian Monte Carlo sampler: a case study, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5398, https://doi.org/10.5194/egusphere-egu22-5398, 2022.

Displays

Display file