申请/专利权人：成都晶石石油科技有限公司

申请日：2021-01-21

公开（公告）日：2024-04-16

公开（公告）号：CN112907458B

主分类号：G06T5/70

分类号：G06T5/70;G06F17/14;G01V1/28

优先权：

专利状态码：有效-授权

法律状态：2024.04.16#授权;2021.06.22#实质审查的生效;2021.06.04#公开

摘要：本发明公开了用于地震勘探的F‑XY域改进非局部均值去噪方法及设备，采集的地震时间域T‑XY数据转换为F‑XY数据，并进行三维体频率域切片预处理，得到二维频率切片；系统进行频率域切片改进NLM去噪处理，获得所有切片的去噪结果；完成对待处理的所有切片的去噪结果后，系统进行Fourier反变换获得T‑XY数据体，利用基于F‑XY域的改进非局部均值去噪计算流程不仅加快了空间滤波算法的计算速度，同时还提高了滤波效果；一方面频率切片数据能有效的保证改进的NLM算法有较高的计算效率，另一方面，采用改进的NLM滤波算法通过高斯加权欧氏距离计算相邻窗口的相关性。

主权项：1.用于地震勘探的F-XY域改进非局部均值去噪方法，其特征在于，包括以下步骤：S1：将采集的地震时间域T-XY数据沿时间方向转换为频率域F-XY数据，并进行三维体频率域切片预处理，得到二维频率域切片；S2：系统对频率域切片进行改进NLM去噪处理，获得所有频率域切片的去噪结果；其中，所述步骤S2的具体过程包括：S21：对频率域切片fi的实部地震体数据和虚部地震体数据进行改进NLM去噪处理，得到去噪处理后的频率域切片；S22：选择去噪后频率域切片，将频率域切片分为N个D大小的窗口，在D窗口中设置一个大小为d*d的方形搜索窗口和一个大小为d*d的方形中心窗口；S23：设置D窗口的中心窗口大小为d*d，中心点为y；设置搜索窗口大小为d*d，中心点为x，所述搜索窗口yd*d在D窗口中进行移动并与中心窗口xd*d进行相似程度计算，得到计算后的搜索窗口和中心窗口的距离数据；S24：根据得到的距离数据，系统计算区域D中的区域信噪比SNR，得到区域信噪比SNR数据；S25：系统对区域信噪比SNR数据进行归一化处理后，得到归一化后的SNR数据；S26：系统根据得到归一化后的SNR数据对SNR进行自适应计算，得到每个区域的平滑因子h；S27：系统重复进行步骤S21-S26所需要计算的频率域切片，得到所有频率域切片的去噪结果；对频率域切片fi的实部和虚部进行改进NLM算法去噪，所述去噪过程通过公式表示；其中，ω，v分别代表地震数据x和y之间的相似度，采用x，y为中心的邻域距离决定Di，j，所述Di，j为i和j两个区域的高斯加权欧氏距离；选择频率域切片fi，开两个窗口进行数据处理，其中一个窗口是大窗D，其含义为将数据分为N个D大小的窗口，设置一个小的搜索窗口d*d和一个中心窗口d*d，两个窗口的中心分别为y和x，搜索窗口yd*d在D窗口中进行移动，并与中心窗口xd*d进行相似程度的计算其权重计算方法为：其中，vNi,vNj表示矩形邻域，计算两个邻域的二范数，h为平滑因子，Zx为归一化系数；所述对Di,j的计算采用相近两个区域的高斯加权欧氏距离计算，其表达式：Gl表示具有标准差为σ的高斯核，l代表邻域内nl个元素之一；其准确表达式为：其中，x0,y0代表高斯核的中心，x,y代表所计算的l点的位置坐标；通过遍历区域D，利用标准差计算公式，求取标准差σ，从而可以获得高斯核Gl，最后可以得到Di,j；得到平滑因子h的具体步骤包括：对区域D中全局扫描的时，求取区域内的小窗能量和平均能量vij，vmean，ND为区域D中扫描时所有可分移动小窗口总数，从而求取区域内的信噪比再对求取的区域内SNR进行归一化处理，得到归一化后的在对SNR进行自适应计算，从而得到每个区域的平滑因子其中，C为平滑阈值约束参数取50；S3：完成对待处理的所有频率域的去噪结果后，系统进行Fourier反变换获得时间域T-XY数据体，即去噪后的数据体。

全文数据：

权利要求：

百度查询： 成都晶石石油科技有限公司 用于地震勘探的F-XY域改进非局部均值去噪方法及设备

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD