起伏地表二阶弹性波方程差分策略稳定性分析.pdfVIP

第 52卷第5期 石 油 物 探 Vo1．52，No．5 2013年 9月 GEOPHYSICALPR0SPECTING FORPETR0LEUM Sep．，2013 文章编号：1000—1441(2013)05—0457—07 起伏地表二阶弹性波方程差分策略稳定性分析 唐 文 ，王尚旭 ，袁三一。 (1．西安石油大学地球科学与工程学院，陕西西安710065；2．中国石油大学(北京)CNPC物探重 点实验室，北京 102249) 摘要：采用坐标变换的思想将起伏地表情况的曲线网格变换成水平地表情况的矩形 网格，推导出了变换后的二 阶弹性波动方程，研究了基于坐标变换的起伏地表二阶弹性波场有限差分数值模拟方法。重点针对起伏地表变 换成水平地表后的3种不同特点的自由边界条件，包括隐式、单边伪节点和混合 自由地表边界条件的处理，并给 出了具体的差分格式。数值模拟结果表明，起伏地表二阶弹性波方程模拟方法是有效可行 的，但不同的地表处 理方法对模拟结果有不同的影响。在这 3种 自由地表处理方法中，隐式的自由地表边界条件有最好的稳定性。 关键词 ：起伏地表；二阶弹性波方程；有限差分法；边界条件；稳定性 DOI：10．3969／j．issn．1000—1441．2013．05．003 中图分类号：P631．4 文献标识码：A 对起伏地表弹性介质模型的模拟具有现实意 义Ⅲ1≈]，尤其是模拟近地表散射波效应。起伏地表地 1 坐标变换后的二阶弹性波动方程 震波场的数值模拟技术很多，但各有优缺点。有限 元法[4-6]计算成本比较高；边界元法_7不适宜模拟复 考虑弹性波在一个均匀的、各 向同性的非规则 杂的地质模型；伪谱法8【经常会伴随有折叠效应、吉 二维区域 (图la)传播。根据动量守恒、应力一应变 布斯效应；有限差分技术_9]是一种广泛应用的简单 关系和虎克定律，二阶弹性波动方程 可表示为 有效的模拟方法。如果使用常规的矩形网格去逼近 fl9211一一诉；象+十考+十L(；一一礁 不规则地表，经常会导致出现绕射干扰波。为了应 用有限差分技术模拟起伏地表条件下的波场，许多 I【一一学+十 +十L(诉一一诣)—血Ox—3z 学者进行了深入研究。Jih等[1，Hestholm等[11-14]， 其 中，(z，)为笛卡儿坐标系下的坐标；Uzu(x，， 王祥春等I15]以及孙章庆[1。]采用非规则网格进行剖 ￡)为沿 z方向的位移；砌：== (z，z，￡)为沿 方向的 分，对 自由表面采用旋转坐标系和线性插值等技术 位移；Vp一 (z，)为纵波速度；7Js—VS(z，)为横 进行处理。由于采用非规则网格，因此增加了复杂 波速度 。 度并降低了模拟精度。Komatitsch等_1]提出了一 为了在矩形区域矩形网格下实现弹性波模拟， 种直接在曲坐标系下进行有限差分的张量方法，该 我们采用了坐标变换将物理域(图la)变换到了计 方法计算效率高，但需要较多内存，并且对界面和地 算域 (图lb)。采用的映射函数为[1卜] 表的网格剖分相当麻烦，也不易推广到三维。Moczo 等[1。和黄 自萍等[19]对有限元和有限差分方法进行 』z(，)一 (2) 【(，r／)一 0()+7[ 一 o