基于小波域Pre-GSOR迭代的地震波场数值模拟法

本发明涉及一种基于小波域Pre‑GSOR迭代的地震波场数值模拟法，先在背景波数中引入衰减因子并改写Helmholz方程，根据表示定理获得复波数L‑S方程并改写为归一化等效L‑S方程；将经法离散的归一化等效复波数L‑S方程对应的离散方程组变换到小波域中，经阈值处理后，生成新的稀疏线性方程组；利用预条件广义逐次超松弛迭代方法求解归一化等效复波数L‑S方程对应的小波域稀疏线性方程组，基于块ILU分解的ARMS分解方法生成非对角稀疏预条件矩阵，有效改善强散射介质下迭代级数的收敛性。本方法解决了广义逐次超松弛迭法在求解强散射介质下的L‑S方程离散线性方程组时，收敛速度缓慢或停滞的问题。