基于动态褶积模型的动态子波估计

目前常用的地震子波估计方法基于传统的褶积模型，没有考虑子波在传播过程中的动态衰减特性。针对这一问题，对传统褶积模型进行了改进，建立了能够描述子波衰减特性的动态褶积模型，并以动态褶积模型为基础提出了动态子波估计方法。该方法首先对衰减地震道做频谱分析；然后假设反射系数的频谱为白谱，通过对衰减地震道频谱值做平滑，去除反射系数谱值的影响，获取动态子波的振幅谱值；最后对子波相位做一定假设，从而求得动态子波。理论模型分析与实例计算结果表明，该方法是可行的。     

动态地震子波估计

目前常用的地震子波估计方法有确定性子波估计和统计性子波估计,各有其优缺点.2类子波估计方法都是基于传统的褶积理论模型,即子波在传播过程中保持波形不变,没有考虑子波在传播过程中的动态性.实际上由于地层对子波具有吸收作用,子波在传播过程中逐步衰减,子波传播是一个动态衰减过程.改进静态褶积模型提出动态褶积模型,从而估计动态子波.假设地层反射系数为白噪序列,震源子波在地层中传播形成衰减地震道,对此地震道作小波变换时频谱分析,采用箱状平滑对衰减地震道频谱作平滑,去除反射系数的影响,从而获得子波振幅谱,再假设子波为最小相位,通过希尔伯特变换求得子波相位,从而实现动态子波估计.理论模型分析与实例计算验证了方法的可行性.     

动态反褶积研究及效果分析

目前常用的反褶积方法有多种,这些反褶积方法都是基于传统的褶积模型,在做反褶积前,先要对衰减地震道做振幅补偿,因而在提高分辨率时具有局限性。在传统的褶积模型中,假设子波在传播过程中振幅不发生变化,实际上子波在传播过程中振幅是逐渐衰减的,因此需要建立能够描述子波动态传播的褶积模型。以动态褶积模型为理论基础,不需要对衰减地震道做振幅补偿,可直接实现分辨率处理。将对脉冲反褶积与动态反褶积结果进行比较,效果表明,动态反褶积方法具有双重功能:具有提高地震分辨率的功能;具有由于吸收引起振幅衰减的补偿功能。理论模型分析与实例计算验证了动态反褶积方法的正确性。     

动态反褶积方法研究

目前常用的反褶积方法有许多,但大都基于传统的褶积模型。传统褶积模型假设子波波形保持不变,这与子波在传播中振幅逐渐衰减的情况不符,因此需要建立能够描述实际子波传播的动态褶积模型。以动态褶积模型为理论基础的反褶积方法不需要对衰减地震道做吸收衰减补偿,可直接实现提高分辨率处理。文中对脉冲反褶积和动态反褶积结果进行了比较,理论模型分析与实例验证了动态反褶积方法的正确性。     

动态反褶积研究及效果分析

目前常用的反褶积方法有多种，这些反褶积方法都是基于传统的褶积模型，在做反褶积前，先要对衰减地震道做振幅补偿，因而在提高分辨率时具有局限性。在传统的褶积模型中，假设子波在传播过程中振幅不发生变化，实际上子波在传播过程中振幅是逐渐衰减的，因此需要建立能够描述子波动态传播的褶积模型。以动态褶积模型为理论基础，不需要对衰减地震道做振幅补偿，可直接实现分辩率处理。将对脉冲反褶积与动态反褶积结果进行比较，效果表明，动态反褶积方法具有双重功能：具有提高地震分辨率的功能；具有由于吸收引起振幅衰减的补偿功能。理论模型分析与实例计算验证了动态反褶积方法的正确性。     

动态反褶积中的反射系数序列时频特征研究

基于动态褶积模型的反褶积方法假设反射系数序列是白噪的,而实际地震资料的反射系数序列通常是非白噪的,因而会影响反褶积效果。基于动态褶积模型研究了反射系数白噪和非白噪情况下地震道信息、反射系数和子波三者之间的关系,指出地震信号的对数时频谱是子波的对数时频谱和反射系数的对数时频谱之和。通过白噪和非白噪反射系数序列模型试验,研究了利用合成地震道的时频谱和对数时频谱求取子波时频谱的效果。试验结果表明:用对数时频谱求取的子波谱在频率方向精度更高;在反射系数序列非白噪的情况下,通过平滑合成地震道的对数时频谱获取子波时频谱,动态反褶积的效果更好。某油田实际地震资料处理结果验证了上述结论。     

