Applying the perturbation method, the nonlinear evolution equations are expanded as multi-order approximate equations. And based on Lame equation and Lame function, these multi-order approximate equations can be solved by Jacobi elliptic function expansion method, where multi-order exact solutions of nonlinear evolution equations are derived. Then the invariants of the multi-order exact solutions are found among different nonlinear evolution equations.