Summary: Due to the non-locality of fractional derivatives, fractional partial differential equations were better to describe anomalous diffusion phenomena than other methods. However, while enjoying the convenience from mathematical modeling, it also caused lots of trouble especially in solving multidimensional cases. An efficient numerical algorithm was proposed for solving the three-dimensional space fractional advection diffusion equation (SFADE) by generalizing the Douglas-Gunn scheme. Stability and convergence of the method were proved by the matrix method. The derived alternating direction implicit (ADI) finite difference scheme had second order accuracy in both time and space directions, respectively. The efficiency and convergence orders were finally demonstrated by some numerical examples.