泥岩盖层封闭性的预测方法及装置

摘要

本发明公开了一种泥岩盖层封闭性的预测方法及装置，其中，该方法包括：获得对泥岩盖层样品执行压力测试得到的泥岩盖层样品的突破压力；根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定粉碎后的泥岩盖层样品的微观特征参数；根据突破压力和微观特征参数，确定突破压力与微观特征参数之间的关联关系；根据待分析区泥岩盖层的微观特征参数，基于突破压力与微观特征参数之间的关联关系，确定待分析区泥岩盖层的突破压力；根据待分析区泥岩盖层的突破压力，预测待分析区泥岩盖层封闭性，本发明可以提高泥岩盖层封闭性的预测效率和精度。

权利要求

1.一种泥岩盖层封闭性的预测方法，其特征在于，包括：获得对泥岩盖层样品执行压力测试得到的泥岩盖层样品的突破压力；

根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定所述粉碎后的泥岩盖层样品的微观特征参数；

根据所述突破压力和所述微观特征参数，确定突破压力与微观特征参数之间的关联关系；

根据待分析区泥岩盖层的微观特征参数，基于所述突破压力与微观特征参数之间的关联关系，确定待分析区泥岩盖层的突破压力；

根据所述待分析区泥岩盖层的突破压力，预测待分析区泥岩盖层封闭性；

其中，根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定所述粉碎后的泥岩盖层样品的微观特征参数，包括：根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附等温曲线和氮气脱附等温曲线，确定所述粉碎后的泥岩盖层样品的微观特征参数；

其中，根据所述突破压力和所述微观特征参数，确定突破压力与微观特征参数之间的关联关系，包括：确定所述突破压力和所述微观特征参数之间的相关性；

根据所述相关性，对所述突破压力和所述微观特征参数进行线性拟合，确定突破压力与微观特征参数之间的关联关系；

其中，所述突破压力与微观特征参数之间的关联关系包括：P＝25.8508‑1.42166×Φ‑0.36571×Ravg‑0.000803696×S+18.5794×Φmid+24.0194×Φbig‑1449.79×Φmic；

其中，P为突破压力；Φ为有效孔隙度；Ravg为平均孔隙半径；S为比表面积；Φmid为中孔孔隙体积；Φbig为大孔孔隙体积；Φmic为微孔孔隙体积。

2.如权利要求1所述的方法，其特征在于，根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附等温曲线和氮气脱附等温曲线，确定所述粉碎后的泥岩盖层样品的微观特征参数，包括：根据氮气吸附等温曲线和氮气脱附等温曲线，基于多点比表面积测试模型和线性回归，确定所述粉碎后的泥岩盖层样品的比表面积；

根据氮气吸附等温曲线和氮气脱附等温曲线，基于密度泛函模型，确定所述粉碎后的泥岩盖层样品的孔径分布和孔隙度。

3.如权利要求1所述的方法，其特征在于，所述突破压力和所述微观特征参数之间的相关性包括：所述突破压力与比表面积、微孔孔隙体积和中孔孔隙体积之间的相关性均为正相关；

所述突破压力与大孔孔隙体积、平均孔径和有效孔隙度之间的相关性均为负相关。

4.一种泥岩盖层封闭性的预测装置，其特征在于，包括：样品突破压力确定模块，用于获得对泥岩盖层样品执行压力测试得到的泥岩盖层样品的突破压力；

碎样品微观参数确定模块，用于根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定所述粉碎后的泥岩盖层样品的微观特征参数；

关联关系确定模块，用于根据所述突破压力和所述微观特征参数，确定突破压力与微观特征参数之间的关联关系；

区域突破压力确定模块，用于根据待分析区泥岩盖层的微观特征参数，基于所述突破压力与微观特征参数之间的关联关系，确定待分析区泥岩盖层的突破压力；

泥岩盖层封闭性预测模块，用于根据所述待分析区泥岩盖层的突破压力，预测待分析区泥岩盖层封闭性；

