[发明专利]一种叠后声波阻抗反演方法

权利要求书

1.一种叠后声波阻抗反演方法，其特征在于：所述方法包括以下步骤：

(1)对叠后地震数据进行处理，使得叠后地震数据的频谱与由测井数据得到的声波阻抗的频谱相匹配，从而获得相对声波阻抗数据；

(2)基于步骤(1)得到的相对声波阻抗数据进行构造解释，以构建用于基于模型反演的初始声波阻抗模型或者变差函数，如果构建的是变差函数，则再由变差函数进行随机模拟得到初始声波阻抗模型；

(3)对步骤(2)得到的初始声波阻抗模型进行正演计算得到合成地震数据；

(4)将合成地震数据与实际叠后地震数据进行比较，求得声波阻抗模型的修改量，利用该修改量对初始声波阻抗模型进行修改更新得到更新后的声波阻抗模型，再对更新后的声波阻抗模型进行正演计算得到合成地震数据；然后将更新后的声波阻抗模型作为初始声波阻抗模型，重复步骤(4)直到合成地震数据与实际叠后地震数据在最小平方意义下最接近时，此时的声波阻抗模型即为最终的反演结果；

所述步骤(2)中构建用于基于模型反演的初始声波阻抗模型是这样实现的：以构造解释结果为框架进行约束，然后在该框架内填入所述相对声波阻抗数据，最后进行网格化，形成网格化的初始声波阻抗模型；

所述步骤(2)中构建变差函数是这样实现的：

对于每一个网格点，首先根据储层的平面展布以及厚度分别确定横向变程范围和垂向变程范围，然后分别在横向和垂向上将此范围内的所述相对声波阻抗值数据进行统计，获得该点纵向和垂向实验变差函数曲线，最后利用变差函数模型对该曲线进行拟合，以确定变差函数模型中的参数，从而最终得到该点的变差函数。

2.根据权利要求1所述的叠后声波阻抗反演方法，其特征在于：所述步骤(1)具体包括：

(11)由测井数据计算得到所有井上的声波阻抗：将声波测井数据与密度测井数据相乘；

(12)在以对数频率为横坐标、声波阻抗为纵坐标的坐标轴的图上进行线性拟合，得到声波阻抗频谱；

(13)利用所述声波阻抗频谱对地震数据的频谱进行整形：设计一个滤波器，滤波器的脉冲响应为声波阻抗频谱，利用该滤波器对地震数据进行滤波处理；

(14)对经过整形后的地震数据进行-90度相位旋转，获得相对声波阻抗数据。

3.根据权利要求2所述的叠后声波阻抗反演方法，其特征在于：所述步骤(2)中进行构造解释是指解释地质层位与断层。

4.根据权利要求1所述的叠后声波阻抗反演方法，其特征在于：所述步骤(3)是这样实现的：

对步骤(2)得到的初始声波阻抗模型，由公式(1)得到反射系数模型，然后由公式(2)进行正演计算得到合成地震数据：

Ri=Zi+1-ZiZi+1+Zi---(1)]]>

e＝D-w·R(2)

公式(1)中的Z表示所述相对声波阻抗数据，R表示反射系数，公式(2)中的D表示实际叠后地震数据，w表示子波，i表示第i层。

5.根据权利要求4所述的叠后声波阻抗反演方法，其特征在于：所述步骤(4)中所述求得声波阻抗模型的修改量是利用公式(5)实现的：

g≈w^T(D-wR)(5)

公式(5)中，w^T表示子波矩阵的转置。