[发明专利]一种地震速度反演方法

说明书

本发明公开了一种地震速度反演方法，包括：输入初始的观测地震数据、初始的速度模型与极性编码矩阵；对观测炮集利用极性编码矩阵合成编码超级炮集；利用编码超级炮集与极性编码矩阵生成加权编码的参考道；利用参考道和目标子波对超级炮记录采用维纳滤波生成目标子波对应的超级炮；利用目标子波和编码函数与速度场生成正演超级炮；利用正演超级炮与生成的超级炮记录做差反传，得到反演梯度；利用抛物插值计算最佳步长，并对速度场进行更新；重复S2到S7直到满足迭代次数或者速度场要求，输出最终的速度场。本发明对维纳滤波采用的参考道是一种加权的多震源地震道，其不需要估计子波，减少子波估计的计算量。

技术领域

本发明涉及一种地震速度反演方法，属于油气物探工程领域。

背景技术

在当前的勘探领域中，地震速度场的反演对于地震数据的偏移成像和解释具有重要的意义，全波形反演能够充分利用地震数据中所包含的振幅、相位、走时、波形等信息，对地下模型参数进行精细刻画。

在当前，通过维纳滤波实现多尺度的全波形反演是一种常用的手段。但是在这种方式下，维纳滤波必须使用原始子波，存在原始子波依赖性强、计算量巨大、反演稳定性等诸多问题，由于采集成本、坏道、噪声、地形等诸多因素的影响，实际采集到的数据往往提取的原始子波不准确，给地震数据的速度反演造成了一定的困扰。

因此，发展一种不依赖子波的快速地震速度反演方法是十分必要的。

发明内容

本发明的目的在于，提供一种地震速度反演方法，它可以解决当前技术中存在的问题，对不知道地震子波的地震数据，提高地震数据反演速度场的精度。

为解决上述技术问题，本发明提供一种地震速度反演方法，包括以下步骤：

S1，输入初始的观测地震数据、初始的速度模型与极性编码矩阵；

S2，对观测炮集利用极性编码矩阵合成编码超级炮集；

S3，利用编码超级炮集与极性编码矩阵生成加权编码的参考道；

S4，利用参考道和目标子波对超级炮记录采用维纳滤波生成目标子波对应的超级炮；

S5，利用目标子波和编码函数与速度场生成正演超级炮；

S6，利用正演超级炮与S4生成的超级炮记录做差反传，互相关得到反演梯度；

S7，利用抛物插值计算最佳步长，并对速度场进行更新；

S8，重复S2到S7直到满足迭代次数或者速度场要求，输出最终的速度场。

与现有技术相比，本发明对维纳滤波采用的参考道是一种加权的多震源地震道，其不需要估计子波，减少子波估计的计算量。

附图说明

附图是用来提供对本发明的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明，但并不构成对本发明的限制。在附图中：

图1为本发明实施例的流程示意图；

图2为本发明的编码矩阵示意图；

图3为本发明使用的速度场模型图；

图4为本发明使用的单炮记录；

图5为本发明实施例的初始速度场示意图；

图6为本发明的观测数据中的合成超级炮记录示意图；

图7为使用本发明方法生成的参考道；

图8为使用本发明得到目标子波对应的超级炮记录；

图9为使用本发明得到反演速度场；

图10为不准确子波得到反演速度场；