An energy-preserving finite difference method is first presented for solving the nonlinear coupled space-fractional Klein-Gordon (KG) equations. The discrete conservation law, boundedness of the numerical solutions and convergence of the numerical schemes are obtained. These results are proved by the recent developed fractional Sobolev inequalities, the matrix analytical methods and so on. Numerical experiments are carried out to confirm the theoretical findings.