I present a finite-difference method for modeling P-SV wave propagation in heterogeneous media This is an extension of the method I previously proposed for modeling SH-wave propagation by using velocity and stress in a discrete grid The two components of the velocity cannot be defined at the same node for a complete staggered grid: the stability condition and the P-wave phase velocity dispersion curve do not depend on the Poisson's ratio, while the S-wave phase velocity dispersion curve behavior is rather insensitive to the Poisson's ratio Therefore, the same code used for elastic media can be used for liquid media, where S-wave velocity goes to zero, and no special treatment is needed for a liquid-solid interface Typical physical phenomena arising with P-SV modeling, such as surface waves, are in agreement with analytical results The weathered-layer and corner-edge models show in seismograms the same converted phases obtained by previous authors This method gives stable results for step discontinuities, as shown for a liquid layer above an elastic half-space The head wave preserves the correct amplitude Finally, the corner-edge model illustrates a more complex geometry for the liquid-solid interface As the Poisson's ratio v increases from 025 to 05, the shear converted phases are removed from seismograms and from the time section of the wave field

/pdf/p-sv-wave-propagation-in-heterogeneous-media-velocity-stress-1srlln7ms7.pdf

P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method

A large and growing number of original papers on both the experimental and the theoretical aspects of stress wave propagation is appearing in the scientific literature, and two international conferences solely concerned with the subject have been held during the last five years. The purpose of this paper is to review recent experimental and theoretical advances in the propagation of deformation waves of arbitrary shape through elastic and anelastic solids, and also to attempt to outline the problems on which present efforts are being directed and to predict probable lines of future development.

Stress Waves in Solids

A gross Earth datum is a single measurable number describing some property of the whole Earth, such as mass, moment of interia, or the frequency of oscillation of some identified elastic-gravitational normal mode. We show how to determine whether a given finite set of gross Earth data can be used to specify an Earth structure uniquely except for fine-scale detail; and how to determine the shortest length scale which the given data can resolve at any particular depth. We apply the general theory to the linear problem of finding the depth-variation of a frequency-independent local Q from the observed quality factors Q of a finite number of normal modes. We also apply the theory to the non-linear problem of finding density vs depth from the total mass, moment, and normal-mode frequencies, in case the compressional and shear velocities are known.

/pdf/the-resolving-power-of-gross-earth-data-31h1wjivgj.pdf

The Resolving Power of Gross Earth Data

Publisher Summary This chapter focuses on elastic wave propagation in stratified media. The development of the theory of elastic wave propagation in stratified media has been strongly influenced by the problems of seismic wave propagation and the nature of the seismograms recorded from earthquakes. For purely analytic developments of elastic wave propagation, the level of manageable algebraic complexity is reached in a model with one or two uniform layers overlying a uniform half space. This chapter shows how the excitation of elastic waves, within a horizontally stratified structure, can be conveniently developed in terms of reflection and transmission matrices. This procedure has allowed the construction of the full response of the medium or approximations with desired properties so that theoretical seismograms may be calculated for realistic distributions of elastic parameters. Although this development has been for isotropic media, nearly all the results apply directly to the case of full anisotropy if 3 × 3 reflection and transmission matrices allowing coupling between all wave types are employed. This development of the wavefield for both source and receiver within the stratification may be used for other classes of wave propagation.

Seismic Wave Propagation in Stratified Media

Research into ultrasonic NDE techniques for the inspection of multilayered structures relies strongly on the use of modeling tools which calculate dispersion curves and reflection and transmission spectra. These predictions are essential to enable the best inspection strategies to be identified and their sensitivities to be evaluated. General purpose multilayer modeling tools may be developed from a number of matrix formulations which have evolved in the latter half of this century and there is now a formidable number of publications on the subject. This paper presents a review of the main developments of the matrix techniques, and their use in response and modal models, with emphasis on ultrasonics applications. >

W. M. Ewing

Papers

Elastic Waves in Layered Media