Fourth‐order finite‐difference P-SV seismograms

Geophysics ◽  
1988 ◽  
Vol 53 (11) ◽  
pp. 1425-1436 ◽  
Author(s):  
Alan R. Levander
Keyword(s):  
Free Surface ◽  
Finite Difference ◽  
Fourth Order ◽  
Equations Of Motion ◽  
Dispersion Analysis ◽  
Numerical Dispersion ◽  
Staggered Grid ◽  
Order System ◽  
Surface Boundary ◽  
Lamb’S Problem

I describe the properties of a fourth‐order accurate space, second‐order accurate time, two‐dimensional P-SV finite‐difference scheme based on the Madariaga‐Virieux staggered‐grid formulation. The numerical scheme is developed from the first‐order system of hyperbolic elastic equations of motion and constitutive laws expressed in particle velocities and stresses. The Madariaga‐Virieux staggered‐grid scheme has the desirable quality that it can correctly model any variation in material properties, including both large and small Poisson’s ratio materials, with minimal numerical dispersion and numerical anisotropy. Dispersion analysis indicates that the shortest wavelengths in the model need to be sampled at 5 gridpoints/wavelength. The scheme can be used to accurately simulate wave propagation in mixed acoustic‐elastic media, making it ideal for modeling marine problems. Explicitly calculating both velocities and stresses makes it relatively simple to initiate a source at the free‐surface or within a layer and to satisfy free‐surface boundary conditions. Benchmark comparisons of finite‐difference and analytical solutions to Lamb’s problem are almost identical, as are comparisons of finite‐difference and reflectivity solutions for elastic‐elastic and acoustic‐elastic layered models.

Lowrank seismic-wave extrapolation on a staggered grid

Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. T157-T168 ◽  
Author(s):  
Gang Fang ◽  
Sergey Fomel ◽  
Qizhen Du ◽  
Jingwei Hu
Keyword(s):  
Seismic Wave ◽  
Wave Equations ◽  
Variable Density ◽  
Dispersion Analysis ◽  
Numerical Dispersion ◽  
Staggered Grid ◽  
Order System ◽  
Constant Density ◽  
Reverse Time ◽  
Time Migration

We evaluated a new spectral method and a new finite-difference (FD) method for seismic-wave extrapolation in time. Using staggered temporal and spatial grids, we derived a wave-extrapolation operator using a lowrank decomposition for a first-order system of wave equations and designed the corresponding FD scheme. The proposed methods extend previously proposed lowrank and lowrank FD wave extrapolation methods from the cases of constant density to those of variable density. Dispersion analysis demonstrated that the proposed methods have high accuracy for a wide wavenumber range and significantly reduce the numerical dispersion. The method of manufactured solutions coupled with mesh refinement was used to verify each method and to compare numerical errors. Tests on 2D synthetic examples demonstrated that the proposed method is highly accurate and stable. The proposed methods can be used for seismic modeling or reverse-time migration.

Accuracy of heterogeneous staggered-grid finite-difference modeling of Rayleigh waves

Geophysics ◽  
2006 ◽  
Vol 71 (4) ◽  
pp. T109-T115 ◽  
Author(s):  
Thomas Bohlen ◽  
Erik H. Saenger
Keyword(s):  
Free Surface ◽  
Finite Difference ◽  
Surface Topography ◽  
Rayleigh Waves ◽  
Complex Surface ◽  
Numerical Dispersion ◽  
Staggered Grid ◽  
Image Method ◽  
The Earth ◽  
Grid Points

Heterogeneous finite-difference (FD) modeling assumes that the boundary conditions of the elastic wavefield between material discontinuities are implicitly fulfilled by the distribution of the elastic parameters on the numerical grid. It is widely applied to weak elastic contrasts between geologic formations inside the earth. We test the accuracy at the free surface of the earth. The accuracy for modeling Rayleigh waves using the conventional standard staggered-grid (SSG) and the rotated staggered grid (RSG) is investigated. The accuracy tests reveal that one cannot rely on conventional numerical dispersion discretization criteria. A higher sampling is necessary to obtain acceptable accuracy. In the case of planar free surfaces aligned with the grid, 15 to 30 grid points per minimum wavelength of the Rayleigh wave are required. The widely used explicit boundary condition, the so-called image method, produces similar accuracy and requires approximately half the sampling of the wavefield compared to heterogeneous free-surface modeling. For a free-surface not aligned with the grid (surface topography), the error increases significantly and varies with the dip angle of the interface. For an irregular interface, the RSG scheme is more accurate than the SSG scheme. The RSG scheme, however, requires 60 grid points per minimum wavelength to achieve good accuracy for all dip angles. The high computation requirements for 3D simulations on such fine grids limit the application of heterogenous modeling in the presence of complex surface topography.

3D finite-difference modeling of elastic wave propagation in the Laplace-Fourier domain

Geophysics ◽  
2012 ◽  
Vol 77 (4) ◽  
pp. T137-T155 ◽  
Author(s):  
Petr V. Petrov ◽  
Gregory A. Newman
Keyword(s):  
Wave Propagation ◽  
Finite Difference ◽  
Fourth Order ◽  
Order Approximation ◽  
Staggered Grid ◽  
Order System ◽  
Finite Difference Schemes ◽  
Fourier Domain ◽  
Krylov Methods ◽  
Complex Objects

With the recent interest in the Laplace-Fourier domain full waveform inversion, we have developed new heterogeneous 3D fourth- and second-order staggered-grid finite-difference schemes for modeling seismic wave propagation in the Laplace-Fourier domain. Our approach is based on the integro-interpolation technique for the velocity-stress formulation in the Cartesian coordinate system. Five averaging elastic coefficients and three averaging densities are necessary to describe the heterogeneous medium, with harmonic averaging of the bulk and shear moduli, and arithmetic averaging of density. In the fourth-order approximation, we improved the accuracy of the scheme using a combination of integral identities for two elementary volumes — “small” and “large” around spatial gridpoints where the wave variables are defined. Two solution approaches are provided, both of which are solved with iterative Krylov methods. In the first approach, the stress variables are eliminated and a linear system for the velocity components is solved. In the second approach, we worked directly with the first-order system of velocity and stress variables. This reduced the computer memory required to store the complex matrix, along with reducing (by 30%) the number of arithmetic operations needed for the solution compared to the fourth-order scheme for velocity only. Numerical examples show that our finite-difference formulations for elastic wavefield simulations can achieve more accurate solutions with fewer grid points than those from previously published second and fourth-order frequency-domain schemes. We applied our simulator to the investigation of wavefields from the SEG/EAGE model in the Laplace-Fourier domain. The calculation is sensitive to the heterogeneity of the medium and capable of describing the structures of complex objects. Our technique can also be extended to 3D elastic modeling within the time domain.

Finite‐difference modeling of wave propagation in a fluid–solid configuration

Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 618-624 ◽  
Author(s):  
Robbert van Vossen ◽  
Johan O. A. Robertsson ◽  
Chris H. Chapman
Keyword(s):  
Wave Propagation ◽  
Finite Difference ◽  
Fourth Order ◽  
Spatial Sampling ◽  
Numerical Dispersion ◽  
Staggered Grid ◽  
Imaging Method ◽  
Complex Structures ◽  
Scholte Waves ◽  
Scholte Wave

Finite‐difference (FD) techniques are widely used to model wave propagation through complex structures. Two main sources of error can be identified: (1) from numerical dispersion and numerical anisotropy and (2) by modeling the response of internal grid boundaries. Conventional discretization criteria to reduce the effects of numerical dispersion and numerical anisotropy have long been established (5‐8 gridpoints per wavelength for a fourth‐order accurate FD scheme). We analyze the second source of errors, comparing different staggered‐grid FD solutions to the Cagniard‐de Hoop solution in models with fluid–solid contacts. Our results confirm that it is sufficient to rely on conventional discretization criteria if the fluid–solid interface is aligned with the grid. If accurate modeling of the Scholte wave is required, then a new imaging method we propose should be used to allow for conventional sampling of the wavefield to minimize numerical dispersion. However, for an interface not aligned with the grid (irregular interfaces), a spatial sampling of at least 15 gridpoints per minimum wavelength is required to obtain acceptable results, particularly in seismic seabed applications where Scholte waves may need to be modeled more accurately.

Acoustic and elastic modeling by optimal time-space-domain staggered-grid finite-difference schemes

Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. T17-T40 ◽  
Author(s):  
Zhiming Ren ◽  
Yang Liu
Keyword(s):  
Finite Difference ◽  
Stability Condition ◽  
Computational Cost ◽  
Numerical Dispersion ◽  
Staggered Grid ◽  
Equation Modeling ◽  
Optimal Time ◽  
Finite Difference Schemes ◽  
Space Domain ◽  
Time Space

Staggered-grid finite-difference (SFD) methods are widely used in modeling seismic-wave propagation, and the coefficients of finite-difference (FD) operators can be estimated by minimizing dispersion errors using Taylor-series expansion (TE) or optimization. We developed novel optimal time-space-domain SFD schemes for acoustic- and elastic-wave-equation modeling. In our schemes, a fourth-order multiextreme value objective function with respect to FD coefficients was involved. To yield the globally optimal solution with low computational cost, we first used variable substitution to turn our optimization problem into a quadratic convex one and then used least-squares (LS) to derive the optimal SFD coefficients by minimizing the relative error of time-space-domain dispersion relations over a given frequency range. To ensure the robustness of our schemes, a constraint condition was imposed that the dispersion error at each frequency point did not exceed a given threshold. Moreover, the hybrid absorbing boundary condition was applied to remove artificial boundary reflections. We compared our optimal SFD with the conventional, TE-based time-space-domain, and LS-based SFD schemes. Dispersion analysis and numerical simulation results suggested that the new SFD schemes had a smaller numerical dispersion than the other three schemes when the same operator lengths were adopted. In addition, our LS-based time-space-domain SFD can obtain the same modeling accuracy with shorter spatial operator lengths. We also derived the stability condition of our schemes. The experiment results revealed that our new LS-based SFD schemes needed a slightly stricter stability condition.

Sign in / Sign up

Export Citation Format

Share Document