A more practical and economic optimization mathematical model of load distribution for industrial boilers is presented. The optimum load rates of boilers are advanced by the method of static boiler model. The target function is built up, and in particular the restrictions are established. The theory of load assignment is checked through a project case. For various industrial boilers, once the definite characteristic functions of fuel consumption are given, the optimal load assigmnent can be obtained using this optimization model.