瑞雷波数值模拟中的边界条件及模拟实例分析-中南大学学报.pdf

39 4 () Vol.39 No .4 2008 8 J. Cent. South Univ. (Science and Technology) Aug. 2008 1 2 2 1 (1. 410083 2. 430074) 2 ×12 PML 2 5%2 3 P631 A 1672−7207(2008)04−0824−07 Boundary conditions and case analysis of numerical modeling of Rayleigh wave XIONG Zhang-qiang1 2 2 1 , ZHANG Da-zhou , QIN Zheng , ZHOU Wen-bin (1. School of Info-Physics and Geomatics Engineering, Central South University, Changsha 410083, China; 2. Institute of Geophysics and Geomatics, China University of Geosciences, Wuhan 430074, China) Abstract: A free surface boundary condition for numerical simulation of Rayleigh wave was built by using 2 ×12 order high precision staggered grids finite difference when the source was at the free surface. PML absorbing boundary condition was applied to eliminate boundary reflections. Wave-field modeling of homogeneous isotropic medium and two layer models prove the effectiveness and correctness of using the free surface boundary condition and absorbing boundary condition. The results of model simulating show that relative error between numerical and analytic solutions is within 5%. At the same accurate level, improving difference order has less computation time compared to decreasing grids space. On the base of two-layer model, a seismic record is gotten more practically from modeling the soft interlayer model. The dispersion curve of seismic record has many modes and exhibits the shape