黏弹TTI介质中旋转交错网格高阶有限差分数值模拟

以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.     

TTI介质弹性波FCT有限差分数值模拟

TTI(Tilted Transversely Isotropic)各向异性是对地下岩石中广泛存在的规则发育的裂缝和层理的一种有效的弹性近似,基于TTI介质的地震波数值模拟技术是分析地震波在复杂各向异性介质中的传播机理的有效工具.同时,高精度的数值模拟算法也能为后续的逆时偏移技术提供重要的技术支撑.由于TTI介质中地震波方程的弹性参数众多且变化复杂,常规有限差分技术在解决TTI介质正演模拟问题时往往会产生严重的数值频散现象,降低了数值模拟精度.通量校正传输(FluxCorrected Transport,FCT)技术能够有效地压制由空间离散产生的数值频散.本文将FCT技术用于TTI介质中弹性波方程的交错网格高阶精度差分正演,在数值模拟过程中通过对波场进行漫射和反漫射校正实现了空间网格频散的压制.模型模拟结果表明,与常规有限差分算法相比,本文算法能够有效的压制大网格条件下的数值频散,提高模拟精度.     

二维三分量TTI 介质Lebedev 网格正演模拟

鉴于三维各向异性介质(TTI、单斜等)正演模拟在计算量与内存上的巨大消耗以及标准交错网格机制波场插值带来的数值频散,本文采用二维三分量Lebedev交错网格有限差分方法对TTI介质进行波场模拟,利用二维介质便可得到3个相互垂直分量的弹性波场,并利用余弦相似度将其与完全三维正演波场进行对比,分析了该方法的模拟精度。对比测试结果表明,本文方法避免了插值误差,能够精确反映地震波在二维观测平面内的运动学特征,并且平面内质点的偏振速度、振幅能量与三维结果具有较高的相似度,而模拟占用的计算机资源却只相当于三维模拟中的一个二维剖面,是一种高效、准确的各向异性介质数值模拟方法。     

求解3D弹性波方程的并行WNAD方法及其TI介质中的波场模拟

准确模拟TTI介质中弹性波的传播是研究地震各向异性、AVO反演的基础. 在二维加权近似解析离散化(WNAD)算法的基础上, 本文发展的并行WNAD算法是一种研究三维横向各向同性(TI)介质中弹性波传播的、快速高效的数值模拟方法. 我们首先介绍三维WNAD方法的构造过程, 然后与经典的差分格式--交错网格(SG)算法进行了比较. 理论分析和数值算例表明, WNAD算法比交错网格算法更适合在高性能计算机上进行大规模弹性波场模拟. 同时, 本文利用并行的WNAD方法研究了弹性波在TTI介质中的传播规律, 观测了TI介质中弹性波传播的重要特征:横波分离、体波耦合和速度各向异性等. 在TTI介质分界面处, 弹性波产生更加复杂的折射、反射和波型转化, 使得波场非常复杂, 研究和辨别不同类型的波能够加深我们对由裂隙诱导的各向异性介质的认识.     

三维TTI介质中的纯准P波方程及求解方法（英文）

以具有垂直对称轴的横向各向同性(Vertical Transversely Isotropic,简称VTI)介质中纯准P波方程为基础,通过投影变换推导出倾斜横向各向同性(Titled Transversely Isotropic,简称TTI)介质中的纯准P波方程,该方程从根本上消除伪横波的影响。方程形式简单并且在δε的情况下仍能稳定求解。文章还给出TTI介质中纯准P波方程数值求解的高阶有限差分格式并推导对应的稳定性条件,依据完全匹配层(Perfectly Matched Layer,简称PML)的分裂思路给出TTI介质中纯准P波方程的PML吸收边界条件。理论分析和数值计算结果均表明:本文推导的纯准P波方程能更好地描述TTI介质中的准P波,文章给出的正演算法也能够准确模拟TTI介质中准P波的传播过程。     

二维弹性及粘弹性TTI介质中地震波场数值模拟:四种不同网格高阶有限差分算法研究

