横向各向同性介质和斜方介质各向异性参数的约束条件

Thomsen提出的横向各向同性(TI)介质各向异性参数(α0、β0、ε、δ、γ)是各向异性理论研究和实际资料处理中的常用参数,Thomsen参数的取值必须符合物理学定律和实际地学情况,随意的取值可能导致无意义乃至错误的结果.本文根据热力学定律和弹性常数的物理意义,结合大量的实测数据,提出常见TI介质的Thomsen参数需要满足以下约束条件:1/4f=1-β02/α02ε>-f/2; 1/2f-1δf-1); -1/2γε)/4(1-f)-1/2.Thomsen参数ε、δ、γ的取值区间主要受P、S波参考速度比β0/α0的约束,其值可正可负,实测参数中ε和γ正多负少;δ、γ、β0/α0的取值范围都有上下界,而ε没有上限;γ与ε正相关且取值上限受到ε的限制;ε、γ与δ之间不存在明显的约束.鉴于椭圆各向异性介质和薄互层等效TI介质在实际应用中的普遍性,专门给出了这两种介质各向异性参数的附加约束条件.一些在实际资料中可观测到的P波特殊偏振方向和SV波三叉现象也能为各向异性参数提供额外的约束.之后,我们将TI介质弹性常数和各向异性参数的约束条件扩展到对正交各向异性介质弹性常数和各向异性参数的约束.本文提出的各向异性参数约束条件简单实用,为各向异性理论研究和数值模拟中参数的选择提供依据,为实际资料的各向异性参数反演提供约束,既能避免出现无物理意义的研究结果,又能加速反演的搜索过程,提高生产效率,具有理论指导意义和实际应用价值.     

On shear-wave triplication in transversely isotropic media

The exact solution to the problem of qSV triplication in homogenous transversely isotropic media has been long known, but the result is algebraically complex and is seldom applied in practice. We present an appropriate approximation (not assuming weak qSV-anisotropy) that simplifies the conditions for the onset of off-axis triplication as anisotropy is increased, identifying the anisotropy parameter σ as the controlling parameter. It follows that commonly reported surface-seismic P-wave move-out measurements imply that many formations in the earth's sedimentary crust support off-axis qSV triplications. For typical Vp/Vs velocity ratios and a horizontally stratified earth, however, off-axis qSV triplications appear to only occur for shear-wave incidence angles too far from the vertical to be sampled by surface-seismic converted-wave survey geometries.     

Pure P- and S-wave equations in transversely isotropic media

Pure-mode wave propagation is important for applications ranging from imaging to avoiding parameter tradeoff in waveform inversion. Although seismic anisotropy is an elastic phenomenon, pseudo-acoustic approximations are routinely used to avoid the high computational cost and difficulty in decoupling wave modes to obtain interpretable seismic images. However, such approximations may result in inaccuracies in characterizing anisotropic wave propagation. We propose new pure-mode equations for P- and S-waves resulting in an artefact-free solution in transversely isotropic medium with a vertical symmetry axis. Our approximations are more accurate than other known approximations as they are not based on weak anisotropy assumptions. Therefore, the S-wave approximation can reproduce the group velocity triplications in strongly anisotropic media. The proposed approximations can be used for accurate modelling and imaging of pure P- and S-waves in transversely isotropic media.     

The impact of reservoir scale on amplitude variation with offset

To investigate the vertical and horizontal impact of reservoir scale on the amplitude‐versus‐offset characteristics, we conduct seismic numerical simulations on models containing spatially confined lithologic units with different scales. We find that the reservoir scale has a nonlinear effect on the amplitude‐versus‐offset intercepts and gradients. As the reservoir width increases, amplitude‐versus‐offset intercept and gradient both first increase, then decrease, and finally remain stable. The amplitude‐versus‐offset intercept is maximum when the reservoir width is 80% larger than the Fresnel zone radius, whereas the amplitude‐versus‐offset gradient peaks at 1.5 times the Fresnel zone radius. Both amplitude‐versus‐offset intercept and gradient are approximately proportional to the reservoir width prior to reaching their maxima. When the lateral extent of the reservoir is more than three times the Fresnel zone radius, the amplitude‐versus‐offset attributes are constant. Modelling the reservoir thickness shows that intercept and gradient behave in a manner similar to that of tuning of thin beds. Both the amplitude‐versus‐offset intercept and gradient first increase and then decrease with the thickness, peaking at the tuning thickness. The thickness contribution to amplitude‐versus‐offset variations is negligible when the thickness is larger than 1.6 times of the tuning thickness. Considering the magnitude of the changes in amplitude‐versus‐offset intercept and gradient caused by reservoir scale, the width causes a maximum 433% intercept increase and a 344% gradient increase, whereas the thickness causes a maximum 100% intercept increase and a 73% gradient increase. Cross‐plotting the amplitude‐versus‐offset intercept and gradient shows the reservoir scale change gives rise to an anti‐clockwise spiraling effect. In conclusion, the lateral and vertical extents of the reservoir both play an important role in amplitude variation with offset. Our analysis shows that the lateral reservoir extent has a larger impact on the amplitude variation with offset than the vertical tuning effect.     

Three-parameter Radon transform in layered transversely isotropic media

Radon transform is a powerful tool with many applications in different stages of seismic data processing, because of its capability to focus seismic events in the transform domain. Three-parameter Radon transform can optimally focus and separate different seismic events, if its basis functions accurately match the events. In anisotropic media, the conventional hyperbolic or shifted hyperbolic basis functions lose their accuracy and cannot preserve data fidelity, especially at large offsets. To address this issue, we propose an accurate traveltime approximation for transversely isotropic media with vertical symmetry axis, and derive two versions of Radon basis functions, time-variant and time-invariant. A time-variant basis function can be used in time domain Radon transform algorithms while a time-invariant version can be used in, generally more efficient, frequency domain algorithms. Comparing the time-variant and time-invariant Radon transform by the proposed basis functions, the time-invariant version can better focus different seismic events; it is also more accurate, especially in presence of vertical heterogeneity. However, the proposed time-invariant basis functions are suitable for a specific type of layered anisotropic media, known as factorized media. We test the proposed methods and illustrate successful applications of them for trace interpolation and coherent noise attenuation.     

P-wave complex-valued traveltimes in homogeneous attenuating transversely isotropic media

Computation of complex-valued traveltimes provides an efficient approach to describe the seismic wave attenuation for applications like attenuation tomography, inverse Q filtering and Kirchhoff migration with absorption compensation. Attenuating acoustic transverse isotropy can be used to describe the directional variation of velocity and attenuation of P-waves in thin-bedding geological structures. We present an approximate method to solve the acoustic eikonal equation for an attenuating transversely isotropic medium with a vertical symmetry axis. We take into account two similar parameterizations of an attenuating vertical symmetry axis medium. The first parameterization uses the normal moveout velocity, whereas the second parameterization uses the horizontal velocity. For each parameterization, we combine perturbation theory and the Shanks transform in different ways to derive analytic solutions. Numerical examples show that the analytic solutions derived from the second parameterization yield better accuracy. The Shanks transform solution with respect to only the anellipticity parameter from the second parameterization is demonstrated numerically to be the most accurate among all the analytic solutions.     

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.     

多层垂直对称轴横向各向同性介质精确走时计算

给出了计算多层垂直对称轴横向各向同性（VTI）介质精确射线路径和走时的方法,所用的体波相速度公式、群速度公式和Snell定律都是严格的显式解析公式. 任意基本波的射线路径和走时计算问题都可以转化成一个等效的透射问题,再用文中的公式来计算,具体实现方法用一个多次波和一个首波的实例给出. 最后分别用精确公式和Thomsen近似公式计算了相同模型相同基本波的走时曲线. 比较两者计算结果可发现, 近似公式反复使用会使误差积累,同时揭示了近似公式适用范围的局限性,强调了使用近似公式需要注意其适用范围的重要性.      

