Enhancing Internal Multiple Prediction by Using the Inverse Scattering Series

Author(s):  
J. Wu ◽  
Z. James Wu ◽  
F. Xavier de Melo ◽  
C. Lapilli ◽  
C. Kostov
Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. Q27-Q40 ◽  
Author(s):  
Katrin Löer ◽  
Andrew Curtis ◽  
Giovanni Angelo Meles

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.


Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. V255-V269 ◽  
Author(s):  
Jian Sun ◽  
Kristopher A. Innanen

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.


Geophysics ◽  
2021 ◽  
pp. 1-94
Author(s):  
Ole Edvard Aaker ◽  
Adriana Citlali Ramírez ◽  
Emin Sadikhov

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.


Geophysics ◽  
2021 ◽  
pp. 1-52
Author(s):  
Ole Edvard Aaker ◽  
Adriana Citlali Ramírez ◽  
Emin Sadikhov

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.


Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. V73-V82 ◽  
Author(s):  
Jian Sun ◽  
Kristopher A. Innanen

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.


Geophysics ◽  
2021 ◽  
pp. 1-70
Author(s):  
Jing Wu ◽  
Zhiming Wu ◽  
Frederico Xavier de Melo ◽  
Cintia Mariela Lapilli ◽  
Clément Kostov ◽  
...  

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.


2020 ◽  
Author(s):  
Jing Wu ◽  
Frederico Xavier de Melo ◽  
Cintia Mariela Lapilli ◽  
Clement Kostov ◽  
Zhiming James Wu

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. S459-S478
Author(s):  
Chao Ma ◽  
Qiang Fu ◽  
Arthur B. Weglein

The industry-standard surface-related multiple elimination (SRME) method provides an approximate predictor of the amplitude and phase of free-surface multiples. This approximate predictor then calls upon an energy-minimization adaptive subtraction step to bridge the difference between the SRME prediction and the actual free-surface multiple. For free-surface multiples that are proximal to other events, the criteria behind energy-minimization adaptive subtraction can be invalid. When applied under these circumstances, a proximal primary can often be damaged. To reduce the dependence on the adaptive process, a more accurate free-surface multiple prediction is required. The inverse scattering series (ISS) free-surface multiple elimination (FSME) method predicts free-surface multiples with accurate time and accurate amplitude of free-surface multiples for a multidimensional earth, directly and without any subsurface information. To quantify these differences, a comparison with analytic data was carried out, confirming that when a free-surface multiple interferes with a primary, applying SRME with adaptive subtraction can and will damage the primary, whereas ISS free-surface elimination will precisely remove the free-surface multiple without damaging the interfering primary. On the other hand, if the free-surface multiple is isolated, then SRME with adaptive subtraction can be a cost-effective toolbox choice. SRME and ISS FSME each have an important and distinct role to play in the seismic toolbox, and each method is the indicated choice under different circumstances.


Sign in / Sign up

Export Citation Format

Share Document