首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The three dimensional scattering of near-field, from a point source, is studied for acceleration in the time domain. The perturbation method is applied to define the acceleration for the first order scattering from a weak inhomogeneity in a homogeneous surrounding. A body force, arising from the interaction between the primary waves and the inhomogeneity, acts as the source generating the scattered motion. The acceleration of scattered waves is related to the velocity and density fluctuations of the inhomogeneity. No restrictions are placed on the inhomogeneity size or locations of the source and receiver. Decoupling of scattered motion enables the identification of different phases. Integral expressions are derived for the scattering acceleration due to the incidence of near-field wave (from an impulsive point force) at a radially inhomogeneous volume element. These integrals are solved further for scattering from an inhomogeneous spherical shell. The accelerations for back scattering are obtained as a special case. These accelerations are simple analytically solvable expressions in closed form. Only spherical asymmetry ofP wave velocity inhomogeneity can affect the scatteredS acceleration. ScatteredP acceleration is affected by the gradient ofS wave velocity inhomogeneity. The back scattering of near-field from a spherical shell, is independent of radial inhomogeneity ofP wave velocity. Inhomogeneity with smoothly perturbedS wave velocity does not back-scatter any acceleration. Accelerations are computed numerically for scattering from a part of inhomogeneous spherical shell. Hypothetical models are considered to study the effects of the distances of spherical shell from source, receiver, its thickness and its position relative to the direction of impulsive force.  相似文献   

2.
Three dimensional scattering of near-field is studied for dilatation and rotation in the time domain. The perturbation method is applied to solve the equation of motion for the first order scattering from a weak inhomogeneity in an otherwise homogeneous medium. The inhomogeneity is assumed close enough to the point source so that the near-field intermediate wave is dominating over the far-field sphericalP andS pulses. The integral expressions are derived to relate dilatation and rotation of scattering to the radial fluctuations of velocities and density in the inhomogeneity. These integrals are solved to calculate the strains of scattering from (a part of) an inhomogeneous spherical shell of arbitrary curvature. Variable curvature may allow the shape of inhomogeneity volume element to change uniformly from spherical to rectangular. Rotation of scattering from a spherical shell is independent ofP wave velocity inhomogeneity. Dilatation of scattering does not involveS wave velocity inhomogeneity but its gradient. The back scattering results are obtained as a special case. Strains are computed numerically, for hypothetical models to study the effects of various parameters viz., velocity inhomogeneity, distance of source from inhomogeneity and from receiver, and thickness of inhomogeneity. The curvature of the spherical shell is varied to study the effects of the shape of inhomogeneous volume element on scattering.  相似文献   

3.
用数值模拟方法研究了地下异质体对波散射,分析了散射波对表面波相速度影响,结果表明:由相速度曲线扰动难以有效确定异质体方位。在表面波场中,散射波能量相对较小,不易辨别,以振源为中心,在测线上对称布置一组测点,通过两侧测点质点响应曲线比较,得到表面散射波场。在表面散射波场中,散射瑞利波是主要的,由散射瑞利波走时随测点位置变化可判断地下异质体的方位。最后,给出地下异质体方位测试方案。  相似文献   

4.
As Morozov [Morozov, I. B. (2004). Crustal scattering and some artefacts in receiver function images. Bull. Seismol. Soc. Am., 94 (4), 1492–1499.] suggested, for a teleseismic array targeting subducting crust in a zone of active subduction, scattering from the strong horizontal velocity heterogeneity beneath the trench zone itself produces subhorizontally-propagating waves that should be observed as coherent dipping events in receiver functions (RF). Due to similar RF delay times and moveouts, these events could be difficult to distinguish from backscattered P- and S-wave modes. To further verify this suggestion, we performed a full-waveform, 3-D visco-elastic finite-difference modelling of teleseismic wave propagation within a simplified model of a subduction zone. The synthetics show strong scattering from the area beneath the trench, dominated by the mantle and crustal P-waves propagating at 6.2–8.1 km/s and slower. These scattered waves occupy the same time and moveout intervals as the backscattered converted modes, and also have similar amplitudes. Although their amplitude decay characters are different, the uncertainty in the knowledge of the velocity and density structure of the subduction zone could make distinguishing between these modes difficult. However, under minimal assumptions, recent observations of receiver function amplitudes decreasing away from the trench support the interpretation of (sub-) trench-zone scattering.Although still limited in its representation of crustal heterogeneity, 3-D modelling suggests that scattering from near-Moho crustal structures plays a key role in the formation of teleseismic wavefields. Recognition of scattered noise in teleseismic records could help to constrain major crustal structures, particularly those with strong horizontal velocity contrasts at near-Moho depths, such as crustal sutures, subduction fault zones, and mountain roots. Matching of the observed arrivals with wavefield synthetics could help constrain the locations and parameters of such structures and also help substantiate the interpretations.  相似文献   