Azimuth moveout correction for transversely isotropic media

The azimuth moveout (AMO) operator in homogeneous transversely isotropic media with a vertical symmetry axis (VTI), as in isotropic media, has an overall skewed saddle shape. However, the AMO operator in anisotropic media is complicated; it includes, among other things, triplications at low angles. Even in weaker anisotropies, with the anisotropy parameter η= 0.1 (10% anisotropy), the AMO operator is considerably different from the isotropic operator, although free of triplications. The structure of the operator in VTI media (positive η) is stretched (has a wider aperture) compared with operators in isotropic media, with the amount of stretch being dependent on the strength of anisotropy. If the medium is both vertically inhomogeneous, i.e. the vertical velocity is a function of depth (v(z)), and anisotropic, which is a common combination in practical problems, the shape of the operator again differs from that for isotropic media. However, the difference in the AMO operator between the homogeneous and the v(z) cases, even for anisotropic media, is small. Stated simply, anisotropy influences the shape and aperture of the AMO operator far more than vertical inhomogeneity does.     

Angle gathers in wave‐equation imaging for transversely isotropic media

In recent years, wave‐equation imaged data are often presented in common‐image angle‐domain gathers as a decomposition in the scattering angle at the reflector, which provide a natural access to analysing migration velocities and amplitudes. In the case of anisotropic media, the importance of angle gathers is enhanced by the need to properly estimate multiple anisotropic parameters for a proper representation of the medium. We extract angle gathers for each downward‐continuation step from converting offset‐frequency planes into angle‐frequency planes simultaneously with applying the imaging condition in a transversely isotropic with a vertical symmetry axis (VTI) medium. The analytic equations, though cumbersome, are exact within the framework of the acoustic approximation. They are also easily programmable and show that angle gather mapping in the case of anisotropic media differs from its isotropic counterpart, with the difference depending mainly on the strength of anisotropy. Synthetic examples demonstrate the importance of including anisotropy in the angle gather generation as mapping of the energy is negatively altered otherwise. In the case of a titled axis of symmetry (TTI), the same VTI formulation is applicable but requires a rotation of the wavenumbers.     

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.     

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.     

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.     

Exact traveltime computation in multi-layered transversely isotropic media with vertical symmetry axis

An approach to calculate the accurate ray paths and traveltimes in multi-layered VTI media (transversely isotropic media with a vertical symmetry axis) is proposed. The expressions of phase velocity, group velocity and Snell’s law used for computation are all explicit and exact. The calculation of ray paths and traveltimes for a given ele-mentary wave is equivalent to that of a transmission problem which is much easier to be treated with the formulae proposed. In the section of numerical examples, the proce...     

Linearized amplitude variation with offset and azimuth and anisotropic poroelasticity

The normal-to-shear weakness ratio is commonly used as a fracture fluid indicator, but it depends not only on the fluid types but also on the fracture intensity and internal architecture. Amplitude variation with offset and azimuth is commonly used to perform the fluid identification and fracture characterization in fractured porous rocks. We demonstrate a direct inversion approach to utilize the observable azimuthal data to estimate the decoupled fluid (fluid/porosity term) and fracture (normal and shear weaknesses) parameters instead of the calculation of normal-to-shear weakness ratio to help reduce the uncertainties in fracture characterization and fluid identification of a gas-saturated porous medium permeated by a single set of parallel vertical fractures. Based on the anisotropic poroelasticity and perturbation theory, we first derive a linearized amplitude versus offset and azimuth approximation using the scattering function to decouple the fluid indicator and fracture parameters. Incorporating Bayes formula and convolution theory, we propose a feasible direct inversion approach in a Bayesian framework to obtain the direct estimations of model parameters, in which Cauchy and Gaussian distribution are used for the a priori information of model parameters and the likelihood function, respectively. We finally use the non-linear iteratively reweighted least squares to solve the maximum a posteriori solutions of model parameters. The synthetic examples containing a moderate noise demonstrate the feasibility of the proposed approach, and the real data illustrates the stabilities of estimated fluid indicator and dry fracture parameters in gas-saturated fractured porous rocks.     

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.     

