地震声波数值模拟中的透射边界条件.pdfVIP

第9卷 第 l3期 2009年7月 科 学 技 术 与 工 程 Vo1．9 No．13 July2009 1671．-1819(2oo9)13--3607．-05 ScienceTechnologyandEngineering @ 2009 Sci．Tech．Engng． 地球科学 地震声波数值模拟中的透射边界条件 邢 丽 (上海第二工业大学数学系，上海 201209) 摘 要 地震波传播的数值模拟中，在有限区域建立吸收边界条件是一个很重要的问题。在不同类型的吸收边界条件 中，比 较经典地是Clayton和Engquist提 出得Clayton-Engquist吸收边界条件，透射边界条件 (MTF)是通过直接模拟各种单向波动的 共同运动学特征来建立的一种边界条件，算法简单容易。重点讨论 了MTF边界条件并通过数值计算与被传统应用的CE吸 收边界条件作比较。 关键词 地震声波数值模拟 有限元 有限差分 吸收边界条件 CE MTF 中图法分类号 P315．31； 文献标志码 A l 地震声波数学模型简要介绍 一 号旁-o，以及三阶近似：03U一2 + 考虑二维空间均匀各项同性介质中的声波方程 等一手c窘_0，这就是著名的Clayton-En- 一Ot2I(++0y)2U，z：J) (¨1)gquist旁轴近似 (paraxialapproximation)吸收边界条 件，简称CE。除此，本文重点介绍透射边界条件。 计算区域 ：{(，Y)：0≤ ≤口，0≤y≤ ()}，其中 Y=B(x)为地表 曲线 ，当B(x) c(C为常数)时为平 2 透射边界条件 面地表，=0，=a和Y=0为人工边界。对声波方 首先给出穿越人工边界上一点的单向波动的一 程 (1)需加以下几种定解条件： 般表达式 。设 轴在所考虑的边界点上垂直于人 (1)初始条件：当t≤O时，介质处于静止状态， 工边界并指向外部无限域，如图1。则单向波沿 轴 因此声波方程的初始条件为： 的视传播速度(apparentpropagation)的一般表达式为 (，Y，t)I≤0=0。 (2)自然边界(地表)：)，=B()上的边界条件： u(，)=∑ (c一) (2) J uO ： 0 On 人工边界 —— _。’ (3)人工边界 ：x=0，=a和Y=0上的人工边 计算区域 界条件 ，Clayton和 Engquist…提出了一系列具有不 I I 2 cat 1 cAt O 同精度的吸收边界条件 ： 。 … a 一 阶近似：_OU+ COt=o，二 阶近似：O + Ox xOt 图1 多次透射公式的几何关系 2009年3月23日收到 作者简介：邢 丽(1978一)，女，研究方向：偏微分方程数值解。E一