抛物线拉冬变换消除多次波的应用要素分析
|
摘要
抛物线拉冬变换(PRT)在国外已是消除多次波的最有效方法之一,但国内却很少应用。其中一个重要原因就是对抛物线拉冬变换的物理意义、适用条件、算法和参数的选择等关键问题缺乏足够的认识。本文在模型研究的基础上,通过理论模型试验,分别讨论了时间域拉冬变换的特点和存在的问题、频率域拉冬变换的不同实现方法及其应用效果、拉冬域的假频问题、复系数矩阵方程的求解和白噪系数的选取、拉冬域多次波切除和截断效应、自适应滤波切除等应用要素,获得了一些有指导意义的结论1时间域PRT不宜直接使用;2频率域正变换最小平方法有利于消除多次波;3拉冬域的参数取值要满足采样定理;4求解复系数方程组白噪系数宜取0.1～1;5加一巴特沃兹滤波器式的增益函数可实现自适应滤除多次波。应用本文方法对实际资料进行处理,取得了很好的应用效果。
Parabolic Radon transform is one of the effective methods for removing multiple abroad,but unfortunately is less used at home.One of the important reasons is lack of enough knowledge about the key problems such as physical meaning of parabolic Radon transform,its adaptive conditions,algorithm and selection of parameters etc. On the basis of model study and through theoretical model test,the paper separately discussed the characters and existed problems of Radon transform in time domain,different realizing ways and applied effects of Radon transform in frequency domain,alias of Radon domain,solution of complex coefficient matrix equation and selection of white coefficient,mute and truncation effects of multiple in Radon domain and self adaptive filtering mute etc.,acquired some useful conclusions:①PRT in time domain can't be directly used;②the least square method of forward transform in frequency domain is of benefit to removing multiple;③the parameter values in Radon domain should meet the sampling theorem;④the selected white coefficient of solving complex coefficient equation is in a region of 0.1~1.0;⑤to add a gain function similar to Butterworth filter can realize self adaptive filtering of multiple.Using the method presented in the paper for practical data processing achieved a good applied results.
Parabolic Radon transform is one of the effective methods for removing multiple abroad,but unfortunately is less used at home.One of the important reasons is lack of enough knowledge about the key problems such as physical meaning of parabolic Radon transform,its adaptive conditions,algorithm and selection of parameters etc. On the basis of model study and through theoretical model test,the paper separately discussed the characters and existed problems of Radon transform in time domain,different realizing ways and applied effects of Radon transform in frequency domain,alias of Radon domain,solution of complex coefficient matrix equation and selection of white coefficient,mute and truncation effects of multiple in Radon domain and self adaptive filtering mute etc.,acquired some useful conclusions:①PRT in time domain can't be directly used;②the least square method of forward transform in frequency domain is of benefit to removing multiple;③the parameter values in Radon domain should meet the sampling theorem;④the selected white coefficient of solving complex coefficient equation is in a region of 0.1~1.0;⑤to add a gain function similar to Butterworth filter can realize self adaptive filtering of multiple.Using the method presented in the paper for practical data processing achieved a good applied results.
引文
[1]Ham psonD.Inverse velocity stacking and m ultipleelim ination.Journal of theCanadianSociety ofExplorationGeophysicists,1986,22:44～55
[2]FosterJ andMosherC.Suppression of multiplereflections using theRadon transform.Geophysics,1992,57:386～395
[3]SacchiM D andUlrychT J.Im proving resolution ofRadon operators using a model reweighed least-squares procedure.Journal ofSeismicExploration,1995,4:315～328
[4]ZhouB andGreenhalghS A.Multiple suppression by2D filtering in the parabolic domain:a wave- equation-based m ethod.GeophysicalProspecting,1996,44:375～401
[5]KabirM M N andMarfurtK J.Toward trueam plitude m ultiple removal.TheL eadingEdge,1999,18:66～73
[6]L andaE,BelferI andKeydarS.Multiple attenuationin the parabolic domain using wavefrontcharacteristics of multiple generating prim aries.Geophysics,1999,64:1806～1815
[7]沈操,牛滨华,余钦范.拉冬变换的MATL AB实现.物探化探计算技术,2000,22(4):346～350,368
[8]朱生旺,魏修成,李锋等.用抛物线Radon变换稀疏解分离和压制多次波.石油地球物理勘探,2002,37(2):110～115
[9]WangY.Multiple attenuation:coping with the spatialtruncation effect in theRadon transform domain.GeophysicalProspecting,2003,51:75～87
[10]李振春,张军华.地震资料数字处理讲义.山东东营:石油大学出版社,2003
[11]MarfurtK J,SchneiderR V andMuellerM C.Pitfalls of using conventional and discreteRadontransforms on poorly sam pled data.Geophysics,1996,61:1467～1482
[12]牛滨华,孙春岩,张中杰等.多项式Radon变换.地球物理学报,2001,44(2):263～271
[13]TurnerG.Aliasing in the tau- p transform and theremove of spatial aliased coherent noise.Geophysics,1990,55:1496～1503
[14]徐士良.Fortran常用算法程序集,北京:清华大学出版社,1992
[2]FosterJ andMosherC.Suppression of multiplereflections using theRadon transform.Geophysics,1992,57:386～395
[3]SacchiM D andUlrychT J.Im proving resolution ofRadon operators using a model reweighed least-squares procedure.Journal ofSeismicExploration,1995,4:315～328
[4]ZhouB andGreenhalghS A.Multiple suppression by2D filtering in the parabolic domain:a wave- equation-based m ethod.GeophysicalProspecting,1996,44:375～401
[5]KabirM M N andMarfurtK J.Toward trueam plitude m ultiple removal.TheL eadingEdge,1999,18:66～73
[6]L andaE,BelferI andKeydarS.Multiple attenuationin the parabolic domain using wavefrontcharacteristics of multiple generating prim aries.Geophysics,1999,64:1806～1815
[7]沈操,牛滨华,余钦范.拉冬变换的MATL AB实现.物探化探计算技术,2000,22(4):346～350,368
[8]朱生旺,魏修成,李锋等.用抛物线Radon变换稀疏解分离和压制多次波.石油地球物理勘探,2002,37(2):110～115
[9]WangY.Multiple attenuation:coping with the spatialtruncation effect in theRadon transform domain.GeophysicalProspecting,2003,51:75～87
[10]李振春,张军华.地震资料数字处理讲义.山东东营:石油大学出版社,2003
[11]MarfurtK J,SchneiderR V andMuellerM C.Pitfalls of using conventional and discreteRadontransforms on poorly sam pled data.Geophysics,1996,61:1467～1482
[12]牛滨华,孙春岩,张中杰等.多项式Radon变换.地球物理学报,2001,44(2):263～271
[13]TurnerG.Aliasing in the tau- p transform and theremove of spatial aliased coherent noise.Geophysics,1990,55:1496～1503
[14]徐士良.Fortran常用算法程序集,北京:清华大学出版社,1992
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心