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

针对复杂构造地震勘探,常规正演模拟算法存在精度和计算效率低的问题,笔者研究并提出了基于逐元算法的高精度谱元法数值模拟方法。基于波动方程强、弱形式,推导二维波动方程的谱元解法,进而利用逐元算法,无需形成全局矩阵,仅对于独立的小规模单元结构进行计算,提高计算效率,减少存储空间消耗。应用基于逐元算法的高精度谱元法对复杂金属矿模型进行数值模拟,得出本方法计算精度高、数值频散小,计算效率高。

High precision spectral element numerical simulation of seismic exploration based on element by element method

Aiming at solving the problem of low precision and low computation efficiency in conventional simulation 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.