其中，所述碎样品微观参数确定模块具体用于：

根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附等温曲线和氮气脱附等温曲线，确定所述粉碎后的泥岩盖层样品的微观特征参数；

其中，所述关联关系确定模块具体用于：

确定所述突破压力和所述微观特征参数之间的相关性；

根据所述相关性，对所述突破压力和所述微观特征参数进行线性拟合，确定突破压力与微观特征参数之间的关联关系；

其中，所述突破压力与微观特征参数之间的关联关系包括：P＝25.8508‑1.42166×Φ‑0.36571×Ravg‑0.000803696×S+18.5794×Φmid+24.0194×Φbig‑1449.79×Φmic；

其中，P为突破压力；Φ为有效孔隙度；Ravg为平均孔隙半径；S为比表面积；Φmid为中孔孔隙体积；Φbig为大孔孔隙体积；Φmic为微孔孔隙体积。

5.如权利要求4所述的装置，其特征在于，所述碎样品微观参数确定模块进一步用于：根据氮气吸附等温曲线和氮气脱附等温曲线，基于多点比表面积测试模型和线性回归，确定所述粉碎后的泥岩盖层样品的比表面积；

根据氮气吸附等温曲线和氮气脱附等温曲线，基于密度泛函模型，确定所述粉碎后的泥岩盖层样品的孔径分布和孔隙度。

6.如权利要求4所述的装置，其特征在于，所述突破压力和所述微观特征参数之间的相关性包括：所述突破压力与比表面积、微孔孔隙体积和中孔孔隙体积之间的相关性均为正相关；

所述突破压力与大孔孔隙体积、平均孔径和有效孔隙度之间的相关性均为负相关。

7.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至3任一所述方法。

8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至3任一所述方法的计算机程序。

说明书

泥岩盖层封闭性的预测方法及装置

技术领域

[0001] 本发明涉及石油勘探开发技术领域，特别涉及一种泥岩盖层封闭性的预测方法及装置。

背景技术

[0002] 泥岩盖层只有具备较强的微观封闭能力时，才能形成大面积的有效盖层，对成藏系统内的油气进行封盖。大量研究表明，盖层毛细管封闭能力的强弱主要取决于其内突破压力的大小，突破压力越大，毛细管封闭能力越强；反之，则越弱。游离相油气欲通过泥岩盖层向上运移，首先必须克服其下致密层段的毛细管阻力，当油气能量小于或等于其毛细管阻力时，油气不能穿过超压泥岩盖层进行运移，盖层才能对油气起封隔作用，因此突破压力是反应盖层封盖能力最直观的参数。

[0003] 现有技术多采用直接测定的方法测定突破压力值，首先钻取多个泥岩圆柱样品，然后在配置的地层水溶液中使圆柱样品饱和地层水，最后再进行气驱，获得突破压力值，由于泥岩块状岩心样品取心比较困难，尤其是深层泥岩块状样品更难获取，且泥岩圆柱体钻取时成功率低，在饱和地层水过程中也易松散破碎，导致样品破损，使得样品数量少，实验数据少，只能对局部的盖层封闭能力进行预测，难以实现对大面积盖层封闭能力的预测。

[0004] 此外，由于影响盖层突破压力大小的因素较多，现有技术通过对根据深度变化或测井曲线计算泥岩盖层突破压力值时，计算参数比较单一，计算结果的不确定性较大。

[0005] 针对上述问题，目前尚未提出有效的解决方案。

发明内容

[0006] 本发明实施例提供一种泥岩盖层封闭性的预测方法，用于提高泥岩盖层封闭性的预测效率和精度，该方法包括：

[0007] 获得对泥岩盖层样品执行压力测试得到的泥岩盖层样品的突破压力；

[0008] 根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定粉碎后的泥岩盖层样品的微观特征参数；

[0009] 根据突破压力和微观特征参数，确定突破压力与微观特征参数之间的关联关系；

[0010] 根据待分析区泥岩盖层的微观特征参数，基于突破压力与微观特征参数之间的关联关系，确定待分析区泥岩盖层的突破压力；

[0011] 根据待分析区泥岩盖层的突破压力，预测待分析区泥岩盖层封闭性。

[0012] 本发明实施例提供一种泥岩盖层封闭性的预测装置，用于提高泥岩盖层封闭性的预测效率和精度，该装置包括：

[0013] 样品突破压力确定模块，用于获得对泥岩盖层样品执行压力测试得到的泥岩盖层样品的突破压力；

[0014] 碎样品微观参数确定模块，用于根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定粉碎后的泥岩盖层样品的微观特征参数；

[0015] 关联关系确定模块，用于根据突破压力和微观特征参数，确定突破压力与微观特征参数之间的关联关系；

[0016] 区域突破压力确定模块，用于根据待分析区泥岩盖层的微观特征参数，基于突破压力与微观特征参数之间的关联关系，确定待分析区泥岩盖层的突破压力；

[0017] 泥岩盖层封闭性预测模块，用于根据待分析区泥岩盖层的突破压力，预测待分析区泥岩盖层封闭性。

[0018] 本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述泥岩盖层封闭性的预测方法。

[0019] 本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质存储有执行上述泥岩盖层封闭性的预测方法的计算机程序。

[0020] 本发明实施例通过：获得对泥岩盖层样品执行压力测试得到的泥岩盖层样品的突破压力；根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定粉碎后的泥岩盖层样品的微观特征参数；根据突破压力和微观特征参数，确定突破压力与微观特征参数之间的关联关系；根据待分析区泥岩盖层的微观特征参数，基于突破压力与微观特征参数之间的关联关系，确定待分析区泥岩盖层的突破压力；根据待分析区泥岩盖层的突破压力，预测待分析区泥岩盖层封闭性。本发明可以基于碎泥氮气吸附测试和压力测试，建立突破压力与微观特征参数之间的关联关系，考虑了与突破压力相关的多种微观特征参数，提高了预测精度，基于突破压力与微观特征参数之间的关联关系可以实现区域泥岩盖层封闭性的预测，解决了泥岩样品取样困难的问题，提高了预测效率，适用性强。

附图说明

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

[0022] 图1为本发明实施例中泥岩盖层封闭性的预测方法流程的示意图；

[0023] 图2为本发明实施例中不同孔径下孔隙体积变化的示意图；

[0024] 图3为本发明实施例中不同孔径下比表面积变化的示意图；

[0025] 图4为本发明实施例中不同孔径下累积孔隙体积的示意图；

[0026] 图5为本发明实施例中突破压力与微观特征参数相关性的示意图；

[0027] 图6为本发明实施例中泥岩盖层封闭性的预测装置结构的示意图。

具体实施方式

[0028] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

[0029] 为了解决现有技术采用直接测定的方法测定突破压力值时，泥岩样品取心困难，样品易破损，使得样品数量少，预测效率低的技术问题，本发明实施例提供一种泥岩盖层封闭性的预测方法，用于提高泥岩盖层封闭性的预测效率和精度，图1为本发明实施例中泥岩盖层封闭性的预测方法流程的示意图，如图1所示，该方法包括：

[0030] 步骤101：获得对泥岩盖层样品执行压力测试得到的泥岩盖层样品的突破压力；

[0031] 步骤102：根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定粉碎后的泥岩盖层样品的微观特征参数；

[0032] 步骤103：根据突破压力和微观特征参数，确定突破压力与微观特征参数之间的关联关系；

[0033] 步骤104：根据待分析区泥岩盖层的微观特征参数，基于突破压力与微观特征参数之间的关联关系，确定待分析区泥岩盖层的突破压力；

[0034] 步骤105：根据待分析区泥岩盖层的突破压力，预测待分析区泥岩盖层封闭性。

[0035] 如图1所示，本发明实施例通过：获得对泥岩盖层样品执行压力测试得到的泥岩盖层样品的突破压力；根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定粉碎后的泥岩盖层样品的微观特征参数；根据突破压力和微观特征参数，确定突破压力与微观特征参数之间的关联关系；根据待分析区泥岩盖层的微观特征参数，基于突破压力与微观特征参数之间的关联关系，确定待分析区泥岩盖层的突破压力；根据待分析区泥岩盖层的突破压力，预测待分析区泥岩盖层封闭性。本发明可以基于碎泥氮气吸附测试和压力测试，建立突破压力与微观特征参数之间的关联关系，考虑了与突破压力相关的多种微观特征参数，提高了预测精度，基于突破压力与微观特征参数之间的关联关系可以实现区域泥岩盖层封闭性的预测，解决了泥岩样品取样困难的问题，提高了预测效率，适用性强。

[0036] 具体实施时，步骤101中，可以按照如下方式对泥岩盖层样品执行压力测试，包括：

[0037] 第一步：钻取圆柱体：在块状泥岩盖层样品上钻取标准规格圆柱体泥岩盖层样品若干；

[0038] 第二步：岩心流体饱和：对泥岩盖层样品抽真空后，注入模拟地层水并加压，饱和至少24小时，直至泥岩盖层样品完全饱和地层水流体；

[0039] 第三步：突破压力值测试：首先对饱和地层水后完好的泥岩盖层样品测试岩心长度，直径等基础参数，然后将泥岩盖层样品置于岩心夹持器中，上游与高压气源相连，出口处接软管并置于水中以便观察气泡冒出，连接好仪器，并排出管路中的残余气体，为确保不漏气，先将围压加至8～10MPa，气体压力从0.01MPa开始按上一次压力的15％每隔10min逐级往上加压，当出口有连续气泡冒出，说明已达到泥岩盖层样品的突破压力，实验过程中记录泥岩盖层样品的突破时间和突破压力，这样就可以得到泥岩盖层样品的突破压力。

[0040] 在一个实施例中，步骤102可以包括：

[0041] 根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附等温曲线和氮气脱附等温曲线，确定粉碎后的泥岩盖层样品的微观特征参数。

[0042] 在一个实施例中，上述根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附等温曲线和氮气脱附等温曲线，确定粉碎后的泥岩盖层样品的微观特征参数，可以包括：

[0043] 根据氮气吸附等温曲线和氮气脱附等温曲线，基于多点比表面积测试模型和线性回归，确定粉碎后的泥岩盖层样品的比表面积；

[0044] 根据氮气吸附等温曲线和氮气脱附等温曲线，基于密度泛函模型，确定粉碎后的泥岩盖层样品的孔径分布和孔隙度。

[0045] 具体实施时，由于影响盖层突破压力大小的因素较多，现有技术通过对根据深度变化或测井曲线计算泥岩盖层突破压力值时，计算参数比较单一，计算结果的不确定性较大，为了解决上述计算结果不确定性较大的技术问题，步骤102中，可以将成功获得突破压力值的圆柱体泥岩盖层样品粉碎至40‑60目，然后按照如下方式进行氮气吸附测试：在150℃条件下进行真空脱气8h以去除泥碎岩盖层样品中残留的气体，然后通入高纯氮气，在‑195.8℃条件下进行等温物理吸附‑脱附测试，其中，孔径测量范围为0.35‑500nm，吸附测试分为两个阶段，第一阶段相对压力范围(p/p0)为0‑0.011，每个压力点平衡时间为30s；第二阶段吸附测试相对压力(p/p0)范围为0.011‑0.998，每个压力点平衡时间为10s；脱附测试相对压力(p/p0)范围为0.011‑0.998，每个压力点平衡时间为10s，这样可以得到碎泥岩盖层样品的氮气吸附等温曲线和氮气脱附等温曲线。

[0046] 接着，根据碎泥岩盖层样品的氮气吸附测试结果，确定碎泥岩盖层样品的微观特征参数，其中，微观特征参数可以包括：比表面积、孔径分布和孔隙度等参数。根据碎泥岩盖层样品的氮气吸附等温曲线和氮气脱附等温曲线，可以基于多点BET模型和线性回归，确定碎泥岩盖层样品的比表面积，还可以基于DFT模型确定碎泥岩盖层样品的孔径分布和孔隙度，其中，孔径分布可以包括孔径、平均孔隙半径，以及微孔、中孔和大孔孔隙体积等参数。图2为本发明实施例中不同孔径下孔隙体积变化的示意图，图3为本发明实施例中不同孔径下比表面积变化的示意图，图4为本发明实施例中不同孔径下累积孔隙体积的示意图，如图
2、图3和图4所示，可以将碎泥岩盖层样品的各种微观特征参数进行对比分析，得到碎泥岩盖层样品的孔径、孔隙体积以及比表面积之间的相关性，可以用于后续步骤的分析。

[0047] 在一个实施例中，步骤103可以包括：

[0048] 确定突破压力和微观特征参数之间的相关性；

[0049] 根据相关性，对突破压力和微观特征参数进行线性拟合，确定突破压力与微观特征参数之间的关联关系。

[0050] 在一个实施例中，上述突破压力和微观特征参数之间的相关性可以包括：

[0051] 突破压力与比表面积、微孔孔隙体积和中孔孔隙体积之间的相关性均为正相关；

[0052] 突破压力与大孔孔隙体积、平均孔径和有效孔隙度之间的相关性均为负相关。

[0053] 在一个实施例中，上述突破压力与微观特征参数之间的关联关系包括：

[0054] P＝25.8508‑1.42166×Φ‑0.36571×Ravg‑0.000803696×S+18.5794×Φmid+24.0194×Φbig‑1449.79×Φmic

[0055]                                                (1)

[0056] 其中，P为突破压力；Φ为有效孔隙度；Ravg为平均孔隙半径；S为比表面积；Φmid为中孔孔隙体积；Φbig为大孔孔隙体积；Φmic为微孔孔隙体积。

[0057] 具体实施时，步骤103中，可以对步骤101中获得的泥岩盖层样品的突破压力和步骤102中获得的碎泥岩盖层样品的微观特征参数进行相关性分析，如图5所示，通过图5可以发现，突破压力与比表面积、微孔孔隙体积和中孔孔隙体积之间的相关性均为正相关，突破压力与大孔孔隙体积、平均孔径和有效孔隙度之间的相关性均为负相关，说明泥岩盖层样品中微孔隙发育数量越多，比表面积越大，对油气的吸附能力越强，吸附阻力越大；大孔孔隙发育越多，越利于游离态油气通过，不利于油气的封闭。根据上述相关性分析发现，可以对大量实测的突破压力以及微观特征参数进行拟合，建立拟合关系式，确定突破压力的计算公式，如上述公式(1)所示。

[0058] 具体实施时，步骤104、步骤105中，可以根据待分析区的地震数据和测井数据，确定待分析区的泥岩盖层的比表面积、孔径分布和孔隙度等微观特征参数，然后基于突破压力与微观特征参数之间的关联关系，即上述公式(1)，确定待分析区的泥岩盖层的突破压力，根据待分析区的泥岩盖层的突破压力，预测待分析区的泥岩盖层封闭性，突破压力越大，泥岩盖层封闭性能力越强；反之，则越弱，这样就可以实现对大面积盖层封闭能力的预测。

[0059] 基于同一发明构思，本发明实施例中还提供了一种泥岩盖层封闭性的预测装置，如下面的实施例。由于泥岩盖层封闭性的预测装置解决问题的原理与泥岩盖层封闭性的预测方法相似，因此装置的实施可以参见方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

[0060] 本发明实施例提供一种泥岩盖层封闭性的预测装置，用于提高泥岩盖层封闭性的预测效率和精度，图6为本发明实施例中泥岩盖层封闭性的预测装置结构的示意图，如图6所示，该装置包括：

[0061] 样品突破压力确定模块01，用于获得对泥岩盖层样品执行压力测试得到的泥岩盖层样品的突破压力；

[0062] 碎样品微观参数确定模块02，用于根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定粉碎后的泥岩盖层样品的微观特征参数；

[0063] 关联关系确定模块03，用于根据突破压力和微观特征参数，确定突破压力与微观特征参数之间的关联关系；

[0064] 区域突破压力确定模块04，用于根据待分析区泥岩盖层的微观特征参数，基于突破压力与微观特征参数之间的关联关系，确定待分析区泥岩盖层的突破压力；

[0065] 泥岩盖层封闭性预测模块05，用于根据待分析区泥岩盖层的突破压力，预测待分析区泥岩盖层封闭性。

[0066] 在一个实施例中，碎样品微观参数确定模块02具体用于：

[0067] 根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附等温曲线和氮气脱附等温曲线，确定粉碎后的泥岩盖层样品的微观特征参数。

[0068] 在一个实施例中，碎样品微观参数确定模块02进一步用于：

[0069] 根据氮气吸附等温曲线和氮气脱附等温曲线，基于多点比表面积测试模型和线性回归，确定粉碎后的泥岩盖层样品的比表面积；

[0070] 根据氮气吸附等温曲线和氮气脱附等温曲线，基于密度泛函模型，确定粉碎后的泥岩盖层样品的孔径分布和孔隙度。

[0071] 在一个实施例中，关联关系确定模块03具体用于：

[0072] 确定突破压力和微观特征参数之间的相关性；

[0073] 根据相关性，对突破压力和微观特征参数进行线性拟合，确定突破压力与微观特征参数之间的关联关系。

[0074] 在一个实施例中，突破压力和微观特征参数之间的相关性包括：

[0075] 突破压力与比表面积、微孔孔隙体积和中孔孔隙体积之间的相关性均为正相关；

[0076] 突破压力与大孔孔隙体积、平均孔径和有效孔隙度之间的相关性均为负相关。

[0077] 在一个实施例中，突破压力与微观特征参数之间的关联关系包括：

[0078] P＝25.8508‑1.42166×Φ‑0.36571×Ravg‑0.000803696×S+18.5794×Φmid+24.0194×Φbig‑1449.79×Φmic；

[0079] 其中，P为突破压力；Φ为有效孔隙度；Ravg为平均孔隙半径；S为比表面积；Φmid为中孔孔隙体积；Φbig为大孔孔隙体积；Φmic为微孔孔隙体积。

[0080] 本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述泥岩盖层封闭性的预测方法。

[0081] 本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质存储有执行上述泥岩盖层封闭性的预测方法的计算机程序。

[0082] 综上，本发明实施例通过：获得对泥岩盖层样品执行压力测试得到的泥岩盖层样品的突破压力；根据对粉碎后的泥岩盖层样品执行氮气吸附测试得到的氮气吸附特征数据，确定粉碎后的泥岩盖层样品的微观特征参数；根据突破压力和微观特征参数，确定突破压力与微观特征参数之间的关联关系；根据待分析区泥岩盖层的微观特征参数，基于突破压力与微观特征参数之间的关联关系，确定待分析区泥岩盖层的突破压力；根据待分析区泥岩盖层的突破压力，预测待分析区泥岩盖层封闭性。本发明可以基于碎泥氮气吸附测试和压力测试，建立突破压力与微观特征参数之间的关联关系，考虑了与突破压力相关的多种微观特征参数，提高了预测精度，基于突破压力与微观特征参数之间的关联关系可以实现区域泥岩盖层封闭性的预测，解决了泥岩样品取样困难的问题，提高了预测效率，适用性强。

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

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

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

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

[0087] 以上仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。