[发明专利]全波形反演中分时窗提取数据的方法及装置

说明书

本发明提供了一种全波形反演中分时窗提取数据的方法及装置。该方法包括：获取地震数据及预设速度模型，利用声波波动方程及速度模型进行正演模拟，得到正演地震数据；利用初至拾取方式拾取初至时间，估算面波同向轴的平均视速度，确定面波到达检波器时间；利用初至时间，确定第一及第二时间窗口上限值，以及第一及第二时间窗口下限值，以确定第一时窗函数及第二时窗函数；利用面波到达检波器的时间，确定第三时间窗口上限值及第三时窗函数；根据第一时窗函数、第二时窗函数、第三时窗函数及地震数据，确定数据子体，利用数据子体进行全波形反演。本发明利用三个时间窗对地震炮集分时窗提取有效数据，有效去除噪声，为全波形反演提供高信噪比数据。

技术领域

本发明涉及全波形反演技术领域，尤指一种全波形反演中分时窗提取数据的方法及装置。

背景技术

全波形反演技术是定量地震成像的前沿性技术之一，该技术为高精度地震成像提供了可靠的速度模型。

由于近地表环境的影响，陆上采集数据往往受到低频噪声面波的影响。面波在地震炮集上的形态为一个以炮点为中心点的扇形区域。这部分数据在地震预处理阶段通过滤波的方法进行消除。当面波能量很大，而且引起次生强噪声的时候，常规的处理方法不能有效去除该区域的噪声，提高信噪比，因此往往导致全波形反演很难收敛。另外，由于地震炮集中初至波与反射波能量级别不一样，全波形反演中往往把早至波与反射波分开来进行操作。

发明内容

为了解决上述问题，本发明实施例提供一种全波形反演中分时窗提取数据的方法，所述方法包括：

获取地震数据及预设速度模型，利用声波波动方程及所述速度模型进行正演模拟，得到正演地震数据；

利用初至拾取方式拾取初至时间，并估算面波同向轴的平均视速度，利用所述平均视速度确定面波到达检波器的时间；

利用所述初至时间，确定第一时间窗口上限值及第一时间窗口下限值，以及第二时间窗口上限值及第二时间窗口下限值，根据所述第一时间窗口上限值及第一时间窗口下限值确定第一时窗函数，并根据所述第二时间窗口上限值及第二时间窗口下限值确定第二时窗函数；

利用所述面波到达检波器的时间，确定第三时间窗口上限值，根据所述第三时间窗口上限值确定第三时窗函数；

根据所述第一时窗函数、第二时窗函数、第三时窗函数及地震数据，确定对应的数据子体，利用所述数据子体进行全波形反演。

可选的，在本发明一实施例中，所述利用所述初至时间，分别确定第一时间窗口上限值及第一时间窗口下限值，以及第二时间窗口上限值及第二时间窗口下限值包括：估算初至波的最大周期以及面波的最大周期；利用所述初至时间、初至波的最大周期及炮点到检波点的距离，确定第一时间窗口上限值及第一时间窗口下限值，以及第二时间窗口上限值；利用所述初至时间、面波的最大周期及炮点到检波点的距离，确定及第二时间窗口下限值。

可选的，在本发明一实施例中，所述根据所述第一时间窗口上限值及第一时间窗口下限值确定第一时窗函数包括：根据所述第一时间窗口上限值、第一时间窗口下限值及初至波的最大周期，确定所述第一时窗函数。

可选的，在本发明一实施例中，所述根据所述第二时间窗口上限值及第二时间窗口下限值确定第二时窗函数包括：根据所述第二时间窗口上限值、第二时间窗口下限值、初至波的最大周期以及面波的最大周期，确定所述第二时窗函数。

可选的，在本发明一实施例中，所述利用所述面波到达检波器的时间，确定第三时间窗口上限值包括：估算面波的最大周期，利用所述面波到达检波器的时间、面波的最大周期及炮点到检波点的距离，确定第三时间窗口上限值。

可选的，在本发明一实施例中，所述根据所述第三时间窗口上限值确定第三时窗函数包括：根据所述第三时间窗口上限值及面波的最大周期，确定第三时窗函数。