The ideas of noncommutative space originate from the research on time-space coordinate on an extremely small scale. Subsequently, the noncommutative space has gradually attracted some attention. The researchers started to explore noncommutative effect in some other fields. With the establishment of noncommutative quantum mechanics, it becomes significant to explore the noncommutative effect of exactly solvable models. The kinds of harmonic oscillators are very important and fundamental models in physics. But in noncommutative phase space, coordinate and coordinate are noncommutative, and momentum and momentum are also noncommutative. These results in the difficulty in obtaining the energy spectra of oscillators systems. In this paper the quantum properties of a two-dimensional coupling harmonic oscillator in noncommutative phase space are studied. Firstly, the Hamiltonian of the system is constructed, which includes all possible coupling types, namely, coordinate-coordinate coupling, momentum-momentum coupling, and coordinate-momentum cross-coupling. Secondly, the explicit expression of the noncommutative energy spectrum for the Hamiltonian is obtained by using the invariant eigen-operator method. In this work it is shown explicitly that the changes in the energy levels are related to the noncommutative parameters and coupling parameters. Thirdly, the effects of coupling parameters and non-commutative parameters on the energy spectra are analyzed in detail in the form of graphs. The results show that the energy levels under the influence of non-commutative parameters are non-degenerated. As the values of non-commutation parameters $ \theta $ and $ \phi $ increase, some energy levels increase and tend to change linearly, and other energy levels first decrease and then increase. If the limit values of the non-commutative parameters are taken as follows: $ \theta \to 0 $ and $ \phi \to 0 $, then the noncommutative energy spectra will be consistent with the energy spectra of the two-dimensional harmonic oscillator in the commutative space in general. On the other hand, the energy levels will split under the influence of coupling parameters. Moreover, the degree to which the energy levels split can increase as the kinds of couplings in the system increase. It is found that the coordinate coupling parameter $ \eta $ and the momentum coupling parameter $ \sigma $ have the same influence on the energy levels, but the coordinate momentum cross-coupling parameter $ \kappa $ has less influence on the energy levels than $ \eta $ and $ \sigma $. Overall, the above results are completely different from those of two-dimensional oscillator in the usual commutative space, which is degenerated except for the ground state.