This paper focuses on a fast and high-order finite difference method for two-dimensional space-fractional complex Ginzburg-Landau equations. We firstly establish a three-level finite difference scheme for the time variable followed by the linearized technique of the nonlinear term. Then the fourth-order compact finite difference method is employed to discretize the spatial variables. Hence the accuracy of the discretization is $\mathcal{O}$(τ² + $h_1^4$ + $h_2^4$) in L₂-norm, where τ is the temporal step-size, both h₁ and h₂ denote spatial mesh sizes in x- and y- directions, respectively. The rigorous theoretical analysis, including the uniqueness, the almost unconditional stability, and the convergence, is studied via the energy argument. Practically, the discretized system holds the block Toeplitz structure. Therefore, the coefficient Toeplitz-like matrix only requires $\mathcal{O}$(M₁M₂) memory storage, and the matrix-vector multiplication can be carried out in $\mathcal{O}$(M₁M₂(log M₁ + log M₂)) computational complexity by the fast Fourier transformation, where M₁ and M₂ denote the numbers of the spatial grids in two different directions. In order to solve the resulting Toeplitz-like system quickly, an efficient preconditioner with the Krylov subspace method is proposed to speed up the iteration rate. Numerical results are given to demonstrate the well performance of the proposed method.