致密砂岩储层孔隙的微观拓扑结构确定方法及装置

技术领域

本申请涉及油气勘探与开发技术领域，特别涉及致密砂岩储层孔隙的微观拓扑结构确定方法及装置。

背景技术

随着油气勘探开发技术的发展，致密油的比重逐年增加。而致密油的开采难度相比于常规原油要高很多，与常规砂岩储层相比致密储层微观孔隙结构比较复杂，一般表现为孔喉细小，孔隙空间分布不均，孔喉配置关系相对复杂等。致密储层孔喉主要以纳米、微米为主，具有很强的非均质性。致密储层微观孔隙结构对储层储集能力以及渗流能力、流体分布等产生直接影响。因此，更加清晰的了解孔隙结构对于致密储层物性等特征的影响，有利于更加深入的了解致密油藏的形成机理，控制和优化油气藏的开发以及更好的评价油气藏。

储层孔隙结构的分类评价研究具有重要的意义，它可以帮助我们更好的了解储层的质量，更好的寻找优质储层，更好的控制和优化采油工艺，从而提高采油效率和收益。储层孔隙结构的分类评价是指在表征孔隙结构的前提下，选择合适的孔喉结构参数，借助一定的分类评价方法和模型，综合运用各种数据研究储层的孔隙结构并进行分类的过程。

目前已经有多种方法和手段对孔隙结构进行表征，这些方法按照测量方式的角度可以分为直接法和间接法。孔隙结构的直接观测法是基于通过高分辨的分析设备直接观察孔隙结构的三维特征。但致密孔隙中发育大量纳米级孔隙及喉道，低于仪器的检测下限，故致密储层中的很多孔隙以及他们之间的连通关系不会被直接观测法检测到，进而使得测试结果不够准确。而间接测试法可以通过测试方法得当孔喉分布特征，但对于孔隙结构样式不能直接体现，且不同的方法具有不同的测试原理，因此在测量精度和范围上有所差距。

随着勘探开发的深入，用单一方法定性的孔隙结构及储层质量表征已经不能满足现阶段的开发需要，因此需要针对孔隙结构及储层质量进行细致、深入研究。

发明内容

本申请的目的是提供致密砂岩储层孔隙的微观拓扑结构确定方法及装置，以提供一种对于致密砂岩储层孔隙更为细致、深入的研究方法。

为解决上述技术问题，本说明书第一方面提供一种致密砂岩储层孔隙的微观拓扑结构确定方法，包括：确定致密储层的多个岩石样品薄片中的孔隙类型、孔径大小；确定所述多个岩石样品薄片中的喉道类型；确定所述多个岩石样品的孔隙度和渗透率；在对所述多个岩石样品分别进行高压压汞实验和恒速压汞实验时，记录实验数据；根据所述实验数据确定岩石样品中的孔喉大小及分布特征；确定所述多个岩石样品薄片中的孔隙类型与喉道类型的主要组合配置关系；对所述实验数据进行处理，得到全孔径下的分形维数；根据孔隙类型与喉道类型的主要组合配置关系、孔隙度、渗透率、全孔径孔喉的分形维数，确定致密砂岩储层孔隙的微观拓扑结构特征表，所述微观拓扑结构特征表中包括致密岩石的多个微观拓扑结构类别，所述微观拓扑结构特征表用于确定目标地区的油气勘探或开发方案。

在一些实施例中，致密砂岩储层孔隙的微观拓扑结构类别包括以下至少一种：球棍网状结构、星链网状结构、离散树状结构、复杂管束状结构。

在一些实施例中，确定致密砂岩储层孔隙的微观拓扑结构类别，包括：在孔隙类型与喉道类型的主要组合配置关系为大粒间孔-宽片状喉道组合、孔隙度大于10％、渗透率大于0.1mD、压汞曲线具有平台、反映大孔的分形维数分布在2.26～2.4、反映宽喉道的分形维数分布在2.04～2.11的情况下，确定岩石样品的微观拓扑结构类别为球棍网状结构；和/或，在孔隙类型与喉道类型的主要组合配置关系为小中粒间孔-窄片状喉道组合、孔隙度为8％～12％，渗透率为0.05mD～0.1mD、压汞曲线具有局部的平台特征、反映小孔的分形维数分布在2.46～2.81、反映窄喉道的分形维数分布在2.04～2.16的情况下，确定岩石样品的微观拓扑结构类别为星链网状结构；和/或，在孔隙类型与喉道类型的主要组合配置关系为粒内孔—窄片状喉道组合、孔隙度为5％～10％、渗透率为0.02mD～0.06mD、压汞曲线不具有平台，汞饱和度与压力呈指数关系、反映小孔的分形维数分布在2.32～2.81、反映窄喉道的分形维数分布在2.04～2.16的情况下，确定岩石样品的微观拓扑结构类别为离散树状结构；和/或，在孔隙类型与喉道类型的主要组合配置关系为微孔—极窄片状及管束状喉道组合、孔隙度为5％、渗透率小于0.03mD、压汞曲线不具平台特征、反映极窄喉道的分形维数分布在2.08～2.16的情况下，确定岩石样品的微观拓扑结构类别为复杂管束状结构。

在一些实施例中，在确定致密砂岩储层孔隙的微观拓扑结构特征表之后，还包括通过以下方法确定目标地区致密砂岩储层孔隙的微观拓扑结构类别：确定所述目标地区岩石样品薄片中的孔隙类型、孔径大小；确定所述目标地区岩石样品薄片中的喉道类型；确定所述目标地区岩石样品的孔隙度和渗透率；在对所述目标地区岩石样品分别进行高压压汞实验和恒速压汞实验时，记录实验数据；根据所述实验数据确定岩石样品中的孔喉大小及分布特征；确定所述目标地区岩石样品薄片中的孔隙类型与喉道类型的主要组合配置关系；根据实验数据，处理得到全孔径孔喉的分形维数；将所述目标地区岩石样品的孔隙类型与喉道类型的主要组合配置关系、孔隙度、渗透率、全孔径孔喉的分形维数与致密砂岩储层孔隙的微观拓扑结构特征表进行比对，将所述微观拓扑结构特征表中与所述目标地区岩石样品的各特征最相近的目标微观拓扑结构类别作为所述目标地区岩石样品的微观拓扑结构类别，以根据所述目标地区岩石样品的微观拓扑结构类别确定油气勘探或开发方案。

