An efficient probabilistic workflow for estimating induced earthquake parameters in 3D heterogeneous media

We present an efficient probabilistic workflow for the estimation of source parameters of induced seismic events in three-dimensional heterogeneous media. Our workflow exploits a linearized variant of the Hamiltonian Monte Carlo (HMC) algorithm. Compared to traditional Markov chain Monte Carlo (MCMC) algorithms, HMC is highly efficient in sampling high-dimensional model spaces. Through a linearization of the forward problem around the prior mean (i.e., the “best” initial model), this efficiency can be further improved. We show, however, that this linearization leads to a performance in which the output of an HMC chain strongly depends on the quality of the prior, in particular because not all (induced) earthquake model parameters have a linear relationship with the recordings observed at the surface. To mitigate the importance of an accurate prior, we integrate the linearized HMC scheme into a workflow that (i) allows for a weak prior through linearization around various (initial) centroid locations, (ii) is able to converge to the mode containing the model with the (global) minimum misfit by means of an iterative HMC approach, and (iii) uses variance reduction as a criterion to include the output of individual Markov chains in the estimation of the posterior probability. Using a three-dimensional heterogeneous subsurface model of the Groningen gas field, we simulate an induced earthquake to test our workflow. We then demonstrate the virtue of our workflow by estimating the event's centroid (three parameters), moment tensor (six parameters), and the earthquake's origin time. Using the synthetic case, we find that our proposed workflow is able to recover the posterior probability of these source parameters rather well, even when the prior model information is inaccurate, imprecise, or both inaccurate and imprecise.Applied Geophysics and Petrophysic

Improved surface-wave response from ambient noise in Malargüe, Argentina, using seismic interferometry by multidimensional deconvolution

Applied Geophysics and Petrophysic

Seismic monitoring of Nature’s Heat Geothermal Project in Kwintsheul (Netherlands)

In 2018, a geothermal doublet started operating in Kwinstheul, Netherlands, for supplying heat to 64 hectares of greenhouses corresponding to Nature’s Heat joint initiative. This kind of geothermal operation requires extraction, circulation, and reinjection of fluids at a depth of 2.4 km. The reservoir used for the geothermal operation has shown good hydraulic parameters which allow the circulation of the fluid. Several authors agree that this kind of geothermal operation is unlikely to generate felt seismicity, nevertheless, adequate seismic monitoring is critical to guarantee sustainable and safe use of the subsurface. To monitor the operation of Nature’s Heat project, 30 three-component short-period seismic sensors were installed by Delft University of Technology and Seismotech (Greece). A challenge for seismic monitoring in Kwinstheul is the high levels of seismic noise coming from anthropogenic and operational activities. Despite the high background noise levels, a seismic event of Md 0.16 was recorded on July 14, 2019. To understand the relation of the event and improve the safety of the geothermal operation, we are developing an optimized monitoring scheme.Applied Geophysics and Petrophysic

Hydroacoustic Travel Time Variations as a Proxy for Passive Deep-Ocean Thermometry: A Cookbook

We report on the extraction of deep ocean travel time variations from time-lapse cross-correlations between a hydrophone station and a three-component broadband seismometer. The signals we cross-correlate in this study result from repeated activity by the Monowai seamount, one of the most active submarine volcanoes of the Tonga-Kermadec ridge. In particular, we introduce a specific workflow to exploit repetitive hydroacoustic underwater source activity, which we detail to such an extent that it serves as an example (or “cookbook”). For this reason, we have made the source code publicly available. The workflow proposed in this study (a) overcomes differences in instrument sensitivity and sample rates, (b) involves the selection of eligible cross-correlations based on a source activity criterium as well as slowness analysis, and (c) extracts the travel time variations in distinct frequency bands. In our case, the two frequency bands are 3–6 and 6–12 Hz. We find that the estimated travel time variations in both frequency bands consist of a complex periodic pattern superimposed on a robust linear trend. This linear trend is decreasing, which we attribute to increasing water temperatures along the propagation path of the hydroacoustic signals.Applied Geophysics and Petrophysic

