基于Lp伪范数和γ范数稀疏低秩约束的地震去噪方法

基于lp伪范数和
γ
范数稀疏低秩约束的地震去噪方法
技术领域
1.本发明涉及油气勘探及地震资料处理领域，尤其是一种基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法。


背景技术：

2.野外地震记录会受到随机噪声干扰，这对于解释工作带来了不小的麻烦，因此地震记录去噪是油气勘探中非常重要的过程。基于稀疏低秩理论的地震记录去噪技术是地震去噪的重要方法，原理是利用时频域地震记录的稀疏性质和低秩性质进行去噪处理。
3.基于l1范数稀疏约束和核范数低秩约束的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法是基于稀疏低秩理论地震去噪的重要方法。该方法采用l1范数对时频域的地震记录进行稀疏约束，采用核范数对时频域的地震记录进行低秩约束，提高了地震去噪能力。anvari等人将l1范数和核范数用于时频域的地震去噪中，证明了该方法的可行性。但是l1范数并不能完全地挖掘稀疏性，核范数也不能精确的表示数据矩阵的秩，导致去噪结果噪音仍有残余且有效信息存在损失，因此需要一种去噪方法克服以上的问题。


技术实现要素：

4.本发明的目的在于：本发明提供了一种基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法，解决了现有基于l1范数稀疏约束和核范数低秩约束的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法中所造成的去噪结果噪音仍有残余且有效信息存在损失的问题，提高去噪效果。
5.本发明采用的技术方案如下：
6.根据本发明的第一方面，提供一种基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法，所述方法包括：根据地震资料得到时频域的地震记录根据时频域的地震记录构建基于l
p
伪范数稀疏约束和γ范数低秩约束的目标函数；采用交替方向乘子法对目标函数进行迭代计算得到时频域去噪后的地震记录将时频域去噪后的地震记录转换到时空域。
7.进一步，所述根据地震资料得到时频域的地震记录包括：输入时空域含噪的地震记录s和子波数据w；利用子波数据w确定短时傅里叶变换的窗口长度，滑动距离，并对时空域含噪的地震记录s进行短时傅里叶变换，得到时频域的地震记录
8.进一步，所述根据时频域的地震记录构建基于l
p
伪范数稀疏约束和γ范数低秩约束的目标函数，包括：
9.根据时频域的地震记录构建目标函数中的保真项，如公式(1)所示：
10.11.其中，表示时频域含噪的地震记录，表示时频域去噪的地震记录，表示frobenius范数；
12.对时频域去噪的地震记录加上l
p
伪范数进行稀疏约束，如公式(2)所示：
[0013][0014]
其中，λ1表示l
p
伪范数稀疏约束的权重系数，表示l
p
伪范数，其表达式如公式(3)所示：
[0015][0016]
其中a、b、m和n分别表示时频地震记录的第a行、第b列、总行数和总列数；
[0017]
对时频域去噪的地震记录加上γ范数进行低秩约束，如公式(4)所示：
[0018][0019]
其中，λ2表示γ范数低秩约束的权重系数，||||
γ
表示γ范数，其表达式如公式(5)所示；
[0020][0021]
其中，表示时频域去噪的地震记录经过奇异值分解后的第k个特征值；
[0022]
引入拉格朗日乘子项r和对偶项c，得到基于l
p
伪范数和γ范数稀疏低秩约束的地震去噪目标函数，如公式(6)所示：
[0023][0024]
其中λ3表示对偶项权重系数。
[0025]
进一步，所述采用交替方向乘子法对目标函数进行迭代计算得到时频域去噪后的地震记录包括：
[0026]
分别设置稀疏约束项、低秩约束项和对偶项的权重系数λ1、λ2和λ3的值，以及l
p
伪范数中参数p的值和γ范数中参数γ的值；
[0027]
设置拉格朗日乘子项初始值r1＝0及其对偶项初始值c1＝0；
[0028]
根据交替方向乘子算法以及软阈值收缩算法更新时频域待去噪的地震记录计算公式如公式(7)所示：
[0029][0030]
其中，表示矩阵点乘运算，i 1表示第i 1次更新，sign表示符号函数，表达式如公式(8)所示：
[0031][0032]
其中，y
a,b
表示矩阵中的第a行，第b列元素；
[0033]
根据交替方向乘子算法以及差分凸规划更新拉格朗日乘子项r
i 1
，计算公式如公式(9)所示：
[0034]ri 1
＝u
·
diag(σ
*
)
·vt
ꢀꢀꢀ
(27)
[0035]
其中，diag表示对角矩阵，u是m
×
m阶酉矩阵，σ
*
是特征值，v
t
是n
×
n阶酉矩阵，矩阵u和v
t
由矩阵(x
i 1
ci)经过奇异值分解后得到，奇异值分解如公式(10)所示：
[0036]
svd(x
i 1
ci)＝u
·
diag(σ)
·vt
ꢀꢀꢀ
(28)
[0037]
其中，svd表示奇异值分解，通过如下公式(11)对特征值σ进行更新：
[0038][0039]
其中，t 1表示第t 1次迭代，经过多次迭代后，σ
t 1
收敛到最优点σ
*
；
[0040]
根据交替方向乘子算法以及费马引理更新对偶项c
i 1
，如公式(12)所示：
[0041]ci 1
＝ci x
i 1-r
i 1
ꢀꢀꢀ
(30)
[0042]
判断若是，则根据公式(7)更新时频域待去噪的地震记录若否，则结束循环。
[0043]
进一步，所述将时频域去噪后的地震记录转换到时空域，包括：利用子波数据w确定短时傅里叶逆变换的窗口长度，滑动距离；利用短时傅里叶逆变换将时频域去噪的地震记录转换到时空域。
[0044]
根据本发明的第二方面，提供一种基于lp伪范数和γ范数稀疏低秩约束的地震去噪装置，所述装置包括：预处理模块，被配置为根据地震资料得到时频域的地震记录目标函数构建模块，被配置为根据时频域的地震记录构建基于l
p
伪范数稀疏约束和γ范数低秩约束的目标函数；去噪模块，被配置为采用交替方向乘子法对目标函数进行迭代计算得到时频域去噪后的地震记录转换模块，被配置为将时频域去噪后的地震记录转换到时空域。
[0045]
进一步，所述去噪模块被进一步配置为：
[0046]
分别设置稀疏约束项、低秩约束项和对偶项的权重系数λ1、λ2和λ3的值，以及l
p
伪范数中参数p的值和γ范数中参数γ的值；
[0047]
设置拉格朗日乘子项初始值r1＝0及其对偶项初始值c1＝0；
[0048]
根据交替方向乘子算法以及软阈值收缩算法更新时频域待去噪的地震记录计算公式如公式(7)所示：
[0049][0050]
其中，表示矩阵点乘运算，i 1表示第i 1次更新，sign表示符号函数，表达式如公式(8)所示：
[0051][0052]
其中，y
a,b
表示矩阵中的第a行，第b列元素；
[0053]
根据交替方向乘子算法以及差分凸规划更新拉格朗日乘子项r
i 1
，计算公式如公式(9)所示：
[0054]ri 1
＝u
·
diag(σ
*
)
·vt
ꢀꢀꢀ
(33)
[0055]
其中，diag表示对角矩阵，u是m
×
m阶酉矩阵，σ
*
是特征值，v
t
是n
×
n阶酉矩阵，矩阵u和v
t
由矩阵(x
i 1
ci)经过奇异值分解后得到，奇异值分解如公式(10)所示：
[0056]
svd(x
i 1
ci)＝u
·
diag(σ)
·vt
ꢀꢀꢀ
(34)
[0057]
其中，svd表示奇异值分解，通过如下公式(11)对特征值σ进行更新：
[0058][0059]
其中，t 1表示第t 1次迭代，经过多次迭代后，σ
t 1
收敛到最优点σ
*
；
[0060]
根据交替方向乘子算法以及费马引理更新对偶项c
i 1
，如公式(12)所示：
[0061]ci 1
＝ci x
i 1-r
i 1
ꢀꢀꢀ
(36)
[0062]
判断若是，则根据公式(7)更新时频域待去噪的地震记录若否，则结束循环。
[0063]
根据本发明的第三方面，提供一种电子设备，包括：控制器；存储器，用于存储一个或多个程序，当所述一个或多个程序被所述控制器执行时，使得所述控制器实现本发明各个实施例所述的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法。
[0064]
根据本发明的第四方面，提供一种计算机可读存储介质，其上存储有计算机可读指令，当所述计算机可读指令被计算机的处理器执行时，使计算机执行本发明各个实施例所述的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法。
[0065]
根据本发明的第五方面，提供一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存储介质中，其特征在于，计算机设备的处理器从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算机设备执行本发明各个实施例所述的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法。
[0066]
本发明在l1范数稀疏约束和核范数低秩约束的基础上，用l
p
伪范数和γ范数分别替换l1范数和核范数，并结合交替方向乘子法提出了一种基于l
p
伪范数稀疏约束和γ范数低秩约束的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法，对含噪地震记录采用上述中的去噪方法进行去噪处理，通过反复迭代，输出最优去噪结果。本发明解决了现有l1稀疏约束和核范数低秩约束去噪结果噪音仍有残余且有效信息存在损失的问题，提高去噪效果。
附图说明
[0067]
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本技术的实施例，并与说明书一起用于解释本技术的原理。显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域普通技术者来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
[0068]
图1为本发明的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法流程框图；
[0069]
图2为本发明的无噪声的地震剖面；
[0070]
图3为本发明的含噪声的地震剖面；
[0071]
图4为本发明的去噪结果示意图；
[0072]
图5为本发明的所去噪声示意图；
[0073]
图6为本发明的装置结构图；
[0074]
图7为本发明的电子设备的计算机系统的结构示意图。
具体实施方式
[0075]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明，即所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
[0076]
因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0077]
需要说明的是，术语“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0078]
下面结合图1-5对本发明作详细说明。
[0079]
本发明解决的技术问题：解决了l1范数稀疏约束和核范数低秩约束导致的去噪结果噪音仍有残余且有效信息存在损失的问题，提高去噪效果。
[0080]
请参阅图1，是本发明实施例的基于lp伪范数和γ范数稀疏低秩约束的地震去噪
方法的整体流程框图。本发明实施例提供一种基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法，包括如下步骤：
[0081]
步骤1：对地震资料进行预处理并获取时频域的地震记录
[0082]
步骤2：根据时频域的地震记录构建基于l
p
伪范数稀疏约束和γ范数低秩约束的目标函数。
[0083]
步骤3：采用交替方向乘子法对目标函数进行迭代计算得到时频域去噪后的地震记录
[0084]
步骤4：将时频域去噪后的地震记录转换到时空域。
[0085]
在一些实施例中，步骤1包括如下步骤：
[0086]
步骤1.1：输入时空域含噪的地震记录s和子波数据w；
[0087]
步骤1.2：利用子波数据w确定短时傅里叶变换的窗口长度，滑动距离，并对时空域含噪的地震记录s进行短时傅里叶变换，得到时频域的地震记录
[0088]
在一些实施例中，步骤2包括如下步骤：
[0089]
步骤2.1：根据时频域的地震记录构建目标函数中的保真项，如公式1所示：
[0090][0091]
其中，表示时频域含噪的地震记录，表示时频域去噪的地震记录，表示frobenius范数；
[0092]
步骤2.2：在公式(1)的基础上，对时频域去噪的地震记录加上l
p
伪范数进行稀疏约束，如公式2所示：
[0093][0094]
其中，λ1表示l
p
伪范数稀疏约束的权重系数，表示l
p
伪范数，其表达式如公式3所示：
[0095][0096]
其中a、b、m和n分别表示时频地震记录的第a行、第b列、总行数和总列数；
[0097]
步骤2.3：在公式(2)的基础上，对时频域去噪的地震记录加上γ范数进行低秩约束，如公式4所示：
[0098][0099]
其中，λ2表示γ范数低秩约束的权重系数，||||
γ
表示γ范数，其表达式如公式5所示；
[0100][0101]
其中，表示时频域去噪的地震记录经过奇异值分解后的第k个特征值；
[0102]
步骤2.4：在公式(4)的基础上，引入拉格朗日乘子项r和对偶项c，得到基于l
p
伪范数和γ范数稀疏低秩约束的地震去噪目标函数，如公式4所示：
[0103][0104]
其中λ3表示对偶项权重系数。
[0105]
在一些实施例中，步骤3包括如下步骤：
[0106]
步骤3.1：分别设置稀疏约束项、低秩约束项和对偶项的权重系数λ1、λ2和λ3的值，l
p
伪范数和γ范数中p与γ的值。
[0107]
步骤3.2：设置拉格朗日乘子项初始值r1＝0及其对偶项初始值c1＝0。
[0108]
步骤3.3：根据交替方向乘子算法以及软阈值收缩算法更新时频域待去噪的地震记录计算公式如公式5所示：
[0109][0110]
其中，表示矩阵点乘运算，i 1表示第i 1次更新，sign表示符号函数，表达式如公式6所示：
[0111][0112]
其中，y
a,b
表示矩阵中的第a行，第b列元素。
[0113]
步骤3.4：根据交替方向乘子算法以及差分凸规划更新拉格朗日乘子项r
i 1
，计算公式如公式9所示：
[0114]ri 1
＝u
·
diag(σ
*
)
·vt
ꢀꢀꢀ
(45)
[0115]
其中，diag表示对角矩阵，u是m
×
m阶酉矩阵，σ
*
是特征值，v
t
是n
×
n阶酉矩阵，矩阵u和v
t
由矩阵(x
i 1
ci)经过奇异值分解后得到，如公式10所示：
[0116]
svd(x
i 1
ci)＝u
·
diag(σ)
·vt
ꢀꢀꢀ
(46)
[0117]
其中，svd表示奇异值分解。再对特征值σ进行更新，如公式11所示：
[0118][0119]
其中，t 1表示第t 1次迭代，经过几次迭代后，σ
t 1
收敛到最优点σ
*
。
[0120]
步骤3.5：根据交替方向乘子算法以及费马引理更新对偶项c
i 1
，如公式12所示：
[0121]ci 1
＝ci x
i 1-r
i 1
ꢀꢀꢀ
(48)
[0122]
步骤3.6：判断如果是，则返回步骤3.3进行循环，如果否，结束循环。
[0123]
在一些实施例中，步骤4包括如下步骤：
[0124]
步骤4.1：利用子波数据w确定短时傅里叶逆变换的窗口长度，滑动距离；
[0125]
步骤4.2：利用短时傅里叶逆变换将时频域去噪的地震记录转换到时空域。
[0126]
综上所述，本发明在l1范数稀疏约束和核范数低秩约束的基础上，将l1范数和核范数分别用l
p
伪范数和γ范数替换，解决了l1范数稀疏约束和核范数低秩约束造成的去噪结果噪音仍有残余且有效信息存在损失的问题，提高去噪效果。
[0127]
下面本发明实施例将结合具体的实验以进一步说明本发明的可行性和进步性。基于如上所述的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法，本发明实施例对野外地震记录进行试验。所述野外地震记录是在油气勘探中利用相关地震设备所采集到的地震记录。野外地震记录至少包括时空域含噪的地震记录s和子波数据w。
[0128]
在利用本发明实施例所述的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法对野外地震记录进行去噪处理，结果如图4和图5所示。将图4和图2进行对比，其中图4为本发明的去噪结果示意图，图2为本发明的无噪声的地震剖面，可以发现去噪后的地震记录与真实的无噪声的地震记录非常的接近，并结合图3，图3为本发明的含噪声的地震剖面，可以发现地震记录的去噪结果剖面相较于含噪声的地震记录剖面能更清晰的反映地下的层位和构造，证明了本发明的正确性。图5中可以看到所去除的噪声剖面上没有出现层位结构性的特征，说明去噪后地震记录的损失非常小，进一步证明了该方法的正确性。本发明在l1范数稀疏约束和核范数低秩约束的基础上，用l
p
伪范数和γ范数分别替换l1范数和核范数，并结合交替方向乘子法提出了一种基于l
p
伪范数稀疏约束和γ范数低秩约束的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法。对含噪地震记录采用上述中的去噪方法进行去噪处理，通过反复迭代，输出最优去噪结果。本发明成功解决了现有l1稀疏约束和核范数低秩约束去噪结果噪音仍有残余且有效信息存在损失的问题，提高去噪效果。
[0129]
请参阅图6，是本发明的装置结构图。本发明实施例还提供一种基于lp伪范数和γ范数稀疏低秩约束的地震去噪装置，所述装置600包括：
[0130]
预处理模块601，被配置为根据地震资料得到时频域的地震记录
[0131]
目标函数构建模块602，被配置为根据时频域的地震记录构建基于l
p
伪范数稀疏约束和γ范数低秩约束的目标函数；
[0132]
去噪模块603，被配置为采用交替方向乘子法对目标函数进行迭代计算得到时频域去噪后的地震记录
[0133]
转换模块604，被配置为将时频域去噪后的地震记录转换到时空域。
[0134]
在一些实施例中，所述预处理模块被进一步配置为：输入时空域含噪的地震记录s和子波数据w；利用子波数据w确定短时傅里叶变换的窗口长度，滑动距离，并对时空域含噪的地震记录s进行短时傅里叶变换，得到时频域的地震记录
[0135]
在一些实施例中，所述目标函数构建模块被进一步配置为：
[0136]
根据时频域的地震记录构建目标函数中的保真项，如公式(1)所示：
[0137]
[0138]
其中，表示时频域含噪的地震记录，表示时频域去噪的地震记录，表示frobenius范数；
[0139]
对时频域去噪的地震记录加上l
p
伪范数进行稀疏约束，如公式(2)所示：
[0140][0141]
其中，λ1表示l
p
伪范数稀疏约束的权重系数，表示l
p
伪范数，其表达式如公式(3)所示：
[0142][0143]
其中a、b、m和n分别表示时频地震记录的第a行、第b列、总行数和总列数；
[0144]
对时频域去噪的地震记录加上γ范数进行低秩约束，如公式(4)所示：
[0145][0146]
其中，λ2表示γ范数低秩约束的权重系数，||||
γ
表示γ范数，其表达式如公式(5)所示；
[0147][0148]
其中，表示时频域去噪的地震记录经过奇异值分解后的第k个特征值；
[0149]
引入拉格朗日乘子项r和对偶项c，得到基于l
p
伪范数和γ范数稀疏低秩约束的地震去噪目标函数，如公式(6)所示：
[0150][0151]
其中λ3表示对偶项权重系数。
[0152]
在一些实施例中，所述去噪模块被进一步配置为：
[0153]
分别设置稀疏约束项、低秩约束项和对偶项的权重系数λ1、λ2和λ3的值，以及l
p
伪范数中参数p的值和γ范数中参数γ的值；
[0154]
设置拉格朗日乘子项初始值r1＝0及其对偶项初始值c1＝0；
[0155]
根据交替方向乘子算法以及软阈值收缩算法更新时频域待去噪的地震记录计算公式如公式(7)所示：
[0156][0157]
其中，表示矩阵点乘运算，i 1表示第i 1次更新，sign表示符号函数，表达式如公式(8)所示：
[0158][0159]
其中，y
a,b
表示矩阵中的第a行，第b列元素；
[0160]
根据交替方向乘子算法以及差分凸规划更新拉格朗日乘子项r
i 1
，计算公式如公式(9)所示：
[0161]ri 1
＝u
·
diag(σ
*
)
·vt
ꢀꢀꢀ
(57)
[0162]
其中，diag表示对角矩阵，u是m
×
m阶酉矩阵，σ
*
是特征值，v
t
是n
×
n阶酉矩阵，矩阵u和v
t
由矩阵(x
i 1
ci)经过奇异值分解后得到，奇异值分解如公式(10)所示：
[0163]
svd(x
i 1
ci)＝u
·
diag(σ)
·vt
ꢀꢀꢀ
(58)
[0164]
其中，svd表示奇异值分解，通过如下公式(11)对特征值σ进行更新：
[0165][0166]
其中，t 1表示第t 1次迭代，经过多次迭代后，σ
t 1
收敛到最优点σ
*
；
[0167]
根据交替方向乘子算法以及费马引理更新对偶项c
i 1
，如公式(12)所示：
[0168]ci 1
＝ci x
i 1-r
i 1
ꢀꢀꢀ
(60)
[0169]
判断若是，则根据公式(7)更新时频域待去噪的地震记录若否，则结束循环。
[0170]
在一些实施例中，所述转换模块被进一步配置为：利用子波数据w确定短时傅里叶逆变换的窗口长度，滑动距离；利用短时傅里叶逆变换将时频域去噪的地震记录转换到时空域。
[0171]
需要说明的是，上述实施例所提供的地震去噪装置与前述实施例所提供的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法属于同一构思，其中各个模块和单元执行操作的具体方式已经在方法实施例中进行了详细描述，这里不再赘述。
[0172]
本发明实施例提供了一种电子设备，包括：控制器；存储器，用于存储一个或多个程序，当一个或多个程序被控制器执行时，以执行上述各个实施例中所述的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法。
[0173]
请参阅图7，图7是本发明的一示例性实施例示出的电子设备的计算机系统的结构示意图，其示出了适于用来实现本技术实施例的电子设备的计算机系统的结构示意图。
[0174]
需要说明的是，图7示出的电子设备的计算机系统700仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。
[0175]
如图7所示，计算机系统700包括中央处理单元(central processing unit，cpu)701，其可以根据存储在只读存储器(read-only memory，rom)702中的程序或者从存储部分708加载到随机访问存储器(random access memory，ram)703中的程序而执行各种适当的动作和处理，例如执行上述实施例中的方法。在ram 703中，还存储有系统操作所需的各种
程序和数据。cpu 701、rom702以及ram 703通过总线704彼此相连。输入/输出(input/output，i/o)接口705也连接至总线704。
[0176]
以下部件连接至i/o接口705：包括键盘、鼠标等的输入部分706；包括诸如阴极射线管(cathode ray tube，crt)、液晶显示器(liquid crystal display，lcd)等以及扬声器等的输出部分707；包括硬盘等的存储部分708；以及包括诸如lan(local area network，局域网)卡、调制解调器等的网络接口卡的通信部分709。通信部分709经由诸如因特网的网络执行通信处理。驱动器710也根据需要连接至i/o接口705。可拆卸介质711，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器710上，以便于从其上读出的计算机程序根据需要被安装入存储部分708。
[0177]
特别地，根据本发明的实施例，上文所述的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法中描述的各个步骤过程可以被实现为计算机软件程序。例如，本发明实施例包括一种计算机程序产品，其包括承载在计算机可读介质上的计算机程序，该计算机程序包含用于执行流程图所示的方法的计算机程序。在这样的实施例中，该计算机程序可以通过通信部分709从网络上被下载和安装，和/或从可拆卸介质711被安装。在该计算机程序被中央处理单元(cpu)701执行时，执行本发明的系统中限定的各种功能。
[0178]
需要说明的是，本发明实施例所示的计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质或者是上述两者的任意组合。计算机可读存储介质例如可以是电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子可以包括但不限于：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(ram)、只读存储器(rom)、可擦式可编程只读存储器(erasable programmable read only memory，eprom)、闪存、光纤、便携式紧凑磁盘只读存储器(compact disc read-only memory，cd-rom)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本发明中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。而在本发明中，计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的计算机程序。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的计算机程序可以用任何适当的介质传输，包括但不限于：无线、有线等等，或者上述的任意合适的组合。
[0179]
描述于本发明实施例中所涉及到的模块可以通过软件的方式实现，也可以通过硬件的方式来实现，所描述的模块也可以设置在处理器中。其中，这些模块的名称在某种情况下并不构成对该模块本身的限定。
[0180]
本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如前的识别方法。该计算机可读存储介质可以是上述实施例中描述的电子设备中所包含的，也可以是单独存在，而未装配入该电子设备中。
[0181]
本发明实施例还提供了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算
机设备执行上述各个实施例中提供的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法。
[0182]
本发明实施例还提供了一种计算机系统，包括中央处理单元(central processing unit，cpu)，其可以根据存储在只读存储器(read-only memory，rom)中的程序或者从存储部分加载到随机访问存储器(random access memory，ram)中的程序而执行各种适当的动作和处理，例如执行上述实施例中的基于lp伪范数和γ范数稀疏低秩约束的地震去噪方法。在ram中，还存储有系统操作所需的各种程序和数据。cpu、rom以及ram通过总线彼此相连。输入/输出(input/output，i/o)接口也连接至总线。
[0183]
以下部件连接至i/o接口：包括键盘、鼠标等的输入部分；包括诸如阴极射线管(cathode ray tube，crt)、液晶显示器(liquid crystal display，lcd)等以及扬声器等的输出部分；包括硬盘等的存储部分；以及包括诸如lan(local area network，局域网)卡、调制解调器等的网络接口卡的通信部分。通信部分经由诸如因特网的网络执行通信处理。驱动器也根据需要连接至i/o接口。可拆卸介质，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器上，以便于从其上读出的计算机程序根据需要被安装入存储部分。
[0184]
上述内容，仅为本发明的较佳示例性实施例，并非用于限制本发明的实施方案，本领域普通技术人员根据本发明的主要构思和精神，可以十分方便地进行相应的变通或修改，故本发明的保护范围应以权利要求书所要求的保护范围为准。