[发明专利]一种数据驱动的面波预测方法、装置及电子设备和系统

说明书

本发明提供了一种数据驱动的面波预测方法、装置、计算机可读存储介质和系统。该方法构建十字排列的炮‑检地震观测系统，其中，由若干竖向排列的炮点组成的震源线与由若干横向排列的检波点组成的检波线呈十字相交；获取十字排列的炮‑检地震观测系统的地震数据，并将其由时间域转换为频率域，确定震源线上的参与面波预测的炮点个数以及检波线上参与面波预测的检波点的个数，将其组成若干炮检对；针对每一组炮检对，计算频率域的地震数据中的每一个频率成分下的面波预测结果；将所有炮检对的所有频率成分的面波预测结果从频率域转换为时间域，以获得该炮‑检地震观测系统的面波预测结果。本发明不依赖于近地表模型，无需人工交互，可广泛应用于地面地震资料处理。

技术领域

本发明涉及地面地震数据处理技术领域，尤其涉及一种数据驱动的面波预测方法、装置、计算机可读存储介质和系统。

背景技术

近年来，随着高精度地震勘探的需求不断发展，对地面地震数据处理的要求也在不断提升。瑞雷面波是地面地震资料中存在的主要干扰，面波在勘探地震中能提供的有用信息一般很少，因为它屏蔽了其它更加有用的体波波至，因此被认为是噪声。常用的面波分离方法有F-K滤波法、Radon变换法、频散曲线反演法等，这些方法需要工人干预，耗费大量人力。

发明内容

针对上述问题，本发明的实施例提供了一种数据驱动的面波预测方法、装置、计算机可读存储介质和系统。

第一方面，本发明实施例提供了一种数据驱动的面波预测方法，所述方法包括：

构建十字排列的炮-检地震观测系统，其中，由若干竖向排列的炮点组成的震源线与由若干横向排列的检波点组成的检波线呈十字相交；

获取十字排列的炮-检地震观测系统的地震数据，并将其由时间域转换为频率域，

确定震源线上的参与面波预测的炮点个数以及检波线上参与面波预测的检波点的个数，将其组成若干炮检对；

针对每一组炮检对，计算频率域的地震数据中的每一个频率成分下的面波预测结果；

将所有炮检对的所有频率成分的面波预测结果从频率域转换为时间域，以获得该炮-检地震观测系统的面波预测结果。

根据本发明的实施例，上述针对每一组炮检对，计算频率域的地震数据中的每一个频率成分下的面波预测结果的步骤，包括：

针对每一组炮检对，估算炮点的位置以及检波器的位置；

基于炮点的位置以及检波器的位置，计算频率域的地震数据中的每一个频率成分下的面波预测结果。

根据本发明的实施例，基于炮点的位置以及接收器的位置，根据下式计算频率域的地震数据中的每一个频率成分下的面波预测结果：

式中，G是格林函数，表示由后一坐标位置处的震源引起的前一坐标位置处的压力，ω表示角频率，ρ表示震源位置处的密度，k为常数，炮点位置为x₁,检波器位置为x₂，两个面S′和S分别为检波器边界和震源边界，分别由检波线上参与面波预测的检波点的个数和震源线上的参与面波预测的炮点个数决定。

根据本发明的实施例，上述检波线上参与面波预测的检波点的个数为3至7个。

根据本发明的实施例，上述震源线上的参与面波预测的炮点个数为3至7个。

根据本发明的实施例，上述检波线上参与面波预测的检波点的个数和震源线上的参与面波预测的炮点个数不同。

根据本发明的实施例，在十字排列的炮-检地震观测系统中，所述检波线上的检波点分别位于震源线的两侧。

第二方面，本发明实施例提供了一种数据驱动的面波预测装置，包括：