[发明专利]一种用于减少纵波压制的叠加自相关滤波方法及装置

说明书

本发明公开了一种用于减少纵波压制的叠加自相关滤波方法及装置，方法包括对各道之间的自相关函数进行叠加；通过变分模态分解对修饰后的叠加自相关函数进行分解；分别对纵波频率分量及横波频率分量进行标准化处理，并得到重构叠加自相关函数；得到降噪后的各道的信号。本发明的有益效果是：通过使用变分模态分解对修饰后的叠加自相关函数进行分解，将其分解成纵波频率分量及横波频率分量，最后将纵波频率分量及横波频率分量进行标准化后进行重构得到重构叠加自相关函数，再进行修饰后得到最终的滤波函数，通过使用最终的滤波函数与各道进行卷积得到降噪后的各道的信号，由于本方法中的纵波频率能量得到增强，从而可避免滤波函数对纵波的压制。

技术领域

本发明涉及地震数据处理技术领域，尤其是涉及一种用于减少纵波压制的叠加自相关滤波方法及装置。

背景技术

微震去噪是微震监测采集数据处理前的重要步骤，对水力压裂效果评价有很大的影响(如申请号为CN202010760856.4的中国发明专利)。

目前，在自相关叠加滤波对微震进行去噪时，通常对各道间的自相关函数进行叠加，然后对得到的叠加自相关函数进行开窗截断相关函数两边的信息，同时去除中心点的突变，最后将得到最终的自相关函数与各道进行卷积，得到最终的结果。该方法目的是通过自相关消除各道间的时差影响，通过自相关突出自身的主要频率，再通过叠加，得到道间的相似频率，从而自适应的得到滤波器。但是由于接收到的纵横波的频率和能量差异(纵波频率高于横波频率，纵波能量远低于横波能量)，其自相关会对能量较弱的纵波频率部分进行压制，与此同时进行卷积滤波时会对能量较弱的纵波进一步压制。

在微震监测中，信号的初至是比较重要的信息，尤其是纵波初至，自相关叠加滤波由于会对纵波进行压制，这使得该方法并不能完全应用于微震中。

发明内容

有鉴于此，有必要提供一种用于减少纵波压制的叠加自相关滤波方法及装置，用以解决现有的自相关叠加滤波方法对微震进行去噪时会对纵波频率部分进行压制、导致去噪后的信号丢失重要信息的技术问题。

为了实现上述目的，本发明提供了一种用于减少纵波压制的叠加自相关滤波方法，包括：

对微震采集数据中各道之间的自相关函数进行叠加，获取叠加自相关函数；

对叠加自相关函数进行修饰；

通过变分模态分解对修饰后的叠加自相关函数进行分解，得到纵波频率分量及横波频率分量；

分别对纵波频率分量及横波频率分量进行标准化处理，并通过标准化处理后的纵波频率分量及横波频率分量，得到重构叠加自相关函数；

对重构叠加自相关函数进行修饰；

将修饰后的重构叠加自相关函数与各道进行卷积，得到降噪后的各道的信号。

在一些实施例中，对微震采集数据中各道之间的自相关函数进行叠加，获取叠加自相关函数；

其中，所述叠加自相关函数的具体计算公式如下：

其中，N是道数，acf_i是第i道的自相关函数，f是叠加自相关函数，i是大于0的自然数。

在一些实施例中，对叠加自相关函数进行修饰，具体为：通过三角截断的权值公式对叠加自相关函数进行修饰；

其中，三角截断的权值公式为：

其中，w是三角截断的权值公式，d是三角截断的边界阈值，t是采样点；

通过三角截断的权值公式对叠加自相关函数进行修饰的具体公式如下：

f₁＝f·w