Seismic scalar wave equation modeling by a convolutional differentiator

By means of the inverse Fourier transform, an easy way to design convolutional differentiators is introduced. The tapered differentiators have the same form as a high-order finite difference (FD) operator, and share some of the advantages of both the conventional FD and Fourier methods. The differentiators have been implemented for 2-D zero-offset and common-source seismic modeling. Results using filters with as few as five points are acceptable. Accuracy of the convolutional differentiator approach is appreciably better than that of the fourth-order FD method. The method can readily be extended to elastic wave field modeling. -from Authors