Cauchy稀疏约束Bayesian估计地震盲反褶积框架与算法研究

以二阶统计学方法为基础,从Canadas等提出的非最小相位子波和非白噪反射系数地震盲反褶积框架出发,给出了Cauchy稀疏约束Bayesian估计地震盲反褶积框架。基于反射系数与子波相互独立（或弱相关）的假设,分别构建了反射系数和子波最优估计方程,并采用预条件共轭梯度法迭代反演实现反射系数和子波的同时估计。在方法具体实现时,以传统脉冲反褶积结果作为迭代初值,通过迭代得到反射系数和任意相位子波;然后再对子波进行最小相位化,通过反演得到反子波;最后将反子波与地震道褶积,得到反褶积结果。利用理论模型和实际数据对算法进行了试算,结果表明,给出的地震盲反褶积理论框架是正确的;与直接（共轭梯度求解正则方程）稀疏同时迭代反演法的对比显示,预条件共轭梯度算法稳定,精度高,收敛快。     

用统计法提取地震传播子波

一个地震道模型常常给定为地震子波与地下反射系数序列相褶积的结果,再在该道上加上白噪声。因而,了解传播子波是我们从地震道中估算反射系数序列的立足点。假设地震子波是唯一的,本文讨论了从几个地震道提取子波的统计方法。该方法使我们在不考虑相位谱的传统限定假设的前提下获得地震子波。此外,在这篇文章中,我们提出了相位展开方法,并考察了一些记录道中反射系数相位谱的统计特性。     

一种新的混合相位反褶积方法

在地震资料处理中,常用预测反褶积方法提高地震资料的分辨率,它具有两个基本的假设条件:①地震子波是最小相位;②反射系数序列是随机过程。上述假设条件与实际地震子波的混合相位存在较大差别。为此,笔者在借鉴前人研究的基础上,提出了一种新的混合相位反褶积方法。其实现过程为:①根据地震记录的自相关和地震子波自相关之间的关系,由最小二乘法拟合出地震子波的振幅谱,进而计算出地震子波的自相关;②利用计算出来的地震子波自相关设计一个最小相位滤波器,然后加一个延迟量于自相关函数的主对角线上,将最小相位反滤波器分解为一个长度很短的滤波器,再与一个长度较长的滤波器进行褶积,可得到一系列的混合相位子波和混合相位滤波器;③将得到的混合相位滤波器和原始地震道进行褶积,选取一定的判别准则,以选出最佳的混合相位滤波器以及相应的混合相位子波。通过模型试算和实际资料处理可知,分时窗反褶积的效果要好于单时窗反褶积,混合相位反褶积效果好于最小相位反褶积。在分时窗求取地震子波时,为了得到更好的反褶积效果,要注意时窗宽度与算子长度之间的关系。     

基于谱模拟技术的混合相位地震子波估计方法

针对现有的地震子波振幅谱模型的不足提出了新的振幅谱模型及其估计方法。介绍了子波的Z变换在单位园上无零点和子波自相关函数已知的条件下混合相位地震子波估计的方法原理和判别准则,并在此基础上提出了基于地震子波振幅谱模拟技术的混合相位地震子波估计方法。合成和实际数据处理结果均表明,基于地震子波振幅谱模拟技术的混合相位地震子波估计方法具有较好的性能。     

应用时变子波的盲反射系数反演

由于地震数据的非稳态和全盲特性,实际地震子波通常是未知的且其波形会随传播时间发生变化,故常规基于稳态子波假设的反射系数反演精度将难以保证。为此,利用广义S变换将非稳态地震数据变换到时频域,基于自相关理论逐点提取子波并重新构建时变子波矩阵,而不再是从一道地震记录中仅提取一个时不变子波,将其应用于反演模型以实现非稳态地震资料反射系数的盲反演。模型测试和实际资料处理结果表明,相较于传统时不变子波反演方法,所提方法考虑了地震信号的时变特征,且是数据驱动的,更适用于实际地震资料,能得到更精确的时变子波结果和较高分辨率的反射系数剖面。     

基于分形方法的地震子波提取及应用

