LAUSR

Abstract The dynamic stress strain simulation is one of the several digital rock physics methods that calculate elastic property of digital rocks. Simulations for calculating P- and S-wave property of one digital rock includes one compressional simulation and three shear simulations. Moreover, the simulation frequency must be low enough to meet the long wavelength assumption. Hence, heavy computational burden is unavoidable. To improve the computational efficiency, in addition to parallel computation, one possible way is to optimize the boundary condition applied on the digital rocks. The idea is to integrate the four kinds of boundary condition which correspond to the four separated simulations in the existed method into one boundary condition, accordingly the number of simulation reduces to one. The computational efficiency of the new simulation method improves approximately four times, and the numerical accuracy is proved to be unaltered.