[发明专利]一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法

权利要求书

1.一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法，其特征在于，包括以下几个步骤：

步骤101，使用固定时长的移动时窗截取短时地震信号，对每一个移动时窗内的短时地震信号做傅里叶变换以获得其复数信号，分别提取其实数部和虚数部，通过反傅里叶变换分别处理实数部和虚数部，对应的得到移动时窗内原始信号在时间域上的偶函数和奇函数部分；

步骤102，提取井旁的地震子波，计算所述移动时窗内信号奇、偶函数部分的振幅谱，同时计算已提取的地震子波的振幅谱，用奇、偶函数部分的振幅谱除以地震子波的振幅谱，得到标准化后的奇、偶函数部振幅谱；分别计算奇、偶函数部的峰值振幅属性；

步骤103，使用获得的奇、偶函数峰值振幅属性，与地震相关属性特征共同组成多属性数据集，采用多属性分析并与井下实测的储层孔隙度数据拟合，以获得设定范围内储层孔隙度预测结果。

2.根据权利要求1所述的基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法，其特征在于，步骤101具体的方法如下：

对地震尺度上的薄储层而言，由于入射地震波在顶底反射叠加效应，将叠加形成的复合反射作为一对脉冲信号；定义所述薄储层的地震反射系数脉冲对g为：

g(t)＝r1δ(t-T/2)+r2δ(t+T/2)， (1)

其中，r1，r2分别为地层顶、底反射系数；T是薄储层的时间厚度，δ表示狄拉克δ函数，t代表时间；

获得移动时窗内原始信号的偶函数部分和奇函数部分的方法为：

定义反射脉冲对的奇函数部分ge和偶函数部分go分别为：

ge(t)＝reδ(t-T/2)+reδ(t+T/2)， (2a)

go(t)＝roδ(t-T/2)+roδ(t+T/2)， (2b)

其中，re、ro分别代表薄层反射脉冲对g的偶数和奇数部反射系数；

101a)对反射系数脉冲对g做傅里叶变换，将其转换至频率域为G(f)：

G(f)＝(2re)cos(πfT)+i(2ro)sin(πfT)； (3)

其中，re，ro分别为反射系数脉冲对偶数部、奇数部的反射系数大小；f为频率；

101b)分别取其实数部和虚数部，分别获得：

Re[G(f)]＝(2re)cos(πfT) (4a)

Im[G(f)]＝(2ro)sin(πfT) (4b)

即反射率频谱的实数部即等于其偶函数的频谱，反射率频谱的虚数部即等于其奇函数的频谱；

101c)则反射系数脉冲对g在时间域的偶函数和奇函数分别由取反射系数脉冲对频谱的实数部和虚数部并做反傅里叶变换得到：

g e = 1 2 π ∫ Re [ G ( ω ) e i ω t ] d ω - - - ( 5 a ) ]]>

g o = 1 2 π ∫ Im [ G ( ω ) e i ω t ] d ω - - - ( 5 b ) ]]>

其中，公式右侧ω、e分别为角频率和自然常数。