H51M-1566
Relocation of Groningen seismicity using refracted waves

Friday, 18 December 2015
Poster Hall (Moscone South)
Jeannot Trampert, Utrecht University, Utrecht, 3584, Netherlands, Elmer Ruigrok, Utrecht University, Utrecht, Netherlands, Hanneke Paulssen, Utrecht University, Dept. of Earth Sciences, Utrecht, Netherlands and Bernard Dost, Royal Netherlands Meteorological Institute, De Bilt, 3730, Netherlands
Abstract:
The Groningen gas field is a giant natural gas accumulation in the Northeast of the Netherlands. The gas is in a reservoir at a depth of about 3 km. The naturally-fractured gas-filled sandstone extends roughly 45 by 25 km laterally and 140 m vertically. Decades of production have led to significant compaction of the sandstone. The (differential) compaction is thought to have reactivated existing faults and being the main driver of induced seismicity.

Precise earthquake location is difficult due to a complicated subsurface, and that is the likely reason, the current hypocentre estimates do not clearly correlate with the well-known fault network. The seismic velocity model down to reservoir depth is quite well known from extensive seismic surveys and borehole data. Most to date earthquake detections, however, were made with a sparse pre-2015 seismic network. For shallow seismicity (<5 km depth) horizontal source-receiver distances tend to be much larger than vertical distances. Consequently, preferred source-receiver travel paths are refractions over high-velocity layers below the reservoir. However, the seismic velocities of layers below the reservoir are poorly known.

We estimated an effective velocity model of the main refracting layer below the reservoir and use this for relocating past seismicity. We took advantage of vertical-borehole recordings for estimating precise P-wave (refraction) onset times and used a tomographic approach to find the laterally varying velocity field of the refracting layer. This refracting layer is then added to the known velocity model, and the combined model is used to relocate the past seismicity. From the resulting relocations we assess which of the faults are being reactivated.