Many seismic imaging methods use wavefield extrapolation operators to redatum sources and receivers from the surface into the subsurface. We discuss wavefield extrapolation operators that account for internal multiple reflections, in particular propagator matrices, transfer matrices and Marchenko focusing functions. A propagator matrix is a square matrix that 'propagates' a wavefield vector from one depth level to another. It accounts for primaries and multiples and holds for propagating and evanescent waves. A Marchenko focusing function is a wavefield that focuses at a designated point in space at zero time. Marchenko focusing functions are useful for retrieving the wavefield inside a heterogeneous medium from the reflection response at its surface. By expressing these focusing functions in terms of the propagator matrix, the usual approximations (such as ignoring evanescent waves) are avoided. While a propagator matrix acts on the full wavefield vector, a transfer matrix (according to the definition used in this paper) 'transfers' a decomposed wavefield vector (containing downgoing and upgoing waves) from one depth level to another. It can be expressed in terms of decomposed Marchenko focusing functions. We present propagator matrices, transfer matrices and Marchenko focusing functions in a consistent way and discuss their mutual relations. In the main text we consider the acoustic situation and in the appendices we discuss other wave phenomena. Understanding these mutual connections may lead to new developments of Marchenko theory and its applications in wavefield focusing, Green's function retrieval and imaging.

@en

Much of recent progress in geophysics can be attributed to the adaptation of heterogeneous high-performance computing architectures. It is projected that the next major leap in many areas of science, and hence hopefully in geophysics too, will be due to the emergence of quantum computers. Finding a right combination of hardware, algorithms, and a use case, however, proves to be a very challenging task - especially when looking for a relevant application that scales efficiently on a quantum computer and is difficult to solve using classical means. We find that maximizing stack power for residual statics correction, an NP-hard combinatorial optimization problem, appears to naturally fit a particular type of quantum computing known as quantum annealing. We express the underlying objective function as a quadratic unconstrained binary optimization, which is a quantum-native formulation of the problem. We choose some solution space and define a proper encoding to translate the problem variables into qubit states. We find that these choices can have a significant impact on the maximum problem size that can fit on the quantum annealer and on the fidelity of the final result. To improve the latter, we embed the quantum optimization step in a hybrid classical-quantum workflow, which aims to increase the frequency of finding the global, rather than some local, optimum of the objective function. Finally, we find that a generic, black-box, hybrid classical-quantum solver also could be used to solve stack-power maximization problems proximal to industrial relevance and capable of surpassing deterministic solvers prone to cycle skipping. A custom-built workflow capable of solving larger problems with an even higher robustness and greater control of the user appears to be within reach in the very near future.

@en

Marchenko equation-based methods promise data-driven, true-amplitude internal multiple elimination. The method is exact in 1-D acoustic media, however it needs to be expanded to account for the presence of 2- and 3-D elastodynamic wave-field phenomena, such as compressional (P) to shear (S) mode conversions, total reflections or evanescent waves. Mastering high waveform-fidelity methods such as this, could further advance amplitude vs offset analysis and lead to improved reservoir characterization. This method-expansion may comprise of re-evaluating the underlying assumptions and/or appending the scheme with additional constraints (e.g. minimum phase). To do that, one may need to better understand the construction of the Marchenko equation solutions, the so-called focusing functions, in a mathematically simple and numerically stable fashion. The latter could be a challenge at large angles of incidence where the elastodynamic effects and evanescent waves start playing a dominant role. We demonstrate that the elastodynamic focusing functions are the bridge between the Marchenko equation theory and the transfer matrix formalism. Using the latter, we show how we can try to gain further insights into how time-reversal (correlations) behaves when either of the elastic modes becomes evanescent. We also show how this construction allows us to shed light on into the mathematical properties of elastodynamic inverse transmissions, which takes us a step closer towards understanding the elastodynamic minimum phase reconstruction.@en

