[发明专利]一种基于稀疏反演的滑动扫描谐波压制方法

权利要求书

1.一种基于稀疏反演的滑动扫描谐波压制方法，特点是通过以下具体步骤实现：

1)可控震源采集地震数据，滑动扫描相关前数据和力信号；

2)以一道为例，从相关后记录中取绝对值最大x_max，从序列中搜寻|x₀(t)|＞αx_max，形成以x_max为中心的一段区域，确定可控震源的强初至；

其中：x_max为绝对值最大，α为＜1的常数；

3)利用从力信号中分离出的基波分量进行自相关得到零相位子波，利用该子波在相关后记录中进行匹配追踪，得到初至的稀疏化反射系数，得到相关前记录强初至的基波分量；

4)从相关前记录中减去强初至的基波分量得到剩余记录将其与力信号的1阶谐波s₁(t)进行互相关，得到相关后记录x₁(t)；

5)对相关后记录x₁(t)用1阶谐波自相关得到的子波w₁(t)在强初至区域内采用匹配追踪方法，得到一阶谐波的初至系数，从而得到强初至在相关前记录中的1阶谐波分量；

6)从步骤4)中得到的剩余记录中减去步骤5)得到的强初至一阶谐波分量，得到另一个剩余记录，按照从步骤4)到步骤5)提取一阶谐波分量的方法，提取出强初至在相关前记录中的2阶、3阶或3阶以上高阶谐波分量；

7)采用步骤3)得到的相关前记录强初至的基波分量和步骤5)和步骤6)得到的1阶、2阶或3阶以上高阶谐波分量，按照采集的滑动时间，从前炮中在相应位置减去这些谐波分量，得到没有谐波污染的相关前记录；

8)把步骤7)得到的没有谐波污染的相关前记录和从前炮力信号中提取出的基波分量进行互相关，就得到前炮压制谐波后的相关后正常单炮记录。

2.根据权利要求1所述的方法，步骤1)所述的滑动扫描相关前数据和力信号是首先从最后一炮开始处理，对最后一炮的力信号利用短时傅立叶变换在时频域进行分离，分离出基波、二次谐波、三次谐波成分，将分离出的基波分量与相关前记录作互相关，得到相关后记录x₀(t)。