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

权利要求书

1.一种基于三维地震数据的砂泥岩孔隙类型反演方法，其特征在于，包括：

设定当前采样点的线号、道号和时间采样点号，提取当前采样点的孔隙度φ′、基质矿物的等效体积模量K′₀、基质矿物的等效剪切模量μ′₀、砂岩的百分比含量v′_s、泥岩的百分比含量或泥质含量v′_c、密度ρ′和横波速度V′_s；

设定当前采样点的砂岩孔隙扁率α′_s的初始值公式、搜索范围和搜索间隔；并设定每个采样点的泥岩孔隙的扁率α′_c的初始值、搜索范围和搜索间隔；

设定完当前采样点的砂岩孔隙扁率α′_s和泥岩孔隙扁率α′_c的初始值后，计算孔隙度为φ′时的干岩石剪切模量μ′_dry；

计算当前采样点的横波速度V′_{s_0}；

进行当前采样点的计算的两种横波速度V′_{s_0}和V′_s的误差分析，并保存当前砂岩孔隙的扁率α′_s和泥岩孔隙的扁率α′_c的值及其所对应的误差值σ。

2.如权利要求1所述的基于三维地震数据的砂泥岩孔隙类型反演方法，其特征在于，还包括：输入砂泥岩研究区的测井参数：纵波速度V_p、横波速度V_s、孔隙度φ和密度ρ；统计得到V_p与ρ的关系式，拟合得到系数α和β；再利用双波法统计φ与V_p和V_s的关系式，拟合得到系数A、B、C。

3.如权利要求2所述的基于三维地震数据的砂泥岩孔隙类型反演方法，其特征在于，还包括：输入砂泥岩研究区的基于叠后反演得到的伽玛GR三维地震数据，并采计算泥质含量即泥岩的百分比含量v_c的三维地震数据。

4.如权利要求3所述的基于三维地震数据的砂泥岩孔隙类型反演方法，其特征在于，还包括：输入基于叠前弹性反演得到纵波阻抗I_P、横波阻抗I_s的三维地震数据；计算V_p、V_s和ρ的三维地震数据；计算孔隙度φ的三维地震数据。