P波广角入射的透射P波特征研究（英文）

P波入射到地下反射界面时,如果入射角超过临界角,透射P波可能会出现极性反转。本文首先从Zoeppritz方程出发,研究了广角入射条件下的反射系数和透射系数特征,发现满足一定条件时,随着入射角的增加,透射P波系数的特征曲线从第一象限经过原点进入第三象限,透射P波系数发生极性反转,它的辐角发生突变。文中给出了透射P波系数为零时入射角的表达式及其近似公式,通过对单界面模型进行有限差分正演,证明了透射P波系数为零的存在性。储层从饱含水到饱含气变化时,储层的P波速度和密度明显降低,通过单界面储层模型的透射P波的合成记录,说明了储层饱含气时透射P波系数可能发生极性反转。经理论分析和数值模拟得到了一些结论,这些结论为井间地震勘探的应用奠定了理论基础。     

甚低纬哨声低电离层透射过程的数值模拟

利用全波解算法模拟哨声波束在甚低纬地区黎明前低电离层透射的三维能量分布，依据波场能量和偏振分布及其对波参量和电子浓度剖面的依赖特征，分析了哨声透射、反射及与大地－电离层波导耦合过程．结果表明，哨声模波存在于９０ｋｍ以上高度，吸收、反射、波束扩展及波模转换主要发生于电离层底部８０－９０ｋｍ区间；到达地面的透射能量密度衰减２０ｄＢ以上，透射衰减随频率变化不大，但随波入射角呈不对称变化；距透射区１５０ｋｍ以外区域的测向方位角有很大偏差；入射波能量的很少一部分（对５ｋＨｚ约为－２５ｄＢ）被反射并激发起哨声模波，反射波束能量集中于入射波束附近，并随频率下降而迅速增强．计算也表明，地面接收到的甚低纬哨声回波可能与使回波向极侧偏移的电离纬向梯度有关．     

用斜入射超声波评价水泥胶结质量的方法研究(英文)

超声波水泥胶结质量评价测井是目前较为广泛使用的固井质量评价方法。将传统的脉冲回波技术与超声波斜入射方法相结合,可以形成评价固井质量的新技术。本文使用有限差分方法模拟研究了超声波斜入射时所激发出的波形特征与水泥性质以及胶结状况之间的关系。通过对三种常用水泥的数值模拟,得到了AO模式和SO模式Lamb波在不同胶结水泥及胶结情况下的声衰减系数;通过对不同厚度胶结水泥的数值模拟,得到了次生AO模式Lamb波和AO模式Lamb波的到时差与水泥性质及厚度的关系,并理论推导了水泥层厚度与时间差之间的关系,且从有限差分时域波形中提取的时间差与由理论公式计算结果能够较好的吻合。因此利用AO模式Lamb波到时及衰减信息可以在一定程度上评价水泥胶结状况及水泥厚度。     

转换波和非转换波与射线参数p的关系

本文对Aki固体-固体分界面的反射和透射公式进行适当变换,使之成为一组适合不同坐标系下,固体、液体、气体和真空之间所有可能形成的分界面上,入射P、SV和SH波所形成的反射和透射公式,从而把不同边界条件的反射和透射方程统一起来,并且阐明了转换波和非转换波动力学特征的差异.指出:在不同弹性介质分界面上入射平面弹性波,转换波的反射系数和透射系数是射线参数p的奇函数;非转换波的反射系数和透射系数是射线参数p的偶函数.根据奇函数和偶函数的性质可见:垂直入射时不存在转换波.同时也为简化Zoeppricz方程,开展AVO分析打开新思路.     

表面多次波衰减的研究现状与进展

表面多次波是一种相干干扰,在海洋地球物理勘探中表面多次波尤为明显.本文给出了表面多次波衰减方法的研究现状的系统阐述,并分析了表面多次波衰减研究的发展趋势.表面多次波衰减主要有滤波法和波动理论方法.对于滤波法给出了广泛应用的抛物Radon变换法的实例,指出了抛物Radon 变换法存在滤波时能量泄漏的问题.同时本文分析了基于波动理论的处理优势和存在的问题,该方法可以处理任意介质的情形,但是对于数据的采集要求比较严格,另外在进行表面多次波衰减处理以前,要进行一系列的预处理,这些预处理的误差直接决定着多次波衰减的效果.最后,本文给出了表面多次波衰减的发展趋势——将滤波法和波动理论法结合来实现多次波的衰减.两类方法如何有效地结合是该类方法一个重要的研究方向.     

地下洞室群对地面运动影响问题的级数解答-P波入射.

采用波函数展开法，给出了平面P波入射下半空间中洞室群对地面运动影响问题的一个级数解答. 数值结果表明，地下洞室群对附近地面运动具有显著的放大作用．建议地铁等地下工程在规划和设计时，考虑工程建设后对沿线地面运动的影响．      

地形对地震烈度衰减的影响

本文采用阶跃截面固体的一维弹性波传播模型,研究了地形对地震烈度衰减的影响。采用固体入射纵波来模拟无地形影响的地震烈度衰减;采用固体透射波来模拟有地形影响的地震烈度衰减。利用阶跃截面处的放松连续条件,得出了烈度与高度差的近似定量关系。同时,以新疆地区逆冲走滑型地震为例,讨论了烈度与高度差的近似定量关系描述地形对地震烈度衰减影响的有效性。     

海洋地震资料多次波衰减方法综述

多次波衰减方法的研究一直都是海洋地震资料处理中的重点和难点.为了能够选择合理有效的多次波衰减方法,提高海洋地震资料处理的信噪比和成像质量,本文着重论述了当前较为成熟的两大类多次波衰减方法——滤波法和波动理论法,并就多次波新技术的发展进行了总结.本文首先就滤波法和波动理论法的基本原理和研究思路进行了分析与探讨,对两大类方法中较为经典的代表性方法进行了详细的公式推导和原理解释,并且在此基础上比较分析了不同多次波衰减方法各自的适用性和优缺点.其次,在论述经典多次波衰减方法的基础上,对近些年来多次波衰减新技术的发展进行了归纳和总结,包括综合法、稀疏反演法和多次波成像技术等,指出了经典方法存在的局限性和新技术发展的必要性.最后,在总结近年来国内外多次波衰减问题研究进展的基础上,指出了该领域研究中当前存在的问题和不足,并且针对这些问题对未来海洋多次波衰减技术的发展做出了展望.     

基于混合边界条件的有限单元法GPR正演模拟

从Maxwell方程组出发,推导了探地雷达(GPR)有限元波动方程.阐述了透射边界条件和Sarma边界条件的原理,推导了这两种边界条件的理论公式;通过在衰减层内加入过渡带优化了Sarma边界条件的加载方法,压制了介质区和衰减层交界面处的人为反射.考虑到透射边界条件与Sarma边界条件不同的理论机制,提出了一种结合透射边界条件和Sarma边界条件的混合边界条件,它利用Sarma边界条件对到达边界区域的GPR波能量衰减功能和透射边界对GPR波能量的透射功能,使GPR波经过Sarma边界条件的衰减吸收后,再通过透射边界条件将剩余能量透射出去,集成了二者的优势.并以二维均匀模型中的中心脉冲激励源方式为例,通过Matlab程序实现,以GPR的全波场快照的直观方式,对比了有、无边界条件及不同边界条件对人工截断边界的处理效果,说明了该混合边界条件对到达截断边界处的GPR波的处理优于单一边界条件.最后,以基于混合边界条件的有限单元法对两个典型的GPR地电模型进行了正演模拟,指导了GPR数据处理与工程实践.     

双相各向异性介质中弹性波传播有限元方程及数值模拟

基于Biot双相各向异性介质理论和动态问题的哈密顿原理，推导出任意双相各向异性介质中弹性波传播的有限元方程，并给出双相各向异性介质中弹性波有限元方程的数值解法．最后进行有限元法的数值模拟，对双相各向异性介质中弹性波传播特征进行了模拟与分析．       

