波动方程有限差分法频散衰减方法研究.pdfVIP

第 34卷第 12期 刘洋等 ：波动方程有限差分法频散衰减方法研究 据柯 朗一弗里德里希斯一列维条件 ，在三维正演模拟中，每个时问步长必须满足下式 ： ‘ At≤ 一——：： 一 (6) ． K／ ／ __ IliaX △ 。△v。 Az。 式中，△ 为时间采样 间隔，S；77 为最大模型速度 ，m／s；Am、Ay分别为x、y方向的空间采样间隔，m。 由于有式 (6)的影响，在二维波动方程式正演下，当网格变小 ，有限差分总计算量将会成 3次方 增长；在二维波动方程式正演下，有限差分总计算量将会成 4次方增长 。图1是使用标准波动方程式在 各向同性介质 中250ms处 的波场快照，震源使用的均是雷克子波 ，其 中图 1(a)、(b)、(c)使用 的计 算网格是 201×201。图 1(d)使用的是 401×401，但是 图1(d)在 X方 向和Z方向的采样间隔是图 1 (a)～ (c)的 1／2，网格大小是图 1(a)～ (c)所使用的 1／4。因此 ，尽管 图 1(d)比图 1(a)的效 果要好得多，但是图 1(d)需要的计算量大约是 图 1(a)的8倍。如果能够在 2阶精度的有限差分法 的基础上压制频散 ，那么就能够使用 比较低 的计算量而得到更好效果 的结果 。 2 声波波动方程式优化 笔者将式 (1)的解用下式表示 ： “(￡，x)一fp(￡，X，，K)dK (7) 式中，X一 (，，z)表示空间位置 ；为式(1)的一个解 ；K一 (忌 志， )表示空间波数 ；是、走、尼。分别为 X、y、Z方向的波数 。 上面已经证明了在用有限差分法模拟波场时，大波数的相速度与小波速的相速度不同，而且波数越大 的，相 比滞后群速度得越多。因此，给出下式 ： “(t，X)一 lg(K)P(t，X，∞，K)dK (8) 式 中，g(K)为阻尼 因子 。相 比式(7)，式(8)多出一个阻尼因子 g(K)，其是平面波解的阻尼 函数 ，g(K)将 在大波数时压制频散噪声 。下面给出具有式 (1)形式的解 ： 一(z，，)等at。一一筹a32。—+十。一等一Y。—+。碧～a。—+。—v急———(x———，——，——一)l(＼——刁——z——●——a—一t—+r一刁at—+r’——a—z●——a—一t)I／ ㈤I、J 7d Y 式 中，r(x，，)为随空间变化的非负实数 。表面上看 ，式(9)比式 (1)要复杂得多；然而 ，式(9)等号右边多 出的部分只是对式(1)简单 的修改 。式(9)能够有效地利用低阶(2阶或者 4阶)有限差分法 。 图2验证了不 同空间精度有 限差分法计算的单炮正演模型 ，图2 (a)与图 2 (b)均使用的是空间 2阶有限差分法 ，区别在于图2 (a)使用 的是式 (1)，而图2 (b)使用 的是式 (9)，图2 (b)中的频 散 比图2 (a)中的频散现象明显少得多。图2 (c)和图2 (d)均使用式 (1)计算，但是前者使用的是 空间8阶精度有限差分 ，而后者使用的是空间2阶精度有 限差分，且后者使用的网格在 x、Z方 向的采 样间隔均是前者的 1／2。图 2 (b)的效果几乎可以达到图 2 (c)、图 2 (d)的效果 ，在单炮正演记录上 基本观测不到频散的影响；而图2 (b)的计算量 比图2 (d)的计算量要小得多。 为 了保证式 (9)的解有意义，r(x，Y，z)必须满足下式 ： r( ， ，)< — ：二= ：= == (10) ／＼／／1十 1十．一1 · A—z2 还有一个简单的处理r(x，，z)值的方法是把r(x，，)设置为一个常数R，R的值与差分方法的阶数 和差分法计算所用的网格大小有关 ，在该次正演模拟中