A field data example of Marchenko multiple elimination

Internal multiple reflections have been widely considered as coherent noise in measured seismic data, and many approaches have been developed for their attenuation. The Marchenko multiple elimination (MME) scheme eliminates internal multiple reflections without model information or adaptive subtraction. This schemewas originally derived from coupled Marchenko equations, but it was modified to make it model independent. It filters primary reflections with their two-way traveltimes and physical amplitudes from measured seismic data. The MME scheme is applied to a deepwater field data set from the Norwegian North Sea to evaluate its success in removing internal multiple reflections. The result indicates that most internal multiple reflections are successfully removed and primary reflections masked by overlapping internal multiple reflections are recovered. INTRODUCTION Standard migration schemes map all of the reflections in the data into reflectors in the model domain if they are all primary reflections. These schemes assume that all events in the measured data are reflected only once in the subsurface. Because of this assumption, the measured seismic data should be preprocessed before being migrated by standard migration schemes. Therefore, multiple reflection elimination schemes play a crucial role for standard migration schemes. To date, several approaches have been developed to deal with multiple reflections. Some focus on free-surface-multiple reflections, whereas others focus on internal multiple reflections. Free-surface-multiple reflections can be strong enough to cause artifacts in the image from marine and land data such that much attention has been attracted from industry and academia. Free-surfacemultiple elimination (SRME) (Verschuur et al., 1992) and estimation of primaries by sparse inversion (EPSI) (van Groenestijn and Verschuur, 2009) are the two schemes that have been widely accepted as robust tools for free-surface-multiple attenuation in industry. For SRME, all orders of free-surface-multiple reflections are predicted and a minimum-energy criterion is used to subtract predicted events from the measured data. The EPSI scheme replaces the two-stage processing of SRME, prediction, and adaptive subtraction by an inversion scheme based on the full-waveform inversion approach (van Groenestijn and Verschuur, 2009). Both have achieved success on field data sets. Another strategy is to image the primary and freesurface multiple reflections simultaneously (Brown and Guitton, 2005; Whitmore et al., 2010; Verschuur and Berkhout, 2011; Wang et al., 2014, 2017; Lu et al., 2015), where free-surface multiple reflections give extended illumination of the subsurface. However, crosstalk is present in the resulting image as coherent noise. Less effort has been devoted to deal with internal multiple reflections. As pioneers, Araújo et al. (1994) derive an internal multiple attenuation scheme from the inverse scattering series (ISS). This is the first data-driven scheme that was developed by Weglein et al. (1997) and modified by Ten Kroode (2002) and Löer et al. (2016). Internal multiple elimination (IME) is a layer-related scheme extended from SRME (Berkhout and Verschuur, 1997). The IME scheme downward extrapolates shot records to a virtual surface and attenuates internal multiple reflections related to that surface. Therefore, velocity information is required for its implementation. The ISS and IME schemes have been demonstrated on numerical and field data sets (Matson et al., 1999; Verschuur and Berkhout, 2005; Luo et al., 2011). Adaptive subtraction is needed for both schemes to achieve a multiple-attenuated data set because of the approximate nature of the predicted events. Using internal multiple reflections in imaging is done via full wavefield migration, a dataconsistent closed-loop scheme (Berkhout, 2014). Davydenko and Verschuur (2018) present a field data application. Recently, Marchenko redatuming schemes have been proposed to remove internal multiple reflections and create images free from artifacts (Slob et al., 2014; Wapenaar et al., 2014). Meles et al. (2015) combine convolutional interferometry with the Marchenko scheme to give an internal multiple reflection attenuation scheme. Manuscript received by the Editor 26 May 2019; revised manuscript received 5 September 2019; published ahead of production 30 October 2019; published online 9 January 2020. Delft University of Technology, 2628 CN Delft, The Netherlands. E-mail: l.zhang-1@tudelft.nl (corresponding author); e.c.slob@tudelft.nl. © 2020 Society of Exploration Geophysicists. All rights reserved. S65 GEOPHYSICS, VOL. 85, NO. 2 (MARCH-APRIL 2020); P. S65–S70, 6 FIGS. 10.1190/GEO2019-0327.1 D ow nl oa de d 01 /0 9/ 20 to 1 45 .9 4. 67 .7 1. R ed is tr ib ut io n su bj ec t t o SE G li ce ns e or c op yr ig ht ; s ee T er m s of U se a t h ttp :// lib ra ry .s eg .o rg / Staring et al. (2018) propose to attenuate the first-order internal multiple reflections using an adaptive Marchenko double-focusing method. Model information and adaptive subtraction are required for the implementation of these schemes. Zhang and Staring (2018) modify a Marchenko multiple elimination (MME) scheme (van der Neut and Wapenaar, 2016), which in theory removes all orders of internal multiple reflections without model information or adaptive subtraction. The MME scheme has been extended to also account for transmission loss in primary reflections and free-surface multiple reflections (Zhang and Slob, 2019). Thus, free-surface and internal multiple reflections can be removed and transmission loss in primary reflections can be compensated for in one step without model information or adaptive subtraction. In this paper, the MME scheme is applied to a deepwater field data set from the Norwegian North Sea. It is the first field data example to validate its capabilities for removal of internal multiple reflections without model information or adaptive subtraction. The paper is organized as follows. In the “Theory” section, we give a brief overview of the theory of the MME scheme. The detailed theory can be found in Zhang and Staring (2018). In the “Field example” section, we apply the MME scheme to a field data set for internal multiple reflection elimination. The performance of the MME scheme is analyzed in the “Discussion” section, and we end with our conclusions.

