[发明专利]一种水陆检数据标定的方法

说明书

本发明提供了一种水陆检数据标定的方法，包括：根据时间‑空间域水、陆检数据计算频率‑波数域水、陆检数据，并切除其中的直达波场数据；根据时间采样间隔、空间采样间隔、海水速度计算频率‑波数耦合系数序列，并结合海水密度计算频率‑波数耦合因子序列；基于频率‑波数耦合因子序列、频率‑波数域陆检数据、频率‑波数域水检数据以及频率‑波数域切除函数经耦合自相关/耦合互相关变换处理、直达波场数据切除计算标定频率‑波数域陆检数据及标定时间‑空间域陆检数据；根据时间‑空间域水检数据和标定时间‑空间域陆检数据分别计算上行波场数据及下行波场数据。根据本发明提供了更优的上、下行波场数据，计算精确、省时、快捷，计算效率高。

技术领域

本发明涉及资源勘探技术领域，尤其涉及地震勘探数据处理技术领域，具体涉及水陆检地震数据处理技术领域，尤其适用于海上采集的OBC/OBN水陆检地震数据。

背景技术

海洋OBC/OBN采集的数据，提供了同一位置水检数据和陆检数据两种数据。这两种数据分别使用水中检波器和陆地检波器记录。水中检波器是一种压力检波器，记录的是地震波产生的压力变化，是一个标量；陆地检波器是一种质点速度检波器，记录的是质点速度变化，是一个矢量。由于这两种检波器记录机理不同，对于同一位置处虚反射和海水鸣震多次波干扰，表现出不同特征。与水中检波器记录的海水鸣震和虚反射多次波干扰相比，陆地检波器记录的海水鸣震和虚反射多次波干扰表现出极性和振幅特征差异。两种检波器记录的海水鸣震和虚反射多次波干扰，其振幅不同，极性相反，且相差一个与海底反射系数成比例的常数，这个常数值就是标定因子。基于这一机理，一些文献中分别提出了水陆检数据常数标定方法，例如参见文献ZL201310736942.1“一种水陆检波器地震数据标定因子反演的方法”，或“OBC水陆检数据标定因子估算方法”，石油地球物理勘探，2016,51(1)。在这些方法基础上，另一些文献中提出了在水陆检数据常相位校正基础上使用时间域匹配滤波器方法，实现水陆检数据匹配，例如参见文献ZL201310531139.4“一种水陆检波器地震数据匹配方法”，或“OBC水陆检数据匹配技术”，石油地球物理勘探，2015,50(1)。

常规的水陆检数据常数标定法和时间域匹配滤波法，在实际数据处理中取得了一定的成功。但是由于陆检数据是由陆地检波器记录的，记录的是质点速度变化，是一个矢量，因此不但考虑时间方向变化，还要考虑空间方向变化，即时间-空间(2D)变化。而常规方法仅仅考虑时间方向变化(时间域匹配滤波器)，是1D的，忽略了空间方向变化，使用1D多道方法计算，计算简单，但是精度低。随着油气勘探的发展，对勘探精度要求越来越高，从理论上实现高精度勘探，即2D/3D计算，是一个迫切需求。因为随着时间方向越深，空间偏移距方向越大，实际地震数据能量都在衰减。

因此，针对现有技术中的上述缺点、问题，需要提出一种优化的水陆检数据标定的方法，降低计算成本、提高计算精确性、减小时间消耗，从而提高为后续上下行波场联合反褶积和镜像偏移成像处理提供更优的上行波场和下行波场数据，进而提高整体的效率。

发明内容

有鉴于此，本发明的目的在于提出一种改进的用于标定水陆检数据的方法，从而解决现有技术中存在的上述问题。

基于上述目的，本发明提供了一种水陆检数据标定的方法，其中该方法包括以下步骤：

根据时间-空间域水检数据h[j,i]、时间-空间域陆检数据g[j,i]计算频率-波数域水检数据H[k,m]及频率-波数域陆检数据G[k,m]，并切除其中的直达波场数据；

根据时间采样间隔△t、空间采样间隔△x、海水速度c_w计算频率-波数耦合系数序列d_k,m，并结合海水密度ρ_w计算频率-波数耦合因子序列A[k,m]；