Numerical modeling in seismic survey was performed to simulate wave propagation and to understand the behavior of seismic waves. Compared with seismic survey in field, it has no limits in space and time and can acquire data for various types of model in a short time. To perform numerical modeling, many numerical methods are introduced such as finite difference method(FDM), staggered-grid FDM, spectral element method, and finite element method. Especially, finite difference method is used widely among the various numerical methods. Because it is relatively easy to implement and can accurately simulate the wave propagation for complex model. However,it requires high performance computing resources. To overcome this limitation, more efficient algorithms are needed.
In this study, seismic numerical modeling algorithm was developed. To develop the algorithm, two techniques are applied in time domain. One is logarithmic grid, and the other is expanding domain method. In order to verify the proposed algorithm, several numerical tests are conducted.
Based on numerical test results, reduction of the computation and memory usage is confirmed. Finally, this algorithm can be applied more efficiently to 3D numerical modeling in time domain.