多松弛时间格子Boltzmann方法的黏滞吸收边界

当介质的几何结构复杂或者内部存在强物性界面时，传统的地震波正演模拟方法的计算结果往往难以满足实际精细波场计算的要求。多松弛时间格子Boltzmann方法（MRT-LBM）是一种新兴数值模拟方法，具有稳定性好、计算精度高、边界处理灵活等优点。针对MRT-LBM数值模拟时面临的人工截断边界，提出了一种基于多松弛参数的黏滞吸收边界方案。由于边界反射的压制效果对衰减参数很敏感，因此通过大量数值模拟实验选择了最优参数组合，获得了适用性强的吸收边界条件。该吸收边界条件的算法简单，且可扩展性强。最后应用均匀模型和简单非均质模型验证了其吸收效果，并用复杂的BP模型验证了其适用性。