In this paper, a massive MIMO enabled backhaul model for two-tier heterogeneous networks with in-band full-duplex transmission and imperfect CSI is provided. Because of a massive number of antennas at the macro base station and densification of small cells, circuit power consumption increases in the system. It is requisite to study energy efficiency of such backhaul networks. This paper aims at maximizing EE for optimal user association, spectrum allocation, and power allocation while accounting for QoS and backhaul constraints. A joint optimization problem is formulated as a non-convex mixed integer non-linear problem which provides an un-affordable computational complexity. To crack the problem efficiently, it is progressively divided into sub-problems until each one is convex and is solved separately to obtain the optimal solution. Furthermore, complexity of the proposed algorithm is evaluated and Jain’s fairness index is measured. Simulation results of the proposed distributive algorithm are shown and compared with that of the existing algorithm (Max SINR) to manifest the strength of the proposed algorithm. The impact of self-interference cancellation factor and channel estimation error on EE is clearly reflected by the numerical results. In addition, it is observed that using the proposed algorithm results in better EE than increasing the number of antennas.