Extension of the sampling theorem
![]() | |
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 |
Contents
Problem
Explain why sampling at 4 ms is sufficient to reproduce exactly a signal whose spectrum is within the range to , but not if the bandwidth is shifted upward?
Background
To convert an analog (continuous) signal to the digital form, values of the signal are measured at a fixed interval , called the sampling interval. The result is a series of numbers representing the amplitudes and polarities of the signal at the times , where is an integer (often we omit and specify the time by giving only). We write for the digital function corresponding to .
Sampling is equivalent to multiplying the analog signal by a comb, a function consisting of an infinite series of unit impulses spaced at a fixed interval . The equation of comb(t) is
( )
The transform of comb(t) is another comb with impulses in the frequency domain at intervals :
( )
where [see Sheriff and Geldart, 1995, equation (15.155)].
The Nyquist frequency, , is half the frequency of sampling , that is,
( )
While a signal that does not contain frequencies higher than will be recorded accurately, a frequency higher than by the amount (i.e., a frequency that is not sampled at least twice per cycle) will produce an alias frequency equal to (Sheriff and Geldart, 1995, section 9.2.2c).
The convolution theorem, equation (9.3f), has a converse [see Sheriff and Geldart, 1995, equation (15.146)]:
( )
Solution
We examine first the case where the signal frequency spectrum extends from to 125 Hz, next consider modifications when the spectrum is shifted upward 125 Hz for 4-ms sampling), and finally, the case where the signal is shifted upward by an amount .
Part (i) of Figure 9.4a shows the signal and its frequency spectrum, the latter being confined to the range Hz [the negative frequencies arise when we use Euler’s equations (Sheriff and Geldart, 1995, problem 15.12a) to transform equation (9.1a) into (9.1c). Part (ii) shows comb(t), with elements 4 ms apart, and its transform, a comb with elements spaced Hz apart. Digitizing the signal is equivalent to multiplying it by comb(t) [part (iii)], which is equivalent to convolving the spectrum with comb(f) [see equation (9.4b)]; this results in a repetition of the spectrum for each impulse in comb(f) [see equation (9.3j)]. To eliminate the repeated spectra, we multiply by a boxcar [part (iv)], thus getting back the original spectrum [part (v)]. Multiplication in the frequency domain is equivalent to convolution in the time domain, so we convolve with a sinc function in part (V) and recover the original time-domain signal.
If the signal is shifted upward by , where is an integer, then the sampling interval is larger than 1/2 cycle for all frequencies and the resulting spectrum is that of the alias function rather than of the signal, as illustrated in Figure 9.4b. Because we no longer have the signal spectrum after sampling, we cannot recover the signal.
If the signal is shifted upward so that it overlaps the Nyquist frequency (Figure 9.4c), the sampling causes the frequencies above to overlap and distort the signal spectrum. Now we cannot separate the signal spectrum from the distorted spectrum and hence cannot recover the signal.
Continue reading
Previous section | Next section |
---|---|
Fourier transforms of the unit impulse and boxcar | Alias filters |
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