基于致密储层的微裂缝有效开度确定方法及装置

技术领域

本发明涉及石油勘探领域，尤其是致密油田的开发技术领域，具体涉及一种基于致密储层的微裂缝有效开度确定方法及装置。

背景技术

针对题述技术内容，申请号为201310706032.9(一种测定致密油临界充注孔喉半径门限的方法及装置)，以及申请号为201410713670.8(一种致密油气储层物性下限的确定方法)，公布了一下内容：本发明涉及一种致密油气储层物性下限的确定方法，应用于致密油气储层物性下限的确定。首先建立平均毛细管压力的求取模型，其次建立最小流动孔喉半径的求取模型，最后求储层物性下限。以及申请号为201710249629.3(致密油藏有效储层物性下限值的确定方法和装置)，本发明提供了一种致密油藏有效储层物性下限值的确定方法和装置，提出了一种定量确定所述临界孔喉半径的计算方式，从而可量得到所述待测深度处的物性下限值。申请号为CN201410855862.2(致密油核磁共振T2截止值及流体饱和度确定方法、装置)，对待测样品进行碳含量测定，生成碳含量测定结果；根据碳含量测定结果确定待测样品的孔喉半径下限值；根据孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值。本方案利用实际样品测定的致密油可动流体孔喉半径下限值求取T2截止值，由于实际样品测定的致密油可动流体孔喉半径下限值采用由粗变细孔喉测定碳含量逐渐逼近的方法，更加具有真实性。

上述方法存在的问题主要有：(1)现有方法和技术均针对单一孔隙介质，对于含微裂缝的双重介质致密油储层不适应；(2)现有方法主要是确定单一孔隙介质致密油储层的孔喉半径下限，没有将孔隙与喉道分开。事实上致密油储层的孔隙与喉道存在明显差异，而且其各自表述的含义和对储层的影响也完全不同；(3)没有提出针对双重介质致密油储层有效微裂缝开度下限的确定方法，而微裂缝对致密油储层的渗透性和有效性起到极其重要的作用。

发明内容

针对现有技术中的问题，本发明的基于致密储层的微裂缝有效开度确定方法，针对致密储层微裂缝发育，孔隙喉道细小，且非均质性强的特点，将决定储层渗流能力的喉道半径和微裂缝开度作为有效性的关键评价参数。并以此为基础准确的计算出致密储层的微裂缝有效开度(下限)。

为解决上述技术问题，本发明提供以下技术方案：

第一方面，本发明提供一种基于致密储层的微裂缝有效开度确定方法，所述方法包括：

对所述致密储层的岩心进行岩心分析，以得到所述岩心不同喉道对渗透率的贡献率；

根据所述贡献率确定所述岩心的有效渗透率；

根据所述有效渗透率确定所述岩心的微裂缝有效开度。

一实施例中，所述对所述致密储层的岩心进行岩心分析，以得到所述岩心不同喉道对渗透率的贡献率，包括：

对所岩心进行压汞实验，以生成所述岩心的全尺度喉道分布曲线；

根据所述全尺度喉道分布曲线计算所述贡献率。

一实施例中，所述根据所述贡献率确定所述岩心的有效渗透率，包括：

对所述全尺度喉道分布曲线进行归一化处理，并计算平均喉道半径；

确定所述平均喉道半径与所述渗透率之间的关系式；

根据所述关系式确定所述岩心的有效喉道半径；

根据所述有效喉道半径确定所述有效渗透率。

一实施例中，所述根据所述有效渗透率确定所述岩心的微裂缝有效开度，包括：

根据所述有效渗透率确定间距最大的微裂缝所对应的有效渗透率；

根据所述间距最大的微裂缝所对应的有效渗透率确定所述微裂缝有效开度。

一实施例中，所述对所岩心进行压汞实验，以生成所述岩心的全尺度喉道分布曲线，包括：

对所述岩心进行恒速压汞实验，以生成恒速压汞曲线；

对所述岩心进行高压压汞实验，以生成高压压汞曲线；

对所述恒速压汞曲线与所述高压压汞曲线进行拼接，以生成所述全尺度喉道分布曲线。

