[发明专利]一种岩性圈闭有效性识别的方法

权利要求书

1.一种岩性圈闭有效性识别的方法，其特征在于，包括以下步骤：

获取目标地层内不同地层深度对应的突破压力和排驱压力；

根据所述目标地层内已有钻井的纵波声波时差测井曲线和密度测井曲线，获得不同地层深度的纵波阻抗数值；建立纵波阻抗数值与排驱压力的关系式、纵波阻抗数值与突破压力的关系式，并获得所述目标地层中盖层的封堵性表征关系式；

根据所述目标地层中盖层的岩性分布，获得不同岩性的盖层的纵波阻抗数值分布范围；

根据所述目标地层的三维地震数据体得到所述目标地层的纵波阻抗三维数据体，并根据所述目标地层中盖层的所述纵波阻抗数值的分布范围筛选不同岩性的盖层，得到所述目标地层中的盖层厚度分布图；

根据所述封堵性表征关系式和所述三维地震数据体获得封堵性参数三维数据体，并根据所述封堵性参数三维数据体，筛选出合格的封堵层，得到所述目标地层的封堵层厚度有效界限值和所述目标地层的封堵层厚度分布图；

根据所述目标地层的三维地震数据体，计算所述目标地层的曲率有效界限值和所述目标地层的曲率分布图；

建立所述目标地层的岩性圈闭有效性参数计算式，并根据所述目标地层的盖层厚度有效界限值、封堵层厚度有效界限值、曲率有效界限值得到所述目标地层的岩性圈闭有效性阈值，得到所述目标地层的有效圈闭发育区分布图；

所述岩性圈闭有效性参数计算式为：

其中，f(C)为盖层厚度；f(S)为封堵层厚度，f(F)为曲率。

2.根据权利要求1所述的岩性圈闭有效性识别的方法，其特征在于，所述封堵层厚度有效界限值为5m；

所述盖层分为泥页岩层、膏岩层和致密灰岩层；当所述盖层为泥页岩层时，所述盖层厚度有效界限值为5m；当所述盖层为膏岩层时，所述盖层有效界限值为5m，当所述盖层为致密灰岩层时，所述盖层有效界限值为100m；

所述曲率有效界限值为所述目标地层内断层处曲率的最小值。

3.根据权利要求1所述的岩性圈闭有效性识别的方法，其特征在于，所述目标地层的有效圈闭发育区分布图获得步骤如下：

根据所述岩性圈闭有效性参数计算式，并带入所述目标地层的所述盖层厚度、所述封堵层厚度、所述曲率，获得所述目标地层内的岩性圈闭有效性参数值分布图；

再通过所述目标地层的盖层厚度有效界限值、封堵层厚度有效界限值、曲率有效界限值计算得到所述目标地层的岩性圈闭有效性阈值；

判断所述岩性圈闭有效性参数值与所述岩性圈闭有效性阈值的大小，得到所述目标地层的有效圈闭发育区分布图。

4.根据权利要求1所述的岩性圈闭有效性识别的方法，其特征在于，所述纵波阻抗数值通过以下公式计算得到：

IMP＝1000000/AC×DEN

其中：IMP为纵波阻抗数值，g/cm³*m/s；AC为纵波声波时差，μs/m；DEN为补偿密度，g/cm³。

5.根据权利要求1或2所述的岩性圈闭有效性识别的方法，其特征在于，建立纵波阻抗数值与排驱压力的关系式、纵波阻抗数值与突破压力的关系式，筛选出相关性系数最高的关系式作为封堵性表征关系式，包括以下步骤：

获取不同所述地层深度的所述纵波阻抗数值和对应的排驱压力、获取不同所述地层深度的所述纵波阻抗数值IMP和对应的突破压力；

分别采用不同的函数拟合方法拟合所述纵波阻抗数值和所述排驱压力、所述纵波阻抗数值和所述突破压力，根据拟合得到的函数关系式的相关性系数筛选出相关性系数最高的关系式，作为封堵性表征关系式；

所述函数拟合方法包括：线性拟合、多项式拟合、指数拟合、对数拟合、幂函数拟合。

6.根据权利要求2所述的岩性圈闭有效性识别的方法，其特征在于，所述盖层厚度为所述泥页岩层、膏岩层、致密灰岩层的厚度之和；所述目标地层内不同岩性的盖层分布根据所述目标地层的录井岩屑、钻井岩芯资料获取。