基于小波变换与奇异值分解的地震资料去噪新方法.pdfVIP

石油天然气学报 (江汉石油学院学报) 2010年2月 第32卷 第i期 JournalofOilandGasTechnology(J．JPI) Feb．2010 Vo1．32 No．1 基于小波变换与奇异值分解的地震资料去噪新方法 王小品，贺振华 (成都理工大学油气藏地质与开发工程国家重点实验室，四)lI成都61o059) 熊8尧军 (成都理工大学地球探测与信息技术教育部重点实验室，四JJl成都610059) [摘要]为 了有效地去除地震资料中的随机噪声，充分利用小波变换 (wT)去噪和奇异值分解 (SVD) 去噪方法的优点，提 出了一种新 的基于小波变换和奇异值分解 (WT-SVD) 的地震资料去噪方法。该方 法首先进行小波软 阚值去噪，有效地 降低噪声的方差；然后进行基于倾角扫描 的奇异值分解去噪，识别 噪声点， 自动追踪 同相轴 ，并进行同相轴拉平处理，充分利用 了奇异值分解方法处理水平 同相轴噪声效 果好 的优点。理论模型和实际资料 的去噪结果表 明，该研究提 出的wT—SVD方法简单易行，比单一的 SVD方法和wT方法的去噪效果更显著，有效地消除了地震资料 中的随机噪声，显著地提高了地震资料 的信噪 比。 [关键词]小波变换；奇异值分解 ；倾角扫描 ；地震资料去噪 [中图分类号]P631．44 [文献标识码]A [文章编号]1000—9752 (2O1。】01—0232—05 常规的基于小波变换 (WT)的去噪方法采用多尺度分解较好地保留了信号的细节信息，可以有效 地消除地震资料中的随机噪声 ，对于精确地震勘探十分重要 ，但是该方法的应用难点在于重建小波系数 时的阈值的选取 ，当阈值选取过大时，有效信号细节被去除过多；当阈值选取过小时，达不到去噪的 目 的。虽然 目前国内外的学者提出了很多阈值的选取方法，如 Donoho和Johnston提出的软硬阈值法 ， sure阈值等 ，但是这些阈值方法只能给出粗略阈值，处理具体信号时必须依靠经验给出精细阈值 ]。 另外 ，常规的基于奇异值分解 (SVD)的去噪方法_4]，从矩阵的角度出发，将包含信号信息的矩阵分解 到一系列奇异值和奇异值矢量对应的信号子空间中，然后通过去除奇异值然后用有效奇异值重构矩阵达 到去噪的目的，其可以有效地恢复受噪声污染的水平地震同相轴 。对于受噪声污染的弯曲同相轴 ，如斜 坡、穹隆或透镜体等，其在去噪的同时会极大地损伤有效信号。即使采用局部同相轴拉平技术 ]，由于 在很小的时窗内少量水平同相轴很难体现地震道相关性强的优势 ，SVD方法的去噪效果也会大打折扣。 为了有效地去除地震资料中的随机噪声 ，笔者充分利用了wT方法和 SVD方法去噪方法的优点，提 出 了一种新的基于小波变换和奇异值分解 (wT—SVD)的地震资料去噪方法。该方法首先采用常规 的小 波软阈值去噪方法 进行第 1步去噪处理，其可以有效地消除部分随机噪声，并降低噪声的方差，显著 增强地震道有效信号的相关性m；然后采用改进的SVD方法——基于倾角扫描 ]的 SVD方法进行第 2 步去噪处理。提出的基于倾角扫描的SVD方法通过倾角扫描确定有效同相轴的方向，其方向 (同相轴 斜率)只取整数提高了此方法的可操作性 ，并根据奇异值曲线的波动有效地进行纯噪声点的识别，有效 地克服了常规的局部同相轴拉平技术错误地把一些噪声点当作有效信号点的缺陷，可以有效地消除地震 资料中的随机噪声。 1 方法原理 1．1 小波软阈值去噪原理 设 ，(￡)为有限能量的信号或平方可积信号，则其连续小波变换可定义为： [收稿日期]2009—12—20 [基金项 目]国家自然科学基金项目 。 [作者简介]王小品 (1983一)，男，2006年大学毕业，硕士生，现主要从事物探资料处理新方法方面的研究工作。 第 32卷第 1期 王小品等：基于小波变换与奇异值分解的地震资料去噪新方法 7．Uf 6)一 一 ． ( )dz (1) 式中，。为尺度参数；6为位移因子；函数 ，。( )为子波；r(n，6)为小波系数；￡表示时间。然后在细 节的尺度上设置阈值去噪即可，即： 一 sgn(w~ ’ — l，当 f f