首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
Summary. The reflection and refraction of general (homogeneous or inhomo-geneous) plane P and type-I S ( SV ) body waves incident on plane boundaries are considered for general linear viscoelastic solids. Reflection—refraction laws, physical characteristics of the waves, and the nature of critical angles are examined in detail at welded boundaries and a free surface. General visco-elasticity with no low-loss approximations predicts that contrasts in intrinsic absorption at boundaries give rise to inhomogeneous reflected and refracted waves with elliptical particle motions, velocities and maximum attenuations that vary with frequency and angle of incidence, energy propagation at speeds and directions different from phase propagation, phase propagation that in general is parallel to the boundary for at most one angle of incidence, and reflection—transmission coefficients dependent on energy flow due to wave interaction. None of these physical characteristics are predicted for waves incident on boundaries that respond instantaneously.  相似文献   

2.
Summary. The propagation of surface waves in a laterally varying medium can be described by representing the wavetrain as a superposition of modal contributions for a reference structure. As the guided waves propagate through a heterogeneous zone the modal coefficients needed to describe the wavetrain vary with position, leading to interconversions between modes and reflection into backward travelling modes. The evolution of the modal terms may be described by a set of first-order differential equations which allow for coupling to both forward and backward travelling waves; the coefficients in these equations depend on the differences between the actual structure and the reference structure. This system is established using the orthogonality properties of the modal eigenfunctions and is valid for SH -waves, P - SV -waves and full anisotropy.
The reflected and transmitted wavefields for a region of heterogeneity can be related to the incident wave by introducing reflection and transmission matrices which connect the modal coefficients in these fields to those in the incident wavetrain. By considering a sequence of models with increasing width of heterogeneity we are able to derive a set of Ricatti equations for the reflection and transmission matrices which may be solved by initial value techniques. This avoids an awkward two-point boundary value problem for a large number of coupled equations. The method is demonstrated for 1 Hz Lg - and Sn -waves in a multilayered model for which there are 19 coupled modes.
The method is applicable to three-dimensional heterogeneity, and we are able to show that the interconversion between Love and Rayleigh waves, in the presence of gradients in seismic properties transverse to the propagation path, leads to a net rate of increase of the transverse components of the seismogram at the expense of the other components.  相似文献   

3.
Seismic waves in a stratified half space   总被引:5,自引:0,他引:5  
Summary. For a buried source in a stratified elastic half space, the surface displacements are calculated by numerical integration of the Fourier–Bessel transform of the response. In the transform space this response is conveniently represented in terms of the reflection and transmission properties of the half space. For a layered medium this procedure avoids all problems associated with growing exponential terms in the evanescent regime. A slightly attenuative medium is assumed, so that the surface wave poles are shifted off the real slowness axis and thus a contour of integration along this axis may be employed. A general point source is represented by an arbitrary moment tensor.
The procedure is illustrated by calculations of three component seismograms including all P , SV and SH contributions for body and surface waves at moderate ranges. For local earthquakes we illustrate the striking effect of focal depth and also show the effect of sedimentary cover on strong ground motion.  相似文献   

4.
Seismic body waves in anisotropic media: synthetic seismograms   总被引:5,自引:0,他引:5  
Summary. Synthetic seismograms and particle motion diagrams are computed for simple, layered Earth models containing an anisotropic layer. The presence of anisotropy couples the P, SV and SH wave motion so that P waves incident on the anisotropic layer from below produce P, SV and small-amplitude SH waves at the surface both the P velocity and the amplitudes of the converted phases vary with azimuth. Significant SH amplitudes may be generated even when the wavelength of the P wave is much greater than the thickness of the anisotropic layer. Incident SV or SH waves may each generate large amplitudes of both SV and SH motion. This strong coupling is largely independent of the degree of velocity anisotropy of the medium. The arrivals from short-period S waves exhibit S-wave splitting, but arrivals from longer period S waves superpose into a modified waveform. This strong coupling does not allow the arrival of separate phases with pure SV and SH polarization except along directions of symmetry where the motion decouples.  相似文献   

