Joost van der Neut

Recently, the focus of reflection seismologists has shifted to applications for which a high-resolution image of the subsurface is required. Least-squares reverse time migration (LSRTM) is a common tool used to compute such images. Even so, its high computational costs have led seismologists to use target-oriented LSRTM for imaging only a small target of interest within a larger subsurface block. Redatuming the data to the upper boundary of the target of interest is one approach to target-oriented LSRTM, but many redatuming methods cannot account for multiple scattering within the overburden. We apply a target-oriented LSRTM algorithm that integrates Marchenko redatuming and double focusing to a field data set. This redatuming method accounts for all orders of multiple scattering in the overburden, thus improving the accuracy of target-oriented LSRTM. Moreover, we determine the effectiveness of a double-focusing algorithm in reducing the data size by decreasing spatial and temporal dimensions of the model and the data. The algorithm’s performance is evaluated using field data acquired in the Norwegian Sea. The numerical results indicate that the target-oriented LSRTM algorithm with Marchenko double-focusing can reduce the internal multiple effects and increase the resolution of the resulting image.

Geophysical monitoring of subsurface reservoirs relies on detecting small changes in the seismic response between a baseline and monitor study. However, internal multiples, related to the over- and underburden, can obstruct the view of the target response, hence complicating the time-lapse analysis. In order to retrieve a response that is free from over- and underburden effects, the data-driven Marchenko method is used. This method effectively isolates the target response, which can then be used to extract more precise time-lapse changes. Additionally, the method also reveals target-related multiples that probe the reservoir more than once, which further define the changes in the reservoir. To verify the effectiveness of the method, a numerical example is constructed. This test shows that when using the isolated target response, the observed time differences resemble the expected time differences in the reservoir. Moreover, the results obtained with target-related multiples also benefit from the Marchenko-based isolation of the reservoir. It is, therefore, concluded that this method has potential to observe dynamic changes in the subsurface with increased accuracy.

Least-Squares Reverse-Time Migration (LSRTM) is a method that seismologists utilize to compute a high-resolution subsurface image. Nevertheless, LSRTM is a computationally demanding problem. One way to reduce the computational costs of the LSRTM is to choose a small region of interest and compute the image of that region. However, finding representations that account for the wavefields entering the target region from the surrounding boundaries is necessary. This paper confines the region of interest between two boundaries above and below this region. The acoustic reciprocity theorem is employed to derive representations for the wavefields at the upper and lower boundaries of the target region. With the help of these representations, a target-enclosed LSRTM algorithm is developed to compute a high-resolution image of the region of interest. Moreover, the possibility of using virtual receivers created by Marchenko redatuming is investigated.

Least-Squares Reverse-Time Migration (LSRTM) is a method that seismologists utilize to compute a high-resolution subsurface image. Nevertheless, LSRTM is a computationally demanding problem. One way to reduce the computational costs of the LSRTM is to choose a small region of interest and compute the image of that region. However, finding representations that account for the wavefields entering the target region from the surrounding boundaries is necessary. This paper confines the region of interest between two boundaries above and below this region. The acoustic reciprocity theorem is employed to derive representations for the wavefields at the upper and lower boundaries of the target region. With the help of these representations, a target-enclosed LSRTM algorithm is developed to compute a high-resolution image of the region of interest. Moreover, the possibility of using virtual receivers created by Marchenko redatuming is investigated.

Geophysical monitoring of subsurface reservoirs relies on detecting small changes in the seismic response between a baseline and monitor study. However, internal multiples, related to the over- and underburden, can obstruct the view of the target response, hence complicating the time-lapse analysis. To retrieve a response that is free from the over- and underburden effects, the data-driven Marchenko method is used. This method effectively isolates the target response, which can then be used to extract more precise time-lapse changes. In addition, the method also reveals target-related multiples that probe the reservoir more than once, which further defines the changes in the reservoir. To verify the effectiveness of the method, a numerical example is constructed. This test finds that, when using the isolated target response, the observed time differences resemble the expected time differences in the reservoir. Moreover, the results obtained with target-related multiples also benefit from the Marchenko-based isolation of the reservoir. It is, therefore, concluded that this method has the potential to observe dynamic changes in the subsurface with increased accuracy.

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.

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 to be collected at the top of the medium. So far , an exact solution could only be obtained if the medium obeys stringent monotonicity conditions and if all forward - scattered (non- converted and converted ) transmissions between the acquisition level and the inner depth level are 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 to be 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 .

SPIE Digital Library Proceedings

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.

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 the artifacts 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 have systematically examined the different approaches to Marchenko redatuming, imaging, and multiple elimination using a common mathematical framework.