Abstract: A D1Q3 evolution model of lattice Boltzmann method (LBM) is established to numerically solve a class of spatial fractional convection-diffusion equation in Riesz sense. By discretizing the integral term of fractional order operator, the fractional convection-diffusion equation is transformed into a standard one with Riesz derivative. With Taylor expansion, Chapman-Enskog and multi-scales expansion, equilibrium distribution functions of the model are derived in all directions. Furthermore, the macroscopic equation to be solved is recovered correctly. Finally, numerical experiments are carried out to verify the model.