[发明专利]一种全波形速度场约束下的盐下碳酸盐岩储层反演方法

权利要求书

1.一种全波形速度场约束下的盐下碳酸盐岩储层反演方法，其特征在于，对于目标工区基于全波形反演和岩石物理分析的结果进行叠前反演得到IP-IS弹性参数体，实现盐下碳酸盐岩储层精细预测；包括以下步骤：

1)目标工区岩石物理分析及井震标定；

2)全波形速度场约束下的目标工区低频模型的构建；

3)基于部分叠加地震数据进行叠前同时反演；

4)将步骤3)得到的反演结果进一步处理，得到IP-IS弹性参数体，作为反演的最终结果，对储层进行精细预测；

所述步骤2)中采用全波形速度场约束下的多井插值法进行低频模型构建，具体过程包括：

21)首先，速度场应为全波形反演得到的速度场，并对其进行去野值、平滑、插值处理；

22)提取目的层段测井纵波阻抗曲线，对其低通滤波，保留10hz以下有效成分，分别统计目的层段全波形速度以及低通滤波后测井纵波阻抗的均值和方差，利用正态分布参数估计公式，将全波形速度场转化为测井纵波阻抗数据体；

23)应用简单克里格法消除井点位置的纵波阻抗预测误差，得到纵波阻抗低频模型；

24)分别对目的层段测井纵波阻抗与密度，纵波阻抗与横波阻抗进行交会分析，并进行线性回归，利用线性回归公式将步骤23)得到的纵波阻抗低频模型转换为横波阻抗低频模型和密度低频模型，做为叠前地震反演所需的低频模型；

所述步骤3)中基于部分叠加地震数据进行叠前同时反演具体包括：

31)基于部分叠加地震数据体，采用Aki-Richard近似公式进行叠前同时反演计算，得到缺乏低频成分的弹性参数体，包括纵波阻抗、横波阻抗和密度；

32)采用频率域合并的方式，将步骤2)得到的低频模型频带范围内的值，合并到步骤31)得到的反演结果中，补充原始地震缺失的低频成分；在所述步骤31)中，Aki-Richard近似公式如下：

式中，R_pp(θ)为纵波的反射系数，v_p和v_s分别为上层介质的纵、横波速度，Δv_p和Δv_s分别为上、下层介质的纵、横波速度差，θ为地震波的入射角。

2.根据权利要求1所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法，其特征在于，所述步骤1)中还包括，在测井资料一致性处理及横波预测的基础上，进行岩石物理分析，明确能够区分储层及非储层的反演敏感弹性参数；根据测井数据、部分叠加地震数据进行井震标定，提取多井子波。

3.根据权利要求1所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法，其特征在于，所述步骤22)中，正态分布参数估计公式如下：

result＝input*(a/b)+(c-d)*(a/b) (1)

式中，a为低通滤波后测井纵波阻抗的标准差，b为目的层段全波形速度的标准差，c为低通滤波后测井纵波阻抗的均值，d为目的层段全波形速度的均值。

4.根据权利要求1所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法，其特征在于，所述步骤4)中将步骤32)得到的反演结果进一步处理具体为，将步骤3)得到的补充低频成分后的纵波阻抗减去横波阻抗，得到IP-IS弹性参数体，作为反演的最终结果。