5.
采用波函数展开法结合Graf加法定理,利用桩体与土体之间的位移和应力连续边界条件,给出了全空间单排桩对柱面SH波散射解析解,并采用傅里叶逆变换,求得时域结果。该方法考虑了入射波曲率的影响,在频域内分析了柱面SH波入射时排桩散射的频谱规律,给出了时域排桩柱面SH波散射位移云图,讨论了桩数与桩间距对柱面SH波入射时排桩散射的影响。研究表明:在振源距排桩较近(d/a=10)时,低频段(η=0~1.0)排桩对柱面SH波的隔离作用显著;振源距排桩较远时,各频段排桩对柱面SH波的隔离作用均较为显著。桩间距固定时,增加桩数,排桩后方首次出现的最大位移响应相应减小,排桩散射的影响范围随之增大;对柱面SH波进行阻隔时,为提高阻隔效率的同时节约成本,不仅需要考虑桩数,还应考虑波源与桩体、桩间空隙相对位置的影响;排桩分布宽度固定时,由于入射波曲率对排桩散射的影响,减小桩间距,排桩后方首次出现的最大位移响应有可能会出现放大的情况,故应采用合理的桩间距对柱面SH波进行阻隔。  相似文献   

6.
The heterogeneous upper mantle low velocity zone   总被引:2,自引:1,他引:2  
The upper mantle low velocity zone (LVZ) is a depth interval with slightly reduced seismic velocity compared to the surrounding depth intervals. The zone is present below a relatively constant depth of 100 km in most continental parts of the world, both in cratonic areas with high average velocity and tectonically active areas with low average velocity. Evidence for the low velocity zone arises from controlled and natural source seismology, including studies of surface waves and of primary and multiple reflections of body waves from the bounding interfaces, calculations of receiver functions, and absolute velocity tomography. The available data indicates a more pronounced reduction in seismic velocity and Q-value for S-waves than P-waves as well as high electrical conductivity in the LVZ. Seismic waves are strongly scattered by the zone, which demonstrates the existence of small-scale heterogeneity. The depth to the base of the LVZ is systematically shallower in cold, stable cratonic areas than in hot, active regions of the world. Because of its global occurrence below a relative constant depth of 100 km, the LVZ cannot be explained by metamorphic or compositional variation and rheological changes. Calculated upper mantle temperatures indicate that the rocks are close to the solidus in an interval with variable thickness below 100 km depth, provided that the rocks contain water and carbon dioxide. The presence of, even small amounts of such fluids in the mantle rocks will lower the solidus by several hundred degrees and introduce a characteristic kink on the solidus curve around 80–100 km depth. The seismic velocities and Q-values are significantly reduced of rocks, which are close to the solidus or contain small amounts of partial melt. Hence, the LVZ may be explained by upper mantle temperatures being close to the solidus in a depth interval below 100 km. Assuming that the rocks contain only limited amounts of fluids, this mechanism may explain the low velocities, Q-values, and resistivity, as well as the intrinsic scattering, and the characteristic variation in thickness of the low velocity zone.  相似文献   

