Internal multiple prediction in the time and offset domains

Andrew Iverson, Kristopher A. H. Innanen, Daniel O. Trad

Varying strategies have been implemented to predict and attenuate internal multiples (Xio et al. 2003). One data driven method to predict internal multiples uses the inverse scattering series (Weglein et al., 1997). The method predicts internal multiples where the only required inputs for the algorithm are the data itself and a search limiting parameter epsilon. It has been displayed that reformulating the calculation domain such that the input and output domains are equivalent allows for the use of a nonstationary epsilon (Innanen, 2015). The version of the algorithm in offset-time is outlined where epsilon can vary in both spatial and temporal dimensions. With the domain of calculation in offset-time this allows for the testing and determining of an epsilon schedule based on the input seismic dataset. For some geologic models in this domain a single epsilon value is insufficient for the prediction of multiples while also minimizing artifacts. Displayed is a nonstationary epsilon implementation that assists in artifact minimization.