Tilted Transversely Isotropic Research Articles

Elastic anisotropic finite-difference full-wave modeling and imaging of 2D tilted transversely isotropic (TTI) media

Imaging below the tilted transversely isotropic (TTI) media pose serious problems and distortions of the subsurface geological targets of interest for hydrocarbon exploration, which is very common in fold-thrust belts or active tectonic areas like the Rocky Mountain foothills of the Canada or foothills of the Himalaya. To model and image below the TTI media having symmetry axis orthogonal to dipping anisotropic layers, a 2D TTI thrust sheet is considered for elastic anisotropic finite-difference full-wave modeling and generation of synthetic seismic data using the robust staggered-grid scheme of wave propagation. Free-surface boundary conditions at the top and absorbing boundary conditions for the other three sides of the model have been imposed to reduce undesirable edge effects and suppress dispersions. The synthetic seismic data generated for the 2D TTI thrust model are migrated using the anisotropic Kirchhoff pre-stack time (PSTM) and depth migration (PSDM) technique followed by the migration velocity analysis (MVA) algorithm. The traveltime computations are made using robust eikonal equations with fourth-order elastic anisotropic finite-difference solutions to handle large tilts (ν) along with anisotropic ray-tracing to image below the TTI thrust sheet. The MVA algorithm adopted also estimates the anisotropic parameters ε and δ accurately within the acceptable error limit through successive iterations of inversion and parameter update, which results in sufficient flattening of the reflection events in the common image gathers (CIGs). This is a powerful and stable methodology to handle large tilt (ν = 60°) of the steeply dipping reflectors in complex geological structure like the TTI thrust model without much distortion of CIGs for which the symmetry axis of each TTI block is set orthogonal to the bottom of the reflectors having different tilts.

Q-compensated least-squares reverse time migration in TTI media using the visco-acoustic TTI wave equation based on the SLS model

The anelasticity and anisotropy widely exist in real subsurface media. Strong anelasticity will lead to phase dispersion and amplitude attenuation during seismic wave propagation. Anisotropy cause seismic waves to have obviously different kinematic characteristics from that of isotropy. For seismic migration, ignoring the anelasticity and anisotropy of subsurface media will significantly reduce the quality of migration images, even cause imaging failure. We propose a Q-compensated least-squares reverse time migration (Q-LSRTM) in tilted transversely isotropic (TTI) media to correct these effects. According to the Born approximation of seismic wave equation, a linearized visco-acoustic TTI pure qP-wave modeling operator is derived using a new visco-acoustic TTI wave equation for one standard linear solid (SLS) model, which can deal with the anelasticity and can simulate pure qP-wave steadily in attenuating anisotropic media without qSV-wave artifacts. Then, the corresponding adjoint equation is formulated using the adjoint-state method to calculate the gradient sensitivity kernel for the visco-acoustic TTI media. Because of the least-squares inversion, the Q compensation can be achieved during the iterations, so that the over-amplification of noises can be avoided naturally. In addition, compared with conventional LSRTM, the proposed method compensates for the anelasticity and corrects the anisotropy, so as to produce images with better spatial resolution and amplitude fidelity. Numerical examples demonstrate the feasibility and advantages of the proposed method for the data including strong attenuation effects over conventional LSRTM.

Modeling of pure visco-qP-wave propagation in attenuating tilted transversely isotropic media based on decoupled fractional Laplacians

The pseudoviscoacoustic anisotropic wave equation is widely used in the oil and gas industry for modeling wavefields in attenuating anisotropic media. Compared to the full viscoelastic anisotropic wave equation, it can greatly reduce the computational cost of wavefield modeling while maintaining the visco-qP-wave kinematics very well. However, even if we place the source in a thin isotropic layer, there will be some unwanted S-wave artifacts in the qP wavefield simulated by the pseudoviscoacoustic anisotropic wave equation due to the stepped approximation of inclined layer interfaces. Furthermore, the wavefield simulated by the pseudoviscoacoustic anisotropic wave equation may suffer from numerical instabilities when the anisotropy parameter epsilon is less than delta. To overcome these problems, we derive a pure-viscoacoustic tilted transversely isotropic (TTI) wave equation in media with anisotropy in velocity and attenuation based on the exact complex-valued phase velocity formula. The pure-viscoacoustic TTI wave equation has decoupled amplitude dissipation and phase dispersion terms, which is suitable for further reverse time migration with Q compensation. For numerical simulations, we adopt the second-order Taylor series expansion to replace the mixed-domain spatially variable fractional Laplacian operator, which guarantees the decoupling of the wavenumber from the space-related fractional order. Then, we use an efficient and stable hybrid finite-difference and pseudospectral method (HFDPSM) to solve the pure-viscoacoustic TTI wave equation. Numerical tests indicate that the simulation results of the newly derived pure-viscoacoustic TTI wave equation are stable, free from S-wave artifacts, and accurate. We further demonstrate that HFDPSM outperforms the pseudospectral method in terms of numerical simulation stability and computing efficiency.

Detection of Natural Tilted Fractures From Azimuthal Seismic Amplitude Data Based on Linear-Slip Theory

Tilted transverse isotropy (TTI) is common for naturally fractured reservoir when the fractures exist a tilted axis of symmetry. Seismic characterization of natural tilted fractures from observable azimuthal reflected amplitude data, however, is quite difficult. We focus on the detection of natural tilted fractures from azimuthal seismic amplitude data, especially for low- and high-angle tilted fractures. Using the linear-slip theory, we first express the effective elastic stiffness matrix of a TTI medium as a function of background elastic moduli, fracture density, and three angles, including the incident angle, the azimuth angle, and the tilt angle. Based on the scattering function and the first-order perturbations in stiffness components across a weak-contrast reflection interface separating two weak-anisotropy TTI media, we then derive a weak-anisotropy and linearized PP-wave reflection coefficient containing the tilt angle and fracture density. For low- and high-angle aligned fracture sets, we formulate the low- and high-angle approximate PP-wave reflection coefficient equations, respectively, and propose a Bayesian seismic inversion approach to estimate the tilt angle and the corresponding fracture densities using the azimuthal differences in seismic reflected amplitude data. Synthetic and real datasets are used to demonstrate the feasibility and reliability of the proposed inversion approach. Test results make us believe that our proposed inversion approach allow us to obtain the tilted fracture properties of hydrocarbon reservoirs in a more accurate manner than previous cases of vertically transverse isotropy (VTI) or horizontally transverse isotropy (HTI).

Fracture Detection With Azimuthal Seismic Amplitude Difference Inversion in Weakly Monoclinic Medium

The rock with a vertical group of aligned fractures in the tilted transverse isotropy (TTI) background medium is regarded as an effective monoclinic medium. Thomsen parameters and fracture weaknesses can reflect the development degree of TTI background anisotropy and vertical fractures, respectively. Thus, this paper first establishes the approximate stiffness matrix of the monocline medium through the Bond transformation from Thomsen weak anisotropy theory and effective medium theory. Moreover, considering the Born scattering theory, a linearized approximate formula of PP wave reflection coefficient is further deduced. This approximate formula contains the parameters of incident angle, azimuth, tilted angle, compressional and shear moduli, density, fracture weaknesses, and Thomsen parameters. Then, the contribution of each parameter on the reflection coefficient by numerical simulation is analyzed respectively. Finally, from the perspective of Bayesian inference, an azimuthal seismic amplitude difference inversion approach incorporating Cauchy prior information and low-frequency regularization of model parameters is proposed to realize the simultaneous estimation of fracture weaknesses and Thomsen parameters. The synthetic test demonstrates that the inversion approach proposed in this article can obtain information on vertical fractures and TTI background medium and can be applied to azimuthal seismic data with moderate noise well. Furthermore, the inversion results of the synthetic test and field data illustrate the stability and feasibility of the inversion approach and also confirm the fracture weaknesses and Thomsen parameters of the monoclinic medium can be reasonably estimated by the inversion equation derived in this paper.

An iterative factored topography-dependent eikonal solver for anisotropic media

Accurate and efficient eikonal solvers for heterogeneous media play an important role in many areas of seismology, such as seismic tomography, migration, and earthquake localization. Incorporating seismic anisotropy and complex topography remain a computational challenge for finite-difference eikonal solvers. In recent years, the topography-dependent eikonal equation (TDEE) has been proposed as an effective way to calculate seismic traveltimes for isotropic and anisotropic media with irregular topography. However, the Lax-Friedrichs sweeping method used in previous studies to approximate the viscosity solution of TDEE for anisotropic media is more dissipative and needs a much higher number of iterations to converge. In addition, the TDEE solution for the initial point source has an upwind source singularity, which makes all TDEE solvers, even the high-order ones, exhibit polluted convergence and relatively large errors that propagate from the point source to the entire computational domain. To solve these problems, we have formulated the factored topography-dependent anisotropic eikonal (FTDAE) equation in tilted transversely isotropic (TTI) media using the factorization principle. Then, the resulting quartic equation can be numerically solved by using a fixed-point iteration technique based on the simpler elliptical anisotropic eikonal (EAE) equation with a high-order source term due to anelliptical anisotropy introduced by TTI media. At each iteration, the unknown traveltime in the EAE equation is factored into two functions: One of the functions is specified analytically to capture the source singularity, such that the unknown factor is differentiable in the source neighborhood and could be solved by the fast sweeping method. Numerical examples indicate that our FTDAE solver can treat source singularity successfully and achieve high accuracy after just a few iterations, independently of the mesh size, which could provide a more efficient and robust tool for traveltime calculation in the presence of seismic anisotropy and complex surfaces.

Anisotropic Born scattering for the qP scalar wavefield using a low-rank symbol approximation

We have developed a procedure to derive low-rank evolution operators in the mixed space-wavenumber domain for modeling the qP Born-scattered wavefield at perturbations of an anisotropic medium under the pseudoacoustic approximation. To approximate the full wavefield, this scattered field is then added to the reference wavefield obtained with the corresponding low-rank evolution operator in the background medium. Being built upon a Hamiltonian formulation using the dispersion relation for qP-waves, this procedure avoids pseudo-S-wave artifacts and provides a unified approach for linearizing anisotropic pseudoacoustic evolution operators. Therefore, it is immediately applicable to any arbitrary class of anisotropy. As an additional asset, the scattering operators explicitly contain the sensitivity kernels of the Born-scattered wavefield with respect to the anisotropic medium parameters. This enables direct access to important information such as its offset dependence or directional characteristics as a function of the individual parameter perturbations. For our numerical tests, we specify the operators for a mildly anisotropic tilted transversely isotropic (TTI) medium. We validate our implementation in a simple model with weak contrasts and simulate reflection data in the BP TTI model to indicate that the procedure works in a more realistic scenario. The Born-scattering results indicate that our procedure is applicable to strongly heterogeneous anisotropic media. Moreover, we use the analytical capabilities of the kernels by means of sensitivity tests to demonstrate that using two different medium parameterizations leads to different results. The mathematical formulation of the method is such that it allows for an immediate application to least-squares migration in pseudoacoustic anisotropic media.

Decoupled approximation and separate extrapolation of P- and SV-waves in transversely isotropic media

Characterizing the kinematics of seismic waves in elastic vertical transversely isotropic (VTI) media involves four independent parameters. To reduce the complexity, the acoustic approximation for P-waves reduces the number of required parameters to three by setting the vertical S-wave velocity to zero. However, because only the SV-wave phase velocities parallel or perpendicular to the symmetry axis are indirectly set to zero, the acoustic approximation leads to coupled P-wave components and SV-wave artifacts. The new acoustic approximation suggests setting the vertical S-wave velocity as a phase angle-dependent variable so that the SV-wave phase velocity is zero at all phase angles. We find that manipulating this parameter is a valid approach for P-wave approximation but doing so inevitably leads to zero- or nonzero-valued spurious SV-wave components. Thus, we have developed a novel approach to efficiently approximate and thoroughly separate the two wave modes in VTI media. First, the exact P- and SV-wave phase velocity expressions are rewritten by introducing an auxiliary function. After confirming the insensitivity of this function, we construct a new expression for it and obtain simplified P- and SV-wave phase velocity expressions, which are three and four parameters, respectively. This approximation process leads to the same reasonable error for both wave modes. Accuracy analysis indicates that, for the P-wave, the overall accuracy performance of our approach is comparable to that of some existing three-parameter approximations. We then derive the corresponding P- and SV-wave equations in tilted transversely isotropic (TTI) media and provide two available solutions, the hybrid finite-difference/pseudospectral scheme and the low-rank approach. Numerical examples illustrate the separability and high accuracy of the proposed P- and SV-wave simulation methods in TTI media.

An iterative fast sweeping method for the eikonal equation in 2D anisotropic media on unstructured triangular meshes

Computation of traveltimes and raypaths is important for anisotropic tomography inversions. The eikonal-equation-based method outperforms traditional ray methods by producing more accurate results. However, most existing eikonal solvers are formulated on structured regular meshes, which are no longer accurate for models with the presence of irregular topography and subsurface interfaces. To solve the eikonal equation in vertically transversely isotropic (VTI) or tilted transversely isotropic (TTI) models with irregular geometry, we have formulated a new iterative fast sweeping method (FSM) on unstructured triangular meshes. The fixed-point iteration is implemented to capture the high-order nonlinear terms therein, and an FSM on unstructured triangular meshes is implemented to solve the resulting elliptically anisotropic eikonal equation at every iteration. We test the new algorithm for direct arrivals and reflected arrivals, and then we use the calculated traveltimes to track the raypath in VTI/TTI media. Numerical tests demonstrate the validity and accuracy of the new method for models with rough topography and subsurface interface.

Azimuthal amplitude variation with offset parameterization and inversion for fracture weaknesses in tilted transversely isotropic media

The characterization of fracture-induced tilted transverse isotropy (TTI) seems to be more suitable to actual scenarios of geophysical exploration for fractured reservoirs. Fracture weaknesses enable us to describe fracture-induced anisotropy. With the incident and reflected PP-wave in TTI media, we have adopted a robust method of azimuthal amplitude variation with offset (AVO) parameterization and inversion for fracture weaknesses in a fracture-induced reservoir with TTI symmetry. Combining the linear-slip model with the Bond transformation, we have derived the stiffness matrix of a dipping-fracture-induced TTI medium characterized by normal and tangential fracture weaknesses and a tilt angle. Integrating the first-order perturbations in the stiffness matrix of a TTI medium and scattering theory, we adopt a method of azimuthal AVO parameterization for PP-wave reflection coefficient for the case of a weak-contrast interface separating two homogeneous weakly anisotropic TTI layers. We then adopt an iterative inversion method by using the partially incidence-angle-stacked seismic data with different azimuths to estimate the fracture weaknesses of a TTI medium when the tilt angle is estimated based on the image well logs prior to the seismic inversion. Synthetic examples confirm that the fracture weaknesses of a TTI medium are stably estimated from the azimuthal seismic reflected amplitudes for the case of moderate noise. A field data example demonstrates that geologically reasonable results of fracture weaknesses can be determined when the tilt angle is treated as the prior information. We determine that the azimuthal AVO inversion approach provides an available tool for fracture prediction in a dipping-fracture-induced TTI reservoir.

A fast-marching eikonal solver for tilted transversely isotropic media

Fast and accurate traveltime computation for quasi-P waves in anisotropic media is an essential ingredient of many seismic processing and interpretation applications such as Kirchhoff modeling and migration, microseismic source localization, and traveltime tomography. Fast-sweeping methods are widely used for solving the anisotropic eikonal equation due to their flexibility in solving general equations compared to the fast-marching method. However, it has been observed that fast sweeping can be much less efficient than fast marching for models with curved characteristics and practical grid sizes. By representing a tilted transversely isotropic (TTI) equation as a sequence of elliptically isotropic (EI) eikonal equations, we determine that the fast-marching algorithm can be used to compute fast and accurate traveltimes for TTI media. The tilt angle is absorbed into the description of the effective EI model; therefore, the adopted approach does not compromise on the solution accuracy. Through tests on benchmark synthetic models, we test our fast-marching algorithm and discover considerable improvement in accuracy by using factorization and a second-order finite-difference stencil. The adopted methodology opens the door to the possibility of using the fast-marching algorithm for a wider class of anisotropic eikonal equations.

Q-compensated reverse time migration in tilted transversely isotropic media

Anisotropy and absorption are critical to the modeling and analysis of seismic amplitude,phase, and traveltime data. To neglect any of these phenomena, which are often bothoperating simultaneously, degrades the resolution and interpretability of migrated images.However, a full accounting of anisotropy and anelasticity is computationally complex andexpensive. One strategy for accommodating these aspects of wave propagation, while keepingcost and complexity under control, is to do so within an acoustic approximation. Weset up a procedure for solving the time-domain viscoacoustic wave equation for tilted transverselyisotropic (TTI) media, based on a standard linear solid model and, from this, developa viscoacoustic reverse time migration (Q-RTM) algorithm. In this approach, amplitudecompensation occurs within the migration process through a manipulation of attenuationand phase dispersion terms in the time domain differential equations. Specifically, theback-propagation operator is constructed by reversing the sign only of the amplitude lossoperators, but not the dispersion-related operators, a step made possible by reformulatingthe absorptive TTI equations such that the loss and dispersion operators appear separately.The scheme is tested on synthetic examples to examine the capacity of viscoacoustic RTM to correct for attenuation, and the overall stability of the procedure.

Broadband finite-difference Q-compensated reverse time migration algorithm for tilted transverse isotropic media

[Formula: see text]-compensated reverse time migration (Q-RTM) is effective for improving the seismic imaging quality degraded by low-[Formula: see text] anomalies. However, it is difficult to apply existing pseudospectral-based Q-RTM methods to large-scale problems due to the obstacles to high-efficient parallelization posed by the global pseudospectral operators. On the other hand, finite-difference Q-RTM is intrinsically appropriate for domain decomposition and parallel computation, thus being suitable for industrial-sized problems but facing a twofold challenge: (1) to effectively compensate the phase in a broad-bandwidth sense during the wave back-propagation process and (2) to accurately handle the tilted transverse isotropic (TTI) medium with attenuation. We have developed a new framework of finite-difference Q-RTM algorithm by expanding the linear viscoacoustic constitutive relation to a series of integer-order differential terms and a unique integral term that can decouple the amplitude and the phase to allow accurate compensation in a broader frequency range. This framework has two typical implementations: (1) optimizing the frequency-dependent phase velocity (while fixing the negative constant [Formula: see text]) and (2) optimizing the [Formula: see text] value (while fixing the frequency-dependent phase velocity). We generalize this broadband finite-difference Q-RTM algorithm to TTI media, where an artificial [Formula: see text] is applied to suppress the shear-wave artifacts induced by the acoustic TTI approximation. Numerical examples demonstrate that this Q-RTM method accurately compensates the phase and amplitude in a broad frequency range of 5–70 Hz and produces high-quality images. Due to the local nature of finite-difference operators, this algorithm is expected to outperform the existing pseudospectral-based Q-RTM methods in terms of computational efficiency and implementation convenience for real-world Q-RTM projects.

Multistage elastic full-waveform inversion for tilted transverse isotropic media

