[发明专利]地震数据的处理方法和装置

权利要求书

1.一种地震数据的处理方法，其特征在于，包括：

将目标地震数据表示为高斯束叠加形式的地震数据，其中，所述高斯束叠加形式的地震数据中包括所述目标地震数据的波形信息；

将所述高斯束叠加形式的地震数据进行稀疏处理，处理得到波形函数矩阵；

基于处理之后得到的所述波形函数矩阵进行地震数据的重构，得到重构之后的地震数据；

其中，将所述高斯束叠加形式的地震数据进行稀疏处理，处理得到波形函数矩阵包括：

基于波形信息矩阵建立稀疏矩阵优化模型，其中，所述波形信息矩阵中包括所述波形信息；

对所述稀疏矩阵优化模型中的波形信息矩阵进行稀疏处理，得到处理之后的所述稀疏矩阵优化模型，以使稀疏处理之后的所述波形信息矩阵中的每个数值均不等于预设数值；

对处理之后的所述稀疏矩阵优化模型进行稀疏求解，求解得到所述波形函数矩阵；

其中，对处理之后的所述稀疏矩阵优化模型进行稀疏求解，求解得到所述波形函数矩阵包括：

反复通过第一公式和第二公式对稀疏矩阵优化模型进行迭代计算，求解得到所述波形函数矩阵，直至满足第三公式时，确定所述第二公式中的输出结果为所述波形函数矩阵，当ρ_k≤∈时，终止迭代，输出波形矩阵∈为大于零的常数，即为上述波形函数矩阵，为波形函数矩阵的第k行向量；

其中，α₁，α₂和α_K为对角线矩阵Λ中的数值，G_k为迭代变量，p为稀疏度参数，k为迭代指数，迭代指数k的初值等于0，其中，k＝0时，对应着为矩阵算子，D为数据矩阵，为稀疏处理之后的所述波形信息矩阵，I_m表示m维恒等算子，为矩阵的第k行向量，

2.根据权利要求1所述的方法，其特征在于，对所述稀疏矩阵优化模型中的波形信息矩阵进行稀疏处理包括：

对所述目标地震数据中的高斯束束中心相邻道上进行相关性度量分析，得到地震子波，其中，所述地震子波中包含所述地震子波所处的时间和倾角方向；

根据所述地震子波所处的时间和所述倾角方向确定所述波形信息矩阵中数值等于所述预设数值的向量；

在所述波形信息矩阵中删除数值等于所述预设数值的向量，得到处理之后的所述稀疏矩阵优化模型。

3.根据权利要求1所述的方法，其特征在于，基于处理之后得到的所述波形函数矩阵进行地震数据的重构，得到重构之后的地震数据包括：

通过公式d(x,t)＝∑_γc_γψ_γ(x,t)进行地震数据的重构，得到重构之后的地震数据，其中，为所述波形函数矩阵，d(x,t)为所述重构之后的地震数据，ψ_γ(x,t)为基函数ψ(x,t；x_m,t_k,p_x)的简写，γ表示参数(x_m,t_k,p_x)的集合。

4.根据权利要求1所述的方法，其特征在于，将目标地震数据表示为高斯束叠加形式的地震数据包括：

获取高斯束基函数的半宽w和所述高斯束中任意相邻的两个束中心之间的束中心间距a；

根据所述半宽和所述束中心间距计算所述高斯束的束中心位置x_m，其中，x_m＝am，m＝1,2,…,M，M为束中心的个数；

通过公式确定所述目标地震数据中每道地震数据与所述束中心相对应的基函数，其中，t_k为所述高斯束基函数的束中心所处的时间位置，p_x为水平慢度，ω_r为相关低频；

通过所述目标地震数据中每道地震数据与所述束中心相对应的基函数，将所述目标地震数据表示为所述高斯束叠加形式的地震数据。