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

权利要求书

1.一种用于减少纵波压制的叠加自相关滤波方法，其特征在于，包括：

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

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

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

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

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

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

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

其中，VMD(f₁, 2)是变分模态分解将f₁分解为两个模态分量的函数，f₁为修饰后的叠加自相关函数，f_s为横波频率分量，f_p为纵波频率分量；

重构叠加自相关函数的计算公式如下：

其中，f₂为重构叠加自相关函数，f_s为横波频率分量，f_p为纵波频率分量，std为标准化算法；

对重构叠加自相关函数进行修饰，具体包括：

对重构叠加自相关函数进行压制中值突变以及镜像平均处理；

对经压制中值突变以及镜像平均处理后的重构叠加自相关函数的两边的冗余信息进行消除，得到最终的滤波函数；

对重构叠加自相关函数进行压制中值突变以及镜像平均处理的具体公式如下：

其中，f_m是经压制中值突变以及镜像平均处理后的重构叠加自相关函数，f₂是重构叠加自相关函数，i是自然数；

最终的滤波函数的计算公式如下：

其中，w是三角截断的权值公式，d是三角截断的边界阈值，t为采样点，f_m是经压制中值突变以及镜像平均处理后的重构叠加自相关函数，f_new是最终的滤波函数。

2.根据权利要求1所述的用于减少纵波压制的叠加自相关滤波方法，其特征在于，对微震采集数据中各道之间的自相关函数进行叠加，获取叠加自相关函数；

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

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

3.根据权利要求1所述的用于减少纵波压制的叠加自相关滤波方法，其特征在于，对叠加自相关函数进行修饰，具体为：通过三角截断的权值公式对叠加自相关函数进行修饰；

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

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

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

其中，w是三角截断的权值公式，f是叠加自相关函数，f₁是修饰后的叠加自相关函数。