孔喉特征参数的获取方法、装置、计算机设备及存储介质

摘要

本申请是关于一种孔喉特征参数的获取方法、装置、计算机设备及存储介质，涉及油气勘探技术领域。该方法包括：通过岩石样本的特征信息和所述岩石样本的第一类型孔洞的特征数据，获得所述第一类型孔洞的进汞饱和度；获取指示岩石样本中第二类型孔洞的体积占比的第一毛管压力曲线；通过所述第一毛管压力曲线、所述进汞饱和度以及所述第一类型孔洞的特征数据，获得第二毛管压力曲线；再根据所述第二毛管压力曲线获取所述岩石样本的孔喉特征参数。本方案提升了毛管压力曲线评价岩石孔喉结构的准确性。

说明书

技术领域

本申请涉及油气勘探技术领域，特别涉及一种孔喉特征参数的获取方法、装置、计算机设备及存储介质。

背景技术

毛管压力曲线可以反映岩石中的不同孔隙的分布情况，是评价岩石孔喉结构和储层岩石渗流能力的重要依据。

在相关的技术中，由于核磁共振弛豫时间T2与孔喉的半径呈正相关关系，技术人员通过核磁共振弛豫时间T2的分布获取毛管压力曲线。

目前相关的技术中，核磁共振受限于顺磁离子扩散和岩石内部磁场梯度等影响因素，无法评价岩石表面未饱和水的溶洞性质，导致得到的孔喉特征参数的准确性较低。

发明内容

本申请实施例提供了一种孔喉特征参数的获取方法、装置、计算机设备及存储介质，可以通过获得的第二毛管压力曲线，有效地评价岩石孔喉结构，得到准确的孔喉特征参数，该技术方案如下：

一方面，提供了一种孔喉特征参数的获取方法，所述方法包括：

获取岩石样本的特征信息，所述特征信息包括所述岩石样本的长度、直径、孔隙度和渗透率；

对所述岩石样本进行特征提取处理，获得所述岩石样本的第一类型孔洞的特征数据，所述第一类型孔洞包括表面溶洞和尺寸大于尺寸阈值的孔；

根据所述特征信息和所述第一类型孔洞的特征数据，获取所述第一类型孔洞的进汞饱和度；

获取所述岩石样本对应的第一毛管压力曲线，所述第一毛管压力曲线用于指示岩石样本中第二类型孔洞的体积占比，所述第二类型孔洞的体积占比是所述第二类型孔洞的体积占所述岩石样本的总孔隙体积的比例，所述第二类型孔洞是尺寸小于尺寸阈值的孔隙；

根据所述第一毛管压力曲线、所述进汞饱和度以及所述第一类型孔洞的特征数据，获取第二毛管压力曲线；

根据所述第二毛管压力曲线获取所述岩石样本的孔喉特征参数，所述孔喉特征参数包括最大孔隙半径、平均孔隙半径、中值、孔喉半径均值、分选系数、歪度、峰态、均质系数。

又一方面，提供了一种孔喉特征参数的获取装置，所述装置包括：

特征信息获取模块，用于获取岩石样本的特征信息，所述特征信息包括所述岩石样本的长度、直径、孔隙度和渗透率；

特征数据获取模块，用于对所述岩石样本进行特征提取处理，获得所述岩石样本的第一类型孔洞的特征数据，所述第一类型孔洞包括表面溶洞和尺寸大于尺寸阈值的孔；

进汞饱和度获取模块，用于根据所述特征信息和所述第一类型孔洞的特征数据，获取所述第一类型孔洞的进汞饱和度；

第一毛管压力曲线获取模块，用于获取所述岩石样本对应的第一毛管压力曲线，所述第一毛管压力曲线用于指示岩石样本中第二类型孔洞的体积占比，所述第二类型孔洞的体积占比是所述第二类型孔洞的体积占所述岩石样本的总孔隙体积的比例，所述第二类型孔洞是尺寸小于尺寸阈值的孔隙；

第二毛管压力曲线获取模块，用于根据所述第一毛管压力曲线、所述进汞饱和度以及所述第一类型孔洞的特征信息，获取第二毛管压力曲线；

孔喉特征参数获取模块，用于根据所述第二毛管压力曲线获取所述岩石样本的孔喉特征参数，所述孔喉特征参数包括最大孔隙半径、平均孔隙半径、中值、孔喉半径均值、分选系数、歪度、峰态、均质系数。

可选的，所述第二毛管压力曲线获取模块，用于通过所述进汞饱和度和所述第一类型孔洞的特征数据，修正所述第一毛管压力曲线，获取所述第二毛管压力曲线。

可选的，所述进汞饱和度获取模块，包括：

总孔隙体积获取子模块，用于根据所述特征信息，获得所述岩石样本的总孔隙体积；

数据处理子模块，用于根据所述第一类型孔洞的特征数据和所述岩石样本的总孔隙体积，获得所述进汞饱和度。

可选的，所述进汞饱和度处理子模块，包括：

体积获取单元，用于根据所述第一类型孔洞的特征数据，获得所述第一类型孔洞的体积；

进汞饱和度计算单元，用于将所述第一类型孔洞的体积，除以所述岩石样本的总孔隙体积，作为所述进汞饱和度。

可选的，所述特征数据获取模块，包括：

三维数字岩心获取子模块，用于对所述岩石样本进行微CT扫描，获得所述岩石样本的三维数字岩心；

三维数字岩心处理子模块，用于根据所述三维数字岩心，获得所述岩石样本的第一类型孔洞的特征数据，所述第一类型孔洞的特征数据包括所述第一类型孔洞的直径和体积。

可选的，所述孔喉特征参数获取模块，用于对所述第二毛管压力曲线中累积的所述进汞饱和度进行归一化处理，获得所述岩石样本的孔喉特征参数。

可选的，所述第一毛管压力曲线获取模块，包括：

高压压汞处理子模块，通过高压压汞实验获得所述岩石样本的毛管压力曲线；

麻皮效应校正子模块，用于对所述毛管压力曲线进行麻皮效应校正，获得所述第一毛管压力曲线。

再一方面，提供了一种计算机设备，所述计算机设备包含处理器和存储器，所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集由所述处理器加载并执行以实现上述的孔喉特征参数的获取方法。

又一方面，提供了一种计算机可读存储介质，所述存储介质中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集由处理器加载并执行以实现上述孔喉特征参数的获取方法。

又一方面，提供了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算机设备执行上述孔喉特征参数的获取方法。

本申请提供的技术方案可以包括以下有益效果：

通过岩石样本的特征信息和该岩石样本中尺寸大于尺寸阈值的第一类型孔洞的特征数据，获得该第一类型孔洞的进汞饱和度；并且获取指示岩石样本中第二类型孔洞的体积占比的第一毛管压力曲线；然后再通过该第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据，获得第二毛管压力曲线；再根据该第二毛管压力曲线获取该岩石样本的孔喉特征参数。由于第一毛管压力曲线能够体现该第二类型孔洞的特征，而进汞饱和度以及该第一类型孔洞的特征数据能够体现该第一类型孔洞的特征，因此，通过第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据获得的第二毛管压力曲线，能够同时体现岩石样本中的该第一类型孔洞和该第二类型孔洞的特征，因此，本方案能够提升通过毛管压力曲线评价岩石孔喉结构的准确性。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。

附图说明

此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本申请的实施例，并与说明书一起用于解释本申请的原理。

图1是根据一示例性实施例示出的一种孔喉特征参数的获取系统的系统构成图；

图2是根据一示例性实施例示出的一种孔喉特征参数的获取方法的流程示意图；

图3是图2所示实施例涉及的碳酸盐岩气藏开发地质特征描述SY/T6110-2002的示意图；

图4是根据一示例性实施例示出的一种孔喉特征参数的获取方法的流程示意图；

图5是图4所示实施例涉及的一种微CT扫描得到的三维数字岩心的示意图；

图6是图4所示实施例涉及的一种高压压汞实验得到的毛管压力曲线的示意图；

图7是图4所示实施例涉及的一种修正后得到的毛管压力曲线的示意图；

图8是根据一示例性实施例示出的孔喉特征参数的获取预测装置的结构方框图；

图9是根据一示例性实施例示出的一种计算机设备的结构框图。

具体实施方式

这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。

应当理解的是，在本文中提及的“若干个”是指一个或者多个，“多个”是指两个或两个以上。“和/或”，描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。

在对本申请所示的各个实施例进行说明之前，首先对本申请涉及到的几个概念进行介绍：

1)孔隙度(porosity)

孔隙度是岩石的孔隙体积与其外表体积之比，以百分数表示。如果孔隙体积由岩石体积和颗粒体积之差获得，其包括连通的和不连通的全部孔隙，则称绝对孔隙度；反之，如果孔隙体积只限于连通的，即有效孔隙，则称有效孔隙度。岩石中流动液体在渗流过程中，所涉及到的孔隙称流动孔隙，其孔隙度称流动孔隙度。常用的是有效孔隙度，它是油藏开发计算以及评价储层性质必不可少的岩石物性参数。对于带有裂隙的岩石，将该岩石的裂缝孔隙度和粒间孔隙的孔隙度相加，即可以得到该岩石的总孔隙度。

2)渗透率(permeability)

渗透率是衡量多孔介质允许流体通过能力的一种量度。根据达西定律，粘度为μ的液体在△p压差作用下通过截面积为A、长度为L的多孔介质，其通过的流量Q与A、△p成正比，与μ、L成反比，即Q＝KA△p/μL，K为比例常数，称渗透率，K值大，多孔介质允许流体通过的能力也大，反之则小。渗透率是多孔介质的自身性质，与所通过的流体性质无关。渗透率的单位通常采用μm

3)毛管压力曲线(capillary curve)

岩心毛管压力与孔隙中润湿相流体饱和度的关系曲线被称为毛管压力曲线。不同的毛管压力曲线的形态特征反映了岩石中不同的孔隙结构特征。

4)微CT扫描(micro-CT scanning)

微CT扫描作为一种无损检测物体内部结构的技术，是当前建立三维数字岩心最直接和最准确的方法，其原理是根据岩石中不同密度的成分对X射线吸收系数不同，来区分岩石的孔隙和骨架。

5)归一化(normalize)

归一化又称规格化；归一化处理是指将数据调整到规定的标准形式的过程。

6)高压压汞法(Mercury Intrusiong Capillary Porosimetry，MICP)

高压压汞法是依靠外加压力使汞克服表面张力进入岩石的孔隙，来测量孔隙的孔径和分布的方法。

7)麻皮效应(pockmarks-effect)

汞在岩石表面的粗糙坑凹处、拐角处和棱线上形成一些不贴和的汞蒸气空腔，这种现象称作麻皮效应。

8)核磁共振(Nuclear Magnetic Resonance，NMR)

在静磁场中，核磁共振是指在某些特征频率的射频场吸收能量的物理现象。电磁学将应用于核磁矩的共振方法称为核磁共振。质子犹如小磁铁，具有磁矩。各种原子核由于结构不同具有不同的磁矩，这种磁矩被称为核磁矩。质子在磁场中所产生的旋进有一定的频率，其大小由质子的磁矩及所在磁场的强度而定。如在垂直原来固定磁场方向另加一可变磁场，并使其频率正好等于质子的旋进频率，便会产生了共振现象，这就是所谓的核磁共振。

9)排驱压力(displacement pressure)

排驱压力指某一岩样中的润湿相流体被非润湿相流体开始排替所需的最低压力。

图1是根据一示例性实施例示出的一种孔喉特征参数的获取系统的系统构成图。该系统包括：特征信息采集设备110、终端120。

特征信息采集设备110是具有采集岩石特征信息能力的特征信息采集设备，比如，特征信息采集设备可以是共聚焦激光扫描显微镜(Laser Scanning ConfocalMicroscopy，LSCM)、X射线计算机层析成像仪(X-Computerized Tomography，X-CT)、聚集离子束-扫描电镜(Focused Ion Beam-Scanning Electron Microscope，FIB-SEM)等高精度仪器等。

终端120可以是具有处理特征信息功能的终端设备，比如，终端可以是手机、平板电脑、智能眼镜、智能手表、膝上型便携计算机和台式计算机等等。终端120的数量不做限定。

其中，特征信息采集设备110与终端120之间通过通信网络相连，能够进行信息的传递。

或者，特征信息采集设备110和终端120之间也可以通过互联网进行信息的传递，比如，特征信息采集设备110和终端120可以分别安装具有信息传递功能的应用程序，并以该应用程序对应在互联网中的服务器为中继进行特征信息的传递。

特征信息采集设备110与终端120之间通过通信网络相连。可选的，通信网络是有线网络或无线网络。

在本申请实施例中，特征信息采集设备110可以将处理后的特征信息数据发送给终端120，由终端120对特征信息做进一步处理。

可选的，上述的无线网络或有线网络使用标准通信技术和/或协议。网络通常为因特网、但也可以是任何网络，包括但不限于局域网(Local Area Network，LAN)、城域网(Metropolitan Area Network，MAN)、广域网(Wide Area Network，WAN)、移动、有线或者无线网络、专用网络或者虚拟专用网络的任何组合)。在一些实施例中，使用包括超文本标记语言(Hyper Text Mark-up Language，HTML)、可扩展标记语言(Extensible MarkupLanguage，XML)等的技术和/或格式来代表通过网络交换的数据。此外还可以使用诸如安全套接字层(Secure Socket Layer，SSL)、传输层安全(Transport Layer Security，TLS)、虚拟专用网络(Virtual Private Network，VPN)、网际协议安全(Internet ProtocolSecurity，IPsec)等常规加密技术来加密所有或者一些链路。在另一些实施例中，还可以使用定制和/或专用数据通信技术取代或者补充上述数据通信技术。

请参考图2，其是根据一示例性实施例示出的一种孔喉特征参数的获取方法的流程示意图。该方法可以由计算机设备执行，该计算机设备可以是上述图1所示实施例中的终端，或者，该计算机设备可以包括上述图1所示实施例中的特征信息采集设备和终端。如图2所示，该孔喉特征参数的获取方法可以包括如下步骤：

步骤201，获取岩石样本的特征信息，该特征信息包括该岩石样本的长度、直径、孔隙度和渗透率；

步骤202，对该岩石样本进行特征提取处理，获得该岩石样本的第一类型孔洞的特征数据，该第一类型孔洞包括表面溶洞和尺寸大于尺寸阈值的孔；

步骤203，根据该特征信息和该第一类型孔洞的特征数据，获取该第一类型孔洞的进汞饱和度；

步骤204，获取该岩石样本对应的第一毛管压力曲线，该第一毛管压力曲线用于指示岩石样本中第二类型孔洞的体积占比，该第二类型孔洞的体积占比是该第二类型孔洞的体积占该岩石样本的总孔隙体积的比例，该第二类型孔洞是尺寸小于尺寸阈值的孔隙；

步骤205，根据该第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据，获取第二毛管压力曲线；

步骤206，根据该第二毛管压力曲线获取该岩石样本的孔喉特征参数，该孔喉特征参数包括最大孔隙半径、平均孔隙半径、中值、孔喉半径均值、分选系数、歪度、峰态、均质系数。

其中，岩石样本可以是碳酸盐岩样本，碳酸盐岩除了有微米级孔喉和裂缝外，毫米级、厘米级溶洞也较为发育，而岩石储集空间具有显著的多尺度性，这也导致了不同类型碳酸盐岩孔喉分布很难被完整的获取；表面溶洞可以是岩石表面未饱和水的溶洞，表面溶洞可以是碳酸盐岩气藏开发地质特征描述SY/T 6110-2002上描述的溶洞，例如该溶洞可以是溶洞半径在1厘米到5厘米之间的溶洞。

上述尺寸阈值可以是碳酸盐岩气藏开发地质特征描述SY/T 6110-2002上孔的尺寸阈值，该第一类型孔洞可以是该碳酸盐岩气藏开发地质特征描述SY/T 6110-2002上的大孔和溶洞，该第二类型孔洞可以是该碳酸盐岩气藏开发地质特征描述SY/T 6110-2002上的中孔、小孔和微孔。

上述尺寸大于尺寸阈值的孔(即上述第一类型孔洞)可以是孔隙半径大于0.4毫米的孔隙(也可以称为大孔)，上述尺寸小于尺寸阈值的孔隙(即上述第二类型孔洞)可以是孔隙半径小于0.4毫米的孔隙(也可称为中小孔)。

上述第一类型孔洞的特征数据可以是用来描述表面溶洞和尺寸大于尺寸阈值的孔的结构特征的数据，例如，该特征数据可以是第一类型孔洞的孔洞数、孔洞密度，孔洞长度、孔洞宽度，也可以是通过这些物理参数获得的第一类型孔洞的体积、面积百分比、孔隙度和分形维数。

请参考图3，其示出了本申请实施例涉及的碳酸盐岩气藏开发地质特征描述SY/T6110-2002的示意图，该图中的数值单位为毫米。

本申请实施例中提出的该第二毛管压力曲线可以是通过该第一毛管压力曲线和能够指示该岩石样本上第一类型孔洞分布的参数得到的，而该第一毛管压力曲线又能指示该岩石样本中第二类型孔洞的分布情况，因此，通过该第二毛管压力曲线能够获得该岩石样本中完整的微米至厘米级孔喉结构特征参数；该第二毛管压力曲线可以用于指示该岩石样本中所有孔隙的分布情况。

