[发明专利]一种基于二阶水平多重同步挤压变换的油气储层识别方法

权利要求书

1.一种基于二阶水平多重同步挤压变换的油气储层识别方法，其特征在于，包括以下步骤：

S1、输入待分析的地震信号x(t)，t为时间；

S2、建立高斯调制线性chirp模型，用于描述具有缓慢或快速变化线性群时延的地震信号，其信号模型的具体表达式为：

其中为实数，以下公式中的均为该含义。

S3、通过类似于同步挤压变换(SST)中瞬时频率估计的方法，在时频域推导出基于上述模型的二阶局部群时延估计

其中，η为频率，为群时延估计算子，为局部频率估计算子，为局部复调制算子，其具体形式为：

其中

S4、根据多重同步挤压原理，在时频域构建一种以信号二阶群时延为中心的水平多重同步挤压算子用于挤压原时频谱群时延曲线处的时频系数，得到二阶水平多重同步挤压变换值Ts^N(t，η)；

S5、对Ts^N(t，η)求模，得到二阶水平多重同步挤压变换的地震时频谱，通过对比高低频能量衰减达到识别油气储层的目的。

2.根据权利要求1所述的基于二阶水平多重同步挤压变换的油气储层识别方法，其特征在于，所述步骤S2中高斯调制线性chirp模型的明确表达式为：

其中

信号x(t)的傅里叶变换为X(η)，可被定义为X(η)＝A(η)e^i2πθ(η)。X(η)在每个频率点处局部都可近似为一个高斯调制线性chirp模型，其中调幅部分A(η)和调频部分θ(η)分别表示为：

其中

3.根据权利要求1所述的基于二阶水平多重同步挤压变换的油气储层识别方法，其特征在于，所述步骤S3通过类似于同步挤压变换(SST)中瞬时频率估计的方法，在时频域推导出基于上述模型的二阶局部群时延估计

其中，为群时延估计算子，为局部频率估计算子，为局部复调制算子，其具体形式为：

其中

4.根据权利要求1所述的基于二阶水平多重同步挤压变换的油气储层识别方法，其特征在于，所述步骤S4根据多重同步挤压原理，在时频域构建一种以信号二阶群时延为中心的水平多重同步挤压估计用于挤压原时频谱群时延曲线处的时频系数，从而获得新的时频系数；

其中，水平多重同步挤压估计的具体表达式如下：

此时水平多重同步挤压估计有如下的收敛序列：

当N→+∞时，收敛于TF不动点。

经过二阶水平多重同步挤压变换的地震时频谱为：

其中δ(t)为单位脉冲函数，表示SHMSST时频域上的二阶多重局部群延迟估计算子。