Viscoelastic time-reversal imaging

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. A45-A50 ◽  
Author(s):  
Tieyuan Zhu
Keyword(s):  
Wave Equation ◽  
Time Reversal ◽  
Wave Equations ◽  
Homogeneous Model ◽  
S Wave ◽  
The Third ◽  
Complex Source ◽  
Imaging Approach ◽  
Time Invariance ◽  
Time Invariant

The time invariance of wave equations, an essential precondition for time-reversal (TR) imaging, is no longer valid when introducing attenuation. I evaluated a viscoelastic (VE) TR imaging algorithm based on a novel VE wave equation. By reversing the sign of the P- and S-wave loss operators, the VE wave equation became time invariant for the TR operation. Attenuation effects were thus compensated for during TR wave propagation. I developed the formulations of VE forward modeling and TR imaging. I tested my imaging approach in three numerical experiments. The first experiment used a 2D homogeneous model with full-aperture receivers to examine the time invariance of the VE TR imaging equation. Using the same model, the second experiment was used to demonstrate the method’s ability to characterize a point source. In the third experiment, I applied this method to characterize a complex source using borehole geophones. Numerical results illustrated that the VE TR imaging improved our knowledge of the source location, radiation pattern, and amplitude.

True-amplitude linearized waveform inversion with the quasi-elastic wave equation

Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. R827-R844 ◽  
Author(s):  
Zongcai Feng ◽  
Gerard Schuster
Keyword(s):  
Wave Equation ◽  
Elastic Wave ◽  
Wave Equations ◽  
Waveform Inversion ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Elastic Wave Equation ◽  
P Waves ◽  
Time Migration

We present a quasi-elastic wave equation as a function of the pressure variable, which can accurately model PP reflections with elastic amplitude variation with offset effects under the first-order Born approximation. The kinematic part of the quasi-elastic wave equation accurately models the propagation of P waves, whereas the virtual-source part, which models the amplitudes of reflections, is a function of the perturbations of density and Lamé parameters [Formula: see text] and [Formula: see text]. The quasi-elastic wave equation generates a scattering radiation pattern that is exactly the same as that for the elastic wave equation, and only requires the solution of two acoustic wave equations for each shot gather. This means that the quasi-elastic wave equation can be used for true-amplitude linearized waveform inversion (also known as least-squares reverse time migration) of elastic PP reflections, in which the corresponding misfit gradients are with respect to the perturbations of density and the P- and S-wave impedances. The perturbations of elastic parameters are iteratively updated by minimizing the [Formula: see text]-norm of the difference between the recorded PP reflections and the predicted pressure data modeled from the quasi-elastic wave equation. Numerical tests on synthetic and field data indicate that true-amplitude linearized waveform inversion using the quasi-elastic wave equation can account for the elastic PP amplitudes and provide a robust estimate of the perturbations of P- and S-wave impedances and, in some cases, the density. In addition, true-amplitude linearized waveform inversion provides images with a wider bandwidth and fewer artifacts because the PP amplitudes are accurately explained. We also determine the 2D scalar quasi-elastic wave equation for P-SV reflections and the 3D vector equation for PS reflections.

Approximation of pure acoustic seismic wave propagation in TTI media

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB97-WB107 ◽  
Author(s):  
Chunlei Chu ◽  
Brian K. Macy ◽  
Phil D. Anno
Keyword(s):  
Wave Propagation ◽  
Wave Equation ◽  
Elastic Wave ◽  
Wave Equations ◽  
Computational Cost ◽  
P Wave ◽  
Anisotropy Axis ◽  
S Wave ◽  
Time Migration ◽  
Tti Media

Pseudoacoustic anisotropic wave equations are simplified elastic wave equations obtained by setting the S-wave velocity to zero along the anisotropy axis of symmetry. These pseudoacoustic wave equations greatly reduce the computational cost of modeling and imaging compared to the full elastic wave equation while preserving P-wave kinematics very well. For this reason, they are widely used in reverse time migration (RTM) to account for anisotropic effects. One fundamental shortcoming of this pseudoacoustic approximation is that it only prevents S-wave propagation along the symmetry axis and not in other directions. This problem leads to the presence of unwanted S-waves in P-wave simulation results and brings artifacts into P-wave RTM images. More significantly, the pseudoacoustic wave equations become unstable for anisotropy parameters [Formula: see text] and for heterogeneous models with highly varying dip and azimuth angles in tilted transversely isotropic (TTI) media. Pure acoustic anisotropic wave equations completely decouple the P-wave response from the elastic wavefield and naturally solve all the above-mentioned problems of the pseudoacoustic wave equations without significantly increasing the computational cost. In this work, we propose new pure acoustic TTI wave equations and compare them with the conventional coupled pseudoacoustic wave equations. Our equations can be directly solved using either the finite-difference method or the pseudospectral method. We give two approaches to derive these equations. One employs Taylor series expansion to approximate the pseudodifferential operator in the decoupled P-wave equation, and the other uses isotropic and elliptically anisotropic dispersion relations to reduce the temporal frequency order of the P-SV dispersion equation. We use several numerical examples to demonstrate that the newly derived pure acoustic wave equations produce highly accurate P-wave results, very close to results produced by coupled pseudoacoustic wave equations, but completely free from S-wave artifacts and instabilities.