Standard Marchenko redatuming and imaging schemes neglect evanescent waves and are based on the assumption that decomposition into downgoing and upgoing waves is possible in the subsurface. Recently we have shown that propagator matrices, which circumvent these assumptions, can be expressed in terms of Marchenko focusing functions. In this paper we generalize the relation between the propagator matrix and the Marchenko focusing functions for a 3D inhomogeneous dissipative medium. Moreover, for the same type of medium we discuss a relation between the transfer matrix and the Marchenko focusing functions.@en

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 workflow that addresses these problems. We apply robust estimation of primaries by sparse inversion (R-EPSI) to solve simultaneously for the surface-related primaries Green’s function and the source wavelet. A significant advantage of the inversion approach of the R-EPSI method is that it does not rely on an adaptive subtraction step that typically limits other demultiple methods such as surface-related multiple elimination. The resulting Green’s function is used as the input to a Marchenko equation-based approach to predict the complex interference pattern of all overburden-generated internal multiples at once. In this approach, no a priori information about the subsurface is needed. In theory, the interbed multiples can be predicted with correct amplitude and phase and, again, no adaptive filters are required. We illustrate this workflow by applying it on an Arabian Gulf field data example. It is crucial that all preprocessing steps are performed in an amplitude-preserving way to restrict any impact on the accuracy of the multiple prediction. In practice, some minor inaccuracies in the processing flow may end up as prediction errors for which corrections will be needed. Hence, we conclude that the use of conservative adaptive filters were necessary to obtain the best results after interbed multiple removal. The obtained results indicate promising suppression of surface-related and interbed multiples.@en

With the Marchenko method it is possible to retrieve Green's functions between virtual sources in the subsurface and receivers at the surface from reflection data at the surface and focusing functions. A macro model of the subsurface is needed to estimate the first arrival; the internal multiples are retrieved entirely from the reflection data. The retrieved Green's functions form the input for redatuming by multidimensional deconvolution (MDD). The redatumed reflection response is free of internal multiples related to the overburden. Alternatively, the redatumed response can be obtained by applying a second focusing function to the retrieved Green's functions. This process is called Marchenko redatuming by double focusing. It is more stable and better suited for an adaptive implementation than Marchenko redatuming by MDD, but it does not eliminate the multiples between the target and the overburden. An attractive efficient alternative is plane-wave Marchenko redatuming, which retrieves the responses to a limited number of plane-wave sources at the redatuming level. In all cases, an image of the subsurface can be obtained from the redatumed data, free of artefacts caused by internal multiples. Another class of Marchenko methods aims at eliminating the internal multiples from the reflection data, while keeping the sources and receivers at the surface. A specific characteristic of this form of multiple elimination is that it predicts and subtracts all orders of internal multiples with the correct amplitude, without needing a macro subsurface model. Like Marchenko redatuming, Marchenko multiple elimination can be implemented as an MDD process, a double dereverberation process, or an efficient plane-wave oriented process. We systematically discuss the different approaches to Marchenko redatuming, imaging and multiple elimination, using a common mathematical framework.@en

Correct handling of strong elastic, internal, multiples remains a challenge for seismic imaging. Methods aimed at eliminating them are currently limited by monotonicity violations, a lack of a-priori knowledge about mode conversions, or unavailability of multi-component sources and receivers for not only particle velocities but also the traction vector. Most of these challenges vanish in acoustic media such that Marchenko-equation-based methods are able in theory to remove multiples exactly (within a certain wavenumber-frequency band). In practice, however, when applied to (elastic) field data, mode conversions are unaccounted for. Aiming to support a recently published marine field data study, we build a representative synthetic model. For this setting, we demonstrate that mode conversions can have a substantial impact on the recovered multiple-free reflection response. Nevertheless, the images are significantly improved by acoustic multiple elimination. Moreover, after migration the imprint of elastic effects is considerably weaker and unlikely to alter the seismic interpretation.@en