[发明专利]一种岩性识别方法及识别系统

权利要求书

1.一种岩性识别方法，其特征在于，包括如下步骤：

S1、获取测井资料；

S2、获得岩性曲线，以及根据所述测井资料计算获得波阻抗曲线A；

S3、根据所述岩性曲线获得至少一个岩性的岩性敏感曲线B；

S4、对所述波阻抗曲线A进行低通滤波处理，通过公式(1)计算得到低频趋势波阻抗A_lpf，以及对所述敏感曲线B进行高通滤波处理，通过公式(2)计算得到高频敏感曲线B_hpf；所述公式(1)为：A_lpf＝F(A)，所述公式(2)为：B_hpf＝F(B)；其中，F(A)、F(B)均为傅里叶变换函数；

以及S5、根据所述低频趋势波阻抗A_lpf与高频敏感曲线B_hpf，首先，将所述高频敏感曲线Bhpf转化到高频波阻抗Ahpf值域范围内，得到高频敏感波阻抗Chpf；然后，在频率域内根据公式(5)对所述高频敏感波阻抗Chpf与高频波阻抗Ahpf进行加权计算，计算得到岩性识别敏感因子；最后，通过所述岩性识别敏感因子进行地震反演，以对岩性进行识别；所述公式(5)为：

S＝f[(F(Ahpf)+F(Chpf)]；

其中，F(x)为傅里叶变换函数，f(x)为傅里叶逆变换函数。

2.根据权利要求1所述的岩性识别方法，其特征在于，步骤S1中，所述测井资料包括：声波时差、密度、自然伽马、自然电位中的一种或几种。

3.根据权利要求1所述的岩性识别方法，其特征在于，所述步骤S5包括：

S51、统计高频波阻抗A_hpf值域变化范围，得到其最小值A_min和最大值A_max，以及统计高频敏感曲线B_hpf值域变化范围，得到其最小值B_min和最大值B_max；

S52、将所述高频敏感曲线B_hpf与所述高频波阻抗A_hpf进行交汇分析，根据公式(3)拟合求取相关函数Fx(B_hpf)，所述公式(3)为：Fx(B_hpf)＝a*(B_hpf)²+b*(B_hpf)+c；其中，所述a、b、c均为相关常数；

S53、根据公式(4)将所述高频敏感曲线B_hpf转化到所述高频波阻抗A_hpf值域范围内，得到高频敏感波阻抗C_hpf；所述公式(4)为：

以及S54、在频率域内根据公式(5)对所述高频敏感波阻抗C_hpf与高频波阻抗A_hpf进行加权计算，由此得到岩性敏感因子S。

4.根据权利要求3所述的岩性识别方法，其特征在于，所述步骤S52中，所述高频敏感曲线B_hpf与高频波阻抗A_hpf的相关系数大于或等于0.7。

5.一种实现权利要求1-4任一项所述岩性识别方法的岩性识别系统，其特征在于，包括：

第一曲线获取单元，其用于获取岩性曲线，以及根据所述测井资料计算获得波阻抗曲线A；

第二曲线获取单元，其用于根据所述岩性曲线获得至少一个岩性的岩性敏感曲线B；

第一曲线处理单元，其用于对所述波阻抗曲线A进行低通滤波处理，计算得到低频趋势波阻抗A_lpf；

第二曲线处理单元，其用于对所述敏感曲线B进行高通滤波处理，计算得到高频敏感曲线B_hpf；

以及敏感因子获取单元，其用于根据所述低频趋势波阻抗A_lpf与高频敏感曲线B_hpf计算得到岩性识别敏感因子。

6.根据权利要求5所述的岩性识别系统，其特征在于，还包括：数据输出单元，其用于根据所述岩性识别敏感因子输出岩性统计数据。