In this paper, an integrated approach to optimize electrical, natural gas, and district heating networks simultaneously is studied. Several interdependencies between these infrastructures are considered in details including a nonlinear part-load performance for boilers and CHPs besides the valve-point effect for generators. A novel approach based on selecting an appropriate set of state-variables for the problem is proposed that eliminates the addition of any new variable to convert irregular equations into a regular set while the optimization problem is still solvable. As a large optimization problem, the optimal solution cannot be achieved by conventional mathematical techniques. Hence, it is better to use evolutionary algorithms instead. In this paper, the well-known modified teaching-learning based optimization algorithm is utilized to solve the multi-period optimal power flow problem of multi-carrier energy networks. The proposed scheme is implemented and applied to a typical multi-carrier energy network. Results are compared with some other conventional heuristic algorithms and the applicability and superiority of the proposed methodology is verified. (C) 2015 Elsevier Ltd. All rights reserved.