Depth-domain wavelet estimation using the subspace-constrained Huber norm
ZHANG Jie1,2, CHEN Xuehua1,2, JIANG Wei1,2, DAN Zhiwei3, XIAO Wei3
1. State Key Laboratory of Oil & Gas Reservoir Geology and Exploitation, Chengdu University of Technology, Chengdu, Sichuan 610059, China; 2. Key Laboratory of Earth Exploration and Information Technology of Ministry of Education, Chengdu University of Technology, Chengdu, Sichuan 610059, China; 3. Data Processing Company, Geophysical Branch, China Oilfield Services Limited, CNOOC, Zhanjiang, Guangdong 524057, China
Abstract:Generally,in the constant-velocity depth domain,only a few hundred meters of logging information is available,and is equivalent to 2~5 times the length of a constant depth-domain seismic wavelet.It is difficult to estimate a reliable seismic wavelet from such a short data segment.To address this issue,we proposed a method for estimating the constant-velocity depth-domain wavelet using the subspace-constrained Huber norm.The corresponding procedure can be divided into three steps.First,the reflectivity and the seismic trace near the well location are transformed from the real depth domain to the constant-velocity depth domain.Next,set a threshold and generate the initial synthetic seismogram by convoluting the constant-velocity depth-domain reflectivity and the initial seismic wavelet.Finally,the seismic wavelet is updated by iterative least square method according to the residual error of the synthetic seismogram and the seismic trace until the iteration termination condition is reached.We compare the proposed method and the conventional least square-based methods through the synthetic and field seismic data.The results show that the high performance of our method for estimating the reliable wavelet from limited seismic and well-log data segment.The proposed method can be further extended to estiamte the depth-variant seismic wavelets.
李振春.地震偏移成像技术研究现状与发展趋势[J].石油地球物理勘探,2014,49(1):1-21.LI Zhenchun.Research status and development trends for seismic migration technology[J].Oil Geophysical Prospecting,2014,49(1):1-21.
[2]
罗红梅,王长江,刘书会,等.深度域高精度井震动态匹配方法[J].石油地球物理勘探,2018,53(5):997-1005.LUO Hongmei,WANG Changjiang,LIU Shuhui,et al.Well-to-seismic calibration in the depth domain using dynamic depth warping[J].Oil Geophysical Prospecting,2018,53(5):997-1005.
[3]
Fletcher R P,Archer S,Nichols D,et al.Inversion after depth imaging[C].SEG Technical Program Expanded Abstracts,2012,31:1-5.
[4]
周赏,汪关妹,张万福,等.深度域地震资料解释技术应用及效果[J].石油地球物理勘探,2017,52(增刊1):92-98.ZHOU Shang,WANG Guanmei,ZHANG Wanfu,et al.Depth-domain seismic data interpretation[J].Oil Geophysical Prospecting,2017,52(S1):92-98.
[5]
彭军,周家雄,王宇,等.基于褶积理论的深度域地震资料反演[J].科学技术与工程,2014,14(19):40-45.PENG Jun,ZHOU Jiaxiong,WANG Yu,et al.Depth domain seismic inversion base on convolution theory[J].Science Technology and Engineering,2014,14(19):40-45.
[6]
Letki L,Darke K,Borges Y.A comparison between time domain and depth domain inversion to acoustic impedance[C].SEG Technical Program Expanded Abstracts,2015,34:3376-3380.
[7]
Singh Y.Deterministic inversion of seismic data in the depth domain[J].The Leading Edge,2012,31(5):538-545.
何惺华.深度域地震资料若干问题初探[J].石油物探,2004,43(4):353-358.HE Xinghua.Discussion on seismic data of depth domain[J].Geophysical Prospecting for Petroleum,2004,43(4):353-358.
[10]
林伯香,胡中平,薛诗桂.利用变换深度域速度函数制作深度域合成地震记录[J].石油地球物理勘探,2006,41(6):640-643.LIN Boxiang,HU Zhongping,XUE Shigui.Using velocity function in transformed depth domain to make synthetic seismograph in depth domain[J].Oil Geophysical Prospecting,2006,41(6):640-643.
[11]
Hu Z,Zhu H,Lin B,et al.Wavelet analysis and convolution in depth domain[C].SEG Technical Program Expanded Abstracts,2007,26:2733-2737.
[12]
Lazear G.Mixed-phase wavelet estimation using fourth-order cumulants[J].Geophysics,1993,58(7):1042-1051.
[13]
Sacchi M D,Ulrych T J.Nonminimum-phase wavelet estimation using higher order statistics[J].The Lea-ding Edge,2000,19(1):80-83.
Mirko V D B,Pham D T.Robust wavelet estimation and blind deconvolution of noisy surface seismics[J].Geophysics,2008,73(5):V37-V46.
[17]
Zhang R,Deng Z.A depth variant seismic wavelet extraction method for inversion of poststack depth-domain seismic data[J].Geophysics,2017,82(6):R569-R579.
[18]
Danielson V,Karlsson T V.Extraction of signatures from seismic and well data[J].First Break,1984,2(4):15-21.
[19]
Lines L R,Treitel S.Wavelets,well logs and Wiener filters[J].First Break,1985,3(8):9-14.
[20]
梁光河.地震子波提取方法研究[J].石油物探,1998,37(1):31-39.LIANG Guanghe.On the methods of seismic wavelet extraction[J].Geophysical Prospecting for Petro-leum,1998,37(1):31-39.
[21]
Buland A,Omre H.Bayesian wavelet estimation from seismic and well data[J].Geophysics,2003,68(6):2000-2009.
[22]
Gunning J,Glinsky M E.Wavelet extractor:A Baye-sian well-tie and wavelet extraction program[J].Computers and Geosciences,2006,32(5):681-695.
[23]
Bianco E.Tutorial:Wavelet estimation for well ties[J].The Leading Edge,2016,35(6):541-543.
[24]
Bunch A W H,White R E.Least-squares filters wi-thout transient errors:an examination of the errors in least-squares filter design[J].Geophysical Prospecting,1985,33(5):657-673.
[25]
Ji J.Robust inversion using biweight norm and its application to seismic inversion[J].Exploration Geophysics,2012,43(2):70-76.
[26]
Huber P J.Robust regression: Asymptotics,conjec-tures,and Monte Carlo[J].The Annals of Statistics,1973,1(5):799-821.
[27]
Scales J A,Gersztenkorn A.Robust methods in inverse theory[J].Inverse Problems,1988,4(4):1071-1091.
[28]
Daubchies I,DeVore R,Fornasier M,et al.Iteratively re-weighted least squares minimization for sparse recovery[J].Communications on Pure and Applied Mathematics,2010,63(1):1-38.
[29]
Kazemi N,Sacchi M S.Sparse multichannel blind deconvolution[J].Geophysics,2014,79(5):V143-V152.