本发明涉及地质勘测领域，具体涉及一种岩性识别方法及识别系统。
背景技术：
随着对岩性油气藏的深入研究，如何实现复杂岩性储层的精准识别成为勘探开发过程中亟待解决的主要问题。基于地震属性分析、地震反演等储层预测技术在预测岩性及物性研究中发挥了重要的作用，波阻抗反演成为较为常用的储层预测手段。但由于复杂岩性带的不同岩性的波阻抗存在较大重叠段，因此造成直接应用波阻抗识别岩性边界较为困难，不能达到精准预测储层要求。
技术实现要素：
针对上述缺陷，本发明提供了一种岩性识别方法及识别系统，其可以对原始波阻抗进行岩性识别能力改造，提高其识别岩性精度，最大发挥地震储层预测能力。
本发明为解决上述技术问题所提供的技术方案如下：。
与现有技术相比，本发明的有益效果是：本发明可以对原始波阻抗进行岩性识别能力改造，提高其识别岩性边界的精度，最大程度发挥地震储层预测能力。
附图说明
图1为实施例一中岩性识别方法的步骤流程图；
图2为现有技术中使用原始波阻抗曲线产生的岩性统计图；
图3是实施例一中岩性识别方法中输出的岩性统计图；
图4是实施例二中岩性识别系统的结构示意图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
实施例一：
图1示出了本发明的1、一种岩性识别方法，其特征在于，包括如下步骤：
s1、获取测井资料；本实施例中，测井流程以及所用设备均为现有技术，在此不在赘述；且所述测井资料包括：声波时差、密度、自然伽马、自然电位中的一种或几种；
s2、获取岩性曲线，以及根据所述测井资料计算获得波阻抗曲线a(如通过声波时差、密度计算获得)；
s3、根据所述岩性曲线获得(如通过统计直方图和交汇图分析等方式)至少一个岩性的岩性敏感曲线b；
s4、对所述波阻抗曲线a进行低通滤波处理，计算得到低频趋势波阻抗alpf，以及对所述敏感曲线b进行高通滤波处理，计算得到高频敏感曲线bhpf；
s5、根据所述低频趋势波阻抗alpf与高频敏感曲线bhpf计算得到岩性识别敏感因子，并通过所述岩性识别敏感因子进行地震反演，以对岩性进行识别；
以及s6、输出岩性统计图。
具体的，步骤s4中，通过公式(1)计算获得所述低频趋势波阻抗alpf，所述公式(1)为：alpf＝f(a)；通过公式(2)计算获得所述高频敏感曲线bhpf，所述公式(2)为：bhpf＝f(b)。
此外，所述步骤s5包括：
s51、统计高频波阻抗ahpf值域变化范围，得到其最小值amin和最大值amax，以及统计高频敏感曲线bhpf值域变化范围，得到其最小值bmin和最大值bmax；本实施例中，高频波阻抗ahpf值域变化范围为(amin,amax)＝(7275,13433)，高频敏感曲线bhpf值域变化范围为(bmin,bmax)＝(39,120)；
s52、将所述高频敏感曲线bhpf与所述高频波阻抗ahpf进行交汇分析，根据公式(3)拟合求取相关函数fx(bhpf)，所述公式(3)为：fx(bhpf)＝a*(bhpf)2+b(bhpf)+c；其中，所述a、b、c均为相关常数，可通过关系函数拟合求得，本实施例中，a＝-1.14572,b＝169.313,c＝5155.07；此外，为保证数据的精确度，所述高频敏感曲线bhpf与高频波阻抗ahpf的相关系数大于或等于0.7(本实施例中优选为0.7)；
s53、根据公式(4)将所述高频敏感曲线bhpf转化到所述高频波阻抗ahpf值域范围内，得到高频敏感波阻抗chpf；所述公式(4)为：
以及s54、在频率域内根据公式(5)对所述高频敏感波阻抗chpf与高频波阻抗ahpf进行加权计算，由此得到岩性敏感因子s；所述公式(5)为：
s＝f[(f(alpf)+f(chpf)]；
其中，f(x)为傅里叶变换函数，f(x)为傅里叶逆变换函数。
图2为现有技术中，使用原始波阻抗曲线产生的岩性统计图，图3是经本发明的岩性识别方法输出的岩性统计图。从中可以看出，由于不同岩性的波阻抗分布区域存在差别，图2的砂岩与泥岩波阻抗分布区域存在较大重叠区，其严重影响岩性边界识别的精度；而本发明创造性的引入了岩性识别敏感因子s，其继承了波阻抗的趋势特征，又在细节上加入了高频敏感曲线，由此可更加精准地识别岩性边界，从图3中可以直观看出，砂岩与泥岩波阻抗分布较为独立，重叠区域几乎没有，由此可极大提高波阻抗反演分辨岩性边界的能力，修正岩性识别不清的缺陷。
实施例二：
本实施例还提供一种实现上述岩性识别方法的岩性识别系统，如图4所述，其包括：
第一曲线获取单元1，其用于获取岩性曲线，以及根据所述测井资料计算获得波阻抗曲线a；
第二曲线获取单元2，其用于根据所述岩性曲线获得至少一个岩性的岩性敏感曲线b；
第一曲线处理单元3，其用于对所述波阻抗曲线a进行低通滤波处理，计算得到低频趋势波阻抗alpf；
第二曲线处理单元4，其用于对所述敏感曲线b进行高通滤波处理，计算得到高频敏感曲线bhpf；
敏感因子获取单元5，其用于根据所述低频趋势波阻抗alpf与高频敏感曲线bhpf计算得到岩性识别敏感因子；
以及数据输出单元6，其用于根据所述岩性识别敏感因子输出岩性统计数据。
需要说明的是，上述实施例一、二中的技术特征可进行任意组合，且组合而成的技术方案均属于本发明的保护范围。
与现有技术相比，本发明的有益效果是：本发明创造性的引入了岩性识别敏感因子s，其继承了波阻抗的趋势特征，又在细节上加入了高频敏感曲线，由此可更加精准地识别岩性边界，从图3中可以直观看出，砂岩与泥岩波阻抗分布较为独立，重叠区域几乎没有，由此可极大提高波阻抗反演分辨岩性边界的能力，修正岩性识别不清的缺陷，最大程度提高地震储层预测能力。
需要说明的是，上述实施例一、二的技术特征可进行任意组合，且组合而成的技术方案均属于本发明的保护范围。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。