[发明专利]致密滩坝砂储层评价方法

权利要求书

1.一种致密滩坝砂储层评价方法，其特征在于，该方法包括以下步骤：

利用压汞实验数据和润湿相模型计算单个致密滩坝砂储集层样品点的孔隙分形维数；

在目的层坝砂与滩砂储集层分类基础上，研究坝砂与滩砂测井参数特点，分析孔隙分形维数与各测井参数相关性，并优选出相关度较高的参数；

将所述测井参数归一化后回归拟合分别建立致密坝砂与致密滩砂的分形维数测井预测模型；

根据已建立的所述致密坝砂与致密滩砂的分形维数测井预测模型对单井纵向上孔隙分形维数值进行预测；

根据所述单井纵向上的孔隙分形维数值，分小层加权平均计算孔隙分形维数值，获得该井所对应的平面孔隙分形维数值；

选取已动用区块的开发效果，分析其与孔隙分形维数之间相关性，建立基于孔隙分形维数的致密滩坝砂储层评价标准；

根据所述孔隙分形维数的致密滩坝砂储层评价标准，对致密滩坝砂的储层进行分类，确定开发动用次序。

2.如权利要求1所述的致密滩坝砂储层评价方法，其特征在于，所述步骤利用压汞实验数据和润湿相模型计算单个致密滩坝砂储集层样品点的孔隙分形维数，具体包括：根据毛管压力和湿相饱和度关系，获得储层毛管压力曲线的分形几何公式为：

式中，V为进汞过程中润湿相饱和度，Pmin为最大孔喉半径对应的毛管压力，单位为MPa，Pc为毛管压力，单位为MPa，结合拉普拉斯方程，公式为：

r为孔隙半径，单位为μm；对其两边取对数，得

lgV＝lg(1-SHg)＝(3-D)lgr-(3-D)lgrmax，

式中，SHg为累积进汞饱和度；

通过绘制lg(1-SHg)和lgr分形曲线，拟合出曲线的斜率K，即可根据D＝3-K求取孔喉的分形维数，计算过程中，分形维数出现两段式，此时，需要根据各孔隙空间的孔隙度加权平均得到整个孔喉的总分形维数D，公式为：

式中，表示第一段孔隙对应的孔隙度，D1表示第一段孔隙对应的分形维数；表示第二段孔隙对应的孔隙度，D2表示第二段孔隙对应的分形维数。

3.如权利要求1或2所述的致密滩坝砂储层评价方法，其特征在于，所述步骤在目的层坝砂与滩砂储集层分类基础上，研究坝砂与滩砂测井参数特点，分析孔隙分形维数与各测井参数相关性，并优选出相关度较高的参数，具体包括：研究坝砂与滩砂的测井参数特点，将归一化后的所述测井参数投影到雷达图中，分析孔隙分形维数与各测井参数相关性。

4.如权利要求3所述的致密滩坝砂储层评价方法，其特征在于，所述步骤将所述测井参数归一化后回归拟合分别建立致密坝砂与致密滩砂的分形维数测井预测模型，其中，所述致密坝砂分形维数测井预测模型为：

D_坝砂＝2.90+0.13GR-0.02DEN-0.44RT-0.14K，

致密滩砂分形维数测井预测模型为：

D_滩砂＝2.95+0.06GR-0.30AC-0.39RT-0.01K，

其中，D_坝砂、D_滩砂为分形维数，GR为自然伽马值，DEN为岩性密度值，单位为g/cm³；RT为电阻率值，单位为Ω·m；K为微电极曲线差值绝对值，单位为Ω·m；AC为声波时差值，单位为μs/ft。

5.如权利要求4所述的一种致密滩坝砂储层评价方法，其特征在于，所述步骤根据所述单井纵向上的孔隙分形维数值，分小层加权平均计算孔隙分形维数值，获得该井所对应的平面孔隙分形维数值，具体包括：确定单井同一层位致密滩坝砂储集层平均孔隙分形维数值，对段内的坝砂与滩砂分形维数分别进行计算，然后加权平均获得对应平面的孔隙分形维数值。所述平均孔隙分形维数计算公式为：

其中，D_坝为坝砂的分形维数值，H_坝为坝砂厚度，D_滩为滩砂的分形维数值，H_滩为滩砂厚度，H为坝砂和滩砂的总厚度，D为该井所对应平面上的分形维数值。