Super-Resolution Inversion of Non-Stationary Seismic Traces

Abstract

In reflection seismology, the inversion of subsurface reflectivity from the observed seismic traces (super-resolution inversion) plays a crucial role in target detection. Since the seismic wavelet in reflection seismic data varies with the travel time, the reflection seismic trace is non-stationary. In this case, a relative amplitude-preserving super-resolution inversion has been a challenging problem. In this paper, we propose a super-resolution inversion method for the non-stationary reflection seismic traces. We assume that the amplitude spectrum of seismic wavelet is a smooth and unimodal function, and the reflection coefficient is an arbitrary random sequence with sparsity. The proposed method can obtain not only the relative amplitude-preserving reflectivity but also the seismic wavelet. In addition, as a by-product, a special Q field can be obtained.

The proposed method consists of two steps. The first step devotes to making an approximate stabilization of non-stationary seismic traces. The key points include: firstly, dividing non-stationary seismic traces into several stationary segments, then extracting wavelet amplitude spectrum from each segment and calculating Q value by the wavelet amplitude spectrum between adjacent segments; secondly, using the estimated Q field to compensate for the attenuation of seismic signals in sparse domain to obtain approximate stationary seismic traces. The second step is the super-resolution inversion of stationary seismic traces. The key points include: firstly, constructing the objective function, where the approximation error is measured in L₂ space, and adding some constraints into reflectivity and seismic wavelet to solve ill-conditioned problems; secondly, applying a Hadamard product parametrization (HPP) to transform the non-convex problem based on the L_p(0<p<1) constraint into a series of convex optimization problems in L₂ space, where the convex optimization problems are solved by the singular value decomposition (SVD) method and the regularization parameters are determined by the L-curve method in the case of single-variable inversion. In this paper, the effectiveness of the proposed method is demonstrated by both synthetic data and field data.

Keywords

Non-stationary seismic trace / reflectivity inversion / super-resolution / attenuation compensation / L_p regularization

References

[1] E. J. Candès and C. Fernandez-Granda, Towards a mathematical theory of super-resolution, Communications on Pure and Applied Mathematics, 67 (2014), 906-956.

[2] M. D. Sacchi and D. R. Velis and A. H. Comínguez, Minimum entropy deconvolution with frequency-domain constraints, Geophysics, 59 (1994), 938-945.

[3] M. D. Sacchi, Reweighting strategies in seismic deconvolution, Geophysical Journal International, 129 (1997), 651-656.

[4] M. Sajid and D. Ghosh, A fast and simple method of spectral enhancement, Geophysics, 79 (2014), V75-V80.

[5] S. A. M. Oliveira, and W. M. Lupinacci, L₁ norm inversion method for deconvolution in attenuating media, Geophysical Prospecting, 61 (2013), 771-777.

[6] S. Y. Yuan, S. X. Wang, M. Ma, Y. Z, Ji and L. Deng, Sparse Bayesian learning-based timevariant deconvolution, IEEE Transactions on Geoscience and Remote Sensing, 55 (2017), 6182-6194.

[7] H. L. Chen, J. H. Gao, N. H. Liu and Y. Yang, Multitrace semiblind nonstationary deconvolution, IEEE Geoscience and Remote Sensing Letters, 16 (2019), 1195-1199.

[8] N. D. Hargreaves and A. J. Calvert, Inverse Q filtering by Fourier transform, Geophysics, 56 (1991), 519-527.

[9] Y. H. Wang, A stable and efficient approach of inverse Q filtering, Geophysics, 67 (2002), 657-663.

[10] Y. H. Wang, Quantifying the effectiveness of stabilized inverse Q filtering, Geophysics, 68 (2003), 337-345.

[11] S. D. Stainsby and M. H. Worthington, Q estimation from vertical seismic profile data and anomalous variations in the central North Sea, Geophysics, 50 (1985), 615-626.

[12] S. W. Patton, Robust and least-squares estimation of acoustic attenuation from well-log data, Geophysics, 53 (1988), 1225-1232.

[13] E. Blias, Accurate interval Q-factor estimation from VSP data, Geophysics, 77 (2012), WA149WA156.

[14] C. L. Hackert and J. O. Parra, Improving Q estimates from seismic reflection data using well-log-based localized spectral correction, Geophysics, 69 (2004), 1521-1529.

[15] J. H. Gao, S. L. Yang, D. X. Wang and R. S. Wu, Estimation of quality of factor Q from the instantaneous frequency at the envelope peak of a seismic signal, Journal of Computational Acoustics, 19 (2011), 155-179.

[16] L. L. Wang, J. H. Gao, W. Zhao and X. D. Jiang, Enhancing resolution of nonstationary seismic data by molecular-Gabor transform, Geophysics, 78 (2012), V31-V41.

[17] Y. H. Wang, Inverse Q-filter for seismic resolution enhancement, Geophysics, 71 (2006), V51V60.

[18] I. L. S. Braga and F. S. Moraes, High-resolution gathers by inverse Q filtering in the wavelet domain, Geophysics, 78 (2013), V53-V61.

