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

权利要求书

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

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。

2.根据权利要求1所述的一种基于插值算子的非规则网格混叠地震数据分离方法，其特征在于，步骤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。

3.根据权利要求2所述的一种基于插值算子的非规则网格混叠地震数据分离方法，其特征在于，W_j通过下式获得：

其中，x为规则网格点j到非规则网格点r_i的距离，I₀(·)为第一类改进的零阶贝塞尔函数，a为超参数。

4.一种基于插值算子的非规则网格混叠地震数据分离系统，其特征在于，该系统包括：

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

插值算子确定模块：基于规则网格位置以及非规则网格位置设计基于Sinc插值的插值算子W及共轭算子W^T；

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

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。