1. Wave equation and the finite difference method - 2023.1 English

Vitis Libraries

Release Date
2023-12-20
Version
2023.1 English

Equation (1) is the basic mathematical equation for elastic wave propagation.

\[\frac{1}{v(\mathbf{x})^2}\frac{\partial^2p(\mathbf{x}, t)}{\partial t^2} = \nabla^2 p(\mathbf{x}, t)\hspace{50mm}(1)\]

where \(\nabla^2\) is the Laplace operator, \(p(\mathbf{x}, t)\) is the wave field and \(v(\mathbf{x})\) is the wave velocity. In 2D-RTM, \(\mathbf{x} = \{x, z\}\), which is the length and depth of the imaging area respectively. In 3D-RTM, \(\mathbf{x} = \{x, y, z\}\), which is the length, width and depth of the imaging area respectively.

Due to the non-uniform structures in the subsurface, there is no theoretical solution for Equation 1. Hence, RTM is based on a numerical solution called a finite difference method. All spacial dimensions and time domains are discretized properly. The partial derivative with respect to time and the Laplace operator also have corresponding transformation according to the discretization scheme.

For instance, the Pluto synthetic seismic velocity model which was used to evaluation our 2D-RTM implementation has been discretized into a regular grid both spatially and temporally. The number of partiitions are listed in Table 1.

Table 124 Table 1: Pluto model discretization
Parameter Length (\(N_{x}\)) Depth (\(N_{z}\)) Time (\(N_{t}\))
No. grids 6960 1201 12860

where \(N_x\), \(N_z\), \(N_t\) denote the number of grids along \(x\), \(z\), \(t\) respectively.

Similar to the approaches presented in many research publications, we also adopted Order-2 temporal scheme and Order-8 spatial scheme. Following these discretization schemes, the 2D wave equation can be re-interpreted by Equation (2) and (3) shown below. Equation (2*) and (3) are for 3D RTM.

\[\nabla^2p = \sum_{n=-4}^4 c_n(\frac{p(x+n\Delta x)}{\Delta x^2} + \frac{p(z+n\Delta z)}{\Delta z^2})\hspace{20mm}(2)\]
\[\nabla^2p = \sum_{n=-4}^4 c_n(\frac{p(x+n\Delta x)}{\Delta x^2} + \frac{p(y+n\Delta y)}{\Delta y^2}+ \frac{p(z+n\Delta z)}{\Delta z^2})\hspace{20mm}(2*)\]
\[p(t+1) = 2p(t) - p(t-1) + v^2 \Delta t^2 \nabla^2p \label{eq:temperal}\hspace{30mm}(3)\]

Where all \(c_n\) are the coeffients shown in Table 2 below.

Table 125 Table 2: Finete Difference Method Order-8 Coefficients
\(c_{4}\), \(c_{-4}\) \(c_{3}\), \(c_{-3}\) \(c_{2}\), \(c_{-2}\) \(c_{1}\), \(c_{-1}\) \(c_{0}\)
\(-\frac{1}{560}\) \(\frac{8}{315}\) \(-\frac{1}{5}\) \(\frac{8}{5}\) \(-\frac{1}{560}\)