3-D Wave simulation in anelastic media using the Kelvin-Voigt constitutive equation

We design a numerical algorithm for wave simulation in anelastic media in the presence of free surface, which can be used to model seismic waves at the Earth's surface and ultrasonic waves in heterogeneous materials. The stress-strain relation is based on the Kelvin-Voigt mechanical model, which has the advantage of not requiring additional field variables. The model requires two anelastic parameters and twice the spatial derivatives of the lossless case. The high-frequency components of the wave field are more attenuated than the low-frequency components, with the attenuation factors being approximately proportional to the square of the frequency. The modeling simulates 3-D waves by using the Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. We stretch the mesh in the vertical direction to increase the minimum grid spacing and reduce the computational cost. Instabilities of the Chebyshev differential operator due to the implementation of the free-surface boundary conditions are solved with a characteristic approach, where the characteristic variables are evaluated at the source central frequency. The results of the modeling are verified by comparisons to the analytical solutions for Lamb's problem and propagation in unbounded homogeneous media. Examples illustrating the propagation of Rayleigh and Love waves are presented.

[1]  Hervé Guillard,et al.  Adaptive spectral methods with application to mixing layer computations , 1992 .

[2]  B. Gurevich,et al.  Seismic attenuation in finely layered porous rocks : Effects of fluid flow and scattering , 1997 .

[3]  J. Carcione.  Constitutive model and wave equations for linear, viscoelastic, anisotropic media , 1995 .

[4]  Ari Ben-Menahem,et al.  Seismic waves and sources , 1981 .

[5]  Harold M. Mooney,et al.  Some numerical solutions for Lamb's problem , 1974, Bulletin of the Seismological Society of America.

[6]  Edward G. S. Paige,et al.  Rayleigh-Wave Theory and Application , 1985 .

[7]  M. Korn,et al.  Incorporation of attenuation into time-domain computations of seismic wave fields , 1987 .

[8]  C. Canuto.  Spectral methods in fluid dynamics , 1991 .

[9]  F. Schubert,et al.  THREE-DIMENSIONAL TIME DOMAIN MODELING OF ULTRASONIC WAVE PROPAGATION IN CONCRETE IN EXPLICIT CONSIDERATION OF AGGREGATES AND POROSITY , 2001 .

[10]  Jianghai Xia,et al.  Estimation of near‐surface shear‐wave velocity by inversion of Rayleigh waves , 1999 .

[11]  Jochen Fröhlich,et al.  A Pseudospectral Chebychev method for the 2D Wave Equation with Domain Stretching and Absorbing Boun , 1995 .

[12]  Moshe Reshef,et al.  Elastic wave calculations by the Fourier method , 1984 .

[13]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[14]  Géza Seriani,et al.  Numerical simulation of interface waves by high‐order spectral modeling techniques , 1992 .

[15]  José M. Carcione,et al.  SOME ASPECTS OF THE PHYSICS AND NUMERICAL MODELING OF BIOT COMPRESSIONAL WAVES , 1995 .

[16]  W. Pilant.  Elastic waves in the earth , 1979 .

[17]  R. Kosloff,et al.  Absorbing boundaries for wave propagation problems , 1986 .

[18]  J. Carcione.  Time-dependent boundary conditions for the 2-D linear anisotropic-viscoelastic wave equation , 1994 .

[19]  Dan Kosloff,et al.  A modified Chebyshev pseudospectral method with an O(N –1 ) time step restriction , 1993 .

[20]  C. Pekeris,et al.  THE SEISMIC SURFACE PULSE. , 1955, Proceedings of the National Academy of Sciences of the United States of America.

[21]  Bengt Fornberg,et al.  The pseudospectral method; accurate representation of interfaces in elastic wave calculations , 1988 .

[22]  D. Kosloff,et al.  Solution of the equations of dynamic elasticity by a Chebychev spectral method , 1990 .

[23]  Alan R. Levander,et al.  Finite-difference forward modeling in seismology , 1989 .

[24]  J. Carcione.  Modeling anelastic singular surface waves in the Earth , 1992 .

[25]  Bernard J. Matkowsky,et al.  Adaptive approximation of solutions to problems with multiple layers by Chebyshev pseudo-spectral methods , 1995 .

[26]  Bengt Fornberg,et al.  A practical guide to pseudospectral methods: Introduction , 1996 .

[27]  Eli Turkel,et al.  A fourth-order accurate finite-difference scheme for the computation of elastic waves , 1986 .

[28]  S. Fagin.  Seismic Modeling of Geologic Structures: Applications to Exploration Problems , 1991 .

[29]  Géza Seriani,et al.  A spectral scheme for wave propagation simulation in 3-D elastic-anisotropic media , 1992 .

[30]  Attenuation properties of viscoelastic material , 1989 .

[31]  H. Levy,et al.  Numerical Solutions of Differential Equations , 1984 .

[32]  Ekkehart Tessmer,et al.  3-D elastic modeling with surface topography by a Chebychev spectral method , 1994 .

[33]  José M. Carcione,et al.  Time-domain Modeling of Constant-Q Seismic Waves Using Fractional Derivatives , 2002 .

[34]  José M. Carcione,et al.  Seismic modeling in viscoelastic media , 1993 .

[35]  J. Augenbaum,et al.  An adaptive pseudospectral method for discontinuous problems , 1989 .

[36]  David E. James,et al.  Encyclopedia of Solid Earth Geophysics , 2020, Encyclopedia of Earth Sciences Series.