基于地质统计学的叠前振幅随偏移距变化油气预测方法技术

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

Oil and gas prediction method based on Geostatistics for pre stack amplitude variation with offset

The invention provides an oil and gas prediction method based on geostatistical analysis of pre stack amplitude and offset variation. The oil and gas prediction method comprises: according to logging data established for indicating oil and gas in the oil curve, indicating curve of gas, oil and gas section is provided with a first oil gas indication, non oil and gas segment has second oil and gas indicator; the calculated AVO attribute data; analysis of oil and gas are indicative curve and near wells AVO attribute data the best correlation optimal AVO attribute data; the indicating oil and gas oil and gas section of oil and gas are part of the histogram analysis and optimal AVO attribute data curve, non oil and gas and oil and gas section part of the membership function between the calculated curve indicating the oil and gas; oil and gas index curve geological statistical analysis obtained the variogram function; the membership functions and the relationship of oil and gas prediction results of simulation analysis of oil and gas are difference function. The method can improve the accuracy and resolution of the fluid prediction according to the method of the present invention.

【技术实现步骤摘要】
基于地质统计学的叠前振幅随偏移距变化油气预测方法
本专利技术属于地球物理勘探技术中的流体预测方法，更具体地讲，涉及一种能够高精度、高分辨地预测薄储层流体的地质统计学AVO分析方法。
技术介绍
AVO(Amplitudevariationoffset，振幅随偏移距变化)技术是一种研究地震反射振幅随炮检距变化的技术，振幅随偏移距的变化是地下岩石及孔隙流体弹性参数的函数，因此根据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中，采用克里金插值算法中的克里金方程计算克里金权值系数，计算未知网格节点的克里金概率值。与现有技术相比，本专利技术的有益效果包括：提高流体预测的精度及分辨率。附图说明图1示出了根据本专利技术的示例性实施例的基于地质统计学的叠前AVO油气预测方法油气模拟步骤的流程图。图2示出了根据本专利技术的示例性实施例的基于地质统计学的叠前AVO油气预测方法点模拟的流程图。具体实施方式在下文中，将结合示例性实施例和附图来详细说明本专利技术的基于地质统计学的叠前AVO油气预测方法。根据本专利技术示例性实施例的基于地质统计学的叠前AVO油气预测方法可以包括如下步骤：A、根据测井数据建立油气指示曲线，在所述油气指示曲线上，油气段具有第一油气指示值(例如1)，非油气段具有第二油气指示值(例如0)。B、计算得到AVO属性数据体。AVO属性数据体包括L1(一范数)、L2(二范数)、P+G(截距加梯度)、P-G(截距减梯度)、P*G(截距乘梯度)、L1P+G、L1P-G、L1P*G、L2P+G、L2P-G、L2P*G等属性体。其中，P代表截距，G代表梯度。C、分析得到油气指示曲线与井旁道AVO属性数据相关性最好的最优AVO属性数据体。D、将油气指示曲线中非油气段部分和油气段部分分别与所述最优AVO属性数据体进行直方图分析，计算得到油气指示曲线中非油气段部分和油气段部分的隶属度函数关系。E、对油气指示曲线进行地质统计学分析得到变差函数。F、利用所述隶属度函数关系和所述变差函数进行油气模拟分析，得到油气预测结果。所述步骤F具体包括：F01、在振幅随偏移距变化属性体上创建数据网格。F02、按照随机路径序贯地游走数据网格查找网格节点，并判断网格节点是否需要模拟，其中，需要模拟的网格节点为不确定油气指示值的未知网格节点，不需要进行模拟的网格节点为已确定油气指示值的已知网格节点。F03、以螺旋搜索的方式查找需要模拟的未知网格节点领域内的已知网格节点，直到已知网格节点达到预定数量。F04、通过克里金插值算法，利用所述预定数量的已知网格节点计算未知网格节点的克里金值，然后利用克里金值和所述变差函数计算得到克里金概率值。F05、利用未知网格节点的克里金概率值和对应的振幅随偏移距变化属性数据体的值，联合比对所述隶属度函数关系，确定未知网格节点的油气指示值，并将确定了油气指示值的未知网格节点作为已知网格节点。F06、重复步骤F02至F05，直至确定出所有的未知网格节点的油气指示值。本专利技术的主要技术思路是将AVO计算的属性体，通过一系列计算，转换为0和1构成的数据体，1指示油气，0指示非油气，实现油气检测。在本专利技术的一个示例性实施例中，基于地质统计学的叠前AVO油气预测方法可以包括如下步骤：(1)首先通过测井数据分析，建立一条HCI指示曲线，HCI曲线是一条油气段HCI值为1、非油气段值为0的指示值曲线。但本专利技术不限于此，油气段和非油气的指示值可以为其它任何能够区分二者不同的值。这里，HCI是HydrocarbonIndicator的缩写，HCI曲线就是油气指示曲线，准确性是由测井曲线控制，HCI是一条由0和1组成的指示油气和水的曲线。(2)通过计算传统AVO分析的结果，得到一系列相关的AVO属性体，比如L1、L2、P+G、P-G、P*G、L1P+G、L1P-G、L1P*G、L2P+G、L2P-G和L2P*G等属性体，其中，P为截距，G为梯度，L1为一范数，L2为二范数。(3)分析HCI曲线与井旁AVO属性道数据之间的相关性，分析并得到HCI-AVO属性相关性最高的HCI-AVO属性曲线，并对其进行保存。(4)对保存的HCI-AVO属性曲线对HCI-0和HCI-1曲线进行直方图分析，并计算得到HCI-0及HCI-1的隶属度函数关系。这里，HCI-AVO属性曲线是通过HCI曲线与井旁道AVO属性数据相关性最好的AVO属性，并将其0和1所对应的AVO属性数据分别保存的到HCI本文档来自技高网...

【技术保护点】
一种基于地质统计学的叠前振幅随偏移距变化油气预测方法，其特征在于，所述油气预测方法包括以下步骤：A、根据测井数据建立油气指示曲线，在所述油气指示曲线上，油气段具有第一油气指示值，非油气段具有第二油气指示值；B、计算得到振幅随偏移距变化属性数据体；C、分析得到油气指示曲线与井旁道振幅随偏移距变化属性数据相关性最好的最优振幅随偏移距变化属性数据体；D、将油气指示曲线中非油气段部分和油气段部分分别与所述最优振幅随偏移距变化属性数据体进行直方图分析，计算得到油气指示曲线中非油气段部分和油气段部分的隶属度函数关系；E、对油气指示曲线进行地质统计学分析得到变差函数；F、利用所述隶属度函数关系和所述变差函数进行油气模拟分析，得到油气预测结果。

【技术特征摘要】
1.一种基于地质统计学的叠前振幅随偏移距变化油气预测方法，其特征在于，所述油气预测方法包括以下步骤：A、根据测井数据建立油气指示曲线，在所述油气指示曲线上，油气段具有第一油气指示值，非油气段具有第二油气指示值；B、计算得到振幅随偏移距变化属性数据体；C、分析得到油气指示曲线与井旁道振幅随偏移距变化属性数据相关性最好的最优振幅随偏移距变化属性数据体；D、将油气指示曲线中非油气段部分和油气段部分分别与所述最优振幅随偏移距变化属性数据体进行直方图分析，计算得到油气指示曲线中非油气段部分和油气段部分的隶属度函数关系；E、对油气指示曲线进行地质统计学分析得到变差函数；F、利用所述隶属度函数关系和所述变差函数进行油气模拟分析，得到油气预测结果。2.根据权利要求1所述的基于地质统计学的叠前振幅随偏移距变化油气预测方法，其特征在于，所述振幅随偏移距变化属性体包括L1、L2、P+G、P-G、P*G、L1P+G、L1P-G、L1P*G、L2P+G、L2P-G和L2P*G，其中，P为截距，G为梯度，L1为一范数，L2为二范数。3.根据权利要求1所述的基于地质统计学的叠前振幅随偏移距变化油气预测方法，其特征在于，所述步骤F包括：F01、在振幅随偏移距变化属性体上创建数据网格；F02、查找数据网格上的网格节点，...

【专利技术属性】
技术研发人员：刘开元，范晓，陈小二，邹文，吴秋波，康昆，巫骏，
申请(专利权)人：中国石油集团川庆钻探工程有限公司地球物理勘探公司，
类型：发明
国别省市：四川,51