[发明专利]一种噪声压制的方法、装置及介质在审
申请号: | 202211102537.X | 申请日: | 2022-09-09 |
公开(公告)号: | CN115639607A | 公开(公告)日: | 2023-01-24 |
发明(设计)人: | 黄炜霖;汪发亮;高菲;李洁丽;刘伟杰;李季东;刘科 | 申请(专利权)人: | 中国石油大学(北京) |
主分类号: | G01V1/36 | 分类号: | G01V1/36;G01V1/30;G01V1/28 |
代理公司: | 北京集佳知识产权代理有限公司 11227 | 代理人: | 丁曼曼 |
地址: | 102200*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 噪声 压制 方法 装置 介质 | ||
本申请公开了一种噪声压制的方法、装置及介质,涉及数字信号处理领域。包括:获取检波器对地震信号进行观测得到的原始观测数据;引入阻尼算子对原始观测数据进行降秩操作并获取降秩后的矩阵;根据降秩后的矩阵获取低秩滤波数据;采用滤波器对低秩滤波数据中的预设频率的信号进行滤波以便获取第一去噪数据。该方法中,在对原始数据进行降秩操作的过程中引入了阻尼算子,由于阻尼算子可以削弱由强扰动噪声和随机噪声的低秩投影引入的人为干扰,因此实现了对地震数据中的非相干噪声的压制。此外,本申请还提供一种噪声压制的装置以及计算机可读存储介质,与上述的噪声压制的方法具有相同或相对应的技术特征,效果同上。
技术领域
本申请涉及数字信号处理领域,特别是涉及一种噪声压制的方法、装置及介质。
背景技术
在记录地震数据的过程中,由于地震数据中可能会掺杂随机噪声(随机噪声属于非相干干扰波),而有的强度较大的噪声,会对地震数据的记录形成很大的干扰,导致记录的地震数据不准确,故而,需要对地震数据中的非相干噪声进行压制。
目前,通常采用奇异谱分析(Singular Spectrum Analysis,SSA)来压制地震非相干噪声。SSA方法利用地震信号的低秩特征,将噪声压制视为低秩重建问题。然而,在采用SSA方法的过程中,由于SSA中的低秩投影被迫引入显著的人为干扰,以拟合最小二乘意义下的高振幅不稳定干扰,因此,采用SSA方法不能有效地对地震数据中的非相干噪声进行压制。
由此可见,如何对地震数据中的非相干噪声进行压制是本领域人员亟需解决的技术问题。
发明内容
本申请的目的是提供一种噪声压制的方法、装置及介质,用于压制地震数据中的非相干噪声。
为解决上述技术问题,本申请提供一种噪声压制的方法,包括:
获取检波器对地震信号进行观测得到的原始观测数据;
引入阻尼算子对所述原始观测数据进行降秩操作并获取降秩后的矩阵;
根据所述降秩后的矩阵获取低秩滤波数据;
采用滤波器对所述低秩滤波数据中的预设频率的信号进行滤波以便获取第一去噪数据。
优选地,在所述采用滤波器对所述低秩滤波数据中的预设频率的信号进行滤波以便获取第一去噪数据之后,所述方法还包括:
根据所述低秩滤波数据、所述第一去噪数据获取加权矩阵;
采用所述滤波器对所述低秩滤波数据、所述加权矩阵中的所述预设频率的信号进行滤波以便获取第二去噪数据。
优选地,所述阻尼算子根据阻尼因子确定,获取所述阻尼因子包括:
获取预设的阻尼因子的上限值、下限值、最大迭代次数、当前迭代次数;
根据所述上限值、所述下限值、所述最大迭代次数、所述当前迭代次数获取所述阻尼因子。
优选地,所述引入阻尼算子对所述原始观测数据进行降秩操作并获取降秩后的矩阵包括:
对所述原始观测数据进行频率域变换并获取变换后的原始观测数据;
对所述变换后的原始观测数据重新排列并获取汉克尔化后的矩阵;
根据所述阻尼因子、截断秩对所述汉克尔化后的矩阵进行降秩操作并获取所述降秩后的矩阵。
优选地,所述根据所述降秩后的矩阵获取低秩滤波数据包括:
对所述降秩后的矩阵进行反对角线平均以便获取所述低秩滤波数据。
优选地,在所述采用所述滤波器对所述低秩滤波数据、所述加权矩阵中的所述预设频率的信号进行滤波以便获取第二去噪数据之后,所述方法还包括:
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油大学(北京),未经中国石油大学(北京)许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202211102537.X/2.html,转载请声明来源钻瓜专利网。