Homotopy scattering series for seismic forward modelling with variable density and velocity

We have derived a convergent scattering series solution for the frequency-domain wave equation in acoustic media with variable density and velocity. The convergent scattering series solution is based on the homotopy analysis of a vectorial integral equation of the Lippmann–Schwinger type. By using the Green's function and partial integration, we have derived the vectorial integral equation of the Lippmann–Schwinger type that involves the pressure gradient field as well as the pressure field from the wave equation. The vectorial Lippmann–Schwinger equation can in principle be solved via matrix inversion, but the computational cost of matrix inversion scales like $N^3$, where $N$ is the number of grid blocks. The computational cost can be significantly reduced if one solves the vectorial Lippmann–Schwinger equation iteratively. A simple iterative solution is the Born series, but it is only convergent when the scattering potential is sufficiently small. In this study, we have used the so-called homotopy analysis method to derive an iterative solution for the vectorial Lippmann–Schwinger equation which can be made convergent even in strongly scattering media. The computational cost of our convergent scattering series scales as $N^2$. Our algorithm, which is based on the homotopy analysis method, involves a convergence control operator that we select using hierarchical matrices. We use a three-layer model and a resampled version of the SEG/EAGE salt model to show the performance of the developed convergent scattering series.     

基于波动方程理论的地震波场数值模拟方法综述

地震波场数值模拟的重要性日益凸显,其在地震勘探、各向异性介质属性研究、强地面振动预测、理论地震图合成等方面的应用日渐增多.基于波动方程理论的模拟方法有很多,如:有限差分法、伪谱法、有限元法、谱元法、边界元法等,而现有的综述性文献仅仅是关注某一种或是几种算法,还未有融入最新的文献后系统论述以上五种方法的综述性文献.鉴于此...     

近断层地震动破裂向前方向性与滑冲效应对隔震建筑结构抗震性能的影响

选择台湾集集地震和美国北岭地震的近断层地震动记录作为输入,考察了近断层地震动破裂向前方向性与滑冲效应引起的两种不同速度脉冲运动对单自由度体系和长周期橡胶支座隔震建筑结构抗震性能的影响.反应谱分析表明,破裂向前方向性与滑冲效应对工程结构地震响应的影响是随结构周期变化的.在中短周期段,含破裂向前方向性效应地震动的谱加速度值大于含滑冲效应地震动的谱加速度值;而在长周期段,含滑冲效应地震动的谱加速度大于含破裂向前方向性效应的谱加速度值.并且,与无脉冲地震动作用相比,含破裂向前方向性与滑冲效应脉冲的近断层地震动作用下隔震建筑的地震响应显著增大.滑冲效应引起的速度脉冲使隔震建筑底部的层间变形和楼层剪力明显增大,这意味着滑冲效应脉冲比向前方向性效应脉冲对长周期建筑结构的破坏更具危害性.     

地震波正演模拟中无单元Galerkin法初探

无单元Galerkin法作为一种新型有限元开始应用于地震勘探领域.地震波正演模拟结果的好坏取决于基函数类型、节点分布、节点编号方式、稳定条件以及边界条件等因素,通过对各个因素的分析可知:基函数的阶数越高,结果精度越高,但数据量和运算量就越大;不同的基函数联合使用时,会带来数值计算上的误差;节点可以不均匀分布,但必须满足空间采样要求;合理的节点编号有利于节约数据量和减少运算量;采用衰减边界条件,能有效克服无单元Galerkin法在边界上处理的难点.以上分析结果,对无单元Galerkin法在地震波正演模拟中的应用具有一定的指导意义.     

地震干涉技术被动源地震成像

被动源地震成像是基于地下反射波响应和透射波响应之间的关系,通过在地表接收由地下非人工震源发出的透射波响应,利用互相关运算来合成反射波响应的方法.合成的反射波响应,又称虚炮集记录等价于地表地震剖面记录.本文研究了由地下随机分布的被动源发出的透射波记录来合成虚炮集记录的基本原理,推导了反射波记录和透射波记录的关系公式;并通过数值模拟方法,分析了被动源成像的可行性.结果表明:利用被动源地震信息进行成像是可行的.     

T-matrix approach to seismic forward modelling in the acoustic approximation

