[发明专利]油气分析方法及装置

权利要求书

1.一种油气分析方法，其特征在于，所述方法包括：

获取待分析区域多道地震信号；

针对每道地震信号，对所述地震信号进行预处理获得该地震信号的多个短时功率谱，所述多个短时功率谱是通过不同的时间窗口获得；

将所述短时功率谱映射到Bark域，获得所述短时功率谱对应的多段临界频带；

针对每段所述临界频带，计算该临界频带的临界频带谱；

利用等响度曲线对所述临界频带谱的幅值进行缩放处理，获得预加重系数；

根据全极点模型对所述预加重系数进行变换，获得线性预测参数；

提取每道地震信号的所述线性预测参数的谱峰值并判断待分析区域的油气分布。

2.根据权利要求1所述的方法，其特征在于，所述对所述地震信号进行预处理获得该地震信号的多个短时功率谱的步骤包括：

去除所述地震信号中的噪声，获得第一中间信号；

对所述第一中间信号进行加窗处理，获得每个时间窗口内的第二中间信号；

针对每个所述第二中间信号进行离散傅里叶变换，获得该第二中间信号的短时功率谱。

3.根据权利要求1所述的方法，其特征在于，所述将所述短时功率谱映射到Bark域的计算公式如下：

z₁(w)＝6ln{w/1200π+sqrt[(w/1200π)²+1]}

其中，w代表频率。

4.根据权利要求1所述的方法，其特征在于，所述利用等响度曲线对所述临界频带谱进行加重处理，获得预加重系数的步骤为：

计算所述临界频带谱与该临界频带谱对应的等响度曲线之积，获得预加重系数；

计算等响度曲线的公式如下：

E₁(w)＝[(w²+56.8*10⁶)w⁴]/[(w²+6.3*10⁶)²*(w²+0.38*10⁹)]

其中，E₁(w)为等响度曲线，w代表频率。

5.根据权利要求1所述的方法，其特征在于，所述根据全极点模型对所述预加重系数进行变换，获得线性预测参数的步骤前，所述方法还包括：

对所述预加重系数进行压缩处理，获得新的预加重系数；

根据全极点模型对所述新的预加重系数进行变换，获得所述线性预测参数。

6.一种油气分析装置，其特征在于，所述装置包括第一获取模块、处理模块、映射模块、计算模块、加重模块、第二获取模块和判断模块；

所述第一获取模块，用于获取待分析区域多道地震信号；

所述处理模块，用于针对每道地震信号，对所述地震信号进行预处理获得该地震信号的多个短时功率谱，所述多个短时功率谱是通过不同的时间窗口获得；

所述映射模块，用于将所述短时功率谱映射到Bark域，获得所述短时功率谱对应的多段临界频带；

所述计算模块，用于针对每段所述临界频带计算该临界频带的临界频带谱；

所述加重模块，用于利用等响度曲线对所述临界频带谱的幅值进行缩放处理，获得预加重系数；

所述第二获取模块，用于根据全极点模型对所述预加重系数进行变换，获得线性预测参数；

所述判断模块，用于提取每道地震信号的所述线性预测参数的谱峰值并判断待分析区域的油气分布。

7.根据权利要求6所述的装置，其特征在于，所述处理模块具体用于：

去除所述地震信号中的噪声，获得第一中间信号；

对所述第一中间信号进行加窗处理，获得每个时间窗口内的第二中间信号；

针对每个所述第二中间信号进行离散傅里叶变换，获得该第二中间信号的短时功率谱。