在一些实施例中，对所述实验数据进行处理，得到全孔径下的分形维数，包括：确定高压压汞和恒速压汞实验的各压力数据对应的孔喉大小及孔喉类型；将高压压汞实验的实验数据中的压力数据、汞饱和度数据分别取对数，构造第一数据集，所述第一数据集中的实验数据按照孔喉大小及类型划分至以下四个类别：大孔隙、宽喉道、小孔隙、窄喉道；对于所述第一数据集中的各个类别的实验数据，分别确定实验数据的拟合直线，得到高压压汞实验对应的多条拟合直线；根据拟合直线的斜率确定反映对应孔喉类别的分形维数；将恒压压汞实验的实验数据中的压力数据、汞饱和度数据分别取对数，构造第二数据集，将所述第二数据集中的数据按照孔喉大小及类别划分至以下四个类别：大孔隙、宽喉道、小孔隙、窄喉道；对于所述第二数据集中的各个类别的实验数据，分别确定实验数据的拟合直线，得到恒压压汞实验对应的多条拟合直线；根据拟合直线的斜率确定反映对应孔喉类别的分形维数；将所述高压压汞实验对应的多条拟合直线与所述恒压压汞实验对应的多条拟合直线进行拼接，得到全孔径下的分形维数。

在一些实施例中，将所述高压压汞实验对应的多条拟合直线与所述恒压压汞实验对应的多条拟合直线进行拼接，得到全孔径下的分形维数，包括：将所述高压压汞实验对应的多条拟合直线与所述恒压压汞实验对应的多条拟合直线进行拼接，拼接值采用所述高压压汞实验和恒速压汞实验对应的多条拟合直线中小孔喉和大孔喉分界点rg2及rh2的平均值。

在一些实施例中，根据拟合直线的斜率确定反映对应孔喉类别的分形维数，包括：计算目标拟合直线的斜率与2的和值；将所述和值作为反映目标孔喉类别的分形维数，所述目标孔喉类别是所述目标拟合直线关联的实验数据所对应的孔喉类别。

本说明书第二方面提供一种致密砂岩储层孔隙的微观拓扑结构确定装置，包括：第一确定单元，用于确定致密储层的多个岩石样品薄片中的孔隙类型、孔径大小；第二确定单元，用于确定所述多个岩石样品薄片中的喉道类型；第三确定单元，用于确定所述多个岩石样品的孔隙度和渗透率；记录单元，用于在对所述多个岩石样品分别进行高压压汞实验和恒速压汞实验时记录实验数据；第四确定单元，用于根据所述实验数据确定岩石样品中的孔喉大小及分布特征；第五确定单元，用于确定所述多个岩石样品薄片中的孔隙类型与喉道类型的主要组合配置关系；处理单元，用于对所述实验数据进行处理，得到全孔径下的分形维数；第六确定单元，用于根据孔隙类型与喉道类型的主要组合配置关系、孔隙度、渗透率、全孔径孔喉的分形维数，确定致密砂岩储层孔隙的微观拓扑结构特征表，所述微观拓扑结构特征表中包括致密岩石的多个微观拓扑结构类别，所述微观拓扑结构特征表用于确定目标地区的油气勘探或开发方案。

本说明书第三方面提供一种电子设备，包括：存储器和处理器，所述处理器和所述存储器之间互相通信连接，所述存储器中存储有计算机指令，所述处理器通过执行所述计算机指令，从而实现第一方面任一项所述方法的步骤。

本说明书第四方面提供一种计算机存储介质，所述计算机存储介质存储有计算机程序指令，所述计算机程序指令被执行时实现第一方面任一项所述方法的步骤。

本说明书提供的致密砂岩储层孔隙的微观拓扑结构确定方法及装置，提出了基于拓扑结构及孔喉配置关系的致密砂岩储层孔隙微观拓扑结构分类及表征方法。不同类型的孔隙网络模型是基于拓扑结构，将致密砂岩储层中不同类型孔隙及喉道的组合关系直观特征与孔隙结构的间接测试结果结合。该孔隙结构模型(也即微观拓扑结构类别)可以直接的反映孔隙-喉道特征，确定孔隙和喉道类型、大小、配置关系、孔喉分布特征等；并且基于拓扑结构的不同孔隙结构模型(也即微观拓扑结构类别)能够与不同方法的测试结果均有可对比性(如压汞)，不同的孔隙结构模型(也即微观拓扑结构类别)可以通过测试数据进行验证。

通过多类型孔喉拓扑结构孔喉模型，更全面的认识储层的孔隙结构特征，在孔隙结构的尺度、精度、样式等方面的精细表征，明确孔隙结构类型、模式及不同类型孔隙结构下储层质量的差异及流体渗透特征的差异，可以对储层进行分类，为下一步研究储层质量差异机理及差异分布特征奠定基础，对致密储层的表征和评价意义重大，可为油田后期勘探及开发提供可靠的地质依据。

附图说明

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

图1示出了本说明书提供的一种致密砂岩储层孔隙的微观拓扑结构确定方法的流程图；

图2示出了某区块(西233区块)长7段储层孔隙类型图版；

图3示出了某区块(西233区块)长7段储层孔隙类型分布直方图；

图4示出了某区块(西233区块)长7段储层喉道类型图版；

图5示出了某区块(西233区块)长7段储层致密砂岩高压压汞毛细管压力图；

图6示出了某区块(西233区块)长7段储层致密砂岩四种类型高压压汞毛细管压力图；

图7示出了某区块(西233区块)7段储层致密砂岩高压压汞孔喉半径分布图；

图8示出了Ⅰ型恒速压汞曲线及喉道半径分布直方图；

图9示出了II型恒速压汞曲线及喉道半径分布直方图；

图10示出了研究区恒速压汞所测孔隙半径分布直方图；

图11示出了研究区恒速压汞所测喉道半径分布直方图；

图12示出了某区块(西233区块)长7段储层致密砂岩孔喉组合类型图；

图13示出了研究区长7段致密砂岩样品高压压汞分形维数拟合曲线；

图14示出了研究区长7段致密砂岩样品恒速压汞分形维数拟合曲线；

图15示出了研究区长7段致密砂岩样品孔隙度与高压压汞分形维数相关性；

图16示出了研究区长7段致密砂岩样品孔隙度与恒速压汞分形维数相关性；

图17示出了研究区长7段致密砂岩样品高压压汞与恒速压汞全孔径分布特征；

图18示出了某区块(西233区块)长7段储层致密砂岩孔喉拓扑结构；

图19示出了某区块(西233区块)7段储层致密砂岩镜下孔喉拓扑结构模式图；

图20示出了致密砂岩储层孔隙的微观拓扑结构特征表的示意图。

具体实施方式

