[发明专利]一种气泡压制的方法和装置

说明书

本发明提供了一种气泡压制的方法和装置，其中，该方法包括：根据待测工区中的多个单炮地震数据，确定待测工区中的最小相位子波；基于实际勘探过程中气泡产生时间，对最小相位子波进行切除处理，得到切除处理后的最小相位子波；根据切除处理后的最小相位子波，计算得到待测工区中的气泡压制算子；利用气泡压制算子对每道地震数据进行卷积处理，得到气泡压制后的地震数据。在本发明实施例中，通过实际数据来求取子波，充分考虑了实际因素对子波的影响，能有效提高气泡的压制程度，消除气泡效应，提高数据的信噪比，从而提高海上地质勘探的分辨率。

技术领域

本发明涉及地质勘探技术领域，特别涉及一种气泡压制的方法和装置。

背景技术

海上地震勘探时，一般可以采用气枪作为震源。当气枪震源激发时，高压空气释放到水中，迅速形成球形气泡，球形气泡在水中不断膨胀压缩，形成多个气泡脉冲，该过程称为气泡效应。在气泡效应的过程中，气泡膨胀会产生地震波，由于多次气泡胀缩，气枪子波会产生多次振荡，大大影响了所采集的地震数据的精度，降低了地震数据的分辨率。因此，有必要对海上地震数据进行气泡压制处理，从而达到去除海上地震数据中由于气泡效应所产生的数据噪声的目的。

目前，一般可以基于气枪震源所产生的远场子波，来确定海上地震勘探时所采集到的数据进行气泡压制后的效果。然而，由于远场子波是采用理论方法得到的模拟子波，模拟远场子波时要作很多假设，例如：海面是水平的，海水是静止的等等。这些假设条件与实际勘探情况相差甚远，使得模拟的远场子波与实际勘探中产生的子波存在明显区别，这就导致采用远场子波的气泡压制效果很不理想，有时很难满足实际生产的需要。

发明内容

本发明提供了一种气泡压制的方法和装置，以达到获取与实际情况吻合度较高的海上气泡压制效果的目的。

本发明实施例提供了一种气泡压制的方法，可以包括：根据待测工区中的多个单炮地震数据，确定所述待测工区中的最小相位子波；基于实际勘探过程中气泡产生时间，对所述最小相位子波进行切除处理，得到切除处理后的最小相位子波；根据所述切除处理后的最小相位子波，计算得到待测工区中的气泡压制算子；利用所述气泡压制算子对所述多个单炮地震数据中各个单炮地震数据的每道地震数据进行卷积处理，得到气泡压制后的地震数据。

在一个实施例中，确定所述待测工区中的最小相位子波，可以包括：对所述多个单炮地震数据中满足预设要求的单炮地震数据进行自相关处理；将多个单炮地震数据中的各道地震数据的自相关结果进行叠加，得到一道自相关叠加数据；计算所述一道自相关叠加数据对应的零相位子波；基于所述一道自相关叠加数据对应的零相位子波，计算得到最小相位子波。

在一个实施例中，基于所述一道自相关叠加数据的零相位子波，计算得到最小相位子波，可以包括：以所述零相位子波作为输入，脉冲作为期望，计算得到所述零相位子波的反子波；以所述零相位子波的反子波作为输入，脉冲作为期望，计算得到所述零相位子波对应的最小相位子波。

在一个实施例中，所述实际勘探过程中气泡产生时间的范围可以为：大于等于40ms，并且小于等于120ms。

在一个实施例中，可以按照以下公式对所述最小相位子波进行切除处理，得到切除处理后的最小相位子波：

c(n)＝a(n)h(n)

上式中，

上式中，c(n)表示所述切除处理后的最小相位子波，a(n)表示所述最小相位子波，h(n)表示切除函数，n表示采样点数，n₀表示气泡产生时间对应的采样点个数。

在一个实施例中，所述切除函数还可以包括：汉明窗函数、钟形窗函数中的一种。

在一个实施例中，可以按照以下公式根据所述切除处理后的最小相位子波，计算得到待测工区中的气泡压制算子：