[go: up one dir, main page]

CN108710150A - 一种基于稳健奇异谱分析的地震不规则噪声去除方法 - Google Patents

一种基于稳健奇异谱分析的地震不规则噪声去除方法 Download PDF

Info

Publication number
CN108710150A
CN108710150A CN201810493096.8A CN201810493096A CN108710150A CN 108710150 A CN108710150 A CN 108710150A CN 201810493096 A CN201810493096 A CN 201810493096A CN 108710150 A CN108710150 A CN 108710150A
Authority
CN
China
Prior art keywords
matrix
seismic data
hankel matrix
hankel
factor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810493096.8A
Other languages
English (en)
Other versions
CN108710150B (zh
Inventor
刘春成
李勇
张益明
王彬
牛聪
周钰邦
黄饶
叶云飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Chengdu Univeristy of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China National Offshore Oil Corp CNOOC, CNOOC Research Institute Co Ltd, Chengdu Univeristy of Technology filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201810493096.8A priority Critical patent/CN108710150B/zh
Publication of CN108710150A publication Critical patent/CN108710150A/zh
Application granted granted Critical
Publication of CN108710150B publication Critical patent/CN108710150B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364—Seismic filtering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Acoustics & Sound (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种基于稳健奇异谱分析的地震不规则噪声去除方法，其特征在于，包括以下步骤：将单频地震数据嵌入Hankel矩阵，构建单频地震数据的Hankel矩阵；基于L1和L2混合范数，得到单频地震数据的Hankel矩阵的目标函数；根据预先设定的所述目标函数中的正则化参数，对所述目标函数进行求解，得到基于L1和L2混合范数的Hankel矩阵的奇异谱分析；对所述奇异谱分析依次进行交替最小化和加权最小二乘法最小化，求解得到单频地震数据的Hankel矩阵的因子矩阵，进而得到单频地震数据的Hankel矩阵的估计，完成单频地震数据的降噪，本发明可以广泛应用于地震数据处理领域中。

Description

一种基于稳健奇异谱分析的地震不规则噪声去除方法
技术领域
本发明是关于一种基于稳健奇异谱分析的地震不规则噪声去除方法，属于地震数据处理领域。
背景技术
地震数据的噪声压制是地震数据处理中的一个重要流程.针对不同的地震干扰类型有不同的去噪方法，但最大的挑战仍然是如何做到保幅去噪，能否做到保幅去噪主要取决于对混杂在一起的信号和干扰进行有效分离的程度，调和分析方法通过诸如旋转、伸缩、平移等变换处理，将地震数据从原始的记录空间投影到不同的正交子空间，为研究人员提供了另一种分析和处理地震数据的视角，因此，如何针对地震数据的特点设计有针对性的特征函数或特征记录成为一个重要的研究方向。
不规则/或相干噪声衰减主要通过f-x和t-x域中的预测误差滤波器来进行，当然不规则噪声衰减也可以通过降秩方法实现，降秩方法认为地震波在地下介质中传播遵循波动规律，某一频率切片上的地震数据应呈现线性相关性，其对应的轨迹矩阵为低秩矩阵，而地震数据的高斯噪声和不规则噪声破坏了地震数据的相关性，导致地震轨迹矩阵的秩增加，因此，在高保真的条件下，地震数据降噪的过程就是对轨迹矩阵降秩的过程。降秩方法可以分为不同的类别，例如通过Karhunen-Loève变换过滤的特征图像滤波可以直接在t-x或f-x-y域的地震数据上进行降秩。
最近，对地震数据的噪声压制处理引入了奇异频谱分析方法(SSA)，也称为Cadzow滤波，以减弱非相干噪声和重建地震数据。奇异频谱分析方法首先将地震数据嵌入Hankel矩阵(汉克尔矩阵)，导致Hankel矩阵的秩增加；然后通过降秩算法进行降噪，秩降低通常通过奇异值分解(SVD)来实现，SVD是一种非鲁棒矩阵分解方法，当存在不规则噪声时，常常导致次优结果。上述不规则噪声是指非高斯噪声，由已知或未知分布的大型孤立事件引入。采用稳健范数之后，基于奇异频谱分析的降秩滤波也用于抑制相干噪声，归纳起来，现有的基于降秩方法的地震不规则噪声去除方法主要为通过对Hankel矩阵进行基于L2范数的截断SVD方法，或采用基于稳健范数(例如turkey范数)的截断SVD方法，其核心问题是基于L2范数的截断SVD方法无法真正有效地去除不规则噪声，因为L2范数的核心本质就是求平均。基于turkey范数的截断SVD方法中，turkey范数是凹函数，其优化问题本质上就是非凸优化，无法找到全局最优解以及智能获得局部最优解，且该方法依赖初始值，从而会影响降噪的性能。
发明内容
针对上述问题，本发明的目的是提供一种能够智能获取全局最优解且不会影响降噪性能的基于稳健奇异谱分析的地震不规则噪声去除方法。
为实现上述目的，本发明采取以下技术方案：一种基于稳健奇异谱分析的地震不规则噪声去除方法，其特征在于，包括以下步骤：将单频地震数据嵌入Hankel矩阵，构建单频地震数据的Hankel矩阵；基于L1和L2混合范数，得到单频地震数据的Hankel矩阵的目标函数；根据预先设定的所述目标函数中的正则化参数，对所述目标函数进行求解，得到基于L1和L2混合范数的Hankel矩阵的奇异谱分析；对所述奇异谱分析依次进行交替最小化和加权最小二乘法最小化，求解得到单频地震数据的Hankel矩阵的因子矩阵，进而得到单频地震数据的Hankel矩阵的估计，完成单频地震数据的降噪。
进一步，所述目标函数为：
其中，λ1和λ2均表示正则化参数，λ1用于调节L1范数的权值大小，λ2用于调节L2范数的权值大小；M表示Hankel矩阵；表示Hankel矩阵的估计。
进一步，当正则化参数λ1＝0时，所述目标函数退化为：
所述目标函数即为L1范数下的拟合解，用于压制非高斯噪声；当正则化参数λ2＝0时，所述目标函数退化为：
所述目标函数即为最小二乘解，用于压制地震数据的高斯背景噪声。
进一步，所述奇异谱分析为：
其中，argminJ表示退化后的目标函数；m表示Hankel矩阵的行向量；U和V分别表示奇异谱分析的因子矩阵。
进一步，对所述奇异谱分析依次进行交替最小化和加权最小二乘法最小化，求解得到单频地震数据的Hankel矩阵的因子矩阵，进而得到单频地震数据的Hankel矩阵的估计，完成单频地震数据的降噪，具体为：对基于L1和L2混合范数下的Hankel矩阵的奇异谱分析进行交替最小化：
其中，Ew(V)表示自变量为因子矩阵V的目标函数，W表示权值矩阵，mj表示Hankel矩阵的第j列，且M＝(m1，m2，…，mn)＝(m1，m2，…，mn)H，mn表示Hankel矩阵的第n行的共轭转置；vj表示因子矩阵V的列向量；Ew(U)表示自变量为因子矩阵U的目标函数；uj表示因子矩阵U的列向量；因子矩阵U和V为：
U＝(u1，u2，…，uK)＝(u1，u2，…，um)H
V＝(v1，v2，…，vK)＝(v1，v2，…，vn)H
其中，uK表示因子矩阵U的列向量，um表示因子矩阵U的行向量，vK表示因子矩阵V的列向量，vn表示因子矩阵V的行向量；对进行交替最小化后的奇异谱分析进行加权最小二乘法最小化：
其中，mi表示Hankel矩阵的第i行的共轭转置；采用QR分解，迭代求解上述公式中的加权最小二乘法最小化问题，得到单频地震数据的Hankel矩阵的因子矩阵U和V，进而得到单频地震数据的Hankel矩阵的估计完成单频地震数据的降噪。
本发明由于采取以上技术方案，其具有以下优点：1、本发明基于L1和L2混合范数，代替截断SVD采用的二次误差准则函数，对单频地震数据进行降噪，其中，L2范数能够对高斯噪声进行降噪，L1范数能够对非高斯规则噪声进行降噪，从而能够有效地同时去除高斯噪声和非高斯噪声。2、由于L1范数和L2范数都是凸函数，采用L1和L2混合范数得到的目标函数也是凸函数，而最优化问题本质上是凸优化问题，容易获得全局最优解，因此，本发明的方法能够获得更好的降噪性能，可以广泛应用于地震资料处理领域中。
附图说明
图1是采用本发明方法对地震数据的Hankel矩阵进行降噪的对比示意图，其中，图1(a)表示干净地震数据的Hankel矩阵，图1(b)表示含噪地震数据的Hankel矩阵，图1(c)表示降噪后地震数据的Hankel矩阵；
图2是原始的叠后地震剖面示意图；
图3是采用本发明方法对图2中叠后地震剖面进行降噪后的示意图；
图4是采用本发明方法对图2中叠后地震剖面去除的干扰的示意图；
图5是原始的叠后地震剖面示意图；
图6是采用本发明方法对图5中叠后地震剖面进行降噪后的示意图；
图7是采用本发明方法对图5中叠后地震剖面去除的干扰的示意图。
具体实施方式
以下结合附图来对本发明进行详细的描绘。然而应当理解，附图的提供仅为了更好地理解本发明，它们不应该理解成对本发明的限制。
小窗口内的地震数据Dj(ω)通过频率-空间域内的平面波叠加得到：
其中，i和j表示空间轴上的地震道集索引值，且N表示地震道数目；ω表示时间频率，在上述公式中，假设地震数据Dj(ω)由具有不同射线参数Pk的k个线性事件(即平面波)组成，Ak(ω)表示第k个平面波的复振幅，Δx表示地震剖面之间的空间间隔，K表示线性事件的序数。
奇异谱分析方法是将地震数据Dj(ω)变换后取单独频率的数据得到单频地震数据D(ω)＝(D1(ω)，D2(ω），…，DN(ω))T，并将单频地震数据嵌入信号，构建轨迹矩阵，即Hankel矩阵M(ω)：
其中，表示Hankel算子，为方便起见，选择Hankel矩阵的行向量数目使Hankel矩阵逼近为方阵 表示复数域。省略符号ω，对所有频率进行分析，并对k个平面波进行叠加得到Hankel矩阵的秩rank(M)＝K。外加噪声D将增加矩阵M的秩K，需要通过降秩来减小外加噪声。
因此，奇异谱分析通过以下表达式表示：
其中，表示反对角平均算子；表示降秩操作符，通过秩k逼近Hankel矩阵M；表示Hankel算子；运算符通过平均对角线将Hankel形式转换成向量形式，对于多维信号则采用块Hankel矩阵和阻尼对角平均算子。如图1所示，降秩操作可以通过截断SVD、随机SVD或采用Lanczos双对角化方法，以及高效矩阵向量乘法的快速算法来实现，其中，快速算法采用快速傅里叶变换。
Hankel矩阵M的秩K可以通过求解下式近似得到：
其中，表示Frobenius范数；为矩阵m和n表示矩阵X的大小，xij表示矩阵X的元素；公式(4)具有全局唯一的局部最小值，这个解具有解析表达式，由截断奇异值分解(TSVD)给出：
其中，Uk和Vk表示包含奇异值向量的矩阵与秩k矩阵最大奇异值σq相关的向量，且q表示特征值或奇异值序号，q＝1，…，K；∑k表示特征值矩阵(即奇异值)，这些奇异值也是矩阵(表示实数域)的对角元素，称为Eckart-Young定理。即使公式(4)的解决方案是简单的，二次匹配函数使得截断奇异值分解对非高斯噪声相当敏感，这个缺点限制了奇异谱分析方法在其中包含异常值的应用。
基于上述现有奇异谱分析存在的问题，本发明提供的基于稳健奇异谱分析的地震不规则噪声去除方法，包括以下步骤：
1)将单频地震数据嵌入Hankel矩阵，构建单频地震数据的Hankel矩阵M。
2)基于L1和L2混合范数，代替截断SVD采用的二次误差准则函数，得到单频地震数据的Hankel矩阵M的目标函数，具体为：
如图1所示，传统的基于L2范数的奇异谱分析只能压制噪声中的高斯部分，而基于L1范数的奇异谱分析对呈现出厚拖尾的非高斯部分具有更好的鲁棒性，因此，结合L1范数与L2范数在不同噪声环境下的优势，得到基于L1和L2混合范数的Hankel矩阵M的目标函数J：
其中，λ1和λ2均表示正则化参数，λ1用于调节L1范数的权值大小，λ2用于调节L2范数的权值大小，表示Hankel矩阵M的估计。
3)设定单频地震数据的Hankel矩阵M的目标函数中的正则化参数λ1和λ2，其中，上述正则化参数λ1和λ2可以根据实际情况进行设定。
当λ1＝0时，上述目标函数(6)退化为：
公式(7)即为L1范数下的拟合解，用于压制非高斯噪声，噪声的概率密度函数可用拉普拉斯分布描述。
当λ2＝0时，上述目标函数(6)退化为：
公式(8)即为最小二乘解，用于压制地震数据的高斯背景噪声。
4)根据设定的正则化参数λ1和λ2，对单频地震数据的Hankel矩阵M的目标函数进行优化求解，得到基于L1和L2混合范数的Hankel矩阵M的奇异谱分析，即特征向量矩阵(U，V)：
其中，argminJ表示退化后的目标函数；m表示Hankel矩阵M的行向量；U和V分别表示特征向量矩阵的因子矩阵，且
5)对基L1和L2混合范数下的Hankel矩阵M的奇异谱分析依次进行交替最小化(ADMM)和加权最小二乘法最小化(IRLS)，求解得到单频地震数据的Hankel矩阵的因子矩阵U和V，进而得到单频地震数据的Hankel矩阵的稳健估计，完成单频地震数据的降噪，具体为：
5.1)对基于L1和L2混合范数下的Hankel矩阵M的奇异谱分析进行交替最小化：
其中，Ew(V)表示自变量为因子矩阵V的目标函数，W表示权值矩阵，mj表示Hankel矩阵M的第j列，且M＝(m1，m2，…，mn)＝(m1，m2，…，mn)H，mn表示Hankel矩阵M的第n行的共轭转置；vj表示因子矩阵V的列向量；Ew(U)表示自变量为因子矩阵U的目标函数；uj表示因子矩阵U的列向量；因子矩阵U和V可以表示为：
U＝(u1，u2，…，uK)＝(u1，u2，…，um)H (12)
V＝(v1，v2，…，vK)＝(v1，v2，…，vn)H (13)
其中，uK表示因子矩阵U的列向量，um表示因子矩阵U的行向量，vK表示因子矩阵V的列向量，vn表示因子矩阵V的行向量。
5.2)对进行交替最小化后的奇异谱分析进行加权最小二乘法最小化，采用QR分解，迭代求解直到收敛或到达最大迭代步数，得到单频地震数据的Hankel矩阵M的因子矩阵U和V，进而得到单频地震数据的Hankel矩阵的稳健估计，完成单频地震数据的降噪。
公式(10)和(11)可以分解成更小的加权最小二乘法最小化问题，即将公式(10)和(11)进行加权最小二乘法最小化：
其中，mi表示Hankel矩阵M的第i行的共轭转置。
采用QR分解，迭代求解上述公式(14)和(15)中的加权最小二乘法最小化问题，得到单频地震数据的Hankel矩阵M的因子矩阵U和V，进而得到单频地震数据的Hankel矩阵的稳健估计完成单频地震数据的降噪。
如图2～7所示，通过对比两组原始的叠后地震剖面，采用本发明的地震不规则噪声去除方法进行降噪后的叠后地震剖面，以及采用本发明的地震不规则噪声去除方法对叠后地震剖面去除的干扰，可以看出本发明的地震不规则噪声去除方法能够有效去除地震数据中的高斯和不规则噪声，且地震数据本身没有受到损害。
上述各实施例仅用于说明本发明，其中各部件的结构、连接方式和制作工艺等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。

