[发明专利]近地表层间多次波预测方法及装置

权利要求书

1.一种近地表层间多次波预测方法，其特征位于，包括：

采集近地表研究区的二维叠前地震数据；

对于选定的每一个目标地震道，从所述二维叠前地震数据中，提取目标地震道对应的共炮点道集在接收点侧预测孔径内和共接收点道集在炮点侧预测孔径内的所有地震道；

对于炮点和接收点分别位于炮点侧预测孔径和接收点侧预测孔径内的每一炮检对，从对应共炮点道集和共接收点道集中分别提取出一个地震道；

对于炮点和接收点分别位于炮点侧预测孔径和接收点侧预测孔径内的每一炮检对，在层间多次波界面上给定的搜索孔径内搜索最小履行时的向下反射点，并确定该炮检对的格林函数；

对于炮点和接收点分别位于炮点侧预测孔径和接收点侧预测孔径内的每一炮检对，根据每一炮检对的格林函数和提取的两个地震道，预测出该炮检对的层间多次波；

将炮点和接收点分别位于炮点侧预测孔径和接收点侧预测孔径内的所有炮检对预测的层间多次波相加，得到预测的目标地震道的层间多次波。

2.如权利要求1所述的方法，其特征在于，所述方法还包括：

获取层间多次波向下反射点的横坐标、深度坐标和等效速度；

对于每一炮检对，根据层间多次波向下反射点的横坐标、深度坐标和等效速度，确定层间多次波界面的平均深度和等效速度。

3.如权利要求2所述的方法，其特征在于，通过如下公式分别计算层间多次波界面的平均深度和等效速度：

其中，z_av表示层间多次波界面的平均深度；v_av表示层间多次波界面的等效速度；x_S表示共炮点的横坐标；x_R表示共接收点的横坐标；z_q表示层间多次波向下反射点的深度坐标；x_q表示层间多次波向下反射点的横坐标；v_q表示层间多次波向下反射点的等效速度。

4.如权利要求2所述的方法，其特征在于，对于选定的每一个目标地震道，从所述二维叠前地震数据中，提取目标地震道对应的共炮点道集在接收点侧预测孔径内和共接收点道集在炮点侧预测孔径内的所有地震道，包括：

对于每一个目标地震道，获取给定的炮点侧预测孔径和接收点侧预测孔径；

根据给定的接收点侧预测孔径、层间多次波界面的平均深度和等效速度，确定目标地震道对应共炮点道集的切除函数；根据给定的炮点侧预测孔径、层间多次波界面的平均深度和等效速度，确定目标地震道对应共接收点道集的切除函数；

根据目标地震道对应共炮点道集的切除函数，从所述二维叠前地震数据中，提取目标地震道对应共炮点道集在接收点侧预测孔径内的地震道；根据目标地震道对应共接收点道集的切除函数，从所述二维叠前地震数据中，提取目标地震道对应共接收点道集在炮点侧预测孔径内的地震道。

5.如权利要求4所述的方法，其特征在于，通过如下公式提取目标地震道对应共炮点道集在接收点侧预测孔径内的地震道和共接收点道集在炮点侧预测孔径内的地震道：

其中，d_csp表示提取的共炮点道集；d_crp表示提取的共接收点道集；C_S表示共炮点道集的切除函数；C_R表示共接收点道集的切除函数；z_av表示层间多次波界面的平均深度；v_av表示层间多次波界面的等效速度；x_S表示共炮点的横坐标；x_R表示共接收点的横坐标；x_P4表示接收点侧预测孔径内接收点P4的横坐标；x_P3表示炮点侧预测孔径内炮点P3的横坐标。