本申请公开了一种各向异性扩散滤波方法、装置、电子设备及介质。该方法可以包括：针对叠后地震数据提取归一化的相干值；根据叠后地震数据计算特征值与特征向量；根据特征值、特征向量与相干值，构建扩散张量；通过一致性增强扩散方法进行扩散，重复上述步骤直至达到设定的扩散时间，对地震图像进行平滑。本发明专利技术通过使用相干属性对各向异性扩散滤波的过程加以控制，在去噪的同时保持并增强对边界的刻画。在去噪的同时保持并增强对边界的刻画。在去噪的同时保持并增强对边界的刻画。
【技术实现步骤摘要】
各向异性扩散滤波方法、装置、电子设备及介质
[0001]本专利技术涉及地球物理勘探
，更具体地，涉及一种各向异性扩散滤波方法、装置、电子设备及介质。
技术介绍
[0002]在地震数据的采集过程中，由于野外的各种干扰噪声会在地震记录中产生影响，在地震数据预处理过程中并不能将噪声全部压制掉，所以导致地震数据的信噪比低，所获得的偏移叠加剖面存在噪声影响，剖面中的构造信息不清晰，不利于之后的地震资料解释工作。
[0003]因此国内外相关的学者展开了大量的工作，提出了中值滤波、时频域滤波、Kuwahara滤波等方法，然而这些方法都没办法在压制噪声的同时保护边界。
[0004]为此，Fehmers等提出构造导向的滤波方法，杨培杰等提出方向性边界保持的断层增强方法，Dave采用了双边滤波的方法进行地震图像增强。这些方法都没有利用叠后地震剖面的相干体属性进行约束。
[0005]因此，有必要开发一种基于相干属性的各向异性扩散滤波方法、装置、电子设备及介质，可以用于对叠后地震数据剖面的噪声压制同时保持构造的边界信息，对断层等构造的刻画进行增强。
[0006]公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
[0007]本专利技术提出了一种各向异性扩散滤波方法、装置、电子设备及介质，其能够通过使用相干属性对各向异性扩散滤波的过程加以控制，在去噪的同时保持并增强对边界的刻画。
[0008]第一方面，本公开实施例提供了一种各向异性扩散滤波方法，包括：
[0009]针对叠后地震数据提取归一化的相干值；
[0010]根据所述叠后地震数据计算特征值与特征向量；
[0011]根据所述特征值、所述特征向量与所述相干值，构建扩散张量；
[0012]通过一致性增强扩散方法进行扩散，重复上述步骤直至达到设定的扩散时间，对地震图像进行平滑。
[0013]优选地，根据所述叠后地震数据计算特征值与特征向量包括：
[0014]设定高斯滤波的参数值，针对叠后地震数据进行空间域高斯核函数滤波，获得滤波数据；
[0015]根据所述滤波数据，计算梯度结构张量；
[0016]根据所述梯度结构张量进行特征值分解，计算所述特征值与所述特征向量。
[0017]优选地，通过公式(1)进行空间域高斯核函数滤波：
[0018]g＝I*G
σ
ꢀꢀꢀ
(1)
[0019]其中，I为叠后地震数据，G
σ
为尺度为σ的高斯核函数，*为卷积运算符，g为滤波后地震数据。
[0020]优选地，通过公式(2)计算梯度结构张量：
[0021][0022]其中，为梯度运算符，G
ρ
为尺度为ρ的高斯核函数，T为梯度结构张量。优选地，通过公式(3)构建扩散张量：
[0023]D＝ε(λ1v1v
1T
+λ2v2v
2T
)
ꢀꢀꢀ
(3)
[0024]其中，ε为归一化的相干值，λ1、λ2为特征值，v1、v2为特征向量，D为扩散张量。
[0025]优选地，通过公式(4)进行扩散：
[0026][0027]优选地，使用公式(4)进行非线性的各向异性扩散，将扩散方程进行差分离散，给定时间步长和总扩散时间进行扩散。
[0028]作为本公开实施例的一种具体实现方式，
[0029]第二方面，本公开实施例还提供了一种各向异性扩散滤波装置，包括：
[0030]相干值提取模块，针对叠后地震数据提取归一化的相干值；
[0031]计算模块，根据所述叠后地震数据计算特征值与特征向量；
[0032]扩散张量构建模块，根据所述特征值、所述特征向量与所述相干值，构建扩散张量；
[0033]扩散模块，通过一致性增强扩散方法进行扩散，重复上述步骤直至达到设定的扩散时间，对地震图像进行平滑。
[0034]优选地，根据所述叠后地震数据计算特征值与特征向量包括：
[0035]设定高斯滤波的参数值，针对叠后地震数据进行空间域高斯核函数滤波，获得滤波数据；
[0036]根据所述滤波数据，计算梯度结构张量；
[0037]根据所述梯度结构张量进行特征值分解，计算所述特征值与所述特征向量。
[0038]优选地，通过公式(1)进行空间域高斯核函数滤波：
[0039]g＝I*G
σ
ꢀꢀꢀ
(1)
[0040]其中，I为叠后地震数据，G
σ
为尺度为σ的高斯核函数，*为卷积运算符，g为滤波后地震数据。
[0041]优选地，通过公式(2)计算梯度结构张量：
[0042][0043]其中，为梯度运算符，G
ρ
为尺度为ρ的高斯核函数，T为梯度结构张量。优选地，通过公式(3)构建扩散张量：
[0044]D＝ε(λ1v1v
1T
+λ2v2v
2T
)
ꢀꢀꢀ
(3)
[0045]其中，ε为归一化的相干值，λ1、λ2为特征值，v1、v2为特征向量，D为扩散张量。
[0046]优选地，通过公式(4)进行扩散：
[0047][0048]优选地，使用公式(4)进行非线性的各向异性扩散，将扩散方程进行差分离散，给定时间步长和总扩散时间进行扩散。
[0049]第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
[0050]存储器，存储有可执行指令；
[0051]处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的各向异性扩散滤波方法。
[0052]第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的各向异性扩散滤波方法。
[0053]其有益效果在于：本专利技术采用梯度结构张量的方法获得了方向信息，以此方向信息计算特征值特征向量构建出扩散张量，考虑到了叠后剖面的方向信息，保护了构造边界；将叠后地震数据的相干值计算出来作为横向不连续因子加以控制扩散张量，更好的考虑到了断层等构造的特征信息，增强了对断层等构造边界的刻画。
[0054]本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。
附图说明
[0055]通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。
[0056]图1示出了根据本专利技术的一个实施例的各向异性扩散滤本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种各向异性扩散滤波方法，其特征在于，包括：针对叠后地震数据提取归一化的相干值；根据所述叠后地震数据计算特征值与特征向量；根据所述特征值、所述特征向量与所述相干值，构建扩散张量；通过一致性增强扩散方法进行扩散，重复上述步骤直至达到设定的扩散时间，对地震图像进行平滑。2.根据权利要求1所述的各向异性扩散滤波方法，其中，根据所述叠后地震数据计算特征值与特征向量包括：设定高斯滤波的参数值，针对叠后地震数据进行空间域高斯核函数滤波，获得滤波数据；根据所述滤波数据，计算梯度结构张量；根据所述梯度结构张量进行特征值分解，计算所述特征值与所述特征向量。3.根据权利要求2所述的各向异性扩散滤波方法，其中，通过公式(1)进行空间域高斯核函数滤波：g＝I*G
σ
ꢀꢀꢀꢀꢀ
(1)其中，I为叠后地震数据，G
σ
为尺度为σ的高斯核函数，*为卷积运算符，g为滤波后地震数据。4.根据权利要求2所述的各向异性扩散滤波方法，其中，通过公式(2)计算梯度结构张量：其中，为梯度运算符，G
ρ
为尺度为ρ的高斯核函数，T为梯度结构张量。5.根据权利要求1所述的各向异性扩散滤波方法，其中，通过公式(3)构建扩散张量：D＝ε(λ1v1v
1T
+λ2...
【专利技术属性】
技术研发人员：王鸣谦，许璐，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。