This work proposes a novel approach for the simultaneous synthesis of Heat Exchanger Networks (HEN) and Utility Systems of chemical processes and energy systems. Given a set of hot and cold process streams and a set of available utility systems, the method determines the optimal selection, arrangement and design of utility systems and the heat exchanger network aiming to rigorously consider the trade-off between efficiency and capital costs. The mathematical formulation uses the SYNHEAT superstructure for the HEN, and ad hoc superstructures and nonlinear models to represent the utility systems. The challenging nonconvex MINLP is solved with a two-stage algorithm. A sequential synthesis algorithm is specifically developed to generate a good starting solution. The algorithm is tested on a literature test problem and two industrial problems, the optimization of the Heat Recovery Steam Cycle of a Natural Gas Combined Cycle and the heat recovery system of an Integrated Gasification Combined Cycle. (C) 2017 Elsevier Ltd. All rights reserved.