[发明专利]页岩储层含气性分析方法及装置

权利要求书

1.一种页岩储层含气性分析方法，其特征在于，包括：

依据横波测井曲线计算得到页岩储层随深度变化的多个弹性参数曲线，所述多个弹性参数曲线包括纵波阻抗曲线、横波阻抗曲线、杨式模量曲线、泊松比曲线和纵横波速度比曲线；

依据积差方法求取通过岩心测试岩心实测散点数据与弹性参数曲线上对应弹性参数之间的关联系数，所述岩心实测散点数据包括有机含碳量、孔隙度和含气量；

选取与岩心实测散点数据之间的关联系数大于预设阈值的弹性参数作为敏感弹性参数；

以敏感弹性参数为自变量、与所述敏感弹性参数对应的岩心实测散点数据为因变量，通过最小二乘法确定与所述敏感弹性参数对应的岩心实测散点数据最佳匹配的多元回归方程；

基于横波测井曲线和三维叠前道集数据通过叠前同时反演技术得到弹性参数体数据，所述弹性参数数据体包括纵波阻抗三维数据体、横波阻抗三维数据体和密度三维数据体；

依据所述弹性参数数据体得到对应的敏感弹性参数数据体；

依据所述敏感弹性参数数据体和所述多元回归方程，得到页岩储层的有机含碳量数据体、孔隙度数据体和含气量数据体；

利用测得的泥质含量对页岩储层的声波速度曲线进行过滤，并依据页岩储层随深度变化的声波速度趋势，建立单井正常压实地层的声波速度趋势线；

依据所述声波速度趋势线得到单井地层孔隙压力，并依据实际测得的地层孔隙压力数据得到区域伊顿指数；

依据所述声波速度趋势线和所述区域伊顿指数，得到地层孔隙压力数据体；

依据地层孔隙压力数据体和静水压力数据体得到地层压力系数体；

采用极值法对有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行归一化处理；

对归一化处理后的有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行加权运算，得到页岩储层含气性综合评价指数。

2.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，所述岩心实测散点数据与弹性参数曲线之间的关联系数为其中，R为关联系数，n为进行关联分析的样本个数，x_i为有机含碳量、孔隙度和含气量中的其中一个，yi为弹性参数曲线中的对应样本值。

3.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，所述多元回归方程为Y＝β₀+β₁χ₁+β₂χ₂+...+β_pχ_p+ε，其中，Y为有机含碳量、孔隙度和含气量中的其中一个，χ₁、χ₂、χ_p为敏感弹性参数，β₀、β₁、β₂、β_p为建立回归方程的回归系数，ε为随机误差项，当ε为0时得到最佳匹配的多元回归方程。

4.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，所述声波速度趋势线为其中，V_P为泥岩地层的声波速度，V_PMatrix为页岩储层骨架的声波速度，V_PTop为深度为0时声波的初始速度，b为页岩储层实压系数，Zml为垂直深度。

5.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，所述单井地层孔隙压力为P_p＝P_o-(P_o-P_w)×(ν_inst/ν_n)^N，其中，P_p为预测的页岩储层孔隙压力，P_o为静岩压力，P_w为静水压力，ν_n为正常压实页岩储层的声波速度，ν_inst为页岩储层实测的声波速度，N为幂指数。

6.根据权利要求5所述的页岩储层含气性分析方法，其特征在于，所述地层压力系数体为