为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施方式中的附图，对本申请实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本申请一部分实施方式，而不是全部的实施方式。基于本申请中的实施方式，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式，都应当属于本申请保护的范围。

本说明书提供一种致密砂岩储层孔隙的微观拓扑结构确定方法，如图1所示，包括如下步骤：

S10：确定致密储层的多个岩石样品薄片中的孔隙类型、孔径大小。

例如，通过对研究区45口取心井铸体薄片观察和扫描电镜图像分析，储层孔隙类型较多，主要以粒间孔、粒内孔以及微孔隙为主，同时还有由于颗粒破裂作用形成的微裂缝。孔隙半径范围在2μm～115μm，平均孔隙半径在35μm，整体处于小孔隙范围内，为典型的致密储层。

常规铸体薄片由于放大倍数有限只能对较大的孔隙以及喉道类型、大小及产状进行定性到半定量表征，但是很难对致密储层的微小孔隙和喉道进行准确地测量。因此结合扫描电镜能够更加准确的识别微观孔候特征和大小以及组合关系。

通过对样品铸体薄片观察和扫描电镜图像分析，明确储层孔隙类型，厘定主要的孔隙类型，明确孔隙半径范围及平均孔隙半径。

(1)粒间孔

研究区长7段储层粒间孔主要为原生残余粒间孔和混合成因粒间孔，根据孔隙半径大小可以细分为大粒间孔和小粒间孔。粒间孔在研究区发育最为普遍，但是孔隙分布不均一，非均质性很强。原生残余粒间孔是指在埋藏过程中受到压实作用和胶结作用的影响未被充填的孔隙，形态呈倒三角形。混合成因粒间孔大多数都与溶蚀作用有关，碎屑颗粒边缘接触流体发生溶蚀或者胶结物溶蚀生成一些孔隙。粒间孔隙半径相对较大，一般介于15μm～115μm，连通性较好容易形成有利储集空间(如图2中的a、b所示)。

(2)粒内孔

研究区粒内孔主要是指粒内溶孔，由于长石、岩屑和少量云母等塑性成分容易溶蚀，导致颗粒内部形成次生溶孔，孔径主要在5μm～60μm。而发生溶蚀的碎屑颗粒边缘也会由于流体的经过导致局部部位溶蚀而呈不规则形状。粒内溶孔发育程度和分布特征与碎屑成分相关，导致其分布不均，部分长石颗粒可完全溶蚀形成铸模孔(如图2中的c、d所示)。

(3)微孔隙

微孔隙发育在填隙物内部，其孔径较小，大多介于2μm～10μm，在普通电子显微镜下难以识别，通常利用扫描电镜识别。微孔隙大多成群分布形成微孔隙群，孔隙与喉道之间连通性较差(如图2中的e所示)。

(4)微裂缝

研究区微裂缝发育较少，在局部地区发育，微裂缝宽度一般小于10μm，与周围孔隙连通性较好(如图2中的f所示)。

对研究区孔隙类型进行统计发现混合粒间孔含量最高，原生残余粒间孔和粒内溶孔次之，微孔隙和微裂缝含量较低，如图3所示。图3中纵坐标的频率也就是含量的意思。

S20：确定所述多个岩石样品薄片中的喉道类型。

喉道的大小和形态影响着流体在孔隙系统中的流动，从而影响岩石的渗透率，进而影响储集层的储集性能。通过铸体薄片和扫描电镜分析，明确储层发育的喉道类型及主要的喉道类型。

喉道的大小和形态影响着流体在孔隙系统中的流动，从而影响岩石的渗透率，进而影响储集层的储集性能。通过铸体薄片和扫描电镜发现，研究区长7段储层主要发育片状喉道，管束状喉道次之，缩颈型喉道发育较少。

片状喉道在形态上可以表现为片状或者片弯状(如图4中的a、c所示)，根据喉道的宽度可分为宽片状喉道、窄片状喉道、极窄片状喉道。此种喉道多与压实作用有关，碎屑颗粒在受到较强的压力作用下会呈一定的方向排列在储层中而形成片状喉道。

管束状喉道(如图4中的d所示)主要发育在胶结物和杂基中，在埋藏过程中原生孔隙被完全充填而导致孔隙损失殆尽，而杂基和胶结物中的微孔隙既可以是孔隙也可以是喉道，这些微孔隙群通过细小喉道相互连接，这些喉道共同组成管束状喉道。

缩颈型喉道(如图4中的b所示)一般发育在溶蚀作用和压实作用相对较强部位，颗粒之间收缩形成。当颗粒在沉积过程中受到压实作用而使得颗粒紧密接触，而颗粒间孔隙可以部分保留，但是喉道却被压缩变窄。此类喉道多见于颗粒之间点接触和线接触的储集层中。

S30：确定所述多个岩石样品的孔隙度和渗透率。

孔隙度和渗透率是通过现有方法可以得到的，这里不再赘述。

S40：在对所述多个岩石样品分别进行高压压汞实验和恒速压汞实验时，记录实验数据。

S50：根据所述实验数据确定岩石样品中的孔喉大小及分布特征。

(1)高压压汞表征孔喉大小及分布特征

通过分析致密砂岩样品高压压汞进汞及退汞曲线形态以及排驱压力、各参数特征，可将压汞曲线划分为不同类型。通过高压压汞数据，绘制孔喉半径分布范围，分析样品孔喉半径呈现出分布曲线峰值特征，明确不同类型的孔喉半径分布区间及峰值。

通过分析研究区长7段致密砂岩样品高压压汞曲线形态以及各参数特征，可将压汞曲线划分为4类，为I型、II型、III型和IV型。从I型到IV型排驱压力(即图5中汞开始大量进入汞时的压力)平均值显著增大，从0.67Mpa到5.37Mpa。而最大孔喉半径呈减小趋势，从1.097μm减小至0.267μm。I型最大进汞饱和度最大为95.54％，孔喉半径中值平均为0.263μm，而IV型最大进汞饱和度最小为48.57％，孔喉半径中值平均为0.03μm。图5示出了储层致密砂岩高压压汞毛细血管压力图的集合，图6分别示出了I型到IV型的高压压汞毛细血管压力图。

基于压力数据与孔喉的对应关系、高压压汞数据，绘制孔喉半径分布范围，可以直观的看出研究区长7段样品孔喉半径呈现出单峰分布特点(如图7所示)。半径主要分布在0.0036μm～1.566μm。I型曲线孔喉半径分布在0.267μm～1.566μm，峰值为0.541μm；II型曲线孔喉半径分布在0.015μm～0.541μm，峰值为0.134μm；III型曲线孔喉半径分布在0.009μm～0.179μm，峰值为0.067μm；IV型曲线孔喉半径分布在0.0036μm～0.134μm，峰值为0.036μm。与常规铸体薄片下所分析测量的孔喉大小相比，高压压汞能够更准确更精细的识别孔喉分布情况。但是由于最大进汞饱和度的影响，部分较小的孔隙难以在外界压力的影响下注入汞从而导致半径较小的孔喉高压压汞无法进行表征。

