共查询到20条相似文献，搜索用时 468 毫秒
1.
Arthur Frankel 《Pure and Applied Geophysics》1989,131(4):639-685
This paper reviews applications of the finite-difference and finite-element methods to the study of seismic wave scattering in both simple and complex velocity models. These numerical simulations have improved our understanding of seismic scattering in portions of the earth where there is significant lateral heterogeneity, such as the crust. The methods propagate complete seismic wavefields through highly complex media and include multiply scattered waves and converted phases (e.g.,P toSV, SV toP, body wave to surface wave). The numerical methods have been especially useful in cases of moderate and strong scattering in complex media where multiple scattering becomes important. Progress has been made with numerical methods in understanding how near-surface, low-velocity basin structures scatter surface waves and vertically-incident body waves. The numerical methods have proven useful in evaluating scattering of surface waves and body waves from topography of both the free surface and interfaces buried at depth. Numerical studies have demonstrated the importance of conversions from body waves to surface waves (andvice versa) when lateral heterogeneities and topographic relief are present in the uppermost crust. Recently, several investigations have applied numerical methods to study seismic wave propagation in velocity models which vary randomly in space. This stochastic approach seeks to understand the effects of small-scale complexity in the earth which cannot be resolved deterministically. These experiments have quantified the relationships between the statistical properties of the random heterogeneity and the measurable properties of high-frequency (1 Hz) seismograms. These simulations have been applied to the study of many features observed in actual high-frequency seismic waves, including: the amplitude and time decay of seismic coda, the apparent attenuation from scattering, the dispersion of waveforms, and the travel time and waveform variations across arrays of receivers. 相似文献
2.
3.
随机弹性介质中地震波散射衰减分析(英文) 总被引：2，自引：0，他引：2
地震波衰减一直是许多学科研究的热点,因为可以反映介质的特性。导致地震波衰减的因素很多,如:传播过程中由于能量扩散导致的几何衰减,固体岩石内部晶粒间相对滑移导致的摩擦衰减,岩石结构不均匀引起的地震波散射衰减。本文主要从统计的观点出发,通过多次数值模拟的方法研究纵波散射在随机弹性介质中所引发的衰减。首先用随机理论建立了二维空间随机弹性介质模型,然后用错格伪谱法的数值方法模拟了波在随机介质中的传播,再通过波场中虚拟检波器的记录,用谱比法估计了弹性波在随机介质中的散射衰减。不同非均匀程度随机弹性介质中的数值结果表明:介质不均匀程度越高,散射衰减越大;在散射体尺寸小于波长的前提下,不同散射体尺寸的计算结果说明:散射体尺寸越大,弹性波衰减越明显。最后提出了一种不均匀孔隙介质中流体流动衰减的方法。通过对随机孔隙介质中地震波的总衰减和散射衰减分别进行了计算,并定量得出了随机孔隙介质中流体流动衰减,结果表明:在实际地震频段下,当介质不均匀尺度101米量级时,散射衰减比流体流动衰减要大,散射衰减是地震波在实际不均匀岩石孔隙介质中衰减的主要原因。 相似文献
4.
在阿波罗月震记录中普遍存在着强烈持久的尾波信号,这样的波形特征无法用均匀分层月球模型解释.一个普遍被接受的解释是月震尾波由月球浅层结构对月震波的散射引起.我们采用基于交错网格的伪谱和有限差分混合方法模拟研究非均匀上月壳对月震波的散射效应,在此基础上解释月震尾波的形成机制,并估计出上月壳速度扰动的强度.我们发现,在均匀分层模型基础上,进一步考虑上月壳中的非均匀结构对月震波的散射效应,能有效地解释月震信号中强烈持久的尾波.我们认为月震尾波可能是由上月壳中的低波速、低衰减和散射这三个因素的共同作用所引起.采用不同的扰动标准差模拟上月壳的非均匀性,并比较模拟波形与真实月震图的相似程度,我们发现上月壳中速度扰动的标准差应该在3%到5%之间,很可能接近于3%. 相似文献
5.
6.
地下岩石结构/构造和矿物组成在小尺度上存在的不均匀性和各向异性会导致地震波速度的分形分布，进而衍生出自组织介质模型. 本文针对地球内部岩浆垂向侵入及横向溢流为特征的岩浆作用，开展自相似型各向异性自组织介质的地震波动力学响应研究. 基于二维自组织介质模型和声波方程，本文利用有限差分法模拟了不同强度岩浆侵入和溢流作用形成的自组织介质中地震波场，并进行了波场特征分析. 研究结果表明：（1）在横向溢流作用为主的自组织介质中，横向相关长度小于地震波波长的条件下，聚焦效应随着横向相关长度的增大而增强，地震波能量随之增强；反之，在横向相关长度大于地震波波长的条件下，散射效应随横向相关长度的增大而增强，地震波能量随之减弱. （2）在垂向侵入作用为主的自组织介质中，随着垂向相关长度增大，散射效应增强，地震波的能量也随之迅速减弱. （3）岩浆作用的自组织强度增强，地震波的能量增大，而中心频率无变化. （4）自组织介质具有相同的小尺度扰动，其速度梯度背景越小，地震波能量在炮点附近越集中. （5）当多层自组织介质共存，特征层位的波场特征仍以各自的自组织特性为主，因此利用实际资料中不同层位的地震反射特征可推测地下介质的自组织性质. 相似文献
7.
The perturbation method in elastic wave scattering 总被引：3，自引：0，他引：3
Ru-Shan Wu 《Pure and Applied Geophysics》1989,131(4):605-637
Methods of theoretical study in seismic wave scattering are reviewed with the emphasis on the perturbation method. Detailed analysis for weak scattering using Born approximation is given. For elastic random media, the mean square amplitudes of scattered waves are derived using a new approach by working directly in the spectrum domain. The conditions for the scalar wave approximation are obtained. The problem of sensitivity of fore- and backscattering to heterogeneities with different scales and properties (velocity or impedance) is discussed. 相似文献
8.
The concept of attenuation operators and complex velocities is applied to scattering attenuation in two and three dimensions, using the minimum-phase assumption for the attenuation operator. Acoustic 2D finite-difference computations of synthetic seismograms show, that the attenuation operator describes well the decay and lowpass filtering of the averaged wave form, which follows from averaging travel-time-corrected wave forms along the wave front. In the case of exponential random media, analytical forms of the attenuation operators and complex velocities are available. The complex velocities are incorporated into the reflectivity method. As an application, synthetic seismograms are presented for theS
n
wave, attenuated by lithospheric velocity and density fluctuations. The limitations of attenuation operators and complex velocities for scattering are also discussed. With these quantities it is not possible to model phenomena related to the scattered waves themselves, such as amplitude and travel-time fluctuations along the wave front, codas and precursors. 相似文献
9.
10.
地球深部圈层及沉积盆地是一种分区非均匀介质系统,其中不规则地层边界(含起伏地表)对地震波的主要特征有显著影响,而地层的随机非均匀性则主要影响地震波的散射和衰减特征.为了精确刻画不规则地层边界对地震波的反射、透射效应以及非均质体散射引起的地震波衰减效应,全局广义R/T递推传播矩阵法(GGRTM)被提出并逐步发展成为继有限元和有限差分方法之后的另一种复杂介质高精度地震波传播半解析求解方法.在已有的此类方法中,不规则边界均匀地层GGRTM法的优势在于对不规则地层边界的反射和透射效应的准确模拟,而非均质地层薄板化GGRTM法则能准确描述非均质体散射对地震波衰减的影响.本文吸收这两种已有方法的优势,提出了一种考虑非均匀介质、不规则边界的全局广义R/T递推传播矩阵混合方法,并将其用于对边界不规则、层内非均质的复杂模型的二维SH波场模拟.随后在本文方法与边界元法对比研究的基础上讨论了方法的模拟精度.研究结果表明本文提出的混合法是一种解决复杂模型高精度地震模拟的有效方法. 相似文献
11.
F. F. Pollitz 《Pure and Applied Geophysics》2002,159(9):2085-2112
-- I present a new algorithm for calculating seismic wave propagation through a three-dimensional heterogeneous medium using the framework of mode coupling theory originally developed to perform very low frequency (f < ~0.01т.05 Hz) seismic wavefield computation. It is a Greens function approach for multiple scattering within a defined volume and employs a truncated traveling wave basis set using the locked mode approximation. Interactions between incident and scattered wavefields are prescribed by mode coupling theory and account for the coupling among surface waves, body waves, and evanescent waves. The described algorithm is, in principle, applicable to global and regional wave propagation problems, but I focus on higher frequency (typically f S ~0.25 Hz) applications at regional and local distances where the locked mode approximation is best utilized and which involve wavefields strongly shaped by propagation through a highly heterogeneous crust. Synthetic examples are shown for P-SV-wave propagation through a semi-ellipsoidal basin and SH-wave propagation through a fault zone. 相似文献
12.
在给出真实模型和相应光滑背景模型的情况下,如何计算扰动模型(散射体)产生的散射波场是一个有实际意义的正演问题.在Gabor变换域描述散射体,且入射波场为短时宽带信号时,散射波场可以在频率域用高斯束或时间域用高斯波包描述.相对于波动方程方法,高斯束和高斯波包的计算效率更高;背景模型光滑时,高斯束和高斯波包方法的精度也接近波动方程方法.文中导出了声波假设下应用高斯束和高斯波包计算散射波的方法.测试分析了高斯波包的计算精度.给出了一般散射体的散射波模拟策略.同时针对一个理论模型完成了本文方法计算散射波的实验,实验结果表明高斯波包散射波计算方法是有效可行的. 相似文献
13.
In this work, we develop the indirect boundary element method (IBEM) to simulate the seismic site response in a realistic, large-scale 3-D sedimentary basin. Most previous applications of boundary element method have used full-space Green’s functions for wave propagation between element points. We use half-space Green’s functions, which include the seismic wavefield interactions at the free surface and require only the boundary elements of the basin interface. In this way, the size of the matrix equation for solution in the IBEM can be reduced to approximately a quarter of that using full-space Green’s functions. The site response modeling of the Granada basin in southern Spain using the IBEM shows that the basin-induced scattering waves were identified as propagating back and forth inside the basin. The scattered waves also generate surface waves that are weakly propagated outside of the basin. The wave propagation inside and outside of the basin shows different patterns. We observe that the scattered wave is locally amplified, and its propagation direction deviates from that of the incident waves propagation direction. Therefore, the computed seismic response in the basin could provide us with good estimates of the seismic motion. 相似文献
14.
针对山地地震勘探数据低信噪比问题,近地表散射波分离意义显得尤为突出,地震干涉测量法为此提供了一种技术手段.本文将地震干涉测量理论和散射理论结合起来,导出了近地表散射波地震干涉测量表达式,分为互相关型和褶积型表达式,它们由实际波场和背景波场干涉测量构成.根据近地表散射波分离理论,结合陆上地震勘探实际观测系统,采用褶积和反褶积混合型地震干涉测量配置,用实际地震资料展示了近地表散射波分离技术的应用效果.经过理论分析和砾石区实际资料试验,表明地震干涉测量不仅能分离测线上散射源产生的散射波,而且能分离部分侧面散射波.该技术的优点在于它适应于起伏地形和不均匀近地表结构,并且不需要起伏地形和近地表速度信息.为了从实际资料中消除近地表散射波,本文采用多道匹配滤波自适应减法,在砾石区见到较好效果. 相似文献
15.
Seismic wave field fluctuations below Lanzhou Comprehensive Test Ban Treaty Organization seismic array are determined and
interpreted to describe the scattering characteristics of the crust beneath the northeastern margin of Qinghai–Tibet plateau
in China. The frequency-dependent intensities of the mean and fluctuation wave fields are analyzed for 21 deep (206–632 km
depth) teleseismic events. The observed wave field fluctuations in the frequency range 0.5–2.5 Hz can be explained by the
scattering of the teleseismic P wave front at random media-type structures. The stable and narrow range of the fluctuation
parameters from these teleseismic events indicates evidence for scattered seismic phases which are generated inside the crust.
The reasonable structural models contain heterogeneities with 1–3% velocity fluctuations and 2.4–8.2 km correlation lengths
in the crust with 53 km thickness. 相似文献
16.
本义利用1979-1981年中国科学院地球物理研究所在随县-马鞍山所做的人工源地震探测的资料,采用二维射线追踪来拟合地震波的走时和振幅,取得了该地带的地壳与上地慢速度结构.结果表明,该地带地壳厚度为32-37km,莫霍面起伏变化显著.地壳速度的横向变化甚为强烈,地壳中部存在低速层.在巢湖以西宽数10km的范围内,地壳结构与其它地方明显不同,上地壳速度明显偏高,莫霍面上隆并有间断;在其两侧,明显可见有两条深断裂.结合其它地球物理场特征及大地构造概况,对结果进行分析,认为郯庐构造带南段被两条深断裂所夹,构造带宽度可达数10km,并在地质历史时期曾有过地幔物质上涌. 相似文献
17.
随县—马鞍山地带地壳与上地幔结构及郯庐构造带南段的某些特征 总被引：14，自引：2，他引：14
本义利用1979—1981年中国科学院地球物理研究所在随县—马鞍山所做的人工源地震探测的资料,采用二维射线追踪来拟合地震波的走时和振幅,取得了该地带的地壳与上地慢速度结构.结果表明,该地带地壳厚度为32—37km,莫霍面起伏变化显著.地壳速度的横向变化甚为强烈,地壳中部存在低速层.在巢湖以西宽数10km的范围内,地壳结构与其它地方明显不同,上地壳速度明显偏高,莫霍面上隆并有间断;在其两侧,明显可见有两条深断裂.结合其它地球物理场特征及大地构造概况,对结果进行分析,认为郯庐构造带南段被两条深断裂所夹,构造带宽度可达数10km,并在地质历史时期曾有过地幔物质上涌. 相似文献
18.
Sensitivity of seismic waves to structure 总被引：2，自引：0，他引：2
Luděk Klime? 《Studia Geophysica et Geodaetica》2012,56(2):483-520
We study how the perturbations of a generally heterogeneous isotropic or anisotropic structure manifest themselves in the wavefield, and which perturbations can be detected within a limited aperture and a limited frequency band. A short-duration broad-band incident wavefield with a smooth frequency spectrum is considered. In-finitesimally small perturbations of elastic moduli and density are decomposed into Gabor functions. The wavefield scattered by the perturbations is then composed of waves scattered by the individual Gabor functions. The scattered waves are estimated using the first-order Born approximation with the paraxial ray approximation. For each incident wave, each Gabor function generates at most 5 scattered waves, propagating in specific directions and having specific polarisations. A Gabor function corresponding to a low wavenumber may generate a single broad-band unconverted wave scattered in forward or narrow-angle directions. A Gabor function corresponding to a high wavenumber usually generates 0 to 5 narrow-band Gaussian packets scattered in wide angles, but may also occasionally generate a narrow-band P to S or S to P converted Gaussian packet scattered in a forward direction, or a broad-band S to P (and even S to S in a strongly anisotropic background) converted wave scattered in wide angles. In this paper, we concentrate on the Gaussian packets caused by narrow-band scattering. For a particular source, each Gaussian packet scattered by a Gabor function at a given spatial location is sensitive to just a single linear combination of 22 values of the elastic moduli and density corresponding to the Gabor function. This information about the Gabor function is lost if the scattered wave does not fall into the aperture covered by the receivers and into the legible frequency band. 相似文献
19.
裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征. 相似文献
20.
Numerical modelling ofSH wave seismograms in media whose material properties are prescribed by a random distribution of many perfectly elastic cavities and by intrinsic absorption of seismic energy (anelasticity) demonstrates that the main characteristics of the coda waves, namely amplitude decay and duration, are well described by singly scattered waves in anelastic media rather than by multiply scattered waves in either elastic or anelastic media. We use the Boundary Integral scheme developed byBenites
et al. (1992) to compute the complete wave field and measure the values of the direct waveQ and coda wavesQ in a wide range of frequencies, determining the spatial decay of the direct wave log-amplitude relation and the temporal decay of the coda envelope, respectively. The effects of both intrinsic absorption and pure scattering on the overall attenuation can be quantified separately by computing theQ values for corresponding models with (anelastic) and without (elastic) absorption. For the models considered in this study, the values of codaQ
–1 in anelastic media are in good agreement with the sum of the corresponding scatteringQ
–1 and intrinsicQ
–1 values, as established by the single-scattering model ofAki andChouet (1975). Also, for the same random model with intrinsic absorption it appears that the singly scattered waves propagate without significant loss of energy as compared with the multiply scattered waves, which are strongly affected by absorption, suggesting its dominant role in the attenuation of coda waves. 相似文献