Three-dimensional Elastic Wave Forward Simulation GPU Implementation

Seismic forward simulation technology,as a basic tool for seismology and seismic exploration technology,can determine the relationship between mineral reservoirs and seismic attributes and reveal the seismic response characteristics of different physical properties media.Therefore,in all fields,seismic forward modeling technology has been A wide range of practical applications.In seismic exploration,fully isotropic media and transverse anisotropy media are the most common symmetry models.In this paper,based on the theory of elastic wave fluctuations,forward modeling algorithms for completely isotropic media and transversely anisotropic media are implemented by the finite difference method,and the P-wave component of the forward result of the fully isotropic media is realized.The main research content of this article is as follows:(1)A GPU parallel computing algorithm for numerical simulation of elastic wave equations is implemented,and a general optimization method for GPU parallel algorithms is discussed.For this algorithm,with P2 P technology,the algorithm can be multi-GPU extended,can calculate a larger model and further improve the computational efficiency,and for the repeated calculation of the displacement derivative will encounter the problem of limited access,the introduction of the use of shared memory To improve algorithm performance.(2)The simulation methods for isotropic media and transversely anisotropic media are studied.Through the Lame parameter of the medium,the propagation of the wave field in the isotropic medium can be simulated directly.Through the Thomson parameter representation of the medium,the stiffness matrix coefficients of the medium can be calculated,and the propagation of the seismic waves in the VTI medium and the HTI medium can be positive.Simulation.(3)In the isotropic medium,the rotation and divergence operations of the forward modeling result can separate P-wave and P-wave.However,the calculated results have changed from the original phase and amplitude of the seismic record.The Hilbert transform after the result is multiplied by the normalization factor to restore the original phase and amplitude.