基于多用户峰度准则的海洋强噪声衰减方法(英文)

海洋地震勘探过程中,由于采集设备的老化或电源的不稳定而造成的漏电,在地震记录表现为强噪音干扰,利用常规噪音衰减方法处理此类强噪音效果并不理想。鉴于强噪音在统计学上具有相同的特性,本文在基于峰度的盲分离(blind source separation,BSS)算法研究基础上,推导出一种基于多用户峰度(multiuser kurtosis,MUK)准则的噪音衰减算法来估计地震记录中具有相同统计特性的强噪音,并将其从地震记录中分离,从而达到衰减强噪音的目的。模型试验与实际资料的处理表明:该方法能够在好的衰减海洋地震勘探记录中的强噪音,保留了更多的有效信息,提高海洋地震数据的信噪比,具有可行性和应用前景。     

基于深度卷积神经网络的地震数据随机噪声压制

本文以提高地震数据的成像质量为目标,提出一种智能的卷积神经网络降噪框架,从带有噪声的地震数据中自适应地学习地震信号。为了加速网络训练和避免训练时出现梯度消失现象,我们在网络中加入残差学习和批标准化的方法,并采用了ReLU激活函数和Adam优化算法优化网络。此外,Marmousi和F3数据集被用来对网络进行训练和测试,经过充分训练的网络不仅能在学习中保留地震数据特征,而且能去除随机噪声。首先充分地训练网络,从中提取出随机噪声,并保留学习到的地震数据特征,之后通过重建地震数据估算测试集中的波形特征。合成记录和实际数据的处理结果显示了深度卷积神经网络在随机噪声压制任务中的潜力,并通过实验验证表明了深度卷积神经网络框架有很好的去噪效果。      

一种神经网络改进小波的地震数据随机噪声去除方法

地震资料的有效信号反射弱,且易受多次波的影响,不可避免地存在随机噪声干扰。提出一种基于神经网络改进小波的地震数据随机噪声去除方法,采用神经网络模型,识别出随机噪声信号,对该信号进行小波包分解,获取多类别随机噪声信号,采用级联BP神经网络模型提取出多类别随机噪声信号,实现地震数据的随机信号压制。实验结果显示,这种改进小波方法对地震数据随机噪声信号的去噪效果较好,在复杂沉积地质结构被探测介质的地震数据随机噪声压制方面具有较强的适用性。     

基于数据增广和CNN的地震随机噪声压制

卷积神经网络(Convolutional Neural Network,CNN)是一种基于数据驱动的学习算法,简化了传统从特征提取到分类的两阶段式处理任务,被广泛应用于计算机科学的各个领域.在标注数据不足的地震数据去噪领域,CNN的推广应用受到限制.针对这一问题,本文提出了一种基于数据生成和增广的地震数据CNN去噪框架.对于合成数据,本文对无噪地震数据添加不同方差的高斯噪声,增广后构成训练集,实现基于小样本的CNN训练.对于实际地震数据,由于无法获得真实的干净数据和噪声来生成训练样本集,本文提出一种直接从无标签实际有噪数据生成标签数据集的方法.在所提出的方法中,我们利用目前已有的去噪方法从实际地震数据中分别获得估计干净数据和估计噪声,前者与未知的干净数据具有相似纹理,后者与实际噪声具有相似的概率分布.人工合成数据和实际数据实验结果表明,相较于F-X反褶积,BM3D和自适应频域滤波算法,本文方法能更好地压制随机噪声和保护有效信号.最后,本文采用神经网络可视化方法对去噪CNN的机理进行了探索,一定程度上解释了网络每一层的学习内容.     

Seismic ground-roll noise attenuation using deep learning

We propose to adopt a deep learning based framework using generative adversarial networks for ground-roll attenuation in land seismic data. Accounting for the non-stationary properties of seismic data and the associated ground-roll noise, we create training labels using local time–frequency transform and regularized non-stationary regression. The basic idea is to train the network using a few shot gathers such that the network can learn the weights associated with noise attenuation for the training shot gathers. We then apply the learned weights to test ground-roll attenuation on shot gathers, that are not a part of training input to obtain the desired signal. This approach gives results similar to local time–frequency transform and regularized non-stationary regression but at a significantly reduced computational cost. The proposed approach automates the ground-roll attenuation process without requiring any manual input in picking the parameters for each shot gather other than in the training data. Tests on field-data examples verify the effectiveness of the proposed approach.     

