基于二维希尔伯特变换的地震倾角求取方法及其在构造导向滤波中的应用（英文）

在地震勘探数据采集中,随机噪声严重影响地震资料质量,给后期解释工作带来很大困难。如何在不损失剖面有效信息的前提下压制随机噪声,有效地提高地震资料的信噪比和保真度,是本文的研究目标。构造导向滤波技术的核心是构造方向表征的求取以及如何实现非平稳滤波,来达到提高地震数据信噪比和保真度的目的。本文首先通过分析函数二维导数与希尔伯特变换的频率响应关系,推导出了基于二维希尔伯特变换的非迭代地震同相轴倾角求取算子,进而达到了构造方向表征的求取;其次选取多项式拟合作为构造导向滤波中的非平稳滤波方法,扩展了非平稳多项式拟合的应用范围;最后沿构造倾角方向进行变振幅同相轴的非平稳多项式拟合,实现和构建了新的自适应构造导向滤波方法。理论模型和实际地震资料处理的结果表明,所提出的方法实现了既保护构造信息又有效地压制了随机噪声的目的。     

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

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

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

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

基于SW统计量的自适应时频峰值滤波压制地震勘探随机噪声研究

由于金属矿区地震记录中随机噪声性质复杂且信噪比低,常规降噪方法难以达到预期的滤波效果.时频峰值滤波(TFPF)方法是实现低信噪比地震勘探记录中随机噪声压制的有效方法,但其在复杂地震勘探随机噪声下时窗参数优化问题仍难以解决.本文充分利用地震勘探噪声的统计特性,结合Shapiro-Wilk(SW)统计量辨识地震勘探记录中的微弱有效信号,提出基于SW统计量的自适应时频峰值滤波降噪方法(S-TFPF).在S-TFPF方案中,对于有效信号集中区,S-TFPF方法根据信号频率特征,选择有利于信号保持的较短时窗长度;对于噪声集中区,按噪声方差自适应增加时窗长度,增强随机噪声压制能力.S-TFPF应用于合成记录和共炮点记录的滤波结果表明,与传统时频峰值滤波方法相比,S-TFPF方法可以有效抑制低信噪比地震勘探记录中的随机噪声,更好地恢复出同相轴.     

基于结构自适应中值滤波器的随机噪声衰减方法

本文提出一种保护断层、裂缝等地层边缘特征的结构自适应中值滤波器,用于衰减地震资料中的随机噪声.基于地震反射同相轴局部呈线型结构的假设,采用梯度结构张量估计地层倾向,分析地层结构的规则程度,在此基础上引入地震剖面中线型和横向不连续性两种结构特征的置信度量.结构自适应中值滤波器根据这两种置信度量调整滤波器窗函数的尺度和形状,根据地层倾角调整滤波器窗函数的方向,从而使得滤波操作窗能够最佳匹配信号的局部结构特征.将本文方法用于合成和实际数据的处理,并与两种常用中值滤波方法进行对比,结果表明,该方法能够更好地解决地震剖面的随机噪声衰减和有效信号保真的问题,在增强反射同相轴的横向一致性的同时有效保持了剖面内的地层边缘和细节特征,显著改善了地震资料的品质.     

基于复数域混合SVD滤波方法及在随机噪声压制中的应用

在低信噪比地震资料处理中,压制随机噪声是其中的关键处理环节.传统的频率空间域预测滤波方法,容易使得去噪后高频段的有效信号严重畸变,不利于进一步提高分辨率.为此,提出了一种基于复数域混合SVD滤波法压制三维地震数据中随机干扰的方法,该方法首先在时空域对地震数据作傅里叶变换,再依据Eigenimage滤波法与Cadzow滤波法建立混合Hankel矩阵,然后对其运行奇异值分解,最后通过秩约化的方法来压制随机干扰.理论模型和实际地震数据的应用表明:该方法可有效地去除随机噪声,保护有效波,明显地改善了叠后三维地震资料的信噪比.     

一种优化预测滤波压制随机噪声的方法

