Multigrid deconvolution of seismic data

John Millar, John C. Bancroft

A multigrid method for deconvolution is developed. Tests on synthetic seismic data show the method to be stable, and convergent, and comparisons to a time-domain Gauss- Seidel iterative approach are made to emphasize the long wavelength correction properties of multigrid. The main advantage to the time-domain methods outlined is a straightforward extension to using a nonstationary wavelet. As well, a method for imaging reflectors in-between sample points was found. This may be achieved by interpolating the data to a higher resolution, prior to a multigrid deconvolution.