基于Seislet-TV双正则化约束的地震随机噪声压制方法

人工地震数据总是受到随机噪声的干扰,地震数据时-空变的特性使得常规去噪方法处理效果并不理想,容易导致有效信号的损失.目前广泛应用的预测滤波类方法存在处理时变数据能力不足的问题.随着压缩感知理论的不断完善,稀疏变换阈值算法能够解决时变地震数据噪声压制问题,但是常规的稀疏变换方法,如傅里叶变换,小波变换等,并不是特殊针对地震数据设计的,很难提供地震数据最佳的压缩特征,同时,常规阈值算法容易导致去噪结果过于平滑.因此开发更加有效的时-空变地震数据信噪分离方法具有重要的工业价值.本文将地震数据信噪分离问题归纳为数学基追踪问题,在压缩感知理论框架下,利用特殊针对地震数据设计的VD-seislet稀疏变换方法,结合全变差(TV)算法,构建seislet-TV双正则化条件,并利用分裂Bregman迭代算法求解约束最优化问题,实现地震数据的有效信噪分离.通过理论模型和实际数据测试本文方法,并且与工业标准FXdecon方法进行比较,结果表明基于seislet-TV双正则化约束条件的迭代方法能够更加有效地保护时-空变地震信号,压制地震数据中的强随机噪声.     

Random noise suppression for seismic data using a non-local Bayes algorithm

For random noise suppression of seismic data, we present a non-local Bayes (NLBayes) filtering algorithm. The NL-Bayes algorithm uses the Gaussian model instead of the weighted average of all similar patches in the NL-means algorithm to reduce the fuzzy of structural details, thereby improving the denoising performance. In the denoising process of seismic data, the size and the number of patches in the Gaussian model are adaptively calculated according to the standard deviation of noise. The NL-Bayes algorithm requires two iterations to complete seismic data denoising, but the second iteration makes use of denoised seismic data from the first iteration to calculate the better mean and covariance of the patch Gaussian model for improving the similarity of patches and achieving the purpose of denoising. Tests with synthetic and real data sets demonstrate that the NL-Bayes algorithm can effectively improve the SNR and preserve the fidelity of seismic data.     

Sub-basalt imaging problems and the application of Artificial Neural Networks

In many areas of the world, the presence of shallow high velocity, highly heterogeneous layers complicate seismic imaging of deeper reflectors. Of particular economic interest are areas where potentially hydrocarbon-bearing strata are obscured by layers of basalt. Basalt layers are highly reflective and heterogeneous. Using reflection seismic, top basalt is typified by a high-amplitude, coherent reflector with poor resolution of reflectors below the basalt, and even bottom basalt. Here, we present a new approach to the imaging problem using the pattern recognition abilities of a back-propagation Artificial Neural Network (ANN). ANNs are computational systems that attempt to mimic natural biological neural networks. They have the ability to recognize patterns and develop their own generalizations about a given data set. Back-propagation neural networks are trained on data sets for which the solution is known and tested on the data that are not previously presented to the ANN in order to validate the network result. We show that Artificial Neural Networks, due to their pattern recognition capabilities, can invert the medium statistics based on the seismic character. We produce statistically defined models involving a basalt analogous layer, and calculate full wavefield finite difference synthetic seismograms. We vary basalt layer thickness and source frequency to generate a synthetic model that produces seismic that is similar to real sub-basalt seismic, i.e. high amplitude top basalt reflector and the absence of base basalt and sub-basalt events. Using synthetic shot gathers, generated in a synthetic representation of the sub-basalt case, we can invert the velocity medium standard deviation by using an ANN. By inverting the velocity medium standard deviation, we successfully identified the transition from basalt to sub-basalt on the synthetic shot gathers. We also show that ANNs are capable of identifying the basalt to sub-basalt transition in the presence of incoherent noise. This is important for any future applications of this technique to the real-world seismic data, as this data is never completely noise-free. There is always a certain level of residual (noise remaining after initial noise filtering) environmental/ambient noise present on the recorded seismics, hence, neural network training with noise-free synthetic seismic is less than optimal.     

时频峰值滤波去噪技术及其应用