7.
Near-surface diffractors are one of the problems in land seismic exploration. They can scatter the surface wave energy emanating from the seismic source and contaminate the signal received by seismic receivers. The scattered energy from the near-surface diffractors manifests itself on seismic shot gathers as strong hyperbolic events, called diffractions, masking the weakly reflected body waves. Diffractions present complications to most of surface-wave suppression schemes, especially when they have been scattered by scatterers outside the receiver line. Different methods have been used to eliminate diffractions from seismic data, including geophone arrays, filtering, and inverse scattering. Each of those methods has its own limitations. In this study, we present processing algorithms to map and attenuate near-surface diffractors of surface waves in seismic shot gathers. The mapping algorithm is based on semblance measurements and time–offset relations, while the attenuation algorithm is based on the least-square fitting of a source wavelet. The algorithms are applied on synthetic data from two different models. The first model has three near-surface diffractors, while the second model has three clusters of near-surface diffractors. Each cluster consists of three near-surface diffractors with a different geometry for each cluster. The results show that the proposed algorithms are successful in locating and attenuating most near-surface diffractors, except when the separation between individual diffractors is below the wavelength of the diffracted surface wave.  相似文献   

8.
《Comptes Rendus Geoscience》2014,346(11-12):307-316
We analyze the sensitivity of cross correlations to the anisotropy of the incident field in the context of seismic ambient noise monitoring of small velocity changes. Numerical simulations of elastic waves are performed in a 2D scattering plate with a focus on the comparative character of the direct and coda waves in the cross-correlation. We show that coda waves reconstructed from cross-correlations are far more robust than direct waves in the presence of azimuthal anisotropy of the incident field. We observe similar behavior with real data recorded on Erebus volcano, where a database of impulsive icequakes is used to simulate an anisotropic source field. We propose a simplified approach to evaluate the sensitivity of scattered waves to the anisotropy of the incoming noise field. We rely on previous results obtained for direct waves and on intrinsic properties of scattered waves to predict the errors produced by strong source anisotropy with numerical experiments. These results also yield realistic values for monitoring the accuracy to be expected with real data at crustal scales. Our analysis shows that high-precision noise-based monitoring could be performed with coda waves in the correlation functions, even in the presence of variations in the azimuthal distribution of the ambient noise field.  相似文献   

9.
The requirements for the performance of a ground-penetrating radar (GPR) system for detecting subsurface cavities are analyzed by numerical modeling of the GPR problem. The algorithm used to solve the forward GPR problem is approximated to a real experiment with regard to the design of the GPR system, the parameters of the source and receiver, and their position relative to the medium under study and its inhomogeneity. We calculated the spatiotemporal distribution of the field of the detected signal from a pulse source located at the interface between the medium and a cavity anomaly of a given geometry. The results were used to estimate the dynamic range of the GPR system necessary for determining the anomaly. We also performed GPR surveys of low-contrast inhomogeneities (cavities in mines) using GROT 12 GPR systems and analyzed the survey results by numerical modeling. It is shown that the GPR performance required to detect and locate inhomogeneities of interest at a certain sounding depth can be estimated in the experiment design phase.  相似文献   

10.
地震照明分析能够定量分析影响逆时偏移成像质量的因素。文章基于Poynting矢量方法进行地震定向照明分析及成像补偿方法研究,利用Poynting矢量进行波场传播角度分解,通过计算不同入射及散射波方向的部分照明建立局部照明矩阵。该矩阵包含了照明分析中波与地质结构相互作用所涉及的所有角度信息,可以方便地分析目标附近所有可能的入射波和散射波对局部照明的贡献。根据局部照明矩阵可以获得对地下结构的照明描述、采集倾角响应,并依此进行成像补偿计算。以SEG/EAGE模型作为数值计算实例,进行基于Poynting矢量的照明分析及成像补偿,有效地提高了成像质量。该方法为角度域逆时偏移成像校正及分辨率分析等提供了高效的处理手段。  相似文献   

11.
We study some effects that produce 3D topographical irregularities under incoming elastic waves from point sources using the indirect boundary element method. This technique is based on the representation of elastic waves in terms of single-layer boundary sources. In this way reflected and diffracted waves are constructed at the boundaries from where they are radiated by means of boundary sources. The field emited from the point source is computed with the moment tensor for a shear dislocation and analytical expressions of the elastodynamic 3D Green's functions. In this way we can construct a double couple with variable orientation. We compare the technique with that of Bouchon, who used the discrete wave number method for a shear dislocation in a halfspace with a triangular source function. We subsequently apply our method to simulate the seismic response of a mountain of ellipsoidal geometry. The source function that we have considered is a triangular pulse and we show results corresponding to the velocity registered over the surface of the irregularity and that of the halfspace. This motion is presented by means of snapshots showing the evolution of the wavefields that are present in the problem.  相似文献   

