Bob Paap
Bob Paap has a total of 92 co-authors. Here are the top co-authors by number of shared publications:
- Brecht B.T. Wassing
Applied Geosciences, Energy & Materials Transition, TNO, Princetonlaan 6, 3584 CB Utrecht, The Netherlands
2 shared publications Dirk Kraaijpoel
TNO, Applied Geosciences, Utrecht, The Netherlands
2 shared publicationsGuido C.N. van Yperen
2 shared publicationsHans Veldkamp
2 shared publicationsHolger Cremer
2 shared publications- Jan ter Heege
TNO -Geological Survey of the Netherlands Princetonlaan 6 3584 CB Utrecht the Netherlands
2 shared publications
Latest Publications
- Open Access
The UK and Ireland Geophysical Array -- Concept and Design
Scientific exploration of the UK and Ireland's subsurface has made important contributions to scholarship and prosperity for people and the planet, including economic growth, sustainable use of natural resources, storage of greenhouse gases, and inspiring curiosity about the Earth beneath our feet. This article outlines a vision for an array of seismological instruments spanning the UK and Ireland, UKI Array, augmented by other types of geophysical sensors, to maximise the value offered by existing equipment pools. The mission is to research natural phenomena and structure in the deep and shallow Earth, to solve problems concerning hazards and resources, to connect scientists to schools and the broader public, and thus to inspire a new generation to learn about geophysics. The vision was created through a community driven process of engagement and participation. This paper describes the concept and design of the UKI-Array; a companion paper discusses related opportunities and potential applications.
Open AccessThe UK and Ireland Geophysical Array -- Concept and Design
Scientific exploration of the UK and Ireland's subsurface has made important contributions to scholarship and prosperity for people and the planet, including economic growth, sustainable use of natural resources, storage of greenhouse gases, and inspiring curiosity about the Earth beneath our feet. This article outlines a vision for an array of seismological instruments spanning the UK and Ireland, UKI Array, augmented by other types of geophysical sensors, to maximise the value offered by existing equipment pools. The mission is to research natural phenomena and structure in the deep and shallow Earth, to solve problems concerning hazards and resources, to connect scientists to schools and the broader public, and thus to inspire a new generation to learn about geophysics. The vision was created through a community driven process of engagement and participation. This paper describes the concept and design of the UKI-Array; a companion paper discusses related opportunities and potential applications.
- Open Access
- 15 January 2025
- Published by Elsevier in Applied Ocean Research
- Vol. 154
Leveraging Distributed Acoustic Sensing for monitoring vessels using submarine fiber-optic cables
- 15 January 2025
- Published by Elsevier in Applied Ocean Research
- Vol. 154
Open AccessLeveraging Distributed Acoustic Sensing for monitoring vessels using submarine fiber-optic cables
- 5 January 2025
- Published by Society of Exploration Geophysicists in Geophysics
- Vol. 90 (3) , L75-L84
Distributed acoustic sensing, transfer functions, and moment magnitude determination: Application to high-frequency downhole data
Distributed acoustic sensing (DAS) enables sampling seismic wavefields along optical fibers at a spatial resolution of less than 1 m, over distances beyond several tens of kilometers. This makes DAS a powerful tool to record seismic events densely along 2D directions, whether horizontally along the earth’s surface or vertically in boreholes. Compared with traditional seismic sensors measuring ground motion units, DAS provides uniaxial strain measurements along the fiber with often imperfectly known transfer functions between the measurements and true ground motion. This can generate uncertainties in the derivation of seismic source parameters, such as the magnitude, that require an absolute measurement of the ground motion and a known instrument response. In this study, we examine the DAS transfer function, mapping DAS data to reference velocity records obtained from multiple colocated accelerometers. Our investigation makes use of downhole recordings from the Frontier Observatory for Research in Geothermal Energy field site situated in Utah, USA. Overall, we find that the DAS response estimated at different depth positions follows a consistent trend and deviates significantly from a flat response only below 80 Hz. An average site-specific DAS system response is then used to convert microseismic event recordings into calibrated velocity records with improved amplitude accuracies. Subsequently, moment magnitudes MW are derived from the P-wave records with results matching the independent accelerometer-based estimations with high fidelity for events with MW>−1.0.
- 5 January 2025
- Published by Society of Exploration Geophysicists in Geophysics
- Vol. 90 (3) , L75-L84
Distributed acoustic sensing, transfer functions, and moment magnitude determination: Application to high-frequency downhole data
Distributed acoustic sensing (DAS) enables sampling seismic wavefields along optical fibers at a spatial resolution of less than 1 m, over distances beyond several tens of kilometers. This makes DAS a powerful tool to record seismic events densely along 2D directions, whether horizontally along the earth’s surface or vertically in boreholes. Compared with traditional seismic sensors measuring ground motion units, DAS provides uniaxial strain measurements along the fiber with often imperfectly known transfer functions between the measurements and true ground motion. This can generate uncertainties in the derivation of seismic source parameters, such as the magnitude, that require an absolute measurement of the ground motion and a known instrument response. In this study, we examine the DAS transfer function, mapping DAS data to reference velocity records obtained from multiple colocated accelerometers. Our investigation makes use of downhole recordings from the Frontier Observatory for Research in Geothermal Energy field site situated in Utah, USA. Overall, we find that the DAS response estimated at different depth positions follows a consistent trend and deviates significantly from a flat response only below 80 Hz. An average site-specific DAS system response is then used to convert microseismic event recordings into calibrated velocity records with improved amplitude accuracies. Subsequently, moment magnitudes MW are derived from the P-wave records with results matching the independent accelerometer-based estimations with high fidelity for events with MW>−1.0.
- 1 January 2022
- Published by EAGE Publications bv
- Vol. 2022 (1) , 1-5
Modeling Approach for Evaluating Time-Lapse Effects of CO2 Storage on Particle Velocity and Strain Rate Data
- 1 January 2022
- Published by EAGE Publications bv
- Vol. 2022 (1) , 1-5
Modeling Approach for Evaluating Time-Lapse Effects of CO2 Storage on Particle Velocity and Strain Rate Data
- Open Access
A Neural Network Approach for Improved Seismic Event Detection in the Groningen Gas Field, The Netherlands
Over the past decades, the Groningen Gas Field (GGF) has been increasingly faced by induced earthquakes resulting from gas production. The seismic monitoring network at Groningen has been recently densified to improve the seismic network performance, resulting in increasing amounts of seismic data. Although traditional automated event detection techniques generally are successful in detecting events from continuous data, its detection success is challenged in cases of lower signal-to-noise ratios. The data stream coming from these networks has initiated specific interest in neural networks for automated classification and interpretation. Here, we explore the feasibility of neural networks in detecting the occurrence of seismic events. For this purpose, a three-layered feedforward neural network was trained using public data of a seismic event in the GGF obtained from the Royal Netherlands Meteorological Institute (KNMI) data portal. The first arrival times and duration of earthquake waveforms determined by KNMI for a subset of the station data, were used to detect the arrival times and event duration for the other uninterpreted station data. Subsequently various attributes were used as input for the neural network, that were based on different short term averaging/long term averaging (STA/LTA) and frequency sub-band settings. Using these input data, the network's parameters were iteratively improved to maximize its capability in successfully discriminating seismic events from noise and determine the event duration. Results show an increase of 65 % in accurately detecting seismic events and determining their duration as compared to the reference method. This clears the way for improved interpretation of signal waveforms and automated seismic event classification in the Groningen area.
Open AccessA Neural Network Approach for Improved Seismic Event Detection in the Groningen Gas Field, The Netherlands
Over the past decades, the Groningen Gas Field (GGF) has been increasingly faced by induced earthquakes resulting from gas production. The seismic monitoring network at Groningen has been recently densified to improve the seismic network performance, resulting in increasing amounts of seismic data. Although traditional automated event detection techniques generally are successful in detecting events from continuous data, its detection success is challenged in cases of lower signal-to-noise ratios. The data stream coming from these networks has initiated specific interest in neural networks for automated classification and interpretation. Here, we explore the feasibility of neural networks in detecting the occurrence of seismic events. For this purpose, a three-layered feedforward neural network was trained using public data of a seismic event in the GGF obtained from the Royal Netherlands Meteorological Institute (KNMI) data portal. The first arrival times and duration of earthquake waveforms determined by KNMI for a subset of the station data, were used to detect the arrival times and event duration for the other uninterpreted station data. Subsequently various attributes were used as input for the neural network, that were based on different short term averaging/long term averaging (STA/LTA) and frequency sub-band settings. Using these input data, the network's parameters were iteratively improved to maximize its capability in successfully discriminating seismic events from noise and determine the event duration. Results show an increase of 65 % in accurately detecting seismic events and determining their duration as compared to the reference method. This clears the way for improved interpretation of signal waveforms and automated seismic event classification in the Groningen area.
- 11 November 2019
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 220 (2) , 1284-1299
Simulation of induced seismic ground motions using coupled geomechanical and seismic wave propagation models
Summary: Numerical simulations of seismic wave propagation usually rely on a simple source model consisting of an idealized point location and a moment tensor. In general, this is a valid approximation when the source dimensions are small relative to the distance of points at which the seismic wave motions are to be evaluated. Otherwise, a more realistic spatio-temporal source representation is required to accurately calculate ground motions at the position of monitoring stations. Here, we present a generic approach to couple geomechanical simulations to seismic wave propagation models using the concept of the equivalent force field. This approach allows the simulation of seismic wave propagation resulting from the spatio-temporal dependent earthquake nucleation and rupture processes. Within the geomechanical package two separate geomechanics codes are used to simulate both the slow loading stage leading to earthquake nucleation as well as the successive dynamic rupture stage. We demonstrate the approach to a case of induced seismicity, where fault reactivation occurs due to production from a natural gas reservoir.
- 11 November 2019
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 220 (2) , 1284-1299
Simulation of induced seismic ground motions using coupled geomechanical and seismic wave propagation models
Summary: Numerical simulations of seismic wave propagation usually rely on a simple source model consisting of an idealized point location and a moment tensor. In general, this is a valid approximation when the source dimensions are small relative to the distance of points at which the seismic wave motions are to be evaluated. Otherwise, a more realistic spatio-temporal source representation is required to accurately calculate ground motions at the position of monitoring stations. Here, we present a generic approach to couple geomechanical simulations to seismic wave propagation models using the concept of the equivalent force field. This approach allows the simulation of seismic wave propagation resulting from the spatio-temporal dependent earthquake nucleation and rupture processes. Within the geomechanical package two separate geomechanics codes are used to simulate both the slow loading stage leading to earthquake nucleation as well as the successive dynamic rupture stage. We demonstrate the approach to a case of induced seismicity, where fault reactivation occurs due to production from a natural gas reservoir.
- Open Access
- 1 January 2019
- Published by Netherlands Journal of Geosciences Foundation in Netherlands Journal of Geosciences
- Vol. 98
Review of induced seismicity in geothermal systems worldwide and implications for geothermal systems in the Netherlands
Geothermal energy is a viable alternative to gas for the heating of buildings, industrial areas and greenhouses, and can thus play an important role in making the transition to sustainable energy in the Netherlands. Heat is currently produced from the Dutch subsurface through circulation of water between two wells in deep (1.5–3 km) geothermal formations with temperature of up to ∼100 °C. As the number of these so-called doublets is expected to increase significantly over the next decades, and targeted depths and temperatures increase, it is important to assess potential show-stoppers related to geothermal operations. One of these potential hazards is the possibility of the occurrence of felt seismic events, which could potentially damage infrastructure and housing, and affect public support. Such events have been observed in several geothermal systems in other countries. Here we review the occurrence (or the lack) of felt seismic events in geothermal systems worldwide and identify key factors influencing the occurrence and magnitude of these events. Based on this review, we project the findings for seismicity in geothermal systems to typical geothermal formations and future geothermal developments in the Netherlands. The case study review shows that doublets that circulate fluids through relatively shallow, porous, sedimentary aquifers far from the crystalline basement are unlikely to generate felt seismic events. On the other hand, stimulations or circulations in or near competent, fractured, basement rocks and production and reinjection operations in high-temperature geothermal fields are more prone to induce felt events, occasionally with magnitudes ofM> 5.0. Many of these operations are situated in tectonically active areas, and stress and temperature changes may be large. The presence of large, optimally oriented and critically stressed faults increases the potential for induced seismicity. The insights from the case study review suggest that the potential for the occurrence ofM> 2.0 seismicity for geothermal operations in several of the sandstone target formations in the Netherlands is low, especially if faults can be avoided. The potential for induced seismicity may be moderate for operations in faulted carbonate rocks. Induced seismicity always remains a complex and site-specific process with large unknowns, and can never be excluded entirely. However, assessing the potential for inducing felt seismic events can be improved by considering the relevant (site-specific) geological and operational key factors discussed in this article.
- 1 January 2019
- Published by Netherlands Journal of Geosciences Foundation in Netherlands Journal of Geosciences
- Vol. 98
Open AccessReview of induced seismicity in geothermal systems worldwide and implications for geothermal systems in the Netherlands
Geothermal energy is a viable alternative to gas for the heating of buildings, industrial areas and greenhouses, and can thus play an important role in making the transition to sustainable energy in the Netherlands. Heat is currently produced from the Dutch subsurface through circulation of water between two wells in deep (1.5–3 km) geothermal formations with temperature of up to ∼100 °C. As the number of these so-called doublets is expected to increase significantly over the next decades, and targeted depths and temperatures increase, it is important to assess potential show-stoppers related to geothermal operations. One of these potential hazards is the possibility of the occurrence of felt seismic events, which could potentially damage infrastructure and housing, and affect public support. Such events have been observed in several geothermal systems in other countries. Here we review the occurrence (or the lack) of felt seismic events in geothermal systems worldwide and identify key factors influencing the occurrence and magnitude of these events. Based on this review, we project the findings for seismicity in geothermal systems to typical geothermal formations and future geothermal developments in the Netherlands. The case study review shows that doublets that circulate fluids through relatively shallow, porous, sedimentary aquifers far from the crystalline basement are unlikely to generate felt seismic events. On the other hand, stimulations or circulations in or near competent, fractured, basement rocks and production and reinjection operations in high-temperature geothermal fields are more prone to induce felt events, occasionally with magnitudes ofM> 5.0. Many of these operations are situated in tectonically active areas, and stress and temperature changes may be large. The presence of large, optimally oriented and critically stressed faults increases the potential for induced seismicity. The insights from the case study review suggest that the potential for the occurrence ofM> 2.0 seismicity for geothermal operations in several of the sandstone target formations in the Netherlands is low, especially if faults can be avoided. The potential for induced seismicity may be moderate for operations in faulted carbonate rocks. Induced seismicity always remains a complex and site-specific process with large unknowns, and can never be excluded entirely. However, assessing the potential for inducing felt seismic events can be improved by considering the relevant (site-specific) geological and operational key factors discussed in this article.
- 7 June 2018
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 214 (3) , 1947-1960
Wave propagation modelling of induced earthquakes at the Groningen gas production site
Gas extraction from the Groningen natural gas field, situated in the Netherlands, frequently induces earthquakes in the reservoir that cause damage to buildings and pose a safety hazard and a nuisance to the local population. Due to the dependence of the national heating infrastructure on Groningen gas, the short-term mitigation measures are mostly limited to a combination of spatiotemporal redistribution of gas production and strengthening measures for buildings. All options become more effective with a better understanding of both source processes and seismic wave propagation. Detailed wave propagation simulations improve both the inference of source processes from observed ground motions and the forecast of ground motions as input for hazard studies and seismic network design. The velocity structure at the Groningen site is relatively complex, including both deep high-velocity and shallow low-velocity deposits showing significant thickness variations over relatively small spatial extents. We performed a detailed three-dimensional wave propagation modelling study for an induced earthquake in the Groningen natural gas field using the spectral-element method. We considered an earthquake that nucleated along a normal fault with local magnitude of $${{\rm{M}}_{\rm{L}}} = 3$$. We created a dense mesh with element size varying from 12 to 96 m, and used a source frequency of 7 Hz, such that frequencies generated during the simulation were accurately sampled up to 10 Hz. The velocity/density model is constructed using a three-dimensional geological model of the area, including both deep high-velocity salt deposits overlying the source region and shallow low-velocity sediments present in a deep but narrow tunnel valley. The results show that the three-dimensional density/velocity model in the Groningen area clearly play a large role in the wave propagation and resulting surface ground motions. The 3d structure results in significant lateral variations in site response. The high-velocity salt deposits have a dispersive effect on the radiated wavefield, reducing the seismic energy reaching the surface near the epicentre. In turn, the presence of low-velocity tunnel valley deposits can locally cause a significant increase in peak ground acceleration. Here we study induced seismicity on a local scale and use SPECFEM3D to conduct full waveform simulations and show how local velocity variations can affect seismic records.- 7 June 2018
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 214 (3) , 1947-1960
Wave propagation modelling of induced earthquakes at the Groningen gas production site
Gas extraction from the Groningen natural gas field, situated in the Netherlands, frequently induces earthquakes in the reservoir that cause damage to buildings and pose a safety hazard and a nuisance to the local population. Due to the dependence of the national heating infrastructure on Groningen gas, the short-term mitigation measures are mostly limited to a combination of spatiotemporal redistribution of gas production and strengthening measures for buildings. All options become more effective with a better understanding of both source processes and seismic wave propagation. Detailed wave propagation simulations improve both the inference of source processes from observed ground motions and the forecast of ground motions as input for hazard studies and seismic network design. The velocity structure at the Groningen site is relatively complex, including both deep high-velocity and shallow low-velocity deposits showing significant thickness variations over relatively small spatial extents. We performed a detailed three-dimensional wave propagation modelling study for an induced earthquake in the Groningen natural gas field using the spectral-element method. We considered an earthquake that nucleated along a normal fault with local magnitude of $${{\rm{M}}_{\rm{L}}} = 3$$. We created a dense mesh with element size varying from 12 to 96 m, and used a source frequency of 7 Hz, such that frequencies generated during the simulation were accurately sampled up to 10 Hz. The velocity/density model is constructed using a three-dimensional geological model of the area, including both deep high-velocity salt deposits overlying the source region and shallow low-velocity sediments present in a deep but narrow tunnel valley. The results show that the three-dimensional density/velocity model in the Groningen area clearly play a large role in the wave propagation and resulting surface ground motions. The 3d structure results in significant lateral variations in site response. The high-velocity salt deposits have a dispersive effect on the radiated wavefield, reducing the seismic energy reaching the surface near the epicentre. In turn, the presence of low-velocity tunnel valley deposits can locally cause a significant increase in peak ground acceleration. Here we study induced seismicity on a local scale and use SPECFEM3D to conduct full waveform simulations and show how local velocity variations can affect seismic records.- 1 March 2016
- Published by Society of Exploration Geophysicists in Geophysics
- Vol. 81 (2) , KS61-KS70
Calibration of a local magnitude relationship for microseismic events using earthquake data
Calibration of local seismic arrays is crucial to quantify the magnitude of local microseismic events. This is traditionally done by performing a site survey consisting of calibration shots at known locations with known strength. However, when a calibration survey cannot be conducted, alternative calibration methods are required. We have demonstrated the feasibility of magnitude estimation of microseismic events recorded at the storage site at Ketzin (Germany), by analyzing earthquake signals that have been recorded on the noncalibrated geophone array and on standardized seismometers of a regional seismological network. We estimated a linear calibration relation by fitting three different signal attributes for 10 distant earthquakes. Through extrapolation of the linear fit toward lower magnitudes, we estimated the local magnitude of weak local seismic events. Reported magnitudes of the distant earthquakes show significant variation because there are different methods and constants used to calculate a magnitude estimate from recordings on multiple seismometers. As a consequence, there was a considerable spread in our magnitude estimates for the Ketzin events. We found magnitude values in the range between and for two local events recorded at the Ketzin site. We expect that by incorporating additional distant earthquake data, the uncertainty in this estimate can be further reduced.- 1 March 2016
- Published by Society of Exploration Geophysicists in Geophysics
- Vol. 81 (2) , KS61-KS70
Calibration of a local magnitude relationship for microseismic events using earthquake data
Calibration of local seismic arrays is crucial to quantify the magnitude of local microseismic events. This is traditionally done by performing a site survey consisting of calibration shots at known locations with known strength. However, when a calibration survey cannot be conducted, alternative calibration methods are required. We have demonstrated the feasibility of magnitude estimation of microseismic events recorded at the storage site at Ketzin (Germany), by analyzing earthquake signals that have been recorded on the noncalibrated geophone array and on standardized seismometers of a regional seismological network. We estimated a linear calibration relation by fitting three different signal attributes for 10 distant earthquakes. Through extrapolation of the linear fit toward lower magnitudes, we estimated the local magnitude of weak local seismic events. Reported magnitudes of the distant earthquakes show significant variation because there are different methods and constants used to calculate a magnitude estimate from recordings on multiple seismometers. As a consequence, there was a considerable spread in our magnitude estimates for the Ketzin events. We found magnitude values in the range between and for two local events recorded at the Ketzin site. We expect that by incorporating additional distant earthquake data, the uncertainty in this estimate can be further reduced.