基于自适应多分辨率奇异值分解的大地电磁数据处理

针对强电磁干扰极易掩盖微弱的大地电磁有用信号,本文结合奇异值分解在去噪方面的优越性,提出基于自适应多分辨率奇异值分解(Adaptive Multi-Resolution Singular Value Decomposition, AMRSVD)的大地电磁数据处理方法.首先对大地电磁数据构建Hankel矩阵,利用MRSVD得到不同分辨率的近似信号和细节信号;然后选用近似信号和细节信号的标准差差值,对大地电磁数据进行信噪辨识;接着结合MRSVD和相邻细节信号的标准差差值,提出先验信息未知情况下的AMRSVD法;最后对辨识出的强干扰运用AMRSVD去除噪声,重构有用信号.实验结果表明,该方法的处理效率高,能有效分离出相关性较强的噪声,时间序列和视电阻率-相位曲线均得到有效改善.     

基于奇异谱分析的重磁位场分离方法

奇异谱分析是一种近年兴起的时间序列分析方法,它利用降秩原理实现信号分离.该方法将数据空间投影到不同特征的子空间中,并用奇异值来表征这些子空间的性质,最后通过截取奇异值实现数据的重构.重磁位场分离可以看成一种多信号叠加的分离问题.不同特征的重磁异常具有不同特征的奇异谱,这是奇异谱分析用于解决位场分离问题的应用基础.本文通过建立理论模型,分析重磁异常的奇异谱特征,得出适用于重磁位场分离的最优参数选择方法,并与传统方法进行比较.对比发现,无论是横向叠加模型、垂向叠加模型还是斜向叠加模型,奇异谱分析都具有很好的分离效果.最后,将奇异谱分析用于鄂东南某矿区的重力资料处理中,实现弱异常的识别和分离.     

分频段排干扰方法的建立与应用

在处理连续观测数据过程中,有时信号呈现出周期变化,如固体潮观测值中的潮汐信号等.直接对信号进行处理,如拟合、回归、调和分析等,是最常用的方法.当低频噪声背景和高频噪声均较强时,信号频率域(可称之为有用频段)内的信号被淹没了.直接提取信号较为困难.在这种情况下,本文试图利用分频段排干扰的方法,先将低频噪声和高频噪声分别排除掉,以便于进一步分析处理.本文实例计算表明,通过这样处理的应变固体潮观测资料,突出地显示了固体潮曲线的变化特征,而处理前的资料所对应的观测曲线,其固体潮汐的变化特征被噪声淹没了.     

SVD与EMD联合去噪方法在地震勘探数据处理中的研究与应用

将基于倾角扫描的奇异值分解与经验模式分解法相结合应用到地震资料随机噪声压制中。首先利用经验模式分解法消除部分噪声,增强地震道有效信号的相关性,再利用奇异值分解对地震信号进行同相轴自动追踪,截取小时窗数据体,并进行同相轴拉平处理,经SVD计算小时窗数据中心点的值来代替计算样点的值,最终实现随机噪声的压制。理论模型试算和实际资料处理表明,本文提出的EMD-SVD方法简单易行,比单一的SVD方法去噪效果更显著有效地消除了地震资料中的随机噪声,提高了地震资料的信噪比,并改善了叠加剖面的质量。     

基于压缩感知重构算法的大地电磁强干扰分离

为压制大地电磁信号中的强人文干扰,提出一种基于压缩感知重构算法的大地电磁信号去噪方法.通过构建与常见典型强干扰相匹配而对有用信号不敏感的冗余字典原子,利用改进的正交匹配追踪算法,分离出大地电磁信号中的强干扰成分.为了验证所述方法的强干扰分离效果,首先通过在实测大地电磁信号中加入理想的强干扰信号进行了仿真分离实验,然后从大量实测数据中选取三种含有不同类型强干扰的时间域片段,用所述方法对实测数据中的强干扰进行分离,最后将所述方法应用于青海试验点以及庐枞矿集区某测点实测数据的综合处理.仿真实验结果表明,该方法在分离出强干扰的同时,能够较好地保留有用信号.实测数据处理结果表明,该方法能够有效压制强干扰,改善强干扰区大地电磁数据的质量.     

改进的正交多项式变换在地震资料去噪中的应用

水平叠加虽然在很大程度上压制了噪声,提高了地震剖面的信噪比,但CMP遭集上还存在不少不是一次波的规则干扰和随机噪声,不利于叠前资料的岩性反演和叠后资料的波阻抗反演.本文提出了基于改进的正交多项式变换压制地震资料中随机噪声的方法,其优势在于:通过对不同时间信号的奇异值分解,确定有效信号正交多项式系数谱的阶数;再利用小波变换,改善有效信号和噪声在低阶上的混叠.文中给出了具体处理的过程,数据试验和实际资料的处理结果表明该方法不仅能有效地压制噪声,而且还能较好地保护地震数据中AVO变化特征.     

基于分数阶傅里叶变换的混合Cadzow滤波法

常规频率域SVD滤波法在随机噪声压制处理时,容易造成弯曲同相轴损伤。针对这一问题,本文提出一种基于分数阶傅里叶变换的混合Cadzow滤波法压制三维地震数据中随机噪声的方法。首先应用分数阶傅里叶变换,将地震数据变换到时频面,再依据Eigenimage滤波法与Cadzow滤波法建立混合的高维Hankel矩阵,然后对其运行奇异值分解,最后通过秩约化的方法来压制随机噪声。此方法用于四川某地区的地震数据处理,结果表明该方法可有效地去除随机噪声,保护有效信号,提高叠后地震数据的信噪比。     

基于信噪辨识的矿集区大地电磁噪声压制

为了避免形态滤波方法在大地电磁强干扰分离中的"过处理"、进一步保留大地电磁低频段的有用信息,提出基于信噪辨识的矿集区大地电磁噪声压制方法.首先,从信号处理的角度剖析矿集区典型强干扰与天然大地电磁微弱信号之间的定量辨识关系,利用形态分形维数和形态膨胀谱熵对大地电磁信号与强干扰进行信噪辨识.然后,结合形态滤波技术和阈值法,仅对辨识出明显不是天然大地电磁信号的异常波形进行噪声压制.最后,重构大地电磁有用信号,并对算法进行性能评价.仿真结果表明,形态分形维数和形态膨胀谱熵能较好地定量辨识大地电磁信号与强干扰,大地电磁信号中一些缓变化的低频信息得到了更为精细的保留;与形态滤波整体处理相比,本文所提方法获得的卡尼亚电阻率曲线更为光滑、连续,视电阻率值相对稳定,其结果更为真实地反映了测点本身所固有的大地电磁深部构造信息.     

地震随机噪声压缩感知迭代压制方法

复杂地表和复杂介质条件下,随机噪声往往严重影响着复杂地震信号的信噪比,同时深层地球物理目标探查中弱地震信号总是被随机噪声所掩盖,如何有效地压制随机噪声干扰、恢复有效地震信号仍然是高精度地震勘探中的关键问题.压缩感知理论突破了奈奎斯特采样定理的限制,利用有效地震信号的可压缩性和稀疏性,提供了从不可压缩随机噪声中进行有效信号分离的数据原理.本文系统分析压缩感知框架下地震随机噪声压制的稀疏优化反问题,提出了基于迭代软阈值算法的"采集-重建-修复"方案对该问题进行求解.在实现高度稀疏表征的基础上进行地震数据的压缩感知随机观测,通过迭代反演对有效地震信号进行重构,有效提高复杂地震数据的信噪比,同时,当求解稀疏优化问题时,如果出现正则化项引起重构信号衰减现象,可以匹配除偏对衰减的有效信号进行修复.通过与工业标准 f-x预测滤波方法进行比较,理论模型和实际数据处理的结果表明,压缩感知迭代噪声压制方法对复杂地震数据中的随机噪声有较好的压制效果,可以有效恢复出被较强非平稳随机噪声干扰的时空变同相轴信息.     

基于曲波变换的航空电磁数据调平方法研究

航空电磁法由于高效和高精度的特点广泛应用于地质填图、矿产资源、地下水、及环境与工程等勘查.然而,航空电磁系统处于动态环境,噪声影响严重,航空电磁数据处理至关重要.航空电磁数据噪声除随机成分外,还包括有各种效应引起的畸变,数据去噪需要依据噪声特征进行处理.航空电磁数据调平是航空电磁数据处理中至关重要的步骤,它能有效去除数据中由飞机飞行条件变化导致系统状态变化而产生的异常.传统的调平方法由于效率较低、易产生数据畸变等受到限制.为了克服这些局限性,我们提出一种基于曲波变换的数据调平方法.该方法得益于曲波变换多尺度和多方向性特征,可以有效地提取数据中的调平误差并予以去除.与此同时,利用该方法我们可以对非规则测区数据进行直接调平,无需进行测区分割,显著提高调平效率和普适性.为了检验本文曲波变换调平方法的有效性,我们将其应用于理论数据以及在爱尔兰Waterford地区实测的航电数据调平.实验结果表明该方法有效地去除调平误差的同时很好地保留有用信号.     

Noise suppression for microseismic data by non‐subsampled shearlet transform based on singular value decomposition

The existence of strong random noise in surface microseismic data may decrease the utility of these data. Non‐subsampled shearlet transform can effectively suppress noise by properly setting a threshold to the non‐subsampled shearlet transform coefficients. However, when the signal‐to‐noise ratio of data is low, the coefficients related to the noise are very close to the coefficients associated with signals in the non‐subsampled shearlet transform domain that the coefficients related to the noise will be retained and be treated as signals. Therefore, we need to minimise the overlapping coefficients before thresholding. In this paper, a singular value decomposition algorithm is introduced to the non‐subsampled shearlet transform coefficients, and low‐rank approximation reconstructs each non‐subsampled shearlet transform coefficient matrix in the singular value decomposition domain. The non‐subsampled shearlet transform coefficients of signals have bigger singular values than those of the random noise, which implies that the non‐subsampled shearlet transform coefficients can be well estimated by taking only a few largest singular values. Therefore, those properties of singular value decomposition may significantly help minimise overlapping of noise and signals coefficients in the non‐subsampled shearlet transform domain. Finally, the denoised microseismic data are obtained easily by giving a simple threshold to the reconstructed coefficient matrix. The performance of the proposed method is evaluated on both synthetic and field microseismic data. The experimental results illustrate that the proposed method can eliminate random noise and preserve signals of interest more effectively.     

瞬变电磁波场转换算法的改进

地球物理电磁场数据与虚拟地震波场数据之间存在数学上的等效转换关系,通过这种等效转换,可有效提高地球物理电磁法对地下目标体分界面的辨识度.但是这种转换在数学上属于不适定问题,可采用奇异值分解法处理.由于大奇异值控制计算矩阵的主要信息,小的奇异值控制计算矩阵的次要信息,传统的截断奇异值分解法只保留大奇异值,而忽略小的奇异值,导致数值解不够精确.本文提出一种新的修正方案——改进截断奇异值法,采用岭估计方法计算由小的奇异值引起的虚拟波场.模型计算结果表明：改进截断奇异值法比传统的奇异值分解法得到的波场转换结果更好,对某煤矿采空区探测数据进行了处理,成功分辨出采空区分界面.     

Random noise attenuation via the randomized canonical polyadic decomposition

Tensor algebra provides a robust framework for multi-dimensional seismic data processing. A low-rank tensor can represent a noise-free seismic data volume. Additive random noise will increase the rank of the tensor. Hence, tensor rank-reduction techniques can be used to filter random noise. Our filtering method adopts the Candecomp/Parafac decomposition to approximates a N-dimensional seismic data volume via the superposition of rank-one tensors. Similar to the singular value decomposition for matrices, a low-rank Candecomp/Parafac decomposition can capture the signal and exclude random noise in situations where a low-rank tensor can represent the ideal noise-free seismic volume. The alternating least squares method is adopted to compute the Candecomp/Parafac decomposition with a provided target rank. This method involves solving a series of highly over-determined linear least-squares subproblems. To improve the efficiency of the alternating least squares algorithm, we uniformly randomly sample equations of the linear least-squares subproblems to reduce the size of the problem significantly. The computational overhead is further reduced by avoiding unfolding and folding large dense tensors. We investigate the applicability of the randomized Candecomp/Parafac decomposition for incoherent noise attenuation via experiments conducted on a synthetic dataset and field data seismic volumes. We also compare the proposed algorithm (randomized Candecomp/Parafac decomposition) against multi-dimensional singular spectrum analysis and classical prediction filtering. We conclude the proposed approach can achieve slightly better denoising performance in terms of signal-to-noise ratio enhancement than traditional methods, but with a less computational cost.     

基于经验模态分解的分数维地震随机噪声衰减方法

经验模态分解算法(EMD)是一种基于有效波和噪声尺度差异进行波场分离的随机噪声压制方法,但由于实际地震数据波场复杂,导致模态混叠较严重,仅凭该方法进行去噪很难达到理想效果.本文基于EMD算法对信号多尺度的分解特性,结合Hausdorff维数约束条件,提出一种用于地震随机噪声衰减的新方法.首先对地震数据进行EMD自适应分解,得到一系列具有不同尺度的、分形自相似性的固有模态分量(IMF);在此基础上,基于有效信号和随机噪声的Hausdorff维数差异,识别混有随机噪声的IMF分量,对该分量进行相关的阈值滤波处理,从而实现有效信号和随机噪声的有效分离.文中从仿真信号试验出发,到模型地震数据和实际地震数据的测试处理,同时与传统的EMD处理结果相对比.结果表明,本文方法对地震随机噪声的衰减有更佳的压制效果.     

Cancelling directional EM noise in magnetoteilurics1

The prospecting of densely urbanized areas by the measurement of magnetic and electric natural fields is severely hampered by electromagnetic (EM) noise. Active man-made EM noise sources can generally be considered fixed in space, thus affecting the magnetotelluric (MT) signals of a measuring site mainly along their polarization directions. Taking advantage of the impulsive nature of polarized EM noise, a time-domain directional noise cancelling (DNC) technique is proposed. The comparison of noisy data with data predicted, using a low noise reference signal or with data interpolated whenever no reference is available, allows the detection of the most likely noise sources with prevailing directional patterns using a Bayes's criterion. The DNC approach is general and can be adapted, depending on the reference signal used (single-site or remote-reference). In field data, hodograms of the prediction residuals basically confirm the directional noise model assumed in DNC. An example is presented in which the DNC technique has been applied to a single-site MT survey carried out in northern Italy, where the signal was heavily corrupted by noise with prevailing directional properties due to the densely urbanized area. MT apparent resistivities and phases obtained at the site of the survey before and after DNC are presented and discussed.     

3D targeted multiple attenuation

Short-period multiple reflections pose a particular problem in the North Sea where predictive deconvolution is often only partially successful. The targeted multiple attenuation (TMA) algorithm comprises computation of the covariance matrix of preflattened prestack or post-stack seismic data, the determination of the dominating eigenvectors of the covariance matrix, and subtraction of the related eigenimages followed by reverse flattening. The main assumption made is that the flattened multiple reflections may be represented by the first eigenimage(s) which implies that the spatial amplitude variations of primaries and associated multiples are similar. This assumption usually limits the method to short-period multiple reflections. TMA is applicable post-stack or prestack to common-offset gathers. It is computationally fast, robust towards random noise, irregular geometry and spatial aliasing, and it preserves the amplitudes of primaries provided they are not parallel to the targeted multiples. Application of TMA to 3D wavefields is preferable because this allows a better discrimination between primaries and multiples. Real data examples show that the danger of partially removing primary energy can be reduced by improving the raw multiple model that is based on eigenimages, for example by prediction filtering.     

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

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

4D signal enhancement using singular-value decomposition: application to mapping oil–water contact movement across the Nelson Field

A new method for time‐lapse signal separation and enhancement using singular‐value decomposition is presented. Singular‐value decomposition is used to separate a 4D signal into its constituent parts: common geology, time‐lapse response and noise. Synthetic tests which demonstrate the advantages of the singular‐value decomposition technique over traditional differencing methods are also presented. This signal separation and enhancement technique is used to map out both the original and moved oil–water contacts across the Nelson Field. The singular‐value decomposition technique allows the oil–water contact to be mapped across regions which would have been missed using traditional differencing methods. In particular, areas toward the edges of the field are highlighted by the technique. The oil–water contact is observed to move upwards across the field, with the largest movements being associated, as anticipated, with natural production. The results obtained are broadly consistent with those predicted by the reservoir simulator model. Singular‐value decomposition is demonstrated to be a useful tool for enhancing the time‐lapse signal and for gaining confidence in areas where traditional differencing fails.     

Dip‐adaptive singular‐value decomposition filtering for seismic reflection enhancement

We present a singular value decomposition (SVD) filtering method for the enhancement of coherent reflections and for attenuation of noise. The method is applied in two steps. First normal move‐out (NMO) correction is applied to shot or CMP records, with the purpose of flattening the reflections. We use a spatial SVD filter with a short sliding window to enhance coherent horizontal events. Then the data are sorted in common‐offset panels and the local dip is estimated for each panel. The next SVD filtering is performed on a small number of traces and a small number of time samples centred around the output sample position. Data in a local window are corrected for linear moveout corresponding to the dips before SVD. At the central time sample position, we sum over the dominant eigenimages of a few traces, corresponding to SVD dip filtering. We illustrate the method using land seismic data from the Tacutu basin, located in the north‐east of Brazil. The results show that the proposed method is effective and is able to reveal reflections masked by ground‐roll and other types of noise.