有限差分的pml1.pdf

CHINESE JOURNAL OF GEOPHYSICS Vol.54, No.3, 2011, pp: 393∼402 A FINITE-DIFFERENCE SCHEME BASED ON PML BOUNDARY CONDITION WITH HIGH POWER GRID STEP VARIATION SUN Lin-Jie, YIN Xing-Yao School of Geosciences, China University of Petroleum, Qingdao 266555, China Abstract Heterogeneous media is a kind of multi-scale media. The scale of the fracture reservoir could reach up to millimeter. Simulating a fractured carbonate model with finite-difference method is calling for the same scale of sample steps. A finite-difference method with high power of grid steps variation is presented here. Spatial and temporal sampling steps can be optimized according to material properties and model geometry. The ratio of grid steps could reach up to hundreds or more. We improve the method of numerical grid division in order to add more flexibility and efficiency to finite-difference method with grid variation. The numerical results support that the method can describe heterogeneous media on millimeter scale accurately with efficiency, and the property made the method suitable for modeling cavity or fractured carbonate reservoirs. Key words Finite-difference, PML, Grid variation, Heterogeneous media 1 INTRODUCTION Seismic wave simulation is to generate seismograms with given subsurface structure and parameters by theoretical calculation[1]. Finite difference method is commonly used in seismic simulation, by using difference scheme to replace the spatial and temporal derivatives in wave equation[2]; the simulation method could achieve satisfactory accuracy and efficiency[3]. Numerical simulation with finite-difference methods usually works with global spatial and temporal grid steps. In situations of complex structures with high veloc