[发明专利]一种突出含气储层低频地震异常的地震信号时频分析方法

权利要求书

1.一种突出含气储层低频地震异常的地震信号时频分析方法，其包括以下主要步骤：

(1)输入地震数据集，对数据集中的每一道地震记录x(t)进行傅里叶变换，得到其频谱x(f)，其中，x(t)包含N个采样点，t为时间，单位是秒，f为频率，单位是赫兹；

(2)给定低频控制参数σ_L和高频控制参数σ_H，按照下式构建第j个频率f_j对应的分析窗：

式中，j＝0,1,…,N-1，

(3)按照下式对w(f_j,τ_a,τ_b)进行能量归一化处理：

式中，δ(j-f)是单位脉冲函数，*表示褶积运算；

(4)计算能量归一化分析窗与x(f)的哈达玛积，即：

式中，⊙表示两个向量的哈达玛积运算，[·]^T表示列向量；

(5)重复步骤(2)至(4)，直至对所有频率值都进行了相应的计算，最终得到N阶矩阵X(f′,f)＝[X(f₀,f),…,X(f_j,f),…,X(f_N-1,f)]；

(6)沿着矩阵X(f′,f)的第一频率轴f′的方向进行反傅里叶变换，得到x(t)的时频分析结果X(t,f)；

(7)基于X(t,f)，在选定频率范围[f_a,f_b]内，对该频率范围内的每个频率值f_j，计算其对应的流度曲线y_j(t)：

其中，2＜a＜j＜b＜N/2；

(8)分别抽取每一道地震记录x(t)的低、中、高频流度曲线y_L(t)、y_M(t)、y_H(t)，形成低、中、高频流度数据集Y_L、Y_M、Y_H；

(9)根据需求，分别从Y_L、Y_M、Y_H中抽取低、中、高频流度剖面，并对比这些剖面中振幅能量异常随频率的变化，若振幅能量异常随频率的增加而消失，则将剖面上对应的位置标记为油气储层。