发明名称 一种基于局部方差分析的自适应曲波阈值地震去噪方法
摘要 本发明提供了一种基于局部方差分析的自适应曲波阈值地震去噪方法,该方法在深入分析曲波系数在不同尺度、不同角度下特征的基础上,提出了一种自适应阈值技术实现每个角度下有效信号与噪声的分离,并且借助曲波域的局部方差信息增强对非噪声成分的保护,从而减小对有效信号(特别是弱信号)的伤害。在合成地震记录与陆上地震数据的实践证明,阈值迭代法的实际去噪表现并不理想,而本发明去噪效果更加真实有效,在压制噪声与保护细节方面具有独特优势。
申请公布号 CN103713324B 申请公布日期 2015.05.13
申请号 CN201410004693.1 申请日期 2014.01.06
申请人 孙赞东 发明人 孙学凯;孙赞东;王招明;韩剑发
分类号 G01V1/36(2006.01)I 主分类号 G01V1/36(2006.01)I
代理机构 北京国林贸知识产权代理有限公司 11001 代理人 李桂玲;许文娟
主权项 一种基于局部方差分析的自适应曲波阈值地震去噪方法,其特征在于,所述方法包括以下步骤:步骤1:对地震数据进行曲波变换获得曲波变换系数;步骤2:利用最精细尺度下曲波系数的标准差来确定随机噪声在时间域的标准差,采用Monta‑Carlo分析获得曲波域各尺度及角度的变换权重,并以此计算曲波域中每一尺度及角度下的噪声标准差;步骤3:利用含噪信号在每一尺度及角度下的噪声标准差,结合局部方差分析,推算有效信号的标准差;所述的有效信号的标准差的计算公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mi>max</mi><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mi>y</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>j</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>&sigma;</mi><mi>a</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>j</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>,</mo><mn>0</mn><mo>)</mo></mrow></msqrt></mrow>]]></math><img file="FDA0000673475890000011.GIF" wi="794" he="120" /></maths>式中<img file="FDA0000673475890000012.GIF" wi="78" he="89" />代表含噪信号的总体方差,j,l分别代表尺度与角度,σ<sub>y</sub>(j,l)代表输入含噪信号的总体标准差,σ<sub>a</sub>(j,l)代表噪声标准差;步骤4:利用上述获得的噪声与有效信号的标准差建立阈值,进行软阈值处理;步骤5:进行曲波逆变换,获得时间域的去噪结果。
地址 102249 北京市昌平区府学路18号中国石油大学(北京)