张春枫,杜启振,张富源,符力耘,魏国华
ZHANG Chunfeng,DU Qizhen,ZHANG Fuyuan,FU Liyun,WEI Guohua

• 1:中国石油大学(华东)山东省深层油气重点实验室
• 2:青岛海洋科学与技术试点国家实验室海洋矿产资源评价与探测技术功能实验室
• 3:中国石化胜利油田分公司物探研究院

摘要(Abstract)：

由于地表条件的限制，野外采集的地震数据难以做到规则采样，不规则采样地震数据影响了地震数据处理与解释。针对这一问题，采用匹配追踪傅里叶插值算法对不规则地震数据进行重构，并在该插值算法中的迭代部分作出改进，使其只在频率-波数域中计算傅里叶系数谱，减少了迭代中大量的正反傅里叶变换，以降低匹配追踪傅里叶插值算法的计算时间，并将改进后的方法应用到随机欠采样凹陷正演地震记录、Marmousi正演地震记录以及实际CMP地震资料中进行规则化重构。结果表明：改进后的匹配追踪傅里叶插值算法可以将随机欠采样地震数据规则化，重构后其同相轴变得连续，频波谱能量收敛，有效信号得以恢复，且改进后该算法的重构时间明显降低，其加速比在30倍以上，表明了该方法的正确性与可行性。改进后的匹配追踪算法为不规则地震数据的重构提供了新的理论依据。
Due to the limitation of surface conditions,it is difficult to conduct regular sampling for seismic data collected in the field,and irregular sampling of seismic data affects seismic data processing and interpretation. To solve this problem,matching pursuit Fourier interpolation algorithm is used to reconstruct irregular seismic data, and iterative part of interpolation algorithm is improved,so that it only calculates Fourier coefficient spectrum in frequency wavenumber domain,reducing a large amount of forward and inverse Fourier transforms in iteration,so as to reduce calculation time of matching pursuit Fourier interpolation algorithm. The improved method is applied to forward seismic records of random undersampling graben,Marmousi forward seismic records and actual CMP seismic data for regular reconstruction. The results show that improved matching pursuit Fourier interpolation algorithm can regularize random undersampling seismic data. After reconstruction,its events become continuous,frequency spectrum energy converges,and effective signal are recovered. Moreover,the reconstruction time of improved algorithm is significantly reduced,and its acceleration ratio is more than 30 times,indicating the correctness and feasibility of the method. Improved matching pursuit algorithm provides new theoretical basis for reconstruction of irregular seismic data.

关键词(KeyWords)： 地震数据规则化;地震数据重构;地震数据插值;匹配追踪;傅里叶变换
regularization of seismic data;seismic data reconstruction;seismic data interpolation;matching pursuit;Fourier transform

Abstract:

Keywords:

基金项目(Foundation): 国家自然科学基金重点项目“弹性全波形反演与逆时偏移成像”（41930429）

作者(Author): 张春枫,杜启振,张富源,符力耘,魏国华
ZHANG Chunfeng,DU Qizhen,ZHANG Fuyuan,FU Liyun,WEI Guohua

DOI: 10.19597/j.issn.1000-3754.202210057

参考文献(References)：

• [1]霍志周，熊登，张剑锋.地震数据重建方法综述[J].地球物理学进展，2013,28(4):1749-1756.HUO Zhizhou,XIONG Deng,ZHANG Jianfeng.The overview of seismic data reconstruction methods[J].Progress in Geophysics,2013,28(4):1749-1756.
• [2] SPITZ S. Seismic trace interpolation in the F-X domain[J].Geophysics,1991,56(6):785-794.
• [3] GULUNAY N, CHAMBERS R E. Unaliased f-k domain trace interpolation(UFKI)[C]. Denver:SEG Annual Meeting,1996.
• [4]梁晨曦，张固澜，李磊，等.能量归一块匹配三维协同滤波及地震去噪应用[J].石油物探，2022,61(3):499-511.LIANG Chenxi, ZHANG Gulan, LI Lei, et al. An energy-normalized block-matching three-dimensional collaborative filtering and its application in seismic denoising[J]. Geophysical Prospecting for Petroleum,2022.
• [5] PORSANI M J.Seismic trace interpolation using half-step prediction filters[J].Geophysics,1999,64(5):1461-1467.
• [6] RONEN J. Wave-equation trace interpolation[J]. Geophysics,1987,52(7):973-984.
• [7] TRICKEET S,BURROUGHS L,MILTON A,et al.Rank-reduction-based trace interpolation[C].Denver:SEG Annual Meeting,2010.
• [8]李言言，李志明，王琴.基于非凸张量秩最小化的三维地震数据插值[J].工程地球物理学报，2022,19(3):410-418.LI Yanyan, LI Zhiming, WANG Qin. A nonconvex tensor rank minimization-based approach for three-dimensional seismic data interpolation[J].Journal of Engineering Geophysics,2022,19(3):410-418.
• [9] LARNER K.Trace interpolation and the design of seismic surveys[J].Geophysics,1981,46(4):407-415.
• [10]王本锋，韩东，李家阔.基于扩张卷积的智能化规则缺失炮插值重建方法[J].地球物理学报，2022,65(6):2226-2243.WANG Benfeng, HAN Dong, LI Jiakuo. Intelligent regularly missing shots reconstruction based on dilated convolution[J].Chinese Journal of Geophysics,2022,65(6):2226-2243.
• [11]金成玫，陈生昌.基于压缩感知和深度学习的地震数据重建[J].石油物探，2022,61(5):782-792.JIN Chengmei,CHEN Shengchang. Seismic data reconstruction based on compressed sensing and deep learning[J].Geophysical Prospecting for Petroleum,2022,61(5):782-792.
• [12] DUIJNDAM A J W,SCHONEWILLE M A,HINDRIKS C O H.Reconstruction of band-limited signals,irregularly sampled along one spatial direction[J]. Geophysics, 1999, 64(2):524-538.
• [13] HINDRIKS K,DUIJNDAM A J W.Reconstruction of 3-D seismic signals irregularly sampled along two spatial coordinates[J].Geophysics,2000,65(1):253-263.
• [14] LIU B, SACCHI M D. Minimum weighted norm interpolation of seismic records[J].Geophysics,2004,69(6):1560-1568.
• [15] XU S,ZHANG Y,PHAM D,et al.Antileakage Fourier transform for seismic data regularization[J].Geophysics,2005,70(4):87-95.
• [16] XU S,ZHANG Y,LAMBARéG.Antileakage Fourier transform for seismic data regularization in higher dimensions[J].Geophysics,2010,75(6):113-120.
• [17]刘仕友，宋维琪，闫安菊.基于稀疏重建分析的道集优化方法[J].石油物探，2023,62(2):297-304.LIU Shiyou, SONG Weiqi, YAN Anju. A gather optimization method based on a sparse reconstruction analysis[J]. Geophysical Prospecting for Petroleum,2023,62(2):297-304.
• [18]王天野.基于傅里叶变换的反泄漏地震数据重建方法[D].大庆：东北石油大学，2017.WANG Tianye.Research on seismic data reconstruction based on Fourier transform[D].Daqing:Northeast Petroleum University,2017.
• [19]陈建友，王晓凯，杨长春.一种均匀网格反泄露傅里叶变换的频率域高效实现方法[J].石油物探，2018,57(1):45-49.CHEN Jianyou,WANG Xiaokai,YANG Changchun.Efficient implementation of regular grid antileakage Fourier transform in the frequency domain[J]. Geophysical Prospecting for Petroleum,2018,57(1):45-49.
• [20] MALLAT S G,ZHANG Z.Matching pursuits with time-frequency dictionaries[J].IEEE Transactions on signal processing,1993,41(12):3397-3415.
• [21] NGUYEN T, WINNETT R. Seismic interpolation by optimally matched Fourier components[C]. SAN Antonio:SEG Annual Meeting,2011.
• [22]段文胜，王鹏，党青宁，等.应用匹配追踪傅里叶插值技术实现OVT域连片处理[J].石油地球物理勘探，2017,52(4):669-677,622.DUAN Wensheng, WANG Peng, DANG Qingning, et al. 5D data regularization based on matching pursuit Fourier interpolational for the OVT domain data merging processing[J].Oil Geophysical Prospecting,2017,52(4):669-677,622.