[发明专利]基于变子波反射系数反演的地震高分辨处理方法

说明书

技术领域

本发明涉及石油天然气地球物理勘探中的地震资料优化方法，尤其涉及基于变子波反射系数反演的地震高分辨处理方法。

背景技术

随着油气勘探的不断深入，很多储集体，特别是地层的储集体体积较小，识别上述储集体则必须提高地震的分辨率。

常用的提高地震时间分辨率的方法有反褶积处理，地震谱白化处理，以及频率域拓频处理技术。上述方法都是在地震资料的有限频带内增加地震的频率成分，而地震频带的带宽却没有拓宽。

近年来出现的地震谱反演方法是基于地震记录的褶积模型反演反射系数，现有的反射系数反演大多基于静态的褶积的模型，实际上大地是个滤波器，地震波在传播的过程中会经受吸收衰减的影响，不同介质的吸收作用又是不同的，因此基于地震波吸收衰减的地震波传播模型更符合实际，这种模型比静态褶积模型更精确地反应地震波在地下的传播规律。

因此，在如何通过拓宽地震频带来提高地震分辨率的问题上，仍然存在很大的研发空间。

发明内容

本发明为解决现有技术中，地震分辨率的表示精度较低的问题，提供了一种基于变子波反射系数反演的地震高分辨处理方法。

本发明的技术方案如下；

基于变子波反射系数反演的地震高分辨处理方法，所述处理方法为拓展地震记录的频带的方法，用于提高地震记录的时间分辨率；处理方法包括先根据地震记录提取浅层地震子波、子波时间域长度；设定时窗并提取每个时窗内的瞬时子波；再通过时间域的子波褶积模型建立频率域的变子波褶积模型；对每一道地震记录，将所述瞬时子波和变子波褶积模型建立模型参数，反演出反射系数；最后将初始提取的浅层地震子波与反射系数褶积处理，获取频带拓宽的地震记录。

所述处理方法的具体步骤为，

步骤1，获取浅层地震资料；

步骤2，获取浅层地震子波以及浅层地震子波时间域长度：

步骤3，设置时窗：

以所述步骤2的浅层地震子波时间域长度的1.5～2倍设置时窗，将地震记录按所述时窗长度均匀分段；

步骤4：提取瞬时子波：

设各所述时窗内地震记录的平滑振幅谱为利用公式（1）提取所述步骤3中每个所述时窗内的变子波，即提取每个所述时窗内中心时刻的瞬时子波，提取过程如公式1所示；

wTi(t)=ifft(wT(f))---1;]]>

其中，ifft为反傅立叶变换算子，T_i为所述地震道上的第i个时窗；

时窗内的变子波通常为最小相位子波，零相位子波以及混合相位子波，常用的是零相位子波和最小相位子波。

步骤5，建立频率域的变子波褶积模型：

步骤5-1，分解反射系数脉冲序列：

设反射系数序列为a(k),k=0,1,2...K,将所述反射系数序列a(k)分解为K个脉冲之和，分解过程如公式2所示：