SUMMARYSeismic anisotropy is an important physical phenomenon that significantly affects wave propagation in complex sedimentary basins. When geological structures exhibit steep dips or severe folding, the symmetry axis of the transversely isotropic (TI) representation of the region can be rotated, leading to tilted transversely isotropic (TTI) media. We seek to find the optimal full-waveform inversion (FWI) strategy to estimate both the seismic velocities and the anisotropic parameters, including the tilt angle, in the presence of elastic TTI media. We first formulate the forward and inverse problems for elastic TTI media and analyse the radiation patterns of the model parameters. Based on the analyses of the radiation patterns, we propose two similar multistage FWI strategies that add inversion parameters over three stages, beginning with the isotropic parameters (horizontal P- and vertical S-wave velocity) and moving to the anisotropic parameters; the tilt angle is directly inverted in the last stage. Since diving waves, which are useful for providing long-wavelength updates, are mainly controlled by horizontal motion in anisotropic media, it is reasonable to choose the horizontal P-wave velocity rather than the vertical P-wave velocity. Then, the anisotropic parameters are inverted mainly using the reflected waves based on the isotropic background model built in the first stage. The main difference between the two multistage FWI strategies is whether the anisotropic parameter η is inverted. Comparing the two multistage FWI strategies with the simultaneous inversion strategy for a downsized version of the synthetic BP TTI model, we confirm that the multistage FWI strategies yield better inversion results than the simultaneous inversion strategy. When we compare the two multistage FWI strategies with each other for surface seismic data, ignoring η during the FWI process (focused multistage FWI) yields better inversion results for the tilt angle than those obtained with the inversion of η because η has less influence on the FWI than the other parameters and is not recovered well, which plays a role in degrading the tilt angle. Numerical examples support our conclusions that the focused multistage FWI strategy (neglecting η) is the optimal FWI strategy for TTI media and achieves computational efficiency for surface seismic data.

Least-squares reverse time migration in TTI media using a pure qP-wave equation

Anisotropy is a common phenomenon in subsurface strata and should be considered in seismic imaging and inversion. Seismic imaging in a vertical transversely isotropic (VTI) medium does not take into account the effects of the tilt angles, which can lead to degraded migrated images in areas with strong anisotropy. To correct such waveform distortion, reduce related image artifacts, and improve migration resolution, a tilted transversely isotropic (TTI) least-squares reverse time migration (LSRTM) method is presented. In the LSRTM, a pure qP-wave equation is used and solved with the finite-difference method. We have analyzed the stability condition for the pure qP-wave equation using the matrix method, which is used to ensure the stability of wave propagation in the TTI medium. Based on this wave equation, we derive a corresponding demigration (Born modeling) and adjoint migration operators to implement TTI LSRTM. Numerical tests on the synthetic data show the advantages of TTI LSRTM over VTI RTM and VTI LSRTM when the recorded data contain strong effects caused by large tilt angles. Our numerical experiments illustrate that the sensitivity of the adopted TTI LSRTM to the migration velocity errors is much higher than that to the anisotropic parameters (including epsilon, delta, and tilted angle parameters), and its sensitivity to the epsilon model and tilt angle is higher than that to the delta model.

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.

Efficient mixed-domain pure P-wave equations for 2D and 3D tilted transversely isotropic media

In transversely isotropic (TI) media, accurate and efficient pure P-wave extrapolation is the basis of seismic imaging and inversion algorithms. For pure P-wave equations in TI media, combining the finite-difference (FD) method with the pseudo-spectral (PS) method is an effective solution. However, if the mixed-domain pure P-wave equation contains multiple wavenumbers, using the hybrid approach will involve multiple Fourier transforms, which will result in high computational costs. Referring to the weak anisotropy approximation, we propose a new approximate P-wave phase velocity expression. We then use an optimization strategy to reduce the number of wavenumber terms in the corresponding dispersion relations and derive the mixed-domain pure P-wave equations in 2D and 3D tilted transversely isotropic (TTI) media. Through numerical experiments in 2D and 3D TTI models, we verify the feasibility and efficiency of the proposed mixed-domain P-wave equations.

