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

说明书

技术领域

本发明属于油气及煤层气地震勘探与开发领域，具体涉及一种叠后声波阻抗反演方法。

背景技术

随着油气勘探开发的需要，储层预测和精细描述越来越引起重视。围绕这个目的而开展的探索和研究也越来越多，地震反演就是最重要的一点。当反演的目标是波阻抗时，称之为波阻抗反演。目前的地震波阻抗反演按照所使用的资料分为叠后反演和叠前反演两大类，按照反演方法又可以分为直接反演和间接反演。直接反演就是直接从地震资料出发，进行运算得到声波阻抗。间接反演是从一个初始模型出发，合成地震记录，然后构造目标函数，最后求得目标函数在一定的范数意义下极值。无论哪种反演，其实质都是想去除子波的影响，从而将地震剖面转化为能与钻井、地质等资料直接对比的形式，因此反演在许多情况下提高了常规地震的分辨率并提高了油藏参数研究的水平。

叠后声波阻抗反演概括起来不外乎有两大类：基于反射系数逆公式的直接反演和基于正演模型的迭代反演。

基于反射系数逆公式的直接反演有道积分和递推反演。地震道积分近似等于对数波阻抗，但是该方法无法求得地层绝对波阻抗，而且使用时无法用地质或测井资料进行约束。以稀疏脉冲反演为代表的递推反演虽然更加精确和稳定，但是在实现和应用上都更复杂。

基于模型的反演可以分为确定性反演和随机反演两种。前者通常都是由测井、地质以及地震资料给出一个确定性的初始模型，然后迭代反演，得到与地震资料最佳匹配的地震波阻抗模型。这种反演方法，在井越多并且地震资料品质越好的情况下，给出的初始模型更接近实际，其结果也就越可靠。但是直接由叠后地震资料解释层位，一方面分辨率低，另一方面也受地震资料品质的影响。虽然目前广泛使用的随机反演可以通过地质统计的手段提供多个符合统计参数要求的模型，但是其横向变差函数还是受叠后地震资料品质的影响，即用于计算纵向变差函数的测井以及地质资料的分辨率与用于计算横向变差函数的地震资料的分辨率不匹配。

发明内容

本发明的目的在于解决上述现有技术中存在的难题，提供一种叠后声波阻抗反演方法，在进行基于模型的反演之前，对地震资料进行处理，能快速而稳健地得到分辨率更高的剖面，从而在进行初始模型的建立或者变差函数的计算时，能生成一个既满足测井和地质资料分辨率又满足地震资料的模型，然后再进行基于模型的反演。

本发明是通过以下技术方案实现的：

一种叠后声波阻抗反演方法，包括以下步骤：

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

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

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

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

所述步骤(1)具体包括：

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

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

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

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

所述步骤(2)中进行构造解释是指解释地质层位与断层。

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

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