Constraints on the anisotropic parameters for pseudoelastic vertical transverse isotropy wave equations and the applications on imaging carbonate reservoirs

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. C85-C94 ◽  
Author(s):  
Houzhu (James) Zhang ◽  
Hongwei Liu ◽  
Yang Zhao
Keyword(s):  
Wave Equation ◽  
Elastic Wave ◽  
Wave Velocity ◽  
Wave Equations ◽  
Transverse Isotropy ◽  
S Wave ◽  
Elastic Wave Equation ◽  
Computation Efficiency ◽  
The Stability ◽  
Vertical Transverse Isotropy

Seismic anisotropy is an intrinsic elastic property. Appropriate accounting of anisotropy is critical for correct and accurate positioning seismic events in reverse time migration. Although the full elastic wave equation may serve as the ultimate solution for modeling and imaging, pseudoelastic and pseudoacoustic wave equations are more preferable due to their computation efficiency and simplicity in practice. The anisotropic parameters and their relations are not arbitrary because they are constrained by the energy principle. Based on the investigation of the stability condition of the pseudoelastic wave equations, we have developed a set of explicit formulations for determining the S-wave velocity from given Thomsen’s parameters [Formula: see text] and [Formula: see text] for vertical transverse isotropy and tilted transverse isotropy media. The estimated S-wave velocity ensures that the wave equations are stable and well-posed in the cases of [Formula: see text] and [Formula: see text]. In the case of [Formula: see text], a common situation in carbonate, a positive value of S-wave velocity is needed to avoid the wavefield instability. Comparing the stability constraints of the pseudoelastic- with the full-elastic wave equation, we conclude that the feasible range of [Formula: see text] and [Formula: see text] was slightly larger for the pseudoelastic assumption. The success of achieving high-accuracy images and high-quality angle gathers using the proposed constraints is demonstrated in a synthetic example and a field example from Saudi Arabia.

scholarly journals Comparison of Numerical Methods for SWW Equations

2019 ◽  
Vol 7 (09) ◽  
Author(s):  
Shkelqim Hajrulla ◽  
Leonard Bezati ◽  
Besiana Hamzallari
Keyword(s):  
Wave Equation ◽  
Numerical Methods ◽  
Finite Volume Method ◽  
Water Wave ◽  
Decomposition Method ◽  
Wave Equations ◽  
Kdv Equation ◽  
The Third ◽  
Laplace Decomposition Method ◽  
Volume Method

Abstract: In this paper we consider three methods of approximation for the nonlinear water wave equation. In particular we are interested of KdV equation as a stationary water wave. The first is the method of approximation with a polynomial, the second method is the finite–volume method and the third method is Laplace decomposition method (LDM). A comparison between the methods is mentioned in this article. We treat the considered methods comparing the obtained solutions with the exact ones. We give in particular the numerical results compared with the analytical results. We show that the used methods are effective and convenient for solving the water wave equations. We can propose and sure that the method of approximation with a polynomial gives accurate results.

Extended F-Expansion Method for Solving the modified Korteweg-de Vries (mKdV) Equation

2020 ◽  
Vol 11 (1) ◽  
pp. 93-100
Author(s):  
Vina Apriliani ◽  
Ikhsan Maulidi ◽  
Budi Azhari
Keyword(s):  
Wave Equation ◽  
Exact Solutions ◽  
Internal Waves ◽  
Water Waves ◽  
Wave Equations ◽  
Elliptic Functions ◽  
Expansion Method ◽  
Mkdv Equation ◽  
Innovative Methods ◽  
De Vries

