基于逐元算法的地震勘探高精度谱元法数值模拟

摘要

针对复杂构造地震勘探,常规正演模拟算法存在精度和计算效率低的问题,笔者研究并提出了基于逐元算法的高精度谱元法数值模拟方法.基于波动方程强、弱形式,推导二维波动方程的谱元解法,进而利用逐元算法,无需形成全局矩阵,仅对于独立的小规模单元结构进行计算,提高计算效率,减少存储空间消耗.应用基于逐元算法的高精度谱元法对复杂金属矿模型进行数值模拟,得出本方法计算精度高、数值频散小,计算效率高.%Aiming at solving the problem of low precision and low computation efficiency in conventional sim-ulation algorithm of the complex structure seismic exploration,the authors propose a high precision spectral element simulation based on element by element algorithm.Based on strong and weak forms of wave equation,the authors derivate the spectral element solution process of second-order seismic wave equation, followed by calculating the separate small-scale element structure using element by element algorithm without global matrix, which improves the computational efficiency and reduces storage space consumption.Numerical simulation of a complicated orebody model proves that this method gets high precision,small numerical dispersion and high computation efficiency.