Morlet小波匹配追踪方法的时频原子参数分析

匹配追踪方法是目前时频分辨率最高的时频分析方法.但传统的匹配追踪方法在寻找与信号匹配最佳的时频原子时,需在较大的参数搜索范围内对时频原子参数进行多重循环寻优迭代,计算效率较低.针对匹配追踪算法实现方式,本文主要工作及结论如下:(1)以Morlet小波为时频原子,以频率、相位、时延和尺度为时频原子参数,分析了信号在时频原子上的投影振幅对这四项时频原子参数的灵敏度,结果表明灵敏度顺序为:时延、频率、尺度和相位,据此结论可在算法实现时根据灵敏度顺序调整不同参数寻优迭代时的参数搜索范围,改善算法的计算效率;(2)针对频率、相位和时延参数已知尺度参数变化的时频原子,分析了利用Hilbert变换得到的瞬时参数表征时频原子参数初值的准确度,结果表明准确度顺序为:时延、频率和相位.据此结论可在算法实现时确定参数寻优的优先级顺序,改善算法的计算效率.与连续小波变换、广义S变换等时频分析方法的对比及在地震数据去噪中的应用效果表明匹配追踪方法效果较好.     

谱图重排的谱分解理论及其在储层探测中的应用

谱分解理论是把单道地震记录分解为连续的时频谱平面,是地震资料处理和解释的重要技术之一.由于谱分解方法的多解性,所以同一道地震记录因为分解方法不同,得到的时频谱是不一样的.短时傅里叶变换,小波变换,S变换和匹配追踪算法都是对信号开窗分析,这些方法都受到不确定性原理的限制.Wigner-Ville变换避开了不确定性原理的限制,但是交叉项的存在限制了本方法的使用.本文利用谱图重排的时频分析方法(RSPWV)对合成的单道地震记录和实际的地震资料进行了分析.与短时傅里叶变换和匹配追踪算法的比较得出:此方法具有较高的时频分辨率,能够很好地识别气层.     

基于结构化匹配追踪的地震资料去噪技术研究

地震数据往往受噪声干扰,影响有效反射波的识别,因此提高地震资料信噪比在数据处理中尤为重要,我国在此方面研究已经取得进展.结构化匹配追踪算法是一种采用原子字典分组式及树状结构搜索的全局结构化算法,具有运算效率高、速度快的特点.本文在分析结构化匹配追踪算法的基础上提出基于该算法的地震资料去噪技术,根据反射波信号与噪声信号的频段差异,对分解得到的时频原子设置阈值频段,去除含噪成分,从而达到去噪目的.与常规去噪技术相比,该技术精准、快速地分解地震反射波信号与噪声信号,最大限度地减小对地震反射波信号的损害.     

匹配追踪算法及其在煤层气富集区预测中的应用

为了更准确地刻画地震信号的时频信息,实现高效率的时频谱分解,本文引入一种基于复数道地震记录的匹配追踪算法,阐述了该算法的基本原理和具体算法流程,基于合成记录进行了匹配追踪算法的抗噪能力和端点效应分析,并对实际煤田地震资料进行了匹配追踪的计算和分析.结果表明:1)匹配追踪算法更适用于信噪比为10dB以上的信号;2)端点平滑方法的引入可有效提高算法的准确性,解决了因端点振幅值过大导致的算法陷入局部匹配的问题;3)实际煤田地震资料中,煤层气富集区下部具有低频能量强、高频能量弱的特点,可将其作为煤层气富集区预测的一个标志.     

基于正交时频原子的地震信号快速匹配追踪

匹配追踪能够根据地震信号自身的特点进行自适应分解,但因计算量巨大使其不能被广泛应用.本文提出了一种正交时频原子快速匹配追踪方法,该方法在迭代分解过程中,充分利用地震信号的局部性特征作为先验信息点,在先验信息点附近采用动态搜索策略寻找最佳匹配时频原子,同时对时频原子进行正交变换,消除了时频原子库中的冗余分量,最终将信号分解为一系列正交时频原子的线性叠加.测试结果表明,该方法不仅保持了匹配追踪的分解精度,而且使计算效率有了质的提高.     

地震反演驱动的改进匹配追踪煤层识别方法

