Multiples attenuation is a key process in seismic data processing, and the quality of multiples elimination will affect the final imaging directly. In this paper, we present a 1.5D MATLAB implementation of the inverse scattering series internal multiple prediction algorithm developed by Weglein and collaborators in the 1990s. This method does not require any subsurface information. However, near offset traces information will be needed for applying this method. We discuss the whole prediction operation, and illustrate the procedure with a synthetic example. Effects of various epsilon values chosen will reveal a more efficient method of choosing the epsilon value. Usefulness of our 1.5D internal multiple prediction algorithm in situations where primaries are mixed together with internal multiples, and dipping interface exists are also demonstrated.
View full article as PDF (1.49 Mb)