Front. Earth Sci. Frontiers in Earth Science Front. Earth Sci. 2296-6463 Frontiers Media S.A. 1103522 10.3389/feart.2023.1103522 Earth Science Original Research Reconstruction of seismic data based on SFISTA and curvelet transform Tian and Qin 10.3389/feart.2023.1103522 Tian Lin 1 2 * Qin Si 1 1 School of Electronic Engineering, Yili Normal University, Yining, China 2 Key Laboratory of Signal Detection and Control Technology, Yili Normal University, Yining, China

Edited by: Peng Zhenming, University of Electronic Science and Technology of China, China

Reviewed by: Shu Li, Jishou University, China

Bibo Yue, Southwest Petroleum University, China

Gulan Zhang, Southwest Petroleum University, China

 *Correspondence: Lin Tian, 1610356358@qq.com

This article was submitted to Environmental Informatics and Remote Sensing, a section of the journal Frontiers in Earth Science

 23 02 2023 2023 11 1103522 20 11 2022 03 02 2023 Copyright © 2023 Tian and Qin. 2023 Tian and Qin

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

In seismic data processing, the reconstruction and interpolation of missing traces are essential tasks. To overcome the limitations of irregularly sampled seismic data, this paper proposes a seismic data interpolation method combining the smoothing fast iterative soft threshold algorithm (SFISTA) and the curvelet transform; this method uses the curvelet domain as the sparse domain. For comparison, the contourlet transform is used for different sparse domains, and the fast iterative shrinkage-thresholding algorithm (FISTA) is used for different optimization algorithms. Numerical modeling and real data show that the SFISTA method in the curvelet domain can give better reconstruction effects and higher reconstruction accuracy than those in the contourlet domain with the FISTA method; in addition, the SFISTA method in the curvelet domain can be used to reconstruct the missing traces of three-dimensional seismic data.

 seismic data reconstruction sparse constraint curvelet transform contourlet transform SFISTA 1 Introduction

In complex environments, seismic data reconstruction has great significance as a recovery technique. Under external disturbance, irregularly sampled seismic data will affect further geological data processing such as migration imaging and data interpretation. In order to obtain high-quality seismic data, interpolation reconstruction is needed to approximate the original data. In recent years, under the compressive sensing theory, seismic data reconstruction methods based on sparse constraints have become more and more popular. It mainly consists of the sparse transform, measurement matrix, and reconstruction algorithm. The sparse transforms that are often used include the Fourier transform (Zhang et al., 2013; Ciabarri et al., 2014), curvelet transform (Hennenfent et al., 2010; Liu et al., 2015; Zhang et al., 2017; Zhang et al., 2019; Tian and Qin, 2022), contourlet transform (Eslami and Radha, 2006), and seislet transform (Liu W et al., 2016). Because the curvelet transform undergoes multi-scale and multi-direction analysis and can perform the optimal local decomposition of seismic data (Yang et al., 2017), the curvelet transform is employed in this paper as a sparse transform, and the contourlet transform is also used in this paper for comparison analysis.

A classical sparse recovery problem usually requires minimizing the L 0 norm, which is NP-hard. The L 1 norm is the approximate of the L 0 norm, which is a convex function, and can be solved by the convex optimization algorithms or tools; so, the L 0 norm is replaced by the L 1 norm for simplicity and effectiveness. The iterative soft threshold algorithm (ISTA) shows great advantages (Daubechies et al., 2004; Mohsin et al., 2015) for a convex optimization algorithm. Gradually, a faster ISTA algorithm (FISTA) has been developed. FISTA is more suitable the synthesis approach to sparse recovery. And for the analysis approach to sparse recovery, Tan et al. (2014) proposed a monotone version of the fast iterative shrinkage-thresholding algorithm, which is the smoothing fast iterative soft threshold algorithm (SFISTA). In this paper, we introduce the SFISTA method-based curvelet transform to the seismic data interpolation reconstruction problem (Tan et al., 2014; Liu Y et al., 2016; Pokala and Seelamantula, 2020; Shen et al., 2020). The theory is given in section 2, and experimental results are given in section 3.

 2 Theory

Assuming that the observed seismic data is y and the downsampling matrix is U , the irregular missing seismic data can be modeled as follows: y = U x + n , where n is a randomly generated noise and x is the original seismic data.

The interpolation problem can be expressed as follows: min x 1 2 y − U x 2 2 + λ Ψ x 1 , where ∙ 2 is the L 2 norm, ∙ 1 is the L 1 norm, λ is the regularization parameter, and Ψ is the analysis operator. Moreover, Ψ * Ψ = I , where Ψ * denotes the adjoint of the operator Ψ and I is the identity matrix. The sparse transform can be expressed as x = Ψ α , where Ψ is the transform operator and α is the sparse coefficient. Let us denote the following equation: f x = 1 2 U x − y 2 2 , where f x is the smooth part. g Ψ x = λ Ψ x 1 , where g Ψ x is the non-smooth part that needs to be replaced by the Moreau envelope with an approximately smooth form g μ Ψ x , where μ is the smooth approximate parameter. ∇ g μ Ψ x = 1 μ Ψ * Ψ x − T λ μ Ψ x , where ∇ is the gradient and T λ μ is the soft threshold operator; the threshold of T λ μ is the multiplication of λ and μ .

In this paper, the process of the smoothing fast iterative soft threshold algorithm (Shen et al., 2020) for optimization is expressed as follows: ∇ f x ^ k = U T U x ^ k − y , where k is the number of iterations, and x ^ 0 = x 0   a n d   t 0 = 1 . ∇ g μ Ψ x k = 1 μ Ψ * Ψ x k − T λ μ Ψ x k .

The core iterations of the SFISTA are as follows: x k + 1 = x ^ k − 1 L ∇ f x ^ k + ∇ g μ Ψ x k , where L is the Lipschitz constant. At 0 < γ ≤ 1 , L = 1 / γ . Substituting (6) and (7) in (8), the aforementioned equation is equivalent to the following equation: x k + 1 = x ^ k − γ μ x k − Ψ * T λ μ Ψ x k + γ U T y − U x ^ k , where γ is the iteration step size. t k + 1 is expressed as follows: t k + 1 = 1 + 1 + 4 t k 2 2 .

Then, we derive the following equation: x ^ k + 1 = x k + 1 + t k − 1 t k + 1 x k + 1 − x k .

The specific algorithm steps are as follows:

SFISTA for seismic data interpolation

Parameters: λ = 10 − 3 , μ = 10 − 3 / λ = 1 , γ = 1 / 1 + 1 / μ

Initialization: x 0 , x ^ 0 , t 0 = 1

When not convergent, the following equations are used to calculate:

x k + 1 = x ^ k − γ μ x k − Ψ * T λ μ Ψ x k + γ U T y − U x ^ k

t k + 1 = 1 + 1 + 4 t k 2 2

x ^ k + 1 = x k + 1 + t k − 1 t k + 1 x k + 1 − x k

R L N E = ∥ x ^ − x ∥ 2 ∥ x ∥ 2

end if R L N E < 10 − 6 or the maximum number of iterations is reached.

end

Output: x

FISTA for seismic data interpolation

Parameters: λ = 10 − 3 , γ = 1

Initialization: α 0 , α ^ 0 , t 0 = 1

When not convergent, the following equations are used to calculate:

α k + 1 = T γ λ γ Ψ U T y − U Ψ * α k + α ^ k

t k + 1 = 1 + 1 + 4 t k 2 2

α ^ k + 1 = α k + 1 + t k − 1 t k + 1 α k + 1 − α k

R L N E = ∥ α ^ − α ∥ 2 ∥ α ∥ 2

end if R L N E < 10 − 6 or the maximum number of iterations is reached.

end

Output: α

3 Examples

In this section, we conduct numerical experiments with different seismic data to demonstrate the reconstruction performance of the SFISTA method in the curvelet domain. The algorithm performance is evaluated by interpolation results, the average amplitude spectrum, single-channel interpolation effect, reconstruction error, and signal-to-noise ratio. Numerical experiments are used to test the method. At last, the paper continues to test the 3D interpolation effect of the interpolation method of the proposed method. The experiments are conducted on a Millet computer running on the Windows 10 operating system and Inter Core m3-6Y30.

3.1 Synthetic examples 3.1.1 Seismic data interpolation in the contourlet domain

The regularization parameters of the FISTA and SFISTA methods are set to 10 − 3 . The step size of the FISTA is 1, and the step size of the SFISTA is set to 0.5 (Liu Y et al., 2016); the number of iterations is set as 500. The two-dimensional data tests have the same parameters.

Part data on the Marmousi2 model (Martin et al., 2006) are selected as test data; the data are shown in Figure 1A. Figure 1B shows the irregular data with 50% of the traces removed randomly, and this part’s records are padded with zeros. The sparse domain is the contourlet domain. Figures 1C,D are interpolation results of Figure 1B using FISTA and SFISTA methods based on the contourlet transform. Figures 1E,F show the reconstruction errors that correspond to the FISTA and SFISTA methods, in which the reconstruction residuals of the SFISTA method have smaller amplitude ranges than those of the FISTA method. It is obvious that the reconstruction method of SFISTA based on the contourlet transform is more effective.

Reconstruction results of synthetic data in the contourlet domain. (A) Original data; (B) 50% of randomly missing data; (C) interpolated data obtained by the FISTA method; (D) interpolated data obtained by the SFISTA method; (E) reconstructed errors obtained by the FISTA method; (F) reconstructed errors obtained by the SFISTA method.

The performance of the proposed method in seismic data interpolation could be evaluated using different qualitative and quantitative analyzing tools. The reconstruction error is the general quantitative evaluation tool used in seismic data interpolation; the reconstruction error is defined to be ∑ x − x ^ , where x denotes the original data and x ^ denotes the reconstructed seismic data. If the reconstruction error is smaller, the reconstructed seismic data will be closer to the original data. SNR is the signal-to-noise ratio, which is defined as S N R = 101 o g 10 x 2 2 x − x ^ 2 2 . A higher SNR value means that the data have better reconstruction quality. The reconstruction error and signal-to-noise ratio are illustrated in Table 1. Table 1 shows that the SFISTA method has a smaller reconstruction error value and higher SNR. Compared to the two methods, the SFISTA method based on the contourlet transform has a lower reconstruction error and higher SNR, and therefore, the SFISTA method shows better performance.

Reconstruction error and SNR of different methods.

FISTA SFISTA
Reconstruction error 784.34 450.45
SNR 10.12 12.41
3.1.2 Seismic data interpolation in the curvelet domain

The part of the Marmousi2 model is shown in Figure 2A. Figure 2B shows the irregular data with 50% of the traces missing randomly. Figures 2C, D show interpolation effects using FISTA and SFISTA methods based on the curvelet transform. Figures 2E, F show the reconstruction errors that correspond to the FISTA and SFISTA methods, in which the reconstruction residuals of the SFISTA method have a smaller magnitude than those of the FISTA method. It is obvious that the reconstruction method of SFISTA in the curvelet domain is more effective.

Reconstruction results of synthetic data in the curvelet domain. (A) Original data; (B) 50% of randomly missing data; (C) interpolated data obtained by the FISTA method; (D) interpolated data obtained by the SFISTA method; (E) reconstructed errors obtained by the FISTA method; (F) reconstructed errors obtained by the SFISTA method.

Comparing Figures 1, 2, the SFISTA method in the curvelet domain has lower reconstruction errors when using the qualitative analyzing method. Comparing Tables 1, 2, the SFISTA method in the curvelet domain has lower reconstruction errors and higher SNR when using the quantitative analyzing method. The SFISTA method shows better performance when using the curvelet domain as the sparse domain.

Reconstruction error and SNR of different methods.

FISTA SFISTA
Reconstruction error 223.16 193.62
SNR 17.75 20.67
3.1.3 Coral reef model tests in the contourlet domain

The selected coral reef synthetic data are illustrated in Figure 3A. Figure 3B shows corrupted data with 50% of traces missing randomly. Figures 3C, D show interpolation effects using FISTA and SFISTA methods in the contourlet domain. Figures 3E, F show the reconstructed errors between the reconstruction and original data; Figure 3F shows a smaller reconstructed error. Table 3 shows the reconstruction error and SNR, which demonstrates the validity of the SFISTA method in the contourlet domain.

Reconstruction results of the coral reef model in the contourlet domain. (A) Original data; (B) 50% of randomly missing data; (C) interpolated data obtained by the FISTA method; (D) interpolated data obtained by the SFISTA method; (E) reconstructed errors obtained by the FISTA method; (F) reconstructed errors obtained by the SFISTA method.