本文将时频峰值滤波(TFPF)去噪技术应用于共炮点地震资料的随机噪声压制.时频峰值滤波技术是通过频率调制将信号调制成解析信号的瞬时频率,利用解析信号的Wigner-Ville分布的峰值进行瞬时频率估计,恢复有效信号,与其它去噪方法相比,TFPF具有在较少的约束条件下压制强随机噪声的优点.本文针对实际地震资料的非线性特性,利用加窗的Wigner-Ville分布实现TFPF,使得地震信号在一个窗长内近似满足线性瞬时频率条件,减小由地震信号非线性引起的偏差.本文对共炮点地震记录做时频峰值滤波处理,滤波结果表明在地震勘探资料中存在强随机噪声的情况下,利用局部线性化处理的时频峰值滤波技术可以有效地压制地震资料中的随机噪声,恢复出湮没在随机噪声中的地震反射信号.信噪比提高3～6 dB.     

A seismic interpolation and denoising method with curvelet transform matching filter

A new seismic interpolation and denoising method with a curvelet transform matching filter, employing the fast iterative shrinkage thresholding algorithm (FISTA), is proposed. The approach treats the matching filter, seismic interpolation, and denoising all as the same inverse problem using an inversion iteration algorithm. The curvelet transform has a high sparseness and is useful for separating signal from noise, meaning that it can accurately solve the matching problem using FISTA. When applying the new method to a synthetic noisy data sets and a data sets with missing traces, the optimum matching result is obtained, noise is greatly suppressed, missing seismic data are filled by interpolation, and the waveform is highly consistent. We then verified the method by applying it to real data, yielding satisfactory results. The results show that the method can reconstruct missing traces in the case of low SNR (signal-to-noise ratio). The above three problems can be simultaneously solved via FISTA algorithm, and it will not only increase the processing efficiency but also improve SNR of the seismic data.     

基于分块字典学习理论的地震数据去噪

随着油气勘探观测环境愈发复杂,采集的地震数据常常掺杂各种噪声信号,导致勘探目标引起的有效微弱信号被覆盖,严重影响高精度的地震勘探数据解译,因而有效的压制地震勘探数据噪声显得越发重要。本文采用字典学习策略,将复杂地震数据进行分块,通过分块数据的字典学习获取字典原子,构建高精度的字典学习地震数据稀疏表示,通过两次迭代更新字典原子,进行数据去噪。将本文的字典学习算法应用于含随机噪声的模拟数据和实测地震勘探数据处理,验证该算法的可行性及有效性。结果表明,本文算法有效去除了随机噪声,保留了有效信号同相轴,提高了信噪比,可为复杂含噪地震数据的去噪处理提供新的技术手段。     

Microseismic events enhancement and detection in sensor arrays using autocorrelation‐based filtering

Passive microseismic data are commonly buried in noise, which presents a significant challenge for signal detection and recovery. For recordings from a surface sensor array where each trace contains a time‐delayed arrival from the event, we propose an autocorrelation‐based stacking method that designs a denoising filter from all the traces, as well as a multi‐channel detection scheme. This approach circumvents the issue of time aligning the traces prior to stacking because every trace's autocorrelation is centred at zero in the lag domain. The effect of white noise is concentrated near zero lag; thus, the filter design requires a predictable adjustment of the zero‐lag value. Truncation of the autocorrelation is employed to smooth the impulse response of the denoising filter. In order to extend the applicability of the algorithm, we also propose a noise prewhitening scheme that addresses cases with coloured noise. The simplicity and robustness of this method are validated with synthetic and real seismic traces.     

基于全变分原理的多震源混合数据直接偏移方法

多震源混合地震采集技术,即将多个震源以一定编码方式连续地激发,得到多炮混合的地震数据.该技术能减少地震采集时间,节约采集成本,但是混合数据的直接偏移会在成像剖面中引入严重的串扰噪声,影响成像效果.从数学上看,地震成像属于典型的数学物理反问题,可以采用线性反演方法求解一个正则化约束的最小二乘(LS)优化问题,获得更高质量的成像结果.全变分(TV)正则化方法是图像去噪和复原领域中广泛应用的热点技术,其能在去除噪声的过程中保留图像的边缘信息和不连续性.在对TV图像去噪复原方法原理分析的基础上,本文将多震源混合数据直接偏移成像问题转换成图像复原的极小化能量泛函问题,用TV正则化代替传统最小二乘偏移(LSM)中的L2范数正则化,提出基于全变分原理的混合数据直接偏移方法.该方法使用基于梯度的快速迭代收缩阈值与快速梯度投影组合算法——FISTA/FGP求解最优化问题,能有效压制串扰噪声,增强同相轴连续性,提高成像分辨率.理论模型测试结果表明:将本方法应用于混合数据,无论是去噪效果还是成像精度都得到显著改善.     

最小二乘反演迭代算法在压制海上变深度缆采集数据虚反射中的应用

海水面的虚反射(鬼波)引起海上拖缆采集数据陷波,导致地震记录频带变窄,而近年发展的变深度缆采集技术,具有多样的陷波特征,通过专门的去虚反射处理方法可获得宽频数据.本文基于已有研究成果,将最小二乘反演迭代压制虚反射算法应用于某海上变深度缆宽频处理.基于频率波数域镜像记录生成方法获得镜像炮集记录,并采用最小二乘解从变深度缆原始和镜像炮集记录中提取上行波.针对镜像炮集记录生成受初始速度模型精度的影响,使得某深度缆接收的上行波和下行波之间的实际延迟时间存在误差,采用最小二乘反演迭代算法最优化计算下行波与上行波之间的平均延迟时间和上行波记录,并采用时空数据窗口滑动克服延迟时间随炮检距和目的层深度变化问题.合成数据及某海上实际变深度缆数据处理测试结果表明,该方法能较好地压制变深度缆由海水面产生的虚反射,能达到拓宽地震记录频带目的.     

Local outlier factor as part of a workflow for detecting and attenuating blending noise in simultaneously acquired data

A number of deblending methods and workflows have been reported in the past decades to eliminate the source interference noise recorded during a simultaneous shooting acquisition. It is common that denoising algorithms focusing on optimizing coherency and weighting down/ignoring outliers can be considered as deblending tools. Such algorithms are not only enforcing coherency but also handling outliers either explicitly or implicitly. In this paper, we present a novel approach based on detecting amplitude outliers and its application on deblending based on a local outlier factor that assigns an outlier-ness (i.e. a degree of being an outlier) to each sample of the data. A local outlier factor algorithm quantifies outlier-ness for an object in a data set based on the degree of isolation compared with its locally neighbouring objects. Assuming that the seismic pre-stack data acquired by simultaneous shooting are composed of a set of non-outliers and outliers, the local outlier factor algorithm evaluates the outlier-ness of each object. Therefore, we can separate the data set into blending noise (i.e. outlier) and signal (i.e. non-outlier) components. By applying a proper threshold, objects having high local outlier factors are labelled as outlier/blending noise, and the corresponding data sample could be replaced by zero or a statistically adequate value. Beginning with an explanation of parameter definitions and properties of local outlier factor, we investigate the feasibility of a local outlier factor application on seismic deblending by analysing the parameters of local outlier factor and suggesting specific deblending strategies. Field data examples recorded during simultaneous shooting acquisition show that the local outlier factor algorithm combined with a thresholding can detect and attenuate blending noise. Although the local outlier factor application on deblending shows a few shortcomings, it is consequently noted that the local outlier factor application in this paper obviously achieves benefits in terms of detecting and attenuating blending noise and paves the way for further geophysical applications.     

Detection and identification of seismic signals recorded at Krakatau volcano (Indonesia) using artificial neural networks

The Anak Krakatau volcano (Indonesia) has been monitored by a multi-parametric system since 2005. A variety of signal types can be observed in the records of the seismic stations installed on the island volcano. These include volcano-induced signals such as LP, VT, and tremor-type events as well as signals not originating from the volcano such as regional tectonic earthquakes and transient noise signals. The work presented here aims at the realization of a system that automatically detects and identifies the signals in order to estimate and monitor current activity states of the volcano. An artificial neural network approach was chosen for the identification task. A set of parameters was defined, describing waveform and spectrogram properties of events detected by an amplitude-ratio-based (STA/LTA) algorithm. The parameters are fed into a neural network which is, after a training phase, able to generalize input data and identify corresponding event types. The success of the identification depends on the network architecture and training strategy. Several tests have been performed in order to determine appropriate network layout and training for the given problem. The performance of the final system is found to be well suited to get an overview of the seismic activity recorded at the volcano. The reliability of the network classifier, as well as general drawbacks of the methods used, are discussed.     

自适应非局部均值地震随机噪声压制(英文)

非局部均值滤波是一种基于图像信息冗余的去噪方法,其认为图像自身的有效结构具有一定的重复性,而随机噪声则不具备这一特点,通过利用图像本身的自相似性来达到压制随机噪声的目的,是一种全局的去噪方法。本文把这一思想引入地震数据随机噪声压制中,针对传统非局部均值滤波计算量过大的问题,文章采用分块非局部均值的方式来减少计算量;针对滤波参数选取会影响非局部均值滤波效果的问题,提出一种简单的自适应滤波参数地震数据分块非局部均值算法。模型和实际数据处理结果表明:相对于传统的去噪算法(如f-x反褶积),该方法在压制随机噪声的同时对有效信号保护地更好,具有更高的保真度,更有利于后续的处理和解释工作。     

