Propagator and transfer matrices, Marchenko focusing functions and their mutual relations

Abstract

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.