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

This Article Full Text (PDF) References Citation Map Services Email this article to a friend Similar articles in this journal Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Day, S. M. Search for Related Content GeoRef GeoRef Citation

Efficient simulation of constant Q using coarse-grained memory variables

Department of Geological Sciences San Diego State University, San Diego, California 92182

Abstract

Improvements in computing speed have progressively increased the usable bandwidth of seismic wave-field simulations computed with time-stepped numerical schemes (e.g., finite difference, finite element, pseudospectral). As computational bandwidth increases, anelastic losses become increasingly significant for some important applications such as earthquake ground-motion modeling, whole earth seismogram simulation, and exploration seismic profile modeling, and these losses need to be included in the simulations. As bandwidth increases, however, the memory variables necessary to incorporate realistic anelastic losses account for an increasing proportion of total computational storage requirements, a consequence of the broad relaxation spectrum of typical earth materials. To reduce these storage requirements, we introduce a new method in which the memory variables are coarse grained, that is, redistributed in such a way that only a single relaxation time is represented at each node point (and therefore a single memory variable per stress component is required). Guided by a perturbation analysis, we effect this redistribution in such a way that spatial variability of this single relaxation time simulates the full relaxation spectrum. Such coarse graining reduces memory-variable storage requirements by a factor of 8 for 3D problems or a factor of 4 for 2D problems.

In fourth-order finite-difference computations for the 3D acoustic-wave equation, the method simulates frequency-independent Q within a 3% tolerance over 2 decades in frequency, and it is highly accurate and free of artifacts over the entire usable bandwidth of the underlying finite-difference scheme. These results should also hold for the elastodynamic equations. The method is readily generalized to approximate specific frequency-dependent Q models such as power laws or to further reduce memory requirements. In its present implementation, the main limitation of the method is that it generates artifacts at wavelengths equal to 4 grid cell dimensions and shorter, which may, in some limited circumstances, overlap the usable bandwidth of very high-order finite-difference and/or pseudospectral schemes.

This article has been cited by other articles:

				A. Askan, V. Akcelik, J. Bielak, and O. Ghattas Full Waveform Inversion for Seismic Velocity and Anelastic Losses in Heterogeneous Structures Bulletin of the Seismological Society of America, December 1, 2007; 97(6): 1990 - 2008. [Abstract] [Full Text] [PDF]

				S. Ma and P. Liu Modeling of the Perfectly Matched Layer Absorbing Boundaries and Intrinsic Attenuation in Explicit Finite-Element Methods Bulletin of the Seismological Society of America, October 1, 2006; 96(5): 1779 - 1794. [Abstract] [Full Text] [PDF]

				P. Liu and R. J. Archuleta Efficient Modeling of Q for 3D Numerical Simulation of Wave Propagation Bulletin of the Seismological Society of America, August 1, 2006; 96(4A): 1352 - 1358. [Abstract] [Full Text] [PDF]

				K. B. Olsen, A. Akinci, A. Rovelli, F. Marra, and L. Malagnini 3D Ground-Motion Estimation in Rome, Italy Bulletin of the Seismological Society of America, February 1, 2006; 96(1): 133 - 146. [Abstract] [Full Text] [PDF]

				R. Benites and K. B. Olsen Modeling Strong Ground Motion in the Wellington Metropolitan Area, New Zealand Bulletin of the Seismological Society of America, December 1, 2005; 95(6): 2180 - 2196. [Abstract] [Full Text] [PDF]

				Prediction of Nonlinear Soil Effects Bulletin of the Seismological Society of America, October 1, 2004; 94(5): 1609 - 1629.

				Seismic-Wave Propagation in Viscoelastic Media with Material Discontinuities: A 3D Fourth-Order Staggered-Grid Finite-Difference Modeling Bulletin of the Seismological Society of America, October 1, 2003; 93(5): 2273 - 2280.

				Estimation of Q for Long-Period (>2 sec) Waves in the Los Angeles Basin Bulletin of the Seismological Society of America, April 1, 2003; 93(2): 627 - 638.

				Stability and Accuracy Analysis of Coarse-Grain Viscoelastic Simulations Bulletin of the Seismological Society of America, February 1, 2003; 93(1): 283 - 300.

				Strong Motion in Central Mexico: A Model Based on Data Analysis and Simple Modeling Bulletin of the Seismological Society of America, December 1, 2002; 92(8): 3087 - 3101.

				Simulation of Broadband Ground Motion Including Nonlinear Soil Effects for a Magnitude 6.5 Earthquake on the Seattle Fault, Seattle, Washington Bulletin of the Seismological Society of America, March 1, 2002; 92(2): 831 - 853.

				Viscoelastic Wave Simulation in Basins by a Variable-Grid Finite-Difference Method Bulletin of the Seismological Society of America, December 1, 2001; 91(6): 1741 - 1749.

				Discontinuous-Grid Finite-Difference Seismic Modeling Including Surface Topography Bulletin of the Seismological Society of America, December 1, 2001; 91(6): 1750 - 1764.

				Memory-Efficient Simulation of Anelastic Wave Propagation Bulletin of the Seismological Society of America, June 1, 2001; 91(3): 520 - 531.