(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211122313.5 (22)申请日 2022.09.15 (71)申请人 长安大学 地址 710064 陕西省西安市南 二环路中段 长安大学 (72)发明人 左琛　郭晨　闫正林　 (74)专利代理 机构 西安铭泽知识产权代理事务 所(普通合伙) 61223 专利代理师 田甜 (51)Int.Cl. G06F 30/23(2020.01) G06F 30/27(2020.01) G01R 27/02(2006.01) G06F 111/08(2020.01) (54)发明名称 一种微观储层岩石电导 率预测方法 (57)摘要 本发明提供了一种微观储层岩石电导率预 测方法， 属于油气储集层导电机理研究领域， 包 括利用二维闵可夫斯基函数描述三维数字岩心 模型M的几何结构特征； 对其二维闵可夫斯基函 数和电导率进行典型相关分析， 得到典型协变量 XC和YC； 用核密度法估计 XC和YC的联合概率密度 分布情况， 得到先验概率分布模型； 用二维闵可 夫斯基函数描述待预测二维岩心样本的几何结 构特征； 对其二维闵可夫斯基函数进行典型相关 分析， 得到典型协变量X； 将典型协变量X输入到 训练好的先验概率分布模型中， 得到后验概率分 布模型； 根据后验概率分布模型预测其电导率。 该方法能够对微观储层岩石的结构和其电导率 之间的关系进行 预测。 权利要求书2页 说明书5页 附图5页 CN 115455772 A 2022.12.09 CN 115455772 A 1.一种微观储层岩石电导 率预测方法， 其特 征在于， 包括： 利用有限元仿真软件， 计算 三维数字岩心模型M的电导 率； 对用于描述三维数字岩心模型M的二维闵可夫斯基函数和电导率进行典型相关分析， 得到典型协变量XC和Y C； 用核密度估计法估计XC和Y C的联合概率密度分布情况， 得到先验概 率分布模型； 利用二维闵可 夫斯基函数描述待预测二维岩心样本图像的几何结构特 征； 对待预测二维岩心样本图像的二维闵可夫斯基函数进行典型相关分析， 得到典型协变 量X； 将典型协变量X输入到训练好的先验概率分布模型中， 得到后验概率分布模型； 其中， 后验概率分布模型中概率最高点对应的电导率即为待预测 二维岩心样本的电导率可能性 最大的值。 2.根据权利要求1所述的一种微观储层岩石电导率预测方法， 其特征在于， 所述三维数 字岩心模型M的生成方法包括： 将二维岩石铸体薄片图像作为多点统计算法的训练图像， 利用多点统计算法生成多个 三维数字岩心模型I； 将多个三维数字岩心模型I作为生成对抗网络的训练图像， 生成多个三维数字岩心模 型M。 3.根据权利要求1所述的一种微观储层岩石电导 率预测方法， 其特 征在于， 还 包括： 对待预测二维岩心样本的二维闵可夫斯基函数进行典型相关分析之前， 先利用单分类 支持向量机模型， 对待预测岩心样本的二维闵可夫斯基函数进行异常值检测， 筛选掉不能 用于预测的异常样本 。 4.根据权利要求1所述的一种微观储层岩石电导率预测方法， 其特征在于， 所述二维闵 可夫斯基函数包括三个闵可夫斯基参数， 三个所述闵可夫斯基参数分别是孔隙度、 周长面 积比和欧拉数。 5.根据权利要求4所述的一种微观储层岩石电导率预测方法， 其特征在于， 所述孔隙度 的计算公式为： 其中， φ为 孔隙度， S P为孔隙的面积总和， SR为 二维岩心样本总面积； 所述周长面积比的计算公式为： 其中， W为周长面积比， CP为样本孔隙总周长， S P为孔隙的面积总和； 所述欧拉数的计算公式为： E＝NP‑HP 其中， E为欧拉数， NP为 孔隙连接体数， HP为 孔洞数。 6.根据权利要求1所述的一种微观储层岩石电导率预测方法， 其特征在于， 所述三维数 字岩心模型的二维闵可 夫斯基函数的计算方式包括：权　利　要　求　书 1/2 页 2 CN 115455772 A 2将三维数字岩心模型分为多层切片， 计算每层切片的二维闵可 夫斯基函数后取均值。 7.根据权利要求1所述的一种微观储层岩石电导 率预测方法， 其特 征在于， 还 包括： 对后验概 率分布结果设置 置信区间以衡量预测电导 率的合理性。 8.根据权利要求7所述的一种微观储层岩石电导率预测方法， 其特征在于， 所述置信区 间的计算公式如下 所示： 其中， CI指代置信区间， 为样本的均值， z为置信水平值， n为样本数量， z值 取0.95， S为 样本的标准差 。权　利　要　求　书 2/2 页 3 CN 115455772 A 3