Inverse filtering

From SEG Wiki
Jump to: navigation, search

If a filter operator f(t) were defined such that convolution of f(t) with the known seismogram x(t) yields an estimate of the earth’s impulse response e(t), then


(4)

By substituting equation (4) into equation (3a), we get


(3a)
(5)

When x(t) is eliminated from both sides of the equation, the following expression results:


(6)

where δ(t) represents the Kronecker delta function:


(7)

By solving equation (6) for the filter operator f(t), we obtain


(8)

Therefore, the filter operator f(t) needed to compute the earth’s impulse response from the recorded seismogram turns out to be the mathematical inverse of the seismic wavelet w(t). Equation (8) implies that the inverse filter converts the basic wavelet to a spike at t = 0. Likewise, the inverse filter converts the seismogram to a series of spikes that defines the earth’s impulse response. Therefore, inverse filtering is a method of deconvolution, provided the source waveform is known (deterministic deconvolution). The procedure for inverse filtering is described in Figure 2.2-1.

Figure 2.2-1  A flowchart for inverse filtering.

See also

External links

find literature about
Inverse filtering
SEG button search.png Datapages button.png GeoScienceWorld button.png OnePetro button.png Schlumberger button.png Google button.png AGI button.png