准致密储层快速评价方法及准致密储层多指标评价方法

准致密储层快速评价方法及准致密储层多指标评价方法
【专利摘要】本发明公开了一种准致密储层快速评价方法及准致密储层多指标评价方法。所述准致密储层快速评价方法包括：根据岩心资料和地质资料，得出不同储层深度下的孔隙度、含水饱和度、冲洗带饱和度以及厚度；根据每个储层深度下的含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中可动烃体积与水体积的比值；根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度以及岩石孔隙中可动烃体积与水体积的比值，确定该储层深度下的储层流体识别指标的大小；以及根据每个储层深度下的储层流体识别指标的大小，识别该储层深度下的流体性质，以对所述储层的流体性质进行评价。由此，可以提高储层评价的快速性、准确性和可靠性。
【专利说明】
准致密储层快速评价方法及准致密储层多指标评价方法
技术领域
[0001] 本发明涉及储层评价领域，具体地，涉及一种准致密储层快速评价方法及准致密 储层多指标评价方法。
【背景技术】
[0002] 准致密碎肩岩储层具有低孔渗、高束缚水、连通性差、非均质性强的特征，储层评 价水平是有效开发该储层的关键因素。测井评价作为三种核心储层评价手段之一，其定性 评价的可靠性和定量评价的精度，是决定储层评价水平的基础和关键。但准致密储层测井 评价一直是储层测井评价面临的技术难点。由于准致密储层超低的储集空间，孔隙和孔隙 流体对测井响应贡献小，储层非均性强，测井响应差异大，因此，传统、常规的储层测井评价 方法、评价指标已经不能满足准致密储层测井评价可靠性和精度要求，制约了储层评价水 平的提高和发展，并且制约该类储层的有效开发。

【发明内容】

[0003] 本发明的目的是提供一种准致密储层快速评价方法及准致密储层多指标评价方 法，所述评价方法能够基于新构建的评价指标来对储层进行评价，从而提高储层评价的快 速性、准确性和可靠性。
[0004] 为了实现上述目的，本发明提供一种准致密储层快速评价方法，该方法包括：根 据岩心资料和地质资料，得出不同储层深度下的孔隙度、含水饱和度、冲洗带饱和度以及厚 度；根据每个储层深度下的含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中 可动烃体积与水体积的比值；根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、 厚度以及岩石孔隙中可动烃体积与水体积的比值，确定该储层深度下的储层流体识别指标 的大小；以及根据每个储层深度下的储层流体识别指标的大小，识别该储层深度下的流体 性质，以对所述储层的流体性质进行评价。
[0005] 优选地，通过以下方式确定所述岩石孔隙中可动烃体积与水体积的比值：
[0007] 通过以下方式确定所述储层流体识别指标大小：
[0009] 其中，GW&表示第i个储层深度下的岩石孔隙中可动烃体积与水体积的比值；Sw i 表示第i个储层深度下的含水饱和度；Sx〇1表示第i个储层深度下的冲洗带饱和度；FLQI i 表示第i个储层深度下的储层流体识别指标；PHIEi表示第i个储层深度下的孔隙度；以及 呒表示第i个储层深度下的厚度。
[0010] 优选地，根据每个储层深度下的储层流体识别指标的大小，识别该储层深度下的 流体性质的步骤包括：在FLQI1= 0. 00的情况下，识别第i个储层深度下的流体性质为干 层；在FLQIf 0. 50的情况下，识别第i个储层深度下的流体性质为水层；在FLQI 1. 00 的情况下，识别第i个储层深度下的流体性质为含水含气层；在FLQI1= 1. 50的情况下，识 别第i个储层深度下的流体性质为含气层；在FLQI1= 2. 00的情况下，识别第i个储层深 度下的流体性质为气水同层；在FLQI1= 2. 50的情况下，识别第i个储层深度下的流体性 质为含水气层；在FLQI1= 3. 00的情况下，识别第i个储层深度下的流体性质为差气层；以 及在FLQIi= 3. 50的情况下，识别第i个储层深度下的流体性质为气层。
[0011] 本发明还提供一种准致密储层多指标评价方法，该方法包括：根据本发明提供的 所述准致密储层快速评价方法，确定每个储层深度下的储层流体识别指标的大小；根据岩 心资料和地质资料，得出每个储层深度下的泥质含量和渗透率；对测井资料进行环境校正， 得出所述储层的测井曲线，其中，该测井曲线包括由每个储层深度下的深电阻率测井值形 成的深电阻率曲线；根据每个储层深度下的孔隙度、含水饱和度和冲洗带饱和度，确定该储 层深度下的岩石孔隙中可动烃的体积、以及可动烃体积与岩石孔隙体积的比值；根据每个 储层深度下的孔隙度和深电阻率测井值，确定该储层深度下的储层流体性质指标的大小； 根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度、泥质含量、渗透率、岩石 孔隙中可动烃的体积、可动烃体积与岩石孔隙体积的比值、岩石孔隙中可动烃体积与水体 积的比值、以及储层流体性质指标，构建该储层深度下的多指标识别图版；以及根据每个储 层深度下的储层流体识别指标的大小、以及每个储层深度下的多指标识别图版，共同识别 该储层深度下的流体性质，以对所述储层的流体性质进行评价。
[0012] 优选地，通过以下方式确定所述岩石孔隙中可动烃的体积、以及所述可动烃体积 与岩石孔隙体积的比值：
[0013] BVWI^ PHIE ^Sxo-PHIE^Sw,
[0015] 其中，BVWL表示第i个储层深度下的岩石孔隙中可动烃的体积；SGM i表示第i个 储层深度下的可动烃体积与岩石孔隙体积的比值；PHIEi表示第i个储层深度下的孔隙度； SWl表示第i个储层深度下的含水饱和度；以及Sxo i表示第i个储层深度下的冲洗带饱和 度。
[0016] 优选地，根据每个储层深度下的孔隙度和深电阻率测井值，确定该储层深度下的 储层流体性质指标的大小的步骤包括：构建所述储层的孔隙度-电阻率交会图版；将每个 储层深度下的孔隙度和深电阻率测井值映射到所述孔隙度-电阻率交会图版中，形成每个 储层深度下的孔隙度-电阻率坐标点；以及根据每个储层深度下的孔隙度-电阻率坐标 点在所述孔隙度-电阻率交会图版中的位置，确定该储层深度下的储层流体性质指标的大 小。
[0017] 优选地，所述测井曲线还包括由每个储层深度下的自然伽玛测井值形成的自然伽 玛曲线和由每个储层深度下的声波测井值形成的声波曲线；以及所述准致密储层多指标评 价方法还包括：对所述自然伽玛曲线和所述声波曲线进行归一化；判断所述测井曲线中是 否包括由每个储层深度下的中子测井值形成的中子曲线和由每个储层深度下的密度测井 值形成的密度曲线；在判定所述测井曲线中未包括所述中子曲线的情况下，利用归一化的 自然伽玛曲线和声波曲线，重构归一化的中子曲线；在判定所述测井曲线中未包括所述密 度曲线的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的密度曲线；在判定 所述测井曲线中包括所述中子曲线的情况下，判断该中子曲线是否合格，并在判定所述中 子曲线不合格的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的中子曲线， 以及在判定所述中子曲线合格的情况下，对所述中子曲线进行归一化；在判定所述测井曲 线中包括所述密度曲线的情况下，判断该密度曲线是否合格，并在判定所述密度曲线不合 格的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的密度曲线，以及在判定 所述密度曲线合格的情况下，对所述密度曲线进行归一化。
[0018] 优选地，通过以下方式重构归一化的中子曲线和归一化的密度曲线：
[0019] NCNLri= 0. 5707*NAC ,+0. 4114^6^-0. 001
[0020] NDENri= 0. 1092-0. 6922*NAC ,+0. 1254^6?
[0021] 其中，NAQ表示归一化的声波曲线中第i个储层深度下的声波测井值；NGR i表示 归一化的自然伽玛曲线中第i个储层深度下的自然伽玛测井值；NCNL"表示重构的归一化 的中子曲线中第i个储层深度下的中子测井值；以及NDEN"表示重构的归一化的密度曲线 中第i个储层深度下的密度测井值。
[0022] 优选地，该准致密储层多指标评价方法还包括：在重构归一化的中子曲线之后，利 用所重构的归一化的中子曲线重构中子曲线；和/或在重构归一化的密度曲线之后，利用 所重构的归一化的密度曲线重构密度曲线。
[0023] 优选地，通过以下方式重构所述中子曲线和所述密度曲线：
[0024] CNLri= 0. 1566+0. 0639*NCNL ri
[0025] DENri= 2· 4929+0. 1305*NDEN ri
[0026] 其中，CNL"表示重构的中子曲线中第i个储层深度下的中子测井值；DEN"表示重 构的密度曲线中第i个储层深度下的密度测井值；NCNL"表示重构的归一化的中子曲线中 第i个储层深度下的中子测井值；以及NDEN"表示重构的归一化的密度曲线中第i个储层 深度下的密度测井值。
[0027] 优选地，所述测井曲线还包括由每个储层深度下的浅电阻率测井值形成的浅电阻 率曲线；以及所述准致密储层多指标评价方法还包括：对所述深电阻率曲线和所述浅电阻 率曲线进行归一化；根据归一化的自然伽玛曲线、归一化的声波曲线、归一化的深电阻率曲 线、以及每个储层深度下的泥质含量，确定该储层深度下的高角度缝指标的大小；根据归一 化的深电阻率曲线、归一化的浅电阻率曲线、未经归一化的深电阻率曲线、未经归一化的浅 电阻率曲线、以及每个储层深度下的泥质含量，确定该储层深度下的低角度缝指标的大小； 以及根据每个储层深度下的高角度缝指标的大小和低角度缝指标的大小，对该储层深度下 的裂缝发育情况进行评价，以对所述储层的裂缝发育情况进行评价。
[0028] 优选地，通过以下方式确定所述高角度缝指标的大小：
[0029]
[0030] 通过以下方式确定所述低角度缝指标的大小：
[0031]
[0032] 其中，VFRL表示第i个储层深度下的高角度缝指标；VCL i表示第i个储层深度下 的泥质含量；DMaXl表示第i个储层深度下，归一化的自然伽玛曲线的一阶导数、归一化的 声波曲线的一阶导数、归一化的深电阻率曲线的一阶导数中的最大值；DMi ni表示第i个储 层深度下，归一化的自然伽玛曲线的一阶导数、归一化的声波曲线的一阶导数、以及归一化 的深电阻率曲线的一阶导数中导数间差值的最小值；HFRL表示第i个储层深度下的低角 度缝指标；NRDi表示归一化的深电阻率曲线中第i个储层深度下的深电阻率测井值；NRS i 表示归一化的浅电阻率曲线中第i个储层深度下的浅电阻率测井值;AvLRDi表示未经归一 化的深电阻率曲线中所有深电阻率测井值的对数均值；以及AvLRS^示未经归一化的浅电 阻率曲线中所有浅电阻率测井值的对数均值。
[0033] 优选地，该准致密储层多指标评价方法还包括：根据归一化的声波曲线、归一化 的中子曲线和每个储层深度下的泥质含量，确定该储层深度下的天然气挖掘效应指标的大 小；以及根据每个储层深度下的天然气挖掘效应指标的大小，对该储层深度下的流体性质 进行定性识别，以对所述储层的流体性质进行评价。
[0034] 优选地，通过以下方式确定所述天然气挖掘效应指标的大小：
[0035]
[0036] 其中，CAVL表示第i个储层深度下的天然气挖掘效应指标；VCL i表示第i个储层 深度下的泥质含量；NAQ表示归一化的声波曲线中第i个储层深度下的声波测井值；以及 NCNQ表示归一化的中子曲线中第i个储层深度下的中子测井值。
[0037] 在上述技术方案中，通过构建储层流体识别指标，并基于该储层流体识别指标的 大小来定量评价储层的流体性质，能够实现对准致密储层的快速评价，并且能够提高准致 密储层评价的可靠性和精度。此外，还可以通过构建的多种指标来构建多指标识别图版，并 基于所述储层流体识别指标的大小和所述多指标识别图版，来对所述储层的流体性质进行 综合定量评价，这样，可以排除准致密储层评价中的多解性，进一步提高准致密储层评价的 可靠性和精度。通过实例验证，本发明提供的储层评价方法所得到的评价准确度能够达到 90%~95%。由此，本发明提供的准致密储层评价方法能够满足准致密储层测井评价可靠 性和精度的要求，有利于储层评价水平的提高和发展，保证该类储层的有效开发，并且能够 为后续地质、地球物理综合研究奠定坚实基础。
[0038] 本发明的其他特征和优点将在随后的【具体实施方式】部分予以详细说明。
【附图说明】
[0039] 附图是用来提供对本发明的进一步理解，并且构成说明书的一部分，与下面的具 体实施方式一起用于解释本发明，但并不构成对本发明的限制。在附图中：
[0040] 图1示出了根据本发明的实施方式的准致密储层快速评价方法的流程图；
[0041] 图2示出了根据本发明的实施方式的准致密储层多指标评价方法的流程图；
[0042] 图3示出了根据本发明的实施方式的孔隙度-电阻率交会图版的示例；
[0043] 图4示出了根据本发明的实施方式的多指标识别图版的示例；以及
[0044] 图5和图6分别示出了利用本发明提供的准致密储层多指标评价方法对四川盆地 某井区进行实际评价所得的结果。
【具体实施方式】
[0045] 以下结合附图对本发明的【具体实施方式】进行详细说明。应当理解的是，此处所描 述的【具体实施方式】仅用于说明和解释本发明，并不用于限制本发明。
[0046] 图1示出了根据本发明的实施方式的准致密储层快速评价方法的流程图。如图 1所示，该方法可以包括：步骤S11，根据岩心资料和地质资料，得出不同储层深度下的孔隙 度、含水饱和度、冲洗带饱和度以及厚度；步骤S12,根据每个储层深度下的含水饱和度和 冲洗带饱和度，确定该储层深度下的岩石孔隙中可动烃体积与水体积的比值；步骤S13,根 据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度以及岩石孔隙中可动烃体积 与水体积的比值，确定该储层深度下的储层流体识别指标的大小；以及步骤S14,根据每个 储层深度下的储层流体识别指标的大小，识别该储层深度下的流体性质，以对所述储层的 流体性质进行评价。
[0047] 具体地，首先在步骤S11，根据岩心资料和地质资料，得出不同储层深度下的孔隙 度PHIE、含水饱和度Sw、冲洗带饱和度Sxo以及厚度H。该步骤的实现是本领域技术人员公 知的，对此，本发明在此不进行详细描述。
[0048] 在得出不同储层深度下的孔隙度PHIE、含水饱和度Sw、冲洗带饱和度Sxo以及厚 度Η之后，可以根据每个储层深度下的含水饱和度Sw和冲洗带饱和度Sxo，确定出该储层深 度下的岩石孔隙中可动烃体积与水体积的比值GWR。该比值是本发明新构建的一项指标，和 下文描述的岩石孔隙中可动烃的体积BVWI和可动烃体积与岩石孔隙体积的比值SGM这两 项指标共同构成"相对体积指标"，该组指标用于后续的储层评价。
[0049] 可以通过以下等式（1)来确定所述岩石孔隙中可动烃体积与水体积的比值：
[0051] 其中，GW&表示第i个储层深度下的岩石孔隙中可动烃体积与水体积的比值； SWl表示第i个储层深度下的含水饱和度；Sxo i表示第i个储层深度下的冲洗带饱和度； 1 < i < N，N表示储层的采样点总数，其中，每个采样点对应一个储层深度，第i个采样点 也等同表示第i个储层深度。
[0052] 之后，在确定出每个储层深度下的岩石孔隙中可动烃体积与水体积的比值之后， 就可以结合该比值、以及每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度和厚度，来 共同确定出每个储层深度下的储层流体识别指标的大小。具体地，可通过以下等式（2)来 确定所述储层流体识别指标的大小：
[0053]
[0054] 其中，FLQL表示第i个储层深度下的储层流体识别指标；PHIE i表示第i个储层 深度下的孔隙度；以及私表示第i个储层深度下的厚度。
[0055] 储层流体识别指标FLQI同样是本发明新构建的一项指标，其用于对储层进行快 速、初步自动评价，特别是针对储层的流体性质进行快速、初步自动评价，如下面进一步描 述的。
[0056] 在FLQI1= 0. 00的情况下，识别第i个储层深度下的流体性质为干层；在FLQI 1 = 〇. 50的情况下，识别第i个储层深度下的流体性质为水层；在FLQI1= 1. 00的情况下，识别 第i个储层深度下的流体性质为含水含气层；在FLQI1= 1. 50的情况下，识别第i个储层 深度下的流体性质为含气层；在FLQI1= 2. 00的情况下，识别第i个储层深度下的流体性 质为气水同层；在FLQI1= 2. 50的情况下，识别第i个储层深度下的流体性质为含水气层； 在FLQI1= 3. 00的情况下，识别第i个储层深度下的流体性质为差气层；以及在FLQI 1 = 3. 50的情况下，识别第i个储层深度下的流体性质为气层。
[0057] 通过上述方式，可以快速识别出每个储层深度下的流体性质，通过将这些识别结 果进行组合，可以对所述储层的流体性质做出整体评价，分析出为干层、水层、含水含气层、 含气层、气水同层、含水气层、差气层和气层的各储层段，以为日后的储层开发、地质研究奠 定基础。
[0058] 通过本发明构建的储层流体识别指标来快速评价储层的流体性质，可以实现对准 致密储层的流体性质的自动、快速地定量评价，具有较高的准确性和可靠性，填补了现有技 术中缺乏对准致密储层的有效评价方法的空白。该指标可以用于对储层的流体性质的预评 价，有利于快速、准确评价储层的流体性质。此外，该指标还可以与其他定量评价指标一起， 共同对储层的流体性质进行定量评价，以实现基于多指标的准致密储层评价方法，进一步 提高评价的可靠性和准确性。下文将具体描述这一准致密储层多指标评价方法。
[0059] 图2示出了根据本发明的实施方式的准致密储层多指标评价方法的流程图。如图 2所示，该方法可以包括：步骤S21，根据本发明提供的上述准致密储层快速评价方法，确定 每个储层深度下的储层流体识别指标的大小；步骤S22,根据岩心资料和地质资料，得出每 个储层深度下的泥质含量和渗透率；步骤S23,对测井资料进行环境校正，得出所述储层的 测井曲线，其中，该测井曲线包括由每个储层深度下的深电阻率测井值形成的深电阻率曲 线；步骤S24,根据每个储层深度下的孔隙度、含水饱和度和冲洗带饱和度，确定该储层深 度下的岩石孔隙中可动烃的体积、以及可动烃体积与岩石孔隙体积的比值；步骤S25,根据 每个储层深度下的孔隙度和深电阻率测井值，确定该储层深度下的储层流体性质指标的大 小；步骤S26,根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度、泥质含量、 渗透率、岩石孔隙中可动烃的体积、可动烃体积与岩石孔隙体积的比值、岩石孔隙中可动烃 体积与水体积的比值、以及储层流体性质指标的大小，构建该储层深度下的多指标识别图 版；以及步骤S27,根据每个储层深度下的储层流体识别指标的大小、以及每个储层深度下 的多指标识别图版，共同识别该储层深度下的流体性质，以对所述储层的流体性质进行评 价。
[0060] 具体地，首先，根据上面结合图1所描述的准致密储层快速评价方法，确定出每个 储层深度下的储层流体识别指标的大小，即FLQL。此外，除了根据岩心资料和地质资料，得 出不同储层深度下的孔隙度、含水饱和度、冲洗带饱和度以及厚度之外，还要根据该岩心资 料和地质资料，得出每个储层深度下的泥质含量和渗透率。该步骤的实现是本领域技术人 员公知的，对此，本发明在此不进行详细描述。
[0061] 接下来，在步骤S23,对测井资料进行环境校正，得出所述储层的测井曲线，其中， 该测井曲线可以包括由每个储层深度下的深电阻率测井值形成的深电阻率曲线。对测井资 料进行环境校正，是为了排除井径、泥浆等非地层因素对测井资料的影响。如何对测井资料 进行环境校正的具体方式是本领域的技术人员公知的，对此，本发明在此不做具体描述。
[0062] 接下来，在步骤S24,根据每个储层深度下的孔隙度、含水饱和度和冲洗带饱和度， 确定该储层深度下的岩石孔隙中可动烃的体积、以及可动烃体积与岩石孔隙体积的比值， 即上文提到的"相对体积指标"中的另外两项指标BVWI和SGM。可以通过以下等式（3)和 (4)来分别确定这两项指标：
[0063] BVWI^ PHIE ^Sxo-PHIE^Sw, (3)
[0065] 其中，BVWL表示第i个储层深度下的岩石孔隙中可动烃的体积；SGM i表示第i个 储层深度下的可动烃体积与岩石孔隙体积的比值；PHIEi表示第i个储层深度下的孔隙度； SWl表示第i个储层深度下的含水饱和度；以及Sxo i表示第i个储层深度下的冲洗带饱和 度。如此，上述等式（1)还可以被扩展成以下等式（5):
[0067] 所述"相对体积指标"中的每个指标作为后面构建多指标识别图版中的三项定量 指标。通过新增这三项定量指标来进行储层评价，可以克服准致密储层评价定量参数误差 大的难题。
[0068] 虽然使用了 S22、S23和S24这样的表述，但是并不意味着这三个步骤是有先后次 序的，这三个步骤可以是同时并行的。
[0069] 接下来，在步骤S25,根据每个储层深度下的孔隙度和深电阻率测井值，确定该储 层深度下的储层流体性质指标的大小。具体地，该步骤S25可以包括如下几个步骤：构建所 述储层的孔隙度-电阻率交会图版；将每个储层深度下的孔隙度和深电阻率测井值映射到 所述孔隙度-电阻率交会图版中，形成每个储层深度下的孔隙度-电阻率坐标点。图3示 出了所述孔隙度-电阻率交会图版的示例。在该交会图版中，深电阻率测井值RD作为纵坐 标，孔隙度PHIE作为横坐标。每个储层深度下的孔隙度和深电阻率测井值被映射到该孔隙 度-电阻率交会图版中，这样，就形成每个储层深度下的孔隙度-电阻率坐标点。如上所述， 总共有N个储层深度，因此，可以形成N个孔隙度-电阻率坐标点。
[0070] 可以根据先验数据或经验，在该交会图版中预先划分出若干区域。在图3所示的 交会图版中，通过5条线将交会图版划分为6个非线性区域，分别为区域I、区域II、区域 III、区域IV、区域V以及区域VI。每个区域代表一个储层流体性质指标IGAS的值。这样， 就可以根据每个储层深度下的孔隙度-电阻率坐标点在所述孔隙度-电阻率交会图版中的 位置（即，落到哪个区域），来确定该储层深度下的储层流体性质指标的大小。
[0071] 例如，针对第i个储层深度，如果其孔隙度-电阻率坐标点落在区域I，则确定该第 i个储层深度下，IGAS1= 4. 5 ;如果其孔隙度-电阻率坐标点落在区域II，则确定该第i个 储层深度下，IGAS1= 3. 5 ;如果其孔隙度-电阻率坐标点落在区域III，则确定该第i个储 层深度下，IGAS1= 2. 5 ;如果其孔隙度-电阻率坐标点落在区域IV，则确定该第i个储层深 度下，IGAS1= 1. 5 ;以及如果其孔隙度-电阻率坐标点落在区域V或区域VI，则确定该第i 个储层深度下，IGASi= 0.5。通过这种方式，可以确定出每个储层深度下的储层流体性质 指标的大小。该储层流体性质指标也作为后面构建多指标识别图版中的一项定量指标。
[0072] 接下来，在步骤S26中，构建每个储层深度下的多指标识别图版。具体地，在本 发明中，通过10种指标来构建多指标识别图版。例如，如图4所示，这10种指标分别是： 含水饱和度（Sw);孔隙度（PHIE);孔隙度与厚度的乘积（PHIE*H);孔隙度、冲洗带饱和度 与厚度的乘积（PHIE*Swo*H);泥质含量（VCL);渗透率（PERM);岩石孔隙中可动烃的体积 (BVMI);可动烃体积与岩石孔隙体积的比值（SGM);岩石孔隙中可动烃体积与水体积的比 值（GWR);以及储层流体性质指标（IGAS)。
[0073] 在本发明中，储层的流体性质可以分为八种，分别是：干层、水层、含水含气层、含 气层、气水同层、含水气层、差气层和气层。可以根据先验数据构建出每种流体性质对应的 标准多指标识别图版。之后，将新构建的一储层深度下的多指标识别图版与这八种标准多 指标识别图版进行匹配，确定与之最接近的标准多指标识别图版，并最终识别出该储层深 度下的流体性质。
[0074] 在本发明中，可以结合使用图1中描述的方法确定出的每个储层深度下的储层流 体识别指标的大小、以及每个储层深度下的多指标识别图版，来共同识别该储层深度下的 流体性质，以对所述储层的流体性质做出综合评价。
[0075] 通过上述准致密储层多指标评价方法，可以融合多种定量指标来对准致密储层进 行综合评价，相比于仅采用一种指标进行评价，可以有效消除准致密储层评价的多解性，从 而可以进一步提高准致密储层评价结果的可靠性和准确性。
[0076] 此外，在上面通过步骤S23获得的测井曲线中，不仅可以包括所述深电阻率曲线， 还可以包括由每个储层深度下的自然伽玛测井值形成的自然伽玛曲线和由每个储层深度 下的声波测井值形成的声波曲线。并且，所述准致密储层多指标评价方法还可以包括：对所 述自然伽玛曲线和所述声波曲线进行归一化；判断所述测井曲线中是否包括由每个储层深 度下的中子测井值形成的中子曲线和由每个储层深度下的密度测井值形成的密度曲线；在 判定所述测井曲线中未包括所述中子曲线的情况下，利用归一化的自然伽玛曲线和声波曲 线，重构归一化的中子曲线；在判定所述测井曲线中未包括所述密度曲线的情况下，利用归 一化的自然伽玛曲线和声波曲线，重构归一化的密度曲线；在判定所述测井曲线中包括所 述中子曲线的情况下，判断该中子曲线是否合格，并在判定所述中子曲线不合格的情况下， 利用归一化的自然伽玛曲线和声波曲线，重构归一化的中子曲线，以及在判定所述中子曲 线合格的情况下，对所述中子曲线进行归一化；在判定所述测井曲线中包括所述密度曲线 的情况下，判断该密度曲线是否合格，并在判定所述密度曲线不合格的情况下，利用归一化 的自然伽玛曲线和声波曲线，重构归一化的密度曲线，以及在判定所述密度曲线合格的情 况下，对所述密度曲线进行归一化。
[0077] 通常情况下，测井曲线中要包括所述中子曲线和密度曲线。中子曲线和密度曲线 可用于辅助泥质含量、孔隙度、渗透率等参数的计算。然而，密度曲线和中子曲线数据的采 集对井况条件要求较高，有时可能采集不到密度曲线和中子曲线，或者采集到的密度曲线 和中子曲线质量差，会影响后期的储层评价。这就导致测井曲线中可能缺失中子曲线和/ 或密度曲线，或者虽然包括中子曲线和/或密度曲线，但是所包括的中子曲线和/或密度曲 线不合格。在这种情况下，为了不影响其他参数的计算及后续的储层评价，可以按照上述方 法，对缺失的或不合格的中子曲线和/或密度曲线进行重构。
[0078] 具体地，可以通过以下等式（6)来重构归一化的中子曲线，以及通过以下等式（7) 来重构归一化的密度曲线：
[0079] NCNLri= 0· 5707*NAC，0· 4114*吣民-0· 001 (6)
[0080] NDENri= 0. 1092-0. 6922*NAC ,+0. 1254^6? (7)
[0081] 其中，NAQ表示归一化的声波曲线中第i个储层深度下的声波测井值；NGR i表示 归一化的自然伽玛曲线中第i个储层深度下的自然伽玛测井值；NCNL"表示重构的归一化 的中子曲线中第i个储层深度下的中子测井值；以及NDEN"表示重构的归一化的密度曲线 中第i个储层深度下的密度测井值。
[0082] 由于采用归一化的自然伽玛曲线和归一化的声波曲线来进行重构，因此，通过上 面等式（6)和（7)重构出的中子曲线和密度曲线也是归一化的中子曲线和归一化的密度曲 线。
[0083] 如上所述，中子曲线和密度曲线可用于其他定量指标的计算。因此，在重构出归一 化的中子曲线和/或归一化的密度曲线之后，还可以利用所重构的归一化的中子曲线重构 中子曲线；和/或利用所重构的归一化的密度曲线重构密度曲线。
[0084] 具体地，可通过以下等式（8)重构所述中子曲线，以及通过等式（9)重构所述密度 曲线：
[0085] CNLri= 0· 1566+0. 0639*NCNLri (8)
[0086] DENri= 2· 4929+0. 1305*NDENri (9)
[0087] 其中，CNL"表示重构的中子曲线中第i个储层深度下的中子测井值；以及DEN"表 示重构的密度曲线中第i个储层深度下的密度测井值。
[0088] 这样，就可以在无法采集中子曲线和密度曲线、或者采集到的中子曲线和密度曲 线质量不合格的情况下，重构出中子曲线和密度曲线。由此，可以满足大斜度、水平井以及 多井综合评价的资料需求。
[0089] 如何判断中子曲线和密度曲线是否合格的方法是本领域的技术人员公知的，例 如，通过确定密度曲线、中子曲线分别与自然伽玛曲线、深电阻率曲线、声波曲线等之间的 相关性差，来判断密度曲线、中子曲线是否合格（相关性差过大，判断曲线不合格）。因此， 对于曲线是否合格的判断方法本发明不做具体描述。
[0090] 另外，在所述中子曲线和密度曲线未缺失、并且合格的情况下，对所述中子曲线和 密度曲线进行归一化。对自然伽玛曲线、声波曲线、中子曲线和密度曲线进行归一化的方法 有很多，并且均是本领域的技术人员公知的。在本发明的一个示例实施方式中，可以采用以 下等式（10)来进行归一化：
[0092] 其中，logi表示曲线上第i个测井值；Nlog ;表示该第i个测井值的归一化值；Αν log表示曲线上所有测井值的平均值；以及Std log表示标准差。
[0093] 上述示例归一化方法不用于限制本发明，本领域的技术人员也可以采用其他归一 化方法来对自然伽玛曲线、声波曲线、中子曲线和密度曲线进行归一化，对此，本发明在此 不--具体阐述。
[0094] 在本发明的另一实施方式中，所述测井曲线还可以包括由每个储层深度下的浅电 阻率测井值形成的浅电阻率曲线；以及所述准致密储层多指标评价方法还可以包括：对所 述深电阻率曲线和所述浅电阻率曲线进行归一化；根据归一化的自然伽玛曲线、归一化的 声波曲线、归一化的深电阻率曲线、以及每个储层深度下的泥质含量，确定该储层深度下的 高角度缝指标的大小；根据归一化的深电阻率曲线、归一化的浅电阻率曲线、未经归一化的 深电阻率曲线、未经归一化的浅电阻率曲线、以及每个储层深度下的泥质含量，确定该储层 深度下的低角度缝指标的大小；以及根据每个储层深度下的高角度缝指标的大小和低角度 缝指标的大小，对该储层深度下的裂缝发育情况进行评价，以对所述储层的裂缝发育情况 进行评价。
[0095] 在对准致密储层进行评价时，储层的裂缝发育情况也是一项评价指标。在本发明 中，通过构造高角度缝指标和低角度缝指标，可以快速且准确地定性分析出储层的裂缝发 育情况。
[0096] 具体地，可通过以下等式（11)确定所述高角度缝指标的大小：
[0097]
[0098] 其中，VFRL表示第i个储层深度下的高角度缝指标；VCL i表示第i个储层深度下 的泥质含量；DMaXl表示第i个储层深度下，归一化的自然伽玛曲线的一阶导数、归一化的 声波曲线的一阶导数、归一化的深电阻率曲线的一阶导数中的最大值；DMi ni表示第i个储 层深度下，归一化的自然伽玛曲线的一阶导数、归一化的声波曲线的一阶导数、以及归一化 的深电阻率曲线的一阶导数中导数间差值的最小值。
[0099] 此外，可通过以下等式（12)确定所述低角度缝指标的大小：
[0101] 其中，HFRL表示第i个储层深度下的低角度缝指标；NRD 1表示归一化的深电阻率 曲线中第i个储层深度下的深电阻率测井值；NRSi表示归一化的浅电阻率曲线中第i个储 层深度下的浅电阻率测井值；AvLRDi表示未经归一化的深电阻率曲线中所有深电阻率测井 值的对数均值；以及AvLRSi表示未经归一化的浅电阻率曲线中所有浅电阻率测井值的对数 均值。
[0102] 之后，就可以根据所得到的高角度缝指标的大小和低角度缝指标的大小，对每个 储层深度下的裂缝发育情况进行定性分析。例如，在VFRI 1= 1. 0和/或HFRI 1. 0的情 况下，可以定性分析出该第i个储层深度下的裂缝发育；而在VFRI1= 0. 1且HFRI 0. 1 的情况下，可以定性分析出该第i个储层深度下的裂缝未发育。由此，可以分析出整个准致 密储层的裂缝发育情况，以为日后的储层开发、地质研究奠定基础。
[0103] 在本发明的另一实施方式中，所述准致密储层多指标评价方法还可以包括：根据 归一化的声波曲线、归一化的中子曲线和每个储层深度下的泥质含量，确定该储层深度下 的天然气挖掘效应指标的大小；以及根据每个储层深度下的天然气挖掘效应指标的大小， 对该储层深度下的流体性质进行定性识别，以对所述储层的流体性质进行评价。
[0104] 具体地，可通过以下等式（13)确定所述天然气挖掘效应指标的大小：
[0105]
[0106] 其中，CAVL表示第i个储层深度下的天然气挖掘效应指标；VCL i表示第i个储层 深度下的泥质含量；NAQ表示归一化的声波曲线中第i个储层深度下的声波测井值；以及 NCNQ表示归一化的中子曲线中第i个储层深度下的中子测井值。
[0107] 如果是经重构得出的归一化的中子曲线，则NCNQ等同于NCNL"。如果不是经重构 得出的归一化的中子曲线，那么NCNQ则代表对测井曲线中包括的原本中子曲线进行归一 化后的结果。
[0108] 区别于上面通过储层流体识别指标FLQI以及多指标识别图版来对储层的流体性 质进行定量识别，在此种实施方式中，通过每个储层深度下的天然气挖掘效应指标CAVI的 大小，可以实现对每个储层深度下的流体性质的定性识别。在储层评价体系中，不仅需要定 量识别，也需要定性识别。通过本发明构造出的天然气挖掘效应指标，可以提高储层流体性 质的定性识别的可靠性。例如，在CAVI 1= 1. 0的情况下，可以定性识别出该第i个储层深 度下含气；而在CAVI1= 0. 1的情况下，可以定性识别出该第i个储层深度下不含气。
[0109] 图5和图6分别示出了利用本发明提供的准致密储层多指标评价方法对四川盆地 某井区进行实际评价所得的结果。经过与实际测试结果进行比较，通过本发明提供的准致 密储层快速评价方法及准致密储层多指标评价方法，能够实现对准致密储层的准确评价。 [0110] 综上所述，在本发明提供的准致密储层快速评价方法及准致密储层多指标评价方 法中，通过构建储层流体识别指标，并基于该储层流体识别指标的大小来定量评价储层的 流体性质，能够实现对准致密储层的快速评价，并且能够提高准致密储层评价的可靠性和 精度。此外，还可以通过构建的多种指标来构建多指标识别图版，并基于所述储层流体识别 指标的大小和所述多指标识别图版，来对所述储层的流体性质进行综合定量评价，这样，可 以排除准致密储层评价中的多解性，进一步提高准致密储层评价的可靠性和精度。通过实 例验证，本发明提供的储层评价方法所得到的评价准确度能够达到90%~95%。由此，本 发明提供的准致密储层评价方法能够满足准致密储层测井评价可靠性和精度的要求，有利 于储层评价水平的提高和发展，保证该类储层的有效开发，并且能够为后续地质、地球物理 综合研究奠定坚实基础。
[0111] 以上结合附图详细描述了本发明的优选实施方式，但是，本发明并不限于上述实 施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简 单变型，这些简单变型均属于本发明的保护范围。
[0112] 另外需要说明的是，在上述【具体实施方式】中所描述的各个具体技术特征，在不矛 盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明对各种可 能的组合方式不再另行说明。
[0113] 此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本 发明的思想，其同样应当视为本发明所公开的内容。
【主权项】
1. 一种准致密储层快速评价方法，其特征在于，该方法包括： 根据岩屯、资料和地质资料，得出不同储层深度下的孔隙度、含水饱和度、冲洗带饱和度 W及厚度； 根据每个储层深度下的含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中 可动控体积与水体积的比值； 根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度W及岩石孔隙中可动 控体积与水体积的比值，确定该储层深度下的储层流体识别指标的大小；W及 根据每个储层深度下的储层流体识别指标的大小，识别该储层深度下的流体性质，W 对所述储层的流体性质进行评价。2. 根据权利要求1所述的方法，其特征在于， 通过W下方式确定所述岩石孔隙中可动控体积与水体积的比值：通过W下方式确定所述储层流体识别指标大小：兵bVVKi巧/J、采1 'n甩/岩休议1'的石CfLl银T Kl A刀肛'伴中只3/]\'伴中只的&1且； SWi表示第i个储层深度下的含水饱和度； SXOi表示第i个储层深度下的冲洗带饱和度； 化QIi表示第i个储层深度下的储层流体识别指标； PHIEi表示第i个储层深度下的孔隙度；W及 Hi表示第i个储层深度下的厚度。3. 根据权利要求2所述的方法，其特征在于，根据每个储层深度下的储层流体识别指 标的大小，识别该储层深度下的流体性质的步骤包括： 在化QIi= 0. 00的情况下，识别第i个储层深度下的流体性质为干层； 在化QIi= 0. 50的情况下，识别第i个储层深度下的流体性质为水层； 在化QIi= 1. 00的情况下，识别第i个储层深度下的流体性质为含水含气层； 在化QIi= 1. 50的情况下，识别第i个储层深度下的流体性质为含气层； 在化QIi= 2. 00的情况下，识别第i个储层深度下的流体性质为气水同层； 在化QIi= 2. 50的情况下，识别第i个储层深度下的流体性质为含水气层； 在化QIi= 3. 00的情况下，识别第i个储层深度下的流体性质为差气层；W及 在化QIi= 3. 50的情况下，识别第i个储层深度下的流体性质为气层。4. 一种准致密储层多指标评价方法，其特征在于，该方法包括： 根据权利要求1-3中任一权利要求所述的准致密储层快速评价方法，确定每个储层深 度下的储层流体识别指标的大小； 根据岩屯、资料和地质资料，得出每个储层深度下的泥质含量和渗透率； 对测井资料进行环境校正，得出所述储层的测井曲线，其中，该测井曲线包括由每个储 层深度下的深电阻率测井值形成的深电阻率曲线； 根据每个储层深度下的孔隙度、含水饱和度和冲洗带饱和度，确定该储层深度下的岩 石孔隙中可动控的体积、W及可动控体积与岩石孔隙体积的比值； 根据每个储层深度下的孔隙度和深电阻率测井值，确定该储层深度下的储层流体性质 指标的大小； 根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度、泥质含量、渗透率、 岩石孔隙中可动控的体积、可动控体积与岩石孔隙体积的比值、岩石孔隙中可动控体积与 水体积的比值、W及储层流体性质指标，构建该储层深度下的多指标识别图版；W及 根据每个储层深度下的储层流体识别指标的大小、W及每个储层深度下的多指标识别 图版，共同识别该储层深度下的流体性质，W对所述储层的流体性质进行评价。5. 根据权利要求4所述的准致密储层多指标评价方法，其特征在于，通过W下方式确 定所述岩石孔隙中可动控的体积、W及所述可动控体积与岩石孔隙体积的比值： BVWIi= PHIE i*Sx〇i-PHIEi*SWi其中，BVWIi表示第i个储层深度下的岩石孔隙中可动控的体积； SGMi表示第i个储层深度下的可动控体积与岩石孔隙体积的比值； PHIEi表示第i个储层深度下的孔隙度； SWi表示第i个储层深度下的含水饱和度；W及 SXOi表示第i个储层深度下的冲洗带饱和度。6. 根据权利要求4所述的准致密储层多指标评价方法，其特征在于，根据每个储层深 度下的孔隙度和深电阻率测井值，确定该储层深度下的储层流体性质指标的大小的步骤包 括： 构建所述储层的孔隙度-电阻率交会图版； 将每个储层深度下的孔隙度和深电阻率测井值映射到所述孔隙度-电阻率交会图版 中，形成每个储层深度下的孔隙度-电阻率坐标点；W及 根据每个储层深度下的孔隙度-电阻率坐标点在所述孔隙度-电阻率交会图版中的位 置，确定该储层深度下的储层流体性质指标的大小。7. 根据权利要求4所述的准致密储层多指标评价方法，其特征在于，所述测井曲线还 包括由每个储层深度下的自然伽玛测井值形成的自然伽玛曲线和由每个储层深度下的声 波测井值形成的声波曲线；W及所述准致密储层多指标评价方法还包括： 对所述自然伽玛曲线和所述声波曲线进行归一化； 判断所述测井曲线中是否包括由每个储层深度下的中子测井值形成的中子曲线和由 每个储层深度下的密度测井值形成的密度曲线； 在判定所述测井曲线中未包括所述中子曲线的情况下，利用归一化的自然伽玛曲线和 声波曲线，重构归一化的中子曲线； 在判定所述测井曲线中未包括所述密度曲线的情况下，利用归一化的自然伽玛曲线和 声波曲线，重构归一化的密度曲线； 在判定所述测井曲线中包括所述中子曲线的情况下，判断该中子曲线是否合格，并在 判定所述中子曲线不合格的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化 的中子曲线，W及在判定所述中子曲线合格的情况下，对所述中子曲线进行归一化； 在判定所述测井曲线中包括所述密度曲线的情况下，判断该密度曲线是否合格，并在 判定所述密度曲线不合格的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化 的密度曲线，化及在判定所述密度曲线合格的情况下，对所述密度曲线进行归一化。8. 根据权利要求7所述的准致密储层多指标评价方法，其特征在于，通过W下方式重 构归一化的中子曲线和归一化的密度曲线： NCNLfi= 0. 5707*NAC 1+0. 4114*NGRi-0.0 Ol NDEN"= 0. 1092-0. 6922*NAC 1+0. 1254*NGRi 其中，NACi表示归一化的声波曲线中第i个储层深度下的声波测井值； NGI^i表示归一化的自然伽玛曲线中第i个储层深度下的自然伽玛测井值； NCNLfi表示重构的归一化的中子曲线中第i个储层深度下的中子测井值；W及 NDENfi表示重构的归一化的密度曲线中第i个储层深度下的密度测井值。9. 根据权利要求7或8所述的准致密储层多指标评价方法，其特征在于，该准致密储层 多指标评价方法还包括： 在重构归一化的中子曲线之后，利用所重构的归一化的中子曲线重构中子曲线；和/ 或 在重构归一化的密度曲线之后，利用所重构的归一化的密度曲线重构密度曲线。10. 根据权利要求9所述的准致密储层多指标评价方法，其特征在于，通过W下方式重 构所述中子曲线和所述密度曲线： CNLfi= 0. 1566+0. 0639蝴CNL。 DENfi= 2. 4929+0. 1305*NDEN。 其中，CNLfi表示重构的中子曲线中第i个储层深度下的中子测井值； DENfi表示重构的密度曲线中第i个储层深度下的密度测井值； NCNLfi表示重构的归一化的中子曲线中第i个储层深度下的中子测井值；W及 NDENfi表示重构的归一化的密度曲线中第i个储层深度下的密度测井值。11. 根据权利要求7或8所述的准致密储层多指标评价方法，其特征在于，所述测井曲 线还包括由每个储层深度下的浅电阻率测井值形成的浅电阻率曲线；W及所述准致密储层 多指标评价方法还包括： 对所述深电阻率曲线和所述浅电阻率曲线进行归一化； 根据归一化的自然伽玛曲线、归一化的声波曲线、归一化的深电阻率曲线、W及每个储 层深度下的泥质含量，确定该储层深度下的高角度缝指标的大小； 根据归一化的深电阻率曲线、归一化的浅电阻率曲线、未经归一化的深电阻率曲线、未 经归一化的浅电阻率曲线、W及每个储层深度下的泥质含量，确定该储层深度下的低角度 缝指标的大小；W及 根据每个储层深度下的高角度缝指标的大小和低角度缝指标的大小，对该储层深度下 的裂缝发育情况进行评价，W对所述储层的裂缝发育情况进行评价。12. 根据权利要求11所述的准致密储层多指标评价方法，其特征在于，通过W下方式 确定所述高角度缝指标的大小：通过W下方式确定所述低角度缝指标的大小：其中，VFRIi表示第i个储层深度下的高角度缝指标； V化1表示第i个储层深度下的泥质含量； DMaXi表示第i个储层深度下，归一化的自然伽玛曲线的一阶导数、归一化的声波曲线 的一阶导数、归一化的深电阻率曲线的一阶导数中的最大值； DMini表示第i个储层深度下，归一化的自然伽玛曲线的一阶导数、归一化的声波曲线 的一阶导数、W及归一化的深电阻率曲线的一阶导数中导数间差值的最小值； HFRIi表示第i个储层深度下的低角度缝指标； NRDi表示归一化的深电阻率曲线中第i个储层深度下的深电阻率测井值； NRSi表示归一化的浅电阻率曲线中第i个储层深度下的浅电阻率测井值； A化RD康示未经归一化的深电阻率曲线中所有深电阻率测井值的对数均值；化及 A化RSi表示未经归一化的浅电阻率曲线中所有浅电阻率测井值的对数均值。13. 根据权利要求7或8所述的准致密储层多指标评价方法，其特征在于，该准致密储 层多指标评价方法还包括： 根据归一化的声波曲线、归一化的中子曲线和每个储层深度下的泥质含量，确定该储 层深度下的天然气挖掘效应指标的大小；W及 根据每个储层深度下的天然气挖掘效应指标的大小，对该储层深度下的流体性质进行 定性识别，W对所述储层的流体性质进行评价。14. 根据权利要求13所述的准致密储层多指标评价方法，其特征在于，通过W下方式 确定所述天然气挖掘效应指标的大小：其中，CAVIi表示第i个储层深度下的天然气挖掘效应指标； V化1表示第i个储层深度下的泥质含量； NACi表示归一化的声波曲线中第i个储层深度下的声波测井值；W及 NCNLi表示归一化的中子曲线中第i个储层深度下的中子测井值。
【文档编号】E21B49/00GK105986813SQ201510079776
【公开日】2016年10月5日
【申请日】2015年2月13日
【发明人】胡华伟, 程洪亮, 李定军, 周维娜, 颜晓, 杨杰, 黎静容, 李涛, 陈俊, 李强, 肖维德, 姜镭
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司西南油气分公司