A novel explicit optimized scheme for numerical simulation of elastic-wavefield separation

Abstract Numerical simulation of elastic-wave equation helps us better understand the information of underground structures and elastic-wave imaging has attracted the widespread attention of researchers. Using elastic-wave imaging requires separating the compressional and shear wavefields. Therefore, we develop a novel explicit optimized scheme to simulate the separated elastic wavefield. We construct a kind of 1-norm objective function directly utilizing the dispersion error and employ the simulated annealing algorithm to acquire improved finite-difference operators, whose optimal coefficients can effectively suppress spatial numerical dispersion. Meanwhile, we introduce a rotated staggered-grid (RSG) approach to enhance computational stability. Then, our proposed scheme, called the optimized RSG approach, is applied to the elastic-wave equations and decoupled elastic-wave equations to simulate the decoupled compressional and shear wavefield propagation. Numerical dispersion analysis is consistent with numerical results. The waveform comparison shows that the optimized RSG approach possesses higher accuracy, and several complex models are used to validate the applicability and effectiveness of the presented scheme.