地震波散射:理论与应用

传统的球面对称(或层状构造)地球模型正经历着一场革命。地球被揭示从地壳、地幔到地核到处都呈现多尺度的横向非均匀性。这些具有不同尺度的非均匀体对地震波具有不同的效应。速度和密度的非均匀体能改变波形,引起走时和振幅的起伏以及产生直达波的视衰减。地球岩石层的非均匀体还能产生P尾波、S尾波和Lg尾波等。核-幔边界附近的非均匀体能产生对PKP波的散射而成为PKIKP的前驱波,PKKP波的散射波可成为主震相自己的前驱波。近源和近台站的复杂构造可通过共振、散射来改变地震波形。粗糙地形或粗糙界面能造成体波和面波的耦合。地壳内规则排列的裂隙可产生有效各向异性而使S波分裂。由三维非均匀体所引起的地震波的变化,在广义上被称为地震波散射。近十多年来,由于高质量的高频数字地震资料的逐渐增多,对地震波散射的研究在急速发展并引起了越来越多的地震学家、工程学家和勘探地球物理学家的兴趣。本文综述了地震波散射的基本理论和在这一领域各方面的最新进展。其内容大致为:一、地球横向非均匀性的谱及各种散射态式二、地震波散射的研究方法1.理论研究(1)不连续非均匀介质的边界匹配方法(2)弱散射的微扰法(3)高频近似法(4)随机方法和非均匀体的统计特性2.数值模拟和物理模型试验3.野外观测三、弹性波散射的基本特征和标量波近似1.弹性波瑞雷散射2.弹性波瑞雷-甘斯散射3.随机介质的弹性波散射四.地震波散射的表现和应用1.透射起伏2.尾波产生及包络消减3.散射衰减4.核-幔边界附近的散射5.地表地形、近地表结构和深部构造引起的散射6.裂缝散射和有效各向异性7.散射和非线性。     

井孔声场中的弹性波散射效应及其应用

地下岩石普遍存在着非均匀性,当井孔周围岩石的非均匀体的尺度(如压裂产生的裂缝系统)与声波波长相当时,会产生很强的散射波,因此可以通过声波测井的散射效应评价地层的非均匀性.文章基于弹性介质散射体在统计意义上的分布规律,建立了三维井孔非均匀介质模型,数值模拟了非均匀地层偶极声波测井的散射效应,并分析了散射波的衰减特性.模拟结果表明:与均匀介质相比,井孔周围非均匀体产生了明显的散射,一个显著的特征是在直达波之后出现了较强的尾波,且尾波的频率随着时间逐渐降低.在此基础之上,根据尾波的衰减特性提出了一种利用散射效应评价地层非均匀性的方法,并将该方法应用于致密储层的压裂效果评价.应用结果表明压裂后偶极声波测井散射效应产生的尾波预示着井孔周围岩石的体积改造,可以用来评价致密储层压裂效果,这为今后利用声波测井资料评价非常规储层(如页岩储层)的压裂效果提供了一条新的思路.     

随机弹性介质中地震波散射衰减分析(英文)

地震波衰减一直是许多学科研究的热点,因为可以反映介质的特性。导致地震波衰减的因素很多,如:传播过程中由于能量扩散导致的几何衰减,固体岩石内部晶粒间相对滑移导致的摩擦衰减,岩石结构不均匀引起的地震波散射衰减。本文主要从统计的观点出发,通过多次数值模拟的方法研究纵波散射在随机弹性介质中所引发的衰减。首先用随机理论建立了二维空间随机弹性介质模型,然后用错格伪谱法的数值方法模拟了波在随机介质中的传播,再通过波场中虚拟检波器的记录,用谱比法估计了弹性波在随机介质中的散射衰减。不同非均匀程度随机弹性介质中的数值结果表明:介质不均匀程度越高,散射衰减越大;在散射体尺寸小于波长的前提下,不同散射体尺寸的计算结果说明:散射体尺寸越大,弹性波衰减越明显。最后提出了一种不均匀孔隙介质中流体流动衰减的方法。通过对随机孔隙介质中地震波的总衰减和散射衰减分别进行了计算,并定量得出了随机孔隙介质中流体流动衰减,结果表明:在实际地震频段下,当介质不均匀尺度101米量级时,散射衰减比流体流动衰减要大,散射衰减是地震波在实际不均匀岩石孔隙介质中衰减的主要原因。     

地层衰减对地震波速度逆散射反演的影响研究

