[发明专利]基于微观赋存形态预测水合物饱和度的方法及装置

权利要求书

1.一种基于微观赋存形态预测水合物饱和度的方法，其特征在于，包括：

根据电阻率或核磁共振测井解释的水合物饱和度和实际测量的纵横波速度，确定水合物多种微观赋存形态的比例；

以电阻率或核磁共振测井解释的水合物饱和度作为待预测水合物饱和度的初始值，并利用确定的水合物多种微观赋存形态的比例来分配不同赋存形态的水合物饱和度，循环确定每一周期对应的纵横波速度误差值，每个循环周期均执行以下操作：

根据当前周期的水合物饱和度值和确定的多种赋存形态百分比，以及预先建立的统一接触胶结模型，得到当前周期计算的纵横波速度；所述统一接触胶结模型是根据引入的正则化接触胶结半径来改进胶结模型并采用Hashin-Shtrikman边界条件约束而建立的；

确定当前周期计算的纵横波速度与实际测量的纵横波速度之间的误差，得到当前周期对应的纵横波速度误差值；在不满足预设循环结束条件时，给定下一周期的水合物饱和度值；

在满足预设循环结束条件时，比对各个周期对应的纵横波速度误差值，将纵横波速度误差值最小时对应的水合物饱和度值作为预测的水合物饱和度值。

2.如权利要求1所述的基于微观赋存形态预测水合物饱和度的方法，其特征在于，根据当前周期的水合物饱和度值和确定的多种赋存形态百分比，以及预先建立的统一接触胶结模型，得到当前周期计算的纵横波速度，包括：

根据当前周期的水合物饱和度值和确定的多种赋存形态百分比，得到岩石基质模量，流体模量以及基于预先建立的统一接触胶结模型得到干岩石骨架模量；

根据岩石基质模量，流体模量以及干岩石骨架模量，得到当前周期的饱和岩石模量；

根据当前周期的饱和岩石模量，得到当前周期计算的纵横波速度。

3.如权利要求1所述的基于微观赋存形态预测水合物饱和度的方法，其特征在于，在待测储层的孔隙度(φ)低于临界孔隙度(φ_c)时，所述统一接触胶结模型包括：

在待测储层的孔隙度(φ)高于临界孔隙度(φ_c)时，所述统一接触胶结模型包括：

其中，K_dry为有效干岩石体积模量，μ_dry为有效干岩石体积剪切模量，K_ma骨架为基质的体积模量，K_ect和μ_ect为临界孔隙度情况下的岩石体积模量和剪切模量，φ为待测储层的孔隙度，φ_c为临界孔隙度。

4.如权利要求1所述的基于微观赋存形态预测水合物饱和度的方法，其特征在于，根据电阻率或核磁共振测井解释的水合物饱和度和实际测量的纵横波速度，确定水合物多种微观赋存形态的比例，包括：

根据预先建立的岩石物理模型，构建纵横波速度与水合物饱和度的交会图；所述岩石物理模型根据统一接触胶结模型得到；

根据所述交会图，确定水合物多种微观赋存形态的比例。

5.如权利要求4所述的基于微观赋存形态预测水合物饱和度的方法，其特征在于，所述岩石物理模型为：

K_sat＝K_dry+λ²K_fl[φ_e+(λ-φ_e)K_fl/K_ma]^-1；

μ_sat＝μ_dry；

其中，式中Biot系数λ＝1-K_dry/K_ma，有效孔隙度φ_e＝φ(1-S_hms-S_hc)，K_fl为流体的体积模量，S_hc，S_hpf和S_hms为胶结型、孔隙充填型和骨架支撑型水合物的饱和度，K_h和K_w分别为水合物和水的体积模量。