流固耦合地震波动问题的显式谱元模拟方法

摘要

流固耦合地震波动问题主要研究由流体和固体构成的复杂系统中地震波传播特性及其规律.传统模拟方法中一般以声波方程、弹性波方程的数值解分别描述理想流体和弹性固体中的波动,并实时地处理两种不同性质介质之间的相互耦合作用,数值格式复杂且限制数值模拟精度与计算效率.本文采用谱元法结合多次透射公式人工边界条件实现了一种流固耦合地震波动问题的高阶显式数值计算方法.该方法利用了流固耦合问题统一计算框架,可将饱和多孔介质的Biot波动方程分别退化为理想流体的声波方程和弹性固体的弹性波方程.通过P波垂直入射的水平成层理想流体-饱和多孔介质-弹性固体场地模型、P波斜入射的不规则层状界面以及任意形状界面的理想流体-饱和多孔介质-弹性固体场地模型等三个算例,与传递函数法解析解以及集中质量有限元法计算结果进行对比分析,证明了本文方法的正确性与有效性.数值模拟结果表明,本文方法相较传统有限元法可以少得多的节点数量获得更高的数值精度,并且在较宽的频率范围内都能可靠地模拟出流固耦合系统的动力响应,充分体现出本文方法兼顾高精度、计算效率和复杂场地建模灵活的特点.