Seismic data processing is an important aspect to improve the signal to noise ratio. The main work of this paper is to combine the characteristics of seismic data, using wavelet transform method, to eliminate and control such random noise, aiming to improve the signal to noise ratio and the technical methods used in large data systems, so that there can be better promotion and application. In recent years, prestack data denoising of alldigital threedimensional seismic data is the key to data processing. Contrapose the characteristics of alldigital threedimensional seismic data, and, on the basis of previous studies, a new threshold function is proposed. Comparing between conventional hard threshold and soft threshold, this function not only is easy to compute, but also has excellent mathematical properties and a clear physical meaning. The simulation results proved that this method can well remove the random noise. Using this threshold function in actual seismic processing of unconventional lithologic gas reservoir with low porosity, low permeability, low abundance, and strong heterogeneity, the results show that the denoising method can availably improve seismic processing effects and enhance the signal to noise ratio (SNR).
The study area is located in the middle of Sulige gas field, with small effective reservoir thickness, expressed as “thin, more scattered, and miscellaneous” feature on the vertical of reservoir; therefore, to remove noise in seismic data, which improves the reliability of the signal, is particularly important for reservoir prediction in later work, as it is always the goal of geophysical workers.
Random noise cannot seek a unified time pattern, with wide bandwidth, and cannot be determined by its apparent velocity and direction of propagation, thus removing random noise and the key point is to improve its SNR. In order to effectively remove random noise, various denoising methods have been proposed; thereinto, the wavelet transform method works better in recent years. Wavelet transform theory developed in mid and late eighties of last century [
Since desired signal and random noise have different characteristics on different wavelet scales, such as during valid signal occurs within a short time, instantaneous modulus maxima of the wavelet coefficients occurs, and with the decomposition scale increment, wavelet coefficients and dense maxima of random noise in desired signal will decrease while decomposition scale increasing [
But hard and soft threshold methods commonly used in wavelet threshold methods have some shortcomings, such as the following: remodeling factor in hard threshold processing occurs intermittently at the threshold point; the reconstructed signal is prone to PseudoGibbo phenomenon; while although reconstruction coefficient of soft threshold has good overall continuity, there is always a constant bias between reconstructed coefficients and decomposition coefficients [
Consider the following observation signal:
Wavelet threshold method based on the minimum mean square error (MSE) proposed by Donoho has been widely used. The basic flow is shown in Figure
The flow of wavelet threshold.
Hard threshold function is
Soft threshold function is
The value
A weighted average threshold estimation method is given in literatures [
Traditional weighted average methods generally set the weighting factor
The new threshold function has the following characteristics:
When
When
This shows that the new threshold function is an adjustable function between hard and soft threshold function (Figure
The new threshold.
To verify effect of the new threshold denoising function, we use Matlab toolbox containing Gaussian noise Heavisine signals. Using traditional hard and soft threshold function and the new threshold function make denoising tests. Then comparative analysis of results, select sym6 to make wavelet decomposition, the maximum resolution scale
Methods of signal to noise ratio (SNR/dB) and mean square error (MSE).
Threshold  Plus noise signal  Hard threshold  Soft threshold function  The new threshold 


SNR  14.8606  24.7590  25.2749  26.5035 
RMSE  1.0324  0.3303  0.3113  0.2702 
The effect diagram of different threshold function.
By the numerical denoising effect evaluation in 3–1 one can see several improved threshold function denoising effects showing soft, hard threshold method is better than the classical one note. We note that the denoising effects of exponential threshold function are better than other kinds of methods; and the method of the new threshold function in root mean square error slightly inferior to the exponential threshold function, but in the signal to noise ratio of this index, is obviously better than the exponential threshold method [
Typical analysis from the original records, shot records SNR differences, higher SNR from 1.8 s to 2.1 s, the whole region can continuously track mark phase axis, with good quality single shot, shallow, midwave group is complete, 2.5 s or less energy loss is more serious. Poor quality single shot shallow, deep SNR poor record serrated early to beat poor continuity (Figures
Better quality of the original single shot.
Poor quality of the original single shot.
Since the surface area of seismic and geological conditions is complex, the entire work area accounts for about 7% of Salt Lake wetlands; desert area is about 50%; grass and sand grass are about 35%; rock outcrop is about 8%. Overall, qualities of collected data are quite different; SNR and continuous quality display banded regional variation, and desert zone data was poor, which is shown in Figure
Stacked section INLINE551 CDP691–961.
The original single shot shows that work area noise is relatively developed (Figure
Mechanical interference: performance on the record is outliers and abnormal trace, and characteristics of such interference is with strong randomness distribution and different stimulation time; different positions have different distributions, easily recognizable from the amplitude [
Linear ramp interference: it includes multiple refraction, surface wave, acoustic and some secondary oblique interference, multiple refraction, and secondary oblique interference which are with larger apparent velocity on the record, generally more than 1000 m/s, with lowfrequency surface waves, mainly distributed in the gun point neighborhood, single shot on the record as low apparent velocity, generally not more than 800 m/s, with high sound waves frequency and strong energy.
Random interference.
Refracted waves, multiple refraction, direct waves, surface waves, and sound waves.
Combined with the above analysis of raw data, we can see that interference of threedimensional lowfrequency surface wave is strong in the study area, effective reflections on the original record are full inundated by surface wave, with low SNR; in addition, shallow direct wave interference is strong, with broad band distribution, which seriously interferes with shallow reflection signals [
To illustrate the effectiveness of the proposed denoising method, we selected a small crosssectional view of seismic data of a work area to actual test. Where Figure
Actual seismic records.
The premise of denoising is without prejudice to the effective reflection information, for low SNR data; if it cannot well suppress noise, it is difficult to have effectively phase axis. In order to achieve better denoising effect, after several trials comparing, one finds that using this method is superior to conventional denoising methods [
From the results, the data for the effective band are from 10 to 70 Hz; frequency ranges are from 13 to 45 Hz, which reached the processing demands and meet the requirements of seismic data interpretation and reservoir inversion. To the whole cross section, the processing result is ideally effective, and SNR of profile and resolution are moderate, with natural wave, and the main target layer wave groups characteristics are distinct, clear, stable, reliable, and with good profile backgrounds, which achieve good results. This confirms the effectiveness of using data denoising algorithm for the complex seismic and geological conditions blocks. (Figures
Stacked section INLINE1081.
Poststack migration profile (INLINE791 CDP194–476).
To further verify compliance between denoising method and actual seismic exploration, we selected well tie cross section of two representative wells in the area to calibrate and analyze, and we found that synthesis record and main boreholeside well layer are highly consistent. Results show that processing outcome and actual drilling revealing geological features is very accordant (Figures
X Wells calibration results.
Y Wells calibration results.
Based on wavelet threshold methods proposed by D. L. Dohono and I. M. Johnstone, we present a new threshold function, which enables the superiority of threshold wavelet transform method fully realized. In order to test the effectiveness of the proposed method, we first made ideal synthetic seismograms and achieved the desired results. Finally, we used the method for denoising actual seismic data, and the results show that this method can effectively remove random noise in original seismic data and improve SNR of seismic data, which provides high quality data for the later work of stratigraphic interpretation and has some practical significance.
The authors declare that there is no conflict of interests regarding the publication of this paper.
Foundation item: The Joint Specialized Research Fund for the Doctoral Program of Higher Education (20135122110010) and National Natural Science Foundation of China (4127236).