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

Inversion of the multidimensional marchenko equation

Focusing functions are defined as wavefields that focus at a specified location in a heterogeneous subsurface. These functions can be directly related to Green's functions and hence they can be used for seismic imaging of complete wavefields, including not only primary reflections but all orders of internal multiples. Recently, it has been shown...

Considerations about the solution space of a VTI marine CSEM Inversion problem using vertical antennas

We exploit the randomness of a genetic inversion algorithm to map the global minimum of the solution space of Controlled-Source Electromagnetic inversion problems. In this study, we focus on the information content that vertical electric or magnetic receivers could add to solve for anisotropic conductivities of the subsurface. By analyzing the...

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

Modeling of GPR data in a stack of VTI-layers with an analytical code

The electromagnetic response in a layered vertical transverse isotropic medium: A new look at an old problem

We determined that the electromagnetic vertical transverse isotropic response in a layered earth can be obtained by solving two equivalent scalar equations, which were for the vertical electric field and for the vertical magnetic field, involving only a scalar global reflection coefficient. Besides the complete derivation of the full...

Probing the solution space of an EM inversion problem with a genetic algorithm

In an inversion for the subsurface conductivity distribution using frequency-domain Controlled-Source Electromagnetic data, various amounts of horizontal components may be included. We investigate which combination of components are best suited to invert for a vertical transverse isotropic (VTI) subsurface. We do this by probing the...

Marchenko imaging below an overburden with random scatterers

Marchenko imaging is a new way to deal with internal multiple scattering in migration. It has been designed for layered media with smooth interfaces. Here we analyze the performance of the Marchenko scheme for a medium with many point scatterers. Although the conditions for Marchenko imaging are violated, we observe from a numerical experiment...

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

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

Data-driven green's function retrieval from reflection data: Theory and example

Recently we introduced a new approach for retrieving the Green's response to a virtual source in the subsurface from reflection data at the surface. Unlike in seismic interferometry, no receiver is needed at the position of the virtual source. Here we present the theory behind this new method. First we introduce the Green's function G and a so...

ESSEMOD - electroseismic and seismoelectric flux-normalized modeling for horizontally layered, radially symmetric configurations

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

Simulating migrated and inverted seismic data by filtering a geologic model

The simulation of migrated and inverted data is hampered by the high computational cost of generating 3D synthetic data, followed by processes of migration and inversion. For example, simulating the migrated seismic signature of subtle stratigraphic traps demands the expensive exercise of 3D forward modeling, followed by 3D migration of the...

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