A Low Overhead Heterogeneous Parallel Optimization Method Based on Three-Dimensional Elastic Wave Numerical Simulation

Applying the staggered grid finite difference method (SGFDM) for simulating acoustic responses in large-scale, complex, three-dimensional models poses substantial challenges in geophysics, especially in high-resolution stratigraphic model, due to the high computational burden. To address this problem, we have proposed a low overhead parallel optimization method (LOPOM) suitable for heterogeneous architectures, using high-resolution borehole models as examples. The LOPOM enhances computational intensity and optimizes memory bandwidth utilization. This is achieved through the technique of data reuse along the discontinuities of the model and by minimizing such discontinuities within the halo region. LOPOM was implemented and tested on the CUDA platform and the Sunway supercomputer. In each instance, LOPOM demonstrated an optimal acceleration ratio, thus proving its robust performance. Furthermore, the method’s effectiveness was confirmed through its application to fracture-vuggy formation models and digital core models. Finally, the numerical simulation results and the actual logging data are combined to illustrate the application value of LOPOM.