Internal multiple prediction using inverse-scattering series with sparsity promotion — Part 1: Background, analysis, and synthetic modeling

Geophysics ◽  
2021 ◽  
pp. 1-94
Author(s):  
Ole Edvard Aaker ◽  
Adriana Citlali Ramírez ◽  
Emin Sadikhov
Keyword(s):  
Plane Wave ◽  
Inverse Scattering ◽  
Seismic Data ◽  
Data Completeness ◽  
Inverse Scattering Series ◽  
Recent Developments ◽  
Internal Multiples ◽  
Multiple Prediction ◽  
Synthetic Modeling ◽  
The Way

The presence of internal multiples in seismic data can lead to artefacts in subsurface images ob-tained by conventional migration algorithms. This problem can be ameliorated by removing themultiples prior to migration, if they can be reliably estimated. Recent developments have renewedinterest in the plane wave domain formulations of the inverse scattering series (ISS) internal multipleprediction algorithms. We build on this by considering sparsity promoting plane wave transformsto minimize artefacts and in general improve the prediction output. Furthermore, we argue forthe usage of demigration procedures to enable multidimensional internal multiple prediction withmigrated images, which also facilitate compliance with the strict data completeness requirementsof the ISS algorithm. We believe that a combination of these two techniques, sparsity promotingtransforms and demigration, pave the way for a wider application to new and legacy datasets.

scholarly journals Relating source-receiver interferometry to an inverse-scattering series to derive a new method to estimate internal multiples

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. Q27-Q40 ◽  
Author(s):  
Katrin Löer ◽  
Andrew Curtis ◽  
Giovanni Angelo Meles
Keyword(s):  
Inverse Scattering ◽  
Computational Cost ◽  
Synthetic Data ◽  
Data Set ◽  
Inverse Scattering Series ◽  
Internal Multiples ◽  
Multiple Prediction ◽  
New Formula ◽  
3D Applications ◽  
Insight Into

We have evaluated an explicit relationship between the representations of internal multiples by source-receiver interferometry and an inverse-scattering series. This provides a new insight into the interaction of different terms in each of these internal multiple prediction equations and explains why amplitudes of estimated multiples are typically incorrect. A downside of the existing representations is that their computational cost is extremely high, which can be a precluding factor especially in 3D applications. Using our insight from source-receiver interferometry, we have developed an alternative, computationally more efficient way to predict internal multiples. The new formula is based on crosscorrelation and convolution: two operations that are computationally cheap and routinely used in interferometric methods. We have compared the results of the standard and the alternative formulas qualitatively in terms of the constructed wavefields and quantitatively in terms of the computational cost using examples from a synthetic data set.

A plane-wave formulation and numerical analysis of elastic multicomponent inverse scattering series internal multiple prediction

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. V255-V269 ◽  
Author(s):  
Jian Sun ◽  
Kristopher A. Innanen
Keyword(s):  
Numerical Analysis ◽  
Plane Wave ◽  
Inverse Scattering ◽  
Simulated Data ◽  
Velocity Model ◽  
Prediction Algorithm ◽  
Inverse Scattering Series ◽  
Elastic Data ◽  
Reference Medium ◽  
Multiple Prediction

Internal multiple prediction and removal is a critical component of seismic data processing prior to imaging, inversion, and quantitative interpretation. Inverse scattering series methods predict multiples without identification of generators, and without requiring a velocity model. Land environments present several challenges to the inverse scattering series prediction process. This is particularly true for algorithm versions that explicitly account for elastic conversions and incorporate multicomponent data. The theory for elastic reference medium inverse scattering series internal multiple prediction was introduced several decades ago, but no numerical analysis or practical discussion of how to prepare data for it currently exists. We have focused our efforts on addressing this gap. We extend the theory from 2D to 3D, analyze the properties of the input data required by the existing algorithm, and, motivated by earlier research results, reformulate the algorithm in the plane-wave domain. The success of the prediction process relies on the ordering of events in either pseudodepth or vertical traveltime being the same as the ordering of reflecting interfaces in true depth. In elastic-multicomponent cases, it is difficult to ensure that this holds true because the events to be combined may have undergone multiple conversions as they were created. Several variants of the elastic-multicomponent prediction algorithm are introduced and examined for their tendency to violate ordering requirements (and create artifacts). A plane-wave domain prediction, based on elastic data that have been prepared (1) using variable, “best-fit” velocities as reference velocities, and (2) with an analytically determined vertical traveltime stretching formula, is identified as being optimal in the sense of generating artifact-free predictions with relatively small values of the search parameter [Formula: see text], while remaining fully data driven. These analyses are confirmed with simulated data from a layered model; these are the first numerical examples of elastic-multicomponent inverse scattering series internal multiple prediction.