Multi-parameter reflection waveform inversion for acoustic transversely isotropic media with a vertical symmetry axis

Full waveform inversion in transversely isotropic media with a vertical symmetry axis provides an opportunity to better match the data at the near and far offsets. However, multi-parameter full waveform inversion, in general, suffers from serious cycle-skipping and trade-off problems. Reflection waveform inversion can help us recover a background model by projecting the residuals of the reflected wavefield along the reflection wavepath. Thus, we extend reflection waveform inversion to acoustic transversely isotropic media with a vertical symmetry axis utilizing the proper parameterization for reduced parameter trade-off. From a radiation patterns analysis, an acoustic transversely isotropic media with a vertical symmetry axis is better described by a combination of the normal-moveout velocity and the anisotropic parameters η and δ for reflection waveform inversion applications. We design a three-stage inversion strategy to construct the optimal resulting model. In the first stage, we only invert for the background by matching the simulated reflected wavefield from the perturbations of and δ with the observed reflected wavefield. In the second stage, the background and η are optimized simultaneously and the far-offset reflected wavefield mainly contribute to their updates. We perform Born modelling to compute the reflected wavefield for the two stages of reflection waveform inversion. In the third stage, we perform full waveform inversion for the acoustic transversely isotropic media with a vertical symmetry axis to delineate the high-wavenumber structures. For this stage, the medium is described by a combination of the horizontal velocity , η and ε instead of , η and δ. The acoustic multi-parameter full waveform inversion utilizes the diving waves to improve the background as well as utilizes reflection for high-resolution information. Finally, we test our inversion algorithm on the modified Sigsbee 2A model (a salt free part) and a two-dimensional line from a three-dimensional ocean bottom cable dataset. The results demonstrate that the proposed reflection waveform inversion approach can recover the background model for acoustic transversely isotropic media with a vertical symmetry axis starting from an isotropic model. This recovered background model can mitigate the cycle skipping of full waveform inversion and help the inversion recover higher resolution structures.     

Modelling of viscoacoustic wave propagation in transversely isotropic media using decoupled fractional Laplacians

A new wave equation is derived for modelling viscoacoustic wave propagation in transversely isotropic media under acoustic transverse isotropy approximation. The formulas expressed by fractional Laplacian operators can well model the constant-Q (i.e. frequency-independent quality factor) attenuation, anisotropic attenuation, decoupled amplitude loss and velocity dispersion behaviours. The proposed viscoacoustic anisotropic equation can keep consistent velocity and attenuation anisotropy effects with that of qP-wave in the constant-Q viscoelastic anisotropic theory. For numerical simulations, the staggered-grid pseudo-spectral method is implemented to solve the velocity–stress formulation of wave equation in the time domain. The constant fractional-order Laplacian approximation method is used to cope with spatial variable-order fractional Laplacians for efficient modelling in heterogeneous velocity and Q media. Simulation results for a homogeneous model show the decoupling of velocity dispersion and amplitude loss effects of the constant-Q equation, and illustrate the influence of anisotropic attenuation on seismic wavefields. The modelling example of a layered model illustrates the accuracy of the constant fractional-order Laplacian approximation method. Finally, the Hess vertical transversely isotropic model is used to validate the applicability of the formulation and algorithm for heterogeneous media.     

各向异性介质中的弹性阻抗及其反演

地震反演已成为油藏描述中的重要组成部分.绝大多数的常规地震反演是叠后地震数据体上进行,很少考虑各向异性存在的情况.随着勘探开发的发展,地震各向异性和叠前地震波阻抗反演引起了人们极大关注.本文在各向同性介质中弹性阻抗研究基础上,推导出了各向异性介质中的弹性阻抗方程,提出了地震各向异性介质中用弹性阻抗进行储层参数描述的技术路线和框架,并对反演过程中存在的问题进行了有益探讨.