ARTICLE

Spectral decomposition with sparsity constraint and its application

HAN XU1 XINWEN WANG1 LI GAO2
Show Less
1 School of Earth Sciences and Recourses, China University of Geosciences, Beijing 100083, P.R. China,
2 Institute of Geomechanics, Chinese Academy of Geological Sciences, Beijing 100081, P.R. China,
JSE 2018, 27(1), 89–103;
Submitted: 2 December 2016 | Accepted: 10 November 2017 | Published: 1 February 2018
© 2018 by the Authors. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution -Noncommercial 4.0 International License (CC-by the license) ( https://creativecommons.org/licenses/by-nc/4.0/ )
Abstract

Xu, H., Wang, X. And Gao, L., 2018. Spectral decomposition with sparsity constraint and its application. Journal of Seismic Exploration, 27: 89-101. Spectral decomposition has been widely used in seismic signal processing and interpretation at present because it can reveal lots of valuable information hidden in the broadband seismic response. Unlike the conventional frequency-based methods, spectral decomposition is able to estimate the frequency contents of the signal at the specific time. How to seek an optimal solution is the most significant aspect for Spectral Decomposition with Sparsity Constraints (SDSC). In this paper, the Li regularized L2-norm is employed as the objective function, the Ricker wavelet is chosen to construct a wavelet library and the optimal solution is obtained by the Iterative Soft Thresholding Algorithm (ISTA). We apply SDSC to the synthetic and field examples. The former show that the SDSC method does a better job in both time and frequency resolution than the traditional spectral decomposition technique, namely the continuous wavelet transform (CWT) method, which always suffers from the conflict between time resolution and frequency resolution. Applications to field data further indicate the potential of the SDSC method in identifying the strong anomalies related to hydrocarbon, and detecting the variations in amplitude associated with faulting.

Keywords
spectral decomposition
sparsity constraint
continuous wavelet transform
hydrocarbon detection
fault identification
References
  1. Allen, J. B., 1977, Short term spectral analysis, synthetic and modification by discretefourier transform. IEEE Transactions on Acoustics, Speech, and Signal Processing,25, 235-238.
  2. Beck, A. and Teboulle, M., 2009, A fast iterative shrinkage-thresholding algorithm forlinear inverse problems. Siam J. Imaging Sci., 2: 183-202.
  3. Bonar, D. and Sacchi, M., 2010, Complex spectral decomposition via inversionstrategies. 80th Annual International Meeting, SEG, Expanded Abstracts:1408-1412.
  4. Bonar, D. and Sacchi, M., 2013, Spectral decomposition with f-x-y preconditioning.Geophys. Prospect., 61(suppl.1):152-165.
  5. Castagna, J. P., Sun, S. and Siegfried, R. W., 2003, Instantaneous spectral analysis:detection of low-frequency shadows associated with hydrocarbons. The LeadingEdge, 22: 120-127.
  6. Chen, S., Li, X. and Wu, X., 2014, Application of frequency-dependent AVO inversionto hydrocarbon detection. J. Seismic Explor., 23: 241-264.
  7. Chen, Y., Liu, T., Chen, X., Li, J. and Wang, E., 2014. Time-frequency analysis ofseismic data using synchrosqueezing wavelet transform. J. Seismic Explor., 23:303-312.
  8. Cohen, L., 1966, Generalized phase-space distribution functions. J. Math. Phys., 7:781-786.
  9. Daubechies, I., Defrise, M. and Mol, C. D., 2004, An iterative thresholding algorithmfor linear inverse problems with a sparsity constraint. Commun. Pur. Appl. Math.,57: 1413-1541.
  10. Donoho, D., 2006, For most large underdetermined systems of linear equations, theminimalell-1 norm solution is also the sparsest solution. Commun. Pur. Appl.Math., 59: 907-934.
  11. Farfour, M., Wang, J., Jo, Y. and Wan, K., 2013, Spectral decomposition and reservoirengineering data in mapping thin bed reservoir, Stratton field, South Texas. J.Seismic Explor., 22: 77-91.
  12. Jeffrey, C. and William, J., 1996, On the existence of discrete Wigner distributions.IEEE Signal Proc. Let., 6:304-306.
  13. Landweber, L., 1951, An iteration formula for Fredholm integral equations of the firstkind. American Journal of Mathematics, 73: 615-624.
  14. Liu, J. and Marfurt, K., 2007, Instantaneous spectral attributes to detect channels.Geophysics, 72: P23-P31.
  15. Liu, W., Cao, S., Liu, Y. and Chen, Y., 2016. Synchrosqueezing transform and itsapplications in seismic data analysis. J. Seismic Explor., 25: 27-44.
  16. Liu, Y. and Fomel, S., 2013, Seismic data analysis using local time-frequencydecomposition. Geophys. Prospect., 61:516-525.
  17. Mallat, S. and Zhang, Z., 1993, Matching pursuits with time frequency dictionaries.IEEE T. Signal Proces., 41: 3397-3415.
  18. Marfurt, K. and Kirlin, R. L., 2001, Narrow-band spectral analysis and thin-bed tuning.Geophysics, 66: 1274-1283.
  19. Odebeatu, E., Zhang, J. and Chapman, M., 2006, Application of spectral decompositionto detection of dispersion anomalies associated with gas saturation. The LeadingEdge, 25: 206-210.
  20. Partyka, G., Gridley, J. and Lopez, J., 1999, Interpretational application of spectraldecomposition in reservoir characterization. The Leading Edge, 22: 353-360.
  21. Portniaguine, O. and Castagna, J., 2004, Inverse spectral decomposition. 74th Annual
  22. International Meeting, SEG, Expanded Abstracts: 1786-1789.
  23. Radad, M., Gholami, A. and Siahkoohi, H., 2016, A fast method for generating highresolution single-frequency seismic attributes. J. Seismic Explor., 25: 11-25.
  24. Reine, C., M. van der Baan and R. Clark., 2009, The robustness of seismic attenuationmeasurements using fixed- and variable-window time-frequency transforms.Geophysics, 74: WA23-WA135.
  25. Rioul, O. and Vetterlil. M., 1991, Wavelet and signal processing. IEEE Signal Proc.Mag., 8: 14-37.
  26. Robinson, E.A. and Treitel, S., 1980, Geophysical Signal Analysis. Prentice-Hall, Inc.
  27. Selesnick, I. Sparse Signal Restoration.' Connexions (2009).
  28. Sinha, S., Routh, P., Anno, P. and Castagna, J., 2005, Spectral decomposition of seismicdata with continuous-wavelet transform. Geophysics, 70: 19-25.
  29. Stockwell, R., Mansinha, L. and Lowe. R., 1996, Localization of the complex spectrum:the S transform. IEEE Trans. Signal Process., 44: 998-1001.
  30. Wu, X., Chapman, M., Wilson, A. and Li, X., Estimating seismic dispersion fromprestack data suing frequency-dependent AVO analysis. J. Seismic Explor., 23:425-429.
  31. Wu, X. and Liu, T., 2010. Analysis of seismic spectral attenuation based on
  32. Wigner-Wille distribution for sandstone reservoir characterization- a case studyfrom West Sichuan Depression, China. J. Seismic Explor., 19: 69-85.JOURNAL OFSEISMIC EXPLORATIONVolume 27, 2018EDITORS-IN-CHIEF :M. Nafi Tokséz (Cambridge, MA, U.S.A.)Arthur B. Weglein (Houston, TX, U.S.A.)EDITORIAL BOARD
  33. L. Amundsen (Trondheim, Norway) E. Liu (Houston, TX, U.S.A.)
  34. J.M. Carcione (Trieste, Italy) R. Marschall (Hannover, Germany)
  35. P. Cary (Calgary, AB, Canada) G.A. McMechan (Richardson, TX, U.S.A.)
  36. Y.K. Chen (Oak Ridge, TN, U.S.A.) M. Meier (Houston, TX, U.S.A.)
  37. W.K. Chung (Busan, South Korea) F. Poletto (Trieste, Italy)
  38. V. Cheverda (Novosibirsk, Russia) M. Sacchi (Edmonton, AB, Canada)
  39. P.F. Daley (Calgary, AB, Canada) M. Sen (Austin, TX, U.S.A.)
  40. E. Filpo (Rio de Janeiro, Brazil) G. Rubino (London, ON, Canada)
  41. D. Foster (Houston, TX, U.S.A.) S. Shapiro (Berlin, Germany)
  42. B. Gurevich (Perth, Australia) J. Sheiman (Houston, TX, U.S.A.)
  43. W. Ha (Busan, South Korea) C.S. Shin (Seoul, South Korea)
  44. M. Hall (Calgary, AB, Canada) W. Séllner (Oslo, Norway)
  45. B.A. Hardage (Austin, TX, U.S.A.) P.L. Stoffa (Austin, TX, U.S.A.)
  46. T. Hertweck (Karlsruhe, Germany) L. Thomsen (Houston, TX, U.S.A.)
  47. L.T. Ikelle (Bryan, TX, U.S.A.) I. Tsvankin (Golden, CO, U.S.A.)
  48. M. Jacobsen (Bergen, Norway) C. Vanelle (Hamburg, Germany)
  49. G. Lambaré (Paris, France) A.R. Verdel (Utrecht, Netherlands)
  50. E. Landa (Tel Aviv, Israel) D.J. Verschuur (Delft, Netherlands)
  51. L. Lines (Calgary, AB, Canada) D.H. Yang (Beijing, P.R. China)
  52. J.E. Lira (Rio de Janeiro, Brazil) J. Zhang (Houston, TX, U.S.A.)Scope
  53. The Journal of Seismic Exploration is an international medium for the publication of research inseiusmic modeling, processing, inversion, interpretation, field techniques , borehole techniques,tomography, instrumentation and software.Policy
  54. Short publication time, high scientific level, easily readable and concise papers are the mainobjectives. This will be achieved by fast review (4 months) and fast publication (4 months).
  55. THERE IS NO PAGE CHARGE. Authors receive a PDF of their paper, instead of reprints.Subscriptions
  56. The 2018 subscription price is £ 775, inclusive of postage for the print edition. In 2018 Vol. 27will be published, containing around 600 pages, in 6 issues. Online subscription price is £1,165,composed of the print + online version. Subscriptions will be continued automatically, unlesscancelled before December 1, 2018.
  57. Subscription orders should be sent to: Geophysical Press Ltd., Urb. Los Flamingos — Tee 6,29679 Benahavis, Spain. E-mail: jse@geophysical-press.com ‧ www.geophysical-press.comGUIDE FOR AUTHORS
  58. Authors are invited to submit manuscripts of no more than 20 pp. to:
  59. Journal of Seismic Exploration, Editorial Office, Urb. Los Flamingos, 29679 Benahavis,
  60. Spain. As PDF file or MS Word file to: jse@geophysical-press.comLanguage
  61. The journal will only contain articles in the English language.Preparation of the texta) The manuscript should be typewritten with double spacing and wide margins and includeat the beginning of the paper an abstract of not more than 500 words. Words to be printedin italics should be underlined. The metric system should be used throughout.b) The title page should include: the title, the name(s) of the author(s), and their affiliations.Referencesa} References in the text with the name of the author(s) without initial{s), followed by thepublication date in brackets.b) {f reference is made in the text to publications written by more than two authors, the nameof the first author should be used, followed by “et al.'. This indication, however, shouldnever be used in the list of references. Here the names of all authors and co-authors shouldbe mentioned.c) The reference list should be in alphabetical order and on sheets separate from the text.Tables
  62. Tables should be compiled as separate sheets. A title should be provided for each table andthey should be referred to in the text.Illustrationsa) All illustrations should be numbered consecutively and referred to in the text.b) Drawings should be completely lettered, the size of the lettering being appropriate to thatof the drawings, but taking into account the possible need for reduction in size (of not morethan 50%). The page format of the journal should be considered in designing the drawings.c) Photographs must be of good quality, printed on glossy paper.d) Figure captions should appear on the same page as the corresponding illustrations.Proofs
  63. One set of proofs wiil be sent to the author, to be checked for printer’s errors.Reprints50 reprints of each article are supplied free of charge. Additional reprints can be orderedtogether with the return of the proofs.Submission of manuscripts
  64. Papers should be submitted to the Editorial Office of this journal, as a PDF or MS Word file.
  65. Submission of an article is understood to imply that the article is original and unpublished andis not being considered for publication elsewhere.THERE IS NO PAGE CHARGEJOURNAL OFSEISMIC EXPLORATIONVolume 27, Number 1O. Giireli and T. KayiranM. Bai, J. Wu, J. Xie andD. ZhangA. Dahroug, S. Mahmoud,E. Morsy, H. Nada andH. KarsliG. Li, Y. Li and X. LuY.T. Zhou and W.X. HanH. Xu, X.W. Wang andL. GaoFebruary 2018CONTENTSOn the non-uniqueness of therefraction solution.............Least-squares reverse timemigration of blended data withlow-rank constraint alongstructural direction.............Contribution of inverted seismicdata in extracting attributes tocharacterize Pliocene agedchannel located off-shore theNile Delta, Egypt.............Seismic events detection in stronglow-frequency backgroundnoise by complex shock filter. . . .Multiples attenuation in thepresence of blending noise. .... .Spectral decomposition withsparsity constraint and itsapplication ...................40 TVNXNnOTxa回圖區x>S之(8I0Z) ZOI-T ‘LZ TOAVolume 27, Number 1, February 2018ISSN 0963-0651JOURNAL OFSEISMIC EXPLORATIONIMAGING TECHNOLOGY IN APPLIED GEOSCIENCES
  66. EDITORS-IN-CHIEF : M. NAFI TOKSOZ and ARTHUR B. WEGLEINGEOPHYSICAL PRESS
Share
Back to top
Journal of Seismic Exploration, Electronic ISSN: 0963-0651 Print ISSN: 0963-0651, Published by AccScience Publishing