The overburden structures often can distort the responses of the target region in seismic data, especially in land datasets. Ideally, all effects of the overburden and underburden structures should be removed, leaving only the responses of the target region. This can be achieved using the Marchenko method. The Marchenko method is capable of estimating Green's functions between the surface of the Earth and arbitrary locations in the subsurface. These Green's functions can then be used to redatum wavefields to a level in the subsurface. As a result, the Marchenko method enables the isolation of the response of a specific layer or package of layers, free from the influence of the overburden and underburden. In this study, we apply the Marchenko-based isolation technique to land S-wave seismic data acquired in the Groningen province, the Netherlands. We apply the technique for combined removal of the overburden and underburden, which leaves the isolated response of the target region, which is selected between 30 and 270 m depth. Our results indicate that this approach enhances the resolution of reflection data. These enhanced reflections can be utilised for imaging and monitoring applications.

High-resolution seismic reflections are essential for imaging and monitoring applications. In seismic land surveys using sources and receivers at the surface, surface waves often dominate, masking the reflections. In this study, we demonstrate the efficacy of a two-step procedure to suppress surface waves in an active-source reflection seismic data set. First, we apply seismic interferometry (SI) by cross-correlation, turning receivers into virtual sources to estimate the dominant surface waves. Then, we perform adaptive subtraction to minimize the difference between the surface waves in the original data and the result of SI. We propose a new approach where the initial suppression results are used for further iterations, followed by adaptive subtraction. This technique aims to enhance the efficacy of data-driven surface-wave suppression through an iterative process. We use a 2-D seismic reflection data set from Scheemda, situated in the Groningen province of the Netherlands, to illustrate the technique’s efficiency. A comparison between the data after recursive interferometric surface-wave suppression and the original data across time and frequency–wavenumber domains shows significant suppression of the surface waves, enhancing visualization of the reflections for subsequent subsurface imaging and monitoring studies.

Seismic interferometry (SI) retrieves the Green function between two receiver locations using their recordings from a boundary of sources. When using sources and receivers only at the surface, the virtual-source gathers retrieved by SI contain pseudo-physical reflections as well as ghost (non-physical) reflections. These ghost reflections are the results of the cross-correlation or auto-correlation (AC) of primary reflections from two different depth levels, and they contain information about the seismic properties of specific layers in the subsurface. We investigated the application of ghost reflections for layer-specific characterization of the shallow subsurface using SI by AC. First, we showed the technique's potential using synthetic data for a subsurface model with a lateral change in velocity, a gradient in depth for velocity, a thickness change and a velocity change of the target layer. Then, we applied the technique to shallow subsurface field data. We also focused on improving the retrieval of ghost reflections by removing the free-surface multiples and muting undesired events in active-source gathers before applying SI. Our results demonstrate that the ghost reflections can be used advantageously to characterize the layer that causes them to appear in the results of SI. Consequently, they can also provide valuable information for imaging and monitoring shallow subsurface structures.

The Azambuja fault is a NNE trending structure located 50 km north of Lisbon, the capital and most populous city of Portugal. The fault has been considered as a possible source for the historical, large earthquakes. Understanding this fault is a priority in seismic hazard evaluation of this region. The fault has a clear morphological signature. Miocene and Pliocene sediments are tilted eastward and cut by steeply dipping mesoscale fault segments, presenting reverse and normal offsets with a net downthrow to the east. Neotectonic studies indicate a Quaternary slip on the fault of 0.05–0.06 mm/year. However, no direct evidence of the Azambuja fault affecting the Pleistocene or Holocene sediments was found so far. Here, we present the findings from high-resolution seismic reflection studies using both P- and S-waves over the Holocene deposits. The detection of small-throw faulting in ductile sediments is a challenging task. We show that multiple signatures, like perturbations in the reflection hyperbolae visible in shot and CMP gathers, interruptions of reflectors in stacked sections, lateral seismic velocity variations obtained by horizon velocity analysis, all at coincident locations, strongly suggest that the activity of the Azambuja fault has affected the Holocene sediments in the study area. The lateral velocity variations are corroborated by wavepath eikonal traveltime tomography and velocity analysis supported by seismic modeling. By means of 2D viscoelastic modeling, we explain the absence of fault-related diffractions and negligible back-scattered energy from the fault. Using data from nearby boreholes, we find that the 15 ka old alluvium cover has indeed been disturbed by the presence of shallow fault strands. Considering the estimated vertical throws and the empirical relationships between fault length, co-seismic rupture and magnitude, a slip rate of 0.07 mm/y, slightly larger than previously thought, is expected for this fault.