综上所述，在本申请实施例所示的方案中，通过岩石样本的特征信息和该岩石样本的第一类型孔洞的特征数据，获得该第一类型孔洞的进汞饱和度；获取指示岩石样本中第二类型孔洞的体积占比的第一毛管压力曲线；通过该第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据，获得第二毛管压力曲线；再根据该第二毛管压力曲线获取该岩石样本的孔喉特征参数。由于第一毛管压力曲线能够体现该第二类型孔洞的特征，而进汞饱和度以及该第一类型孔洞的特征数据能够体现该第一类型孔洞的特征，因此，通过第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据获得的第二毛管压力曲线，能够同时体现岩石样本中的该第一类型孔洞和该第二类型孔洞的特征，因此，本方案能够提升通过毛管压力曲线评价岩石孔喉结构的准确性。

请参考图4，其是根据一示例性实施例提供的一种孔喉特征参数的获取方法的方法流程图。该方法可以由计算机设备执行，该计算机设备可以是上述图1所示实施例中的终端，或者，该计算机设备可以包括上述图1所示实施例中的特征信息采集设备和终端。如图4所示，该孔喉特征参数的获取方法可以包括如下步骤：

步骤401，获取岩石样本的特征信息，该特征信息包括该岩石样本的长度、直径、孔隙度和渗透率。

在一种可选的实现方式中，该岩石样本是碳酸盐岩岩样，该碳酸盐岩岩样含有微米至厘米级别的孔喉。

在一种可选的实现方式中，该岩石样本是通过切割工具切割得到的岩石柱塞样本。

示例性的，该岩石样本可以被切割成直径为2.5厘米、长为2厘米的岩石柱塞样本。

在一种可选的实现方式中，对岩石样本进行核磁共振处理，测量岩石孔隙中含氢流体的核磁弛豫信号，通过不同大小孔隙中氢离子的响应，最终得出岩石特征信息。

在一种可选的实现方式中，使用游标卡尺和气测渗透率仪来获得相关的物理参数，根据该物理参数计算得到岩石特征信息。

步骤402，对该岩石样本进行特征提取处理，获得该岩石样本的第一类型孔洞的特征数据，该第一类型孔洞包括表面溶洞和尺寸大于尺寸阈值的孔。

在一种可选的实现方式中，通过采样电子光束对该岩石样本的岩心进行抛光，再用高像素的拍照仪器对该岩心进行拍照，得到高分辨率的图像；每当岩心进行一次抛光，便用高像素的拍照设备拍照，于是得到一张二维图像；经过多次抛光后，得到若干数量的二维图像；将这若干张连续的二维图像进行分离组合，得到该岩石样本的三维数字岩心；通过该三维数字岩心得到该岩石样本的第一类型孔洞的特征信息。

在一种可选的实现方式中，先将岩石样本切片，得到岩石薄片；再利用电子显微镜等装置获取岩石薄片的特征信息；通过数据处理软件，将连续的岩石薄片上的特征信息组合，便得到该岩石样本的三维数字岩心；再通过该岩石样本的三维数字岩心获得该岩石样本的第一类型孔洞的特征数据。

在一种可选的实现方式中，对该岩石样本进行微CT扫描，再将该微CT扫描得到的数据传输给数据处理软件；由该数据处理软件处理该数据，得到三维数字岩心，然后根据该三维数字岩心获取岩石样本的第一类型孔洞的特征数据。

例如，该微CT扫描通常包括三个步骤：(1)对岩样进行预处理后开展CT实验获得投影数据；(2)通过三维图像重建方法，利用该投影数据重建岩心灰度图像；(3)采用图像二值分割方法分离该灰度图像中的孔隙空间和岩石骨架，从而构建三维数字岩心，并根据该三维数字岩心获取岩石样本的特征数据。

其中，对岩石样本进行微CT扫描时，可以使用微米级别的CT扫描仪；而数据处理软件可以是具备三维可视化能力的软件；数字岩心是岩心资料经过数据化处理得到的；通过数字岩心，可以保证岩心原始图像信息资料的完整性和真实性。

请参考图5其示出了本申请实施例涉及的微CT扫描得到的三维数字岩心的示意图，其中，上层图片和下层图片在位置上是一一对应关系；其中，图501是孔隙型岩石样本的三维数字岩心模型图像、图502是裂缝孔隙型岩石样本的三维数字岩心模型图像、图503是孔洞型岩石样本的三维数字岩心模型图像、图504是裂缝孔洞型岩石样本的三维数字岩心模型图像、图505是孔隙型岩石样本的岩心灰度图像、图506是裂缝孔隙型岩石样本的岩心灰度图像、图507是孔洞型岩石样本的岩心灰度图像、图508是裂缝孔洞型岩石样本的岩心灰度图像。