针对低信噪比地震资料中的随机噪声干扰,提出了一种优化的预测滤波方法.通过对AR模型数学表达式进行改写,然后引入反褶积代价函数,最后推导得到压制随机噪声的递归算法.该算法实现了单道压制随机噪声,保护了有效信号的细节部分,比基于常规最小二乘代价函数的滤波算法收敛更快.通过模拟和实际地震资料验证表明,该方法能较好地压制低信噪比资料中的随机噪声干扰,同时能较好地保护有效地震信号.     

微地震信号的变速FK滤波-自适应极化滤波方法

针对微地震信号的特点,研究了变速FK滤波-自适应极化滤波方法.在设计FK滤波因子时,通过计算相邻道滑动时窗内信号的互相关系数,确定时窗内信号的视速度范围,实现变速FK滤波;对时窗内的信号进行二维FFT变换时,分别在时窗的纵横向上补零充位,消除信号在二维变换时造成的时空域的混跌和泄漏;在设计自适应极化滤波因子时,针对复杂波场中波的偏振方向的不确定性问题,本文通过计算相邻道信号偏振投影的最大互相关系数,确定波的跟踪分量;把波的跟踪分量作为极化滤波因子里的期望方向,改进常规滤波因子,实现自适应极化滤波.对理论模型和实际资料的处理结果表明,该方法理论正确,实际资料的处理效果也得到了明显的改善.     

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

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

局部相关加权中值滤波技术及其在叠后随机噪声衰减中的应用

随机噪声的衰减和同相轴连续性的提高可以极大地改善地震资料解释的精度.本文提出一种新的滤波技术,既能够有效地衰减随机噪声又可以很好地保护地震资料中的断层等信息不被破坏,增强地震剖面中弯曲、倾斜同相轴的连续性.该方法结合新的加权中值滤波技术和两种构造信息保护滤波策略,实现基于预测数据体和基于倾角走向的加权中值滤波.通过设计...     

Suppressing random noise to broaden the useful bandwidth of seismic reflection data

The signal-to-noise ratio （SNR） of seismic reflection data in many areas is rather poor and conventional two-dimensional filters designed to suppress noise with different moveout from the signal tend to generate artifacts. We have extended a method of multichannel filtering, based on the hypothesis that signals on adjacent channels are similar, for enhancing the SNR on stacked sections. Using only the mid-range frequencies where the SNR is highest, the event trend is found for overlapping windows on the section and the average signal vector is calculated. Then the data from the full bandwidth section are projected onto the spatially varying unit similarity vectors and the results are merged for the overlapping windows. Application of the method to synthetic data containing steeply dipping events and to a stacked section for a marine 2D line has produced good results. The modifications we have introduced carry a small overhead in computing time but they should enable the method to be used effectively even on sections containing steep dips.     

Further improvement of temporal resolution of seismic data by autoregressive (AR) spectral extrapolation

Seismic data have still no enough temporal resolution because of band-limited nature of available data even if it is deconvolved. However, lower and higher frequency information belonging to seismic data is missing and it is not directly recovered from seismic data. In this paper, a method originally applied by Honarvar et al. [Honarvar, F., Sheikhzadeh, H., Moles, M., Sinclair, A.N., 2004. Improving the time-resolution and signal–noise ratio of ultrasonic NDE signals. Ultrasonics 41, 755–763.] which is the combination of the most widely used Wiener deconvolution and AR spectral extrapolation in frequency domain is briefly reviewed and is applied to seismic data to improve temporal resolution further. The missing frequency information is optimally recovered by forward and backward extrapolation based on the selection of a high signal–noise ratio (SNR) of signal spectrum deconvolved in signal processing technique. The combination of the two methods is firstly tested on a variety of synthetic examples and then applied to a stacked real trace. The selection of necessary parameters in Wiener filtering and in extrapolation are discussed in detail. It is used an optimum frequency windows between 3 and 10 dB drops by comparing results from these drops, while frequency windows are used as standard between 2.8 and 3.2 dB drops in study of Honarvar et al. [Honarvar, F., Sheikhzadeh, H., Moles, M., Sinclair, A.N., 2004. Improving the time-resolution and signal–noise ratio of ultrasonic NDE signals. Ultrasonics 41, 755–763.]. The results obtained show that the application of the purposed signal processing technique considerably improves temporal resolution of seismic data when compared with the original seismic data. Furthermore, AR based spectral extrapolated data can be almost considered as reflectivity sequence of layered medium. Consequently, the combination of Wiener deconvolution and AR spectral extrapolation can reveal some details of seismic data that cannot be observed in raw signal or which lost during the previous processing.     

