各向异性扩散滤波方法、装置、电子设备及介质.pdf

(19)国家知识产权局(12)发明专利申请(10)申请公布号CN115993652A(43)申请公布日2023.04.21(21)申请号202111217473.3(22)申请日2021.10.19(71)申请人中国石油化工股份有限公司地址100027北京市朝阳区朝阳门北大街22号申请人中国石油化工股份有限公司石油物探技术研究院(72)发明人王鸣谦许璐(74)专利代理机构北京思创毕升专利事务所11218专利代理师孙向民廉莉莉(51)Int.Cl.G01V1/36(2006.01)G01V1/30(2006.01)权利要求书2页说明书8页附图3页(54)发明名称各向异性扩散滤波方法、装置、电子设备及介质(57)摘要本申请公开了一种各向异性扩散滤波方法、装置、电子设备及介质。该方法可以包括：针对叠后地震数据提取归一化的相干值；根据叠后地震数据计算特征值与特征向量；根据特征值、特征向量与相干值，构建扩散张量；通过一致性增强扩散方法进行扩散，重复上述步骤直至达到设定的扩散时间，对地震图像进行平滑。本发明通过使用相干属性对各向异性扩散滤波的过程加以控制，在去噪的同时保持并增强对边界的刻画。CN115993652ACN115993652A权利要求书1/2页1.一种各向异性扩散滤波方法，其特征在于，包括：针对叠后地震数据提取归一化的相干值；根据所述叠后地震数据计算特征值与特征向量；根据所述特征值、所述特征向量与所述相干值，构建扩散张量；通过一致性增强扩散方法进行扩散，重复上述步骤直至达到设定的扩散时间，对地震图像进行平滑。2.根据权利要求1所述的各向异性扩散滤波方法，其中，根据所述叠后地震数据计算特征值与特征向量包括：设定高斯滤波的参数值，针对叠后地震数据进行空间域高斯核函数滤波，获得滤波数据；根据所述滤波数据，计算梯度结构张量；根据所述梯度结构张量进行特征值分解，计算所述特征值与所述特征向量。3.根据权利要求2所述的各向异性扩散滤波方法，其中，通过公式(1)进行空间域高斯核函数滤波：g＝I*Gσ(1)其中，I为叠后地震数据，Gσ为尺度为σ的高斯核函数，*为卷积运算符，g为滤波后地震数据。4.根据权利要求2所述的各向异性扩散滤波方法，其中，通过公式(2)计算梯度结构张量：其中，为梯度运算符，Gρ为尺度为ρ的高斯核函数，T为梯度结构张量。5.根据权利要求1所述的各向异性扩散滤波方法，其中，通过公式(3)构建扩散张量：TTD＝ε(λ1v1v1+λ2v2v2)(3)其中，ε为归一化的相干值，λ1、λ2为特征值，v1、v2为特征向量，D为扩散张量。6.根据权利要求1所述的各向异性扩散滤波方法，其中，通过公式(4)进行扩散：7.根据权利要求6所述的各向异性扩散滤波方法，其中，使用公式(4)进行非线性的各向异性扩散，将扩散方程进行差分离散，给定时间步长和总扩散时间进行扩散。8.一种各向异性扩散滤波装置，其特征在于，包括：相干值提取模块，针对叠后地震数据提取归一化的相干值；计算模块，根据所述叠后地震数据计算特征值与特征向量；扩散张量构建模块，根据所述特征值、所述特征向量与所述相干值，构建扩散张量；扩散模块，通过一致性增强扩散方法进行扩散，重复上述步骤直至达到设定的扩散时间，对地震图像进行平滑。9.一种电子设备，其特征在于，所述电子设备包括：存储器，存储有可执行指令；处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1‑7中任一项所述的各向异性扩散滤波方法。2CN115993652A权利要求书2/2页10.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1‑7中任一项所述的各向异性扩散滤波方法。3CN115993652A说明书1/8页各向异性扩散滤波方法、装置、电子设备及介质技术领域[0001]本发明涉及地球物理勘探技术领域，更具体地，涉及一种各向异性扩散滤波方法、装置、电子设备及介质。背景技术[0002]在地震数据的采集过程中，由于野外的各种干扰噪声会在地震记录中产生影响，在地震数据预处理过程中并不能将噪声全部压制掉，所以导致地震数据的信噪比低，所获得的偏移叠加剖面存在噪声影响，剖面中的构造信息不清晰，不利于之后的地震资料解释工作。[0003]因此国内外相关的学者展开了大量的工作，提出了中值滤波、时频域滤波、Kuwahara滤波等方法，然而这些方法都没办法在压制噪声的同时保护边界。[0004]为此，Fehmers等提出构造导向的滤波方法，杨培杰等提出方向性边界保持的断层增强方法，Dave采用了双边滤波的方法进行地震图像增强。这些方法都没有利用叠后地震剖面的相干体属性进行约束。[0005]因此，有必要开发一种基于相干属性的各向异性扩散滤波方法、装置、电子设备及介质，可以用于对叠