The intelligence algorithm used in the backbone grid search is of weak robustness and has uncertain results. To tackle the drawback, this paper proposes a Mixed Integer Linear Programming (MILP) model for generating the power system backbone grid. With Kirchhoff's current law, the power balance of the power system, operation constraints for generators and lines are rigorously derived as analytical mathematical expressions. The distinguishing feature of the method is the connectivity of the backbone grid is expressed as some linear constraints with integer variables based on graph theory. Finally, optimal power flow is used to refine MILP results. The simulations of IEEE-39 bus and IEEE-118 bus show that the proposed method is effective and has high reliability, which provides a new idea for the development of this research.