火山岩微观结构孔隙表征方法、装置、设备及存储介质与流程

1.本技术实施例涉及岩石结构表征技术领域，特别涉及一种火山岩微观结构孔隙表征方法、装置、设备及存储介质。


背景技术：

2.近年来，随着石油与天然气勘探程度的进一步发展，常规油气藏的勘探愈发困难，油气资源的勘探领域逐渐转向非常规油气。
3.迄今，作为非常规油气资源之一的火山岩储层已在我国勘探开发方面取得了重大进展。相对于常规储层，火山岩储层地质特征特殊，进而导致储层岩石微观孔喉匹配关系更加复杂，且难以准确表征，这给开采开发带来众多难题。为了准确表征油气藏储层地质特征，更加高效、经济地开采与改造油气藏，需要对火山岩储层微观非均质性进行剖析研究，精细、准确对火山岩储层进行岩石微观结构表征。然而，如何对火山岩储层进行准确地微观结构表征仍是待解决的问题。


技术实现要素：

4.本技术实施例提供了一种火山岩微观结构孔隙表征方法、装置、设备及存储介质，可以得到更准确的孔隙半径，从而可以实现对火山岩致密岩样进行精确的微观结构表征。技术方案如下：
5.一方面，本技术实施例提供一种火山岩微观结构孔隙表征方法，所述方法包括：
6.获取火山岩岩心的第一岩样的至少一个孔隙图像；
7.确定至少一个所述孔隙图像中的孔隙的孔隙轮廓；
8.基于至少一个所述孔隙的孔隙轮廓，确定至少一个所述孔隙的孔隙形状；
9.基于至少一个所述孔隙的孔隙形状，确定至少一个所述孔隙的初始孔隙半径；
10.基于至少一个所述孔隙的初始孔隙半径，确定所述第一岩样的孔隙半径范围；
11.基于所述第一岩样的孔隙半径范围，确定与所述孔隙半径范围对应的压力范围；
12.从关系曲线中确定与所述压力范围对应的目标关系曲线，所述关系曲线是指对所述火山岩岩心的第二岩样进行恒速压汞实验得到的压力与体积之间的关系曲线；
13.基于所述目标关系曲线，得到压力与体积之间的标准关系曲线；
14.从所述标准关系曲线中确定至少一个所述孔隙的目标体积值；
15.基于至少一个所述孔隙的目标体积值，确定至少一个所述孔隙的等效孔隙半径；
16.基于至少一个所述孔隙的孔隙形状和所述等效孔隙半径，确定至少一个所述孔隙的最终孔隙半径。
17.另一方面，本技术实施例提供一种火山岩微观结构孔隙表征装置，所述装置包括：
18.图像获取模块，用于获取火山岩岩心的第一岩样的至少一个孔隙图像；
19.轮廓确定模块，用于确定至少一个所述孔隙图像中的孔隙的孔隙轮廓；
20.形状确定模块，用于基于至少一个所述孔隙的孔隙轮廓，确定至少一个所述孔隙
的孔隙形状；
21.半径确定模块，用于基于至少一个所述孔隙的孔隙形状，确定至少一个所述孔隙的初始孔隙半径；
22.范围确定模块，用于基于至少一个所述孔隙的初始孔隙半径，确定所述第一岩样的孔隙半径范围；
23.所述范围确定模块，还用于基于所述第一岩样的孔隙半径范围，确定与所述孔隙半径范围对应的压力范围；
24.曲线确定模块，用于从关系曲线中确定与所述压力范围对应的目标关系曲线，所述关系曲线是指对所述火山岩岩心的第二岩样进行恒速压汞实验得到的压力与体积之间的关系曲线；
25.所述曲线确定模块，还用于基于所述目标关系曲线，得到压力与体积之间的标准关系曲线；
26.体积确定模块，用于从所述标准关系曲线中确定至少一个所述孔隙的目标体积值；
27.所述半径确定模块，还用于基于至少一个所述孔隙的目标体积值，确定至少一个所述孔隙的等效孔隙半径；
28.所述半径确定模块，还用于基于至少一个所述孔隙的孔隙形状和所述等效孔隙半径，确定至少一个所述孔隙的最终孔隙半径。
29.另一方面，本技术实施例提供一种计算机设备，所述计算机设备包括处理器和存储器，所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集由所述处理器加载并执行以实现如上述方面所述的火山岩微观结构孔隙表征方法。
30.再一方面，本技术实施例提供一种计算机可读存储介质，所述存储介质中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集由处理器加载并执行以实现如上述方面所述的火山岩微观结构孔隙表征方法。
31.本技术实施例提供的技术方案带来的有益效果可以包括：
32.通过确定每个孔隙的孔隙形状，基于孔隙形状确定初始孔隙半径；基于初始孔隙半径确定最终标准的恒速压汞曲线；最后基于恒速压汞曲线和孔隙形状确定孔隙的最终孔隙半径，本技术实施例利用扫描电镜对恒速压汞结果中的孔隙大小进行修正，得到更准确的孔隙半径，从而可以实现对火山岩致密岩样进行精确的微观结构表征。
附图说明
33.为了更清楚地说明本技术实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
34.图1是本技术一个实施例提供的火山岩微观结构孔隙表征方法的流程图；
35.图2是本技术一个实施例提供的确定标准关系曲线的示意图；
36.图3是本技术一个实施例提供的确定形状分布频率的示意图；
37.图4是相关技术提供的孔隙大小分布频率的示意图；
38.图5是本技术实施例提供的孔隙大小分布频率的示意图；
39.图6是本技术一个实施例提供的火山岩微观结构孔隙表征装置的示意图；
40.图7是本技术一个实施例提供的计算机设备的结构示意图。
具体实施方式
41.为使本技术的目的、技术方案和优点更加清楚，下面将结合附图对本技术实施方式作进一步地详细描述。
42.在相关技术中，岩石微观结构表征的技术手段众多，根据其性质可分为定性表征手段和定量表征手段。定性表征主要研究孔隙和吼道的类型、形态分布和连通性，使用的主要技术手段有铸体薄片、场发射扫描电镜(field emission scanning electron microscope，fesem)和环境扫描电镜(environment scanning electron microscope，esem)等。扫描电镜和铸体薄片可以直接观察致密储层的孔喉结构，主要用于定性评价，同时也可提供定量表征的结果；然而，铸体薄片对于微观尺度孔隙分辨率较低，扫描电镜随着分辨率的提高会导致样品分析视域的减小，在后续的图像处理后对二维孔喉结构信息的提取有限，需要大量的实验数据来保证结果具有代表性。定量表征是对孔喉大小和分布的定量分析与参数表征，现阶段使用的主要技术手段有恒速压泵、高压压泵及核磁共振等。基于汞对大多数造岩矿物成分均为非润湿相的认识，压汞法以恒定压力或者恒定速度向岩样中注入汞，根据得到的进汞量与进汞压力之间的实验数据，计算获取岩样的孔喉大小及其分布曲线，其中恒速压汞法可以区分岩样的孔隙和喉道，实现岩样孔隙和喉道的大小与分布表征。但是在数据分析处理中，一般将岩样孔隙等效为球体进而确定出孔喉半径，这会引起计算结果与实际情况间表现出较大的差异。核磁共振技术通过测定饱和岩样的t2谱曲线，进而计算得到岩样的孔隙度、孔径大小与分布等参数，相较其他方法更加快速、高效且无损，只是该方法不能直接获取孔喉半径分布信息，需要借助其他定量表征手段。
43.请参考图1，其示出了本技术一个实施例提供的火山岩微观结构孔隙表征方法的流程图。该方法可应用于诸如pc(personal computer，个人计算机)之类的具有数据处理和存储能力的计算机设备中，例如，计算机设备可以是终端或计算机设备。上述方法可以包括以下几个步骤。
44.步骤101，获取火山岩岩心的第一岩样的至少一个孔隙图像。
45.火山岩是火山爆发后形成的多孔形石材，在可能的实现方式中，本技术实施例中的火山岩还可以称之为火山岩储层。岩心是指根据地质工作或工程的需要，使用岩心环状钻头及其他取心工具，从孔内取出的圆柱状岩石样品，岩心是研究和了解地下地质和矿产情况的重要实物物质资料。
46.示例性地，孔隙图像是指包含孔隙的图像。孔隙是指被岩石颗粒包围的较大储集空间，是流体的基本储集空间。
47.示例性地，技术人员取一块火山岩储层的柱塞岩样，依据《gb/t 29172-2012岩心分析方法》并进行洗盐和烘干，测量岩心的孔隙度及渗透率等基本物性(基本物性对岩心的选取以及后续处理流程有个参考作用)。技术人员将岩心切割成两段，分别记为第一岩样和
第二岩样，然后再次进行烘干处理。技术人员选取第一岩样开展扫描电镜实验，依据《sy/t 5362-2014岩石样品扫描电子显微镜分析方法》中的规定利用扫描电子显微镜观察新鲜断面样品(第一岩样)：技术人员首先在较低倍数(例如，500倍)下扫描获取第一岩样的全局图像；然后，将上述全局图像划分为不同区域，并提高显示镜的倍数(例如，从500倍提升到1000倍)，拍摄第一岩样中的每一处孔隙图像，得到火山岩岩心的第一岩样的至少一个孔隙图像。技术人员可以将上述火山岩岩心的第一岩样的至少一个孔隙图像上传至计算机设备，以使得计算机设备对该至少一个孔隙图像进行处理。
48.步骤102，确定至少一个孔隙图像中的孔隙的孔隙轮廓。
49.在可能的实现方式中，通过图像处理软件确定至少一个孔隙图像中的孔隙的孔隙轮廓。例如，通过photoshop软件提取出该至少一个孔隙图像中的孔隙的孔隙轮廓。
50.步骤103，基于至少一个孔隙的孔隙轮廓，确定至少一个孔隙的孔隙形状。
51.孔隙的孔隙形状可以是指孔隙的截面形状。
52.在可能的实现方式中，计算机设备可以基于孔隙的孔隙轮廓特征，确定出该孔隙的孔隙形状，例如，计算机设备可以基于孔隙轮廓，确定孔隙的孔隙因子，然后再基于孔隙的孔隙因子确定出孔隙形状。有关计算机设备如何确定孔隙的孔隙形状，此处先不作介绍说明。
53.步骤104，基于至少一个孔隙的孔隙形状，确定至少一个孔隙的初始孔隙半径。
54.初始孔隙半径是指基于孔隙形状确定的初始的半径。
55.在可能的实现方式中，孔隙形状包括以下任意一种：第一三边形(三边形1)、第二三边形(三边形2)、第三三边形(三边形3)、第一四边形(四边形1)、第二四边形(四边形2)、第一五边形(五边形1)和第二五边形(五边形2)。第一三边形、第二三边形和第三三边形是互不相同的三边形，第一三边形和第二三边形是不规则的三边形，第三三边形是规则的三角形，第一四边形和第二四边形是两个不相同的四边形，第二四边形是规则的四边形，第一五边形和第二五边形是不相同的五边形，第二五边形是指规则的五边形(需要说明的是，本技术实施例中的三边形是指具有三个顶点的图形，四边形是指具有四个顶点的图形，五边形是指具有五个顶点的图形)。
56.在孔隙形状为第一三边形的情况下，初始孔隙半径可以为第一三边形的宽度的一半；在孔隙形状为第二三边形的情况下，初始孔隙半径可以为第二三边形的高度的一半；在孔隙形状为第三三边形的情况下，初始孔隙半径可以通过如下公式确定：初始孔隙半径＝2*第三三边形的面积/第三三边形的周长；在孔隙形状为第一四边形的情况下，初始孔隙半径可以为第一四边形的高度的一半；在孔隙形状为第二四边形的情况下，初始孔隙半径可以通过如下公式确定：初始孔隙半径＝2*第二四边形的面积/第二四边形的周长；在孔隙形状为第一五边形的情况下，初始孔隙半径可以为第一五边形的宽度的一半；在孔隙形状为第二五边形的情况下，初始孔隙半径可以为第二五边形的面积/第二五边形的周长。
57.步骤105，基于至少一个孔隙的初始孔隙半径，确定第一岩样的孔隙半径范围。
58.示例性地，将至少一个孔隙的初始孔隙半径中最小的初始孔隙半径确定为第一岩样的孔隙半径范围的下限，将至少一个孔隙的初始孔隙半径中最大的初始孔隙半径确定为第一岩样的孔隙半径范围的上限。
59.步骤106，基于第一岩样的孔隙半径范围，确定与孔隙半径范围对应的压力范围。
60.在可能的实现方式中，通过如下公式确定与孔隙半径r对应的压力值pc：
[0061][0062]
示例性地，计算机设备可以计算最小的初始孔隙半径对应的压力值，以及最大的初始孔隙半径对应的压力值，将最小的初始孔隙半径对应的压力值确定为压力范围的下限，将最大的初始孔隙半径对应的压力值确定为压力范围的上限。
[0063]
步骤107，从关系曲线中确定与压力范围对应的目标关系曲线。
[0064]
在本技术实施例中，关系曲线是指对火山岩岩心的第二岩样进行恒速压汞实验得到的压力与体积之间的关系曲线。示例性地，技术人员对第二岩样开展恒速压汞实验，得到第二岩样的注汞压力(简称压力)和注汞体积(简称体积)的关系曲线。示例性地，第一岩样的长度为2厘米，第二岩样的长度为3厘米。
[0065]
汞对绝大多数造岩矿物来说都是非润湿的。如果对汞施加压力，当注入汞的压力达到孔隙喉道的毛管压力时，汞就会克服毛管阻力进入孔隙内，根据不断注入汞的孔隙体积百分数和对应压力，便能绘制出关系曲线。恒速压汞法是以较低的恒定流速将汞注入到岩石孔隙中，进汞速度低，可认为是准静态的过程。汞在注入过程中，压力是变化的，通过监测压力升降情况可将岩石内部孔隙和喉道分开。汞先进入喉道1，压力上升到一定值后，汞突破喉道1进入孔隙1，此时，压力迅速降低，当孔隙1充满汞后，压力开始上升，此时，汞正进入喉道2中，压力上升到一定值后，汞突破喉道2进入孔隙2中，压力再次迅速降低，同理，汞依次充满其它喉道和其它孔隙，直至与喉道1连通的所有喉道及控制的孔隙充满后，汞才会流入比喉道1半径小的喉道所控制的孔隙单元。重复以上过程，当压力达到测试仪器的压力上限时，实验结束。据进汞过程中压力涨落可获得孔隙结构方面信息，定量获取喉道和孔隙数目，得到孔隙半径分布、喉道半径分布等孔隙结构特征参数。
[0066]
示例性地，关系曲线的横坐标为体积，纵坐标为压力。将关系曲线中与压力范围对应的曲线确定为目标关系曲线。
[0067]
步骤108，基于目标关系曲线，得到压力与体积之间的标准关系曲线。
[0068]
将目标关系曲线进行放大，结合放大后的目标关系曲线，剔除目标关系曲线中的异常点，得到优化后的点，由于目标关系曲线其实是由多个点组成的曲线，所以对目标关系曲线进行放大后，可以得到多个点的起伏情况，剔除目标关系曲线中起伏情况异常的点，得到优化后的点。基于优化后的点绘制压力和体积的标准关系曲线，标准关系曲线也可以称之为标准的“m”曲线。如图2所示，其示出了本技术一个实施例提供的确定标准关系曲线的示意图。计算机设备基于关系曲线210和压力范围，确定目标关系曲线，假设孔隙半径范围为0.226微米-0.533微米，基于该孔隙半径范围确定出压力范围，然后在关系曲线210中截取目标关系曲线220，放大目标关系曲线220，并剔除目标关系曲线220中的异常点，得到优化后的点，基于优化后的点绘制出压力与体积的标准关系曲线230。(需要说明的是，图2仅是示例性地，在可能的实现方式中，关系曲线还可以是其它类型的曲线，本技术实施例对此不作限定)。
[0069]
步骤109，从标准关系曲线中确定至少一个孔隙的目标体积值。
[0070]
示例性地，如图2所示，v1代表一个孔隙的目标体积值，v2代表一个孔隙的目标体积值，v3代表一个孔隙的目标体积值，以此类推，得到至少一个孔隙的目标体积值。
[0071]
在可能的实现方式中，利用插入排序算法从标准关系曲线中确定出至少一个孔隙的目标体积值。
[0072]
步骤110，基于至少一个孔隙的目标体积值，确定至少一个孔隙的等效孔隙半径。
[0073]
示例性地，将孔隙等效为球体，基于球体体积公式和目标体积值，确定至少一个孔隙的等效孔隙半径。
[0074]
等效孔隙半径是指基于等效孔隙确定的半径。
[0075]
球体体积公式为其中，v表示体积，r表示半径。在已知上述公式中的v时，可以反推出r，因此，基于至少一个孔隙的目标体积值，可以确定至少一个孔隙的等效孔隙半径。
[0076]
步骤111，基于至少一个孔隙的孔隙形状和等效孔隙半径，确定至少一个孔隙的最终孔隙半径。
[0077]
最终孔隙半径是指最终确定的孔隙的半径。
[0078]
示例性地，计算机设备可以通过如下方式确定至少一个孔隙的最终孔隙半径：
[0079]
1、确定与至少一个孔隙的孔隙形状对应的半径确定方式。
[0080]
半径确定方式是指确定最终孔隙半径的方式。
[0081]
由于不同的孔隙形状，其确定最终孔隙半径的方式不一样，因此，需要先基于孔隙形状确定半径确定方式。
[0082]
2、基于半径确定方式与等效孔隙半径，确定至少一个孔隙的最终孔隙半径。
[0083]
在可能的实现方式中，在孔隙的孔隙形状为三边形的情况下，将第一数值的等效孔隙半径确定为最终孔隙半径。
[0084]
在孔隙的孔隙形状为四边形的情况下，将第二数值的等效孔隙半径确定为最终孔隙半径。
[0085]
在孔隙的孔隙形状为五边形的情况下，将第三数值的等效孔隙半径确定为最终孔隙半径。
[0086]
示例性地，第一数值、第二数值和第三数值是互不相同的数值。
[0087]
在可能的实现方式中，在孔隙形状为三边形的情况下，计算机设备可以通过如下方式确定最终孔隙半径：
[0088][0089]
其中，r
三边形内切球半径
表示孔隙形状为三边形的孔隙的最终孔隙半径，r
等效孔隙半径
表示等效孔隙半径。
[0090]
在可能的实现方式中，在孔隙形状为四边形的情况下，计算机设备可以通过如下方式确定最终孔隙半径：
[0091][0092]
其中，r
四边形内切球半径
表示孔隙形状为四边形的孔隙的最终孔隙半径，r
等效孔隙半径
表示等效孔隙半径。
[0093]
在可能的实现方式中，在孔隙形状为五边形的情况下，计算机设备可以通过如下
方式确定最终孔隙半径：
[0094][0095]
其中，r
五边形内切球半径
表示孔隙形状为五边形的孔隙的最终孔隙半径，r
等效孔隙半径
表示等效孔隙半径。
[0096]
火山岩储层一般都是经历了强烈压实和胶结(硅质胶结、黏土胶结)作用，这导致其孔隙形状和大小变化特征异常复杂；如果要准确地表征火山岩微观结构，则需结合恒速压汞实验和其他方法的优点。本技术实施例提供了一种结合恒速压汞实验的火山岩微观孔隙结构表征方法，针对火山岩致密岩样孔隙结构，利用扫描电镜对恒速压汞结构中孔隙及喉道大小进行修正，从而对火山岩致密岩样进行精确的微观结构表征。
[0097]
综上所述，本技术实施例提供的技术方案中，通过确定每个孔隙的孔隙形状，基于孔隙形状确定初始孔隙半径；基于初始孔隙半径确定最终标准的恒速压汞曲线；最后基于恒速压汞曲线和孔隙形状确定孔隙的最终孔隙半径，本技术实施例利用扫描电镜对恒速压汞结果中的孔隙大小进行修正，得到更准确的孔隙半径，从而可以实现对火山岩致密岩样进行精确的微观结构表征。
[0098]
在示意性实施例中，通过如下方式确定孔隙的孔隙形状：
[0099]
第一、基于至少一个孔隙的孔隙轮廓，确定至少一个孔隙的截面积和至少一个孔隙的周长。
[0100]
在可能的实现方式中，通过imagej软件对孔隙轮廓进行处理，得到至少一个孔隙的截面积和周长。
[0101]
第二、基于至少一个孔隙的截面积和至少一个孔隙的周长，确定至少一个孔隙的形状因子。
[0102]
在可能的实现方式中，计算机设备通过如下公式确定形状因子g：
[0103][0104]
其中，a表示孔隙的截面积，p表示孔隙的周长，示例性地，公式中的截面积的单位可以是cm2，周长的单位可以是cm。
[0105]
第三、基于至少一个孔隙的形状因子，确定至少一个孔隙的孔隙形状。
[0106]
如图3所示，通过对第一岩样的全局图像310划分不同区域，拍摄得到第一岩样中的每一处孔隙图像320；确定孔隙图像320中的孔隙的孔隙轮廓330；基于孔隙轮廓330确定孔隙的孔隙因子；计算机设备可以通过形状因子分布图340(形状因子分布图340用于指示孔隙形状与形状因子的对应关系)确定孔隙的孔隙形状。示例性地，若形状因子位于三边形2和四边形1的交集处，则选取与形状因子更接近的下限对应的形状作为孔隙的孔隙形状。例如，形状因子为0.0218，则将四边形1确定为孔隙的孔隙形状，因为0.0281更接近0.0217(即，四边形1对应的形状因子的下限)。
[0107]
在示意性实施例中，计算机设备确定至少一个孔隙的孔隙形状之后，还可以执行下述步骤：
[0108]
第一、基于至少一个孔隙的孔隙形状，确定孔隙形状为同一类型的孔隙的个数。
[0109]
计算机设备确定出各个孔隙的孔隙形状后，可以孔隙形状为同一类型的孔隙的个
数。例如，计算机设备可以确定孔隙形状为三边形1的孔隙的个数、孔隙形状为四边形1的孔隙的个数、孔隙形状为五边形1的孔隙的个数等。
[0110]
第二、基于孔隙形状为同一类型的孔隙的个数以及至少一个孔隙的总个数，确定孔隙形状为不同类型的孔隙各自对应的形状分布频率，形状分布频率用于指示孔隙形状为同一类型的孔隙的个数与至少一个孔隙的总个数的比例。
[0111]
示例性地，计算机设备确定出孔隙形状为某一类型的孔隙的个数后，可以将该值与孔隙的总个数的比值，确定为该孔隙形状为某一类型的孔隙的形状分布频率。如图3所示，其示出了本技术一个实施例提供的形状分布频率的示意图350。
[0112]
在示意性实施例中，计算机设备确定出孔隙的最终孔隙半径之后，还可以下述步骤：基于至少一个孔隙的孔隙形状、最终孔隙半径，确定岩心的最终等效半径。
[0113]
示例性地，计算机设备可以基于孔隙形状为三边形的孔隙的最终孔隙半径、孔隙形状为四边形的孔隙的最终孔隙半径、孔隙形状为五边形的孔隙的最终孔隙半径，确定岩心的最终等效半径。
[0114]
在可能的实现方式中，计算机设备可以通过如下方式确定岩心的最终等效半径：基于不同类型的孔隙形状的孔隙各自对应的形状分布频率与不同类型的孔隙形状的孔隙各自对应的最终孔隙半径，确定岩心的最终等效半径。
[0115]
在可能的实现方式中，岩心的最终等效半径可以通过如下方式确定：
[0116]
1、基于至少一个孔隙的最终孔隙半径，确定孔隙形状为同一类型的孔隙的目标孔隙半径。
[0117]
示例性地，将孔隙形状为同一类型的孔隙各自的最终孔隙半径相加求平均，得到孔隙形状为同一类型的孔隙的目标孔隙半径。
[0118]
2、将各个孔隙形状为同一类型的孔隙对应的形状分布频率与目标孔隙半径的乘积之和，确定为岩心的最终等效半径。
[0119]
以图3所示的形状分布频率的示意图为例进行介绍说明，假设第一岩样中包括孔隙形状为三边形3、四边形2、五边形1和五边形2的孔隙，孔隙形状为三边形3的孔隙的形状分布频率为0.3，孔隙形状为四边形2的孔隙的形状分布频率为0.4，孔隙形状为五边形1的孔隙的形状分布频率为0.2，孔隙形状为五边形2的孔隙的形状分布频率为0.1，假设孔隙形状为三边形3的孔隙的目标孔隙半径为半径1，孔隙形状为四边形2的孔隙的目标孔隙半径为半径2，孔隙形状为五边形1的孔隙的目标孔隙半径为半径3，孔隙形状为五边形2的孔隙的目标孔隙半径为半径4。则岩心的最终等效半径为(0.3*半径1 0.4*半径2 0.2*半径3 0.1*半径4)。
[0120]
针对某一盆地二叠系火山岩储层某一岩样，基于扫描电镜实验和恒速压汞实验的结果，依据本技术实施例提出的流程得到修正后的火山岩孔隙半径(即最终孔隙半径)与分布频率，孔隙半径与分布频率可以称之为孔隙大小分布频率。对比分析发现，从扫描电镜实验中直观获取的孔隙大小分布频率与基于本技术实施例修正后的孔隙大小分布频率一致，这说明本技术实施例提出的结合恒速压汞实验的火山岩微观孔隙结构表征方法是可行的，且可操作性较强。图4示出了相关技术中基于恒速压汞实验未修正的孔隙大小分布频率的示意图400，图5示出了基于本技术实施例修正后的孔隙大小分布频率的示意图500。
[0121]
下述为本技术装置实施例，可以用于执行本技术方法实施例。对于本技术装置实
施例中未披露的细节，请参照本技术方法实施例。
[0122]
请参考图6其示出了本技术一个实施例提供的一种火山岩微观结构孔隙表征装置的框图。该装置具有实现上述方法示例的功能，所述功能可以由硬件实现，也可以由硬件执行相应的软件实现。该装置600可以包括：图像获取模块610、轮廓确定模块620、形状确定模块630、半径确定模块640、范围确定模块650、曲线确定模块660、体积确定模块670。
[0123]
图像获取模块610，用于获取火山岩岩心的第一岩样的至少一个孔隙图像；
[0124]
轮廓确定模块620，用于确定至少一个所述孔隙图像中的孔隙的孔隙轮廓；
[0125]
形状确定模块630，用于基于至少一个所述孔隙的孔隙轮廓，确定至少一个所述孔隙的孔隙形状；
[0126]
半径确定模块640，用于基于至少一个所述孔隙的孔隙形状，确定至少一个所述孔隙的初始孔隙半径；
[0127]
范围确定模块650，用于基于至少一个所述孔隙的初始孔隙半径，确定所述第一岩样的孔隙半径范围；
[0128]
所述范围确定模块650，还用于基于所述第一岩样的孔隙半径范围，确定与所述孔隙半径范围对应的压力范围；
[0129]
曲线确定模块660，用于从关系曲线中确定与所述压力范围对应的目标关系曲线，所述关系曲线是指对所述火山岩岩心的第二岩样进行恒速压汞实验得到的压力与体积之间的关系曲线；
[0130]
所述曲线确定模块660，还用于基于所述目标关系曲线，得到压力与体积之间的标准关系曲线；
[0131]
体积确定模块670，用于从所述标准关系曲线中确定至少一个所述孔隙的目标体积值；
[0132]
所述半径确定模块640，还用于基于至少一个所述孔隙的目标体积值，确定至少一个所述孔隙的等效孔隙半径；
[0133]
所述半径确定模块640，还用于基于至少一个所述孔隙的孔隙形状和所述等效孔隙半径，确定至少一个所述孔隙的最终孔隙半径。
[0134]
综上所述，本技术实施例提供的技术方案中，通过确定每个孔隙的孔隙形状，基于孔隙形状确定初始孔隙半径；基于初始孔隙半径确定最终标准的恒速压汞曲线；最后基于恒速压汞曲线和孔隙形状确定孔隙的最终孔隙半径，本技术实施例利用扫描电镜对恒速压汞结果中的孔隙大小进行修正，得到更准确的孔隙半径，从而可以实现对火山岩致密岩样进行精确的微观结构表征。
[0135]
在示意性实施例中，所述形状确定模块630，用于：
[0136]
基于至少一个所述孔隙的孔隙轮廓，确定至少一个所述孔隙的截面积和至少一个所述孔隙的周长；
[0137]
基于至少一个所述孔隙的截面积和至少一个所述孔隙的周长，确定至少一个所述孔隙的形状因子；
[0138]
基于至少一个所述孔隙的形状因子，确定至少一个所述孔隙的孔隙形状。
[0139]
在示意性实施例中，所述半径确定模块640，包括：方式确定单元、半径确定单元(图中未示出)。
[0140]
方式确定单元，用于确定与至少一个所述孔隙的孔隙形状对应的半径确定方式；
[0141]
半径确定单元，用于基于所述半径确定方式与所述等效孔隙半径，确定至少一个所述孔隙的最终孔隙半径。
[0142]
在示意性实施例中，所述半径确定单元，用于：
[0143]
在所述孔隙的孔隙形状为三边形的情况下，将第一数值的所述等效孔隙半径确定为所述最终孔隙半径；
[0144]
在所述孔隙的孔隙形状为四边形的情况下，将第二数值的所述等效孔隙半径确定为所述最终孔隙半径；
[0145]
在所述孔隙的孔隙形状为五边形的情况下，将第三数值的所述等效孔隙半径确定为所述最终孔隙半径。
[0146]
在示意性实施例中，所述装置还包括：频率确定模块(图中未示出)。
[0147]
基于至少一个所述孔隙的孔隙形状，确定孔隙形状为同一类型的孔隙的个数；
[0148]
基于所述孔隙形状为同一类型的孔隙的个数以及至少一个所述孔隙的总个数，确定孔隙形状为不同类型的孔隙各自对应的形状分布频率，所述形状分布频率用于指示所述孔隙形状为同一类型的孔隙的个数与至少一个所述孔隙的总个数的比例。
[0149]
在示意性实施例中，所述半径确定模块640，还用于：
[0150]
基于至少一个所述孔隙的孔隙形状、所述最终孔隙半径，确定所述岩心的最终等效半径。
[0151]
在示意性实施例中，所述半径确定模块640，用于：
[0152]
基于至少一个所述孔隙的最终孔隙半径，确定孔隙形状为同一类型的孔隙的目标孔隙半径；
[0153]
将各个孔隙形状为同一类型的孔隙对应的形状分布频率与所述目标孔隙半径的乘积之和，确定为所述岩心的最终等效半径。
[0154]
需要说明的一点是，上述实施例提供的装置在实现其功能时，仅以上述各个功能模块的划分进行举例说明，实际应用中，可以根据实际需要而将上述功能分配由不同的功能模块完成，即将设备的内容结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。
[0155]
关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
[0156]
请参考图7，其示出了本技术一个实施例提供的计算机设备700的结构示意图。该计算机设备700可用于实施上述实施例中提供的计算机设备侧的火山岩微观结构孔隙表征方法。具体来讲：
[0157]
所述计算机设备700包括中央处理单元(central processing unit，cpu)701、包括ram(random access memory，随机存取存储器)702和rom(read-only memory，只读存储器)703的系统存储器704，以及连接系统存储器704和中央处理单元701的系统总线705。所述计算机设备700还包括帮助计算机内的各个器件之间传输信息的基本输入/输出系统(i/o系统，input/output系统)706，和用于存储操作系统713、应用程序714和其他程序模块715的大容量存储设备707。
[0158]
所述基本输入/输出系统706包括有用于显示信息的显示器708和用于用户输入信
息的诸如鼠标、键盘之类的输入设备709。其中所述显示器708和输入设备709都通过连接到系统总线705的输入输出控制器710连接到中央处理单元701。所述基本输入/输出系统706还可以包括输入输出控制器710以用于接收和处理来自键盘、鼠标、或电子触控笔等多个其他设备的输入。类似地，输入输出控制器710还提供输出到显示屏、打印机或其他类型的输出设备。
[0159]
所述大容量存储设备707通过连接到系统总线705的大容量存储控制器(未示出)连接到中央处理单元701。所述大容量存储设备707及其相关联的计算机可读介质为计算机设备700提供非易失性存储。也就是说，所述大容量存储设备707可以包括诸如硬盘或者cd-rom(compact disc read-only memory，只读光盘)驱动器之类的计算机可读介质(未示出)。
[0160]
不失一般性，所述计算机可读介质可以包括计算机存储介质和通信介质。计算机存储介质包括以用于存储诸如计算机可读指令、数据结构、程序模块或其他数据等信息的任何方法或技术实现的易失性和非易失性、可移动和不可移动介质。计算机存储介质包括ram、rom、eprom(erasable programmable read-only memory，可擦除可编程只读存储器)、eeprom(electrically erasable programmable read-only memory，电可擦可编程只读存储器)、闪存(flash memory)或其他固态存储设备，cd-rom、dvd(digital versatile disc，数字通用光盘)或其他光学存储、磁带盒、磁带、磁盘存储或其他磁性存储设备。当然，本领域技术人员可知所述计算机存储介质不局限于上述几种。上述的系统存储器704和大容量存储设备707可以统称为存储器。
[0161]
根据本技术的各种实施例，所述计算机设备700还可以通过诸如因特网等网络连接到网络上的远程计算机运行。也即计算机设备700可以通过连接在所述系统总线705上的网络接口单元711连接到网络712，或者说，也可以使用网络接口单元711来连接到其他类型的网络或远程计算机系统(未示出)。
[0162]
所述存储器还包括一个或者一个以上的程序，所述一个或者一个以上程序存储于存储器中，且经配置以由一个或者一个以上处理器执行。上述一个或者一个以上程序包含用于实现上述计算机设备侧的火山岩微观结构孔隙表征方法的指令。
[0163]
在示意性实施例中，还提供了一种计算机设备，所述计算机设备包括处理器和存储器，所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集。所述至少一条指令、至少一段程序、代码集或指令集经配置以由一个或者一个以上处理器执行，以实现上述计算机设备侧的火山岩微观结构孔隙表征方法。
[0164]
在示意性实施例中，还提供了一种计算机可读存储介质，所述存储介质中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或所述指令集在被计算机设备的处理器执行时实现上述计算机设备侧的火山岩微观结构孔隙表征方法。
[0165]
可选地，上述计算机可读存储介质可以是rom(read-only memory，只读存储器)、ram(random access memory，随机存取存储器)、cd-rom(compact disc read-only memory，只读光盘)、磁带、软盘和光数据存储设备等。
[0166]
在示例性实施例中，还提供了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存储介质中。计算机
设备的处理器从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算机设备执行上述计算机设备侧的火山岩微观结构孔隙表征方法。
[0167]
以上所述仅为本技术的示例性实施例，并不用以限制本技术，凡在本技术的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本技术的保护范围之内。