Matching Demodulation Synchrosqueezing S Transform and its Application in Seismic Time–Frequency Analysis

Synchrosqueezing transform (SST) benefits from an instantaneous frequency (IF) estimator in the time-frequency domain, providing an energy-concentrated time-frequency representation to describe the time-varying frequency of seismic signals. To enhance the concentration performance of SST, this paper theoretically proposes a statistical SST (SSST) by constructing a spectrum-weighted IF estimator in the short-time Fourier transform (STFT) domain. In this SSST, a linear chirp signal is introduced to better capture its chirp rate, and then the quadratic STFT spectra are weighted by a window-related weighting function. On this basis, two equations related to the IF and chirp rate are constructed to derive the spectrum-weighted IF estimator. Finally, the STFT coefficients are squeezed to the estimated IF trajectories by a frequency fixed-point iterative algorithm, thereby providing a more concentrated time-frequency representation for multi-component signals than existing advanced methods, while enabling to retrieve each component. Two synthetic examples and one field seismic data on thin interbeds are utilized to demonstrate the effectiveness of the proposed SSST and show its ability to highlight the time-varying frequency features of seismic signals, which is a promising seismic data analysis tool, such as characterizing thin interbed thickness variations.

Synchrosqueezing transform (SST) is an effective time-frequency analysis (TFA) approach for the processing of nonstationary signals. The SST shows a satisfactory ability of the TF localization of the nonlinear signal with a slowly time-varying instantaneous frequency (IF). However, for the signal of which ridge curves in the TF domain are fast varying, or even almost parallel to the frequency axis, the SST will provide a blurred TF representation (TFR). To solve this issue, the transient-extracting transform (TET) was recently put forward. The TET can effectively characterize and extract transient features in the much concentrated TFR for the strongly frequency-modulated (FM) signal, especially the impulse-like signal. However, contrary to the SST, it is not suitable for weak FM modes. In this study, we propose a TFA method called the time-synchroextracting general chirplet transform (TEGCT). The TEGCT can achieve a highly concentrated TFR for strong FM signals as well as weak FM ones. Quantized indicators, the concentration measurement and the peak signal-to-noise ratio, are used to analyze the performances of the proposed method compared with those of other methods. The comparisons show that the TEGCT can provide a result with better TF localization. Then, the proposed method was applied to the spectrum analysis of the seismic data for oil reservoir characteristics. The horizontal slices of the offshore 3-D seismic data show that the TEGCT delineates more distinct and continuous subsurface channels in a fluvial-delta deposition system. All the results illustrate that our proposed method is a good potential tool for seismic processing and interpretation in the geoscience.

The use of linear time-frequency (TF) representa-tions as instantaneous frequency (IF) estimators arises in various fundamental disciplines, mainly thanks to their simplicity and immunity against interfering cross-terms. In this paper, we derive the optimal window width for IF estimation of noisy multicomponent signals based on a family of linear TF transforms that use Gaussian windowing functions and the Fourier oscillatory kernel. Closed-form formulas concerning the estimation bias and the variance are presented, which, thanks to their generality, describe the statistical performance of IF estimators based on transforms with fixed, time-adaptive, frequency-adaptive, or time-frequency adaptive windows. The optimal window is dependent on the unknown first derivative of the IF; therefore, we employ a low-complexity procedure to infer the derivative and optimize the width accordingly. Two adaptive and fully automated TF representations (TFRs) are developed; the first employs a time-adaptive window that minimizes the sum of the mean-squared errors (MSEs) of the IF estimates at each time instant, while in the second TFR, the window is adaptive over time and frequency, minimizing the estimation MSE at each location in the TF domain. Examples using synthetic and real-world signals demonstrate that the proposed algorithms may outperform many popular state-of-the-art techniques, including those that are signal-adaptive, in terms of IF estimation.

Synchronous Spline-kernelled Chirplet Squeezing Transform and Its Application for Seismic Data Analysis

Synchrosqueezing transform (SST) is a high resolution time frequency representation technology for nonstationary signal analysis. The short time Fourier transform-based synchrosqueezing transform (FSST) and the S transform-based synchrosqueezing transform (SSST) time frequency methods are effective tools for bearing fault signal analysis. The fault signals belong to a non-Gaussian and nonstationary alpha (α) stable distribution with 1<α<2 and even the noises being also α stable distribution. The conventional FSST and SSST methods degenerate and even fail under α stable distribution noisy environment. Motivated by the fact that fractional low order STFT and fractional low order S-transform work better than the traditional STFT and S-transform methods under α stable distribution noise environment, we propose in this paper the fractional lower order FSST (FLOFSST) and the fractional lower order SSST (FLOSSST). In addition, we derive the corresponding inverse FLOSST and inverse FLOSSST. The simulation results show that both FLOFSST and FLOSSST perform better than the conventional FSSST and SSST under α stable distribution noise in instantaneous frequency estimation and signal reconstruction. Finally, FLOFSST and FLOSSST are applied to analyze the time frequency distribution of the outer race fault signal. Our results show that FLOFSST and FLOSSST extract the fault features well under symmetric stable (SαS) distribution noise.

A robust parameter estimation of FHSS signals using time–frequency analysis in a non-cooperative environment

