Abstract
The compressed sensing (CS) method, commonly utilized for restructuring sparse signals, has been extensively used to attenuate the random noise in seismic data. An important basis of CS-based methods is the sparsity of sparse coefficients. In this method, the sparse coefficient vector is acquired by minimizing the norm as a substitute for the norm. Many efforts have been made to minimize the norm (0 < p < 1) to obtain a more desirable sparse coefficient representation. Despite the improved performance that is achieved by minimizing the norm with 0 < p < 1, the related sparse coefficient vector is still suboptimal since the parameter p is greater than 0 rather than infinitely approaching 0 . Therefore, the CS method with the limit is proposed to enhance the sparse performance and thus generate better denoised results in this paper. Our proposed method is referred to as the CS-LHR method because the solving process for minimizing is the log-sum heuristic recovery (LHR). Furthermore, to improve the computational efficiency, we incorporate the majorization-minimization (MM) algorithm in this CS-LHR method. Experimental results of synthetic and real seismic records demonstrate the remarkable performance of CS-LHR in random noise suppression.
1 Introduction
Random noise is frequently present in raw seismic data, which disrupts the continuity of seismic events and reduces the signal-to-noise ratio (SNR) of seismic data. Low SNR and discontinuous seismic events can blur the stratigraphic information in seismic profiles, reduce the interpretability of seismic data, and lead to incorrect identification of subsurface targets. Hence, it is essential to perform seismic noise separation and attenuation during both prestack and poststack seismic data processing (Wu et al., 2019; Dong et al., 2022a; 2022b; Liu et al., 2022a; Liu et al., 2022b; Wu B Y et al., 2022; Zhong et al., 2022; Zhong et al., 2023).
In recent years, numerous signal processing methods have emerged for the separation and suppression of seismic noise (Yuan et al., 2012; Li et al., 2014; 2022; Zhang et al., 2021; Ni et al., 2022; Sun et al., 2022; Wu H et al., 2022). These methods include singular spectrum analysis (Oropeza and Sacchi, 2011), empirical mode decomposition-based techniques (Bekara and Baan, 2009), wavelet transform (Yang et al., 2018), and curvelet transform (Qu et al., 2016). Most of these methods are typically developed based on the distinguishing characteristics of seismic signals and specific types of noise in transform domains. Notably, sparse representation-based techniques have gained significant popularity (Candès et al., 2006; Chen et al., 2017; Wu B Y et al., 2022). While seismic data is not inherently sparse, it can be effectively transformed into a sparse signal by sparse transformation (Siahsar et al., 2016). Random noise cannot be transformed into a sparse signal due to lacking sparsity. Then, during sparse transformation, the noisy seismic signal is separated into a sparse signal and random noise. Subsequently, the denoised seismic signal is reconstructed using the sparse signal, thereby the separation of the seismic signal and random noise is achieved by sparse transformation and sparse signal reconstruction. In practical applications, the denoising effectiveness of sparse transformation is linked to the sparsity of the resulting sparse signal. Greater sparsity leads to improved denoising performance. Thus, enhancing the sparsity of sparse transformation is crucial for its denoising applications (Wu H et al., 2022).
Compressed sensing (CS) is a well-established method that combines sparse transformation and signal reconstruction (Donoho, 2006). In contrast to the conventional Nyquist–Shannon sampling theory, the CS method can reconstruct signals without higher sampling rates and has received significant attention and been widely applied in separating random noise. Regrettably, obtaining sparse signals through the norm minimization in the CS method is an NP-hard problem (Candès and Wakin, 2008; Yang et al., 2022). As such, the minimization of norm and norm (0<p<1) are often adopted as replacements for the minimization of norm in some improved CS methods (Yang et al., 2009; Liu et al., 2023a). And the norm (0<p<1) minimization has been demonstrated as having superior sparsity capabilities compared to norm minimization (Wu B Y et al., 2022; Liu et al., 2023b). Although the aforementioned methods for tackling NP-hard problems can enhance the convergence and effectiveness of the solution process, they also diminish the sparsity of the sparse signal; their sparse performance can still be further enhanced by incorporating the limiting form .
Thus, a novel algorithm leveraging norm minimization with in the CS method is proposed in this paper; it can enhance the sparsity of sparse signals, achieve proficient signal reconstruction, and effectively suppress seismic noise. The minimization utilizing the limiting form is referred to as log-sum heuristic recovery (LHR) because the expansion of the limiting form norm is a logarithmic sum (Zou and Hastie, 2015). Therefore, we also call the proposed algorithm the CS-LHR method. In our approach, the minimization with poses a non-convex problem, making its solution process more intricate than that of convex problems. Encouragingly, significant progress has been achieved in addressing non-convex problems, and the majorization-minimization (MM) iterative optimization algorithm is one effective method for solving such problems. The MM algorithm substitutes a complex optimization problem with a series of simpler ones, thus approximating the objective function that encompasses non-differentiable and non-convex traits with a differentiable and convex surrogate function to facilitate optimal solution retrieval (Fazel et al., 2003; Foo et al., 2009). To ensure good convergence rates, our workflow incorporates the majorization-minimization (MM) algorithm.
In the subsequent sections of this article, we provide a detailed description of the proposed workflow. Subsequently, a synthetic dataset and a field dataset containing noise are utilized to demonstrate the effectiveness of the approach. The results show that our method can suppress noise from seismic reflections effectively and results in a seismic profile with good continuity of seismic events and high resolution.
2 Compressed sensing with the limit form (CS-LHR)
Compressed sensing (CS), which challenges the traditional Nyquist–Shannon sampling theory, has emerged as a hot topic in the field of signal processing (Donoho, 2006.). Although many studies on the applications of CS have been conducted, there is still value in exploring how to enhance its performance (Candès and Wakin, 2008; Yang et al., 2009). In this paper, we will explore how to enhance the sparsity of sparse signals in the CS method and apply the related research to seismic signal denoising.
If a signal can be sparsely represented, it can be written aswhere the matrix consists of a set of sparse basis vectors, and the vector represents a sparse signal within the space defined by these basis vectors in matrix . In the CS method, the sparse signal can be obtained by minimizing the norm as and should satisfy the constraint as
In Eq. 2 is referred to as the sensing matrix, and is the sampled data of obtained by the sensing matrix .
Then, the CS theory can be described by
In practice, the matrix and the sensing matrix are pre-determined. As the description by Eq. 1 and Eq. 3, the sparse signal , when acquired, allows for the reconstruction of X utilizing .
Eq. 3 presented above is an NP-hard problem that is challenging to solve. To overcome this NP-hard problem, the optimization of norm in Eq. 3 can be replaced with a convex optimization of norm which is the convex approximation of the norm, and easier to solve. Then, the related CS method with the norm is (Candès and Wakin, 2008)where represents the norm.
Although the convex relaxation in Eq. 4 reduces the complexity of the original NP-hard problem, it unfortunately yields a solution with suboptimal sparsity due to the norm deviating significantly from the norm. To address this issue, the norm is introduced (Candès et al., 2006). Subsequently, the optimization problem described by Eq. 4 with the norm can be written as
In which and p (0,1]. For , is equivalent toin Eq. 6 N is the length of the sparse signal and denotes the element of (Caiafa and Cichocki, 2013).
Presented above corresponds to a non-convex optimization and exhibits superior sparsity performance compared to the norm minimization. Although the advantages of CS method with the p (0,1] norm shown as Eq. 5 have been demonstrated, the CS method with limit has not been studied. It is important that the norm minimization based on differs from other p values in p (0,1]; it possesses greater sparse capability (Deng et al., 2012). Additionally, the norm minimization based on also differs from (p=0) that yields an NP-hard problem; it is solvable. Thus, in order to acquire a sparse signal with high sparsity, we propose a novel approach that combines the optimization of the limit norm with the CS method to enhance the sparsity of , described as
According to L’Hôspital’s rule (Caiafa and Cichocki, 2013), in Eq. 7 can be expressed aswhere is a small positive number to guarantee the stability of the algorithm. In practice, should be set to a value slightly smaller than the expected non-zero element . Typically, the solve process of Eq. 8 is robust enough to tolerate different choices of . Therefore, combined with Eq. 8, 7 can be rewritten asin which the logarithmic sum is denoted as , and is the log-sum heuristic recovery (LHR) model. Therefore, the improved method described by Eq. 9 is named as the CS-LHR method by us because it is the composition of the CS and the LHR. In contrast to the traditional CS methods, the CS-LHR method can attain a best sparse signal that exhibits the optimal sparsity.
Note that Eq. 9 is non-convex due to the non-convexity of its log-sum. According to recent progress in non-convex optimization, the non-convex problem can be solved efficiently. In this paper, we incorporate the alternating direction method of the majorization-minimization (MM) algorithm into our workflow to ensure faster convergence (Fazel et al., 2003; Foo et al., 2009). The MM algorithm transforms the original non-differentiable, non-convex function into a differentiable and convex surrogate function, facilitating the retrieval of optimal solutions. Then, Eq. 9 can be equivalently expressed as Eq. 10 based on the MM algorithm aswhere represents the vector of weighted parameters with each element . Eq. 10 demonstrates that the log-sum penalty function performs the re-weighted minimization, which promotes sparsity more effectively compared to the norm (0<p<1) minimization. Moreover, each iteration of the MM algorithm for solving Eq. 10 corresponds to a convex optimization that can be easily solved.
Eq. 10 can also be rewritten asin Eq. 11 is the positive weighting parameter. Once the solution of Eq. 10 is obtained, the signal can be recovered.
3 Seismic denoising by the CS-LHR
Section 2 suggests that the CS-LHR can achieve the optimal sparse signal through the limit norm minimization. The resulting optimal sparse signal complies with the constraints and is well-suited for signal reconstruction. This paper focuses on the application of the CS-LHR method to seismic signal denoising. A general form of an observed seismic signal that is contaminated by noise can be expressed asin Eq. 12 is a sparse signal, represents the noise term that can either be stochastic or deterministic, and represents the index of time sampling point. Assuming that and represent the sensing and sparse basis matrices, and are independent as the CS theory, it becomes feasible to separate the noise from the observed seismic signal . Subsequently, the denoised result can be reconstructed by . Moreover, the effectiveness of denoising is associated with the sparsity of . Greater sparsity leads to improved noise separation. Based on the foregoing analysis, the CS-LHR, utilizing limit norm optimization, yields an optimally sparse signal .
According to the CS theory, the sparse basis matrix and the sensing matrix should be irrelevant (Donoho, 2006). However, achieving complete independence between and in practical applications is challenging. Then, some special matrices are chosen as to ensure a certain degree of independence with the sparse basis matrix . In this study, a Gaussian random matrix is chosen as sensing matrix due to its excellent characteristic of having minimal correlation with other matrices. Additionally, since a seismic signal is non-stationary, the sparse basis matrix should be provided by an algorithm which facilitates the analysis of non-stationary signals. To obtain a sparse basis matrix for non-stationary signals, the sparse S-transform is introduced (Wang et al., 2016). Furthermore, as the log-sum penalty term in Eq. 8 is non-convex, a suitable starting point for iterative computation is necessary. Consequently, we initialize with the solution of Eq. 3 with the . The proposed workflow is summarized in Algorithm 1 (Table 1):
TABLE 1
| Algorithm 1 Workflow of seismic denoising by CS-LHR |
|---|
| Input: Observed seismic data , the sensing matrix , the sparse basis , the positive parameter ; |
| Initialization: Initialize from Eq. 3. Determine each through ; |
| Repeat: Update and determine each by until convergence; |
| Output: The sparse coefficient vector ; |
| End: Recover the free-noise . |
The workflow of seismic denoising by CS-LHR.
4 Synthetic and real data examples
In order to illustrate the effectiveness of the proposed CS-LHR method, we initially apply it to synthetic seismic data with different levels of signal noise ratios (SNRs). Then, the proposed method is utilized for field data denoising. Figure 1 displays a 2-D synthetic seismic trace without noise. Figure 2 shows the 2-D synthetic trace with different SNRs (5dB, 5dB, -3dB, and 3 dB). Figure 3 exhibits a 3-D noisy field data acquired over the Scotian shelf, offshore Canada, and termed Penobscot. For comparison, the traditional CS method based on norm (0<p<1) is utilized as an alternative method.
FIGURE 1
FIGURE 2
FIGURE 3
4.1 Seismic signal enhancement with different SNRs
The denoising results for noisy 2-D synthetic data using the traditional CS method are depicted in Figures 4A,B, while those obtained from the CS-LHR method are presented in Figures 5A,B. It is clear that the traditional CS method can effectively attenuate noise in smooth areas of noisy synthetic data; however, it introduces artifacts and exhibits a poor denoising effect in the oscillatory areas marked by black ellipses. Although the traditional CS method based on norm (0<p<1) exhibits greater effect on noise attenuation compared to that based on norm, its sparse signal exhibits varying reconstruction capabilities across different regions of the signals. During noise separation using sparse transformation in the traditional CS method, the suboptimal sparsity of the sparse signal leads to the inclusion of some noise characteristics in the sparse signal, which become apparent in the reconstructed original signal. While the sparse signal effectively captures the primary information within smooth areas of the signal, it also incorporates some noise features in oscillatory regions. To enhance the denoising performance across all areas of the noisy signal, it is crucial that the sparsity of the sparse signal is increased to enhance the reconstruction effectiveness of the original signal. Consequently, this paper introduces the CS-LHR method which can achieve the best sparse reconstruction ability due to the norm minimization based on . The denoising results shown in Figure 5 illustrate the random noises are successfully removed while the seismic events are preserved well. Notably, the proposed method demonstrates exceptional noise filtering capabilities, even under low signal-to-noise ratio (SNR) conditions.
FIGURE 4
FIGURE 5
4.2 Field data applications
To verify the effectiveness of the proposed method, 3-D noise-contaminated field data obtained from the Scotian shelf, offshore Canada, referred to as Penobscot, are shown in Figure 3. The 3-D field data comprise 401 inlines and 401 crosslines, with a time sampling interval of 4 ms.
In Figure 3, black arrows indicate discontinuous seismic events, while green arrows represent seismic artifacts caused by random noise. The green lines correspond to the location of X-line 1273, as depicted in Figure 6A. Obviously, this seismic volume contains significant random noises which hinder subsequent seismic data processing and interpretation. Our method’s denoising result is shown in Figure 6B, where the improved resolution and well-preserved reflection events are evident. Regions marked by the black ellipses demonstrate efficient attenuation of random noise, enhanced continuity, and resolution of seismic events. Additionally, the seismic fault structures indicated by black arrows are preserved well. Further, Figure 6C shows no useful information in the difference profile.
FIGURE 6
To further demonstrate the effectiveness of our method, we compare it with the traditional CS method with 0<p<1 on the same field data. The corresponding results are presented in Figure 7. Figure 7A shows the denoising result, and Figure 7B represents the related difference profile. We can observe that valid seismic events are generally preserved in Figure 7A. However, compared to the CS-LHR result, the fault structures and seismic events, indicated by black arrows and black ellipses in Figure 7A, respectively, are less subtle. Additionally, some valuable information that can improve the resolution of the denoising result is contained in the difference profile Figure 7B.
FIGURE 7
5 Conclusion
This paper proposes the CS-LHR method, a novel method for seismic noise attenuation. Compared to the traditional CS methods with 0<p 1, the CS-LHR method with the limit provides enhanced sparse representation ability and denoising performance. Testing results on field data demonstrate that our workflow efficiently recovers noise-free signals. Additionally, we implement the MM algorithm to improve calculation efficiency.
The CS method can be used for denoising, but its primary contribution to the scientific domain lies in accomplishing the compression and reconstruction of original signals via sparse signal representation. This process facilitates the reduction of data acquisition and transmission costs while preserving data quality, essential for diverse applications including medical imaging, remote diagnosis, earth observation, and wireless transmission. The CS-LHR method introduced in this paper can achieve the optimal sparsity of the sparse signal, leading to additional reductions in storage and transmission costs. This holds particular significance for industrial applications driven by cost considerations.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
FS: Conceptualization, Methodology, Writing–original draft. QZ: Formal Analysis, Writing–review and editing. ZW: Validation, Writing–review and editing. WH: Validation, Improving the quality, Writing–review and editing, Formal Analysis.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the Natural Science Foundation of Guangxi (2021GXNSFBA196071 and 2021JJA170177), Guangxi Key Laboratory of Precision Navigation Technology and Application, and Key Laboratory of Microwave and Optical Wave Application Technology, Guilin University of Electronic Technology.
Acknowledgments
The authors would like to thank the dGB Earth Sciences, Nova Scotia Department of Energy, Canada-Nova Scotia Offshore Petroleum Board, for providing the 3-D field data used in this study.
Conflict of interest
Author ZW was employed by The china state shipbuilding corporation limited.
Author WH was employed by Xi’an Engineering Investigation and Design Research Institute of China National Nonferrous Metals Industry Co., Ltd.
The remaining 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
1
BekaraM.BaanM. V. D. (2009). Random and coherent noise attenuation by empirical mode decomposition. Geophysics74 (5), 89–98. 10.1190/1.3157244
2
CaiafaC. F.CichockiA. (2013). Multidimensional compressed sensing and their applications. Wires Data Min. Knowl. Discov.3, 355–380. 10.1002/widm.1108
3
CandèsE. J.RombergJ.TaoT. (2006). Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information. IEEE Trans. Inf. Theory52, 489–509. 10.1109/TIT.2005.862083
4
CandèsE. J.WakinM. B. (2008). An introduction to compressive sampling. IEEE Signal Process. Mag.25, 21–30. 10.1109/MSP.2007.914731
5
ChenY. K.ZhouY. T.ChenW.ZuS.HuangW.ZhangD. (2017). Empirical low-rank approximation for seismic noise attenuation. IEEE Trans. Geoscience Remote Sens.55 (8), 4696–4711. 10.1109/TGRS.2017.2698342
6
DengY.LiuY. B.DaiQ. H.ZhangZ.WangY. (2012). Noisy depth maps fusion for multiview stereo via matrix completion. IEEE J. Sel. Top. Signal Process.6, 566–582. 10.1109/JSTSP.2012.2195472
7
DongX. T.LinJ.LuS. P.HuangX.WangH.LiY. (2022a). Seismic shot gather denoising by using a supervised-deep-learning method with weak dependence on real noise data: a solution to the lack of real noise data. Surv. Geophys.43 (5), 1363–1394. 10.1007/s10712-022-09702-7
8
DongX. T.LinJ.LuS. P.WangH.LiY. (2022b). Multiscale spatial attention network for seismic data denoising. IEEE Trans. Geoscience Remote Sens.60, 1–17. 10.1109/TGRS.2022.3178212
9
DonohoD. L. (2006). Compressed sensing. IEEE Trans. Inf. Theory52, 1289–1306. 10.1109/TIT.2006.871582
10
FazelM.HindiH.BoydS. P. (2003). “Log-det heuristic for matrix rank minimization with applications to Hankel and Euclidean distance matrices,” in The 2003 American Control Conference, Denver, CO, USA, 04-06 June 2003, 2156–2172.
11
FooC. S.DoC. B.NgA. Y. (2009). “A majorization-minimization algorithm for (multiple) hyper-parameter learning,” in The 26th annual international conference on machine learning (Montreal, Quebec, CA: Association for Computing Machinery), 321–328. 10.1145/1553374.1553415
12
LiF. Y.SunF. Y.LiuN. H.XieR. (2022). Denoising seismic signal via resampling local applicability functions. IEEE Geoscience Remote Sens. Lett.19, 1–5. 10.1109/LGRS.2020.3048110
13
LiF. Y.ZhangB.MarfurtK. J.HallI. (2014). “Random noise suppression using normalized convolution filter,” in The SEG technical program expanded abstracts 2014 (Denver, USA: SEG). 10.1190/segam2014-1478.1
14
LiuN. H.LeiY. B.LiuR. C.YangY.WeiT.GaoJ. (2023a). Sparse time-frequency analysis of seismic data: sparse representation to unrolled optimization. IEEE Trans. Geoscience Remote Sens.61, 1–10. 10.1109/TGRS.2023.3300578
15
LiuN. H.LeiY. B.YangY.WangZ.LiuR.GaoJ.et al (2023b). Sparse time-frequency analysis of seismic data via convolutional neural network. Interpretation12 (1), T47–T62. 10.1190/int-2023-0020.1
16
LiuN. H.WangJ. L.GaoJ. H.ChangS.LouY. (2022a). Similarity-informed self-learning and its application on seismic image denoising. IEEE Trans. Geoscience Remote Sens.60, 1–13. 10.1109/tgrs.2022.3210217
17
LiuN. H.WangJ. L.GaoJ. H.YuK.LouY.PuY.et al (2022b). NS2NS: self-learning for seismic image denoising. IEEE Trans. Geoscience Remote Sens.60, 1–11. 10.1109/TGRS.2022.3217289
18
NiW. D.LouY. H.XuX. L.XueA.WangW. (2022). “Seismic random noise attenuation via noise assisted-multivariate EMD based MSSA,” in The 2022 International Conference on Automation, Robotics and Computer Engineering (ICARCE), Wuhan, China, 16-17 December 2022. 10.1109/ICARCE55724.2022.10046637
19
OropezaV.SacchiM. (2011). Simultaneous seismic data denoising and reconstruction via multichannel singular spectrum analysis. Geophysics76 (3), 25–32. 10.1190/1.3552706
20
QuS.ZhouH.LiuR. W.ChenY.ZuS.YuS.et al (2016). Deblending of Simultaneous-source seismic data using fast iterative shrinkage-thresholding algorithm with firm-thresholding. Acta Geophys.64 (4), 1064–1092. 10.1515/acgeo-2016-0043
21
SiahsarM. A. N.GholtashiS.KahooA. R.MarviH.AhmadifardA. (2016). Sparse time-frequency representation for seismic noise reduction using low-rank and sparse decomposition. Geophysics81 (2), 117–124. 10.1190/geo2015-0341.1
22
SunF. Y.LiaoG. S.LouY. H.JiangX. (2022). Seismic data denoising with correlation feature optimization via S-mean. IEEE Geoscience Remote Sens. Lett.19, 1–5. 10.1109/LGRS.2021.3117965
23
WangY. Q.PengZ. M.HeY. M. (2016). “Time-frequency representation for seismic data using sparse S transform,” in The2nd IEEE International Conference on Computer and Communications (ICCC), Chengdu, China, 14-17 October 2016. 10.1109/CompComm.2016.7925036
24
WuB. Y.YuJ. Y.RenH. R.LouY. H.LiuN. H. (2022). Seismic traffic noise attenuation using $l_{p}$ -norm robust PCA. IEEE Geoscience Remote Sens. Lett.17, 1998–2001. 10.1109/lgrs.2019.2955737
25
WuH.ZhangB.LinT. F.LiF.LiuN. (2019). White noise attenuation of seismic trace by integrating variational mode decomposition with convolutional neural network. Geophysics84 (5), 307–317. 10.1190/geo2018-0635.1
26
Wu HH.ZhangB.LiuN. H. (2022). Self-adaptive denoising net: self-supervised learning for seismic migration artifacts and random noise attenuation. J. Petroleum Sci. Eng.214, 110431. 10.1016/j.petrol.2022.110431
27
YangJ. Y.PengY. G.XuW. L.DaiQ. H. (2009). Ways to sparse representation: an overview. Sci. China (Series F:Information Sci.52, 695–703. 10.1007/s11432-009-0045-5
28
YangY.LeiY. B.LiuN. H.WangZ.GaoJ.DingJ. (2022). Sparse TFNet: a physically informed auto encoder for sparse time–frequency analysis of seismic data. IEEE Trans. Geoscience Remote Sens.60, 1–12. 10.1109/TGRS.2022.3213851
29
YangY.LiD. Q.TongT. G.ZhangD.ZhouY.ChenY. (2018). Denoising controlled-source electromagnetic data using least-squares inversion. Geophysics83 (4), 229–244. 10.1190/geo2016-0659.1
30
YuanS. Y.WangS. X.LiG. F. (2012). Random noise reduction using Bayesian inversion. J. Geophys. Eng.9, 60–68. 10.1088/1742-2132/9/1/007
31
ZhangY. J.ZhangH. R.YangY. (2021). Seismic random noise separation and attenuation in based on MVMD and MSSA. IEEE Trans. Geoscience Remote Sens.60, 5908916. 10.1109/TGRS.2021.3131655
32
ZhongT.ChengM.DongX. T. (2023). RMCHN: a residual modular cascaded heterogeneous network for noise suppression in DAS-VSP records. IEEE Geoscience Remote Sens. Lett.20, 7500205. 10.1109/LGRS.2022.3229556
33
ZhongT.ChengM.DongX. T.LiY.WuN. (2022). Seismic random noise suppression by using deep residual U-Net. J. Petroleum Sci. Eng.209, 109901. 10.1016/j.petrol.2021.109901
34
ZouH.HastieT. (2015). Regularization and variable selection via the elastic net. J. R. Stat. Soc. Ser. B67, 301–320. 10.1111/j.1467-9868.2005.00503.x
Summary
Keywords
compressed sensing, log-sum heuristic recovery, seismic denoising, lp norm, the log-sum heuristic recovery (LHR)
Citation
Sun F, Zhang Q, Wang Z and Hou W (2024) Compressed sensing with log-sum heuristic recover for seismic denoising. Front. Earth Sci. 11:1285622. doi: 10.3389/feart.2023.1285622
Received
30 August 2023
Accepted
18 December 2023
Published
09 January 2024
Volume
11 - 2023
Edited by
Xintong Dong, Jilin University, China
Reviewed by
Hao Wu, China University of Geosciences Wuhan, China
Yihuai Lou, Zhejiang University, China
Updates
Copyright
© 2024 Sun, Zhang, Wang and Hou.
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.
*Correspondence: Fengyuan Sun, fysun@guet.edu.cn
Disclaimer
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.