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

说明书

本发明公开了一种基于微观赋存形态预测水合物饱和度的方法及装置，方法包括：根据电阻率或核磁共振测井解释的水合物饱和度和实际测量的纵横波速度，确定水合物多种微观赋存形态的比例；以测井解释的水合物饱和度作为初始值，结合该比例，循环确定每一周期的纵横波速度误差值：根据当前周期的水合物饱和度值，以及预先建立的统一接触胶结模型，得到当前周期计算的纵横波速度；确定当前周期计算与实际测量的纵横波速度的误差；在不满足循环结束条件时，给定下一周期的水合物饱和度值；在满足条件时，比对各个周期的纵横波速度误差值，得到预测的水合物饱和度值。本发明可以提高预测天然气水合物饱和度的精度，对石油天然气开采具有科学的指导意义。

技术领域

本发明涉及石油地球物理勘探技术领域，尤其涉及基于微观赋存形态预测水合物饱和度的方法及装置。

背景技术

本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。

近年来，天然气水合物因被视为潜在的替代资源和对气候和环境的影响而越来越受到广泛关注。因此，从地球物理测量数据中准确地预测天然气水合物的饱和度十分重要。然而，目前预测天然气水合物饱和度的方案精度低。

发明内容

本发明实施例提供一种基于微观赋存形态预测水合物饱和度的方法，用以提高预测水合物饱和度的精度，该方法包括：

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

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

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

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

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

本发明实施例还提供一种基于微观赋存形态预测水合物饱和度的装置，用以提高预测水合物饱和度的精度，该装置包括：

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

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

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

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

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