CN113418848B - 岩心的渗透参数的确定方法、装置、设备及存储介质 - Google Patents

Abstract

本公开提供了一种岩心的渗透参数的确定方法、装置、设备及存储介质，属于油气勘探技术领域。该方法包括：获取岩心对应的三维数字图像，三维数字图像为通过对岩心进行扫描得到的多个体素对应的灰度值；根据多个体素对应的灰度值，确定岩心中的基质孔隙和裂缝；确定基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径，以及，确定裂缝的第二孔隙度；根据基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径，确定基质孔隙的第一渗透参数；根据第二孔隙度和第一渗透参数，确定岩心的第二渗透参数。由于三维数字图像和上述参数都是通过计算机设备模拟获取的，而不需要对岩心进行真实的实验，所以提高了确定岩心的第二渗透参数的效率。

Description

岩心的渗透参数的确定方法、装置、设备及存储介质

技术领域

本申请涉及油气勘探技术领域，特别涉及一种岩心的渗透参数的确定方法、装置、设备及存储介质。

背景技术

随着常规油气产量的日渐降低，大力勘探和开发致密油气藏等非常规油气资源已经是必然趋势。在对致密油气藏等非常规油气资源开发过程中，为了确定致密油气藏储层的采收率，需要确定储层对应的岩心的渗透率。

相关技术中，通过实验法确定岩心的渗透率；该方法为：将岩心放入密封的测试管内，且岩心与测试管的管壁贴合；从测试管的一端压入测试液，按照预设压力梯度增加测试液的压力，直至测试液从测试管的另一端流出；确定测试液从测试管的另一端流出时对应的测试液的目标压力，确定与目标压力匹配的岩心的渗透率。

但是，在上述相关技术中，当岩心为致密性岩心时，岩心的渗透率较低，测试液从测试管的另一端流出时对应的测试液的目标压力较高，而测试液体从测试管的一端压入时的初始压力是固定的，这样测试液体从初始压力按照预设压力梯度增加至目标压力所花费的时间较长，导致通过实验法确定岩心的渗透率的效率低。

发明内容

本申请实施例提供了一种岩心的渗透参数的确定方法、装置、设备及存储介质，可以提高确定岩心的渗透参数的效率。所述技术方案如下：

一方面，本申请提供了一种岩心的渗透参数的确定方法，所述方法包括：

获取岩心对应的三维数字图像，所述三维数字图像为通过对所述岩心进行扫描得到的多个体素对应的灰度值；

根据所述多个体素对应的灰度值，确定所述岩心中的基质孔隙和裂缝；

确定所述基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径，以及，确定所述裂缝的第二孔隙度；

根据所述基质孔隙的第一孔隙度、所述第一孔隙半径和所述第一喉道半径，确定所述基质孔隙的第一渗透参数；

根据所述第二孔隙度和所述第一渗透参数，确定所述岩心的第二渗透参数。

在一种可能的实现方式中，所述根据所述多个体素对应的灰度值，确定所述岩心中的基质孔隙和裂缝，包括：

根据所述多个体素对应的灰度值，对所述三维数字图像中的孔隙进行图像提取，得到所述岩心中的孔隙对应的孔隙三维图像；

确定所述孔隙三维图像中的每个孔隙的体积；

确定所述体积大于预设体积的孔隙为所述裂缝，确定所述体积小于所述预设体积的孔隙为所述基质孔隙。

在另一种可能的实现方式中，所述确定所述裂缝的第二孔隙度，包括：

确定所述裂缝的长度、开度和数量以及确定所述岩心的体积；

根据所述长度、所述开度、所述数量和所述岩心的体积，通过以下公式一，确定所述裂缝的第二孔隙度；

公式一：φ_f＝(L_f·W_f·N)/V

其中，φ_f表示所述第二孔隙度，L_f表示所述长度，W_f表示所述开度，N表示所述数量，V表示所述岩心的体积。

在另一种可能的实现方式中，所述基质孔隙包括多个子单元体，所述根据所述基质孔隙的第一孔隙度、所述第一孔隙半径和所述第一喉道半径，确定所述基质孔隙的第一渗透参数，包括：

从所述基质孔隙中选取多个第一目标子单元体；

对于每个第一目标子单元体，通过单相流模拟确定所述第一目标子单元体的第三渗透参数，以及，根据所述第一目标子单元体的孔隙分布，确定所述第一目标子单元体的第三孔隙度、第二孔隙半径和第二喉道半径，所述第一目标子单元体的个数大于三个；

根据所述每个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和所述第二喉道半径，确定第一关系数据，所述第一关系数据用于表示所述第一渗透参数、所述第一孔隙度、所述第一孔隙半径和所述第一喉道半径四者之间的关系；

根据所述第一孔隙度、所述第一孔隙半径、所述第一喉道半径和所述第一关系数据，确定所述基质孔隙的第一渗透参数。

在另一种可能的实现方式中，所述根据所述每个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径，确定第一关系数据，包括：

将所述每个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径输入以下公式二，确定所述公式二中的第一参数、第二参数和第三参数，得到所述第一关系数据；

公式二：k_a＝φ_e(ar_p+br_t)+c

其中，k_a表示所述基质孔隙的第一渗透参数，φ_e表示所述基质孔隙的第一孔隙度，r_p表示所述第一孔隙半径，r_t表示所述第一喉道半径，a表示所述第一参数，b表示所述第二参数，c表示所述第三参数。

在另一种可能的实现方式中，所述基质孔隙包括多个子单元体；所述根据所述第二孔隙度和所述第一渗透参数，确定所述岩心的第二渗透参数，包括：

从所述基质孔隙中选取多个第二目标子单元体；

对于每个第二目标子单元体，通过单相流模拟确定所述第二目标子单元体的第四渗透参数，以及，根据所述第二目标子单元体的裂缝分布，确定所述第二目标子单元体中裂缝的第四孔隙度，所述多个第二目标子单元体的个数大于三个；

根据所述第一渗透参数、所述每个第二目标子单元体的第四渗透参数和第四孔隙度，确定第二关系数据，所述第二关系数据用于表示所述第二渗透参数、所述第二孔隙度和所述第一渗透参数三者之间的关系；

根据所述第二孔隙度、所述第一渗透参数和所述第二关系数据，得到所述岩心的第二渗透参数。

在另一种可能的实现方式中，所述根据所述第一渗透参数、所述每个第二目标子单元体的第四渗透参数和第四孔隙度，确定所述第二关系数据，包括：

将所述第一渗透参数、所述每个第二目标子单元体的第四渗透参数和第四孔隙度输入以下公式三，确定所述公式三中的第四参数、第五参数和第六参数，得到所述第二关系数据；

公式三：k＝k_a(d·k_a·φ_f+eφ_f+f)

其中，k表示所述第二渗透参数，k_a表示所述第一渗透参数，φ_f表示所述裂缝的第二孔隙度，d表示所述第四参数，e表示所述第五参数，f表示所述第六参数。

另一方面，本申请提供了一种岩心的渗透参数的确定装置，所述装置包括：

获取模块，用于获取岩心对应的三维数字图像，所述三维数字图像为通过对所述岩心进行扫描得到的多个体素对应的灰度值；

第一确定模块，用于根据所述多个体素对应的灰度值，确定所述岩心中的基质孔隙和裂缝；

第二确定模块，用于确定所述基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径，以及，确定所述裂缝的第二孔隙度；

第三确定模块，用于根据所述基质孔隙的第一孔隙度、所述第一孔隙半径和所述第一喉道半径，确定所述基质孔隙的第一渗透参数；

第四确定模块，用于根据所述第二孔隙度和所述第一渗透参数，确定所述岩心的第二渗透参数。

在一种可能的实现方式中，所述第一确定模块，用于根据所述多个体素对应的灰度值，对所述三维数字图像中的孔隙进行图像提取，得到所述岩心中的孔隙对应的孔隙三维图像；确定所述孔隙三维图像中的每个孔隙的体积；确定所述体积大于预设体积的孔隙为所述裂缝，确定所述体积小于所述预设体积的孔隙为所述基质孔隙。

在另一种可能的实现方式中，所述第二确定模块，用于确定所述裂缝的长度、开度和数量以及确定所述岩心的体积；根据所述长度、所述开度、所述数量和所述岩心的体积，通过以下公式一，确定所述裂缝的第二孔隙度；

公式一：φ_f＝(L_f·W_f·N)/V

其中，φ_f表示所述第二孔隙度，L_f表示所述长度，W_f表示所述开度，N表示所述数量，V表示所述岩心的体积。

在另一种可能的实现方式中，所述基质孔隙包括多个子单元体；所述第三确定模块，用于从所述基质孔隙中选取多个第一目标子单元体；对于每个第一目标子单元体，通过单相流模拟确定所述第一目标子单元体的第三渗透参数，以及，根据所述第一目标子单元体的孔隙分布，确定所述第一目标子单元体的第三孔隙度、第二孔隙半径和第二喉道半径，所述第一目标子单元体的个数大于三个；根据所述每个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径，确定第一关系数据，所述第一关系数据用于表示所述第一渗透参数、所述第一孔隙度、所述第一孔隙半径和所述第一喉道半径四者之间的关系；根据所述第一孔隙度、所述第一孔隙半径、所述第一喉道半径和所述第一关系数据，确定所述基质孔隙的第一渗透参数。

在另一种可能的实现方式中，所述第三确定模块，用于将所述每个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径输入以下公式二，确定所述公式二中的第一参数、第二参数和第三参数，得到所述第一关系数据；

公式二：k_a＝φ_e(ar_p+br_t)+c

其中，k_a表示所述基质孔隙的第一渗透参数，φ_e表示所述基质孔隙的第一孔隙度，r_p表示所述第一孔隙半径，r_t表示所述第一喉道半径，a表示所述第一参数，b表示所述第二参数，c表示所述第三参数。

在另一种可能的实现方式中，所述基质孔隙包括多个子单元体；所述第四确定模块，用于从所述基质孔隙中选取多个第二目标子单元体；对于每个第二目标子单元体，通过单相流模拟确定所述第二目标子单元体的第四渗透参数，以及，根据所述第二目标子单元体的裂缝分布，确定所述第二目标子单元体中裂缝的第四孔隙度，所述多个第二目标子单元体的个数大于三个；根据所述第一渗透参数、所述每个第二目标子单元体的第四渗透参数和第四孔隙度，确定第二关系数据，所述第二关系数据用于表示所述第二渗透参数、所述第二孔隙度和所述第一渗透参数三者之间的关系；根据所述第二孔隙度、所述第一渗透参数和所述第二关系数据，得到所述岩心的第二渗透参数。

在另一种可能的实现方式中，所述第四确定模块，用于将所述第一渗透参数、所述每个第二目标子单元体的第四渗透参数和第四孔隙度输入以下公式三，确定所述公式三中的第四参数、第五参数和第六参数，得到所述第二关系数据；

公式三：k＝k_a(d·k_a·φ_f+eφ_f+f)

其中，k表示所述第二渗透参数，k_a表示所述第一渗透参数，φ_f表示所述裂缝的第二孔隙度，d表示所述第四参数，e表示所述第五参数，f表示所述第六参数。

另一方面，本申请实施例提供了一种计算机设备，所述计算机设备包括：处理器和存储器，所述存储器中存储有至少一条程序代码，所述至少一条程序代码由所述处理器加载并执行以实现上述任一可能实现方式所述的岩心的渗透参数的确定方法中所执行的操作。

另一方面，本申请实施例提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有至少一条程序代码，所述至少一条程序代码由处理器加载并执行以实现上述任一可能实现方式所述的岩心的渗透参数的确定方法中所执行的操作。

本申请实施例提供的技术方案的有益效果至少包括：

本申请实施例提供了一种岩心的渗透参数的确定方法，由于根据岩心对应的三维数字图像，确定岩心中基质孔隙的第一渗透参数以及岩心中裂缝的第二孔隙度，在通过基质孔隙的第一渗透参数和裂缝的第二孔隙度，确定岩心的第二渗透参数的过程中，三维数字图像、第一渗透参数和第二孔隙度都是通过计算机设备模拟获取的，不需要对岩心进行真实的实验，所以提高了确定岩心的第二渗透参数的效率。

附图说明

为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是根据本申请实施例提供的一种岩心的渗透参数的确定方法流程图；

图2是根据本申请实施例提供的一种岩心对应的二维切片灰度图像；

图3是根据本申请实施例提供的一种岩心中的孔隙对应的孔隙图像；

图4是根据本申请实施例提供的一种岩心中的裂缝对应的三维数字图像；

图5是根据本申请实施例提供的一种多个第一目标子单元体的模拟示意图；

图6是根据本申请实施例提供的一种第二目标子单元体对应的三维数字图像；

图7是根据本申请实施例提供的一种岩心的渗透参数的确定装置的框图；

图8是根据本申请实施例提供的一种计算机设备的结构框图。

具体实施方式

为使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请实施方式作进一步地详细描述。

图1是根据本申请实施例提供的一种岩心的渗透参数的确定方法流程图。

参见图1，该方法包括：

101、计算机设备获取岩心对应的三维数字图像，三维数字图像为通过对岩心进行扫描得到的多个体素对应的灰度值。

在本步骤中，岩心可以是任一形状的岩心。可选的，岩心的形状是规则形状，例如，圆柱形。岩心的体积可以是1立方厘米至1立方米之间的任一数值；例如，100立方厘米。在本申请实施例中，对岩心的体积和形状不作具体限定，可以根据需要进行设定并修改。

在一种可能的实现方式中，计算机设备获取岩心对应的三维数字图像的步骤为：计算机设备对岩心进行CT(Computed Tomography，电子计算机断层扫描)扫描，得到岩心对应的多个二维切片灰度图像，将多个二维切片灰度图像进行组合，得到岩心对应的多个体素，确定每个体素对应的灰度值，根据岩心对应的多个体素以及每个体素对应的灰度值，得到该岩心对应的三维数字图像。

在一种可能的实现方式中，参见图2，计算机设备对岩心进行CT扫描，得到岩心对应的二维切片灰度图像如图2所示。

102、计算机设备根据多个体素对应的灰度值，确定岩心中的基质孔隙和裂缝。

岩心包括骨架和孔隙，孔隙包括基质孔隙和裂缝，岩心的三维数字图像包括骨架相的体素和孔隙相的体素。其中，骨架相的体素对应的灰度值大于孔隙相的体素对应的灰度值。

在一种可能的实现方式中，计算机设备根据体素对应的灰度值，确定孔隙，进一步，通过孔隙的体积，将孔隙划分为基质孔隙和裂缝。相应的，本步骤为：计算机设备根据多个体素对应的灰度值，对三维数字图像中的孔隙进行图像提取，得到岩心中的孔隙对应的孔隙三维图像；确定孔隙三维图像中的每个孔隙的体积；确定体积大于预设体积的孔隙为裂缝，确定体积小于预设体积的孔隙为基质孔隙。

需要说明的一点是，为了使骨架相和孔隙相之间的边界更加清晰，计算机设备对三维数字图像中的孔隙进行图像提取之前，先对三维数字图像进行平滑处理。在一种可能的实现方式中，参见图3，计算机设备先对三维数字图像进行平滑处理，然后根据多个体素对应的灰度值，对三维数字图像中的孔隙进行图像提取，得到岩心中的孔隙对应的孔隙三维图像。该孔隙对应的孔隙三维图像如图3所示。

在一种可能的实现方式中，体积大于预设体积的孔隙为裂缝，体积小于预设体积的孔隙为基质孔隙，孔隙的体积与孔隙的体素数量相关。相应的，计算机设备确定孔隙三维图像中的每个孔隙的体积，确定体积大于预设体积的孔隙为裂缝，确定体积小于预设体积的孔隙为基质孔隙的步骤为：计算机设备确定孔隙三维图像中的每个孔隙的体素数量，确定体素数量大于预设数量的孔隙为裂缝，确定体素数量小于预设数量的孔隙为基质孔隙。其中，预设数量可以为10至150之间的任一数值，例如，50、60、70等。在本申请实施例中，对预设数量的数值不作具体限定，可以根据需要进行设定并修改。

在另一种可能的实现方式中，裂缝与基质孔隙的形状不同，计算机设备可以根据形状对孔隙进行划分。相应的，计算机设备根据多个体素对应的灰度值，确定岩心中的基质孔隙和裂缝的步骤为：计算机设备根据多个体素对应的灰度值，对三维数字图像中的孔隙进行图像提取，得到岩心中的孔隙对应的孔隙三维图像；确定孔隙三维图像中的每个孔隙的长宽比；确定长宽比大于预设比值的孔隙为裂缝，确定长宽比不大于预设比值的孔隙为基质孔隙。

其中，预设比值可以为3至10之间的任一数值，例如，3、4、5等。在本申请实施例中，对预设比值不作具体限定，可以根据需要进行设定并修改。

在另一种可能的实现方式中，计算机设备还可以根据形状和体积，对孔隙进行划分。相应的，计算机设备根据多个体素对应的灰度值，确定岩心中的基质孔隙和裂缝的步骤为：计算机设备根据多个体素对应的灰度值，对三维数字图像中的孔隙进行图像提取，得到岩心中的孔隙对应的孔隙三维图像；确定孔隙三维图像中的每个孔隙的长宽比和体积；确定长宽比大于预设比值且体积大于预设体积的孔隙为裂缝，确定长宽比不大于预设比值且体积不大于预设体积的孔隙为基质孔隙。

在一种可能的实现方式中，参见图4，计算机设备确定体素数量大于预设数量的孔隙为裂缝之后，对孔隙三维图像中的裂缝进行图像提取，得到裂缝三维图像；该裂缝三维图像如图4所示。

103、计算机设备确定基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径，以及，确定裂缝的第二孔隙度。

在一种可能的实现方式中，计算机设备确定基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径的步骤为：计算机设备确定基质孔隙的体积以及岩心的体积，确定基质孔隙的体积与岩心的体积的比值为基质孔隙的第一孔隙度；以及，计算机设备确定基质孔隙中每个孔隙的半径，确定半径的平均值为基质孔隙的第一孔隙半径；以及，计算机设备确定基质孔隙中每个孔隙的喉道半径，确定喉道半径的平均值为基质孔隙的第一喉道半径。

在一种可能的实现方式中，计算机设备根据体素数量，确定基质孔隙的体积以及岩心的体积。相应的，计算机设备确定基质孔隙的体积以及岩心的体积，确定基质孔隙的体积与岩心的体积的比值为基质孔隙的第一孔隙度的步骤为：计算机设备确定基质孔隙的体素数量以及岩心的体素数量，将基质孔隙的体素数量与岩心的体素数量的比值作为基质孔隙的第一孔隙度。

在一种可能的实现方式中，计算机设备确定裂缝的第二孔隙度的步骤为：计算机设备确定裂缝的长度、开度和数量以及确定岩心的体积；根据裂缝的长度、开度、数量和岩心的体积，通过以下公式一，确定裂缝的第二孔隙度；

公式一：φ_f＝(L_f·W_f·N)/V

其中，φ_f表示第二孔隙度，L_f表示长度，W_f表示开度，N表示数量，V岩心的体积。

在另一种可能的实现方式中，计算机设备根据裂缝的体素数量，确定裂缝的第二孔隙度。相应的，计算机设备确定裂缝的第二孔隙度的步骤为：计算机设备确定裂缝的体素数量和岩心的体素数量，将裂缝的体素数量和岩心的体素数量的比值作为裂缝的第二孔隙度。

104、计算机设备根据基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径，确定基质孔隙的第一渗透参数。

在一种可能的实现方式中，基质孔隙包括多个子单元体。相应的，本步骤为：计算机设备从基质孔隙中选取多个第一目标子单元体；对于每个第一目标子单元体，通过单相流模拟确定第一目标子单元体的第三渗透参数，以及，根据第一目标子单元体的孔隙分布，确定第一目标子单元体的第三孔隙度、第二孔隙半径和第二喉道半径，第一目标子单元体的个数大于三个；根据每个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径，确定第一关系数据，第一关系数据用于表示第一渗透参数、第一孔隙度、第一孔隙半径和第一喉道半径四者之间的关系；根据第一孔隙度、第一孔隙半径、第一喉道半径和第一关系数据，确定基质孔隙的第一渗透参数。

需要说明的一点是，在确定基质孔隙的第一渗透参数的过程中，为了降低缝隙孔隙的干扰，第一目标子单元体可以为只包含基质孔隙的子单元体。在一种可能的实现方式中，计算机设备从基质孔隙中选取多个第一目标子单元体的步骤为：计算机设备确定每个子单元体内包含的孔隙类型，从多个子单元体中选取孔隙类型只包含基质孔隙的多个第一目标子单元体。

其中，第一目标子单元体的形状可以为圆柱体、也可以是球体，还可以是长方体。在一种可能的实现方式中，参见图5，第一目标子单元体的形状为正方体，也即，图5中虚线标出的正方体为第一目标子单元体。正方体的边长可以是100体素至500体素之间的任一数值；例如，100体素、200体素、300体素等。在本申请实施例中，对第一目标子单元体的尺寸不作具体限定，可以根据需要进行设定并修改。例如，每个第一目标子单元体的尺寸为200体素*200体素*200体素。

需要说明的一点是，计算机设备通过对多个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径，进行回归分析得到第一关系数据。其中，第一目标子单元体的个数越多，根据每个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径，确定的第一关系数据的精确性越高，计算机设备对应的运算量越大。

在一种可能的实现方式中，计算机设备通过单相流模拟确定第一目标子单元体的第三渗透参数的步骤为：计算机设备通过LBM(Lattice Boltzmann Method，格子玻尔兹曼方法)单相流模拟，确定第一目标子单元体中基质孔隙对应的速度场，根据第一目标子单元体中基质孔隙对应的速度场，确定第一目标子单元体的第三渗透参数。

在一种可能的实现方式中，计算机设备根据多个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径，确定第一关系数据的步骤为：计算机设备将每个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径输入以下公式二，确定公式二中的第一参数、第二参数和第三参数，得到第一关系数据；

公式二：k_a＝φ_e(ar_p+br_t)+c

其中，k_a表示基质孔隙的第一渗透参数，φ_e表示基质孔隙的第一孔隙度，r_p表示第一孔隙半径，r_t表示第一喉道半径，a表示第一参数，b表示第二参数，c表示第三参数。

需要说明的一点是，第一关系数据为k_a＝φ_e(ar_p+br_t)+c；其中，第一参数a、第二参数b和第三参数c为已知参数，基质孔隙的第一孔隙度φ_e、第一孔隙半径r_p和第一喉道半径r_t为自变量，基质孔隙的第一渗透参数k_a为因变量。

105、计算机设备根据第二孔隙度和第一渗透参数，确定岩心的第二渗透参数。

在一种可能的实现方式中，基质孔隙包括多个子单元体；相应的，本步骤为：计算机设备从基质孔隙中选取多个第二目标子单元体；对于每个第二目标子单元体，通过单相流模拟确定第二目标子单元体的第四渗透参数，以及，根据第二目标子单元体的裂缝分布，确定第二目标子单元体中裂缝的第四孔隙度，多个第二目标子单元体的个数大于三个；根据第一渗透参数、每个第二目标子单元体的第四渗透参数和第四孔隙度，确定第二关系数据，第二关系数据用于表示第二渗透参数、第二孔隙度和第一渗透参数三者之间的关系；根据第二孔隙度、第一渗透参数和第二关系数据，得到岩心的第二渗透参数。

需要说明的一点是，在确定岩心的第二渗透参数的过程中，为了综合基质孔隙和缝隙孔隙两者对该岩心的第二渗透参数的影响，第二目标子单元体为包含基质孔隙和缝隙孔隙的子单元体。

在一种可能的实现方式中，计算机设备从基质孔隙中选取多个第二目标子单元体的步骤为：计算机设备确定每个子单元体内包含的孔隙类型，从多个子单元体中选取孔隙类型包含基质孔隙和缝隙孔隙的多个第二目标子单元体。

在一种可能的实现方式中，计算机设备通过单相流模拟确定第二目标子单元体的第四渗透参数的步骤为：计算机设备通过LBM单相流模拟，确定第二目标子单元体中缝隙对应的速度场，根据第二目标子单元体中缝隙对应的速度场，确定第二目标子单元体的第四渗透参数。

在一种可能的实现方式中，计算机设备根据第二目标子单元体的裂缝分布，确定第二目标子单元体中裂缝的第四孔隙度的步骤为：计算机设备根据第二目标子单元体的裂缝分布，确定第二子单元体中裂缝的长度、开度和数量，以及，确定第二单元体的体积；根据第二目标子单元体中裂缝长度、开度、数量和第二单元体的体积，通过以下公式二，确定第二目标子单元体中裂缝的第四孔隙度，得到每个第二目标子单元体中裂缝的第四孔隙度；

公式二：φ_f＝(L_f·W_f·N)/V

其中，φ_f表示第四孔隙度，L_f表示裂缝长度，W_f表示裂缝开度，N表示裂缝的数量，V表示第二目标子单元体的体积。

需要说明的另一点是，多个第二目标子单元体也可以是从多个第一目标子单元体中选取的，相应的，计算机设备从基质孔隙中选取多个第二目标子单元体的步骤为：计算机设备从多个第二目标子单元体中选取的多个第二目标子单元体，该第二目标子单元体数量至少为三个；计算机设备向每个第一目标子单元体中加入模拟裂缝，得到多个第二目标子单元体。例如，参见图6，计算机设备向第一目标子单元体对应的三维数字图像中加入模拟裂缝，得到第二目标子单元体对应的三维数字图像。

在本步骤中，计算机设备根据第一渗透参数、每个第二目标子单元体的第四渗透参数和第四孔隙度，确定第二关系数据的步骤为：计算机设备将第一渗透参数、每个第二目标子单元体的第四渗透参数和第四孔隙度输入以下公式三，确定公式三中的第四参数、第五参数和第六参数，得到第二关系数据；

公式三：k＝k_a(d·k_a·φ_f+eφ_f+f)

其中，k表示第二渗透参数，k_a表示第一渗透参数，φ_f表示裂缝的第二孔隙度，d表示第四参数，e表示第五参数，f表示第六参数。

需要说明的一点是，第二关系数据为k＝k_a(d·k_a·φ_f+eφ_f+f)；其中，第四参数d、第五参数e和第六参数f为已知参数，第一渗透参数k_a和第二孔隙度φ_f为自变量，岩心的第二渗透参数k为因变量。在确定公式三中的第四参数、第五参数和第六参数的过程中，自变量φ_f为第四孔隙度，因变量k为第四渗透参数。也即，通过第二目标子单元体的第四渗透参数和第四孔隙度，确定公式三中的参数，进而得到第二渗透参数、第二孔隙度和第一渗透参数三者之间的具体关系。

在一种可能的实现方式中，计算机设备根据第二孔隙度、第一渗透参数和第二关系数据，得到岩心的第二渗透参数的步骤为：计算机设备将第二孔隙度、第一渗透参数输入公式三，得到岩心的第二渗透参数。

需要说明的另一点是，计算机设备确定岩心的第二渗透参数之后，还可以根据第一渗透参数和第二渗透参数，确定裂缝对整个岩心系统的渗透率的贡献百分比。贡献百分比能够表征裂缝对岩心的渗透率的影响程度。

在一种可能的实现方式中，计算机设备根据第一渗透参数和第二渗透参数，确定裂缝对整个岩心系统的渗透率的贡献百分比的步骤为：计算机设备根据第一渗透参数和第二渗透参数，通过以下公式四，确定裂缝对整个岩心系统的渗透率的贡献百分比；

公式四：η＝(k-k_a)/k×100％

其中，η表示贡献百分比，k表示第二渗透参数，k_a表示第一渗透参数。

本申请实施例提供了一种岩心的渗透参数的确定方法，由于根据岩心对应的三维数字图像，确定岩心中基质孔隙的第一渗透参数以及岩心中裂缝的第二孔隙度，在通过基质孔隙的第一渗透参数和裂缝的第二孔隙度，确定岩心的第二渗透参数的过程中，三维数字图像、第一渗透参数和第二孔隙度都是通过计算机设备模拟获取的，不需要对岩心进行真实的实验，所以提高了确定岩心的第二渗透参数的效率。

图7是根据一示例性实施例示出的一种岩心的渗透参数的确定装置的框图。参见图7，该装置包括：

获取模块701，用于获取岩心对应的三维数字图像，三维数字图像为通过对岩心进行扫描得到的多个体素对应的灰度值；

第一确定模块702，用于根据多个体素对应的灰度值，确定岩心中的基质孔隙和裂缝；

第二确定模块703，用于确定基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径，以及，确定裂缝的第二孔隙度；

第三确定模块704，用于根据基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径，确定基质孔隙的第一渗透参数；

第四确定模块705，用于根据第二孔隙度和第一渗透参数，确定岩心的第二渗透参数。

在一种可能的实现方式中，第一确定模块702，用于根据多个体素对应的灰度值，对三维数字图像中的孔隙进行图像提取，得到岩心中的孔隙对应的孔隙三维图像；确定孔隙三维图像中的每个孔隙的体积；确定体积大于预设体积的孔隙为裂缝，确定体积小于预设体积的孔隙为基质孔隙。

在另一种可能的实现方式中，第二确定模块703，用于确定裂缝的长度、开度和数量以及确定岩心的体积；根据长度、开度、数量和岩心的体积，通过以下公式一，确定裂缝的第二孔隙度；

公式一：φ_f＝(L_f·W_f·N)/V

其中，φ_f表示第二孔隙度，L_f表示长度，W_f表示开度，N表示数量，V表示岩心的体积。

在另一种可能的实现方式中，基质孔隙包括多个子单元体；第三确定模块704，用于从基质孔隙中选取多个第一目标子单元体；对于每个第一目标子单元体，通过单相流模拟确定第一目标子单元体的第三渗透参数，以及，根据第一目标子单元体的孔隙分布，确定第一目标子单元体的第三孔隙度、第二孔隙半径和第二喉道半径，第一目标子单元体的个数大于三个；根据每个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径，确定第一关系数据，第一关系数据用于表示第一渗透参数、第一孔隙度、第一孔隙半径和第一喉道半径四者之间的关系；根据第一孔隙度、第一孔隙半径、第一喉道半径和第一关系数据，确定基质孔隙的第一渗透参数。

在另一种可能的实现方式中，第三确定模块704，用于将每个第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径输入以下公式二，确定公式二中的第一参数、第二参数和第三参数，得到第一关系数据；

公式二：k_a＝φ_e(ar_p+br_t)+c

其中，k_a表示基质孔隙的第一渗透参数，φ_e表示基质孔隙的第一孔隙度，r_p表示第一孔隙半径，r_t表示第一喉道半径，a表示第一参数，b表示第二参数，c表示第三参数。

在另一种可能的实现方式中，基质孔隙包括多个子单元体；第四确定模块705，用于从基质孔隙中选取多个第二目标子单元体；对于每个第二目标子单元体，通过单相流模拟确定第二目标子单元体的第四渗透参数，以及，根据第二目标子单元体的裂缝分布，确定第二目标子单元体中裂缝的第四孔隙度，多个第二目标子单元体的个数大于三个；根据第一渗透参数、每个第二目标子单元体的第四渗透参数和第四孔隙度，确定第二关系数据，第二关系数据用于表示第二渗透参数、第二孔隙度和第一渗透参数三者之间的关系；根据第二孔隙度、第一渗透参数和第二关系数据，得到岩心的第二渗透参数。

在另一种可能的实现方式中，第四确定模块705，用于将第一渗透参数、每个第二目标子单元体的第四渗透参数和第四孔隙度输入以下公式三，确定公式三中的第四参数、第五参数和第六参数，得到第二关系数据；

公式三：k＝k_a(d·k_a·φ_f+eφ_f+f)

其中，k表示第二渗透参数，k_a表示第一渗透参数，φ_f表示裂缝的第二孔隙度，d表示第四参数，e表示第五参数，f表示第六参数。

本申请实施例提供了一种岩心的渗透参数的确定方法，由于根据岩心对应的三维数字图像，确定岩心中基质孔隙的第一渗透参数以及岩心中裂缝的第二孔隙度，在通过基质孔隙的第一渗透参数和裂缝的第二孔隙度，确定岩心的第二渗透参数的过程中，三维数字图像、第一渗透参数和第二孔隙度都是通过计算机设备模拟获取的，不需要对岩心进行真实的实验，所以提高了确定岩心的第二渗透参数的效率。

