一种地震波数值模拟的方法

摘要

本发明涉及地震波波场数值模拟技术领域，本发明提出的针对复杂地表、复杂构造中的地震波数值模拟的方法，针对原算法中容易出现震源噪音等假波问题的原因，在原算法中强制采用更严格的数值流方式进行相邻单元间的波场交换，如局部Lax‑Friedrich数值流，同时通过引入惩罚因子来降低强制采用局部Lax‑Friedrich数值流带来的时间步长降低等影响，既保持基于蛙跳时间DGFEM地震数值模拟算法的大步长、低计算代价的优点，又可以实现高精度无假波数值模拟。