• journal_title：Geophysics
• Contributor：Rui Zhang ; Mrinal K Sen ; Sanjay Srinivasan
• Publisher：Society of Exploration Geophysicists
• Date：2013-01-01
• Format：text/html
• Language：en
• Identifier：10.1190/geo2011-0502.1
• journal_abbrev：Geophysics
• issn：0016-8033
• volume：78
• issue：1
• firstpage：R1
• section：Seismic Inversion

Resolving thin layers and clearly delineating layer boundaries in inverted seismic sections are very important goals for exploration and production. Many seismic inversion methods based on a least-squares optimization approach with Tikhonov-type regularization can lead to unfocused transitions between adjacent layers. A basis pursuit inversion (BPI) algorithm based on the <mml:math display="inline"><mml:msub><mml:mi>L</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:math>L1 norm optimization method can, however, resolve sharp boundaries between layers. We have formulated a BPI algorithm for amplitude-versus-angle inversion and investigated its potential to improve contrasts between layers. Like the BPI for poststack case, the sparse layer constraint, rather than the sparse spike constraint, is used to construct the model space as a wedge dictionary. All the elements of the dictionary are bed reflectivities, which include solutions consisting of thin beds as well. With this dictionary, we use an <mml:math display="inline"><mml:msub><mml:mi>L</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:math>L1 norm optimization framework to derive three reflectivities, namely, <mml:math display="inline"><mml:msub><mml:mi>R</mml:mi><mml:mi>p</mml:mi></mml:msub></mml:math>Rp, <mml:math display="inline"><mml:msub><mml:mi>R</mml:mi><mml:mi>s</mml:mi></mml:msub></mml:math>Rs, and <mml:math display="inline"><mml:msub><mml:mi>R</mml:mi><mml:mi>ρ</mml:mi></mml:msub></mml:math>Rρ. Although BPI does not require a starting model, high-resolution absolute velocities (<mml:math display="inline"><mml:mrow><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="normal"><mml:mrow><mml:mi mathvariant="normal">P</mml:mi></mml:mrow></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:math>VP, <mml:math display="inline"><mml:mrow><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="normal"><mml:mrow><mml:mi mathvariant="normal">S</mml:mi></mml:mrow></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:math>VS) and density (<mml:math display="inline"><mml:mi>ρ</mml:mi></mml:math>ρ) can be obtained by incorporating initial models in the BPI derived reflectivities. Tests on synthetic and field data show that the BPI algorithm can indeed detect and enhance layer boundaries by effectively removing the wavelet interference.