[发明专利]基于波场分解的绕射多次波压制方法

权利要求书

1.一种基于波场分解的绕射多次波压制方法，其特征在于，包括以下步骤：

a、将原始地震数据排列成共炮点记录的数据矩阵；

b、对各个共炮点记录进行线性拉冬变换；

c、以反射系数为准重新排列，获得不同角度入射的共p剖面平面波数据，其中在共p剖面中，反射波呈现近似线性的连续，而绕射波则表现为近似于双曲线形状；

d、依选择合适大小的时窗，将抽取剖面的地震数据通过快速傅里叶变换将数据映射到f-x域，得到其f-x域信号的实部和虚部；

e、对于每一个频率切片；

e1、对于实部数据进行EMD分解，抽取合适的IMF分量，为前两个IMF分量，重构作为绕射波场实部，剩余部分作为反射波场实部，其中EMD分解即频率域经验模态分解，IMF分量可通过计算信号与局部极大值点三次样条连线和局部极小值点连线构成的上下包络线平均值之差获取；

e2、对于虚部数据进行EMD分解，抽取合适的IMF分量，为前两个IMF分量，重构作为绕射波场虚部，剩余部分作为反射波场虚部；

e3、将分解好的实部和虚部数据重构成复数数据；

f、将重构的数据通过反傅里叶变换到t-x域，生成绕射波场和反射波场两组数据；

g、滑动时窗，对下一个时窗重复以上操作，直到全部数据完成；

h、将共p剖面的已分解的反射波场和绕射波场重新排列成平面波数据；

i、做反拉冬变换，将平面波数据变成共炮点道集形式的反射波场P_r和绕射波场数据P_s；

j、将分解后的反射波场P_r和绕射波场数据P_s进行分级多次波预测，获得预测模型；

展开反射模态多次波、绕射模态多次波、混合模态多次波模型；

k、对比绕射多次波与四种多次波预测模型的曲率，从中选择相似度最高模型进行预测减去。