|
|
||||||||
Article |
URS Corporation
566 El Dorado Street
Pasadena, California 91101
robert_graves{at}urscorp.com
(R.W.G.)
Department of Geological Sciences
San Diego State University
San Diego, California 92182
day{at}moho.sdus.edu
(S.M.D.)
Manuscript received 21 March 2002.
We analyze the stability and accuracy of the coarse-grain memory variable technique used for viscoelastic wave-field simulations. Our analysis shows that the general behavior of the coarse-grain system is well described by effective parameters (ME and QE) that are derived from the harmonic average of the moduli over the volume of the coarse-grain cell. The use of these effective parameters proves essential for analyzing the performance and accuracy of the coarse-grain system for Q values less than about 20. A necessary stability condition for the coarse-grain system requires that the weighting coefficients be bounded between 0 and 1. Specifying the weights using the approach of Day and Bradley (2001) satisfies this condition for Q values of about 3 and larger; however, using unconstrained optimization techniques will often produce weights that violate this condition at much higher Q values. We also derive a variation of the original coarse-grain methodology called the element-specific modulus (ESM) formulation in which each element of the coarse-grain cell uses a different unrelaxed modulus. We demonstrate that the accuracy of the coarse-grain system for Q values lower than about 20 is significantly improved with the ESM formulation without increasing the computational cost. Finally, we present a technique for optimizing the accuracy of the coarse-grain system for very low Q based on the use of the effective quality factor (QE). We demonstrate that using conventional optimization techniques that do not employ the effective parameter QE will actually degrade the accuracy of the coarse-grain system.
This article has been cited by other articles:
![]() |
S. R. Ford, D. S. Dreger, K. Mayeda, W. R. Walter, L. Malagnini, and W. S. Phillips Regional Attenuation in Northern California: A Comparison of Five 1D Q Methods Bulletin of the Seismological Society of America, August 1, 2008; 98(4): 2033 - 2046. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. W. Graves The Seismic Response of the San Bernardino Basin Region during the 2001 Big Bear Lake Earthquake Bulletin of the Seismological Society of America, February 1, 2008; 98(1): 241 - 252. [Abstract] [Full Text] [PDF] |
||||
![]() |
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] |
||||
![]() |
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. |
||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |