[发明专利]一种广义W变换地震含气储层表征方法

权利要求书

1.一种广义W变换地震含气储层表征方法，其特征在于，包括以下步骤：

S1、输入待分析的过井地震道信号s(t)；其中t表示时间；

S2、对S1输入的过井地震道信号s(t)，选择合适的参数集Γ＝(k,p)范围，估计各个参数集Γ下的广义高斯窗口g(t,f；Γ)，计算方法如下：

其中，t表示时间中心，f表示频率中心，Γ＝(k,p)为一组可调节的参数集，k表示尺度因子，调节广义高斯窗口的宽度，p表示趋势因子，调节广义高斯窗口的变化率所述合适的参数集范围为Γ_range＝{k∈(1，3),p∈(1，5]}；σ(f₀，f；Γ)为标准差函数：

其中，f₀为过井地震道信号s(t)的瞬时主频，根据Hilbert变换估计获得；

S3、根据上述获得不同参数集下Γ_range的广义高斯窗口g(t,f；Γ)，对过井地震道信号s(t)进行分解，获取时频变换结果GWT(t,f；Γ)，计算方法如下：

S4、采用聚焦策略(CM)，对S2中广义高斯窗口函数g(t,f；Γ)的参数集范围Γ_range＝{k∈(1,3),p∈(1,5]}内进行自适应优化，得到广义W变换的CM：

其中表示广义W变换的归一化形式，并选取CM值最大的参数集作为优化参数集Γ_opt：

最终得到优化参数集Γ_opt下的广义W变换的时频谱；

S5、采用傅里叶变换计算过井地震道信号的有效频带范围，将有效频带范围中的频率最大值的75％、60％、50％、33％设置为对比分析的频率值。通过对比分析S4得到的时频谱中上述四个频率值的衰减情况来表征该井是否存在含气储层。

2.根据权利要求1所述的广义W变换地震含气储层表征方法，其特征在于，所述步骤S2中过井地震道信号s(t)在各个参数集下Γ＝(k,p)的广义高斯窗口为：

其中，t表示时间中心，f表示频率中心，Γ＝(k，p)为一组可调节的参数集，k表示尺度因子，调节广义高斯窗口的宽度，p表示趋势因子，调节广义高斯窗的变化率，合适的参数集范围为：Γ_range＝{k∈(1,3),p∈(1,5]}；σ(f₀，f；Γ)为标准差函数：

其中，f₀为过井地震道信号s(t)的瞬时主频。

3.根据权利要求1所述的广义W变换地震含气储层表征方法，其特征在于，所述步骤S3中，根据上述获得不同参数集Γ_range下的广义高斯窗口g(t,f；Γ)，对过井地震道信号s(t)进行分解，获取时频变换结果GWT(t,f；Γ)：

4.根据权利要求1所述的广义W变换地震含气储层表征方法，其特征在于，所述步骤S4可采用下式对(8)中的GWT(t,f；Γ)进行逆变换，重构地震信号s(t)：

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