论文部分内容阅读
[摘要]随着勘探的逐步深入,勘探目标越来越复杂,对地震成像也提出了越来越高的要求,采用波动方程积分法的克希霍夫偏移的成像精度难以满足地质要求,逆时偏移以其精确的算法和良好的模型测试效果,越来越多的在生产中得到应用。逆时偏移基于双程波波动方程进行波场逆时外推,并应用成像条件提取成像值,可以得到逆时偏移数据体。通过在潜山勘探中的应用,逆时偏移对潜山内幕的成像明显好于克希霍夫偏移,成像信噪比有了很大提高,反映的地质信息更加丰富。
[关键词]逆时偏移 潜山成像 双程波
[中图分类号] P62 [文献码] B [文章编号] 1000-405X(2014)-11-117-1
目前,常用的叠前深度偏移成像方法主要有单程波波动方程偏移和克希霍夫积分偏移两种。单程波波动方程偏移基于双向波方程的单向波分解,其模型假设是水平层状地层或横向变化不太剧烈的地层。始于20世纪80年代的基于双程波方程的逆时偏移技术,受计算机计算能力的限制,因其计算量巨大而发展较为缓慢。近年来,随着计算机技术的快速发展和计算能力的大幅提高,逆时偏移愈来愈受到人们的重视。
1逆时偏移原理
逆时偏移是双程波动方程偏移,它可以对多次波、回转波等对陡倾角及反转构造进行准确成像,充分利用了地震采集的地震信号。逆时偏移是通过对双程波动方程的二阶导数进行差分离散,以差分代替微分,求解波动方程。通过有限差分法来构建合适的波场传播算子,实现波场外推。逆时偏移实现过程包括基于双程波方程逆时波场外推和应用成像条件两个步骤。对于单炮道集的成像,首先对震源波场利用双程波方程进行正向外推,并保存所有时间步的外推波场。然后对接收波场利用双程波方程进行逆时外推,在时间上每逆时外推一步之后应用成像条件,得到该时刻的成像结果。所有时刻的成像结果叠加,得到该炮集的成像结果。所有炮集的逆时偏移结果叠加即可得到最终叠前深度偏移成像结果。
在三维或二维正演模拟及逆时深度偏移中,当利用截断误差为 O=(△x2,△y2,△z2,△t2)的差分格式时,为保证频散较小及递推过程稳定,差分网格要求取得非常小,这样计算需要的计算机内存及运算时间会大大增加。Dablain利用高阶差分方模拟和偏移过程。利用高阶差分方程时,网格值可以取得大些,而计算精度并不降低。这里认为截断误差高于四阶的差分方程为高阶差分方程。三维声波方程的高阶差分方程可以用统一的方式推导出来,即
其中:u=(x,y,z,t)为地表记录的压力波场; v=(x,y,z)为纵横向可变的介质速度。为导出方程(1)的离散差分格式,需把观测对应的地下介质分布区域或要对其进行地震波模拟的模型区域离散化,即把它们剖分成一个个的小方块。
令
为导出高阶差分方程,把波场进行Taylor展开,得到截断误差为O=(△xM△yM,△zM,△t2),对应的三维逆时深度偏移高阶差分程为
2逆时偏移效果分析
在辽河凹陷西部凹陷北部高升潜山,逆时偏移处理取得较好的成像效果,潜山成像精度提高,能量聚焦,潜山内幕信噪比、分辨率明显提高(图1)。与克希霍夫深度偏移相比,逆时偏移在陡坡部位、低信噪比部位成像都有了提高,反映出来的地质信息更丰富,对潜山顶部小断层的刻画更清楚,潜山内幕的形态也更清楚。
3 小结
目前,逆时偏移技术在大面积的生产当中应用较少,多数是以试验性为目的。本文充分发挥逆时偏移技术在复杂构造条件下的成像优势,成功的实现了潜山的准确成像。由于目前针对逆时偏移的速度建模技术还不成熟,我们建模的手段仍然是基于克希霍夫偏移的,所以逆时偏移技术仍有很大的发挥空间。
参考文献
[1] 底青云,王妙月
[关键词]逆时偏移 潜山成像 双程波
[中图分类号] P62 [文献码] B [文章编号] 1000-405X(2014)-11-117-1
目前,常用的叠前深度偏移成像方法主要有单程波波动方程偏移和克希霍夫积分偏移两种。单程波波动方程偏移基于双向波方程的单向波分解,其模型假设是水平层状地层或横向变化不太剧烈的地层。始于20世纪80年代的基于双程波方程的逆时偏移技术,受计算机计算能力的限制,因其计算量巨大而发展较为缓慢。近年来,随着计算机技术的快速发展和计算能力的大幅提高,逆时偏移愈来愈受到人们的重视。
1逆时偏移原理
逆时偏移是双程波动方程偏移,它可以对多次波、回转波等对陡倾角及反转构造进行准确成像,充分利用了地震采集的地震信号。逆时偏移是通过对双程波动方程的二阶导数进行差分离散,以差分代替微分,求解波动方程。通过有限差分法来构建合适的波场传播算子,实现波场外推。逆时偏移实现过程包括基于双程波方程逆时波场外推和应用成像条件两个步骤。对于单炮道集的成像,首先对震源波场利用双程波方程进行正向外推,并保存所有时间步的外推波场。然后对接收波场利用双程波方程进行逆时外推,在时间上每逆时外推一步之后应用成像条件,得到该时刻的成像结果。所有时刻的成像结果叠加,得到该炮集的成像结果。所有炮集的逆时偏移结果叠加即可得到最终叠前深度偏移成像结果。
在三维或二维正演模拟及逆时深度偏移中,当利用截断误差为 O=(△x2,△y2,△z2,△t2)的差分格式时,为保证频散较小及递推过程稳定,差分网格要求取得非常小,这样计算需要的计算机内存及运算时间会大大增加。Dablain利用高阶差分方模拟和偏移过程。利用高阶差分方程时,网格值可以取得大些,而计算精度并不降低。这里认为截断误差高于四阶的差分方程为高阶差分方程。三维声波方程的高阶差分方程可以用统一的方式推导出来,即
其中:u=(x,y,z,t)为地表记录的压力波场; v=(x,y,z)为纵横向可变的介质速度。为导出方程(1)的离散差分格式,需把观测对应的地下介质分布区域或要对其进行地震波模拟的模型区域离散化,即把它们剖分成一个个的小方块。
令
为导出高阶差分方程,把波场进行Taylor展开,得到截断误差为O=(△xM△yM,△zM,△t2),对应的三维逆时深度偏移高阶差分程为
2逆时偏移效果分析
在辽河凹陷西部凹陷北部高升潜山,逆时偏移处理取得较好的成像效果,潜山成像精度提高,能量聚焦,潜山内幕信噪比、分辨率明显提高(图1)。与克希霍夫深度偏移相比,逆时偏移在陡坡部位、低信噪比部位成像都有了提高,反映出来的地质信息更丰富,对潜山顶部小断层的刻画更清楚,潜山内幕的形态也更清楚。
3 小结
目前,逆时偏移技术在大面积的生产当中应用较少,多数是以试验性为目的。本文充分发挥逆时偏移技术在复杂构造条件下的成像优势,成功的实现了潜山的准确成像。由于目前针对逆时偏移的速度建模技术还不成熟,我们建模的手段仍然是基于克希霍夫偏移的,所以逆时偏移技术仍有很大的发挥空间。
参考文献
[1] 底青云,王妙月