本文利用交错网格、辅助网格、旋转交错网格、同位网格有限差分方法分别模拟了二维弹性TTI介质和二维黏弹性TTI介质中的地震波传播.在稳定性条件内,选用不同的网格间距及时间间隔,通过波场快照、合成理论地震图较为系统分析对比了这四种不同网格有限差分数值模拟在计算精度、CPU时间、相移、频散、以及保幅方面的优缺点.数值模拟结果表明:1)这四种不同网格有限差分算法都是很好的波场数值模拟算法;2)就CPU计算时间而言,旋转交错网格有限差分算法的计算效率最高;3)从计算精度来看,同位网格有限差分的计算精度最高;4)从振幅保护方面来看,四种网格的保护振幅的能力相当;5)相移方面,当网格间距增大时,交错网格和旋转交错网格有可能出现相移现象;6)频散方面,同位网格的频散现象不明显.     

基于LS-RSGFD方法优化的横向各向同性(TI)介质一阶qP波高精度数值模拟

为克服各向异性介质弹性波数值模拟中存在着计算量大和波场分离困难等局限,研究了声学近似的VTI介质和TTI介质一阶qP波数值模拟方法.首先对VTI介质弹性波方程进行声学近似,推导了VTI介质一阶qP波方程;然后基于精确的TTI介质频散关系,引入一个包含各向异性控制参数σ的新辅助波场,推导了稳定的TTI介质二阶耦合qP波波动方程,并通过引入波场的伪速度分量,推导了等价的一阶应力-速度形式.结合旋转交错网格有限差分(RSGFD)和基于最小二乘优化的有限差分(LS-FD)两种各具优势的方法,研究了最小二乘旋转交错网格有限差分(LS-RSGFD)方法,并用其数值求解VTI和TTI介质一阶qP波方程,然后通过构造其LS-RSGFD格式,实现了高精度的各向异性介质qP波波场数值模拟.数值模拟结果表明:TI介质一阶qP波方程能够准确地模拟各向异性介质中qP波的运动学特征,引入控制参数σ能够有效地减弱不稳定性问题,保证非均匀TTI介质中qP波场的稳定传播;利用优化的LS-RSGFD方法可以得到高精度的合成地震记录,同时还可以相对地提高计算效率.     

稳定的TTI介质拟声波逆时偏移及伪横波的联合压制策略

为克服各向同性和VTI介质逆时偏移方法对复杂地质构造成像的局限,研究了TTI介质拟声波逆时偏移方法.首先从精确的TTI介质频散关系出发,引入一个各向异性控制参数σ,推导了新的二阶耦合TTI介质拟声波方程,以保证波场延拓的稳定性;然后引入波场的伪速度分量,推导了等价的一阶拟声波方程.相比于规则网格有限差分法,交错网格有限差分(SGFD)法能够有效地压制数值频散,模拟精度更高;因此利用高阶SGFD法求解TTI介质一阶拟声波方程,构建逆时偏移所需的正向和逆时波场延拓算子,并应用归一化互相关成像条件实现精确的TTI介质逆时偏移成像.最后,简单讨论了伪横波的产生机制,并给出了伪横波的联合压制策略.模型试验结果验证了方法的有效性和稳定性.     

交错网格高阶差分法三维弹性波数值模拟

本文应用交错网格高阶有限差分方法模拟弹性波在三维各向同性介质中的传播。采用时间上二阶、空间上高阶近似的交错网格高阶差分公式求解三维弹性波位移-应力方程,并在计算边界处应用基于傍轴近似法得到的三维弹性波方程吸收边界条件。在此基础上进行了三维盐丘地质模型的地震波传播数值模拟试算。试算结果表明该方法模拟精度高,在很大程度上减小了数值频散,绕射波更加丰富,而且适用于介质速度具有纵向变化和横向变化的情况。     

二维各向异性介质中地震波场的高阶同位网格有限差分模拟

本文将DRP/opt MacCormack有限差分格式用于模拟二维各向异性介质中的地震波传播.DRP/opt MacCormack是一种同位网格下的差分格式,避免了传统的交错网格在计算各向异性问题时由于变量插值而导致的误差.而且相对于低阶同位网格差分格式,它具有低色散、低耗散的优点.此格式将中心差分算子分成前向和后向两个空间单边差分,然后在4-6步Runge-Kutta时间积分中使用单边差分组合.在具有垂直对称轴的横向各向同性（VTI）模型下,通过对比DRP/opt MacCormack有限差分和谱元方法的模拟结果,验证了前者具有很高的精度和稳定性.由于实际地质条件下TI介质的对称轴通常是倾斜的(TTI),本文在二维三分量框架下模拟TTI介质中的地震波场.结果显示横波分裂和切平面/反平面运动耦合的特征.数值实验表明DRP/opt MacCormack是一种有效的研究各向异性介质中地震波传播规律的差分格式.     

A hybrid absorbing boundary condition for elastic staggered‐grid modelling

We recently proposed an efficient hybrid scheme to absorb boundary reflections for acoustic wave modelling that could attain nearly perfect absorptions. This scheme uses weighted averaging of wavefields in a transition area, between the inner area and the model boundaries. In this paper we report on the extension of this scheme to 2D elastic wave modelling with displacement‐stress formulations on staggered grids using explicit finite‐difference, pseudo‐implicit finite‐difference and pseudo‐spectral methods. Numerical modelling results of elastic wave equations with hybrid absorbing boundary conditions show great improvement for modelling stability and significant absorption for boundary reflections, compared with the conventional Higdon absorbing boundary conditions, demonstrating the effectiveness of this scheme for elastic wave modelling. The modelling results also show that the hybrid scheme works well in 2D rotated staggered‐grid modelling for isotropic medium, 2D staggered‐grid modelling for vertically transversely isotropic medium and 2D rotated staggered‐grid modelling for tilted transversely isotropic medium.     

Numerical simulation of seismic wave in elastic and viscoelastic TTI media*

Based on the two-dimensional (2D) three-component first-order velocity-stress equation, the high order staggered mesh finite difference numerical simulation method was used to simulate the elastic and viscoelastic tilted transversely isotropic (TTI) media. The perfect matched layer (PML) absorption boundary condition was selected to eliminate the boundary effect. The results show that: ① Under the condition of fixed elastic parameters of elastic TTI medium, when the polarization angle and azimuth are 60° and 45° respectively, the degree of shear wave splitting is significantly greater than the angle of 0°; ② The influence of viscoelasticity on TTI medium is mainly reflected in the amplitude. If the quality factor decreases, the attenuation of the seismic wave amplitude increases, causing the waveform to become wider and distorted. If the quality factor increases, the viscoelastic medium becomes closer to elastic medium; ③ For TTI medium with different polarization angle and azimuth angle in the upper and lower layers, the shear wave can multiple splits at the interface of medium. The symmetry of seismograms is affected by the polarization angle and azimuth angle of TTI medium; ④ Viscoelasticity has a great influence on reflected wave, transmitted wave and converted wave in the low-velocity model. When the viscoelasticity is strong, the weaker waves may not be shown.     

Anisotropic reverse-time migration for tilted TI media

Seismic anisotropy in dipping shales results in imaging and positioning problems for underlying structures. We develop an anisotropic reverse‐time depth migration approach for P‐wave and SV‐wave seismic data in transversely isotropic (TI) media with a tilted axis of symmetry normal to bedding. Based on an accurate phase velocity formula and dispersion relationships for weak anisotropy, we derive the wave equation for P‐wave and SV‐wave propagation in tilted transversely isotropic (TTI) media. The accuracy of the P‐wave equation and the SV‐wave equation is analyzed and compared with other acoustic wave equations for TTI media. Using this analysis and the pseudo‐spectral method, we apply reverse‐time migration to numerical and physical‐model data. According to the comparison between the isotropic and anisotropic migration results, the anisotropic reverse‐time depth migration offers significant improvements in positioning and reflector continuity over those obtained using isotropic algorithms.     

Pseudo‐spectral method using rotated staggered grid for elastic wave propagation in 3D arbitrary anisotropic media

Staggering grid is a very effective way to reduce the Nyquist errors and to suppress the non‐causal ringing artefacts in the pseudo‐spectral solution of first‐order elastic wave equations. However, the straightforward use of a staggered‐grid pseudo‐spectral method is problematic for simulating wave propagation when the anisotropy level is greater than orthorhombic or when the anisotropic symmetries are not aligned with the computational grids. Inspired by the idea of rotated staggered‐grid finite‐difference method, we propose a modified pseudo‐spectral method for wave propagation in arbitrary anisotropic media. Compared with an existing remedy of staggered‐grid pseudo‐spectral method based on stiffness matrix decomposition and a possible alternative using the Lebedev grids, the rotated staggered‐grid‐based pseudo‐spectral method possesses the best balance between the mitigation of artefacts and efficiency. A 2D example on a transversely isotropic model with tilted symmetry axis verifies its effectiveness to suppress the ringing artefacts. Two 3D examples of increasing anisotropy levels demonstrate that the rotated staggered‐grid‐based pseudo‐spectral method can successfully simulate complex wavefields in such anisotropic formations.     