5.
Summary. We develop a méthod of reconstructing the elastic paraméters as functions of depth, for a horizontally stratified, isotropic elastic half-space. Unlike previous schemes, which have been able to retrieve the shear wave speed and density from SH seismograms slant stacked at two angles, our méthod makes use of P - SV data at a single stacking paraméter to obtain all three elastic constants. The data required are the elements of the full reflection matrix at the surface, corresponding to measurements of two separate components of the response to two independent sources, one explosive, the other generating shear waves.
In developing this inverse scheme fundamental differences emerge between the acoustic or SH problem, and the coupled P - SV case, the most important being in the nature of the interfacial scattering matrix. We show that it is not possible to make use of the downward reflection data for an interface to determine directly the remaining reflection and transmission coefficients, but that the scattering data may be completed by applying a simple iterative procedure at each interface.
We show the result of applying our inverse scheme to seismograms generated for a six-layered model, including a low-velocity layer. We are able to reconstruct both wave speeds and the density as functions of depth, all quantities being in close agreement with the original model.  相似文献   

6.
Seismic body waves in anisotropic media: propagation through a layer   总被引:3,自引:0,他引:3  
Summary. The square-root energy ratios and pulse shapes are presented for P, SV and SH waves transmitted through a layer of orthorhombic olivine between two isotropic half-spaces. Off incident planes of symmetry, incident P waves generate two small amplitude SH waves (one from each interface), whose amplitudes decrease slowly with increasing period. Incident SV (or SH ) waves can generate large amplitude SH (or SV ) waves which decrease rapidly with increasing period. For incident S waves, many pulses not present in isotropic models are generated, often of large relative amplitude, with many of the transmitted S pulses showing evidence of double arrivals, either in the form of S-wave splitting, or a modification of the shape of the input waveform.  相似文献   

7.
The diffraction of P, S and Rayleigh waves by 3-D topographies in an elastic half-space is studied using a simplified indirect boundary element method (IBEM). This technique is based on the integral representation of the diffracted elastic fields in terms of single-layer boundary sources. It can be seen as a numerical realization of Huygens principle because diffracted waves are constructed at the boundaries from where they are radiated by means of boundary sources. A Fredholm integral equation of the second kind for such sources is obtained from the stress-free boundary conditions. A simplified discretization scheme for the numerical and analytical integration of the exact Green's functions, which employs circles of various sizes to cover most of the boundary surface, is used.
The incidence of elastic waves on 3-D topographical profiles is studied. We analyse the displacement amplitudes in the frequency, space and time domains. The results show that the vertical walls of a cylindrical cavity are strong diffractors producing emission of energy in all directions. In the case of a mountain and incident P, SV and SH waves the results show a great variability of the surface ground motion. These spatial variations are due to the interference between locally generated diffracted waves. A polarization analysis of the surface displacement at different locations shows that the diffracted waves are mostly surface and creeping waves.  相似文献   

8.
Summary. A formulation is derived for calculating the energy division among waves generated by plane waves incident on a boundary between generally anisotropic media. A comprehensive account is presented for P, SV and SH waves incident from an isotropic half-space on an orthorhombic olivine half-space, where the interface is parallel to a plane of elastic symmetry. For comparison, a less anisotropic medium having transverse isotropy with a horizontal axis of symmetry is also considered. The particle motion polarizations of waves in anisotropic medium differ greatly from the polarizations in isotropic media, and are an important diagnostic of the presence of anisotropy. Incident P and SV waves generate quasi- SH waves, and incident SH waves generate quasi- P and quasi- SV waves, often of considerable relative magnitude. The direction of energy transport diverges from the propagation direction.  相似文献   

9.
Summary. A ray-theoretical approach is successfully employed for obtaining the relations between the intrinsic dissipation factors of P - and S -waves, Q α and Q β, and the apparent dissipation factors of toroidal and spheroidal oscillations, Q Tor and Q SphThese are given in quite simple forms and are expressed, to zero-order approximatiom which neglects the effect of reflection of waves at an interface and/or at a free surface), as where α and β are the P - and S -wave velocities and the integrals are evaluated long the surface to surface P - and S -rays.
Some numerical comparisons support the validity of these formulae for higher modes of toroidal and radial oscillations of a realistic earth model.
A first-order approximation reveals that the existence of a discontinuity in an earth induces a systematic fluctuation in normal-mode Q values for any given phase velocity as a function of frequency (or radial mode number).  相似文献   

10.
An explicit analytical formula for the complete elastodynamic Green tensor for homogeneous unbounded weak transversely isotropic media is presented. The formula was derived by analytical calculations of higher-order approximations of the ray series. The ray series is finite and consists of seven non-zero terms. The formula for the Green tensor is complete and correct for the whole frequency range, thus it describes correctly the wavefield at all distances and at all directions including the shear-wave singularity direction. The Green tensor consists of P, SV and SH far-field waves and four coupling waves. Three of them couple P and SV waves, and the fourth wave couples the SV and SH waves. The P-SV coupling waves behave similarly to the near-field waves in isotropy. However, the SV-SH coupling wave, which is called 'shear-wave coupling', behaves exceptionally and it has no analogy in the Green tensor for isotropy. The formula for the elastostatic Green tensor is also derived.  相似文献   

11.
We propose a vertical array analysis method that decomposes complex seismograms into body and surface wave time histories by using a velocity structure at the vertical array site. We assume that the vertical array records are the sum of vertically incident plane P and S waves, and laterally incident Love and Rayleigh waves. Each phase at the surface is related to that at a certain depth by the transfer function in the frequency domain; the transfer function is obtained by Haskell's matrix method, assuming a 1-D velocity structure. Decomposed P , S and surface waves at the surface are estimated from the vertical array records and the transfer functions by using a least-squares method in the frequency domain; their time histories are obtained by the inverse Fourier transform. We carried out numerical tests of this method based on synthetic vertical array records consisting of vertically incident plane P and S waves and laterally incident plane Love and Rayleigh waves. Perfect results of the decomposed P , S , Love and Rayleigh waves were obtained for synthetic records without noise. A test of the synthetic records in which a small amount of white noise was added yielded a reasonable result for the decomposed P , S and surface waves. We applied this method to real vertical array records from the Ashigara valley, a moderate-sized sedimentary valley. The array records from two earthquakes occurring at depths of 123 and 148 km near the array (epicentral distance of about 31 km) exhibited long-duration later phases. The analysis showed that duration of the decomposed S waves was a few seconds and that the decomposed surface waves appeared a few seconds after the direct S -wave arrival and had very long duration. This result indicated that the long-duration later phases were generated not by multireflected S waves, but by basin-induced surface waves.  相似文献   

12.
Summary. A set of recurrence relations which are computationally more efficient than those of the reflection matrix method of Kennett & Kerry is presented for P - and SV -wave generation in a ( n + 1) layered medium. The recurrence relations contain no growing terms and thus provide a stable algorithm for computing complete P and SV synthetic seismograms. Our algorithm requires a fewer algebraic operations for computing the reflectivity and transmissivity coefficients, ranging from 15 per cent less for a source in the half-space to 30 per cent less for a source in the top layer, than the reflection matrix method.  相似文献   

13.
We derive asymptotic formulae for the toroidal and spheroidal eigenfrequencies of a SNREI earth model with two discontinuities, by considering the constructive interference of propagating SH and P-SV body waves. For a model with a smooth solid inner core, fluid outer core and mantle, there are four SH and 10 P-SV ray parameters regimes, each of which must be examined separately. The asymptotic eigenfrequency equations in each of these regimes depend only on the intercept times of the propagating wave types and the reflection and transmission coefficients of the waves at the free surface and the two discontinuities. If the classical geometrical plane-wave reflection and transmission coefficients are used, the final eigenfrequency equations are all real. In general, the asymptotic eigenfrequencies agree extremely well with the exact numerical eigenfrequencies; to illustrate this, we present comparisons for a crustless version of earth model 1066A.  相似文献   