Reconstruction error and SNR of different methods.

FISTA SFISTA
Reconstruction error 1695.09 407.17
SNR 9.11 25.07
3.1.4 Coral reef model tests in the curvelet domain

The selected Coral reef synthetic data are tested in the curvelet domain; the original data are shown in Figure 4A. Figure 4B shows corrupted data with 50% of traces missing randomly. Figures 4C, D show interpolation results using FISTA and SFISTA methods in the curvelet domain. Figures 4E, F show the reconstructed errors between the reconstruction and original data; Figure 4F shows a smaller reconstructed error. Table 4 shows the reconstruction error and SNR, which demonstrates the validity of the proposed method.

Reconstruction results of the coral reef model in the curvelet domain. (A) Original data; (B) 50% of randomly missing data; (C) interpolated data obtained by the FISTA method; (D) interpolated data obtained by the SFISTA method; (E) reconstructed errors obtained by the FISTA method; (F) reconstructed errors obtained by the SFISTA method.

Reconstruction error and SNR of different methods.

FISTA SFISTA
Reconstruction error 902.74 323.55
SNR 15.45 27.72

Comparing the detailed values in Tables 3, 4, the reconstruction performance of the SFISTA method in the curvelet domain is better than that in the contourlet domain.

 3.2 Field examples 3.2.1 Seismic data interpolation in the contourlet domain

The field data are illustrated in Figure 5A. Figure 5B shows corrupted data with 50% of traces missing randomly. Figures 5C, D show interpolation effects using FISTA and SFISTA methods in the contourlet domain. By comparing the interpolation results of the ellipse regions, we learn that the reconstruction results of the SFISTA method in the contourlet domain are similar to the original data. Figures 5E, F show the reconstructed errors between the reconstruction and original data; Figure 5F shows a smaller reconstructed error. Table 5 shows the reconstruction error and SNR, which demonstrates the validity of the SFISTA method in the contourlet domain.

Reconstruction results of field data in the contourlet domain. (A) Original data; (B) 50% of randomly missing data; (C) interpolated data obtained by the FISTA method; (D) interpolated data obtained by the SFISTA method; (E) reconstructed errors obtained by the FISTA method; (F) reconstructed errors obtained by the SFISTA method.

Reconstruction error and SNR of different methods.

FISTA SFISTA
Reconstruction error 3904.69 3633.50
SNR 6.43 7.16
3.2.2 Seismic data interpolation in the curvelet domain

The field data are illustrated in Figure 6A. Figure 6B shows incomplete seismic data with 50% of traces missing randomly. Figures 6C, D show interpolation effects using FISTA and SFISTA methods in the curvelet domain. By comparing the interpolation results of the ellipse regions, we learn that the reconstruction data on the SFISTA method in the curvelet domain are similar to the original data. Figures 6E, F show the reconstructed errors between the reconstruction and original data; Figure 6F shows a smaller reconstructed error. Table 6 shows the reconstruction error and SNR results of the two algorithms in the curvelet domain. Comparing the two methods, SFISTA performs better with a lower error and higher SNR.

Reconstruction results of field data in the curvelet domain. (A) Original data; (B) 50% of randomly missing data; (C) interpolated data obtained by the FISTA method; (D) interpolated data obtained by the SFISTA method; (E) reconstructed errors obtained by the FISTA method; (F) reconstructed errors obtained by the SFISTA method.

Reconstruction error and SNR of different methods.

FISTA SFISTA
Reconstruction error 2968.24 2220.78
SNR 9.18 12.91

Figure 7A plots the interpolation single trace of the missing trace with zero values. In Figures 7, the black, pink, and blue lines correspond to the original data and those obtained by FISTA and SFISTA methods, respectively. Figure 7B shows the zoomed traces from the transparent gray window of Figure 7A. A detailed comparison reveals that the blue line is similar to the black line. We can, therefore, affirm that the proposed SFISTA method can better restore the significant features of the useful signal than the FISTA method.

(A) Comparison of a missing single-trace amplitude interpolated with FISTA and SFISTA methods; (B) zoomed traces from the transparent gray window of (A).

