Sparse spectral attributes to detect hydrocarbons: A case study from deep dolomite reservoirs

2019 ◽  
Vol 7 (3) ◽  
pp. T701-T711
Author(s):  
Jianhu Gao ◽  
Bingyang Liu ◽  
Shengjun Li ◽  
Hongqiu Wang

Hydrocarbon detection is always one of the most critical sections in geophysical exploration, which plays an important role in subsequent hydrocarbon production. However, due to the low signal-to-noise ratio and weak reflection amplitude of deep seismic data, some conventional methods do not always provide favorable hydrocarbon prediction results. The interesting dolomite reservoirs in Central Sichuan are buried over an average depth of 4500 m, and the dolomite rocks have a low porosity below approximately 4%, which is measured by well-logging data. Furthermore, the dominant system of pores and fractures as well as strong heterogeneity along the lateral and vertical directions lead to some difficulties in describing the reservoir distribution. Spectral decomposition (SD) has become successful in illuminating subsurface features and can also be used to identify potential hydrocarbon reservoirs by detecting low-frequency shadows. However, the current applications for hydrocarbon detection always suffer from low resolution for thin reservoirs, probably due to the influence of the window function and without a prior constraint. To address this issue, we developed sparse inverse SD (SISD) based on the wavelet transform, which involves a sparse constraint of time-frequency spectra. We focus on investigating the applications of sparse spectral attributes derived from SISD to deep marine dolomite hydrocarbon detection from a 3D real seismic data set with an area of approximately [Formula: see text]. We predict and evaluate gas-bearing zones in two target reservoir segments by analyzing and comparing the spectral amplitude responses of relatively high- and low-frequency components. The predicted results indicate that most favorable gas-bearing areas are located near the northeast fault zone in the upper reservoir segment and at the relatively high structural positions in the lower reservoir segment, which are in good agreement with the gas-testing results of three wells in the study area.

Geophysics ◽  
2021 ◽  
Vol 86 (3) ◽  
pp. V245-V254
Author(s):  
Yangkang Chen

Time-frequency analysis is a fundamental approach to many seismic problems. Time-frequency decomposition transforms input seismic data from the time domain to the time-frequency domain, offering a new dimension to probe the hidden information inside the data. Considering the nonstationary nature of seismic data, time-frequency spectra can be obtained by applying a local time-frequency transform (LTFT) method that matches the input data by fitting the Fourier basis with nonstationary Fourier coefficients in the shaping regularization framework. The key part of LTFT is the temporal smoother with a fixed smoothing radius that guarantees the stability of the nonstationary least-squares fitting. We have developed a new LTFT method to handle the nonstationarity in all time, frequency, and space ( x and y) directions of the input seismic data by extending fixed-radius temporal smoothing to nonstationary smoothing with a variable radius in all physical dimensions. The resulting time-frequency transform is referred to as the nonstationary LTFT method, which could significantly increase the resolution and antinoise ability of time-frequency transformation. There are two meanings of nonstationarity, i.e., coping with the nonstationarity in the data by LTFT and dealing with the nonstationarity in the model by nonstationary smoothing. We evaluate the performance of our nonstationary LTFT method in several standard seismic applications via synthetic and field data sets, e.g., arrival picking, quality factor estimation, low-frequency shadow detection, channel detection, and multicomponent data registration, and we benchmark the results with the traditional stationary LTFT method.


2019 ◽  
Vol 38 (4) ◽  
pp. 298-305
Author(s):  
Prashant Kumar Mishra ◽  
Sanjai Kumar Singh ◽  
Pradip Kumar Chaudhuri

The resolution limit of seismic data is an intricate issue that depends not only on frequency and data quality (signal-to-noise ratio) but also on the tools and technology used to analyze seismic response. In this context, the subject of thin-bed delineation is extremely significant for coal-laminated (causing large acoustic impedance contrasts) clastic sequences of the Western Onshore Basin, India. Most of the clastic reservoirs in the area are of subseismic resolution (below 10 m in thickness) due to the low dominant frequency available in seismic data (19–35 Hz). This is where improving seismic resolution is essential for a detailed structural and stratigraphic interpretation. We have implemented a modified workflow with which, by using state-of-the-art techniques of time-frequency decomposition and cepstral analysis, significant seismic bandwidth extension has been achieved. This in turn yields improved vertical resolution of the seismic data with better geologic interpretability. Our approach is named the “syn-cepstral method” after its two integral constituents — synchrosqueezing transform and cepstral analysis. Applying the syn-cepstral method produces better well-to-seismic ties and resolves additional events in comparison to the original seismic data. The validity of syn-cepstral methodology has been demonstrated by 1D and 2D modeling studies followed by application to a 3D seismic data set from the Western Onshore Basin of India. The improvement in thin-bed delineation arising from the increased bandwidth of the resultant data has been validated by well-to-seismic ties and amplitude map interpretation. Thus, while thin clastic reservoir beds in the logs show no discernible presence in the original seismic data, upon application of the syn-cepstral method, the resultant seismic data show improved interpretability of these units.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. WA227-WA240 ◽  
Author(s):  
Guoyin Zhang ◽  
Chengyan Lin ◽  
Yangkang Chen

Microseismic data have a low signal-to-noise ratio (S/N). Existing waveform classification and arrival-picking methods are not effective enough for noisy microseismic data with low S/N. We have adopted a novel antinoise classifier for waveform classification and arrival picking by combining the continuous wavelet transform (CWT) and the convolutional neural network (CNN). The proposed CWT-CNN classifier is applied to synthetic and field microseismic data sets. Results show that CWT-CNN classifier has much better performance than the basic deep feedforward neural network (DNN), especially for microseismic data with low S/N. The CWT-CNN classifier has a shallow network architecture and small learning data set, and it can be trained quickly for different data sets. We have determined why CWT-CNN has better performance for noisy microseismic data. CWT can decompose the microseismic data into time-frequency spectra, where effective signals and interfering noise are easier to distinguish. With the help of CWT, CNN can focus on the specific frequency components to extract useful features and build a more effective classifier.


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.


Geophysics ◽  
2021 ◽  
pp. 1-62
Author(s):  
Wencheng Yang ◽  
Xiao Li ◽  
Yibo Wang ◽  
Yue Zheng ◽  
Peng Guo

As a key monitoring method, the acoustic emission (AE) technique has played a critical role in characterizing the fracturing process of laboratory rock mechanics experiments. However, this method is limited by low signal-to-noise ratio (SNR) because of a large amount of noise in the measurement and environment and inaccurate AE location. Furthermore, it is difficult to distinguish two or more hits because their arrival times are very close when AE signals are mixed with the strong background noise. Thus, we propose a new method for detecting weak AE signals using the mathematical morphology character correlation of the time-frequency spectrum. The character in all hits of an AE event can be extracted from time-frequency spectra based on the theory of mathematical morphology. Through synthetic and real data experiments, we determined that this method accurately identifies weak AE signals. Compared with conventional methods, the proposed approach can detect AE signals with a lower SNR.


Sensors ◽  
2018 ◽  
Vol 18 (11) ◽  
pp. 3985 ◽  
Author(s):  
Siyu Chen ◽  
Yanzhang Wang ◽  
Jun Lin

Residence time difference (RTD) fluxgate sensor is a potential device to measure the DC or low-frequency magnetic field in the time domain. Nevertheless, jitter noise and magnetic noise severely affect the detection result. A novel post-processing algorithm for jitter noise reduction of RTD fluxgate output strategy based on the single-frequency time difference (SFTD) method is proposed in this study to boost the performance of the RTD system. This algorithm extracts the signal that has a fixed frequency and preserves its time-domain information via a time–frequency transformation method. Thereby, the single-frequency signal without jitter noise, which still contains the ambient field information in its time difference, is yielded. Consequently, compared with the traditional comparator RTD method (CRTD), the stability of the RTD estimation (in other words, the signal-to-noise ratio of residence time difference) has been significantly boosted with sensitivity of 4.3 μs/nT. Furthermore, the experimental results reveal that the RTD fluxgate is comparable to harmonic fluxgate sensors, in terms of noise floor.


Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. R989-R1001 ◽  
Author(s):  
Oleg Ovcharenko ◽  
Vladimir Kazei ◽  
Mahesh Kalita ◽  
Daniel Peter ◽  
Tariq Alkhalifah

Low-frequency seismic data are crucial for convergence of full-waveform inversion (FWI) to reliable subsurface properties. However, it is challenging to acquire field data with an appropriate signal-to-noise ratio in the low-frequency part of the spectrum. We have extrapolated low-frequency data from the respective higher frequency components of the seismic wavefield by using deep learning. Through wavenumber analysis, we find that extrapolation per shot gather has broader applicability than per-trace extrapolation. We numerically simulate marine seismic surveys for random subsurface models and train a deep convolutional neural network to derive a mapping between high and low frequencies. The trained network is then tested on sections from the BP and SEAM Phase I benchmark models. Our results indicate that we are able to recover 0.25 Hz data from the 2 to 4.5 Hz frequencies. We also determine that the extrapolated data are accurate enough for FWI application.


Geophysics ◽  
2013 ◽  
Vol 78 (6) ◽  
pp. V229-V237 ◽  
Author(s):  
Hongbo Lin ◽  
Yue Li ◽  
Baojun Yang ◽  
Haitao Ma

Time-frequency peak filtering (TFPF) may efficiently suppress random noise and hence improve the signal-to-noise ratio. However, the errors are not always satisfactory when applying the TFPF to fast-varying seismic signals. We begin with an error analysis for the TFPF by using the spread factor of the phase and cumulants of noise. This analysis shows that the nonlinear signal component and non-Gaussian random noise lead to the deviation of the pseudo-Wigner-Ville distribution (PWVD) peaks from the instantaneous frequency. The deviation introduces the signal distortion and random oscillations in the result of the TFPF. We propose a weighted reassigned smoothed PWVD with less deviation than PWVD. The proposed method adopts a frequency window to smooth away the residual oscillations in the PWVD, and incorporates a weight function in the reassignment which sharpens the time-frequency distribution for reducing the deviation. Because the weight function is determined by the lateral coherence of seismic data, the smoothed PWVD is assigned to the accurate instantaneous frequency for desired signal components by weighted frequency reassignment. As a result, the TFPF based on the weighted reassigned PWVD (TFPF_WR) can be more effective in suppressing random noise and preserving signal as compared with the TFPF using the PWVD. We test the proposed method on synthetic and field seismic data, and compare it with a wavelet-transform method and [Formula: see text] prediction filter. The results show that the proposed method provides better performance over the other methods in signal preserving under low signal-to-noise ratio.


2020 ◽  
Vol 8 (2) ◽  
pp. T217-T229
Author(s):  
Yang Mu ◽  
John Castagna ◽  
Gabriel Gil

Sparse-layer reflectivity inversion decomposes a seismic trace into a limited number of simple layer responses and their corresponding reflection coefficients for top and base reflections. In contrast to sparse-spike inversion, the applied sparsity constraint is less biased against layer thickness and can thus better resolve thin subtuning layers. Application to a 3D seismic data set in Southern Alberta produces inverted impedances that have better temporal resolution and lateral stability and a less blocky appearance than sparse-spike inversion. Bandwidth extension harmonically extrapolated the frequency spectra of the inverted layers and nearly doubled the usable bandwidth. Although the prospective glauconitic sand tunes at approximately 37 m, bandwidth extension reduced the tuning thickness to 22 m. Bandwidth-extended data indicate a higher correlation with synthetic traces than the original seismic data and reveal features below the original tuning thickness. After bandwidth extension, the channel top and base are more evident on inline and crossline profiles. Lateral facies changes interpreted from the inverted acoustic impedance of the bandwidth-extended data are consistent with observations in wells.


Geophysics ◽  
1989 ◽  
Vol 54 (2) ◽  
pp. 181-190 ◽  
Author(s):  
Jakob B. U. Haldorsen ◽  
Paul A. Farmer

Occasionally, seismic data contain transient noise that can range from being a nuisance to becoming intolerable when several seismic vessels try simultaneously to collect data in an area. The traditional approach to solving this problem has been to allocate time slots to the different acquisition crews; the procedure, although effective, is very expensive. In this paper a statistical method called “trimmed mean stack” is evaluated as a tool for reducing the detrimental effects of noise from interfering seismic crews. Synthetic data, as well as field data, are used to illustrate the efficacy of the technique. Although a conventional stack gives a marginally better signal‐to‐noise ratio (S/N) for data without interference noise, typical usage of the trimmed mean stack gives a reduced S/N equivalent to a fold reduction of about 1 or 2 percent. On the other hand, for a data set containing high‐energy transient noise, trimming produces stacked sections without visible high‐amplitude contaminating energy. Equivalent sections produced with conventional processing techniques would be totally unacceptable. The application of a trimming procedure could mean a significant reduction in the costs of data acquisition by allowing several seismic crews to work simultaneously.


Sign in / Sign up

Export Citation Format

Share Document