[发明专利]核范数与广义全变差联合约束的地震随机噪声压制方法

权利要求书

1.核范数与广义全变差联合约束的地震随机噪声压制方法，其特征在于，包括以下步骤：

(1)输入三维地震数据体，将地震数据按Inline线进行输入，再对每条Inline线地震数据按时间采样点与CMP号排列成原始含躁地震数据矩阵D，其中时间采样点按行排列，CMP按列排列；

(2)构建目标函数；

构建核范数与广义全变差联合约束的地震随机噪声压制方法的目标函数，即

其中，F(X)即为目标函数，D为原始含躁地震数据矩阵，X为噪声压制后的地震数据矩阵，为矩阵X-D的Frobenius范数，||X||_*为X的核范数，||X||_GTV为X的广义全变差，λ与μ分别为核范数约束系数和广义全变差约束系数；

(3)采用分裂Bergman迭代策略，将目标函数分解为两个子问题进行迭代，针对具体子问题采用相应算法进行求解；

(4)将噪声压制后的地震数据矩阵X视为已知量，输入核范数约束系数λ与辅助参数β，采用奇异值软阈值化算法求解子问题一，求取辅助变量矩阵Θ；

(5)输入广义全变差约束系数μ以及上一次Bergman迭代得到的地震数据矩阵，以及步骤四得到的辅助变量矩阵Θ，采用梯度投影算法求解子问题二，并求取更新后的地震数据矩阵；

(6)判断Bergman迭代是否停止，若停止迭代，则输出随机噪声压制后的地震数据矩阵；

(7)对每条Inline线地震数据都进行随机噪声压制处理，得到整个随机噪声压制后的地震数据体。