Chapter 10 - Instantaneous Frequency Estimation and Localization

This study proposes a high-resolution processing technique for seismic data based on the improved synchrosqueezing transform (SST). The SST rearranges the complex spectrum of the wavelet transform along the frequency axis. However, the energy is not concentrated in the position with the faster frequency change rate. To overcome these problems, the proposed method first transformed the seismic signal into an analytical one via Hilbert transformation and then determined the phase correction value (x0(t)) before and after the transformation. This was achieved using the instantaneous frequency attribute (ω(t)) of the analytical signal and a specific frequency (ω0) selected as the dominant frequency of the input data to adjust the change rate of the signal frequency before the SST falls close to zero, which can improve the time–frequency resolution after compression transformation. The essence of this method only adjusts the phase of the signal before and after transformation, without affecting the inverse expression of the SST. The one- and two-dimensional model data results show that the proposed method can better identify the time–frequency distribution characteristics of seismic signals and improve their resolution and the focusing effect in the time–frequency domain. The proposed method shows good prospects for application in gas reservoir detection and identification.

Classical time-frequency analysis is based on the amplitude responses of bandpass filters, discarding phase information. Instantaneous frequency analysis, in contrast, is based on the derivatives of these phases. This method of frequency calculation is of interest for its high precision and for reasons of similarity to cochlear encoding of sound. This article describes a methodology for high resolution analysis of sparse sounds, based on instantaneous frequencies. In this method, a comparison between tonotopic and instantaneous frequency information is introduced to select filter positions that are well matched to the signal. Second, a cross-check that compares frequency estimates from neighboring channels is used to optimize filter bandwidth, and to signal the quality of the analysis. These cross-checks lead to an optimal time-frequency representation without requiring any prior information about the signal. When applied to a signal that is sufficiently sparse, the method decomposes the signal into separate time-frequency contours that are tracked with high precision. Alternatively, if the signal is spectrally too dense, neighboring channels generate inconsistent estimates-a feature that allows the method to assess its own validity in particular contexts. Similar optimization principles may be present in cochlear encoding.

Slope synchronous chirplet transform and its application to tacho-less order tracking of rotating machineries

The Wigner distribution (WD) produces highly concentrated time-frequency (TF) representation of nonstationary signals. It may be used as an efficient signal analysis tool, including the cases of frequency modulated signals corrupted with the Gaussian noise. In some applications, a significant amount of impulse noise is present. Then, the WD fails to produce satisfactory results. The robust periodogram has been introduced for spectral estimation of this kind of noisy signals. It can produce good concentration for pure harmonic signals. However, it is not so efficient in the cases of signals with rapidly varying frequency. This is the motivation for introducing the robust WD. It is a reliable TF representation tool for wide class of nonstationary signals corrupted with impulse noise. This distribution produces good accuracy of the instantaneous frequency (IF) estimation. Using the Huber (1981) loss function, a generalization of the WD is presented. It includes both the standard and the robust WD as special cases. This distribution can be used for TF analysis of signals corrupted with a mixture of impulse and Gaussian noise. The presented theory is illustrated on examples, including applications on the IF estimation and time-varying filtering of signals corrupted with a mixture of the Gaussian and impulse noise. The case study analysis of the IF estimators' accuracy, based on the standard and the robust WD forms, is performed. In order to improve the IF estimation, a median filter is applied on the obtained IF estimate.

XWD-algorithm for the instantaneous frequency estimation revisited: Statistical analysis

This paper explores three groups of time–frequency distributions: the Cohen’s, affine, and reassigned classes of time–frequency representations (TFRs). This study provides detailed insight into the theory behind the selected TFRs belonging to these classes. Extensive numerical simulations were performed with examples that illustrate the behavior of the analyzed TFR classes in the joint time–frequency domain. The methods were applied both on synthetic and real-life non-stationary signals. The obtained results were assessed with respect to time–frequency concentration (measured by the Rényi entropy), instantaneous frequency (IF) estimation accuracy, cross-term presence in the TFRs, and the computational cost of the TFRs. This study gives valuable insight into the advantages and limitations of the analyzed TFRs and assists in selecting the proper distribution when analyzing given non-stationary signals in the time–frequency domain.

Time-frequency(T-F) analysis has proved to be a powerful tools in seismic data processing and interpretation. One of the most important application of T-F analysis in seismic data processing is the spectral decomposition technique. Using the spectral magnitude, we can identify wave filed attenuation due to absorption and lateral changes in stratigraphy. The principal approaches include the short time Fourier transform and wavelet transform, which have limitations such as trade-offs between time and frequency resolution. Synchrosqueezing wave packet transform(SSWPT) has recently been proposed as a alternative method of the general empirical mode decomposition and has a consolidated mathematical background. It is an extension of the wavelet transform incorporating frequency reassignment techniques, this new tools produces high-resolution time-frequency representation. We introduce the SSWPT to seismic data analysis and validate its effectiveness with nonstationary synthetic example. Field examples illustrate that the SSWPT can be potentially used to detect frequency shadows caused by hydrocarbons and to identify subtle stratigraphic features for reservoir characterization.

Self-paced learning for instantaneous frequencies estimation in heavy noise environments