Multi—transmitting formula for attenuating waves

The MTF is extended to case of attenuating incident wave by introducing an attenuation coefficient.The reflection coefficients of this modified MTF and MTF are evaluated and compared when an attenuating wave impinges on the boundary,and the results demonstrate that MTF can be used to absorb slightly attenuating waves and the modi-fied MTF is more capable of absorbing heavily attenuating waves than MTF.The accuracy of modified MTF is also tested by numerical examples of fluid saturated porous media.     

多人工波速优化透射边界在谱元法地震波动模拟中的应用

将作者最近发展的多人工波速优化透射边界(记为ca j-MTF)应用于高精度谱元法的地震波动模拟中,并与经典的廖氏透射(MTF)边界、完美匹配层(PML)边界、黏弹性边界以及一阶旁轴近似边界进行了比较分析.结果显示:①ca j-MTF边界与MTF边界在形式上非常接近,它继承了后者公式简单、易于实现、精度可控、计算量低以及...     

透射边界高频耦合失稳机理及稳定实现——SH波动

就大型近场波动的高效数值模拟而言,稳定实现高阶人工边界是一个尚未圆满解决的问题.本文针对使用多次透射公式的SH波动集中质量有限元模拟,依据GKS定理的群速度解释,进一步阐明了人工边界与内域离散格式耦合所导致高频失稳的机理,即两者支持群速度指向内域的外行高频平面谐波,波动能量自发地从人工边界进入內域,从而导致失稳,而这类谐波是由集中质量有限元离散引入的.本文提出了消除此种耦合失稳的一种方法:通过修改有限元刚度阵来改变内域离散格式,并保证修改格式的精度不低于原有格式的精度.理论分析和数值实验表明此法能稳定实现透射边界.本文研究结果具有推广应用前景.     

热弹性介质中波传播特征

深层-超深层油气地震勘探涉及高温介质地震波传播问题,热弹介质参数对地震波传播有重要影响.含弛豫时间修正项的Lord-Shulman双曲型耦合热弹波动方程从理论上预测了热弹性介质中存在快纵波、慢纵波(一种准静态慢纵波,简称热波)和横波的传播,两个纵波为热耗散衰减波而横波不受介质热特性的影响.本文结合平面波频散分析和格林函数法数值模拟,详细研究两个热耗散衰减波的频散和衰减特征,着重分析热导率、热膨胀系数及比热的变化对波速和衰减的影响.研究表明热导率作为主要参数决定了波速与衰减的临界变化,热膨胀系数对波速和衰减的幅度有明显影响,比热则兼顾了前两个热弹系数的影响特征.最后,利用热弹性动力学频率域的二阶格林函数进行波场快照数值模拟,展示热弹性介质中纵波、横波和热波的传播行为.     

Arbitrary-order Taylor series expansion-based viscoacoustic wavefield simulation in 3D vertical transversely isotropic media

Taking the anisotropy of velocity and attenuation into account, we investigate the wavefield simulation of viscoacoustic waves in 3D vertical transversely isotropic attenuating media. The viscoacoustic wave equations with the decoupled amplitude attenuation and phase dispersion are derived from the fractional Laplacian operator and using the acoustic approximation. With respect to the spatially variable fractional Laplacian operator in the formulation, we develop an effective algorithm to realize the viscoacoustic wavefield extrapolation by using the arbitrary-order Taylor series expansion. Based on the approximation, the mixed-domain fractional Laplacian operators are decoupled from the wavenumbers and fractional orders. Thus, the viscoacoustic wave propagation can be conveniently implemented by using a generalized pseudospectral method. In addition, we perform the accuracy and efficiency analyses among first-, second- and third-order Taylor series expansion pseudospectral methods with different quality factors. Considering both the accuracy and computational cost, the second-order Taylor series expansion pseudospectral method can generally satisfy the requirements for most attenuating media. Numerical modelling examples not only illustrate that our decoupled viscoacoustic wave equations can effectively describe the attenuating property of the medium, but also demonstrate the accuracy and the high robustness of our proposed schemes.     

