This paper focuses on the hierarchical power distribution optimization of multiple local energy network (LEN) systems that are formed in three levels and can be operated in six typical modes. The decentralized optimal model for each LEN (the first level) and LENs (the second level) as well as the concentrated optimal model for the top level of the system are built, respectively. For each LEN, all the basic unities such as power generated by wind turbines and photovoltaic, and their upper nodes are considered. For LENs, the aggregated results (e.g., supply-demand requirements) from each LEN are dispatched. Furthermore, in the concentrated optimal control model (the third level), the ultimate supply-demand requirements of networked LENs together with other resources such as electric vehicles are considered. Due to the large number of control resources, the whole system is formulated as a large-scale global optimization (LSGO) problem. The self-adaptive differential evolution with neighborhood search method (SaNSDE) modified with the Lagrange multiplier method is used to solve the problem. The algorithm is firstly examined on 10 constrained benchmark functions, then it is applied to our problem. Experimental results show that the proposed model and algorithm are effective and efficient.