[发明专利]一种数据驱动的地震信号零相位反褶积方法

说明书

本发明公开了一种数据驱动的地震信号零相位反褶积方法，包括：1)从地震记录中估计地震子波振幅谱；2)基于估计的振幅谱，对地震记录在频率域进行维纳滤波；3)对维纳滤波器做逆傅里叶变换，得到反褶积结果。本发明利用子波振幅谱估计方法，不依赖于模型参数，能够较好地对子波振幅谱进行有效地估计，具有好的鲁棒性。用该方法对地震记录进行零相位反褶积，可以收到了很好的效果，同时该方法也可以有效地处理非平稳地震数据。

技术领域

本发明属于地震勘探技术领域，涉及一种数据驱动的地震信号零相位反褶积方法。

背景技术

对高分辨率地震资料的要求除了与常规资料相比具有高的纵向分辨率、宽的频带外，还应该具有较高的保真度，即要求处理过的地震资料不仅能解释精细的地层构造，还必须满足后续岩性预测及储层反演的需要，这就需要较为精确地描述出地震子波在地层中的变化。所以在较短的数据集里、较低的信噪比下有效估计子波振幅谱是解决上述问题的关键。

理论上估计子波振幅谱的方法为：对一段平稳的地震信号，如果反射系数统计独立，具有有限的二阶矩，地震信号相关函数的傅里叶变换即为子波振幅谱的平方。但是地层反射系数概率分布所呈现出来的复杂性和地震信号本身的特点，使得上述方法在应用中受到了一些限制：地层的反射系数有时并不表现出独立特性，即并不是“白的”；反射系数的概率分布更接近于α-stable分布，该分布二阶以上的统计量不存在；对于低信噪比的地震资料，该方法也难以达到满意的结果。

谱模拟方法可以克服以上缺点，是一种非常有效的估计子波振幅谱方法。谱模拟是零相位反褶积的常用方法，模拟指的是从地震信号振幅谱中拟合出子波的振幅谱。

一般情况下，由于地震信号的非平稳性，地震子波不仅是时变而且是空变的，要获得准确的相位信息非常困难。在实际操作中，子波相位的不准确会引起反褶积后的剖面同相轴无规律的时移，处理后的结果会形成很多假象。地震信号零相位褶积后只能使振幅谱发生了变化，而相位谱则与原来相同，反褶积剖面上的相位特征得以保持，同相轴无时移现象。零相位反褶积所要求的就是准确地估计出子波的振幅谱，目前通常采用的方法主要有以下两种：多道振幅谱平均法和谱模拟方法。

通过将多道地震记录对数谱平均的方法，可以得到对子波振幅谱的估计。理论上讲，该方法不足之处在于：子波是空变的，各道的子波存在差异；其理论基础只是针对无噪的地震记录，未考虑到噪声对上述渐进特性的影响；地下介质的属性可能表现出不同范围的连续性，为保证效果，用于叠加的地震道应该尽量不相关，所以随机选取地震道的准则对估计质量也有影响。

Rosa和Ulrych(1991)提出一种谱模拟(spectral modeling)方法，该方法将地震子波谱视为频率的连续、光滑的函数，将其建模为一个参数待定的函数，该公式被广泛应用于子波振幅的估计，其主要的不足是估计结果对模型参数比较敏感。

Jinghuai Gao等(2017)通过对子波振幅谱的研究，提出了谱模拟的一个新的迭代算法，该算法不依赖于模型参数，能够数据驱动地、较好地对子波振幅谱进行估计，对噪声具有好的鲁棒性，该方法较使用特定的函数对地震记录谱进行拟合相比，减少了人为因素的影响，能够更客观地反映谱的变化；同时，用该方法处理短时窗的地震数据，也能收到了好的估计效果，这意味着该方法可以对非平稳地震记录中子波的变化进行很好地刻画。

由于反射系数的概率分布非常复杂、类型多样，很多统计性反褶积算法对噪声较为敏感，这使得基于谱模拟的零相位反褶积技术成为较稳健的提高分辨率的方法，该方法要求准确地估计出子波的振幅谱。传统的方法是基于多项式函数拟合子波的振幅谱，模型参数较多、估计结果受参数影响很大。

发明内容

本发明的目的在于克服上述现有技术的缺点，使用新的谱模拟技术，对地震信号进行零相位反褶积，提供一种数据驱动的地震信号零相位反褶积方法。

为达到上述目的，本发明采用以下技术方案予以实现：