This study attempts to validate a mathematical formalism of introducing attenuation into Schoenberg’s linear slip model. This formalism is based on replacing the real-valued weaknesses by complex-valued ones. During an ultrasonic experiment, performed at a central frequency of [Formula: see text] on a plate-stack model with [Formula: see text]-thick Plexiglas™ plates, the velocity and attenuation (inverse of the quality factor [Formula: see text]) of P-, SH-, and SV-waves are measured in directions from 25° to 90° with the symmetry axis for dry and oil-saturated models and loading uniaxial pressures of 2 and [Formula: see text]. The velocity and attenuation data are fitted by the derived theoretical functions. The values of the real and imaginary parts of the complex-valued weaknesses are estimated. Thereal parts of the weaknesses, which have a clear physical meaning (they affect the weakening of the material), are three times larger for the dry model than for the oil-saturated one. The imaginary parts of the weaknesses are responsible for attenuation; their values are an order of magnitude smaller than the real parts. The derived expressions for angle-dependent velocities and attenuations can be used to distinguish between dry and oil-saturated fractures. In particular, the P-wave attenuation function in the symmetry-axis direction (normal to fracture planes) is different in dry and saturated media. The experiment shows that the plate-stack model is inhomogeneous because of the nonuniform pressure distribution, which degrades the experimental results and creates difficulties in the inversion for the complex-valued weaknesses — particularly in joint inversion of P- and S-wave data.

A wave equation, derived using the acoustic medium assumption for P-waves in transversely isotropic (TI) media with a vertical symmetry axis (VTI media), yields a good kinematic approximation to the familiar elastic wave equation for VTI media. The wavefield solutions obtained using this VTI acoustic wave equation are free of shear waves, which significantly reduces the computation time compared to the elastic wavefield solutions for exploding‐reflector type applications. From this VTI acoustic wave equation, the eikonal and transport equations that describe the ray theoretical aspects of wave propagation in a TI medium are derived. These equations, based on the acoustic assumption (shear wave velocity = 0), are much simpler than their elastic counterparts, yet they yield an accurate description of traveltimes and geometrical amplitudes. Numerical examples prove the usefulness of this acoustic equation in simulating the kinematic aspects of wave propagation in complex TI models.

When transversely isotropic (VTI) media with vertical symmetry axes are characterized using the zero‐dip normal moveout (NMO) velocity [[Formula: see text]] and the anisotropy parameter ηinstead of Thomsen’s parameters, time‐related processing [moveout correction, dip moveout (DMO), and time migration] become nearly independent of the vertical P- and S-wave velocities ([Formula: see text] and [Formula: see text], respectively). The independence on [Formula: see text] and [Formula: see text] is well within the limits of seismic accuracy, even for relatively strong anisotropy. The dependency on [Formula: see text] and [Formula: see text] reduces even further as the ratio [Formula: see text] decreases. In fact, for [Formula: see text], all time‐related processing depends exactly on only [Formula: see text] and η. This fortunate dependence on two parameters is demonstrated here through analytical derivations of time‐related processing equations in terms of [Formula: see text] and η. The time‐migration dispersion relation, the NMO velocity for dipping events, and the ray‐tracing equations extracted by setting [Formula: see text] (i.e., by considering VTI as acoustic) not only depend solely on [Formula: see text] and η but are much simpler than the counterpart expressions for elastic media. Errors attributed to this use of the acoustic assumption are small and may be neglected. Therefore, as in isotropic media, the acoustic model arising from setting [Formula: see text], although not exactly true for VTI media, can serve as a useful approximation to the elastic model for the kinematics of P-wave data. This approximation can boost the efficiency of imaging and DMO programs for VTI media as well as simplify their description.

在各向异性地震波场中,qP波与qS波常常是耦合在一起的.多分量地震数据处理中一个关键环节就是波型分离(即模式解耦),以纵波成分为主的常规单分量地震数据的成像则需要合理描述标量qP波的传播算子.本文作者曾构建了在运动学上同弹性波动方程等价,动力学上突出标量qP波的伪纯模式波动方程.为了彻底消除qS波残余,本文根据波矢量与qP波偏振矢量之间的偏差,提出从伪纯模式波场提取纯模式标量qP波的方法.数值分析展示了投影偏差算子在波数域和空间域的特征.基于不同复杂程度理论模型的试验结果表明,联合"伪纯模式传播算子"与"投影偏差校正"可为各向异性介质分离模式波场传播过程提供一种简便的描述工具.

We present an approach for P-wave modeling in inhomogeneous transversely isotropic media with tilted symmetry axis (TTI media), suitable for anisotropic reverse-time migration. The proposed approach is based on wave equations derived from first principles — the equations of motion and Hooke’s law — under the acoustic TI approximation. Consequently, no assumptions are made about the spatial variation of medium parameters. A rotation of the stress and strain tensors to a local coordinate system, aligned with the TI-symmetry axis, makes it possible to benefit from the simple and sparse form of the TI-elastic tensor in that system. The resulting wave equations can be formulated either as a set of five first-order or as a set of two second-order partial differential equations. For the constant-density case, the second-order TTI wave equations involve mixed and nonmixed second-order spatial derivatives with respect to global, nonrotated coordinates. We propose a numerical implementation of these equations using high-order centered finite differences. To minimize modeling artifacts related to the use of centered first-derivative operators, we use discrete second-derivative operators for the nonmixed second-order spatial derivatives and repeated discrete first-derivative operators for the mixed derivatives. Such a combination of finite-difference operators leads to a stable wave propagator, provided that the operators are designed properly. In practice, stability is achieved by slightly weighting down terms that contain mixed derivatives. This has a minor, practically negligible, effect on the kinematics of wave propagation. The stability of the presented scheme in inhomogeneous TTI models with rapidly varying anisotropic symmetry axis direction is demonstrated with numerical examples.

Acoustic transversely isotropic (TI) media are defined by artificially setting the shear‐wave velocity in the direction of symmetry axis, VS0, to zero. Contrary to conventional wisdom that equating VS0= 0 eliminates shear waves, we demonstrate their presence and examine their properties. Specifically, we show that SV‐waves generally have finite nonzero phase and group velocities in acoustic TI media. In fact, these waves have been observed in full waveform modeling, but apparently they were not understood and labeled as numerical artifacts.

逆时偏移作为一种高精度偏移方法已成为复杂构造成像的重要技术，描述纵波独立传播的延拓方程是各向异性介质逆时偏移的一个关键问题.在对VTI介质几个经典相速度近似公式回顾的基础上，针对常用于描述纯P波的Harlan近似公式在各向异性参数ε较大情况下近似精度较低的问题，本文对Harlan公式中的非椭圆项进行了修正，在非椭圆项前添加了一个与各向异性参数ε有关的修正系数，得到了三种改进型Harlan公式，并以近似精度最高的改进式为基础，推导了TTI介质纯P波方程.针对该伪微分方程，本文利用伪谱法和有限差分法联合实现波场延拓，对于常密度二阶方程，基于中心网格实现；对于一阶应力-速度方程则基于旋转交错网格实现.通过数值试验分析了TTI介质纯P波一阶应力-速度方程的近似精度，并以一阶纯P波方程为基础进行了TTI介质逆时偏移数值模拟试验.结果表明，本文给出的方法能够较准确地描述TTI介质纯P波波场特征，可以应用至各向异性介质逆时偏移.

Modeling and reverse time migration based on the tilted transverse isotropic (TTI) acoustic wave equation suffers from instability in media of general inhomogeniety, especially in areas where the tilt abruptly changes. We develop a stable TTI acoustic wave equation implementation based on the original elastic anisotropic wave equation. We, specifically, derive a vertical transversely isotropic wave system of equations that is equivalent to their elastic counterpart and introduce the self-adjoint differential operators in rotated coordinates to stabilize the TTI acoustic wave equations. Compared to the conventional formulations, the new system of equations does not add numerical complexity; a stable solution can be found by either a pseudospectral method or a high-order explicit finite difference scheme. We demonstrate by examples that our method provides stable and high-quality TTI reverse time migration images.

Conventional modeling and migration for tilted transversely isotropic (TTI) media may suffer from numerical instabilities and shear wave artifacts due to the coupling of the P-wave and SV-wave modes in the TTI coupled equations. Starting with the separated P- and SV-phase velocity expressions for vertical transversely isotropic (VTI) media, we extend these decoupled equations for modeling and reverse time migration (RTM) in acoustic TTI media. Compared with the TTI coupled equations published in the geophysical literature, the new TTI decoupled equations provide a more stable solution due to the complete separation of the P-wave and SV-wave modes. The pseudospectral method is the most convenient method to implement these equations due to the form of wavenumber expressions and has the added benefit of being highly accurate and thus avoiding numerical dispersion. The rapid expansion method (REM) in time is employed to produce a broad band numerically stable time evolution of the wavefields. Synthetic results validate the proposed TTI decoupled equations and show that modeling and RTM in TTI media with the decoupled equations remain numerically stable even for models with strong anisotropy and sharp contrasts.

Reverse time migration (RTM) in complex anisotropic media requires calculation of the propagation of a single-mode wave, the quasi-P-wave. This was conventionally realized by solving a [Formula: see text] system of second-order partial differential equations. The implementation of this [Formula: see text] system required at least twice the computational resources as compared with the acoustic wave equation. The S-waves, an introduced auxiliary function in this system, were treated as artifacts in the RTM. Furthermore, the [Formula: see text] system suffered numerical stability problems at the places in which abrupt changes of symmetric axis of anisotropy exist, which brings more challenges to real data implementation. On the other hand, the Alkhalifah’s equation, which governs the pure quasi-P-wave propagation, was hard to solve because it was a pseudodifferential equation. We proposed a pure quasi-P-wave equation that can be easily implemented within current imaging framework. Our new equation was obtained by decomposing the original pseudodifferential operator into two numerical solvable operators: a differential operator and a scalar operator. The combination of these two operators yielded an accurate phase of quasi-P-wave propagation. Our solution was S-wave free and numerically stable for very complicated models. Because only one equation was required to resolve numerically, the new proposed scheme was more efficient than those conventional schemes that solve the [Formula: see text] second-order differential equations system. For tilted transverse isotropy (TTI) RTM implementation, the required increase of numerical cost was minimal, and we could expect at least a factor of two of improvement of efficiency. We showed the effectiveness and robustness of our method with numerical examples with simple and very complicated TTI models, the SEG Advanced Modeling (SEAM) model. Further extension of our approach to orthorhombic anisotropy or tilted orthorhombic anisotropy was straightforward.