Interpolated Multichannel Singular Spectrum Analysis (I-MSSA): a reconstruction method that honors true trace coordinates

Geophysics ◽  
2020 ◽  
pp. 1-65
Author(s):  
Fernanda Carozzi ◽  
Mauricio D. Sacchi

The Multichannel Singular Spectrum Analysis (MSSA) reconstruction algorithm denoises and reconstructs seismic traces on a regular grid. We present a modified version of MSSA that can cope with denoising and reconstruction of traces with irregular coordinates. The proposed method, Interpolated Multichannel Singular Spectrum Analysis (I-MSSA), connects off-the-grid observations to the desired gridded data via a non-invertible bilinear interpolation operator. The algorithm consists of two steps. In the first step, we use the steepest descent method to estimate the gridded data that honors off-the-grid observations. The second step guarantees convergence to a solution by applying the MSSA filter to the gridded data. The final solution is the reconstructed volume that honors off-the-grid observations. We apply the algorithm to synthetic and field data. We also provide an application where 3D prestack data corresponding to an orthogonal survey is fully reconstructed using cross-spread gathers. We use I-MSSA to reconstruct each subset individually. The output is a complete seismic volume described in a regular CMP grid.

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. V25-V32 ◽  
Author(s):  
Vicente Oropeza ◽  
Mauricio Sacchi

We present a rank reduction algorithm that permits simultaneous reconstruction and random noise attenuation of seismic records. We based our technique on multichannel singular spectrum analysis (MSSA). The technique entails organizing spatial data at a given temporal frequency into a block Hankel matrix that in ideal conditions is a matrix of rank [Formula: see text], where [Formula: see text] is the number of plane waves in the window of analysis. Additive noise and missing samples will increase the rank of the block Hankel matrix of the data. Consequently, rank reduction is proposed as a means to attenuate noise and recover missing traces. We present an iterative algorithm that resembles seismic data reconstruction with the method of projection onto convex sets. In addition, we propose to adopt a randomized singular value decomposition to accelerate the rank reduction stage of the algorithm. We apply MSSA reconstruction to synthetic examples and a field data set. Synthetic examples were used to assess the performance of the method in two reconstruction scenarios: a noise-free case and data contaminated with noise. In both cases, we found extremely low reconstructions errors that are indicative of an optimal recovery. The field data example consists of a 2D prestack volume that depends on common midpoint and offset. We use the MSSA reconstruction method to complete missing offsets and, at the same time, increase the signal-to-noise ratio of the seismic volume.


Sensors ◽  
2021 ◽  
Vol 21 (4) ◽  
pp. 1403
Author(s):  
Xin Jin ◽  
Xin Liu ◽  
Jinyun Guo ◽  
Yi Shen

Geocenter is the center of the mass of the Earth system including the solid Earth, ocean, and atmosphere. The time-varying characteristics of geocenter motion (GCM) reflect the redistribution of the Earth’s mass and the interaction between solid Earth and mass loading. Multi-channel singular spectrum analysis (MSSA) was introduced to analyze the GCM products determined from satellite laser ranging data released by the Center for Space Research through January 1993 to February 2017 for extracting the periods and the long-term trend of GCM. The results show that the GCM has obvious seasonal characteristics of the annual, semiannual, quasi-0.6-year, and quasi-1.5-year in the X, Y, and Z directions, the annual characteristics make great domination, and its amplitudes are 1.7, 2.8, and 4.4 mm, respectively. It also shows long-period terms of 6.09 years as well as the non-linear trends of 0.05, 0.04, and –0.10 mm/yr in the three directions, respectively. To obtain real-time GCM parameters, the MSSA method combining a linear model (LM) and autoregressive moving average model (ARMA) was applied to predict GCM for 2 years into the future. The precision of predictions made using the proposed model was evaluated by the root mean squared error (RMSE). The results show that the proposed method can effectively predict GCM parameters, and the prediction precision in the three directions is 1.53, 1.08, and 3.46 mm, respectively.


2020 ◽  
Vol 14 (3) ◽  
pp. 295-302
Author(s):  
Chuandong Zhu ◽  
Wei Zhan ◽  
Jinzhao Liu ◽  
Ming Chen

AbstractThe mixture effect of the long-term variations is a main challenge in single channel singular spectrum analysis (SSA) for the reconstruction of the annual signal from GRACE data. In this paper, a nonlinear long-term variations deduction method is used to improve the accuracy of annual signal reconstructed from GRACE data using SSA. Our method can identify and eliminate the nonlinear long-term variations of the equivalent water height time series recovered from GRACE. Therefore the mixture effect of the long-term variations can be avoided in the annual modes of SSA. For the global terrestrial water recovered from GRACE, the peak to peak value of the annual signal is between 1.4 cm and 126.9 cm, with an average of 11.7 cm. After the long-term and the annual term have been deducted, the standard deviation of residual time series is between 0.9 cm and 9.9 cm, with an average of 2.1 cm. Compared with the traditional least squares fitting method, our method can reflect the dynamic change of the annual signal in global terrestrial water, more accurately with an uncertainty of between 0.3 cm and 2.9 cm.


Sign in / Sign up

Export Citation Format

Share Document