基于反射地震记录变子波模型提高地震记录分辨率
|
摘要
本文给出了地震记录变子波模型的一种近似数学表达式.基于该表达式研究了反射系数序列不满足白噪假设和子波在地下传播时发生变化这两种情况下地震道谱的组成及结构,讨论了谱白化及反褶积方法在这两种情况下效果不佳的原因.然后基于变子波模型,提出了一种新的提高地震记录分辨率的方法:第一步,用自适应于地震记录的Gabor分子窗把地震记录恰当地划分成若干片断,每段内信号近似平稳,然后将地震记录变换到时间—频率域;第二步,在变换域对每个分子窗内信号的振幅谱进行处理以拓宽频带;最后把处理后的时间—频率域函数反变换回时间域得到提高分辨率后的结果.本文提出的方法具有能较好地适用于反射系数不满足白噪假设的情况及提高分辨率后的地震记录能较好地保持原地震记录的相对能量关系等优点,模型和实际资料算例结果均表明,本文方法在拓宽地震资料频带及保持地震记录局部能量相对关系方面均明显优于谱白化方法.
We present a mathematic expression of the changing wavelet model of the seismogram (CWMS).Based on this expression,we study the content and structure of spectra of seismic traces produced in two cases,i.e.the reflectivity does not meet the white assumption and seismic wavelets change during travel in subsurface.And we analyze why the spectral-whitening and deconvolution methods cannot yield good results in these two cases.Then based on CWMS,we propose a new method to enhance the resolution of seismic traces.Firstly,we divide the seismic signal into proper segments by Gabor molecule-windows adaptive to the signal,where each segment is regarded stationary,and then transform it into the Gabor domain.Secondly,we deal with the amplitude spectrum for each segment to expand the frequency band of the spectrum in the Gabor domain.Finally,taking inverse Gabor transform on modified time-frequency coefficients,we get the high-resolution signal.Compared with the spectral-whitening and deconvolution methods,our method have the following advantages:(1) it can work well even when the reflectivity doesn't satisfy white presumption;(2) it can preserve local energy relationship with original seismic trace which may be important for seismic data interpretation. Tests on synthetic and real data show that our method is much better than the spectral-whitening method in expanding frequency bandwidth and preserving local energy relationship of seismograms.
We present a mathematic expression of the changing wavelet model of the seismogram (CWMS).Based on this expression,we study the content and structure of spectra of seismic traces produced in two cases,i.e.the reflectivity does not meet the white assumption and seismic wavelets change during travel in subsurface.And we analyze why the spectral-whitening and deconvolution methods cannot yield good results in these two cases.Then based on CWMS,we propose a new method to enhance the resolution of seismic traces.Firstly,we divide the seismic signal into proper segments by Gabor molecule-windows adaptive to the signal,where each segment is regarded stationary,and then transform it into the Gabor domain.Secondly,we deal with the amplitude spectrum for each segment to expand the frequency band of the spectrum in the Gabor domain.Finally,taking inverse Gabor transform on modified time-frequency coefficients,we get the high-resolution signal.Compared with the spectral-whitening and deconvolution methods,our method have the following advantages:(1) it can work well even when the reflectivity doesn't satisfy white presumption;(2) it can preserve local energy relationship with original seismic trace which may be important for seismic data interpretation. Tests on synthetic and real data show that our method is much better than the spectral-whitening method in expanding frequency bandwidth and preserving local energy relationship of seismograms.
引文
[1] Yilmaz O.Seismic data processing.Tulsa:Society of Exploration Geophysicists,1987． 49～61,85～94
[2] Ricker N.The form and nature of seismic wavelets and the structure of seismograms.Geophysics,1940,5:348～366
[3] Dobrin M B,Savit C H.Introduction to geophysical prospecting,4~(th)ed.New York:McGraw Hill Book Co,1988
[4] Ziolkowski A.Why don't we measure seismic signatures?Geophysics,1991,56:190～201
[5] Grchenig K.Foundations of Time-Frequency Analysis.Boston:Birkhauser,2000． 301～307
[6] Margrave G F.Theory of nonstationary linear filtering in the Fourier domain with application to time variant filtering.Geophysics,1998,63(1) :244～259
[ 7 ] Kjartannson E. Constant Q-wave propagation and attenuation. J. Geophys. Res. , 1979 ,84:4737 ～ 4748
[ 8 ] Hale T D. Q-adaptive deconvolution. Stanford Exploration Project Report Number 30,1982． 133 ～ 158
[ 9 ] Gibson B, Larner K L. Comparison of spectral flattening techniques. unpublished technical document. Western Geophysical Company, 1982
[10] Margrave G F, Lamoureux M P. Gabor deconvolution. CREWES Research Report, 2001,13:241～276, 252～253
[11] Robinson E A, Treitel S. Principles of digital Wiener filtering. Geophys. Prosp. ,1967,15:311～333
[12] Margrave G F, Henley D C, Lamoureux M P, et al. An update on Gabor deconvolution. CREWES Research Report, 2002,14:1 ～ 27
[13] Margrave G F, Gibson P C, Grossman J P, et al. The Gabor transform, pseudodifferential operators, and seismic deconvolution. Integrated Computer-Aided Engineering, 2005,12:43 ～ 45
[14] Goupillaud P L. An approach to inverse filtering of near surface layer effects from seismic records. Geo physics, 1961,26:754 ～ 760
[15] Girod B. Signal and Systems. Chichester: John Wiley & Sons Ltd, 2001． 261～309
[16] Barnes A E. Instantaneous frequency and amplitude at the envelope peak of a constant-phase wavelet. Geo physics ,1991,56:1058 ～ 1060
[17] Coppens F. Spectrum equavilente; A method to enhance the quality of seismic data. Geophysical Prospecting , 1984, 32(2) :157 ～ 350
[18] Grossman J P, Margrave G F, Lamoureux M P.Constructing adaptive, nonuniform Gabor frames from partitions of unity. CREWES Research Report ,2002,14:1 ～ 10
[19] Grossman J P. Theory of adaptive, nonstationary filtering in the Gabor domain with applications to seismic inversion [Ph. D. thesis]. Calgary: Department of Geology and Geophysics and Department of Mathematics and Statistics of University of Calgary, 2005
[20] Margrave G F, Lamoureux M P, Grossman J P, et al. Gabor deconvolution of seismic data for source waveform and Q correction. 72nd Annual International Meeting, Soc. of Expl. Geophys. , 2002． 2190 ～ 2193
[2] Ricker N.The form and nature of seismic wavelets and the structure of seismograms.Geophysics,1940,5:348～366
[3] Dobrin M B,Savit C H.Introduction to geophysical prospecting,4~(th)ed.New York:McGraw Hill Book Co,1988
[4] Ziolkowski A.Why don't we measure seismic signatures?Geophysics,1991,56:190～201
[5] Grchenig K.Foundations of Time-Frequency Analysis.Boston:Birkhauser,2000． 301～307
[6] Margrave G F.Theory of nonstationary linear filtering in the Fourier domain with application to time variant filtering.Geophysics,1998,63(1) :244～259
[ 7 ] Kjartannson E. Constant Q-wave propagation and attenuation. J. Geophys. Res. , 1979 ,84:4737 ～ 4748
[ 8 ] Hale T D. Q-adaptive deconvolution. Stanford Exploration Project Report Number 30,1982． 133 ～ 158
[ 9 ] Gibson B, Larner K L. Comparison of spectral flattening techniques. unpublished technical document. Western Geophysical Company, 1982
[10] Margrave G F, Lamoureux M P. Gabor deconvolution. CREWES Research Report, 2001,13:241～276, 252～253
[11] Robinson E A, Treitel S. Principles of digital Wiener filtering. Geophys. Prosp. ,1967,15:311～333
[12] Margrave G F, Henley D C, Lamoureux M P, et al. An update on Gabor deconvolution. CREWES Research Report, 2002,14:1 ～ 27
[13] Margrave G F, Gibson P C, Grossman J P, et al. The Gabor transform, pseudodifferential operators, and seismic deconvolution. Integrated Computer-Aided Engineering, 2005,12:43 ～ 45
[14] Goupillaud P L. An approach to inverse filtering of near surface layer effects from seismic records. Geo physics, 1961,26:754 ～ 760
[15] Girod B. Signal and Systems. Chichester: John Wiley & Sons Ltd, 2001． 261～309
[16] Barnes A E. Instantaneous frequency and amplitude at the envelope peak of a constant-phase wavelet. Geo physics ,1991,56:1058 ～ 1060
[17] Coppens F. Spectrum equavilente; A method to enhance the quality of seismic data. Geophysical Prospecting , 1984, 32(2) :157 ～ 350
[18] Grossman J P, Margrave G F, Lamoureux M P.Constructing adaptive, nonuniform Gabor frames from partitions of unity. CREWES Research Report ,2002,14:1 ～ 10
[19] Grossman J P. Theory of adaptive, nonstationary filtering in the Gabor domain with applications to seismic inversion [Ph. D. thesis]. Calgary: Department of Geology and Geophysics and Department of Mathematics and Statistics of University of Calgary, 2005
[20] Margrave G F, Lamoureux M P, Grossman J P, et al. Gabor deconvolution of seismic data for source waveform and Q correction. 72nd Annual International Meeting, Soc. of Expl. Geophys. , 2002． 2190 ～ 2193
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心