利用谱元法计算弹性波场的若干理论问题.pdfVIP

中国科学 G 辑: 物理学 力学 天文学 2007 年 第 37 卷 第 1 期: 41~59 收稿日期: 2006-01-10; 接受日期: 2006-08-28 国家自然科学基金 (批准号:、UNESCO 国际理论物理中心基金、澳大利亚教育和科技部国际科学合作计 划以及中国科学院百人计 划联合基金资助项目 * E-mail: wangxm@ 《中国科学》杂志社 SCIENCE IN CHINA PRESS 利用谱元法计算弹性波场的若干理论问题 王秀明①②* G. SERIANI③ 林伟军① (① 中国科学院声学研究所, 北京 100080; ② CSIRO Petroleum, ARRC, POBOX 1130, Technology Park, Bentley, WA 6102, Australia; ③ Instituto Nazonnale di Oceangrafia e di Geofisica Sperimentale Borga Grotta Gigante 42/c-Sgonico, I-34016, Trieste, Italy) 摘要 近来在解决弹性动力学问题中新发展起来的谱元法是一种基于组成整 个研究区域的单元 Galerkin 算法, 并利用弹性波动力学方程的弱形式进行数值求 解的方法. 针对求解弹性波动方程的谱元法中的若干数学问题以及相关的算法原 理进行了系统研究, 介绍了Legendre 和Chebyshev多项式构造的基函数及其用在 参考元上数值积分的 Gauss-Lobbatto 配置点 , 推导了在利用 Legendre 和 Chebyshev 多项式展开时的单元积分的具体表达式. 在声阻抗差别很大的非均匀 介质中, 通过引入空间域中在预先条件下的共轭梯度的元到元算法和时间域中时 间积分的交错网格的预期/多次校正算法, 谱元法不需要形成有限元中的全局矩 阵和有效载荷矢量, 能够在很大程度上提高计算精度和计算效率. 另外, 在某些 情况下, 如果再利用单元积分的解析式, 也不需要形成单元矩阵, 由于采用解析 式, 计算中可以节约大量内存而不会损失太多计算效率. 元到元算法中使用了最 优的张量乘积技术, 使得该方法比有限元法在内存需求量和计算时间等方面更为 有效. 关键词 弱形式 弹性波模拟 谱元 元到元算法 目前, 尽管计算机的运算能力已有相当大的提高并在不断的迅速发展, 但在计算声学和 数值模拟波的传播方面 , 精确模拟弹性动力学问题仍面临挑战 , 例如 , 在石油勘探领域中, 常 常需要模拟地震波在三维带有不规则自由界面的大尺度复杂地质构造中的传播, 同时, 在理 解波在含有多相流体的孔隙地层中的传播特性时, 我们还需要考虑流变理论. 在地震学和超 声检测方面, 实验观测点往往放置在自由表面上或是自由表面附近, 而自由表面可能会对接 收到的信号有很强的影响. 因此, 模拟和理解这些不规则界面对弹性波传播的影响在很多领 域已经引起了人们的高度关注, 如地震勘探[1~3]、天然地震[4]以及无损检测等[5]. 42 中国科学 G 辑 物理学 力学 天文学 第 37 卷 目前, 已有几种较为成熟的数值方法用于模拟弹性波在具有自由起伏界面的非均匀介质 中的传播, 包括有限元方法、边界元法、有限差分算法、伪谱法以及谱元法. 有限元法可以处理带有自由起伏界面(气-固交界)的复杂介质模型. 该方法基于波动方程 弱形式, 能够自动满足自由边界条件 , 同时在解决非均匀复杂弹性介质模型时也很有优势[6,7]. 然而与显式的有限差分相比, 有限元法要花费更多的计算时间[8], 这可能是有限元方法没有在 石油地球物理勘探领域被广泛应用的主要原因. 另外, 低阶有限元法会出现令人不满意的频 散现象[9], 而传统的高阶有限元法会产生虚假波[4]. 对于有限差分方法而言, 它处理有规则的边界模型及介质速度或密度差别不是很大时的 模型有优势[10~12]. 处理不规则自由表面的有限差分算法或者改进数值精度的方法已被很多学 者做过研究并报道[1~3,13], 尽管提出的算法容易实现, 但这里需要说明的是, 在交错网格的有 限差分中无论采用多么高阶的精度, 即使是平面自由表面, 其处理自由界面的精度也不超过 二阶. 在流体力学中的伪谱法已用于解决弹性动力学问题[14~17], 该方法基于快速 Fourier 变换或 Chebyshev变换 , 精度较高, 同时在 Nyquist波长取很少的空间采样点时空间频散小 , 使其成为 一种目前被