[发明专利]一种地震数据重建方法

权利要求书

1.一种地震数据重建方法，其特征在于，包括以下步骤：

步骤A：获取待重建的地震数据及其采集矩阵；

步骤B：对地震数据和采集矩阵进行分块处理；

步骤C：建立每个小块的地震数据及其对应的采集矩阵的双曲正切函数正则化的稀疏反演模型，其中，所述双曲正切函数正则化的稀疏反演模型为：

式中，min表示最小化符号，J_σ(x)为目标函数，F^-1为傅里叶变换的逆，x为地震数据在傅里叶变换域的系数，:＝表示将minJ_σ(x)定义为s.t.表示前面目标函数最小化的约束条件为SF^-1x＝d，其中σ为超参数，d为待重建的地震数据，S为采集矩阵，表示地震数据在傅里叶变换域系数的第i个分量的一种双曲正切函数；N为地震数据在傅里叶变换域的系数x的分量个数，x_i表示地震数据在傅里叶变换域的系数x的第i个分量；

步骤D：通过每个小块数据求取该稀疏反演模型的最优解，得到重建后的地震数据；

步骤D21：设定循环次数IterN、初始超参数σ₀和终止迭代超参数σ_IterN，并将观测得到的地震数据作为初始解d₀；

步骤D22：计算当前第k次迭代的迭代解在变换域的系数x_k，其表达式为：

x_k＝Fd_k (2)

式中，k表示迭代序号；F表示傅里叶变换；当k＝0时，表示初始迭代，d_k表示第k次迭代时重建所得的地震数据；

计算当前第k次迭代的修正牛顿方向g_k，其表达式为：

式中，σ_k为第k次迭代的超参数；

求出更新解：

式中，λ为步长；

步骤D23：更新解在时间空间域的投影：

进行k的自增1运算，如果k＜IterN，重复所述步骤D22到所述步骤D23的迭代过程；如果k≥IterN，则迭代结束；

步骤D24：输出小块数据的最后一次迭代时重建所得的地震数据d_IterN；

步骤E：对所有分块数据进行数据重建，并将所有重建的数据拼接成完整的重建数据。