Modeling of nonmigrating tides in the middle atmosphere

On the basis of calculations using the general circulation model of the middle and upper atmosphere, the relative role of sources of nonmigrating tides distributed in atmosphere has been investigated. It is shown that in winter, when planetary waves in stratosphere are well developed, the main contribution to the generation of nonmigrating tides is caused by nonlinear interaction between migrating tides and a quasi-stationary planetary wave with zonal wave number 1 (SPW1). Taking into account the longitudinal ozone inhomogeneities in the model leads to the occurrence of additional sources of nonmigrating tides caused by longitudinally inhomogeneous heating of the atmosphere, the contribution of which can be comparable to that from nonlinear interaction under an attenuating amplitude of SPW1 in the stratosphere.     

基于混合有限元格式的完美匹配层与多次透射公式人工边界比较研究

介绍了完美匹配层(PML)人工边界可以吸收不同频率和任意角度入射波的原理以及PML人工边界的构造方法. 在此基础上,将PML人工边界应用于地震波动数值模拟的速度应力混合有限元格式中,探讨了PML应用的可行性,并通过数值试验研究了PML人工边界的反射率,比较了PML人工边界与多次透射公式(MTF)人工边界应用于体波和面波模拟中数值反射的差异,对两种边界的透射效果进行了分析. 结果表明, 尽管数值离散后PML人工边界不再保持完美匹配特性,但PML人工边界在近场波动数值模拟中可获得比MTF人工边界更为理想的吸收效果,在角点透射、大角度掠射情形下尤为明显;PML人工边界在混合有限元格式的数值算法中,未见失稳等不良反应,比MTF人工边界有更好的稳定性;在合理选择参数的情况下,PML人工边界的运算量可接受.      

Comparison of perfectly matched layer and multi-transmitting formula artificial boundary condition based on hybrid finite element formulation

The theory of perfectly matched layer (PML) artificial boundary condition (ABC), which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construc- tion process of PML boundary based on elastodynamic partial differential equation (PDE) system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investi- gated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula (MTF) boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in corner and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computa- tional efficiency of PML boundary is only a little lower than MTF boundary.     

Research Note: Overpressure prediction from C‐wave seismic data

Recently, mode converted shear waves (C‐waves) have been shown to enable overpressure prediction in media where primary wave acquisition is inhibited by gas and fluid effects – C‐wave moveout is analysed and a long standing relationship between differential stress and primary‐wave (P‐wave) velocity is modified and employed. Though pore‐pressure prediction based on C‐waves is supported by empirical evidence from laboratory and field experiments, a theoretical justification has yet to be developed. In this research note, we provide a supporting algebra for the original relationship between pore pressure and C‐wave velocity.     

Thermal Generation of Alfvén Waves in Oscillatory Magnetoconvection: Diffusively Modified Modes

Thermal instabilities in the form of oscillatory magnetoconvection representing diffusively modified Alfvén waves in an electrically-conducting Bénard fluid layer of rigid walls in the presence of a vertical magnetic field are investigated. Emphasis of the article is on the transition from a nearly undamped Alfvén wave to diffusively modified Alfvén waves, and on the effect of physically realisable magnetic field boundary conditions on magnetoconvection. It is found that the extra magnetic dissipation in the magnetic Hartmann boundary layers can enhance oscillatory magnetoconvection in the form of strongly modified Alfvén waves. Oscillatory magnetoconvection produced solely by the Alfvén wave mechanism can be the most unstable mode even in the presence of a strong viscous effect. This article also represents the first study on the effect of an electrically conducting wall on magnetoconvection which is associated with a nonlinear eigenvalue problem. We find that the electrically perfectly conducting condition does not yield a good approximation for magnetoconvection with an electrically highly conducting wall. The size of oscillation frequency with an electrically highly conducting wall can be more than a factor of 2 larger than that obtained using the perfectly conducting condition.