Estimating tilted fracture weaknesses from azimuthal differences in seismic amplitude data

Tilted transverse isotropy (TTI) provides a useful model for the elastic response of a medium containing aligned fractures with a symmetry axis oriented obliquely in the vertical and horizontal coordinate directions. Robust methods for determining the TTI properties of a medium from seismic observations to characterize fractures are sought. Azimuthal differencing of seismic amplitude data produces quantities that are particularly sensitive to TTI properties. Based on the linear slip fracture model, we express the TTI stiffness matrix in terms of the normal and tangential fracture weaknesses. Perturbing stiffness parameters to simulate an interface separating an isotropic medium and a TTI medium, we derive a linearized P-to-P reflection coefficient expression in which the influence of tilt angle and fracture weaknesses separately emerge. We formulate a Bayesian inversion approach in which amplitude differences between seismic data along two azimuths, interpreted in terms of the reflection coefficient approximation, are used to determine fracture weaknesses and tilt angle. Tests with simulated data confirm that the unknown parameter vector involving fracture weakness and tilted fracture weaknesses is stably estimated from seismic data containing a moderate degree of additive Gaussian noise. The inversion approach is applied to a field surface seismic data acquired over a fractured reservoir; from it, interpretable tilted fracture weaknesses, consistent with expected reservoir geology, are obtained. We determine that our inversion approach and the established inversion workflow can produce the properties of systems of tilted fractures stably using azimuthal seismic amplitude differences, which may add important information for characterization of fractured reservoirs.

Modeling of pure qP- and qSV-waves in tilted transversely isotropic media with the optimal quadratic approximation

Seismic forward modeling in tilted transverse isotropic (TTI) media is crucial for the application of reverse time migration and full-waveform inversion. Modeling based on conventional coupled pseudoacoustic wave equations not only generates SV-wave artifacts, but it also suffers from instabilities in which the anisotropy parameter [Formula: see text]. To address these issues, we have started with the exact vertical transversely isotropic phase velocity formula and developed novel pure qP- and qSV-wave governing equations in TTI media by using the optimal quadratic approximation. For the convenience of using finite-difference (FD) method to solve the new pure qP- and qSV-wave wave equations, we decompose the equations into a combination of a time-space-domain wave equation that can be solved by the FD method and a Poisson equation that can be solved by the pseudospectral method. We find that the high-frequency errors caused by the pseudospectral method and the usual truncation errors in FD schemes may be responsible for the instability of the numerical simulations. To stabilize the computation, we design a 2D low-pass filtering operator to eliminate severe high-frequency numerical noise. Several numerical examples demonstrate that modeling using the new pure qP-wave equations does not have qSV-wave artifacts interference and is stable for [Formula: see text]. Our results indicate that our method can achieve highly accurate and stable modeling results even in extremely complex TTI media.

Tutorial — time conversion of depth migrated data: part II, TTI preSDM

Contemporary depth imaging projects often require final pre-Stack Depth-Migrations (preSDM) to be converted to time for comparison to vintage pre-Stack Time-Migrations (preSTM), or to facilitate conversion to ‘geological’ depth through calibration to well and check-shot data. Here, I consider the situation of having performed an anisotropic Tilted Transverse Isotropy (TTI) preSDM and wanting to convert it to time via vertical stretch in order to compare it, say, to an anisotropic preSTM. Such a comparison is inherently invalid, as time-migration will explicitly treat any anisotropy as if it were Vertical Transverse Isotropy (VTI), and, in addition, the lateral positioning error inherent in preSTM will render such comparisons questionable on steeply dipping structures. I show here that the most appropriate type of ‘velocity’ to use for conversion to time of TTI preSDM reflection events should be the vertical component of the phase velocity. Conversely, if we are considering point-to-point measurements, such as the direct arrival travel time, a down-hole or check-shot measurement, then the group velocity should be used, as it is with this speed that energy travels. In addition, subsequent depth conversion of any time product for interpretational purposes would best be accomplished using a velocity calibrated to well check-shots.