PML boundary three-dimensional seismic wave propagation simulation method utilizing CUDA

Abstract

The present invention provides a PML (Perfectly Matched Layer) boundary three-dimensional seismic wave propagation simulation method utilizing CUDA (Compute Unified Device Architecture), comprising the steps of: reading in an initial parameter and a Ricker wavelet; decomposing a stress in a first acoustic wave fluctuation equation of a three-dimensional seismic wave into stress components in three directions which are mutually perpendicular, and adding a PML boundary condition to form a second acoustic wave fluctuation equation; using a finite difference method to discretize the second acoustic wave fluctuation equation, so as to form a finite difference acoustic wave fluctuation equation; bordering a first three-dimensional speed model to form a second three-dimensional speed model; calculating attenuation coefficients of edges, angles and surfaces of the second three-dimensional speed model in a CPU (Central Processing Unit) according to the finite difference acoustic wave fluctuation equation; calculating stresses of net points of the second three-dimensional speed model in a GPU (Graphics Processing Unit) according to the attenuation coefficients; and outputting a three-dimensional seismic wave simulation record at the last moment according to the stresses. The PML boundary three-dimensional seismic wave propagation simulation method of the present invention can achieve high speed-up ratio and shorten simulation time.

Description

technical field [0001] The present invention relates to seismic wave simulation technology, in particular to a three-dimensional GPU (Graphics Processing Unit) seismic wave propagation simulation technology with PML (Perfect Matching Layer) boundaries, specifically a CUDA (Compute Unified Device Architecture, computing unified device architecture) The PML Boundary 3D Seismic Wave Propagation Simulation Method. Background technique [0002] Seismic exploration is the upstream business of the oil industry, and has made significant contributions to the discovery of Shengli Oilfield, Dagang Oilfield and Huabei Oilfield. Nowadays, it is very necessary to develop a high-precision seismic exploration method under the situation that it is increasingly difficult to find oil reserves in Chinese territory of. [0003] Seismic exploration includes seismic acquisition and data processing. How to extract effective information from the collected data and then obtain the distribution of un...

AI Technical Summary

Problems solved by technology

[0005] The conventional implementation of seismic wave simulation is to use the central processing unit (CPU) to calculate, but even with the help of parallel technologies such as OpenMP and MPI (multi-point interface), the calculation is still very time-consuming and cannot meet the requirements of use, and the reverse time migration technology It also needs to use the seismic wave simulation process many times. Therefore, the time-consuming seismic wave simulation makes the reverse time migration technology with a huge amount of calculation impossible to be practically applied.

Images