[发明专利]基于地质统计学的叠前振幅随偏移距变化油气预测方法

说明书

本发明提供了一种基于地质统计学分析的叠前振幅随偏移距变化油气预测方法。所述油气预测方法包括：根据测井数据建立油气指示曲线，在油气指示曲线上，油气段具有第一油气指示值，非油气段具有第二油气指示值；计算得到AVO属性数据体；分析得到油气指示曲线与井旁道AVO属性数据相关性最好的最优AVO属性数据体；将油气指示曲线中非油气段部分和油气段部分分别与最优AVO属性数据体进行直方图分析，计算得到油气指示曲线中非油气段部分和油气段部分的隶属度函数关系；对油气指示曲线进行地质统计学分析得到变差函数；利用隶属度函数关系和变差函数进行油气模拟分析得到油气预测结果。根据本发明的方法能够提高流体预测的精度及分辨率。

技术领域

本发明属于地球物理勘探技术中的流体预测方法，更具体地讲，涉及一种能够高精度、高分辨地预测薄储层流体的地质统计学AVO分析方法。

背景技术

AVO(Amplitude variation offset，振幅随偏移距变化)技术是一种研究地震反射振幅随炮检距变化的技术，振幅随偏移距的变化是地下岩石及孔隙流体弹性参数的函数，因此根据AVO变化规律所反应出来的地下岩性及其孔隙流体性质可以用来直接预测油气和估算储层岩性参数。传统AVO分析就是采用Zoeppritz方程或其近似方程由实际地震道集记录反演估算岩石的地震参数(如密度、纵横波速度)的一种反演方法。常规AVO分析方法通常是把非线性的叠前反演问题线性化，解的稳定性受初始模型的影响较大，且容易陷入局部极小，而AVO分析反演方法所得到的结果分辨率较差不利于储层流体的预测。

发明内容

针对现有技术中存在的问题，本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如，本发明的目的在于解决传统AVO分析分辨率较低、精度不高的问题。

为了实现上述目的，本发明提供了一种基于地质统计学分析的叠前振幅随偏移距变化油气预测方法。所述油气预测方法包括：A、根据测井数据建立油气指示曲线，在所述油气指示曲线上，油气段具有第一油气指示值，非油气段具有第二油气指示值；B、计算得到AVO属性数据体；C、分析得到油气指示曲线与井旁道AVO属性数据相关性最好的最优AVO属性数据体；D、将油气指示曲线中非油气段部分和油气段部分分别与所述最优AVO属性数据体进行直方图分析，计算得到油气指示曲线中非油气段部分和油气段部分的隶属度函数关系；E、对油气指示曲线进行地质统计学分析得到变差函数；F、利用所述隶属度函数关系和所述变差函数进行油气模拟分析，得到油气预测结果。在本发明的一个示例性实施例中，所述振幅随偏移距变化属性体包括L1(一范数)、L2(二范数)、P+G、P-G、P*G、L1P+G、L1P-G、L1P*G、L2P+G、L2P-G、L2P*G，其中，P代表截距，G代表梯度。

在本发明的一个示例性实施例中，所述步骤F可以包括：F01、在振幅随偏移距变化属性体上创建数据网格；F02、查找数据网格上的网格节点，并判断网格节点是否需要模拟，其中，需要模拟的网格节点为不确定油气指示值的未知网格节点，不需要进行模拟的网格节点为已确定油气指示值的已知网格节点；F03、查找需要模拟的未知网格节点领域内的已知网格节点，直到已知网格节点达到预定数量；F04、通过克里金插值算法，利用所述预定数量的已知网格节点计算未知网格节点的克里金值，然后利用克里金值和所述变差函数计算得到克里金概率值；F05、利用未知网格节点的克里金概率值和对应的振幅随偏移距变化属性数据体的值，联合比对所述隶属度函数关系，确定未知网格节点的油气指示值，并将确定了油气指示值的未知网格节点作为已知网格节点；F06、重复步骤F02至F05，直至确定出所有的未知网格节点的油气指示值。

在本发明的一个示例性实施例中，在所述步骤F02中，可以按照随机路径序贯地游走数据网格查找网格节点；

在本发明的一个示例性实施例中，在所述步骤F03中，可以以螺旋搜索的方式进行查找。

在本发明的一个示例性实施例中，在所述步骤F04中，采用克里金插值算法中的克里金方程计算克里金权值系数，计算未知网格节点的克里金概率值。