For differential equations with multiple order spatial derivatives, there are some shortcomings by the classical high order compact (HOC) discretization. At least one of them is to reduce the computational efficiency due to the multiple inverse manipulations of matrixes. This motivates us to design a new kind of compact method what is called combined high order compact (CHOC) methods. The basic idea lying in this kind of method is to solve all the spatial derivatives simultaneously. Then, it is used to solve coupled nonlinear Schrödinger (CNLS) equations which contain both the first and second order derivatives. This scheme is not only more compact and accurate than standard HOC scheme and standard finite difference method with the same order, but also it can construct structure-preserving schemes. It preserves the symplectic structure and mass, and sometimes energy and momentum. Numerical experiments indicate that the new scheme can simulate the CNLS equations very accurately and efficiently. The mass and momentum are exactly preserved. The energy is preserved in some especially cases.