首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 79 毫秒
1.
The blockage of the L g wave by crustal barriers such as continental margins and graben structures has long been recognized as providing a very useful tool for mapping large-scale lateral crustal variations along the propagation path. Numerical simulation of L g -wave propagation in complex anelastic media using the pseudospectral method provides insight into the nature of the propagation process using both snapshots of the wavefield and synthetic seismograms. A variety of 2-D structures have been investigated, including the influence of sediments, crustal thickness and attenuation.
Thick sedimentary basins covering a graben structure can have a major influence, since they remove L g energy by generating P conversion and scattering–the principal mechanisms for strong L g attenuation across a graben. The reduction of the L g energy is reinforced by anelastic attenuation in the sediments as well as the influence of the gradually thinning crustal waveguide associated with an elevated Moho.
The extinction of L g in a sequence of explosions fired across the central graben of the North Sea can be simulated by numerical calculations for the structure derived from refraction experiments.  相似文献   

2.
  相似文献   

3.
The traveltime perturbation equations for the quasi-compressional and the two quasi-shear waves propagating in a factorized anisotropic inhomogeneous (FAI) media are derived. The concept of FAI media simplifies considerably these equations. In the FAI medium, the density normalized elastic parameters a ijkl ( X i ) can be described by the relation a ijkl ( X i) = f 2( x i ) A ijkl, where A ijkl are constants, independent of coordinates x i and f 2( x i) is a continuous smooth function of x i . The types of anisotropy ( A ijkl ) and inhomogeneity [ f ( x i)] are not restricted. The traveltime perturbations of individual seismic body waves ( q P , qS 1 and qS 2) propagating in the FAI medium depend, of course, both on the structural pertubations [δ f 2( x i)] and on the anisotropy perturbations (δ A ijkl ), but both these effects are fully separated. The perturbation equations for the time delay between the two qS -waves propagating in the FAI medium are simplified even more. If the unperturbed (background) medium is isotropic, the perturbation of the time delay does not depend on the structural perturbations (δ f 2( x i) at all. This striking result, valid of course only in the framework of first-order perturbation theory, will simplify considerably the interpretation of the time delay between the two split qS -waves in inhomogeneous anisotropic media. Numerical examples are presented.  相似文献   

4.
Summary. The paper gives the results of a study of the anisotropy of seismic wave velocities within the Ashkhabad test field in Central Asia. The anisotropy was studied by analysing variations in the values of apparent velocities of first arrivals for epicentral distances ranging from 30 to 130 km and by analysing the delays (Δ ts1-s2 ) between the arrival times of shear waves with different polarizations.
The velocities of P -waves vary with azimuth from 5.3 to 6.27 km s-1 and the velocities of S -waves vary from 3.15 to 3.5 km s-1.
The delay times Δ tS1 - S2 depend on the direction of the propagation. The character of the variation of the propagation velocity of the longitudinal wave, the presence of two differently polarized shear waves S 1 and S 2 propagating at different velocities, and the character of the distribution of Δ tS1 - S2 on the stereogram suggest that the symmetry of the anisotropic medium is close to hexagonal with a nearly horizontal symmetry axis coinciding with the direction of maximal velocity. The azimuth of the symmetry axis of the medium is 140° and coincides with the direction of geological faults.  相似文献   

5.
Some comments on the descriptions of the polarization states of waves   总被引：4，自引：0，他引：4  
Summary. Procedures are developed for specifying the polarization characteristics of n -dimensional waves, and in particular three-dimensional waves of geophysical interest. We show that when a wave is in a pure state or is totally polarized, all the polarization information can be represented by a single vector u in an n -dimensional unitary space. Simple measures of the degree of polarization of the wave are constructed from the characteristic equation of the spectral matrix S . These measures are functions only of the scalar invariants of S and consequently S need not be diagonalized. If S represents a purely polarized wave, the unitary vector which contains the polarization information about the wave can be obtained directly from S using any 2 n – 1 equations of n 2 possible equations. By multiplying by a phase-factor this unitary vector can be written in the form u = r 1+ i r 2 where r 1 and r 2 are orthogonal vectors in a real space. For an elliptically polarized wave, r 1 and r 2 locate the major and minor axes of the ellipse, and the ellipticity is given by the ratio of their magnitudes. The polarization parameters of ULF magnetic waves at the Earth's surface are computed from one set of five equations ( n = 3) and compared with parameters calculated using established techniques.  相似文献   

6.
The attenuation mechanism of seismic waves in northwestern Himalayas   总被引：3，自引：0，他引：3  
We analysed local earthquake waveforms recorded on a broad-band seismic network in northwestern Himalayas to compute the intrinsic and scattered attenuation parameters from coda waves. Similar to other tectonically active and heterogeneous regions, attenuation-frequency relation for western Himalaya is   Q −1 c = (113 ± 7)  f (1.01±0.05)  where   Qc   is the coda Q parameter. Intrinsic  ( Q −1 i )  and scattering  ( Q −1 s )  attenuations was separated using   Qc   and direct S -wave Q data  ( Qd )  . It is observed that estimated   Q −1 c   is close to   Q −1 i   and both of them are much larger than   Q −1 s   suggesting that coda decay is predominantly caused by intrinsic attenuation. At higher frequencies, both the attenuation parameters   Qc   and,   Qd   are similar indicating that coda is predominantly composed of back-scattered S waves at these frequencies.  相似文献   

7.
3-D images of P velocity and P - to S -velocity ratio have been produced for the upper crust of the Friuli area (northeastern Italy) using local earthquake tomography. The data consist of 2565 P and 930 S arrival times of high quality. The best-fitting V P and V P / V S 1-D models were computed before the 3-D inversion. V P was measured on two rock samples representative of the investigated upper layers of the Friuli crust. The tomographic V P model was used for modelling the gravity anomalies, by converting the velocity values into densities along three vertical cross-sections. The computed gravity anomalies were optimized with respect to the observed gravity anomalies. The crust investigated is characterized by sharp lateral and deep V P and V P / V S anomalies that are associated with the complex geological structure. High V P / V S values are associated with highly fractured zones related to the main faulting pattern. The relocated seismicity is generally associated with sharp variations in the V P / V S anomalies. The V P images show a high-velocity body below 6 km depth in the central part of the Friuli area, marked also by strong V P / V S heterogeneities, and this is interpreted as a tectonic wedge. Comparison with the distribution of earthquakes supports the hypothesis that the tectonic wedge controls most of the seismicity and can be considered to be the main seismogenic zone in the Friuli area.  相似文献   

8.
We report on calculations of the on-shore run-up of waves that might be generated by the impact of subkilometre asteroids into the deep ocean. The calculations were done with the COULWAVE code, which models the propagation and shore-interaction of non-linear moderate- to long-wavelength waves  ( kh < π)  using the extended Boussinesq approximation. We carried out run-up calculations for several different situations: (1) laboratory-scale monochromatic wave trains onto simple slopes; (2) 10–100 m monochromatic wave trains onto simple slopes; (3) 10–100 m monochromatic wave trains onto a compound slope representing a typical bathymetric profile of the Pacific coast of North America; (4) time-variable scaled trains generated by the collapse of an impact cavity in deep water onto simple slopes and (5) full-amplitude trains onto the Pacific coast profile. For the last case, we also investigated the effects of bottom friction on the run-up. For all cases, we compare our results with the so-called 'Irribaren scaling': The relative run-up   R / H 0=ξ= s ( H 0/ L 0)−1/2  , where the run-up is   R , H 0  is the deep-water waveheight, L 0 is the deep-water wavelength, s is the slope and ξ is a dimensionless quantity known as the Irribaren number. Our results suggest that Irribaren scaling breaks down for shallow slopes   s ≤ 0.01  when  ξ < 0.1 − 0.2  , below which   R / H 0  is approximately constant. This regime corresponds to steep waves and very shallow slopes, which are the most relevant for impact tsunami, but also the most difficult to access experimentally.  相似文献   

9.
We investigate the particle orbits of long-period (about 20 s) P waves observed with the global seismic network. By analysing 84 three-component seismograms recorded at 25 stations from 60 earthquakes occurring beneath 300 km, we quantitatively evaluate the orbits by three sets of eigenvalues and eigenvectors, using a covariance matrix method. The eigenvalues for P waves recorded at stations located on continents are explained by the standard horizontal layered structure model (iasp91). On the other hand, the orbits observed at stations close to island arcs are affected not only by the horizontal layered structure but also by heterogeneity due to subducting plates, mantle diapirs and so on. On the basis of a single-scattering model for a plane P wave, we quantify the heterogeneities by an isotropic scattering coefficient g0. Fitting the theoretical eigenvalues to the observed ones, we estimate g0 for the crust and upper mantle beneath continents to be less than 0.0005 km-1, and the mean g0 for the structure beneath island arcs to be about 0.0015 to 0.003 km-1.  相似文献   

10.
Seismic sources with observable glut moments of spatial degree two   总被引：1，自引：0，他引：1  
Let ζΛ and r Λ. be the hypocentral position and time of an extended indigenous seismic source. Backus showed that the force moment tensors of the source, Γ( m +1, n )(ζΛ, r Λ), determine and are determined by the motion which the source produces. For small m + n , only the long-period motion is relevant. The glut moment tensor Λ( m,n ) (ζΛ, r Λ.) can be calculated uniquely from γ( m +1, n )(ζΛ r Λ) only if m = 0 or m = 1. The tensor G =Λ(2,0)(ζΛ) gives the spatial variance tensor WΛ of the source, and WΛ. roughly describes the size, shape and orientation of the source region. Therefore the failure of the observed F =Γ(3,0)(ζΛ) to determine G uniquely is of seismological interest. In the present paper we show that F determines G uniquely if we assume the source to be a simple straight line source (SSLS) or an ideal fault in an isotropic medium with isotropic prestress (IFIMIP). We give tests on F which determine whether it can come from a SSLS, from an IFIMIP or from a simple plane surface source (SPSS). If we assume the source to be a SPSS then knowing F and the fault plane determines G to within an unknown scalar multiple of a certain tensor tangent to the fault plane. Moreover F determines the fault plane uniquely unless F can come from a SSLS. If it can, then F determines this virtual source line uniquely, and F permits the fault plane to be any plane containing the virtual source line.  相似文献   

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

12.
Summary. The limitations of asymptotic wave theory and its geometrical manifestations are newly formalized and scrutinized. Necessary and sufficient conditions for the existence of acoustic and seismic rays and beams in general inhomogeneous media are expressed in terms of new physical parameters: the threshold frequency ω0 associated with the P/S decoupling condition, the cut-off frequency ωc associated with the radiation-zone condition, the total curvature of the wavefront and the Fresnel-zone radius.
The analysis is facilitated with the introduction of a new ancillary functional – the hypereikonal which is capable of representing ordinary as well as evanescent waves. The hypereikonal is the natural extension of the eikonal theory.
With the aid of the above new parameters, simple conditions are obtained for the decoupled far field, the decoupled near field, two point dynamic ray tracing, paraxial wavefields and Gaussian beams.  相似文献   

13.
Summary. The geopotential is usually expressed as an infinite series of spherical harmonics, and the odd zonal harmonics are the terms independent of longitude and antisymmetric about the equator: they define the 'pear-shape' effect. The coefficients J 3, J 5, J 7, … of these harmonics have been evaluated by analysing the variations in eccentricity of 28 satellite orbits from near-equatorial to polar. Most of the orbits from our previous determination in 1974 are used again, but three new orbits are added, including two at inclinations between 62° and 63°, which have been specially observed for more than five years by the Hewitt cameras. With the help of the new orbits and revised theory, we have obtained sets of J -coefficients with standard deviations about 40 per cent lower than before. A 9-coefficient set is chosen as representative, and is as follows (all × 109): J 3=– 2530 ± 4, J 5=–245 ± 5, J 7=–336 ± 6, J 9=–90 ± 7, J 11= 159 ± 9, J 13=–158 ± 15, J 15=– 20 ± 15, J 17=– 236 ± 14, J 19=– 27 ± 19. With this set of values, the pear-shape asymmetry of the geoid (north polar minus south polar radius) amounts to 45.1 m instead of the previous 44.7 m. The accuracy of the longitude-averaged geoid profile is estimated as 50 cm, except at latitudes above 86°. The geoid profile and predicted amplitude of the oscillation in eccentricity are compared with those from other sources.  相似文献   

14.
Observation of Coriolis coupled modes below 1 mHz   总被引：3，自引：0，他引：3  
We present observations of spectral energy at toroidal mode frequencies in vertical seismic recordings of the 1998 Balleny Islands earthquake. Since toroidal modes on a spherically symmetric, nonrotating Earth have horizontally polarized particle motion these observations call for an explanation. We first rule out local and instrumental effects as being responsible for the verticalcomponent signal of the toroidal modes 0 T 3 (0.59 mHz) and 0 T 4 (0.77 mHz). The global effects that we consider are general heterogeneous mantle structure, ellipticity of figure and rotation. We find that rotation through Coriolis coupling of loworder spheroidal and toroidal oscillations is the dominant mechanism.  相似文献   

15.
b
The spatial distribution of earthquakes is a fractal, which is characterized by a fractal dimension. However, if a spatial distribution has a heterogeneous fractal structure, a single value of fractal dimension [e.g. Do (capacity dimension) or D 2 (correlation dimension)] is not enough to characterize it. From a multifractal viewpoint, we analysed the spatial distribution of microearthquakes in the Kanto region by using a local density function. Generalized dimensions, Dq , of the spatial distribution were calculated from the slopes of generalized correlation integrals, Cq(r) versus distance r , on a log-log plot, examining the self-similarity of the spatial distribution of microearthquakes. Self-similar structures are held well at scales from 1.26 to 12.6 km. Our results suggest that the spatial distribution of microearthquakes in the Kanto region is not a homogeneous fractal structure but a heterogeneous one with generalized dimensions D2 = 2.2 ≤ D 3≤…≤ D ∞= 1.7. The value of D ∞, the lower limit of fractal dimension, is the fractal dimension of the most intensive clustering in the heterogeneous fractal set. The fractal dimension of the most intensive clustering of microearthquakes in the Kanto region is 1.7.  相似文献   

16.
Summary A comparison is made of synthetic and natural magnetites after oxidation and reduction. During these reactions new spinel phases are formed from basaltic magnetites which fall outside the composition diagram FeO-Fe2O3-TiO2. They are solid solutions of MgFe2O4, MnFe2O4, FeAl2O4 and FeFe2O4 depending on the minor components of natural magnetites like MgO, MnO, Al2O3. A irascibility gap makes the situation more complicated.  相似文献   

17.
Summary. Records of 21 earthquakes with ray paths crossing south-eastern China have been examined to establish the nature of propagation of the seismic wave Sn . Sn is seen at stations within 1500km of earthquakes, but not beyond 1500km. No single localized region of attenuation can account for all the missing or weak Sn waves. The inefficient propagation of Sn through south-eastern China is surprising in light of efficient Sn propagation in other continental platforms which, like south-eastern China, appear to have been stable throughout the Cenozoic.  相似文献   

18.
Summary. Using nine IDA records for the Indonesian earthquake of 1977 August 19, we have formed an optimal linear combination of the records and have measured the frequency and Q of 0 S 0 and 1 S 0. The frequency was measured using the moment ratio method. The attenuation was measured by the minimum width method and by the time-lapse method. The frequency and attenuation were measured simultaneously by varying them to obtain a best fit to the data. A 2000-hr stack, the sum of nine individual records, for 0 S 0 gave a frequency of 0.814664 mHz±4 ppm. The values for the Q of 0 S 0 for the three different methods of measurement were 5600,5833 and 5700, respectively. The error in the estimates of Q -1 is about 5 per cent for the minimum power method. For 1 S 0 a 300-hr stack yielded a frequency of 1.63151 mHz±30 ppm. The values of Q for this mode were 1960, 1800 and 1850, respectively, with an error in Q -1 of about 12 per cent for the minimum power method.  相似文献   

19.
Summary. The potential function ø for a magnetic body of susceptibility μ in a medium of susceptibility μ* satisfies the integral equation
Here Φ* is the potential function for the region without the heterogeneity and R is the distance from the point of observation to the point on the surface, s , of the body. δΦ /δn is the normal derivative, in the direction of the outward normal. The equation allows for the effects of demagnetization. For numerical purposes the surfaces can be divided into N facets over which δΦ/δ n is a constant. The unknown quantities δΦ/δnj can be found from the system of equations defined by:
The prime on the summation sign denotes that the summation does not include the i th element. The magnetic field in the direction of the unit vector P( P 1, P 2, P3 ) is given by:   相似文献   

20.
Summary. In cases where directional data, such as palaeomagnetic directions, lie nearly along a great circle, a good approximation to the maximum likelihood estimate of the intermediate concentration parameter k 2 in the Bingham probability distribution is given by: 2( t 2/ N ) – 1 = I 1(1/2 k 2)/ I 0(1/2 k 2), where t 2 is the intermediate eigenvalue, N is the number of samples, and the Ii are the appropriate modified Bessel functions of the first kind. This estimate, the asymptotic limit as the smallest eigenvalue t 1→ 0, corresponds to restricting all points to lie on a great circle. The limit is also useful as an endpoint for interpolation, especially since numerical calculation in this region is difficult.  相似文献   

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

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