CN105986813A - 准致密储层快速评价方法及准致密储层多指标评价方法 - Google Patents

Abstract

本发明公开了一种准致密储层快速评价方法及准致密储层多指标评价方法。所述准致密储层快速评价方法包括：根据岩心资料和地质资料，得出不同储层深度下的孔隙度、含水饱和度、冲洗带饱和度以及厚度；根据每个储层深度下的含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中可动烃体积与水体积的比值；根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度以及岩石孔隙中可动烃体积与水体积的比值，确定该储层深度下的储层流体识别指标的大小；以及根据每个储层深度下的储层流体识别指标的大小，识别该储层深度下的流体性质，以对所述储层的流体性质进行评价。由此，可以提高储层评价的快速性、准确性和可靠性。

Description

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

技术领域

本发明涉及储层评价领域，具体地，涉及一种准致密储层快速评价方法及准致密储层多指标评价方法。

背景技术

准致密碎屑岩储层具有低孔渗、高束缚水、连通性差、非均质性强的特征，储层评价水平是有效开发该储层的关键因素。测井评价作为三种核心储层评价手段之一，其定性评价的可靠性和定量评价的精度，是决定储层评价水平的基础和关键。但准致密储层测井评价一直是储层测井评价面临的技术难点。由于准致密储层超低的储集空间，孔隙和孔隙流体对测井响应贡献小，储层非均性强，测井响应差异大，因此，传统、常规的储层测井评价方法、评价指标已经不能满足准致密储层测井评价可靠性和精度要求，制约了储层评价水平的提高和发展，并且制约该类储层的有效开发。

发明内容

本发明的目的是提供一种准致密储层快速评价方法及准致密储层多指标评价方法，所述评价方法能够基于新构建的评价指标来对储层进行评价，从而提高储层评价的快速性、准确性和可靠性。

为了实现上述目的，本发明提供一种准致密储层快速评价方法，该方法包括：根据岩心资料和地质资料，得出不同储层深度下的孔隙度、含水饱和度、冲洗带饱和度以及厚度；根据每个储层深度下的含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中可动烃体积与水体积的比值；根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度以及岩石孔隙中可动烃体积与水体积的比值，确定该储层深度下的储层流体识别指标的大小；以及根据每个储层深度下的储层流体识别指标的大小，识别该储层深度下的流体性质，以对所述储层的流体性质进行评价。

优选地，通过以下方式确定所述岩石孔隙中可动烃体积与水体积的比值：

${\mathrm{GWR}}_i=\frac{{\mathrm{Sxo}}_i-{\mathrm{Sw}}_i}{{\mathrm{Sw}}_i}$

通过以下方式确定所述储层流体识别指标大小：

