Research On High Performance Computing Method For Numerical Simulation Of Seismic Response

Seismic wave field numerical simulation is the propagation process of seismic wave through the underground medium,and the numerical simulation of seismic wave field plays a fundamental role in seismic acquisition,data processing and seismic interpretation.Nowadays one-way wave method is mainly used in seismic wave field numerical simulation,which simulates a reflected wave through the upward wave of seismic wave in upward propagation process and a diffracted wave as a result of wave impedance mutation point,and it is simple and efficient.One-way wave is widely used in the construction of its operator.Split-Step Fourier Migration is not able to adapt to the dramatic change of lateral velocity,while the Phase shift plus interpolation can nearly realize the speed in the change of the vertical direction and horizontal direction,but it needs a number of reference velocity in each layer.The corresponding numerical simulation of seismic wave field can simulate the seismic model of complex media and obtain accurate wave field information.However,to the abnormal geological seismic wave field numerical simulation,the existing method is to simulate in the way of large grid.It will give rise to small cracks and its development belt and cavity edge information cannot be reflected,which fails to meet the need of seismic interpretation.In view of the above questions,this article analyzes the selection of continuation operator,the processing of geologic abnormal body,increasing the speed of seismic wave field numerical simulation and so on.The specific research works are as follows:1.The main difficulty of the numerical simulation of seismic wavefield is in the constitution of tits operator,which adapts to the rapid and lateral changes in velocity.This article has been studied the advantages and disadvantages of the Split-Step Fourier operator and phase-shift interpolation continuation operator.And then constructing continuation operators that adapt to dramatic changes in lateral velocity.2.In the forward modeling,the general method is to simulate in the way of large grid.regard to a geologic body of the small scale(such as fault fracture zone,fracture development zone),the large grid irregular shape(solution pores)and the change of lateral velocity.combined with the idea of the grid subdivision of the geological model of thought,this paper comes up with an seismic wave field continuation operator based on adaptive grid subdivision that can be adapted to the fast change of lateral velocity.3.With regard to a mass calculation after the small grid,this paper presents different and parallel designs responding to seismic wave field numerical simulation.The poststack seismic wavefield numerical simulation is simpler than the prestack,and it is parallelized by three modes based on MPI,OpenMP,MPI+OpenMP.Prestack seismic wavefield numerical simulation uses a parallel model based on MPI,MPI+OpenMP and MPI+OpenMP+CUDA.4.Based on the VS2012 development tool,a high-performance computing system based on MPI+OpenMP seismic response numerical simulation was implemented.For the differences in pre-stack and post-stack seismic responses,different window settings were made and different modes were selected.Finally,a high-performance computing system based on numerical simulation of seismic response is fully implemented.