1. College of Geophysics, China University of Petroleum (Beijing), Beijing 102249, China; 2. State Key Laboratory of Petroleum Resources and Prospecting, Beijing 102249, China
Abstract:Although acoustic impedance inversion technology is quite mature, there are still some issues such as the ill-posedness of inverse problems, low resolution of inversion results, and the inability to delineate the stratigraphic boundaries. To this end, a seismic “blocky” acoustic impedance inversion method based on L1-2 regularization is proposed. Based on existing research, this paper introduces L1-2 regularization into the model-based acoustic impedance inversion and directly obtains the acoustic impedance inversion results from post-stack seismic data according to the idea of total variation regularization. First, the paper deduces a linearized forward modeling equation of acoustic impedance and analyzes its accuracy. Then, based on Bayesian theory, the paper constructs the objective function of the acoustic impedance inversion by the L1-2 regularization and solves the function through an iterative reweighted least squares (IRLS) algorithm to obtain the acoustic impedance inversion results. Since the acoustic impedance inversion is a single-trace inversion method, when it is applied to multi-trace data inversion, there is a spatial discontinuity. Therefore, an f-x space predictive filtering method is used to alleviate the discontinuity caused by noise and single-trace inversion. The quantitative comparison of correlation coefficients proves that the inversion results obtained by the L1-2 norm are better than those obtained by L1 and L2 norms, and synthetic and field data inversion examples demonstrate the effectiveness and feasibility of the proposed method.
MAURYA S P,SARKAR P.Comparison of post stack seismic inversion methods: A case study from Blackfoot Field, Canada[J]. International Journal of Scientific and Engineering Research, 2016, 7(8): 1091-1101.
[3]
张凌远, 张宏兵, 尚作萍, 等. 基于Zoeppritz方程的叠前和叠后混合多参数非线性地震反演[J].石油地球物理勘探,2021,56(1): 164-171.ZHANG Lingyuan,ZHANG Hongbing,SHANG Zuoping, et al.Nonlinear multi-parameter hybrid inversion of pre-stack and post-stack seismic data based on Zoeppritz equation[J].Oil Geophysical Prospecting,2021,56(1): 164-171.
[4]
周东红,李景叶,陈莉.基于弹性波动方程的叠后地震反演方法[J].石油地球物理勘探, 2018, 53(2): 395-402.ZHOU Donghong,LI Jingye,CHEN Li. Poststack inversion based on the elastic wave equation[J].Oil Geophysical Prospecting,2018,53(2): 395-402.
[5]
LINDSETH R O.Synthetic sonic logs: A process for stratigraphic interpretation[J].Geophysics,1979,44(1): 3-26.
[6]
邹振桓,杨文采.地震道的广义线性反演[J].石油地球物理勘探,1987,22(4): 363-375.ZOU Zhenhuan,YANG Wencai.Generalized linear inversion of seismic traces[J]. Oil Geophysical Prospecting, 1987,22(4): 363-375.
[7]
LI S,PENG Z P.Seismic acoustic impedance inversion with multi-parameter regularization[J].Journal of Geophysics and Engineering,2017,14(3): 520-532.
[8]
冉喜阳,周怀来,张益明,等.一种基于稀疏窗S变换的分频—重构波阻抗反演方法[J].中国海上油气,2019,31(4): 75-84.RAN Xiyang,ZHOU Huailai,ZHANG Yiming,et al.A frequency division-reconstruction wave impedance inversion method based on sparse window S transform[J]. China Offshore Oil and Gas,2019,31(4): 75-84.
[9]
郭朝斌,杨小波,陈红岳,等.约束稀疏脉冲反演在储层预测中的应用[J].石油物探,2006,45(4): 397-400.GUO Chaobin,YANG Xiaobo,CHEN Hongyue,et al. Constrained sparse pulse inversion research in north of Haitongji depression[J].Geophysical Prospecting for Petroleum,2006,45(4): 397-400.
[10]
HAMID H,PIDLISECKY A.Multitrace impedance inversion with lateral constraints[J].Geophysics,2015,80(6): M101-M111.
[11]
WANG L Q,ZHOU H,WANG Y F,et al.Three-parameter prestack seismic inversion based on L1-2 minimization[J].Geophysics,2019,84(5): R753-R766.
[12]
DONOHO D L.Compressed sensing[J]. IEEE Transactions on Information Theory, 2006, 52(4): 1289-1306.
刘国昌,蔡加铭,闫海洋,等.利用稀疏约束非平稳多项式回归去除地震噪声及拾取初至[J].石油地球物理勘探, 2020, 55(3): 548-556.LIU Guochang,CAI Jiaming,YAN Haiyang,et al.Using sparse-constrained nonstationary polynomial regression to remove seismic noises and picking up first arrival[J].Oil Geophysical Prospecting,2020,55(3): 548-556.
[15]
程三,张志勇,周峰,等.非结构化网格渐进自适应正则化反演算法[J].石油地球物理勘探,2022,57(2): 467-477.CHENG San, ZHANG Zhiyong, ZHOU Feng, et al. Study on step-by-step regularization inversion based on adaptive unstructured mesh[J].Oil Geophysical Prospecting, 2022, 57(2): 467-477.
[16]
CLAERBOUT J F,MUIR F.Robust modeling with erratic data[J]. Geophysics,1973,38(5): 826-844.
[17]
TIKHONOV A N, GONCHARSKY A V, STEPANOV V V, et al. Numerical Methods for the Solution of Ill-Posed Problems[M]. Springer, Nether-lands, 1995.
TARANTOLA A.Inverse Problem Theory and Methods for Model Parameter Estimation[M]. Society for Industrial and Applied Mathematics, 2005.
[20]
ZHANG F C,DAI R H,LIU H Q.Seismic inversion based on L1-norm misfit function and total variation regularization[J].Journal of Applied Geophysics,2014,109(10): 111-118.
[21]
SUN J J,LI Y G.Adaptive Lp inversion for simultaneous recovery of both blocky and smooth features in a geophysical model[J].Geophysical Journal International, 2014,197(2): 882-899.
[22]
WANG Y F, MA X, ZHOU H,et al.L1-2 minimization for exact and stable seismic attenuation compensation[J]. Geophysical Journal International,2018,213(3): 1629-1646.
[23]
CANALES L L.Random noise reduction[C].SEG Technical Program Expanded Abstracts,1984,3: 525-527.
[24]
SPITZ S.Seismic trace interpolation in the F-X domain[J]. Geophysics,1991,56(6): 785-794.
[25]
姚建阳.用地震道积分法提高地层的识别能力[J].石油物探,1990,29(1): 40-49.YAO Jianyang.Enhancement of the discernibility on formation by seismic trace integration[J].Geophysical Prospecting for Petroleum,1990, 29(1): 40-49.