[发明专利]一种从地震记录中去除面波的方法和装置

权利要求书

1.一种从地震记录中去除面波的方法，其特征在于，该方法包括：

S1：对待处理的地震记录进行分频扫描，获得面波频段的分布范围；

S2：根据所述面波频段的分布范围，对所述地震记录进行小波分频，获取不同频段上的地震记录，所述不同频段上的地震记录包括包含面波的地震记录和不包含面波段的地震记录；

S3：根据所述面波频段的分布范围，从所述不同频段上的地震记录中识别出包含面波的地震记录，对所述包含面波的地震记录进行径向道变换；

S4：根据所述面波频段的分布范围确定高通滤波的频率，对S3中径向道变换后的地震记录在径向道域进行高通滤波；

S5：对S4中高通滤波后的地震记录进行反径向道变换；

S6：将S5中反径向道变换后的地震记录和S2中所述的不包含面波的地震记录进行小波重构，得到去除面波的地震记录。

2.根据权利要求1所述的从地震记录中去除面波的方法，其特征在于，所述地震记录记录了多道采样点，每道有多个采样点，所述采样点记录了信号的振幅值。

3.根据权利要求1所述的从地震记录中去除面波的方法，其特征在于，所述对所述地震记录进行分频扫描，获得面波频段的分布范围包括：

对所述地震记录进行频谱分析，获得地震记录的频段范围；

将所述地震记录进行频段划分，获得不同的频段的地震记录；

对所述地震记录在每个频段范围内进行频率滤波，获得不同频段的时间域的地震记录；

根据所述不同频段的时间域的地震记录确定面波，获得面波频段的分布范围。

4.根据权利要求1所述的从地震记录中去除面波的方法，其特征在于，所述根据所述面波频段的分布范围，对所述地震记录进行小波分频，获取不同频段上的地震记录包括：

根据面波频段的分布范围设置所述小波分频的分频范围；

采用马拉特(Mallat)算法对所述地震记录进行分频，把地震记录分解成地震记录小波细节和地震记录小波逼近，并将上一次分解的地震记录小波逼近作为下一次输入，获取不同频段上的地震记录，公式如下：

W2j+1f=S2jf*H]]>

S2j+1f=S2jf*L]]>

上式中，High，Low为一对正交镜像滤波器；j是分解次数，j＝0，1，......j-1，j＝0表示第1次分解；2^j是分析尺度，是第j+2次地震记录小波细节，是第j+1次地震记录小波逼近，是第j+2次地震记录小波逼近。

5.根据权利要求4所述的从地震记录中去除面波的方法，其特征在于，所述对所述包含面波的地震记录进行径向道变换包括：

将t-x域的所述包含面波的地震记录变换到τ-v域，并利用三次样条插值来实现径向道变换；从t-x域的变换到τ-v域的映射关系如下：

t→τ＝t-t₀

x→v＝(x-x₀)/(t-t₀)

上式中，t为双程旅行时；x为偏移距；x₀为径向道变换的原点位置；t₀为径向道变换的原点位置的双程旅行时；τ为延迟时；v为速度。