一种基于离散小波矩量法的地震散射波场数值模拟方法

摘要

本发明涉及一种基于离散小波矩量法的地震散射波场数值模拟方法；首先，选取形式较为简单的正交分域基函数，对待求波场进行展开，通过矩量法将LS方程转化为矩阵方程；其次，选取具有较好消失矩、正则性和紧支撑的Daubechie小波函数，利用Mallat塔式分解算法得到Daubechie小波滤波器，进而构造出正交小波变换矩阵对系数矩阵进行稀疏化，即得到关于小波展开系数的线性方程组；接着采用共轭梯度法求解该线性方程组；最后，对线性方程组的解向量进行小波逆变换，即得到待求波场的展开系数，进而求得计算区域上的待求波场值。本方法解决了传统矩量法在解LS方程时，离散系数矩阵满秩导致占用存储空间大、计算效率低的问题。