(2)恒速压汞表征孔喉大小及分布特征

分析致密砂岩样品恒速压汞曲线特征和排驱压力、总喉道进汞饱和度、最终总进汞饱和度、喉道半径平均值、孔隙度、渗透率、最大连通喉道半径等特征。

通过分析研究区长7段致密砂岩样品恒速压汞曲线特征和相关孔隙结构参数，对比发现研究区恒速压汞进汞曲线有两种类型。第一种类型表现为排驱压力较低，总喉道进汞饱和度和最终总进汞饱和度较高，以及喉道半径平均值较高(如图8所示)；第二种类型表现为排驱压力较高，总喉道进汞饱和度和最终总进汞饱和度较低，以及喉道半径平均值较低(如图9所示)。类型I样品孔隙度10.713％，渗透率0.087mD，最大连通喉道半径0.876μm，当进汞压力在1Mpa及以下时，喉道进汞饱和度和最终总进汞饱和度曲线趋势基本一致，但随着压力增加开始分离。最终总进汞饱和度64.103％，总喉道进汞饱和度38.483％。类型II样品孔隙度6.28％，渗透率0.013mD，最大连通喉道半径0.365μm，当进汞压力在3Mpa及以下时，喉道进汞饱和度和最终总进汞饱和度曲线趋势基本一致，但随着压力增加开始分离。最终总进汞饱和度32.357％，总喉道进汞饱和度29.382％。通过恒速压汞实验数据分析可以得出不同的砂岩样品孔渗不同，其孔隙半径和喉道半径也有较大差别(如图10和图11所示)。物性较好的样品，孔隙半径主要分布在130μm～260μm，平均半径约为196.6μm，喉道范围在0.2μm～0.7μm，平均喉道半径0.525μm，最大连通喉道半径0.876μm；物性较差的样品，孔隙半径主要分布在80μm～160μm，平均半径约为124.2μm，喉道范围在0.04μm～0.22μm，平均喉道半径0.18μm，最大连通喉道半径0.365μm。研究表明西233区块长7段致密砂岩储层的物性主要受孔隙、喉道半径以及最大连通喉道半径的影响。

S60：确定所述多个岩石样品薄片中的孔隙类型与喉道类型的主要组合配置关系。

孔隙及喉道的组合配置关系不同会形成不同的孔隙结构特征。根据孔隙和喉道特征，结合铸体薄片、扫描电镜、高压压汞等资料对储层孔隙-喉道配置关系进行分类。例如：大粒间孔—宽片状喉道组合、小-中粒间孔—窄片状喉道组合、粒内孔—窄片状喉道组合、微孔—极窄片状、管束状喉道组合。

研究表明致密砂岩孔隙结构相对于常规砂岩更复杂，在同一深度段同一块岩石样品中也会存在多种孔隙和喉道类型，他们之间的组合配置关系(即主要连接关系)不同就会形成不同的孔隙结构特征。根据孔隙和喉道特征，结合铸体薄片、扫描电镜、高压压汞等资料对研究区长7段储层孔喉组合关系进行分类(表1)。主要划分为4种类型，分别是大粒间孔—宽片状喉道组合、小-中粒间孔—窄片状喉道组合、粒内孔—窄片状喉道组合、微孔—极窄片状、管束状喉道组合。

表1西233区块长7段致密砂岩孔喉组合关系解释标准

大粒间孔—宽片状喉道组合主要表现为孔隙粒径较大，以原生残余粒间孔和粒内溶孔为主，大粒间孔相对含量大于50％，粒内溶孔含量大于30％，孔隙半径介于65μm～115μm。喉道类型主要为宽片状喉道，其含量大于50％，窄片状、管束状喉道较少总和小于50％，偶见缩颈型喉道，喉道半径为0.5μm～2μm。此种孔喉组合类型常发育在块状细砂岩相和交错层理细砂岩相、平行层理细砂岩相中，储集物性好(如图12中的a所示)。

小-中粒间孔—窄片状喉道组合表现为孔隙粒径较小，以小-中原生残余粒间孔为主，含量大于50％，粒内溶孔次之，含量大于20％，大孔径粒间孔和微孔隙含量较少，孔隙半径介于25μm～65μm。喉道类型主要为窄片状喉道，其含量大于50％，宽片状和管束状喉道含量相对较少，喉道半径为0.3μm～1μm。此种孔喉组合类型常发育在含变形构造粉砂岩相、波纹层理粉砂岩相中，平行层理细砂岩相和块状细砂岩相也有少量发育，储集物性相对较好(如图12中的b所示)。

粒内孔—窄片状喉道组合中以粒内溶孔为主含量大于50％，其次是小-中原生残余粒间孔和微孔隙群，含量分别小于20％，孔隙半径较小介于5μm～25μm。喉道类型主要为窄片状喉道，其含量大于50％，其次为管束状喉道，含量相对小于30％，宽片状喉道和颈型喉道较少，喉道平均半径小于0.5μm。此种孔喉组合类型常发育在含泥砾粉砂岩相、波纹层理粉砂岩相、含泥质条带粉砂岩相中，块状细砂岩相也有少量发育，储集物性中等(如图12中的c所示)。

微孔—极窄片状、管束状喉道组合中的孔隙以微孔隙为主，含量高于50％，微孔隙主要为发育在黏土矿物和泥质胶结物中的孔隙，微孔隙半径小于5μm。喉道类型主要为窄片状喉道和管束状喉道以及部分缩颈型喉道，喉道平均半径小于0.5μm。此种孔喉组合类型常发育在水平层理泥质粉砂岩相、含植物碎屑泥质粉砂岩相中，含泥质条带粉砂岩相也有少量发育，储集物性较差(如图12中的d所示)。

可以将不同孔隙-喉道配置关系孔隙度及渗透率投点到孔渗交汇图中，分析不同孔隙-喉道配置关系样品的孔隙度及渗透率特征，明确不同孔隙-喉道配置关系的物性差异。

S70：对所述实验数据进行处理，得到全孔径下的分形维数。