第二方面，本发明提供一种基于致密储层的微裂缝有效开度确定装置，该装置包括：

岩心分析单元，用于对所述致密储层的岩心进行岩心分析，以得到所述岩心不同喉道对渗透率的贡献率；

有效渗透率确定单元，用于根据所述贡献率确定所述岩心的有效渗透率；

有效开度确定单元，用于根据所述有效渗透率确定所述岩心的微裂缝有效开度。

一实施例中，所述岩心分析单元包括：

分布曲线生成模块，用于对所岩心进行压汞实验，以生成所述岩心的全尺度喉道分布曲线；

贡献率计算模块，用于根据所述全尺度喉道分布曲线计算所述贡献率。

一实施例中，所述有效渗透率确定单元包括：

平均喉道半径计算模块，用于对所述全尺度喉道分布曲线进行归一化处理，并计算平均喉道半径；

关系式确定模块，用于确定所述平均喉道半径与所述渗透率之间的关系式；

有效喉道半径确定模块，用于根据所述关系式确定所述岩心的有效喉道半径；

有效渗透率确定模块，用于根据所述有效喉道半径确定所述有效渗透率。

一实施例中，所述有效开度确定单元包括：

有效渗透率确定模块，用于根据所述有效渗透率确定间距最大的微裂缝所对应的有效渗透率；

有效开度确定模块，用于根据所述间距最大的微裂缝所对应的有效渗透率确定所述微裂缝有效开度。

一实施例中，所述分布曲线生成模块包括：

恒速压汞曲线生成模块，用于对所述岩心进行恒速压汞实验，以生成恒速压汞曲线；

高压压汞曲线生成模块，用于对所述岩心进行高压压汞实验，以生成高压压汞曲线；

拼接模块，用于对所述恒速压汞曲线与所述高压压汞曲线进行拼接，以生成所述全尺度喉道分布曲线。

第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行程序时实现基于致密储层的微裂缝有效开度确定方法的步骤。

第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现基于致密储层的微裂缝有效开度确定方法的步骤。

由上述描述可知，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法及装置，针对致密油储层微裂缝发育，孔隙喉道细小，且非均质性强的特点，将决定储层渗流能力的喉道半径和微裂缝开度作为有效性的关键评价参数，将孔隙与喉道区别开来，以平均喉道半径和微裂缝开度作为储层渗流能力的主要表征参数并且对于双重介质致密油储层，首次提出了确定微裂缝开度的方法。具体地，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法根据全尺度喉道分布，重新获取平均喉道半径，并计算不同尺寸喉道对渗透率的贡献，进而确定有效喉道半径下限，依据平均喉道半径与渗透率的相关公式，得到有效渗透率下限，根据微裂缝开度、微裂缝间距和渗透率的关系式计算可得到有效微裂缝开度的下限。从而克服了现有技术中，只针对单一孔隙介质，对于含微裂缝的双重介质致密油储层不适应；只能确定单一孔隙介质致密油储层的孔喉半径下限，没有将孔隙与喉道分开；没有提出针对双重介质致密油储层有效微裂缝开度下限的确定方法的问题。

附图说明

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

图1为本发明的实施例中的基于致密储层的微裂缝有效开度确定方法的流程示意图；

图2为本发明的实施例中步骤100的流程示意图；

图3为本发明的实施例中步骤101的流程示意图；

图4为本发明的实施例中步骤200的流程示意图；

图5为本发明的实施例中步骤300的流程示意图；

图6为本发明的具体应用实例中基于致密储层的微裂缝有效开度确定方法的流程示意图；

图7为本发明的具体应用实例中恒速压汞实验喉道分布曲线示意图；

图8为本发明的具体应用实例中高压压汞实验喉道分布曲线示意图；

图9为本发明的具体应用实例中全尺度喉道分布曲线示意图；

图10为本发明的具体应用实例中平均喉道半径与渗透率的相关关系曲线示意图；

图11为本发明的具体应用实例中喉道对渗透率累计贡献的分布曲线示意图；

图12为本发明的实施例中基于致密储层的微裂缝有效开度确定装置的结构框图；

图13为本发明的实施例中岩心分析单元结构框图；

图14为本发明的实施例中有效渗透率确定单元结构框图；