Internal multiple prediction using inverse scattering series withsparsity promotionPart II: Application strategy and field data examples

Geophysics ◽  
2021 ◽  
pp. 1-52
Author(s):  
Ole Edvard Aaker ◽  
Adriana Citlali Ramírez ◽  
Emin Sadikhov
Keyword(s):  
Inverse Scattering ◽  
Field Data ◽  
Norwegian Sea ◽  
Prediction Quality ◽  
Inverse Scattering Series ◽  
Exploration And Production ◽  
Internal Multiples ◽  
Multiple Prediction ◽  
Seismic Images

Incorrect imaging of internal multiples can lead to substantial imaging artefacts. It is estimatedthat the majority of seismic images available to exploration and production companies have had nodirect attempt at internal multiple removal. In Part I of this article we considered the role of spar-sity promoting transforms for improving practical prediction quality for algorithms derived fromthe inverse scattering series (ISS). Furthermore, we proposed a demigration-migration approach toperform multidimensional internal multiple prediction with migrated data and provided a syntheticproof of concept. In this paper (Part II) we consider application of the demigration-migration approach to field data from the Norwegian Sea, and provide a comparison to a post-stack method (froma previous related work). Beyond application to a wider range of data with the proposed approach,we consider algorithmic and implementational optimizations of the ISS prediction algorithms tofurther improve the applicability of the multidimensional formulations.

Multidimensional inverse-scattering series internal multiple prediction in the coupled plane-wave domain

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. V73-V82 ◽  
Author(s):  
Jian Sun ◽  
Kristopher A. Innanen
Keyword(s):  
Plane Wave ◽  
Inverse Scattering ◽  
Radon Transforms ◽  
Seismic Reflection Data ◽  
Arrival Times ◽  
Computational Expense ◽  
Reflection Data ◽  
Inverse Scattering Series ◽  
Text Prediction ◽  
Multiple Prediction

The inverse-scattering series internal multiple prediction and attenuation algorithm predicts multiples using certain combinations of input seismic reflection data events, which are computed in the wavenumber/pseudodepth or plane-wave/vertical traveltime (i.e., [Formula: see text]) domains. Significant differences can arise in the algorithms’ output and computational expense depending on which domain is used. Many of these are traceable to the response of the algorithm to the users’ choice of the search-limiting parameter [Formula: see text]. The question of which domain is optimal can be addressed with benchmark synthetics. The compactness of the input to the plane-wave domain algorithm leads to the expectation that it will have a reduced computational expense. Also, the lack of increase in the dominant period (i.e., the “width”) of input events as the horizontal slowness increases leads to the expectation that it will respond well to a constant [Formula: see text]. Both of these expectations are borne out with a 1.5D benchmark example. A 2D plane-wave prediction requires the data to be transformed to the [Formula: see text], or coupled plane-wave, domain, involving source- and receiver-side horizontal slownesses. An implementation of this transform leads to the first numerical examples of full 2D inverse series [Formula: see text] prediction. The arrival times, relative amplitudes, and moveout patterns of multiples from dipping horizons are seen in a benchmark synthetic example to be faithfully determined in the plane-wave formulation; waveform mismatches are, however, observed, which are traceable to the numerics of the forward and inverse transforms. High-resolution Radon transforms are a good candidate to improve the match.

Enhancing Internal Multiple Prediction by Using the Inverse Scattering Series

2020 ◽  
Author(s):  
J. Wu ◽  
Z. James Wu ◽  
F. Xavier de Melo ◽  
C. Lapilli ◽  
C. Kostov
Keyword(s):  
Inverse Scattering ◽  
Inverse Scattering Series ◽  
Multiple Prediction

scholarly journals Comparison of monotonicity challenges encountered by the inverse scattering series and the Marchenko demultiple method for elastic waves

Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. Q11-Q26
Author(s):  
C. Reinicke ◽  
M. Dukalski ◽  
K. Wapenaar
Keyword(s):  
Inverse Scattering ◽  
Acoustic Waves ◽  
Prior Information ◽  
Monotonicity Condition ◽  
Scattering Media ◽  
Coupled Modes ◽  
Inverse Scattering Series ◽  
Temporal Event Ordering ◽  
Internal Multiples ◽  
The Time Domain

The reflection response of strongly scattering media often contains complicated interferences between primaries and (internal) multiples, which can lead to imaging artifacts unless handled correctly. Internal multiples can be kinematically predicted, for example by the Jakubowicz method or by the inverse scattering series (ISS), as long as monotonicity, that is, “correct” temporal event ordering, is obeyed. Alternatively, the (conventional) Marchenko method removes all overburden-related wavefield interactions by formulating an inverse problem that can be solved if the Green’s and the so-called focusing functions are separable in the time domain, except for an overlap that must be predicted. For acoustic waves, the assumptions of the aforementioned methods are often satisfied within the recording regimes used for seismic imaging. However, elastic media support wave propagation via coupled modes that travel with distinct velocities. Compared to the acoustic case, not only does the multiple issue become significantly more severe, but also violation of monotonicity becomes much more likely. By quantifying the assumptions of the conventional Marchenko method and the ISS, unexpected similarities as well as differences between the requirements of the two methods come to light. Our analysis demonstrates that the conventional Marchenko method relies on a weaker form of monotonicity. However, this advantage must be compensated by providing more prior information, which in the elastic case is an outstanding challenge. Rewriting, or remixing, the conventional Marchenko scheme removes the need for prior information but leads to a stricter monotonicity condition, which is now almost as strict as for the ISS. Finally, we introduce two strategies on how the remixed Marchenko solutions can be used for imperfect, but achievable, demultiple purposes.

Elimination of land internal multiples based on the inverse scattering series

2011 ◽  
Vol 30 (8) ◽  
pp. 884-889 ◽  
Author(s):  
Yi Luo ◽  
Panos G. Kelamis ◽  
Qiang Fu ◽  
ShouDong Huo ◽  
Ghada Sindi ◽  
...  
Keyword(s):  
Inverse Scattering ◽  
Inverse Scattering Series ◽  
Internal Multiples

Enhancing internal multiple prediction by using the inverse scattering series: methodology and field application

Geophysics ◽  
2021 ◽  
pp. 1-70
Author(s):  
Jing Wu ◽  
Zhiming Wu ◽  
Frederico Xavier de Melo ◽  
Cintia Mariela Lapilli ◽  
Clément Kostov ◽  
...  
Keyword(s):  
Inverse Scattering ◽  
Field Data ◽  
Nearest Neighbor ◽  
Computational Cost ◽  
Field Application ◽  
Nearest Neighbor Search ◽  
Opening Angle ◽  
Model Quality ◽  
Inverse Scattering Series ◽  
Multiple Prediction

We introduce four approaches that dramatically enhance the application of the inverse scattering series method for field data internal multiple prediction. The first approach aims to tackle challenges related to input data conditioning and interpolation. We addressed this through an efficient and fit-for-purpose data regularization strategy, which in this work was a nearest-neighbor search followed by differential moveout to accommodate various acquisition configurations. The second approach addresses cost challenges through applying angle constraints over both the dip angle and opening angle, reducing computational cost without compromising the model’s quality. We also propose an automatic solution for parameterization. The third approach segments the prediction by limiting the range of the multiple’s generator, which can benefit the subsequent adaptive subtraction. The fourth approach works on improving predicted model quality. The strategy includes correctly incorporating the 3D source effect and obliquity factor to enhance the amplitude fidelity of the predicted multiples in terms of frequency spectrum and angle information. We illustrate challenges and report on the improvements in cost, quality or both from the new innovative approaches, using examples from synthetic data and from three field data 2D lines representative of shallow and of deep water environments.

Sign in / Sign up

Export Citation Format

Share Document