Nevertheless, is actually still contaminated with residual random noise. Huang et al. (2016) suggest using damped MSSA (DMSSA) algorithm to attenuate the residual random noise that the conventional MSSA fails to. Here, we further extend the DMSSA algorithm to simultaneous reconstruction and denoising of 3D seismic data. Since the missing data at each frequency slice can be regarded as the random noise, the derivation are similar with Huang et al. (2016) except for the extra reconstruction process based on the weighted POCS-like method. Therefore, we conclude the approximation of as:
The revised TSVD for DMSSA algorithm can be represented in operator notation as follows:
The filtered data is recovered with random noise attenuated and missing traces reconstructed by properly averaging along the anti-diagonals of the low-rank matrix
obtained via revised TSVD in DMSSA:
In this paper, we pay our attention to the reconstruction of noisy data, more specifically, those noisy data with very low SNR. Here, we recall the traditional weighted POCS-like algorithm (Oropeza and Sacchi, 2011) as follows:
Substituting the MSSA filter with the DMSSA filter , the improved DMSSA algorithm for simultaneous seismic data reconstruction and denoising can be formulated as follows:
The MSSA based approach is based on the assumption that seismic data is of low rank at each frequency slice in the frequency-space domain. Since seismic data are locally composed of plane-wave components, it intuitively satisfies such assumption. According to our experience, when the SNR is relatively high, in other words with less noise, the performance of MSSA and DMSSA based approaches is similar to each other. However, when the SNR is very low, the proposed DMSSA based approach will outperform the MSSA based approach a lot.