The practice of exploration and production has proved that explosives are excited in different surrounding rocks and the seismic wavelets collected have different characteristics. In this paper, by establishing a nume...The practice of exploration and production has proved that explosives are excited in different surrounding rocks and the seismic wavelets collected have different characteristics. In this paper, by establishing a numerical model of the explosion in the well, using finite element analysis technology for numerical simulation, the simulation calculated the stress structure in the near-source area of the earthquake excitation, and extracted the seismic wavelet. The results show that the simulation seismic wavelet characteristics of different thin interbedded sand and mudstone structures have changed significantly. Through excitation simulation, the amplitude and spectrum information of seismic wavelets can be compared and analyzed, and the excitation parameters can be optimized. .展开更多
Phase spectrum estimation of the seismic wavelet is an important issue in high-resolution seismic data processing and interpretation. On the basis of two patterns of constant-phase rotation and root transform for wave...Phase spectrum estimation of the seismic wavelet is an important issue in high-resolution seismic data processing and interpretation. On the basis of two patterns of constant-phase rotation and root transform for wavelet phase spectrum variation, we introduce six sparse criteria, including Lu's improved kurtosis criterion, the parsimony criterion, exponential transform criterion, Sech criterion, Cauchy criterion, and the modified Cauchy criterion, to phase spectrum estimation of the seismic wavelet, obtaining an equivalent effect to the kurtosis criterion. Through numerical experiments, we find that when the reflectivity is not a sparse sequence, the estimated phase spectrum of the seismic wavelet based on the criterion function will deviate from the true value. In order to eliminate the influence of non-sparse reflectivity series in a single trace, we apply the method to the multi-trace seismogram, improving the accuracy of seismic wavelet phase spectrum estimation.展开更多
On the assumption that the wavelet is causal and nonminimum phase,an autoregressive moving average(ARMA) model is introduced to fit the seismic trace.Seismic wavelet extraction is converted to parameters estimation of...On the assumption that the wavelet is causal and nonminimum phase,an autoregressive moving average(ARMA) model is introduced to fit the seismic trace.Seismic wavelet extraction is converted to parameters estimation of the ARMA model.Singular value decomposition(SVD) of an appropriate matrix formed by autocorrelation is exploited to determine the autoregressive(AR) order,and the cumulant-based SVD-TLS(total least squares) approach is proposed to obtain the AR parameters.The author proposes a new moving average(MA) model order determination method via combining the information theoretic criteria method and higher-order cumulant method.The cumulant approach is used to achieve the MA parameters.Theoretical analysis and numerical simulations demonstrate the feasibility of the wavelet extraction approach.展开更多
This paper puts forward wavelet transform method to identify P and S phases in three component seismograms using polarization information contained in the wavelet transform coefficients of signal. The P and S wave loc...This paper puts forward wavelet transform method to identify P and S phases in three component seismograms using polarization information contained in the wavelet transform coefficients of signal. The P and S wave locator functions are constructed by using eigenvalue analysis method to wavelet transform coefficient across several scales. Locator functions formed by wavelet transform have stated noise resistance capability, and is proved to be very effective in identifying the P and S arrivals of the test data and actual earthquake data.展开更多
This paper presents a wavelet-based approach for estimating the response of the base-isolated structure under seismic ground motions. The seismic ground motion record is expressed as the multi-scale wavelet coefficien...This paper presents a wavelet-based approach for estimating the response of the base-isolated structure under seismic ground motions. The seismic ground motion record is expressed as the multi-scale wavelet coefficients which presents the time frequency characteristics of the seismic excitation. The wavelet domain governing differential equation between the wavelet coefficients of the excitation and response is derived. Numerical study on a one-storey base isolated structure is performed. The result shows that the wavelet based response computation method is of high precision.展开更多
Based on the characteristics of gradual change style seismic signal onset which has more high frequency signal components but less magnitude, this paper selects Gauss linear frequency modulation wavelet as base functi...Based on the characteristics of gradual change style seismic signal onset which has more high frequency signal components but less magnitude, this paper selects Gauss linear frequency modulation wavelet as base function to study the change characteristics of Gauss linear frequency modulation wavelet transform with difference wavelet and signal parameters, analyzes the error origin of seismic phases identification on the basis of Gauss linear frequency modulation wavelet transform, puts forward a kind of new method identifying gradual change style seismic phases with background noise which is called fixed scale wavelet transform ratio, and presents application examples about simulation digital signal and actual seismic phases recording onsets identification.展开更多
The field seismic data is disturbed by the interferential information, which has low signal to noise ratio (SNR). That is disadvantage for seismic data interpretation. So it is important to remove the noise of seismic...The field seismic data is disturbed by the interferential information, which has low signal to noise ratio (SNR). That is disadvantage for seismic data interpretation. So it is important to remove the noise of seismic data. Independent component analysis (ICA) can remove most of the noise interference. However, ICA has some defects in noise reduction, because it needs some conditions that seismic data is independent reciprocally for denoising. To solve these defects, this paper proposes an improved ICA algorithm to noise reduction. Through simulation experiments, it can be obtained that the best decomposition levels of the new algorithm is 3. At last, the proposed improved ICA is applied to deal with the actual seismic data. The results show that it can effectively eliminate most of seismic noise such as random noise, linear interference, surface waves, and so on. The improved ICA is not only easy to denoising, but also has excellent mathematical theoretical properties.展开更多
In the time-frequency analysis of seismic signals, the matching pursuit algorithm is an effective tool for non-stationary signals, and has high time-frequency resolution and a transient structure with local self-adapt...In the time-frequency analysis of seismic signals, the matching pursuit algorithm is an effective tool for non-stationary signals, and has high time-frequency resolution and a transient structure with local self-adaption. We expand the time-frequency dictionary library with Ricker, Morlet, and mixed phase seismic wavelets, to make the method more suitable for seismic signal time-frequency decomposition. In this paper, we demonstrated the algorithm theory using synthetic seismic data, and tested the method using synthetic data with 25% noise. We compared the matching pursuit results of the time-frequency dictionaries. The results indicated that the dictionary which matched the signal characteristics better would obtain better results, and can reflect the information of seismic data effectively.展开更多
The construction of basic wavelet was discussed and many basic analyzing wavelets was compared. Acomplex analyzing wavelet which is continuous, smoothing, orthogonal and exponential decreasing was presented, andit was...The construction of basic wavelet was discussed and many basic analyzing wavelets was compared. Acomplex analyzing wavelet which is continuous, smoothing, orthogonal and exponential decreasing was presented, andit was used to decompose two blasting seismic signals with the continuous wavelet transforms (CWT). The resultshows that wavelet analysis is the better method to help us determine the essential factors which create damage effectsthan Fourier analysis.展开更多
Non-liner wave equation inversion,wavelet analysis and artificial neural networks were used to obtain stratum parameters and the distribution of thin coal seams.The lithology of the water-bearing/resisting layer in th...Non-liner wave equation inversion,wavelet analysis and artificial neural networks were used to obtain stratum parameters and the distribution of thin coal seams.The lithology of the water-bearing/resisting layer in the Quaternary system was also predicted.The implementation process included calculating the well log parameters,stratum contrasting the seismic data and the well logs,and extracting,studying and predicting seismic attributes.Seismic inversion parameters,including the layer velocity and wave impedance,were calculated and effectively used for prediction and analysis.Prior knowledge and seismic interpretation were used to remedy a dearth of seismic data during the inversion procedure.This enhanced the stability of the inversion method.Non-linear seismic inversion and artificial neural networks were used to interpret coal seismic lithology and to study the water-bearing/resisting layer in the Quaternary system.Interpretation of the 1～2 m thin coal seams,and also of the water-bearing/resisting layer in the Quaternary system,is provided.The upper mining limit can be lifted from 60 m to 45 m.The predictions show that this method can provide reliable data useful for thin coal seam exploitation and for lifting the upper mining limit,which is one of the principles of green mining.展开更多
The coherence cube technology has become an important technology for the seismic attribute interpretation, which extracts the discontinuities of the events through analyzing the similarities of adjacent seismic channe...The coherence cube technology has become an important technology for the seismic attribute interpretation, which extracts the discontinuities of the events through analyzing the similarities of adjacent seismic channels to identify the fault form. The coherence cube technology which uses constant time window lengths can not balance the shallow layers and the deep layers, because the frequency band of seismic data varies with time. When analyzing the shallow layers, the time window will crossover a lot of events, which will lead to weak focusing ability and failure to delineate the details. While the time window will not be long enough for analyzing deep layers, which will lead to low accuracy because the coherences near the zero points of the events are heavily influenced by noise. For solving the problem, we should make a research on the coherence cube technology with self-adaptive time window. This paper determines the sample points' time window lengths in real time by computing the instantaneous frequency bands with Wavelet Transformation, which gives a coherence computing method with the self-adaptive time window lengths. The result shows that the coherence cube technology with self-adaptive time window based on Wavelet Transformation improves the accuracy of fault identification, and supresses the noise effectively. The method combines the advantages of long time window method and short time window method.展开更多
In order to solve the problems of the fine division of sedimentary sequence cycles and their change in two-dimensional space as well as lateral extension contrast, we developed a method of wavelet depth-frequency anal...In order to solve the problems of the fine division of sedimentary sequence cycles and their change in two-dimensional space as well as lateral extension contrast, we developed a method of wavelet depth-frequency analysis. The single signal and composite signal of different Milankovitch cycles are obtained by numerical simulation. The simulated composite signal can be separated into single signals of a single frequency cycle. We also develop a well-seismic calibration insertion technology which helps to realize the calibration from the spectrum characteristics of a single well to the seismic profile. And then we determine the change and distribution characteristics of spectrum cycles in the two-dimensional space. It points out the direction in determining the variations of the regional sedimentary sequence cycles, underground strata structure and the contact relationship.展开更多
文摘The practice of exploration and production has proved that explosives are excited in different surrounding rocks and the seismic wavelets collected have different characteristics. In this paper, by establishing a numerical model of the explosion in the well, using finite element analysis technology for numerical simulation, the simulation calculated the stress structure in the near-source area of the earthquake excitation, and extracted the seismic wavelet. The results show that the simulation seismic wavelet characteristics of different thin interbedded sand and mudstone structures have changed significantly. Through excitation simulation, the amplitude and spectrum information of seismic wavelets can be compared and analyzed, and the excitation parameters can be optimized. .
基金supported by the Major Basic Research Development Program of China (973 Project No. 2007CB209608)
文摘Phase spectrum estimation of the seismic wavelet is an important issue in high-resolution seismic data processing and interpretation. On the basis of two patterns of constant-phase rotation and root transform for wavelet phase spectrum variation, we introduce six sparse criteria, including Lu's improved kurtosis criterion, the parsimony criterion, exponential transform criterion, Sech criterion, Cauchy criterion, and the modified Cauchy criterion, to phase spectrum estimation of the seismic wavelet, obtaining an equivalent effect to the kurtosis criterion. Through numerical experiments, we find that when the reflectivity is not a sparse sequence, the estimated phase spectrum of the seismic wavelet based on the criterion function will deviate from the true value. In order to eliminate the influence of non-sparse reflectivity series in a single trace, we apply the method to the multi-trace seismogram, improving the accuracy of seismic wavelet phase spectrum estimation.
基金supported by the National High Technology Research and Development Program of China (863 Program, No.2007AA09Z301)the Graduate Innovation Fund of China University of Petroleum and National Natural Science Foundation of China (40974072)
文摘On the assumption that the wavelet is causal and nonminimum phase,an autoregressive moving average(ARMA) model is introduced to fit the seismic trace.Seismic wavelet extraction is converted to parameters estimation of the ARMA model.Singular value decomposition(SVD) of an appropriate matrix formed by autocorrelation is exploited to determine the autoregressive(AR) order,and the cumulant-based SVD-TLS(total least squares) approach is proposed to obtain the AR parameters.The author proposes a new moving average(MA) model order determination method via combining the information theoretic criteria method and higher-order cumulant method.The cumulant approach is used to achieve the MA parameters.Theoretical analysis and numerical simulations demonstrate the feasibility of the wavelet extraction approach.
基金supported by National Key Basic Research Development Program (Grant No. 2007CB209600)National Major Science and Technology Program (Grant No. 2008ZX05010-002)
文摘This paper puts forward wavelet transform method to identify P and S phases in three component seismograms using polarization information contained in the wavelet transform coefficients of signal. The P and S wave locator functions are constructed by using eigenvalue analysis method to wavelet transform coefficient across several scales. Locator functions formed by wavelet transform have stated noise resistance capability, and is proved to be very effective in identifying the P and S arrivals of the test data and actual earthquake data.
文摘This paper presents a wavelet-based approach for estimating the response of the base-isolated structure under seismic ground motions. The seismic ground motion record is expressed as the multi-scale wavelet coefficients which presents the time frequency characteristics of the seismic excitation. The wavelet domain governing differential equation between the wavelet coefficients of the excitation and response is derived. Numerical study on a one-storey base isolated structure is performed. The result shows that the wavelet based response computation method is of high precision.
基金State Natural Science Foundation of China (40074007) Science and Technology Key Project during the Ten-Year Plan(2001BA601B02-03-06) and the Natural Science Foundation of Shandong Province (Y2000E08).
文摘Based on the characteristics of gradual change style seismic signal onset which has more high frequency signal components but less magnitude, this paper selects Gauss linear frequency modulation wavelet as base function to study the change characteristics of Gauss linear frequency modulation wavelet transform with difference wavelet and signal parameters, analyzes the error origin of seismic phases identification on the basis of Gauss linear frequency modulation wavelet transform, puts forward a kind of new method identifying gradual change style seismic phases with background noise which is called fixed scale wavelet transform ratio, and presents application examples about simulation digital signal and actual seismic phases recording onsets identification.
基金Funded by the Project of China Geological Survey (No.1212010916040)the Sichuan Science and Technology Program (No.2017JY0051)the Sichuan Science and Technology Program (No.2018GZ0200)
文摘The field seismic data is disturbed by the interferential information, which has low signal to noise ratio (SNR). That is disadvantage for seismic data interpretation. So it is important to remove the noise of seismic data. Independent component analysis (ICA) can remove most of the noise interference. However, ICA has some defects in noise reduction, because it needs some conditions that seismic data is independent reciprocally for denoising. To solve these defects, this paper proposes an improved ICA algorithm to noise reduction. Through simulation experiments, it can be obtained that the best decomposition levels of the new algorithm is 3. At last, the proposed improved ICA is applied to deal with the actual seismic data. The results show that it can effectively eliminate most of seismic noise such as random noise, linear interference, surface waves, and so on. The improved ICA is not only easy to denoising, but also has excellent mathematical theoretical properties.
文摘In the time-frequency analysis of seismic signals, the matching pursuit algorithm is an effective tool for non-stationary signals, and has high time-frequency resolution and a transient structure with local self-adaption. We expand the time-frequency dictionary library with Ricker, Morlet, and mixed phase seismic wavelets, to make the method more suitable for seismic signal time-frequency decomposition. In this paper, we demonstrated the algorithm theory using synthetic seismic data, and tested the method using synthetic data with 25% noise. We compared the matching pursuit results of the time-frequency dictionaries. The results indicated that the dictionary which matched the signal characteristics better would obtain better results, and can reflect the information of seismic data effectively.
文摘The construction of basic wavelet was discussed and many basic analyzing wavelets was compared. Acomplex analyzing wavelet which is continuous, smoothing, orthogonal and exponential decreasing was presented, andit was used to decompose two blasting seismic signals with the continuous wavelet transforms (CWT). The resultshows that wavelet analysis is the better method to help us determine the essential factors which create damage effectsthan Fourier analysis.
基金Projects 40574057 and 40874054 supported by the National Natural Science Foundation of ChinaProjects 2007CB209400 by the National Basic Research Program of ChinaFoundation of China University of Mining and Technology (OF4471)
文摘Non-liner wave equation inversion,wavelet analysis and artificial neural networks were used to obtain stratum parameters and the distribution of thin coal seams.The lithology of the water-bearing/resisting layer in the Quaternary system was also predicted.The implementation process included calculating the well log parameters,stratum contrasting the seismic data and the well logs,and extracting,studying and predicting seismic attributes.Seismic inversion parameters,including the layer velocity and wave impedance,were calculated and effectively used for prediction and analysis.Prior knowledge and seismic interpretation were used to remedy a dearth of seismic data during the inversion procedure.This enhanced the stability of the inversion method.Non-linear seismic inversion and artificial neural networks were used to interpret coal seismic lithology and to study the water-bearing/resisting layer in the Quaternary system.Interpretation of the 1～2 m thin coal seams,and also of the water-bearing/resisting layer in the Quaternary system,is provided.The upper mining limit can be lifted from 60 m to 45 m.The predictions show that this method can provide reliable data useful for thin coal seam exploitation and for lifting the upper mining limit,which is one of the principles of green mining.
文摘The coherence cube technology has become an important technology for the seismic attribute interpretation, which extracts the discontinuities of the events through analyzing the similarities of adjacent seismic channels to identify the fault form. The coherence cube technology which uses constant time window lengths can not balance the shallow layers and the deep layers, because the frequency band of seismic data varies with time. When analyzing the shallow layers, the time window will crossover a lot of events, which will lead to weak focusing ability and failure to delineate the details. While the time window will not be long enough for analyzing deep layers, which will lead to low accuracy because the coherences near the zero points of the events are heavily influenced by noise. For solving the problem, we should make a research on the coherence cube technology with self-adaptive time window. This paper determines the sample points' time window lengths in real time by computing the instantaneous frequency bands with Wavelet Transformation, which gives a coherence computing method with the self-adaptive time window lengths. The result shows that the coherence cube technology with self-adaptive time window based on Wavelet Transformation improves the accuracy of fault identification, and supresses the noise effectively. The method combines the advantages of long time window method and short time window method.
基金supported by the Fundamental Research Funds for the Central Universities(Grant No.2011YYL128)the CNPC Innovation Foundation(GrantNo.2012D-5006-0103)the Ministry of Land and Resources special funds for scientific research on public cause(Grant No.201311107)
文摘In order to solve the problems of the fine division of sedimentary sequence cycles and their change in two-dimensional space as well as lateral extension contrast, we developed a method of wavelet depth-frequency analysis. The single signal and composite signal of different Milankovitch cycles are obtained by numerical simulation. The simulated composite signal can be separated into single signals of a single frequency cycle. We also develop a well-seismic calibration insertion technology which helps to realize the calibration from the spectrum characteristics of a single well to the seismic profile. And then we determine the change and distribution characteristics of spectrum cycles in the two-dimensional space. It points out the direction in determining the variations of the regional sedimentary sequence cycles, underground strata structure and the contact relationship.
