[发明专利]多次波干扰的压制方法、存储介质和计算机设备

权利要求书

1.一种地震数据中多次波干扰的压制方法，其特征在于，包括以下步骤：

获取原始地震数据，对原始地震数据进行预处理，得到共反射点道集数据；

利用逆散射级数法对共反射点道集数据进行多次波预测，得到层间多次波模型数据；

将共反射点道集数据和层间多次波模型数据变换到F-K域，在F-K域，利用层间多次波模型数据，设计自适应滤波器，利用自适应滤波器对共反射点道集数据进行滤波，以将共反射点道集数据分离成含有部分一次波的多次波数据和未包含多次波的一次波数据，并通过F-K反变换将层间多次波模型数据、含有部分一次波的多次波数据和未包含多次波的一次波数据变换到时间空间域；

利用时间空间域的层间多次波模型数据，对时间空间域的含有部分一次波的多次波数据执行曲波域匹配相减，以压制含有部分一次波的多次波数据中的多次波数据，从而得到剩余的一次波数据；

将剩余的一次波数据与未被多次波干扰的一次波数据进行相加融合，得到压制层间多次波干扰后的地震数据。

2.根据权利要求1所述的地震数据中多次波干扰的压制方法，其特征在于，所述预处理包括静校正、去噪、反褶积、速度分析和克希霍夫叠前时间偏移。

3.根据权利要求1所述的地震数据中多次波干扰的压制方法，其特征在于，利用逆散射级数法对共反射点道集数据进行多次波预测，得到层间多次波模型数据，包括：

根据共反射点道集数据，利用逆散射级数法确定层间多次波的预测算子；

根据层间多次波的预测算子，确定层间多次波模型数据。

4.根据权利要求3所述的地震数据中多次波干扰的压制方法，其特征在于，利用逆散射级数法确定层间多次波的预测算子，包括：

利用下式确定层间多次波的预测算子：

其中，b₁为单频平面波场，b₁(k_g,k_s,q_g+q_s)＝-2iq_sA(ω)^-1D(k_g,k_s,ω)，b₃为层间多次波的预测算子，D(k_g，k_s，ω)为共反射点道集数据，ε_s、q_s、k₁和k_s分别为震源的深度、垂向波数、水平波数积分变量和水平波数，ε_g、q_g、k₂和k_g分别为检波器的深度、垂向波数、水平波数积分变量和水平波数，λ为条件参数，q₁和q₂分别表示震源和检波器的垂向波数积分变量，ω为角频率，A(ω)为地震子波，i为虚数单位，z₁、z₂和z₃分别是空间上从左到右3个散射点的深度。

5.根据权利要求4所述的地震数据中多次波干扰的压制方法，其特征在于，根据层间多次波的预测算子，确定层间多次波模型数据，包括：

利用下式确定层间多次波模型数据：

M(k_g,k_s,ω)＝(-2iq_s)^-1A(ω)b₃(k_g,k_s,q_g+q_s)

其中，M(k_g，k_s，ω)表示层间多次波模型数据，q_s和k_s分别为震源的垂向波数和水平波数，q_g和k_g分别为检波器的垂向波数和水平波数，ω为角频率，i为虚数单位，A(ω)为地震子波，b₃为层间多次波的预测算子。