摘要
提出了一种基于W加权核范数最小化的地震信号盲去噪算法。采用主成分分析法估计地震信号噪声水平,借助加权核范数最小化(WNNM)实现去噪。在去噪中通过权值分配控制矩阵奇异值的收缩程度,提升了算法性能。分别对三种地震信号进行去噪,并与双树复小波变换、曲波变换、WNNM算法进行了性能对比。研究结果表明,该算法在噪声水平未知的情况下,能有效去除地震信号所含噪声,去噪效果优于传统去噪算法。
A seismic signal blind denoising algorithm is proposed based on W-weighted nuclear norm minimization.The noise level of seismic signals is estimated by principal component analysis and the denoising is realized by weighted nuclear norm minimization(WNNM).In denoising,the shrinkage degree of singular values of a matrix is controlled by weight assignment,and the performances of the algorithm is improved.Three kinds of seismic signals are denoised,respectively.The performance is compared with double tree complex wavelet transform,curvelet transform and the WNNM algorithm.The research results show that the proposed algorithm can effectively remove the noises contained in seismic signals when the noise level is unknown.Moreover,the denoising effect is superior to those of the traditional denoising algorithms.
引文
[1]Zhang H,Chi Y,Zhou Y T,et al.Three dimensional seismic signal denoising based on fourdimensional block matching cooperative filtering combined with principle component analysis[J].Laser&Optoelectronics Progress,2018,55(4):041007.张欢,池越,周亚同,等.结合主成分分析的四维块匹配协同滤波三维地震信号去噪[J].激光与光电子学进展,2018,55(4):041007.
[2]Zhang M,LüX Q,Wu L,et al.Multiplicative denoising method based on deep residual learning[J].Laser&Optoelectronics Progress,2018,55(3):031004.张明,吕晓琪,吴凉,等.基于深度残差学习的乘性噪声去噪方法[J].激光与光电子学进展,2018,55(3):031004.
[3]Zheng C L,He D D,Fei Q G.Improved digital image correlation method based on gray gradient denoised by regularization method[J].Acta Optica Sinica,2018,38(8):0812002.郑成林,何顶顶,费庆国.基于灰度梯度正则化去噪的改进数字图像相关法[J].光学学报,2018,38(8):0812002.
[4]Zhang J H,LüN,Tian L Y,et al.An overview of the methods and techniques for seismic data noise attenuation[J].Progress in Geophysics,2006,21(2):546-553.张军华,吕宁,田连玉,等.地震资料去噪方法技术综合评述[J].地球物理学进展,2006,21(2):546-553.
[5]Shao J,Sun C Y,Tang J,et al.Micro-seismic data denoising based on sparse representations over learned dictionary in the wavelet domain[J].Oil Geophysical Prospecting,2016,51(2):254-260.邵婕,孙成禹,唐杰,等.基于字典训练的小波域稀疏表示微地震去噪方法[J].石油地球物理勘探,2016,51(2):254-260.
[6]Hennenfent G,Herrmann F J.Seismic denoising with nonuniformly sampled curvelets[J].Computing in Science&Engineering,2006,8(3):16-25.
[7]Li J H,Zhang Y J,Qi R,et al.Wavelet-based higher order correlative stacking for seismic data denoising in the curvelet domain[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2017,10(8):3810-3820.
[8]Liu C M,Wang D L,Wang T,et al.Random seismic noise attenuation based on the Shearlet transform[J].Acta Petrolei Sinica,2014,35(4):692-699.刘成明,王德利,王通,等.基于Shearlet变换的地震随机噪声压制[J].石油学报,2014,35(4):692-699.
[9]Chen Y K,Zhou Y T,Chen W,et al.Empirical lowrank approximation for seismic noise attenuation[J].IEEE Transactions on Geoscience and Remote Sensing,2017,55(8):4696-4711.
[10]Zhou Y T,Shi C J,Chen H M,et al.Spike-like blending noise attenuation using structural low-rank decomposition[J].IEEE Geoscience and Remote Sensing Letters,2017,14(9):1633-1637.
[11]Gu S H,Zhang L,Zuo W M,et al.Weighted nuclear norm minimization with application to image denoising[C]//IEEE Conference on Computer Vision and Pattern Recognition,2014:2862-2869.
[12]Gu S H,Xie Q,Meng D Y,et al.Weighted nuclear norm minimization and its applications to low level vision[J].International Journal of Computer Vision,2017,121(2):183-208.
[13]Wang D X.The research of seismic random noise suppression method based on the WNNM[D].Changchun:Jilin University,2017:63-86.王代香.基于WNNM的地震勘探随机噪声压制方法的研究[D].长春:吉林大学,2017:63-86.
[14]Liu X H,Tanaka M,Okutomi M.Noise level estimation using weak textured patches of a single noisy image[C]//IEEE International Conference on Image Processing,2012:665-668.
[15]Pyatykh S,Hesser J,Zheng L.Image noise level estimation by principal component analysis[J].IEEETransactions on Image Processing,2013,22(2):687-699.
[16]Ke Q F,Kanade T.Robust L1norm factorization in the presence of outliers and missing data by alternative convex programming[C]//IEEE Computer Society Conference on Computer Vision and Pattern Recognition,2005,1:739-746.
[17]Wright J,Ganesh A,Rao S,et al.Robust principal component analysis:exact recovery of corrupted lowrank matrices via convex optimization[C]//Advances in Neural Information Processing Systems,2009,87(4):20-56.
[18]Fazel M,Hindi H,Boyd S P.A rank minimization heuristic with application to minimum order system approximation[C]//American Control Conference,2001:4734-4739.
[19]Blomberg N,Rojas C R,Wahlberg B.Approximate regularization paths for nuclear norm minimization using singular value bounds[C]//IEEE Signal Processing and Signal Processing Education Workshop,2015:190-195.
[20]Au-Yeung E.Dictionary learning for sparse representation based on nuclear norm minimization[C]//International Conference on Sampling Theory and Applications,2017:570-574.
[21]Candès E J,Recht B.Exact matrix completion via convex optimization[J].Foundations of Computational Mathematics,2009,9(6):717-772.
[22]Cai J F,Candès E J,Shen Z W.A singular value thresholding algorithm for matrix completion[J].SIAM Journal on Optimization,2010,20(4):1956-1982.
[23]Dabov K,Foi A,Katkovnik V,et al.Image denoising by sparse 3D transform-domain collaborative filtering[J].IEEE Transactions on Image Processing,2007,16(8):2080-2095.
[24]Chan R H,Ng M K.Conjugate gradientmethods for Toeplitz systems[J].SIAM Review,1996,38(3):427-482.
[25]Bhaumik D K,Kapur K,Gibbons R D.Testing parameters of a gamma distribution for small samples[J].Technometrics,2009,51(3):326-334.
[26]Witkovsk′y V.Computing the distribution of a linear combination of inverted gamma variables[J].IEEETransactions on Geoscience&Remote Sensing,2001,37(1):79-90.
