Generalized internal multiple imaging

Internal multiples deteriorate the image when the imaging procedure assumes only single scattering, especially if the velocity model does not have sharp contrasts to reproduce such scattering in the Green's function through forward modeling. If properly imaged, internal multiples (internally scattered energy) can enhance the seismic image. Conventionally, to image internal multiples, accurate, sharp contrasts in the velocity model are required to construct a Green's function with all the scattered energy. As an alternative, we have developed a generalized internal multiple imaging procedure that images any order internal scattering using the background Green's function (from the surface to each image point), constructed from a smooth velocity model, usually used for conventional imaging. For the first-order internal multiples, the approach consisted of three steps, in which we first back propagated the recorded surface seismic data using the background Green's function, then crosscorrelated the back-propagated data with the recorded data, and finally crosscorrelated the result with the original background Green's function. This procedure images the contribution of the recorded first-order internal multiples, and it is almost free of the single-scattering recorded energy. The cost includes one additional crosscorrelation over the conventional single-scattering imaging application. We generalized this method to image internal multiples of any order separately. The resulting images can be added to the conventional single-scattering image, obtained, e.g., from Kirchhoff or reverse-time migration, to enhance the image. Application to synthetic data with reflectors illuminated by multiple scattering (double scattering) demonstrated the effectiveness of the approach.