Seismic reflection imaging, accounting for primary and multiple reflections

An illustration of adaptive Marchenko imaging

In Marchenko imaging, wavefields are retrieved at specified focal points in the subsurface through an iterative scheme derived from the multidimensional Marchenko equation. The method requires seismic-reflection data at the earth’s surface (after free-surface multiple elimination) and an estimate of the direct wavefield from the surface to each...

On Green’s function retrieval by iterative substitution of the coupled Marchenko equations

Iterative substitution of the coupled Marchenko equations is a novel methodology to retrieve the Green's functions from a source or receiver array at an acquisition surface to an arbitrary location in an acoustic medium. The methodology requires as input the single-sided reflection response at the acquisition surface and an initial focusing...

Marchenko imaging

Traditionally, the Marchenko equation forms a basis for 1D inverse scattering problems. A 3D extension of the Marchenko equation enables the retrieval of the Green’s response to a virtual source in the subsurface from reflection measurements at the earth’s surface. This constitutes an important step beyond seismic interferometry. Whereas seismic...

Green's function retrieval from reflection data, in absence of a receiver at the virtual source position

The methodology of Green’s function retrieval by cross-correlation has led to many interesting applications for passive and controlled-source acoustic measurements. In all applications, a virtual source is created at the position of a receiver. Here a method is discussed for Green’s function retrieval from controlled-source reflection data,...

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

Interferometric reservoir monitoring with a single passive source

Changes in the subsurface can be imaged by subtracting seismic reflection data at two different states, one serving as the initial survey or base, and the second as the monitor survey. Conventionally, the reflection data are acquired by placing active seismic sources at the acquisition surface. Alternatively, these data can be acquired from...

Interferometric redatuming of autofocused primaries and internal multiples

Recently, an iterative scheme has been introduced to retrieve the down- and upgoing Green's functions at an arbitrary level ?F inside an acoustic medium as if there were a source at the surface. This scheme requires as input the reflection response acquired at the surface and the direct arrival of the transmission response from the surface to...

Three-dimensional Marchenko equation for Green's function retrieval “beyond seismic interferometry”

In recent work we showed with heuristic arguments that the Green's response to a virtual source in the subsurface can be obtained from reflection data at the surface. This method is called “Green's function retrieval beyond seismic interferometry”, because, unlike in seismic interferometry, no receiver is needed at the position of the virtual...

Deblending by direct inversion

Deblending of simultaneous-source data is usually considered to be an underdetermined inverse problem, which can be solved by an iterative procedure, assuming additional constraints like sparsity and coherency. By exploiting the fact that seismic data are spatially band-limited, deblending of densely sampled sources can be carried out as a...

Controlled-source interferometric redatuming by crosscorrelation and multidimensional deconvolution in elastic media

Various researchers have shown that accurate redatuming of controlled seismic sources to downhole receiver locations can be achieved without requiring a velocity model. By placing receivers in a horizontal or deviated well and turning them into virtual sources, accurate images can be obtained even below a complex near-subsurface. Examples...

Deghosting, demultiple, and deblurring in controlled-source seismic interferometry

With controlled-source seismic interferometry we aim to redatum sources to downhole receiver locations without requiring a velocity model. Interferometry is generally based on a source integral over cross-correlation (CC) pairs of full, perturbed (time-gated), or decomposed wavefields. We provide an overview of ghosts, multiples, and spatial...

Improved surface?wave retrieval from ambient seismic noise by multi?dimensional deconvolution

The methodology of surface?wave retrieval from ambient seismic noise by crosscorrelation relies on the assumption that the noise field is equipartitioned. Deviations from equipartitioning degrade the accuracy of the retrieved surface?wave Green's function. A point?spread function, derived from the same ambient noise field, quantifies the...

A representation for Green’s function retrieval by multidimensional deconvolution

Green’s function retrieval by crosscorrelation may suffer from irregularities in the source distribution, asymmetric illumination, intrinsic losses, etc. Multidimensional deconvolution (MDD) may overcome these limitations. A unified representation for Green’s function retrieval by MDD is proposed. From this representation, it follows that the...

Passive seismic interferometry by multidimensional deconvolution

We introduce seismic interferometry of passive data by multidimensional deconvolution (MDD) as an alternative to the crosscorrelation method. Interferometry by MDD has the potential to correct for the effects of source irregularity, assuming the first arrival can be separated from the full response. MDD applications can range from reservoir...

Deghosting, demultiple, and deblurring in controlled-source seismic interferometry

With controlled-source seismic interferometry we aim to redatum sources to downhole receiver locations without requiring a velocity model. Interferometry is generally based on a source integral over cross-correlation (CC) pairs of full, perturbed (timegated), or decomposed wavefields. We provide an overview of ghosts, multiples, and spatial...