[发明专利]一种碳酸盐岩储层孔隙度计算方法及装置

权利要求书

1.一种碳酸盐岩储层孔隙度计算方法，其特征在于，包括以下步骤：

开展岩心孔隙度实验数据标定，交会拟合建立孔隙度大小与测井值大小之间的关系，确定碳酸盐岩储层孔隙度测井响应的非线性特征，对于测井响应非线性特征，找出适合于描述该测井曲线-孔隙度散点分布规律的非线性函数关系；

建立碳酸盐岩储层体积模型和测井响应方程，其中，对于非线性特征，基于非线性函数关系，建立非线性方程；

上述的测井响应方程联立形成方程组；

迭代求解方程组，得到储层孔隙度。

2.根据权利要求1所述的碳酸盐岩储层孔隙度计算方法，其特征在于，在所述步骤迭代求解方程组之后，还包括以下步骤：

误差分析，在误差值满足预设要求的情况下，则计算结果中的储层孔隙度满足误差要求；在误差值不满足预设要求的情况下，继续迭代求解方程组，直到计算结果中的储层孔隙度满足误差要求。

3.根据权利要求2所述的碳酸盐岩储层孔隙度计算方法，其特征在于，在所述误差分析中：

将求解结果带入到方程组中，得到正演测井响应值，令测井得到的测井响应值与正演测井响应值之间的差值为求解误差；

所述预设要求为：

根据满足实际要求的误差值标准，当求解误差的绝对值小于误差值标准时，满足预设要求。

4.根据权利要求2所述的碳酸盐岩储层孔隙度计算方法，其特征在于，在所述迭代求解方程组中：

利用牛顿法对方程组进行数值求解。

5.根据权利要求1所述的碳酸盐岩储层孔隙度计算方法，其特征在于，在所述开展岩心孔隙度实验数据标定中：

将测井曲线和已有的岩心孔隙度置于同一坐标系中，调整岩心孔隙度数据的位置，使得岩心孔隙度随地层深度的变化趋势与测井曲线的变化趋势一致或基本一致。

6.根据权利要求5所述的碳酸盐岩储层孔隙度计算方法，其特征在于，在所述交会拟合建立孔隙度大小与测井值大小之间的关系，确定碳酸盐岩储层孔隙度测井响应的非线性特征中：

取同一深度下的孔隙度值和测井曲线值，以孔隙度值为横坐标，测井曲线值为纵坐标，在坐标系中标注出不同深度下的孔隙度值和测井曲线值所对应的散点，若散点的分布规律呈现非线性特点，则该测井曲线即为碳酸盐岩储层孔隙度测井响应的非线性特征。

7.根据权利要求1-6中任意一项所述的碳酸盐岩储层孔隙度计算方法，其特征在于，所述测井响应方程包括自然伽玛测井响应方程、中子测井响应方程和声波测井响应方程。

8.根据权利要求7所述的碳酸盐岩储层孔隙度计算方法，其特征在于，所述自然伽玛测井响应方程为：

所述中子测井响应方程为：

所述声波测井响应方程为：

其中，GR、CNL、AC分别为测井得到的自然伽玛、中值、声波测井值；

V₁、V₂、V₃、分别为泥质、灰岩、白云岩、孔隙体积百分比；

GR_clay、GR_lime、GR_dom、GR_fluid分别为泥质、灰岩、白云岩、孔隙流体自然伽玛值，为已知参数；

CNL_clay、CNL_lime、CNL_dom、CNLf_luid分别为泥质、灰岩、白云岩、孔隙流体中子值，为已知参数；

AC_clay、AC_lime、AC_dom、ACf_luid分别为泥质、灰岩、白云岩、孔隙流体声波值，为已知参数；

f为岩心标定实验得到的孔隙度与测井响应曲线之间的函数关系。