2-D elastic wave forward modeling with frequency-space 25-point finite-difference operators

2018 ◽  
Author(s):  
Hongliang Li ◽  
Shoudong Wang ◽  
Doudou Wang ◽  
Mingxiao Cui ◽  
Kailong Su
2009 ◽  
Vol 6 (3) ◽  
pp. 259-266 ◽  
Author(s):  
Jianping Liao ◽  
Huazhong Wang ◽  
Zaitian Ma

Geophysics ◽  
1992 ◽  
Vol 57 (2) ◽  
pp. 218-232 ◽  
Author(s):  
A. Vafidis ◽  
F. Abramovici ◽  
E. R. Kanasewich

Two finite‐difference schemes for solving the elastic wave equation in heterogeneous two‐dimensional media are implemented on a vector computer. A modified Lax‐Wendroff scheme that is second‐order accurate both in time and space and is a version of the MacCormack scheme that is second‐order accurate in time and fourth‐order in space. The algorithms are based on the matrix times vector by diagonals technique that is fully vectorized and is described using a novel notation for vector supercomputer operations. The technique described can be implemented on a vector processor of modest dimensions and increase the applicability of finite differences. The two difference operators are compared and the programs are tested for a simple case of standing sinusoidal waves for which the exact solution is known and also for a two‐layer model with a line source. A comparison of the results for an actual well‐to‐well experiment verifies the usefulness of the two‐dimensional approach in modeling the results.


2014 ◽  
Vol 2014 ◽  
pp. 1-14 ◽  
Author(s):  
Ke-Yang Chen

Elastic wave equation simulation offers a way to study the wave propagation when creating seismic data. We implement an equivalent dual elastic wave separation equation to simulate the velocity, pressure, divergence, and curl fields in pure P- and S-modes, and apply it in full elastic wave numerical simulation. We give the complete derivations of explicit high-order staggered-grid finite-difference operators, stability condition, dispersion relation, and perfectly matched layer (PML) absorbing boundary condition, and present the resulting discretized formulas for the proposed elastic wave equation. The final numerical results of pure P- and S-modes are completely separated. Storage and computing time requirements are strongly reduced compared to the previous works. Numerical testing is used further to demonstrate the performance of the presented method.


Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. T293-T308
Author(s):  
Longfei Gao ◽  
David Keyes

We present a mechanism to explicitly couple the finite-difference discretizations of 2D acoustic and isotropic elastic-wave systems that are separated by straight interfaces. Such coupled simulations allow for the application of the elastic model to geological regions that are of special interest for seismic exploration studies (e.g., the areas surrounding salt bodies), with the computationally more tractable acoustic model still being applied in the background regions. Specifically, the acoustic wave system is expressed in terms of velocity and pressure while the elastic wave system is expressed in terms of velocity and stress. Both systems are posed in first-order forms and are discretized on staggered grids. Special variants of the standard finite-difference operators, namely, operators that possess the summation-by-parts property, are used for the approximation of spatial derivatives. Penalty terms, which are also referred to as the simultaneous approximation terms, are designed to weakly impose the elastic-acoustic interface conditions in the finite-difference discretizations and couple the elastic and acoustic wave simulations together. With the presented mechanism, we are able to perform the coupled elastic-acoustic wave simulations stably and accurately. Moreover, it is shown that the energy-conserving property in the continuous systems can be preserved in the discretized systems with carefully designed penalty terms.


Sign in / Sign up

Export Citation Format

Share Document