首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 421 毫秒
1.
Shear-coupled PL     
Summary. Observed teleseismic shear-coupled PL -waves ( SPL ) display a variety of waveforms depending on factors such as source depth, source type and velocity structure. Using a WKBJ spectral method for SV -wave propagation, synthetic seismograms of SPL are produced to examine the factors important in SV and SPL excitation. Results show that SPL is preferentially excited by shallow sources compared to deep sources. This is due to large source area reverberations which consequently leak as SV -waves into the mantle. Interaction at the receiver area then sets up the classic prograde elliptical motion by which SPL can be identified. This is in accordance with the teleseismic observations and indicates that most previous models for the propagation of SPL were not appropriate for shallow source since emphasis was placed on wave interactions occurring only near the receiver.  相似文献   

2.
Summary. Moment tensor inversion methods can be applied with success in the determination of source properties of simple earthquakes. However, these methods utilize the assumption of a point source, which is inadequate for modelling many complicated, shallow earthquakes. For complex earthquakes, an inversion using finite faulting models is desirable but the number of parameters involved requires that a good starting model be found or that independent constraints be placed on some of the parameters. A method is presented for low-pass filtering both the data and Green's functions, passing only signals with wavelengths greater than the dimension of the entire fault. The filter tends to smooth complications in the waveforms and allows application of the point source moment tensor inversion. This method is applied to body waves from the 1978 Thessaloniki, Greece, earthquake, the 1971 San Fernando earthquake and to a multiple-point source synthetic model of the San Fernando event. For the Thessaloniki event, although a multiple-source mechanism has been suggested, inversion results before and after filtering were essentially identical, indicating that a point source mechanism is sufficient in modelling the long-period, teleseismic body waves. In the case of the San Fernando earthquake, the point source Green's functions were incapable of simultaneously modelling the P - and SH -waves. Inversion of P -waves alone resulted in extreme parameter resolution problems, but allowed constraint in one axis of the moment tensor and suggested an overall source time function. Inversion of a synthetic San Fernando data set yielded similar results, but allowed an investigation of the shortcomings of the method under controlled circumstances. Although the results may require substantial interpretation, the method presented represents a simple first step in the analysis of complex earthquakes.  相似文献   

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

4.
Summary. The space-time elastic wave motion generated by an impulsive line source in a homogeneous anisotropic medium is calculated with the aid of the Cagniard-de Hoop method. Two types of sources are considered in detail, viz. a line source of expansion (model for an explosive source) and a line force (model for a mechanical vibrator). Numerical results are presented for the radiated particle velocity in the medium including those regions of space where shear-wave triplication occurs. There is a marked difference in the time response observed for the two types of sources and for the different positions of the receiver with respect to the source position. These waveform differences are important when the radiated wave is used to determine experimentally the elastic properties of the medium. As compared with the traditional Fourier-integral transform method to handle this problem, the computation time with the present method is considerably less.  相似文献   

5.
Summary. An algorithm which is part analytical and part numerical is suggested for the computation of complete synthetic seismograms for complex three-dimensional geological structures with radial symmetry. A partial separation of variables based on the combination of a finite Fourier integral transform with respect to the spatial coordinate z together with the finite difference method is the essence of the algorithm. Upon application of the finite transform the problem reduces to solving a system of equations containing only partial derivatives with respect to one spatial coordinate ( r ) and time. As radial symmetry is assumed, there is no functional dependence on φ in the cylindrical system of coordinates ( r , φ, z ). The coefficients of the transformed equations may contain finite Fourier integrals of the z dependence of the elastic parameters. Several examples of synthetic seismograms computed for both SH - and P – SV -waves propagating in complex subsurface geometries are presented and their interpretation discussed.  相似文献   

6.
b
Long-period data of the Global Digital Seismograph Network (GDSN) recorded over the three-year period from 1984 to 1986 were studied for the occurrence of S-P and P-S conversions from the upper mantle transition zone that appear as precursors to teleseismic S arrivals. Conversions of this type were identified on a large number of single-station records. Simple stacking of many records enhanced the appearance of converted phases and demonstrated that no major lateral variations in the nature of the transition zone exist between various tectonic regions. S-P and P-S conversions from the 400 km discontinuity were best observed at distances between 70 and 85 while conversions from the 670 km discontinuity showed up best at distances beyond 87. The analysis of published source mechanisms and comparison with synthetic seismograms suggests that the appearance of converted phases is primarily governed by the earthquake radiation pattern. Phases that have undergone S-P conversions beneath the receiver are best observed from dip-slip events that radiate strong SV - and weak P -waves towards the station. P-S conversions beneath the source area, on the other hand, are frequently observed from events that radiate strong P and little SV energy towards the station, and also from some strike-slip events. Comparison of observed with synthetic seismograms suggests that the PREM model of Dziewonski & Anderson (1981) explains most of the observations. Observed S-P and P-S conversions from the 670 km discontinuity, however, often have larger amplitudes than in the synthetics. Constructive interference of converted waves with the P -wave coda, source radiation effects and a velocity contrast across the 670 km discontinuity which is higher than in PREM may all contribute to the discrepancy.  相似文献   

7.
Summary. Two apparently equivalent methods used in the literature for calculating the mean energy flux and other properties of harmonic plane waves in a viscoelastic medium are shown to give slightly different results when applied to inhomogeneous P/SV -waves. The discrepant results of the two methods have appeared in recent papers and textbooks. Hence, it is important to point out and explain them. It is also shown that the discrepancies occur even for perfect elasticity. The first method makes use of the familiar classical formulae for the energy flux and other quantities, which can be derived from first principles. However, in the second method, the formulae for these quantities are inferred from a specifically derived conservation relation_ The validity of the results of the first method cannot be questioned. It is shown that the errors in the results of the second method are due to the presence of a certain mathematical non-uniqueness.  相似文献   

8.
Summary. An asymptotic procedure for the computation of wave fields in two-dimensional laterally inhomogeneous media is proposed. It is based on the simulation of the wave field by a system of Gaussian beams. Each beam is continued independently through an arbitrary inhomogeneous structure. The complete wave field at a receiver is then obtained as an integral superposition of all Gaussian beams arriving in some neighbourhood of the receiver. The corresponding integral formula is valid even in various singular regions where the ray method fails (the vicinity of caustic, critical point, etc.). Numerical examples are given.  相似文献   

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

10.
Summary. Seismic investigations using shear-wave and converted wave techniques show that very often reflected PS - and SS -waves have anomalous polarizations ( accessory components ). This phenomenon cannot be explained in terms of isotropic models with dipping boundaries. Computations of synthetic seismograms of reflected PS - and SS -waves were made for different models of transversely isotropic media with dipping anisotropic symmetry axes not normal to the boundaries. Synthetic seismograms were computed by ray techniques using an optimization algorithm to construct all rays arriving at a given receiver. These computations indicate that accessory components arise when the medium above the boundary is anisotropic, where they are caused by the constructive interference of qSV - and qSH -waves. If a low-velocity layer is present, displacement vectors of both waves have horizontal projections which are approximately orthogonal. The algorithm for wave separation is presented and some results of its use are given.  相似文献   

11.
Summary. A set of recurrence relations similar to that of Kennett suitable for SH -wave generation in an ( n + l)-layered half-space is presented. The recurrence relations contain no growing terms, thus providing a stable and efficient algorithm for computing complete SH synthetic seismograms. The complete expansion of these recurrence relations gives the explicit form of the transfer function for SH -waves. The transfer function for a point source in layer s of the stratification is a series of 2 n terms in the denominator and a series of 2 n−s +1 terms in the numerator. The result of Wang from ray summation is shown to be a special case of our general result. Numerical comparison of the algorithm of this paper with the propagator matrix method is also made.  相似文献   

12.
A curious observation has been made on radial receiver functions calculated from teleseisms recorded by 29 broad-band seismometers distributed over Iceland. The arrival time of the direct P phase of the radial receiver functions depends critically upon the azimuth of the teleseismic source. For a seismic station in West Iceland, the direct P  phase of the radial receiver function arrives consistently later for easterly source azimuths than for westerly source azimuths. The reverse applies for stations in East Iceland. In the original seismograms, the delayed P phase of the receiver function appears up to 450 ms later on the radial than on the vertical component. The seismometer locations in East and West Iceland are separated by the Neovolcanic Zone, a constructive plate boundary. The delayed P phases occur for seismic rays travelling across this zone. However, it is not obvious how wave propagation across the plate boundary zone could cause the observed delays. The tentative explanation proposed here involves the regional dip of the Icelandic lava sequences towards the Neovolcanic Zone. A dipping interface at shallow depth results in a P–S converted phase arriving shortly after the P phase. These phases cannot be separated in the radial receiver functions, given the bandwidth of the observed signals. However, a calculation of receiver functions from estimates of the P , SV and SH wavefields clearly reveals a P–S converted phase at about 500 ms for easterly source azimuths in West Iceland and for westerly source azimuths in East Iceland. The amplitudes of the direct P phase and the P–S phase converted at a dipping interface would be expected to vary strongly with azimuth in accordance with the observed behaviour.  相似文献   

13.
P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.  相似文献   

14.
Summary. The paper describes some results of experimental seismic prospecting investigations of the anisotropic properties of sedimentary rocks at depths of less than 2.5 km. Shear and converted PS -waves were recorded. Examples of SV and SH velocity distributions and shear-wave polarization are given. The main conclusion is that the medium essentially differs from that usually adopted for sedimentary rocks in that the symmetry axis is not normal to the layering.  相似文献   