基于横向各向同性BISQ方程的弹性波传播数值模拟

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制. 近年来，利用同时处理这两种力学机制的BISQ(Biot-Squirt)模型，弹性波衰减和频散的问题已被广泛研究；然而基于BISQ方程的波场数值模拟尚未见到公开的报道．本文从BISQ方程出发，利用交错网格方法对横向各向同性孔隙介质中不同频率和相界情况，以及双层介质中的弹性波传播进行数值模拟，研究了在同时考虑两种流动机制作用情况下地震波和声波的传播特性及传播过程中出现的各种波动现象．      

带正则化校正的TTI介质qP波方程及其逆时偏移方法和应用

各向异性研究对地下介质精确成像有着重要的意义,在当前计算机硬件迅速发展及宽方位地震数据采集日益普遍的情况下,成像必须考虑介质的各向异性.逆时偏移是基于双程波动方程的较为精确的数值解的成像方法,所以相对于其他地震成像方法,它具有很大的优势,譬如不受反射界面的倾角限制、偏移速度结构合适时能够使回转波及多次波正确成像.在各向同性介质中,可使用标量波方程来模拟波场.而在各向异性介质中,P波和SV波是相互耦合的,即不存在单纯的标量波传播,通常利用能代表耦合波场中P波分量运动学特征的拟声波(qP波)进行偏移成像.本文中,我们推导出了TTI介质下qP波控制方程.该方程可采用显式有限差分格式进行求解.通过声学近似,若沿对称轴方向的剪切波速度为零,对于对称轴方向不变且ε≥δ的模型来说,可得到稳定的数值解.但对于TTI介质来说,由于沿对称轴方向各向异性参数是变化的,声学近似会引起波场传播及数值计算的不稳定.因此,我们提出了正则化有限横波的方法,很好地解决了这一问题.最后,给出了Foothill模型的测试结果及某探区实际资料试算结果,展示了采用这个方程进行复杂TTI模型正演和高质量逆时偏移成像结果,证实了该方法的正确性和实际资料应用中的有效性.     

Numerical simulation of elastic wave propagation based on the transversely isotropic BISQ equation

Introduction More real models are being developed by the modern seismology. As we all know, the earth is not a simple elastic body. Oil and gas reservoir, ground surface, seashore zone, sea bottom layer, etc, are porous solid media with fluids. It has been confirmed that there are two main fluid flow mechanisms in these media (Dvorkin, Nur, 1993), i.e., the Biot flow mechanism (Biot, 1956, 1962) based on the macroscopic property and the Squirt-flow mechanism (Mavko, Nur, 1979) based on the …     

Pure acoustic wave propagation in transversely isotropic media by the pseudospectral method

Seismic wave propagation in transversely isotropic (TI) media is commonly described by a set of coupled partial differential equations, derived from the acoustic approximation. These equations produce pure P‐wave responses in elliptically anisotropic media but generate undesired shear‐wave components for more general TI anisotropy. Furthermore, these equations suffer from instabilities when the anisotropy parameter ε is less than δ. One solution to both problems is to use pure acoustic anisotropic wave equations, which can produce pure P‐waves without any shear‐wave contaminations in both elliptical and anelliptical TI media. In this paper, we propose a new pure acoustic transversely isotropic wave equation, which can be conveniently solved using the pseudospectral method. Like most other pure acoustic anisotropic wave equations, our equation involves complicated pseudo‐differential operators in space which are difficult to handle using the finite difference method. The advantage of our equation is that all of its model parameters are separable from the spatial differential and pseudo‐differential operators; therefore, the pseudospectral method can be directly applied. We use phase velocity analysis to show that our equation, expressed in a summation form, can be properly truncated to achieve the desired accuracy according to anisotropy strength. This flexibility allows us to save computational time by choosing the right number of summation terms for a given model. We use numerical examples to demonstrate that this new pure acoustic wave equation can produce highly accurate results, completely free from shear‐wave artefacts. This equation can be straightforwardly generalized to tilted TI media.