Bulletin of the Seismological Society of America
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Bulletin of the Seismological Society of America; February 2009; v. 99; no. 1; p. 302-313; DOI: 10.1785/0120080002
© 2009 Seismological Society of America
This Article
Right arrow Figures Only
Right arrow Full Text
Right arrow Full Text (PDF)
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Pei, Z.-L.
Right arrow Articles by Zhang, L.-X.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation

A Wavelet-Optimized Adaptive Grid Method for Finite-Difference Simulation of Wave Propagation

Zheng-Lin Pei

China National Petroleum Corporation Key Laboratory of Geophysical Exploration, University of Petroleum, Beijing 100083, China

Li-Yun Fu, Geng-Xin Yu, and Lu-Xin Zhang

Key Laboratory of the Earth’s Deep Interior, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China lfu{at}mail.iggcas.ac.cn

To accurately and efficiently model the complex spatiotemporal localization of wave propagation in heterogeneous media, we present a wavelet-optimized gridding scheme for the staggered-grid finite-difference (FD) numerical method. An average-interpolating (AI) wavelet transform with its wavelet coefficients obtained from linear combinations of samples, rather than from integrals, enables the gridding scheme constructed in the physical space to determine where to refine/coarsen the grid adapted to the complex spatial localization of problems and the relatively irregular sampling of seismic signals. The AI wavelet is a perfect transform with its wavelet coefficients decaying rapidly with increasing resolution. The wavelet scale parameter is related directly to the level of meshes. We can reconstruct the solution with good accuracy using a finite number of terms of wavelet expansion under some threshold sets for refinement/coarsening in a low- or high- degree of regularity. The adaptive meshes, in response to the heterogeneities of media and the resolution scales of wave fields, lead not only to a great saving of computation time and memory but also to an enabling method for solving problems with high accuracy in some local complex zones. Numerical examples demonstrate that the proposed wavelet-optimized adaptive FD scheme is efficient and accurate for wave propagation simulation.







HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2009 by the Seismological Society of America.