基于分形理论，利用高压压汞和恒速压汞所得的进汞饱和度对数Lg(SHg)与毛管压力对数Lg(Pc)之间的分段关系分析孔隙结构的分型维数。并且由于高压压汞和恒速压汞进汞过程相对一致，所以可以结合这两种实验方法来获得致密砂岩全孔径分形维数。

致密砂岩孔隙结构分形研究大多数都是停留在单一压汞技术分析上，但是高压压汞与恒速压汞都存在自身的局限性，难以很好对致密砂岩孔隙结构进行准确表征。然而恒速压汞虽然能够表征孔喉的大小以及分布特征，但是由于最大进汞压力的限制，无法对0.12μm以下的孔喉进行表征。高压压汞能对0.12μm以下孔喉进行测量，但无法测定与窄片状喉道相连的大孔分布，而恒速压汞则可以满足，因此综合这两种方法来对全孔喉特征进行表征。同时由于高压压汞和恒速压汞进汞过程相对一致，所以结合这两种实验方法来获得研究区致密砂岩全孔径分形维数是有效的。

基于分形理论综合应用到本研究区致密砂岩孔隙结构表征，利用高压压汞和恒速压汞所得的进汞饱和度对数Lg(SHg)与毛管压力对数Lg(Pc)之间的分段关系分析孔隙结构的分型维数。同时综合前人相关计算模型对本研究区推导出如下关系式：

Lg(SHg)＝ (D- 2)Lg(Pc)+ C (1)

其中，D为分形维数；S

(1)高压压汞多重分形维数

根据致密砂岩样品高压压汞实验汞饱和度与毛细管压力数据分析高压压汞分形拟合曲线特征，明确分形拟合曲线可划分为几个段，以及不同分形区间所对应的孔隙、喉道类型，厘定分形维数数值范围，以及拟合曲线相关系数。

高压压汞分形拟合曲线往往表现出曲线具有四段分布的特征，分别可以对应大孔隙、宽喉道、小孔隙、窄喉道。各样品拟合曲线相关系数R较大，相关性较好，大于0.80，表明每段划分具有较强的代表性和准确性且具有多维分形特征。

(2)恒速压汞的多重分形维数

根据致密砂岩样品恒速压汞实验汞饱和度与毛细管压力数据分析恒速压汞分形拟合曲线特征，明确分形拟合曲线可划分为几个段，以及不同分形区间所对应的孔隙、喉道类型，厘定分形维数数值范围，以及拟合曲线相关系数。

根据研究区致密砂岩样品高压压汞和恒速压汞的进汞饱和度与毛细管压力数据进行分析。发现高压压汞分形拟合曲线表现出曲线具有四段分布的特征，分别可以对应大孔隙、宽喉道、小孔隙、窄喉道(如图13所示)，各样品拟合曲线相关系数R较大，相关性较好，都分布在0.85～0.99之间，表明每段划分具有较强的代表性和准确性。因此可以证实研究区样品的孔隙结构比较复杂具有多维分形特征。如下表2所示。

在表2中，Dg1表示反应大孔隙的拟合曲线的分形维数，Dg1右侧的R

表2研究区长7段致密砂岩样品高压压汞分形维数

根据研究区长7段致密砂岩样品恒速压汞分形维数拟合曲线特点可以发现，拟合曲线主要为3段、4段式。但是不管3段还是4段式都能有效的区分大孔隙和宽喉道，其差异性是3段式小孔隙和窄喉道在图上分布在一起很难进行区分，而4段式则可以将小孔隙和窄喉道进行有效区分(如图14所示)。与高压压汞相比恒速压汞在表征大孔隙和宽喉道时相关系数R较大，相关性较好，因此能够更有效地对致密砂岩储层中较大的孔隙与喉道进行识别(如下表3所示)。

在表3中，Dh1表示反应大孔隙的拟合曲线的分形维数，Dh1右侧的R

表3研究区长7段致密砂岩样品恒速压汞分形维数

致密砂岩储层孔喉大小和分布可以通过高压压汞和恒速压汞实验来定量表征，高压压汞能够表征大孔隙，但是对小孔隙和喉道区分度有限难以表征。而恒速压汞进汞速度恒定且较低，能够通过压力变化来辨别出孔隙和喉道。所以结合高压压汞和恒速压汞共同表征孔隙和喉道的大小与分布特征。

(3)高压压汞与恒速压汞相结合的全孔径分形维数

通过观察高压压汞及恒速压汞的分形维数拟合曲线特征，分别明确大孔隙分形维数(Dg1)、宽喉道分形维数(Dg2)、小孔隙分形维数(Dg3)及窄喉道分形维数(Dg4)的数值范围。因致密砂岩孔隙结构分形维数一般介于2～3，对于分形维数不在此范围内的孔隙或喉道类型予以去除，明确具有分析特征的孔隙及喉道类型。

对具有分形特征的孔隙、喉道的分形维数与孔隙度及渗透率进行相关性分析，明确与储层孔渗相关较好的分形维数，确定不同的分形维数段对储层物性的影响。

根据实验原理，高压压汞对小孔喉表征较好，而恒速压汞对大孔喉表征更为准确，通过综合对比高压压汞与恒速压汞分型维数特征以及分型维数与储层物性相关性，将两者进行结合表征全孔径孔隙结构。

根据高压压汞与恒速压汞大孔隙、宽喉道、小孔隙、窄喉道分型维数交点处的半径分别对应rg1、rg2、rg3和rh1、rh2、rh3，将高压压汞与恒速压汞进行拼接，拼接值选用高压压汞与恒速压汞所测定的小孔喉rg2与大孔喉rh2的临界处平均值，进而来表征全孔径孔隙结构特征。

通过观察高压压汞分型维数拟合曲线特征，大孔隙分形维数(Dg1)范围在2.26～2.64；宽喉道分形维数(Dg2)范围在2.04～2.11；小孔隙分形维数(Dg3)范围在3.32～3.46；窄喉道分形维数(Dg4)范围在2.04～2.16。由于致密砂岩孔隙结构分形维数一般介于2～3具有分型特征，否则不具有分型特征。因此研究区长7段致密砂岩储层小孔隙不具分型特征，而大孔隙、宽喉道、窄喉道具有分形特征。针对于窄喉道分形维数与孔隙度相关性较好(R为-0.754)，表明Dg4对储层孔渗影响显著(如图15所示)。分析恒速压汞分型维数拟合曲线特征，大孔隙分形维数(Dh1)范围在2.60～2.74；宽喉道分形维数(Dh2)范围在2.10～2.29；小孔隙分形维数(Dh3)范围在3.68～4.30；窄喉道分形维数(Dh4)范围在2.29～2.50。同样小孔隙不具分型特征，而大孔隙、宽喉道、窄喉道具有分形特征。且对于恒速压汞宽喉道分形维数与孔隙度相关性较好(R为0.908)，因此Dh2可以准确反映致密砂岩储层宽喉道特征(如图16所示)。

