Simultaneous seismic data denoising and reconstruction via multichannel singular spectrum analysis
Simultaneous seismic data denoising and reconstruction via multichannel singular spectrum analysis
复制标题
DOI：
10.1190/1.3552706
复制
发表时间：
2011-05-01
期刊：
影响因子：
3.3
通讯作者：
Sacchi, Mauricio
中科院分区：
文献类型：
--
作者：
Oropeza, Vicente;Sacchi, Mauricio
我们提出了一个降秩算法，允许同时重建和随机噪声衰减的地震记录。我们基于我们的技术多通道奇异谱分析（MSSA）。该技术需要将给定时间频率的空间数据组织成块汉克尔矩阵，该块汉克尔矩阵在理想条件下是秩为k的矩阵，其中k是分析窗口中平面波的数量。加性噪声和缺失样本将增加数据的块汉克尔矩阵的秩。因此，降秩被提出作为衰减噪声和恢复丢失的道的手段。本文提出了一种类似于凸集投影法的地震数据重建迭代算法。此外，我们建议采用随机奇异值分解来加速算法的降秩阶段。我们应用MSSA重建合成的例子和现场数据集。合成的例子被用来评估该方法在两个重建方案的性能：一个无噪声的情况下，数据污染的噪声。在这两种情况下，我们发现极低的重建错误，这表明最佳的恢复。野外数据示例由依赖于公共中点和偏移的二维叠前体组成。我们使用MSSA重建方法来完成缺失的偏移距，同时提高地震体的信噪比。
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 k, where k 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.