用于地震勘探的F-XY域改进非局部均值去噪方法及设备技术

本发明专利技术公开了用于地震勘探的F

【技术实现步骤摘要】
用于地震勘探的F
‑
XY域改进非局部均值去噪方法及设备


[0001]本专利技术涉及地震数据去噪领域，具体涉及用于地震勘探的F
‑
XY域改进非局部均值去噪 方法及设备。

技术介绍

[0002]目前，地震记录中包含各种噪音，严重影响后续地震资料处理，因此在常规的子波处理， 速度分析之前需要进行叠前噪音消除；线性切除是最简单，直接的去噪方式之一，对相干噪 音的消除效果好，但存在时间域去噪不可逆的特性，并且切除的数据量不易约束。小波变换 分频可以逐点对相干干扰方向进行扫描，对直线变换的干扰有较好的去噪效果，由于小波变 换分频只能在有限的频带进行去噪，其计算原理是在数据上直接做减法，会导致剖面上能量 不均，因此用小波变换分频去噪还需要在去噪后进行其他修饰处理。
[0003]地震资料是一个庞大的多维数据体，在去噪中不仅要考虑去噪方法的新颖性，同时还需 要考虑庞大地震数据去噪的计算效率问题，一套良好的去噪处理方法，需要充分考虑计算的 时间成本，去噪常常在更高级的地震方法处理之前进行，是地震数据处理的基础；可惜的是， 目前去噪算法并没有很好的解决去噪效果和时间成本之间关系，实际应用中，更多的是在效 果和运算复杂度之间求得一个平衡；现阶段去噪方法往往存在以下几个问题：去噪效果好的 方法往往需要花费巨大时间成本，如何降低方法的计算时间是亟待解决的问题；去噪方法一 般是整个频率段的整体去噪，很难针对每个频率进行分别去噪；去噪中的去噪系数往往都是 人工选择的，没有跟随数据的信噪比进行针对性的调整去噪强度，导致去噪效果不佳。

技术实现思路

