统计法同态反褶积

同态反褶积在地震资料处理中已得到了广泛应用。这种非线性系统,对地震子波不做最小相位假设,只需能分离出地震子波和反射系数,在分离出的地震子波中,选取最合适的地震子波需要花费大量的人工时间,并需要一定的经验。本文采用统计法在许多反子波中按一定的准则选取最佳反子波用于同态反褶积。判断准则是采用最小熵反褶积中的MED准则和PARSIMONY准则,从而解决了人工选取的费时和由于没有统一的选择标准出现的因人而异等问题。通过大量的地震测线试验,按MED准则和PARSIMONY准则求得的最佳反子波相同。求得的地震子波不仅能用于同态反褶积还可用于地震子波的最小相位和零相位转换,转换后的最小相位地震子波能满足对地震子波做最小相位假设反褶积的要求,转换后的零相位地震子波可用于地震剖面的最终输出,符合解释要求。     

可控震源自适应地表一致性反褶积

随着地震勘探资料处理技术的发展，用于地震勘探数据处理的反褶积方法也不断完善。但是，多年来对可控震源地震数据处理主要采用纯相位滤波方法来进行零相位至小相位的子波转换。这种方法对于有Ｑ吸收的实际地震数据存在一定问题，它使可控震源相关波不再是零相位，而是混反褶积方法。该方法能自动提取地表一致Ｑ吸收量，并模拟出地表一致性Ｑ吸收可控震源相关子波，用子波反褶积方法。该方法能自动提取地表一致性Ｑ吸收量，并模拟出     

同态反褶积的改进与应用

在地震记录的褶积模型中,当假设地震子波为最小相位、反射系数为白噪声时,反褶积的不定解问题便变成了可解。由于判定实际地震记录中地震子波的相位特性很难,因此出现了几种混合相位子波的反褶积,如同态反褶积、最小熵反褶积、L_1模反褶积等,这些方法各有优、缺点。统计法同态反褶积使混合相位反褶积达到了实用化的程度,但应用中要人工选定地震子波和反子波。本文将统计法同态反褶积与其他几种混合相位反褶积方法结合起来,即对统计出的多个子波,采用MED准则、L_1模准则和PARSIMONY准则选定最佳子波和反子波,使混合相位反褶积在实用化方面又向前推进了一步。     

����������ȡ�����λ�Ӳ��ķ���

为有效地提高地震资料的分辨率，提出了一种新的混合相位地震子波提取方法。相同振幅谱、不同相位谱的地震子波分别对应不同的最大和最小相位分量。在复赛谱域，遵照一定的规则，子波振幅谱可以分割为两部分，其分别对应于子波的最大和最小相位分量，分割比例的不同决定了子波不同的相位，由此可以得到最大相位子波、最小相位子波、零相位子波和一系列混合相位子波。根据地震记录信噪比谱的情况，设定地震记录的希望输出子波，对地震记录进行确定性子波反褶积。依据反褶积后地震记录方差模大小和分辨率情况，最终从一系列相同振幅谱不同相位谱子谱集合中确定理想的混合相位子波，通过确定性子波反褶积达到提高地震记录分辨率的目标。对理论合成地震记录和实际地震记录进行了子波提取实验，并利用提取的地震子波对地震资料进行了提高分辨率的处理，实验结果表明了该方法的正确性和有效性。     

同态反褶积的子波稳定性方法研究

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

Research on the methods for wavelet stability of homomorphic deconvolution.

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

海上地震资料零相位化处理技术研究

在海洋地震资料采集中,气枪震源子波接近最小相位子波,反射信号的峰值点并不代表地下反射界面,需要进行相位处理。传统的最小相位化处理技术不能彻底消除相位对资料的影响,传统反褶积方法从地震资料中提取零相位化算子并对资料进行零相位化处理技术,由于提取算子的不稳定性,很难得到满意的处理效果。为此,探讨了采用震源模拟子波、远场实测子波和海底提取子波求取零相位化算子并对地震资料进行零相位化处理的技术。首先分析了非零相位处理对后续高分辨率处理和解释标定的影响;然后讨论了海洋地震资料子波的获取方式,阐明了在深水、浅水两种条件下如何利用不同子波进行零相位化算子求取和应用;最后以海洋资料处理为例,证实了深、浅水资料零相位化处理技术的可靠性。     

因式分解串联反褶积

当前使用较多的最小相位反褶积方法均假设子波为最小相位，所以不适应非最小相位的实际地震子波，限制了地震数据处理中分辨率的提高。即使是实际地震子波为最小相位，采用最小相位反褶积方法，也会导致较多的高频噪声，使信噪比降低。本文提出的因式分解串联反褶积方法避开子波的最小相位假设，通过对自相关进行因式分解，得到子波的可能因式，然后对分解所得到的子波因式逐个求反褶积算子，并进行最大、最小相位筛选，再将各个因式串联进行试探和判别，最后实现非最小相位子波的反褶积。对理论数据和实际数据的应用表明，因式分解串联反褶积方法能显著提高分辨率，且在一定程度上提高了信噪比，即使是在最小相位子波情况下，其效果也优于当前使用的最小相位反褶积方法。     

零相位校正和子波反褶积

熵，作为描述一个序列有序程度的量，从热力学引入地震处理中，用以描述反射系数的随机程度，并把它作为一种约束条件应用到反褶积技术中。在陆相薄互层地质结构中，地层在局部不能表现出较强的反射值，但在整体上具有强反射界面的特点，在应用最小熵反褶积时，要有强反射界面作为约束。在子波反褶积和零相位校正中，子波相位求取在约束条件下是个最优化问题。本文研究了最小熵反褶积的应用条件，根据陆上地震勘探中子波延迟较短的特点，改进了前人在应用模拟退火法中对子波长度已知的假设，简化了输入条件，在最小熵的约束下完成了零相位校正和子波反褶积。     

频率域地表一致性反褶积方法及应用效果分析

提高地震资料分辨率的重要手段是采用反褶积技术。反褶积方法很多，如子波反褶积、脉冲反褶积和预测反褶积等。这些方法各有优缺点，在高分辨率资料处理中应用受到许多条件制约，不能有效地提高资料的分辨率，其效果不能满足解释人员的要求。作频率域地表一致性反褶积需要对记录进行频谱分析和频谱分解，同时要设计反褶积算子，然后在共炮点域和共接收点域分两步对资料进行褶积。这种方法能够展宽频谱，压缩地震子波，并能校正地震信号的相位谱，输出零相位子波，较大程度地提高地震资料的分辨率。该方法在胜利油田车22地区三维等资料处理中应用，取得了较好的效果。     

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

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

相位校正S域反褶积方法及应用

反Q滤波和时频（S）域反褶积都可以提高衰减地震记录的分辨率。通常指数增益导致反Q过程不稳定;而时频域反褶积虽然可以很好地适应地震记录的非平稳性,但只能去除子波振幅谱的影响,无法准确提取子波相位谱,从而影响反褶积结果。为此,引入反Q滤波的相位校正方法,使衰减地震记录零相位化,优化反褶积的基础条件;再通过改进的广义S变换在S域实现反褶积,达到高分辨处理的目的。模型数据测试和辽西实际数据应用的结果均表明,相位校正S域反褶积能得到更准确的结果,是识别薄互层和刻画地层尖灭点的有效方法。     

三维叠后零相位同态反褶积

在提高分辨率的处理中,人们比较熟悉叠前统计子波反褶积和地表一致性反褶积,它们在提取地表一致性信息和提高分辨率方面起到了十分有效的作用.然而在叠后处理中往往采用单道反褶积方法,但是,从高保真角度讲,采用二维和三维叠后多边或一体化反褶积将更加必要.本文首先应用自组织神经元聚类和分形果类方法求出对应地层组的三维一体化子波,然后完成叠后三维零相位固态反褶积处理,从而达到保持地质异常作信息、提高分辨率的目的.     

一种基于自然梯度的地震盲反褶积方法

常规反褶积方法在地震子波和反射系数未知的情况下，需要做统计性假设，这样往往会降低地震数据处理的分辨率和可靠性。基于自然梯度的盲反褶积算法无需任何假设，利用自然梯度法估算地震子波，得到更接近于实际情况的混合相位子波，用于反褶积，使地震资料的分辨率得到了提高。首先介绍了源数据体可分离的数学依据；然后阐述了基于自然梯度盲反褶积算法的数学原理，并引入补偿函数，推导出适合地震盲反褶积的迭代公式；最后分别利用合成地震记录和实际地震数据对方法进行了验证，结果表明，基于自然梯度的盲反褶积算法能较好地适用于混合相位子波，最大可能地估算原始反射系数。     

最大熵子波反褶积

最大熵子波反褶积是建立在最大熵反褶积的基础上。最大熵反褶积的优点在于不需要求子波的自相关值和期望输出与子波的互相关值,即不需要解托布里兹方程组求出反滤波因子,而只要求出各i阶偏相关系数ρ_i(i=1,…,M)就可由数据直接递推地求出第M阶向前预测误差和向后预测误差。这种作法可以避免在数据时窗外用补零的办法求相关值而带来的误差。最大熵反褶积是在地震子波为最小相位条件下进行的,其向前预测误差就等于反射系数序列。最大熵子波反褶积与最大熵反褶积的不同点在于,该方法将所得反射系数序列作为初值,用此反射系数序列与地震记录作相关提取单道子波,再采用多道统计平均法求取最终子波。为使其最小相位化,对最终子波要进行整形,然后再作子波反褶积。经模型试验与实际资料处理表明,该方法效果良好。     

纯振幅反褶积

在未知子波的情况下,可由地震剖面直接求取零相位子波,由零相位子波求出反褶积因子,并对地震记录作反褶积处理。纯振幅反褶积与其它类型反褶积不同,它不要求子波为最小相位,它的反褶积因子为零相位的对称波形,褶积结果不改变原剖面的相位特征,即没有时移。文中叙述了纯振幅反褶积因子的求法及实现高分辨率、高信噪比的可行性。     

基于零偏移距VSP的时变子波反褶积方法

提高地面地震数据分辨率的一种有效方法是从VSP数据中提取反褶积算子并将其应用于地面地震数据,然而,传统的VSP子波反褶积方法通常是在不考虑地震子波时变特性的情况下,提取固定反褶积算子,且由于地震子波是时变的,从VSP中提取一个子波反褶积算子应用于地面地震记录是不准确的,由此提出一种利用VSP数据提取时变子波反褶积算子的方法。该方法将零偏移距VSP数据的下行直达波作为单程VSP地震子波,利用地下介质的衰减函数将单程地震子波的振幅谱转换为双程地震子波的振幅谱,在最小相位的假设下,从时变子波振幅谱中获得时变反褶积算子,然后将算子应用于叠后地面地震数据。理论合成数据和实际地震数据证明了该方法可有效补偿振幅能量并可提高地震记录的分辨率。