Multi-refractor imaging with stacked refraction convolution section

Multi‐refractor imaging is a technique for constructing a single two‐dimensional image of a number of refractors by stacking multiple convolved and cross‐correlated reversed shot records. The method is most effective with high‐fold data that have been obtained with roll‐along acquisition programs because the stacking process significantly improves the signal‐to‐noise ratios. The major advantage of the multi‐refractor imaging method is that all the data can be stacked to maximize the signal‐to‐noise ratios before the measurement of any traveltimes. However, the signal‐to‐noise ratios can be further increased if only those traces that have arrivals from the same refractor are used, and if the correct reciprocal times or traces are employed. A field case study shows that multi‐refractor imaging can produce a cross‐section similar to the familiar reflection cross‐section with substantially higher signal‐to‐noise ratios for the equivalent interfaces.     

DECOMPOSED VIBRATOR PATTERNS TO IMPROVE SEISMIC SURVEY RESULTS1

Effective noise reduction of single-sweep recorded data is achieved by application of a velocity filter process on a decomposed vibrator pattern. This technique promises high resolution results with a minimum effect on signal characteristic. A comparison of the stacked section of records vertically stacked in the field with the stacked section of velocity-filtered receiver gathers shows a significant increase in resolution and signal-to-noise ratio.     

利用数字图像处理技术提高地震剖面图像信噪比

提出了利用数字图像处理技术提高地震剖面信噪比的新方法，首先根据数字图像处理要求的格式，对地震剖面数据进行转换，得到地震剖面图像，分析了地震数据特点和初步地震图像的实验结果后，设计了新的预处理方法——“二维沿层滤波”，在此基础上，利用可以计算帧间运动速度及其变化都较大的改进的光流分析技术，计算出多幅地震剖面对应点的偏移量，然后应用图像积累技术对这多幅地震剖面进行积累，实现对三维地震数据体提高信噪比的处理，该方法充分利用了三维地震信息，不但可以提高整个数据体的信噪比，而且可以减少信号能量的损失，并保持原来的信号能量关系，使地震剖面的质量得到明显提高，为地震解释奠定良好的基础。     

S变换谱分解技术在深反射地震弱信号提取中的应用

在深反射地震资料处理中,当来自深部的有效弱信号和噪声干扰频带差异较小且难以区分时,传统滤波方法的应用会受到限制.谱分解方法是一种使用离散傅里叶变换,基于信号的频率-振幅谱等信息生成高分辨率地震图像的方法,通常用来识别介质物性横向分布特征,处理复杂介质内频谱变化和局部相位的不稳定性等问题,包括定位复杂断层和小尺度断裂等.S变换作为一种新的时频分析方法,具有自动调节分辨率的能力,近些年来被广泛应用到勘探地震、大地电磁等数据处理中,逐渐成为地球物理方法中噪声压制的有效方法之一.与常规石油反射地震资料相比,深反射主动源地震为了探测深部结构信息,常采用大药量激发方式、长排列观测系统等,导致深部有效信号基本湮灭在噪声干扰之中.针对深反射数据特点,本文结合谱分解和S变换技术,首先设计了简单的脉冲函数实验数据,证实S变换方法的有效性,同时说明谱分解方法的效果受所用时频分析方法影响较大,而其中决定分辨能力的变换窗函数的选取尤为重要.在此基础上,分别应用到深反射地震资料的单道和叠加剖面实际数据上,对比分析了传统变换谱分解和S变换谱分解的应用效果,单道资料对比结果表明:相比传统谱分解,S变换谱分解方法具有自动调节分辨率的能力,能够精确的标定深反射地震资料中弱信号不同时刻的频率分量;叠加剖面资料应用结果表明:由S变换谱分解得到的剖面结果与其他谱分解方法结果整体上具有较高的一致性,同时清晰地刻画出原叠加剖面上被噪声湮灭的低频细节特征,提高了剖面的分辨率及同相轴连续性;对比结果明显看出,Gabor变换谱分解方法得到的结果同相轴较为破碎,分析原因认为这是由Gabor变换的时频分解方法的定长窗函数所致,窗口大小不会随着信号频率的变化来调节长度,只能在处理的过程中根据一定的记录长度范围选取窗函数参数,而S变换谱分解方法在窗函数的选取时,通过时变信号的局部频率特征自动调节窗口长度,能够更好的刻画各个频段的细节特征,在深反射剖面成像应用中效果尤为明显.本文结果表明S变换谱分解技术在深地震叠加剖面上的应用有效地提高了来自深部弱反射信号的信噪比和分辨率,并刻画出了叠加剖面上所不具有的低频细节特征,在实际深反射地震资料处理中能有效保护低频弱信号获得更好的成像效果.本文为深地震反射资料中弱信号的保护处理找到一种有效的方法.     

