A fast rank-reduction algorithm for three-dimensional seismic data interpolation

2016 ◽  
Vol 132 ◽  
pp. 137-145 ◽  
Author(s):  
Yongna Jia ◽  
Siwei Yu ◽  
Lina Liu ◽  
Jianwei Ma
Geophysics ◽  
2021 ◽  
pp. 1-92
Author(s):  
Yangkang Chen ◽  
Sergey Fomel ◽  
Hang Wang ◽  
shaohuan zu

The prediction error filter (PEF) assumes that the seismic data can be destructed to zero by applying a convolutional operation between the target data and prediction filter in either time-space or frequency-space domain. Here, we extend the commonly known PEF in 2D or 3D problems to its 5D version. To handle the non-stationary property of the seismic data, we formulate the PEF in a non-stationary way, which is called the non-stationary prediction error filter (NPEF). In the NPEF, the coefficients of a fixed-size PEF vary across the whole seismic data. In NPEF, we aim at solving a highly ill-posed inverse problem via the computationally efficient iterative shaping regularization. The NPEF can be used to denoise multi-dimensional seismic data, and more importantly, to restore the highly incomplete aliased 5D seismic data. We compare the proposed NPEF method with the state-of-the-art rank-reduction method for the 5D seismic data interpolation in cases of irregularly and regularly missing traces via several synthetic and real seismic data. Results show that although the proposed NPEF method is less effective than the rank-reduction method in interpolating irregularly missing traces especially in the case of low signal to noise ratio (S/N), it outperforms the rank-reduction method in interpolating aliased 5D dataset with regularly missing traces.


2021 ◽  
Vol 18 (4) ◽  
pp. 529-538
Author(s):  
Liyan Zhang ◽  
Ang Li ◽  
Jianguo Yang ◽  
Shichao Li ◽  
Yulai Yao ◽  
...  

Abstract To improve the imaging quality of wide-azimuth seismic data and enhance the uniformity of the attributes between adjacent bins, we developed a novel interpolation method in the offset-vector tiles (OVT) domain for wide-azimuth data. The orthogonal matching pursuit (OMP) interpolation method based on the Fourier transform is a frequency-domain processing technique based on discrete Fourier interpolation that achieves the goal of anti-aliasing by extracting the weight factor in the effective band from low-frequency data without aliasing. For data reconstruction, the OMP-based data interpolation technique in the OVT domain comprehensively uses the seismic data in five dimensions: the vertical and horizontal coordinates, time, offset and azimuth. Compared with conventional three-dimensional data interpolation, five-dimensional interpolation in the OVT domain is more accurate and achieves better results in practical applications.


2020 ◽  
Vol 1631 ◽  
pp. 012110
Author(s):  
Xiaoguo Xie ◽  
Shuling Pan ◽  
Bing Luo ◽  
Cailing Chen ◽  
Kai Chen

2021 ◽  
Author(s):  
Vladimir Cheverda ◽  
Vadim Lisitsa ◽  
Maksim Protasov ◽  
Galina Reshetova ◽  
Andrey Ledyaev ◽  
...  

Abstract To develop the optimal strategy for developing a hydrocarbon field, one should know in fine detail its geological structure. More and more attention has been paid to cavernous-fractured reservoirs within the carbonate environment in the last decades. This article presents a technology for three-dimensional computing images of such reservoirs using scattered seismic waves. To verify it, we built a particular synthetic model, a digital twin of one of the licensed objects in the north of Eastern Siberia. One distinctive feature of this digital twin is the representation of faults not as some ideal slip surfaces but as three-dimensional geological bodies filled with tectonic breccias. To simulate such breccias and the geometry of these bodies, we performed a series of numerical experiments based on the discrete elements technique. The purpose of these experiments is the simulation of the geomechanical processes of fault formation. For the digital twin constructed, we performed full-scale 3D seismic modeling, which made it possible to conduct fully controlled numerical experiments on the construction of wave images and, on this basis, to propose an optimal seismic data processing graph.


Geophysics ◽  
1972 ◽  
Vol 37 (3) ◽  
pp. 417-430 ◽  
Author(s):  
G. G. Walton

The three‐dimensional seismic method is a different way of gathering and presenting seismic data. Instead of showing the subsurface beneath a profile line, 3-D displays give an, areal picture from the shallowest reflector to the deepest one that can be found seismically. Data are collected in the field with cross‐spreads that provide over 2000 evenly spaced depth points on each reflecting interface. Several variations of the cross‐spread technique give the same subsurface coverage while providing flexibility in data gathering. Because of the dense coverage, the method is best suited for problems requiring great detail, such as production problems. The usual presentation of 3-D data is a visual, moving display of emerging wavefronts covering four sq mi of surface. From this dynamic display, average velocity to each reflector and the dip direction and magnitude can be computed. The method has proved especially useful for the recognition of faults and determination of fault directions.


Geophysics ◽  
2021 ◽  
pp. 1-96
Author(s):  
Yapo Abolé Serge Innocent Oboué ◽  
Yangkang Chen

Noise and missing traces usually influence the quality of multidimensional seismic data. It is, therefore, necessary to e stimate the useful signal from its noisy observation. The damped rank-reduction (DRR) method has emerged as an effective method to reconstruct the useful signal matrix from its noisy observation. However, the higher the noise level and the ratio of missing traces, the weaker the DRR operator becomes. Consequently, the estimated low-rank signal matrix includes a unignorable amount of residual noise that influences the next processing steps. This paper focuses on the problem of estimating a low-rank signal matrix from its noisy observation. To elaborate on the novel algorithm, we formulate an improved proximity function by mixing the moving-average filter and the arctangent penalty function. We first apply the proximity function to the level-4 block Hankel matrix before the singular value decomposition (SVD), and then, to singular values, during the damped truncated SVD process. The relationship between the novel proximity function and the DRR framework leads to an optimization problem, which results in better recovery performance. The proposed algorithm aims at producing an enhanced rank-reduction operator to estimate the useful signal matrix with a higher quality. Experiments are conducted on synthetic and real 5-D seismic data to compare the effectiveness of our approach to the DRR approach. The proposed approach is shown to obtain better performance since the estimated low-rank signal matrix is cleaner and contains less amount of artifacts compared to the DRR algorithm.


Sign in / Sign up

Export Citation Format

Share Document