[发明专利]一种油气储层渗透率预测方法及装置

权利要求书

1.一种油气储层渗透率预测方法，其特征在于，包括：

获取测井信息及叠前地震角度道集，并根据测井信息及叠前地震角度道集生成多条不同角度的弹性阻抗曲线；

对每条弹性阻抗曲线分别进行多种数学变换；

对数学变换后的每任意两条曲线进行角度旋转重组，得到多条新的属性曲线；

获取渗透率测井曲线，并将每条新的属性曲线与渗透率测井曲线分别进行曲线拟合，取拟合相关系数最大时对应的新属性作为目标属性；

确定目标属性的构成元素，并根据构成元素将叠前地震反演得到的叠前地震弹性阻抗反演数据体转化为渗透率数据体，利用渗透率数据体预测油气储层的渗透率。

2.如权利要求1所述的油气储层渗透率预测方法，其特征在于，所述测井信息包括：纵波速度曲线、横波速度曲线和密度曲线；

叠前地震角度道集中的地震波的入射角范围为-90°-90°；

根据测井信息生成多条不同角度的弹性阻抗曲线，包括：

在地震波入射角范围内，根据纵波速度曲线、横波速度曲线和密度曲线，每隔预设角度生成一条弹性阻抗曲线。

3.如权利要求2所述的油气储层渗透率预测方法，其特征在于，所述弹性阻抗曲线X(θ_i)的计算公式如下：

A＝V_pρ；

A₀＝V_p0ρ₀；

其中，θ_i为第i条弹性阻抗曲线对应的地震波入射角的角度；V_p为纵波速度，单位为米/秒；V_s为横波速度，单位为米/秒；ρ为密度，单位为千克/立方米；V_p0为目的层段纵波速度，单位为米/秒；V_s0为目的层段横波速度，单位为米/秒；ρ₀为密度的平均值，单位为千克/立方米；K为横波速度与纵波速度之比，无量纲；M为目的层段平均纵波阻抗A₀的数量级。

4.如权利要求3所述的油气储层渗透率预测方法，其特征在于，所述数学变换的方式包括：

保持不变、以自然指数为底的对数变换、以自然指数为底的幂变换、倒数变换、平方变换、开方变换。

5.如权利要求1所述的油气储层渗透率预测方法，其特征在于，所述新的属性曲线T_pq(χ)的计算公式如下：

其中，χ为旋转角度；记数学变换后的多条新的属性曲线的集合为curve-A，A_p、A_q分别为曲线集curve-A中的第p、q条属性曲线；分别为曲线集curve-A中的第p、q条属性曲线平均值的数量级；T_pq(χ)为曲线集curve-A中的第p、q条属性曲线旋转χ度重组得到的新的属性曲线。

6.如权利要求5所述的油气储层渗透率预测方法，其特征在于，所述将每条新的属性曲线与渗透率测井曲线分别进行曲线拟合所采用的拟合公式为：

P＝aT²_pq(χ)+bT_pq(χ)+c

其中，P表示已知渗透率测井曲线，a，b，c分别表示二项式拟合常系数。

7.如权利要求6所述的油气储层渗透率预测方法，其特征在于，取拟合相关系数最大时对应的新属性作为目标属性，包括：

将每一条新的属性曲线通过二项式拟合，得到一组二项式拟合常系数；

将每组常系数与对应的新的属性曲线带入到二项式公式，计算得到每条新的属性曲线对应的新的渗透率预测曲线；

根据相关系数公式，分别计算每条新的渗透率预测曲线与已知渗透率测井曲线的相关系数，取相关系数最大值对应的新属性作为目标属性。