Abstract:Most of the projection onto convex sets (POCS) algorithms widely used for seismic signal reconstruction use linear or exponential threshold models,which have high computational efficiency,but have poor reconstruction effect due to the difficulty in eliminating the noise caused by the leakage of missing signals. Therefore,this paper proposes a POCS seismic signal reconstruction method based on the regional threshold model,which transforms the numerical threshold into a regional threshold,and the regional filtering window is iteratively updated as the threshold. The core idea is to reserve the transform coefficients of effective signal as much as possible by selecting a rectangular or a sectorial region of fixed size as a threshold according to a certain law in each POCS reconstruction iteration based on the frequency-wavenumber (F-K) spectrum distribution range of missing seismic signal in the spatiotemporal domain,and reserving and zeroing the transform coefficients inside and outside the region respectively. The rectangular and sectorial threshold models for POCS reconstruction of seismic signals are thus constructed. The numerical results demonstrated that compared with the POCS reconstruction of the exponential threshold model in the F-K domain,the regional threshold model in the F-K domain has a higher reconstruction accuracy for continuous missing signals. Compared with the sectorial region threshold model,the reconstruction accuracy and computational efficiency of the rectangular region threshold model are slightly higher. Compared with the exponential threshold model reconstructed by POCS in the curvelet domain,the reconstruction accuracy of the regional threshold model in the F-K domain is similar,but the computational efficiency is increased by about 90%.
GAO J J,CHEN X H,LI J Y,et al. Irregular seismic data reconstruction based on exponential threshold model of POCS method[J]. Applied Geophysics,2010,7(3):229-238.
[2]
WANG B F,WU R S,CHEN X H, et al. Simultaneous seismic data interpolation and denoising with a new adaptive method based on dreamlet transform[J]. Geophysical Journal International,2015,201(2):1182-1194.
[3]
SPITZ S. Seismic trace interpolation in the F-X domain[J]. Geophysics,1991,56(6):785-794.
[4]
周亚同,滕琳琳,李玲玲. 基于高阶扩展快速行进法的缺失地震数据重建[J]. 石油地球物理勘探,2015,50(5):873-880.ZHOU Yatong,TENG Linlin,LI Lingling. Seismic data reconstruction with the high-order expansion fast marching method[J]. Oil Geophysical Prospecting,2015,50(5):873-880.
[5]
RONEN J. Wave-equation trace interpolation[J]. Geophysics,1987,52(7):973-984.
GAO J J,SACCHI M D,CHEN X H. A fast reduced-rank interpolation method for prestack seismic volumes that depend on four spatial dimensions[J]. Geophysics,2013,78(1):V21-V30.
[8]
MA J W. Three-dimensional irregular seismic data reconstruction via low-rank matrix completion[J]. Geophysics,2013,78(5):V181-V192.
[9]
NAGHIZADEH M,SACCHI M D. On sampling functions and Fourier reconstruction methods[J]. Geophysics,2010,75(6):WB137-WB151.
[10]
张良,韩立国,许德鑫,等. 基于压缩感知技术的Shearlet变换重建地震数据[J]. 石油地球物理勘探,2017,52(2):220-225.ZHANG Liang,HAN Liguo,XU Dexin,et al. Seismic data reconstruction with Shearlet transform based on compressed sensing technology[J]. Oil Geophysical Prospecting,2017,52(2):220-225.
[11]
赵子越,李振春,张敏. 利用压缩感知技术的离散正交S变换地震数据重建[J]. 石油地球物理勘探,2020,55(1):29-35.ZHAO Ziyue,LI Zhenchun,ZHANG Min. Seismic data reconstruction using discrete orthonormal S-transform based on compressive sensing[J]. Oil Geophysical Prospecting,2020,55(1):29-35.
[12]
张华,王冬年,李红星,等. 基于非均匀曲波变换的高精度地震数据重建[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.
[13]
孙苗苗,李振春,李志娜,等. 基于压缩感知的加权MCA地震数据重构方法[J]. 地球物理学报,2019,62(3):1007-1021.SUN Miaomiao,LI Zhenchun,LI Zhina,et al. Reconstruction of seismic data with weighted MCA based on compressed sensing[J].Chinese Journal of Geophysics,2019,62(3):1007-1021.
[14]
董烈乾,张慕刚,汪长辉,等. 应用快速POCS算法的非均匀地震数据重建[J]. 石油地球物理勘探,2023,58(2):334-339.DONG Lieqian,ZHANG Mugang,WANG Changhui,et al. Reconstruction of non-uniformly sampled seismic data based on fast POCS algorithm[J]. Oil Geophysical Prospecting,2023,58(2):334-339.
[15]
石光明,刘丹华,高大化,等. 压缩感知理论及其研究进展[J]. 电子学报,2009,37(5):1070-1081.SHI Guangming,LIU Danhua,GAO Dahua,et al. Advances in theory and application of compressed sensing[J]. Acta Electronica Sinica,2009,37(5):1070-1081.
[16]
戴琼海,付长军,季向阳. 压缩感知研究[J]. 计算机学报,2011,34(3):3425-3434.DAI Qionghai,FU Changjun,JI Xiangyang. Research on compressed sensing[J]. Chinese Journal of Computers,2011,34(3):3425-3434.
[17]
GAN S W,WANG S D,CHEN Y K,et al. Compressive sensing for seismic data reconstruction via fast projection onto convex sets based on seislet transform[J]. Journal of Applied Geophysics,2016,130:194-208.
[18]
马泽川,李勇,陈力鑫,等. 地震资料快速两步插值算法[J]. 石油地球物理勘探,2020,55(5):997-1004.MA Zechuan,LI Yong,CHEN Lixin,et al. Fast two-step interpolation algorithm for seismic data[J]. Oil Geophysical Prospecting,2020,55(5):997-1004.
[19]
张华,陈小宏. 基于jitter采样和曲波变换的三维地震数据重建[J]. 地球物理学报,2013,56(5):1637-1649.ZHANG Hua,CHEN Xiaohong. Seismic data reconstruction based on jittered sampling and curvelet transform[J]. Chinese Journal of Geophysics,2013,56(5):1637-1649.
[20]
GAO J J,STANTON A,NAGHIZADEH M,et al. Convergence improvement and noise attenuation considerations for beyond alias projection onto convex sets reconstruction[J]. Geophysical Prospecting, 2013, 61(S1):138-151.
[21]
刁塑,张华,张恒琪,等. 利用抗假频凸集投影算法的规则缺失地震数据重建[J]. 石油地球物理勘探,2020,55(4):716-724.DIAO Su,ZHANG Hua,ZHANG Hengqi,et al. Reconstructing regularly missing seismic data based on anti-aliasing POCS algorithm[J]. Oil Geophysical Prospecting,2020,55(4):716-724.
[22]
闫浩飞, 刘国峰, 薛典军,等. 基于凸集投影方法的重磁数据规则缺失重建[J]. 地球物理学进展,2016,31(5):2192-2197.YAN Haofei,LIU Guofeng,XUE Dianjun,et al. Reconstruction of gravity/magnetic data with the projection-onto-convex-sets methods[J]. Progress in Geophysics,2016,31(5):2192-2197.
[23]
曾小牛,李夕海,刘继昊,等. 基于凸集投影的重力数据扩充下延一体化方法[J]. 石油地球物理勘探,2019,54(5):1166-1173.ZENG Xiaoniu,LI Xihai,LIU Jihao,et al. An integration of interpolation,edge padding and downward continuation for gravity data based on the projection onto convex sets[J]. Oil Geophysical Prospecting, 2019, 54(5):1166-1173.
[24]
ZHANG H,CHEN X H,WU X M. Seismic data reconstruction based on CS and Fourier theory[J]. Applied Geophysics,2013,10(2):170-180.
[25]
WANG Y F,CAO J J,YANG C C. Recovery of seismic wavefields based on compressive sensing by an l1norm constrained trust region method and the piecewise random subsampling[J]. Geophysical Journal International,2011,187(1):199-213.
[26]
SHAHIDI R,TANG G,MA J W,et al. Application of randomized sampling schemes to curvelet-based sparsity-promoting seismic data recovery[J]. Geophysical Prospecting,2013,61(5):973-997.
[27]
王本锋,陆文凯,陈小宏,等. 基于3D Curvelet变换的频率域高效地震数据插值方法研究[J]. 石油物探,2018,57(1):65-71.WANG Benfeng,LU Wenkai,CHEN Xiaohong,et al. Efficient seismic data interpolation using three-dimensional Curvelet transform in the frequency domain[J]. Geophysical Prospecting for Petroleum, 2018, 57(1):65-71.
[28]
王本锋,陈小宏,李景叶,等. POCS联合改进的Jitter采样理论曲波域地震数据重建[J]. 石油地球物理勘探,2015,50(1):20-28.WANG Benfeng,CHEN Xiaohong,LI Jingye,et al. Seismic data reconstruction based on POCS and improved Jittered sampling in the curvelet domain[J]. Oil Geophysical Prospecting,2015,50(1):20-28.
[29]
ABMA R,KABIR N. 3D interpolation of irregular data with a POCS algorithm[J]. Geophysics,2006,71(6):E91-E97.
[30]
黄小刚,王一博,刘伊克,等. 半径-斜率域加权反假频地震数据重建[J]. 地球物理学报, 2014,57(7):2278-2290.HUANG Xiaogang,WANG Yibo,LIU Yike,et al. Weighted anti-alias seismic data reconstruction in R-P domain[J]. Chinese Journal of Geophysics, 2014, 57(7):2278-2290.
[31]
YANG P L,GAO J H,CHEN W C. On analysis-based two-step interpolation methods for randomly sampled seismic data[J]. Computers & Geosciences,2013,51:449-461.
[32]
刘国昌,陈小宏,郭志峰,等. 基于Curvelet变换的缺失地震数据插值方法[J]. 石油地球物理勘探,2011,46(2):237-246.LIU Guochang,CHEN Xiaohong,GUO Zhifeng,et al. Missing seismic data rebuilding by interpolation based on Curvelet transform[J]. Oil Geophysical Prospecting,2011,46(2):237-246.
[33]
CANDES E J,ROMBERG J K. Signal recovery from random projections[C]. Proceedings Volume 5674,Computational Imaging Ⅲ,SPIE,2005,76-86.
[34]
温睿,刘国昌,冉扬. 压缩感知地震数据重建中的三个关键因素分析[J]. 石油地球物理勘探,2018,53(4):682-693.WEN Rui,LIU Guochang,RAN Yang. Three key factors in seismic data reconstruction based on compressive sensing[J]. Oil Geophysical Prospecting,2018,53(4):682-693.
[35]
徐卫,张华,张落毅. 基于复值曲波变换的地震数据重建方法[J]. 物探与化探,2016,40(4):750-756.XU Wei,ZHANG Hua,ZHANG Luoyi. A study of seismic data reconstruction based on complex-valued curvelet transform[J]. Geophysical and Geochemical Exploration,2016,40(4):750-756.
[36]
周亚同,刘志峰,张志伟. 形态分量分析框架下基于DCT与曲波字典组合的地震信号重建[J]. 石油物探,2015,54(5):560-568,581.ZHOU Yatong,LIU Zhifeng,ZHANG Zhiwei. Seismic signal reconstruction under the morphological component analysis framework combined with discrete cosine transform (DCT) and curvelet dictionary[J]. Geophysical Prospecting for Petroleum,2015,54(5):560-568,581.
[37]
张岩,崔淋淇,宋利伟,等. 基于空间-波数域联合深度学习的数值频散压制[J]. 石油地球物理勘探,2022,57(3):510-524.ZHANG Yan,CUI Linqi,SONG Liwei,et al. Numerical dispersion suppression based on joint deep learning in the space and wave number domains[J]. Oil Geophysical Prospecting,2022,57(3):510-524.