各向异性介质中弹性波的数值模拟Ξ.PDFVIP

Vol. 21 　No. 3 　　 　　　　　　 　 　　　　　　　　　 　 第21卷 第3期 固 体 力 学 学 报 2000年 　9 月 ACTA MECHANICA SOLIDA SINICA September 2000 各向异性介质中弹性波的数值模拟 张剑锋 　　刘铁林 (大连理工大学工业装备结构分析国家重点实验室 ,大连 ,116024) 摘 　要 　提出了一种非均匀各向异性介质中弹性波传播的数值模拟算法. 该方法可以灵活 地运用于具有任意地表形状、内部孔洞、固液边界和不规则内部交界面的介质情况 ,另外 ,该方 法自然满足复杂几何边界的自由表面条件. 这种基于三角形和四边形离散网格的算法使用的是 围绕每个节点的积分平衡方程 ,而不是其它有限差分法中使用的各个节点满足的弹性动力学的 微分方程. 该文工作是非均匀各向同性介质中弹性波传播格子法研究的继续. 除了研究各向异 性介质中波的传播以外 ,还给出了一种能够省时的四边形网格的格子法. 关键词 　各向异性 ,数值模拟 ,有限差分 ,格子法 1 　引言 沉积岩通常是具有各向异性的结构 ,例如 ,细小尺度的层状介质、具有一定方向的微小 裂缝 结构以及有优先指向的非球形粒子或各向异性矿物质[1 ] . 各向异性介质中地震波传播 [2 ] 的数值模拟的研究可以使我们更好地解释地震图. 目前的模拟方法有渐进射线理论 ,反射 法[3 ] ,有限差分方法[4 ] . 垂直方向分层的任意各向异性结构运用反射法时 ,构造传播算子的 特征值和特征向量要用数值形式给出 ,此时 ,反射法耗时甚多. 有限差分法功能强且使用灵 活 ,它可以处理复杂的横向和纵向不均匀以及任意各向异性介质中弹性波的传播问题. 如对 于横观各向同性介质的矢量化有限差分算法[5 ] ,模拟二维和三维各向异性介质中弹性波传 [6 ] [7 ] [8 ] 播问题的伪谱法 和快速展开法 . 文献[4 ]运用交错网格差分法 计算了横观各向同性介 质中的地震波传播问题. 然而 ,这些有限差分算法都是基于长方形网格 ,因此 ,很难用它们来 解决具有地表起伏、内部孔洞和曲线型内部交界面时的波传播问题. 另外 ,用阶梯状折线来 模拟原曲线型交界面时 ,会产生人为的绕射现象 ,特别是用高阶有限差分法时更明显[9 ] . 文献[10 ]提出了非均匀介质中 PSV 波传播数值模拟的格子法. 本文将其推广到非均匀 各向异性结构 ,特别是由平行裂缝或细密层形成的各向异性结构. 本文采用一种简化四边形 格子替代原三角形格子 ,它比原三角形格子法的计算效率高. 格子法可以非常容易地处理具 有复杂几何表面和内部交界面固体中弹性波的传播问题. 关于垂直自由表面的轴有对称性 的横观各向同性弹性半无限空间介质 ,用本文方法计算并与 Lamb 解进行了比较 ,结果表 明 ,本文方法具有很高的精度. 研究了均匀横观各向同性介质的辐射状模型 ,此时 ,无需大量 计算即可得到剪切波分离现象. 还研究了表面有半圆形凹陷的弹性各向异性半无限介质中 ( ) 表面波传播以及交界面是斜面的固 各向异性 液混合体中弹性波的传播问题. ( ) 国家自然科学基金资助项目 . 收到第 1 稿 ,199