[发明专利]地震资料中的随机噪音消除方法及系统在审
申请号: | 201810934965.6 | 申请日: | 2018-08-16 |
公开(公告)号: | CN110837121A | 公开(公告)日: | 2020-02-25 |
发明(设计)人: | 宋林;李洪建 | 申请(专利权)人: | 中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院 |
主分类号: | G01V1/36 | 分类号: | G01V1/36 |
代理公司: | 北京思创毕升专利事务所 11218 | 代理人: | 孙向民;廉莉莉 |
地址: | 100728 北*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 地震 资料 中的 随机 噪音 消除 方法 系统 | ||
1.一种地震资料中的随机噪音消除方法,其特征在于,包括:
针对单道地震信号进行连续小波变换,获得小波系数矩阵;
针对所述小波系数矩阵进行SVD分解,获得多个奇异值;
根据所述小波系数矩阵,计算奇异熵;
根据所述奇异熵筛选重构奇异值,对所述重构奇异值进行SVD逆分解,进而获得重构后的小波系数矩阵;
针对所述重构后的小波系数矩阵进行小波逆变换,得到去噪后的地震信号。
2.根据权利要求1所述的地震资料中的随机噪音消除方法,其中,根据所述小波系数矩阵,计算奇异熵包括:
根据所述小波系数矩阵,获得l×l维的对角参数矩阵Λl×l;
根据所述对角参数矩阵Λl×l的主对角线元素,计算奇异谱;
根据所述奇异谱,计算奇异熵在阶次i处的增量;
根据所述奇异熵在阶次i处的增量,计算所述奇异熵。
3.根据权利要求2所述的地震资料中的随机噪音消除方法,其中,通过公式(1)获得l×l维的对角参数矩阵Λl×l:
其中,Am×n为小波系数矩阵,Um×l、Vn×l均为参数矩阵,T为矩阵的转置。
4.根据权利要求2所述的地震资料中的随机噪音消除方法,其中,所述奇异谱为:
其中,σh为第h个奇异谱参量,λj表示Λl×l的主对角线元素。
5.根据权利要求2所述的地震资料中的随机噪音消除方法,其中,所述奇异熵在阶次i处的增量为:
其中,ΔEi为奇异熵在阶次i处的增量,λj表示Λl×l的主对角线元素。
6.根据权利要求2所述的地震资料中的随机噪音消除方法,其中,所述奇异熵为:
其中,Ek为奇异熵,ΔEi为奇异熵在阶次i处的增量。
7.一种地震资料中的随机噪音消除系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
针对单道地震信号进行连续小波变换,获得小波系数矩阵;
针对所述小波系数矩阵进行SVD分解,获得多个奇异值;
根据所述小波系数矩阵,计算奇异熵;
根据所述奇异熵筛选重构奇异值,对所述重构奇异值进行SVD逆分解,进而获得重构后的小波系数矩阵;
针对所述重构后的小波系数矩阵进行小波逆变换,得到去噪后的地震信号。
8.根据权利要求7所述的地震资料中的随机噪音消除系统,其中,根据所述小波系数矩阵,计算奇异熵包括:
根据所述小波系数矩阵,获得l×l维的对角参数矩阵Λl×l;
根据所述对角参数矩阵Λl×l的主对角线元素,计算奇异谱;
根据所述奇异谱,计算奇异熵在阶次i处的增量;
根据所述奇异熵在阶次i处的增量,计算所述奇异熵。
9.根据权利要求8所述的地震资料中的随机噪音消除系统,其中,通过公式(1)获得l×l维的对角参数矩阵Λl×l:
其中,Am×n为小波系数矩阵,Um×l、Vn×l均为参数矩阵,T为矩阵的转置。
10.根据权利要求8所述的地震资料中的随机噪音消除系统,其中,所述奇异熵在阶次i处的增量为:
其中,ΔEi为奇异熵在阶次i处的增量,λj表示Λl×l的主对角线元素;
所述奇异熵为:
其中,Ek为奇异熵。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院,未经中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810934965.6/1.html,转载请声明来源钻瓜专利网。