|
|
Seismic reflectivity inversion based on L1-L1-norm sparse representation |
Zhan-Zhan SHI1,2, Yan-Qing XIA1, Huai-Lai ZHOU2, Yuan-Jun WANG2, Xiang-Rong TANG2 |
1. The Engineering & Technical College of Chengdu University of Technology,Leshan 614000,China; 2. College of Geophysics,Chengdu University of Technology,Chengdu 610059,China |
|
|
Abstract High-resolution seismic inversion is confronted with two problems:First,seismic inversion is an ill-posed problem and has multiplicity of solutions,and second,noise and distortion are generated in the flows of acquisition and processing to reduce the stability of the inversion algorithm.Aimed at solving these two problems, this paper proposes an inversion method of seismic reflectivity based on L1-L1-norm sparse representation.Firstly,the L1-norm regularization term is used to reduce the inversion multiplicity,and then the L1-norm fitting term is used to enhance the noise robustness.The wavelet is extracted by well logging and seismic data to construct the over-complete wedge wavelet dictionary,and then the seismic signal is sparsely decomposed by the L1-L1-norm sparse representation,so as to realize the high-resolution reflectivity inversion.The experimental results of wedge model and actual seismic data show that the inversion algorithm is stable and has good noise robustness,and the inversion results are accurate and credible through logging data calibration.
|
Received: 29 October 2018
Published: 15 August 2019
|
|
|
|
|
|
Wedge geological model a—reflection coefficient model;b—synthetic seismic section
|
|
Comparison of random noise sensitivity a、b、c、d—inversion results of traditional algorithms;e、f、g、h—inversion results of L1-L1 norm sparse representation;random noise intensities from top to bottom are 0%,1%,5% and 10%,respectively
|
|
Comparison of outlier sensitivity a、b、c、d—inversion results of traditional algorithms;e、f、g、h—inversion results of L1-L1 norm sparse representation;the number of outliers from top to bottom is 0%,1%,5%,and 10%,respectively
|
|
Comparative analysis of inverted reflectivity of wedge model a—traditional algorithm inversion results under the condition of 1% random noise and 1% outliers;b—L1-L1 norm sparse representation inversion results under the condition of 1% random noise and 1% outliers;c—traditional algorithm inversion results under the condition of 5% random noise and 5% outliers;d—L1-L1 norm sparse representation inversion results under the condition of 5% random noise and 5% outliers
|
|
Vertical section (a) and its amplitude spectrum (b) extracted from a 3D seismic data set intersect well L1
|
|
Comparison of inverted reflectivity section a—traditional algorithms;b—L1-L1 norm system sparse representation
|
|
Comparison of inverted impedance sections a—traditional algorithms;b—L1-L1 norm system sparse representation
|
|
Zoom of inverted impedance sections a—traditional algorithms;b—L1-L1 norm system sparse representation
|
[1] |
Widess M B . How thin is a thin bed?[J]. Geophysics, 1973,38(6):1176-1180.
|
[2] |
Kallweit R S . The limits of resolution zero-phase wavelets[J]. Geophysics, 1982,47(7):1035-1046.
|
[3] |
Voogd N D, Rooijen H D . Thin-layer response and spectral bandwidth[J]. Geophysics, 1983,48(1):12-18.
|
[4] |
Chung H M, Lawton D C . Amplitude responses of thin beds; sinusoidal approximation versus Ricker approximation[J]. Geophysics, 1995,60(1):223-230.
|
[5] |
Van Riel P . Resolution in seismic trace inversion by parameter estimation[J]. Geophysics, 1985,50(9):1440-1455.
|
[6] |
Russell B H . Introduction to Seismic Inversion Methods[M]. Society of Exploration Geophysicists, 1988.
|
[7] |
Schuster G T . Seismic Inversion[M]. Society of Exploration Geophysicists, 2017.
|
[8] |
Partyka G . Interpretational applications of spectral decomposition in reservoir characterization[J]. Leading Edge, 1999,18(3):173-184.
|
[9] |
Marfurt K J, Kirlin R L . Narrow-band spectral analysis and thin-bed tuning[J]. Geophysics, 2001,66(4):1274-1283.
|
[10] |
Puryear C I, Castagna J P . An algorithm for calculation of bed thickness and reflection coefficients from amplitude spectrum[J]. SEG Technical Program Expanded Abstracts, 2006: 1767-1770.
|
[11] |
Puryear C I, Castagna J P . Layer-thickness determination and stratigraphic interpretation using spectral inversion:Theory and application[J]. Geophysics, 2008,73(2):R37-R48.
|
[12] |
Puryear C I, Portniaguine O N, Cobos C M , et al. Constrained least-squares spectral analysis:Application to seismic data[J]. Geophysics, 2012,77(5):V143-V167.
|
[13] |
Nguyen T, Castagna J . High-resolution reflectivity inversion[J]. Journal of Seismic Exploration, 2010,19(4):303-320.
|
[14] |
Zhang R, Castagna J . Seismic sparse-layer reflectivity inversion using basis pursuit decomposition[J]. Geophysics, 2011,76(6):R147-R158.
|
[15] |
Mallat S G, Zhang Z . Matching pursuits with time-frequency dictionaries[J]. IEEE Transactions on Signal Processing, 1993,41(12):3397-3415.
|
[16] |
Chen S S, Donoho D L, Saunders M A . Atomic decomposition by basis pursuit[J]. Siam Review, 2001,43(1):129-159.
|
[17] |
Beck A, Teboulle M . A fast iterative Shrinkage-Thresholding algorithm for linear inverse problems[J]. Siam Journal on Imaging Sciences, 2009,2(1):183-202.
|
[18] |
Becker S, Bobin J, Candès E J . NESTA:A fast and accurate first-order method for sparse recovery[J]. Siam Journal on Imaging Sciences, 2009,4(1):1-39.
|
[19] |
Zhang Z . Integrating globality and locality for robust representation based classification[J]. Mathematical Problems in Engineering, 2014,2014:1-10.
|
[20] |
Tropp J A, Gilbert A C, Strauss M J . Algorithms for simultaneous sparse approximation.Part I: Greedy pursuit[J]. Signal Processing, 2006,86(3):572-588.
|
[21] |
Tropp J A . Algorithms for simultaneous sparse approximation:part II:Convex relaxation[J]. Signal Processing, 2006,86(3):589-602.
|
[22] |
Zhang Z, Xu Y, Yang J , et al. A survey of sparse representation:algorithms and applications[J]. IEEE Access, 2017,3:490-530.
|
[23] |
Fuchs J J. Fast implementation of a l1-l1 regularized sparse representations algorithm[C]// IEEE International Conference on Acoustics,Speech and Signal Processing IEEE, 2009: 3329-3332.
|
[24] |
Yang J, Zhang Y . Alternating direction algorithms for l1-problems in compressive sensing[J]. Siam Journal on Scientific Computing, 2009,33(1):250-278.
|
[25] |
Claerbout J F, Muir F . Robust modeling with erratic data[J]. Geophysics, 1973,38(5):826-844.
|
[26] |
Morozov I B . Geometrical attenuation,frequency dependence of Q,and the absorption band problem[J]. Geophysical Journal International, 2008,175(1):239-252.
|
[27] |
Cheng H X, Kennett B L N . Frequency dependence of seismic wave attenuation in the upper mantle beneath the Australian region[J]. Geophysical Journal International, 2010,150(1):45-57.
|
[28] |
Zhou H. Practical seismic data processing[M]. Cambridge University Press, 2014.
|
[29] |
Zhang F, Dai R, Liu H . Seismic inversion based on L1-norm misfit function and total variation regularization[J]. Journal of Applied Geophysics, 2014,109:111-118.
|
[30] |
Gabay D, Mercier B . A dual algorithm for the solution of nonlinear variational problems via finite element approximation[J]. Computers & Mathematics with Applications, 1976,2(1):17-40.
|
[31] |
Boyd S, Parikh N, Chu E , et al. Distributed optimization and statistical learning via the alternating direction method of multipliers[J]. Foundations & Trends in Machine Learning, 2011,3(1):1-122.
|
[32] |
He B, Yuan X . On the O(1/n) convergence rate of the douglas-rachford alternating direction method[J]. Siam Journal on Numerical Analysis, 2012,50(2):700-709.
|
[33] |
Chung H, Lawton D C . Frequency characteristics of seismic reflections from thin beds[J]. Canadian Journal of Exploration Geophysics, 1995,31(1):8-32.
|
[1] |
Yong LUO, Hai-Bo MAO, Xiao-Hai YANG, Wen-Jie LI, Wen-Chao CHEN. Seismic random seismic noise attenuation method on basis of the double sparse representation[J]. Geophysical and Geochemical Exploration, 2018, 42(3): 608-615. |
[2] |
GUO Qi, ZENG Zhao-Fa, YU Chen-Xia, ZHANG Si-Meng. Seismic random noise suppression based on the high-precision dictionary learning algorithm[J]. Geophysical and Geochemical Exploration, 2017, 41(5): 907-913. |
|
|
|
|