首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 49 毫秒
1.
Broad-band data from South American earthquakes recorded by Californian seismic networks are analysed using a newly developed seismic wave migration method—the slowness backazimuth weighted migration (SBWM). Using the SBWM, out-of-plane seismic P -wave reflections have been observed. The reflection locations extend throughout the Earth's lower mantle, down to the core–mantle boundary (CMB) and coincide with the edges of tomographically mapped high seismic velocities. Modelling using synthetic seismograms suggests that a narrow (10–15 km) low- or high-velocity lamella with about 2 per cent velocity contrast can reproduce the observed reflected waveforms, but other explanations may exist. Considering the reflection locations and synthetic modelling, the observed out-of-plane energy is well explained by underside reflections off a sharp reflector at the base of the subducted lithosphere. We also detect weaker reflections corresponding to the tomographically mapped top of the slab, which may arise from the boundary between the Nazca plate and the overlying former basaltic oceanic crust. The joint interpretation of the waveform modelling and geodynamic considerations indicate mass flux of the former oceanic lithosphere and basaltic crust across the 660 km discontinuity, linking processes and structure at the top and bottom of the Earth's mantle, supporting the idea of whole mantle convection.  相似文献   

2.
Summary. Conventional migration of deep seismic reflection data often produces disappointingly poor results even when the original unmigrated data are of high quality and are relatively noise free. Surprisingly, deep data often subjectively appear to be best migrated at velocities which are up to 50% less than appropriate interval velocities derived from crustal refraction experiments or directly from stacking velocities. The explanation for this behaviour is that near surface features distort and attenuate the seismic wave field and produce apparent discontinuities in deep reflections. Since these discontinuities are spurious they are not associated with the appropriate large diffractions which real discontinuities at depth would produce. During the process of migration reflections are invented in order to cancel out the missing diffractions thereby producing a "smiley" section which appears overmigrated. Since the lateral extent of individual smiles increases with increasing migration velocity, increasing two-way-time, and increasing seismic wavelength, the effect" is almost unnoticeable on conventional shallow seismic data but is overwhelming for deep crustal data.  相似文献   

3.
Summary. For a smooth earth model, observations of a set of high-frequency toroidal modes at fixed slowness yield only a single piece of information, the tau value for that slowness. In this note, a procedure for obtaining the shear velocity structure from free oscillation data for an earth model with velocity discontinuities is developed, based on the method of tau inversion. The information content of the high-frequency modes is greater in this case, and the nature and depths of the discontinuities may be deduced. It is shown, for the real Earth, that the tau values obtained from free oscillation data are affected significantly by the presence of the Moho, but a simple iterative scheme may be used to remove this contamination. Brune's method of deducing mode frequencies from body wave pulses is shown to produce significant errors for a model with a pronounced Moho discontinuity, and the same iterative scheme may also be employed to correct for this effect.  相似文献   

4.
This paper presents a method to invert underside-reflection ( P d P or S d S arrivals) data for lateral depth variations of upper-mantle discontinuities, combining traveltime and amplitude data. The method greatly improves the resolution of small-scale undulations obtained by existing imaging methods and does not suffer from the long-wavelength biases that are likely to be present in currently available models. Existing inversion methods account for the large size of the Fresnel zone of underside reflections, but not for its complexity, arising from the mini-max traveltime nature of PP- and SS -related waves. This neglect results in long-wavelength artefacts from small-scale undulations of the discontinuities, obscuring true long-wavelength depth variations. The inversion method presented in this paper uses a complex-valued sensitivity kernel, derived from the representation of underside reflections through a Kirchhoff integral formulation. The sensitivity kernel accounts for the varying sensitivity of the waveforms to discontinuity structure over the Fresnel zone. The method is applied to a large, synthetic data set. The data set consists of P d P amplitudes and traveltimes. The results show that the new inversion method resolves depth variations on a lateral scale that is smaller than the size of the Fresnel zone of individual underside reflections (but larger than the dominant wavelength), retaining the resolution of large-scale variations. The results presented here suggest that the discontinuity depth variations induced by slab penetration of the 670 discontinuity could be resolved by current broad-band P 670 P data sets.  相似文献   