Using post-critical reflection data, it is possible to obtain useful information that allows more reliable geological characterization of the subsurface. However, the strong distortion caused by the phase shift in post-critical wavelets makes the use of post-critical reflections rather challenging. For this reason, an approach which is capable of estimating the phase shift of each wavelet of a reflection event in a data-driven manner is desirable. In this vein, in case the frequency spectrum of a wavelet can be correctly estimated, it is possible to estimate the instantaneous phase shift. In this work, we propose an approach which can perform such estimation based on spectral recomposition of seismic data. We design an inversion approach in order to reconstruct the seismic spectrum of the wavelets of a reflection event, which subsequently allows us to estimate the instantaneous phase of each wavelet of the near-surface reflection events without performing prior velocity analysis and/or critical-angle estimation. After finding the instantaneous phase for each wavelet of a reflection event, we show next how one can find the respective phase shifts that can then be corrected.

Using the critical angle information of a reflection event, it is possible to calculate several essential physical parameters that are key to reliable geological characterization of the subsurface. However, estimation of the critical angle usually requires several steps of seismic processing. For this reason, an approach which is capable of estimating the critical angle directly from the data is of interest. Once the critical angle is estimated, it is possible to estimate further the Poisson's ratio and the seismic velocities. In this work, we propose an approach which can perform this estimation, based on spectral recomposition of seismic data. We design an inversion scheme in order to reconstruct the seismic spectrum of wavelets of a reflection event, which subsequently allows us to estimate the critical angle of near-surface reflection events without performing prior velocity analysis. After finding the critical angle, we show next how to estimate the Poisson's ratio and the compressional- and shear-wave velocities of the medium above the reflector. The approach leads to quite accurate values for Poisson's ratio even for noisy data, in case the number of layers is not large.

Bulk-density ( ρ) of soil is an important indicator of soil compaction and type. A knowledge of the spatial variability of in situ soil density is important in geotechnical engineering, hydrology and agriculture. Surface geophysical methods have so far shown limited success in providing an accurate and high-resolution image of 3-D soil-density distribution. In this pursuit, 3-D seismic full-waveform inversion (FWI) is promising, provided the robustness and accuracy of density inversion via this approach can be established in the near-surface scale. Ho wever , simultaneous reconstruction of ρand seismic wave velocities through multiparameter FWI remains a challenging task. Near-surface seismic data are commonly dominated by dispersive surface waves whose velocities are controlled by the value and distribution of shear-wave velocity ( V S ). One major difficulty in estimating reliab ly ρfrom near -surface seismic data is due to the relati vel y low sensiti vity of the seismic w av efield to ρcompared to seismic v elocities. Additionally, the accuracy of the estimated ρdecreases due to error in V S -an issue known as parameter coupling. Parameter coupling makes it difficult to estimate accurately ρwithin the framework of conventional gradient-based FWI. More sophisticated optimization approaches (e.g. truncated Newton) can reduce the effect of parameter coupling, but these approaches are commonl y not af fordab le in near -surface applications due to heavy computational burden. In this research, w e ha v e inv estigated how choosing correctly the force direction of the seismic source can contribute to a higher accuracy of ρestimates through 3-D FWI. Using scattered wavefields, the Hessian, and inversion tests, an in-depth and systematic investigation of data sets corresponding to different force directions has been carried out. A comparison of the scattered wavefields due to a point-localized ρperturbation for different force directions shows the robustness of the horizontal-force data set to noise compared to the vertical-force data set. Fur ther more, for a point-scatterer model, an analysis of the gradients of the misfit function using the Hessian shows that utilizing a horizontal-force source enables one to reconstruct the high-resolution gradient with relati vel y small parameter coupling. Finally, inversion tests for two different subsoil models demonstrate that 3-D FWI on a horizontal-force-source seismic data set is capable of providing a more accurate 3-D ρdistribution in soil compared to a vertical-force-source data set. Our results show that the use of a horizontal-force source might allow avoiding computationally demanding, costly optimization approaches in 3-D FWI.

The structural health of historic quay walls needs to be evaluated well in light of the new conditions they are subjected to. For that, information about their current subsurface structure and condition of their subsurface constructional elements, but also information about the surrounding subsurface structure is crucial. Such information can be supplied by seismic imaging and characterization. We show preliminary results from a high-resolution S-wave survey we performed at a historic quay wall in Overamstel, the Netherlands. We recorded data along four lines – two parallel and two perpendicular to the quay wall. We used a sledge-hammer and a beam as a source and 10-Hz horizontal-component geophones, both oriented in the crossline direction. We show that applying simple processing along the two parallel lines to obtain stacked sections already allowed extracting useful structural information of the subsurface.

Knowing the location and characteristics of shallow subsurface structures like tunnels, cavities, archeological ruins, etc. is of importance for different disciplines and application. To image and/or characterize such objects of interest, different geophysical methods are used. For imaging of a very shallow network of tunnels, the high-resolution seismic method with active sources provide valuable information. We show the results of analysis of an S-wave profile recorded over a network of very shallow tunnels in the Netherlands. The survey used a high-frequency vibratory S-wave source and horizontal particle-velocity geophones, both oriented in the crossline direction, along three lines. We process the reflection data along one of the lines to obtain a stacked section in depth. We also use a method inspired by seismic interferometry to localize a scatterer along the line. We show that both techniques image well the subsurface structures taking into account the 3D ambiguity of processing 2D data.

Intersecting faults are often ignored in the geomechanical simulation of induced seismicity. To investigate the effects of fault intersection and the resulting reservoir geometry on induced seismicity, caused, for instance, by gas extraction, we have developed 3D geomechanical models considering two intersecting normal faults and the surrounding horst structure. We simulate the stress field and the dynamic fault reactivation in a uniformly depleted reservoir. We observe that a smaller intersection angle increases the incremental Coulomb stress at the lower reservoir juxtaposition, thus changing the temporal rupture pattern of the seismic event. In our dynamic simulation, the rupture propagates from the main fault to the secondary fault. We conclude that the fault intersection has important effects on the induced seismicity and should be taken into account when evaluating the seismicity risk in a specific region.

Seismic incoherent noise and waves scattered from objects in the crossline directions can cause 2D elastic full-waveform inversion (FWI) to produce artifacts in the resulting 2D models. We develop a complete workflow that can determine subsurface S-wave velocity (VS) models inverted from 2D near-surface seismic data more stably. We make use of a combination of supervirtual interferometry and a matched filter to accurately retrieve dominant surface waves from the field data, whereas the incoherent noise and 3D scattering events are significantly suppressed. The subsurface structures obtained from inverting the retrieved data can be interpreted together with the sections resulting from FWI of the original data to mitigate the potential misinterpretation of artifacts. Our results demonstrate that it is possible to invert 2D near-surface seismic data even when the data quality is lowered by the presence of strong noise and 3D scattered events caused by objects located in the crossline direction.

Traditional least-squares full-waveform inversion (FWI) suffers from severe local minima problems in case of the presence of strongly dispersive surface waves. Additionally, recorded wavefields are often characterized by amplitude errors due to varying source coupling and incorrect 3D-to-2D geometrical-spreading correction. Thus, least-squares FWI is considered less than suitable for near-surface applications. In this paper, we introduce an amplitude-unbiased coherency measure as a misfit function that can be incorporated into FWI. Such coherency was earlier used in phase-weighted stacking (PWS) to enhance weak but coherent signals. The benefit of this amplitude-unbiased misfit function is that it can extract information uniformly for all seismic signals (surface waves, reflections, and scattered waves). Using the adjoint-state method, we show how to calculate the gradient of this new misfit function. We validate the robustness of the new approach using checkerboard tests and synthetic data contaminated by random noise. We then apply the new FWI approach to a field dataset acquired at an archaeological site located in Ostia, Italy. The goal of this survey was to map the unexcavated archaeological remains with high-resolution. We identify a known tumulus in the FWI results. The instantaneous-phase coherency FWI results also establish that the shallow subsurface under the survey lines is quite heterogeneous. The instantaneous-phase coherency FWI of near-surface data can be a promising tool to image shallow small-scale objects buried under shallow soil covers, as found at archaeological sites.