首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 45 毫秒
1.
Elastic scattered waves from a continuous and heterogeneous layer   总被引：3，自引：0，他引：3  
Elastic scattering from a continuous and laterally unbounded heterogeneous layer has been formulated using the Born approximation. A general solution of the scattered wave equation for the above-stated medium has been given in terms of a Fourier integral over plane waves. Far-field asymptotic expressions for weak elastic scattering by a finite, continuous and inhomogeneous layer have been presented which agree with earlier results. For perturbations of the two elastic parameters and the density having the same form of spatial variation, the spectrum of plane waves scattered from a heterogeneous layer is expressed as a product of an 'elastic scattering factor'and a 'distribution factor'. As in earlier results for small-scale heterogeneity, the scattering pattern depends on various combinations of perturbations of elastic parameters and density. In order to show the general characteristics of the elastic wave scattering, some scattering patterns have been given.  相似文献   

2.
Wave propagation in weakly anisotropic inhomogeneous media is studied by the quasi-isotropic approximation of ray theory. The approach is based on the ray-tracing and dynamic ray-tracing differential equations for an isotropic background medium. In addition, it requires the integration of a system of two complex coupled differential equations along the isotropic ray.
The interference of the qS waves is described by traveltime and polarization corrections of interacting isotropic S waves. For qP waves the approach leads to a correction of the traveltime of the P wave in the isotropic background medium.
Seismograms and particle-motion diagrams obtained from numerical computations are presented for models with different strengths of anisotropy.
The equivalence of the quasi-isotropic approximation and the quasi-shear-wave coupling theory is demonstrated. The quasi-isotropic approximation allows for a consideration of the limit from weak anisotropy to isotropy, especially in the case of qS waves, where the usual ray theory for anisotropic media fails.  相似文献   

3.
Velocity estimation remains one of the main problems when imaging the subsurface with seismic reflection data. Traveltime inversion enables us to obtain large-scale structures of the velocity field and the position of seismic reflectors. However, as the media currently under study are becoming more and more complex, we need to know the finer-scale structures. The problem is that below a certain range of velocity heterogeneities, deterministic methods become difficult to use, so we turn to a probabilistic approach. With this in view, we characterize the velocity field as a random field defined by its first and second statistical moments. Usually, a seismic random medium is defined as a homogeneous velocity background perturbed by a small random field that is assumed to be stationary. Thus, we make a link between such a random velocity medium (together with a simple reflector) and seismic reflection traveltimes. Assuming that the traveltimes are ergodic, we use 2-D seismic reflection geometry to study the decrease in the statistical traveltime fluctuations as a function of the offset (the source–receiver distance). Our formulae are based on the Rytov approximation and the parabolic approximation for acoustic waves. The validity and the limits are established for both of these approximations in statistically anisotropic random media. Finally, theoretical inversion procedures are developed for the horizontal correlation structure of the velocity heterogeneities for the simplest case of a horizontal reflector. Synthetic seismograms are then computed (on particular realizations of random media) by simulating scalar wave propagation via finite difference algorithms. There is good agreement between the theoretical and experimental results.  相似文献   

4.
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.  相似文献   

5.
Finite-frequency sensitivity kernels for head waves   总被引：2，自引：0，他引：2  
Head waves are extremely important in determining the structure of the predominantly layered Earth. While several recent studies have shown the diffractive nature and the 3-D Fréchet kernels of finite-frequency turning waves, analogues of head waves in a continuous velocity structure, the finite-frequency effects and sensitivity kernels of head waves are yet to be carefully examined. We present the results of a numerical study focusing on the finite-frequency effects of head waves. Our model has a low-velocity layer over a high-velocity half-space and a cylindrical-shaped velocity perturbation placed beneath the interface at different locations. A 3-D finite-difference method is used to calculate synthetic waveforms. Traveltime and amplitude anomalies are measured by the cross-correlation of synthetic seismograms from models with and without the velocity perturbation and are compared to the 3-D sensitivity kernels constructed from full waveform simulations. The results show that the head wave arrival-time and amplitude are influenced by the velocity structure surrounding the ray path in a pattern that is consistent with the Fresnel zones. Unlike the 'banana–doughnut' traveltime sensitivity kernels of turning waves, the traveltime sensitivity of the head wave along the ray path below the interface is weak, but non-zero. Below the ray path, the traveltime sensitivity reaches the maximum (absolute value) at a depth that depends on the wavelength and propagation distance. The sensitivity kernels vary with the vertical velocity gradient in the lower layer, but the variation is relatively small at short propagation distances when the vertical velocity gradient is within the range of the commonly accepted values. Finally, the depression or shoaling of the interface results in increased or decreased sensitivities, respectively, beneath the interface topography.  相似文献   

6.
A Gaussian correlation function characterizes smoothly heterogeneous media, while real heterogeneities in the Earth are often non-Gaussian in nature. Using the Born approximation, mean square amplitudes of the scattered waves have been derived for an elastic media characterized by the Von Karman correlation function. Heterogeneities with different power laws can be defined by the Von Karman correlation function. The sensitivity of fore- and backscattering to heterogeneities with different scales and properties (that is velocity and impedance) is discussed in this paper. The analytical expression for total scattered energy for the incident P waves is also derived for a random medium having the Von Karman correlation function. We find that at high frequencies, the scattered power of converted waves is a function of frequency. In the case of codawave excitation by local earthquakes, which must be handled by the full elastic-wave theory, we can define any type of inhomogeneity by the Von Karman correlation function. It also supports the idea that the lithosphere might have multiple-scale inhomogeneities.  相似文献   

7.
In this study, we test the adequacy of 2-D sensitivity kernels for fundamental-mode Rayleigh waves based on the single-scattering (Born) approximation to account for the effects of heterogeneous structure on the wavefield in a regional surface wave study. The calculated phase and amplitude data using the 2-D sensitivity kernels are compared to phase and amplitude data obtained from seismic waveforms synthesized by the pseudo-spectral method for plane Rayleigh waves propagating through heterogeneous structure. We find that the kernels can accurately predict the perturbation of the wavefield even when the size of anomaly is larger than one wavelength. The only exception is a systematic bias in the amplitude within the anomaly itself due to a site response.
An inversion method of surface wave tomography based on the sensitivity kernels is developed and applied to synthesized data obtained from a numerical simulation modelling Rayleigh wave propagation over checkerboard structure. By comparing recovered images to input structure, we illustrate that the method can almost completely recover anomalies within an array of stations when the size of the anomalies is larger than or close to one wavelength of the surface waves. Surface wave amplitude contains important information about Earth structure and should be inverted together with phase data in surface wave tomography.  相似文献   

8.
Summary. Scattering of surface waves by lateral heterogeneities is analysed in the Born approximation. It is assumed that the background medium is either laterally homogeneous, or smoothly varying in the horizontal direction. A dyadic representation of the Green's function simplifies the theory tremendously. Several examples of the theory are presented. The scattering and mode conversion coefficients are shown for scattering of surface waves by the root of an Alpine-like crustal structure. Furthermore a 'great circle theorem'in a plane geometry is derived. A new proof of Snell's law is given for surface wave scattering by a quarter-space. It is shown how a stationary phase approximation can be used to simplify the Fourier synthesis of the scattered wave in the time domain. Finally a procedure is suggested to do 'surface wave holography'.  相似文献   

9.
  相似文献   

10.
An efficient inverse scattering method is developed for imaging near-surface heterogeneities using scattered surface waves. Three dimensional elastodynamic wave propagation and scattering in a laterally invariant embedding medium is considered. The Born Approximation is used and the scattered wavefield is expressed as a domain type integral representation. The computation time of Green's tensor elements is reduced by considering the radial symmetry of the medium. The method is validated by numerical tests. Ultrasonic laboratory data obtained from a scale model experiment are used for imaging the near-surface inhomogeneities caused by an epoxy-filled hole in the surface of an aluminum block. Both synthetic and the scale model tests show that the location, the actual density contrast and the depth of the inhomogeneities are reasonably well estimated.  相似文献   

11.
Polarization anomaly of Love waves caused by lateral heterogeneity   总被引：1，自引：0，他引：1  
We calculate surface waves propagating in a laterally heterogeneous structure beneath the Kuril trench, where significant Love-wave polarization anomalies, called quasi-Love waves, are generated. Since 3-D wave propagation in the two-dimensionally heterogeneous structure can be assumed, we apply the 2.5-D finite difference method to the surface-wave calculations. The calculations show that a velocity contrast of 7 per cent at depths of less than 210 km beneath the Kuril trench cannot generate quasi-Love waves, and that an unlikely contrast of 20 per cent is required to generate clear quasi-Love waves. The possible cause of the quasi-Love waves inferred from previous studies on coupled free oscillations is a lateral variation in azimuthal anisotropy. The lateral variation in azimuthal anisotropy beneath the Kuril trench suggests a change in the mantle flow induced by the subducting slab.  相似文献   

12.
We present the results of Rayleigh wave and Love wave phase velocity tomography in the western United States using ambient seismic noise observed at over 250 broad-band stations from the EarthScope/USArray Transportable Array and regional networks. All available three-component time-series for the 12-month span between 2005 November 1 and 2006 October 31 have been cross-correlated to yield estimated empirical Rayleigh and Love wave Green's functions. The Love wave signals were observed with higher average signal-to-noise ratio (SNR) than Rayleigh wave signals and hence cannot be fully explained by the scattering of Rayleigh waves. Phase velocity dispersion curves for both Rayleigh and Love waves between 5 and 40 speriod were measured for each interstation path by applying frequency–time analysis. The average uncertainty and systematic bias of the measurements are estimated using a method based on analysing thousands of nearly linearly aligned station-triplets. We find that empirical Green's functions can be estimated accurately from the negative time derivative of the symmetric component ambient noise cross-correlation without explicit knowledge of the source distribution. The average traveltime uncertainty is less than 1 s at periods shorter than 24 s. We present Rayleigh and Love wave phase speed maps at periods of 8, 12, 16,and 20 s. The maps show clear correlations with major geological structures and qualitative agreement with previous results based on Rayleigh wave group speeds.  相似文献   

13.
Summary. In order to separate the scattering effect from the intrinsic attenuation, we need a multiple scattering model for seismic wave propagation in random heterogeneous media. In this paper, we apply radiative transfer theory to seismic wave propagation and formulate in the frequency domain the energy density distribution in space for a point source. We consider the cases of isotropic scattering and strong forward scattering. Some numerical examples are shown. It is seen that the energy density–distance curves have quite different shapes depending on the values of medium seismic albedo B 0=ηs/(ηs+ηa) where ηs is the scattering coefficient and ηa is the absorption coefficient of the medium. For a high albedo ( B > 0.5) medium, the energy–distance curve is of arch shape and the position of the peak is a function of the extinction coefficient of the medium ηe=ηs+ηa. Therefore it is possible to separate the scattering effect and the absorption based on the measured energy density distribution curves.  相似文献   

14.
We investigate the sensitivity of finite-frequency body-wave observables to mantle anisotropy based upon kernels calculated by combining adjoint methods and spectral-element modelling of seismic wave propagation. Anisotropy is described by 21 density-normalized elastic parameters naturally involved in asymptotic wave propagation in weakly anisotropic media. In a 1-D reference model, body-wave sensitivity to anisotropy is characterized by 'banana–doughnut' kernels which exhibit large, path-dependent variations and even sign changes. P -wave traveltimes appear much more sensitive to certain azimuthally anisotropic parameters than to the usual isotropic parameters, suggesting that isotropic P -wave tomography could be significantly biased by coherent anisotropic structures, such as slabs. Because of shear-wave splitting, the common cross-correlation traveltime anomaly is not an appropriate observable for S waves propagating in anisotropic media. We propose two new observables for shear waves. The first observable is a generalized cross-correlation traveltime anomaly, and the second a generalized 'splitting intensity'. Like P waves, S waves analysed based upon these observables are generally sensitive to a large number of the 21 anisotropic parameters and show significant path-dependent variations. The specific path-geometry of SKS waves results in favourable properties for imaging based upon the splitting intensity, because it is sensitive to a smaller number of anisotropic parameters, and the region which is sampled is mainly limited to the upper mantle beneath the receiver.  相似文献   

15.
When full 3-D modelling is too costly or cumbersome, computations of 3-D elastic wave propagation in laterally heterogeneous, multilayered 2-D geological structures may enhance considerably our ability to predict strong ground motion for seismological and engineering purposes. Towards this goal, we extend the method based on the combination of the thin-layer finite-element and boundary-element methods (TLFE-BEM) and calculate windowend f - k spectra of the 3-D wavefield. The windowed f - k spectra are spatially localized spectra from which the local properties of the wavefield can be extracted. The TLFE-BEM is particularly suited for calculating the complete wavefield where surface waves are dominant in multilayered media. The computations are performed in the frequency domain, providing the f - k spectra directly. From the results for the 3-D wavefield excited by a point source in a 2-D multilayered, sloped structure, it can be said that the phase velocity of the fundamental-mode Rayleigh wave in a laterally heterogeneous multilayered medium, estimated from the windowed f - k spectra, varies with the location of the point source. For the model calculated in this article, the phase velocity varies between the value for the flat layered structure of the thick-layer side and that for the structure just under the centre of the window. The exact subsurface structure just under the centre of an array in a laterally heterogeneous medium cannot be obtained if we use the f - k spectral analysis assuming a flat layered structure.  相似文献   

16.
Summary. General expressions are derived for the kernels of the set of integral equations that relates the spectral matrix of the six components of a random electromagnetic wave field in a magnetoplasma to the wave distribution function for the field. The dependence of the kernels on wave-normal direction is examined, with particular reference to the propagation of very low-frequency waves in the whistler mode.  相似文献   

17.
Scattering of surface waves modelled by the integral equation method   总被引：1，自引：0，他引：1  
The integral equation method is used to model the propagation of surface waves in 3-D structures. The wavefield is represented by the Fredholm integral equation, and the scattered surface waves are calculated by solving the integral equation numerically. The integration of the Green's function elements is given analytically by treating the singularity of the Hankel function at   R = 0  , based on the proper expression of the Green's function and the addition theorem of the Hankel function. No far-field and Born approximation is made. We investigate the scattering of surface waves propagating in layered reference models imbedding a heterogeneity with different density, as well as Lamé constant contrasts, both in frequency and time domains, for incident plane waves and point sources.  相似文献   

18.
The Kirchhoff (or tangent plane) approximation, derived from the theoretically complete Kirchhoff–Helmholtz integral representation for the seismic wavefield, has been used extensively for the analysis of seismic-wave scattering from irregular interfaces; however, the accuracy of this method for curved interfaces has not been rigorously established. This paper describes an efficient Kirchhoff algorithm to simulate scattered waves from an arbitrarily curved interface in an elastic medium. Synthetic seismograms computed using this algorithm are compared with exact synthetics computed using analytical formulae for scattering of plane P waves by a spherical elastic inclusion. A windowing technique is used to remove strong internal reverberations from the analytical solution. Although the Kirchhoff method tends to underestimate the total scattering intensity, the accuracy of the approximation improves with increasing value of the wavenumber-radius product, kR . The arrival times and pulse shapes of primary reflections from the sphere are well approximated using the Kirchhoff approach regardless of curvature of the scattering surface, but the amplitudes are significantly underestimated for kR ≤ 5. The results of this work provide some new guidelines to assess the accuracy of Kirchhoff-synthetic seismograms for curved interfaces.  相似文献   

19.
We test the feasibility of using Green's functions extracted from records of ambient seismic noise to monitor temporal changes in the Earth crust properties by repeated measurements at regional distances. We use about 11 yr of continuous recordings to extract surface waves between three pairs of stations in California. The correlations are computed in a moving 1-month window and we analyse the temporal evolution of measured interstation traveltimes. The comparison of the arrival times in the positive and negative correlation time of Rayleigh and Love waves allows us to separate time-shifts associated with any form of physical change in the medium, those resulting from clock drift or other instrumental errors, and those due to change in the localization of the noise sources. This separation is based on the principle of time symmetry. When possible, we perform our analysis in two different period bands: 5–10 and 10–20 s. The results indicate that significant instrumental time errors (0.5 s) are present in the data. These time-shifts can be measured and tested by closure relation and finally corrected independently of any velocity model. The traveltime series show a periodic oscillation that we interpret as the signature of the seasonal variation of the region of origin of the seismic noise. Between 1999 and 2005, the final arrival time fluctuations have a variance of the order of 0.01 s. This allows us to measure interstation traveltimes with errors smaller than 0.3 per cent of the interstation traveltime and smaller than 1 per cent of the used wave period. This level of accuracy was not sufficient to detect clear physical variation of crustal velocity during the considered 11 yr between the three stations in California. Such changes may be more easily detectable when considering pairs of stations more closely located to each other and in the vicinity of tectonically active faults or volcanoes.  相似文献   

20.
Summary. Using a single scattering approximation, we derive equations for the scattering attenuation coefficients of P- and S -body waves. We discuss our results in the light of some recent energy renormalization approaches to seismic wave scattering. Practical methods for calculating the scattering attenuation coefficients for various earth models are emphasized. The conversions of P - to S -waves and S- to P -waves are included in the theory. The earth models are assumed to be randomly inhomogeneous, with their properties known only through their average wavenumber power spectra. We approximate the power spectra with piecewise constant functions, each segment of which contributes to the net, frequency-dependent, scattering attenuation coefficient. The smallest and largest wavenumbers of a segment can be plotted along with the wavevectors of the incident and scattered waves on a wavenumber diagram. This diagram gives a geometric interpretation for the frequency behaviour associated with each spectral segment, including a 'transition' peak that is due entirely to the wavenumber limits of the segment. For regions of the earth where the inhomogeneity spectra are concentrated in a band of wavenumbers, it should be possible to observed such a peak in the apparent attenuation of seismic waves. We give both the frequency and distance limits on the accuracy of the theoretical results.  相似文献   

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

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