一种基于曲波变换和全变差的联合衰减随机噪声处理方法
【专利摘要】本发明涉及地震勘探【技术领域】的一种基于曲波变换和全变差的联合衰减随机噪声处理方法，其先获得单炮地震数据，然后将获取的单炮数据或单炮数据的叠加进行曲波变换，再进行多尺度曲波阈值去噪，然后采用全变差去噪技术进行去噪，最后通过图形显示软件输出地震数据。该方法根据随机噪声在曲波域的分布规律来选择最佳阈值，使资料的信噪比达到最高，使用曲波变换得到最优去噪效果的目的；通过全变差最小化技术，调整曲波系数，克服了单独使用曲波变换带来的伪曲线缺点，使显示的地层数据更真实、可靠，以便进一步进行地层分析，获取更为准确的含油、含矿的分析结果。
【专利说明】一种基于曲波变换和全变差的联合衰减随机噪声处理方法
【技术领域】
[0001] 本发明涉及一种地震勘探【技术领域】，特别涉及地震数据处理中随机噪声的压制方 法。
【背景技术】
[0002] 曲波变换（Curvelet变换）是一种较新的多尺度几何变换算法。1999年，Cantos 和Donoho在Ridgelet变换的基础上提出了连续曲波（Curvelet)变换，即第一代Curve let 变换；2002年，Cantos等人提出了第二代Curvelet变换；2005年，Cantos等人提出了两种 基于第二代Curvelet变换理论的快速离散实现方法：1)非均匀空间抽样的二维FFT算法 (Unequally-Spaced Fast Fourier Transform, USFFT) ;2) Wrapping 算法（Wrapping-Based Transform)。由于曲波对频率和方向具有优越的局部性而被广泛的应用于地震数据处理领 域，其中，2004年Herrmann F等最先将Curvelet变换应用到地震数据处理领域。
[0003] 在地震勘探中，常规的随机噪声衰减方法对有效信号的损害比较大，如RPF (径 向预测去噪)，AMC0D(倾角相干叠加），多项式拟合去噪，中值滤波，SVD方法以及工业界普 遍使用的RNA去噪方法等，尽管RNA在诸多方法中算得上是比较好的随机噪声衰减方法， 但此方法对有效信号也有一定的损害，为了较好的去除随机噪声，Neelamani等人在2008 年引入了一种多尺度的变换方法--曲波变换来衰减随机噪声，取得了较好的效果。尽管 Curvelet变换在去噪方面存在诸多优势，但是Curvelet变换也存在固有的缺陷--去噪 时，容易在地震剖面中产生强的能量团，同时在同相轴边缘产生不光滑现象。
【发明内容】
[0004] 本发明的目的在于提供一种基于曲波变换和全变差的联合衰减随机噪声处理方 法，该方法可以有效压制影响叠前单炮数据、叠后数据品质的随机噪声，最大限度的减少对 有效信号的损害，是一种有效和高效的随机噪声衰减技术。
[0005] 为实现上述技术目的，本发明可采用如下技术方案：一种基于曲波变换和全变差 的联合衰减随机噪声处理方法，依次包括如下步骤： 1)地震数据的获取： 首先获取单炮数据：布置好激发点位置和多个接收点位置，然后通过在激发点埋置的 炸药爆炸激发产生地震波，各接收点安放的检波器接收到地下反射界面反射上来的地震 波，所有检波器接收到的同一炮激发的地震波数据，形成单炮数据；
【权利要求】
1. 一种基于曲波变换和全变差的联合衰减随机噪声处理方法，其特征在于依次包括如 下步骤： 1) 地震数据的获取： 首先获取单炮数据：布置好激发点位置和多个接收点位置，然后通过在激发点埋置的 炸药爆炸激发产生地震波，各接收点安放的检波器接收到地下反射界面反射上来的地震 波，所有检波器接收到的同一炮激发的地震波数据，形成单炮数据； 2) 将获取的单炮数据或单炮数据的叠加进行曲波变换 对叠加数据或单炮数据/进行曲波变换得到曲波变换系数其中可 以表示为：
次数，是则结束迭代，否则继续迭代，直到结束迭代； 5)输出地震数据：通过图形显示软件来输出地震数据。
2. 根据权利要求1所述的一种基于曲波变换和全变差的联合衰减随机噪声处理方法， 其特征在于步骤2)中叠加数据的获取方法为：对获取的多个单炮数据分别进行面波和异 常振幅衰减、地表一致性能量补偿、地表一致性反褶积、地表一致性剩余静校正，然后分选 到共中心点数据域，进行速度解释，用解释的速度进行动校正，对动校正后的数据在共中心 点域内进行叠加，获取叠加数据。
3. 根据权利要求1或2所述的一种基于曲波变换和全变差的联合衰减随机噪声处理 方法，其特征在于步骤5)中，图形显示软件包括Omega处理软件中的图形显示工具、Cgg处 理软件中的图形显示工具，Promax处理软件中的图形显示工具，Focus处理软件中的图形 显示工具，Geoeast处理软件中的图形显示工具，Landmark解释软件中的图形显示工具， Geoframe解释软件中的图形显示工具，Matlab中的wigb程序中的一种。
【文档编号】G01V1/36GK104062683SQ201410107764
【公开日】2014年9月24日 申请日期:2014年3月21日 优先权日:2014年3月21日
【发明者】薛永安, 王勇, 王山岭, 陈习峰, 庞全康, 陆树勤, 刘立民, 管文华, 潘成磊, 付波, 陈丹 申请人:中国石油化工股份有限公司, 中国石油化工股份有限公司江苏油田分公司