Multi-trace post-stack seismic data sparse inversion with nuclear norm constraint

2020 ◽  
Author(s):  
Ronghuo Dai ◽  
Fanchang Zhang ◽  
Cheng Yin ◽  
Ying Hu
2018 ◽  
Vol 23 (5) ◽  
pp. 627-635
Author(s):  
Songze Tang ◽  
Heng Li ◽  
Liang Xiao

Geophysics ◽  
2020 ◽  
pp. 1-60
Author(s):  
Ouyang Shao ◽  
Lingling Wang ◽  
Xiangyun Hu ◽  
Zhidan Long

Because there are many similar geological structures underground, seismic profiles have an abundance of self-repeating patterns. Thus, we can divide a seismic profile into groups of blocks with similar seismic structure. The matrix formed by stacking together similar blocks in each group should be of low rank. Hence, we can transfer the seismic denoising problem to a serial of low-rank matrix approximation (LRMA) problem. The LRMA-based model commonly adopts the nuclear norm as a convex substitute of the rank of a matrix. However, the nuclear norm minimization (NNM) shrinks the different rank components equally and may cause some biases in practice. Recently introduced truncated nuclear norm (TNN) has been proven to more accurately approximate the rank of a matrix, which is given by the sum of the set of smallest singular values. Based on this, we propose a novel denoising method using truncated nuclear norm minimization (TNNM). The objective function of this method consists of two terms, the F-norm data fidelity and a truncated nuclear norm regularization. We present an efficient two-step iterative algorithm to solve this objective function. Then, we apply the proposed TNNM algorithm to groups of blocks with similar seismic structure, and aggregate all resulting denoised blocks to get the denoised seismic data. We update the denoised results during each iteration to gradually attenuate the heavy noise. Numerical experiments demonstrate that, compared with FX-Decon, the curvelet, and the NNM-based methods, TNNM not only attenuates noise more effectively even when the SNR is as low as -10 dB and seismic data have complex structures, but also accurately preserves the seismic structures without inducing Gibbs artifacts.


Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. V145-V152 ◽  
Author(s):  
Ketil Hokstad ◽  
Roger Sollie

The basic theory of surface-related multiple elimination (SRME) can be formulated easily for 3D seismic data. However, because standard 3D seismic acquisition geometries violate the requirements of the method, the practical implementation for 3D seismic data is far from trivial. A major problem is to perform the crossline-summation step of 3D SRME, which becomes aliased because of the large separation between receiver cables and between source lines. A solution to this problem, based on hyperbolic sparse inversion, has been presented previously. This method is an alternative to extensive interpolation and extrapolation of data. The hyperbolic sparse inversion is formulated in the time domain and leads to few, but large, systems of equations. In this paper, we propose an alternative formulation using parabolic sparse inversion based on an efficient weighted minimum-norm solution that can be computed in the angular frequency domain. The main advantage of the new method is numerical efficiency because solving many small systems of equations often is faster than solving a few big ones. The method is demonstrated on 3D synthetic and real data with reflected and diffracted multiples. Numerical results show that the proposed method gives improved results compared to 2D SRME. For typical seismic acquisition geometries, the numerical cost running on 50 processors is [Formula: see text] per output trace. This makes production-scale processing of 3D seismic data feasible on current Linux clusters.


Electronics ◽  
2021 ◽  
Vol 10 (21) ◽  
pp. 2667
Author(s):  
Xiaodong Yu ◽  
Rui Ding ◽  
Jingbo Shao ◽  
Xiaohui Li

Due to the high dimensionality and high data redundancy of hyperspectral remote sensing images, it is difficult to maintain the nonlinear structural relationship in the dimensionality reduction representation of hyperspectral data. In this paper, a feature representation method based on high order contractive auto-encoder with nuclear norm constraint (CAE-HNC) is proposed. By introducing Jacobian matrix in the CAE of the nuclear norm constraint, the nuclear norm has better sparsity than the Frobenius norm and can better describe the local low dimension of the data manifold. At the same time, a second-order penalty term is added, which is the Frobenius norm of the Hessian matrix expressed in the hidden layer of the input, encouraging a smoother low-dimensional manifold geometry of the data. The experiment of hyperspectral remote sensing image shows that CAE-HNC proposed in this paper is a compact and robust feature representation method, which provides effective help for the ground object classification and target recognition of hyperspectral remote sensing image.


2020 ◽  
Author(s):  
Mengyao Jiao ◽  
Tianyue Hu ◽  
Weikang Kuang ◽  
Yang Liu ◽  
Shaohuan Zu

2021 ◽  
pp. 1-57
Author(s):  
Chen Liang ◽  
John Castagna ◽  
Marcelo Benabentos

Sparse reflectivity inversion of processed reflection seismic data is intended to produce reflection coefficients that represent boundaries between geological layers. However, the objective function for sparse inversion is usually dominated by large reflection coefficients which may result in unstable inversion for weak events, especially those interfering with strong reflections. We propose that any seismogram can be decomposed according to the characteristics of the inverted reflection coefficients which can be sorted and subset by magnitude, sign, and sequence, and new seismic traces can be created from only reflection coefficients that pass sorting criteria. We call this process reflectivity decomposition. For example, original inverted reflection coefficients can be decomposed by magnitude, large ones removed, the remaining reflection coefficients reconvolved with the wavelet, and this residual reinverted, thereby stabilizing inversions for the remaining weak events. As compared with inverting an original seismic trace, subtle impedance variations occurring in the vicinity of nearby strong reflections can be better revealed and characterized when only the events caused by small reflection coefficients are passed and reinverted. When we apply reflectivity decomposition to a 3D seismic dataset in the Midland Basin, seismic inversion for weak events is stabilized such that previously obscured porous intervals in the original inversion, can be detected and mapped, with good correlation to actual well logs.


Sign in / Sign up

Export Citation Format

Share Document