Propagation of sound pulses in a homogeneous layer under an inhomogeneous half space

1968 ◽  
Vol 68 (3) ◽  
pp. 131-148
Author(s):  
Mani Prasad Sharma
1996 ◽  
Vol 86 (3) ◽  
pp. 750-760
Author(s):  
F. Abramovici ◽  
L. H. T. Le ◽  
E. R. Kanasewich

Abstract This article presents some numerical experiments in using a computer program for calculating the displacements due to a P source in a vertically inhomogeneous structure, based on the Fourier-Bessel representation. The structure may contain homogeneous, inhomogeneous, elastic, or viscoelastic layers. The source may act in any type of sublayer or in the half-space. Synthetic results for the simple case of a homogeneous layer overlaying a homogeneous half-space compare favorably with computations based on the Cagniard method. Numerical seismograms for an elastic layer having velocities and density varying linearly with depth were computed by integrating numerically the governing differential systems and compared with results based on the Haskell model of splitting the linear layer in homogeneous sublayers. Even an adaptive process with a variable step size based on the Haskell model has a poorer performance on the accuracy-cpu time scale than numerical integration.


1982 ◽  
Vol 72 (5) ◽  
pp. 1483-1498
Author(s):  
F. Abramovici ◽  
E. R. Kanasewich ◽  
P. G. Kelamis

abstract The displacement components for a horizontal stress discontinuity along a buried finite fault in an elastic homogeneous layer on top of an elastic half-space are given analytically in terms of generalized rays. For a particular case of a concentrated horizontal force pointing in an arbitrary direction, detailed time-dependent expressions are given. For a simple model of a “crustal” layer over a “mantle” half-space, the numerical seismograms in the near- and intermediate-field show some interesting features. These include a prominent group of compressional waves whose radial component is substantial at distances four times the crustal thickness. All the dominant shear arrivals (s, SS, and sSS) are important and show large variations of amplitude as the source depth and receiver distance are varied. Some of the prominent individual generalized rays are shown, and it is found that they can be grouped naturally into families based on the number of interactions with the boundaries. The subdivision into individual generalized rays is useful for analysis and for checks on the numerical stability of the synthetic seismograms. Since the solution is analytic and the numerical evaluation is complete up to any desired time, the results are useful in comparing other approximate methods for the computation of seismograms.


Geophysics ◽  
1982 ◽  
Vol 47 (1) ◽  
pp. 89-99 ◽  
Author(s):  
D. Kao

Magnetotelluric (MT) response is studied for a vertically inhomogeneous earth, where conductivity (or resistivity) varies exponentially with depth as [Formula: see text]. Horizontal electric and magnetic fields in such an inhomogeneous medium are given in terms of modified Bessel functions. Impedance and apparent resistivity are calculated for (1) an inhomogeneous half‐space having conductivity varying exponentially with depth, (2) an inhomogeneous half‐space overlain by a homogeneous layer, and (3) a three‐layer model with the second layer as an inhomogeneous or transitional layer. Results are presented graphically and are compared with those of homogeneous multilayer models. In the case of resistivity increasing exponentially with depth, the results of the above inhomogeneous models are equivalent to those of Cagniard two‐layer models, with [Formula: see text]. In the case of resistivity decreasing exponentially with depth, the homogeneous multilayer approximation depends upon the number of layers and the layer parameters chosen; |Z/ωμ| as a function of frequency is more useful than the apparent resistivity in determining the values of p and [Formula: see text].


1979 ◽  
Vol 69 (6) ◽  
pp. 1693-1714
Author(s):  
F. Abramovici ◽  
J. Gal-Ezer

abstract The time-dependent solution for a multipolar source in a structure consisting of a homogeneous layer over a homogeneous half-space is obtained as a sum of generalized rays. Numerical seismograms are calculated for a horizontal strikeslip and a horizontal dip-slip for a point-source, a finite line-source, and a finite two-dimensional source in the form of a rectangle. For comparison, the displacements in a homogeneous space and half-space are also calculated. The seismograms for finite sources are similar to those for a point-source but show less conspicuous phases, the arriving pulses being wider and less sharp.


2020 ◽  
Vol 222 (1) ◽  
pp. 582-594
Author(s):  
Thomas Forbriger ◽  
Lingli Gao ◽  
Peter Malischewsky ◽  
Matthias Ohrnberger ◽  
Yudi Pan

SUMMARY Other than commonly assumed in seismology, the phase velocity of Rayleigh waves is not necessarily a single-valued function of frequency. In fact, a single Rayleigh mode can exist with three different values of phase velocity at one frequency. We demonstrate this for the first higher mode on a realistic shallow seismic structure of a homogeneous layer of unconsolidated sediments on top of a half-space of solid rock (LOH). In the case of LOH a significant contrast to the half-space is required to produce the phenomenon. In a simpler structure of a homogeneous layer with fixed (rigid) bottom (LFB) the phenomenon exists for values of Poisson’s ratio between 0.19 and 0.5 and is most pronounced for P-wave velocity being three times S-wave velocity (Poisson’s ratio of 0.4375). A pavement-like structure (PAV) of two layers on top of a half-space produces the multivaluedness for the fundamental mode. Programs for the computation of synthetic dispersion curves are prone to trouble in such cases. Many of them use mode-follower algorithms which loose track of the dispersion curve and miss the multivalued section. We show results for well established programs. Their inability to properly handle these cases might be one reason why the phenomenon of multivaluedness went unnoticed in seismological Rayleigh wave research for so long. For the very same reason methods of dispersion analysis must fail if they imply wave number kl(ω) for the lth Rayleigh mode to be a single-valued function of frequency ω. This applies in particular to deconvolution methods like phase-matched filters. We demonstrate that a slant-stack analysis fails in the multivalued section, while a Fourier–Bessel transformation captures the complete Rayleigh-wave signal. Waves of finite bandwidth in the multivalued section propagate with positive group-velocity and negative phase-velocity. Their eigenfunctions appear conventional and contain no conspicuous feature.


Sign in / Sign up

Export Citation Format

Share Document