中图诲上袖钮-中国海上油气.pdfVIP

第21卷第4期 中国海上油气 V01．21No．4 2009年8月 CHlNAOFFSHORE()ILAND(；AS Aug．2009 标量声波波动方程高阶交错网格有限差分法 陈可洋 (大庆油田有限责任公司勘探开发研究院) 摘 要 为进一步提高标量声波波动方程正演数值模拟精度，减小数值频散和计算量，改善边界吸 收效果，提出将高阶交错网格有限差分法直接用于求解标量声波波动方程的方法，并推导出了时间 和空间导数的高阶交错网格有限差分格式，给出了数值模拟计算所需的稳定性条件和吸收边界条 件。模型试算结果表明，与传统的中心有限差分法相比，本文方法不仅可以快速提高时空导数的差 分精度，而且具有频散小、计算效率高等特点，可推广应用于三维标量声波高精度正演数值模拟和 叠后逆时深度偏移。 关键词 标量声波波动方程 高阶交错网格有限差分格式正演数值模拟 叠后逆时深度偏移 2 目前国内外对于标量声波波动方程的求解已形 的至I团导效，则U(z)在z 成多种精度较高的正演数值模拟方法[1-13]，而如何 2N+1阶Taylor展开式为 进一步提高数值模拟精度、减小数值频散和计算量、 让z-4-(m一1／2)ax3一u(xo)+ 改善边界的吸收效果等一直是有限差分算法研究的 核心问题。Alford和Kelly[14’153提出二阶中心有限 1／2)] +。(△zz卅z) 2∑N+I(优一-+E)073(，}“娑耸， ∑(优一t掣“}，+。(舻№) 差分法计算标量声波波动方程，但数值结果表明低阶 (2) Gazd— 差分算法具有强频散和低精度的缺点；Jenti 161、DanD．Kosloff[173等采用伪谱法计算标量声 agf 止百au(x)b一蚤N口。{“[z。+ 波波动方程，得到了高精度的数值模拟记录，但计算 时问长、存储量非常大，且较难处理吸收边界条件和 (3) 自由表面边界条件，还可能存在Gibbs效应。为了 有效解决上述方法所存在的问题，本文提出将高阶 式(2)、(3)中：“≯代表甜对于变量z的i阶导数，N 交错网格有限差分法直接应用于求解标量声波波动 为有限差分近似阶数。 方程，以期提高局部地震波场数值模拟和逆时深度 将式(2)代入式(3)，有 成像精度，并提高计算效率。 1 方法原理 “。～ amU(x2i+1)(引+ 1．1标量声波波动方程 ∑N N∑--I埃群Ax2件I 惫台(2i+1)! 声波波动方程有如下形式： “。一可2(Ux．r+U。) (1)