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

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

Info

Publication number
CN112326523A
CN112326523A CN202010984074.9A CN202010984074A CN112326523A CN 112326523 A CN112326523 A CN 112326523A CN 202010984074 A CN202010984074 A CN 202010984074A CN 112326523 A CN112326523 A CN 112326523A
Authority
CN
China
Prior art keywords
type
rock sample
hole
capillary pressure
pressure curve
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010984074.9A
Other languages
English (en)
Other versions
CN112326523B (zh
Inventor
李骞
李滔
赵梓寒
张飞
戚涛
李秋
姜艺
张寰一
吴婷婷
江俊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010984074.9A priority Critical patent/CN112326523B/zh
Publication of CN112326523A publication Critical patent/CN112326523A/zh
Application granted granted Critical
Publication of CN112326523B publication Critical patent/CN112326523B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering

Landscapes

  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Immunology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Dispersion Chemistry (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pulmonology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Graphics (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

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

Description

孔喉特征参数的获取方法、装置、计算机设备及存储介质
技术领域
本申请涉及油气勘探技术领域，特别涉及一种孔喉特征参数的获取方法、装置、计算机设备及存储介质。
背景技术
毛管压力曲线可以反映岩石中的不同孔隙的分布情况，是评价岩石孔喉结构和储层岩石渗流能力的重要依据。
在相关的技术中，由于核磁共振弛豫时间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值大，多孔介质允许流体通过的能力也大，反之则小。渗透率是多孔介质的自身性质，与所通过的流体性质无关。渗透率的单位通常采用μm2。
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是裂缝孔洞型岩石样本的岩心灰度图像。
在一种可选的实现方式中，该特征数据是该第一类型孔洞的形态特征参数，该特征参数包括孔洞的数量、孔洞的最大卡规直径、孔洞的面积、孔洞的圆度、孔洞的粗糙度等等。
请参考表一，其示出了本申请实施例涉及的不同类型孔隙的岩石样本的特征参数的参数表。
Figure BDA0002688615440000111
表一
步骤403，根据该特征信息，获得该岩石样本的总孔隙体积。
在一种可选的实现方式中，该特征信息包括岩石样本的体积以及岩石样本的孔隙度，根据该岩石样本的体积和该孔隙度，计算得到岩石样本的总孔隙体积。
示例性的，该岩石样本的孔隙度为10％，该岩石样本的体积为10立方米，则该岩石样本的总孔隙体积为1立方米。
在一种可选的实现方式中，该特征信息包括岩石样本的体积和密度，根据完整的岩石样本的质量与含有孔隙的岩石样本的质量的差值，得到该岩石样本的总孔隙体积。
步骤404，根据该第一类型孔洞的特征数据和该岩石样本的总孔隙体积，获得进汞饱和度。
在一种可选的实现方式中，可以根据该第一类型孔洞的特征数据得到第一类型孔洞的体积；将该第一类型孔洞的体积，除以该总孔隙体积，作为该进汞饱和度。
示例性的，该第一类型孔洞的体积为1立方米，总孔隙体积为5立方米，则该进汞饱和度为20％。
步骤405，通过高压压汞实验获得该岩石样本的毛管压力曲线。
在一种可选的实现方式中，通过高压压汞的方法，对该岩石样本进行处理；通过外加压力使汞克服表面张力进入孔隙，从而获得该岩石样本的原始毛管压力曲线。
请参考图6，其示出了本申请实施例涉及的高压压汞实验得到的毛管压力曲线的示意图。
步骤406，对该毛管压力曲线进行麻皮效应校正，获得该第一毛管压力曲线，该第一毛管压力曲线用于指示岩石样本中第二类型孔洞的体积占比，该第二类型孔洞的体积占比是该第二类型孔洞的体积占该岩石样本中总孔隙体积的比例，该第二类型孔洞是尺寸小于尺寸阈值的孔隙。
在一种可选的实现方式中，依据行业标准SY/T 5346-2005《岩石毛管压力曲线的测定》对毛管压力曲线进行麻皮效应校正，排除麻皮效应数值。
在一种可能的实现方式中，通过作图法求取排驱压力；在单对数坐标中沿毛管压力曲线的平坦部分做切线与纵轴相交的压力点即排驱压力点，此点对应的汞饱和度即为麻皮效应数值；再通过麻皮效应值对毛管压力曲线进行修正得到第一毛管压力曲线。
步骤407，根据该第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据，获取第二毛管压力曲线。
在一种可选的方式中，通过该进汞饱和度和该第一类型孔洞的特征数据，修正该第一毛管压力曲线，获取该第二毛管压力曲线。
在一种可选的实现方式中，该第一类型孔洞的特征数据包括与该第一类型孔洞的体积对应的该第一类型孔洞的直径。
示例性的，该第一类型孔洞的直径为1米；通过微CT扫描得到该第一类型孔洞的体积为5立方米；该直径与该体积存在对应关系。
在一种可选的实现方式中，通过该第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的直径，计算得到该第二毛管压力曲线。
在一种可选的实现方式中，根据该进汞饱和度和该第一类型孔洞的直径，对该第一毛管压力曲线进行修正，得到该第二毛管压力曲线。
请参考图7，其示出了本申请实施例涉及的修正后得到的毛管压力曲线的示意图。
步骤408，根据该第二毛管压力曲线获取该岩石样本的孔喉特征参数，该孔喉特征参数包括最大孔隙半径、平均孔隙半径、中值、孔喉半径均值、分选系数、歪度、峰态、均质系数。
在一种可选的实现方式中，对该第二毛管压力曲线中累积的进汞饱和度进行归一化处理，再定量计算该岩石样本的孔喉特征参数。
在一种可选的实现方式中，该孔喉特征参数用于指示该岩石样本的渗流能力和该岩石样本的湿润性。
在一种可选的实现方式中，该最大孔隙半径的计算公式表示为：
Figure BDA0002688615440000131
其中，Rmax是最大孔隙半径；Pcd是排驱压力。
在一种可选的实现方式中，该平均孔隙半径的计算公式表示为：
Figure BDA0002688615440000132
其中，
Figure BDA0002688615440000133
是平均孔隙半径；Ri是孔隙半径；SHgi是不同孔隙对应的进汞饱和度。
在一种可选的实现方式中，该中值的计算公式表示为：
Figure BDA0002688615440000141
其中，R50是汞饱和度中值；孔隙半径；Pc50是饱和度中值压力。
在一种可选的实现方式中，该孔喉半径均值的计算公式表示为：
Figure BDA0002688615440000142
Figure BDA0002688615440000143
其中，D是孔喉直径；
Figure BDA0002688615440000144
是岩石样本对应的孔隙度；Rm是孔喉半径均值；
Figure BDA0002688615440000145
是汞饱和度为n％时岩石样本的孔隙度。
在一种可选的实现方式中，该分选系数的计算公式表示为：
Figure BDA0002688615440000146
其中，Sp是分选系数；
Figure BDA0002688615440000147
是汞饱和度为n％时岩石样本的孔隙度。
在一种可选的实现方式中，该歪度的计算公式表示为：
Figure BDA0002688615440000148
其中，Skp是歪度；
Figure BDA0002688615440000149
是汞饱和度为n％时岩石样本的孔隙度。
在一种可选的实现方式中，该峰态的计算公式表示为：
Figure BDA00026886154400001410
其中，Kp是峰态；
Figure BDA00026886154400001411
是汞饱和度为n％时岩石样本的孔隙度。
在一种可选的实现方式中，该均质系数的计算公式表示为：
Figure BDA00026886154400001412
其中，α是均质系数；
Figure BDA00026886154400001413
是平均孔隙半径；Rmax是最大孔隙半径。
请参考表二，其示出了本申请实施例涉及的该种数字岩心处理方法校正前后的岩石孔喉特征参数的参数表。
Figure BDA00026886154400001414
Figure BDA0002688615440000151
表二
综上所述，在本申请实施例所示的方案中，通过岩石样本的特征信息和该岩石样本的第一类型孔洞的特征数据，获得该第一类型孔洞的进汞饱和度；获取指示岩石样本中第二类型孔洞的体积占比的第一毛管压力曲线；通过该第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据，获得第二毛管压力曲线；再根据该第二毛管压力曲线获取该岩石样本的孔喉特征参数。由于第一毛管压力曲线能够体现该第二类型孔洞的特征，而进汞饱和度以及该第一类型孔洞的特征数据能够体现该第一类型孔洞的特征，因此，通过第一毛管压力曲线、该进汞饱和度以及该第一类型孔洞的特征数据获得的第二毛管压力曲线，能够同时体现岩石样本中的该第一类型孔洞和该第二类型孔洞的特征，因此，本方案能够提升通过毛管压力曲线评价岩石孔喉结构的准确性。
下述为本公开装置实施例，可以用于执行本公开方法实施例。对于本公开装置实施例中未披露的细节，请参照本公开方法实施例。
图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)、磁带、软盘和光数据存储设备等。
本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本申请的其它实施方案。本申请旨在涵盖本申请的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本申请的一般性原理并包括本申请未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本申请的真正范围和精神由权利要求指出。
应当理解的是，本申请并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本申请的范围仅由所附的权利要求来限制。

Claims (10)

1.一种孔喉特征参数的获取方法，其特征在于，所述方法包括：
获取岩石样本的特征信息，所述特征信息包括所述岩石样本的长度、直径、孔隙度和渗透率；
对所述岩石样本进行特征提取处理，获得所述岩石样本的第一类型孔洞的特征数据，所述第一类型孔洞包括表面溶洞和尺寸大于尺寸阈值的孔；
根据所述特征信息和所述第一类型孔洞的特征数据，获取所述第一类型孔洞的进汞饱和度；
获取所述岩石样本对应的第一毛管压力曲线，所述第一毛管压力曲线用于指示岩石样本中第二类型孔洞的体积占比，所述第二类型孔洞的体积占比是所述第二类型孔洞的体积占所述岩石样本的总孔隙体积的比例，所述第二类型孔洞是尺寸小于尺寸阈值的孔隙；
根据所述第一毛管压力曲线、所述进汞饱和度以及所述第一类型孔洞的特征数据，获取第二毛管压力曲线；
根据所述第二毛管压力曲线获取所述岩石样本的孔喉特征参数，所述孔喉特征参数包括最大孔隙半径、平均孔隙半径、孔喉半径均值、分选系数、歪度、峰态、均质系数。
2.根据权利要求1所述的方法，其特征在于，所述根据所述第一毛管压力曲线、所述进汞饱和度以及所述第一类型孔洞的特征数据，获取第二毛管压力曲线，包括：
通过所述进汞饱和度和所述第一类型孔洞的特征数据，修正所述第一毛管压力曲线，获取所述第二毛管压力曲线。
3.根据权利要求1所述的方法，其特征在于，所述根据所述特征信息和所述第一类型孔洞的特征数据，获取所述第一类型孔洞的进汞饱和度，包括：
根据所述特征信息，获得所述岩石样本的总孔隙体积；
根据所述第一类型孔洞的特征数据和所述岩石样本的总孔隙体积，获得所述进汞饱和度。
4.根据权利要求3所述的方法，其特征在于，所述根据所述第一类型孔洞的特征数据和所述岩石样本的总孔隙体积，获得所述进汞饱和度，包括：
根据所述第一类型孔洞的特征数据，获得所述第一类型孔洞的体积；
将所述第一类型孔洞的体积，除以所述岩石样本的总孔隙体积，作为所述进汞饱和度。
5.根据权力要求1所述的方法，其特征在于，所述对所述岩石样本进行特征提取处理，获得所述岩石样本的第一类型孔洞的特征数据，包括：
对所述岩石样本进行微CT扫描，获得所述岩石样本的三维数字岩心；
根据所述三维数字岩心，获得所述岩石样本的第一类型孔洞的特征数据，所述第一类型孔洞的特征数据包括所述第一类型孔洞的直径和体积。
6.根据权利要求1所述的方法，其特征在于，所述根据所述第二毛管压力曲线获取所述岩石样本的孔喉特征参数，包括：
对所述第二毛管压力曲线中累积的所述进汞饱和度进行归一化处理，获得所述岩石样本的孔喉特征参数。
7.根据权利要求1所述的方法，其特征在于，所述获取所述岩石样本对应的第一毛管压力曲线，包括：
通过高压压汞实验获得所述岩石样本的原始毛管压力曲线；
对所述原始毛管压力曲线进行麻皮效应校正，获得所述第一毛管压力曲线。
8.一种孔喉特征参数的获取装置，其特征在于，所述装置包括：
特征信息获取模块，用于获取岩石样本的特征信息，所述特征信息包括所述岩石样本的长度、直径、孔隙度和渗透率；
特征数据获取模块，用于对所述岩石样本进行特征提取处理，获得所述岩石样本的第一类型孔洞的特征数据，所述第一类型孔洞包括表面溶洞和尺寸大于尺寸阈值的孔；
进汞饱和度获取模块，用于根据所述特征信息和所述第一类型孔洞的特征数据，获取所述第一类型孔洞的进汞饱和度；
第一毛管压力曲线获取模块，用于获取所述岩石样本对应的第一毛管压力曲线，所述第一毛管压力曲线用于指示岩石样本中第二类型孔洞的体积占比，所述第二类型孔洞的体积占比是所述第二类型孔洞的体积占所述岩石样本的总孔隙体积的比例，所述第二类型孔洞是尺寸小于尺寸阈值的孔隙；
第二毛管压力曲线获取模块，用于根据所述第一毛管压力曲线、所述进汞饱和度以及所述第一类型孔洞的特征信息，获取第二毛管压力曲线；
孔喉特征参数获取模块，用于根据所述第二毛管压力曲线获取所述岩石样本的孔喉特征参数，所述孔喉特征参数包括最大孔隙半径、平均孔隙半径、中值、孔喉半径均值、分选系数、歪度、峰态、均质系数。
9.一种计算机设备，其特征在于，所述计算机设备包含处理器和存储器，所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集由所述处理器加载并执行以实现如权利要求1至7任一所述的孔喉特征参数的获取方法。
10.一种计算机可读存储介质，其特征在于，所述存储介质中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集由处理器加载并执行以实现如权利要求1至7任一所述的孔喉特征参数的获取方法。
CN202010984074.9A 2020-09-18 2020-09-18 孔喉特征参数的获取方法、装置、计算机设备及存储介质 Active CN112326523B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010984074.9A CN112326523B (zh) 2020-09-18 2020-09-18 孔喉特征参数的获取方法、装置、计算机设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010984074.9A CN112326523B (zh) 2020-09-18 2020-09-18 孔喉特征参数的获取方法、装置、计算机设备及存储介质

Publications (2)

Publication Number Publication Date
CN112326523A true CN112326523A (zh) 2021-02-05
CN112326523B CN112326523B (zh) 2024-06-25

Family

ID=74304513

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010984074.9A Active CN112326523B (zh) 2020-09-18 2020-09-18 孔喉特征参数的获取方法、装置、计算机设备及存储介质

Country Status (1)

Country Link
CN (1) CN112326523B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115506758A (zh) * 2021-06-04 2022-12-23 中国石油天然气股份有限公司 一种致密储层排驱压力确定方法、装置、设备和存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103884633A (zh) * 2014-03-05 2014-06-25 中国石油天然气股份有限公司 一种确定岩石渗透率的方法及装置
CN103900942A (zh) * 2014-01-10 2014-07-02 中国石油大学(华东) 基于高压压汞分析的储层微观孔喉参数的连续表征方法
CN104458525A (zh) * 2013-09-13 2015-03-25 中国石油天然气股份有限公司 一种常规压汞曲线表征微观孔隙喉道特征的采集处理方法
CN105353419A (zh) * 2015-11-25 2016-02-24 中国地质大学（北京） 一种基于储集层分类的核磁毛管压力曲线构造方法
CN107132171A (zh) * 2017-04-13 2017-09-05 东北石油大学 一种基于压汞‑氮吸附联测数据确定致密储层孔径分布的方法
CN107153733A (zh) * 2017-05-09 2017-09-12 中国石油集团川庆钻探工程有限公司 一种复杂孔隙结构岩石的毛管压力曲线拟合方法
CN107290259A (zh) * 2016-03-30 2017-10-24 中国石油化工股份有限公司 低渗砂岩储层有效渗流孔喉半径的计算方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104458525A (zh) * 2013-09-13 2015-03-25 中国石油天然气股份有限公司 一种常规压汞曲线表征微观孔隙喉道特征的采集处理方法
CN103900942A (zh) * 2014-01-10 2014-07-02 中国石油大学(华东) 基于高压压汞分析的储层微观孔喉参数的连续表征方法
CN103884633A (zh) * 2014-03-05 2014-06-25 中国石油天然气股份有限公司 一种确定岩石渗透率的方法及装置
CN105353419A (zh) * 2015-11-25 2016-02-24 中国地质大学（北京） 一种基于储集层分类的核磁毛管压力曲线构造方法
CN107290259A (zh) * 2016-03-30 2017-10-24 中国石油化工股份有限公司 低渗砂岩储层有效渗流孔喉半径的计算方法
CN107132171A (zh) * 2017-04-13 2017-09-05 东北石油大学 一种基于压汞‑氮吸附联测数据确定致密储层孔径分布的方法
CN107153733A (zh) * 2017-05-09 2017-09-12 中国石油集团川庆钻探工程有限公司 一种复杂孔隙结构岩石的毛管压力曲线拟合方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115506758A (zh) * 2021-06-04 2022-12-23 中国石油天然气股份有限公司 一种致密储层排驱压力确定方法、装置、设备和存储介质
CN115506758B (zh) * 2021-06-04 2024-05-28 中国石油天然气股份有限公司 一种致密储层排驱压力确定方法、装置、设备和存储介质

Also Published As

Publication number Publication date
CN112326523B (zh) 2024-06-25

Similar Documents

Publication Publication Date Title
CN105372168B (zh) 一种混凝土孔隙三维结构的构建方法
US9625601B2 (en) Nuclear magnetic resonance rock sample analysis method and instrument with constant gradient field
US20170115242A1 (en) Measuring connectivity between different pore types in porous media
Bhat et al. Pulsar scintillation and the local bubble
CA3097362C (en) Determining permeability of porous media based on nuclear magnetic resonance measurement
CN105510204A (zh) 一种基于ct图像的渗透率预测方法
CN110231272A (zh) 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统
EP2872741A2 (en) Digital rock analysis systems and methods with reliable multiphase permeability determination
CN112326523B (zh) 孔喉特征参数的获取方法、装置、计算机设备及存储介质
CN106644879B (zh) 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置
CN111699381A (zh) 基于核磁共振弛豫时间的温度灵敏度测量多孔介质的润湿性
Ganju et al. Effect of particle characteristics on the evolution of particle size, particle morphology, and fabric of sands loaded under uniaxial compression
Bourne et al. Effect of formalin fixation on biexponential modeling of diffusion decay in prostate tissue
Payton et al. The influence of grain shape and size on the relationship between porosity and permeability in sandstone: A digital approach
CN109387468A (zh) 页岩储层纳米孔隙结构特征参数测试分析方法及系统
CN112329975A (zh) 油气资源预测方法、装置、计算机设备及存储介质
Nistor et al. Tomographic reconstruction and morphological analysis of rigid polyurethane foams
CN112734143A (zh) 一种孔隙结构分形表征的方法、装置及系统
CN114112856A (zh) 一种含铀砂岩全尺寸孔径结构表征方法及系统
US11566992B2 (en) Method and device for characterizing microbial carbonate pores, and server
CN111693427B (zh) 油气藏流体可动性的分析方法
Gardner et al. The measurement of gas bubble size distributions in a three phase laboratory gassy soil
CN115963125A (zh) 岩石物理模型建立方法、装置、计算机设备及存储介质
CN113945497A (zh) 油气藏储层流体的可动性的评价方法
CN113376069A (zh) 致密储层物性下限的确定方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant