next up previous [pdf]

Next: Applications Up: 2D linear phase approximation Previous: Directional separability

Frequency responses

We show the magnitude responses of the infinite impulse response (IIR) omnidirectional PWD (OPWD) $ 1-H_2(Z_1,Z_2)$ when $ N=2$ in Figure 4d-f. Compared with the ideal responses in Figure 2, it has good approximations for all the three dip angles in most of the frequency band. There are distortions in the high frequency bands, due to phase approximation errors in $ H_2(Z_1,Z_2)$ for high frequencies. They are not significant in practice, because the frequency band of seismic data is temporally limited.

liir20 liir50 liir80 oiir20 oiir50 oiir80
liir20,liir50,liir80,oiir20,oiir50,oiir80
Figure 4.
Magnitude responses of the IIR implementation of LPWD $ 1-Z_2H_1(Z_1)$ (top) and OPWD $ 1-H_2(Z_1,Z_2)$ (bottom): from left to right $ \theta =20^\circ ,50^\circ ,80^\circ $ .
[pdf] [pdf] [pdf] [pdf] [pdf] [pdf] [png] [png] [png] [png] [png] [png] [scons]

When $ \theta=20^\circ$ or $ 50^\circ$ , the omnidirectional PWD has a low-frequency response similar to that of line-interpolating PWD (LPWD) shown in Figure 4a-b. When $ \theta=80^\circ$ , the LPWD has aliasing in low-frequency bands (Figure 4c), while the OPWD exhibits a more desirable low-frequency response (Figure 4f).

Following Fomel (2002), the IIR LPWD $ 1-Z_2H_1(Z_1)$ can be approximated by the following finite impulse response (FIR) filter:

$\displaystyle H_3(Z_1,Z_2)=B(Z_1)-Z_2B(1/Z_1).$ (16)

Similarly, the IIR OPWD filter $ 1-H_2(Z_1,Z_2)$ can also be approximated by an FIR implementation:

$\displaystyle H_4(Z_1,Z_2)=F(Z_1,Z_2)-F(1/Z_1,1/Z_2).$ (17)

lfir20 lfir50 lfir80 ofir20 ofir50 ofir80
lfir20,lfir50,lfir80,ofir20,ofir50,ofir80
Figure 5.
Magnitude response of the FIR approximation of LPWD $ H_3(Z_1,Z_2)$ (top) and OPWD $ H_4(Z_1,Z_2)$ (bottom): from left to right $ \theta =20^\circ ,50^\circ ,80^\circ $ .
[pdf] [pdf] [pdf] [pdf] [pdf] [pdf] [png] [png] [png] [png] [png] [png] [scons]

The frequency responses of these two FIR approximations are shown in Figure 5. Similar to the IIR implementation, the FIR OPWD (Figure 5d-f) can obtain an expected low-frequency response for all three dip angles. While the FIR LPWD cannot obtain a good low-frequency response when $ \theta=80^\circ$ (Figure 5c). Compared with the IIR implementations, the FIR approximations of both LPWD and OPWD have less desirable high-frequency responses.


next up previous [pdf]

Next: Applications Up: 2D linear phase approximation Previous: Directional separability

2013-08-09