The Wolfspar experience with low-frequency seismic source field data: Motivation, processing, and implications

The promise of fully automatic full-waveform inversion (FWI) — a (seismic) data-driven velocity model building process — has proven elusive in complex geologic settings, with impactful examples using field data unavailable until recently. In 2015, success with FWI at the Atlantis Field in the U.S. Gulf of Mexico demonstrated that semiautomatic velocity model building is possible, but it also raised the question of what more might be possible if seismic data tailor-made for FWI were available (e.g., with increased source-receiver offsets and bespoke low-frequency seismic sources). Motivated by the initial value case for FWI in settings such as the Gulf of Mexico, beginning in 2007 and continuing into 2021 BP designed, built, and field tested Wolfspar, an ultralow-frequency seismic source designed to produce seismic data tailor-made for FWI. A 3D field trial of Wolfspar was conducted over the Mad Dog Field in the Gulf of Mexico in 2017–2018. Low-frequency source (LFS) data were shot on a sparse grid (280 m inline, 2 to 4 km crossline) and recorded into ocean-bottom nodes simultaneously with air gun sources shooting on a conventional dense grid (50 m inline, 50 m crossline). Using the LFS data with FWI to improve the velocity model for imaging produced only incremental uplift in the subsalt image of the reservoir, albeit with image improvements at depths greater than 25,000 ft (approximately 7620 m). To better understand this, reprocessing and further analyses were conducted. We found that (1) the LFS achieved its design signal-to-noise ratio (S/N) goals over its frequency range; (2) the wave-extrapolation and imaging operators built into FWI and migration are very effective at suppressing low-frequency noise, so that densely sampled air gun data with a low S/N can still produce useable model updates with low frequencies; and (3) data density becomes less important at wider offsets. These results may have significant implications for future acquisition designs with low-frequency seismic sources going forward.

Seismic hazard analyses of the Metropolitan Water District Emergency Freshwater Pathway, California

The Sacramento-San Joaquin Delta in central California is particularly susceptible to damage in a large earthquake due to the vulnerability of the levees that protect cities, farms, and infrastructure. The Delta is located adjacent to the seismically active San Andreas fault system and is also subject to strong ground shaking from numerous other seismic sources in central California, including faults within the Delta. We performed a probabilistic seismic hazard analysis (PSHA) to provide seismic design ground motions for the Metropolitan Water District (MWD) Emergency Freshwater Pathway. We have evaluated the appropriateness of the Next Generation of Attenuation (NGA)-West2 ground motion models for use in our analyses of the Delta, evaluated shear-wave velocity ( VS) data in the vicinity of the Pathway, and performed site response analyses. The latter was performed to compute the probabilistic hazard at the top of the peat at five sites along the Pathway. The sixth site was located outside the Delta and on firm soil. The probabilistic hazard for the six sites and for a range of return periods of engineering relevance were computed in the PSHA. For a return period of 2500 years, the peak horizontal ground acceleration (PGA) values on peat ranged from 0.40 g to 0.53 g. The seismic sources that control the hazard at these sites vary as a function of return period and spectral frequency, but in general, the closer the sites are to faults within the San Andreas fault system, the higher the hazard.

Characterization of the seismic field at Virgo and improved estimates of Newtonian-noise suppression by recesses

Fluctuations of gravitational forces cause so-called Newtonian noise (NN) in gravitational-wave (GW) detectors which is expected to limit their low-frequency sensitivity in upcoming observing runs. Seismic NN is produced by seismic waves passing near a detector's suspended test masses. It is predicted to be the strongest contribution to NN. Modeling this contribution accurately is a major challenge. Arrays of seismometers were deployed at the Virgo site to characterize the seismic field near the four test masses. In this paper, we present results of a spectral analysis of the array data from one of Virgo's end buildings to identify dominant modes of the seismic field. Some of the modes can be associated with known seismic sources. Analyzing the modes over a range of frequencies, we provide a dispersion curve of Rayleigh waves. We find that the Rayleigh speed in the NN frequency band 10\,Hz--20\,Hz is very low ($\lesssim$100\,m/s), which has important consequences for Virgo's seismic NN. Using the new speed estimate, we find that the recess formed under the suspended test masses by a basement level at the end buildings leads to a 10 fold reduction of seismic NN.

