previous up next print clean
Next: Estimating the spatial prediction Up: ESTIMATION OF SPATIAL PREDICTION Previous: Wavefield composed of one

Wavefield composed of two plane waves

If we have two plane waves, the wavefield (Figure 6) can be expressed by
\begin{displaymath}
w(t,x) = a_1\delta(t-t_1-p_1x)+a_2\delta(t-t_2-p_2x)\end{displaymath} (11)

 
two-plane-wave
Figure 6
A synthetic wavefield composed of two plane waves.
two-plane-wave
view burn build edit restore

In the f-x domain, it has another expression
\begin{displaymath}
W(f,x) = a_1e^{i2\pi f(t_1+p_1x)}+a_2e^{i2\pi f(t_2+p_2x)}\end{displaymath} (12)
Or

 
W(f,x) = E1+E2 (13)

where $E_1 = a_1e^{i2\pi f(t_1+p_1x)}$ and $E_2=a_2e^{i2\pi f(t_2+p_2x)}$.The adjacent trace can be expressed by  
 \begin{displaymath}
W(f,x-\Delta x)=P^{-1}_1E_1+P^{-1}_2E_2\end{displaymath} (14)
Now we cannot predict one trace from another trace, since there are two unknown propagators. Consequently, we need the information from another trace.  
 \begin{displaymath}
W(f,x-2\Delta x)=P^{-2}_1E_1+P^{-2}_2E_2\end{displaymath} (15)
Equation (14) and (15) can be solved for the propagators. But in reality, we do not know p1 and p2. So we cannot get the explicit form of the propagators. We must use another form of parameter to link the the different traces together, which is called a spatial prediction filter.

If we have trace $x-\Delta x$ and trace $x-2\Delta x$, our aim is to predict trace x. Assume there is a relation between the three traces:
\begin{displaymath}
W(f,x) = C_1W(f,x-\Delta x)+C_2W(f,x-2\Delta x)\end{displaymath} (16)
Substituting Ws with equations (13), (14) and (15), we get

 
E1+E2 = C1(P-11E1+P-12E2)+C2(P-21E1+P-22E2) (17)

$E_{\rm i}$ can be any value, as long as it is a plane wave. Therefore, equation (17) holds for all of E1 and E2. In particular, we have
\begin{displaymath}
1 = C_1P^{-1}_1+C_2P^{-2}_1, \hspace*{0.2in}(E_1 = 1, \hspace*{0.1in}E_2 = 0)\end{displaymath} (18)
\begin{displaymath}
1 = C_1P^{-1}_2+C_2P^{-2}_2, \hspace*{0.2in}(E_1 = 0, \hspace*{0.1in}E_2 = 1)\end{displaymath} (19)
Solving this linear equation, we find the relation between $C_{\rm i}$ and $P_{\rm i}$ (i=1,2). That is
\begin{displaymath}
C_1 = P_1 + P_2, \hspace*{0.3in} C_2 = -P_1P_2\end{displaymath} (20)

Like the one plane wave case, all above deductions have nothing to do with the gridding scheme. All the results are effective for both cases. The differences between the two gridding schemes are

We can extend our deduction to N plane waves (see the Appendix) and we can get similar results: for FIG, coefficient $C_{\rm i}$ is frequency-dependent; for FDG, $C_{\rm i}$ is frequency-independent. This is the basis of estimating the spatial prediction filter.


previous up next print clean
Next: Estimating the spatial prediction Up: ESTIMATION OF SPATIAL PREDICTION Previous: Wavefield composed of one
Stanford Exploration Project
11/11/1997