Time-frequency analysis of seismic data using local attributes |

The Fourier transform has found various applications in signal analysis. The classic Fourier transform indicates the presence of different frequencies within the analysis window, but does not show where in that window the particular frequency components reside. Localized frequency information can be obtained by computing the Fourier transform with a temporally shifted window. Such an approach to time-frequency analysis is known as the STFT (Allen, 1977). The window function is commonly parameterized by window size, overlap, and taper. Once the window function has been chosen for the STFT, temporal and spectral resolutions are fixed for the entire time-frequency map.

The S-transform (Stockwell et al., 1996) is similar to the STFT,
but with a Gaussian-shaped window whose width scales inversely
with frequency. The expression of the S-transform is

Consider a signal on . The Fourier series, assuming a
periodic extension of the boundary conditions, can be expressed as

and represent the series coefficients

where , equation 1 can be written as

If frequency is finite, the range of becomes , . In linear notation, can be obtained by solving the least-squares problem

where denotes the squared norm of a function. By allowing coefficients to change with time , we can define the timevarying coefficients via the following least-squares problem

The Fourier coefficient in equation 7 is a function of time and frequency and . The numerical support of frequency f can be between zero and the Nyquist frequency (Cohen, 1995), and the interval of frequency can be . In practical applications, the range of frequencies can also be assigned by the user. In matrix notation, equation 7 can be written as

where denotes a diagonal matrix which is composed from the elements of .

The problem of minimization in equation 7 is mathematically
ill-posed because it is severely underconstrained: There are more
unknown variables than constraints. To solve this ill-posed problem,
we force the coefficients
to have a desired behavior, such as
smoothness. With the addition of a regularization term, equation 7
becomes

where is the Tikhonov regularization operator (roughening operator) and is a scalar regularization parameter.

Shaping regularization (Fomel, 2007b) provides a particularly convenient method of enforcing smoothness in iterative optimization schemes. In the appendix, we review the method of shaping regularization in detail. Fomel (2009) used shaping regularization to constrain the problem of nonstationary regression. In this paper, we use shaping regularization, analogous to the problem of nonstationary regression, to constrain estimated coefficients. We choose our shaping operator to be Gaussian smoothing with an adjustable radius. In shaping regularization, the radius of the Gaussian smoothing operator controls the smoothness of coefficients .

Once we obtain time-varying Fourier coefficients
, the time-frequency map is defined as

Consider a simple signal (Figure 1a), which includes three monophonic components at 10, 20, and 30 Hz and two broad-band spikes at 2 and 2.3 s. Time-frequency maps generated by the S-transform and our method are shown, respectively, in Figure 1b and 1c. Frequency components appear to be represented well from low to high frequencies in the proposed method. In comparison with the S-transform, the proposed method provides superior resolution for monophonic waves. At the lower frequencies, the S-transform is as good as the proposed method on the spectral resolution for monophonic waves. Note that the S-transform represents spikes better at high frequencies. However, because the width of analysis window is large at low frequencies, the S-transform provides poor time resolution at low frequencies for spikes. When we use the shaping operator as a regularization term, there is an edge effect from smoothing, which appears near 0 s and 4 s (Figure 1c). Figure 1d displays the time-frequency map using a different parameter (30-point smoothing radius) to demonstrate adjustable time-frequency representation of the proposed method.

s-0,st-0,proj-0,proj-0-30
(a) Synthetic signal with three constant frequency
components and two spikes. (b) Time-frequency map of the S-transform.
(c) Time-frequency map of the proposed method with smoothing radius of 15 points.
(d) Time-frequency map of the proposed method with smoothing radius of 30 points.
Figure 1. |
---|

To show the effect of varying frequencies, we provide a composite chirp signal (Figure 2a), which includes two parabolic frequencies, each having a constant amplitude. Figure 2b and 2c shows the results of the S-transform and the proposed method with 10-point smoothing radius, respectively. The two frequency components are detected with higher time and frequency resolution by our method. The S-transform has high spectral resolution near low frequencies, but loses resolution at high frequencies. Resolution of the time-frequency map deteriorates in both methods when the curvature of the time-frequency curve becomes large (as indicated by the arrow).

s-1,st-1,proj-1
(a) Synthetic chirp signal with two parabolic frequency changes.
(b) Time-frequency map of the S-transform. (c) Time-frequency
map of the proposed method.
Figure 2. |
---|

In the proposed method, the smoothing radius used in shaping regularization is a parameter which controls the smoothness of the model. It is different from the method of the STFT and the S-transform, in which division into local windows is applied to the data to localize frequency content in time. In contrast to the wavelet-based methods (Chakraborty and Okaya, 1995; Wang, 2007; Sinha et al., 2005), our method, as a straightforward extension of the classic Fourier analysis, is different because of the choice of basis functions. The wavelet-based methods need to decide the wavelets, which can be regarded as patterns of seismic data. The closer the selected wavelets are to the true pattern of the input data, the better the result of the time-frequency analysis. In this paper, we use more general Fourier basis functions to compute the time-frequency map.

Time-frequency analysis of seismic data using local attributes |

2013-07-26