A tour of Marchenko redatuming: Focusing the subsurface wavefield

Marchenko redatuming can retrieve the impulse response to a subsurface virtual source from the single-sided surface reflection data with limited knowledge of the medium. We illustrate the concepts and practical aspects of Marchenko redatuming on a simple 1D acoustic lossless medium in which the coupled Marchenko equations are exact. Defined...

Wavefield focusing using a generalised, potentially asymmetric homogeneous Green's function

Marchenko-type integrals typically relate so-called focusing functions and Green's functions via the reflection response measured on the open surface of a volume of interest. Originating from one dimensional inverse scattering theory, the extension to two and three dimensions set in motion various new developments regarding imaging in complex...

A modified Marchenko method to retrieve the wave field inside layered metamaterial from reflection measurements at the surface

With the Marchenko method, it is possible to retrieve the wave field inside a medium from its reflection response at the surface. To date, this method has predominantly been applied to naturally occurring materials. This study extends the Marchenko method for applications in layered metamaterials with, in the low-frequency limit, effective...

Reciprocity and Representation Theorems for Flux- and Field-Normalised Decomposed Wave Fields

We consider wave propagation problems in which there is a preferred direction of propagation. To account for propagation in preferred directions, the wave equation is decomposed into a set of coupled equations for waves that propagate in opposite directions along the preferred axis. This decomposition is not unique. We discuss flux-normalised...

Velocity-independent Marchenko focusing in time- and depth-imaging domains for media with mild lateral heterogeneity

The Marchenko method retrieves Green's functions between the acquisition surface and any arbitrary point in the medium. The process generally involves solving an inversion starting with an initial focusing function, e.g., a direct-wave Green's function from the desired subsurface position, typically obtained using an approximate velocity...

The Marchenko method for evanescent waves

With the Marchenko method, Green’s functions in the subsurface can be retrieved from seismic reflection data at the surface. State-of-the-art Marchenko methods work well for propagating waves but break down for evanescent waves. This paper discusses a first step towards extending the Marchenko method for evanescent waves and analyses its...

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...

Implementation of the Marchenko multiple elimination algorithm

The Marchenko multiple elimination (MME) and transmission compensation schemes retrieve primary reflections in the two-way traveltime domain without model information or using adaptive subtraction. Both schemes are derived from projected Marchenko equations and are similar to each other, but they use different time-domain truncation operators...

3D Marchenko applications: implementation and examples

We implement the 3D Marchenko equations to retrieve responses to virtual sources inside the subsurface. For this, we require reflection data at the surface of the Earth that contain no free-surface multiples and are densely sampled in space. The required 3D reflection data volume is very large and solving the Marchenko equations requires a...

Marchenko imaging

An interferometric interpretation of Marchenko redatuming

Monitoring changes in the Groningen subsurface using non-physical arrivals in seismic interferometry

Seismic interferometry (SI) refers to the principle of generating seismic responses by crosscorrelating seismic observations at different receiver locations. Theory requires that the boundary sources emit the same energy, have regular spacing and are spaced densely enough. When these assumptions are not met, not only the desired physical...

Monitoring changes inside subsurface layers using non-physical reflections retrieved from seismic interferometry

Seismic interferometry (SI) is a principle for retrieving responses between two receivers using cross-correlation. After the retrieval, one of the receivers acts as a virtual seismic source whose response is retrieved at the second receiver. Correct response retrieval relies on assumptions, among other, of a lossless medium being illuminated...

Application of seismic interferometry by multidimensional deconvolution to earthquakes data recorded in Malargue, Argentina

Marchenko-based target replacement in laterally varying media

Seismic time-lapse studies are generally concerned with variations in a specific target zone, situated inside an otherwise static medium. In seismic monitoring the entire reflection response at the surface needs to be remodeled for every change in the target zone. Ideally, however, only the response of the target zone is remodeled, which is then...

Seismische reflecties