12.
The propagation of plane waves in an anisotropic elastic medium possessing monoclinic symmetry is discussed. The expressions for the phase velocity ofqP andqSV waves propagating in the plane of elastic symmetry are obtained in terms of the direction cosines of the propagation vector. It is shown that, in general,qP waves are not longitudinal andqSV waves are not transverse. Pure longitudinal and pure transverse waves can propagate only in certain specific directions. Closed-form expressions for the reflection coefficients ofqP andqSV waves incident at the free surface of a homogeneous monoclinic elastic half-space are obtained. These expressions are used for studying numerically the variation of the reflection coefficients with the angle of incidence. The present analysis corrects some fundamental errors appearing in recent papers on the subject.  相似文献   

13.
Voronin  E. A.  Nosov  V. N.  Savin  A. S. 《Doklady Earth Sciences》2020,493(1):569-571
Doklady Earth Sciences - The inverse problem of generating surface waves in a liquid is to determine the parameters of the source of perturbations from the waves it creates on the liquid surface. A...  相似文献   

14.
L. Margerin   《Tectonophysics》2006,416(1-4):229
Most theoretical investigations of seismic wave scattering rely on the assumption that the underlying medium is statistically isotropic. However, deep seismic soundings of the crust as well as geological observations often reveal the existence of elongated or preferentially oriented scattering structures. In this paper, we develop mean field and radiative transfer theories to describe the attenuation and multiple scattering of a scalar wavefield in an anisotropic random medium. The scattering mean free path is found to depend strongly on the propagation direction. We derive a radiative transfer equation for statistically anisotropic random media from the Bethe–Salpeter formalism and propose a Monte Carlo method to solve this equation numerically. At longer times, the energy density is shown to obey a tensorial diffusion equation. The components of the diffusion tensor are obtained in closed form and excellent agreement is found between Monte Carlo simulations and analytical solutions of the diffusion equation. The theory has important potential implications for lithospheric models where scatterers are for example flat structures preferentially aligned along the surface. In this simple geometry, analytical expressions of the Coda Q parameter can be given explicitly in the diffusive regime. Our results suggest that pulse broadening and coda decay are controlled by different parameters, related to the eigenvalues of the diffusion tensor. These eigenvalues can differ by more than one order of magnitude. This theory could be applied to probe the anisotropy of length scales in the lithosphere.  相似文献   

15.
Exact analytical expressions for the scattering phase matrix are presented for various models of a turbulent magnetic field having a symmetry axis in the solar atmosphere. The polarization of the scattered radiation in the plane of the scattering is shown to be virtually independent of the predominant direction of the turbulent magnetic field. The Hanle effect does not operate in a longitudinal magnetic field normal to the surface of the atmosphere, since the phase matrix has the form of a resonance phase matrix with zero magnetic field.  相似文献   

16.
In this paper, a stochastic analysis of surface rotations during earthquakes is presented. The surface motion is modeled by a spacetime random field. A method to obtain rotation based on direct differentiation of the random field is proposed with a reference to a second method based on wave decomposition. The first method leads to a spectral density of surface rotations in terms of point translation spectral density, coherence function and apparent waves velocity, while the second method formulates the rotational spectrum in terms of point translation spectral density and wave parameters. The advantages and disadvantages of both methods are discussed in detail.On leave from: Faculty of Civil Engineering, Technical University, ul. Katowicka 48, 45–951 Opole Poland.  相似文献   

17.
An analytical-numerical examination has been made of the behavior of harmonic Rayleigh-type waves in a compressible elastic half-space of constant density and whose shear modulus increases linearly with depth. The dispersion law is presented for values of Poisson's ratio of 1/2, 1/3 and 1/4 for wide ranges of degree of inhomogeneity.  相似文献   