Forward seismic modelling in the acoustic approximation, for variable velocity but constant density, is dealt with. The wave equation and the boundary conditions are represented by a volume integral equation of the Lippmann-Schwinger (LS) or Fredholm type. A T-matrix (or transition operator) approach from quantum mechanical potential scattering theory is used to derive a family of linear and nonlinear approximations (cluster expansions), as well as an exact numerical solution of the LS equation. For models of 4D anomalies involving small or moderate contrasts, the Born approximation gives identical numerical results as the first-order t-matrix approximation, but the predictions of an exact T-matrix solution can be quite different (depending on spatial extention of the perturbations). For models of fluid-saturated cavities involving large or huge contrasts, the first-order t-matrix approximation is much more accurate than the Born approximation, although it does not lead to significantly more time-consuming computations. If the spatial extention of the perturbations is not too large, it is practical to use the exact T-matrix solution which allows for arbitrary contrasts and includes all the effects of multiple scattering.     

海底天然气水合物OBS多分量地震正演模拟（英文）

根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震（OBS）共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。     

线性粘弹体中地震波场伪谱法模拟技术

以积分本构方程为基础,应用对应原理建立线性粘弹体中的波动方程,采用伪谱法进行数值模拟.实际计算结果表明该方法能有效模拟标准线性粘弹体中地震波的传播.     

正演模拟技术在地震采集设计中的应用

随着地震勘探开发的不断深入和发展,地震勘探的主战场逐渐向复杂地区转移.复杂地区既指地表条件复杂的地区,也指地下地质构造和地层条件复杂的地区,这些都对地震勘探提出了新的挑战和更高的要求.地震正演模拟正是开展此类问题研究的一个重要手段和方法.目前市场上具有正演模拟功能的软件大多是根据射线理论采用射线追踪的方法来完成正演模拟的,这种方法不能很好地反映地震波的动力学特征,特别是在复杂地区难以得到正确的结果.本文利用高阶有限差分有效克服了常规有限差分算法求解波动方程的频散问题,并以高效的OpenMP并行计算模式进行了并行优化,较大程度上提高了正演计算的速度和精度;同时实现了二阶Higdon边界条件,改善了边界吸收效果;也在一定程度上提高了计算速度.塔中地区主要目的层埋藏深,逆断层发育,地震反射特征复杂,增加了地震勘探的难度.本文依据该地区的地质模型,利用波动方程正演技术论证了该地区的地震采集观测系统,为该区地震采集观测系统的设计提供了科学依据.     

地震子波处理的二步法反褶积方法研究

针对玛湖斜坡区三块三维地震资料和赛汉塔拉凹陷二块三维地震资料连片处理中的特点,结合地质任务和处理目标要求,提出了地震数据连片处理中的地震子波处理的方法.该方法主要体现了两次反褶积,一次是采用地表一致性反褶积,将不同震源的频带拓宽到一个标准上;再一次采用相位校正反褶积,将不同震源的数据校正到相同相位上.为了保证提取的相位校正反褶积算子稳定,采用叠后地震道提取（主要考虑到叠后地震道信噪比高,算子稳定性强）,然后将该算子应用到叠前地震道,进行相位校正.     

VSP正反演综述

垂直地震剖面法是一种发展很迅速的地球物理方法,由于垂直地震剖面法是在介质内部点上直接观测,因而能够避开或减弱剖面上部低速带的干扰及外界噪声干扰,可以更直接、更有效地研究波的运动学和动力学特征,解决其地质问题,因此在地球物理勘探领域中最为活跃.本文从垂直地震剖面正反演的角度,并结合国内外实际的例子来介绍垂直地震剖面,同时讨论了垂直地震剖面未来发展的方向.     

地震多次波成像

地震资料含有各种类型多次波,而传统成像方法仅利用地震一次反射波成像,在地震成像前需将多次波去除.然而,多次波携带了丰富的地下结构信息,多次波偏移能够提供除反射波外的额外地下照明.修改传统逆时偏移方法,用包含一次反射波和多次波的原始记录代替震源子波,将SRME方法预测的表面多次波代替一次反射波作为输入数据,可将表面多次波成像.多次波成像的挑战和困难在于大量串扰噪声的产生,针对表面多次波成像中的成像噪声问题,将最小二乘逆时偏移方法与多次波分阶思想结合起来,发展可控阶数的表面多次波反演成像方法,有望初步实现高精度的表面多次波成像.在消除原始记录中的表面多次波后,通过逆散射级数方法预测得到层间多次波,将层间多次波作为逆时偏移方法的输入数据可将其准确归位到地下反射位置.数值实验表明,多次波成像能够有效地为地下提供额外照明,而可控阶表面多次波最小二乘逆时偏移成像方法几乎完全避免成像噪声.

     

Performance-based seismic design of nonstructural building components: The next frontier of earthquake engineering

