[发明专利]基于插值算子的非规则网格混叠地震数据分离方法及系统

说明书

本发明涉及一种基于插值算子的非规则网格混叠地震数据分离方法及系统，该方法包括：S1、获取非规则网格混叠地震数据d_bl、混叠算子Γ、规则网格位置以及非规则网格位置；S2、基于规则网格位置以及非规则网格位置设计基于Sinc插值的插值算子W及共轭算子W^T；S3、循环迭代更新获取规则网格上混叠分离数据。与现有技术相比，本发明具有分离效率高、精度高等优点。

技术领域

本发明涉及一种非规则网格混叠地震数据分离方法及系统，尤其是涉及一种基于插值算子的非规则网格混叠地震数据分离方法及系统。

背景技术

常规地震数据采集方法的采集周期长、采集成本高，为了提高采集效率，地震数据混叠采集技术应运而生。但是混叠噪声的存在影响后续传统地震数据偏移、成像的精度，因此混叠数据分离成为必要步骤。常用的混叠数据分离方法包括滤波类方法和反演类方法等，但是传统的混叠分离方法均针对规则网格混叠数据。而实际混叠数据采集过程中炮点一般布置在非规则网格上，通常利用Binning策略将非规则网格数据转化为规则网格数据，但是降低了后续的混叠分离精度。另外，非均匀傅里叶变换也被引入非规则网格地震数据处理中，但是计算效率较低。因此，如何提高非规则网格混叠地震数据的分离效率和精度成为亟待解决的技术问题。

发明内容

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于插值算子的非规则网格混叠地震数据分离方法及系统。

本发明的目的可以通过以下技术方案来实现：

一种基于插值算子的非规则网格混叠地震数据分离方法，该方法包括：

S1、获取非规则网格混叠地震数据d_bl、混叠算子Γ、规则网格位置以及非规则网格位置；

S2、基于规则网格位置以及非规则网格位置设计基于Sinc插值的插值算子W及共轭算子W^T；

S3、采用下式循环迭代更新获取规则网格上混叠分离数据直至迭代结束：

d^k+1＝C^TT_λC(d^k-W^TΓ^T(ΓWd^k-d_bl))

其中，d^k为第k次迭代时规则网格上混叠分离数据，W^T为插值算子W的共轭算子，Γ^T为混叠算子Γ的共轭算子，C为稀疏Curvelet变换矩阵，C^T为逆稀疏Curvelet变换矩阵，T_λ为以λ为阈值的硬阈值函数，k＝1，2，……，N，N为迭代总次数，k＝1时，d¹＝0。

优选地，步骤S2具体为：

对于非规则网格点r_i，确定以r_i为中心、R为半径的邻域Nr_i；

对于邻域Nr_i内的任意一个规则网格点j，计算规则网格点j对非规则网格点r_i的权重W_j，j＝1，2，……NUM，NUM为邻域Nr_i内规则网格点个数；

将邻域Nr_i内所有规则网格点的权重集合作为插值算子，表示为W＝[W₁，W₂……W_NUM]，同时求取W的共轭算子W^T。

优选地，W_j通过下式获得：