在一种可选的实现方式中，该特征数据是该第一类型孔洞的形态特征参数，该特征参数包括孔洞的数量、孔洞的最大卡规直径、孔洞的面积、孔洞的圆度、孔洞的粗糙度等等。

请参考表一，其示出了本申请实施例涉及的不同类型孔隙的岩石样本的特征参数的参数表。

表一

步骤403，根据该特征信息，获得该岩石样本的总孔隙体积。

在一种可选的实现方式中，该特征信息包括岩石样本的体积以及岩石样本的孔隙度，根据该岩石样本的体积和该孔隙度，计算得到岩石样本的总孔隙体积。

示例性的，该岩石样本的孔隙度为10％，该岩石样本的体积为10立方米，则该岩石样本的总孔隙体积为1立方米。

在一种可选的实现方式中，该特征信息包括岩石样本的体积和密度，根据完整的岩石样本的质量与含有孔隙的岩石样本的质量的差值，得到该岩石样本的总孔隙体积。

步骤404，根据该第一类型孔洞的特征数据和该岩石样本的总孔隙体积，获得进汞饱和度。

在一种可选的实现方式中，可以根据该第一类型孔洞的特征数据得到第一类型孔洞的体积；将该第一类型孔洞的体积，除以该总孔隙体积，作为该进汞饱和度。

示例性的，该第一类型孔洞的体积为1立方米，总孔隙体积为5立方米，则该进汞饱和度为20％。

步骤405，通过高压压汞实验获得该岩石样本的毛管压力曲线。

在一种可选的实现方式中，通过高压压汞的方法，对该岩石样本进行处理；通过外加压力使汞克服表面张力进入孔隙，从而获得该岩石样本的原始毛管压力曲线。

请参考图6，其示出了本申请实施例涉及的高压压汞实验得到的毛管压力曲线的示意图。

步骤406，对该毛管压力曲线进行麻皮效应校正，获得该第一毛管压力曲线，该第一毛管压力曲线用于指示岩石样本中第二类型孔洞的体积占比，该第二类型孔洞的体积占比是该第二类型孔洞的体积占该岩石样本中总孔隙体积的比例，该第二类型孔洞是尺寸小于尺寸阈值的孔隙。

在一种可选的实现方式中，依据行业标准SY/T 5346-2005《岩石毛管压力曲线的测定》对毛管压力曲线进行麻皮效应校正，排除麻皮效应数值。

在一种可能的实现方式中，通过作图法求取排驱压力；在单对数坐标中沿毛管压力曲线的平坦部分做切线与纵轴相交的压力点即排驱压力点，此点对应的汞饱和度即为麻皮效应数值；再通过麻皮效应值对毛管压力曲线进行修正得到第一毛管压力曲线。

步骤407，根据该第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据，获取第二毛管压力曲线。

在一种可选的方式中，通过该进汞饱和度和该第一类型孔洞的特征数据，修正该第一毛管压力曲线，获取该第二毛管压力曲线。

在一种可选的实现方式中，该第一类型孔洞的特征数据包括与该第一类型孔洞的体积对应的该第一类型孔洞的直径。

示例性的，该第一类型孔洞的直径为1米；通过微CT扫描得到该第一类型孔洞的体积为5立方米；该直径与该体积存在对应关系。

在一种可选的实现方式中，通过该第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的直径，计算得到该第二毛管压力曲线。

在一种可选的实现方式中，根据该进汞饱和度和该第一类型孔洞的直径，对该第一毛管压力曲线进行修正，得到该第二毛管压力曲线。

请参考图7，其示出了本申请实施例涉及的修正后得到的毛管压力曲线的示意图。

步骤408，根据该第二毛管压力曲线获取该岩石样本的孔喉特征参数，该孔喉特征参数包括最大孔隙半径、平均孔隙半径、中值、孔喉半径均值、分选系数、歪度、峰态、均质系数。

在一种可选的实现方式中，对该第二毛管压力曲线中累积的进汞饱和度进行归一化处理，再定量计算该岩石样本的孔喉特征参数。

在一种可选的实现方式中，该孔喉特征参数用于指示该岩石样本的渗流能力和该岩石样本的湿润性。

在一种可选的实现方式中，该最大孔隙半径的计算公式表示为：

其中，R

在一种可选的实现方式中，该平均孔隙半径的计算公式表示为：

其中，

在一种可选的实现方式中，该中值的计算公式表示为：

其中，R

在一种可选的实现方式中，该孔喉半径均值的计算公式表示为：

其中，D是孔喉直径；

在一种可选的实现方式中，该分选系数的计算公式表示为：

其中，Sp是分选系数；

在一种可选的实现方式中，该歪度的计算公式表示为：

其中，S

在一种可选的实现方式中，该峰态的计算公式表示为：

其中，K

在一种可选的实现方式中，该均质系数的计算公式表示为：

其中，α是均质系数；

请参考表二，其示出了本申请实施例涉及的该种数字岩心处理方法校正前后的岩石孔喉特征参数的参数表。

表二

综上所述，在本申请实施例所示的方案中，通过岩石样本的特征信息和该岩石样本的第一类型孔洞的特征数据，获得该第一类型孔洞的进汞饱和度；获取指示岩石样本中第二类型孔洞的体积占比的第一毛管压力曲线；通过该第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据，获得第二毛管压力曲线；再根据该第二毛管压力曲线获取该岩石样本的孔喉特征参数。由于第一毛管压力曲线能够体现该第二类型孔洞的特征，而进汞饱和度以及该第一类型孔洞的特征数据能够体现该第一类型孔洞的特征，因此，通过第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据获得的第二毛管压力曲线，能够同时体现岩石样本中的该第一类型孔洞和该第二类型孔洞的特征，因此，本方案能够提升通过毛管压力曲线评价岩石孔喉结构的准确性。

下述为本公开装置实施例，可以用于执行本公开方法实施例。对于本公开装置实施例中未披露的细节，请参照本公开方法实施例。

图8是根据一示例性实施例示出的一种孔喉特征参数的获取方法装置的结构方框图。该孔喉特征参数的获取方法装置可以实现图2所示实施例提供的方法中的全部或者部分步骤。该孔喉特征参数的获取方法装置可以包括：

特征信息获取模块801，用于获取岩石样本的特征信息，所述特征信息包括所述岩石样本的长度、直径、孔隙度和渗透率；

特征数据获取模块802，用于对所述岩石样本进行特征提取处理，获得所述岩石样本的第一类型孔洞的特征数据，所述第一类型孔洞包括表面溶洞和尺寸大于尺寸阈值的孔；

进汞饱和度获取模块803，用于根据所述特征信息和所述第一类型孔洞的特征数据，获取所述第一类型孔洞的进汞饱和度；

第一毛管压力曲线获取模块804，用于获取所述岩石样本对应的第一毛管压力曲线，所述第一毛管压力曲线用于指示岩石样本中第二类型孔洞的体积占比，所述第二类型孔洞的体积占比是所述第二类型孔洞的体积占所述岩石样本的总孔隙体积的比例，所述第二类型孔洞是尺寸小于尺寸阈值的孔隙；

第二毛管压力曲线获取模块805，用于根据所述第一毛管压力曲线、所述进汞饱和度以及所述第一类型孔洞的特征信息，获取第二毛管压力曲线；

孔喉特征参数获取模块806，用于根据所述第二毛管压力曲线获取所述岩石样本的孔喉特征参数，所述孔喉特征参数包括最大孔隙半径、平均孔隙半径、中值、孔喉半径均值、分选系数、歪度、峰态、均质系数。

可选的，所述第二毛管压力曲线获取模块805，用于通过所述进汞饱和度和所述第一类型孔洞的特征数据，修正所述第一毛管压力曲线，获取所述第二毛管压力曲线。

可选的，所述进汞饱和度获取模块803，包括：

总孔隙体积获取子模块，用于根据所述特征信息，获得所述岩石样本的总孔隙体积；

数据处理子模块，用于根据所述第一类型孔洞的特征数据和所述岩石样本的总孔隙体积，获得所述进汞饱和度。

可选的，所述进汞饱和度处理子模块，包括：

体积获取单元，用于根据所述第一类型孔洞的特征数据，获得所述第一类型孔洞的体积；

进汞饱和度计算单元，用于将所述第一类型孔洞的体积，除以所述岩石样本的总孔隙体积，作为所述进汞饱和度。

可选的，所述特征数据获取模块802，包括：

三维数字岩心获取子模块，用于对所述岩石样本进行微CT扫描，获得所述岩石样本的三维数字岩心；

三维数字岩心处理子模块，用于根据所述三维数字岩心，获得所述岩石样本的第一类型孔洞的特征数据，所述第一类型孔洞的特征数据包括所述第一类型孔洞的直径和体积。

