[发明专利]一种碳酸盐岩层间岩溶储层预测方法及装置

权利要求书

1.一种碳酸盐岩层间岩溶储层预测方法，其特征在于，包括如下步骤：

S101、获取地震纯波数据和目的层的相对应时间构造层位数据；

S102、通过对地震纯波数据进行相对波阻抗反演，获得相对波阻抗；

S103、通过对地震纯波数据进行地层吸收因子计算，获得地层吸收因子数据体；

S104、构建相对波阻抗与地层吸收因子数据体的融合属性计算公式；

S105、利用最小二乘法优化融合属性计算公式，获得优化后的融合属性；

S106、利用优化后的融合属性对储层进行预测；

步骤S104中所述相对波阻抗与地层吸收因子数据体的融合属性计算公式具体为：

式II中，

F(t)表示融合属性；

Ip(t)表示相对波阻抗；

Q表示吸收因子；

α表示低频数值；

β表示频率滑动视窗；

步骤S105中所述利用最小二乘法优化融合属性计算公式，获得优化后的融合属性，具体包括：

统计实际钻井储层测井的相对波阻抗和地层吸收因子数据体，利用最小二乘法优化融合属性计算公式，获得油气层与水层分界门槛值以及优化后的频率滑动视窗；

根据优化后的频率滑动视窗，获得优化后的地层吸收因子数据体；

将优化后的地层吸收因子数据体结合融合属性计算公式，获得优化后的融合属性。

2.根据权利要求1所述的碳酸盐岩层间岩溶储层预测方法，其特征在于，步骤S101还包括获取其他数据，所述其他数据包括测井声波、密度、电阻率、伽马和孔隙度数据。

3.根据权利要求1所述的碳酸盐岩层间岩溶储层预测方法，其特征在于，步骤S102中所述通过对地震纯波数据进行相对波阻抗反演，获得相对波阻抗，具体为：

基于褶积模型，通过对地震纯波数据利用道积分进行相对波阻抗反演，得到相对波阻抗的计算公式如式I，根据式I计算获得相对波阻抗；

ln Ip(t)＝2∫S(t)dt I；

式I中，

S(t)表示地震信号；

Ip(t)表示相对波阻抗；

t表示时间。

4.根据权利要求1所述的碳酸盐岩层间岩溶储层预测方法，其特征在于，所述利用最小二乘法优化融合属性计算公式，获得优化后的融合属性，具体通过求解下式III进行优化：

式III中，

m表示实际钻井储层测井的数目；

Ip(t)表示相对波阻抗；

Q表示吸收因子；

α表示低频数值；

β表示频率滑动视窗；

ε表示油气层与水层分界门槛值。

5.根据权利要求1所述的碳酸盐岩层间岩溶储层预测方法，其特征在于，步骤S106中所述利用优化后的融合属性对储层进行预测，具体为将优化后的融合属性和油气层与水层分界门槛值进行比较，得出预测结论。

6.根据权利要求5所述的碳酸盐岩层间岩溶储层预测方法，其特征在于，所述优化后的融合属性大于油气层与水层分界门槛值的地方为含油气地区。