通过综合对比高压压汞与恒速压汞分型维数特征以及分型维数与储层物性相关性，可以看出高压压汞对小孔喉表征较好，而恒速压汞对大孔喉表征更为准确，因此将两者进行结合表征全孔径孔隙结构(表4)。根据高压压汞与恒速压汞大孔隙、宽喉道、小孔隙、窄喉道分型维数交点处的半径分别对应r

在表4中，r

表4研究区长7段致密砂岩样品高压压汞与恒速压汞分形维数交点处孔隙半径

由此可见，S70可以包括如下步骤S71至S76。

S71：确定高压压汞和恒速压汞实验的各压力数据对应的孔喉大小及孔喉类型。

S72：将高压压汞实验的实验数据中的压力数据、汞饱和度数据分别取对数，构造第一数据集，所述第一数据集中的实验数据按照孔喉大小及类型划分至以下四个类别：大孔隙、宽喉道、小孔隙、窄喉道。

S73：对于所述第一数据集中的各个类别的实验数据，分别确定实验数据的拟合直线，得到高压压汞实验对应的多条拟合直线；根据拟合直线的斜率确定反映对应孔喉类别的分形维数。

S74：将恒压压汞实验的实验数据中的压力数据、汞饱和度数据分别取对数，构造第二数据集，将所述第二数据集中的数据按照孔喉大小及类别划分至以下四个类别：大孔隙、宽喉道、小孔隙、窄喉道。

S75：对于所述第二数据集中的各个类别的实验数据，分别确定实验数据的拟合直线，得到恒压压汞实验对应的多条拟合直线；根据拟合直线的斜率确定反映对应孔喉类别的分形维数。

S76：将所述高压压汞实验对应的多条拟合直线与所述恒压压汞实验对应的多条拟合直线进行拼接，得到全孔径下的分形维数。

可以将所述高压压汞实验对应的多条拟合直线与所述恒压压汞实验对应的多条拟合直线进行拼接，拼接值采用所述高压压汞实验和恒速压汞实验对应的多条拟合直线中小孔喉和大孔喉分界点rg2及rh2的平均值。

在一些实施例中，依据上述公式(1)，可以根据拟合直线的斜率确定反映对应孔喉类别的分形维数。具体地，可以是先计算目标拟合直线的斜率与2的和值(即将斜率值加2)，然后将所述和值作为反映目标孔喉类别的分形维数，所述目标孔喉类别是所述目标拟合直线关联的实验数据所对应的孔喉类别。

S80：根据孔隙类型与喉道类型的主要组合配置关系、孔隙度、渗透率、全孔径孔喉的分形维数，确定致密砂岩储层孔隙的微观拓扑结构特征表，所述微观拓扑结构特征表中包括致密岩石的多个微观拓扑结构类别，所述微观拓扑结构特征表用于确定目标地区的油气勘探或开发方案。

致密储层基质中的储集空间是储层油气的场所，并构成了油气在储层中流动的路径最末端的部分。基质储层孔隙作为油气储集空间的功能是相似的。但是基质储层中孔隙喉道的发育模式不同会大大影响储层绝对渗透率，并且导致油气的采出程度的巨大差异。

基质储层的孔隙结构是影响油气采收率的重要因素之一。通过退汞后残余汞包含度S

致密砂岩相对于常规砂岩孔隙结构更为复杂，本发明提出在研究其孔隙结构时通过对孔隙类型、喉道类型、孔喉配置关系、压汞曲线特征、多重分形维数等综合研究，建立基于拓扑结构的多类型孔隙结构模型(也即微观拓扑结构类别)来表征不同的孔隙结构特征。不同孔隙结构模型(也即微观拓扑结构类别)且会导致储层质量及油气的最终采收率差异。可以通过该模式进行储层分类评价，并对预测致密油气的产能及最终采收率也具有重要的意义。

本发明基于根据不同的孔隙类型、喉道类型、孔喉配置关系、高压压汞曲线特征、恒速压汞曲线特征、分形维数特征等，建立了四种基于拓扑结构的孔隙结构模式：球棍网状结构、星链网状结构、离散树状结构、复杂管束状结构。

(1)球棍网状结构

图18中的a示出了球棍网络结构的示意图。图20中示出了球棍网络结构对应的特征。

球棍网状结构孔隙以大粒间孔为主，含有少量粒内溶孔，喉道类型主要为宽片状喉道，缩颈型喉道和管束状喉道基本不发育。孔喉配置关系为大粒间孔—宽片状喉道组合相对应(如图19中的a所示)。压汞曲线特征表现为门槛压力较低小于0.9Mpa，平台明显，粗歪度，分选较好。在球棍孔隙结构模型(也即微观拓扑结构类别)中存在多个循环路径，也就是说，从一个孔隙开始，有一条或多条路径会返回到该孔隙，而不会有重新回溯的路径。

粒间孔隙网络，在两个孔隙之间存在很多通道。流体的驱替表现出经典渗流特征，即在较小的毛细管压力变化会导致较大的汞饱和度的变化。

在压汞实验进汞过程中，当压力到达某个值，汞会大量进入样品中喉道相互连接的储集空间，形成整体排驱现象，压汞曲线呈现出一个明显的平台。

在退汞过程中，首先是最小孔喉中的汞优选排出，排出到周围较大的喉道中。然后随着压力的降低，与孔隙相连的喉道中的汞也逐步清空，将汞排到孔隙中。与外部通道不连通的孔隙中的汞便会残留在样品中无法排出，导致退汞后，具有较高的残余汞含量。

反映大孔、宽喉道，且与孔隙度、渗透率呈正相关的分形维数段的分形维数值高。此种类型孔隙结构储层物性较好，孔隙度及渗透率值相对高。

球棍网状结构孔隙对应的压汞曲线为Ⅰ型，门槛压力较低，平台比较明显，歪度较粗，分选较好。反映大孔的分形维数D

(2)星链网状结构

图18中的b示出了星链网络结构的示意图。图20中示出了星链网络结构对应的特征。

星链网状结构孔隙以小-中粒间孔为主，大粒间孔较少，还含有部分粒内溶孔，喉道类型主要为窄片状喉道，缩颈型喉道和管束状喉道较少。孔喉配置关系主要为小-中粒间孔—窄片状喉道组合。压汞曲线特征表现为门槛压力略低一般小于1.5Mpa，压汞曲线具有局部的平台特征，歪度略粗，分选较中等-好，与小-中粒间孔—窄片状喉道组合相对应(如图19中的b所示)。

