  The close integration of increasing operating units into an integrated energy
system (IES) results in complex interconnections between these units. The
strong dynamic interactions create barriers to designing a successful
distributed coordinated controller to achieve synergy between all the units and
unlock the potential for grid response. To address these challenges, we
introduce a directed graph representation of IESs using an augmented Jacobian
matrix to depict their underlying dynamics topology. By utilizing this
representation, a generic subsystem decomposition method is proposed to
partition the entire IES vertically based on the dynamic time scale and
horizontally based on the closeness of interconnections between the operating
units. Exploiting the decomposed subsystems, we develop a cooperative
distributed economic model predictive control (DEMPC) with multiple global
objectives that regulate the generated power at the grid's requests and satisfy
the customers cooling and system economic requirements. In the DEMPC, multiple
local decision-making agents cooperate sequentially and iteratively to leverage
the potential across all the units for system-wide dynamic synergy.
Furthermore, we discuss how subsystem decomposition impacts the design of
distributed cooperation schemes for IESs and provide a control-oriented basic
guideline on the optimal decomposition of complex energy systems. Extensive
simulations demonstrate that the control strategies with different levels of
decomposition and collaboration will lead to marked differences in the overall
performance of IES. The standard control scheme based on the proposed subsystem
configuration outperforms the empirical decomposition-based control benchmark
by about 20%. The DEMPC architecture further improves the overall performance
of the IES by about 55% compared to the benchmark.
