297 research outputs found
Recommended from our members
Numerical treatment of seismic accelerograms and of inelastic seismic structural responses using harmonic wavelets
The harmonic wavelet transform is employed to analyze various kinds of nonstationary signals common in aseismic design. The effectiveness of the harmonic wavelets for capturing the temporal evolution of the frequency content of strong ground motions is demonstrated. In this regard, a detailed study of important earthquake accelerograms is undertaken and smooth joint time-frequency spectra are provided for two near-field and two far-field records; inherent in this analysis is the concept of the mean instantaneous frequency. Furthermore, as a paradigm of usefulness for aseismic structural purposes, a similar analysis is conducted for the response of a 20-story steel frame benchmark building considering one of the four accelerograms scaled by appropriate factors as the excitation to simulate undamaged and severely damaged conditions for the structure. The resulting joint time-frequency representation of the response time histories captures the influence of nonlinearity on the variation of the effective natural frequencies of a structural system during the evolution of a seismic event. In this context, the potential of the harmonic wavelet transform as a detection tool for global structural damage is explored in conjunction with the concept of monitoring the mean instantaneous frequency of records of critical structural responses
Recommended from our members
Time-frequency representation of earthquake accelerograms and inelastic structural response records using the adaptive chirplet decomposition and empirical mode decomposition
In this paper, the adaptive chirplet decomposition combined with the Wigner-Ville transform and the empirical mode decomposition combined with the Hilbert transform are employed to process various non-stationary signals (strong ground motions and structural responses). The efficacy of these two adaptive techniques for capturing the temporal evolution of the frequency content of specific seismic signals is assessed. In this respect, two near-field and two far-field seismic accelerograms are analyzed. Further, a similar analysis is performed for records pertaining to the response of a 20-story steel frame benchmark building excited by one of the four accelerograms scaled by appropriate factors to simulate undamaged and severely damaged conditions for the structure. It is shown that the derived joint time–frequency representations of the response time histories capture quite effectively the influence of non-linearity on the variation of the effective natural frequencies of a structural system during the evolution of a seismic event; in this context, tracing the mean instantaneous frequency of records of critical structural responses is adopted.
The study suggests, overall, that the aforementioned techniques are quite viable tools for detecting and monitoring damage to constructed facilities exposed to seismic excitations
The estimation of geoacoustic properties from broadband acoustic data, focusing on instantaneous frequency techniques
The compressional wave velocity and attenuation of marine sediments are fundamental to marine science. In order to obtain reliable estimates of these parameters it is necessary to examine in situ acoustic data, which is generally broadband. A variety of techniques for estimating the compressional wave velocity and attenuation from broadband acoustic data are reviewed. The application of Instantaneous Frequency (IF) techniques to data collected from a normal-incidence chirp profiler is examined. For the datasets examined the best estimates of IF are obtained by dividing the chirp profile into a series of sections, estimating the IF of each trace in the section using the first moments of the Wigner Ville distribution, and stacking the resulting IF to obtain a composite IF for the section. As the datasets examined cover both gassy and saturated sediments, this is likely to be the optimum technique for chirp datasets collected from all sediment environments
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
Instantaneous Spectral Analysis: Time-Frequency Mapping via Wavelet Matching with Application to Contaminated-Site Characterization by 3D GPR
Spectral decomposition, by which a time series is transformed from the 1D time/amplitude domain to the 2D time/spectrum domain, has become a popular and useful tool in seismic exploration for hydrocarbons. The windowed, or short-time Fourier transform (STFT) was one early approach to computing the time-frequency (t-f) distribution. This method relies on the user selecting a fixed time window, then computing the Fourier spectrum within the time window while sliding the window along the length of the trace. The primary limitation of the STFT is the fixed window which prevents either time localization of high frequency components (if a long window is used) or spectral resolution of the low-frequency components (if a short window is used)
- …