Fourier finite-difference wave propagation

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. T123-T129 ◽  
Author(s):  
Xiaolei Song ◽  
Sergey Fomel
Keyword(s):  
Finite Difference ◽  
Acoustic Wave Equation ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Complex Velocity ◽  
Difference Wave ◽  
Time Migration ◽  
A Chain ◽  
Anisotropic Case ◽  
Accuracy And Stability

We introduce a novel technique for seismic wave extrapolation in time. The technique involves cascading a Fourier transform operator and a finite-difference operator to form a chain operator: Fourier finite differences (FFD). We derive the FFD operator from a pseudoanalytical solution of the acoustic wave equation. Two-dimensional synthetic examples demonstrate that the FFD operator can have high accuracy and stability in complex-velocity media. Applying the FFD method to the anisotropic case overcomes some disadvantages of other methods, such as the coupling of qP-waves and qSV-waves. The FFD method can be applied to enhance accuracy and stability of seismic imaging by reverse time migration.

Implicit interpolation in reverse‐time migration

Geophysics ◽  
1997 ◽  
Vol 62 (3) ◽  
pp. 906-917 ◽  
Author(s):  
Jinming Zhu ◽  
Larry R. Lines
Keyword(s):  
Wave Equation ◽  
Finite Difference ◽  
Real Data ◽  
Seismic Sources ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Difference Wave ◽  
Time Migration ◽  
Seismic Acquisition ◽  
Recorded Data

Reverse‐time migration applies finite‐difference wave equation solutions by using unaliased time‐reversed recorded traces as seismic sources. Recorded data can be sparsely or irregularly sampled relative to a finely spaced finite‐difference mesh because of the nature of seismic acquisition. Fortunately, reliable interpolation of missing traces is implicitly included in the reverse‐time wave equation computations. This implicit interpolation is essentially based on the ability of the wavefield to “heal itself” during propagation. Both synthetic and real data examples demonstrate that reverse‐time migration can often be performed effectively without the need for explicit interpolation of missing traces.

An optimized wave equation for seismic modeling and reverse time migration

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCA153-WCA158 ◽  
Author(s):  
Faqi Liu ◽  
Guanquan Zhang ◽  
Scott A. Morton ◽  
Jacques P. Leveille
Keyword(s):  
Wave Equation ◽  
Finite Difference ◽  
Acoustic Wave ◽  
Computational Cost ◽  
Numerical Dispersion ◽  
New Wave ◽  
Acoustic Wave Equation ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration

The acoustic wave equation has been widely used for the modeling and reverse time migration of seismic data. Numerical implementation of this equation via finite-difference techniques has established itself as a valuable approach and has long been a favored choice in the industry. To ensure quality results, accurate approximations are required for spatial and time derivatives. Traditionally, they are achieved numerically by using either relatively very fine computation grids or very long finite-difference operators. Otherwise, the numerical error, known as numerical dispersion, is present in the data and contaminates the signals. However, either approach will result in a considerable increase in the computational cost. A simple and computationally low-cost modification to the standard acoustic wave equation is presented to suppress numerical dispersion. This dispersion attenuator is one analogy of the antialiasing operator widely applied in Kirchhoff migration. When the new wave equation is solved numerically using finite-difference schemes, numerical dispersion in the original wave equation is attenuated significantly, leading to a much more accurate finite-difference scheme with little additional computational cost. Numerical tests on both synthetic and field data sets in both two and three dimensions demonstrate that the optimized wave equation dramatically improves the image quality by successfully attenuating dispersive noise. The adaptive application of this new wave equation only increases the computational cost slightly.

Angle-Domain Gathers Computation Using Poynting Vector of Two-Way Acoustic Wave Equation

2014 ◽  
Vol 962-965 ◽  
pp. 2984-2987
Author(s):  
Jia Jia Yang ◽  
Bing Shou He ◽  
Ting Chen
Keyword(s):  
Wave Equation ◽  
Acoustic Wave ◽  
Real Data ◽  
Flow Density ◽  
Acoustic Wave Equation ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Wave Fields ◽  
Time Migration ◽  
Source Wave

Based on two-way acoustic wave equation, we present a method for computing angle-domain common-image gathers for reverse time migration. The method calculates the propagation direction of source wave-fields and receiver wave-fields according to expression of energy flow density vectors (Poynting vectors) of acoustic wave equation in space-time domain to obtain the reflection angle, then apply the normalized cross-correlation imaging condition to achieve the angle-domain common-image gathers. The angle gathers obtained can be used for migration velocity analysis, AVA analysis and so on. Numerical examples and real data examples demonstrate the effectiveness of this method.

scholarly journals Local Explosion Detection and Infrasound Localization by Reverse Time Migration Using 3-D Finite-Difference Wave Propagation

2021 ◽  
Vol 9 ◽  
Author(s):  
David Fee ◽  
Liam Toney ◽  
Keehoon Kim ◽  
Richard W. Sanderson ◽  
Alexandra M. Iezzi ◽  
...  
Keyword(s):  
Numerical Modeling ◽  
Finite Difference ◽  
Grid Point ◽  
Reverse Time ◽  
Computationally Efficient ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Straight Line ◽  
Active Vent ◽  
Time Reversal Mirror

Infrasound data are routinely used to detect and locate volcanic and other explosions, using both arrays and single sensor networks. However, at local distances (<15 km) topography often complicates acoustic propagation, resulting in inaccurate acoustic travel times leading to biased source locations when assuming straight-line propagation. Here we present a new method, termed Reverse Time Migration-Finite-Difference Time Domain (RTM-FDTD), that integrates numerical modeling into the standard RTM back-projection process. Travel time information is computed across the entire potential source grid via FDTD modeling to incorporate the effects of topography. The waveforms are then back-projected and stacked at each grid point, with the stack maximum corresponding to the likely source. We apply our method to three volcanoes with different network configurations, source-receiver distances, and topography. At Yasur Volcano, Vanuatu, RTM-FDTD locates explosions within ∼20 m of the source and differentiates between multiple vents. RTM-FDTD produces a more accurate location for the two Yasur subcraters than standard RTM and doubles the number of detected events. At Sakurajima Volcano, Japan, RTM-FDTD locates the source within 50 m of the active vent despite notable topographic blocking. The RTM-FDTD location is similar to that from the Time Reversal Mirror method, but is more computationally efficient. Lastly, at Shishaldin Volcano, Alaska, RTM and RTM-FDTD both produce realistic source locations (<50 m) for ground-coupled airwaves recorded on a four-station seismic network. We show that RTM is an effective method to detect and locate infrasonic sources across a variety of scenarios, and by integrating numerical modeling, RTM-FDTD produces more accurate source locations and increases the detection capability.

Sign in / Sign up

Export Citation Format

Share Document