一种地震数据叠前成像方法
【技术领域】
[0001] 本发明涉及石油地质勘探地震成像领域，尤其涉及到一种地震数据叠前成像方 法。
【背景技术】
[0002] 随着油气勘探程度的不断加深，地震勘探所面临的问题大大增加，新老勘探区块 无论是复杂地表或是近地表速度纵横向剧烈变化还是地下复杂构造都给地球物理工作者 提出了极富挑战性的难题。各种复杂地表和复杂地下构造的地区，如沙漠、山地、山前、黄土 塬、山前地带的逆冲带、南方碳酸岩裸露区、逆掩推覆体构造带、地层褶皱强烈带、断裂异常 复杂和横向非均质引起横向速度变化剧烈的地区，给地震资料处理带来非常大的难度，常 规的地震资料处理方法和思路很难使这些地区地震资料获得精确成像效果。
【发明内容】
[0003] 本发明的目的是提供一种，以解决现有技术的上述不足。
[0004] 本发明的目的是通过以下技术方案来实现：
[0005] 本发明提供了一种地震数据叠前成像方法，包括以下步骤：
[0006] 利用双程波方程对震源波场进行正向外推并保存外推结果；
[0007] 利用逆时双程波方程对接收波场进行反向外推；
[0008] 每反向外推一步就用成像条件进行求和，获得局部成像数据；
[0009] 将所有炮集的逆时偏移结果进行叠加，得到叠前深度偏移成像。
[0010] 优选的，还包括；将双程波方程转化为边界条件式：
[0011]
【主权项】
1. 一种地震数据叠前成像方法，其特征在于，包括以下步骤： 利用双程波方程对震源波场进行正向外推并保存外推结果； 利用逆时双程波方程对接收波场进行反向外推； 每反向外推一步就用成像条件进行求和，获得局部成像数据； 将所有炮集的逆时偏移结果进行叠加，得到叠前深度偏移成像。
2. 根据权利要求1所述的地震数据叠前成像方法，其特征在于，还包括；将双程波方程 转化为边界条件式：
其中，U为地下点的波场矢量；f(x，t)为地表接收到的矢量波；t为时间；x，z分别为水 平距离和垂直深度，z=0表示地表。
3. 根据权利要求2所述的地震数据叠前成像方法，其特征在于，双程波方程中的二维 各向同性介质中弹性波动方程表达式为：
其中，vx，VZ为质点振动速度的水平分量和垂直分量；τ XX、τ ζζ、τ XZ为应力的三个 分量；α β分别为介质中的纵波速度与横波速度；Q为介质密度；X，ζ分别为空间坐标的两 个方向；t为时间。
4. 根据权利要求3所述的地震数据叠前成像方法，其特征在于，还包括：对得到的二维 各向同性介质中弹性波动方程表达式在交错网格空间中采用高阶差分法对式进行差分离 散，可得到二维各向同性介质中弹性波方程逆时延拓的2N阶精度差分格式，其中，N为任意 正整数。
5. 根据权利要求4所述的地震数据叠前成像方法，其特征在于，所述逆时双程波方程 的二维各向同性介质中中程函方程表达式为；
其中，t为时间；s为慢度；X, z为空间坐标当地震源位置。
【专利摘要】本发明涉及一种地震数据叠前成像方法，包括以下步骤：利用双程波方程对震源波场进行正向外推并保存外推结果；利用逆时双程波方程对接收波场进行反向外推；每反向外推一步就用成像条件进行求和，获得局部成像数据；将所有炮集的逆时偏移结果进行叠加，得到叠前深度偏移成像。本发明的有益效果为：解决了振幅问题：保幅效果好，无振幅和相位近似，通过精确的照明补偿和真振幅成像，使岩性预测成为可能；解决了成像精度问题：能对陡倾角成像，回转波成像，成像效果好；解决了多次波成像问题：多次波收敛、聚焦。
【IPC分类】G01V1-28
【公开号】CN104635260
【申请号】CN201310547814
【发明人】孙传文, 曲鸣飞
【申请人】北京瑞拉爱堡地质勘探技术有限公司
【公开日】2015年5月20日
【申请日】2013年11月8日