一种复杂碳酸盐岩储层渗透率的定量表征方法
【技术领域】
[0001] 本发明涉及一种储层渗透率的定量表征方法，尤其涉及一种复杂碳酸盐岩储层渗 透率的定量表征方法。
【背景技术】
[0002] 如图1所示，均质岩石孔隙结构较规则，孔隙度与渗透率之间存在高精度的函数 关系：
[0003]
[0004] 式中，K为渗透率，Φ为孔隙度，f表示函数关系。基于这种高精度的函数关系，可 以采用孔隙度建立储层渗透率的计算公式。
[0005] 另外，基于核磁测井资料，根据自由流体模型可以建立渗透率的计算公式为：
[0006]
[0007] 式中，C为经验常数，FFI为自由流体的孔隙体积，BVI为束缚水的孔隙体积；
[0008] 根据平均T2(即核磁共振测井横向弛豫时间）模型可以建立渗透率的计算公式 为：
[0009]
[0010] 式中，a为经验常数，1^为Τ2分布的几何平均值（以下简称为Τ2几何平均值）。
[0011] 但是，如图2所示，在复杂碳酸盐岩储层中，岩石孔隙结构的复杂性造成储层的强 非均质性，致使孔隙度与渗透率之间不存在高精度的函数关系，如图3所示。因此，采用由 孔隙度建立渗透率计算公式去计算复杂碳酸盐岩储层的渗透率并不可行。另外，在复杂碳 酸盐岩储层中，由于自由流体模型中的经验常数C和平均T2模型中的经验常数a均难以准 确确定，导致式（2)'和式（3)'也难以准确计算复杂碳酸盐岩储层的渗透率。
[0012] 因此，为了准确确定复杂碳酸盐岩储层的渗透率，必须提取能更准确地反映复杂 碳酸盐岩储层渗透性的最佳参数，然后建立该参数与渗透率之间高精度的函数关系，从而 准确确定复杂碳酸盐岩储层的渗透率。
【发明内容】
[0013] 针对上述问题，本发明的目的是提供一种复杂碳酸盐岩储层渗透率的定量表征方 法，能够准确确定复杂碳酸盐岩储层的渗透率。
[0014] 为实现上述目的，本发明采取以下技术方案：一种复杂碳酸盐岩储层渗透率的定 量表征方法，包括以下步骤：
[0015] 1)选取待表征的复杂碳酸盐岩储层位于不同深度的多块岩心样品，分别对各块岩 心样品进行压汞实验，得到不同深度岩心样品的渗透率和毛管压力曲线；然后分别在不同 深度岩心样品的毛管压力曲线上获取各块岩心样品在不同进汞压力下的单位压力进汞饱 和度；
[0016] 2)分别比较不同深度岩心样品在不同进汞压力下的单位压力进汞饱和度，从中选 取出不同深度岩心样品的最大单位压力进汞饱和度；
[0017] 3)根据步骤1)得到的不同深度岩心样品的渗透率和步骤2)得到的不同深度岩心 样品的最大单位压力进汞饱和度，在复杂碳酸盐岩储层中，建立最大单位压力进汞饱和度 与渗透率之间尚精度的函数关系；基于该尚精度的函数关系，建立以最大单位压力进萊饱 和度为参数的复杂碳酸盐岩储层渗透率的计算公式；
[0018] 4)根据待表征的复杂碳酸盐岩储层的测井结果，建立T2几何平均值、孔隙度和最 大单位压力进汞饱和度之间高精度的函数关系；基于该高精度的函数关系，建立全井段不 同深度的最大单位压力进汞饱和度的计算公式；利用步骤3)得到的复杂碳酸盐岩储层渗 透率的计算公式，在全井段计算复杂碳酸盐岩储层的渗透率。
[0019] 所述步骤3)中最大单位压力进汞饱和度与渗透率之间高精度的函数关系为：
[0020] K=f((SHg/Pc)_)
[0021] 式中，K为渗透率，SHg为进汞饱和度，P。为进汞压力，（SHg/h)_为最大单位压力进 汞饱和度，f表示函数关系。
[0022] 所述步骤3)中以最大单位压力进汞饱和度为参数的复杂碳酸盐岩储层渗透率的 计算公式为：
[0023] K=b(SHg/Pc)naxd
[0024] 式中，K为渗透率，（Sh/P。；^为最大单位压力进汞饱和度，b、d均为常数。
[0025] 所述步骤4)中T2几何平均值、孔隙度和最大单位压力进汞饱和度之间高精度的 函数关系为：
[0026]
[0027] 式中，（SHg/Pc)_为最大单位压力进汞饱和度，Τ^为Τ2几何平均值，Φ为孔隙度，f表示函数关系。
[0028] 所述步骤4)中最大单位压力进汞饱和度的计算公式为：
[0029]
[0030] 式中，（SHg/h)_为最大单位压力进汞饱和度，T^为T2几何平均值，为孔隙度，e、g、h、k均为常数。
[0031] 本发明由于采取以上技术方案，其具有以下优点：1、本发明的复杂碳酸盐岩储层 渗透率的定量表征方法，通过提取在非均质岩石中比孔隙度更能准确反映岩石渗透性的参 数一一最大单位压力进汞饱和度（SHg/Pj_，建立新的复杂碳酸盐岩储层渗透率的计算公 式，能够更准确地确定复杂碳酸盐岩储层的渗透率。2、本发明的复杂碳酸盐岩储层渗透率 的定量表征方法，采用最大单位压力进汞饱和度（/匕)_来准确计算复杂碳酸盐岩储层 的渗透率，并结合测井结果，在全井段建立最大单位压力进汞饱和度（/匕)_的计算公 式，避免开展其他大量的岩心实验，在保证复杂碳酸盐岩储层渗透率定量评价结果准确性 的同时，提供了一种简单实用的新方法，能够有效地节约成本，具有较强的经济型。3、本发 明为复杂碳酸盐岩储层渗透率的计算，提供了一种准确率高、通用性强、经济型好的有效方 法，本发明可以广泛应用于各种复杂碳酸盐岩储层渗透率的定量表征。
【附图说明】
[0032] 图1是均质岩石孔隙结构示意图；
[0033] 图2是非均质岩石孔隙结构示意图；
[0034] 图3是复杂碳酸盐岩储层的孔隙度与渗透率的函数关系示意图；
[0035] 图4是复杂碳酸盐岩储层单对数坐标毛管压力曲线示意图；
[0036] 图5是复杂碳酸盐岩储层双对数坐标毛管压力曲线示意图；
[0037] 图6是复杂碳酸盐岩储层的最大单位压力进汞饱和度与渗透率的函数关系示意 图；
[0038] 图7是复杂碳酸盐岩储层的T2几何平均值平方与孔隙度六次方的乘积，和最大单 位压力进汞饱和度的函数关系示意图；
[0039] 图8是复杂碳酸盐岩储层的测井处理成果示意图。
【具体实施方式】
[0040] 下面结合附图和实施例对本发明进行详细的描述。
[0041 ] 如图4~图7所示，本发明提出了 一种复杂碳酸盐岩储层渗透率的定量表征方法， 其包括以下步骤：
[0042] 1)选取待表征的复杂碳酸盐岩储层位于不同深度的多块岩心样品，分别对各块岩 心样品进行压汞实验，得到不同深度岩心样品的渗透率Κ和毛管压力曲线等实验结果；然 后分别在不同深度岩心样品的毛管压力曲线上获取其在不同进汞压力下的单位压力进汞 饱和度sHg/p。，其中，sHgS进汞饱和度，ρ。为进汞压力。
[0043] 2)比较某一深度岩心样品在不同进汞压力下的单位压力进汞饱和度SHg/P。，从中 选取出该深度岩心样品在不同进汞压力下的单位压力进汞饱和度SHg/P。的最大值，即该深 度岩心样品的最大单位压力进汞饱和度（SHg/Pj_;采用上述方法，分别得到不同深度岩心 样品的最大单位压力进汞饱和度（SHg/Pj_。
[0044] 3)根据步骤1)得到的不同深度岩心样品的渗透率K和步骤2)得到的不同深度岩 心样品的最大单位压力进汞饱和度（SHg/P。) _，在复杂碳酸盐岩储层中，建立最大单位压力 进汞饱和度（/匕)_与渗透率K之间高精度的函数关系：
[0045] K=f((SHg/Pc)_) (1)
[0046] 基于该高精度的函数关系，建立以最大单位压力进汞饱和度（SHg/Pj_S参数的 复杂碳酸盐岩储层渗透率K的计算公式。一般情况下，计算公式的主要表现形式为：
[0047] K=b(SHg/Pc)_d (2)
[0048] 式中，b、d均为常数。
[0049] 但需要注意的是，上述复杂碳酸盐岩储层渗透率K的计算公式是以最大单位压力 进汞饱和度（sHg/pj_S参数的，当需要确定某一地层深度的渗透率κ时，首先需要确定该 地层深度处的最大单位压力进汞饱和度（SHg/Pj_，而最大单位压力进汞饱和度（SHg/Pj_ 需要通过对岩心样品进行压汞实验才能得到；对于没有取得岩心样品的地层深度，则得不 到其最大单位压力进汞饱和度（SHg/Pj_，导致无法确定渗透率K。因此，我们需要将上述复 杂碳酸盐岩储层渗透率K的计算公式推广至全井段范围。
[0050] 4)根据待表征的复杂碳酸盐岩储层的测井结果，建立T2几何平均值T2011、孔隙度# 和最大单位压力进汞饱和度（sHg/pj_2间高精度的函数关系：
[0051]
[0052] 基于该高精度的函数关系，建立最大单位压力进汞饱和度（SHg/Pj_的计算公式， 从而可以计算全井段不同深度的最大单位压