[发明专利]一种基伸缩调频同步提取地震储层预测方法

权利要求书

1.一种基伸缩调频同步提取地震储层预测方法，其特征在于，包括以下步骤：

S1、输入待分析的原始二维地震剖面信号s(x；t)；其中x表示地震道，t表示时间；

S2、对S1输入的二维地震剖面信号s(x；t)进行基伸缩自适应调频变换，获取时频变换结果SBCT(x；f,t)，计算方法如下：

其中，t表示时间中心，f表示频率中心，h(x；τ)表示高斯函数，表示相位函数，τ为时间变量；

S3、根据S2得到的时频变换结果SBCT(x；f,t)估计第x道地震信号在各个时频位置(t,f)处的瞬时频率

S4、根据同步提取原理，在该时频域上构建一种以信号瞬时频率曲线为中心的同步提取算子SEO(x；f,t)用于提取原时频谱能量，获得新的时频系数T(x；f,t)，剔除掉模糊的时频能量；

S5、对S4获得的时频系数取模，得到基伸缩调频同步提取变换后的时频谱；

S6、采用傅氏变换计算原始地震信号s(x；t)的主频范围，将主频范围中的频率最大值的85％设置为高频值，将主频范围中的频率最小值设置为低频值，然后从S5得到的时频谱中抽取高频值对应的频谱，得到高频值的共频率剖面；从S5得到的时频谱中抽取低频值对应的频谱，得到低频值的共频率剖面。通过对比两个共频率剖面中地震信号的衰减情况来确定该地震剖面中是否存在储层。

2.根据权利要求1所述的基伸缩调频同步提取地震储层预测方法，其特征在于，所述步骤S2中二维地震剖面信号s(x；t)的基伸缩调频自适应变换的时频变换结果为：

其中，t表示时间中心，f表示频率中心，h(x；τ)表示高斯函数，表示相位函数，τ为时间变量，所述相位函数为：

其中，k＝1,2,…,n，n为相位函数的阶数，(a₁,a₂,…,a_n)是相位函数的自适应参数，可根据峰度理论确定，即：

其中，argmax(·)表示函数取到最大值时的参数集。

3.根据权利要求1所述的基伸缩调频同步提取地震储层预测方法，其特征在于，所述步骤S3中瞬时频率的计算包括根据S2得到的时频变换结果SBCT(x；f,t)估计信号在各个时频位置(t,f)处的瞬时频率

其中，表示窗函数h(x；t)对时间t的偏导，表示在窗函数下的SBCT结果；表示在函数t^kh(x；t)下的SBCT结果。

4.根据权利要求1所述的基伸缩调频同步提取地震储层预测方法，其特征在于，所述步骤S4中获取新的时频系数T(x；f,t)的具体方法为：根据同步提取原理，在该时频域上构建一种以信号瞬时频率为曲线中心的同步提取算子SEO(x；f,t)提取原时频谱能量获得新的时频系数T(x；f,t)，剔除掉模糊的时频能量，同步提取算子SEO(x；f,t)满足：

其中δ(·)为单位脉冲函数，新的时频系数T(x；f,t)为：

T(x；f,t)＝SBCT(x；f,t)SEO(x；f,t) (7) 。

5.根据权利要求1所述的基伸缩同步提取地震储层预测方法，其特征在于，所述步骤S5可采用下式对(7)中的T(x；f,t)进行逆变换，重构地震信号s(x；t)；

其中，本发明适用于重构地震信号。