Claims (5)

1.一种基于稳健奇异谱分析的地震不规则噪声去除方法，其特征在于，包括以下步骤：
将单频地震数据嵌入Hankel矩阵，构建单频地震数据的Hankel矩阵；
基于L1和L2混合范数，得到单频地震数据的Hankel矩阵的目标函数；
根据预先设定的所述目标函数中的正则化参数，对所述目标函数进行求解，得到基于L1和L2混合范数的Hankel矩阵的奇异谱分析；
对所述奇异谱分析依次进行交替最小化和加权最小二乘法最小化，求解得到单频地震数据的Hankel矩阵的因子矩阵，进而得到单频地震数据的Hankel矩阵的估计，完成单频地震数据的降噪。
2.如权利要求1所述的一种基于稳健奇异谱分析的地震不规则噪声去除方法，其特征在于，所述目标函数为：
其中，λ1和λ2均表示正则化参数，λ1用于调节L1范数的权值大小，λ2用于调节L2范数的权值大小；M表示Hankel矩阵；表示Hankel矩阵的估计。
3.如权利要求2所述的一种基于稳健奇异谱分析的地震不规则噪声去除方法，其特征在于，当正则化参数λ1＝0时，所述目标函数退化为：
所述目标函数即为L1范数下的拟合解，用于压制非高斯噪声；
当正则化参数λ2＝0时，所述目标函数退化为：
所述目标函数即为最小二乘解，用于压制地震数据的高斯背景噪声。
4.如权利要求2所述的一种基于稳健奇异谱分析的地震不规则噪声去除方法，其特征在于，所述奇异谱分析为：
其中，argminJ表示退化后的目标函数；m表示Hankel矩阵的行向量；U和V分别表示奇异谱分析的因子矩阵。
5.如权利要求1所述的一种基于稳健奇异谱分析的地震不规则噪声去除方法，其特征在于，对所述奇异谱分析依次进行交替最小化和加权最小二乘法最小化，求解得到单频地震数据的Hankel矩阵的因子矩阵，进而得到单频地震数据的Hankel矩阵的估计，完成单频地震数据的降噪，具体为：
对基于L1和L2混合范数下的Hankel矩阵的奇异谱分析进行交替最小化：
其中，Ew(V)表示自变量为因子矩阵V的目标函数，W表示权值矩阵，mj表示Hankel矩阵的第j列，且M＝(m1，m2，…，mn)＝(m1，m2，…，mn)H，mn表示Hankel矩阵的第n行的共轭转置；vj表示因子矩阵V的列向量；Ew(U)表示自变量为因子矩阵U的目标函数；uj表示因子矩阵U的列向量；因子矩阵U和V为：
U＝(u1，u2，…，uK)＝(u1，u2，…，um)H
V＝(v1，v2，…，vK)＝(v1，v2，…，vn)H
其中，uK表示因子矩阵U的列向量，um表示因子矩阵U的行向量，vK表示因子矩阵V的列向量，vn表示因子矩阵V的行向量；
对进行交替最小化后的奇异谱分析进行加权最小二乘法最小化：
其中，mi表示Hankel矩阵的第i行的共轭转置；
采用QR分解，迭代求解上述公式中的加权最小二乘法最小化问题，得到单频地震数据的Hankel矩阵的因子矩阵U和V，进而得到单频地震数据的Hankel矩阵的估计完成单频地震数据的降噪。
CN201810493096.8A 2018-05-22 2018-05-22 一种基于稳健奇异谱分析的地震不规则噪声去除方法 Active CN108710150B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810493096.8A CN108710150B (zh) 2018-05-22 2018-05-22 一种基于稳健奇异谱分析的地震不规则噪声去除方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810493096.8A CN108710150B (zh) 2018-05-22 2018-05-22 一种基于稳健奇异谱分析的地震不规则噪声去除方法

Publications (2)

Publication Number Publication Date
CN108710150A true CN108710150A (zh) 2018-10-26
CN108710150B CN108710150B (zh) 2019-09-06

Family

ID=63869361

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810493096.8A Active CN108710150B (zh) 2018-05-22 2018-05-22 一种基于稳健奇异谱分析的地震不规则噪声去除方法

Country Status (1)

Country Link
CN (1) CN108710150B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109633751A (zh) * 2018-12-26 2019-04-16 北京化工大学 一种基于多分辨奇异值分解的地震数据降噪方法
CN110728624A (zh) * 2019-09-29 2020-01-24 厦门大学 一种高分辨率扩散加权图像重建方法
CN111257929A (zh) * 2020-02-17 2020-06-09 成都理工大学 一种奇异值衰减的降秩去噪方法
CN111323816A (zh) * 2020-03-20 2020-06-23 中国海洋石油集团有限公司 基于海洋宽频地震数据波形的瞬时相位梯度属性提取方法
CN111474911A (zh) * 2020-04-28 2020-07-31 浙江浙能技术研究院有限公司 面向高端燃煤发电装备非平稳运行的高斯非高斯特征协同解析与监测方法
CN112700389A (zh) * 2021-01-13 2021-04-23 安徽工业大学 一种活性污泥微生物彩色显微图像去噪方法
CN113687421A (zh) * 2021-08-23 2021-11-23 中国石油大学(北京) 地震信号的数据处理方法、装置、电子设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103295250A (zh) * 2013-05-31 2013-09-11 东莞中山大学研究院 一种基于l1混合范数求解的图像彩色化的方法
EP2781936A2 (en) * 2013-03-22 2014-09-24 CGG Services SA System and method for interpolating seismic data
CN105930310A (zh) * 2016-05-06 2016-09-07 南京信息工程大学 一种自适应正则化平滑l0范数方法
CN105958470A (zh) * 2014-10-20 2016-09-21 国家电网公司 一种基于双线性凸优化理论电力系统双线性抗差估计方法
CN106204482A (zh) * 2016-07-08 2016-12-07 桂林电子科技大学 基于加权稀疏的混合噪声去除方法
CN107728211A (zh) * 2017-08-31 2018-02-23 电子科技大学 基于张量核范数正则化的地震信号恢复算法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2781936A2 (en) * 2013-03-22 2014-09-24 CGG Services SA System and method for interpolating seismic data
CN103295250A (zh) * 2013-05-31 2013-09-11 东莞中山大学研究院 一种基于l1混合范数求解的图像彩色化的方法
CN105958470A (zh) * 2014-10-20 2016-09-21 国家电网公司 一种基于双线性凸优化理论电力系统双线性抗差估计方法
CN105930310A (zh) * 2016-05-06 2016-09-07 南京信息工程大学 一种自适应正则化平滑l0范数方法
CN106204482A (zh) * 2016-07-08 2016-12-07 桂林电子科技大学 基于加权稀疏的混合噪声去除方法
CN107728211A (zh) * 2017-08-31 2018-02-23 电子科技大学 基于张量核范数正则化的地震信号恢复算法

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109633751A (zh) * 2018-12-26 2019-04-16 北京化工大学 一种基于多分辨奇异值分解的地震数据降噪方法
CN110728624A (zh) * 2019-09-29 2020-01-24 厦门大学 一种高分辨率扩散加权图像重建方法
CN110728624B (zh) * 2019-09-29 2021-07-23 厦门大学 一种高分辨率扩散加权图像重建方法
CN111257929A (zh) * 2020-02-17 2020-06-09 成都理工大学 一种奇异值衰减的降秩去噪方法
CN111257929B (zh) * 2020-02-17 2021-03-19 成都理工大学 一种奇异值衰减的降秩去噪方法
CN111323816A (zh) * 2020-03-20 2020-06-23 中国海洋石油集团有限公司 基于海洋宽频地震数据波形的瞬时相位梯度属性提取方法
CN111323816B (zh) * 2020-03-20 2021-12-14 中国海洋石油集团有限公司 基于海洋宽频地震数据波形的瞬时相位梯度属性提取方法
CN111474911A (zh) * 2020-04-28 2020-07-31 浙江浙能技术研究院有限公司 面向高端燃煤发电装备非平稳运行的高斯非高斯特征协同解析与监测方法
CN111474911B (zh) * 2020-04-28 2021-03-16 浙江浙能技术研究院有限公司 面向高端燃煤发电装备非平稳运行的高斯非高斯特征协同解析与监测方法
CN112700389A (zh) * 2021-01-13 2021-04-23 安徽工业大学 一种活性污泥微生物彩色显微图像去噪方法
CN113687421A (zh) * 2021-08-23 2021-11-23 中国石油大学(北京) 地震信号的数据处理方法、装置、电子设备及存储介质

Also Published As

Publication number Publication date
CN108710150B (zh) 2019-09-06

Similar Documents

Publication Publication Date Title
CN108710150B (zh) 一种基于稳健奇异谱分析的地震不规则噪声去除方法
Qiu et al. Deep learning prior model for unsupervised seismic data random noise attenuation
Liu et al. A 1D time-varying median filter for seismic random, spike-like noise elimination
Zhou et al. Seismic noise attenuation using an online subspace tracking algorithm
Lari et al. Adaptive singular spectrum analysis for seismic denoising and interpolation
CN106443787B (zh) 叠前地震道集压噪方法及其装置
Van Veen An analysis of several partially adaptive beamformer designs
CN106680876B (zh) 一种地震数据联合去噪方法
CN108828670B (zh) 一种地震数据降噪方法
CN112578471B (zh) 一种探地雷达杂波噪声去除方法
CN110068865B (zh) 一种低秩矩阵逼近的沙漠地震噪声压制方法
Zhou et al. A hybrid method for noise suppression using variational mode decomposition and singular spectrum analysis
Jumah et al. Dimensionality‐reduced estimation of primaries by sparse inversion
Bai et al. Seismic signal enhancement based on the low‐rank methods
CN113204051A (zh) 一种基于变分模态分解的低秩张量地震数据去噪方法
Hald Denoising of cross-spectral matrices using canonical coherence
Zhang et al. Seismic signal enhancement and noise suppression using structure-adaptive nonlinear complex diffusion
Ma et al. Multiple attenuation with 3D high-order high-resolution parabolic Radon transform using lower frequency constraints
Li et al. Seismic noise suppression using weighted nuclear norm minimization method
Liu et al. Trace-wise coherent noise suppression via a self-supervised blind-trace deep-learning scheme
Bahia et al. Widely linear denoising of multicomponent seismic data
Wang et al. Denoising with weak signal preservation by group-sparsity transform learning
Liu et al. A self‐supervised scheme for ground roll suppression
Yang et al. Deep learning with soft attention mechanism for small-scale ground roll attenuation
Chen et al. Robust reduced-rank seismic denoising

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200616

Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee after: CNOOC RESEARCH INSTITUTE Co.,Ltd.

Patentee after: CHINA NATIONAL OFFSHORE OIL Corp.

Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee before: CNOOC RESEARCH INSTITUTE Co.,Ltd.

Patentee before: CHINA NATIONAL OFFSHORE OIL Corp.

Co-patentee before: Chengdu University of Technology