Double sparsity dictionary for seismic noise attenuation |

We train the filters patch by patch of the whole dataset. We apply the maximally overlapping patches (i.e., the two neighbor patches only shifted by one column) for the training. This creates a richness in the training data and leads to a shift-invariance for the dictionary. In the DDTF, for each patch, the training filer is orthogonal, but the full transform matrix (each column is a training filter on patches) is redundant. Because of the base seislet transform, the DDTF trains the dictionary from a multi-scale seislet domain instead of the single scale domain of the original data. Although we follow the same strategy as Ophir et al. (2011), we utilize different tools for the two main components in the multi-scale learning framework. The DDTF used in the proposed method is much faster than the K-SVD method used by Ophir et al. (2011), and the seislet transform used in the proposed method is more appropriate for seismic data than the wavelets used by Ophir et al. (2011). In K-SVD, there exists a large number of patches for training and the components of the dictionary are updating individually (indicating many SVD decompositions). However, the DDTF method updates all the components (columns of ) by one SVD decomposition (Yu et al., 2015), thus is more efficient. Besides, because of the learning efficiency in DDTF, the DDTF based approaches can be applicable to high-dimensional denoising and interpolation problems. The seislet transform, however, is superior than traditional wavelet transform for compressing seismic data because it utilizes local slope to construct the transform. The superior performance of seislet transform over state-of-the-art transforms in denoising, interpolation, and deblending of seismic data have been demonstrated in a number of references (Gan et al., 2015; Chen et al., 2014; Fomel and Liu, 2010; Liu and Fomel, 2010). Thus, the combination between the seislet transform and DDTF through the DSD framework seems very appealing.

The cascaded DSD framework, combined with thresholding-based denoising, can be briefly summarized as follows:

- Apply the seislet transform to noisy seismic data.
- Apply the DDTF to each band in the seislet domain to form DSD domain.
- Perform thresholding in the DSD domain.
- Apply the adjoint DDTF to the thresholded data in DSD domain.
- Apply the inverse seislet transform to recover the denoised data.

There is an inner iterative process in the above steps when applying the DDTF. It refers to applying the DDTF by the data-driven learning process following equations 6 to 8. The iterative process acts as a training step for the DDTF, and the trained dictionary are finally used to transform the seislet subband data into the double-sparsity domain and used for the adjoint DDTF step in step 4. Figure 1 shows a demonstration of denoising processes of the proposed DSD. It also shows the seislet transform domain and DSD transform domain before and after thresholding.

demo2
Demonstration of the DSD-based denoising processes.
Figure 1. |
---|

Double sparsity dictionary for seismic noise attenuation |

2016-02-27