小波域叠前地震资料衰减参数定性估计方法

叠前地震资料中,高频分量和低频分量随传播距离的衰减特性不同.本文给出了一种在小波域定性估计叠前地震资料衰减参数的方法.该方法利用连续小波变换提取共反射点道集的高、低频分量,以低频分量和高频分量之差定性反映地震波的衰减.通过累加不同偏移距的衰减,提高了估计的稳定性;采用幅度归一化方法,降低了信号幅值对衰减参数估计的影响.将本文提出方法与常用的基于叠后地震资料衰减估计方法用于某油田的地震资料处理,结果表明,本文方法得到的衰减估计结果能够更好地反映油气的空间展布.     

可控震源匹配扫描方法研究

常规的可控震源Chirp扫描信号的地震响应存在着相关信号旁瓣大、分辨率低的缺点;用二元伪随机编码信号作为可控震源的驱动信号,地震响应剖面中存在严重的相关噪声,降低了地震记录的信噪比.本文提出的可控震源匹配扫描方法,利用匹配伪随机序列偶调制正弦载波信号作为可控震源驱动信号进行扫描激发,得到两个匹配的原始地震记录,对两个匹配的原始地震记录分别解码后再进行叠加形成最终的综合解码地震剖面.地震数值模拟的结果表明,该方法基本不存在常规Chirp扫描情形的旁瓣效应,可有效避免m-序列编码扫描情形的相关噪声,其最后形成的综合解码地震剖面可以与脉冲震源的炮集记录相媲美.     

CGP叠加柱面波偏移成像方法

针对松辽盆地薄砂层油藏的地震勘探问题,提出了一种共检波器接收点(CGP)叠加柱面波偏移成像方法。该方法对小炮点距、小检波器点距、中间放炮观测系统采集地震资料,经CGP道集叠加组合成柱面波剖面;采用下行波射线法向下延拓和上行波波动方程向下延拓的方法,使柱面波剖面偏移成像。通过模型分析和对松辽盆地TK8157测线资料进行处理证明该方法的地震分辨率和保真度较高,在发现小砂体、小断层、地层尖灭等方面有较好的效果。     

Common reflection surface stack using dip decomposition for rugged surface topography

We present an extension of the Common Reflection Surface （CRS） stack that provides support for an arbitrary top surface topography. CRS stacking can be applied to the original prestack data without the need for any elevation statics. The CRS-stacked zero- offset section can be corrected （redatumed） to a given planar level by kinematic wave field attributes. The seismic processing results indicate that the CRS stacked section for rugged surface topography is better than the conventional stacked section for S/N ratio and better continuity of reflection events. Considering the multiple paths of zero-offset rays, the method deals with reflection information coming from different dips and performs the stack using the method of dip decomposition, which improves the kinematic and dynamic character of CRS stacked sections.