一种地震资料的偏移成像处理方法
【专利摘要】本发明涉及地震资料的偏移成像处理方法，野外采集得到炮记录后，对其进行平面波编码得到若干不同射线参数的平面波记录，使用逆时偏移算子对平面波记录进行偏移成像，并在最小二乘反演框架下，通过共轭梯度法使得成像结果向反射率模型收敛，最终得到信噪比更高、振幅属性更为可靠的成像结果，本发明将地震成像看做最小二乘反演问题，能够较好地抑制逆时偏移算法中的低频噪音，减弱稀疏采集造成的震源效应，补偿深层反射同相轴能量，大大改善了常规逆时偏移的成像质量；此外通过平面波编码压缩数据克服了最小二乘偏移计算效率的缺陷，是一种有效的高精度成像方法。
【专利说明】一种地震资料的偏移成像处理方法
【技术领域】
[0001] 本发明涉及地震资料处理方法，具体是地震资料的偏移成像处理方法。
【背景技术】
[0002] 近年来，基于双程波动方程的逆时偏移（RTM)方法在地震成像领域吸引了越 来越多地球物理学家的关注。逆时偏移方法相对于实际油气勘探过程中被广泛应用的 Kirchhoff偏移方法和单程波波动方程偏移方法，具有原理简单、成像精确且无倾角限制、 适应复杂速度模型等优点。但是RTM方法在成像时会产生成像噪音，在一定程度上湮没了 地下构造。虽然国内外一些研究学者应用Laplace滤波或者其改进方法能够一定程度上压 制RTM成像产生的低频噪音，但是成像结果的保幅性及地下振幅的均衡性得不到保证。这 些弱点使得逆时偏移难以满足高精度成像以及岩性成像的需要。
[0003] 而基于反演思想的最小-乘偏移（LSM)是保幅成像的最佳方法。但是，将最小_. 乘偏移看做数据空间的反演问题，其计算量约为常规偏移的2Xniter倍（niter为迭代次 数）。尤其是当目标区域较大时，LSM的计算成本无法接受。
【发明内容】
[0004] 本发明目的在于提供一种地震资料的偏移成像处理方法，为叠前平面波最小二乘 逆时偏移（P-LSRTM)方法，将地震成像看做最小二乘反演问题，通过共轭梯度法使成像结 果向真实反射率模型收敛，大大改善了常规RTM的成像质量；通过平面波编码压缩偏移数 据体，提高了计算效率，实现了高精度地震成像的方法。
[0005] 本发明包括以下步骤：
[0006] (1)采用叠加前单炮数据，输入初始反射率模型、速度场和平面波编码函数；
[0007] (2)对炮记录进行平面波编码：
【权利要求】
1. 一种地震资料的偏移成像处理方法，其特征在于，该地震资料的偏移成像处理方法 包括以下步骤： (1) 采用采用叠加前单炮数据输入初始反射率模型、速度场和平面波编码函数； (2) 对炮记录进行平面波编码： (3) 通过反偏移算子计算当前反射率模型下的平面波记录， (4) 计算当前反射率模型下不同射线参数的平面波记录与真实记录的残差； (5) 通过逆时偏移算子反传残差计算更新梯度 (6) 计算迭代更新的步长： (7) 更新反射率模型 (8) 重复步骤3到7,直到当前的反射率模型满足精度； (9) 叠加不同射线参数平面波记录的成像结果，输出当前模型。
【文档编号】G01V1/28GK104360381SQ201410587640
【公开日】2015年2月18日 申请日期:2014年10月20日 优先权日:2014年10月20日
【发明者】李闯, 黄建平, 段心彪, 郭书娟, 李庆洋, 李振春, 李继光, 赵胜天 申请人:李闯