Inversion-Based Time-Domain Inverse $Q$ -Filtering for Seismic Resolution Enhancement

Seismic wave will suffer from resolution loss during its propagation in subsurface media, which is caused by the earth filtering effect. Inverse Q-filtering (IQF) is an essential scheme for eliminating earth filtering effect and enhancing temporal resolution. In this letter, IQF is achieved by solving a time-domain inverse problem whose corresponding forward model is established based on the idea of seismic migration, i.e., extrapolating the unfiltered record (the detected signal if there is no earth filtering effect) from surface downward to each time–depth point and then adopting imaging condition to obtain the attenuated record. We express this migration-based forward model as a matrix-vector equation and construct the objective function of inversion-based time-domain IQF (ITIQF). Instead of adopting frequency-domain incomplete amplitude compensation operator for the purpose of antinoise, we propose three transform-domain sparse regularization strategies, including wavelet-based constraint, curvelet-based constraint, and double constraint (a combination of wavelet- and curvelet-based constraints). Solvers of these constrained ITIQF methods are provided according to convex optimization theory. We demonstrate the superior performance of ITIQF scheme via both synthetic and field data examples.

[1]  Einar Kjartansson,et al.  Constant Q-wave propagation and attenuation , 1979 .

[2]  Jianguo Zhao,et al.  Compressional wave dispersion due to rock matrix stiffening by clay squirt flow , 2016 .

[3]  Sanyi Yuan,et al.  Sparse Bayesian Learning-Based Time-Variant Deconvolution , 2017, IEEE Transactions on Geoscience and Remote Sensing.

[4]  Tom Goldstein,et al.  The Split Bregman Method for L1-Regularized Problems , 2009, SIAM J. Imaging Sci..

[5]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[6]  A. C. Grant Multichannel seismic reflection profiles of the continental crust beneath the Newfoundland Ridge , 1977, Nature.

[7]  Tadeusz J. Ulrych,et al.  Seismic absorption compensation: A least squares inverse scheme , 2007 .

[8]  Yanghua Wang,et al.  A stable and efficient approach of inverse Q filtering , 2002 .

[9]  Sanyi Yuan,et al.  Stable inversion-based multitrace deabsorption method for spatial continuity preservation and weak signal compensation , 2016 .

[10]  Ali Gholami Semi-blind nonstationary deconvolution: Joint reflectivity and Q estimation , 2015 .

[11]  D. Donoho,et al.  Sparse MRI: The application of compressed sensing for rapid MR imaging , 2007, Magnetic resonance in medicine.

[12]  Yanghua Wang,et al.  Inverse Q-filter for seismic resolution enhancement , 2006 .

[13]  Abolghasem Kamkar Rouhani,et al.  Seismic resolution enhancement using scale transform in the time-frequency domain , 2018, GEOPHYSICS.

[14]  Xiaotao Wen,et al.  Q estimation of seismic data using the generalized S-transform , 2016 .

[15]  Mauricio D. Sacchi,et al.  A Fast and Automatic Sparse Deconvolution in the Presence of Outliers , 2012, IEEE Transactions on Geoscience and Remote Sensing.

[16]  Libin Liu,et al.  Data-Driven Attenuation Compensation Via a Shaping Regularization Scheme , 2018, IEEE Geoscience and Remote Sensing Letters.

[17]  L. Denham,et al.  Seismic interpretation , 1984, Proceedings of the IEEE.