基于cerjan衰减函数的pml吸收边界条件 - 物探与化探.pdf

　第３６卷第６期 物　探　与　化　探 Ｖｏｌ．３６，Ｎｏ．６　 　２０１２年１２月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ Ｄｅｃ．，２０１２　 基于Ｃｅｒｊａｎ衰减函数的ＰＭＬ吸收边界条件 １ １，２ １ 李飞 ，程久龙 ，杨思通 （１．山东科技大学矿山灾害预防控制教育部重点实验室，山东 青岛　２６６５９０；２．中国矿业大学煤 炭资源与安全开采国家重点实验室，北京　１０００８３） 摘 要：在Ｃｅｒｊａｎ吸收边界和ＰＭＬ吸收边界的基础上提出了基于Ｃｅｒｊａｎ衰减函数的ＰＭＬ吸收边界。通过两种吸收 边界条件下地震波场模拟计算比较表明，后者具有吸收更干净且能够吸收各种角度的边界反射等优点，并且保留 了前者操作简单的特点，计算时间少，适用于二维情形下的地震波数值模拟。 关键词：地震波；数值模拟；有限差分；Ｃｅｒｊａｎ衰减函数；ＰＭＬ 中图分类号：Ｐ６３１．４　　　文献标识码：Ａ　　　文章编号：１０００－８９１８（２０１２）０６－０９９６－０５ 　　实际中地震波是在半无限介质中传播的，而地 和基于Ｃｅｒｊａｎ衰减函数的ＰＭＬ吸收边界进行比较 震波的数值模拟只能在有限的区域内进行，这种有 研究，分析二者的优缺点。 限区域的边界除自由地表外就是人为的边界，即人 １　弹性波动方程交错网格有限差分数值模拟 工边界。在进行地震波数值模拟时，若在人工边界 处不作任何技术处理，则地震波在人工边界处会产 ［８］ 二维波动方程应力与速度关系 为 生很强的反射，使模拟结果受到严重的干扰［１－３］。 ｕ／ｔ＝１／（ ／ｘ＋ ／ｙ），   ρσ  σ   ｘｘ ｘｙ 为了消除或减弱人工边界的影响，人们提出了  ｗ／ｔ＝１／（ ／ｘ＋ ／ｙ），   ρσ  σ  ［４－５］  ｘｙ ｙｙ  各种方法 。主要的方法有：扩展边界法，即将模 ／ｔ＝（ ＋２）ｕ／ｘ＋ ｗ／ｙ， （１）  σ  λ μ   λ  ｘｘ 拟区域选的足够大，减小边界反射对所需结果的影  ／ｔ＝ ｕ／ｘ＋（ ＋２）ｗ／ｙ， σ  λ  λ μ   