多类型油气藏储层受到煤层反射干扰, 有效识别与去除煤层地震响应是提高煤系储层地震描述精度的关键技术之一.常规煤层匹配识别欠缺考虑煤层弹性参数的影响, 引起煤层定位不准确、识别可靠性低、运算效率低等问题, 本文充分考虑了弹性参数对煤层具有较好的指示作用, 首次提出了地震反演驱动的改进匹配追踪煤层识别方法.首先根据岩石物理分析, 构建煤层指示因子, 提出匹配追踪煤层指示参数的地震反演方法, 该参数能够较好地刻画煤层空间展布; 基于此, 在地震反演驱动下获得了煤层的先验位置, 将其引入匹配追踪算法中约束匹配追踪搜索邻域, 为缓解现有匹配追踪算法运行效率缓慢, 将经验模态分解(EMD)引入匹配追踪领域, 同时将连续相位替换瞬时相位加快匹配追踪运算效率, 研发了煤层指示因子约束下基于EMD字典的改进匹配追踪煤层识别算法.模型测试验证了本文方法识别煤层的可行性, 实际资料的计算结果表明, 该方法在保持精度的情况下大大提高了计算效率, 此外, 预测的煤层与实际测井解释的煤层基本一致, 有助于油气储层预测和识别, 具有重要的实际应用价值.

     

S变换的时频分析特性及其改进

地震信号具有非线性、非平稳的特性,用时频分析的方法处理地震信号,可以最大限度的保留原始信息,也可以较为精确地分析数据的时间、频率特征.现阶段应用于地震数据处理的时频分析方法主要有:短时傅里叶变换,小波变换和匹配追踪算法等等.S变换是近几年发展起来的一种新的时频分析方法.S变换结合了短时傅里叶变换和小波变换的优点,具有相...     

基于匹配追踪稀疏分解的高斯束成像方法

高斯束偏移是一种高效、稳健的深度域成像方法,其不仅保持了射线类偏移的高效、灵活性,而且具有接近波动方程偏移的成像精度.匹配追踪是一种基于匹配寻优算法的信号稀疏分解技术,常用于地震信号的时频分析和去噪处理中.本文将建立在Ricker子波原子库的匹配追踪算法应用于常规的高斯束偏移中,并结合Kirchhoff偏移的单输入道成像方式,发展了一种成像精度更高且适用于低信噪比资料的叠前高斯束成像方法.该方法通过合理地控制匹配追踪稀疏分解的迭代次数,可以有效地去除地震信号中的随机干扰,提高成像结果的信噪比;此外,在偏移过程中,本文方法采用了Kirchhoff偏移的单输入道成像方式,解决了常规高斯束方法对浅层小尺度地质体成像不准的问题,提高浅部反射层的成像精度.两个典型的数值算例验证了本文方法的有效性和适应性.     

复数子波匹配追踪算法识别薄层砂体

针对东营凹陷滩坝砂岩薄砂体的预测识别问题,研究了复小波匹配追踪方法.从地震信号morlet小波分解及重构分析入手,讨论了重构地震信号时关于频率、振幅、相位及时移参数的选取方法及其具体的选取变量.并就关键的振幅参数,研究了改善的计算方法.常规的匹配追踪算法扫描地震信号的所有时间点和频率.在本文中,为了改善时频分辨率,采用局部扫描主频和子波的时间延迟.在具体实现过程中,采用计算地震信号的平均瞬时频率和地震信号的瞬时包络处的时间,做为采样频率和时间采样点,得到与地震信号最佳相关的子波.对研究的方法从理论模型和实际资料进行了论证分析,证明理论上正确合理,实际应用效果显著.     

频谱成像技术研究进展

频谱成像技术是近年来发展起来的一项基于频率谱分解的储层特色解释技术,是地震属性分析中重要组成部分.频谱成像技术具有在空间横向上分辨率高的特点,是一种利用三维地震资料的多尺度信息对储层进行高分辨率成像、检测储层时间厚度变化的工具.地震数据振幅谱可以识别地层的时间厚度变化,相位谱可以检测地质体横向不连续性,在确定油藏边界、计算地层厚度方面比传统地震属性研究方法具有更大的优势.频谱成像技术的核心模块是信号的时频分析,本文回顾了频谱成像技术的基本概念和目前在频谱成像中广泛应用的时频分析方法,简单介绍了几种新的时频分析方法:广义S变换和改进的匹配追踪算法.     

自适应多道匹配追踪去强屏蔽方法研究与应用

阻抗差异较大的地层在地震剖面上呈现为强振幅同相轴,会掩盖附近储层有效信息,需要做去除强屏蔽的目标处理。针对常规匹配追踪方法在构造复杂地区匹配精度和空间连续性不佳的问题,本文提出一种基于自适应权值的多道匹配追踪去强屏蔽方法。首先利用层位构造信息对强反射层进行局部拉平处理,减弱地层空间构造对提取强反射层的影响,然后引入相邻地震道与中心道之间的相关系数作为多道平均的权值,以提高匹配结果的稳定性和横向连续性。同时使用解释得到的层位时间作为子波初始时间,有效提高运算效率。模型试算和实际地震资料应用分析表明,改进方法能有效剥离强反射并突出有效储层信息,剥离后的井震吻合度得到明显提升,相较于常规匹配追踪算法匹配精度更高,同时具有更好的空间连续性与强反射去除效果。      

The attenuated Ricker wavelet basis for seismic trace decomposition and attenuation analysis

