Multienergy systems (MES), in which multiple energy vectors are integrated and optimally operated, are key assets in low-carbon energy systems. Multienergy interactions of distributed energy resources via different energy networks generate the so-called distributed MES (DMES). While it is now well recognized that DMES can provide power system flexibility by shifting across different energy vectors, it is essential to have a systematic discussion on the main features of such flexibility. This article presents a comprehensive overview of DMES modeling and characterization of flexibility applications. The concept of “multienergy node” is introduced to extend the power node model, used for electrical flexibility, in the multienergy case. A general definition of DMES flexibility is given, and a general mathematical and graphical modeling framework, based on multidimensional maps, is formulated to describe the operational characteristics of individual MES and aggregate DMES, including the role of multienergy networks in enabling or constraining flexibility. Several tutorial examples are finally presented with illustrative case studies on current and future DMES practical applications.

One of the emergent prospects for active distribution networks (DN) is to establish new roles to the distribution company (DISCO). The DISCO can act as an aggregator of the resources existing in the DN, also when parts of the network are structured and managed as microgrids (MGs). The new roles of the DISCO may open the participation of the DISCO as a player trading energy in the wholesale markets, as well as in local energy markets. In this paper, the decision making aspects involving the DISCO are addressed by proposing a bilevel optimization approach in which the DISCO problem is modeled as the upper-level problem and the MGs problems and day-ahead wholesale market clearing process are modeled as the lower-level problems. To include the uncertainty of renewable energy sources, a risk-based two-stage stochastic problem is formulated, in which the DISCO's risk aversion is modeled by using the conditional value at risk. The resulting nonlinear bilevel model is transformed into a linear single-level one by applying the Karush-Kuhn-Tucker conditions and the duality theory. The effectiveness of the model is shown in the application to the IEEE 33-bus DN connected to the IEEE RTS 24-bus power system.