Simultaneous seismic interpolation and denoising based on sparse inversion with a 3D low redundancy curvelet transform
Jingjie Cao 1 3 Jingtao Zhao 21 Hebei GEO University, Shijiazhuang, Hebei 050031, China.
2 State Key Laboratory of Coal Resources and Safe Mining, China University of Mining and Technology (Beijing), Beijing 100083, China.
3 Corresponding author. Email: cao18601861@163.com
Exploration Geophysics 48(4) 422-429 https://doi.org/10.1071/EG15097
Submitted: 8 September 2015 Accepted: 28 April 2016 Published: 9 June 2016
Abstract
The simultaneous seismic interpolation and denoising problem can be solved as a sparse inversion problem by using the sparseness of seismic data in a transformed domain as the a priori information, where the properties of the sparse transform will significantly influence the numerical results and computational efficiency. Curvelet transform has nearly optimal sparse expression for seismic data, thus seismic signal processing based on this transform tends to result in preferable results. However, the redundancy of this transform can be 24–32 for three dimensional data which is not computationally cost efficient. This paper introduces a low redundancy curvelet transform to simultaneously interpolate and denoise. The redundancy of the proposed transform can be reduced to 10 for three dimensional data, and this property will improve the computational efficiency of the curvelet transform-based signal processing. The iterative soft thresholding method was chosen to solve the sparse inversion problem. Some practical principles on how to choose the regularisation parameters are discussed, due to the crucial nature of the regularisation parameter for simultaneous interpolation and denoising. Numerical experiments on synthetic and field data demonstrate that the low redundancy transform can provide reliable results while at the same time improving the computational efficiency.
Key words: curvelet transform, inverse problems, sparse inversion, wavefield interpolation.
References
Berkhout, A. J., and Verschuur, D. J., 1997, Estimation of multiple scattering by iterative inversion, part I: theoretical considerations: Geophysics, 62, 1586–1595| Estimation of multiple scattering by iterative inversion, part I: theoretical considerations:Crossref | GoogleScholarGoogle Scholar |
Candes, E., and Donoho, D., 2000, Curvelets – a surprisingly effective nonadaptive representation for objects with edges: Vanderbilt University Press.
Candes, E., and Donoho, D., 2004, New tight frames of curvelets and optimal representation of objects with piecewise C2 singularities: Communications on Pure and Applied Mathematics, 57, 219–266
| New tight frames of curvelets and optimal representation of objects with piecewise C2 singularities:Crossref | GoogleScholarGoogle Scholar |
Candes, E. J., Demanet, L., Donoho, D., and Ying, L., 2006, Fast discrete curvelet transforms: SIAM Journal on Multiscale Modeling and Simulation, 5, 861–899
| Fast discrete curvelet transforms:Crossref | GoogleScholarGoogle Scholar |
Cao, J. J., and Wang, B. F., 2015, An improved projection onto convex sets method for simultaneous interpolation and denoising: Chinese Journal of Geophysics, 58, 2935–2947
Cao, J. J., and Zhao, J. T., 2015, 3D seismic interpolation with a low redundancy, fast curvelet transform: Journal of Seismic Exploration, 24, 121–134
Cao, J. J., Wang, Y. F., Zhao, J. T., and Yang, C. C., 2011, A review on restoration of seismic wavefields based on regularization and compressive sensing: Inverse Problems in Science and Engineering, 19, 679–704
| A review on restoration of seismic wavefields based on regularization and compressive sensing:Crossref | GoogleScholarGoogle Scholar |
Cao, J. J., Wang, Y. F., and Yang, C. C., 2012, Seismic data restoration based on compressive sensing using the regularization and zero-norm sparse optimization: Chinese Journal of Geophysics, 55, 239–251
| Seismic data restoration based on compressive sensing using the regularization and zero-norm sparse optimization:Crossref | GoogleScholarGoogle Scholar |
Cao, J. J., Wang, Y. F., and Wang, B. F., 2015a, Accelerating seismic interpolation with a gradient projection method based on tight frame property of curvelet: Exploration Geophysics, 46, 253–260
| Accelerating seismic interpolation with a gradient projection method based on tight frame property of curvelet:Crossref | GoogleScholarGoogle Scholar |
Cao, J. J., Zhao, J. T., and Hu, Z. Y., 2015b, 3D seismic denoising based on a low-redundancy curvelet transform: Journal of Geophysics and Engineering, 12, 566–576
| 3D seismic denoising based on a low-redundancy curvelet transform:Crossref | GoogleScholarGoogle Scholar |
Chauris, H., and Nguyen, T., 2008, Seismic demigration/migration in the curvelet domain: Geophysics, 73, S35–S46
| Seismic demigration/migration in the curvelet domain:Crossref | GoogleScholarGoogle Scholar |
Chen, S., Donoho, D., and Saunders, M., 1998, Atomic decomposition by basis pursuit: SIAM Journal on Scientific Computing, 20, 33–61
| Atomic decomposition by basis pursuit:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DyaK1cXkvV2ltrw%3D&md5=75ca9aad088f4165f7aad8442aa55227CAS |
Daubechies, I., Defrise, M., and Mol, C. D., 2004, An iterative thresholding algorithm for linear inverse problems with a sparsity constraint: Communications on Pure and Applied Mathematics, 57, 1413–1457
| An iterative thresholding algorithm for linear inverse problems with a sparsity constraint:Crossref | GoogleScholarGoogle Scholar |
Fomel, S., and Liu, Y., 2010, Seislet transform and seislet frame: Geophysics, 75, V25–V38
| Seislet transform and seislet frame:Crossref | GoogleScholarGoogle Scholar |
Gao, J. J., Stanton, A., Naghizadeh, M., Sacchi, M. D., and Chen, X. H., 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 |
Geng, Y., Wu, R. S., and Gao, J. H., 2012, Dreamlet compression of seismic data: Chinese Journal of Geophysics, 55, 2705–2715
Hennenfent, G., and Herrmann, F., 2006, Seismic denoising with non-uniformly sampled curvelets: Computing in Science & Engineering, 8, 16–25
| Seismic denoising with non-uniformly sampled curvelets:Crossref | GoogleScholarGoogle Scholar |
Herrmann, F., 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 |
Herrmann, F., Moghaddam, P., and Stolk, C., 2008a, Sparsity- and continuity-promoting seismic image recovery with curvelet frames: Applied and Computational Harmonic Analysis, 24, 150–173
| Sparsity- and continuity-promoting seismic image recovery with curvelet frames:Crossref | GoogleScholarGoogle Scholar |
Herrmann, F., Wang, D., Hennenfent, G., and Moghaddam, P., 2008b, Curvelet-based seismic data processing: a multiscale and nonlinear approach: Geophysics, 73, A1–A5
| Curvelet-based seismic data processing: a multiscale and nonlinear approach:Crossref | GoogleScholarGoogle Scholar |
Kreimer, N., and Sacchi, M. D., 2012, A tensor higher-order singular value decomposition for prestack seismic data reduction and interpolation: Geophysics, 77, V113–V122
| A tensor higher-order singular value decomposition for prestack seismic data reduction and interpolation:Crossref | GoogleScholarGoogle Scholar |
Kumar, V., and Herrmann, F., 2008, Deconvolution with curvelet-domain sparsity: SEG Technical Program, Expanded Abstracts, 1996–2000.
Kumar, V., Oueity, J., Clowes, R. M., and Herrmann, F., 2011, Enhancing crustal reflection data through curvelet denoising: Tectonophysics, 508, 106–116
| Enhancing crustal reflection data through curvelet denoising:Crossref | GoogleScholarGoogle Scholar |
Lin, T., and Herrmann, F., 2013, Robust estimation of primaries by sparse inversion via one-norm minimization: Geophysics, 78, R133–R150
| Robust estimation of primaries by sparse inversion via one-norm minimization:Crossref | GoogleScholarGoogle Scholar |
Liu, B., and Sacchi, M. D., 2004, Minimum weighted norm interpolation of seismic records: Geophysics, 69, 1560–1568
| Minimum weighted norm interpolation of seismic records:Crossref | GoogleScholarGoogle Scholar |
Ma, J. W., 2013, Three-dimensional irregular seismic data reconstruction via low-rank matrix completion: Geophysics, 78, V181–V192
| Three-dimensional irregular seismic data reconstruction via low-rank matrix completion:Crossref | GoogleScholarGoogle Scholar |
Mansour, H., Herrmann, F., and Yilmaz, O., 2013, Improved wavefield reconstruction from randomized sampling via weighted one-norm minimization: Geophysics, 78, V193–V206
| Improved wavefield reconstruction from randomized sampling via weighted one-norm minimization:Crossref | GoogleScholarGoogle Scholar |
Montefusco, L. B., and Papi, S., 2003, A parameter selection method for wavelet shrinkage denoising: BIT Numerical Mathematics, 43, 611–626
| A parameter selection method for wavelet shrinkage denoising:Crossref | GoogleScholarGoogle Scholar |
Neelamani, R., Baumstein, A., Gillard, D., Hadidi, M., and Soroka, W., 2008, Coherent and random noise attenuation using the curvelet transform: The Leading Edge, 27, 240–248
| Coherent and random noise attenuation using the curvelet transform:Crossref | GoogleScholarGoogle Scholar |
Oropeza, V., and Sacchi, M. D., 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 |
Ronen, J., 1987, Wave-equation trace interpolation: Geophysics, 52, 973–984
| Wave-equation trace interpolation:Crossref | GoogleScholarGoogle Scholar |
Sacchi, M., and Liu, B., 2005, Minimum weighted norm wavefield reconstruction for AVA imaging: Geophysical Prospecting, 53, 787–801
| Minimum weighted norm wavefield reconstruction for AVA imaging:Crossref | GoogleScholarGoogle Scholar |
Spitz, S., 1991, Seismic trace interpolation in the f-x domain: Geophysics, 56, 785–794
| Seismic trace interpolation in the f-x domain:Crossref | GoogleScholarGoogle Scholar |
Stolt, R. H., 2002, Seismic data mapping and reconstruction: Geophysics, 67, 890–908
| Seismic data mapping and reconstruction:Crossref | GoogleScholarGoogle Scholar |
Symes, W. W., 2007, Reverse time migration with optimal checkpointing: Geophysics, 72, SM213–SM221
| Reverse time migration with optimal checkpointing:Crossref | GoogleScholarGoogle Scholar |
Trad, D., Ulrych, T., and Sacchi, M., 2002, Accurate interpolation with high-resolution time-variant Radon transforms: Geophysics, 67, 644–656
| Accurate interpolation with high-resolution time-variant Radon transforms:Crossref | GoogleScholarGoogle Scholar |
van den Berg, E., and Friedlander, M. P., 2009, Probing the Pareto frontier for basis pursuit solutions: SIAM Journal on Scientific Computing, 31, 890–912
| Probing the Pareto frontier for basis pursuit solutions:Crossref | GoogleScholarGoogle Scholar |
Wang, Y. F., Cao, J. J., and Yang, C. C., 2011, Recovery of seismic wavefields based on compressive sensing by an l 1-norm constrained trust region method and the piecewise random subsampling: Geophysical Journal International, 187, 199–213
| Recovery of seismic wavefields based on compressive sensing by an l 1-norm constrained trust region method and the piecewise random subsampling:Crossref | GoogleScholarGoogle Scholar |
Woiselle, A., Starck, J., and Fadili, J., 2011, 3D data denoising and inpainting with the low-redundancy fast curvelet transform: Journal of Mathematical Imaging and Vision, 39, 121–139
| 3D data denoising and inpainting with the low-redundancy fast curvelet transform:Crossref | GoogleScholarGoogle Scholar |
Xiao, T. Y., Yu, S. G., and Wang, Y. F., 2003, Numerical methods of inverse problems [in Chinese]: Science Press (Beijing).
Xu, S., Zhang, Y., and Lambaré, G., 2010, Antileakage Fourier transform for seismic data regularization in higher dimensions: Geophysics, 75, WB113–WB120
| Antileakage Fourier transform for seismic data regularization in higher dimensions:Crossref | GoogleScholarGoogle Scholar |
Zhang, R. F., and Ulrych, T., 2003, Physical wavelet frame denoising: Geophysics, 68, 225–231
| Physical wavelet frame denoising:Crossref | GoogleScholarGoogle Scholar |
Zwartjes, P., 2005, Fourier reconstruction with sparse inversion: Ph.D. thesis, Delft University of Technology.