[0004]本专利技术所要解决的技术问题是地震资料去噪时间较长，去噪效果不明显，平滑因子与真 实数据间关系不匹配等技术问题，目的在于提供用于地震勘探的F
‑
XY域改进非局部均值去 噪方法及设备，通过高斯噪声数据求平均，解决了地震资料去噪时间较长，去噪效果不明显， 平滑因子与真实数据间关系不匹配的问题。
[0005]本专利技术通过下述技术方案实现：
[0006]用于地震勘探的F
‑
XY域改进非局部均值去噪方法及设备，包括以下步骤：
[0007]S1：将采集的地震时间域T
‑
XY数据沿时间方向转换为频率域F
‑
XY数据，并进行三维体 频率域切片预处理，得到频率域切片；其中X可以代表炮号或者CMP号，Y可以代表道号或 者炮检距，T代表地震数据的时间深度；
[0008]S2：系统对频率域切片进行改进NLM去噪处理，获得所有频率域切片的去噪结果；
[0009]S3：完成对待处理的所有频率域的去噪结果后，系统进行Fourier反变换获得时间域T
‑
XY 数据体，即去噪后的数据体。
[0010]利用基于F
‑
XY域的改进非局部均值去噪计算流程不仅加快了空间滤波算法的计算速度， 同时还提高了滤波效果。一方面频率域切片数据能有效的保证改进的NLM算法有
G(l)表示具有标准差为σ的高斯核，l代 表邻域内nl个元素之一；其准确表达式为：其中，(x0,y0) 代表高斯核的中心，(x,y)代表所计算的l点的位置坐标；通过遍历区域D，利用标准差计算 公式，求取标准差σ，从而可以获得高斯核G(l)，最后可以得到D(i,j)。进一步，对权重系 数中的区域平滑因子h进行计算，对区域D中全局扫描的时，需要求取区域内的小窗能量和 平均能量v
ij
，v
mean
，N
D
为区域D中扫描时所有可分移动小窗口总数，从而可求取区域内的 信噪再对求取的区域内SNR进行归一化处理，得到归一化后的 在对SNR进行自适应计算，从而得到每个区域的平滑因子 其中，C为平滑阈值约束参数取50，至此获得窗口的权重系数。
[0025]进一步，所述设备包括：
[0026]一个或多个处理器；
[0027]存储器，用于存储一个或多个程序，
[0028]当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器执 行如权利要求1
‑
7任一所述的用于地震勘探的F
‑
XY域的改进非局部均值去噪方法。
[0029]本专利技术与现有技术相比，具有如下的优点和有益效果：
[0030]本专利技术一种用于地震勘探的F
‑
XY域改进非局部均值去噪方法及设备，能够有效地提高地 震数据去噪的效率和效果。改进NLM计算方法需要遍历所有数据体，针对大规模三维数据时， 计算效率非常低，采用频率域F
‑
XY的切片数据能有效的降低每次循环计算的数据量，并且能 针对特定频率域切片进行NLM去噪，从而达到快速高效的获得去噪结果；其次，改进的NLM 计算方法通过区域的高斯加权欧氏距离代替欧式距离能有效的提高去噪效果，且加入的自适 应平滑参数算法能对高信噪比的数值窗口数据选择低去噪强度，而对低信噪比的窗口数据选 择高去噪强度。
附图说明
[0031]此处所说明的附图用来提供对本专利技术实施例的进一步理解，构成本申请的一部分，并不 构成对本专利技术实施例的限定。在附图中：
[0032]图1是本专利技术计算过程流程图；
[0033]图2是去噪效果对比图；
[0034]图3是NLM去噪与改进后的NLM去噪效果对比图；
[0035]图4是本专利技术整体去噪效果图。
具体实施方式
[0036]为使本专利技术的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本专利技术 作进一步的详细说明，本专利技术的示意性实施方式及其说明仅用于解释本专利技术，并不作为对本 专利技术的限定。
[0037]在以下描述中，为了提供对本专利技术的透彻理解阐述了大量特定细节。然而，对于本领域 普通技术人员显而易见的是：不必采用这些特定细节来实行本专利技术。在其他实例中，为了避 免混淆本专利技术，未具体描述公知的结构、电路、材料或方法。
[0038]在整个说明书中，对“一个实施例”、“实施例”、“一个示例”或“示例”的提及意味着： 结合该实施例或示例描述的特定特征、结构或特性被包含在本专利技术至少一个实施例中。因此， 在整个说明书的各个地方出现的短语“一个实施例”、“实施例”、“一个示例”或“示例”不 一定都指同一实施例或示例。此外，可以以任何适当的组合和、或子组合将特定的特征、结 构或特性组合在一个或多个实施例或示例中。此外，本领域普通技术人员应当理解，在此提 供的示图都是为了说明的目的，并且示图不一定是按比例绘制的。这里使用的术语“和/或
”ꢀ
包括一个或多个相关列出的项目的任何和所有组合。
[0039]在本专利技术的描述中，需要理解的是，术语“前”、“后”、“左”、“右”、“上”、
ꢀ“
下”、“竖直”、“水平”、“高”、“低”“内”、“外”等指示的方位或位置关系为 基于附图所示的方位或位置关系，仅是为了便于描述本专利技术和本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.用于地震勘探的F
‑
XY域改进非局部均值去噪方法，其特征在于，包括以下步骤：S1：将采集的地震时间域T
‑
XY数据沿时间方向转换为频率域F
‑
XY数据，并进行三维体频率域切片预处理，得到二维频率域切片；S2：系统对频率域切片进行改进NLM去噪处理，获得所有频率域切片的去噪结果；S3：完成对待处理的所有频率域的去噪结果后，系统进行Fourier反变换获得时间域T
‑
XY数据体，即去噪后的数据体。2.根据权利要求1所述的用于地震勘探的F
‑
XY域改进非局部均值去噪方法，其特征在于，所述步骤S1的具体过程包括：S11：系统采集到地震时间域T
‑
XY数据，对时间域T
‑
XY数据沿时间方向进行Fourier变换，得到频率域二维F
‑
XY数据；S12：对于得到的频率域F
‑
XY数据，设置滤波频段和进行遍历提取频率域切片操作，得到频率域切片f
i
及f
i
实部地震体数据和虚部地震体数据。3.根据权利要求1所述的用于地震勘探的F
‑
XY域改进非局部均值去噪方法，其特征在于，所述步骤S2的具体过程包括：S21：对频率域切片f
i
的实部地震体数据和虚部地震体数据进行改进NLM去噪处理，得到去噪处理后的频率域切片；S22：选择去噪后频率域切片，将频率域切片分为N个D大小的窗口，在D窗口中设置一个大小为d*d的方形搜索窗口和一个大小为d*d的方形中心窗口；S23：设置D窗口的中心窗口大小为d*d，中心点为y；设置搜索窗口大小为d*d，中心点为x，所述搜索窗口y(d*d)在D窗口中进行移动并与中心窗口x(d*d)进行相似程度计算，得到计算后的搜索窗口和中心窗口的距离数据；S24：根据得到的距离数据，系统计算区域D中的区域信噪比SNR，得到区域信噪比SNR数据；S25：系统对区域信噪比SNR数据进行归一化处理后，得到归一化后的SNR数据；S26：系统根据得到归一化后的SNR数据对SNR进行自适应计算，得到每个区域的平滑因子h；S27：系统重复进行步骤S21
‑
S26所需要计算的频率域切片，得到所有频率域切片的去噪结果。4.根据权利要求1所述的用于地震勘探的F
‑
XY域改进非局部均值去噪方法，其特征在于，所述系统遍历所有N个窗口的数据之后，得到频率域切片去噪后的频率域F
‑
XY域数...

【专利技术属性】
技术研发人员：刘力辉，
申请(专利权)人：成都晶石石油科技有限公司，
类型：发明
国别省市：