[发明专利]基于三维地震数据的砂泥岩孔隙类型反演方法

说明书

技术领域

本发明专利涉及测量技术，特别涉及基于三维地震数据的砂泥岩孔隙类型反演方法。

背景技术

在现有的技术中，只有基于测井资料的孔隙扁率反演方法，主要是针对单井的岩石等效孔隙纵横比进行反演，该方法只能获得井点处的孔隙结构信息，不能有效地获得反映研究区孔隙结构信息的三维数据。

发明内容

本发明的目的在于提供基于三维地震数据的砂泥岩孔隙类型反演方法，以有效地获得反映研究区孔隙结构信息的三维数据。

本发明提供一种基于三维地震数据的砂泥岩孔隙类型反演方法，包括：设定当前采样点的线号、道号和时间采样点号，提取当前采样点的孔隙度φ'、基质矿物的等效体积模量K′₀、基质矿物的等效剪切模量μ′₀、砂岩的百分比含量v'_s、泥岩的百分比含量或泥质含量v'_c、密度ρ'和纵波速度V′_s；设定当前采样点的砂岩孔隙扁率α′_s的初始值公式、搜索范围和搜索间隔；并设定每个采样点的泥岩孔隙的扁率α′_c的初始值、搜索范围和搜索间隔；设定完当前采样点的砂岩孔隙扁率α'_s和泥岩孔隙扁率α'_c的初始值后，计算孔隙度为φ'时的干岩石剪切模量μ′_dry；计算当前采样点的横波速度V′_{s_0}；进行当前采样点的计算的两种横波速度V′_{s_0}和V′_s的误差分析，并保存当前砂岩孔隙的扁率α'_s和泥岩孔隙的扁率α'_c的值所对应的误差值σ。

优选的，还包括：修改当前采样点的砂岩孔隙扁率α′_s和泥岩孔隙扁率α′_c的值，直到搜索范围内所有的扁率值计算完毕；寻找砂岩孔隙的扁率α′_s和泥岩孔隙的扁率α′_c的值所对应的最小误差值σ，并输出对应的砂岩孔隙的扁率α′_s、泥岩孔隙的扁率α'_c和φ_s'、φ_c'作为当前采样点的计算结果。

优选的，还包括：修改当前采样点的线号、道号和时间采样点号，直到所有的采样点计算完毕，并输出砂岩孔隙的扁率α_s及其百分比含量φ_s、泥岩孔隙的扁率α_c及其百分比含量φ_c的三维数据体；对α_s、φ_s、α_c、φ_c进行三维中值滤波，去除异常值的影响，并作为最终的结果输出。

优选的，还包括：输入砂泥岩研究区的测井参数：纵波速度V_p、横波速度V_s、孔隙度φ和密度ρ；统计得到V_p与ρ的关系式，拟合得到系数α和β；再利用双波法统计φ与V_p和V_s的关系式，拟合得到系数A、B、C。

本发明的有益效果在于：基于砂泥岩研究区的三维地震资料和测井资料，计算反映砂泥岩研究区的孔隙类型的三维数据体，获得地下每个采样点的孔隙类型信息(包括组成孔隙的各种孔隙类型的扁率及其百分比含量)，从而有效地描述砂泥岩储层的孔隙类型，为砂泥岩储层预测和钻探提供重要的参考数据。

附图说明

图1是基于三维地震数据的砂泥岩孔隙类型反演方法的流程图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下文中将结合实施例进行详细说明。

如图1所示，本发明基于三维地震数据的砂泥岩孔隙类型反演方法具体包括如下步骤：

步骤(1)，输入砂泥岩研究区的测井参数：纵波速度(V_p)、横波速度(V_s)、孔隙度(φ)和密度(ρ)。利用Gardner公式分别统计得到V_p与ρ的关系式(见公式4-1)，拟合得到系数α和β；再利用双波法统计φ与V_p和V_s的关系式(见公式4-2)，拟合得到系数A、B、C；

φ＝AV_p+BV_s+C (4-2)