[1]  N. Chemingui,et al.  Separated-wavefield imaging using primary and multiple energy , 2015 .

[2]  Lele Zhang,et al.  Marchenko scheme based internal multiple reflection elimination in acoustic wavefield , 2018, Journal of Applied Geophysics.

[3]  N. D. Whitmore,et al.  Imaging of primaries and multiples using a dual-sensor towed streamer , 2010 .

[4]  Fons ten Kroode,et al.  Prediction of internal multiples , 2002 .

[5]  D. Verschuur,et al.  Including and using internal multiples in closed-loop imaging — Field data examples , 2018, GEOPHYSICS.

[6]  Yi Luo,et al.  Elimination of land internal multiples based on the inverse scattering series , 2011 .

[7]  Helmut Jakubowicz Wave Equation Prediction And Removal of Interbed Multiples , 1998 .

[8]  Evert Slob,et al.  Data-driven internal multiple elimination and its consequences for imaging : A comparison of strategies , 2019 .

[9]  Arthur B. Weglein,et al.  Inverse Scattering Series For Multiple Attenuation: An Example With Surface And Internal Multiples , 1994 .

[10]  A. Berkhout,et al.  Estimation of Multiple Scattering By Iterative Inversion, Part 1: Theoretical Considerations , 1995 .

[11]  Andrew Curtis,et al.  Relating source-receiver interferometry to an inverse-scattering series to derive a new method to estimate internal multiples , 2016 .

[12]  D. J. Verschuur,et al.  Adaptive surface-related multiple elimination , 1992 .

[13]  Evert Slob,et al.  Free-surface and internal multiple elimination in one step without adaptive subtraction , 2019, GEOPHYSICS.

[14]  Evert Slob,et al.  Accounting for free-surface multiples in Marchenko imaging , 2017 .

[15]  A. J. Berkhout,et al.  Review Paper: An outlook on the future of seismic imaging, Part I: forward and reverse modelling , 2014 .

[16]  Matteo Ravasi,et al.  Rayleigh-Marchenko redatuming for target-oriented, true-amplitude imaging , 2017 .

[17]  D. J. Verschuur,et al.  Estimation of multiple scattering by iterative inversion, Part I: Theoretical considerations , 1997 .

[18]  Matteo Ravasi,et al.  Internal multiple prediction and removal using marchenko autofocusing and seismic interferometry , 2015 .

[19]  Joost van der Neut,et al.  Adaptive overburden elimination with the multidimensional Marchenko equation , 2016 .

[20]  Evert Slob,et al.  Transmission compensated primary reflection retrieval in the data domain and consequences for imaging , 2019, GEOPHYSICS.

[21]  Joost van der Neut,et al.  Source-receiver Marchenko redatuming on field data using an adaptive double-focusing method , 2018 .

[22]  Evert Slob,et al.  Seismic reflector imaging using internal multiples with Marchenko-type equations , 2014 .

[23]  Ken H. Matson,et al.  Inverse Scattering Internal Multiple Attenuation: Results From Complex Synthetic And Field Data Examples , 1999 .

[24]  Yibo Wang,et al.  Reverse time migration of multiples: Reducing migration artifacts using the wavefield decomposition imaging condition , 2017 .

[25]  Evert Slob,et al.  Marchenko imagingMarchenko imaging , 2014 .

[26]  Arthur B. Weglein,et al.  An inverse-scattering series method for attenuating multiples in seismic reflection data , 1997 .

[27]  D. J. Verschuur,et al.  Restoration of missing offsets by parabolic Radon transform1 , 1995 .

[28]  Antoine Guitton,et al.  Least-squares joint imaging of multiples and primaries , 2004 .

[29]  D. J. Verschuur,et al.  Seismic migration of blended shot records with surface-related multiple scattering , 2011 .

[30]  D. J. Verschuur,et al.  Removal of internal multiples with the common-focus-point (CFP) approach: Part 2 — Application strategies and data examples , 2005 .

[31]  D. J. Verschuur,et al.  Estimating primaries by sparse inversion and application to near-offset data reconstruction , 2009 .

[32]  A. J. Berkhout,et al.  Review Paper: An outlook on the future of seismic imaging, Part II: Full‐Wavefield Migration , 2014 .