[发明专利]一种利用叠后地震数据进行油气检测的方法

权利要求书

1.一种利用叠后地震数据进行油气检测的方法，其特征是包括以下步骤:

(1)数据准备：获取工区的地震叠加纯波数据、工区内各井位的坐标及各井目的层段含油气层厚度和解释的目的层顶、底时间数据；

(2)确定地震数据待分析时间范围和滑动时窗长度；

(3)对工区地震叠后纯波数据提取能量系数，构建地震波衰减因子,获得衰减因子数据体：

针对地震叠后纯波数据的每一地震道，在步骤(2)中确定的分析时间范围内对每一时间采样点以步骤(2)中确定的滑动时窗长度做短时窗傅立叶变换，获得该时间采样点的频谱，确定加权平均频率和有效频带的低频端和高频端；以上述步骤确定的加权平均频率为分界点，把该时间采样点的频谱划分为低频和高频两个部分，把频谱的低频部分以频率反向加权的形式求和，获得频谱的低频部分能量；把频谱的高频部分求和，获得频谱的高频部分能量；

利用上述步骤获取的低频部分能量与高频部分能量的比值，获得该时间采样点的能量系数，并通过该时间采样点的能量系数和上述步骤确定的加权平均频率，构建该时间采样点的地震波衰减因子；

(4)重复上述步骤(3)，获得工区地震叠后纯波数据每一道在步骤(2)确定分析时间范围内每一时间采样点的地震波衰减因子，得到工区地震波衰减因子数据体,根据工区地震波衰减因子数据体数值的大小，分析、识别工区内有利含油气目标区；

对工区地震叠后纯波数据提取能量系数，构建地震波衰减因子,获得衰减因子数据体，包括以下步骤：

(3.1)针对地震叠后纯波数据的每一地震道，在确定的分析时间范围内对每一时间采样点以确定的滑动时窗长度做短时窗傅立叶变换，获得该时间采样点的频谱A(i,f)，其中，i为时间采样点的序号，f为频率采样点的序号，由A(i,f)中最大振幅的1/20为下限确定有效频带的低频端f_L和高频端f_H，加权平均频率f_w通过如下方法确定

(3.2)以上述步骤(3.1)中确定的加权平均频率为分界点，把该时间采样点的频谱划分为低频和高频两个部分，把频谱的低频部分以频率反向加权的形式求和，获得频谱的低频部分能量E_L

把频谱的高频部分求和，获得频谱的高频部分能量E_H

(3.3)利用步骤(3.2)获取的低频部分能量与高频部分能量的比值，获得该时间采样点的能量系数C(i)，

通过该时间采样点的能量系数和上述步骤(3.1)确定的加权平均频率，构建该时间采样点的地震波衰减因子Q(i)

其中，b为常数，此处取值为1；

(3.4)重复上述步骤(3.1)-(3.3)，获得工区地震叠后纯波数据每一道在分析时间范围内每一时间采样点的地震波衰减因子，得到工区地震波衰减因子数据体；

(3.5)对上述步骤(3.4)获得的工区地震波衰减因子数据体，根据其数值的大小，分析、识别工区内有利含油气目标区。

2.根据权利要求1所述的一种利用叠后地震数据进行油气检测的方法，其特征是：步骤(2)的确定滑动时窗长度包括以下步骤：

(2.1)根据井位坐标在工区地震叠后纯波数据中抽取井旁地震道x(i)，i为所抽取井旁地震道时间采样点的序号；

(2.2)以不同滑动时窗长度，利用傅立叶变换方法获得地震道x(i)的目的层段每个时间采样点的频谱A(i,f)，即

A(i,f)＜--＞FFT(x'(i))

其中，i为时间采样点的序号，f为频率采样点的序号，x'(i)是地震道x(i)加海明窗的结果，即

海明窗的表达式为

其中，j为海明窗函数的时间采样点序号，N为滑动时窗长度的样点数，π为圆周率；

(2.3)分析不同时窗长度对频率域振幅谱形态及其对油气层反映的能力，确定滑动时窗长度。