Next: Slope estimation
Up: Fomel: Planewave destructors
Previous: Introduction
The mathematical basis of the planewave destructor filters is the
local plane differential equation
 
(1) 
where P(t,x) is the wave field, and s is the local slope, which may
also depend on t and x. In the case of a constant slope,
equation (1) has the simple general solution
 
(2) 
where f(t) is an arbitrary waveform. Equation (2) is
nothing more than a mathematical description of a plane wave.
If the slope s does not depend on the t coordinate, we can
transform equation (1) to the frequency domain, where it
takes the form of the ordinary differential equation
 
(3) 
and has the general solution
 
(4) 
where is the Fourier transform of P. The complex
exponential term in equation (4) simply represents a shift
of a ttrace according to the slope s and the trace separation
x. In the frequency domain, the operator for transforming the trace
at position x1 to the neighboring trace at position x is a
multiplication by . In other words, a plane wave can be
perfectly predicted by a twoterm predictionerror filter in the
FX domain:
 
(5) 
where a_{0} = 1 and . The goal of predicting
several plane waves can be accomplished by cascading several twoterm
filters. In fact, any FX predictionerror filter, represented in
the Ztransform form as
 
(6) 
can be factored into a product of twoterm filters:
 
(7) 
where are the zeroes of
polynomial (6). According to equation (5),
the phase of each zero corresponds to the slope of a local plane wave
multiplied by the frequency. Zeroes that are not on the unit circle
carry an additional amplitude gain not included in
equation (3).
In order to incorporate timevarying slopes, we need to return back to
the time domain and look for an appropriate analog of the phaseshift
operator (4) and the planeprediction
filter (5). An important property of planewave
propagation across different traces is that the total energy of the
transmitted wave stays invariant throughout the process. This property
is assured in the frequencydomain solution (4) by the fact
that the spectrum of the complex exponential is
equal to one. In the time domain, we can reach an equivalent effect
by using an allpass digital filter. In the Ztransform notation,
convolution with an allpass filter takes the form
 
(8) 
where denotes the Ztransform of the corresponding
trace, and the ratio B(Z_{t})/B(1/Z_{t}) is an allpass digital filter,
approximating the timeshift operator (5). In
finitedifference terms, equation (8) represents an
implicit finitedifference scheme for solving equation (1)
with the initial conditions at a constant x. The coefficients of
filter B(Z_{t}) can be determined, for example, by fitting the filter
frequency response at small frequencies to the response of the
phaseshift operator. The Taylor series technique (equating the
coefficients of the Taylor series expansion around zero frequency)
yields the expression
 
(9) 
for a threepoint centered filter B_{3}(Z_{t}) and the expression
 

 
 (10) 
for a fivepoint centered filter B_{5}(Z_{t}). It is easy to generalize
these expressions to longer filters. Figure 1 shows the
phase of the allpass filters B_{3}(Z_{t})/B_{3}(1/Z_{t}) and
B_{5}(Z_{t})/B_{5}(1/Z_{t}) for two values of the slope s in
comparison with the exact linear function of equation (4).
As expected, the phases fit the exact line at low frequencies, and the
accuracy of the approximation increases with the length of the filter.
phase
Figure 1 Phase of the implicit
finitedifference shift operators in comparison with the exact
solution. Left plot corresponds to s=0.5. Right plot; s=0.8.
In two dimensions, equation (8) transforms to the
prediction equation analogous to (5) with the 2D
prediction filter^{}
 
(11) 
In order to characterize several plane waves, we can cascade several
filters of the form (11) in a manner similar to
equation (7). In all examples of this paper, I used a
modified version of the filter A(Z_{t},Z_{x}), namely the filter
 
(12) 
which avoids the need for polynomial division. In case of the 3point
filter (9), the 2D filter (12) has exactly
six coefficients, with the second t column being a reversed copy
of the first column. When filter (12) is used in
interpolation problems, it can occasionally cause undesired
highfrequency oscillations in the solution, resulting from the
nearNyquist zeroes of the polynomial B(Z_{t}). The oscillations are
easily removed in practice with an appropriate lowpass filtering.
In the next section, I address the problem of estimating the local
slope s with the filters of form (12). Estimating the
slope is a necessary step for applying the finitedifference
planewave filters on real data.
Next: Slope estimation
Up: Fomel: Planewave destructors
Previous: Introduction
Stanford Exploration Project
9/5/2000