One of the phenomenon in marine science that is often encountered is the phenomenon of water waves. Waves that occur below the surface of seawater are called internal waves. One of the mathematical models that can represent solitary internal waves is the modified Korteweg-de Vries (mKdV) equation. Many methods can be used to construct the solution of the mKdV wave equation, one of which is the extended F-expansion method. The purpose of this study is to determine the solution of the mKdV wave equation using the extended F-expansion method. The result of solving the mKdV wave equation is the exact solutions. The exact solutions of the mKdV wave equation are expressed in the Jacobi elliptic functions, trigonometric functions, and hyperbolic functions. From this research, it is expected to be able to add insight and knowledge about the implementation of the innovative methods for solving wave equations. 

scholarly journals Fractional wave equations with attenuation

2013 ◽  
Vol 16 (1) ◽  
Author(s):  
Peter Straka ◽  
Mark Meerschaert ◽  
Robert McGough ◽  
Yuzhen Zhou
Keyword(s):  
Wave Equation ◽  
Power Law ◽  
Weak Solutions ◽  
Wave Equations ◽  
Inhomogeneous Media ◽  
Medical Ultrasound ◽  
Sound Waves ◽  
Fractional Wave Equations

AbstractFractional wave equations with attenuation have been proposed by Caputo [5], Szabo [28], Chen and Holm [7], and Kelly et al. [11]. These equations capture the power-law attenuation with frequency observed in many experimental settings when sound waves travel through inhomogeneous media. In particular, these models are useful for medical ultrasound. This paper develops stochastic solutions and weak solutions to the power law wave equation of Kelly et al. [11].

On the existence of compressional MHD oscillations in an inhomogeneous magnetoplasma

1990 ◽  
Vol 44 (2) ◽  
pp. 361-375 ◽  
Author(s):  
Andrew N. Wright
Keyword(s):  
Magnetic Field ◽  
Wave Equation ◽  
Density Distribution ◽  
Cold Plasma ◽  
Wave Equations ◽  
Perturbation Field ◽  
Plasma Density Distribution ◽  
Background Magnetic Field ◽  
Transverse Fields ◽  
The Magnetic Field

In a cold plasma the wave equation for solely compressional magnetic field perturbations appears to decouple in any surface orthogonal to the background magnetic field. However, the compressional fields in any two of these surfaces are related to each other by the condition that the perturbation field b be divergence-free. Hence the wave equations in these surfaces are not truly decoupled from one another. If the two solutions happen to be ‘matched’ (i.e. V.b = 0) then the medium may execute a solely compressional oscillation. If the two solutions are unmatched then transverse fields must evolve. We consider two classes of compressional solutions and derive a set of criteria for when the medium will be able to support pure compressional field oscillations. These criteria relate to the geometry of the magnetic field and the plasma density distribution. We present the conditions in such a manner that it is easy to see if a given magnetoplasma is able to executive either of the compressional solutions we investigate.

Imaging Cracks in Bones Using Acoustic Nonlinearity: A Simulation Study

2011 ◽  
Vol 97 (5) ◽  
pp. 728-733
Author(s):  
Yang Liu ◽  
Xiasheng Guo ◽  
Zhao Da ◽  
Dong Zhang ◽  
Xiufen Gong
Keyword(s):  
Time Reversal ◽  
Three Dimensional ◽  
Harmonic Component ◽  
Ultrasonic Signal ◽  
Long Bone ◽  
Third Harmonic ◽  
Acoustic Nonlinearity ◽  
The Third ◽  
Nonlinear Approach ◽  
Radial Depth

This article proposes an acoustic nonlinear approach combined with the time reversal technique to image cracks in long bones. In this method, the scattered ultrasound generated from the crack is recorded, and the third harmonic nonlinear component of the ultrasonic signal is used to reconstruct an image of the crack by the time reversal process. Numerical simulations are performed to examine the validity of this approach. The fatigue long bone is modeled as a hollow cylinder with a crack of 1, 0.1, and 0.225 mm in axial, radial and circumferential directions respectively. A broadband 500 kHz ultrasonic signal is used as the exciting signal, and the extended three-dimensional Preisach-Mayergoyz model is used to describe the nonclassical nonlinear dynamics of the crack. Time reversal is carried out by using the filtered third harmonic component. The localization capability depends on the radial depth of the crack.

Pure P- and S-Wave Equations in Anisotropic Media

2020 ◽  
Author(s):  
A. Stovas ◽  
T. Alkhalifah ◽  
U. Bin Waheed
Keyword(s):  
Wave Equations ◽  
Anisotropic Media ◽  
S Wave
Sign in / Sign up

Export Citation Format

Share Document