[发明专利]全波形反演中去除波动方程模拟直达波的方法及装置

说明书

技术领域

本发明涉及地震勘探技术领域，尤其涉及全波形反演中去除波动方程模拟直达波的方法及装置。

背景技术

全波形反演是目前地震勘探领域前沿技术之一。该技术需要多次基于双程波动方程的正演过程，并且要求正演波场与实际波场能够匹配。因此，波动方程正演模块波场的精度是影响全波形反演技术的重要因素之一。由于置放在地表的检波器所接收到正演模拟的直达波具有很强振幅，这些强振幅会导致逆时偏移的构造剖面或者全波形反演的弹性参数剖面浅层存在强振幅噪声，从而影响了偏移或反演的质量。常规去除直达波方法采用预先设定的时间窗，然后根据设定的阈值直接切除直达波。虽然方法能够成功应用于野外采集地震炮集，但是并不适用于波动方程模拟的波场，这是因为，基于速度模型，尤其是光滑速度模拟的地震波场，正演模拟的直达波能量比反射波能量高出很多倍，而且其影响的输入时间窗口较大，常规的直接切除直达波的方法只能切除局部的高能量数值，并不能很好的消除模拟直达波的高能量区域，导致反射波的信噪比较低。

发明内容

本发明的目的在于提供一种全波形反演中去除波动方程模拟直达波的方法及装置，以消除模拟直达波，增强反射波的信噪比。

为达到上述目的，一方面，本发明提供了一种全波形反演中去除波动方程模拟直达波的方法，包括以下步骤：

输入速度模型；

基于声波波动方程对所述速度模型进行波动方程正演，得到模拟的地震记录；

计算直达波到达地表接收器的时间T₀(x)，x为检波点位置到震源位置的距离；

根据公式计算输入时间窗口T(x)，α为控制所述输入时间窗口的拉伸系数；

删除所述地表接收器接收到的其接收时间小于所述T(x)时刻的地震记录。

再一方面，本发明还提供了一种全波形反演中去除波动方程模拟直达波的装置，包括：

输入模块，用于输入速度模型；

波动方程正演模块，用于基于声波波动方程对所述速度模型进行波动方程正演，得到模拟的地震记录；

第一计算模块，用于计算直达波到达地表接收器的时间T₀(x)，x为检波点位置到震源位置的距离；

第二计算模块，用于根据公式计算输入时间窗口T(x)，α为控制所述输入时间窗口的拉伸系数；

直达波切除模块，用于删除所述地表接收器接收到的其接收时间小于所述T(x)时刻的地震记录。

与现有常规直达波切除技术不同，本发明不再是根据设定的阈值直接切除直达波，而是设定一个平滑窗口函数（即T(x)），并根据该窗口函数切除直达波，其可更有效地消除直达波的高能区域，使得直达波与反射波能量之间有平缓的过渡，从而能够突出反射波的能量，即增强了反射波的信噪比。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：

图1为本发明实施例的全波形反演中去除波动方程模拟直达波的方法的方法流程图；

图2为本发明实施例的全波形反演中去除波动方程模拟直达波的装置的结构原理示意图；

图3为本发明实施例中在去除波动方程模拟直达波之前的模拟的地震记录示意图；

图4为本发明实施例中在去除波动方程模拟直达波之后的模拟的地震记录示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。

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

请参阅图1所示，本发明实施例的一种全波形反演中去除波动方程模拟直达波的方法，包括以下步骤：

步骤S1，输入二维的速度模型v(x,z)，其中，x为水平传播方向，z为深度传播方向。