【Matlab】地震波传播时域有限差分仿真
一、引言
地震动力学仿真是岩土工程、地震工程、防灾减灾工程的核心研究方向，主要研究地震波在地下岩土介质中的传播规律、场地动力响应与地表振动特征。天然地震发生过程中，弹性地震波在地层内部持续传播、反射与折射，引发地表土层与工程结构产生振动变形，是造成建筑损毁、边坡失稳、工程灾害的核心诱因。准确掌握地震波在不同地层结构中的传播特性，能够为场地地震安全性评价、建筑抗震设计、地下工程防震减灾、地震灾害预判提供核心数据支撑。
传统地震波研究以现场监测与理论解析计算为主，现场监测成本高、数据随机性强、可控性差，无法系统探究地层介质、场地结构、波速参数对地震传播的影响规律。经典解析法仅适用于均匀单层地层、简单平面波入射的理想工况，无法适配分层地层、不均匀介质、复杂地形的真实场地条件，难以精准还原地震波的散射、反射与波形畸变特征。随着数值仿真技术的发展，时域有限差分算法逐步成为地震波传播仿真的主流手段，凭借时域迭代、全域离散、自适应强的优势，可高效求解复杂地层的地震波传播问题。
时域有限差分法简称FDTD，是求解波动方程的经典数值算法，核心思路是将空间求解域与时间域进行双重离散，通过中心差分格式替换波动偏微分方程，通过逐时刻迭代计算实现波场的动态演化，完整还原地震波的传播、反射、折射与衰减过程。该算法原理清晰、迭代逻辑简单、对复杂介质适配性强，无需复杂的网格剖分与矩阵求解，计算稳定性高，极其适合弹性波、地震波的时域动态仿真。
MATLAB具备强大的矩阵迭代运算、时域离散计算与数据可视化能力，能够快速完成地层网格离散、波动方程差分离散、震源加载、边界条件处理、波场迭代更新与波形云图绘制，可高效实现地震波传播的全过程仿真。本文系统阐述地震波波动理论与时域有限差分
