高阶拉普拉斯算子逆时偏移低频噪声去除方法.pdfVIP

第52卷第 6期 石 油 物 探 Vo1．52，No．6 2013年 11月 GEOPHYSICALPR0SPECTING FOR PETROLEUM NOV．，2013 文章编号：1000—1441(2013)06—0642—08 高阶拉普拉斯算子逆时偏移低频噪声去除方法 郭念民 ，冯雪梅 ，李海山 (1．中国石油天然气集团公司塔里木油田分公司，新疆库尔勒841000；2中国石油天然气股份有 限公司勘探开发研究院西北分院，甘肃兰州 730020) 摘要：叠前逆时偏移被认为是 目前最为准确的成像方法 ，但是因为采用了双程波动方程和时间一致性成像准则， 逆时偏移成像会产生很强的低频噪声。逆时偏移的成像噪声具有区别于有效波的低频特性、方向特性和大角度 特性，这些特性是逆时偏移低频噪声去除方法的实现基础。拉普拉斯算子滤波去噪方法可以直接在偏移后的数 据上进行，实现起来比较简单 ，增加的计算量小，在实际中应用较多；但是常规的拉普拉斯算子滤波去噪不彻底， 而且改变了成像数据的振幅和相位信息。为了改进逆时偏移低频噪声去除效果，分析了拉普拉斯算子去噪的基 本原理，研究了改进的高阶拉普拉斯算子去噪方法，修正了常规方法中存在的上述问题，既保证了去噪的效果 ， 又保持了成像数据振幅和相位信息的相对不变。模型测试结果证明了高阶拉普拉斯算子逆时偏移低频噪声去 除方法的可行性和有效性。 关键词 ：逆时偏移；低频噪声；拉普拉斯算子；去噪方法 DOI：lO．3969／j．issn．1000—1441．2013．06．013 中图分类号 ：P631．4 文献标识码：A 逆时偏移在波动方程偏移理论发展的初期就 反射界面的位置和形态，然后据此设计了阻尼因子 已经提出，1983年Whitmore[、Baysal等[、Mc— 系数 。 Mechan[3J和 Loewenthal等_l4用不同的方法实现 2006年，Guitton等_8采用最小二乘预测滤波 了逆时偏移 。但是，由于受计算和存储条件的限 法实现了逆时偏移成像的去噪，该方法有较好的保 制，逆时偏移在 20世纪八九十年代发展缓慢。21 真性，在去噪的同时不会改变数据的相位和振幅强 世纪以来，复杂地区地震勘探的要求不断提高，加 度。2007年，Liu等E9_基于单程波成像条件提出了 上大规模海量计算技术的发展，使得逆时偏移被重 一 种定性的逆时偏移成像条件，即：将波场分解为 新认识并得到了迅速发展，并逐步被应用于实际生 沿单一方向传播的波场分量，经由不同方向的波场 产 中。 分量间的有效相关组合，消除常规的成像条件所产 逆时偏移成像有其 自身的优势，不对方程分 生的低频噪声。2009年，Suh等l】0]在波场分解成 解，无倾角限制，可以实现复杂构造高角度成像、回 像去噪的基础上做了进一步改进，在 F-K域把波 转波成像、反转构造成像；但同样也有 自身缺陷，相 场分为上、下行波场和左、右行波场，然后应用扇形 比其它偏移方法，逆时偏移会产生很强的低频噪 滤波器来实现剩余噪声的压制。Zhang等E ]提 声。随着逆时偏移成像研究的发展，许多学者提出 出了角度域滤波的方法消除逆时偏移低频噪声，首 了多种去除逆时偏移低频噪声的实现方法 。2001 先输出角度域成像道集，然后直接切除大角度道集 年Youn等l5提出了拉普拉斯算子滤波去噪方法， 压制噪声。 可以较好地去除成像噪声，而且是对偏移后的数据 在这些去噪方法中，角度域滤波方法需要输出 进行操作，效率较高。2004年，Yoon等_6]提出了 角度域共成像点道集，必然会增加大量的计算量； 一 种新的PoyntingVector成像条件，在传统互相 波场分解方法由于对波场方向分