High-order Finite-difference 3D Seismic Wave Parallel Modeling In The Time-space Domain

Large-scale three-dimensional seismic wave numerical simulation is an important part of imaging and inversion.The current industrial requirements for wave equation numerical simulation technology are,on the one hand,higher simulation accuracy,and on the other hand,faster calculation speed.The thesis is carried out under the above background.The main research content of the thesis focuses on the study of three-dimensional parallel numerical simulation method of high-order accuracy in the space-time domain.On the basis of the traditional spatial high-order and temporal second-order finite difference method,the temporal fourth-order and sixth-order schemes are further considered in this paper.In parallel processing,CPU parallel clusters are used to achieve coarse-grained process-level parallelism between nodes through MPI,and then through Open MP to achieve fine-grained thread-level parallelism within the nodes.In addition,by analyzing the impact of model decomposition direction and decomposition dimension on parallel efficiency,an optimized domain decomposition strategy is obtained,which maximizes the computing potential and improves the computational efficiency without the need for additional hardware investment.Finally,the proposed high-order accuracy optimized parallel scheme is applied to two numerical examples to further verify the feasibility and efficiency.