一种改进的利用程函方程计算旅行时的方法.pdfVIP

　2008 年 10 月 第 43 卷 　第 5 期 　 ·正演技术 · 一种改进的利用程函方程计算旅行时的方法 李文杰 　魏修成 　宁俊瑞 　张建伟 ( 中国石化石油勘探开发研究院 ,北京 100083) 李文杰 ,魏修成 , 宁俊瑞 , 张建伟. 一种改进的利用程函方程计算旅行时的方法. 石油地球物理勘探 ,2008 ,43( 5) : 589～594 摘要 　地震波旅行时计算是实现叠前深度偏移成像和地震层析成像的核心内容 ,其精度决定着最终成像效果 。 通过利用有限差分法求解程函方程计算地震波旅行时的方法有许多优点 ,一直受到人们的青睐 。本文在前人 提出的已知三点初至时间用程函方程求第四点初至时间的基础上 ,提出一种新的已知两点初至时间用程函方 程求相邻第三点初至时间的计算方法 ,并将扩展波阵面法与局部算法综合应用 ,使算法更加稳健 。 关键词 　程函方程 　有限差分 　初至时间　射线追踪 　纵波 震波传播的因果性关系[2 ] ,使得所计算的时间并非 1 　引言 真正的初至时间。为了克服该法所具有的局限性 , [3 ] Cao 等人 提出了利用交错网格差分波前追踪方法 叠前深度偏移成像和地震层析成像技术的应用 得到程函方程有限差分解使计算出的各个网格点初 [2 ] 越来越广泛 ,地震波旅行时计算技术作为其核心内 至时间更精确 ; Qin 等人 在 Vidale 提出的方法基 容一直是地震资料处理领域里的研究课题 。以前国 础上 ,采用扩张波阵面法追踪波前 ,并通过在上一个 内外许多学者采用解析法和射线追踪法计算地下各 波前面上寻找一个最小旅行时点 ,然后向外扩张 ,从 点的初至时间及相关研究 ,取得了大量的研究成果 , 而使得在计算过程中能够避免因果性问题出现的同 然而利用这些方法求取初至时间都不同程度地存在 时 ,还可以提高计算精度 ; Eiichi 等提出了一种插值 着局限性 。主要表现为 :解析法只能适用于少数特 法的思想来改进算法的精度[4 ] ; 王华忠等为改进 Vidale 提出的方法也做了许多有意义的工作[ 5 ,6 ] 。 殊的速度分布 ,对于一些较复杂的速度分布则难以 而何兵寿 [7 ] 对 Qin 等人提出的通过扩张波阵面法 应用 ;而利用射线追踪法求取初至时间有时会出现 追踪波前的方法做了进一步的改进 。本文在前人提 盲区 ,特别是针对一些较复杂的速度分布区域尤其 出的已知三点初至时间用程函方程求第四点初至时 如此 ;利用解析法和射线追踪法计算出的初至时间 间的基础上 ,提出了一种新的已知两点初至时间用 可能只是局部极小时间 ,并非真正的初至时间;利用 程函方程求相邻第三点初至时间的计算方法 ,并将 解析法和射线追踪法计算初至时间一般说来都比较 扩展波阵面法与