5.
Summary. The asymptotic properties of spheroidal mode dispersion at high frequency for fixed phase velocity are related to the intercept times τβ( p ) for P and S waves. If the mode eigenfrequency and the ratio of horizontal to vertical displacement at the surface for the mode are known τα( p ) and τβ( p ) may be separately estimated. If discontinuities exist in the velocity model then 'solotone' effects occur, in frequency at fixed slowness, and in τα( p ), τβ( p ) estimated from the mode dispersion as a function of slowness. The coupling of P and S waves in the spheroidal modes means that the interaction of P waves with upper-mantle discontinuities affects also the estimates of the S wave τβ( p ) values for which the corresponding turning points lie in the lower mantle. The asymptotic formalism also shows that sharp pulses formed by superposition of spheroidal modes correspond to multiple PS reflections.
A study of τα( p ), τβ( p ) estimates derived from spheroidal modes with periods from 45–50s, calculated for model 1066B, shows that even in the presence of strong upper-mantle discontinuities the errors in intercept time are only about one-tenth of a period. The asymptotic properties may there-for provide a useful means of estimating intercept times from modes with a few seconds period as a supplement to travel-time methods.  相似文献   

6.
Broad-band P - and S -waves from earthquakes in South America recorded at Californian network stations are analysed to image lateral variations of the D"-discontinuity beneath the Cocos plate. We apply two array processing methods to the data set: a simplified migration method to the P -wave data set and a double-array method to both the P - and S -wave data sets, allowing us to compare results from the two methods. The double-array method images a dipping reflector at a depth range from 2650 to 2700 km in the southern part of the study area. We observe a step-like topography of 100 km to a shallower reflector at about 2600 km depth to the north, as well as evidence for a second (deeper) reflector at a depth range from 2700 to 2750 km in the north. Results from the simplified migration agree well with those from the double-array method, similarly locating a large step in reflector depth in a similar location (about 2650 km depth in the south and about 2550 km in the north) as well as the additional deeper reflector at the depth of about 2750 km in the north. Waveform modelling of the reflected waves from both methods suggests a positive velocity contrast for S waves, but a negative velocity contrast for P waves for the upper reflector in agreement with predictions from mineral physical calculations for a post-perovskite phase transition. The data also show some evidence for the existence of another deeper reflector that could indicate a double intersection of the geotherm with the post-perovskite stability field, that is, the back-transformation of post-perovskite to perovskite close to the core–mantle boundary.  相似文献   

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

8.
Summary . Synthetic seismograms represented by integrals generally display signals associated with the limits of integration. Sometimes these 'end-point' contributions are spurious (e.g. in the WKBJ seismogram) and sometimes they are the main physical interest (e.g. the Kirchhoff integral for an edge). The end-point contributions may be asymptotically approximated using integration by parts or Laplace's method and it may then be possible to reduce them if desired. We describe examples in the WKBJ seism ogram for reflected or transmitted waves in homogeneous layers and for turning waves. We also study signals due to discontinuities in reflection coefficients, by partitioning the real slowness integral so that the discontinuities lie at end points. Examples are the head wave, which is a physically correct signal, and spurious diffractions caused by using plane-wave coefficients for grazing rays in the WKBJ seismogram.  相似文献   

9.
Summary. The study of the asymptotic behaviour of eigenfrequencies of torsional modes of the Earth is of interest in the problem of the existence of discontinuities in the Earth's interior and the determination of their depth and scale, since the solotone effect, which is a persistent oscillatory cohponent in the asymptotic overtone structure, is very sensitive to the presence of discontinuities. The asymptotic behaviour of torsional eigenfrequencies of the second order differential equation for the Earth's free oscillations can be compared with that obtained from eigenfrequencies evaluated from synthetic SH-wave seismograms by Brune's phase correlation method, using various earth models. The solotone effect that appears in the former for discontinuous models can be explained in terms of multiple reflections fram the discontinuities, and can be reconstructed from synthetic SH -wave pulses arising from these reflections. Its properties vary systematically with the depth and the scale of discontinuities and can be superposed for several discontinuities.  相似文献   

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

11.
We study properties of the energy-flux vector and other related energy quantities of homogeneous and inhomogeneous time-harmonic P and S plane waves, propagating in unbounded viscoelastic anisotropic media, both analytically and numerically. We propose an algorithm for the computation of the energy-flux vector, which can be used for media of unrestricted anisotropy and viscoelasticity, and for arbitrary homogeneous or inhomogeneous plane waves. Basic part of the algorithm is determination of the slowness vector of a homogeneous or inhomogeneous wave, which satisfies certain constraints following from the equation of motion. Approaches for determination of a slowness vector commonly used in viscoelastic isotropic media are usually difficult to use in viscoelastic anisotropic media. Sometimes they may even lead to non-physical solutions. To avoid these problems, we use the so-called mixed specification of the slowness vector, which requires, in a general case, solution of a complex-valued algebraic equation of the sixth degree. For simpler cases, as for SH waves propagating in symmetry planes, the algorithm yields simple analytic solutions. Once the slowness vector is known, determination of energy flux and of other energy quantities is easy. We present numerical examples illustrating the behaviour of the energy-flux vector and other energy quantities, for homogeneous and inhomogeneous plane P , SV and SH waves.  相似文献   

