A robust power allocation is proposed for downlink non-orthogonal multiple access (NOMA) heterogeneous networks with EH (Energy harvesting) under imperfect channel state information (CSI). In order to achieve green communication, an EH-aided scheme by leveraging energy from macro base station (MBS) signal and interference signal transmitted from other SBSs is proposed, which reduces the power burden and energy consumption of the SBS. In order to conform to the actual communication scenario, we construct an energy efficiency optimization function under imperfect CSI with considering the constraint of the outage probability interference power in macro cell user (MCU). However, the formulated optimization problem is non-convex due to the fractional form of the objective function and the probabilistic constraints of the outage probability limit. To cope with this problem, we propose a robust power allocation scheme. Firstly, the probabilistic problem is converted into a robust non-probabilistic problem by the minimax probability machine (MPM) and robust optimization theory. Then, the robust non-probabilistic problem can be transformed into the convex optimization problem via Dinkelbach method and sequential convex programming. Finally, the optimal transmission powers of the small cell users (SCUs) are obtained by Lagrange dual approach. The simulation results show that the robust power allocation scheme for NOMA heterogeneous networks with EH under imperfect CSI can significantly improve energy efficiency compared with traditional power allocation algorithms.