压制空间非平稳地震勘探随机噪声的ROAD径向时频峰值滤波算法

径向时频峰值滤波算法是一种有效保持低信噪比地震勘探记录中反射同相轴的随机噪声压制方法,但该算法对空间非平稳地震勘探随机噪声压制效果不理想.本文研究空间非平稳地震勘探随机噪声,即各道噪声功率不同的地震勘探随机噪声,其在径向滤波轨线上表征近似脉冲噪声,在径向时频峰值滤波过程中干扰相邻道滤波结果.为了减小空间非平稳随机噪声的影响,本文提出一种基于绝对级差统计量(ROAD)的径向时频峰值滤波随机噪声压制方法.该方法首先根据径向轨线上信号的绝对级差统计量检测空间非平稳地震勘探随机噪声,然后结合局部时频峰值滤波和径向时频峰值滤波压制地震勘探记录中的随机噪声.将ROAD径向时频峰值滤波方法应用于合成记录和实际共炮点地震记录,结果表明ROAD径向时频峰值滤波方法可以压制空间非平稳地震勘探随机噪声且不损害有效信号,有效抑制随机噪声空间非平稳对滤波结果的影响.与径向时频峰值滤波相比,ROAD径向时频峰值滤波方法更适用于空间非平稳地震勘探随机噪声压制.     

Improving Airgun Signal Detection with Small-Aperture Seismic Array in Yunnan

Repeating airgun sources are eco-friendly sources for monitoring the changes in the physical properties of subsurface mediums, but their signals decay quickly and are buried in the noises soon after traveling short distances. Stacking waveforms from different airgun shots recorded by a single seismic station (shot stacking) is the most popular technique to detect weak signals from noisy backgrounds, and has been widely used to process the data of Fixed Airgun Signal Transmission Stations (FASTS) in China. However, shot stacking sacrifices the time resolution in monitoring to recover a qualified airgun signal by stacking many shots at distance stations, and also suffers from persistent local noises. In this paper, we carried out several small-aperture seismic array experiments around the Binchuan FAST Station (BCFASTS) in Yunnan Province,China, and applied the array technique to improve airgun signal detection. The results show that seismic array processing combining with shot stacking can suppress seismic noises more efficiently, and provide better signal-to-noise ratio (SNR) and coherent airgun signals with less airgun shots. This work suggests that the array technique is a feasible and promising tool in FAST to increase the time resolution and reduce noise interference on routine monitoring.     

低信噪比地震资料FDOC-seislet变换阈值消噪方法研究

地震数据本质上是时变的,不仅有效同相轴表现出确定性信号的时变特征,而且复杂地表和构造条件以及深部探测环境总是引入时变的非平稳随机噪声.标准的频率-空间域预测滤波只适合压制平面波信号假设下的平稳随机噪声,而处理非平稳地震随机噪声时,需要将数据体分割为小窗口进行分析,但效果不够理想,而传统非预测类随机噪声压制方法往往适应性不高,因此开发能够保护地震信号时变特征的随机噪声压制方法具有重要的工业价值.压缩感知是近年出现的一个新的采样理论,通过开发信号的稀疏特性,已经在地震数据处理中的数据插值以及噪声压制中得到了应用.本文系统地分析了压缩感知理论框架下的地震随机噪声压制问题,建立了阈值消噪的数学反演目标函数;针对时变有效信息具有的可压缩性,利用有限差分算法求解炮检距连续方程,构建有限差分炮检距连续预测算子(FDOC),在seislet变换框架下,提出一种新的快速稀疏变换域———FDOC-seislet变换,实现地震数据的高度稀疏表征;结合非平稳随机噪声不可压缩的特征,提出了一种整形迭代消噪方法,该方法是一种广义的迭代收缩阈值(IST)算法,在无法计算稀疏变换伴随算子的条件下,仍然能够对强噪声环境中的时变有效信息进行有效恢复.通过对模型数据和实际数据的处理,验证了FDOC-seislet稀疏变换域随机噪声迭代压制方法能够在保护复杂构造地震波信息的前提下,有效地衰减原始数据中的强振幅随机噪声干扰.