[发明专利]一种井中天然气水合物饱和度岩石物理计算方法及系统

权利要求书

1.一种井中天然气水合物饱和度岩石物理计算方法，其特征在于，该方法包括：

获取工区资料；

根据所述工区资料，分别对当前深度点的含水合物层位及不含水合物层位进行井中地层孔隙度计算，得到井中地层孔隙度；

根据所述工区资料及井中地层孔隙度，经过反演得到当前深度点的地层最佳固结系数；

对井中每个深度点进行井中地层孔隙度及固结系数的计算，得到整个井段的地层孔隙度及固结系数；

根据地层孔隙度及固结系数，利用三相Biot模型预测地层纵波速度及地层密度，并以修正系数为约束计算井中每个深度点的水合物饱和度；

通过修改权重因子调整修正系数，利用新的修正系数为约束计算井中不含水合物层位的饱和度平均值，直至所述不含水合物层位的饱和度平均值趋近于0时，得到本井的权重因子；

利用所述本井的权重因子，重新计算井中每个深度点的水合物饱和度，得到整个井段的水合物饱和度值。

2.根据权利要求1所述的井中天然气水合物饱和度岩石物理计算方法，其特征在于，获取工区资料，包括：

通过工区测井和岩心分析，获得测井数据及岩石矿物数据；其中，

所述测井数据包括：研究区的深度、纵波速度、电阻率、密度、自然伽马及声波速度数据；

所述岩石矿物数据包括：井中地层矿物颗粒的体积模量、剪切模量、体积百分比、密度及中子响应值；

根据岩石矿物数据，得到固体成分的体积模量、剪切模量、密度及中子响应值，计算公式如下：

ρ_ma＝∑f_iρ_i；

其中，K_ma、μ_ma、ρ_ma、分别为基质的体积模量、剪切模量、密度、中子响应值；K_i、μ_i、ρ_i、f_i分别为每种矿物颗粒的体积模量、剪切模量、密度、中子响应值、体积百分比。

3.根据权利要求2所述的井中天然气水合物饱和度岩石物理计算方法，其特征在于，该方法还包括：

根据所述工区资料，利用含水合物层在测井曲线上表现为自然伽马值低于第一阈值、电阻率值高于第二阈值及声波速度值高于第三阈值的特征，划分得到含水合物层，并将其它层位划分为不含水合物层。