On the estimation of attenuation from the ambient seismic field: Inferences from distributions of isotropic point scatterers

International audienceCross-correlation of ambient seismic noise recorded by two seismic stations may result in an estimate of the Green’s function between those two receivers. Several authors have recently attempted to measure attenuation based on these interferometric, receiver–receiver surface waves. By now, however, it is well established that the loss of coherence of the cross-correlation as a function of space depends strongly on the excitation of the medium. In fact,in a homogeneous dissipative medium, uniform excitation is required to correctly recover attenuation. Applied to fundamental-mode ambient seismic surface waves, this implies that the cross-correlation will decay at the local attenuation rate only if noise sources are distributeduniformly on the Earth’s surface. In this study we show that this constraint can be relaxed in case the observed loss of coherence is due to multiple scattering instead of dissipation of energy. We describe the scattering medium as an effective medium whose phase velocity and rate of attenuation are a function of the scatterer density and the average strength of the scatterers. We find that the decay of the cross-correlation in the effective medium coincideswith the local attenuation of the effective medium in case the scattering medium is illuminated uniformly from all angles. Consequently, uniform excitation is not a necessary condition for the correct retrieval of scattering attenuation. We exemplify the implications of this finding for studies using the spectrally whitened cross-correlation to infer subsurface attenuation

Reflecting boundary conditions for interferometry by multidimensional deconvolution: invited paper

Seismic interferometry (SI) takes advantage of existing (ambient) wavefield recordings by turning receivers into so-called “virtual-sources.” The medium’s response to these virtual sources can be harnessed to image that medium. Applications of SI include surface-wave imaging of the Earth’s shallow subsurface and medical imaging. Most interferometric applications, however, suffer from the fact that the retrieved virtual-source responses deviate from the true medium responses. The accrued artifacts are often predominantly due to a non-isotropic illumination of the medium of interest, and prohibit accurate interferometric imaging. Recently, it has been shown that illumination-related artifacts can be removed by means of a so-called multidimensional deconvolution (MDD) process. However, The current MDD formulation, and hence method, relies on separation of waves traveling inward and outward through the boundary of the medium of interest. As a consequence, it is predominantly useful when receivers are illuminated from one side only. This puts constraints on the applicability of the current MDD formulation to omnidirectional wavefields. We present a modification of the formulation of the theory underlying SI by MDD. This modification eliminates the requirement to separate inward-and outward propagating wavefields and, consequently, holds promise for the application of MDD to non-isotropic, omnidirectional wavefieldsApplied Geophysics and PetrophysicsApplied Mechanic

Reflecting boundary conditions for interferometry by multidimensional deconvolution

In this work we investigate a modification of the formulation of the theory underlying seismic interferometry (SI) by multidimensional deconvolution (MDD). The current formulation, and hence method, relies on separation of waves traveling inward and outward of a volume bounded by receivers. As a consequence, it is predominantly useful when receivers are illuminated from one side only. This puts constraints on the applicability of SI by MDD to omnidirectional wave fields. The proposed modification eliminates the requirement to separate inward-and outward propagating wave field and, consequently, improves the applicability of MDD to omnidirectional wave fields. We therefore envisage the modified MDD formulation to hold significant promise in the application to ambient-noise surface wave data.Geoscience & EngineeringCivil Engineering and Geoscience

Towards monitoring the englacial fracture state using virtual-reflector seismology

In seismology, coda wave interferometry (CWI) is an effective tool to monitor time-lapse changes using later arriving, multiply scattered coda waves. Typically, CWI relies on an estimate of the medium's impulse response. The latter is retrieved through simple time-averaging of receiver-receiver cross-correlations of the ambient field, that is, seismic interferometry (SI). In general, the coda is induced by heterogeneities in the Earth. Being comparatively homogeneous, however, ice bodies such as glaciers and ice sheets exhibit little scattering. In addition, the temporal stability of the time-averaged cross-correlations suffers from temporal variations in the distribution and amplitude of the passive seismic sources. Consequently, application of CWI to ice bodies is currently limited. Nevertheless, fracturing and changes in the englacial macroscopic water content alter the bulk elastic properties of ice bodies, which can be monitored with cryoseismological measurements. To overcome the current limited applicability of CWI to ice bodies, we therefore introduce virtual-reflector seismology (VRS). VRS relies on a so-called multidimensional deconvolution (MDD) process of the time-averaged crosscorrelations. The technique results in the retrieval of a medium response that includes virtual reflections from a contour of receivers enclosing the region of interest (i.e. the region to be monitored). The virtual reflections can be interpreted as artificial coda replacing the (lacking) natural scattered coda. Hence, this artificial coda might be exploited for the purpose of CWI. From an implementation point of view, VRS is similar to SI by MDD, which, as its name suggests, also relies on a multidimensional deconvolution process. SI by MDD, however, does not generate additional virtual reflections. Advantageously, both techniques mitigate spurious coda changes associated with temporal variations in the distribution and amplitude of the passive seismic sources. In this work, we apply SI by MDD and VRS to synthetic and active seismic surface-wave data. The active seismic data were acquired on Glacier de la Plaine Morte, Switzerland. We successfully retrieve virtual reflections through the application of VRS to this active seismic data. In application to both synthetic and active seismic data, we show the potential of VRS to monitor time-lapse changes. In addition, we find that SI by MDD allows for a more accurate determination of phase velocity.Applied Geophysics and Petrophysic

On the Potential of 3D Transdimensional Surface Wave Tomography for Geothermal Prospecting of the Reykjanes Peninsula

Seismic travel time tomography using surface waves is an effective tool for three-dimensional crustal imaging. Historically, these surface waves are the result of active seismic sources or earthquakes. More recently, however, surface waves retrieved through the application of seismic interferometry have also been exploited. Conventionally, two-step inversion algorithms are employed to solve the tomographic inverse problem. That is, a first inversion results in frequency-dependent, two-dimensional maps of phase velocity, which then serve as input for a series of independent, one-dimensional frequency-to-depth inversions. As such, a set of localized depth-dependent velocity profiles are obtained at the surface points. Stitching these separate profiles together subsequently yields a three-dimensional velocity model. Relatively recently, a one-step three-dimensional non-linear tomographic algorithm has been proposed. The algorithm is rooted in a Bayesian framework using Markov chains with reversible jumps, and is referred to as transdimensional tomography. Specifically, the three-dimensional velocity field is parameterized by means of a polyhedral Voronoi tessellation. In this study, we investigate the potential of this algorithm for the purpose of recovering the three-dimensional surface-wave-velocity structure from ambient noise recorded on and around the Reykjanes Peninsula, southwest Iceland. To that end, we design a number of synthetic tests that take into account the station configuration of the Reykjanes seismic network. We find that the algorithm is able to recover the 3D velocity structure at various scales in areas where station density is high. In addition, we find that the standard deviation of the recovered velocities is low in those regions. At the same time, the velocity structure is less well recovered in parts of the peninsula sampled by fewer stations. This implies that the algorithm successfully adapts model resolution to the density of rays. It also adapts model resolution to the amount of noise in the travel times. Because the algorithm is computationally demanding, we modify the algorithm such that computational costs are reduced while sufficiently preserving non-linearity. We conclude that the algorithm can now be applied adequately to travel times extracted from station–station cross correlations by the Reykjanes seismic network.Applied Geophysics and Petrophysic

Reflecting boundary conditions for interferometry by multidimensional deconvolution.

Session 1pUWc: Underwater Acoustics: Topics in Underwater Acoustics (Poster Session)Applied Geophysics and PetrophysicsImPhys/Acoustical Wavefield ImagingDynamics of Structure