On the instability in second-order systems for acoustic VTI and TTI media

Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. T171-T186 ◽  
Author(s):  
Kenneth P. Bube ◽  
Tamas Nemeth ◽  
Joseph P. Stefani ◽  
Ray Ergas ◽  
Wei Liu ◽  
...  
Keyword(s):  
Wave Equation ◽  
Dispersion Relation ◽  
Differential Equations ◽  
Shear Wave ◽  
Transversely Isotropic ◽  
Second Order ◽  
Wave Speeds ◽  
Tti Media ◽  
Second Order Systems ◽  
Vti Media

We studied second-order wave propagation systems for vertical transversely isotropic (VTI) and tilted transversely isotropic (TTI) acoustic media with variable axes of symmetry that have their shear-wave speeds set to zero. Acoustic TTI systems are commonly used in reverse-time migration, but these second-order systems are susceptible to instablities appearing as nonphysical stationary noise growing linearly in time, particularly in variable-tilt TTI media. We found an explanation of the cause of this phenomenon. The instabilities are not caused only by the numerical schemes; they are inherent to the differential equations. These instabilities are present even in homogeneous VTI media. These instabilities are caused by zero wave speeds at a wide variety of wavenumbers — a direct consequence of setting the shear-wave speeds to zero — coupled with the second time derivative in these systems. Although the second-order isotropic wave equation allows smooth time-growing solutions, a larger class of time-growing solutions exists for the second-order acoustic TI systems, including nonsmooth solutions. Boundary conditions appear to be less effective in controlling these time-growing solutions than they are for the isotropic wave equation. These systems conserve an incomplete energy that does not prevent the instabilities. The corresponding steady-state systems are no longer elliptic differential equations and can have nonsmooth solutions that are related to the instabilities. We started initially with homogeneous VTI media, and then extended these results to heterogeneous variable-tilt TTI media. We also developed a second-order acoustic system for heterogeneous variable-tilt TTI media derived directly from the full-elastic system for heterogeneous variable-tilt TTI media. All second-order systems with a dispersion relation obtained by setting the shear-wave speeds to zero in the elastic dispersion relation allowed these nonphysical time-growing solutions; however, knowing the cause of these instabilities, it may be possible to prevent or control the activation of these solutions.

First-order systems for elastic and acoustic variable-tilt TI media

Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. T157-T170 ◽  
Author(s):  
Kenneth P. Bube ◽  
Tamas Nemeth ◽  
Joseph P. Stefani ◽  
Wei Liu ◽  
Kurt T. Nihei ◽  
...  
Keyword(s):  
Shear Wave ◽  
Lower Order ◽  
Elastic System ◽  
Transversely Isotropic ◽  
Order System ◽  
Wave Speeds ◽  
First Order ◽  
Vti Media ◽  
Ti Media ◽  
Lower Order Terms

We derive and compare first-order wave propagation systems for variable-tilt elastic and acoustic tilted transversely isotropic (TTI) media. Acoustic TTI systems are commonly used in reverse-time migration. Starting initially with homogeneous vertical transversely isotropic (VTI) media, and then extending to heterogeneous variable-tilt TI media, we derive a pseudoacoustic [Formula: see text] first-order system of differential equations by setting the shear-wave speeds to zero and simplifying the full-elastic system accordingly. This [Formula: see text] system conserves a complete energy, but only when the anelliptic anisotropy parameter [Formula: see text]. For [Formula: see text] (including isotropic media), the system allows linearly time-growing and spatially nonpropagating nonphysical solutions frequently taken for numerical noise. We modified this [Formula: see text] acoustic first-order system by changing the stress variables to obtain a system that stays stable for [Formula: see text]. This system for homogeneous VTI media is generalized to heterogeneous variable-tilt TI media by rotating the stress and strain variables in the full elastic system before setting the shear-wave speeds to zero; the system obtained can be greatly simplified by combining the rotational terms, resulting in only one rotation and extra lower-order terms compared to the [Formula: see text] first-order acoustic system for VTI media. This new system can be simplified further by neglecting the lower-order terms. Both systems (with and without lower-order terms) conserve the same complete energy. Finally, the corresponding [Formula: see text] full elastic system for variable-tilt acoustic TI media can be used for the purposes of benchmarking.

Migration of P‐wave reflection data in transversely isotropic media

Geophysics ◽  
1994 ◽  
Vol 59 (4) ◽  
pp. 591-596 ◽  
Author(s):  
Suhas Phadke ◽  
S. Kapotas ◽  
N. Dai ◽  
Ernest R. Kanasewich
Keyword(s):  
Wave Equation ◽  
Dispersion Relation ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Limited Range ◽  
Horizontal Velocity ◽  
P Wave ◽  
Transversely Isotropic Media ◽  
Isotropic Media

Wave propagation in transversely isotropic media is governed by the horizontal and vertical wave velocities. The quasi‐P(qP) wavefront is not an ellipse; therefore, the propagation cannot be described by the wave equation appropriate for elliptically anisotropic media. However, for a limited range of angles from the vertical, the dispersion relation for qP‐waves can be approximated by an ellipse. The horizontal velocity necessary for this approximation is different from the true horizontal velocity and depends upon the physical properties of the media. In the method described here, seismic data is migrated using a 45-degree wave equation for elliptically anisotropic media with the horizontal velocity determined by comparing the 45-degree elliptical dispersion relation and the quasi‐P‐dispersion relation. The method is demonstrated for some synthetic data sets.

scholarly journals Nonlinear Green’s Functions for Wave Equation with Quadratic and Hyperbolic Potentials

2018 ◽  
Vol 2018 ◽  
pp. 1-9 ◽  
Author(s):  
Asatur Zh. Khurshudyan
Keyword(s):  
Wave Equation ◽  
Differential Equations ◽  
Wave Equations ◽  
Separation Of Variables ◽  
Linear Equations ◽  
Second Order ◽  
Nonlinear Wave Equations ◽  
One Dimensional ◽  
Exponential Nonlinearity ◽  
Second Order Differential Equations

The advantageous Green’s function method that originally has been developed for nonhomogeneous linear equations has been recently extended to nonlinear equations by Frasca. This article is devoted to rigorous and numerical analysis of some second-order differential equations new nonlinearities by means of Frasca’s method. More specifically, we consider one-dimensional wave equation with quadratic and hyperbolic nonlinearities. The case of exponential nonlinearity has been reported earlier. Using the method of generalized separation of variables, it is shown that a hierarchy of nonlinear wave equations can be reduced to second-order nonlinear ordinary differential equations, to which Frasca’s method is applicable. Numerical error analysis in both cases of nonlinearity is carried out for various source functions supporting the advantage of the method.

Some second-order systems of partial differential equations of composite type

1987 ◽  
Vol 106 (1-2) ◽  
pp. 73-88 ◽  
Author(s):  
Lin Wei
Keyword(s):  
Partial Differential Equations ◽  
Differential Equations ◽  
Hyperbolic Systems ◽  
Second Order ◽  
Cauchy Type ◽  
Linear Transformations ◽  
Partial Differential ◽  
Composite Systems ◽  
Composite Type ◽  
Second Order Systems

SynopsisThe Cauchy problem and the Dirichlet-Cauchy type problem of some second-order systems of partial differential equations of composite type of two unknown functions are investigated. Such systems possess some of the characteristics not only of elliptic but also of hyperbolic systems in the same domain. Representations of the solutions are found for the upper half plane. To this end, the composite systems are reduced to the canonical form by means of successive applications of three kinds of linear transformations. Function theoretic methods are used to obtain representation formulae. Furthermore, some composite systems of 2m-unknown function are also considered.

Hybrid Fourier finite-difference 3D depth migration for anisotropic media

Geophysics ◽  
2008 ◽  
Vol 73 (2) ◽  
pp. S27-S34 ◽  
Author(s):  
Tong W. Fei ◽  
Christopher L. Liner
Keyword(s):  
Wave Equation ◽  
Finite Difference ◽  
Acoustic Wave ◽  
Seismic Data ◽  
Transversely Isotropic ◽  
Acoustic Wave Equation ◽  
Near Surface ◽  
Depth Migration ◽  
Surface Distortions ◽  
Vti Media

When a subsurface is anisotropic, migration based on the assumption of isotropy will not produce accurate migration images. We develop a hybrid wave-equation migration algorithm for vertical transversely isotropic (VTI) media based on a one-way acoustic wave equation, using a combination of Fourier finite-difference (FFD) and finite-difference (FD) approaches. The hybrid method can suppress an additional solution that exists in the VTI acoustic wave equation, and it offers speed and other advantages over conventional FFD or FD methods alone. The algorithm is tested on a synthetic model involving log data from onshore eastern Saudi Arabia, including estimates of both intrinsic and layer-induced VTI parameters. Results indicate that VTI imaging in this region offers some improvement over isotropic imaging, primarily with respect to subtle structure and stratigraphy and to image continuity. These benefits probably will be overshadowed by perennial land seismic data issues such as near-surface distortions and multiples.

A separable strong-anisotropy approximation for pure qP-wave propagation in transversely isotropic media

Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. C337-C354 ◽  
Author(s):  
Jörg Schleicher ◽  
Jessé C. Costa
Keyword(s):  
Wave Propagation ◽  
Wave Equation ◽  
Dispersion Relation ◽  
Transversely Isotropic ◽  
Wave Dispersion ◽  
Low Rank ◽  
S Wave ◽  
Finite Difference Approximations ◽  
Transversely Isotropic Media ◽  
Isotropic Media

The wave equation can be tailored to describe wave propagation in vertical-symmetry axis transversely isotropic (VTI) media. The qP- and qS-wave eikonal equations derived from the VTI wave equation indicate that in the pseudoacoustic approximation, their dispersion relations degenerate into a single one. Therefore, when using this dispersion relation for wave simulation, for instance, by means of finite-difference approximations, both events are generated. To avoid the occurrence of the pseudo-S-wave, the qP-wave dispersion relation alone needs to be approximated. This can be done with or without the pseudoacoustic approximation. A Padé expansion of the exact qP-wave dispersion relation leads to a very good approximation. Our implementation of a separable version of this equation in the mixed space-wavenumber domain permits it to be compared with a low-rank solution of the exact qP-wave dispersion relation. Our numerical experiments showed that this approximation can provide highly accurate wavefields, even in strongly anisotropic inhomogeneous media.

Sign in / Sign up

Export Citation Format

Share Document