18.
The dynamic processes within the Earth leave their record in geophysical and geochemical variation about the general stratification with depth. A snapshot of current structure is provided by geophysical evidence, whereas geochemical information provides a perspective over the age of the Earth. The combination of information on the distribution of heterogeneity from geophysical and geochemical sources provides enhanced insight into likely geodynamic processes. A variety of techniques can be used to examine crustal structure, but the major source of information on seismic heterogeneity within the Earth comes from tomographic studies, exploiting surface waves for the lithosphere and body waves for the bulk Earth. A powerful tool for examining the character of mantle heterogeneity is the comparison of images of bulk-sound and shear-wave speed extracted in a single inversion, since this isolates the dependencies on the elastic moduli. Such studies are particularly effective when a common path coverage is achieved for P and S as, for example, when common source and receiver pairs are extracted for arrival times of the phases. The relative behaviour of the bulk-sound and shear-wave speeds allows the definition of heterogeneity regimes. For subduction zones, a large part of the imaged structure comes from S-wave speed variations. The narrow segments of fast wave speeds in the lower mantle, in the depth range 900 – 1500 km, are dominated by S variations, with very little bulk-sound contribution, so images of P-wave speed are controlled by shear-wave variability. Deep in the mantle, there are many features with high seismic-wave speed without an obvious association with subduction in the last 100 million years, which suggests long-lived preservation of components of the geodynamic cycle. The base of the Earth's mantle is a complex zone with widespread indications of heterogeneity on many scales, discontinuities of variable character, and shear-wave anisotropy. Discordance between P- and S-wave speed anomalies suggests the presence of chemical heterogeneity rather than just the effect of temperature.  相似文献   

19.
Imag(in)ing the continental lithosphere   总被引:1,自引:0,他引:1  
This paper is primarily concerned with seismically imaging details in the mantle at an intermediate scale length between the large scales of regional and global tomography and the small scales of reflection profiles and outcrops. This range is roughly 0.1–1 km < a < 10–102 km, where a is the scale. We consider the implications of several models for mantle evolution in a convecting mantle, and possible scales present in the non-convecting tectosphere. Reflection seismic evidence shows that the structures preserved from continental accretion within and at the margins of the Archean cratons are subduction related, and we use subduction as an analog for scales left by past events. In modern orogenic belts we expect to find subduction structures, small scale upper mantle convection structures, and basalt extraction structures. We examine some of the scales that are likely formed by orogenic processes.We also examine the seismic velocity and density contrasts expected between various upper mantle constituents, including fertile upper mantle, depleted upper mantle, normal and eclogitized oceanic crust, and fertile mantle with and without partial melt. This leads directly to predicting the size of seismic signals that can be produced by specular conversion, and scattering from layers and objects with these contrasts.We introduce an imaging scheme that makes use of scattered waves in teleseismic receiver functions to make a depth migrated image of a pseudo-scattering coefficient. Image resolution is theoretically at least an order of magnitude better than traveltime tomography. We apply the imaging scheme to three data sets from 1) the Kaapvaal craton, 2) the Cheyenne Belt, a Paleoproterozoic suture between a protocontinent and an island arc, and 3) the Jemez Lineament, a series of aligned modern volcanic structures at the site of a Proterozoic suture zone. The Kaapvaal image, although not defining a unique base of the tectosphere, shows complicated “layered” events in the region defined as the base of the tectosphere in tomography images. The image of the transition zone discontinuities beneath the Kaapvaal craton is remarkable for clarity. The migrated receiver function image of the upper mantle beneath the Cheyenne belt is complicated, more so than the tomography image, and may indicate limitations in the receiver function imaging system. In contrast the Jemez Lineament image shows large-amplitude negative-polarity layered events beneath the Moho to depths of 120 km, that we interpret as melt-containing sills in the upper mantle. These sills presumably feed the Quaternary–Neogene regional basaltic volcanic field.  相似文献   

20.
The conformal mapping method of complex functions and expansion method of wave functions are adopted, the stresses and displacements at the boundary between the cavities and inner fluids are considered as continuous, and then the theoretical solutions are obtained about the scattering problems of incident P waves by several fluid-filled cavities. Two circular cavities are taken as examples, the influences of frequencies of incident P waves, whether there are inner fluids or not and the spacing distances between the two cavities on the scattering amplitudes are studied, and the results reveal that: the scattering amplitudes of two cavities decrease when the distance between the two cavities increase and they reach constants when the spacing distance exceeds 7.5 times of cavities radius, and the inner fluids would absorb part of energy of the scattered waves and the scattering amplitude of the cavities filled with fluids would slightly smaller than those without fluids.  相似文献   

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

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