非局部平均滤波噪声压制方法及其在VSP资料逆时偏移中的应用
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
基于双程波的逆时偏移会产生低频成像噪声,在成像后运用拉普拉斯滤波法可以取得较好的压噪效果,但是,该方法严重依赖于角度参数,使得滤波后的成像剖面上常常存在同相轴不光滑和噪声压制不完全的情况。基于叠后地震资料同相轴结构的特点,为进一步提高逆时偏移叠加数据的信噪比,引入了非局部平均滤波法,针对拉普拉斯滤波后的成像数据进行进一步的压噪处理。非局部平均滤波法将输入的地震数据分解为不含噪声的地震数据和噪声数据两部分,利用不同成像点与其它成像点间的相似系数,实现滤波处理。二维复杂模型VSP正演模拟数据逆时偏移结果的试算表明,应用非局部平均滤波后的剖面信噪比得到进一步提高,同相轴的连续性明显增强。实际VSP资料逆时偏移低频噪声压制试处理也表明非局部平均滤波方法具有计算精度高、算法稳定性好和易于实现的特点。
Low-frequency noise is the inherent property of reverse-time migration(RTM)based on two-way wave equation.Filtering after imaging is widely used to suppress this kind of noise,especially Laplacian filtering method.But the result of Laplacian method depends on angle parameters,the case of rough events and residual noises always exist.In order to improve signal to noise ratio of imaging data of reverse-time migration,Non-local Means(NLM)filtering approach is introduced in this paper,which is used to suppress noise based on the characteristics of event structure of post-stack seismic data,and is applied to the seismic data after Laplacian filtering.By using NLM method,the input seismic data can be divided into noise and data,and the similarity coefficient of different imaging points is applied to realize filtering.For 2Dcomplicated theoretical model,we carried out forward modeling and RTM on VSP data.The results demonstrate that the combination of Laplacian and NLM filtering is better than the result of Laplacian filtering.Moreover,the application of actual VSP seismic data also suggests NLM filtering has advantages in computational accuracy,stability and practicability.
引文
[1]Zhang Y,Duan L,Xie Y.A stable and practical implementation of least-squares reverse time migration[J].Geophysics,2015,80(1):23-31
    [2]刘文卿,王西文,刘洪,等.盐下构造速度建模与逆时偏移成像研究及应用[J].地球物理学报,2013,56(2):616-625Liu W Q,Wang X W,Liu H,et al.Application of velocity modeling and reverse time migration to subsalt structure[J].Chinese Journal of Geophysics(in Chinese),2013,56(2):612-625
    [3]孙赞东.三维三分量VSP方法原理及应用[M].北京:石油工业出版社,2011:1-4Sun Z D.The principle and application of 3Dthreecomponent VSP method[M].Beijing:Petroleum industry press,2011:1-4
    [4]朱金明,严俊华.VSP逆时偏移处理[J].石油地球物理勘探,1991,26(5):564-570Zhu J M,Yan J H.VSP reverse-time migration[J].Oil Geophysical Prospecting,1991,26(5):564-570
    [5]朱金明,董敏煜,李承楚.VSP的双程无反射波动方程逆时偏移[J].石油地球物理勘探,1989,24(3):256-270Zhu J M,Dong M Y,Li C C.VSP reverse-time migration using two-way nonreflection wave equation[J].Oil Geophysical Prospecting,1989,24(3):256-270
    [6]Claerbout J F.Toward a unified theory of reflector mapping[J].Geophysics,1971,36(3):467-481
    [7]刘红伟,刘洪,邹振,等.地震叠前逆时偏移中的去噪与存储[J].地球物理学报,2010,53(9):2171-2180Liu H W,Liu H,Zou Z,et al.The problems of denoise and storage in seismic reverse time migration[J].Chinese Journal of Geophysics(in Chinese),2010,53(9):2171-2180
    [8]Baysal E,Kosloff D,Sherwood J.A two-way nonreflecting wave equation[J].Geophysics,1984,49(2):132-141
    [9]Loewenthal D,Stoffa P A,Faria E L.Suppressing the unwanted reflections of the full wave equation[J].Geophysics,1987,52(7):1007-1012
    [10]Fletcher R F,Fowler,Kitchenside P.Suppressing artifacts in prestack reverse-time migration[J].Expanded Abstracts of 75th Annual Internat SEG Mtg,2005,2049-2051
    [11]Kaelin B,Guitton A.Imaging condition for reverse time migration[J].Expanded Abstracts of 76th Annual Internat SEG Mtg,2006,2594-2598
    [12]Mulder W,Plessix R.One-way and two-way waveequation migration[J].Expanded Abstracts of 73rd Annual Internat SEG Mtg,2003,881-884
    [13]Antoine Guitton,Bruno Kaelin,Biondo Biond.Least-square attenuation of reverse time migration artifacts[J].Expanded Abstracts of 76th Annual Internat SEG Mtg,2006,2348-2352
    [14]Zhang Y,Sun J.Practical issues of reverse time migration:true amplitude gathers,noise removal and harmonic-source encoding[J].First Break,2009,26(1):29-35
    [15]Buades A,Coll B,Morel J M.A review of image denoising algorithms,with a new one[J].Multiscale Modeling and Simulation,2005(4):490-530
    [16]Buades A,Coll B,Morel J M.Image denoising methods:a new nonlocal principle[J].SIAM Review,2010,52(1):113-147
    [17]Sheng B,Li P,Sun H.Image-based material restyling with fast non-local means filtering[J].Proceedings of the 2009 Fifth International Conference on Image and Graphics,2009,841-846
    [18]Mahmoudi M,Sapiro G.Fast image and video denoising via nonlocal means of similar neighborhoods[J].IEEE Signal Processing Letters,2005,12,839-842
    [19]Brox T,Kleinschmidt O,Cremers D.Efficient nonlocal means for denoising of textural patterns[J].IEEE Transactions on Image Processing,2008,17(7):1083-1092
    [20]CoupéP,Yger P,Prima S,et al.An optimized blockwise nonlocal means denoising filter for 3-D magnetic resonance images[J].IEEE Transactions on Medical Imaging,2008,27(4):425-441
    [21]Deledalle C A,Denis L,Tupin F.NL-InSAR:nonlocal interferogram estimation[J].IEEE Transactions on Geoscience and Remote Sensing,2011,49(4):1441-1452
    [22]Zoican S.Speech de-noising system with nonlocal means algorithm[J].9thInternational Symposium on Electronics and Telecommunications(ISETC),2010,315-318
    [23]Wei D Y,Yin C C.An optimized locally adaptive non-local means denoising filter for cryo-electron microscopy data[J].Journal of Structural Biology,2010,172(1):211-218
    [24]Bonar D,Sacchi M.Denoising seismic data using the nonlocal means algorithm[J].Geophysics,2012,77(1):A5-A8

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心