${\mathrm{FLQI}}_i=\left\{\begin{array}{cc}0.00& {\mathrm{PHIE}}_i*H_i\≤0.22\&\&{\mathrm{GWR}}_i>0.40\\ 0.50& {\mathrm{GWR}}_i\≤0.30\&\&{\mathrm{Sw}}_i>0.70\\ 1.00& {\mathrm{PHIE}}_i*H_i>0.22\&\&{\mathrm{Sw}}_i>0.70\&\&{\mathrm{PHIE}}_i*{\mathrm{Sxo}}_i*H_i>0.25\\ 1.50& {\mathrm{PHIE}}_i*H_i>0.22\&\&{\mathrm{Sw}}_i\≤0.70\&\&{\mathrm{PHIE}}_i*{\mathrm{Sxo}}_i*H_i>0.25\\ 2.00& {\mathrm{Sw}}_i>0.53\&\&{\mathrm{Sw}}_i\≤0.70\\ 2.50& {\mathrm{PHIE}}_i*H_i>0.80\&\&{\mathrm{Sw}}_i>0.70\&\&{\mathrm{PHIE}}_i*{\mathrm{Sxo}}_i*H_i>0.40\\ 3.00& {\mathrm{PHIE}}_i*H_i>0.80\&\&{\mathrm{Sw}}_i\≤0.70\&\&{\mathrm{PHIE}}_i*{\mathrm{Sxo}}_i*H_i>0.40\\ 3.50& {\mathrm{PHIE}}_i*H_i>0.80\&\&{\mathrm{Sw}}_i\≤0.53\end{array}\right.$

其中，GWR_i表示第i个储层深度下的岩石孔隙中可动烃体积与水体积的比值；Sw_i表示第i个储层深度下的含水饱和度；Sxo_i表示第i个储层深度下的冲洗带饱和度；FLQI_i表示第i个储层深度下的储层流体识别指标；PHIE_i表示第i个储层深度下的孔隙度；以及H_i表示第i个储层深度下的厚度。

优选地，根据每个储层深度下的储层流体识别指标的大小，识别该储层深度下的流体性质的步骤包括：在FLQI_i＝0.00的情况下，识别第i个储层深度下的流体性质为干层；在FLQI_i＝0.50的情况下，识别第i个储层深度下的流体性质为水层；在FLQI_i＝1.00的情况下，识别第i个储层深度下的流体性质为含水含气层；在FLQI_i＝1.50的情况下，识别第i个储层深度下的流体性质为含气层；在FLQI_i＝2.00的情况下，识别第i个储层深度下的流体性质为气水同层；在FLQI_i＝2.50的情况下，识别第i个储层深度下的流体性质为含水气层；在FLQI_i＝3.00的情况下，识别第i个储层深度下的流体性质为差气层；以及在FLQI_i＝3.50的情况下，识别第i个储层深度下的流体性质为气层。

本发明还提供一种准致密储层多指标评价方法，该方法包括：根据本发明提供的所述准致密储层快速评价方法，确定每个储层深度下的储层流体识别指标的大小；根据岩心资料和地质资料，得出每个储层深度下的泥质含量和渗透率；对测井资料进行环境校正，得出所述储层的测井曲线，其中，该测井曲线包括由每个储层深度下的深电阻率测井值形成的深电阻率曲线；根据每个储层深度下的孔隙度、含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中可动烃的体积、以及可动烃体积与岩石孔隙体积的比值；根据每个储层深度下的孔隙度和深电阻率测井值，确定该储层深度下的储层流体性质指标的大小；根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度、泥质含量、渗透率、岩石孔隙中可动烃的体积、可动烃体积与岩石孔隙体积的比值、岩石孔隙中可动烃体积与水体积的比值、以及储层流体性质指标，构建该储层深度下的多指标识别图版；以及根据每个储层深度下的储层流体识别指标的大小、以及每个储层深度下的多指标识别图版，共同识别该储层深度下的流体性质，以对所述储层的流体性质进行评价。

优选地，通过以下方式确定所述岩石孔隙中可动烃的体积、以及所述可动烃体积与岩石孔隙体积的比值：

BVWI_i＝PHIE_i*Sxo_i-PHIE_i*Sw_i

${\mathrm{SGM}}_i=\frac{{\mathrm{BVWI}}_i}{{\mathrm{PHIE}}_i}$

其中，BVWI_i表示第i个储层深度下的岩石孔隙中可动烃的体积；SGM_i表示第i个储层深度下的可动烃体积与岩石孔隙体积的比值；PHIE_i表示第i个储层深度下的孔隙度；Sw_i表示第i个储层深度下的含水饱和度；以及Sxo_i表示第i个储层深度下的冲洗带饱和度。

优选地，根据每个储层深度下的孔隙度和深电阻率测井值，确定该储层深度下的储层流体性质指标的大小的步骤包括：构建所述储层的孔隙度-电阻率交会图版；将每个储层深度下的孔隙度和深电阻率测井值映射到所述孔隙度-电阻率交会图版中，形成每个储层深度下的孔隙度-电阻率坐标点；以及根据每个储层深度下的孔隙度-电阻率坐标点在所述孔隙度-电阻率交会图版中的位置，确定该储层深度下的储层流体性质指标的大小。

优选地，所述测井曲线还包括由每个储层深度下的自然伽玛测井值形成的自然伽玛曲线和由每个储层深度下的声波测井值形成的声波曲线；以及所述准致密储层多指标评价方法还包括：对所述自然伽玛曲线和所述声波曲线进行归一化；判断所述测井曲线中是否包括由每个储层深度下的中子测井值形成的中子曲线和由每个储层深度下的密度测井值形成的密度曲线；在判定所述测井曲线中未包括所述中子曲线的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的中子曲线；在判定所述测井曲线中未包括所述密度曲线的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的密度曲线；在判定所述测井曲线中包括所述中子曲线的情况下，判断该中子曲线是否合格，并在判定所述中子曲线不合格的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的中子曲线，以及在判定所述中子曲线合格的情况下，对所述中子曲线进行归一化；在判定所述测井曲线中包括所述密度曲线的情况下，判断该密度曲线是否合格，并在判定所述密度曲线不合格的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的密度曲线，以及在判定所述密度曲线合格的情况下，对所述密度曲线进行归一化。

优选地，通过以下方式重构归一化的中子曲线和归一化的密度曲线：

NCNL_ri＝0.5707*NAC_i+0.4114*NGR_i-0.001

NDEN_ri＝0.1092-0.6922*NAC_i+0.1254*NGR_i

其中，NAC_i表示归一化的声波曲线中第i个储层深度下的声波测井值；NGR_i表示归一化的自然伽玛曲线中第i个储层深度下的自然伽玛测井值；NCNL_ri表示重构的归一化的中子曲线中第i个储层深度下的中子测井值；以及NDEN_ri表示重构的归一化的密度曲线中第i个储层深度下的密度测井值。

优选地，该准致密储层多指标评价方法还包括：在重构归一化的中子曲线之后，利用所重构的归一化的中子曲线重构中子曲线；和/或在重构归一化的密度曲线之后，利用所重构的归一化的密度曲线重构密度曲线。

优选地，通过以下方式重构所述中子曲线和所述密度曲线：

CNL_ri＝0.1566+0.0639*NCNL_ri

DEN_ri＝2.4929+0.1305*NDEN_ri

其中，CNL_ri表示重构的中子曲线中第i个储层深度下的中子测井值；DEN_ri表示重构的密度曲线中第i个储层深度下的密度测井值；NCNL_ri表示重构的归一化的中子曲线中第i个储层深度下的中子测井值；以及NDEN_ri表示重构的归一化的密度曲线中第i个储层深度下的密度测井值。

优选地，所述测井曲线还包括由每个储层深度下的浅电阻率测井值形成的浅电阻率曲线；以及所述准致密储层多指标评价方法还包括：对所述深电阻率曲线和所述浅电阻率曲线进行归一化；根据归一化的自然伽玛曲线、归一化的声波曲线、归一化的深电阻率曲线、以及每个储层深度下的泥质含量，确定该储层深度下的高角度缝指标的大小；根据归一化的深电阻率曲线、归一化的浅电阻率曲线、未经归一化的深电阻率曲线、未经归一化的浅电阻率曲线、以及每个储层深度下的泥质含量，确定该储层深度下的低角度缝指标的大小；以及根据每个储层深度下的高角度缝指标的大小和低角度缝指标的大小，对该储层深度下的裂缝发育情况进行评价，以对所述储层的裂缝发育情况进行评价。

优选地，通过以下方式确定所述高角度缝指标的大小：

通过以下方式确定所述低角度缝指标的大小：

其中，VFRI_i表示第i个储层深度下的高角度缝指标；VCL_i表示第i个储层深度下的泥质含量；DMax_i表示第i个储层深度下，归一化的自然伽玛曲线的一阶导数、归一化的声波曲线的一阶导数、归一化的深电阻率曲线的一阶导数中的最大值；DMin_i表示第i个储层深度下，归一化的自然伽玛曲线的一阶导数、归一化的声波曲线的一阶导数、以及归一化的深电阻率曲线的一阶导数中导数间差值的最小值；HFRI_i表示第i个储层深度下的低角度缝指标；NRD_i表示归一化的深电阻率曲线中第i个储层深度下的深电阻率测井值；NRS_i表示归一化的浅电阻率曲线中第i个储层深度下的浅电阻率测井值；AvLRD_i表示未经归一化的深电阻率曲线中所有深电阻率测井值的对数均值；以及AvLRS_i表示未经归一化的浅电阻率曲线中所有浅电阻率测井值的对数均值。

优选地，该准致密储层多指标评价方法还包括：根据归一化的声波曲线、归一化的中子曲线和每个储层深度下的泥质含量，确定该储层深度下的天然气挖掘效应指标的大小；以及根据每个储层深度下的天然气挖掘效应指标的大小，对该储层深度下的流体性质进行定性识别，以对所述储层的流体性质进行评价。

优选地，通过以下方式确定所述天然气挖掘效应指标的大小：

其中，CAVI_i表示第i个储层深度下的天然气挖掘效应指标；VCL_i表示第i个储层深度下的泥质含量；NAC_i表示归一化的声波曲线中第i个储层深度下的声波测井值；以及NCNL_i表示归一化的中子曲线中第i个储层深度下的中子测井值。

在上述技术方案中，通过构建储层流体识别指标，并基于该储层流体识别指标的大小来定量评价储层的流体性质，能够实现对准致密储层的快速评价，并且能够提高准致密储层评价的可靠性和精度。此外，还可以通过构建的多种指标来构建多指标识别图版，并基于所述储层流体识别指标的大小和所述多指标识别图版，来对所述储层的流体性质进行综合定量评价，这样，可以排除准致密储层评价中的多解性，进一步提高准致密储层评价的可靠性和精度。通过实例验证，本发明提供的储层评价方法所得到的评价准确度能够达到90％～95％。由此，本发明提供的准致密储层评价方法能够满足准致密储层测井评价可靠性和精度的要求，有利于储层评价水平的提高和发展，保证该类储层的有效开发，并且能够为后续地质、地球物理综合研究奠定坚实基础。

本发明的其他特征和优点将在随后的具体实施方式部分予以详细说明。

附图说明

附图是用来提供对本发明的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明，但并不构成对本发明的限制。在附图中：

图1示出了根据本发明的实施方式的准致密储层快速评价方法的流程图；

图2示出了根据本发明的实施方式的准致密储层多指标评价方法的流程图；

图3示出了根据本发明的实施方式的孔隙度-电阻率交会图版的示例；

图4示出了根据本发明的实施方式的多指标识别图版的示例；以及

图5和图6分别示出了利用本发明提供的准致密储层多指标评价方法对四川盆地某井区进行实际评价所得的结果。

具体实施方式

以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。

图1示出了根据本发明的实施方式的准致密储层快速评价方法的流程图。如图1所示，该方法可以包括：步骤S11，根据岩心资料和地质资料，得出不同储层深度下的孔隙度、含水饱和度、冲洗带饱和度以及厚度；步骤S12，根据每个储层深度下的含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中可动烃体积与水体积的比值；步骤S13，根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度以及岩石孔隙中可动烃体积与水体积的比值，确定该储层深度下的储层流体识别指标的大小；以及步骤S14，根据每个储层深度下的储层流体识别指标的大小，识别该储层深度下的流体性质，以对所述储层的流体性质进行评价。

具体地，首先在步骤S11，根据岩心资料和地质资料，得出不同储层深度下的孔隙度PHIE、含水饱和度Sw、冲洗带饱和度Sxo以及厚度H。该步骤的实现是本领域技术人员公知的，对此，本发明在此不进行详细描述。

在得出不同储层深度下的孔隙度PHIE、含水饱和度Sw、冲洗带饱和度Sxo以及厚度H之后，可以根据每个储层深度下的含水饱和度Sw和冲洗带饱和度Sxo，确定出该储层深度下的岩石孔隙中可动烃体积与水体积的比值GWR。该比值是本发明新构建的一项指标，和下文描述的岩石孔隙中可动烃的体积BVWI和可动烃体积与岩石孔隙体积的比值SGM这两项指标共同构成“相对体积指标”，该组指标用于后续的储层评价。

可以通过以下等式(1)来确定所述岩石孔隙中可动烃体积与水体积的比值：

${\mathrm{GWR}}_i=\frac{{\mathrm{Sxo}}_i-{\mathrm{Sw}}_i}{{\mathrm{Sw}}_i}---\left(1\right)$

其中，GWR_i表示第i个储层深度下的岩石孔隙中可动烃体积与水体积的比值；Sw_i表示第i个储层深度下的含水饱和度；Sxo_i表示第i个储层深度下的冲洗带饱和度；1≤i≤N，N表示储层的采样点总数，其中，每个采样点对应一个储层深度，第i个采样点也等同表示第i个储层深度。

之后，在确定出每个储层深度下的岩石孔隙中可动烃体积与水体积的比值之后，就可以结合该比值、以及每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度和厚度，来共同确定出每个储层深度下的储层流体识别指标的大小。具体地，可通过以下等式(2)来确定所述储层流体识别指标的大小：

${\mathrm{FLQI}}_i=\left\{\begin{array}{cc}0.00& {\mathrm{PHIE}}_i*H_i\≤0.22\&\&{\mathrm{GWR}}_i>0.40\\ 0.50& {\mathrm{GWR}}_i\≤0.30\&\&{\mathrm{Sw}}_i>0.70\\ 1.00& {\mathrm{PHIE}}_i*H_i>0.22\&\&{\mathrm{Sw}}_i>0.70\&\&{\mathrm{PHIE}}_i*{\mathrm{Sxo}}_i*H_i>0.25\\ 1.50& {\mathrm{PHIE}}_i*H_i>0.22\&\&{\mathrm{Sw}}_i\≤0.70\&\&{\mathrm{PHIE}}_i*{\mathrm{Sxo}}_i*H_i>0.25\\ 2.00& {\mathrm{Sw}}_i>0.53\&\&{\mathrm{Sw}}_i\≤0.70\\ 2.50& {\mathrm{PHIE}}_i*H_i>0.80\&\&{\mathrm{Sw}}_i>0.70\&\&{\mathrm{PHIE}}_i*{\mathrm{Sxo}}_i*H_i>0.40\\ 3.00& {\mathrm{PHIE}}_i*H_i>0.80\&\&{\mathrm{Sw}}_i\≤0.70\&\&{\mathrm{PHIE}}_i*{\mathrm{Sxo}}_i*H_i>0.40\\ 3.50& {\mathrm{PHIE}}_i*H_i>0.80\&\&{\mathrm{Sw}}_i\≤0.53\end{array}\right.---\left(2\right)$

其中，FLQI_i表示第i个储层深度下的储层流体识别指标；PHIE_i表示第i个储层深度下的孔隙度；以及H_i表示第i个储层深度下的厚度。

储层流体识别指标FLQI同样是本发明新构建的一项指标，其用于对储层进行快速、初步自动评价，特别是针对储层的流体性质进行快速、初步自动评价，如下面进一步描述的。

在FLQI_i＝0.00的情况下，识别第i个储层深度下的流体性质为干层；在FLQI_i＝0.50的情况下，识别第i个储层深度下的流体性质为水层；在FLQI_i＝1.00的情况下，识别第i个储层深度下的流体性质为含水含气层；在FLQI_i＝1.50的情况下，识别第i个储层深度下的流体性质为含气层；在FLQI_i＝2.00的情况下，识别第i个储层深度下的流体性质为气水同层；在FLQI_i＝2.50的情况下，识别第i个储层深度下的流体性质为含水气层；在FLQI_i＝3.00的情况下，识别第i个储层深度下的流体性质为差气层；以及在FLQI_i＝3.50的情况下，识别第i个储层深度下的流体性质为气层。

通过上述方式，可以快速识别出每个储层深度下的流体性质，通过将这些识别结果进行组合，可以对所述储层的流体性质做出整体评价，分析出为干层、水层、含水含气层、含气层、气水同层、含水气层、差气层和气层的各储层段，以为日后的储层开发、地质研究奠定基础。

通过本发明构建的储层流体识别指标来快速评价储层的流体性质，可以实现对准致密储层的流体性质的自动、快速地定量评价，具有较高的准确性和可靠性，填补了现有技术中缺乏对准致密储层的有效评价方法的空白。该指标可以用于对储层的流体性质的预评价，有利于快速、准确评价储层的流体性质。此外，该指标还可以与其他定量评价指标一起，共同对储层的流体性质进行定量评价，以实现基于多指标的准致密储层评价方法，进一步提高评价的可靠性和准确性。下文将具体描述这一准致密储层多指标评价方法。

图2示出了根据本发明的实施方式的准致密储层多指标评价方法的流程图。如图2所示，该方法可以包括：步骤S21，根据本发明提供的上述准致密储层快速评价方法，确定每个储层深度下的储层流体识别指标的大小；步骤S22，根据岩心资料和地质资料，得出每个储层深度下的泥质含量和渗透率；步骤S23，对测井资料进行环境校正，得出所述储层的测井曲线，其中，该测井曲线包括由每个储层深度下的深电阻率测井值形成的深电阻率曲线；步骤S24，根据每个储层深度下的孔隙度、含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中可动烃的体积、以及可动烃体积与岩石孔隙体积的比值；步骤S25，根据每个储层深度下的孔隙度和深电阻率测井值，确定该储层深度下的储层流体性质指标的大小；步骤S26，根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度、泥质含量、渗透率、岩石孔隙中可动烃的体积、可动烃体积与岩石孔隙体积的比值、岩石孔隙中可动烃体积与水体积的比值、以及储层流体性质指标的大小，构建该储层深度下的多指标识别图版；以及步骤S27，根据每个储层深度下的储层流体识别指标的大小、以及每个储层深度下的多指标识别图版，共同识别该储层深度下的流体性质，以对所述储层的流体性质进行评价。

具体地，首先，根据上面结合图1所描述的准致密储层快速评价方法，确定出每个储层深度下的储层流体识别指标的大小，即FLQI_i。此外，除了根据岩心资料和地质资料，得出不同储层深度下的孔隙度、含水饱和度、冲洗带饱和度以及厚度之外，还要根据该岩心资料和地质资料，得出每个储层深度下的泥质含量和渗透率。该步骤的实现是本领域技术人员公知的，对此，本发明在此不进行详细描述。

接下来，在步骤S23，对测井资料进行环境校正，得出所述储层的测井曲线，其中，该测井曲线可以包括由每个储层深度下的深电阻率测井值形成的深电阻率曲线。对测井资料进行环境校正，是为了排除井径、泥浆等非地层因素对测井资料的影响。如何对测井资料进行环境校正的具体方式是本领域的技术人员公知的，对此，本发明在此不做具体描述。

接下来，在步骤S24，根据每个储层深度下的孔隙度、含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中可动烃的体积、以及可动烃体积与岩石孔隙体积的比值，即上文提到的“相对体积指标”中的另外两项指标BVWI和SGM。可以通过以下等式(3)和(4)来分别确定这两项指标：

BVWI_i＝PHIE_i*Sxo_i-PHIE_i*Sw_i (3)

${\mathrm{SGM}}_i=\frac{{\mathrm{BVWI}}_i}{{\mathrm{PHIE}}_i}---\left(4\right)$

其中，BVWI_i表示第i个储层深度下的岩石孔隙中可动烃的体积；SGM_i表示第i个储层深度下的可动烃体积与岩石孔隙体积的比值；PHIE_i表示第i个储层深度下的孔隙度；Sw_i表示第i个储层深度下的含水饱和度；以及Sxo_i表示第i个储层深度下的冲洗带饱和度。如此，上述等式(1)还可以被扩展成以下等式(5)：

${\mathrm{GWR}}_i=\frac{{\mathrm{BVWI}}_i}{{\mathrm{PHIE}}_i*{\mathrm{Sw}}_i}=\frac{{\mathrm{Sxo}}_i-{\mathrm{Sw}}_i}{{\mathrm{Sw}}_i}---\left(5\right)$

所述“相对体积指标”中的每个指标作为后面构建多指标识别图版中的三项定量指标。通过新增这三项定量指标来进行储层评价，可以克服准致密储层评价定量参数误差大的难题。

虽然使用了S22、S23和S24这样的表述，但是并不意味着这三个步骤是有先后次序的，这三个步骤可以是同时并行的。

接下来，在步骤S25，根据每个储层深度下的孔隙度和深电阻率测井值，确定该储层深度下的储层流体性质指标的大小。具体地，该步骤S25可以包括如下几个步骤：构建所述储层的孔隙度-电阻率交会图版；将每个储层深度下的孔隙度和深电阻率测井值映射到所述孔隙度-电阻率交会图版中，形成每个储层深度下的孔隙度-电阻率坐标点。图3示出了所述孔隙度-电阻率交会图版的示例。在该交会图版中，深电阻率测井值RD作为纵坐标，孔隙度PHIE作为横坐标。每个储层深度下的孔隙度和深电阻率测井值被映射到该孔隙度-电阻率交会图版中，这样，就形成每个储层深度下的孔隙度-电阻率坐标点。如上所述，总共有N个储层深度，因此，可以形成N个孔隙度-电阻率坐标点。

可以根据先验数据或经验，在该交会图版中预先划分出若干区域。在图3所示的交会图版中，通过5条线将交会图版划分为6个非线性区域，分别为区域I、区域II、区域III、区域IV、区域V以及区域VI。每个区域代表一个储层流体性质指标IGAS的值。这样，就可以根据每个储层深度下的孔隙度-电阻率坐标点在所述孔隙度-电阻率交会图版中的位置(即，落到哪个区域)，来确定该储层深度下的储层流体性质指标的大小。

例如，针对第i个储层深度，如果其孔隙度-电阻率坐标点落在区域I，则确定该第i个储层深度下，IGAS_i＝4.5；如果其孔隙度-电阻率坐标点落在区域II，则确定该第i个储层深度下，IGAS_i＝3.5；如果其孔隙度-电阻率坐标点落在区域III，则确定该第i个储层深度下，IGAS_i＝2.5；如果其孔隙度-电阻率坐标点落在区域IV，则确定该第i个储层深度下，IGAS_i＝1.5；以及如果其孔隙度-电阻率坐标点落在区域V或区域VI，则确定该第i个储层深度下，IGAS_i＝0.5。通过这种方式，可以确定出每个储层深度下的储层流体性质指标的大小。该储层流体性质指标也作为后面构建多指标识别图版中的一项定量指标。

接下来，在步骤S26中，构建每个储层深度下的多指标识别图版。具体地，在本发明中，通过10种指标来构建多指标识别图版。例如，如图4所示，这10种指标分别是：含水饱和度(Sw)；孔隙度(PHIE)；孔隙度与厚度的乘积(PHIE*H)；孔隙度、冲洗带饱和度与厚度的乘积(PHIE*Swo*H)；泥质含量(VCL)；渗透率(PERM)；岩石孔隙中可动烃的体积(BVMI)；可动烃体积与岩石孔隙体积的比值(SGM)；岩石孔隙中可动烃体积与水体积的比值(GWR)；以及储层流体性质指标(IGAS)。

在本发明中，储层的流体性质可以分为八种，分别是：干层、水层、含水含气层、含气层、气水同层、含水气层、差气层和气层。可以根据先验数据构建出每种流体性质对应的标准多指标识别图版。之后，将新构建的一储层深度下的多指标识别图版与这八种标准多指标识别图版进行匹配，确定与之最接近的标准多指标识别图版，并最终识别出该储层深度下的流体性质。

在本发明中，可以结合使用图1中描述的方法确定出的每个储层深度下的储层流体识别指标的大小、以及每个储层深度下的多指标识别图版，来共同识别该储层深度下的流体性质，以对所述储层的流体性质做出综合评价。

通过上述准致密储层多指标评价方法，可以融合多种定量指标来对准致密储层进行综合评价，相比于仅采用一种指标进行评价，可以有效消除准致密储层评价的多解性，从而可以进一步提高准致密储层评价结果的可靠性和准确性。

此外，在上面通过步骤S23获得的测井曲线中，不仅可以包括所述深电阻率曲线，还可以包括由每个储层深度下的自然伽玛测井值形成的自然伽玛曲线和由每个储层深度下的声波测井值形成的声波曲线。并且，所述准致密储层多指标评价方法还可以包括：对所述自然伽玛曲线和所述声波曲线进行归一化；判断所述测井曲线中是否包括由每个储层深度下的中子测井值形成的中子曲线和由每个储层深度下的密度测井值形成的密度曲线；在判定所述测井曲线中未包括所述中子曲线的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的中子曲线；在判定所述测井曲线中未包括所述密度曲线的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的密度曲线；在判定所述测井曲线中包括所述中子曲线的情况下，判断该中子曲线是否合格，并在判定所述中子曲线不合格的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的中子曲线，以及在判定所述中子曲线合格的情况下，对所述中子曲线进行归一化；在判定所述测井曲线中包括所述密度曲线的情况下，判断该密度曲线是否合格，并在判定所述密度曲线不合格的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的密度曲线，以及在判定所述密度曲线合格的情况下，对所述密度曲线进行归一化。

通常情况下，测井曲线中要包括所述中子曲线和密度曲线。中子曲线和密度曲线可用于辅助泥质含量、孔隙度、渗透率等参数的计算。然而，密度曲线和中子曲线数据的采集对井况条件要求较高，有时可能采集不到密度曲线和中子曲线，或者采集到的密度曲线和中子曲线质量差，会影响后期的储层评价。这就导致测井曲线中可能缺失中子曲线和/或密度曲线，或者虽然包括中子曲线和/或密度曲线，但是所包括的中子曲线和/或密度曲线不合格。在这种情况下，为了不影响其他参数的计算及后续的储层评价，可以按照上述方法，对缺失的或不合格的中子曲线和/或密度曲线进行重构。

具体地，可以通过以下等式(6)来重构归一化的中子曲线，以及通过以下等式(7)来重构归一化的密度曲线：

NCNL_ri＝0.5707*NAC_i+0.4114*NGR_i-0.001 (6)

NDEN_ri＝0.1092-0.6922*NAC_i+0.1254*NGR_i (7)

其中，NAC_i表示归一化的声波曲线中第i个储层深度下的声波测井值；NGR_i表示归一化的自然伽玛曲线中第i个储层深度下的自然伽玛测井值；NCNL_ri表示重构的归一化的中子曲线中第i个储层深度下的中子测井值；以及NDEN_ri表示重构的归一化的密度曲线中第i个储层深度下的密度测井值。

由于采用归一化的自然伽玛曲线和归一化的声波曲线来进行重构，因此，通过上面等式(6)和(7)重构出的中子曲线和密度曲线也是归一化的中子曲线和归一化的密度曲线。

如上所述，中子曲线和密度曲线可用于其他定量指标的计算。因此，在重构出归一化的中子曲线和/或归一化的密度曲线之后，还可以利用所重构的归一化的中子曲线重构中子曲线；和/或利用所重构的归一化的密度曲线重构密度曲线。

具体地，可通过以下等式(8)重构所述中子曲线，以及通过等式(9)重构所述密度曲线：

CNL_ri＝0.1566+0.0639*NCNL_ri (8)

DEN_ri＝2.4929+0.1305*NDEN_ri (9)

其中，CNL_ri表示重构的中子曲线中第i个储层深度下的中子测井值；以及DEN_ri表示重构的密度曲线中第i个储层深度下的密度测井值。

这样，就可以在无法采集中子曲线和密度曲线、或者采集到的中子曲线和密度曲线质量不合格的情况下，重构出中子曲线和密度曲线。由此，可以满足大斜度、水平井以及多井综合评价的资料需求。

如何判断中子曲线和密度曲线是否合格的方法是本领域的技术人员公知的，例如，通过确定密度曲线、中子曲线分别与自然伽玛曲线、深电阻率曲线、声波曲线等之间的相关性差，来判断密度曲线、中子曲线是否合格(相关性差过大，判断曲线不合格)。因此，对于曲线是否合格的判断方法本发明不做具体描述。

另外，在所述中子曲线和密度曲线未缺失、并且合格的情况下，对所述中子曲线和密度曲线进行归一化。对自然伽玛曲线、声波曲线、中子曲线和密度曲线进行归一化的方法有很多，并且均是本领域的技术人员公知的。在本发明的一个示例实施方式中，可以采用以下等式(10)来进行归一化：

$N{\log }_i=\frac{{\log }_i-\mathrm{Av}\log }{\mathrm{Std}\log }---\left(10\right)$

其中，log_i表示曲线上第i个测井值；Nlog_i表示该第i个测井值的归一化值；Av log表示曲线上所有测井值的平均值；以及Std log表示标准差。

上述示例归一化方法不用于限制本发明，本领域的技术人员也可以采用其他归一化方法来对自然伽玛曲线、声波曲线、中子曲线和密度曲线进行归一化，对此，本发明在此不一一具体阐述。

在本发明的另一实施方式中，所述测井曲线还可以包括由每个储层深度下的浅电阻率测井值形成的浅电阻率曲线；以及所述准致密储层多指标评价方法还可以包括：对所述深电阻率曲线和所述浅电阻率曲线进行归一化；根据归一化的自然伽玛曲线、归一化的声波曲线、归一化的深电阻率曲线、以及每个储层深度下的泥质含量，确定该储层深度下的高角度缝指标的大小；根据归一化的深电阻率曲线、归一化的浅电阻率曲线、未经归一化的深电阻率曲线、未经归一化的浅电阻率曲线、以及每个储层深度下的泥质含量，确定该储层深度下的低角度缝指标的大小；以及根据每个储层深度下的高角度缝指标的大小和低角度缝指标的大小，对该储层深度下的裂缝发育情况进行评价，以对所述储层的裂缝发育情况进行评价。

在对准致密储层进行评价时，储层的裂缝发育情况也是一项评价指标。在本发明中，通过构造高角度缝指标和低角度缝指标，可以快速且准确地定性分析出储层的裂缝发育情况。

具体地，可通过以下等式(11)确定所述高角度缝指标的大小：

其中，VFRI_i表示第i个储层深度下的高角度缝指标；VCL_i表示第i个储层深度下的泥质含量；DMax_i表示第i个储层深度下，归一化的自然伽玛曲线的一阶导数、归一化的声波曲线的一阶导数、归一化的深电阻率曲线的一阶导数中的最大值；DMin_i表示第i个储层深度下，归一化的自然伽玛曲线的一阶导数、归一化的声波曲线的一阶导数、以及归一化的深电阻率曲线的一阶导数中导数间差值的最小值。

此外，可通过以下等式(12)确定所述低角度缝指标的大小：

其中，HFRI_i表示第i个储层深度下的低角度缝指标；NRD_i表示归一化的深电阻率曲线中第i个储层深度下的深电阻率测井值；NRS_i表示归一化的浅电阻率曲线中第i个储层深度下的浅电阻率测井值；AvLRD_i表示未经归一化的深电阻率曲线中所有深电阻率测井值的对数均值；以及AvLRS_i表示未经归一化的浅电阻率曲线中所有浅电阻率测井值的对数均值。

之后，就可以根据所得到的高角度缝指标的大小和低角度缝指标的大小，对每个储层深度下的裂缝发育情况进行定性分析。例如，在VFRI_i＝1.0和/或HFRI_i＝1.0的情况下，可以定性分析出该第i个储层深度下的裂缝发育；而在VFRI_i＝0.1且HFRI_i＝0.1的情况下，可以定性分析出该第i个储层深度下的裂缝未发育。由此，可以分析出整个准致密储层的裂缝发育情况，以为日后的储层开发、地质研究奠定基础。

在本发明的另一实施方式中，所述准致密储层多指标评价方法还可以包括：根据归一化的声波曲线、归一化的中子曲线和每个储层深度下的泥质含量，确定该储层深度下的天然气挖掘效应指标的大小；以及根据每个储层深度下的天然气挖掘效应指标的大小，对该储层深度下的流体性质进行定性识别，以对所述储层的流体性质进行评价。

具体地，可通过以下等式(13)确定所述天然气挖掘效应指标的大小：

其中，CAVI_i表示第i个储层深度下的天然气挖掘效应指标；VCL_i表示第i个储层深度下的泥质含量；NAC_i表示归一化的声波曲线中第i个储层深度下的声波测井值；以及NCNL_i表示归一化的中子曲线中第i个储层深度下的中子测井值。

如果是经重构得出的归一化的中子曲线，则NCNL_i等同于NCNL_ri。如果不是经重构得出的归一化的中子曲线，那么NCNL_i则代表对测井曲线中包括的原本中子曲线进行归一化后的结果。

区别于上面通过储层流体识别指标FLQI以及多指标识别图版来对储层的流体性质进行定量识别，在此种实施方式中，通过每个储层深度下的天然气挖掘效应指标CAVI的大小，可以实现对每个储层深度下的流体性质的定性识别。在储层评价体系中，不仅需要定量识别，也需要定性识别。通过本发明构造出的天然气挖掘效应指标，可以提高储层流体性质的定性识别的可靠性。例如，在CAVI_i＝1.0的情况下，可以定性识别出该第i个储层深度下含气；而在CAVI_i＝0.1的情况下，可以定性识别出该第i个储层深度下不含气。

图5和图6分别示出了利用本发明提供的准致密储层多指标评价方法对四川盆地某井区进行实际评价所得的结果。经过与实际测试结果进行比较，通过本发明提供的准致密储层快速评价方法及准致密储层多指标评价方法，能够实现对准致密储层的准确评价。

综上所述，在本发明提供的准致密储层快速评价方法及准致密储层多指标评价方法中，通过构建储层流体识别指标，并基于该储层流体识别指标的大小来定量评价储层的流体性质，能够实现对准致密储层的快速评价，并且能够提高准致密储层评价的可靠性和精度。此外，还可以通过构建的多种指标来构建多指标识别图版，并基于所述储层流体识别指标的大小和所述多指标识别图版，来对所述储层的流体性质进行综合定量评价，这样，可以排除准致密储层评价中的多解性，进一步提高准致密储层评价的可靠性和精度。通过实例验证，本发明提供的储层评价方法所得到的评价准确度能够达到90％～95％。由此，本发明提供的准致密储层评价方法能够满足准致密储层测井评价可靠性和精度的要求，有利于储层评价水平的提高和发展，保证该类储层的有效开发，并且能够为后续地质、地球物理综合研究奠定坚实基础。

以上结合附图详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简单变型，这些简单变型均属于本发明的保护范围。

另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明对各种可能的组合方式不再另行说明。

此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明的思想，其同样应当视为本发明所公开的内容。

Claims (14)

1.一种准致密储层快速评价方法，其特征在于，该方法包括：

根据岩心资料和地质资料，得出不同储层深度下的孔隙度、含水饱和度、冲洗带饱和度以及厚度；

根据每个储层深度下的含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中可动烃体积与水体积的比值；

根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度以及岩石孔隙中可动烃体积与水体积的比值，确定该储层深度下的储层流体识别指标的大小；以及

根据每个储层深度下的储层流体识别指标的大小，识别该储层深度下的流体性质，以对所述储层的流体性质进行评价。

2.根据权利要求1所述的方法，其特征在于，

通过以下方式确定所述岩石孔隙中可动烃体积与水体积的比值：

${\mathrm{GWR}}_i=\frac{{\mathrm{Sxo}}_i-{\mathrm{Sw}}_i}{{\mathrm{Sw}}_i}$

通过以下方式确定所述储层流体识别指标大小：

其中，GWR_i表示第i个储层深度下的岩石孔隙中可动烃体积与水体积的比值；

Sw_i表示第i个储层深度下的含水饱和度；

Sxo_i表示第i个储层深度下的冲洗带饱和度；

FLQI_i表示第i个储层深度下的储层流体识别指标；

PHIE_i表示第i个储层深度下的孔隙度；以及

H_i表示第i个储层深度下的厚度。

3.根据权利要求2所述的方法，其特征在于，根据每个储层深度下的储层流体识别指标的大小，识别该储层深度下的流体性质的步骤包括：

在FLQI_i＝0.00的情况下，识别第i个储层深度下的流体性质为干层；

在FLQI_i＝0.50的情况下，识别第i个储层深度下的流体性质为水层；

在FLQI_i＝1.00的情况下，识别第i个储层深度下的流体性质为含水含气层；

在FLQI_i＝1.50的情况下，识别第i个储层深度下的流体性质为含气层；

在FLQI_i＝2.00的情况下，识别第i个储层深度下的流体性质为气水同层；

在FLQI_i＝2.50的情况下，识别第i个储层深度下的流体性质为含水气层；

在FLQI_i＝3.00的情况下，识别第i个储层深度下的流体性质为差气层；以及

在FLQI_i＝3.50的情况下，识别第i个储层深度下的流体性质为气层。

4.一种准致密储层多指标评价方法，其特征在于，该方法包括：

根据权利要求1-3中任一权利要求所述的准致密储层快速评价方法，确定每个储层深度下的储层流体识别指标的大小；

根据岩心资料和地质资料，得出每个储层深度下的泥质含量和渗透率；

对测井资料进行环境校正，得出所述储层的测井曲线，其中，该测井曲线包括由每个储层深度下的深电阻率测井值形成的深电阻率曲线；

根据每个储层深度下的孔隙度、含水饱和度和冲洗带饱和度，确定该储层深度下的岩石孔隙中可动烃的体积、以及可动烃体积与岩石孔隙体积的比值；

根据每个储层深度下的孔隙度和深电阻率测井值，确定该储层深度下的储层流体性质指标的大小；

根据每个储层深度下的孔隙度、含水饱和度、冲洗带饱和度、厚度、泥质含量、渗透率、岩石孔隙中可动烃的体积、可动烃体积与岩石孔隙体积的比值、岩石孔隙中可动烃体积与水体积的比值、以及储层流体性质指标，构建该储层深度下的多指标识别图版；以及

根据每个储层深度下的储层流体识别指标的大小、以及每个储层深度下的多指标识别图版，共同识别该储层深度下的流体性质，以对所述储层的流体性质进行评价。

5.根据权利要求4所述的准致密储层多指标评价方法，其特征在于，通过以下方式确定所述岩石孔隙中可动烃的体积、以及所述可动烃体积与岩石孔隙体积的比值：

BVWI_i＝PHIE_i*Sxo_i-PHIE_i*Sw_i

${\mathrm{SGM}}_i=\frac{{\mathrm{BVWI}}_i}{{\mathrm{PHIE}}_i}$

其中，BVWI_i表示第i个储层深度下的岩石孔隙中可动烃的体积；

SGM_i表示第i个储层深度下的可动烃体积与岩石孔隙体积的比值；

PHIE_i表示第i个储层深度下的孔隙度；

Sw_i表示第i个储层深度下的含水饱和度；以及

Sxo_i表示第i个储层深度下的冲洗带饱和度。

6.根据权利要求4所述的准致密储层多指标评价方法，其特征在于，根据每个储层深度下的孔隙度和深电阻率测井值，确定该储层深度下的储层流体性质指标的大小的步骤包括：

构建所述储层的孔隙度-电阻率交会图版；

将每个储层深度下的孔隙度和深电阻率测井值映射到所述孔隙度-电阻率交会图版中，形成每个储层深度下的孔隙度-电阻率坐标点；以及

根据每个储层深度下的孔隙度-电阻率坐标点在所述孔隙度-电阻率交会图版中的位置，确定该储层深度下的储层流体性质指标的大小。

7.根据权利要求4所述的准致密储层多指标评价方法，其特征在于，所述测井曲线还包括由每个储层深度下的自然伽玛测井值形成的自然伽玛曲线和由每个储层深度下的声波测井值形成的声波曲线；以及所述准致密储层多指标评价方法还包括：

对所述自然伽玛曲线和所述声波曲线进行归一化；

判断所述测井曲线中是否包括由每个储层深度下的中子测井值形成的中子曲线和由每个储层深度下的密度测井值形成的密度曲线；

在判定所述测井曲线中未包括所述中子曲线的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的中子曲线；

在判定所述测井曲线中未包括所述密度曲线的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的密度曲线；

在判定所述测井曲线中包括所述中子曲线的情况下，判断该中子曲线是否合格，并在判定所述中子曲线不合格的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的中子曲线，以及在判定所述中子曲线合格的情况下，对所述中子曲线进行归一化；

在判定所述测井曲线中包括所述密度曲线的情况下，判断该密度曲线是否合格，并在判定所述密度曲线不合格的情况下，利用归一化的自然伽玛曲线和声波曲线，重构归一化的密度曲线，以及在判定所述密度曲线合格的情况下，对所述密度曲线进行归一化。

8.根据权利要求7所述的准致密储层多指标评价方法，其特征在于，通过以下方式重构归一化的中子曲线和归一化的密度曲线：

NCNL_ri＝0.5707*NAC_i+0.4114*NGR_i-0.001

NDEN_ri＝0.1092-0.6922*NAC_i+0.1254*NGR_i

其中，NAC_i表示归一化的声波曲线中第i个储层深度下的声波测井值；

NGR_i表示归一化的自然伽玛曲线中第i个储层深度下的自然伽玛测井值；

NCNL_ri表示重构的归一化的中子曲线中第i个储层深度下的中子测井值；以及

NDEN_ri表示重构的归一化的密度曲线中第i个储层深度下的密度测井值。

9.根据权利要求7或8所述的准致密储层多指标评价方法，其特征在于，该准致密储层多指标评价方法还包括：

在重构归一化的中子曲线之后，利用所重构的归一化的中子曲线重构中子曲线；和/或

在重构归一化的密度曲线之后，利用所重构的归一化的密度曲线重构密度曲线。

10.根据权利要求9所述的准致密储层多指标评价方法，其特征在于，通过以下方式重构所述中子曲线和所述密度曲线：

CNL_ri＝0.1566+0.0639*NCNL_ri

DEN_ri＝2.4929+0.1305*NDEN_ri

其中，CNL_ri表示重构的中子曲线中第i个储层深度下的中子测井值；

DEN_ri表示重构的密度曲线中第i个储层深度下的密度测井值；

NCNL_ri表示重构的归一化的中子曲线中第i个储层深度下的中子测井值；以及

NDEN_ri表示重构的归一化的密度曲线中第i个储层深度下的密度测井值。

11.根据权利要求7或8所述的准致密储层多指标评价方法，其特征在于，所述测井曲线还包括由每个储层深度下的浅电阻率测井值形成的浅电阻率曲线；以及所述准致密储层多指标评价方法还包括：

对所述深电阻率曲线和所述浅电阻率曲线进行归一化；

根据归一化的自然伽玛曲线、归一化的声波曲线、归一化的深电阻率曲线、以及每个储层深度下的泥质含量，确定该储层深度下的高角度缝指标的大小；

根据归一化的深电阻率曲线、归一化的浅电阻率曲线、未经归一化的深电阻率曲线、未经归一化的浅电阻率曲线、以及每个储层深度下的泥质含量，确定该储层深度下的低角度缝指标的大小；以及

根据每个储层深度下的高角度缝指标的大小和低角度缝指标的大小，对该储层深度下的裂缝发育情况进行评价，以对所述储层的裂缝发育情况进行评价。

12.根据权利要求11所述的准致密储层多指标评价方法，其特征在于，通过以下方式确定所述高角度缝指标的大小：

通过以下方式确定所述低角度缝指标的大小：

其中，VFRI_i表示第i个储层深度下的高角度缝指标；

VCL_i表示第i个储层深度下的泥质含量；

DMax_i表示第i个储层深度下，归一化的自然伽玛曲线的一阶导数、归一化的声波曲线的一阶导数、归一化的深电阻率曲线的一阶导数中的最大值；

DMin_i表示第i个储层深度下，归一化的自然伽玛曲线的一阶导数、归一化的声波曲线的一阶导数、以及归一化的深电阻率曲线的一阶导数中导数间差值的最小值；

HFRI_i表示第i个储层深度下的低角度缝指标；

NRD_i表示归一化的深电阻率曲线中第i个储层深度下的深电阻率测井值；

NRS_i表示归一化的浅电阻率曲线中第i个储层深度下的浅电阻率测井值；

AvLRD_i表示未经归一化的深电阻率曲线中所有深电阻率测井值的对数均值；以及

AvLRS_i表示未经归一化的浅电阻率曲线中所有浅电阻率测井值的对数均值。

13.根据权利要求7或8所述的准致密储层多指标评价方法，其特征在于，该准致密储层多指标评价方法还包括：

根据归一化的声波曲线、归一化的中子曲线和每个储层深度下的泥质含量，确定该储层深度下的天然气挖掘效应指标的大小；以及

根据每个储层深度下的天然气挖掘效应指标的大小，对该储层深度下的流体性质进行定性识别，以对所述储层的流体性质进行评价。

14.根据权利要求13所述的准致密储层多指标评价方法，其特征在于，通过以下方式确定所述天然气挖掘效应指标的大小：

其中，CAVI_i表示第i个储层深度下的天然气挖掘效应指标；

VCL_i表示第i个储层深度下的泥质含量；

NAC_i表示归一化的声波曲线中第i个储层深度下的声波测井值；以及

NCNL_i表示归一化的中子曲线中第i个储层深度下的中子测井值。