在利用地震波数据进行地球物理反演时,地层对地震波的吸收衰减效应会对地层物性参数的准确反演产生较大的影响,因此利用黏弹性声波方程进行反演更符合实际情形.本文在考虑地层衰减效应进行频率空间域正演模拟的基础上,提出基于黏弹性声波方程的频率域逆散射反演算法并对地震波传播速度进行反演重建,在反演过程中分别用地震波传播复速度和实速度来表征是否考虑地层吸收衰减效应.基于反演参数总变差的正则化处理使反演更加稳定,在反演中将低频反演速度模型作为高频反演的背景模型进行逐频反演,由于单频反演过程中背景模型保持不变,故该方法不需要在每次迭代中重新构造正演算子,具有较高的反演效率;此外本文在反演过程中采用了基于MPI的并行计算策略,进一步提高了反演计算的效率.在二维算例中分别对是否考虑地层吸收衰减效应进行了地震波速度反演,反演结果表明考虑衰减效应可以得到与真实模型更加接近的速度分布结果,相反则无法得到正确的地震波速度重建结果.本文算法对复杂地质模型中浅层可以反演得到分辨率较高的速度模型,为其他地震数据处理提供比较准确的速度信息,在地层深部由于地震波能量衰减导致反演分辨率不太理想.     

介质非均匀性参数对散射波包络的影响

地球内部不同尺度的非均匀体会引起高频地震波的散射,这些散射携带了大量关于地球内部介质非均匀精细结构的信息,可通过研究地震波散射信息来获取地球内部介质结构非均匀性.非均匀介质速度微扰动是造成地震波散射的最主要起因,也是引起地震波尾波和散射波包络展宽的主要因素.     

黏弹性VTI介质中敏感度核函数：地震波复走时关于弹性模量及Q值偏导数的计算

实际地层中地震波传播普遍存在速度和衰减各向异性现象,研究黏弹各向异性介质中高频地震波传播理论有助于揭示地震波的传播特征.本文针对黏弹性VTI介质,从Christoffel矩阵的解析特征值出发推导出qP、qSV和qSH波的复相速度和复射线速度的解析表达式,并应用实射线追踪方法确定出均匀复射线速度矢量,由此计算出实射线速度和实射线衰减以及实射线品质因子.基于非均匀复相速度和均匀复射线速度的解析表达式,推导了实射线慢度和实射线衰减关于黏弹性模量(包括弹性模量和Q值)的敏感度核函数,该敏感度核函数反映各个黏弹性模量对地震波复走时的影响程度.不同岩石样本的数值计算结果显示,实走时对弹性模量更为敏感,而射线衰减(虚走时)对弹性模量和Q值的敏感程度相当.本研究可为黏弹性VTI介质中地震射线追踪和复走时层析成像提供理论基础.     

地震学百科知识(十二)——地震波散射

正广义地说,由地球介质的任何三维非均匀性引起的地震波场畸变均可称为地震波散射。不过,通常地震波散射是指由地球小尺度(与波长相比)非均匀性引起的、不适合用射线理论(几何射线)研究的地震波场畸变现象。例如,波场传播方向的偏折和分裂,P波和S波的相互转换,地震尾波的产生,地震波的走时和振幅横跨观测台阵的随机起伏变化等。1地球介质的不均匀性及地震波散射随着对地球结构认识的深入,经典的球     

基于常规测井数据的纵波本征衰减估计

地震波本征衰减反映了地层及其所含流体的一些特性,对油气勘探开发有重要意义.已有的理论研究与实验发现,地震频带内的衰减主要与中观尺度(波长与颗粒尺度之间)的斑状部分饱和、完全饱和岩石弹性非均匀性情况下波诱导的局部流体流有关.这种衰减与岩石骨架、孔隙度及充填流体的性质密切相关.本文着重讨论均匀流体分布、斑状或非均匀流体分布两种情况下部分饱和岩石的纵波模量差异.以经典岩石物理理论和衰减机制认识为基础,通过分析低频松弛状态、高频非松弛状态岩石的弹性模量,讨论储层参数(如孔隙度、泥质含量以及含水饱和度等)与纵波衰减之间的确定性关系.上述方法与模型在陆相砂泥岩地层与海相碳酸盐岩地层中的适用性通过常规测井资料得到了初步验证.     

云南地区S波非弹性衰减Qem值研究/

地震波衰减性质的研究是地震学研究的重要课题.地震台站(地震仪)记录到的地震波,包含了地震震源效应、地震波的传播路径效应、台站场地响应及仪器响应.在使用地震波资料研究震源性质时,必须要扣除地震波传播路径效应、台站场地响应及仪器响应的影响.地震波传播的路径效应(地震波衰减),除了随距离存在几何衰减外,还有一个重要的影响因素--即介质的非弹性衰减,用介质品质因子Q值来度量.Q值是地球介质的基本物理参数之一,是对地震进行定量和研究震源性质(如震源参数的测定等)所必需的重要参数,在震源物理和工程地震研究中具有重要应用价值.     

