首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
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'.  相似文献   

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

3.
Generalized Born scattering of elastic waves in 3-D media   总被引：1，自引：0，他引：1  
It is well known that when a seismic wave propagates through an elastic medium with gradients in the parameters which describe it (e.g. slowness and density), energy is scattered from the incident wave generating low-frequency partial reflections. Many approximate solutions to the wave equation, e.g. geometrical ray theory (GRT), Maslov theory and Gaussian beams, do not model these signals. The problem of describing partial reflections in 1-D media has been extensively studied in the seismic literature and considerable progress has been made using iterative techniques based on WKBJ, Airy or Langer type ansätze. In this paper we derive a first-order scattering formalism to describe partial reflections in 3-D media. The correction term describing the scattered energy is developed as a volume integral over terms dependent upon the first spatial derivatives (gradients) of the parameters describing the medium and the solution. The relationship we derive could, in principle, be used as the basis for an iterative scheme but the computational expense, particularly for elastic media, will usually prohibit this approach. The result we obtain is closely related to the usual Born approximation, but differs in that the scattering term is not derived from a perturbation to a background model, but rather from the error in an approximate Green's function. We examine analytically the relationship between the results produced by the new formalism and the usual Born approximation for a medium which has no long-wavelength heterogeneities. We show that in such a case the two methods agree approximately as expected, but that in a media with heterogeneities of all wavelengths the new gradient scattering formalism is superior. We establish analytically the connection between the formalism developed here and the iterative approach based on the WKBJ solution which has been used previously in 1-D media. Numerical examples are shown to illustrate the examples discussed.  相似文献   

4.
Non-linear elastic response of rocks has been widely observed in laboratory, but very few seismic studies are reported in the literature, even though it is the most natural environment where this feature could be observed. Analytic solutions to the non-linear wave propagation phenomena are not readily available, and there is a need to use approximated techniques. It is clear that when a seismic wave propagates through a homogeneous non-linear elastic media, it will be perturbed by the non-linearity. This perturbation can be treated as a source of scattering, spreading the energy of the primary wave in space and time, contributing to the seismic coda. This is in some sense similar to the effect of heterogeneities. The properties of the coda due to the non-linearity depend on the amount of non-linearity and the seismic moment. Using a perturbation approach we calculate the amplitude of the scattered waves, and show that it can describe reasonably well the main features of real seismic codas.  相似文献   

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

6.
Summary. Previous theoretical studies on the transfer function of crustal plane surface layers have been primarily based on applying recursion formulae to the solutions of wave equations. In this way the detailed physical insight of the transfer function is often obscure and it is very difficult to express the transfer function in an explicit form. To show that this situation can be improved by an alternative approach, we demonstrate in this paper with models of one and two surface layers that the explicit transfer function can be derived by summing multiple-reflected rays between interfaces. Since the derived transfer function is related explicitly to some physical parameters, such as wave travel times between interfaces, reflection and transmission coefficients at the interfaces, and attenuation and dispersion of waves in each layer of the model, it is more flexible than the traditional recursion form when applied to different models. This suggests that the ray theory can play an important role in problems of layered media.  相似文献   

7.
Numerical simulation of the propagation of P waves in fractured media   总被引：1，自引：0，他引：1  
We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q −1 of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q −1 of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) −1 in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.  相似文献   

