一种高效的地震数据随机噪声衰减方法
【专利摘要】本发明是一种高效的地震数据随机噪声衰减方法，分别给定各方向的预测算子长度，各空间方向计算窗长度和时间窗大小，将数据分成小计算窗，在当前计算窗内沿时间方向做Fourier变换，计算三维预测算子序列，与频率空间域原始数据进行褶积，得到该计算窗随机噪声衰减结果，计算得到权系数，加权得到该计算窗随机噪声衰减加权结果，对当前窗以及前一计算窗进行均值处理，做反变换，得到当前计算窗时间－空间域随机噪声衰减结果，循环直至最末计算窗结束。本发明保持了现有技术去噪效果的显著性，比现有技术效率提升60%以上。
【专利说明】-种高效的地震数据随机噪声衰减方法
【技术领域】
[0001] 本发明涉及地球物理勘探地震数据处理技术范畴，是一种高效的地震数据随机噪 声衰减方法。
【背景技术】
[0002] 地震数据在采集过程中，不可避免的受到各种噪声的干扰，致使有效信号失真，提 高地震数据的信噪比是地震勘探中一项重要而基础的工作。频率-空间域预测随机噪声衰 减技术是压制H维叠前地震数据随机噪声的一种行之有效的处理方法，对于提高地震数据 的信噪比有着显著的效果。但是，现有的方法计算量巨大，处理效率低下，处理周期很长，严 重限制了该技术在石油物探实际生产中的推广应用。
[0003] 现有频率-空间域预测随机噪声衰减技术的一大弊端是对地震资料各方向做一 致处理。众所周知，H维叠前地震资料是各向异性的，空间各方向反射波复杂程度不同，各 方向间可能存在巨大差异。现有技术忽视了H维叠前地震资料的该一重要特征，没有对不 同方向区分处理，在某些空间方向计算量浪费严重，导致该一随机噪声衰减技术难W推广。 与此同时，由于各方向一致处理的限制，现有技术难W处理空间分布不均的地震数据，例如 宽线数据等。
【发明内容】
[0004] 本发明目地在于提供一种提高地震数据信噪比、提升处理效率、节约处理成本的 一种高效的地震数据随机噪声衰减方法。
[0005] 本发明采用W下步骤实现：
[0006] 1)采集H维叠前地震数据；
[0007] 2)根据地震数据空间各方向反射波复杂程度的不同，分别给定各方向的预测算子 长度，根据预测算子长度给定各空间方向计算窗长度，并给定时间窗大小；
[0008] 所述的预测算子长度如果某方向反射波较复杂，则可适当增加该个方向预测算 子长度；如果某方向反射波较平缓，则可W适当减少该方向预测算子长度。
[0009] 所述的预测算子长度为5 - 8。
[0010] 所述的给定各空间方向计算窗长度选取各空间预测算子长度的四倍加一。
[0011] 时间窗为500毫砂一个计算时间窗。
[0012] 3)根据各方向计算窗大小将数据分成小计算窗；
[0013] 4)在当前计算窗内沿时间方向做化urier变换得到频率空间域原始数据
[0014] {su，k;f}，其中i，j，k分别表示X、Y、0立个空间方向，f表示频率；
[0015] 5)用下式计算H维预测算子序列:
[0016]
【权利要求】
1. 一种高效的地震数据随机噪声衰减方法，特点是采用以下步骤实现： 1) 采集三维叠前地震数据； 2) 根据地震数据空间各方向反射波复杂程度的不同，分别给定各方向的预测算子长 度，根据预测算子长度给定各空间方向计算窗长度，并给定时间窗大小； 3) 根据各方向计算窗大小将数据分成小计算窗； 4) 在当前计算窗内沿时间方向做Fourier变换得到频率空间域原始数据{smf}，其 中i，j，k分别表示Χ、Υ、0三个空间方向，f表示频率； 5) 用下式计算三维预测算子序列Ix1^J :
其中，l，m，n;p，q，r为自然数，变化范围分别为1，"％LX;1，?，LY ;1，···，!!)； 1^、1^、11)分别是乂、￥、0方向预测算子长度； 6) 利用预测算子序列Ix1^J与频率空间域原始数据{Si^k;f}进行褶积，得到该计算窗 随机噪声衰减结果; 7) 利用下式计算得到权系数qf :
8) 将计算窗随机噪声衰减结果进行加权得到该计算窗随机噪声衰减加权结 果WW 9) 对当前窗以及前一计算窗进行均值处理，或者线性插值； 10) 在该窗内做反Fourier变换，得到当前计算窗时间一空间域随机噪声衰减结果； 11) 当前计算窗如果是最末计算窗则完成随机噪声衰减；如不是最末计算窗则跳至下 一计算窗，并返回步骤4)循环，直至最末计算窗结束。
2. 根据权利要求1的方法，特点是步骤1)所述的预测算子长度如果某方向反射波较复 杂，则可适当增加这个方向预测算子长度；如果某方向反射波较平缓，则可以适当减少该 方向预测算子长度。
3. 根据权利要求1的方法，特点是步骤1)所述的预测算子长度为5 - 8。
4. 根据权利要求1的方法，特点是步骤1)所述的给定各空间方向计算窗长度选取各空 间预测算子长度的四倍加一。
5. 根据权利要求1的方法，特点是步骤1)时间窗为500毫秒一个计算时间窗。
【文档编号】G01V1/36GK104375184SQ201310356433
【公开日】2015年2月25日 申请日期:2013年8月15日 优先权日:2013年8月15日
【发明者】贾继伟, 周兴元, 罗国安 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司