[发明专利]一种地下构造反演成像方法

说明书

本发明公开了一种地下构造反演成像方法，包括以下步骤：S1，输入背景速度场、观测数据、震源，获取初始反射系数模型；S2，对所述观测数据进行编码，获得超炮集；对震源进行编码，获得编码后的震源；S3，根据编码后的震源和初始反射系数模型，获得正演模拟数据；S4，获得正演模拟数据与超炮集之间的数据残差；S5，根据数据残差确定梯度，基于梯度进行反射系数模型迭代，获得最终的反射系数模型。本发明通过对观测数据和震源采用分段阶梯编码，充分加大多震源炮数据之间的差异性，减少临炮之间的干扰，在多震源炮数据成像时较为好的压制串扰噪音，提高了多震源炮数据成像反演时的成像精度。

技术领域

本发明涉及一种地下构造反演成像方法，属于油气物探工程领域。

背景技术

在当前的勘探领域中，地下数据采集技术尤为关键，这其中多震源数据采集方式开始得到关注，所谓多震源数据是多震源同时激发，波场在地下介质中传播，在地面接受得到的数据，其是地下介质对多震源的响应。多震源数据采集方式改善了传统采集方式照明不足的问题，也大大降低了三维数据采集成本。但是当前在对多震源数据进行处理时，由于不相关炮的波场进行互相关而产生的串扰噪声，由于不能有效压制这些串扰噪声，最终影响了成像质量。因此，在对这样的多震源数据做处理的时候，压制串扰噪音以提高其成像精度是非常必要的。

发明内容

本发明的目的在于，提供一种地下构造反演成像方法，它可以解决当前技术中存在的问题，有效压制多震源数据中之间的串扰噪声，提高反演成像精度。

为解决上述技术问题，本发明采用如下的技术方案：一种地下构造反演成像方法，包括以下步骤：

S1，输入背景速度场、观测数据和震源，获取初始反射系数模型；

S2，对所述观测数据进行编码，获得超炮集；对震源进行编码，获得编码后的震源；

S3，根据编码后的震源和初始反射系数模型，获得正演模拟数据；

S4，获得正演模拟数据与超炮集之间的数据残差；

S5，根据数据残差确定梯度，基于梯度进行反射系数模型迭代，获得最终的反射系数模型。

与现有技术相比，本发明对观测数据和震源采用分段阶梯编码，充分加大多震源炮数据之间的差异性，减少临炮之间的干扰，在多震源数据成像时较为好的压制串扰噪音，提高了多震源数据成像反演时的成像精度。此外，利用分段阶梯编码可以直接对多震源数据进行偏移，而不需要对每炮数据进行单独偏移成像再叠加，从而提高了反演的效率。

附图说明

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

图1为本发明的编码示意图；

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

图3为本发明实施例的背景速度场示意图；

图4为本发明实施例观测数据中的单炮观测记录示意图；

图5为本发明实施例中超炮集的记录示意图；

图6为本发明实施例第5次迭代得到的反射系数模型示意图；

图7为本发明实施例最终的反射系数模型示意图；

图8为常规方法得到最终的反射系数模型示意图；

图9为图7和图8的局部放大后的对比，其中图(a)和(c)为常规方法所得，图(b)(d)为本发明所用的方法所得；

图10为常规方法与本发明方法的迭代误差曲线对比；

下面结合附图和具体实施方式对本发明作进一步的说明。