用于地震资料去噪的结构导向方向广义全变差正则化方法技术

本发明专利技术公开了一种用于地震资料去噪的结构自适应方向广义全变差正则化方法，包括以下步骤：(1)基于梯度结构张量，提出了一种新的地震同相轴空变倾角的计算方法；(2)将l

Generalized total variation regularization method of structural guidance direction for seismic data denoising

【技术实现步骤摘要】
用于地震资料去噪的结构导向方向广义全变差正则化方法
：本专利技术属于地震勘探
，涉及一种用于地震资料去噪的结构导向方向广义全变差(Structure-AdaptiveDirectionalTotalGeneralizedVariation,SADTGV)正则化方法。
技术介绍
：噪声压制在地震资料处理中起着非常重要的作用，是提高地震资料解释精度的首要任务。随机噪声在地震数据的时空域随机分布，覆盖了地震数据的有效频带，因此，随机噪声去除是地震数据处理中的一个重要问题。工业界常用的地震资料去噪方法有小波变换方法、中值滤波、KL变换、奇异值分解等。近年来，用于求解反问题的变分正则化方法已被应用于地震资料随机噪声压制、并取得了一定的效果。其中，全变差(TotalVariation,TV)正则化方法是一个开创性工作，其模型定义如下：式中，BV(Ω)表示定义在Ω上的有界变差函数空间，f表示含噪数据，u表示去噪后的数据，λ＞0为正则化参数。此外，TV(u)定义为式中，v表示以原点为中心的欧几里得闭单位球B2(0)本文档来自技高网...

【技术保护点】
1.用于地震资料去噪的结构导向方向广义全变差正则化方法，其特征在于包括以下步骤：/n(1)估计地震剖面在任意一点处同相轴的倾角：/n第一步，对含噪数据进行高斯滤波，然后计算其梯度，记为/n

【技术特征摘要】
1.用于地震资料去噪的结构导向方向广义全变差正则化方法，其特征在于包括以下步骤：
(1)估计地震剖面在任意一点处同相轴的倾角：
第一步，对含噪数据进行高斯滤波，然后计算其梯度，记为



式中，u表示含噪数据，σ表示高斯函数的标准差。对于二维数据，g的分量可以写为



第二步，将梯度g通过并矢积运算映射成如下张量



然后，将T与高斯核G(x；ρ)做卷积，得到u的梯度结构张量



第三步，计算地震数据在某点处的梯度结构张量Tρ的特征值与特征向量，而相应于最大特征值的特征向量为该点处反射同相轴的法向量，将地震同相轴的倾角限制到区间[0,π)内，在任意一点处的倾角定义为：



式中，e表示在x点处反射同相轴的法向量，e1、e2分别表示沿x1和x2方向的单位向量，<·,·>表示向量的内积运算。
(2)构建l2-SADTGV正则化去噪模型
将



中方向广义全变差中...

【专利技术属性】
技术研发人员：王德华，高静怀，丁小丽，张丽丽，周宏安，
申请(专利权)人：西安工业大学，
类型：发明
国别省市：陕西;61