期刊
IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING
卷 60, 期 -, 页码 -出版社
IEEE-INST ELECTRICAL ELECTRONICS ENGINEERS INC
DOI: 10.1109/TGRS.2022.3223280
关键词
Alternating iteration; compact smoothness; reflectivity inversion; relative sparsity; wavelet inversion
类别
资金
- National Key Research and Development Program of China [2018YFC0603501, 2020YFA0713403]
- National Natural Science Foundation of China [42174137, 41974137]
- Natural Science Foundation of Anhui Province [1908085MF218]
This paper proposes a new algorithm called compact smoothness and relative sparsity (CSRS) algorithm for wavelet and reflectivity inversion (WRI) in seismic data processing. The algorithm introduces normalized constraints and regularization to overcome the ill-posedness of WRI inversion and improve the lateral continuity of reflectivities. It also presents an efficient algorithm for solving the optimization problem of WRI by dividing it into reflectivity inversion and wavelet inversion subproblems.
Wavelet and reflectivity inversion (WRI) is an important issue in seismic data processing. To overcome the ill-posedness of WRI inversion with more efficient parameter selection and better lateral continuity of reflectivities, we propose a new WRI algorithm named compact smoothness and relative sparsity (CSRS) algorithm, where a normalized compact constraint and a normalized smooth regularization are proposed for the wavelet inversion, and a relative sparsity constraint is proposed for the reflectivity inversion. The proposed constraints and regularization make the parameters of WRI easy to be selected. The proposed relative sparsity constraint can lead to a reflectivity profile with good lateral continuity, as it can be suitable for various seismic data with a fixed sparsity parameter. We also propose an efficient algorithm for solving corresponding WRI optimization problem. The whole WRI problem is divided into reflectivity inversion subproblem and wavelet inversion subproblem by using alternating iterative method, where the initial wavelet is estimated by smoothing the absolute amplitude spectrum of averaged seismic data. The proximal algorithm is applied to solve both reflectivity inversion subproblem and wavelet inversion subproblem. By replacing Toeplitz matrix multiplication with the fast Fourier transform (FFT) and using compact wavelet, our algorithm can be efficient for 3-D seismic data. The numerical examples on 2-D synthetic data, 2-D offshore field data, and 3-D onshore field data demonstrate that, compared to the Toeplitz-sparse matrix factorization (TSMF) algorithm, the CSRS algorithm with fixed default parameters can get high-resolution reflectivities with better lateral continuity, and requires much less computational time.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据