12.
We present a split-step complex Padé-Fourier migration method based on the one-way wave equation. The downward-continuation operator is split into two downward-continuation operators: one operator is a phase-shift operator and the other operator is a finite-difference operator. A complex treatment of the propagation operator is applied to mitigate inaccuracies and instabilities due to evanescent waves. It produces high-quality images of complex structures with fewer numerical artefacts than those obtained using a real approximation of a square-root operator in the one-way wave equation. Tests on zero-offset data from the SEG/EAGE salt data show that the method improves the image quality at the cost of an additional 10 per cent computational time compared to the conventional Fourier finite-difference method.  相似文献   

13.
High noise levels hamper teleseismic shear wave splitting measurements, which bandpass filtering does not always help. To investigate how robust splitting measurements are to noise, we analysed a set of synthetic records with known splitting parameters and added fixed levels of noise. In the presence of weak anisotropy, single-waveform splitting measurements are unreliable when operating with noisy data sets. A practical rule in terms of S/N ratio and splitting delay time parameters is that splitting is confidently detectable at S/N > 8, regardless of the wave's original polarization orientation. However, for the evidence of weak anisotropy to be detectable and measurable at an S/N value of 4, the backazimuth separation of the phases from the fast polarization direction needs to be higher than 20°. Stacks of individual measurements consistently yield reliable results down to S/N values of 4. Applying stacking to data from DSB (Dublin, Ireland), the fast polarization direction φ and lag time δt are 58° and 0.95  s. This orientation reflects surface trends of deformation in the area, as found elsewhere in the UK. Our result thus reinforces the proposed model that the detected anisotropy in the British Isles originates from lithospheric coherent deformation preserved from the last main tectonic episode.  相似文献   

14.
当前用来计算分形维数的原始数据通常是栅格图像,但是现有的土地利用数据大多以矢量形式存放,若将其转为栅格形式后计算分维数,则将面临基于栅格图像所得到的分形维数的计算精度与像元尺寸以及图像自身大小有较大关系,并且在对栅格图像进行不规则形状提取时容易造成数据丢失等问题。鉴于此,本文利用Visual C#完成了基于Windows平台的矢量数据分形维数计算程序设计,并利用该程序对20世纪90年代三个时期的秦岭中段和鄂东南地区土地利用类型的分维数进行了求算。结果表明该程序是一种方便、实用而且精确的分维数计算方法。此外,从三期各类土地利用类型的稳定性指数来看,在秦岭中段地区表现为：居民点及工矿用地>未利用土地>水田>旱地>草地>林地>水域,而鄂东南地区则表现为：水田>未利用土地>居民点及工矿用地>旱地>林地>草地>水域,这说明在上述两个地区,林地、水域、草地的稳定性都相对较差,它们的潜在变化趋势也较大。  相似文献   

15.
A crustal seismic velocity model for the UK, Ireland and surrounding seas   总被引：1，自引：0，他引：1  
A regional model of the 3-D variation in seismic P -wave velocity structure in the crust of NW Europe has been compiled from wide-angle reflection/refraction profiles. Along each 2-D profile a velocity–depth function has been digitised at 5 km intervals. These 1-D velocity functions were mapped into three dimensions using ordinary kriging with weights determined to minimise the difference between digitised and interpolated values. An analysis of variograms of the digitised data suggested a radial isotropic weighting scheme was most appropriate. Horizontal dimensions of the model cells are optimised at 40 × 40 km and the vertical dimension at 1 km. The resulting model provides a higher resolution image of the 3-D variation in seismic velocity structure of the UK, Ireland and surrounding areas than existing models. The construction of the model through kriging allows the uncertainty in the velocity structure to be assessed. This uncertainty indicates the high density of data required to confidently interpolate the crustal velocity structure, and shows that for this region the velocity is poorly constrained for large areas away from the input data.  相似文献   

16.
There are three types of surfaces which are used for studying wave propagation in anisotropic media: normal surfaces, slowness surfaces and wave surfaces. Normal surfaces and slowness surfaces have been researched in detail. Wave surfaces are the most complicated and comparatively poorly known compared with the other two. Areas of complicated geometrical structure of the wave surfaces are located in the vicinity of conical acoustic axes. There is an elliptical hole on the quick shear wave surface and complicated folds and cusps on the slow shear wave surface. Decomposition of the slow shear wave surface into smooth sheets is used for the study of its geometrical structure. Complexity of shear wave surfaces can be expressed by the number of waves corresponding to a fixed ray. An original approach to the calculation of wave normals depending on ray direction is presented.  相似文献   

17.
A massive vapour cloud explosion occurred at the Buncefield fuel depot near Hemel Hempstead, UK, in the morning of 2005 December 11. The explosion was the result of an overflow from one of the storage tanks with the release of over 300 tons of petrol and generating a vapour cloud that spread over an area of 80 000 m2, before being ignited. Considerable damage was caused in the vicinity of the explosion and a total of 43 people were injured. The explosion was detected by seismograph stations in the UK and the Netherlands and by infrasound arrays in the Netherlands. We analysed the seismic recordings to determine the origin time of 06:01:31.45 ±0.5 s (UTC) from P -wave arrival times. Uncertainties in determination of origin time from acoustic arrival times alone were less than 10 s. Amplitudes of P -, Lg and primary acoustic waves were measured to derive decay relationships as function of distance. From the seismic amplitudes we estimated a yield of 2–10 tons equivalent to a buried explosion. Most seismic stations recorded primary and secondary acoustic waves. We used atmospheric ray tracing to identify the various travel paths, which depend on temperature and wind speed as function of altitude, leading to directional variation. Refracted waves were observed from the troposphere, stratosphere and thermosphere with a good match between observed and calculated traveltimes. The various wave types were also identified through array processing, which provides backazimuth and slowness, of recordings from an infrasound array in the Netherlands. The amplitude of stratospheric refracted acoustic waves recorded by the array microbarometers was used to estimate a yield of 21.6 (±5) tons TNT equivalent. We have demonstrated through joint seismo-acoustic analysis of the explosion that both the seismic velocity model and the atmospheric model are sufficient to explain the observed traveltimes.  相似文献   

18.
The zero-lag cross-correlation technique, used for array analysis in the hypothesis of plane waves, has been modified to allow the wave front to be circular. Synthetic tests have been performed to check the capability of the method, which returns the input test data when the source–array distances are not greater than two or three times the array aperture. For this distance range the method furnishes an estimate of the apparent velocity and the epicentral coordinates of the source. For more distant sources the method becomes equivalent to that based on the planar-wave approximation, which gives an estimate of the backazimuth to the source and the apparent velocity. The method has been applied to seismic data recorded at the active volcano located at Deception Island, Antarctica. 35 volcanic long-period events occurring in a small swarm were selected. Results show that the epicentres are close to the array (between 0.4 and 2 km) and aligned in a SW direction, in agreement with one of the main directions of the fracture system of Deception volcano.  相似文献   

19.
Array techniques are particularly well‐suited for detecting and quantifying the complex seismic wavefields associated with volcanic activity such as volcanic tremor and long‐period events. The methods based on the analysis of the signal in the frequency domain, or spectral methods, have the main advantages of both resolving closely spaced sources and reducing the necessary computer time, but may severely fail in the analysis of monochromatic, non‐stationary signals. Conversely, the time‐domain methods, based on the maximization of a multichannel coherence estimate, can be applied even for short‐duration pulses. However, for both the time and the frequency domain approaches, an exhaustive definition of the errors associated with the slowness vector estimate is not yet available. Such a definition become crucial once the slowness vector estimates are used to infer source location and extent. In this work we develop a method based on a probabilistic formalism, which allows for a complete definition of the uncertainties associated with the estimate of frequency–slowness power spectra from measurement of the zero‐lag cross‐correlation. The method is based on the estimate of the theoretical frequency–slowness power spectrum, which is expressed as the convolution of the true signal slowness with the array response pattern. Using a Bayesian formalism, the a posteriori probability density function for signal slowness is expressed as the difference, in the least‐squares sense, between the model spectrum and that derived from application of the zero‐lag cross‐correlation technique. The method is tested using synthetic waveforms resembling the quasi‐monochromatic signals often associated with the volcanic activity. Examples of application to data from Stromboli volcano, Italy, allow for the estimate of source location and extent of the explosive activity.  相似文献   

20.
  相似文献   

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

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