1. State Key Laboratory of Nuclear Resources and Environment, East China University of Technology, Nanchang, Jiangxi 330013, China; 2. School of Geophysics and Measurement-Control Technology, East China University of Technology, Nanchang, Jiangxi 330013, China
Abstract:The conventional total variation (TV) regularization model only considers the first-order derivative information in the horizontal and vertical directions. When dealing with prestack seismic data with curved reflection events,it can severely damage the amplitude information and cause “staircase effects” by suppressing the lateral gradient characteristics of the amplitude. The local dip information of seismic data is often applied to improve the amplitude-preserving ability of the TV model. However,the calculation of local dip information itself will be impacted by noise. To address this issue,this paper proposes a high-order TV regularization model to suppress random noise in prestack seismic data in the domain of normal moveout(NMO). This method first transforms the prestack seismic data into the NMO domain,NMO is robust to noise and avoids the calculation of the local dip angle. In the NMO domain,the curved event is flattened,and then high-order TV denoising is performed. Finally,the prestack seismic data are restored through inverse NMO. Taking the second-order derivative as an example,a high-order TV regularization inversion denoising objective function is constructed,and a fast optimization method is derived under the split Bregman optimization framework. The processing results of synthetic seismic data and actual seismic data show that this method can not only effectively suppress random noise but also eliminate amplitude distortion caused by curved reflection events and “staircase effects”, improving the amplitude preservation performance of the TV denoising method.
基金资助:本项研究受国家自然科学基金项目“组稀疏结构和等效衰减模型双重约束下的叠前Q值反演方法研究”(42004114)、江西省自然科学基金项目“基于压缩感知的地震数据自适应压缩及反射系数快速反演”(20202BAB211010)、“基于人工智能的江西地区天然地震和非天然地震事件识别方法研究”(20224BAB213047)、江西省教育厅科学技术研究项目“定向高阶导数 TV 正则化保幅地震噪声压制算法研究”(GJJ2200746)、核资源与环境国家重点实验室开放基金项目“致密层系井震结合计算三维TOC实现油铀兼探方法研究”(2020NRE27)、“井震结合下基于谱蓝化—有色反演的松辽盆地南部姚家组砂岩型铀矿预测方法研究”(2022NRE16)和东华理工大学研究生创新专项资金项目“基于同相轴拉平技术的高阶TV正则化地震资料保幅去噪算法研究”(DHYC-202314)联合资助。
WU N,LI Y,YANG B. Noise attenuation for 2-D seismic data by radial-trace time-frequency peak filtering[J]. IEEE Geoscience and Remote Sensing Letters,2011,8(5):874-878.
[4]
董烈乾,汪长辉,李长芬,等. 利用自适应中值滤波方法压制混叠噪声[J]. 地球物理学进展,2018,33(4):1475-1479.DONG Lieqian, WANG Changhui, LI Changfen,et al. Blending noise removal utilizing an adaptive median filter[J]. Progress in Geophysics,2018,33(4):1475-1479.
[5]
国胧予,刘财,刘洋.滤波类方法衰减地震数据噪声[J]. 地球物理学进展,2018,33(5):1890-1896.GUO Longyu,LIU Cai,LIU Yang. Filtering methods attenuate seismic data noise[J]. Progress in Geophysics,2018,33(5):1890-1896.
[6]
TARY J B,HERRERA R H,HAN J,et al. Spectral estimation:What is new? What is next?[J]. Reviews of Geophysics,2014,52(4):723-749.
[7]
周怀来. 基于小波变换的地震信号去噪方法研究与应用[D]. 四川成都:成都理工大学,2006.ZHOU Huailai. Research and Application of Seismic Signal Denoising Method Based on Wavelet Transform[D]. Chengdu University of Technology,Chengdu,Sichuan,2006.
[8]
蔡剑华,王先春,胡惟文. 基于经验模态分解与小波阈值的MT信号去噪方法[J]. 石油地球物理勘探,2013,48(2):303-307.CAI Jianhua,WANG Xianchun,HU Weiwen. A method for MT data denoising based on empirical mode decomposition and wavelet threshold[J]. Oil Geophysical Prospecting,2013,48(2):303-307.
[9]
CANDÈS E,DONOHO D. Continuous curvelet transform Ⅱ: Discretization and frames[J]. Applied and Computational Harmonic Analysis,2005,19(2):198-222.
韩卫雪,周亚同,池越. 基于深度学习卷积神经网络的地震数据随机噪声去除[J]. 石油物探,2018,57(6):862-869,877.HAN Weixue,ZHOU Yatong,CHI Yue. Deep learning convolutional neural networks for random noise attenuation in seismic data[J]. Geophysical Prospecting for Petroleum,2018,57(6):862-869,877.
[12]
张岩,李新月,王斌,等. 基于深度学习的鲁棒地震数据去噪[J]. 石油地球物理勘探,2022,57(1):12-25.ZHANG Yan,LI Xinyue,WANG Bin,et al. Robust seismic data denoising based on deep learning[J]. Oil Geophysical Prospecting,2022,57(1):12-25.
[13]
李佳,王维波,盛立,等. 应用双向长短时记忆神经网络的微地震信号降噪方法[J]. 石油地球物理勘探,2023,58(2):285-294.LI Jia,WANG Weibo,SHENG Li,et al. Denoising of microseismic signal based on bidirectional long short-term memory neural network[J]. Oil Geophysical Prospecting,2023,58(2):285-294.
[14]
时磊,刘俊州,韦婉婉,等. 基于贝叶斯反演的叠前数据同时插值和保幅去噪方法研究[J]. 地球物理学进展,2019,34(6):2309-2314.SHI Lei,LIU Junzhou,WEI Wanwan,et al. Research on simultaneous interpolation and amplitude-preserving denoising based on Bayesian inversion of prestack data[J]. Progress in Geophysics,2019,34(6):2309-2314.
[15]
任浩,李宗杰,薛姣,等. 基于稀疏反演的多道匹配追踪地震信号去噪方法及其应用[J]. 石油物探,2019,58(2):199-207.REN Hao,LI Zongjie,XUE Jiao,et al. Multichannel matching pursuit based on sparse inversion for seismic data denoising and its application[J]. Geophysical Prospecting for Petroleum,2019,58(2):199-207.
[16]
RUDIN L I,OSHER S,FATEMI E. Nonlinear total variation based noise removal algorithms[J]. Physica D: Nonlinear Phenomena,1992,60(1):259-268.
[17]
OSHER S, BURGER M, GOLDFARB D,et al. An iterative regularization method for total variation-based image restoration[J]. Multiscale Modeling & Simulation,2005,4(2):460-489.
[18]
陈勇,韩波,肖龙,等. 多尺度全变分法及其在时移地震中的应用[J].地球物理学报,2010,53(8):1883-1892.CHEN Yong,HAN Bo,XIAO Long, et al. Multiscale total variation method and its application on time-lapse seismic[J]. Chinese Journal of Geophysics,2010,53(8):1883-1892.
[19]
CHAN T F,ESEDOGLU S,PARK F E. Image decomposition combining staircase reduction and texture extraction[J]. Journal of Visual Communication and Image Representation,2007,18(6):464-486.
[20]
LEFKIMMIATIS S,BOURQUARD A,UNSER M. Hessian-based norm regularization for image restoration with biomedical applications[J]. IEEE Transactions on Image Processing,2012,21(3):983-995.
[21]
Kong D H,Peng Z M. Seismic random noise attenuation using shearlet and total generalized variation[J]. Journal of Geophysics and Engineering,2015,12(6):1024-1035.
[22]
BAYRAM İ,KAMASAK M E. Directional total variation[J]. IEEE Signal Processing Letters,2012,19(12):781-784.
[23]
ZHANG H,WANG Y Q. Edge adaptive directional total variation[J]. The Journal of Engineering,2013,2013(11):61-62.
[24]
WANG D H,GAO J H,LIU N H,et al. Structure-oriented DTGV regularization for random noise attenuation in seismic data[J].IEEE Transactions on Geoscience and Remote Sensing,2021,59(2):1757-1771.
[25]
LYSAKER M,LUNDERVOLD A,TAI X C. Noise removal using fourth-order partial differential equation with applications to medical magnetic resonance images in space and time[J]. IEEE Transactions on Image Processing,2003,12(12):1579-1590.
[26]
STEIDL G. A note on the dual treatment of higher-order regularization functionals[J]. Computing,2006,76(1):135-148.
[27]
LYSAKER M,TAI X. Iterative image restoration combining total variation minimization and a second-order functional[J]. International Journal of Computer Vision,2006,66(1):5-18.
[28]
LIU X Y,LI Q,YUAN C,et al. High-order directional total variation for seismic noise attenuation[J]. IEEE Transactions on Geoscience and Remote Sensing,2022,60:1-13.
[29]
FOMEL S. Applications of plane-wave destruction filters[J]. Geophysics,2002,67(6):1946-1960.
[30]
CHEN Y K,HUANG W L,ZHOU Y T,et al. Plane-wave orthogonal polynomial transform for amplitude-preserving noise attenuation[J]. Geophysical Journal International,2018,214(3):2207-2223.
[31]
郝亚炬,张鹏,文晓涛,等. 横向二阶导数TV正则化三维反射系数反演[J].石油地球物理勘探, 2023,58(3):680-689.HAO Yaju,ZHANG Peng,WEN Xiaotao,et al. 3D reflectivity inversion method based on TV regularization of lateral second-order derivatives[J]. Oil Geophysical Prospecting,2023,58(3):680-689.
[32]
GOLDSTEIN T,OSHER S. The split Bregman method for L1-regularized problems[J]. SIAM Journal on Imaging Sciences,2009,2(2):323-343.
[33]
CHEN Y P,PENG Z M,LI M H,et al. Seismic signal denoising using total generalized variation with overlapping group sparsity in the accelerated ADMM framework[J]. Journal of Geophysics and Engineering,2019,16(1):30-51.