Comparing the detailed values in Tables 5, 6, the SFISTA method in the curvelet domain is better than that in the contourlet domain.

The performance of these methods could be investigated more by comparing the average amplitude spectrum (Mafakheri et al., 2022). The average amplitude spectrum presents the original data and those obtained by FISTA and SFISTA methods, as shown in Figure 8 by the black, pink, and blue lines, respectively. The SFISTA method in the curvelet domain gives a closer spectrum to that of the original data, particularly in the range of 5–30 Hz. In this case, our method has better performance.

Average amplitude spectrum from original data, interpolated with FISTA and SFISTA methods.

 3.2.3 Three-dimensional seismic data tests

The experimental results of three sets of 2D data show that SFISTA based on the curvelet transform shows good performance. Next, this method will be tested for 3D seismic data reconstruction. The 3D data (size: 64 × 64 × 64) are from the software package of MathGeo 2020 (https://gitee.com/sevenysw/MathGeo2020). The 3D discrete curvelet transform method comes from the reference of Ying et al. (2005).

The original data are shown in Figure 9A. Figure 9B shows corrupted data with 51% of the traces removed randomly; the iteration number of the FISTA and SFISTA methods is 1000. Figures 9C, D show the interpolation results using the FISTA and SFISTA methods in the curvelet domain. Figures 9E, F show reconstructed errors obtained by the FISTA and SFISTA methods. The quantitatively recovered reconstruction error and SNR are shown in Table 7; Table 7 illustrates that the SFISTA method is better than the FISTA method.

Reconstruction results of three-dimensional field data in the curvelet domain. (A) Original data; (B) 51% of randomly missing data; (C) interpolated data obtained by the FISTA method; (D) interpolated data obtained by the SFISTA method; (E) reconstructed errors obtained by the FISTA method; (F) reconstructed errors obtained by the SFISTA method.

Reconstruction error and SNR of different methods.

FISTA SFISTA
Reconstruction error 27781.58 18396.99
SNR 7.32 10.79

Figure 10A plots the interpolation single trace from Figures 9C, D. The black, pink, and blue lines represent the original data and those obtained by FISTA and SFISTA methods, respectively. Figure 10B shows the zoomed traces from the transparent gray window of Figure 10A. A detailed comparison reveals that the blue line is similar to the black line. We can, therefore, affirm that the proposed SFISTA method can better restore the significant features of the useful signal than the FISTA approach.

(A) Comparison of a missing single-trace amplitude interpolated with FISTA and SFISTA methods; (B) zoomed traces from the transparent gray window of (A).

The performance of these methods could be investigated more by comparing the average amplitude spectrum. The average amplitude spectrum of the original data and FISTA and SFISTA methods are shown in Figure 11 by the black, pink, and blue lines, respectively. Interestingly, in this case, the reconstruction data obtained by our method are very similar to the original data.

Average amplitude spectrum from original data, interpolated with FISTA and SFISTA methods.

 4 Conclusion

In this paper, we proposed a new seismic data reconstruction method combining the curvelet transform and the SFISTA method. Comparing the results obtained when using the contourlet and curvelet domains as the sparse domain, it can be concluded that the optimization algorithm in the curvelet domain has better performance. In the same sparse domain, comparing the FISTA and SFISTA methods, it can be concluded that SFISTA shows better performance. The seismic data reconstruction effects of SFISTA based on the curvelet transform have been demonstrated by quantitative and qualitative comparisons with several sets of data. The proposed method can be used for 2D and 3D seismic data reconstruction.

 Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

 Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

 Funding

This work was supported by the Yili Normal University Foundation under Grant 22XKZZ22, the Talent Position Project of Yili Normal University under Grant YSXSGG22006, and the National Natural Science Foundation of China under Grant 61761043.

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 References Ciabarri F. Mazzotti A. Stucchi E. Bienati N. (2014). Appraisal problem in the 3D least squares Fourier seismic data reconstruction. Geophys. Prospect. 63 (2), 296–314. 10.1111/1365-2478.12192 Daubechies I. Defrise M. Mol C. D. (2004). An iterative thresholding algorithm for linear inverse problems with a sparsity constraint. Commun. Pure Appl. Math. 57 (11), 1413–1457. 10.1002/cpa.20042 Eslami R. Radha H. (2006). Translation-invariant contourlet transform and its application to image denoising. IEEE Trans. Image Process 15 (11), 3362–3374. 10.1109/tip.2006.881992 Hennenfent G. Fenelon L. Herrmann F. J. (2010). Nonequispaced curvelet transform for seismic data reconstruction: A sparsity-promoting approach. GEOPHYSICS 75 (6), 203–210. 10.1190/1.3494032 Liu W. Cao S. Gan S. Chen Y. Zu S. Jin Z. (2016). One-step slope estimation for dealiased seismic data reconstruction via iterative seislet thresholding. IEEE Geoscience Remote Sens. Lett. 13 (10), 1462–1466. 10.1109/lgrs.2016.2591939 Liu W. Cao S. Li G. He Y. (2015). Reconstruction of seismic data with missing traces based on local random sampling and curvelet transform. J. Appl. Geophys. 115, 129–139. 10.1016/j.jappgeo.2015.02.009 Liu Y. Zhan Z. Cai J.-F. Guo D. Chen Z. Qu X. (2016). Projected iterative soft-thresholding algorithm for tight frames in compressed sensing magnetic resonance imaging. IEEE Trans. Med. Imaging 35 (9), 2130–2140. 10.1109/TMI.2016.2550080 Mafakheri J. Kahoo A. R. Anvari R. Mohammadi M. Radad M. Monfared M. S. (2022). Expand dimensional of seismic data and random noise attenuation using Low-Rank estimation. IEEE J. Sel. Top. Appl. Earth Observations Remote Sens. 15, 2773–2781. 10.1109/jstars.2022.3162763 Martin G. S. Wiley R. Marfurt K. J. (2006). Marmousi2: An elastic upgrade for Marmousi. Lead. Edge 25 (2), 156–166. 10.1190/1.2172306 Mohsin Y. Q. Ongie G. Jacob M. (2015). Iterative shrinkage algorithm for patch-smoothness regularized medical image recovery. IEEE Trans. Med. Imaging 34 (12), 2417–2428. 10.1109/tmi.2015.2398466 Pokala P. Seelamantula C. (2020). “Accelerated weighted ℓ1-minimization for MRI reconstruction under tight frames in complex domain,” in International Conference on Signal Processing and Communications (SPCOM), Bangalore, India, 19-24 July 2020. 10.1109/SPCOM50965.2020.9179611 Shen L. Chu Z. Tan L. Chen D. Ye F. (2020). Improving the sound source identification performance of sparsity constrained deconvolution beamforming utilizing SFISTA. Shock Vib. 2020, 1–9. 10.1155/2020/1482812 Tan Z. Eldar Y. C. Beck A. Nehorai A. (2014). Smoothing and decomposition for analysis sparse recovery. IEEE Trans. Signal Process. 62 (7), 1762–1774. 10.1109/tsp.2014.2304932 Tian L. Qin S. (2022). “Seismic data interpolation by the projected iterative soft-threshold algorithm for tight frame,” in 2022 4th International Conference on Image Processing and Machine Vision (IPMV), Hong Kong China, March 25 - 27, 2022. 10.1145/3529446.3529460 Yang H. Long Y. Lin J. Zhang F. Chen Z. (2017). A seismic interpolation and denoising method with curvelet transform matching filter. Acta Geophys. 65 (5), 1029–1042. 10.1007/s11600-017-0078-x Ying L. Demanet L. Candes E. (2005). 3D discrete curvelet transform. Proc. SPIE - Int. Soc. Opt. Eng. 50 (3), 351–361. 10.1117/12.616205 Zhang H. Chen X.-H. Wu X.-M. (2013). Seismic data reconstruction based on CS and Fourier theory. Appl. Geophys. 10 (2), 170–180. 10.1007/s11770-013-0375-3 Zhang H. Chen X.-H. Zhang L.-Y. (2017). 3D simultaneous seismic data reconstruction and noise suppression based on the curvelet transform. Appl. Geophys. 14 (1), 87–95. 10.1007/s11770-017-0607-z Zhang H. Diao S. Chen W. Huang G. Li H. Bai M. (2019). Curvelet reconstruction of non-uniformly sampled seismic data using the linearized Bregman method. Geophys. Prospect. 67 (5), 1201–1218. 10.1111/1365-2478.12762