摘要
有限差分法算法是声波数值模拟算法中最广泛的数值计算方法。该方法具有计算速度快、占用内存相对较小、易于编程实现及模拟精度高等优点。规则网格下的有限差分算法相比于近年来提出的新的有限差分算法计算效率最高,消耗内存最少。然而,引入PML边界,导致常规网格的有限差分离散公式复杂,计算过程需要对计算区域和边界区域进行判断,导致计算效率低。针对这个问题,笔者新推导的离散公式,形式简单,整个模拟区域计算代码一致,可以很好地解决这个问题。OpenMP并行算法,语言简洁和可移植性高,通过结合OpenMP并行算法,对模拟算法进一步优化,可以较大提高数值模拟的计算效率。
Finite difference method is the most widely used method in acoustic numerical simulation,which has the advantages of fast computation,relatively small memory consumption,easy programming and high simulati`on accuracy.Compared with some new improved finite difference algorithms in recent years,finite difference algorithm based on regular grid is the most efficient one with the lowest memory consumption.However,the traditional discrete formula becomes complicated since the introduction of PML boundary.The calculation process needs to distinguish the calculation area from boundary area,which results in low calculation efficiency.To solve this problem,a new discrete formula is derived by the authors,which is simple in form and consistent in the whole simulation area.OpenMP parallel algorithm uses simple language and is with high portability.By combining with OpenMP parallel algorithm and further optimizing the simulation algorithm,the computational efficiency of numerical simulation is greatly improved.
作者
郑如秋
王波涛
冯永照
余卫江
ZHENG Ru-qiu;WANG Bo-tao;FENG Yong-zhao;YU Wei-jiang(Data Processing Company,China Oilfield Services Limited-COSL,Zhanjiang 542051,Guangdong,China)
出处
《世界地质》
CAS
2020年第4期896-904,共9页
World Geology
基金
国家科技重大专项项目(2016ZX05026-002-004)。
关键词
有限差分
规则网格
离散公式
OpenMP并行算法
计算效率
finite difference
regular grid
discrete formula
OpenMP parallel program
computational efficiency
