scholarly journals Seismic Reflection Coefficient Inversion Using Basis Pursuit Denoising in the Joint Time-Frequency Domain

Energies ◽  
2020 ◽  
Vol 13 (19) ◽  
pp. 5025
Author(s):  
Xuegong Zhao ◽  
Hao Wu ◽  
Xinyan Li ◽  
Zhenming Peng ◽  
Yalin Li

Seismic reflection coefficient inversion in the joint time-frequency domain is a method for inverting reflection coefficients using time domain and frequency domain information simultaneously. It can effectively improve the time-frequency resolution of seismic data. However, existing research lacks an analysis of the factors that affect the resolution of inversion results. In this paper, we analyze the influence of parameters, such as the length of the time window, the size of the sliding step, the dominant frequency band, and the regularization factor of the objective function on inversion results. The SPGL1 algorithm for basis pursuit denoising was used to solve our proposed objective function. The applied geological model and experimental field results show that our method can obtain a high-resolution seismic reflection coefficient section, thus providing a potential avenue for high-resolution seismic data processing and seismic inversion, especially for thin reservoir inversion and prediction.

Author(s):  
Shulin Zheng ◽  
Zijun Shen

Complex geological characteristics and deepening of the mining depth are the difficulties of oil and gas exploration at this stage, so high-resolution processing of seismic data is needed to obtain more effective information. Starting from the time-frequency analysis method, we propose a time-frequency domain dynamic deconvolution based on the Synchrosqueezing generalized S transform (SSGST). Combined with spectrum simulation to estimate the wavelet amplitude spectrum, the dynamic convolution model is used to eliminate the influence of dynamic wavelet on seismic records, and the seismic signal with higher time-frequency resolution can be obtained. Through the verification of synthetic signals and actual signals, it is concluded that the time-frequency domain dynamic deconvolution based on the SSGST algorithm has a good effect in improving the resolution and vertical resolution of the thin layer of seismic data.


2021 ◽  
pp. 1-81
Author(s):  
Xiaokai Wang ◽  
Zhizhou Huo ◽  
Dawei Liu ◽  
Weiwei Xu ◽  
Wenchao Chen

Common-reflection-point (CRP) gather is one extensive-used prestack seismic data type. However, CRP suffers more noise than poststack seismic dataset. The events in the CRP gather are always flat, and the effective signals from neighboring traces in the CRP gather have similar forms not only in the time domain but also in the time-frequency domain. Therefore, we firstly use the synchrosqueezing wavelet transform (SSWT) to decompose seismic traces to the time-frequency domain, as the SSWT has better time-frequency resolution and reconstruction properties. Then we propose to use the similarity of neighboring traces to smooth and threshold the SSWT coefficients in the time-frequency domain. Finally, we used the modified SSWT coefficients to reconstruct the denoised traces for the CRP gather. Synthetic and field data examples show that our proposed method can effectively attenuate random noise with a better attenuation performance than the commonly-used principal component analysis, FX filter, and the continuous wavelet transform method.


2017 ◽  
Vol 5 (1) ◽  
pp. T75-T85 ◽  
Author(s):  
Naihao Liu ◽  
Jinghuai Gao ◽  
Zhuosheng Zhang ◽  
Xiudi Jiang ◽  
Qi Lv

The main factors responsible for the nonstationarity of seismic signals are the nonstationarity of the geologic structural sequences and the complex pore structure. Time-frequency analysis can identify various frequency components of seismic data and reveal their time-variant features. Choosing a proper time-frequency decomposition algorithm is the key to analyze these nonstationarity signals and reveal the geologic information contained in the seismic data. According to the Heisenberg uncertainty principle, we cannot obtain the finest time location and the best frequency resolution at the same time, which results in the trade-off between the time resolution and the frequency resolution. For instance, the most commonly used approach is the short-time Fourier transform, in which the predefined window length limits the flexibility to adjust the temporal and spectral resolution at the same time. The continuous wavelet transform (CWT) produces an “adjustable” resolution of time-frequency map using dilation and translation of a basic wavelet. However, the CWT has limitations in dealing with fast varying instantaneous frequencies. The synchrosqueezing transform (SST) can improve the quality and readability of the time-frequency representation. We have developed a high-resolution and effective time-frequency analysis method to characterize geologic bodies contained in the seismic data. We named this method the SST, and the basic wavelet is the three-parameter wavelet (SST-TPW). The TPW is superior in time-frequency resolution than those of the Morlet and Ricker wavelets. Experiments on synthetic and field data determined its validity and effectiveness, which can be used in assisting in oil/gas reservoir identification.


2019 ◽  
Vol 16 (6) ◽  
pp. 1017-1031 ◽  
Author(s):  
Yong Hu ◽  
Liguo Han ◽  
Rushan Wu ◽  
Yongzhong Xu

Abstract Full Waveform Inversion (FWI) is based on the least squares algorithm to minimize the difference between the synthetic and observed data, which is a promising technique for high-resolution velocity inversion. However, the FWI method is characterized by strong model dependence, because the ultra-low-frequency components in the field seismic data are usually not available. In this work, to reduce the model dependence of the FWI method, we introduce a Weighted Local Correlation-phase based FWI method (WLCFWI), which emphasizes the correlation phase between the synthetic and observed data in the time-frequency domain. The local correlation-phase misfit function combines the advantages of phase and normalized correlation function, and has an enormous potential for reducing the model dependence and improving FWI results. Besides, in the correlation-phase misfit function, the amplitude information is treated as a weighting factor, which emphasizes the phase similarity between synthetic and observed data. Numerical examples and the analysis of the misfit function show that the WLCFWI method has a strong ability to reduce model dependence, even if the seismic data are devoid of low-frequency components and contain strong Gaussian noise.


2021 ◽  
Author(s):  
◽  
Jiawen Chua

<p>In most real-time systems, particularly for applications involving system identification, latency is a critical issue. These applications include, but are not limited to, blind source separation (BSS), beamforming, speech dereverberation, acoustic echo cancellation and channel equalization. The system latency consists of an algorithmic delay and an estimation computational time. The latter can be avoided by using a multi-thread system, which runs the estimation process and the processing procedure simultaneously. The former, which consists of a delay of one window length, is usually unavoidable for the frequency-domain approaches. For frequency-domain approaches, a block of data is acquired by using a window, transformed and processed in the frequency domain, and recovered back to the time domain by using an overlap-add technique.  In the frequency domain, the convolutive model, which is usually used to describe the process of a linear time-invariant (LTI) system, can be represented by a series of multiplicative models to facilitate estimation. To implement frequency-domain approaches in real-time applications, the short-time Fourier transform (STFT) is commonly used. The window used in the STFT must be at least twice the room impulse response which is long, so that the multiplicative model is sufficiently accurate. The delay constraint caused by the associated blockwise processing window length makes most the frequency-domain approaches inapplicable for real-time systems.  This thesis aims to design a BSS system that can be used in a real-time scenario with minimal latency. Existing BSS approaches can be integrated into our system to perform source separation with low delay without affecting the separation performance. The second goal is to design a BSS system that can perform source separation in a non-stationary environment.  We first introduce a subspace approach to directly estimate the separation parameters in the low-frequency-resolution time-frequency (LFRTF) domain. In the LFRTF domain, a shorter window is used to reduce the algorithmic delay of the system during the signal acquisition, e.g., the window length is shorter than the room impulse response. The subspace method facilitates the deconvolution of a convolutive mixture to a new instantaneous mixture and simplifies the estimation process.  Second, we propose an alternative approach to address the algorithmic latency problem. The alternative method enables us to obtain the separation parameters in the LFRTF domain based on parameters estimated in the high-frequency-resolution time-frequency (HFRTF) domain, where the window length is longer than the room impulse response, without affecting the separation performance.  The thesis also provides a solution to address the BSS problem in a non-stationary environment. We utilize the ``meta-information" that is obtained from previous BSS operations to facilitate the separation in the future without performing the entire BSS process again. Repeating a BSS process can be computationally expensive. Most conventional BSS algorithms require sufficient signal samples to perform analysis and this prolongs the estimation delay. By utilizing information from the entire spectrum, our method enables us to update the separation parameters with only a single snapshot of observation data. Hence, our method minimizes the estimation period, reduces the redundancy and improves the efficacy of the system.  The final contribution of the thesis is a non-iterative method for impulse response shortening. This method allows us to use a shorter representation to approximate the long impulse response. It further improves the computational efficiency of the algorithm and yet achieves satisfactory performance.</p>


2021 ◽  
Author(s):  
Piotr Krzywiec ◽  
Łukasz Słonka ◽  
Quang Nguyen ◽  
Michał Malinowski ◽  
Mateusz Kufrasa ◽  
...  

&lt;p&gt;In 2016, approximately 850 km of high-resolution multichannel seismic reflection data of the BALTEC survey have been acquired offshore Poland within the transition zone between the East European Craton and the Paleozoic Platform. Data processing, focused on removal of multiples, strongly overprinting geological information at shallower intervals, included SRME, TAU-P domain deconvolution, high resolution parabolic Radon demultiple and SWDM (Shallow Water De-Multiple). Entire dataset was Kirchhoff pre-stack time migrated. Additionally, legacy shallow high-resolution multichannel seismic reflection data acquired in this zone in 1997 was also used. All this data provided new information on various aspects of the Phanerozoic evolution of this area, including Late Cretaceous to Cenozoic tectonics and sedimentation. This phase of geological evolution could be until now hardly resolved by analysis of industry seismic data as, due to limited shallow seismic imaging and very strong overprint of multiples, essentially no information could have been retrieved from this data for first 200-300 m. Western part of the BALTEC dataset is located above the offshore segment of the Mid-Polish Swell (MPS) &amp;#8211; large anticlinorium formed due to inversion of the axial part of the Polish Basin. BALTEC seismic data proved that Late Cretaceous inversion of the Koszalin &amp;#8211; Chojnice fault zone located along the NE border of the MPS was thick-skinned in nature and was associated with substantial syn-inversion sedimentation. Subtle thickness variations and progressive unconformities imaged by BALTEC seismic data within the Upper Cretaceous succession in vicinity of the Kamie&amp;#324;-Adler and the Trzebiat&amp;#243;w fault zones located within the MPS documented complex interplay of Late Cretaceous basin inversion, erosion and re-deposition. Precambrian basement of the Eastern, cratonic part of the study area is overlain by Cambro-Silurian sedimentary cover. It is dissected by a system of steep, mostly reverse faults rooted in most cases in the deep basement. This fault system has been regarded so far as having been formed mostly in Paleozoic times, due to the Caledonian orogeny. As a consequence, Upper Cretaceous succession, locally present in this area, has been vaguely defined as a post-tectonic cover, locally onlapping uplifted Paleozoic blocks. New seismic data, because of its reliable imaging of the shallowest substratum, confirmed that at least some of these deeply-rooted faults were active as a reverse faults in latest Cretaceous &amp;#8211; earliest Paleogene. Consequently, it can be unequivocally proved that large offshore blocks of Silurian and older rocks presently located directly beneath the Cenozoic veneer must have been at least partly covered by the Upper Cretaceous succession; then, they were uplifted during the widespread inversion that affected most of Europe. Ensuing regional erosion might have at least partly provided sediments that formed Upper Cretaceous progradational wedges recently imaged within the onshore Baltic Basin by high-end PolandSPAN regional seismic data. New seismic data imaged also Paleogene and younger post-inversion cover. All these results prove that Late Cretaceous tectonics substantially affected large areas located much farther towards the East than previously assumed.&lt;/p&gt;&lt;p&gt;This study was funded by the Polish National Science Centre (NCN) grant no UMO-2017/27/B/ST10/02316.&lt;/p&gt;


Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. V235-V247 ◽  
Author(s):  
Duan Li ◽  
John Castagna ◽  
Gennady Goloshubin

The frequency-dependent width of the Gaussian window function used in the S-transform may not be ideal for all applications. In particular, in seismic reflection prospecting, the temporal resolution of the resulting S-transform time-frequency spectrum at low frequencies may not be sufficient for certain seismic interpretation purposes. A simple parameterization of the generalized S-transform overcomes the drawback of poor temporal resolution at low frequencies inherent in the S-transform, at the necessary expense of reduced frequency resolution. This is accomplished by replacing the frequency variable in the Gaussian window with a linear function containing two coefficients that control resolution variation with frequency. The linear coefficients can be directly calculated by selecting desired temporal resolution at two frequencies. The resulting transform conserves energy and is readily invertible by an inverse Fourier transform. This modification of the S-transform, when applied to synthetic and real seismic data, exhibits improved temporal resolution relative to the S-transform and improved resolution control as compared with other generalized S-transform window functions.


2007 ◽  
Vol 13 ◽  
pp. 17-20 ◽  
Author(s):  
Erik S. Rasmussen ◽  
Thomas Vangkilde-Pedersen ◽  
Peter Scharling

Intense investigations of deep aquifers in Jylland, western Denmark, during the last seven years have resulted in de tailed mapping of Miocene sand-rich deposits laid down in fluvial channels, delta lobes, shoreface and spit complexes (Fig. 1; Rasmussen 2004). Detailed sedimentological and paly nol ogical studies of outcrops and cores, and interpretation of high-resolution seismic data, have resulted in a well-founded sequence-stratigraphic and lithostratigraphic scheme (Fig. 1) suitable for prediction of the distribution of sand. The Miocene succession onshore Denmark is divided into three sand-rich deltaic units: the Ribe and Bastrup sands and the Odderup Formation (Fig. 2). Prodeltaic clayey deposits of the Vejle Fjord and Arnum Formations interfinger with the sand-rich deposits. Most of the middle and upper Mio- cene in Denmark is composed of clayey sediments referred to the Hodde and Gram Formations (Fig. 2). This paper presents examples of seismic reflection patterns that have proved to correlate with sand-rich deposits from lower Miocene deltaic deposits and that could be applied in future exploration for aquifers and as analogues for oil- and gas-bearing sands in wave-dominated deltas.


Author(s):  
Jordi Burriel-Valencia ◽  
Ruben Puche-Panadero ◽  
Javier Martinez-Roman ◽  
Angel Sapena-Bano ◽  
Martin Riera-Guasp ◽  
...  

Induction machines drive many industrial processes, and their unexpected failure can cause heavy production losses. The analysis of the current spectrum can identify online the characteristic fault signatures at an early stage, avoiding unexpected breakdowns. Nevertheless, frequency domain analysis requires stable working conditions, which is not the case for wind generators, motors driving varying loads, etc. In these cases an analysis in the time-frequency domain -such as a spectrogram- is required for detecting faults signatures. The spectrogram is built using the short frequency Fourier transform, but its resolution depends critically on the time window used to generate it: short windows provide good time resolution, but poor frequency resolution, just the opposite than long windows. Therefore, the window must be adapted at each time to the shape of the expected fault harmonics, by highly skilled maintenance personnel. In this paper, this problem is solved with the design of a new multi-band window, which generates simultaneously many different narrow-band current spectrograms, and combines them into a single, high resolution one, without the need of manual adjustments. The proposed method is validated with the diagnosis of bar breakages during the start-up of a commercial induction motor.


Sign in / Sign up

Export Citation Format

Share Document