PmL吸收边界条件影响因素分析 - 物探与化探.PDF

　第３６卷第４期 物　探　与　化　探 Ｖｏｌ．３６，Ｎｏ．４　 　２０１２年８月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ Ａｕｇ．，２０１２　 ＰＭＬ吸收边界条件影响因素分析 丁　科 （中南林业科技大学土木工程与力学学院，湖南长沙　４１０００４） 摘 要：影响完全匹配层方法吸收效果的主要因素有吸收层厚度和衰减系数，笔者通过地震波场数值模拟，讨论了 吸收层厚度和衰减系数对吸收效果的影响。研究表明，在衰减系数一定时，吸收层厚度越大，边界反射越弱。吸收 层的厚度一般在１２～２０道较为适当，在吸收层的厚度比较小时，随着衰减系数的增大，边界反射逐渐减弱甚至完 全没有边界反射，但是若再进一步增大衰减系数则又会逐步出现边界反射。因此，在实际应用中应该注意衰减系 数的选择，衰减系数一般在５００～２０００之间，衰减系数较大反而影响其效果。 关键词：ＰＭＬ吸收边界；数值模拟；影响因素；波动方程 中图分类号：Ｐ６３１．４　　　文献标识码：Ａ　　　文章编号：１０００－８９１８（２０１２）０４－０６２３－０５ 　　波场数值模拟技术对于研究波在声介质、弹性 界，在模拟区域外增加一定厚度的吸收层，使入射波 介质、粘弹性介质等各种介质中的传播规律有着非 的能量逐渐变为零，这类边界条件目前应用比较广 常重要的作用，也是用来检验各种数据处理方法有 泛。常规吸收边界法是先导出吸收边界条件方程， 效性的常用工具。通过正演模拟，可以正确认识波 然后让方程与计算区域内的波动方程联立求解，使 在复杂介质中传播的运动学和动力学特征，准确分 从人工边界向计算区域反射回来的波能够全部或部 析地下地质构造所产生的反射波场特征。通过对已 分被吸收，而ＰＭＬ法则是在所研究区域的边界上引 知模型的正演模拟研究得到的规律，可以加强人们 入吸收层，波由研究区域边界传到吸收层时不产生 对未知模型的认识能力，有助于解决实际问题。 任何反射，在吸收层内按传播距离的指数规律衰减， 目前波动方程数值解方法主要有有限差分法、 不产生反射，从而达到吸收边界的效果。目前，ＰＭＬ 频率—波数域法、Ｋｉｒｃｈｈｏｆｆ积分法、有限元方法、普 方法已在电磁波、声波、弹性波等多个领域得到广泛 元法、边界元法、广义屏法、格子法等。由于波在地 的应用，是一种效果很好的吸收边界的方法。影响 下介质中的传播是在无限区域内传播，而进行数值 ＰＭＬ方法吸收效果的主要因素有吸收层厚度和衰 模拟时的计算区域必定是有限的，也就是在计算区 减系数，笔者对它们的影响效果进行了分析。 域内介质的各种参数不为零，而在计算区域之外各 １　基本原理 参数均为零，这就人为地设置了边界，波在通过这些 边界时，就会产生反射，干扰了有效波的信息，因此 ＰＭＬ法的主要思想是在计算区域的周围加上 不管哪一种数值模拟方法都会遇到人为的边界反射 完全匹配层，其作用是吸收边界反射波。对于二维 问题。在进行模拟时，必须消除这种人工界面所产 常密度情况下的纵波方程，可以表示为 ２ ２ ２ 生的反射波等的干扰，这样就要有适当的边界处理 １ ｕ（ｘ，ｚ，ｔ） ｕ（ｘ，ｚ，ｔ） ｕ（ｘ，ｚ，ｔ）    ２ ２ ＝ ２ ＋