星链网络结构对应的压汞曲线为II型，门槛压力较低，平台相对明显，歪度略粗，分选中等-好。反映小孔的分形维数D

在该孔隙结构模型(也即微观拓扑结构类别)中，既存在孔隙连通性较好、具有多个循环路径的粒间孔隙，也存在部分孔隙连通性较差、不具有循环孔喉通道的粒内孔隙。

在压汞实验进汞过程中，在中等压力下有局部的“排驱”现象，即压汞曲线具有局部的平台特征。样品中粒间孔隙网络部分，会发生这种在某个窄的压力下大量汞进入的现象。然后在高压下，毛管压力与汞饱和度呈指数关系，这种毛管压力曲线代表汞进入粒内微孔部分。粒间孔所占的比例越大，毛管压力曲线中的“平台”特征越发育。

反映大孔、宽喉道，且与孔隙度、渗透率呈正相关的分形维数段的分形维数值相对较高。此种类型孔隙结构储层物性较好，孔隙度及渗透率值相对较高。

(3)离散树状结构

图18中的c示出了离散树状结构的示意图。图20中示出了离散树状结构对应的特征。

离散树状结构孔隙以粒内孔为主，含有少量小粒间孔，喉道类型主要为窄片状喉道，缩颈型喉道次之。压汞曲线特征表现为门槛压力增大至3.2Mpa，平台不太明显，歪度略偏细，分选一般，与粒内孔—窄片状喉道组合相对应(如图19中的c所示)。

如图20所示，离散树状结构以粒内孔和小粒间孔为主，喉道主要为窄片状喉道，缩颈型喉道次之。对应的压汞曲线为Ⅲ型，门槛压力较大，平台不太明显，歪度略偏细，分选一般。反映小孔的分形维数D

离散树状结构孔隙以粒内孔为主，含有少量小粒间孔，喉道类型主要为窄片状喉道，缩颈型喉道次之。孔喉配置关系主要为粒内孔—窄片状喉道组合。压汞曲线特征表现为门槛压力增大至3.2Mpa，压汞曲线不具有平台，汞饱和度与压力呈指数关系，孔喉分选差。

离散树形孔隙结构中主要为粒内孔隙，具有独特的空间分布特征。孔隙结构呈树形分叉，可分出更加细小的喉道。与球棍网状结构相比，树形孔隙结构不具有通道的循环。也就是说，一个孔隙通过树形孔隙结构网络通道不等再回到该孔隙，除非回溯过去。没有通道循环意味着在进行压汞实验的时候，汞到达每个孔隙只有一条途径。而流体若要进较窄的喉道只能从较宽的喉道进入。这些特征导致在进汞过程中，随着压力的增大，不会出现在某一压力值下大量汞进入样品的特征，也就是毛细管曲线没有明显的平台。相反，随着毛管压力的不断增加，导致进汞饱和度的不断增加。树形孔隙网络模下，毛管压力随着进汞饱和的增加呈指数形式增加。

随着毛管压力的降低，汞首先从最窄的喉道排出，然后随着毛细压力的不断降低，较宽的喉道中的汞也逐步排出。因为树形孔隙结构中，孔喉中没有循环的路径，大孔喉不断分叉成小孔喉，理论上树形孔隙结构中的汞都能够被排出，就是进汞曲线与退汞曲线会保持一致，所以离散树状结构的样品往往具有较高的退汞效率。这种现象对应油气的开采是有利的，因为这种孔隙结构具有较高的采收率。

反映小孔、窄喉道，且与孔隙度、渗透率呈负相关的分形维数段的分形维数值较高。此种类型孔隙结构储层物性较差，孔隙度及渗透率值相对较低。

(4)复杂管束状结构

图18中的d示出了离散树状结构的示意图。图20中示出了复杂管束状结构对应的特征。

复杂管束状结构孔隙以微孔为主，含有少量小粒内溶孔，喉道类型主要为极窄片状喉道和管束状喉道，缩颈型喉道较少。孔喉配置关系主要为微孔—极窄片状/管束状喉道组合。压汞曲线特征表现为门槛压力较大约为9.6Mpa，平台不明显，歪度较细，分选中等-差，与微孔—极窄片状、管束状喉道组合相对应(如图19中的d所示)。

复杂管束状结构，其储集空间主要由微小孔隙构成，汞要进入则需要较大的压力，所以排驱压力较大。随着压力进一步增大，汞可以进入更小的孔隙，压力曲线呈上拱型，具有细歪度。并且由于孔隙小，喉道复杂，汞无法大量的进入他们。所以同等压力下，汞饱和度要低于其他孔隙结构模型(也即微观拓扑结构类别)的样品。

反映小孔、窄喉道，且与孔隙度、渗透率呈负相关的分形维数段的分形维数高。此种类型孔隙结构储层物性差，孔隙度及渗透率值低。IV类孔隙结构以微孔为主，含有少量小粒内溶孔，喉道主要为极窄片状喉道和管束状喉道，缩颈型喉道较少。对应的压汞曲线为IV型，门槛压力较大，平台不明显，歪度较细，分选中等-差。反映极窄喉道的分形维数D

相应地，步骤S80可以包括如下步骤S81、S82、S83和S84。

S81：在孔隙类型与喉道类型的主要组合配置关系为大粒间孔-宽片状喉道组合、孔隙度大于10％、渗透率大于0.1mD、反映大孔的分形维数分布在2.26～2.4、反映宽喉道的分形维数分布在2.04～2.11的情况下，确定岩石样品的微观拓扑结构类别为球棍网状结构。

S82：在孔隙类型与喉道类型的主要组合配置关系为小中粒间孔-窄片状喉道组合、孔隙度为8％～12％，渗透率为0.05mD～0.1mD、反映小孔的分形维数分布在2.46～2.81、反映窄喉道的分形维数分布在2.04～2.16的情况下，确定岩石样品的微观拓扑结构类别为星链网状结构。

S83：在孔隙类型与喉道类型的主要组合配置关系为粒内孔—窄片状喉道组合、孔隙度为5％～10％、渗透率为0.02mD～0.06mD、反映小孔的分形维数分布在2.32～2.81、反映窄喉道的分形维数分布在2.04～2.16的情况下，确定岩石样品的微观拓扑结构类别为离散树状结构。

S84：在孔隙类型与喉道类型的主要组合配置关系为微孔—极窄片状及管束状喉道组合、孔隙度为5％、渗透率小于0.03mD、反映极窄喉道的分形维数分布在2.08～2.16的情况下，确定岩石样品的微观拓扑结构类别为复杂管束状结构。

