Lowrank seismic wave extrapolation on a staggered grid |

We consider first-order acoustic wave equations for a medium of variable velocity and density. For a lossless 2-D medium,

where is acoustic particle velocity with components and ; is the acoustic pressure; is density of the medium; is seismic wave velocity of the medium; and denotes the space location in vector coordinate.

The second-order wave equation corresponding to equation 1 is

In the case of homogeneous velocity and density, equation 2 can be written in the spatial-frequency domain as

where is the 2-D spatial Fourier transform of . Equation 3 has an analytical solution,

Applying the second-order time-marching scheme leads to the k-space scheme (Tabei et al., 2002),

where sinc . In general, velocity and density vary in space. When both the gradient of velocity and the time step are small, replacing with in equation 5 provides a new approximation. Applying inverse Fourier transform to this approximation leads to the scheme

where denotes a spatial Fourier transform. The operator on the right-hand side of equation 6 depends on both and . Following Tabei et al. (2002) we call it second-order -space operator,

where the operator is analogous to the standard gradient operator, but it is a function of parameter . Similar to the definition of the standard gradient operator, we can define

Tabei et al. (2002) suggested a factorization, which can factor their second-order -space operator into parts associated with each spatial direction. This factorization can also be applied to -space operator by replacing constant velocity with variable velocity. The factored operators are called first-order -space operators:

The spatial frequency components and are defined so that . Application of the exponential coefficient in equation 9 requires the corresponding wavefield to be evaluated on grid points staggered by distance of along the positive or negative direction and along the positive or negative direction. The spatial staggering in equation 9 is implicitly incorporated into the spatial derivative by the shift property of the Fourier transform. Using operators in equation 9 within equation 1 enables a new construction of spectral method. The first-order coupled equations for acoustic-wave extrapolation in variable velocity and density media with staggered spatial and temporal grids are therefore:

where , , , .

The partial derivative operators in equation 10 are defined by equation 9. Note that the ordering of and is arbitrary depending on the configuration of staggered grid. However, these operators should be used in pairs, such that the spatial shifting cancel out over any temporal interval length . The spatial and temporal staggered grids used in equation 10 are analogous to staggered scheme employed in previous finite-difference methods(Virieux, 1986; Madariaga, 1976; Virieux, 1984), which are know to increase accuracy and stability by halving spatial and time interval without increasing the number of computational points. Equation 9 can be solved with a localized Fourier transform (Wards et al., 2008). However, this kind of solution has a high computational cost. Song et al. (2012) proposed to apply FFD method to calculate the first-order -space operators in equation 9, which can handle the variable velocity and density accurately and efficiently. Another possible ways to speed up the computation is to represent the extrapolation operator with a lowrank matrix.

Lowrank seismic wave extrapolation on a staggered grid |

2014-06-02