图8示出了本发明一个示例性实施例提供的终端800的结构框图。该终端800可以是：智能手机、平板电脑、MP3播放器(Moving Picture Experts Group Audio Layer III，动态影像专家压缩标准音频层面3)、MP4(Moving Picture Experts Group Audio LayerIV，动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。终端800还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。

通常，终端800包括有：处理器801和存储器802。

处理器801可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器801可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field－Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器801也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称CPU(Central ProcessingUnit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器801可以在集成有GPU(Graphics Processing Unit，图像处理器)，GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器801还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关机器学习的计算操作。

存储器802可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器802还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器802中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器801所执行以实现本申请中方法实施例提供的岩心的渗透参数的确定方法。

在一些实施例中，终端800还可选包括有：外围设备接口803和至少一个外围设备。处理器801、存储器802和外围设备接口803之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口803相连。具体地，外围设备包括：射频电路804、显示屏805、摄像头806、音频电路807、定位组件808和电源809中的至少一种。

外围设备接口803可被用于将I/O(Input/Output，输入/输出)相关的至少一个外围设备连接到处理器801和存储器802。在一些实施例中，处理器801、存储器802和外围设备接口803被集成在同一芯片或电路板上；在一些其他实施例中，处理器801、存储器802和外围设备接口803中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。

射频电路804用于接收和发射RF(Radio Frequency，射频)信号，也称电磁信号。射频电路804通过电磁信号与通信网络以及其他通信设备进行通信。射频电路804将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路804包括：天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路804可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：城域网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity，无线保真)网络。在一些实施例中，射频电路804还可以包括NFC(Near Field Communication，近距离无线通信)有关的电路，本申请对此不加以限定。

显示屏805用于显示UI(User Interface，用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏805是触摸显示屏时，显示屏805还具有采集在显示屏805的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器801进行处理。此时，显示屏805还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏805可以为一个，设置终端800的前面板；在另一些实施例中，显示屏805可以为至少两个，分别设置在终端800的不同表面或呈折叠设计；在再一些实施例中，显示屏805可以是柔性显示屏，设置在终端800的弯曲表面上或折叠面上。甚至，显示屏805还可以设置成非矩形的不规则图形，也即异形屏。显示屏805可以采用LCD(LiquidCrystal Display，液晶显示屏)、OLED(Organic Light-Emitting Diode,有机发光二极管)等材质制备。

摄像头组件806用于采集图像或视频。可选地，摄像头组件806包括前置摄像头和后置摄像头。通常，前置摄像头设置在终端的前面板，后置摄像头设置在终端的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件806还可以包括闪光灯。闪光灯可以是单色温闪光灯，也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同色温下的光线补偿。

音频电路807可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器801进行处理，或者输入至射频电路804以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在终端800的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器801或射频电路804的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路807还可以包括耳机插孔。

定位组件808用于定位终端800的当前地理位置，以实现导航或LBS(LocationBased Service，基于位置的服务)。定位组件808可以是基于美国的GPS(GlobalPositioning System，全球定位系统)、中国的北斗系统、俄罗斯的格雷纳斯系统或欧盟的伽利略系统的定位组件。

电源809用于为终端800中的各个组件进行供电。电源809可以是交流电、直流电、一次性电池或可充电电池。当电源809包括可充电电池时，该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。

在一些实施例中，终端800还包括有一个或多个传感器810。该一个或多个传感器810包括但不限于：加速度传感器811、陀螺仪传感器812、压力传感器813、指纹传感器814、光学传感器815以及接近传感器816。

加速度传感器811可以检测以终端800建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器811可以用于检测重力加速度在三个坐标轴上的分量。处理器801可以根据加速度传感器811采集的重力加速度信号，控制显示屏805以横向视图或纵向视图进行用户界面的显示。加速度传感器811还可以用于游戏或者用户的运动数据的采集。

陀螺仪传感器812可以检测终端800的机体方向及转动角度，陀螺仪传感器812可以与加速度传感器811协同采集用户对终端800的3D动作。处理器801根据陀螺仪传感器812采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。

压力传感器813可以设置在终端800的侧边框和/或显示屏805的下层。当压力传感器813设置在终端800的侧边框时，可以检测用户对终端800的握持信号，由处理器801根据压力传感器813采集的握持信号进行左右手识别或快捷操作。当压力传感器813设置在显示屏805的下层时，由处理器801根据用户对显示屏805的压力操作，实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。

指纹传感器814用于采集用户的指纹，由处理器801根据指纹传感器814采集到的指纹识别用户的身份，或者，由指纹传感器814根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器801授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器814可以被设置终端800的正面、背面或侧面。当终端800上设置有物理按键或厂商Logo时，指纹传感器814可以与物理按键或厂商Logo集成在一起。

光学传感器815用于采集环境光强度。在一个实施例中，处理器801可以根据光学传感器815采集的环境光强度，控制显示屏805的显示亮度。具体地，当环境光强度较高时，调高显示屏805的显示亮度；当环境光强度较低时，调低显示屏805的显示亮度。在另一个实施例中，处理器801还可以根据光学传感器815采集的环境光强度，动态调整摄像头组件806的拍摄参数。

接近传感器816，也称距离传感器，通常设置在终端800的前面板。接近传感器816用于采集用户与终端800的正面之间的距离。在一个实施例中，当接近传感器816检测到用户与终端800的正面之间的距离逐渐变小时，由处理器801控制显示屏805从亮屏状态切换为息屏状态；当接近传感器816检测到用户与终端800的正面之间的距离逐渐变大时，由处理器801控制显示屏805从息屏状态切换为亮屏状态。

本领域技术人员可以理解，图8中示出的结构并不构成对终端800的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。

在示例性实施例中，还提供了一种包括程序代码的存储介质，例如包括程序代码的存储器，上述程序代码可由装置的处理器执行以完成上述方法。可选地，存储介质可以是非临时性计算机可读存储介质，例如，所述非临时性计算机可读存储介质可以是ROM(Read-Only Memory，只读存储器)、RAM(Random Access Memory，随机存取存储器)、CD-ROM(Compact Disc Read-Only Memory，只读光盘)、磁带、软盘和光数据存储设备等。

以上所述仅为本申请的可选实施例，并不用以限制本申请，凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (6)

1.一种岩心的渗透参数的确定方法，其特征在于，所述方法包括：

获取岩心对应的三维数字图像，所述三维数字图像为通过对所述岩心进行扫描得到的多个体素对应的灰度值；

根据所述多个体素对应的灰度值，确定所述岩心中的基质孔隙和裂缝；所述基质孔隙包括多个子单元体；

确定所述基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径，以及，确定所述裂缝的第二孔隙度；

从所述多个子单元体中选取多个第一目标子单元体；

对于每个第一目标子单元体，通过单相流模拟确定所述第一目标子单元体的第三渗透参数，以及，根据所述第一目标子单元体的孔隙分布，确定所述第一目标子单元体的第三孔隙度、第二孔隙半径和第二喉道半径，所述第一目标子单元体的个数大于三个；

将所述第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径输入以下公式二，确定所述公式二中的第一参数、第二参数和第三参数，得到第一关系数据；所述第一关系数据用于表示所述基质孔隙的第一渗透参数、所述第一孔隙度、所述第一孔隙半径和所述第一喉道半径四者之间的关系；

公式二：k_a＝φ_e(ar_p+br_t)+c

其中，k_a表示所述第一渗透参数，φ_e表示所述第一孔隙度，r_p表示所述第一孔隙半径，r_t表示所述第一喉道半径，a表示所述第一参数，b表示所述第二参数，c表示所述第三参数；

根据所述第一孔隙度、所述第一孔隙半径、所述第一喉道半径和所述第一关系数据，确定所述第一渗透参数；

向所述每个第一目标子单元体中加入模拟裂缝，得到多个第二目标子单元体；

对于每个第二目标子单元体，通过单相流模拟确定所述第二目标子单元体的第四渗透参数，以及，根据所述第二目标子单元体的裂缝分布，确定所述第二目标子单元体中裂缝的第四孔隙度；

将所述第一渗透参数、所述第二目标子单元体的第四渗透参数和第四孔隙度输入以下公式三，确定所述公式三中的第四参数、第五参数和第六参数，得到第二关系数据；所述第二关系数据用于表示所述岩心的第二渗透参数、所述第二孔隙度和所述第一渗透参数三者之间的关系；

公式三：k＝k_a(d·k_a·φ_f+eφ_f+f)

其中，k表示所述第二渗透参数，k_a表示所述第一渗透参数，φ_f表示所述第二孔隙度，d表示所述第四参数，e表示所述第五参数，f表示所述第六参数；

将所述第二孔隙度和所述第一渗透参数代入所述第二关系数据中，得到所述第二渗透参数。

2.根据权利要求1所述的方法，其特征在于，所述根据所述多个体素对应的灰度值，确定所述岩心中的基质孔隙和裂缝，包括：

根据所述多个体素对应的灰度值，对所述三维数字图像中的孔隙进行图像提取，得到所述岩心中的孔隙对应的孔隙三维图像；

确定所述孔隙三维图像中的每个孔隙的体积；

确定所述体积大于预设体积的孔隙为所述裂缝，确定所述体积小于所述预设体积的孔隙为所述基质孔隙。

3.根据权利要求1所述的方法，其特征在于，所述确定所述裂缝的第二孔隙度，包括：

确定所述裂缝的长度、开度和数量以及确定所述岩心的体积；

根据所述长度、所述开度、所述数量和所述岩心的体积，通过以下公式一，确定所述裂缝的第二孔隙度；

公式一：φ_f＝(L_f·W_f·N)/V

其中，φ_f表示所述第二孔隙度，L_f表示所述长度，W_f表示所述开度，N表示所述数量，V表示所述岩心的体积。

4.一种岩心的渗透参数的确定装置，其特征在于，所述装置包括：

获取模块，用于获取岩心对应的三维数字图像，所述三维数字图像为通过对所述岩心进行扫描得到的多个体素对应的灰度值；

第一确定模块，用于根据所述多个体素对应的灰度值，确定所述岩心中的基质孔隙和裂缝；所述基质孔隙包括多个子单元体；

第二确定模块，用于确定所述基质孔隙的第一孔隙度、第一孔隙半径和第一喉道半径，以及，确定所述裂缝的第二孔隙度；

第三确定模块，用于从所述多个子单元体中选取多个第一目标子单元体；对于每个第一目标子单元体，通过单相流模拟确定所述第一目标子单元体的第三渗透参数，以及，根据所述第一目标子单元体的孔隙分布，确定所述第一目标子单元体的第三孔隙度、第二孔隙半径和第二喉道半径，所述第一目标子单元体的个数大于三个；将所述第一目标子单元体的第三渗透参数、第三孔隙度、第二孔隙半径和第二喉道半径输入以下公式二，确定所述公式二中的第一参数、第二参数和第三参数，得到第一关系数据；所述第一关系数据用于表示所述基质孔隙的第一渗透参数、所述第一孔隙度、所述第一孔隙半径和所述第一喉道半径四者之间的关系；公式二：k_a＝φ_e(ar_p+br_t)+c；其中，k_a表示所述第一渗透参数，φ_e表示所述第一孔隙度，r_p表示所述第一孔隙半径，r_t表示所述第一喉道半径，a表示所述第一参数，b表示所述第二参数，c表示所述第三参数；根据所述第一孔隙度、所述第一孔隙半径、所述第一喉道半径和所述第一关系数据，确定所述第一渗透参数；

第四确定模块，用于向所述每个第一目标子单元体中加入模拟裂缝，得到多个第二目标子单元体；对于每个第二目标子单元体，通过单相流模拟确定所述第二目标子单元体的第四渗透参数，以及，根据所述第二目标子单元体的裂缝分布，确定所述第二目标子单元体中裂缝的第四孔隙度；将所述第一渗透参数、所述第二目标子单元体的第四渗透参数和第四孔隙度输入以下公式三，确定所述公式三中的第四参数、第五参数和第六参数，得到第二关系数据；所述第二关系数据用于表示所述岩心的第二渗透参数、所述第二孔隙度和所述第一渗透参数三者之间的关系；公式三：k＝k_a(d·k_a·φ_f+eφ_f+f)；其中，k表示所述第二渗透参数，k_a表示所述第一渗透参数，φ_f表示所述第二孔隙度，d表示所述第四参数，e表示所述第五参数，f表示所述第六参数；将所述第二孔隙度和所述第一渗透参数代入所述第二关系数据中，得到所述第二渗透参数。

5.一种计算机设备，其特征在于，所述计算机设备包括：

处理器和存储器，所述存储器中存储有至少一条程序代码，所述至少一条程序代码由所述处理器加载并执行以实现权利要求1至3任一项所述的岩心的渗透参数的确定方法中所执行的操作。

6.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有至少一条程序代码，所述至少一条程序代码由处理器加载并执行以实现如权利要求1至3任一项所述的岩心的渗透参数的确定方法中所执行的操作。

Images