14.
The characteristics of a reflected spherical wave at a free surface are investigated by numerical methods; in particular, the polarization angles and amplitude coefficients of a reflected spherical wave are studied. The classical case of the reflection of a plane P wave from a free surface is revisited in order to establish our terminology, and the classical results are recast in a way which is more suited for the study undertaken. The polarization angle of a plane P wave, for a given angle of incidence, is shown to be 90° minus twice the angle of reflection of the reflected S wave. For a Poisson's ratio less than 1/3, there is a non-normal incident angle for which both amplification coefficients are 2 precisely; for this incident angle the direction of the particle motion at the free surface is also the direction of the incident wave. For a wave emanating from a spherical source, the polarization angle, for all angles of incidence, is always less than, or equal to, the polarization angle of a plane P wave. The vector amplification coefficient of a spherical wave, for all angles of incidence, is always greater than the vector amplification coefficient of a plane P wave. As expected, the results for a spherical wave approach the results for a plane P wave in the far field. Furthermore, there was a good agreement between the theoretical modelling and the numerical modelling using the dynamic finite element method (DFEM).  相似文献   

15.
The Born approximation is applied to the modelling of the propagation of deeply turning longperiod body waves through heterogeneities in the lowermost mantle. We use an exact Green's function for a spherically symmetric earth model that also satisfies the appropriate boundary conditions at internal boundaries and the surface of the earth. The scattered displacement field is obtained by a numerical quadrature of the product of the Green's function, the exciting wavefield and structural perturbations. We study three examples: scattering of longperiod P waves from a plume rising from the coremantle boundary (CMB), generation of longperiod precursors to PKIKP by strong, localized scatterers at the CMB, and propagation of corediffracted P waves through largescale heterogeneities in D". The main results are as follows: (1) the signals scattered from a realistic plume are small with relative amplitudes of less than 2 per cent at a period of 20 s, rendering plume detection a fairly difficult task; (2) strong heterogeneities at the CMB of appropriate size may produce observable longperiod precursors to PKIKP in spite of the presence of a diffraction from the PKP B caustic; (3) corediffracted P  waves ( P diff) are sensitive to structure in D" far off the geometrical ray path and also far beyond the entry and exit points of the ray into and out of D"; sensitivity kernels exhibit ringshaped patterns of alternating sign reminiscent of Fresnel zones; (4) P diff also shows a nonnegligible sensitivity to shear wave velocity in D"; (5) down to periods of 40 s, the Born approximation is sufficiently accurate to allow waveform modelling of P diff through largescale heterogeneities in D" of up to 5 per cent.  相似文献   

16.
Summary. A class of elastic transition zones are modelled by considering a homogeneous half space overlying an inhomogeneous half space with a bounded and monotonically increasing profile for the rigidity modulus and constant Poisson's ratio and density. Reflected P waves due to a compressional point source in the upper half space are studied in the frequency and time domains by means of numerical contour integration in the complex k plane and the Fast Fourier Transform (FFT). Results from the exact fourth-order elasticity theory are compared with those from the approximate decoupled equations for P and SV waves. Agreement is observed between the two theories at high frequencies beyond the caustic range.  相似文献   

17.
18.
Experimental study of shear-waves from shots in anisotropic media   总被引:1,自引:0,他引:1  
Summary. A complex array of vertical-seismic-profiles and reflection surveys in the Taman Peninsula, Krasnodar, using geophones with inclined axes, is used to investigate shear-wave splitting in thick Maikopian clays characterised by pronounced diapirism. The shear waves split into SH (faster) and SV (slower) components and display transverse isotropy with a vertical axis of symmetry with a maximum shear-wave delay of up to 0.5 s . Within the accuracy of the observations there is no azimuthal anisotropy.  相似文献   

19.
Summary. A special class of p th-order autoregressive stochastic models is proposed for seismic records of explosions and earthquakes. These models are based on deterministic formulations proposed by Haskell, Von Seggern & Blandford and Aki. Statistical properties of the estimation procedure are derived and the model applied to seismic records of P waves. The method is shown to have applications in depth and reflection parameter estimation and in event discrimination.  相似文献   

20.
Summary. Attempts are made to estimate the distribution of damping with depth, based on the observed damping of P waves, of free vibrations, and of surface waves and partly on a discussion by D. L. Anderson and R. S. Hart. Results are consistent in order of magnitude, but could easily be wrong by a factor of 2. Comments are added on the need for further work. Stress is laid on the damping of the free nutation, which indicates that α in the modified Lomnitz law is about 0.2 on an average.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号