A STUDY OF SEISMIC DATA RECOVERY BASED ON SPARSE TRANSFORM
LU Jiao-tong1, CAO Si-yuan2, DONG Jian-hua3, ZHANG Yan2
1. Geophysical Engineering Department,Sinopec International Petroleum Service Corporation,Beijing 100029,China;
2. College of Geophysics & Information Engineering,China University of Petroleum,Beijing 102249,China;
3. CNOOC Research Institute,Beijing 100027,China
The seismic data recovery from data with missing traces plays an important role in the later stage seismic processing.The authors studied the sparse transform (F-K transform and Curvelet transform) and popular compressed sensing theory,and then combined the two methods together to build the seismic data recovery model which is based on sparse transform.The F-K transform changes the seismic data from the t-x (time-space) domain into the f-k (frequency-wavenumber) domain.Because of the favorable directionality and locality and multidimensionality,the curvelet transform can represent the seismic data in a more compressible way.On the basis of the recovery model,the missed seismic data are recovered by the two sparse transforms and the recovery results are compared and analyzed.The recovery results prove that the Curvelet transform recovery can get the better reconstruction effect than the F-K transform.Finally the Marmousi2 model and practical seismic data are processed,and the result shows that the seismic data recovery model is correct and effective.
路交通, 曹思远, 董建华, 张. 基于稀疏变换的地震数据重构方法[J]. 物探与化探, 2013, 37(1): 175-179.
LU Jiao-tong, CAO Si-yuan, DONG Jian-hua, ZHANG Yan. A STUDY OF SEISMIC DATA RECOVERY BASED ON SPARSE TRANSFORM. Geophysical and Geochemical Exploration, 2013, 37(1): 175-179.
[1] Spitz S.Seismic trace interpolation in the f-x domain[J].Geophysics,1991,56(6):785-794.[2] Gulunay N, Chambers R E.Unaliased f-k domain trace interpolation(UFKI)//66th Ann. Internat. Mtg.,Soc Expl. Geophys,Expanded Abstracts,1996:1461.[3] Ronen J.Wave equation trace interpolation[J].Geophysics,1987,52(7):973-984.[4] 高建军,陈小宏,李景叶,等.基于非均匀Fourier变换的地震数据重建方法研究[J].地球物理学进展,2009,24(5):1741-1747.[5] Naghizadeh M,Sacchi M.Beyond alias hierarchical scale curvelet interpolation of regularly and irregularly sampled seismic data[J].Geophysics,2010,75(6):189-202.[6] Naghizadeh M,Sacchi M.Hierarchical scale curvelet interpolation of aliased seismic data[J].SEG International Exposition and 80th Annual Meeting,2010.[7] Herrmann J,Hennenfent G.Non-parametric seismic data recovery with curvelet frames[J].Geophysical Journal International,2008,173(1):A1-A5.[8] Donoho D.Compressed sensing[J].IEEE Trans. On Information Theory,2006,52(4):1289-1306.[9] Donoho D,Tsaig Y.Extensions of compressed sensing[J].Signal Processing,2006,86(3):533-548.[10] Berg E D,Friedlander M P.Probing the Pareto frontier for basis pursuit solutions[J].SIAM J. On Scientific Computing,2008,31(2):890-912.