With the development and implementation of performance-based earthquake engineering, harmonization of performance levels between structural and nonstructural components becomes vital. Even if the structural components of a building achieve a continuous or immediate occupancy performance level after a seismic event, failure of architectural, mechanical or electrical components can lower the performance level of the entire building system. This reduction in performance caused by the vulnerability of nonstructural components has been observed during recent earthquakes worldwide. Moreover, nonstructural damage has limited the functionality of critical facilities, such as hospitals, following major seismic events. The investment in nonstructural components and building contents is far greater than that of structural components and framing. Therefore, it is not surprising that in many past earthquakes, losses from damage to nonstructural components have exceeded losses from structural damage. Furthermore, the failure of nonstructural components can become a safety hazard or can hamper the safe movement of occupants evacuating buildings, or of rescue workers entering buildings. In comparison to structural components and systems, there is relatively limited information on the seismic design of nonstructural components. Basic research work in this area has been sparse, and the available codes and guidelines are usually, for the most part, based on past experiences, engineering judgment and intuition, rather than on objective experimental and analytical results. Often, design engineers are forced to start almost from square one after each earthquake event: to observe what went wrong and to try to prevent repetitions. This is a consequence of the empirical nature of current seismic regulations and guidelines for nonstructural components. This review paper summarizes current knowledge on the seismic design and analysis of nonstructural building components, identifying major knowledge gaps that will need to be filled by future research. Furthermore, considering recent trends in earthquake engineering, the paper explores how performance-based seismic design might be conceived for nonstructural components, drawing on recent developments made in the field of seismic design and hinting at the specific considerations required for nonstructural components.     

Performance-based seismic design of nonstructural building components:The next frontier of earthquake engineering

With the development and implementation of performance-based earthquake engineering,harmonization of performance levels between structural and nonstructural components becomes vital. Even if the structural components of a building achieve a continuous or immediate occupancy performance level after a seismic event,failure of architectural,mechanical or electrical components can lower the performance level of the entire building system. This reduction in performance caused by the vulnerability of nonstructural components has been observed during recent earthquakes worldwide. Moreover,nonstructural damage has limited the functionality of critical facilities,such as hospitals,following major seismic events. The investment in nonstructural components and building contents is far greater than that of structural components and framing. Therefore,it is not surprising that in many past earthquakes,losses from damage to nonstructural components have exceeded losses from structural damage. Furthermore,the failure of nonstructural components can become a safety hazard or can hamper the safe movement of occupants evacuating buildings,or of rescue workers entering buildings. In comparison to structural components and systems,there is relatively limited information on the seismic design of nonstructural components. Basic research work in this area has been sparse,and the available codes and guidelines are usually,for the most part,based on past experiences,engineering judgment and intuition,rather than on objective experimental and analytical results. Often,design engineers are forced to start almost from square one after each earthquake event: to observe what went wrong and to try to prevent repetitions. This is a consequence of the empirical nature of current seismic regulations and guidelines for nonstructural components. This review paper summarizes current knowledge on the seismic design and analysis of nonstructural building components,identifying major knowledge gaps that will need to be filled by future research. Furthermore,considering recent trends in earthquake engineering,the paper explores how performance-based seismic design might be conceived for nonstructural components,drawing on recent developments made in the field of seismic design and hinting at the specific considerations required for nonstructural components.     

井间地震波动方程正演数值模拟及其应用

井间地震正演模拟是在给定地层参数及边界条件的情况下求波动方程数值解的过程.本文详细推导了井间二维声波方程四阶差分方程;根据假设的包含楔形体和丘状体的理论地层模型,对其波场进行了模拟计算分析;据胜利油田垦71区块两口井的实际声波测井资料建立了精细速度模型,对其正演模拟的合成记录波场进行分析,并把合成记录反射波成像与实际的过井三维地震剖面进行了对比,验证了所设计观测系统的合理性.     

Study of crustal seismic velocity in the Weihe fault depression basin by seismic tomographic imaging

ＳｔｕｄｙｏｆｃｒｕｓｔａｌｓｅｉｓｍｉｃｖｅｌｏｃｉｔｙｉｎｔｈｅＷｅｉｈｅｆａｕｌｔｄｅｐｒｅｓｉｏｎｂａｓｉｎｂｙｓｅｉｓｍｉｃｔｏｍｏｇｒａｐｈｉｃｉｍａｇｉｎｇＧＵＡＮＧＹＩＮＧＸＵＥ１）（薛广盈）ＹＵＮＹＵＤＩＮＧ２）（丁韫玉）ＺＨ...