A pseudospectral‐finite difference hybrid approach for large‐scale seismic modeling and RTM on parallel computers

We present numerical comparisons between the finite difference method and the Fourier pseudospectral method on the 3D SEG/EAGE salt model to investigate the potentials of the two methods for large scale seismic modeling and reverse-time migration problems on distributed-memory parallel systems. We compare the accuracy of the seismic modeling results by the two methods. We discuss the parallel implementations and present the runtime and speedup rate data on a Linux cluster. Based on these comparison results, we propose a pseudospectral-finite difference hybrid method for largescale seismic modeling and RTM problems on massively distributed-memory parallel systems.