Enhanced Small Reflections Sparse-Spike Seismic Inversion with Iterative Hybrid Thresholding Algorithm
Seismic inversion is a process of imaging or predicting the spatial and physical properties of underground strata. The most commonly used one is sparse-spike seismic inversion with sparse regularization. There are many effective methods to solve sparse regularization, such as L0-norm, L1-norm, weigh...
Saved in:
Main Authors: | , , |
---|---|
Format: | Article |
Language: | English |
Published: |
MDPI AG
2024-12-01
|
Series: | Mathematics |
Subjects: | |
Online Access: | https://www.mdpi.com/2227-7390/13/1/37 |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
_version_ | 1841549179844820992 |
---|---|
author | Yue Feng Ronghuo Dai Zidan Fan |
author_facet | Yue Feng Ronghuo Dai Zidan Fan |
author_sort | Yue Feng |
collection | DOAJ |
description | Seismic inversion is a process of imaging or predicting the spatial and physical properties of underground strata. The most commonly used one is sparse-spike seismic inversion with sparse regularization. There are many effective methods to solve sparse regularization, such as L0-norm, L1-norm, weighted L1-norm, etc. This paper studies the sparse-spike inversion with L0-norm. It is usually solved by the iterative hard thresholding algorithm (IHTA) or its faster variants. However, hard thresholding algorithms often lead to a sharp increase or numerical oscillation of the residual, which will affect the inversion results. In order to deal with this issue, this paper attempts the idea of the relaxed optimal thresholding algorithm (ROTA). In the solution process, due to the particularity of the sparse constraints in this convex relaxation model, this model can be considered as a L1-norm problem when dealt with the location of non-zero elements. We use a modified iterative soft thresholding algorithm (MISTA) to solve it. Hence, it forms a new algorithm called the iterative hybrid thresholding algorithm (IHyTA), which combines IHTA and MISTA. The synthetic and real seismic data tests show that, compared with IHTA, the results of IHyTA are more accurate with the same SNR. IHyTA improves the noise resistance. |
format | Article |
id | doaj-art-d822cb7b720148f5a9cf0220fb5a9d8b |
institution | Kabale University |
issn | 2227-7390 |
language | English |
publishDate | 2024-12-01 |
publisher | MDPI AG |
record_format | Article |
series | Mathematics |
spelling | doaj-art-d822cb7b720148f5a9cf0220fb5a9d8b2025-01-10T13:18:02ZengMDPI AGMathematics2227-73902024-12-011313710.3390/math13010037Enhanced Small Reflections Sparse-Spike Seismic Inversion with Iterative Hybrid Thresholding AlgorithmYue Feng0Ronghuo Dai1Zidan Fan2School of Mathematics and Information, China West Normal University, Nanchong 637009, ChinaSchool of Mathematics and Information, China West Normal University, Nanchong 637009, ChinaSchool of Mathematics and Information, China West Normal University, Nanchong 637009, ChinaSeismic inversion is a process of imaging or predicting the spatial and physical properties of underground strata. The most commonly used one is sparse-spike seismic inversion with sparse regularization. There are many effective methods to solve sparse regularization, such as L0-norm, L1-norm, weighted L1-norm, etc. This paper studies the sparse-spike inversion with L0-norm. It is usually solved by the iterative hard thresholding algorithm (IHTA) or its faster variants. However, hard thresholding algorithms often lead to a sharp increase or numerical oscillation of the residual, which will affect the inversion results. In order to deal with this issue, this paper attempts the idea of the relaxed optimal thresholding algorithm (ROTA). In the solution process, due to the particularity of the sparse constraints in this convex relaxation model, this model can be considered as a L1-norm problem when dealt with the location of non-zero elements. We use a modified iterative soft thresholding algorithm (MISTA) to solve it. Hence, it forms a new algorithm called the iterative hybrid thresholding algorithm (IHyTA), which combines IHTA and MISTA. The synthetic and real seismic data tests show that, compared with IHTA, the results of IHyTA are more accurate with the same SNR. IHyTA improves the noise resistance.https://www.mdpi.com/2227-7390/13/1/37sparse regularizationiterative hybrid thresholdingoptimal thresholdingL0-normsparse spike inversion |
spellingShingle | Yue Feng Ronghuo Dai Zidan Fan Enhanced Small Reflections Sparse-Spike Seismic Inversion with Iterative Hybrid Thresholding Algorithm Mathematics sparse regularization iterative hybrid thresholding optimal thresholding L0-norm sparse spike inversion |
title | Enhanced Small Reflections Sparse-Spike Seismic Inversion with Iterative Hybrid Thresholding Algorithm |
title_full | Enhanced Small Reflections Sparse-Spike Seismic Inversion with Iterative Hybrid Thresholding Algorithm |
title_fullStr | Enhanced Small Reflections Sparse-Spike Seismic Inversion with Iterative Hybrid Thresholding Algorithm |
title_full_unstemmed | Enhanced Small Reflections Sparse-Spike Seismic Inversion with Iterative Hybrid Thresholding Algorithm |
title_short | Enhanced Small Reflections Sparse-Spike Seismic Inversion with Iterative Hybrid Thresholding Algorithm |
title_sort | enhanced small reflections sparse spike seismic inversion with iterative hybrid thresholding algorithm |
topic | sparse regularization iterative hybrid thresholding optimal thresholding L0-norm sparse spike inversion |
url | https://www.mdpi.com/2227-7390/13/1/37 |
work_keys_str_mv | AT yuefeng enhancedsmallreflectionssparsespikeseismicinversionwithiterativehybridthresholdingalgorithm AT ronghuodai enhancedsmallreflectionssparsespikeseismicinversionwithiterativehybridthresholdingalgorithm AT zidanfan enhancedsmallreflectionssparsespikeseismicinversionwithiterativehybridthresholdingalgorithm |