8.
Simulation of SH- and P-SV-wave propagation in fault zones   总被引：3，自引：0，他引：3  
Seismic fault-zone (FZ) trapped waves provide a potentially high-resolution means for investigating FZ and earthquake properties. Seismic waves emitted within and travelling along low-velocity FZ layers may propagate many kilometres within the low-velocity structure associated with the fault. Waveform observation of FZ trapped waves can be modelled in terms of FZ layer velocities, thicknesses and attenuation coefficients. This can greatly improve the resolution of imaged FZ structure and microearthquake locations. At present, broad-band theoretical seismograms are restricted to plane-parallel layers of uniform properties. However, it is not clear how realistic these models are compared with actual fault structures which could, for example, flare outwards near the surface, have irregular boundaries, interior heterogeneities, etc. To address these interpretational uncertainties, we perform finite-difference simulations for irregular FZ geometries and non-uniform material properties within the layers. The accuracy of the numerical solutions are verified by comparison with the analytical solution of Ben-Zion & Aki (1990) for plane-parallel structures. Our main findings are: (1) FZs can widen at the ctustal surface only slightly modifying the trapped waves; (2) velocity variations with depth destroy trapped wave propagation at all wavelengths; (3) FZ trapped waves can be obscured by the presence of a low-velocity surface layer; (4) models with short-scale random structures suggest that trapped waves average out irregular FZ geometries, and hence can be effectively modelled by average-property plane-layered media for the observed range of wavelengths.  相似文献   

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

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

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

12.
We consider the two coupled differential equations of the two radial functions appearing in the displacement components of spheroidal oscillations for a transversely isotropic (TI) medium in spherical coordinates. Elements of the layer matrix have been explicitly written—perhaps for the first time—to extend the use of the Thomson-Haskell matrix method to the derivation of the dispersion function of Rayleigh waves in a transversely isotropic spherical layered earth. Furthermore, an earth-flattening transformation (EFT) is found and effectively used for spheroidal oscillations. The exponential function solutions obtained for each layer give the dispersion function for TI spherical media the same form as that on a flat earth. This has been achieved by assuming that the five elastic parameters involved vary as r p and that the density varies as r p-2, where p is an arbitrary constant and r is the radial distance. A numerical illustration with p = - 2 shows that, in spite of the inhomogeneity assumed within layers, the results for spherical harmonic degree n , versus time period T , obtained here for the Primary Reference Earth Model (PREM), agree well with those obtained earlier by other authors using numerical integration or variational methods. The results for isotropic media derived here are also in agreement with previous results. The effect of transverse isotropy on phase velocity for the first two modes of Rayleigh waves in the period range 20 to 240 s is calculated and discussed for continental and oceanic models.  相似文献   

13.
We investigate the effects of lateral heterogeneities in the upper mantle on the calculation of postglacial land uplift. For the model calculations we use a commercial finite-element code, which enables us to solve the equations governing a layered, isotropic. incompressible, Maxwell-viscoelastic half-space with laterally varying layer thicknesses and physical properties. Following previous investigations performed by Sabadini. Yuen & Portney (1986) and Gasperini & Sabadini (1989), we extend their results using a more realistic loading history and different earth models. We then focus our attention on the question whether lateral heterogeneities in the upper mantle can be modelled correctly using a set of homogeneous earth models. To this end, a comparison of model calculations using both laterally homogeneous and heterogeneous earth models is performed.
We find that lateral heterogeneities in the upper mantle significantly influence the calculated postglacial land uplift. The resolving power of relative sea-level observations for the prescribed lateral heterogeneities used in this study is mainly focused on observations around the load margin and outside the glaciated areas, where differences in predicted land uplift between individual models are large enough to be resolved by observations.
We can qualitatively determine lateral heterogeneities in the upper mantle using a set of laterally homogeneous earth models, if the geological structure, for example a continental margin, is known. However, in order to infer the correct values of lithospheric thickness and asthenospheric viscosity, we need to use laterally heterogeneous models.  相似文献   

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

15.
It has been demonstrated both theoretically and experimentally that the Green's function between two receivers can be retrieved from the cross-correlation of isotropic noise records. Since surface waves dominate noise records in geophysics, tomographic inversion using noise correlation techniques have been performed from Rayleigh waves so far. However, very few numerical studies implying surface waves have been conducted to confirm the extraction of the true dispersion curves from noise correlation in a complicated soil structure. In this paper, synthetic noise has been generated in a small-scale (<1 km) numerical realistic environment and classical processing techniques are applied to retrieve the phase velocity dispersion curves, first step toward an inversion. We compare results obtained from spatial autocorrelation method (SPAC), high-resolution frequency-wavenumber method (HRFK) and noise correlation slantstack techniques on a 10-sensor array. Two cases are presented in the (1–20 Hz) frequency band that corresponds to an isotropic or a directional noise wavefield. Results show that noise correlation slantstack provides very accurate phase velocity estimates of Rayleigh waves within a wider frequency band than classical techniques and is also suitable for accurately retrieving Love waves dispersion curves.  相似文献   

