[发明专利]共炮点道集多级面波压制方法和装置

权利要求书

1.一种共炮点道集多级面波压制方法，其特征在于，包括：

在共炮点道集上确定出面波区域和不含面波的反射波区域；

确定所述面波区域中面波的起始频率、截至频率和最低视速度；

在面波的起始频率至截至频率范围内计算所述反射波区域中的每个频率上振幅的多道统计值；

在面波的起始频率至截至频率范围内对所述面波区域的数据，以所述多道统计值为期望振幅谱，进行零相位滤波处理，得到初步压制面波的结果；

利用面波的截至频率、最低视速度、空间采样间隔和炮检距计算线性相移算子；

利用所述线性相移算子对所述初步压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到压制面波的结果；

其中，利用面波的截至频率、最小视速度、空间采样间隔和炮检距计算线性相移算子，包括：

利用面波的截至频率和空间采样间隔计算满足无面波假频的最小视速度；

确定所述满足无面波假频的最小视速度与面波的最小视速度之间的大小关系；

如果所述满足无面波假频的最小视速度小于面波的最小视速度，则确定线性相移系数为零；

如果所述满足无面波假频的最小视速度大于面波的最小视速度，则将所述满足无面波假频的最小视速度的倒数与面波最小视速度的倒数之间的差值作为线性相移系数；

将所述线性相移系数、频率和炮检距三者相乘得到线性相移算子。

2.如权利要求1所述的方法，其特征在于，在面波的起始频率至截至频率范围内计算所述反射波区域中的每个频率上振幅的多道统计值，包括：

将所述反射波区域的数据变换至频率域，得到频率空间域的数据体；

在面波的起始频率至截至频率范围内计算所述频率空间域数据体的每个频率上振幅的多道统计几何平均值和标准差；

将所述几何平均值和所述标准差相加，得到面波的起始频率至截至频率范围内的每个频率上振幅的多道统计值。

3.如权利要求2所述的方法，其特征在于，在面波的起始频率至截至频率范围内计算所述频率空间域数据体的每个频率上振幅的多道统计几何平均值，包括：

对所述频率空间域数据体取自然对数；

在面波的起始频率至截至频率范围内的每个频率上，对取自然对数后的结果计算几何平均值；

对所述几何平均值取做e指数处理，得到在面波的起始频率至截至频率范围内的所述频率空间域数据体的每个频率上振幅的多道统计几何平均值。

4.如权利要求1所述的方法，其特征在于，利用所述线性相移算子将所述初步压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到压制面波的结果，包括：

计算所述初步压制面波结果的相位谱；

对所述初步压制面波结果的相位谱加上所述线性相移算子得到线性相移后的相位谱；

对所述线性相移后的相位谱，取e的复指数运算得到包含相位信息的信号；

对所述包含相位信息的信号进行FKK域的滤波处理，消除相位信息中的面波信息，并做去e的复指数处理，得到相位滤波后的相位谱；

对所述相位滤波后的相位谱减去线性相移因子，得到反线性相移后的相位谱；

对所述反线性相移后的相位谱与所述初步压制面波后的振幅谱结合反变换到时间域得到压制面波的结果。

5.如权利要求1至4中任一项所述的方法，其特征在于，将所述初步压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到压制面波的结果之后，所述方法还包括：

确定得到的所述压制面波后的结果中面波是否被完全压制；

如果没有，则重复执行以下操作，直至面波被完全压制：将所述压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到再次面波压制的结果。