Monitoring subsurface changes by tracking direct-wave amplitudes and traveltimes in continuous DAS VSP data

Instrumenting wells with distributed acoustic sensors (DAS) and illuminating them with passive or active seismic sources allows precise tracking of temporal variations of direct-wave traveltimes and amplitudes, which can be used to monitor variations in formation stiffness and density. This approach has been tested by tracking direct-wave amplitudes and traveltimes as part of a CCS project where a 15 kt supercritical CO2 injection was monitored with continuous offset VSPs using nine permanently mounted surface orbital vibrators (SOVs) acting as seismic sources and several wells instrumented with DAS cables cemented behind the casing. The results show a significant (from 15 to 30%) increase of strain amplitudes within the CO2 injection interval, and travetime shifts of 0.3 to 0.4 ms below this interval, consistent with full-wave 1.5D numerical simulations and theoretical predictions. The results give independent estimates of the CO2 plume thickness and P-wave velocity reduction within it.

Permanent, seasonal, and episodic seismic sources around Vatnajökull, Iceland, from the analysis of correlograms

In this study, we locate and characterise the main seismic noise sources in the region of the Vatnajökull icecap (Iceland). Vatnajökull is the largest Icelandic icecap, covering several active volcanoes. The seismic context is very complex, with glacial and volcanic events occurring simultaneously and the classification between the two can become cumbersome. Using seismic interferometry on continuous seismic data (2011–2019), we calculate the propagation velocities and locate the main seismic sources by using hyperbolic geometry and a grid-search method. We identify and characterise permanent oceanic sources, seasonal glacial-related sources, and episodic volcanic sources. These results give a better understanding of the background seismic noise sources in this region and could allow the identification of seismic sources associated with potentially threatening events in real-time.

Seismic and Geodetic Crustal Moment-Rates Comparison: New Insights on the Seismic Hazard of Egypt

A comparative analysis of geodetic versus seismic moment-rate estimations makes it possible to distinguish between seismic and aseismic deformation, define the style of deformation, and also to reveal potential seismic gaps. This analysis has been performed for Egypt where the present-day tectonics and seismicity result from the long-lasting interaction between the Nubian, Eurasian, and Arabian plates. The data used comprises all available geological and tectonic information, an updated Poissonian earthquake catalog (2200 B.C.–2020 A.D.) including historical and instrumental datasets, a focal-mechanism solutions catalog (1951–2019), and crustal geodetic strains from Global Navigation Satellite System (GNSS) data. The studied region was divided into ten (EG-01 to EG-10) crustal seismic sources based mainly on seismicity, focal mechanisms, and geodetic strain characteristics. The delimited seismic sources cover the Gulf of Aqaba–Dead Sea Transform Fault system, the Gulf of Suez­–Red Sea Rift, besides some potential seismic active regions along the Nile River and its delta. For each seismic source, the estimation of seismic and geodetic moment-rates has been performed. Although the obtained results cannot be considered to be definitive, among the delimited sources, four of them (EG-05, EG-06, EG-08, and EG-10) are characterized by low seismic-geodetic moment-rate ratios (<20%), reflecting a prevailing aseismic behavior. Intermediate moment-rate ratios (from 20% to 60%) have been obtained in four additional zones (EG-01, EG-04, EG-07, and EG-09), evidencing how the seismicity accounts for a minor to a moderate fraction of the total deformational budget. In the other two sources (EG-02 and EG-03), high seismic-geodetic moment-rates ratios (>60%) have been observed, reflecting a fully seismic deformation.