Wavefield finite time focusing with reduced spatial exposure

Marchenko Green’s Function Retrieval in Layered Elastic Media from Two-Sided Reflection and Transmission Data

By solving a Marchenko equation, Green’s functions at an arbitrary (inner) depth level inside an unknown elastic layered medium can be retrieved from single-sided reflection data, which are collected at the top of the medium. To date, it has only been possible to obtain an exact solution if the medium obeyed stringent monotonicity conditions and if all forward-scattered (non-converted and converted) transmissions between the acquisition level and the inner depth level were known a priori. We introduce an alternative Marchenko equation by revising the window operators that are applied in its derivation. We also introduce an auxiliary equation for transmission data, which are collected at the bottom of the medium, and a coupled equation, which is based on both reflection and transmission data. We show that the joint system of the Marchenko equation, the auxiliary equation and the coupled equation can be succesfully inverted when broadband reflection and transmission data are available. This results in a novel methodology for elastodynamic Green’s function retrieval from two-sided data. Apart from these data, our approach requires P- and S-wave transmission times between the inner depth level and the top of the medium, as well as two angle-dependent amplitude scaling factors, which can be estimated from the data by enforcing energy conservation.

On the Retrieval of Forward-Scattered Waveforms From Acoustic Reflection and Transmission Data With the Marchenko Equation

A Green's function in an acoustic medium can be retrieved from reflection data by solving a multidimensional Marchenko equation. This procedure requires a priori knowledge of the initial focusing function, which can be interpreted as the inverse of a transmitted wavefield as it would propagate through the medium, excluding (multiply) reflected waveforms. In practice, the initial focusing function is often replaced by a time-reversed direct wave, which is computed with help of a macro velocity model. Green's functions that are retrieved under this (direct-wave) approximation typically lack forward-scattered waveforms and their associated multiple reflections. We examine whether this problem can be mitigated by incorporating transmission data. Based on these transmission data, we derive an auxiliary equation for the forward-scattered components of the initial focusing function. We demonstrate that this equation can be solved in an acoustic medium with mass density contrast and constant propagation velocity. By solving the auxiliary and Marchenko equation successively, we can include forward-scattered waveforms in our Green's function estimates, as we demonstrate with a numerical example.

Data-driven retrieval of primary plane-wave responses

Seismic images provided by reverse time migration can be contaminated by artefacts associated with the migration of multiples. Multiples can corrupt seismic images, producing both false positives, that is by focusing energy at unphysical interfaces, and false negatives, that is by destructively interfering with primaries. Multiple prediction/primary synthesis methods are usually designed to operate on point source gathers and can therefore be computationally demanding when large problems are considered. A computationally attractive scheme that operates on plane-wave datasets is derived by adapting a data-driven point source gathers method, based on convolutions and cross-correlations of the reflection response with itself, to include plane-wave concepts. As a result, the presented algorithm allows fully data-driven synthesis of primary reflections associated with plane-wave source responses. Once primary plane-wave responses are estimated, they are used for multiple-free imaging via plane-wave reverse time migration. Numerical tests of increasing complexity demonstrate the potential of the proposed algorithm to produce multiple-free images from only a small number of plane-wave datasets.