有限元_有限差分法二维波动逆时偏移初探.pdfVIP

( ) 2003 年 　第 27 卷 　　　　　　　　　石油大学学报 自然科学版 　　　　　　　　　　Vol . 27 　No . 6 　第 6 期 　　　　　　　　Journal of t he University of Petroleum , China 　　　　　　　　Dec . 2003 　　文章编号 (2003) 有限元 有限差分法二维波动逆时偏移初探 1 2 3 董 　渊 , 杨慧珠 , 杜启振 ( 1. 清华大学计算机科学与技术系 ,北京 100084 ; 2 . 清华大学工程力学系 ,北京 100084 ; 3 . 石油大学资源与信息学院 ,山东东营 25706 1) 　　摘要 :提出了采用有限元 有限差分实现二维波动方程的逆时偏移算法 。该方法在空间上 ,联合采用有限元法和 ( ) ( ) 有限差分法 ;对于地表 水平 方向 ,使用有限元法进行离散 ,将原方程转化为一个一维 深度和时间 问题的方程组 ; ( ) 在深度和时间方向上 ,采用有限差分法来求解 。介绍了算法的基本原理 ,给出了计算实例并与使用 FK 频率 波数 域相移法 、频率 空间域有限差分法的结果进行了比较 。与采用有限元的偏移方法相比 ,本方法可以节省大量内存 ; 与采用有限差分的偏移方法相比 ,可以在一定程度上提高计算精度 。本算法有可能在地震勘探数据处理中发挥一 定的作用 。 关键词 :波动方程 ;逆时偏移 ;有限元法 ;有限差分法 ;计算精度 中图分类号 :P 63 1. 443 　　　文献标识码 :A 有限差分法用于二维波动逆时偏移 。 引　言 1 　原理概述 偏移技术是地震勘探数据处理的关键技术之 一[ 1 ] 。当今的偏移技术主要有 以波动方程为基础 以二维声波方程反演为例, 波场 P 是水平方向 的克希霍夫积分法 、频率 - 波数法 、有限差分法以及 x 、深度方向 z 和时间t 的函数, a ( x , z) 为地下介质 有限元法等[2 ] 。其中的有限差分方法使用最为广 速度 。有限元 有限差分法对 x 方向采用有限元离 泛 ,该算法大多针对波动方程的近似方程进行差分 散, 得到偏微分方程组之后对 z 、t 使用有限差分方 离散 ,能够适应任意复杂的地下结构和速度变化 ,计 法求解 。 2 2 2 算效率高 ,但存在计算精度低 、频散严重的现象 ,计 u u 1 u ( ) 　 2 + 2 = 2 2 , 1 x z