A noise-resistant Wigner-Vile spectrum analysis method based on cyclostationarity and its application in fault diagnosis of rotating

Rolling element bearing and gear are the most common used rotating parts in rotating machinery and they are also the fragile mechanical part. Studying the effective method of timely diagnosis of them is very necessary. The Wigner-Vile spectrum (WVS) is an effective time-frequency analysis and common used method for diagnosis of rotating machinery. However, it would not work effectively when the impulsion characteristic fault signal of rotating machinery is buried by strong background noise. To solve the above problem, the property of cyclostationarity of the rotating machinery signal is used, and the cyclic spectral density basing on second order cyclostationarity statistic is combined with the WVS, and the cyclic spectral density Wigner Vile spectrum (CSDWVS) time-frequency method is proposed in the paper. Through the analysis results of simulation and experiment, the CSDWVS method has the advantages of much more noise-resistant than traditional WVS method, and it could extract the fault feature of the vibration signal of rotating machinery buried in strong background noise. Besides, it also has better time frequency aggregation effect

Seismic characterisation based on time-frequency spectral analysis

We present high-resolution time-frequency spectral analysis schemes to better resolve seismic images for the purpose of seismic and petroleum reservoir characterisation. Seismic characterisation is based on the physical properties of the Earth's subsurface media, and these properties are represented implicitly by seismic attributes. Because seismic traces originally presented in the time domain are non-stationary signals, for which the properties vary with time, we characterise those signals by obtaining seismic attributes which are also varying with time. Among the widely used attributes are spectral attributes calculated through time-frequency decomposition. Time-frequency spectral decomposition methods are employed to capture variations of a signal within the time-frequency domain. These decomposition methods generate a frequency vector at each time sample, referred to as the spectral component. The computed spectral component enables us to explore the additional frequency dimension which exists jointly with the original time dimension enabling localisation and characterisation of patterns within the seismic section. Conventional time-frequency decomposition methods include the continuous wavelet transform and the Wigner-Ville distribution. These methods suffer from challenges that hinder accurate interpretation when used for seismic interpretation. Continuous wavelet transform aims to decompose signals on a basis of elementary signals which have to be localised in time and frequency, but this method suffers from resolution and localisation limitations in the time-frequency spectrum. In addition to smearing, it often emerges from ill-localisation. The Wigner-Ville distribution distributes the energy of the signal over the two variables time and frequency and results in highly localised signal components. Yet, the method suffers from spurious cross-term interference due to its quadratic nature. This interference is misleading when the spectrum is used for interpretation purposes. For the specific application on seismic data the interference obscures geological features and distorts geophysical details. This thesis focuses on developing high fidelity and high-resolution time-frequency spectral decomposition methods as an extension to the existing conventional methods. These methods are then adopted as means to resolve seismic images for petroleum reservoirs. These methods are validated in terms of physics, robustness, and accurate energy localisation, using an extensive set of synthetic and real data sets including both carbonate and clastic reservoir settings. The novel contributions achieved in this thesis include developing time-frequency analysis algorithms for seismic data, allowing improved interpretation and accurate characterisation of petroleum reservoirs. The first algorithm established in this thesis is the Wigner-Ville distribution (WVD) with an additional masking filter. The standard WVD spectrum has high resolution but suffers the cross-term interference caused by multiple components in the signal. To suppress the cross-term interference, I designed a masking filter based on the spectrum of the smoothed-pseudo WVD (SP-WVD). The original SP-WVD incorporates smoothing filters in both time and frequency directions to suppress the cross-term interference, which reduces the resolution of the time-frequency spectrum. In order to overcome this side-effect, I used the SP-WVD spectrum as a reference to design a masking filter, and apply it to the standard WVD spectrum. Therefore, the mask-filtered WVD (MF-WVD) can preserve the high-resolution feature of the standard WVD while suppressing the cross-term interference as effectively as the SP-WVD. The second developed algorithm in this thesis is the synchrosqueezing wavelet transform (SWT) equipped with a directional filter. A transformation algorithm such as the continuous wavelet transform (CWT) might cause smearing in the time-frequency spectrum, i.e. the lack of localisation. The SWT attempts to improve the localisation of the time-frequency spectrum generated by the CWT. The real part of the complex SWT spectrum, after directional filtering, is capable to resolve the stratigraphic boundaries of thin layers within target reservoirs. In terms of seismic characterisation, I tested the high-resolution spectral results on a complex clastic reservoir interbedded with coal seams from the Ordos basin, northern China. I used the spectral results generated using the MF-WVD method to facilitate the interpretation of the sand distribution within the dataset. In another implementation I used the SWT spectral data results and the original seismic data together as the input to a deep convolutional neural network (dCNN), to track the horizons within a 3D volume. Using these application-based procedures, I have effectively extracted the spatial variation and the thickness of thinly layered sandstone in a coal-bearing reservoir. I also test the algorithm on a carbonate reservoir from the Tarim basin, western China. I used the spectrum generated by the synchrosqueezing wavelet transform equipped with directional filtering to characterise faults, karsts, and direct hydrocarbon indicators within the reservoir. Finally, I investigated pore-pressure prediction in carbonate layers. Pore-pressure variation generates subtle changes in the P-wave velocity of carbonate rocks. This suggests that existing empirical relations capable of predicting pore-pressure in clastic rocks are unsuitable for the prediction in carbonate rocks. I implemented the prediction based on the P-wave velocity and the wavelet transform multi-resolution analysis (WT-MRA). The WT-MRA method can unfold information within the frequency domain via decomposing the P-wave velocity. This enables us to extract and amplify hidden information embedded in the signal. Using Biot's theory, WT-MRA decomposition results can be divided into contributions from the pore-fluid and the rock framework. Therefore, I proposed a pore-pressure prediction model which is based on the pore-fluid contribution, calculated through WT-MRA, to the P-wave velocity.Open Acces

Radon spectrogram-based approach for automatic IFs separation

The separation of overlapping components is a well-known and difficult problem in multicomponent signals analysis and it is shared by applications dealing with radar, biosonar, seismic, and audio signals. In order to estimate the instantaneous frequencies of a multicomponent signal, it is necessary to disentangle signal modes in a proper domain. Unfortunately, if signal modes supports overlap both in time and frequency, separation is only possible through a parametric approach whenever the signal class is a priori fixed. In this work, time-frequency analysis and Radon transform are jointly used for the unsupervised separation of modes of a generic frequency modulated signal in noisy environment. The proposed method takes advantage of the ability of the Radon transform of a proper time-frequency distribution in separating overlapping modes. It consists of a blind segmentation of signal components in Radon domain by means of a near-to-optimal threshold operation. The inversion of the Radon transform on each detected region allows us to isolate the instantaneous frequency curves of each single mode in the time-frequency domain. Experimental results performed on constant amplitudes chirp signals confirm the effectiveness of the proposed method, opening the way for its extension to more complex frequency modulated signals

Region-enhanced passive radar imaging

The authors adapt and apply a recently-developed region-enhanced synthetic aperture radar (SAR) image reconstruction technique to the problem of passive radar imaging. One goal in passive radar imaging is to form images of aircraft using signals transmitted by commercial radio and television stations that are reflected from the objects of interest. This involves reconstructing an image from sparse samples of its Fourier transform. Owing to the sparse nature of the aperture, a conventional image formation approach based on direct Fourier transformation results in quite dramatic artefacts in the image, as compared with the case of active SAR imaging. The regionenhanced image formation method considered is based on an explicit mathematical model of the observation process; hence, information about the nature of the aperture is explicitly taken into account in image formation. Furthermore, this framework allows the incorporation of prior information or constraints about the scene being imaged, which makes it possible to compensate for the limitations of the sparse apertures involved in passive radar imaging. As a result, conventional imaging artefacts, such as sidelobes, can be alleviated. Experimental results using data based on electromagnetic simulations demonstrate that this is a promising strategy for passive radar imaging, exhibiting significant suppression of artefacts, preservation of imaged object features, and robustness to measurement noise

Quantitative comparison of analysis methods for spectroscopic optical coherence tomography

Spectroscopic optical coherence tomography (sOCT) enables the mapping of chromophore concentrations and image contrast enhancement in tissue. Acquisition of depth resolved spectra by sOCT requires analysis methods with optimal spectral/spatial resolution and spectral recovery. In this article, we quantitatively compare the available methods, i.e. the short time Fourier transform (STFT), wavelet transforms, the Wigner-Ville distribution and the dual window method through simulations in tissue-like media. We conclude that all methods suffer from the trade-off in spectral/spatial resolution, and that the STFT is the optimal method for the specific application of the localized quantification of hemoglobin concentration and oxygen saturation

An Efficient Algorithm for Instantaneous Frequency Estimation of Nonstationary Multicomponent Signals in Low SNR

A method for components instantaneous frequency (IF) estimation of multicomponent signals in low signal-to-noise ratio (SNR) is proposed. The method combines a new proposed modification of a blind source separation (BSS) algorithm for components separation, with the improved adaptive IF estimation procedure based on the modified sliding pairwise intersection of confidence intervals (ICI) rule. The obtained results are compared to the multicomponent signal ICI-based IF estimation method for various window types and SNRs, showing the estimation accuracy improvement in terms of the mean squared error (MSE) by up to 23%. Furthermore, the highest improvement is achieved for low SNRs values, when many of the existing methods fail.Scopu

Statistical Properties and Applications of Empirical Mode Decomposition

Signal analysis is key to extracting information buried in noise. The decomposition of signal is a data analysis tool for determining the underlying physical components of a processed data set. However, conventional signal decomposition approaches such as wavelet analysis, Wagner-Ville, and various short-time Fourier spectrograms are inadequate to process real world signals. Moreover, most of the given techniques require \emph{a prior} knowledge of the processed signal, to select the proper decomposition basis, which makes them improper for a wide range of practical applications. Empirical Mode Decomposition (EMD) is a non-parametric and adaptive basis driver that is capable of breaking-down non-linear, non-stationary signals into an intrinsic and finite components called Intrinsic Mode Functions (IMF). In addition, EMD approximates a dyadic filter that isolates high frequency components, e.g. noise, in higher index IMFs. Despite of being widely used in different applications, EMD is an ad hoc solution. The adaptive performance of EMD comes at the expense of formulating a theoretical base. Therefore, numerical analysis is usually adopted in literature to interpret the behavior. This dissertation involves investigating statistical properties of EMD and utilizing the outcome to enhance the performance of signal de-noising and spectrum sensing systems. The novel contributions can be broadly summarized in three categories: a statistical analysis of the probability distributions of the IMFs and a suggestion of Generalized Gaussian distribution (GGD) as a best fit distribution; a de-noising scheme based on a null-hypothesis of IMFs utilizing the unique filter behavior of EMD; and a novel noise estimation approach that is used to shift semi-blind spectrum sensing techniques into fully-blind ones based on the first IMF. These contributions are justified statistically and analytically and include comparison with other state of art techniques

Detection and classification of vibrating objects in SAR images

The vibratory response of buildings and machines contains key information that can be exploited to infer their operating conditions and to diagnose failures. Furthermore, since vibration signatures observed from the exterior surfaces of structures are intrinsically linked to the type of machinery operating inside of them, the ability to monitor vibrations remotely can enable the detection and identification of the machinery. This dissertation focuses on developing novel techniques for the detection and M-ary classification of vibrating objects in SAR images. The work performed in this dissertation is conducted around three central claims. First, the non-linear transformation that the micro-Doppler return of a vibrating object suffers through SAR sensing does not destroy its information. Second, the instantaneous frequency (IF) of the SAR signal has sufficient information to characterize vibrating objects. Third, it is possible to develop a detection model that encompasses multiple scenarios including both mono-component and multi-component vibrating objects immersed in noise and clutter. In order to cement these claims, two different detection and classification methodologies are investigated. The first methodology is data-driven and utilizes features extracted with the help of the discrete fractional Fourier transform (DFRFT) to feed machine-learning algorithms (MLAs). Specifically, the DFRFT is applied to the IF of the slow-time SAR data, which is reconstructed using techniques of time-frequency analysis. The second methodology is model-based and employs a probabilistic model of the SAR slow-time signal, the Karhunen-Loève transform (KLT), and a likelihood-based decision function. The performance of the two proposed methodologies is characterized using simulated data as well as real SAR data. The suitability of SAR for sensing vibrations is demonstrated by showing that the separability of different classes of vibrating objects is preserved even after non-linear SAR processing Finally, the proposed algorithms are studied when the range-compressed phase-history data is contaminated with noise and clutter. The results show that the proposed methodologies yields reliable results for signal-to-noise ratios (SNRs) and signal-to-clutter ratios (SCRs) greater than -5 dB. This requirement is relaxed to SNRs and SCRs greater than -10 dB when the range-compressed phase-history data is pre-processed with the Hankel rank reduction (HRR) clutter-suppression technique