核范数与广义全变差联合约束的地震随机噪声压制方法技术

本发明专利技术提供一种核范数与广义全变差联合约束的地震随机噪声压制方法，包括(1)输入三维地震数据体；(2)构建目标函数；(3)采用分裂Bergman迭代策略，将目标函数分解为两个子问题；(4)采用奇异值软阈值化算法求解子问题一，求取辅助变量矩阵；(5)采用梯度投影算法求解子问题二，求取地震数据矩阵；(6)判断Bergman迭代是否停止，若停止迭代，输出地震数据矩阵；(7)将整个地震数据体计算完毕，得到最终的压噪后地震数据体。本发明专利技术将图像处理中的广义全变差约束、矩阵分析中的核范数约束、最优化理论的分裂Bergman迭代有机结合，适用于地震数据随机噪声压制处理等，有效提高地震资料的信噪比与保真度。

【技术实现步骤摘要】
核范数与广义全变差联合约束的地震随机噪声压制方法
本专利技术属于石油地球物理勘探
，具体技术为核范数与广义全变差联合约束的地震随机噪声压制方法，可以有效地压制地震数据中的随机噪声，提高地震资料信噪比，同时保持地震数据的局部光滑性，提高地震资料保真度。
技术介绍
石油勘探发展到今天，勘探开发逐渐转向复杂型、隐蔽型、深层和非常规油气藏，对地震资料质量要求越来越高。受复杂采集环境影响，野外采集的地震资料常含有很强的随机噪声，严重的甚至将有效信息完全淹没。地震随机噪声压制是地震资料处理中需要面对的重要问题。地震随机噪声压制方法有很多，这些方法通常不直接在时间域内对地震资料进行去噪处理，而是在某种变换域内对其变换系数进行相应处理，之后进行反变换来达到去噪目的。用于地震数据随机噪声压制的变换有很多种，如Fourier变换、Radon变换、Wavelet变换、Curvelet变换、Seislet变换、Ridgelet变换、Shearlet变换等。Fourier变换是一种全局变换，基于Fourier变换的阈值去噪方法可能会产生吉布斯现象；Radon变换虽然能较好地识别直线特征，但地震波前通常并不是直线；Wavelet变换在地震去噪方面表现出一定的优越性，但它主要用于表示各向同性的点奇异特征，对于地震波前这种曲线奇异特征并不能进行有效的稀疏表示，因此基于Wavelet变换的地震数据压缩和去噪等方法，都不可避免地在边界和波前奇异附近存在不光滑和模糊现象。近年来发展起来的Curvelet变换、Seislet变换、Ridgelet变换、Shearlet变换等则弥补了传统小波变换的不足，具有良好的多尺度性、多方向性和各向异性，更加适合表示地震波前特征，该变换一经提出，便在地震数据处理领域获得广泛应用。基于变换域的地震随机噪声压制方法都可统一在信号稀疏性表示理论中，并在地震数据处理领域引起广泛关注。例如Curvelet变换得到广泛应用正是由于其具有很好的稀疏表示能力。但在许多情况下，地震数据的稀疏性未必能有效地表现地震数据的内部结构特征，特别是相邻地震道间的波形相似性。而基于相邻地震道波形相似性的噪声压制方法，如核范数约束的噪声压制等，很少考虑到地震数据的局部光滑性，使得随机噪声压制地震数据在地质体边缘部位产生强烈的“阶梯抖动”现象。
技术实现思路
本专利技术的目的在于提供核范数与广义全变差联合约束的地震数据随机噪声压制方法，综合地震数据处理、图像处理、最优化理论、矩阵分析等多个学科，将图像处理中的广义全变差约束、矩阵分析中的核范数约束、最优化理论的分裂Bergman迭代有机结合，适用于地震数据随机噪声压制处理等多个环节，有效提高地震资料的信噪比与保真度。为达到以上技术目的，本专利技术提供以下技术方案。核范数与广义全变差联合约束的地震随机噪声压制方法，其核心是地震数据矩阵核范数约束、广义全变差约束、分裂Bergman迭代求解。该方法依次包括以下步骤：(1)输入三维地震数据体，将地震数据按Inline线进行输入，再对每条Inline线地震数据按时间采样点与CMP号排列成矩阵，其中时间采样点按行排列，CMP按列排列；(2)构建目标函数；(3)(3)采用分裂Bergman迭代策略，将目标函数分解为两个子问题进行迭代，针对具体子问题采用相应算法进行求解；(4)采用奇异值软阈值化算法求解子问题一，求取辅助变量矩阵；(5)采用梯度投影算法求解子问题二，求取地震数据矩阵；(6)判断Bergman迭代是否停止，若停止迭代则输出随机噪声压制后的地震数据矩阵；(7)对每条Inline线地震数据都进行随机噪声压制处理，得到整个随机噪声压制后的地震数据体。其中，所述的步骤(2)具体过程为：构建核范数与广义全变差联合约束随机噪声压制方法的目标函数；所述的步骤(4)具体过程为：将地震数据矩阵视为已知量，输入核范数约束系数与辅助参数，采用奇异值软阈值化算法求解子问题一，求取辅助变量矩阵；所述的步骤(5)具体过程为：输入广义全变差约束系数、上一次Bergman迭代得到的地震数据矩阵，以及步骤(4)得到的辅助变量矩阵，采用梯度投影算法求解子问题二，求取更新后的地震数据矩阵。本专利技术提供的核范数与广义全变差联合约束的地震随机噪声压制方法，优点在于：有机结合核范数和广义全变差，对地震随机噪声的压制处理进行联合约束；通过多道地震数据矩阵的核范数约束有效地表征相邻地震道之间的波形相似性，通过广义全变差约束约束地震数据的局部光滑性，从而更有效地压制地震数据中的随机噪声，提高地震资料信噪比，并保持地震数据的局部光滑性，提高地震资料保真度。附图说明图1为本专利技术的流程图。具体实施方式结合实施例说明本专利技术的具体技术方案。核范数与广义全变差联合约束的地震随机噪声压制方法，如图1的流程：步骤一：输入含躁地震数据体。将地震数据按Inline线进行输入，再对每条Inline线地震数据按时间采样点与CMP号排列成矩阵D，其中时间采样点按行排列，CMP按列排列。步骤二：构建核范数与广义全变差联合约束的地震随机噪声压制方法的目标函数，即其中，F(X)即为目标函数，D为原始含躁地震数据矩阵，X为噪声压制后地震数据矩阵，为矩阵X-D的Frobenius范数，||X||*为X的核范数，||X||GTV为X的广义全变差，λ与μ分别为相应约束系数。目标函数(1)中，矩阵X的核范数为其奇异值之和，即其中Λ(X)是矩阵X的奇异值组成的矢量，σi是矩阵X的第i个奇异值，||Λ(X)||1为矢量Λ(X)的1范数，而矩阵X的广义全变差为其中为导数算子，为二阶导数算子，vec为矩阵矢量化转换算子。因此，目标函数(1)可等价地写为步骤三：采用分裂Bergman迭代策略，将目标函数(4)分解为两个子问题迭代，针对具体子问题采用相应算法进行求解。分裂Bergman迭代策略的基本思想是引进一个辅助变量矩阵Θ，将目标函数(4)改写为其中F(X,Θ)为改写后的目标函数，β为Bergman迭代辅助参数。由分裂Bergman迭代原理可得，随着辅助参数β逐渐增大，目标函数(5)会渐近收敛于目标函数(4)。由此，将目标函数(4)分解为以下两个子问题进行迭代：子问题一：将地震数据矩阵X视为已知量，求解下式(6)以求取辅助变量矩阵Θ，其中F(Θ)为子问题一目标函数。子问题二：将辅助变量矩阵Θ视为已知量，求解下式(7)以求取迭代更新后的地震数据矩阵X，其中上式中的F(X)为子问题二目标函数。不失一般性，假设已进行到第l次Bergman迭代。步骤四：将地震数据矩阵X视为已知量，输入核范数约束系数λ与辅助参数β，采用奇异值软阈值化算法求解子问题一，求取辅助变量矩阵Θ。一般地，在第一次Bergman迭代时，可将地本文档来自技高网...

【技术保护点】
1.核范数与广义全变差联合约束的地震随机噪声压制方法，其特征在于，包括以下步骤：/n(1)输入三维地震数据体，将地震数据按Inline线进行输入，再对每条Inline线地震数据按时间采样点与CMP号排列成矩阵，其中时间采样点按行排列，CMP按列排列；/n(2)构建目标函数；/n(3)(3)采用分裂Bergman迭代策略，将目标函数分解为两个子问题进行迭代，针对具体子问题采用相应算法进行求解；/n(4)采用奇异值软阈值化算法求解子问题一，求取辅助变量矩阵；/n(5)采用梯度投影算法求解子问题二，求取地震数据矩阵；/n(6)判断Bergman迭代是否停止，若停止迭代则输出随机噪声压制后的地震数据矩阵；/n(7)对每条Inline线地震数据都进行随机噪声压制处理，得到整个随机噪声压制后的地震数据体。/n

【技术特征摘要】
1.核范数与广义全变差联合约束的地震随机噪声压制方法，其特征在于，包括以下步骤：
(1)输入三维地震数据体，将地震数据按Inline线进行输入，再对每条Inline线地震数据按时间采样点与CMP号排列成矩阵，其中时间采样点按行排列，CMP按列排列；
(2)构建目标函数；
(3)(3)采用分裂Bergman迭代策略，将目标函数分解为两个子问题进行迭代，针对具体子问题采用相应算法进行求解；
(4)采用奇异值软阈值化算法求解子问题一，求取辅助变量矩阵；
(5)采用梯度投影算法求解子问题二，求取地震数据矩阵；
(6)判断Bergman迭代是否停止，若停止迭代则输出随机噪声压制后的地震数据矩阵；
(7)对每条Inline线地震数据都进行随机噪声压制处理，得到整个随机噪声压制后的地震数据体。

【专利技术属性】
技术研发人员：代荣获，
申请(专利权)人：西华师范大学，
类型：发明
国别省市：四川;51