Abstract: To accurately describe the impact of the fluctuation and randomness of renewable energy output on the operation of regional integrated energy system, in this paper we propose a robust optimization of integrated energy systems based on correlated polyhedral sets. Firstly, a polyhedral set modeling is performed on the uncertain parameters of renewable energy output based on the traditional interval set. Subsequently, to address the issue of inability of the polyhedral set in accurately describing the spatiotemporal correlation between the uncertain parameters of renewable energy output, a correlated polyhedral set model is constructed by bending the polyhedral set boundary. In this way, a robust optimization model is established for the integrated energy system based on the correlated polyhedral set, and the column and constraint generation (C&CG) algorithm is employed to solve the model. Finally, the simulation results indicate that the robust optimization method for integrated energy system based on correlated polyhedral sets can effectively mitigate conservatism and enhance the robustness of optimization results by considering the correlation of renewable energy output.