We have developed an adaptive double-focusing method that is specifically designed for the field-data application of source-receiver Marchenko redatuming. Typically, the single-focusing Marchenko method is combined with a multidimensional deconvolution (MDD) to achieve redatuming. Our method replaces the MDD step by a second focusing step that naturally complements the single-focusing Marchenko method. Instead of performing the MDD method with the directionally decomposed Green’s functions that result from single-focusing, we now use the retrieved upgoing Green’s function and the retrieved downgoing focusing function to obtain a redatumed reflection response in the physical medium. Consequently, we only remove the strongest overburden effects instead of removing all of the overburden effects. However, the gain is a robust method that is less sensitive to imperfections in the data and a sparse acquisition geometry than the MDD method. In addition, it is computationally much cheaper, more straightforward to implement, and it can be parallelized over pairs of focal points, which makes it suitable for application to large data volumes. We evaluate the successful application of our method to 2D field data of the Santos Basin.

No abstract

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 propose a two-step, amplitude- and phase-preserving, inversion-based workflow, which addresses these problems. We apply Robust Estimation of Primaries by Sparse Inversion (R-EPSI) to suppress the surface-related multiples and solve for 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 de-multiple methods such as SRME. The resulting Green's function is used as input to a Marchenko equation-based approach to predict the complex interference pattern of all overburden-generated internal multiples at once, without a priori subsurface information. 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 pre-processing 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 that need to be corrected for. Hence, we decided that the use of conservative adaptive filters is necessary to obtain the best results after interbed multiple removal. The obtained results show promising suppression of both surface-related and interbed multiples.

In recent years, a variety of Marchenko methods for the attenuation of internal multiples has been developed. These methods have been extensively tested on two‐dimensional synthetic data and applied to two‐dimensional field data, but only little is known about their behaviour on three‐dimensional synthetic data and three‐dimensional field data. Particularly, it is not known whether Marchenko methods are sufficiently robust for sparse acquisition geometries that are found in practice. Therefore, we start by performing a series of synthetic tests to identify the key acquisition parameters and limitations that affect the result of three‐dimensional Marchenko internal multiple prediction and subtraction using an adaptive double‐focusing method. Based on these tests, we define an interpolation strategy and use it for the field data application. Starting from a wide azimuth dense grid of sources and receivers, a series of decimation tests are performed until a narrow azimuth streamer geometry remains. We evaluate the effect of the removal of sail lines, near offsets, far offsets and outer cables on the result of the adaptive double‐focusing method. These tests show that our method is most sensitive to the limited aperture in the crossline direction and the sail line spacing when applying it to synthetic narrow azimuth streamer data. The sail line spacing can be interpolated, but the aperture in the crossline direction is a limitation of the acquisition. Next, we apply the adaptive Marchenko double‐focusing method to the narrow azimuth streamer field data from the Santos Basin, Brazil. Internal multiples are predicted and adaptively subtracted, thereby improving the geological interpretation of the target area. These results imply that our adaptive double‐focusing method is sufficiently robust for the application to three‐dimensional field data, although the key acquisition parameters and limitations will naturally differ in other geological settings and for other types of acquisition.