专利名称：一种混合的地层吸收补偿方法
技术领域：
本发明涉及地球物理勘探技术，是一种地震资料处理中混合的地层吸收补偿方
法。
背景技术：
在地球物理勘探领域中，当激发的地震波在地下介质中传播时，地层的Q吸收效应可以使地震波的能量衰减，并且频率越高衰减越快，致使接收到的地震信号频带变窄，分辨率降低；同时由于吸收造成的地震波频散效应使地震波的相位发生畸变，导致地震剖面与测井数据不能有效匹配，给解释人员带来不准确的信息。因此补偿地层的吸收效应对提高地震资料分辨率、薄地层的解释准确性及更好的井_地数据匹配有很大的意义。目前采用的Q吸收补偿方法存在两个大的问题稳定性和效率。时间域的Q吸收补偿方法是通过串联反Q滤波来实现，这种方法对于比较长的地震道来说，其稳定性不能得到有效保证，方法有很大的局限性。频率域的Q吸收补偿方法通过波场延拓来实现，通过在延拓算子中引入地层吸收因子来实现对振幅和频散补偿，这种方法需要逐个数据样点进行计算，对于大的三维叠前地震数据来说，其效率不能满足生产需要。

发明内容
本发明目的在于提供一种稳定性好和效率提高的混合的地层吸收补偿方法。本发明采用以下技术方案实现，具体实施步骤是1)采集地震数据，按照炮集或共中心点道集记录方式排列；2)用傅立叶变换方法把当前道接收到的地震波场信号从时间域变换到频率域；3)把步骤2)得到的频率域波场从所处的频率网格映射到新的频率网格中；所述的映射关系为
权利要求
1.一种混合的地层吸收补偿方法，其特征是采用以下具体实施步骤1)采集地震数据，按照炮集或共中心点道集记录方式排列；2)用傅立叶变换方法把当前道接收到的地震波场信号从时间域变换到频率域；3)把步骤2)得到的频率域波场从所处的频率网格映射到新的频率网格中； 所述的映射关系为Jpeak7tQ上式中Q是新的频率网格变量，f是当前所处的频率网格变量，fpeak是当前地震数据的主频，Y是与地层吸收因子有关的变量，Q是表示当前地层的品质因子，由测量得到；4)把步骤3)得到的新频率网格下的波场进行与频率相关的比例缩放；5)把步骤4)缩放得到的新频率网格下的波场从频率域用傅立叶变换方法变换到时间域得到校正的地层吸收所造成的地震波频散效应的波场信号；6)采用时间域增益限制滤波器为/'/ < fclip8{t)~-IFFT(\ “-^i-LJf-U > “W. ifJ nyq J clip式中IFFT()表示对括号里的表达式用傅立叶方法进行反变换，f为频率，f。lip是由增益限制求取的截止频率，f_为奈奎斯特频率，Glim为增益限制其值为30分贝，t为当前层底的时间，Q为当前地层的品质因子；7)用步骤5)得到的波场信号与步骤6)得到的增益限制滤波器进行循环卷积得到新的波场信号，是当前层吸收补偿后的结果；8)把步骤2)得到的地震波场向下延拓到下一层的顶部；9)重复步骤3)到步骤8)进行当前层的地层吸收补偿，直到所有层都补偿完成；10)重复步骤2)到步骤9)进行下一个地震道的地层吸收补偿，直到所有地震道都补偿完成。
2.根据权利要求1所述的方法，特征是步骤4)所述的频率相关的比例缩放的比例缩放因子为丄(丄)、"丄fpeak^Q上式中f是旧的频率网格变量，fpeak是当前地震数据的主频，Y是与地层吸收因子有关的变量，Q是表示当前地层的品质因子。
3.根据权利要求1所述的方法，特征是步骤7)所述的循环卷积计算为州)=x(t) + b(t)[x(t) * g(i)] + ^b(tf[x{t) * g(t) * +L L其中冲上式中x(t)是步骤5)得到的波场信号， b(t)是与地层吸收因子有关的算子， Q是表示当前地层的品质因子， g(t)是步骤6)得到的增益限制滤波器， y(t)是输出信号。
4.根据权利要求1所述的方法，特征是步骤8)所述的波场向下延拓因子为
全文摘要
本发明涉及地球物理勘探技术的地震资料处理中混合的地层吸收补偿方法。把当前道接收到的地震波场信号从时间域变换到频率域，从所处的频率网格映射到新的频率网格中，进行与频率相关的比例缩放，变换到时间域得到校正的地层吸收所造成的地震波频散效应的波场信号，采用时间域增益限制滤波器进行循环卷积得到新的波场信号。本发明在频率域进行快速的频散校正，在时间域进行稳定的振幅补偿，保证了吸收补偿过程的稳定高效，提高了地震资料分辨率及薄地层的解释准确性。
文档编号G01V1/36GK102269822SQ20101019706
公开日2011年12月7日 申请日期2010年6月2日 优先权日2010年6月2日
发明者李合群, 王文闯 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司