[发明专利]一种基于叠前谱反演优化的高分辨率中深层储层预测方法

权利要求书

1.一种基于叠前谱反演优化的高分辨率中深层储层预测方法，其特征在于：包括以下步骤：

S1：对靶区内中深层所钻遇的储层及非储层的测井曲线进行岩石物理分析，确定能够用于识别储层的敏感弹性参数；

S2：地震叠前角度道集预处理，包括根据井中实测的纵横波速度和密度曲线进行的AVO正演模拟，分析井旁地震的AVO与正演的AVO趋势变化的差异性，对实际道集进行角度域转化，选取吻合度高的角度段进行分角度叠加，得到小、中、大三个角度的部分叠加数据体；

S3：对步骤S2中三个数据体的部分角度叠加分别进行子波提取，应用提取的子波进行谱反演，提高角度叠加地震资料的分辨率；

S4：对步骤S3得到的反射率体进行宽频带俞氏子波滤波整形，得到高分辨率的角度叠加数据体，并对三个角度的数据体与原始数据体进行质控，确保振幅变化的稳定性；

S5：将步骤S4所做的结果作为输入，进行三个角度的叠后反演，目标曲线为对应的弹性阻抗曲线，得到三个角度的弹性阻抗体；

S6：根据Connoly提出的弹性阻抗与弹性参数的关系，构建已知量和未知量的参数矩阵，采用奇异值分解方式最终得到稳定的弹性参数数据体，其中：

所述步骤S3中，通过对部分角度叠加的地震资料提取相对应的子波，并在子波约束下进行谱反演，提高叠前反演的纵向分辨率，在无噪声情况下，叠后地震资料在时间域表示为反射系数与地震子波的褶积形式：

s(t)＝w(t)*r(t) (1)

在水平层状介质的假设下，入射角在θ₁～θ₂度范围内的地震资料的叠加看做是平均角度位置的反射系数谱，任意角度的偏移后剖面表示为该角度的反射系数与该角度对应的地震子波的褶积形式：

s_θ(t)＝w_θ(t)*r_θ(t) (2)

对地震数据进行角度域的道集划分及部分角度的叠加处理，得到三个角度的地震资料，表示为三个角度的反射系数与地震子波的褶积形式，通过传统的谱反演方法理论，将公式(2)变换到频率域：

S_θ(f)＝W_θ(f)×R_θ(f) (3)

对公式(2)采用复数表示，任意角度对应的谱反演目标泛函表示为：

2.根据权利要求1所述的基于叠前谱反演优化的高分辨率中深层储层预测方法，其特征在于，所述步骤S1中，分别对储层和非储层的测井曲线进行岩石物理分析，得到能够区分泥岩与其他岩性的纵波阻抗参数，进一步对其他岩性中的砂岩和特殊岩性进行分析，得到区分二者的纵横波速度比参数，确立了采用纵波阻抗和纵横比双参数进行优质储层预测的前提条件。

3.根据权利要求1所述的基于叠前谱反演优化的高分辨率中深层储层预测方法，其特征在于，所述步骤S4中，宽频带俞氏子波是一种宽频带子波，其时间域可以写为如下形式：

4.根据权利要求1所述的基于叠前谱反演优化的高分辨率中深层储层预测方法，其特征在于，所述步骤S5中，对三个角度的宽频带数据进行弹性阻抗反演，得到各个角度对应的绝对弹性阻抗，将其写成矩阵的形式，输入三个角度的弹性波阻抗，得到储层的敏感参数：