15.
We use teleseismic three-component digital data from the Trabzon, Turkey broadband seismic station TBZ to model the crustal structure by the receiver function method. The station is located at a structural transition from continental northeastern Anatolia to the oceanic Black Sea basin. Rocks in the region are of volcanic origin covered by young sediments. By forward modelling the radial receiver functions, we construct 1-D crustal shear velocity models that include a lower crustal low-velocity zone, indicating a partial melt mechanism which may be the source of surfacing magmatic rocks and regional volcanism. Within the top 5 km, velocities increase sharply from about 1.5 to 3.5 km s−1. Such near-surface low velocities are caused by sedimentation, extending from the Black Sea basin. Velocities at around 20 km depth have mantle-like values (about 4.25 km s−1 ), which easily correlate to magmatic rocks cropping out on the surface. At 25 km depth there is a thin low-velocity layer of about 4.0 km s−1. The average Moho velocity is about 4.6 km s−1, and its depth changes from 32 to 40 km. Arrivals on the tangential components indicate that the Moho discontinuity dips approximately southwards, in agreement with the crustal thickening to the south. We searched for the solution of receiver functions around the regional surface wave group velocity inversion results, which helped alleviate the multiple solution problem frequently encountered in receiver function modelling.
Station TBZ is a recently deployed broadband seismic station, and the aim of this study is to report on the analysis of new receiver function data. The analysis of new data in such a structurally complex region provides constraining starting models for future structural studies in the region.  相似文献   

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

17.
Wavepath traveltime tomography   总被引：1，自引：1，他引：1  
The elastic-wave equation is used to construct sensitivity kernels relating perturbations in elastic parameters to traveltime deviations. Computation of the functions requires a correlation of the forward-propagating seismic wavefield with a backward propagation of the residual wavefield. The computation of the wavefields is accomplished using a finite difference algorithm and is efficiently executed on a CM-2 parallel processor. The source and receiver locations have maximum sensitivity to velocity structure. The sensitivity kernels or wavepaths are well suited for transmission traveltime inversion such as cross-borehole tomography and vertical seismic profiling. Conventional ray tomography and wavepath tomography are applied to a set of P -wave arrival times, from a cross-borehole experiment at Kesterson, California. Because the wavepaths have increased sensitivity near the source and receiver there are differences in resolution of the velocity structure. Both techniques recover the same relative variations in velocity where the coverage is adequate. The wavepath solution is more laterally continuous and the dominant variation is vertical, as is expected for the layered sediments in this region.  相似文献   

18.
Summary. The presence of anisotropy has been established and its parameters have been estimated by the combined use of compressional and shear, SH -and SV -waves. Anisotropy has been studied in argillaceous and carbonate rocks to a depth of 2 km by combined surface and borehole observations along longitudinal and non-longitudinal vertical lines with linear shot lines. Anisotropy in carbonate rocks has been found to be a function of their heterogeneity, ordered fracturing and the amount of clay. The highest anisotropy has been found in argillaceous rocks irrespective of their depth and degree of heterogeneity.  相似文献   

19.
A numerical method is presented for calculating complete theoretical seismograms, under the assumption that the earth models have velocity, density and attenuation profiles which are arbitrary piece-wise continuous functions of depth only. Solutions for the stress-displacement vectors in the medium are expanded in terms of orthogonal cylindrical functions. Our method for solving the resulting two-point boundary value problems differs from that of other investigators in three ways. First, collocation is used in traditionally troublesome situations, e.g. for highly evanescent waves, at turning points, and in regions having large gradient in material properties. Second, in some situations (high frequencies and small gradients) P and S -waves decouple and we use a different solution method for each wave type, instead of trying to force a single method to find all solutions. For example, above the P - and S -waves turning points an approximate fundamental matrix may be used for each wave type. At the P -wave turning point, the fundamental matrix may be used for the S -wave components but collocation is used for the P -wave. Between the P - and S -wave turning points collocation is used for the evanescent P -wave and the fundamental matrix is used for the S -wave. At the S -wave turning point and below, collocation is used for both. Third, the computational algorithm chooses the appropriate solution method and depth domain upon which it is employed based upon a specified error tolerance and the known inaccuracies of the various approximations employed. Once solutions of the boundary value problems are obtained, a Fourier—Bessel transform is then applied to get back into the space-time domain.  相似文献   

20.
Use of Kirchhoff s formula for body wave calculations in the Earth   总被引：1，自引：0，他引：1  
Summary. Kirchhoff's time-dependent surface integral representation of a scalar wavefield is applied to the problem of computing synthetic seismograms for P -waves in the Earth. By means of an appropriate parameterization, the Kirchhoff integral is transformed into a convolution of a weight function with the derivative of the source function in the time domain. The weight function is calculated using simple ray theory. The method extends the applicability of simple ray theory to caustics and other diffraction phenomena and allows certain kinds of departures from spherical symmetry to be taken into account. The method is illustrated in detail by application to the PKP -wavefield in the Earth.  相似文献   

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

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