非均匀介质中地震波应力场的wnad方法及其数值模拟.pdfVIP

第 卷 第 期 地 球 物 理 学 报 ， ５２ ６ Ｖｏｌ．５２ Ｎｏ．６ 　 　 　 　 　 年 月 ， ２００９ ６ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ｊｕｎｅ ２００９ 　 　 　 王 磊，杨顶辉，邓小英 非均匀介质中地震波应力场的 方法及其数值模拟 地球物理学报， ， （）： ， ． ＷＮＡＤ ． ２００９５２６ １５２６ １５３５ 　 ～ ： ／ ＤＯＩ１０．３９６９ ．ｉｓｓｎ．０００１５７３３．２００９．０６．０１４ ｊ ， ， Ｗａｎ Ｌ Ｙａｎ ＤＨ Ｄｅｎ ＸＹ．ＡＷＮＡＤｍｅｔｈｏｄｆｏｒｓｅｉｓｍｉｃｓｔｒｅｓｓｆｉｅｌｄｍｏｄｅｌｉｎ ｉｎｈｅｔｅｒｏｅｎｅｏｕｓｍｅｄｉａ．犆犺犻狀犲狊犲犑． ｇ ｇ ｇ ｇ ｇ （ ）， ， （）： ， ： ／ 犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ２００９ ５２６ １５２６ １５３５ＤＯＩ１０．３９６９ ．ｉｓｓｎ．０００１５７３３．２００９．０６．０１４ 狆狔 ～ ｊ 非均匀介质中地震波应力场的犠犖犃犇方法 及其数值模拟  王 磊，杨顶辉 ，邓小英 　 清华大学数学科学系，北京 １０００８４ 　 摘 要 通过对近似解析离散化（ ）方法的分析，给出了一种求解声波和弹性波方程的带权重的近似解析离散 ＮＡＤ 　 　 化（ ）方法，并用 方法、 修正格式（ ）和二阶中心差分方法计算了二维波动方程初 ＷＮＡＤ ＷＮＡＤ ＬａｘＷｅｎｄｒｏｆｆ ＬＷＣ 值问题的应力场数值误差 结果表明 方法具有更高的数值精度 用 方法、 和四阶交错网格法 ． ＷＮＡＤ ． ＷＮＡＤ ＬＷＣ 对二维非均匀介质中弹性波传播的应力场进行了数值模拟 应力场快照和地表地震记录表明，即使是在粗网格条 ． 件下 方法的模拟结果仍无可见的数值频散和源噪声 另一方面，由于 方法同时计算了地震位移和 ＷＮＡＤ ． ＷＮＡＤ 梯度场，使得应力的计算更为便捷和精确，而且 ＷＮＡＤ方法中波位移梯度局部连接关系的使用使得应力在间断处 能够自动近似地满足应力连续性． 关键词 ＷＮＡＤ，波动方程，数值频散，应力场模拟 　 ：