复杂介质地震波传播模拟中边界元法与有限差分法的比较研究.pdfVIP

第28卷第2期 地球物理学进展 V01．28。No．2 2013年4月(页码：0664—0678)PR()GRESSINGEOPHYSICS Apr．，2013 0678，doi：10．6038／p920130215． for HE of elementmethodand“nite_differencemethod Yan—feng。SUNWei—jia，FUI。i—yun．Comparisonboundary simulating seismicwave in in propagationcomplexmedia．ProgressGeophys．(in 3021 p92015． 复杂介质地震波传播模拟中边界元法 与有限差分法的比较研究 何彦锋1’2， 孙伟家1， 符力耘1 (1．中国科学院地质与地球物理研究所，北京100029；2．中国科学院大学，北京100049) 摘要地震波场数值模拟是理论地球物理学和勘探地球物理学的重要研究手段．在众多数值模拟方法中边界元 法和有限差分法是两种典型的地震波传播模拟计算方法．边界元法是一种半解析一半数值的边界型方法，它显式地 利用边界连续条件，沿着地层边界进行离散，具有降维、高精度和自动满足远场辐射条件的优点；有限差分法是一 种典型的基于微分的区域型方法，它隐式地使用边界连续条件，以空间网格形式进行离散和数值逼近，具有高效、 实用和容易数值实现的优点．本文以一个半圆形均匀Valley模型和两个非均匀断裂／断层模型为例，从计算精度、 计算效率、频散特性以及适用性等方面对这两种方法进行了比较研究．数值计算结果表明：边界元法可以精确地几 何描述有内部断点、断面的复杂构造，能够精确地模拟内部不规则界面之间波的反射／传播；有限差分法不能以足 够的精度描述几何断点和内部不规则边界．边界元法在高频时计算量大于有限差分法，有限差分法则需要更小的 网格间距以压制数值频散．因此，在处理内部非均质和高频计算时，有限差分法更有效；在处理内部不规则边界、断 点、大尺度等问题时，边界元法比有限差分法更有优势． SH波 关键词边界元法，有限差分法，数值模拟，地震波传播，2一D doi：10．6038／p920130215中图分类号P315 文献标识码 A of elementmethodandfinite—differencemethod Comparisonboundary for seismicwave in media simulating propagationcomplex HE SUN ial。 FU Yan-fen91“．Wei-j Li—yunl (1．Institute and 100029，China； ofGeologyGeophysics，ChineseAcademyofSciences，Beijing oChinese o 100049，China) 2．UniversitYf AcademYfSciences，Beiiing an inbo