Alias filters
|
| |
| Series | Geophysical References Series |
|---|---|
| Title | Problems in Exploration Seismology and their Solutions |
| Author | Lloyd P. Geldart and Robert E. Sheriff |
| Chapter | 9 |
| Pages | 295 - 366 |
| DOI | http://dx.doi.org/10.1190/1.9781560801733 |
| ISBN | ISBN 9781560801153 |
| Store | SEG Online Store |
Problem 9.5a
The standard alias filter such as that shown in Figure 9.5a has a 3-dB point at about half the Nyquist frequency $ f_{N} $ and a very steep slope, so that noise above Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): f_{N} is highly attenuated relative to the passband of the system. Assuming an original flat spectrum, alias filtering with a 125-Hz, 72-dB/octave filter, and subsequent resampling from 2 to 4 ms (without additional alias filtering), graph the resulting alias noise versus frequency.
Background
The passband of a system is the range of frequencies that are unattenuated by passage through the system. The limits are usually taken as the frequencies for which the attenuation is 3 dB.

Solution
With the 3-dB point on the high-cut filter at 125 Hz and a 72 dB/octave slope, the attenuation at Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): f_{N} =250 Hz is 75 dB. Frequencies Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \Delta{f} higher than the Nyquist frequency, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): f_{N} +\Delta f , alias to appear as the frequencies Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): f_{N} -\Delta{f} , i.e., they fold back about Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): f_{N} , as shown in Figure 9.5a.
A 65-Hz 72-dB/octave filter is usually applied before resampling to prevent aliasing. If resampling to 4 ms is done without this additional filtering, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): f_{N} for the resampled data is Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): 1/(2\times 0.004)=125 Hz; the aliased noise is shown by the foldback curves in Figure 9.5a.
Problem 9.5b
Some believe that standard alias filters may be unnecessarily restrictive. The standard alias filter for 4-ms sampling is about 65 Hz, 72-dB/octave. Graph the alias noise versus frequency for a 90-Hz, 72-dB/octave filter for 4-ms sampling and draw conclusions.
Solution
The 90-Hz filter is shown in Figure 9.5a. Because the sampling interval is 4 ms, $ f_{N} $ is 125 Hz and the alias noise is given by the foldback curves shown. The 90-Hz filter gives a broader passband than the standard 65-Hz alias filter if signal frequencies above about 80 Hz and signals attenuated more than 60 dB are not important.
Continue reading
| Previous section | Next section |
|---|---|
| Extension of the sampling theorem | The convolutional model |
| Previous chapter | Next chapter |
| Reflection field methods | Geologic interpretation of reflection data |
Also in this chapter
- Fourier series
- Space-domain convolution and vibroseis acquisition
- Fourier transforms of the unit impulse and boxcar
- Extension of the sampling theorem
- Alias filters
- The convolutional model
- Water reverberation filter
- Calculating crosscorrelation and autocorrelation
- Digital calculations
- Semblance
- Convolution and correlation calculations
- Properties of minimum-phase wavelets
- Phase of composite wavelets
- Tuning and waveshape
- Making a wavelet minimum-phase
- Zero-phase filtering of a minimum-phase wavelet
- Deconvolution methods
- Calculation of inverse filters
- Inverse filter to remove ghosting and recursive filtering
- Ghosting as a notch filter
- Autocorrelation
- Wiener (least-squares) inverse filters
- Interpreting stacking velocity
- Effect of local high-velocity body
- Apparent-velocity filtering
- Complex-trace analysis
- Kirchhoff migration
- Using an upward-traveling coordinate system
- Finite-difference migration
- Effect of migration on fault interpretation
- Derivative and integral operators
- Effects of normal-moveout (NMO) removal
- Weighted least-squares