在得到图20所示的微观拓扑结构特征表之后，还可以通过以下方法确定目标地区致密砂岩储层孔隙的微观拓扑结构类别：

S91：确定所述目标地区岩石样品薄片中的孔隙类型、孔径大小。

S92：确定所述目标地区岩石样品薄片中的喉道类型。

S93：确定所述目标地区岩石样品的孔隙度和渗透率。

S94：在对所述目标地区岩石样品分别进行高压压汞实验和恒速压汞实验时，记录实验数据。

S95：根据所述实验数据确定岩石样品中的孔喉大小及分布特征。

S96：确定所述目标地区岩石样品薄片中的孔隙类型与喉道类型的主要组合配置关系。

S97：根据实验数据，处理得到全孔径孔喉的分形维数。

S98：将所述目标地区岩石样品的孔隙类型与喉道类型的主要组合配置关系、孔隙度、渗透率、全孔径孔喉的分形维数与致密砂岩储层孔隙的微观拓扑结构特征表进行比对，将所述微观拓扑结构特征表中与所述目标地区岩石样品的各特征最相近的目标微观拓扑结构类别作为所述目标地区岩石样品的微观拓扑结构类别，以根据所述目标地区岩石样品的微观拓扑结构类别确定油气勘探或开发方案。

本说明书提供的微观拓扑结构特征表，确定了影响油气产能及最终采收率的微观结构特征，根据该微观拓扑结构特征表能够从微观结构的角度更为直观、准确地预测油气的产能及最终采收率。现有的预测方法通常是直接根据某几种测试数据(例如一种测试数据为孔隙度)进行预测的，由于地层的结构的复杂特性，个别种类的测试数据难以准确地表征地层的油气采收率及难度，若对于测试数据的处理不合理、或者采用的测试数据种类不足则会影响预测结果。本方案直接根据指标确定微观结构特征，并基于微观结构特征表来进行预测，预测结果更为准确。

在球棍孔隙结构模型中存在多个循环路径，也就是说，从一个孔隙开始，有一条或多条路径会返回到该孔隙，而不会有重新回溯的路径。在油气采出过程中，首先是最小孔喉中的油气优选排出，排出到周围较大的喉道中。然后排到孔隙中。与外部通道不连通的孔隙中的油气便会残留在样品中无法排出，导致具有较高的残余油气含量。该类型孔隙结构的储层质量较好，孔隙、喉道较大，油气饱和度较高，但油气的最终采收率有限。

离散树形孔隙结构中主要为粒内孔隙，具有独特的空间分布特征。孔隙结构呈树形分叉，可分出更加细小的喉道。与球棍网状结构相比，树形孔隙结构不具有通道的循环。也就是说，一个孔隙通过树形孔隙结构网络通道不能再回到该孔隙，除非回溯过去。没有通道循环意味着在油气采出的时候，到达每个孔隙只有一条途径。流体从较窄的喉道到较宽的喉道再排出，理论上树形孔隙结构中的油气都能够被排出，所以离散树状结构的样品往往具有较高的退汞效率。这种现象对应油气的开采是有利的，因此这种孔隙结构具有较高的油气最终采收率。

星链网状结构模型中，既存在孔隙连通性较好、具有多个循环路径的粒间孔隙，也存在部分孔隙连通性较差、不具有循环孔喉通道的粒内孔隙。故储层质量及油气饱和度介于球棍孔隙结构模型与离散树形孔隙结构模型之间，油气最终采收率也介于球棍孔隙结构模型与离散树形孔隙结构模型之间。

复杂管束状结构模型由于孔隙小，喉道复杂，油气难以大量的进入孔隙及喉道。油气进去后，开采过程中由难以采出，所以储层质量及最终油气采收率都较差。

本说明书还提供一种致密砂岩储层孔隙的微观拓扑结构确定装置，可以用于实现上述致密砂岩储层孔隙的微观拓扑结构确定方法。该装置包括第一确定单元、第二确定单元、第三确定单元、记录单元、第四确定单元、第五确定单元、处理单元和第六确定单元。

第一确定单元用于确定致密储层的多个岩石样品薄片中的孔隙类型、孔径大小。

第二确定单元用于确定所述多个岩石样品薄片中的喉道类型。

第三确定单元用于确定所述多个岩石样品的孔隙度和渗透率。

记录单元用于在对所述多个岩石样品分别进行高压压汞实验和恒速压汞实验时记录实验数据。

第四确定单元用于根据所述实验数据确定岩石样品中的孔喉大小及分布特征。

第五确定单元，用于确定所述多个岩石样品薄片中的孔隙类型与喉道类型的主要组合配置关系。

处理单元用于对所述实验数据进行处理，得到全孔径下的分形维数。

第六确定单元用于根据孔隙类型与喉道类型的主要组合配置关系、孔隙度、渗透率、全孔径孔喉的分形维数，确定致密砂岩储层孔隙的微观拓扑结构特征表，所述微观拓扑结构特征表中包括致密岩石的多个微观拓扑结构类别，所述微观拓扑结构特征表用于确定目标地区的油气勘探或开发方案。

上述装置的具体描述及有益效果请参阅对应的方法实施例，不再赘述。

本发明实施例还提供了一种电子设备，该电子设备可以包括处理器和存储器，其中处理器和存储器可以通过总线或者其他方式连接，例如通过总线连接。

处理器可以为中央处理器(Central Processing Unit，CPU)。处理器还可以为其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片，或者上述各类芯片的组合。

存储器作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态计算机可执行程序以及模块，如本发明实施例中的致密砂岩储层孔隙的微观拓扑结构确定方法对应的程序指令/模块(例如，第一确定单元、第二确定单元、第三确定单元、记录单元、第四确定单元、第五确定单元、处理单元和第六确定单元)。处理器通过运行存储在存储器中的非暂态软件程序、指令以及模块，从而执行处理器的各种功能应用以及数据分类，即实现上述方法实施例中的致密砂岩储层孔隙的微观拓扑结构确定方法。

存储器可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储处理器所创建的数据等。此外，存储器可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中，存储器可选包括相对于处理器远程设置的存储器，这些远程存储器可以通过网络连接至处理器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。

所述一个或者多个模块存储在所述存储器中，当被所述处理器执行时，执行上述致密砂岩储层孔隙的微观拓扑结构确定方法。

上述电子设备具体细节可以参阅对应实施例中的相关描述和效果进行理解，此处不再赘述。

上述实施方式阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。

为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。

以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。