Abstract:Seismic data regularization and uneven subsurface illumination are two of key factors hindering amplitude-preserve imaging.Based on information of known subsurface model,original data are mapped to the image space with Kirchhoff migration operator.Then features are extracted in the image space.Finally the data are demigrated to the data space.Since the migration operator eliminates propagating effects to some extent,it is easier to extract features in the image space.In this paper,we propose an improved seismic regularizationmethod using an iterative algorithm based on least-squares inversion and analyze its dependence on velocity model.Considering the computational cost and flexibility,Kirchhoff migration and demigration are used as the mapping operators.Synthetic and real data examples show that the proposed method is able to reconstruct missing data and improve final images even in complex geological settings.
Verschuur D J,Berkhout A J.Estimation of multiple scattering by iterative inversion,part Ⅱ:Practical aspects and examples[J].Geophysics,1997,62(5):1596-1611.
[2]
Biondi B L.3-D Seismic Imaging[R].Stanford Exploration Project,2004.
[3]
Morice S,Ronen S,Canter P,et al.The impact of positioning differences on 4D repeatability[C].SEG Technical Program Expanded Abstracts,2000,19:1611-1614.
[4]
Trad D,Ulrych T,Sacchi D.Latest views of the sparse Radon transform[J].Geophysics,2003,68(1):386-399.
[5]
Liu B,Sacchi D.Minimum weighted norm interpola-tion of seismic records[J].Geophysics,2004,69(4):1560-1568.
[6]
Xu S,Zhang Y,Pham D,et al.Antileakage Fourier transform for seismic data regularization[J].Geophysics,2005,70(4):V87-V95.
[7]
王伟,陈双廷,王宝彬,等.五维规则化技术研究与应用[J].石油地球物理勘探,2017,52(增刊1):28-33.WANG Wei,CHEN Shuangting,WANG Baobin,et al.Application of 5D regularization in seismic data processing[J].Oil Geophysical Prospecting,2017,52(S1):28-33.
[8]
段文胜,王鹏,党青宁,等.应用匹配追踪傅里叶插值技术实现OVT域连片处理[J].石油地球物理勘探,2017,52(4):669-677.DUAN Wensheng,WANG Peng,DANG Qingning,et al.5D data regularization based on matching pursuit Fourier interpolation for the OVT domain data merging processing[J].Oil Geophysical Prospecting,2017,52(4):669-677.
[9]
崔兴福,刘东奇,张关泉.小波变换实现地震道内插[J].石油地球物理勘探,2003,38(增刊1):93-97.CUI Xingfu,LIU Dongqi,ZHANG Guanquan.Seismic interpolation via wavelet transform[J].Oil Geophysical Prospecting,2003,38(S1):93-97.
[10]
刘国昌,陈小宏,郭志峰,等.基于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.
[11]
王本锋,陆文凯,陈小宏,等.基于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.
[12]
Spitz S.Seismic trace interpolation in the f-x domain[J].Geophysics,1991,56(6):785-794.
[13]
Fomel S.Three Dimension Seismic Data Regulariza-tion[D].Stanford University,2001.
[14]
刘洋,张鹏,刘财,等.地震数据Bregman整形迭代插值方法[J].地球物理学报,2018,61(4):1400-1412.LIU Yang,ZHANG Peng,LIU Cai,et al.Seismic data interpolation based on Bregman shaping iteration[J].Chinese Journal of Geophysics,2018,61(4):1400-1412.
[15]
陈双全,王尚旭,季敏.基于信号保真的地震数据插值[J].石油地球物理勘探,2005,40(5):515-517.CHEN Shuangquan,WANG Shangxu,JI Min.Seismic data interpolation based on signal hi-fi[J].Oil Geophysical Prospecting,2005,40(5):515-517.
[16]
李海山,吴国忱,印兴耀.形态分量分析在地震数据重建中的应用[J].石油地球物理勘探,2012,47(2):236-243.LI Haishan,WU Guochen,YIN Xingyao.Morphological component analysis in seismic data reconstruction[J].Oil Geophysical Prospecting,2012,47(2):236-243.
[17]
刘玉金,李振春.局部平面波模型约束下的迭代加权最小二乘反演三维地震数据规则化[J].石油地球物理勘探,2012,47(3):418-424.LIU Yujin,LI Zhenchun.Three dimensional data re-gularization by iteratively re-weighted least-squares inversion based on the local plane wave model[J].Oil Geophysical Prospecting,2012,47(3):418-424.
[18]
祖绍环,周辉,陈阳康,等.不规则采样的多震源数据整形正则化分离方法[J].石油地球物理勘探,2016,51(2):247-253.ZU Shaohuan,ZHOU Hui,CHEN Yangkang,et al.Irregularly sampled multi-source data deblending based on shaping regularization[J].Oil Geophysical Prospecting,2016,51(2):247-253.
Bagaini C,Spagnolini U.2-D continuation operators and their applications[J].Geophysics,1996,61(6):1846-1858.
[21]
Biondi B,Fomel S and Chemingui N.Azimuth moveout for 3-D prestack imaging[J].Geophysics,1998,63(2):574-588.
[22]
Claerbout J F.Earth Soundings Analysis:Processing Versus Inversion[M].Blackwell Scientific Publications, 2004.
[23]
Wu D,Wang Y,Liu Y.Least-squares seismic data mapping with prestack time migration and demigration[C].SEG Technical Program Expanded Abstracts,2015,34:3931-3935.
[24]
刘玉金,李振春,吴丹,等.局部倾角约束最小二乘偏移方法研究[J].地球物理学报,2013,56(3):1003-1011.LIU Yujin,LI Zhenchun,WU Dan,et al.The research on local slope constrained least-squares migration[J].Chinese Journal of Geophysics,56(3):1003-1011.
[25]
Dai W,Boonyasiriwat C,Schuster G T.3D multi-source least-squares reverse time migration[C].SEG Technical Program Expanded Abstracts,2010,29:3120-3124.
[26]
Wong M,Ronen S,Biondi B.Least-squares reversetime migration/inversion for ocean bottom data:A case study[C].SEG Technical Program Expanded Abstracts,2011,30:3578-3583.
[27]
Symes W.Migration velocity analysis and waveform inversion[J].Geophysical Prospecting,2008,56(6):765-790.
[28]
Bourgeois A,Bourget M,Lailly P,et al.Marmousi model and data[C].EAGE Workshop of Practical Aspects of Seismic Data Inversion,1991,5-16.