平面P、SV和Rayleigh波入射下沉积河谷地震响应IBEM宽频模拟

沉积河谷对地震波的高频散射规律尚不明确。采用一种高精度间接边界元法(IBEM)用于沉积河谷对平面P、SV和Rayleigh波的二维散射宽频求解分析。数值结果表明:IBEM可高效精确地模拟地震波的宽频散射;中低频段(无量纲频率η5.0),河谷地形对地震波的放大效应显著;高频波段(η10.0),河谷放大效应减弱甚至出现缩幅效应。随河谷深度增加,位移放大频段的带宽逐渐减小,第一峰值频率降低,且在低频段频谱曲线振荡剧烈。阻尼比对沉积河谷地表位移幅值具有显著影响,尤其是在高频段衰减作用明显。实际沉积河谷场地地震反应分析,需充分考虑波型、入射波频率和角度、河谷深宽比等因素,以更科学地进行震害解释和地震安全性评价。     

复杂地表边界元-体积元波动方程数值模拟

复杂近地表引起来自深部构造的地震反射信号振幅和相位的异常变化,是影响复杂近地表地区地震资料品质的主要原因.本文采用边界元-体积元方法,通过求解含复杂地表的波动积分方程,来模拟地震波在复杂近地表构造中的传播.其中,边界元法模拟地形起伏和表层地质结构对地震波传播的影响;体积元法模拟起伏地表下非均质低降速层的影响.与其他数值...     

Seismic Wave Propagation and Excitation in Multi-layered Media with Irregular Interfaces Part (Ⅰ): SH Case

Study of seismic wave excitation and propagation in laterally heterogeneous media was an active and important subject in seismology in the past two decades, numerous analytical and numerical efforts have been made in this research field. In this article, I have, first, made a brief review on those developments and then introduced and summarized a unified and efficient method, global generalized reflection-transmission (abbreviated to R/T thereafter) matrices method, for synthetic seismograms in multi-layered media with irregular interfaces developed by the author [24~26]. As demonstrated in this article, this method could be regarded as an extension of the generalized R/T coefficients method for the horizontally layered case [2,5] to the layered media with irregularly shaped interfaces by incorporating the T matrices technique [27,28]. Because of the use of a recursive scheme in computing the global generalized R/T matrices, this method is efficient, particularly for the case with a large number of irre     

THE APPLICATION OF PML BOUNDARY CODITIONS IN THE SIMULATION OF SEISMIC WAVE BY THE HIGH-ORDER STAGGERED-GRID FINITE DIFFERENCES AT THE TRANSVERSELY ISOTROPIC MEDIA

In the realm of the numerical simulation, finite difference method and finite element method are more intuitive and effective than other simulation methods. In the process of simulating seismic wave propagation, the finite differences method is widely used because of its high computational efficiency and the advantage of the algorithm is more efficient. With the demand of precision, more and more researchers have proposed more effective methods of finite differences, such as the high-order staggered-grid finite differences method, which can restore the actual process of wave propagation on the premise of ensuring accuracy and improving the efficiency of operation. In the past numerical simulation of seismic wave field, different models of isotropic medium are mostly used, but it is difficult to reflect the true layer situation. With the research demand of natural seismology and seismic exploration, the research on anisotropic media is more and more extensive. Transversely isotropic(TI)media can well simulate the seismic wave propagation in the formation medium, such as gas-bearing sandstone, mudstone, shale et al., the character of TI media is reflected by introducing the Thomsen parameters to reflect its weak anisotropy of vertical direction by using Thomson parameter. Therefore, studying the process of seismic wave propagation in TI media can restore the true information of the formation to the greatest extent, and provide a more reliable simulation basis for the numerical simulation of seismic wave propagation. In the geodynamic simulation and the numerical simulation of the seismic wave field, under the limited influence of the calculation area, if no boundary conditions are added, a strong artificial boundary reflection will be generated, which greatly reduces the validity of the simulation. In order to minimize the influence of model boundaries on the reflection of seismic waves, it is often necessary to introduce absorbing boundary conditions. At present, there are three types of absorption boundary conditions: one-way wave absorption boundary, attenuation absorption boundary, and perfectly matched layer(PML)absorption boundary. In terms of numerical simulation of seismic waves, the boundary absorption effect of PML is stronger than the first two, which is currently the most commonly used method, and it also represents the cutting-edge development direction of absorption boundary technology. The perfectly matched layer absorbing boundary is effectively applied to eliminating the reflective waves from model boundaries, but for transversely isotropic medium, the effect of the absorbing is not very well. For this reason, the elastic dynamic wave equations in transversely isotropic media are derived, and we describe a second-order accurate time, tenth-order accurate space, formulation of the Madariaga-Virieux staggered-grid finite difference methods with the perfectly matched layer(PML)are given. In addition, we have established vertical transversely isotropic(VTI)media and arbitrary inclined tilted transversely isotropic(TTI)media models, using a uniform half-space velocity model and a two-layer velocity model, respectively. By combining the actual geoscience background, we set the corresponding parameters and simulation conditions in order to make our model more research-oriented. When setting model parameters, different PML thickness, incident angle, source frequency and velocity layer models were transformed to verify the inhibition of boundary reflection effect by PML absorption boundary layer. The implementations of this simulation show that the formula is correct and for the transversely isotropic(TI)media of any angular symmetry axis, when the thickness of the PML layer reaches a certain value, the seismic wave reflection effect generated by the artificial boundary can be well suppressed, and the absorption effect of PML is not subject to changes in incident angle and wave frequency. Therefore, the results of our study indicate that our research method can be used to simulate the propagation process of seismic waves in the transversely isotropic(TI)media without being affected by the reflected waves at the model boundary to restore the actual formation information and more valuable geological research.     

基于多辛结构谱元法的保结构地震波场模拟

近年来构造高精度、高效且具有长时程跟踪能力的保结构算法已逐渐成为地震波模拟算法发展的重要方向之一. 本文基于谱元法(SEM)进行空间域离散结合新推导的三阶辛算法(NTSTO)进行时间域离散,构造了一种具有时-空保结构特性的新算法. 本文给出的多组数值试验对比结果表明,本算法无论在内存消耗、稳定性及计算耗时,还是长时程跟踪能力方面都有上佳的表现; 另外,本文给出的起伏地表多层介质模型的数值算例验证了该算法处理复杂几何形状和复杂介质时的有效性. 该多辛结构谱元法的发展将为长时程地震波传播的计算及模拟提供更为广泛而有效的选择.     

分段光滑曲线边界波动方程数值模拟研究

矩形网格有限差分法在地震波传播数值模拟方面具有计算速度快的显著优势,但该方法在处理复杂边界问题上存在着效率低的严重缺陷.本文针对分段光滑曲线边界定义了尖点处的一种正则导数,给出了矩形网格情形分段光滑曲线网格边界点法向导数的一种插值计算方法.采用矩形网格有限差分法对复杂边界地球介质模型进行地震波场数值模拟,并采用波场系列快照技术揭示地震波在起伏地表和复杂介质中的传播规律.模拟结果表明：法向导数插值计算方法为矩形网格有限差分法处理复杂边界提供了有效途径,采用波场系列快照技术可以清晰地展现地震波在反射界面的反射和透射规律、在尖点的绕射规律以及在自由表面的直达波和多次反射规律.     

横向非均匀介质中弹性地震波传播的数值模拟

应用混合变量弹性动力学方程和线性常微分方程组的矩阵指数解法，将层状介质中广泛应用的弹性波传播矩阵解法推广至横向非均匀介质，给出了一种可计算复杂地质体中弹性波传播的广义传播矩阵数值解法。该方法可模拟任意震源及所产生的各种体波、面波，数值结果表明具有很高的计算精度。     

三维层状孔隙介质中弹性波的一种积分表达式I：理论

孔隙介质弹性波传播理论在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用.而孔隙介质中的弹性波受孔隙度、渗透率、流体黏滞系数等参数的影响,因此研究波场的传播特征将有助于分析和提取这些信息.本文在Biot理论的基础上,针对三维层状孔隙介质模型,利用在合成理论地震图的研究中已经被证实具有稳定、高效且适用范围较广的Luco-Apsel-Chen（LAC）广义反透射方法,给出了弹性波场的一种积分形式的半解析解,可通过数值方法高效、准确地计算层状孔隙介质中的理论波场,所以该积分形式的半解析解可为三维层状孔隙介质波场传播特征的理论数值模拟研究提供一种新的途径和手段.     

Seismic wavefield modeling in media with fluid-filled fractures and surface topography

We present a finite difference (FD) method for the simulation of seismic wave fields in fractured medium with an irregular (non-flat) free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations (in a curvilinear coordinate system) in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry (HTI). Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures.