Simultaneous reconstruction and denoising of seismic data using multi-channel singular spectrum analysis based on hierarchical clustering
CAO Jingjie1,2, XU Changhao1,3, ZHU Yuefei1,2,4
1. Key Laboratory of Intelligent Detection and Equipment for Underground Space of Beijing-Tianjin-Hebei Urban Agglomeration, Ministry of Natural Resources, Shijiazhuang, Hebei 050031, China; 2. College of Mathematics and Physics, Hebei GEO University, Shijiazhuang, Hebei 050031, China; 3. Hebei Key Laboratory of Strategic Critical Mineral Resources, Shijiazhuang, Hebei 050031, China; 4. College of Geoscience and Surveying Engineering, China University of Mining and Technology-Beijing, Beijing 100083, China
Abstract:High signal-to-noise ratio, high fidelity, and high-resolution seismic data are prerequisites for clear imaging of subsurface structures. During seismic data acquisition, factors such as unfavorable topography and bad channels lead to the phenomenon that seismic data cannot satisfy the sampling theorem, and reconstruction is required to obtain complete seismic data. The multi-channel singular spectrum analysis method is a common seismic data reconstruction and denoising method, and its key is to determine the number of effective singular values for each divided data block. This parameter needs to be determined based on different data characteristics, but current manual selection methods require a lot of labor and computational resources. For each data block, based on the framework of multi-channel singular spectrum analysis in the frequency domain, the discrete point curves and spectral patterns of singular values are analyzed, and a hierarchical clustering method is proposed to automatically identify the number of singular values corresponding to the effective signals, which improves the denoising effect and reconstruction quality of seismic data. Then, the hierarchical clustering method is adopted to cluster the singular value sequence obtained by singular value decomposition of the block Hankel matrix and acquire the number of effective singular values of each frequency. Additionally, the number of singular values taken within the effective frequency is maximized to obtain the effective singular values for data blocks. In the framework of simultaneous reconstruction and denoising by damped multi-channel singular spectrum analysis, an improved multi-channel singular spectrum analysis (MSSA) method is put forward to realize simultaneous reconstruction and denoising of seismic data. Tests conducted on simulated and actual seismic data show that the MSSA method based on hierarchical clustering is superior in simultaneous reconstruction and denoising, and the accuracy of the obtained singular values is verified. The method can avoid the manual selection of the effective number of singular values and reduce the workload of seismic data processing, which is of practical significance for the reconstruction and denoising of large-scale seismic data.
付欣, 张小龙, 孙云强, 等. 基于双差反演策略的时移地震AVO反演[J]. 石油地球物理勘探, 2022, 57(4):888-896.FU Xin, ZHANG Xiaolong, SUN Yunqiang, et al. Time-lapse seismic AVO inversion based on a doubledifference inversion strategy[J]. Oil Geophysical Prospecting, 2022, 57(4):888-896.
[2]
陈生昌, 刘亚楠, 李代光. 基于波动方程偏移的地震数据地层成像[J]. 石油地球物理勘探, 2022, 57(5):1105-1113.CHEN Shengchang, LIU Yanan, LI Daiguang. Stratigraphic imaging with seismic data based on wave equation migration[J]. Oil Geophysical Prospecting, 2022, 57(5):1105-1113.
[3]
RONEN J. Wave-equation trace interpolation[J]. Geophysics, 1987, 52(7):973-984
[4]
CANNING A, GARDNER G. Regularizing 3-D data sets with DMO[J]. Geophysics, 1996, 61(4):1103-1114.
[5]
SACCHI M D, ULRYCH T J. High-resolution velocity gather and offset space reconstruction[J]. Geophysics, 1995, 60(4):1169-1177.
[6]
薛亚茹, 唐欢欢, 陈小宏. 高阶高分辨率Radon变换地震数据重建方法[J]. 石油地球物理勘探, 2014, 49(1):95-100, 131.XUE Yaru, TANG Huanhuan, CHEN Xiaohong. Seismic data reconstruction based on high order high resolution Radon transform[J]. Oil Geophysical Prospecting, 2014, 49(1):95-100, 131.
[7]
唐欢欢, 毛伟建, 詹毅. 3D高阶抛物Radon变换在不规则地震数据保幅重建中的应用[J]. 地球物理学报, 2020, 63(9):3452-3464.TANG Huanhuan, MAO Weijian, ZHAN Yi. Reconstruction of 3D irregular seismic data with amplitude preserved by high-order parabolic Radon transform[J]. Chinese Journal of Geophysics, 2020, 63(9):3452-3464.
[8]
HENNENFENT G, HERRMANN F J. Seismic denoising with nonuniformly sampled curvelets[J]. Computing in Science & Engineering, 2006, 8(3):16-25.
[9]
曹静杰, 王彦飞, 杨长春. 地震数据压缩重构的正则化与零范数稀疏最优化方法[J]. 地球物理学报, 2012, 55(2):596-607.CAO Jingjie, WANG Yanfei, YANG Changchun. Seismic data restoration based on compressive sensing using the regularization and zero-norm sparse optimization[J]. Chinese Journal of Geophysics, 2012, 55(2):596-607.
[10]
张华, 王冬年, 李红星, 等. 基于非均匀曲波变换的高精度地震数据重建[J]. 地球物理学报, 2017, 60(11):4480-4490.ZHANG Hua, WANG Dongnian, LI Hongxing, et al. High accurate seismic data reconstruction based on non-uniform curvelet transform[J]. Chinese Journal of Geophysics, 2017, 60(11):4480-4490.
[11]
刘洋, FOMEL S, 刘财, 等. 高阶seislet变换及其在随机噪声消除中的应用[J]. 地球物理学报, 2009, 52(8):2142-2151.LIU Yang, FOMEL S, LIU Cai, et al. High-order seislet transform and its application of random noise attenuation[J]. Chinese Journal of Geophysics, 2009, 52(8):2142-2151.
[12]
LIU Y, FOMEL S. OC-seislet:Seislet transform construction with differential offset continuation[J]. Geophysics, 2010, 75(6):WB235-WB245.
[13]
仝中飞. Curvelet阈值迭代法在地震数据去噪和插值中的应用研究[D]. 吉林长春:吉林大学, 2009.TONG Zhongfei. The Study on Seismic Data Denoising and Interpolation with Curvelet Thresholding Iterative Method[D]. Jilin University, Changchun, Jilin, 2009.
[14]
CAO J, ZHAO J. 3D seismic interpolation with a low redundancy, fast curvelet transform[J]. Journal of Seismic Exploration, 2015, 24(2):121-134.
[15]
葛子建. 基于傅里叶变换的POCS插值方法研究[D]. 四川成都:西南石油大学, 2014.GE Zijian. Research on POCS Interpolation Method Based on Fourier Transform[D]. Southwest Petroleum University, Chengdu, Sichuan, 2014.
[16]
曹静杰, 王本锋. 基于一种改进凸集投影方法的地震数据同时插值和去噪[J]. 地球物理学报, 2015, 58(8):2935-2947.CAO Jingjie, WANG Benfeng. An improved projection onto convex sets method for simultaneous interpolation and denoising[J]. Chinese Journal of Geophysics, 2015, 58(8):2935-2947.
[17]
OROPEZA V E. The Singular Spectrum Analysis Method and Its Application to Seismic Data Denoising and Reconstruction[D]. University of Alberta, Edmonton, 2010.
[18]
OROPEZA V, SACCHI M. Simultaneous seismic data denoising and reconstruction via multichannel singular spectrum analysis[J]. Geophysics, 2011, 76(3):V25-V32.
[19]
XUE Z G, HUANG J, LI Z C, et al. De-noising method for 3D field data based on the MSSA[C]. Extended Abstracts of 75th EAGE Conference & Exhibition, 2013, cp-348-00832.
[20]
CADZOW J A. Singal enhancement-a composite property mapping algorithm[J]. IEEE Transactions on Acoustics, Speech, and Signal Processing, 1988, 36(1):49-62.
[21]
TRICKETT S. F-xy Cadzow noise suppression[C]. SEG Technical Program Expanded Abstracts, 2008, 27:2586-2590..
[22]
GAO J, SACCHI M D, CHEN X. A fast rank reduction method for the reconstruction of 5D seismic volumes[C]. SEG Technical Program Expanded Abstracts, 2011, 30:3622-3626.
[23]
黄建平, 李闯, 李国磊, 等. 基于奇异谱分析的联合去噪及规则化方法[J]. 地球物理学进展, 2014, 29(4):1666-1671.HUANG Jianping, LI Chuang, LI Guolei, et al. Simultaneous seismic data de-noising and regularizationmethod based on singular spectrum analysis[J]. Progress in Geophysics, 2014, 29(4):1666-1671.
[24]
HUANG W, WANG R, CHEN Y, et al. Damped multichannel singular spectrum analysis for 3D random noise attenuation[J]. Geophysics, 2016, 81(4):V261-V270.
[25]
马继涛, 王建花, 刘国昌. 基于频率域奇异值分解的地震数据插值去噪方法研究[J]. 石油物探, 2016, 55(2):205-213.MA Jitao, WANG Jianhua, LIU Guochang. Seismic data noise attenuation and interpolation using singular value decomposition in frequency domain[J]. Geophysical Prospecting for Petroleum, 2016, 55(2):205-213.
[26]
CHEN Y, ZHANG D, JIN Z, et al. Simultaneous denoising and reconstruction of 5-D seismic data via damped rank-reduction method[J]. Geophysical Journal International, 2016, 206(3):1695-1717.
[27]
朱跃飞, 曹静杰, 殷晗钧. 一种自动判定保留的奇异值个数的地震随机噪声压制算法[J]. 石油地球物理勘探, 2022, 57(3):570-581.ZHU Yuefei, CAO Jingjie, YIN Hanjun. Seismic random noise suppression algorithm with automatic determination of the number of retained singular values[J]. Oil Geophysical Prospecting, 2022, 57(3):570-581.
[28]
BECKOUCHE S, MA J. Simultaneous dictionary learning and denoising for seismic data[J]. Geophysics, 2014, 79(3):A27-A31.
[29]
LIU L, MA J. Structured graph dictionary learning and application on the seismic denoising[J]. IEEE Transactions on Geoscience and Remote Sensing, 2019, 57(4):1883-1893..
[30]
LECUN Y, BOTTOU L, BENGIO Y, et al. Gradientbased learning applied to document recognition[J]. Proceedings of the IEEE, 1998, 86(11):2278-2324.
[31]
MASCI J, MEIER U, CIRESAN D, et al. Stacked convolutional auto-encoders for hierarchical feature extraction[C]. Artificial Neural Networks and Machine Learning-ICANN 2011, 2011, 52-59.
[32]
HE K, ZHANG X, REN S, et al. Deep residual learning for image recognition[C]. 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 2016, 770-778.
LIU G, REDA F A, SHIH J, et al. Image inpanting for irregular holes using partial convolution[C]. Computer Vision-ECCV 2018, 2018, 89-105.
[35]
SIAHKOOHI A, KUMAR R, HERRMANN F J. Deep-learning based ocean bottom seismic wavefield recovery[C]. SEG Technical Program Expanded Abstracts, 2019, 38:2232-2237.
[36]
王峰. 基于深度学习的地震数据去噪和重建方法的研究[D]. 浙江杭州:浙江大学, 2020.
[37]
WANG B, ZHANG N, LU W, et al. Deep-learning based seismic data interpolation:A preliminary result[J]. Geophysics, 2019, 84(1):V11-V20.
[38]
OLIVEIRA D A B, FERREIRA R S, SILVA R, et al. Interpolating seismic data with conditional generative adversarial networks[J]. IEEE Geoscience and Remote Sensing Letters, 2018, 15(12):1952-1956.
[39]
CHAI X, TANG G, WANG S, et al. Deep learning for irregularly and regularly missing 3-D data reconstruction[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020, 59(7):6244-6265.
[40]
FANG W, FU L, ZHANG M, et al. Seismic data interpolation based on U-net with texture loss U-net with texture loss for interpolation[J]. Geophysics, 2021, 86(1):V41-V54.