Robust estimation of primaries by sparse inversion and Marchenko equation-based workflow for multiple suppression in the case of a shallow water layer and a complex overburden: A 2D case study in the Arabian Gulf

Suppression of surface-related and internal multiples is an outstanding challenge in seismic data processing. The former is particularly difficult in shallow water, whereas the latter is problematic for targets buried under complex, highly scattering overburdens. We have developed a two-step, amplitude- and phase-preserving, inversion-based...

Marchenko without up/down decomposition on the Marmousi model and retrieval of the refracted waves: Are they caused by the Marchenko algorithm?

Marchenko algorithms retrieve the Green’s function for arbitrary subsurface locations, and the retrieved Green’s function includes the primary and multiple reflected waves. The Marchenko algorithms require the estimate of the direct arrivals and the reflected waves; however, most previous Marchenko algorithms also require the up/down components...

A lossless earth Green's function representation between any two subsurface points from surface reflection GPR data

We present a three-dimensional scheme that can be used to compute the electromagnetic impulse response between any two subsurface points from surface reflection data measured at a single surface of a lossless medium. The scheme first computes a virtual vertical radar profile using the Marchenko scheme from which focusing wavefields are...

Marchenko equations for acoustic Green's function retrieval and imaging in dissipative media

We present a scheme for Marchenko imaging in a dissipative heterogeneous medium. The scheme requires measured reflection and transmission data at two sides of the dissipative medium. The effectual medium is the same as the dissipative medium, but with negative dissipation. We show how the measured double-sided data can be combined to obtain the...

Beyond Marchenko: Obtaining virtual receivers and virtual sources in the subsurface

By solving the Marchenko equations, the Green’s function can be retrieved between a virtual receiver in the subsurface to points at the surface (no physical receiver is required at the virtual location). We extend the idea of these equations to retrieve the Green’s function between any two points in the subsurface; i.e, between a virtual source...

Marchenko wavefield redatuming, imaging conditions, and the effect of model errors

Recently, a novel method to redatum the wavefield in the sub-surface from a reflection response measured at the surface has gained interest for imaging primaries in the presence of strong internal multiples. A prerequisite for the algorithm is an accurate and correct estimate of the direct-wave Green's function. However, usually we use an...

An interferometric interpretation of Marchenko redatuming including free-surface multiples

We present an interferometric interpretation of the iterative Marchenko scheme including both free-surface multiples and internal multiples. Cross-correlations are used to illustrate the combination of causal and acausal events that are essential for the process of multiple removal. The first 4 steps in the scheme are discussed in detail, where...

Creating virtual vertical radar profiles from surface reflection ground penetrating radar data

Autofocusing imaging: Imaging with primaries, internal multiples and free-surface multiples

Recent work on autofocusing with the Marchenko equation has shown how the Green's function for a virtual source in the subsurface can be obtained from reflection data. The response to the virtual source is the Green's function from the location of the virtual source to the surface. The Green's function is retrieved using only the reflection...

Internal multiple suppression by adaptive Marchenko redatuming

Recently, a novel iterative scheme was proposed to retrieve Green's functions in an unknown medium from its single-sided reflection response and an estimate of the propagation velocity. In Marchenko imaging, these Green's functions are used for seismic imaging with complete wavefields, including internal multiple reflections. In this way, common...

Data-driven Green's function retrieval and application to imaging with multidimensional deconvolution

An iterative method is presented that allows one to retrieve the Green's function originating from a virtual source located inside a medium using reflection data measured only at the acquisition surface. In addition to the reflection response, an estimate of the travel times corresponding to the direct arrivals is required. However, no detailed...