Wavelet transforms have been successfully used in seismic data processing with their ability for local time - frequency analysis. However, identification of directionality is limited because wavelet transform coeffici...Wavelet transforms have been successfully used in seismic data processing with their ability for local time - frequency analysis. However, identification of directionality is limited because wavelet transform coefficients reveal only three spatial orientations. Whereas the ridgelet transform has a superior capability for direction detection and the ability to process signals with linearly changing characteristics. In this paper, we present the issue of low signal-to-noise ratio (SNR) seismic data processing based on the ridgelet transform. Actual seismic data with low SNR from south China has been processed using ridgelet transforms to improve the SNR and the continuity of seismic events. The results show that the ridgelet transform is better than the wavelet transform for these tasks.展开更多
Conventional f-x prediction filtering methods are based on an autoregressive model. The error section is first computed as a source noise but is removed as additive noise to obtain the signal, which results in an assu...Conventional f-x prediction filtering methods are based on an autoregressive model. The error section is first computed as a source noise but is removed as additive noise to obtain the signal, which results in an assumption inconsistency before and after filtering. In this paper, an autoregressive, moving-average model is employed to avoid the model inconsistency. Based on the ARMA model, a noncasual prediction filter is computed and a self-deconvolved projection filter is used for estimating additive noise in order to suppress random noise. The 1-D ARMA model is also extended to the 2-D spatial domain, which is the basis for noncasual spatial prediction filtering for random noise attenuation on 3-D seismic data. Synthetic and field data processing indicate this method can suppress random noise more effectively and preserve the signal simultaneously and does much better than other conventional prediction filtering methods.展开更多
In seismic exploration for coal, seismic waves are very difficult to transmit downward because of high velocity protective layers, making the reflection information very hard to receive above ground. Based on the Snel...In seismic exploration for coal, seismic waves are very difficult to transmit downward because of high velocity protective layers, making the reflection information very hard to receive above ground. Based on the Snell law and the Zoeppritz equation, we studied the relationship between the incidence angle and reflection seismic wave energy using a forward model of level media. The result shows that the seismic wave energy has a sudden increase at the critical angle. Based on the energy propagation rule, using big offset to receive the seismic wave energy under a protective layer can effectively reduce its protection effect.展开更多
Converted waves have slow velocity and low signal-to-noise ratio. It is also difficult to pick first-breaks and bin the common-conversion-points (CCP). Some statics methods, which work well for P-wave data, can't b...Converted waves have slow velocity and low signal-to-noise ratio. It is also difficult to pick first-breaks and bin the common-conversion-points (CCP). Some statics methods, which work well for P-wave data, can't be effectively used for solving convertedwave statics problems. This has become the main obstacle to breakthroughs in convertedwave data processing. To improve converted-wave static corrections, first, a statics method based on the common-receiver-point (CRP) stack is used for the initial receiver static correction to enhance the coherency of the CRP stack. Second, a stack-power-maximization static correction which improves the continuity of the CCP stack is used for detailed receiver statics. Finally, a non-surface-consistent residual moveout correction of the CCP gathers is used to enhance the stack power of reflection signals from different depths. Converted-wave statics are solved by the joint use of the three correction methods.展开更多
The modification to the matrix method for constructing the displacement field on the free surface of an anisotropic layered medium was presented. The source of seismic waves was modelled by a randomly oriented force a...The modification to the matrix method for constructing the displacement field on the free surface of an anisotropic layered medium was presented. The source of seismic waves was modelled by a randomly oriented force and seismic tensor. A trial and error method was presented for solving the inverse problem of determining parameters of the earthquake source. A number of analytical and numerical approaches to determining the earthquake source parameters, based on the direct problem solutions, were proposed. The focal mechanisms for the events in the Carpathian region of Ukraine are determined by the graphical method. The theory of determinating the angles of orientation of the fault plane and the earthquake's focal mechanism was presented. The focal mechanisms obtained by two different methods were compared.展开更多
Based on the earthquake catalog reported by the Chinese digital seismic network in recent years,we select the earthquakes with both surface wave magnitude and local magnitude and fit them into a relationship between t...Based on the earthquake catalog reported by the Chinese digital seismic network in recent years,we select the earthquakes with both surface wave magnitude and local magnitude and fit them into a relationship between the two magnitudes.The systematic difference is found from the formula which has been used for 30 years.Because of a large dynamic range and wide frequency range of the current digital observation system,in addition to a larger number of stations and earthquakes being used compared to before,the relation obtained in this paper seems more reliable.Our calculation shows that there is no significant difference before and after magnitude conversion so we suggest the abandonment of magnitude conversion.The site response of a station consists of amplification at different frequencies.The amplification is equal to about 1 and changes little with frequency at stations located on basement rock,and it is greater than 1 at low frequency ranges and less than 1 at high frequency ranges at stations located on sediment layers.The difference between magnitudes from single station located on sediment layer and the average magnitude from the whole network increases from negative to positive with period.It seems that there is no fixed station correction factor and the station correction method does not work to improve the accuracy and magnitude estimates.展开更多
In process of seismic exploration,the noise of seismic signals produces serious interference. Conventional methods of wavelet threshold denoising cannot fully use the characteristics of seismic signals due to its limi...In process of seismic exploration,the noise of seismic signals produces serious interference. Conventional methods of wavelet threshold denoising cannot fully use the characteristics of seismic signals due to its limitations. There is always a certain degree of deviation between estimated value and actual value. In this study,a method of seismic data denoising is proposed,the authors use the current coefficients,the parent coefficients and the neighborhood coefficients based on dual-tree complex wavelet transform( DTCWT) and related sub-band denoising model( TrivaS hrink) to achieve the optimal estimation of shrinking factor and get the noise reduction of seismic records. It is found that the method is better than conventional methods of wavelet threshold denoising in removing random noise.展开更多
The Quality factor is the parameter that can be used to describe the energy attenuation on seismic wave. In theory, we can obtain the relationship between the change of the coda wave quality factor with time and the s...The Quality factor is the parameter that can be used to describe the energy attenuation on seismic wave. In theory, we can obtain the relationship between the change of the coda wave quality factor with time and the strong earthquake preparation process on the basis of the quality factor of a coda wave in a same ray path. However, in reality the coda wave quality factor measured by different seismic coda waves corresponds to different seismic wave ray paths. The change of the quality factor with time is related to non-elastic characteristics of the medium and the volume of scattering ellipsoid constrained by scattered wave phase fronts, besides the change of regional stress field. This paper discusses the relationship between quality factor, epicenter distance and different lapse time, and then discusses the relationship between quality factor and frequency. Furthermore the determination method of the coda wave quality factor is put forward. The improved determination method of the quality factor, which removes the influence of different earthquakes or propagation depth of scattered waves, may increase measurement precision, thus information pertaining to abnormal changes in quality factor and the relationship between the quality factor and earthquake preparation process can be acquired.展开更多
基金This paper is supported by China Petrochemical Key Project in the"11th Five-Year"Plan Technology and the Doctorate Fund of Ministry of Education of China (No.20050491504)
文摘Wavelet transforms have been successfully used in seismic data processing with their ability for local time - frequency analysis. However, identification of directionality is limited because wavelet transform coefficients reveal only three spatial orientations. Whereas the ridgelet transform has a superior capability for direction detection and the ability to process signals with linearly changing characteristics. In this paper, we present the issue of low signal-to-noise ratio (SNR) seismic data processing based on the ridgelet transform. Actual seismic data with low SNR from south China has been processed using ridgelet transforms to improve the SNR and the continuity of seismic events. The results show that the ridgelet transform is better than the wavelet transform for these tasks.
基金This research was financially supported by National Natural Science Foundation of China (Grant No. 40604016) and the National Hi-Tech Research and Development Program (863 Program) (Grants No. 2006AA09A102-09 and No. 2007AA06Z229).
文摘Conventional f-x prediction filtering methods are based on an autoregressive model. The error section is first computed as a source noise but is removed as additive noise to obtain the signal, which results in an assumption inconsistency before and after filtering. In this paper, an autoregressive, moving-average model is employed to avoid the model inconsistency. Based on the ARMA model, a noncasual prediction filter is computed and a self-deconvolved projection filter is used for estimating additive noise in order to suppress random noise. The 1-D ARMA model is also extended to the 2-D spatial domain, which is the basis for noncasual spatial prediction filtering for random noise attenuation on 3-D seismic data. Synthetic and field data processing indicate this method can suppress random noise more effectively and preserve the signal simultaneously and does much better than other conventional prediction filtering methods.
基金Projects 40574058 supported by the National Natural Science Foundation of China2005cb221500 by the National Basic Research Program of China 03(2007) by the Scientific and Technological Project about Geology and Mineral Resources of Henan Land Resources Department
文摘In seismic exploration for coal, seismic waves are very difficult to transmit downward because of high velocity protective layers, making the reflection information very hard to receive above ground. Based on the Snell law and the Zoeppritz equation, we studied the relationship between the incidence angle and reflection seismic wave energy using a forward model of level media. The result shows that the seismic wave energy has a sudden increase at the critical angle. Based on the energy propagation rule, using big offset to receive the seismic wave energy under a protective layer can effectively reduce its protection effect.
文摘Converted waves have slow velocity and low signal-to-noise ratio. It is also difficult to pick first-breaks and bin the common-conversion-points (CCP). Some statics methods, which work well for P-wave data, can't be effectively used for solving convertedwave statics problems. This has become the main obstacle to breakthroughs in convertedwave data processing. To improve converted-wave static corrections, first, a statics method based on the common-receiver-point (CRP) stack is used for the initial receiver static correction to enhance the coherency of the CRP stack. Second, a stack-power-maximization static correction which improves the continuity of the CCP stack is used for detailed receiver statics. Finally, a non-surface-consistent residual moveout correction of the CCP gathers is used to enhance the stack power of reflection signals from different depths. Converted-wave statics are solved by the joint use of the three correction methods.
文摘The modification to the matrix method for constructing the displacement field on the free surface of an anisotropic layered medium was presented. The source of seismic waves was modelled by a randomly oriented force and seismic tensor. A trial and error method was presented for solving the inverse problem of determining parameters of the earthquake source. A number of analytical and numerical approaches to determining the earthquake source parameters, based on the direct problem solutions, were proposed. The focal mechanisms for the events in the Carpathian region of Ukraine are determined by the graphical method. The theory of determinating the angles of orientation of the fault plane and the earthquake's focal mechanism was presented. The focal mechanisms obtained by two different methods were compared.
基金This project was fund by the International Cooperation Programof the Ministry of Science and Technology of China(2005DFA20980)
文摘Based on the earthquake catalog reported by the Chinese digital seismic network in recent years,we select the earthquakes with both surface wave magnitude and local magnitude and fit them into a relationship between the two magnitudes.The systematic difference is found from the formula which has been used for 30 years.Because of a large dynamic range and wide frequency range of the current digital observation system,in addition to a larger number of stations and earthquakes being used compared to before,the relation obtained in this paper seems more reliable.Our calculation shows that there is no significant difference before and after magnitude conversion so we suggest the abandonment of magnitude conversion.The site response of a station consists of amplification at different frequencies.The amplification is equal to about 1 and changes little with frequency at stations located on basement rock,and it is greater than 1 at low frequency ranges and less than 1 at high frequency ranges at stations located on sediment layers.The difference between magnitudes from single station located on sediment layer and the average magnitude from the whole network increases from negative to positive with period.It seems that there is no fixed station correction factor and the station correction method does not work to improve the accuracy and magnitude estimates.
基金Supported by the National "863" Project(No.2014AA06A605)
文摘In process of seismic exploration,the noise of seismic signals produces serious interference. Conventional methods of wavelet threshold denoising cannot fully use the characteristics of seismic signals due to its limitations. There is always a certain degree of deviation between estimated value and actual value. In this study,a method of seismic data denoising is proposed,the authors use the current coefficients,the parent coefficients and the neighborhood coefficients based on dual-tree complex wavelet transform( DTCWT) and related sub-band denoising model( TrivaS hrink) to achieve the optimal estimation of shrinking factor and get the noise reduction of seismic records. It is found that the method is better than conventional methods of wavelet threshold denoising in removing random noise.
基金sponsored by the Natural Science Foundation of Shandong Province (Y2007E09)Joint Earthquake Science Foundation (C08028)Special Application Research of Digital Seismic Wave Data ,Shangdong,China
文摘The Quality factor is the parameter that can be used to describe the energy attenuation on seismic wave. In theory, we can obtain the relationship between the change of the coda wave quality factor with time and the strong earthquake preparation process on the basis of the quality factor of a coda wave in a same ray path. However, in reality the coda wave quality factor measured by different seismic coda waves corresponds to different seismic wave ray paths. The change of the quality factor with time is related to non-elastic characteristics of the medium and the volume of scattering ellipsoid constrained by scattered wave phase fronts, besides the change of regional stress field. This paper discusses the relationship between quality factor, epicenter distance and different lapse time, and then discusses the relationship between quality factor and frequency. Furthermore the determination method of the coda wave quality factor is put forward. The improved determination method of the quality factor, which removes the influence of different earthquakes or propagation depth of scattered waves, may increase measurement precision, thus information pertaining to abnormal changes in quality factor and the relationship between the quality factor and earthquake preparation process can be acquired.