常规自相关子波求取方法假设反射系数功率谱为白噪谱，因而在实际应用中受到限制。以分形理论为基础的地震子波功率谱估计方法，可突破反射系数功率谱为白噪的假设约束，适应反射系数序列具有一定相关性的地震资料处理。理论模型和实际地震资料处理表明，分形技术应用于复赛谱求解地震子波的方法中，可适应非白噪随机反射系数序列地震记录的子波求取。用基于分形技术的复赛谱法求取的子波对地震记录进行反褶积处理，分辨率明显高于自相关子波求取方法。     

混合相位地震子波提取及应用

地震子波在地震正演和反演中有着重要的作用,目前的子波提取主要基于最小相位假设,而实际的地震子波往往表现为混合相位特征。为了 提取具有真实相位的地震子波,从地震信号的统计性出发,利用高阶统计量中的双谱,从反射地震记录中提取出地震子波的真实相位,进一步重构 出地震子波。经过模型仿真和实际资料处理,证明了该地震子波提取方法的有效性和实用性。     

Research on the methods for wavelet stability of homomorphic deconvolution.

同态反褶积受噪声及反射系数非随机性的影响,计算出的子波相位稳定性差,难以用于实际生产。使用小时窗、多道数据计算平均复赛谱和平均自相关函数复赛谱,利用二者计算混合相位子波,从而增加子波相位的稳定性,使同态反褶积的实用性大大增强。在实际资料处理过程中,使用了多种方法增加子波稳定性：①对地震数据进行指数衰减,使地震子波变为最小相位;②对求出的子波进行能量定位,选择能量集中的部分,对计算出的子波进行窗口整形,从而减少子波旁瓣的影响;③调整计算子波与期望子波峰值时差;④在三维叠后地震数据中使用Inline和Crossline双方向道计算中央道算子。理论模型试验结果表明,该方法适用于混合相位子波地震数据;实际数据应用结果表明,采用该方法计算出的子波稳定性强,反褶积后地震数据的分辨率得到了提高。     

基于S变换和变分法的品质因子Q估计方法

品质因子Q是定量描述黏弹性衰减的重要参数,准确估计Q值有利于储层识别、烃类检测,还可用于反Q滤波来提高地震资料的分辨率。传统的Q值估计方法包括对数谱比法、中心频率偏移法和峰值频率偏移法等,其中对数谱比法的抗噪性较差,而中心频率偏移法和峰值频率偏移法则存在依赖地震子波类型的问题。针对这些问题,提出一种基于S变换和变分法提取品质因子Q的方法。首先,通过研究非平稳褶积模型,推导出非平稳地震数据在S变换域的近似表示形式;其次,在近似表示的基础上,建立了关于品质因子Q和地震子波的目标函数,并基于变分法最小化目标函数,得到Q值估计表达式;最后,为了提高该方法的准确性与抗噪性,设计了一种自适应积分区间选择的方案,该方案可根据地震数据的时频谱自动计算积分区域的频率参数。模型算例测试结果表明,所提出的Q值估计方法不依赖于地震子波的类型和窗函数的长度,同时对噪声具有较好的鲁棒性。实际数据计算结果进一步验证了该方法在Q值估计中的有效性。     

基于时频谱模拟的时变混合相位子波提取

为解决地层衰减作用造成传统子波提取方法提取结果不准确的问题,提出一种基于时频谱模拟的时变混合相位子波提取方法.该方法在对非平稳地震记录进行改进的广义S变换的基础上,首次引入时频滤波处理,然后采用谱模拟方法估计每一时刻的子波振幅谱,在子波相位时不变假设下采用高阶累积量的双谱法重构子波的混合相位谱,最后将振幅谱和相位谱结合实现时变混合相位子波的提取.与常规分段提取方法相比,该方法不需要子波分段平稳的假设,当相邻层段子波属性存在较大变化时也能准确估计出时变子波.理论分析和仿真实验结果表明,该方法有效可行.     

改进的统计子波反褶积

地震子波的小相位性，反射系数的白噪声以及理想的期望子波是地震资料反褶积处理需要考虑的三个主要因素。当地震子波为非最小相位时，可利用指数加权方法使之相位化；对于反射系数非白噪声问题，本文提出了利用对数功率谱来估算地震子波的以及反褶积算子的方法；对于理想的期望子波，本文采用“宽带雷克子波”作为期望输出子波。