16.
The deployment of temporary arrays of broadband seismological stations over dedicated targets is common practice. Measurement of surface wave phase velocity across a small array and its depth-inversion gives us information about the structure below the array which is complementary to the information obtained from body-wave analysis. The question is however: what do we actually measure when the array is much smaller than the wave length, and how does the measured phase velocity relates to the real structure below the array? We quantify this relationship by performing a series of numerical simulations of surface wave propagation in 3-D structures and by measuring the apparent phase velocity across the array on the synthetics. A principal conclusion is that heterogeneities located outside the array can map in a complex way onto the phase velocities measured by the array. In order to minimize this effect, it is necessary to have a large number of events and to average measurements from events well-distributed in backazimuth. A second observation is that the period of the wave has a remarkably small influence on the lateral resolution of the measurement, which is dominantly controlled by the size of the array. We analyse if the artefacts created by heterogeneities can be mistaken for azimuthal variations caused by anisotropy. We also show that if the amplitude of the surface waves can be measured precisely enough, phase velocities can be corrected and the artefacts which occur due to reflections and diffractions in 3-D structures greatly reduced.  相似文献   

17.
  相似文献   

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

19.
Summary. This paper discusses several aspects of the calculation of theoretical seismograms for two-dimensional inhomogeneous media with the method of Gaussian beams. The most important steps of this method, kinematic and dynamic ray tracing, can be performed very efficiently, if the model cross-section is subdivided into triangles with linear velocity laws. Each Gaussian beam is characterized by a complex beam constant ε which determines its width and phase-front curvature. Various possibilities to choose ε are discussed, including cases where beam properties at the beam endpoint (and not at the beginning) are prescribed; for instance, the beam width at the endpoint can be specified. In such cases the beam constant is a function of the radiation angle at the source, and the decomposition of a cylindrical wave into beams has to take this into account by weighting the beams differently, at least in principle. The exact weight function is derived and shown to be reasonably well approximated by the weight function, corresponding to angle-independent ε Theoretical seismograms are presented for a laterally heterogeneous model of the crust–mantle transition which is characterized by complications in the reflection from the transition and in the refraction from below. These complications are modelled by and large with success. The seismograms, however, depend to a certain extent on the choice of the beam constant. Moreover, according to the reciprocity principle calculations with source and receiver interchanged should have the same results as calculations for the original configuration. In practice this is not so, and the difference increases with the strength of lateral heterogeneities. Hence, for a successful application of Gaussian beams the model should not vary too strongly in lateral direction.  相似文献   

20.
With the help of complexity indices, we quantitatively studied multifractals, frequency distributions, and linear and nonlinear characteristics of geochemical data for exploration of the Daijiazhuang Pb–Zn deposit. Furthermore, we derived productivity differentiation models of elements from thermodynamics and self-organized criticality of metallogenic systems. With respect to frequency distributions and multifractals, only Zn in rocks and most elements except Sb in secondary media, which had been derived mainly from weathering and alluviation, exhibit nonlinear distributions. The relations of productivity to concentrations of metallogenic elements and paragenic elements in rocks and those of elements strongly leached in secondary media can be seen as linear addition of exponential functions with a characteristic weak chaos. The relations of associated elements such as Mo, Sb, and Hg in rocks and other elements in secondary media can be expressed as an exponential function, and the relations of one-phase self-organized geological or metallogenic processes can be represented by a power function, each representing secondary chaos or strong chaos. For secondary media, exploration data of most elements should be processed using nonlinear mathematical methods or should be transformed to linear distributions before processing using linear mathematical methods.  相似文献   

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

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