[发明专利]一种微观储层岩石电导率预测方法

权利要求书

1.一种微观储层岩石电导率预测方法，其特征在于，包括：

利用有限元仿真软件，计算三维数字岩心模型M的电导率；

对用于描述三维数字岩心模型M的二维闵可夫斯基函数和电导率进行典型相关分析，得到典型协变量XC和YC；并得到三维数字岩心模型M的二维闵可夫斯基函数到XC的变换方式；

用核密度估计法估计XC和YC的联合概率密度分布情况，得到先验概率分布模型；

利用二维闵可夫斯基函数描述待预测二维岩心样本图像的几何结构特征；

对待预测二维岩心样本图像的二维闵可夫斯基函数进行典型相关分析，按照三维数字岩心模型M的二维闵可夫斯基函数到XC的变换方式将待预测二维岩心样本图像的二维闵可夫斯基函数变换为典型协变量X；

将典型协变量X输入到训练好的先验概率分布模型中，得到后验概率分布模型；其中，后验概率分布模型中概率最高点对应的电导率即为待预测二维岩心样本的电导率可能性最大的值。

2.根据权利要求1所述的一种微观储层岩石电导率预测方法，其特征在于，所述三维数字岩心模型M的生成方法包括：

将二维岩石铸体薄片图像作为多点统计算法的训练图像，利用多点统计算法生成多个三维数字岩心模型I；

将多个三维数字岩心模型I作为生成对抗网络的训练图像，生成多个三维数字岩心模型M。

3.根据权利要求1所述的一种微观储层岩石电导率预测方法，其特征在于，还包括：

对待预测二维岩心样本的二维闵可夫斯基函数进行典型相关分析之前，先利用单分类支持向量机模型，对待预测岩心样本的二维闵可夫斯基函数进行异常值检测，筛选掉不能用于预测的异常样本。

4.根据权利要求1所述的一种微观储层岩石电导率预测方法，其特征在于，所述二维闵可夫斯基函数包括三个闵可夫斯基参数，三个所述闵可夫斯基参数分别是孔隙度、周长面积比和欧拉数。

5.根据权利要求4所述的一种微观储层岩石电导率预测方法，其特征在于，所述孔隙度的计算公式为：

其中，φ为孔隙度，SP为孔隙的面积总和，SR为二维岩心样本总面积；

所述周长面积比的计算公式为：

其中，W为周长面积比，CP为样本孔隙总周长，SP为孔隙的面积总和；

所述欧拉数的计算公式为：

E＝NP-HP

其中，E为欧拉数，NP为孔隙连接体数，HP为孔洞数。

6.根据权利要求1所述的一种微观储层岩石电导率预测方法，其特征在于，所述三维数字岩心模型的二维闵可夫斯基函数的计算方式包括：

将三维数字岩心模型分为多层切片，计算每层切片的二维闵可夫斯基函数后取均值。

7.根据权利要求1所述的一种微观储层岩石电导率预测方法，其特征在于，还包括：

对后验概率分布结果设置置信区间以衡量预测电导率的合理性。

8.根据权利要求7所述的一种微观储层岩石电导率预测方法，其特征在于，所述置信区间的计算公式如下所示：

其中，CI指代置信区间，为样本的均值，z为置信水平值，n为样本数量，z值取0.95，S为样本的标准差。