[19] I. Morozov, M. Haiba and W. B. Deng, Inverse attenuation filtering, Geophysics, 83 (2018), V135-V147.

[20] G. W. Zhang and J. H. Gao, Inversion-driven attenuation compensation using synchrosqueezing transform, IEEE Geoscience and Remote Sensing Letters, 15 (2018), 132-136.

[21] E. J. Candès and M. B. Wakin, Enhancing sparsity by reweighted L₁ minimization, Journal of Fourier Analysis and Applications, 14 (2008), 877-905.

[22] Z. B. Xu, H. Zhang, Y. Wang and X. Y. Chang, ${\mathbf{L}}_{1/2}$ regularization, Science China Information Sciences, 53 (2010), 1159-1169.

[23] L. L. Wang, Q. Zhao, J. H. Gao, Z. B. Xu, M. Fehler and X. D. Jiang, Seismic sparse-spike deconvolution via Toeplitz-sparse matrix factorization, Geophysics, 81 (2016), V169-V182.

[24] R. Chartrand, Exact reconstruction of sparse signals via nonconvex minimization, IEEE Signal Processing Letters, 14 (2007), 707-710.

[25] M. J. Lai, Y. Y. Xu and W. T. Yin, Improved iteratively reweighted least squares for unconstrained smoothed l_q minimization, SIAM Journal on Numerical Analysis, 51 (2013), 927957.

[26] J. Woodworth and R. Chartrand, Compressed sensing recovery via nonconvex shrinkage penalties, Inverse Problems, 32 (2016), 075004.

[27] H. L. Chen, J. H. Gao and B. Zhang, An adaptive time-varying seismic super-resolution inversion based on Lp regularization, IEEE Geoscience and Remote Sensing Letters, 2020, doi: 10.1109/LGRS.2020.3000339.

[28] M. Elad, Sparse and Redundant Representations: From Theory to Applications in Signal and Image Processing, Springer Publishing Company, Incorporated, 2010.

[29] P. D. Hoff, Lasso, fractional norm and structured sparse estimation using a Hadamard product parametrization, Computational Statistics and Data Analysis, 115 (2017), 186-198.

[30] P. C. Hansen, REGULARIZATION TOOLS: A Matlab package for analysis and solution of discrete ill-posed problems, Numerical Algorithms, 6 (1994), 1-35.

[31] P. C. Hansen, Analysis of discrete ill-posed problems by means of the L-curve, SIAM Review, 34 (1992), 561-580.

[32] G. H. Golub, M. Heath and G. Wahba, Generalized cross-validation as a method for choosing a good ridge parameter, Technometrics, 21 (1979), 215-223.

[33] J. H. Gao, B. Zhang, W. M. Han, J. G. Peng and Z. B. Xu, A new approach for extracting the amplitude spectrum of the seismic wavelet from the seismic traces, Inverse Problems, 33 (2017), 085005.

[34] J. P. Grossman, G. F. Margrave and M. P. Lamoureux, Adaptive, nonuniform Gabor frames from partitions of unity, CREWES Research Report, 14 (2002), 1-10.

[35] Y. L. Quan and J. M. Harris, Seismic attenuation tomography using frequency shift method, Geophysics, 62 (1997), 895-905.

[36] Y. H. Wang and J. Guo, Modified Kolsky model for seismic attenuation and dispersion, Journal of Geophysics and Engineering, 1 (2004), 187.

[37] R. Saab, R. Chartrand and O. Yilmaz, Stable sparse approximations via nonconvex optimization, IEEE International Conference on Acoustics, Speech and Signal Processing, 2008, 3885-3888.

[38] J. K. Pant, W. Lu and A. Antoniou, New improved algorithms for compressive sensing based on l_p norm, IEEE Transactions on Circuits and Systems II: Express Briefs, 61 (2014), 198-202.

[39] R. Chartrand and W. T. Yin, Iteratively reweighted algorithms for compressive sensing, IEEE International Conference on Acoustics, Speech and Signal Processing, 2008, 3869-3872.

[40] Z. B. Xu, H. L. Guo, Y. Wang and H. Zhang, Representative of L_1/2 regularization among L_q(0&q≤1) regularizations: An experimental study based on phase diagram, Acta Automatica Sinica, 38 (2012), 1225-1228.

[41] H. L. Chen, S. Y. Cao, S. Y. Yuan, X. Pan, S. Q. Chen and S. A. Shen, Nonstationary Sparse Reflectivity Inversion with EADTV Regularization, Eage Conference and Exhibition, 2017.

[42] S. Arridge, P. Maass, O. Öktem and C. B. Schönlieb, Solving inverse problems using datadriven models, Acta Numerica, 28 (2019), 1-174.

[43] X. Xu, J. H. Gao, B. Zhang, H. L. Chen and Y. Yang, Multichannel reflectivity inversion with sparse group regularization based on HPPSG algorithm, IEEE Geoscience and Remote Sensing Letters, 17 (2020), 784-788.

[44] C. A. Richard, B. Brian and H. T. Clifford, Parameter Estimation and Inverse Problems (Second Edition), Academic Press, 2011.