next up previous print clean
Next: Finite-difference datuming Up: Phase-shift datuming Previous: Phase-shift datuming

Phase-shift datuming with lateral velocity variation

For laterally-variant media, phase-shift datuming can be performed by using an adaptation of Gazdag and Sguazzero's (1984) phase shift plus interpolation (PSPI) method. In this method, the wavefield is extrapolated at multiple velocities and a single upward-propagated wavefield is obtained by interpolation. Linear interpolation of two wavefields P1 and P2 is an operation of the type

\begin{displaymath}
\left[
 \begin{array}
{c}
 \tilde{P} \\  \end{array}\right]
...
 ...
\left[
 \begin{array}
{c}
 P_1 \\  P_2 \\  \end{array}\right],\end{displaymath}

where the adjoint operator spreads a value with two weights w1 and w2:

\begin{displaymath}
\left[
 \begin{array}
{c}
 \tilde{P}_1 \\  \tilde{P}_2 \\  \...
 ...ay}\right]
\left[
 \begin{array}
{c}
 P \\  \end{array}\right].\end{displaymath}

The PSPI upward datuming operator for the geometry of Figure [*] can be written as  
 \begin{displaymath}
\displaystyle{
\left[
 \begin{array}
{ccc}
U_3 U_2 U_1 & U_3...
 ...]
\left[
\begin{array}
{c}
P(x,z_s,\omega)\end{array}\right],
}\end{displaymath} (18)
where the matrices Ui represent the extrapolation operators for laterally-variant velocity. The matrix Ui can be further decomposed into the sequence  
 \begin{displaymath}
\displaystyle{
\left[
 \begin{array}
{cc}
 w_1 & w_2 \\  \en...
 ...}\right]
\left[
 \begin{array}
{c}
 F \\  \end{array}\right]
}.\end{displaymath} (19)
The physical interpretation of equation ([*]) is that the wavefield, after Fourier transformation, is split and upward continued with two different velocities. The two wavefields are independently inverse Fourier transformed and then interpolated. This sequence is repeated for each depth level.

The adjoint algorithm is found by transposing each matrix and reversing the multiplication order, as follows:  
 \begin{displaymath}
\displaystyle{
\left[
 \begin{array}
{ccc}
 A & B & C 
 \end...
 ...
\begin{array}
{c}
P(x,z_{\rm dat},\omega)\end{array}\right].
}\end{displaymath} (20)
The matrices U*i represent the operator for downward continuation of the wavefield to the depth level i. The matrix U*i is obtained by taking the adjoint of equation  ([*]):  
 \begin{displaymath}
\displaystyle{
\left[
 \begin{array}
{c}
 F^* \\  \end{array...
 ...left[
 \begin{array}
{c}
 w_1 \\  w_2 \\  \end{array}\right]
}.\end{displaymath} (21)

Popovici (1992) shows that the Split-Step formulation is similar, and that the only difference is in the Ui matrices.


next up previous print clean
Next: Finite-difference datuming Up: Phase-shift datuming Previous: Phase-shift datuming
Stanford Exploration Project
2/12/2001