一种地震资料成像处理方法技术

本发明专利技术提供了一种地震资料成像处理方法，获取地震资料后，首先根据地震资料的主频和频宽选择参考频率与最高频率，进而根据探区速度场大小确定束宽，从而可以给定初始射线参数间隔；然后基于地表倾角信息在震源和接收点处直接出射高斯束，并进行运动学射线追踪和动力学射线追踪，以求取高斯束的旅行时和振幅信息；最后，根据高斯束的旅行时和振幅信息，应用反褶积成像条件计算对应单道成像值，进而对每道循环求得单炮成像值，将所有炮的成像值进行叠加便获得最终的成像剖面。

【技术实现步骤摘要】

本专利技术属于地震资料处理领域，具体为
技术介绍
研宄复杂地表和复杂地下地质体双重复杂条件下的处理技术是地震勘探领域的 热点。复杂地表对地震数据的采集与处理都带来了不小的困难，尤其是在中国西部和南方 山前带探区尤为突出。在近地表速度远小于地下介质速度且地形变化较平缓的条件下，简 单的静校正可以满足要求；但当地表高程起伏较大且近地表速度变化剧烈时，静校正的假 设条件不成立，采用简单的静校正对波场造成的畸变会影响偏移成像的质量。为此，发展双 复杂条件下的地震偏移成像方法成为解决上述问题的关键。 目前，针对起伏地表条件下的地震成像问题，国内外学者提出并发展了一系列有 效的处理方法：第一类是基于波动方程的偏移成像方法。Berryhill (1979,1984)较早提 出了波动方程基准面校正的方法。随后，Yilmaz和Lucas (1986)，Schneider和Pillip等 (1995)，Bevc (1997)以及Yang等（1999)在此基础上结合层替代思想进一步发展了此方 法。虽然波动方程基准面校正能够有效地消除起伏地表对同相轴造成的扭曲。但是，近地 表速度的横向变化会严重影响基准面校正的精度。针对上述缺陷，Reshef (1991)，Beasley 和 Lynn (1992)，何英等（2002)，Shragge (2005)，程玖兵（2006)，叶月明等（2008)以及吕 彬等（2011)提出了多种改进的基于波动理论的偏移成像方法来解决复杂地表问题。然 而，在复杂的地表条件下，波场延拓很难适应不规则的观测系统，并且运算量较大。第二类 是基于射线的偏移成像方法。射线类方法具有更加灵活，易于实现，计算效率高等优点。 Wiggins (1984)提出了适用于起伏地表条件的Kirchhoff偏移方法，Gray等（1995)实现了 直接在起伏地表进行Kirchhoff偏移的方法且成像效果较好，Jager等（2003)提出了起伏 地表条件下的保幅Kirchhoff偏移方法，而董春晖等（2009)及刘国峰等（2010)先后基于 Kirchhoff算子提出了起伏地表直接时间偏移的流程与应用。Gray (2005)发展了一种适用 于起伏地表的局部静校正的束偏移方法，但是当地表高程及速度变化剧烈时，静校正方法 会造成地震波场发生畸变从而影响成像质量，Yue等（2010)提出了一种适用于起伏地表条 件的保幅高斯束偏移方法，该方法通过考虑起伏地表平面波的传播特点，直接在起伏地表 进行局部平面波的分解。然而，由于束中心出射位置与接收点位置之间有水平距离差和高 程差，需要引入相位校正因子以近似表征格林函数，当水平距离差较大时，这种近似所带来 的振幅误差不能忽略，尤其当地形起伏较大且近地表速度变化剧烈时，上述近似处理会极 大地降低偏移成像的质量。
技术实现思路
本专利技术的目的在于针对中国西部和南方山前带探区地震资料成像处理的难题，提 出了，是一种双复杂条件下基于单道处理的保幅高斯束深度域 成像方法，以纵波地震资料和纵波速度作为输入，利用地表倾角信息，通过在每个接收点处 出射高斯束来表征反向延拓的地震波场，并采用反褶积成像条件求取成像值。 本专利技术包括以下步骤； (1)根据地震资料的主频和频带宽度选择参考频率以及最高频率ω h; (2)计算高斯束初始宽度：【主权项】1. ，其特征在于，该地震资料成像处理方法包括以下步 骤： (1) 根据地震资料的主频和频带宽度选择参考频率以及最高频率； (2) 计算高斯束初始宽度； (3) 设置初始射线参数间隔； (4) 根据地表高程、倾角以及纵横波速度信息，结合上述1、2、3步骤所选取并计算的参 数，利用单道处理方式，在震源和接收点处直接出射高斯束，并通过求解运动学射线追踪方 程和动力学射线追踪方程计算得到高斯束的旅行时和振幅信息，然后求得复杂地表条件下 的反向延拓地震波场； (5) 根据基于反褶积成像条件推导出的复杂地表高斯束偏移成像公式，利用步骤4所 求得的高斯束旅行时和振幅信息计算对应单道的成像值； (6) 重复4、5步骤的计算流程，计算出单炮成像值，最终将所有炮的成像值进行叠加得 到所需要的成像剖面。【专利摘要】本专利技术提供了，获取地震资料后，首先根据地震资料的主频和频宽选择参考频率与最高频率，进而根据探区速度场大小确定束宽，从而可以给定初始射线参数间隔；然后基于地表倾角信息在震源和接收点处直接出射高斯束，并进行运动学射线追踪和动力学射线追踪，以求取高斯束的旅行时和振幅信息；最后，根据高斯束的旅行时和振幅信息，应用反褶积成像条件计算对应单道成像值，进而对每道循环求得单炮成像值，将所有炮的成像值进行叠加便获得最终的成像剖面。【IPC分类】G01V1-28【公开号】CN104808242【申请号】CN201410587072【专利技术人】袁茂林, 黄建平, 李继光, 赵胜天, 杨继东, 李振春, 段心彪, 郭书娟 【申请人】袁茂林【公开日】2015年7月29日【申请日】2014年10月20日本文档来自技高网...

【技术保护点】
一种地震资料成像处理方法，其特征在于，该地震资料成像处理方法包括以下步骤：(1)根据地震资料的主频和频带宽度选择参考频率以及最高频率；(2)计算高斯束初始宽度；(3)设置初始射线参数间隔；(4)根据地表高程、倾角以及纵横波速度信息，结合上述1、2、3步骤所选取并计算的参数，利用单道处理方式，在震源和接收点处直接出射高斯束，并通过求解运动学射线追踪方程和动力学射线追踪方程计算得到高斯束的旅行时和振幅信息，然后求得复杂地表条件下的反向延拓地震波场；(5)根据基于反褶积成像条件推导出的复杂地表高斯束偏移成像公式，利用步骤4所求得的高斯束旅行时和振幅信息计算对应单道的成像值；(6)重复4、5步骤的计算流程，计算出单炮成像值，最终将所有炮的成像值进行叠加得到所需要的成像剖面。

【技术特征摘要】

【专利技术属性】
技术研发人员：袁茂林，黄建平，李继光，赵胜天，杨继东，李振春，段心彪，郭书娟，
申请(专利权)人：袁茂林，
类型：发明
国别省市：山东;37