[发明专利]一种地震纵、横波波场分离与去噪的方法

权利要求书

1.一种地震纵、横波波场分离与去噪的方法，其特征在于，所述方法包括：

根据含有M₀道，每道有N₀个采样点的地震数据生成二维数组X₁；

依据纵波叠加速度和P-P反射波时距规律，将所述二维数组X₁中的地震数据进行正常时差校正，生成校正数据L₁；

对所述校正数据L₁进行奇异值分解，生成初始奇异值矩阵∑₁；

根据所述初始奇异值矩阵∑₁的奇异值谱获取所述初始奇异值矩阵∑₁的畸变位置；

将所述畸变位置前的奇异值构成低通滤波奇异值矩阵∑₂，将所述畸变位置后的奇异值构成高通滤波奇异值矩阵∑₃；

奇异值分解所述高通滤波奇异值矩阵∑₃，生成高通滤波重构信号L_Hp；

依据纵波叠加速度和P-P反射波时距规律，将所述高通滤波重构信号L_Hp进行正常时差反校正处理，生成高通反校正数据X₃；

根据所述M₀、N₀的数值判断是否需要对所述高通反校正数据X₃进行扩边处理，如果是，对所述高通反校正数据X₃进行扩边处理，生成扩边矩阵X₄；

对所述扩边矩阵X₄进行二维傅里叶正变换得到f-k谱Y₁；

根据所述f-k谱Y₁提取目标信号的f-k谱，生成目标f-k谱Y₂；

对所述目标f-k谱Y₂进行二维傅里叶反变换，生成反变换数据X₅；

对所述反变换数据X₅进行缩边处理，生成缩边后数据X₆；

将所述缩边后数据X₆转换为与所述采样点的地震数据相同的格式并输出。

2.根据权利要求1所述的地震纵、横波波场分离与去噪的方法，其特征在于，在将所述畸变位置前的奇异值构成低通滤波奇异值矩阵∑₂之后，所述方法还包括：

奇异值分解所述低通滤波奇异值矩阵∑₂，生成低通滤波重构信号L_Lp；

依据纵波叠加速度和P-P反射波时距规律，将所述低通滤波重构信号L_Lp进行正常时差反校正处理，生成低通反校正数据X₂；

将所述反校正数据X₂转换为与所述采样点的地震数据相同的格式并输出。

3.根据权利要求1或2所述的地震纵、横波波场分离与去噪的方法，其特征在于，所述根据所述M₀、N₀的数值判断是否需要对X₃进行扩边处理，包括：

当所述M₀、N₀的数值为2的整数幂时，不需要对所述高通反校正数据X₃进行扩边处理；

当所述M₀、N₀的数值不是2的整数幂时，需要对所述高通反校正数据X₃进行扩边处理。

4.根据权利要求3所述的地震纵、横波波场分离与去噪的方法，其特征在于，如果不需要对所述高通反校正数据X₃进行扩边处理，所述方法还包括：

对所述高通反校正数据X₃进行二维傅里叶正变换得到f-k谱Y₁；

根据所述f-k谱Y₁提取目标信号的f-k谱，生成目标f-k谱Y₂；

对所述目标f-k谱Y₂进行二维傅里叶反变换，生成反变换数据X₅；

将所述反变换数据X₅转换为与所述采样点的地震数据相同的格式并输出。

5.根据权利要求4所述的地震纵、横波波场分离与去噪的方法，其特征在于，所述目标信号为P-S转换波或噪声信号。