Abstract: In this study a selective filtering non-staggered finite difference method, called SFFD, is introduced to simulate seismic wave propagation in 2D media. SFFD utilizes optimized 11-point DRP (dispersion relation preserving) non-staggered finite difference scheme to discretize first-order spatial derivatives. In addition, selective filtering is applied to removing grid-to-grid oscillation caused by non-staggered algorithm. The selective filtering enhances the numerical accuracy and makes the simulation stable to implement. Test result demonstrates that SFFD can achieve the same accuracy as O(Delta;x8, Delta;t4) order staggered finite difference scheme. Moreover, the proposed algorithm is able to handle the media with high Poissonrsquo;s ratio. The media with strong inhomogeneity can also be treated by SFFD. As a non-staggered method, SFFD has potential application to seismic wave simulation in curvilinear coordinate system and general anisotropic media, in which complex interpolation must be performed for staggered scheme.