S-wave in 2D acoustic transversely isotropic media with a tilted symmetry axis

In an acoustic transversely isotropic medium, there are two waves that propagate. One is the P-wave and another one is the S-wave (also known as S-wave artefact). This paper is devoted to analyse the S-wave in two-dimensional acoustic transversely isotropic media with a tilted symmetry axis. We derive the S-wave slowness surface and traveltime function in a homogeneous acoustic transversely isotropic medium with a tilted symmetry axis. The S-wave traveltime approximations in acoustic transversely isotropic media with a tilted symmetry axis can be mapped from the counterparts for acoustic transversely isotropic media with a vertical symmetry axis. We consider a layered two-dimensional acoustic transversely isotropic medium with a tilted symmetry axis to analyse the S-wave moveout. We also illustrate the behaviour of the moveout for reflected S-wave and converted waves.     

二维横各向同性弹性随机介质中的波场特征

本文通过交错网格有限差分正演．模拟了平面地震波在二维横各向同性弹性随机介质模型中的传播及其自激自收时间记录．为研究横各向同性弹性随机介质模型中的波场特征,我们在五个不同的时间区段上,分别计算剖面的三个统计特征(横向中心频率、纵向中心频率、波场能量相对值)．这样,对应每一个横各向同性弹性随机介质模型．均可计算得到15个不同的波场特征量．我们通过在二维横各向同性弹性随机介质中的正演模拟．研究当自相关长度以及介质的各向异性系数变化时,对应的上述波场特征量的变化特点．证实了在随机介质模型中．各向异性系数的变化会引起波场记录上的某些统计特征的变化,归纳得出了若干结论．     

考虑各向异性的层土-盾构隧道地震反应数值模拟

在层状各向异性土体-盾构隧道地震反应分析中,引入了横观各向同性弹塑性模型理论,建立了适合于横观各向同性介质的双渐近-多向透射边界条件。针对地铁区间盾构隧道抗震设计的特点,基于横观各向同性弹塑性模型,研制了考虑层状土体各向异性和施工开挖效应,适合于盾构隧道动力计算的各向异性弹塑性动力有限元程序。在程序中对于不同的材料采用了不同的本构关系和不同的单元形式,并采用了关联流动法则和多种屈服准则,可同时进行各向异性土体与地下结构的二维平面应力、平面应变和轴对称问题的静力和动力数值分析。最后利用所研制的程序进行了上海地铁二号线石门一路站附近区间隧道在不同超越概率地震动输入下的隧道反应计算。结果表明,在层状土体-地铁区间隧道的抗震分析中考虑土体各向异性的影响是必要的,所提出的计算模型是可行的。     

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.     

An efficient wave extrapolation method for anisotropic media with tilt

Wavefield extrapolation operators for elliptically anisotropic media offer significant cost reduction compared with that for the transversely isotropic case, particularly when the axis of symmetry exhibits tilt (from the vertical). However, elliptical anisotropy does not provide accurate wavefield representation or imaging for transversely isotropic media. Therefore, we propose effective elliptically anisotropic models that correctly capture the kinematic behaviour of wavefields for transversely isotropic media. Specifically, we compute source‐dependent effective velocities for the elliptic medium using kinematic high‐frequency representation of the transversely isotropic wavefield. The effective model allows us to use cheaper elliptic wave extrapolation operators. Despite the fact that the effective models are obtained by matching kinematics using high‐frequency asymptotic, the resulting wavefield contains most of the critical wavefield components, including frequency dependency and caustics, if present, with reasonable accuracy. The methodology developed here offers a much better cost versus accuracy trade‐off for wavefield computations in transversely isotropic media, particularly for media of low to moderate complexity. In addition, the wavefield solution is free from shear‐wave artefacts as opposed to the conventional finite‐difference‐based transversely isotropic wave extrapolation scheme. We demonstrate these assertions through numerical tests on synthetic tilted transversely isotropic models.     

S‐wave kinematics in acoustic transversely isotropic media with a vertical symmetry axis

Acoustic transversely isotropic models are widely used in seismic exploration for P‐wave processing and analysis. In isotropic acoustic media only P‐wave can propagate, while in an acoustic transversely isotropic medium both P and S waves propagate. In this paper, we focus on kinematic properties of S‐wave in acoustic transversely isotropic media. We define new parameters better suited for S‐wave kinematics analysis. We also establish the travel time and relative geometrical spreading equations and analyse their properties. To illustrate the behaviour of the S‐wave in multi‐layered acoustic transversely isotropic media, we define the Dix‐type equations that are different from the ones widely used for the P‐wave propagation.     

Depth migration of shot records in heterogeneous, transversely isotropic media using optimum explicit operators

A space–frequency domain 2D depth-migration scheme is generalized for imaging in the presence of anisotropy. The anisotropy model used is that of a transversely isotropic (TI) medium with a symmetry axis that can be either vertical or tilted. In the proposed scheme the anisotropy is described in terms of Thomsen parameters; however, the scheme can accommodate a wide range of anisotropy rather than only weak anisotropy. Short spatial convolution operators are used to extrapolate the wavefields recursively in the space–frequency domain for both qP- and qSV-waves. The weighted least-squares method for designing isotropic optimum operators is extended to asymmetric optimum explicit extrapolation operators in the presence of TI media with a tilted symmetry axis. Additionally, an efficient weighted quadratic-programming design method is developed. The short spatial length of the derived operators makes it possible for the proposed scheme to handle lateral inhomogeneities. The performance of the operators, designed by combining the weighted least-squares and weighted quadratic-programming methods, is demonstrated by migration impulse responses of qP and qSV propagation modes for the weak and strong TI models with both vertical and tilted symmetry axes. Finally, a table-driven shot-record depth-migration scheme is proposed, which is illustrated for finite-difference modelled shot records in TI media.     

横向各向同性介质紧致交错网格有限差分波场模拟(英文)

针对有限差分数值模拟的频散问题,本文将交错网格技术和紧致差分格式相结合,推导了横向各向同性介质一阶速度一应力波动方程的紧致交错网格差分格式;对比分析了紧致交错网格差分格式、交错网格差分格式以及紧致差分格式的截断误差主项,并利用Fourier误差分析方法分析了上述三种差分格式的近似精度;在此基础上,分别采用上述三种差分格式进行了波场数值模拟。结果表明,当差分方程阶数相同时,紧致交错网格差分格式截断误差最小,数值频散最弱,差分精度最高,证实了该方法的有效性。     

Seismic wavefield modelling in two-phase media including undulating topography with the modified Biot/squirt model by a curvilinear-grid finite difference method

In this paper, we deduced the corresponding first-order velocity–stress equation for curvilinear coordinates from the first-order velocity–stress equation based on the modified Biot/squirt model for a two-dimensional two-phase medium. The equations are then numerically solved by an optimized high-order non-staggered finite difference scheme, that is, the dispersion relation preserving/optimization MacCormack scheme. To implement undulating free-surface topography, we derive an analytical relationship between the derivatives of the particle velocity components and use the compact finite-difference scheme plus a traction-image method. In the undulating free surface and the undulating subsurface interface of two-phase medium, the complex reflected wave and transmitted wave can be clearly recognized in the numerical simulation results. The simulation results show that the curvilinear-grid finite-difference method, which uses a body-conforming grid to describe the undulating surface, can accurately reduce the numerical scattering effect of seismic wave propagation caused by the use of ladder-shaped grid to fit the surfaces when undulating topography is present in a two-phase isotropic medium.     

Moveout approximation for horizontal transversely isotropic and vertical transversely isotropic layered medium. Part II: effective model‡

We use residual moveouts measured along continuous full azimuth reflection angle gathers, in order to obtain effective horizontal transversely isotropic model parameters. The angle gathers are generated through a special angle domain imaging system, for a wide range of reflection angles and full range of phase velocity azimuths. The estimation of the effective model parameters is performed in two stages. First, the background horizontal transversely isotropic (HTI)/vertical transversely isotropic (VTI) layered model is used, along with the values of reflection angles, for converting the measured residual moveouts (or traveltime errors) into azimuthally dependent normal moveout (NMO) velocities. Then we apply a digital Fourier transform to convert the NMO velocities into azimuthal wavenumber domain, in order to obtain the effective HTI model parameters: vertical time, vertical compression velocity, Thomsen parameter delta and the azimuth of the medium axis of symmetry. The method also provides a reliability criterion of the HTI assumption. The criterion shows whether the medium possesses the HTI type of symmetry, or whether the azimuthal dependence of the residual traveltime indicates to a more complex azimuthal anisotropy. The effective model used in this approach is defined for a 1D structure with a set of HTI, VTI and isotropic layers (with at least one HTI layer). We describe and analyse the reduction of a multi‐layer structure into an equivalent effective HTI model. The equivalent model yields the same NMO velocity and the same offset azimuth on the Earth's surface as the original layered structure, for any azimuth of the phase velocity. The effective model approximates the kinematics of an HTI/VTI layered structure using only a few parameters. Under the hyperbolic approximation, the proposed effective model is exact.     

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

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

横向各向同性介质中地震波场谱元法数值模拟

横向各向同性介质是地球内部广泛存在的一种各向异性介质,因此为了能够更好地认识地震波在这种介质中的传播特征,用数值方法进行地震波模拟显得十分必要.本文采用谱元法对横向同性介质中的地震波进行模拟,该方法基于弹性力学方程弱形式基础之上,具有有限元适应任意复杂介质模型的韧性和伪谱法的精度.文中阐述了基于Legendre多项式的谱元法的理论和推导过程,该方法可以形成全局对角质量矩阵,在时间域使用显式的差分算法,提高运算效率,最后通过横向各向同性介质的数值计算,模拟结果表明该方法是一种有效的数值模拟方法.     

Determination of the reference symmetry axis of a generally anisotropic medium which is approximately transversely isotropic

For a given stiffness tensor (tensor of elastic moduli) of a generally anisotropic medium, we estimate to what extent the medium is transversely isotropic (uniaxial) and determine the direction of its reference symmetry axis expressed in terms of the unit reference symmetry vector. If the medium is exactly transversely isotropic (exactly uniaxial), we obtain the direction of its symmetry axis. We can also calculate the first–order and second–order spatial derivatives of the reference symmetry vector which may be useful in tracing the reference rays for the coupling ray theory. The proposed method is tested using various transversely isotropic (uniaxial) and approximately transversely isotropic (approximately uniaxial) media.     

双相各向同性介质弹性波高精度波场分离数值模拟方法

为了便于研究双相介质固流相混合弹性波场中纵横波波场的传播规律,提出了基于交错网格的Biot双相各向同性介质弹性波动方程高精度波场分离正演数值模拟方法.采用高阶交错网格有限差分法来构建一阶双曲型双相各向同性介质弹性波动方程正演算子实现波场正演,并在每一步递推过程中,分别计算出同相和流相分量相应的散度场(纯纵波场)和旋度场...     

Triplications for the converted wave in transversely isotropic media with a tilted symmetry axis

In seismic data processing, serious problems could be caused by the existence of triplication and need to be treated properly for tomography and other inversion methods. The triplication in transversely isotropic medium with a vertical symmetry axis has been well studied and concluded that the triplicated traveltime only occurs for S wave and there is no triplication for P and converted PS waves since the P wave convexity slowness always compensates the S wave slowness concavity. Compared with the vertical symmetry axis model, the research of the triplication in transversely isotropic medium with a tilted symmetry axis is still keeping blank. In order to analyse the triplication for the converted wave in the tilted symmetry axis model, we examine the traveltime of the triplication from the curvature of averaged P and S wave slowness. Three models are defined and tested in the numerical examples to illustrate the behaviour of the tilted symmetry axis model for the triplicated traveltime with the change of the rotation angle. Since the orientation of an interface is related to the orientation of the symmetry axis, the triplicated traveltime is encountered for the converted wave in the tilted symmetry axis model assuming interfaces to be planar and horizontal. The triplicated region is influenced by the place and level of the concave curvature of the P and S wave slowness.     

ELASTIC WAVE PROPAGATION IN TRANSVERSELY ISOTROPIC MEDIA USING FINITE DIFFERENCES1

When treating the forward full waveform case, a fast and accurate algorithm for modelling seismic wave propagation in anisotropic inhomogeneous media is of considerable value in current exploration seismology. Synthetic seismograms were computed for P-SV wave propagation in transversely isotropic media. Among the various techniques available for seismic modelling, the finite-difference method possesses both the power and flexibility to model wave propagation accurately in anisotropic inhomogeneous media bounded by irregular interfaces. We have developed a fast high-order vectorized finite-difference algorithm adapted for the vector supercomputer. The algorithm is based on the fourth-order accurate MacCormack-type splitting scheme. Solving the equivalent first-order hyperbolic system of equations, instead of the second-order wave equation, avoids computation of the spatial derivatives of the medium's anisotropic elastic parameters. Examples indicate that anisotropy plays an important role in modelling the kinematic and the dynamic properties of the wave propagation and should be taken into account when necessary.     

Finite-difference elastic wave propagation in 2D heterogeneous transversely isotropic media1

The velocity-stress formulation for propagation of elastic seismic waves through 2D heterogeneous transversely isotropic media of arbitrary orientation is presented. The equations are recast into a finite-difference scheme and solved numerically using fourth-order spatial operators and a second-order temporal operator on a staggered grid. Absorbing, free-surface and symmetry boundary conditions have been implemented. Test cases compare well with other published solutions. Synthetic seismograms are calculated over two idealized models: (i) vertical fractures in granite with a dolerite sill reflector and (ii) a dipping anisotropic shale. Comparisons with the isotropic counterparts show significant differences which may have to be accounted for in seismic processing in the future.     

起伏地表下基于改进BISQ模型双相介质中曲线网格有限差分法波场模拟

本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础,推导出了曲线坐标系下对应的方程,然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程,并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件,实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面,地表的网格线紧贴地形,避免了台阶近似造成的数值散射.数值模拟结果表明,在双相介质起伏自由地表和分界面处,各类波型复杂的反射透射规律可以清晰展现,曲线网格有限差分法可以精确地解决地震波在含起伏地表的双相各向同性介质中的传播问题.     

On possible fluid detection in equivalent transversely isotropic media

We consider a transversely isotropic medium that is long-wave equivalent to a stack of thin, parallel, isotropic layers and is obtained using the Backus average. In such media, we analyse the relations among anisotropy parameters; Thomsen parameters, ε and δ, and a new parameter ϕ. We discuss the last parameter and show its essential properties; it equals 0 in the case of isotropy of equivalent medium and/or constant Lamé coefficient λ in layers. The second property occurs to make ϕ sensitive to variations of λ in thin-bedded sequences. According to Gassmann, in isotropic media the variation of fluid content affects only the Lamé coefficient λ, not μ; thus, the sensitivity to changes of λ is an essential property in the context of possible detection of fluids. We show algebraically and numerically that ϕ is more sensitive to these variations than ε or δ. Nevertheless, each of these parameters is dependent on the changes of μ; to understand this influence, we exhibit comprehensive tables that illustrate the behaviour of anisotropy parameters with respect to specific variations of λ and μ. The changes of μ in layers can be presented by the Thomsen parameter γ that depends on them solely. Hence, knowing the values of elasticity coefficients of equivalent transversely isotropic medium, we may compute ϕ and γ, and based on the aforementioned tables, we predict the expected variation of λ; in this way, we propose a new method of possible fluid detection. Also, we show that the prior approach of possible detection of fluids, proposed by Berryman et al., may be unreliable in specific cases. To establish our results, we use the Monte Carlo method; for the range and chosen variations of Lamé coefficients λ and μ – relevant to sandstones – we generate these coefficients in thin layers and, after the averaging process, we obtain an equivalent transversely isotropic medium. We repeat that process numerous times to get many equivalent transversely isotropic media, and – for each of them  – we compute their anisotropy parameters. We illustrate ϕ, ε and δ in the form of cross-plots that are relevant to the chosen variations of λ and μ. Additionally, we present a table with the computed ranges of anisotropy parameters that correspond to different variations of Lamé coefficients.     

Transmitting boundary for water-saturated transversely isotropic strata based on the u–U formulation

A new transmitting boundary in a cylindrical coordinate system has been developed for modeling the elastic waves radiating out to an infinite boundary in water-saturated transversely isotropic soil strata over a rigid bedrock. The saturated soil strata are assumed to consist of a porous material and modeled as a transversely isotropic two-phase medium, based on the u−U formulation. The newly developed transmitting boundary is combined with the finite-elements model of the near-field region, using the same u−U formulation, and applied to the study of the dynamics of a rigid circular foundation in porous isotropic or transversely isotropic layered strata, either fully or partly saturated with water. The verification and application examples give valuable insights into new and interesting aspects of the dynamic behavior of rigid circular foundations in fully or partly saturated two-phase ground in terms of permeability, transverse anisotropy, and ground-water table level.