图15为本发明的实施例中有效开度确定单元结构框图；

图16为本发明的实施例中曲线生成模块结构框图；

图17为本发明的实施例中的电子设备的结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

基于现有技术中的问题，本发明的实施例提供一种基于致密储层的微裂缝有效开度确定方法的具体实施方式，参见图1，该方法具体包括如下内容：

步骤100：对所述致密储层的岩心进行岩心分析，以得到所述岩心不同喉道对渗透率的贡献率。

可以理解的是，步骤100中的岩心分析可以包括：

显微技术：偏光显微镜、阴极发光显微镜、荧光显微镜、激光显微镜、电子显微镜、显微镜图象分析；

分光技术：紫外光谱、红外光谱、X射线荧光光谱，穆式鲍尔光谱；

其他技术：X射线衍射、电子探针，差热及热重分析、中子活化、核磁共振、核伽马共振、薄片染色微化分析等。

可以理解的是，步骤100将孔隙与喉道分开，考虑到致密油储层的孔隙与喉道存在明显差异，从而更符合致密储层的微裂缝有效开度确定。

步骤200：根据所述贡献率确定所述岩心的有效渗透率。

进一步地，步骤200在实施时，具体为：根据贡献率确定所述岩心的有效渗透率下线。其中有效渗透率是大于能够允许流体渗流的最小喉道称为有效喉道下限，大于有效喉道下限的所有喉道集合对应的渗透率。

步骤300：根据所述有效渗透率确定所述岩心的微裂缝有效开度。

可以理解的是，微裂缝是按裂缝大小(主要依据宽度)而作的一种分类。裂缝宽度在100微米或150微米以下的称微裂缝。微裂缝的开度即为微裂缝的宽度。

由上述描述可知，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法，针对致密油储层微裂缝发育，孔隙喉道细小，且非均质性强的特点，将决定储层渗流能力的喉道半径和微裂缝开度作为有效性的关键评价参数，将孔隙与喉道区别开来，以平均喉道半径和微裂缝开度作为储层渗流能力的主要表征参数并且对于双重介质致密油储层，首次提出了确定微裂缝开度的方法。具体地，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法根据全尺度喉道分布，重新获取平均喉道半径，并计算不同尺寸喉道对渗透率的贡献，进而确定有效喉道半径下限，依据平均喉道半径与渗透率的相关公式，得到有效渗透率下限，根据微裂缝开度、微裂缝间距和渗透率的关系式计算可得到有效微裂缝开度的下限。从而克服了现有技术中，只针对单一孔隙介质，对于含微裂缝的双重介质致密油储层不适应；只能确定单一孔隙介质致密油储层的孔喉半径下限，没有将孔隙与喉道分开；没有提出针对双重介质致密油储层有效微裂缝开度下限的确定方法的问题。

一实施例中，参见图2，步骤100具体包括：

步骤101：对所岩心进行压汞实验，以生成所述岩心的全尺度喉道分布曲线。

步骤102：根据所述全尺度喉道分布曲线计算所述贡献率。

可以理解的是，压汞实验的原理为：汞对一般固体不润湿，欲使汞进入孔需施加外压，外压越大，汞能进入的孔半径越小。测量不同外压下进入孔中汞的量即可知相应孔大小的孔体积。

一实施例中，参见图3，步骤101进一步又可分为：

步骤1011：对所述岩心进行恒速压汞实验，以生成恒速压汞曲线。

步骤1012：对所述岩心进行高压压汞实验，以生成高压压汞曲线。

步骤1013：对所述恒速压汞曲线与所述高压压汞曲线进行拼接，以生成所述全尺度喉道分布曲线。

步骤1011至步骤1013，将步骤1011中的恒速压汞实验数据以及步骤1012中的高压压汞实验实验数据整理成两条以喉道半径为横坐标，喉道体积为纵坐标的分布曲线，并选取一数据点为拼接点对上述两条分布曲线进行拼接。

一实施例中，参见图4，步骤200具体包括：

步骤201：对所述全尺度喉道分布曲线进行归一化处理，并计算平均喉道半径。

步骤202：确定所述平均喉道半径与所述渗透率之间的关系式。

步骤203：根据所述关系式确定所述岩心的有效喉道半径。

步骤204：根据所述有效喉道半径确定所述有效渗透率。

一实施例中，参见图5，步骤300具体包括：

步骤301：根据所述有效渗透率确定间距最大的微裂缝所对应的有效渗透率。

具体地，对岩心进行铸体薄片和扫描电镜测量，得到微裂缝间距最大值，并根据有效渗透率查找该最大值所对应的有效渗透率(下限)。

步骤302：根据所述间距最大的微裂缝所对应的有效渗透率确定所述微裂缝有效开度。

由上述描述可知，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法，针对致密油储层微裂缝发育，孔隙喉道细小，且非均质性强的特点，将决定储层渗流能力的喉道半径和微裂缝开度作为有效性的关键评价参数，将孔隙与喉道区别开来，以平均喉道半径和微裂缝开度作为储层渗流能力的主要表征参数并且对于双重介质致密油储层，首次提出了确定微裂缝开度的方法。具体地，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法根据全尺度喉道分布，重新获取平均喉道半径，并计算不同尺寸喉道对渗透率的贡献，进而确定有效喉道半径下限，依据平均喉道半径与渗透率的相关公式，得到有效渗透率下限，根据微裂缝开度、微裂缝间距和渗透率的关系式计算可得到有效微裂缝开度的下限。从而克服了现有技术中，只针对单一孔隙介质，对于含微裂缝的双重介质致密油储层不适应；只能确定单一孔隙介质致密油储层的孔喉半径下限，没有将孔隙与喉道分开；没有提出针对双重介质致密油储层有效微裂缝开度下限的确定方法的问题。

为进一步地说明本方案，以下以某油田A区块为例，提供基于致密储层的微裂缝有效开度确定方法的具体应用实例，该具体应用实例具体包括如下内容，参见图6。

S0：对所述致密储层的岩心进行岩心分析。

具体地，在全直径岩心上钻取直径为2.5cm，长度为6cm的岩心7块，洗油烘干后测孔隙度分别为10.7％、10.8％、11.1％、9.3％、11.2％、5.8％、2.8％，渗透率分别为0.73mD、0.47mD、0.21mD、0.15mD、0.11mD、0.019mD、0.001mD。接着，将岩心分为四部分，长度分别为2.0m、2.0m、1.0cm和1.0cm；对两个2.0cm长度的岩心分别开展恒速压汞和高压压汞实验，对两个1.0cm长度的岩心分别开展扫描电镜和铸体薄片实验。

S1：生成全尺度喉道分布曲线。

根据恒速压汞实验和高压压汞实验结果，将数据处理成横坐标为喉道半径(r)，纵坐标单位为喉道体积(v)的分布曲线，参见图7以及图8。在0.12μm处将恒速压汞和高压压汞曲线衔接，得到全尺度喉道分布曲线，得到图9(全尺度喉道分布曲线)。

S2：确定平均喉道半径与渗透率之间的关系式。

对全尺度喉道分布进行归一化处理，利用式(1)计算平均喉道半径；

式中：

r

n表示数据点个数；

α

重复步骤四至步骤六得到所有样品的平均喉道半径；做平均喉道半径与渗透率的相关关系曲线(散点图)，得到图10，并得到平均喉道半径与渗透率的相关关系函数为

S3：计算不同喉道对渗透率的贡献率。

具体地，利用式(2)计算得到不同半径喉道对渗透率的相对贡献；

式中：ΔK

S4：计算有效渗透率下限。

作横坐标为喉道半径，纵坐标为渗透率累计贡献的分布曲线，如图11所示；确定7块样品渗透率贡献达到99.99％时的喉道半径，即为有效喉道半径下限(r

S5：确定微裂缝有效开度下限值。

具体地，对所有样品的铸体薄片和扫描电镜进行测量，得到微裂缝间距最大值(D

由上述描述可知，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法，针对致密油储层微裂缝发育，孔隙喉道细小，且非均质性强的特点，将决定储层渗流能力的喉道半径和微裂缝开度作为有效性的关键评价参数，将孔隙与喉道区别开来，以平均喉道半径和微裂缝开度作为储层渗流能力的主要表征参数并且对于双重介质致密油储层，首次提出了确定微裂缝开度的方法。具体地，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法根据全尺度喉道分布，重新获取平均喉道半径，并计算不同尺寸喉道对渗透率的贡献，进而确定有效喉道半径下限，依据平均喉道半径与渗透率的相关公式，得到有效渗透率下限，根据微裂缝开度、微裂缝间距和渗透率的关系式计算可得到有效微裂缝开度的下限。从而克服了现有技术中，只针对单一孔隙介质，对于含微裂缝的双重介质致密油储层不适应；只能确定单一孔隙介质致密油储层的孔喉半径下限，没有将孔隙与喉道分开；没有提出针对双重介质致密油储层有效微裂缝开度下限的确定方法的问题。

基于相似的发明构思，本申请实施例还提供了基于致密储层的微裂缝有效开度确定装置，优选地用于实现上述实施例所描述的方法。由于基于致密储层的微裂缝有效开度确定装置解决问题的原理与基于致密储层的微裂缝有效开度确定方法相似，因此基于致密储层的微裂缝有效开度确定装置的实施可以参见基于致密储层的微裂缝有效开度确定方法实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的系统较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

本发明的实施例提供一种能够实现基于致密储层的微裂缝有效开度确定方法的基于致密储层的微裂缝有效开度确定装置的具体实施方式，参见图12，基于致密储层的微裂缝有效开度确定装置具体包括如下内容：

岩心分析单元10，用于对所述致密储层的岩心进行岩心分析，以得到所述岩心不同喉道对渗透率的贡献率；

有效渗透率确定单元20，用于根据所述贡献率确定所述岩心的有效渗透率；

有效开度确定单元30，用于根据所述有效渗透率确定所述岩心的微裂缝有效开度。

一实施例中，参见图13，所述岩心分析单元10包括：

分布曲线生成模块101，用于对所岩心进行压汞实验，以生成所述岩心的全尺度喉道分布曲线；

贡献率计算模块102，用于根据所述全尺度喉道分布曲线计算所述贡献率。

一实施例中，参见图14，所述有效渗透率确定单元20包括：

平均喉道半径计算模块201，用于对所述全尺度喉道分布曲线进行归一化处理，并计算平均喉道半径；

关系式确定模块202，用于确定所述平均喉道半径与所述渗透率之间的关系式；

有效喉道半径确定模块203，用于根据所述关系式确定所述岩心的有效喉道半径；

有效渗透率确定模块204，用于根据所述有效喉道半径确定所述有效渗透率。

一实施例中，参见图15，所述有效开度确定单元30包括：

有效渗透率确定模块301，用于根据所述有效渗透率确定间距最大的微裂缝所对应的有效渗透率；

有效开度确定模块302，用于根据所述间距最大的微裂缝所对应的有效渗透率确定所述微裂缝有效开度。

一实施例中，参见图16，所述分布曲线生成模块101包括：

恒速压汞曲线生成模块1011，用于对所述岩心进行恒速压汞实验，以生成恒速压汞曲线；

高压压汞曲线生成模块1012，用于对所述岩心进行高压压汞实验，以生成高压压汞曲线；

拼接模块1013，用于对所述恒速压汞曲线与所述高压压汞曲线进行拼接，以生成所述全尺度喉道分布曲线。

由上述描述可知，本发明实施例提供的基于致密储层的微裂缝有效开度确定装置，针对致密油储层微裂缝发育，孔隙喉道细小，且非均质性强的特点，将决定储层渗流能力的喉道半径和微裂缝开度作为有效性的关键评价参数，将孔隙与喉道区别开来，以平均喉道半径和微裂缝开度作为储层渗流能力的主要表征参数并且对于双重介质致密油储层，首次提出了确定微裂缝开度的方法。具体地，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法根据全尺度喉道分布，重新获取平均喉道半径，并计算不同尺寸喉道对渗透率的贡献，进而确定有效喉道半径下限，依据平均喉道半径与渗透率的相关公式，得到有效渗透率下限，根据微裂缝开度、微裂缝间距和渗透率的关系式计算可得到有效微裂缝开度的下限。从而克服了现有技术中，只针对单一孔隙介质，对于含微裂缝的双重介质致密油储层不适应；只能确定单一孔隙介质致密油储层的孔喉半径下限，没有将孔隙与喉道分开；没有提出针对双重介质致密油储层有效微裂缝开度下限的确定方法的问题。

本申请的实施例还提供能够实现上述实施例中的基于致密储层的微裂缝有效开度确定方法中全部步骤的一种电子设备的具体实施方式，参见图17，电子设备具体包括如下内容：

处理器(processor)1201、存储器(memory)1202、通信接口(CommunicationsInterface)1203和总线1204；

其中，处理器1201、存储器1202、通信接口1203通过总线1204完成相互间的通信；通信接口1203用于实现服务器端设备、渗透率计算设备以及用户端设备等相关设备之间的信息传输。

处理器1201用于调用存储器1202中的计算机程序，处理器执行计算机程序时实现上述实施例中的基于致密储层的微裂缝有效开度确定方法中的全部步骤，例如，处理器执行计算机程序时实现下述步骤：

步骤100：对所述致密储层的岩心进行岩心分析，以得到所述岩心不同喉道对渗透率的贡献率。

步骤200：根据所述贡献率确定所述岩心的有效渗透率。

步骤300：根据所述有效渗透率确定所述岩心的微裂缝有效开度。

由上述描述可知，本申请实施例中的电子设备，针对致密油储层微裂缝发育，孔隙喉道细小，且非均质性强的特点，将决定储层渗流能力的喉道半径和微裂缝开度作为有效性的关键评价参数，将孔隙与喉道区别开来，以平均喉道半径和微裂缝开度作为储层渗流能力的主要表征参数并且对于双重介质致密油储层，首次提出了确定微裂缝开度的方法。具体地，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法根据全尺度喉道分布，重新获取平均喉道半径，并计算不同尺寸喉道对渗透率的贡献，进而确定有效喉道半径下限，依据平均喉道半径与渗透率的相关公式，得到有效渗透率下限，根据微裂缝开度、微裂缝间距和渗透率的关系式计算可得到有效微裂缝开度的下限。从而克服了现有技术中，只针对单一孔隙介质，对于含微裂缝的双重介质致密油储层不适应；只能确定单一孔隙介质致密油储层的孔喉半径下限，没有将孔隙与喉道分开；没有提出针对双重介质致密油储层有效微裂缝开度下限的确定方法的问题。

本申请的实施例还提供能够实现上述实施例中的基于致密储层的微裂缝有效开度确定方法中全部步骤的一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的基于致密储层的微裂缝有效开度确定方法的全部步骤，例如，处理器执行计算机程序时实现下述步骤：

步骤100：对所述致密储层的岩心进行岩心分析，以得到所述岩心不同喉道对渗透率的贡献率。

步骤200：根据所述贡献率确定所述岩心的有效渗透率。

步骤300：根据所述有效渗透率确定所述岩心的微裂缝有效开度。

由上述描述可知，本申请实施例中的计算机可读存储介质，针对致密油储层微裂缝发育，孔隙喉道细小，且非均质性强的特点，将决定储层渗流能力的喉道半径和微裂缝开度作为有效性的关键评价参数，将孔隙与喉道区别开来，以平均喉道半径和微裂缝开度作为储层渗流能力的主要表征参数并且对于双重介质致密油储层，首次提出了确定微裂缝开度的方法。具体地，本发明实施例提供的基于致密储层的微裂缝有效开度确定方法根据全尺度喉道分布，重新获取平均喉道半径，并计算不同尺寸喉道对渗透率的贡献，进而确定有效喉道半径下限，依据平均喉道半径与渗透率的相关公式，得到有效渗透率下限，根据微裂缝开度、微裂缝间距和渗透率的关系式计算可得到有效微裂缝开度的下限。从而克服了现有技术中，只针对单一孔隙介质，对于含微裂缝的双重介质致密油储层不适应；只能确定单一孔隙介质致密油储层的孔喉半径下限，没有将孔隙与喉道分开；没有提出针对双重介质致密油储层有效微裂缝开度下限的确定方法的问题。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。

虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。

虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。

为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。

内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。

计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。

本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。

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