基于低频约束的高分辨率Radon变换多次波压制方法研究
|
摘要
抛物线Radon变换是地震资料处理中常用的一种多次波压制方法,常规Radon变换由于采用最小二乘算法的原因,对不同的曲率值采用相同的加权值进行求解,分辨率较低。为提高Radon域数据的分辨率,应采用对不同曲率值应用不同的加权值的方法。一般的思路是利用前一次迭代的结果得到加权矩阵,该算法需要迭代进行,运算量较大。这里拟采用低频约束的方式提高Radon变换的分辨率,即利用前一个频率的运算结果对下一个频率的计算进行约束,该算法分辨率高,计算速度快。模拟数据和实际数据的测试表明,本文方法在多次波压制处理中,可以快速有效地完成多次波去除。
Parabolic Radon transform is one of the common used method of multiple attenuation in seismic data processing. Conventional Radon transform with the least squares algorithm, adopts the same weight for different q values, thus generate low resolution data in Radon domain. In order to increase the resolution in the Radon domain, the method of giving different weight for different q value should be adopted. The general idea is to use the result of previous iteration to get the weighted matrix. The method needs to be solved iteratively, and requires large computation. We propose a method based low frequency constraint to improve the resolution in Radon domain, which use the result of previous frequency to constrain the computation of the next frequency. The method can get a high-resolution result, and can get a good result very quickly. Simulated data and real data tests demonstrated the feasibility and efficiency of the method in multiple suppression.
Parabolic Radon transform is one of the common used method of multiple attenuation in seismic data processing. Conventional Radon transform with the least squares algorithm, adopts the same weight for different q values, thus generate low resolution data in Radon domain. In order to increase the resolution in the Radon domain, the method of giving different weight for different q value should be adopted. The general idea is to use the result of previous iteration to get the weighted matrix. The method needs to be solved iteratively, and requires large computation. We propose a method based low frequency constraint to improve the resolution in Radon domain, which use the result of previous frequency to constrain the computation of the next frequency. The method can get a high-resolution result, and can get a good result very quickly. Simulated data and real data tests demonstrated the feasibility and efficiency of the method in multiple suppression.
引文
[1] HAMPSON D.Inverse velocity stacking for multiple elimination[J].Journal of the Canadian Society of Exploration Geophysics,1986(22):44-45.
[2] BEYLKIN G.The discrete Radon transform[J].IEEE Transactions of Acoustics,Speech,and Signal Processing,1987(35):162-712.
[3] SACCHI M D,ULRYCH T J.High-resolution velocity gathers and offset space reconstruction[J].Geophysics,1995,60(4):1169-1177.
[4] HERRMANN,P.,T.MOJESKY,M.MAGESAN et al.high-resolution Radon transforms:70th Annual International Meeting[C].SEG,2000:1953-1956.
[5] 吴律.Tau-p变换及应用[M].北京:石油工业出版社,1993.WU L.Tau-p transform and it’s application[M].Beijing:Petroleum Industry Press,1993.(In Chinese)
[6] 牛滨华,孙春岩.多项式Radon变换[J].地球物理学报,2001,44 (2):263-271.NIU B H,SUN C Y.Polynomial radon transform[J].Chinese Journal of Geophysics,44(2):263-271.(In Chinese)
[7] 王维红,首皓.线性同相轴波场分离的高分辨率Tau-p变换法[J].地球物理学进展,2006,21(1) :74-78.WANG W H,SHOU H.High resolution Tau-p transform using in linear events wavefield separation[J].Progress in Geophysics,2006,21(1):74-78.(In Chinese)
[8] 沈操,牛滨华,余钦范.Radon变换的Matlab实现[J].物探化探计算技术,2000,22 (4):346-351.SHEN C,NIU B H,YU Q F.Implementing Radon transform in Matlab[J].Techniques for Geophysical and Geochemical Exploration,2000,22(4):346-351.(In Chinese)
[9] 黄新武,孙春岩,牛滨华,等.基于数据一致性预测与压制自由表面多次波[J].地球物理学报,2005,48(1):173-180.HUANG X W,SUN C Y,NIU B H ,et al.Surface-related multiple prediction and suppression based on data-consistence:A theoretical study and test[J].Chinese Journal Of Geophysics,2005,48(1):173-180.(In Chinese)
[10] 刘喜武,刘洪.高分辨率Radon变换方法及其在地震信号处理中的应用[J].地球物理学进展,2004,19 (1):8-15.LIU X W,LIU H,et al.High resolution radon transform and its application in seismic signal processing.Progress[J] in Geophysics,2004,19(1):8-15.(In Chinese)
[11] 朱生旺,魏修成,李峰,等.用抛物线Radon变换稀疏解分离和压制多次波[J].石油地球物理勘探,2002,37 (2):110-115.ZHU S W,WEI X C,LI F,et al.Separating and suppressing multiples by sparse solution of parabolic Radon transform[J].Oil Geophysical Prospecting,2002,37(2):110-115.(In Chinese)
[12] 刘保童,朱光明.一种频率域提高Radon变换分辨率的方法[J].西安科技大学学报,2006,26(1):112-116.LIU B T,ZHU G M.An approach to improve the resolution of Radon transform in frequency domain[J].Journal of Xi’an University of Science and Technology,2006,26(1):112-116.(In Chinese)
[13] 巩向博.高分辨率Radon变换及其应用[D],长春:吉林大学,2008.GONG X B.High Resolution Radon transform and Application[D].Changchun :University,2008.(In Chinese)
[2] BEYLKIN G.The discrete Radon transform[J].IEEE Transactions of Acoustics,Speech,and Signal Processing,1987(35):162-712.
[3] SACCHI M D,ULRYCH T J.High-resolution velocity gathers and offset space reconstruction[J].Geophysics,1995,60(4):1169-1177.
[4] HERRMANN,P.,T.MOJESKY,M.MAGESAN et al.high-resolution Radon transforms:70th Annual International Meeting[C].SEG,2000:1953-1956.
[5] 吴律.Tau-p变换及应用[M].北京:石油工业出版社,1993.WU L.Tau-p transform and it’s application[M].Beijing:Petroleum Industry Press,1993.(In Chinese)
[6] 牛滨华,孙春岩.多项式Radon变换[J].地球物理学报,2001,44 (2):263-271.NIU B H,SUN C Y.Polynomial radon transform[J].Chinese Journal of Geophysics,44(2):263-271.(In Chinese)
[7] 王维红,首皓.线性同相轴波场分离的高分辨率Tau-p变换法[J].地球物理学进展,2006,21(1) :74-78.WANG W H,SHOU H.High resolution Tau-p transform using in linear events wavefield separation[J].Progress in Geophysics,2006,21(1):74-78.(In Chinese)
[8] 沈操,牛滨华,余钦范.Radon变换的Matlab实现[J].物探化探计算技术,2000,22 (4):346-351.SHEN C,NIU B H,YU Q F.Implementing Radon transform in Matlab[J].Techniques for Geophysical and Geochemical Exploration,2000,22(4):346-351.(In Chinese)
[9] 黄新武,孙春岩,牛滨华,等.基于数据一致性预测与压制自由表面多次波[J].地球物理学报,2005,48(1):173-180.HUANG X W,SUN C Y,NIU B H ,et al.Surface-related multiple prediction and suppression based on data-consistence:A theoretical study and test[J].Chinese Journal Of Geophysics,2005,48(1):173-180.(In Chinese)
[10] 刘喜武,刘洪.高分辨率Radon变换方法及其在地震信号处理中的应用[J].地球物理学进展,2004,19 (1):8-15.LIU X W,LIU H,et al.High resolution radon transform and its application in seismic signal processing.Progress[J] in Geophysics,2004,19(1):8-15.(In Chinese)
[11] 朱生旺,魏修成,李峰,等.用抛物线Radon变换稀疏解分离和压制多次波[J].石油地球物理勘探,2002,37 (2):110-115.ZHU S W,WEI X C,LI F,et al.Separating and suppressing multiples by sparse solution of parabolic Radon transform[J].Oil Geophysical Prospecting,2002,37(2):110-115.(In Chinese)
[12] 刘保童,朱光明.一种频率域提高Radon变换分辨率的方法[J].西安科技大学学报,2006,26(1):112-116.LIU B T,ZHU G M.An approach to improve the resolution of Radon transform in frequency domain[J].Journal of Xi’an University of Science and Technology,2006,26(1):112-116.(In Chinese)
[13] 巩向博.高分辨率Radon变换及其应用[D],长春:吉林大学,2008.GONG X B.High Resolution Radon transform and Application[D].Changchun :University,2008.(In Chinese)
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |