储层裂缝渗吸质量预测方法及系统

摘要

本发明提供一种储层裂缝渗吸质量预测方法及系统。该方法包括：根据各个裂缝的毛细管力与上一渗吸时刻的重力确定各个裂缝在当前渗吸时刻的裂缝渗吸高度；当裂缝渗吸高度小于或等于岩心高度时，根据该裂缝的裂缝长度、该裂缝在当前渗吸时刻的裂缝渗吸高度、该裂缝长度的裂缝数量和湿相流体密度得到该裂缝在当前渗吸时刻的渗吸质量；否则根据该裂缝的裂缝长度、岩心高度、该裂缝长度的裂缝数量、湿相流体密度和平均迂曲度得到该裂缝在当前渗吸时刻的渗吸质量；根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案。本发明可以有效指导油气资源的开发过程。

权利要求

1.一种储层裂缝渗吸质量预测方法，其特征在于，包括：

获取各个裂缝的毛细管力和各个裂缝在上一渗吸时刻的重力；

根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度；

创建第一渗吸质量预测模型和第二渗吸质量预测模型；

当裂缝在当前渗吸时刻的裂缝渗吸高度小于或等于岩心高度时，输入该裂缝的裂缝长度、该裂缝在当前渗吸时刻的裂缝渗吸高度、该裂缝长度的裂缝数量和湿相流体密度至所述第一渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量；

当裂缝在当前渗吸时刻的裂缝渗吸高度大于岩心高度时，输入该裂缝的裂缝长度、所述岩心高度、该裂缝长度的裂缝数量、所述湿相流体密度和平均迂曲度至所述第二渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量；

根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案。

2.根据权利要求1所述的储层裂缝渗吸质量预测方法，其特征在于，还包括：根据各个裂缝的裂缝长度、表面张力和接触角确定各个裂缝的毛细管力；

根据各个裂缝在上一渗吸时刻的裂缝渗吸高度和所述湿相流体密度确定各个裂缝在上一渗吸时刻的重力。

3.根据权利要求2所述的储层裂缝渗吸质量预测方法，其特征在于，确定各个裂缝在当前渗吸时刻的裂缝渗吸高度包括：当裂缝的毛细管力小于裂缝在上一渗吸时刻的重力时，将该裂缝在上一渗吸时刻的裂缝渗吸高度作为该裂缝在当前渗吸时刻的裂缝渗吸高度；

当裂缝的毛细管力等于裂缝在上一渗吸时刻的重力时，根据获取的当前渗吸时刻、第一物性参数和各个裂缝的裂缝长度确定各个裂缝在当前渗吸时刻的裂缝渗吸高度；

其中，所述第一物性参数包括表面张力、接触角、湿相流体粘度和所述平均迂曲度。

4.根据权利要求1所述的储层裂缝渗吸质量预测方法，其特征在于，还包括：创建裂缝数量模型；

获取第二物性参数；其中，所述第二物性参数包括：岩心横截面、分形维数、最大裂缝长度、最小裂缝长度、裂缝孔隙度和裂缝长度步长；

输入各个裂缝长度和所述第二物性参数至所述裂缝数量模型中，得到各个裂缝长度的裂缝数量。

5.一种储层裂缝渗吸质量预测系统，其特征在于，包括：

第一获取单元，用于获取各个裂缝的毛细管力和各个裂缝在上一渗吸时刻的重力；

裂缝渗吸高度单元，用于根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度；

渗吸质量预测模型创建单元，用于创建第一渗吸质量预测模型和第二渗吸质量预测模型；

裂缝渗吸质量单元，用于当裂缝在当前渗吸时刻的裂缝渗吸高度小于或等于岩心高度时，输入该裂缝的裂缝长度、该裂缝在当前渗吸时刻的裂缝渗吸高度、该裂缝长度的裂缝数量和湿相流体密度至所述第一渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量；当裂缝在当前渗吸时刻的裂缝渗吸高度大于岩心高度时，输入该裂缝的裂缝长度、所述岩心高度、该裂缝长度的裂缝数量、所述湿相流体密度和平均迂曲度至所述第二渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量；

预测单元，用于根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案。

6.根据权利要求5所述的储层裂缝渗吸质量预测系统，其特征在于，还包括：毛细管力确定单元，用于根据各个裂缝的裂缝长度、表面张力和接触角确定各个裂缝的毛细管力；

重力确定单元，用于根据各个裂缝在上一渗吸时刻的裂缝渗吸高度和所述湿相流体密度确定各个裂缝在上一渗吸时刻的重力。

7.根据权利要求6所述的储层裂缝渗吸质量预测系统，其特征在于，所述裂缝渗吸高度单元具体用于：当裂缝的毛细管力小于裂缝在上一渗吸时刻的重力时，将该裂缝在上一渗吸时刻的裂缝渗吸高度作为该裂缝在当前渗吸时刻的裂缝渗吸高度；

当裂缝的毛细管力等于裂缝在上一渗吸时刻的重力时，根据获取的当前渗吸时刻、第一物性参数和各个裂缝的裂缝长度确定各个裂缝在当前渗吸时刻的裂缝渗吸高度；

其中，所述第一物性参数包括表面张力、接触角、湿相流体粘度和所述平均迂曲度。

8.根据权利要求5所述的储层裂缝渗吸质量预测系统，其特征在于，还包括：裂缝数量模型创建单元，用于创建裂缝数量模型；

第二获取单元，用于获取第二物性参数；其中，所述第二物性参数包括：岩心横截面、分形维数、最大裂缝长度、最小裂缝长度、裂缝孔隙度和裂缝长度步长；

裂缝数量单元，用于输入各个裂缝长度和所述第二物性参数至所述裂缝数量模型中，得到各个裂缝长度的裂缝数量。

9.一种计算机设备，包括存储器、处理器及存储在存储器上并在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一项所述的储层裂缝渗吸质量预测方法的步骤。

10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至4任一项所述的储层裂缝渗吸质量预测方法的步骤。

说明书

储层裂缝渗吸质量预测方法及系统

技术领域

[0001] 本发明涉及油藏勘探开发技术领域，具体地，涉及一种储层裂缝渗吸质量预测方法及系统。

背景技术

[0002] 致密气、页岩气等非常规油气资源的开发对保障国家能源安全具有重要意义。与常规油气储层不同，致密油气储层具有低孔低渗的特点，孔喉结构复杂，毛管力作用强，存在明显的自发渗吸现象。大规模水力压裂是一种普遍采用的提高致密储层产能的技术手段，大量的压裂液被注入储层来产生人造裂缝。然而在矿场应用中，超过50％的压裂液滞留在地下储层中，自发渗吸是造成压裂液返排量低的主要原因之一。研究表明致密储层中发育有大量天然微纳米裂缝，这些天然裂缝与压裂产生的人造裂缝交汇形成复杂的缝网结构，因此研究压裂液在天然/人造裂缝中的自发渗吸机理对致密气、页岩气的开发具有重要意义。

[0003] 目前对于裂缝自发渗吸机理的研究有称重法、体积法、CT扫描、核磁共振和中子照相技术等，但该类方法存在以下问题：(1)致密砂岩孔隙度一般小于10％，这意味着渗吸过程中岩心的质量变化很小，即对测量工具的精度要求很高，常规测量工具难以达到让人满意的精度，因此称重法和体积法虽然可以研究岩心渗吸质量随时间的动态变化，但往往误差较大。(2)CT扫描和核磁共振等方法虽然具有较高的精度，但其成本高昂且无法研究渗吸质量随时间的动态变化；中子照相技术虽然可以研究渗吸过程中的动态变化，但其对实验设备具有非常高的要求，实验成本也极其高昂，很难普及应用，不利于油气资源的开发。

发明内容

[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] 图7是本发明实施例中储层裂缝渗吸质量预测系统的结构框图；

[0029] 图8是本发明实施例中计算机设备的结构框图。

具体实施方式

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

[0031] 本领域技术人员知道，本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。

[0032] 鉴于现有技术精度低，成本高，不利于油气资源的开发，本发明实施例提供一种储层裂缝渗吸质量预测方法及系统，可以预测精确的裂缝渗吸质量，降低预测成本，进一步有效指导了油气资源的开发过程。以下结合附图对本发明进行详细说明。

[0033] 图1是本发明实施例中储层裂缝渗吸质量预测方法的流程图。如图1所示，储层裂缝渗吸质量预测方法包括：

[0034] S101：获取各个裂缝的毛细管力和各个裂缝在上一渗吸时刻的重力。

[0035] 在执行S101之前还包括：

[0036] 根据各个裂缝的裂缝长度、表面张力和接触角确定各个裂缝的毛细管力。

[0037] 可以通过以下公式确定各个裂缝的毛细管力：

[0038]

[0039] 其中，p为裂缝长度为l的裂缝的毛细管力，单位为N；l为裂缝长度，n为裂缝的开度与长度之比，σ为表面张力，θ为接触角。

[0040] 根据各个裂缝在上一渗吸时刻的裂缝渗吸高度和湿相流体密度确定各个裂缝在上一渗吸时刻的重力。

[0041] 可以通过以下公式确定各个裂缝在上一渗吸时刻的重力：

[0042]

[0043] 其中，Gt-Δt为裂缝长度为l的裂缝在上一渗吸时刻(t-Δt渗吸时刻)的重力，单位为N；ρ为湿相流体密度，g为重力加速度， 为裂缝长度为l的裂缝在上一渗吸时刻(t-Δt渗吸时刻)的裂缝渗吸高度，单位为m；t为渗吸时间，单位可以为s或min；Δt为时间步长，单位可以为s或min。裂缝长度为l的裂缝在初始时刻的裂缝渗吸高度 等于0。

[0044] S102：根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度。

[0045] 一实施例中，S102包括：

[0046] 当裂缝的毛细管力小于裂缝在上一渗吸时刻的重力时，将该裂缝在上一渗吸时刻的裂缝渗吸高度作为该裂缝在当前渗吸时刻的裂缝渗吸高度：

[0047] 当p＜Gt-Δt时， 为裂缝长度为l的裂缝在当前渗吸时刻(t时刻)的裂缝渗吸高度，单位为m。

[0048] 第一物性参数包括表面张力、接触角、湿相流体粘度和平均迂曲度。当裂缝的毛细管力等于裂缝在上一渗吸时刻的重力时，根据获取的当前渗吸时刻、第一物性参数和各个裂缝的裂缝长度确定各个裂缝在当前渗吸时刻的裂缝渗吸高度：

[0049] 当p＝Gt-Δt时，

[0050] 其中， 为裂缝长度为l的裂缝在当前渗吸时刻(t时刻)的裂缝渗吸高度，[0051]μ为湿相流体粘度，τ为平均迂曲度。

[0052] S103：创建第一渗吸质量预测模型和第二渗吸质量预测模型。

[0053] S104：当裂缝在当前渗吸时刻的裂缝渗吸高度小于或等于岩心高度时，输入该裂缝的裂缝长度、该裂缝在当前渗吸时刻的裂缝渗吸高度、该裂缝长度的裂缝数量和湿相流体密度至第一渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量。

[0054] 一实施例中，第一渗吸质量预测模型如下：

[0055]

[0056] 其中，MG(l)为裂缝长度为l的裂缝( 时)的当前渗吸时刻的渗吸质量，单位为g或kg，ρ为湿相流体密度，n(l)为裂缝长度l的裂缝数量。

[0057] S105：当裂缝在当前渗吸时刻的裂缝渗吸高度大于岩心高度时，输入该裂缝的裂缝长度、岩心高度、该裂缝长度的裂缝数量、湿相流体密度和平均迂曲度至第二渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量。

[0058] 一实施例中，第二渗吸质量预测模型如下：

[0059]

[0060] 其中， 为裂缝长度为l的裂缝( 时)的当前渗吸时刻的渗吸质量，单位为g或kg；H为岩心高度，单位为cm。

[0061] S106：根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案。

[0062] 图1所示的储层裂缝渗吸质量预测方法适用于致密砂岩和页岩，其中的裂缝渗吸高度为考虑重力时的裂缝渗吸高度，渗吸质量为考虑重力时的渗吸质量。由图1所示的流程可知，本发明实施例的储层裂缝渗吸质量预测方法先根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度，再根据裂缝渗吸高度和岩心高度的比较结果确定各个裂缝在当前渗吸时刻的渗吸质量，最后根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案，可以预测精确的裂缝渗吸质量，降低预测成本，进一步有效指导了油气资源的开发过程。

[0063] 在执行S101之前还包括：创建裂缝数量模型；获取第二物性参数；输入各个裂缝长度和第二物性参数至裂缝数量模型中，得到各个裂缝长度的裂缝数量。

[0064] 一实施例中，第二物性参数包括：岩心横截面、分形维数、最大裂缝长度、最小裂缝长度、裂缝孔隙度和裂缝长度步长。裂缝数量模型如下：

[0065]

[0066] 其中，n(l)为裂缝长度l的裂缝数量，l为裂缝长度，Af为岩心横截面，n为裂缝的开度与长度之比，Df为分形维数，lmax为最大裂缝长度，lmin为最小裂缝长度，φf为裂缝孔隙度，Δl为裂缝长度步长，单位为μm。

[0067] 表1是物性参数表。如表1所示，一实施例中的各个物性参数如下：

[0068] 表1

[0069]

[0070]

[0071] 其中，岩心横截面为

[0072] 一实施例中，还可以通过岩心在当前渗吸时刻的渗吸质量确定岩心在当前渗吸时刻的渗吸饱和度，渗吸饱和度可以有效指导油气资源的开发过程：

[0073]

[0074] 其中，RG为岩心在当前渗吸时刻的渗吸饱和度，单位为％； 为岩心在当前渗吸时刻的渗吸质量，单位为g或kg；d为岩心直径。

[0075] 图2是本发明实施例中湿相流体在气饱和的裂缝型岩心内的自发渗吸示意图。图3是本发明实施例中一簇迂曲平板裂缝组成的裂缝型岩心的理论模型示意图。图4是本发明实施例中湿相流体在单根裂缝内的自发渗吸示意图，也是图3方框中的放大图。

[0076] 如图2所示，裂缝2内的湿相流体3的流动非常类似于平板间的流动，因此将岩心1内的大量天然裂缝假设为一簇具有分形分布且互不交叉的迂曲平板裂缝，如图3所示。

[0077] 如图4所示，与湿相流体3相比，天然气4的粘度和密度可以忽略，根据立方定律，单根裂缝内的自发渗吸流量可以表示为：

[0078]

[0079] 其中q是单根裂缝内的自发渗吸流量，单位为mm3/s；a为裂缝开度，单位为μm；l为裂缝长度，单位为μm；Lf为渗吸前缘至岩心底面的迂曲长度，单位为m；Ls为渗吸前缘至岩心底面的直线长度(裂缝渗吸高度)，单位为m；μ为湿相流体粘度，单位为mPa·s；σ为天然气与3
湿相流体的表面张力，单位为mN/m；；θ为接触角，单位为°；ρ为湿相流体密度，单位为kg/m ；
g为重力加速度，单位为m/s2。式(1)可以写为如下微分形式：

[0080]

[0081] 渗吸前缘至岩心底面的迂曲长度和裂缝渗吸高度之间的关系可以表示为：

[0082] Lf＝τLs。  (3)

[0083] 其中τ为裂缝平均迂曲度，无量纲。将式(3)代入式(2)可得：

[0084]

[0085] 式(4)可以写为如下形式：

[0086]

[0087] 式(5)的解为：

[0088]

[0089] 其中，a＝nl，n为裂缝的开度与长度之比，无量纲，因此为单根裂缝的裂缝渗吸高度；引入朗伯函数(Lambert函数)W(x)，
可以得到式(6)的显式解：

[0090]

[0091]

[0092] 单根裂缝的t时刻渗吸质量mG为：

[0093]

[0094] 其中，mG的单位为g或kg。

[0095] 图5是本发明实施例中含有大量裂缝的岩心截面示意图。如图5所示，岩心中含有大量天然裂缝，首先在岩心截面选取一个代表性单元Au来研究，假设代表性单元Au内的裂缝长度符合分形分布：

[0096]

[0097] 其中，N(l≥ζ)为长度大于ζ的裂缝数量，lmax为最大裂缝长度，单位为μm；Df为分形维数，无量纲。对式(10)两端关于l求导可以得到下式：

[0098]

[0099] 因此代表性单元Au的总孔隙面积Ap可以被计算得到：

[0100]

[0101] 进而代表性单元Au的面积可以被计算得到：

[0102]

[0103] 其中φf为裂缝孔隙度，单位为％；lmin为最小裂缝长度，单位为μm。岩心横截面Af中的裂缝总数量可以被计算得到：

[0104]

[0105] 其中，Af为岩心横截面，单位为cm2，Au和Ap的单位均为cm2；Nf(≥l)为岩心横截面中长度大于l的裂缝数量。对式(14)关于l求导可得：

[0106]

[0107] 根据式(14)，裂缝长度为l的裂缝数量n(l)可以被计算得到：

[0108]

[0109] 对于整个岩心，长度为l的裂缝渗吸质量为单根裂缝渗吸质量与其数量的乘积。岩心的裂缝渗吸质量由停止渗吸和未停止渗吸两部分组成。渗吸停止时的裂缝渗吸质量为：

[0110]

[0111] 其中，Le为渗吸平衡高度，单位为m，其表达式如下：

[0112]

[0113] 因此，

[0114] 渗吸未停止时的裂缝渗吸质量为：

[0115]

[0116] 因此，

[0117] 则岩心的渗吸质量可以表示为：

[0118]

[0119] 其中，lc为临界裂缝长度，表示t时刻渗吸刚停止的裂缝的长度，单位为μm，其表达式如下：

[0120]

[0121] 储层裂缝渗吸质量预测方法的具体流程如下：

[0122] 1、创建裂缝数量模型；获取第二物性参数；输入各个裂缝长度和第二物性参数至裂缝数量模型中，得到各个裂缝长度的裂缝数量。

[0123] 2、根据各个裂缝的裂缝长度、表面张力和接触角确定各个裂缝的毛细管力。根据各个裂缝在上一渗吸时刻的裂缝渗吸高度和湿相流体密度确定各个裂缝在上一渗吸时刻的重力。

[0124] 3、获取各个裂缝的毛细管力和各个裂缝在上一渗吸时刻的重力。

[0125] 4、当裂缝的毛细管力小于裂缝在上一渗吸时刻的重力时，将该裂缝在上一渗吸时刻的裂缝渗吸高度作为该裂缝在当前渗吸时刻的裂缝渗吸高度。当裂缝的毛细管力等于裂缝在上一渗吸时刻的重力时，根据获取的当前渗吸时刻、第一物性参数和各个裂缝的裂缝长度确定各个裂缝在当前渗吸时刻的裂缝渗吸高度。

[0126] 5、创建第一渗吸质量预测模型和第二渗吸质量预测模型。

[0127] 6、当裂缝在当前渗吸时刻的裂缝渗吸高度小于或等于岩心高度时，输入该裂缝的裂缝长度、该裂缝在当前渗吸时刻的裂缝渗吸高度、该裂缝长度的裂缝数量和湿相流体密度至第一渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量。

[0128] 7、当裂缝在当前渗吸时刻的裂缝渗吸高度大于岩心高度时，输入该裂缝的裂缝长度、岩心高度、该裂缝长度的裂缝数量、湿相流体密度和平均迂曲度至第二渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量。

[0129] 8、根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案。

[0130] 9、通过岩心在当前渗吸时刻的渗吸质量确定岩心在当前渗吸时刻的渗吸饱和度以创建或调整储层的油气资源开发方案。

[0131] 图6是本发明实施例中岩心的渗吸质量和渗吸饱和度随时间的变化曲线示意图。如图6所示，图6中的横坐标为时间，单位为min；纵坐标为渗吸质量和渗吸饱和度，单位分别为g和％。

[0132] 综上，本发明实施例的储层裂缝渗吸质量预测方法先根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度，再根据裂缝渗吸高度和岩心高度的比较结果确定各个裂缝在当前渗吸时刻的渗吸质量，最后根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案，可以预测精确的裂缝渗吸质量，降低预测成本，进一步有效指导了油气资源的开发过程。

[0133] 基于同一发明构思，本发明实施例还提供了一种储层裂缝渗吸质量预测系统，由于该系统解决问题的原理与储层裂缝渗吸质量预测方法相似，因此该系统的实施可以参见方法的实施，重复之处不再赘述。

[0134] 图7是本发明实施例中储层裂缝渗吸质量预测系统的结构框图。如图7所示，储层裂缝渗吸质量预测系统包括：

[0135] 第一获取单元，用于获取各个裂缝的毛细管力和各个裂缝在上一渗吸时刻的重力；

[0136] 裂缝渗吸高度单元，用于根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度；

[0137] 渗吸质量预测模型创建单元，用于创建第一渗吸质量预测模型和第二渗吸质量预测模型；

[0138] 裂缝渗吸质量单元，用于当裂缝在当前渗吸时刻的裂缝渗吸高度小于或等于岩心高度时，输入该裂缝的裂缝长度、该裂缝在当前渗吸时刻的裂缝渗吸高度、该裂缝长度的裂缝数量和湿相流体密度至第一渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量；当裂缝在当前渗吸时刻的裂缝渗吸高度大于岩心高度时，输入该裂缝的裂缝长度、岩心高度、该裂缝长度的裂缝数量、湿相流体密度和平均迂曲度至第二渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量；

[0139] 预测单元，用于根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案。

[0140] 在其中一种实施例中，还包括：

[0141] 毛细管力确定单元，用于根据各个裂缝的裂缝长度、表面张力和接触角确定各个裂缝的毛细管力；

[0142] 重力确定单元，用于根据各个裂缝在上一渗吸时刻的裂缝渗吸高度和湿相流体密度确定各个裂缝在上一渗吸时刻的重力。

[0143] 在其中一种实施例中，裂缝渗吸高度单元具体用于：

[0144] 当裂缝的毛细管力小于裂缝在上一渗吸时刻的重力时，将该裂缝在上一渗吸时刻的裂缝渗吸高度作为该裂缝在当前渗吸时刻的裂缝渗吸高度；

[0145] 当裂缝的毛细管力等于裂缝在上一渗吸时刻的重力时，根据获取的当前渗吸时刻、第一物性参数和各个裂缝的裂缝长度确定各个裂缝在当前渗吸时刻的裂缝渗吸高度。

[0146] 在其中一种实施例中，还包括：

[0147] 裂缝数量模型创建单元，用于创建裂缝数量模型；

[0148] 第二获取单元，用于获取第二物性参数；

[0149] 裂缝数量单元，用于输入各个裂缝长度和第二物性参数至裂缝数量模型中，得到各个裂缝长度的裂缝数量。

[0150] 综上，本发明实施例的储层裂缝渗吸质量预测系统先根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度，再根据裂缝渗吸高度和岩心高度的比较结果确定各个裂缝在当前渗吸时刻的渗吸质量，最后根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案，可以预测精确的裂缝渗吸质量，降低预测成本，进一步有效指导了油气资源的开发过程。

[0151] 本发明实施例还提供能够实现上述实施例中的储层裂缝渗吸质量预测方法中全部步骤的一种计算机设备的具体实施方式。图8是本发明实施例中计算机设备的结构框图，参见图8，所述计算机设备具体包括如下内容：

[0152] 处理器(processor)801和存储器(memory)802。

[0153] 所述处理器801用于调用所述存储器802中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的储层裂缝渗吸质量预测方法中的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：

[0154] 获取各个裂缝的毛细管力和各个裂缝在上一渗吸时刻的重力；

[0155] 根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度；

[0156] 创建第一渗吸质量预测模型和第二渗吸质量预测模型；

[0157] 当裂缝在当前渗吸时刻的裂缝渗吸高度小于或等于岩心高度时，输入该裂缝的裂缝长度、该裂缝在当前渗吸时刻的裂缝渗吸高度、该裂缝长度的裂缝数量和湿相流体密度至第一渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量；

[0158] 当裂缝在当前渗吸时刻的裂缝渗吸高度大于岩心高度时，输入该裂缝的裂缝长度、岩心高度、该裂缝长度的裂缝数量、湿相流体密度和平均迂曲度至第二渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量；

[0159] 根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案。

[0160] 综上，本发明实施例的计算机设备先根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度，再根据裂缝渗吸高度和岩心高度的比较结果确定各个裂缝在当前渗吸时刻的渗吸质量，最后根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案，可以预测精确的裂缝渗吸质量，降低预测成本，进一步有效指导了油气资源的开发过程。

[0161] 本发明实施例还提供能够实现上述实施例中的储层裂缝渗吸质量预测方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的储层裂缝渗吸质量预测方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：

[0162] 获取各个裂缝的毛细管力和各个裂缝在上一渗吸时刻的重力；

[0163] 根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度；

[0164] 创建第一渗吸质量预测模型和第二渗吸质量预测模型；

[0165] 当裂缝在当前渗吸时刻的裂缝渗吸高度小于或等于岩心高度时，输入该裂缝的裂缝长度、该裂缝在当前渗吸时刻的裂缝渗吸高度、该裂缝长度的裂缝数量和湿相流体密度至第一渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量；

[0166] 当裂缝在当前渗吸时刻的裂缝渗吸高度大于岩心高度时，输入该裂缝的裂缝长度、岩心高度、该裂缝长度的裂缝数量、湿相流体密度和平均迂曲度至第二渗吸质量预测模型中，得到该裂缝在当前渗吸时刻的渗吸质量；

[0167] 根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案。

[0168] 综上，本发明实施例的计算机可读存储介质先根据各个裂缝的毛细管力与上一渗吸时刻的重力的比较结果确定各个裂缝在当前渗吸时刻的裂缝渗吸高度，再根据裂缝渗吸高度和岩心高度的比较结果确定各个裂缝在当前渗吸时刻的渗吸质量，最后根据各个裂缝在当前渗吸时刻的渗吸质量预测岩心在当前渗吸时刻的渗吸质量以创建或调整储层的油气资源开发方案，可以预测精确的裂缝渗吸质量，降低预测成本，进一步有效指导了油气资源的开发过程。

[0169] 以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

[0170] 本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrative logical block)，单元，和步骤可以通过电子硬件、电脑软件，或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability)，上述的各种说明性部件(illustrative components)，单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为超出本发明实施例保护的范围。

[0171] 本发明实施例中所描述的各种说明性的逻辑块，或单元，或装置都可以通过通用处理器，数字信号处理器，专用集成电路(ASIC)，现场可编程门阵列或其它可编程逻辑装置，离散门或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置来实现。

[0172] 本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用户终端中的不同的部件中。

[0173] 在一个或多个示例性的设计中，本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。