Comparison of the projection onto convex sets and iterative hard thresholding methods for seismic data interpolation and denoising*
Benfeng Wang 1 3 Chenglong Gu 21 School of Ocean and Earth Science, Institute for Advanced Study, Tongji University, Shanghai, 200092, China.
2 Exploration and Development Research Institute of Tarim Oilfield, Xinjiang Korla, 841000, China.
3 Corresponding author. Email: wbf1232007@126.com
Exploration Geophysics 49(6) 825-832 https://doi.org/10.1071/EG17046
Submitted: 11 March 2017 Accepted: 18 November 2017 Published: 31 January 2018
Abstract
Because of the environment limitations, irregularity appears in the observed seismic data. In addition, the observed seismic data contains random noise from the acquisition equipment and surrounding environment, which affects the performances of multi-channel techniques, such as surface related multiple elimination (SRME) and amplitude variation with offset (AVO) analysis. The projection onto convex sets (POCS) method, known as an efficient interpolation method, is suitable for high signal-to-noise ratio (SNR) situations; however, the existing random noise may affect its final performance. In our previously published paper, the POCS formula was deduced in the view of iterative hard thresholding (IHT) method using a projection operator. In this paper, more physical illustrations about its detailed deduction are provided to show the differences between IHT and POCS in noise-free and noisy situations with easy understanding for readers. Then, performances of the POCS and IHT methods are compared in both noise-free and noisy situations, in terms of seismograms, frequency wavenumber (FK) spectra and single traces. For noise-free data, both the POCS and IHT methods can achieve good interpolation results. For noisy data, the POCS method is unsuitable because of the observed noisy data insertion, while the IHT performance is satisfactory because it uses a thresholding operator to eliminate random noise. Numerical examples on noise-free datasets demonstrate the validities of the POCS and IHT methods for interpolation. Tests on noisy data contaminated with additive white Gaussian noise prove the ability and superiority of the IHT method with anti-noise property compared with the POCS method.
Key words: interpolation, iterative hard thresholding (IHT), projection onto convex sets (POCS), random noise.
References
Abma, R., and Kabir, N., 2006, 3D interpolation of irregular data with a POCS algorithm: Geophysics, 71, E91–E97| 3D interpolation of irregular data with a POCS algorithm:Crossref | GoogleScholarGoogle Scholar |
Beckouche, S., and Ma, J., 2014, Simultaneous dictionary learning and denoising for seismic data: Geophysics, 79, A27–A31
| Simultaneous dictionary learning and denoising for seismic data:Crossref | GoogleScholarGoogle Scholar |
Blumensath, T., and Davies, M. E., 2008, Iterative thresholding for sparse approximations: The Journal of Fourier Analysis and Applications, 14, 629–654
| Iterative thresholding for sparse approximations:Crossref | GoogleScholarGoogle Scholar |
Blumensath, T., and Davies, M. E., 2009, Iterative hard thresholding for compressed sensing: Applied and Computational Harmonic Analysis, 27, 265–274
| Iterative hard thresholding for compressed sensing:Crossref | GoogleScholarGoogle Scholar |
Bregman, L., 1967, The relaxation method of finding the common point of convex sets and its application to the solution of problems in convex programming: USSR Computational Mathematics and Mathematical Physics, 7, 200–217
| The relaxation method of finding the common point of convex sets and its application to the solution of problems in convex programming:Crossref | GoogleScholarGoogle Scholar |
Candès, E., Demanet, L., Donoho, D., and Ying, L., 2006, Fast discrete curvelet transforms: Multiscale Modeling & Simulation, 5, 861–899
| Fast discrete curvelet transforms:Crossref | GoogleScholarGoogle Scholar |
Chen, Y., and Ma, J., 2014, Random noise attenuation by f-x empirical-mode decomposition predictive filtering: Geophysics, 79, V81–V91
| Random noise attenuation by f-x empirical-mode decomposition predictive filtering:Crossref | GoogleScholarGoogle Scholar |
Gao, J., Chen, X., Li, J., Liu, G., and Ma, J., 2010, Irregular seismic data reconstruction based on exponential threshold model of POCS method: Applied Geophysics, 7, 229–238
| Irregular seismic data reconstruction based on exponential threshold model of POCS method:Crossref | GoogleScholarGoogle Scholar |
Gao, J., Stanton, A., Naghizadeh, M., Sacchi, M. D., and Chen, X., 2013, Convergence improvement and noise attenuation considerations for beyond alias projection onto convex sets reconstruction: Geophysical Prospecting, 61, 138–151
| Convergence improvement and noise attenuation considerations for beyond alias projection onto convex sets reconstruction:Crossref | GoogleScholarGoogle Scholar |
Hennenfent, G., and Herrmann, F. J., 2008, Simply denoise: wavefield reconstruction via jittered undersampling: Geophysics, 73, V19–V28
| Simply denoise: wavefield reconstruction via jittered undersampling:Crossref | GoogleScholarGoogle Scholar |
Herrmann, F. J., and Hennenfent, G., 2008, Non-parametric seismic data recovery with curvelet frames: Geophysical Journal International, 173, 233–248
| Non-parametric seismic data recovery with curvelet frames:Crossref | GoogleScholarGoogle Scholar |
Liu, G., and Chen, X., 2013, Noncausal f–x–y regularized nonstationary prediction filtering for random noise attenuation on 3D seismic data: Journal of Applied Geophysics, 93, 60–66
| Noncausal f–x–y regularized nonstationary prediction filtering for random noise attenuation on 3D seismic data:Crossref | GoogleScholarGoogle Scholar |
Naghizadeh, M., 2012, Seismic data interpolation and denoising in the frequency-wavenumber domain: Geophysics, 77, V71–V80
| Seismic data interpolation and denoising in the frequency-wavenumber domain:Crossref | GoogleScholarGoogle Scholar |
Oropeza, V., and Sacchi, M., 2011, Simultaneous seismic data denoising and reconstruction via multichannel singular spectrum analysis: Geophysics, 76, V25–V32
| Simultaneous seismic data denoising and reconstruction via multichannel singular spectrum analysis:Crossref | GoogleScholarGoogle Scholar |
Stark, H., and Oskoui, P., 1989, High-resolution image recovery from image-plane arrays, using convex projections: Journal of the Optical Society of America A, 6, 1715–1726
| High-resolution image recovery from image-plane arrays, using convex projections:Crossref | GoogleScholarGoogle Scholar | 1:STN:280:DyaK3c%2Fmtlyquw%3D%3D&md5=424d2b827427c853d57792acf9aefb47CAS |
Wang, B., 2016, An efficient POCS interpolation method in the frequency-space domain: IEEE Geoscience and Remote Sensing Letters, 13, 1384–1387
| An efficient POCS interpolation method in the frequency-space domain:Crossref | GoogleScholarGoogle Scholar |
Wang, B., Wu, R., Geng, Y., and Chen, X., 2014, Dreamlet-based interpolation using POCS method: Journal of Applied Geophysics, 109, 256–265
| Dreamlet-based interpolation using POCS method:Crossref | GoogleScholarGoogle Scholar |
Wang, B., Wu, R.-S., Chen, X., and Li, J., 2015, Simultaneous seismic data interpolation and denoising with a new adaptive method based on dreamlet transform: Geophysical Journal International, 201, 1182–1194
| Simultaneous seismic data interpolation and denoising with a new adaptive method based on dreamlet transform:Crossref | GoogleScholarGoogle Scholar |
Wang, B., Chen, X., Lu, W., Guo, Z., and Guo, K., 2016, Comparison of the POCS and IHT methods for seismic data interpolation and denoising: SEG International Exposition and 86th Annual Meeting, 4845–4849.
Yang, P., Gao, J., and Chen, W., 2012, Curvelet-based POCS interpolation of nonuniformly sampled seismic records: Journal of Applied Geophysics, 79, 90–99
| Curvelet-based POCS interpolation of nonuniformly sampled seismic records:Crossref | GoogleScholarGoogle Scholar |
Zhang, H., Chen, X., and Li, H., 2015, 3D seismic data reconstruction based on complex-valued curvelet transform in frequency domain: Journal of Applied Geophysics, 113, 64–73
| 3D seismic data reconstruction based on complex-valued curvelet transform in frequency domain:Crossref | GoogleScholarGoogle Scholar |