HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Figures Only Full Text Full Text (PDF) Services Email this article to a friend Similar articles in this journal Similar articles in Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Web of Science (2) Citing Articles via Google Scholar Google Scholar Articles by Yang, D. Articles by Lu, M. Search for Related Content GeoRef GeoRef Citation

Optimally Accurate Nearly Analytic Discrete Scheme for Wave-Field Simulation in 3D Anisotropic Media

¹ Department of Mathematical Sciences, Tsinghua University
Beijing 100084, People’s Republic of China
dhyang{at}math.tsinghua.edu.cn
 (D.Y., G.S., M.L.)
² College of Sciences, Southwest Petroleum University
Chengdu 610500, People’s Republic of China
 (G.S.)

We present a new numerical method for elastic wave modeling in 3D isotropic and anisotropic media, which is called the 3D optimal nearly analytic discrete method (ONADM) in this article. This work is an extension of the 2D ONADM (Yang et al., 2004) that models acoustic and elastic waves propagating in 2D media. The formulation is derived by using a multivariable truncated Taylor series expansion and high-order interpolation approximations. Our 3D ONADM enables wave propagation to be simulated in three dimensions through isotropic and anisotropic models. Promising numerical results show that the error of the ONADM for the 3D case is less than those of the conventional finite-difference (FD) method and the so-called Lax–Wendroff correction (LWC) schemes, measured quantitatively by the root- mean-square deviation from analytical solution. The seismic wave fields in the 3D isotropic and anisotropic media are simulated and compared with those obtained by using the fourth-order LWC method and exact solutions for the acoustic wave case. We show that, compared with the conventional FD method and the LWC schemes, the ONADM for the 3D case can reduce significantly the storage space and computational costs. Numerical experiments illustrate that the ONADM provides a useful tool for the 3D large-scale isotropic and anisotropic problems and it can suppress effectively numerical dispersions caused by discretizing the 3D wave equations when too-coarse grids are used, which is the same as the 2D ONADM. Numerical modeling also implies that simultaneously using both the wave displacement and its gradients to approximate the high-order derivatives is important for both decreasing the numerical dispersion caused by the discretization of wave equations and compensating the important wave-field information included in wave-displacement gradients.