一种混合的地层吸收补偿方法

摘要

本发明涉及地球物理勘探技术的地震资料处理中混合的地层吸收补偿方法。把当前道接收到的地震波场信号从时间域变换到频率域，从所处的频率网格映射到新的频率网格中，进行与频率相关的比例缩放，变换到时间域得到校正的地层吸收所造成的地震波频散效应的波场信号，采用时间域增益限制滤波器进行循环卷积得到新的波场信号。本发明在频率域进行快速的频散校正，在时间域进行稳定的振幅补偿，保证了吸收补偿过程的稳定高效，提高了地震资料分辨率及薄地层的解释准确性。

主权项

一种混合的地层吸收补偿方法，其特征是采用以下具体实施步骤：1)采集地震数据，按照炮集或共中心点道集记录方式排列；2)用傅立叶变换方法把当前道接收到的地震波场信号从时间域变换到频率域；3)把步骤2)得到的频率域波场从所处的频率网格映射到新的频率网格中；所述的映射关系为： Ω = 2 πf ( f f peak ) γ , γ = 1 πQ 上式中：Q是新的频率网格变量，f是当前所处的频率网格变量，fpeak是当前地震数据的主频，γ是与地层吸收因子有关的变量，Q是表示当前地层的品质因子，由测量得到；4)把步骤3)得到的新频率网格下的波场进行与频率相关的比例缩放；5)把步骤4)缩放得到的新频率网格下的波场从频率域用傅立叶变换方法变换到时间域得到校正的地层吸收所造成的地震波频散效应的波场信号；6)采用时间域增益限制滤波器为： g ( t ) = IFFT ( f , f < f clip f clip · cos 2 ( f - f clip f nyq - f clip ) , f > f clip ) , f clip = ln 10 · Q πt · G lim 20 式中：IFFT()表示对括号里的表达式用傅立叶方法进行反变换，f为频率，fclip是由增益限制求取的截止频率，fnyq为奈奎斯特频率，Glim为增益限制其值为30分贝，t为当前层底的时间，Q为当前地层的品质因子；7)用步骤5)得到的波场信号与步骤6)得到的增益限制滤波器进行循环卷积得到新的波场信号，是当前层吸收补偿后的结果；8)把步骤2)得到的地震波场向下延拓到下一层的顶部；9)重复步骤3)到步骤8)进行当前层的地层吸收补偿，直到所有层都补偿完成；10)重复步骤2)到步骤9)进行下一个地震道的地层吸收补偿，直到所有地震道都补偿完成。