Parallel large-scale variable grid 3D seismic wave modeling
Cai Zhicheng1,2, Gu Hanming1,2, Cheng Jingwang3, Liu Chuncheng4, Liu Zhibin4, Liu Shaoyong1,2
1. Institute of Geophysics and Geomatics, China University of Geosciences (Wuhan), Wuhan, Hubei 430074, China;
2. Hubei Key Laboratory of Subsurface Multi-scale Imaging, Wuhan, Hubei 430074, China;
3. College of Geophysics and Oil Resource, Yangzte University, Wuhan, Hubei 430100, China;
4. CNOOC Research Institute, Beijing 100027, China
Abstract:3D wave equation numerical modeling is one of the key approaches in the analysis of complex reservoir wavefield characteristics and seismic response recognition model building. This numerical modeling however,needs huge consuming time and large memory. We use MPI parallel to deal with the problems. First the whole calculation task is divided into a few parts and distributed to a number of CPUs to accomplish the modeling. Then the variable grid algorithm is applied in wave equation,and smaller grids are used in geologically complex regions in order to effectively and accurately simulate wave propagation. Therefore,our special MPI parallel strategy is developed,and it is more appropriate to the large scale variable grid 3D seismic modeling algorithm. The new strategy reduces computation consuming time and memory size,and further improves the efficiency of the normal large MPI seismic wave modeling algorithm. Model tests prove that the proposed strategy of parallel large-scale variable grid 3D seismic wave modeling algorithm is more efficient and accurate.
Dong Liangguo,Ma Zaitian,Cao Jingzhong et al. A staggered-grid high-order difference method of one-order elastic wave equation. Chinese Journal of Geophysics,2000,43(3):411-419.
Pei Zhenglin. Three-dimensional numerical simulation of elastic wave propagation in 3-D anisotropic media with staggered-grid high-order difference method. Journal of University of Petroleum(Natural Science Edition),2004,28(5):23-29.
Pei Zhenglin. Numerical simulation of elastic wave equation in 3-D isotropic media with staggered-grid high-order difference method. GPP,2005,44(4):308-315.
[6]
Moczo P. Finite-difference technique for SH-waves in 2-D media using irregular grids-application to the.response problem. Geophysical Journal International,1989,99(2):321-329.
[7]
Jastram C,Behle A. Acoustic modelling on a grid of vertically varying spacing. Geophysical Prospecting,1992,40(2):157-169.
[8]
Jastram C,Tessmer E.Elastic modelling on a grid with vertically varying spacing. Geophysical Prospecting,1994,42(4):357-370.
[9]
Kang T S,Baag C E. Finite-difference.simulation combining discontinuous grids with locally variable time steps. Bulletin of the Seismological Society of America,2004,94(1):207-219.
[10]
Kang T S,Baag C E.An efficient finite-difference method for simulating 3D.response of locali-zed basin structures. Bulletin of the SeismologicalSociety of America,2004,94(5):1690-1705.
Sun Chengyu,Li Shengjun,Ni Changkuan et al. Wave equation numerical modeling by finite difference method with varying grid spacing. GPP,2008,47(2):123-127.
Sun Chengyu,Ding Yucai. Accuracy analysis of wave equation finite difference with dual-variable grid algorithm. OGP,2012,47(4):545-551.
[13]
Hayashi K,Burns D R,Toks z M N. Discontinuous-grid finite-difference.modeling including surface topography. Bulletin of the Seismological Society of America,2001,91(6):1750-1764.
Li Zhenchun,Li Qingyang,Huang Jianping et al. A stable and high-precision dual-variable grid forward modeling and reverse time migration method. GPP,2014,53(2):127-136.
Micikevicius P.3D finite difference computation on GPUs using CUDA. Proceedings of 2nd Workshop on General Purpose Processing on Graphics Processing Units,Association for Computing Machinery,2009,79-84.
He Bingshou,Zhang Huixing,Han Linghe. Forward modeling of elastic wave equation with coarse-grained parallel algorithm. Progress in Geophysics,2010,25(2):650-656.