The widely used wavelets in the context of the matching pursuit are mostly focused on the time–frequency attributes of seismic traces. We propose a new type of wavelet basis based on the classic Ricker wavelet, where the quality factor Q is introduced. We develop a new scheme for seismic trace decomposition by applying the multi-channel orthogonal matching pursuit based on the proposed wavelet basis. Compared with the decomposition by the Ricker wavelets, the proposed method could use fewer wavelets to represent the seismic signal with fewer iterations. Besides, the quality factor of the subsurface media could be extracted from the decomposition results, and the seismic attenuation could be compensated expediently. We test the availability of the proposed methods on both synthetic seismic record and field post-stack data.     

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.     

AVO信息约束匹配追踪技术在富煤储层刻画中的应用

针对西湖凹陷富煤环境下储层刻画精度低问题,本文结合煤层AVO截距、梯度特征,提出一种基于AVO信息约束的匹配追踪技术,压制煤层强反射引起的岩性假象,凸显储层真实、有效信号。该方法首先利用煤层4类AVO负强截距P、正强梯度G特点,构建煤层地震敏感因子P - G,放大煤层地震响应,并压制非煤层强振幅影响,实现煤层位置精细定位;在此基础上,将该煤层地震信息作为匹配追踪需要分解、重构的原始信号,利用复地震道分析技术提高信号快速匹配分解的效率,完成煤层强反射解耦。模型试算及实际资料应用表明:匹配追踪技术在精细定位煤层地震响应基础上,提高了匹配追踪算法去煤层强振幅效率;煤层解耦后地震数据较好地凸显储层横向展布变化,提高主力气层的纵向刻画精度。      

多尺度快速匹配追踪多域联合地震反演方法

多尺度快速匹配追踪多域联合地震反演是一种通过地震数据多尺度分解的迭代反演方法.与此同时,在快速匹配追踪算法中引入低频模型约束,有效提高了收敛精度,使反演结果具有丰富的高低频信息.首先通过对大尺度地震资料进行反演得到低频背景.在此基础上,采用中尺度与小尺度地震数据进行逐级迭代用以获得高频数据,因而有效缓解了常规反演方法对于初始模型精度的依赖.最后利用理论模型与实际地震数据进行测试,通过与常规时间域反演方法的反演结果进行对比可以看出,本文方法在地层连续变化处依然可以对变化地层进行精确刻画,且在纵向分辨率提升的同时保持了较好的横向连续性.

     

多尺度快速匹配追踪多域联合地震反演方法

多尺度快速匹配追踪多域联合地震反演是一种通过地震数据多尺度分解的迭代反演方法.与此同时,在快速匹配追踪算法中引入低频模型约束,有效提高了收敛精度,使反演结果具有丰富的高低频信息.首先通过对大尺度地震资料进行反演得到低频背景.在此基础上,采用中尺度与小尺度地震数据进行逐级迭代用以获得高频数据,因而有效缓解了常规反演方法对于初始模型精度的依赖.最后利用理论模型与实际地震数据进行测试,通过与常规时间域反演方法的反演结果进行对比可以看出,本文方法在地层连续变化处依然可以对变化地层进行精确刻画,且在纵向分辨率提升的同时保持了较好的横向连续性.     

High‐definition frequency decomposition

Spectral decomposition is a widely used technique in analysis and interpretation of seismic data. According to the uncertainty principle, there exists a lower bound for the joint time–frequency resolution of seismic signals. The highest temporal resolution is achieved by a matching pursuit approach which uses waveforms from a dictionary of functions (atoms). This method, in its pure mathematical form can result in atoms whose shape and phase have no relation to the seismic trace. The high‐definition frequency decomposition algorithm presented in this paper interleaves iterations of atom matching and optimization. It divides the seismic trace into independent sections delineated by envelope troughs, and simultaneously matches atoms to all peaks. Co‐optimization of overlapping atoms ensures that the effects of interference between them are minimized. Finally, a second atom matching and optimization phase is performed in order to minimize the difference between the original and the reconstructed trace. The fully reconstructed traces can be used as inputs for a frequency‐based reconstruction and red–green–blue colour blending. Comparison with the results of the original matching pursuit frequency decomposition illustrates that high‐definition frequency decomposition based colour blends provide a very high temporal resolution, even in the low‐energy parts of the seismic data, enabling a precise analysis of geometrical variations of geological features.     

基于约束插值算法的井资料外推

Based on a detailed analysis of differences between seismic data and well logs, we discuss the problem of matching seismic traces and well logs and present a new matching method based on event search in instantaneous phase which greatly improves seismic resolution. The method is based on flattening events in instantaneous phase to compare the seismic traces to the well log traces with the same phase. We calculate the coefficients using the singular value decomposition method to extrapolate the well logs. As a result, the events in the seismic profile are continuous and match well with well logs. We apply this method to the Mao-2 well in Daqing Oilfield with good results.