可选的，所述孔喉特征参数获取模块806，用于对所述第二毛管压力曲线中累积的所述进汞饱和度进行归一化处理，获得所述岩石样本的孔喉特征参数。

可选的，所述第一毛管压力曲线获取模块804，包括：

高压压汞处理子模块，通过高压压汞实验获得所述岩石样本的原始毛管压力曲线；

麻皮效应校正子模块，用于对所述毛管压力曲线进行麻皮效应校正，获得所述第一毛管压力曲线。

综上所述，在本申请实施例所示的方案中，通过岩石样本的特征信息和该岩石样本的第一类型孔洞的特征数据，获得该第一类型孔洞的进汞饱和度；获取指示岩石样本中第二类型孔洞的体积占比的第一毛管压力曲线；通过该第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据，获得第二毛管压力曲线；再根据该第二毛管压力曲线获取该岩石样本的孔喉特征参数。由于第一毛管压力曲线能够体现该第二类型孔洞的特征，而进汞饱和度以及该第一类型孔洞的特征数据能够体现该第一类型孔洞的特征，因此，通过第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据获得的第二毛管压力曲线，能够同时体现岩石样本中的该第一类型孔洞和该第二类型孔洞的特征，因此，本方案能够提升通过毛管压力曲线评价岩石孔喉结构的准确性。

需要说明的一点是，上述实施例提供的装置在实现其功能时，仅以上述各个功能模块的划分进行举例说明，实际应用中，可以根据实际需要而将上述功能分配由不同的功能模块完成，即将设备的内容结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。

关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。

图9示出了本申请一个示例性实施例示出的计算机设备900的结构框图。该计算机设备900可以是图1所示系统中的用户终端或者服务器。

通常，计算机设备900包括有：处理器901和存储器902。

处理器901可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器901可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field－Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器901也可以包括主处理器和协处理器。在一些实施例中，处理器901可以在集成有GPU(Graphics Processing Unit，图像处理器)。

存储器902可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器902还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器902中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器901所执行以实现本申请上述方法实施例中的全部或者部分步骤。

在一些实施例中，计算机设备实现为用户终端时，该计算机设备900还可选包括有：外围设备接口903和至少一个外围设备。处理器901、存储器902和外围设备接口903之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口903相连。可选的，外围设备包括：射频电路904、显示屏905、图像采集组件906、音频电路907、定位组件908和电源909中的至少一种。

外围设备接口903可被用于将I/O(Input/Output，输入/输出)相关的至少一个外围设备连接到处理器901和存储器902。

射频电路904用于接收和发射RF(Radio Frequency，射频)信号，也称电磁信号。可选地，射频电路904包括：天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路904可以通过至少一种无线通信协议来与其它计算机设备进行通信。在一些实施例中，射频电路904还可以包括NFC(NearField Communication，近距离无线通信)有关的电路，本申请对此不加以限定。

显示屏905用于显示UI(UserInterface，用户界面)。当显示屏905是触摸显示屏时，显示屏905还具有采集在显示屏905的表面或表面上方的触摸信号的能力。

图像采集组件906用于采集图像或视频。在一些实施例中，图像采集组件906还可以包括闪光灯。

音频电路907可以包括麦克风和扬声器。在一些实施例中，音频电路907还可以包括耳机插孔。

定位组件908用于定位计算机设备900的当前地理位置，以实现导航或LBS(Location Based Service，基于位置的服务)。

电源909于为计算机设备900中的各个组件进行供电。

在一些实施例中，计算机设备900还包括有一个或多个传感器910。该一个或多个传感器910包括但不限于：加速度传感器911、陀螺仪传感器912、压力传感器913、指纹传感器914、光学传感器915以及接近传感器916。

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

在示例性实施例中，还提供了一种包括指令的非临时性计算机可读存储介质，例如包括计算机程序(指令)的存储器，上述程序(指令)可由计算机设备的处理器执行以完成本申请各个实施例所示的方法中，由服务器或者用户终端执行的方法。例如，所述非临时性计算机可读存储介质可以是只读存储器(Read-Only Memory，ROM)、随机存取存储器(Random Access Memory，RAM)、只读光盘(Compact Disc Read-Only Memory，CD-ROM)、磁带、软盘和光数据存储设备等。

本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本申请的其它实施方案。本申请旨在涵盖本申请的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本申请的一般性原理并包括本申请未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本申请的真正范围和精神由权利要求指出。

应当理解的是，本申请并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本申请的范围仅由所附的权利要求来限制。