Next Article in Journal
Operational Emissions in Prosuming Dwellings: A Study Comparing Different Sources of Grid CO2 Intensity Values in South Wales, UK
Next Article in Special Issue
Efficient Operations of Micro-Grids with Meshed Topology and Under Uncertainty through Exact Satisfaction of AC-PF, Droop Control and Tap-Changer Constraints
Previous Article in Journal
Optimized Configuration and Operating Plan for Hydrogen Refueling Station with On-Site Electrolytic Production
Previous Article in Special Issue
Robust Scheduling of Networked Microgrids for Economics and Resilience Improvement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Operation of a Power Grid with Embedded Networked Microgrids and Onsite Renewable Technologies

1
Department of Marketing and Business Analytics, San Jose State University, San Jose, CA 95192, USA
2
Department of Systems and Industrial Engineering, University of Arizona, Tucson, AZ 85721, USA
*
Author to whom correspondence should be addressed.
Energies 2022, 15(7), 2350; https://doi.org/10.3390/en15072350
Submission received: 23 February 2022 / Revised: 15 March 2022 / Accepted: 18 March 2022 / Published: 23 March 2022
(This article belongs to the Special Issue Operational Optimization of Networked Microgrids)

Abstract

:
The international community has set ambitious targets to replace the use of fossil fuels for electricity generation with renewable energy sources. The use of large-scale (e.g., solar farms) and small-scale solutions (e.g., onsite green technologies) represents one way to achieve these goals. This paper presents a mathematical optimization framework to coordinate the energy decisions between the distribution network and the networked microgrids embedded within it. Utility-scale renewable and conventional generators are considered in the distribution network, while the microgrids include onsite renewable generation and energy storage. The distribution network operator utilizes demand-side management policies to improve the network’s efficiency, and the microgrids operate under these programs by reducing their energy usage, scheduling the electricity usage under dynamic tariffs, and supplying energy to the grid. The uncertainty of renewable energy sources is addressed by robust optimization. The decisions of the distribution network and the microgrids are made independently, whereas the proposed collaboration scheme allows for the alignment of the systems’ objectives. A case study is analyzed to show the capability of the model to assess multiple configurations, eliminating the necessity of load shedding, and increasing the power supplied by the microgrids (22.3 MW) and the renewable energy share by up to 5.03%.

1. Introduction

According to the World Bank [1], between 1960 and 2018 the world’s total greenhouse gas emissions increased from 9.46 to 34.04 Gt of CO 2 and the emissions per capita increased from 3.12 to 4.48 metric tons of CO 2 . Additionally, electricity consumption increased from 1200.15 to 3128.30 kWh per capita between 1971 and 2014 [1]. To satisfy the growing energy demand, efficient electricity supply systems with a higher capacity are necessary. Nevertheless, a linear increment in the current power generation share, which is mainly composed of fossil fuel sources, will result in an accelerated increase in greenhouse gas concentration. For example, according to the U.S. Energy Information Administration (EIA) [2], the U.S. power grid in 2020 was composed of 60.6% fossil fuels, 19.7% nuclear, and 19.8% renewables. The use of renewable energy sources (RES) to generate electricity is an option to reduce greenhouse gas emissions. As a result, there is a growing environmental awareness in the international community, and in some cases, specific goals have been derived to address this climate change, e.g., the United Nations Sustainable Development Goals [3]. Other international efforts, such as the Paris Agreement, provide quantifiable targets to reduce greenhouse gas emissions or to increase the renewable share of the generation mix. For example, the European Union targets are to generate 27% of the electricity from renewables by 2030, while the U.S. state of California’s targets are 50% by 2030 [4].
Reducing the environmental footprint from global-scale systems requires the participation of both large-scale and small-scale entities. In this context, the utilization of renewable energy to satisfy the electricity demand is aligned with the international community’s objectives. Increasing the renewables share in the electricity generation mix can be achieved by two different approaches: (1) By the utilization of large-scale generators, such as wind farms, hydropower, solar farms, and solar concentration; (2) By using multiple small-scale generators owned by the consumers for onsite electricity generation. Each of these approaches presents different advantages and disadvantages. Large-scale generators owned by utility companies allow some control in their operations and more precise forecasts, while from the consumer point of view there will be more to pay in utility bills due to the investment required for these projects and the cost associated with the intermittency. In the case of small-scale solutions, the utility company reduces its control and forecast capabilities over these sources, but it transfers the investment costs to the consumers, while the owners can actively participate in the electricity market due to demand-side management (DSM) policies, the capacity to form microgrids—discussed later in this paper—and reduced electricity bills, even if there is an initial investment. In any case, one of the main challenges of renewable sources compared to conventional generation systems is that they present a higher variability which can jeopardize the grid stability and reliability [5,6], and more advanced control techniques are required to efficiently manage this uncertainty.
On the energy policy side, DSM is a series of activities that distribution network operators (DNOs) undertake to change the users’ energy consumption pattern to achieve power grid efficiency. For a comprehensive review of DSM policies, see [7,8,9]. DSM policies can be divided into two categories: energy efficiency and demand response (DR). While energy efficiency proposes policies to reduce the energy consumption, e.g., the replacement of older equipment with efficient equipment, DR utilizes activities in which the consumers can participate to change the load pattern [9]. Laitsos et al. [10] mention that the main advantages that DSM policies provide to energy systems include (i) peak shaving, (ii) reduction in operational and capacity costs, (iii) improved system reliability, and (iv) emission reduction via the utilization of greener technologies. These advantages are the result of the consumer’s load curve control, carried out in six modalities [9,10]:
  • Peak shaving: This aims to reduce the power usage during high-demand periods (peak), either via direct load control or through pricing contracts.
  • Valley filling: The goal of this strategy is to encourage consumers to increase their load during periods with low demand (off-peak) to secure the system’s stability.
  • Load shifting: This consists of shifting the consumers’ loads from peak periods to off-peak periods with minimum or no change in the total energy consumption.
  • Load growth or load building: The overall consumption increases during the planning horizon.
  • Load reduction or energy conservation: The overall consumption during the planning horizon is reduced.
  • Flexible load curve: The loads can be redistributed during the planning horizon flexibly. The consumers receive benefits for participating in this strategy, while the DNOs have the flexibility to improve the reliability of the system or to respond to emergencies.
In general, consumers are interested in two types of DR policies: incentive-based programs, in which the users receive incentives to participate in specific events defined by the grid operator, and price-based programs, which present different electricity tariffs at certain periods depending on the demand level of the grid and the generation cost. The basic price-based program is the Time-of-Use, in which different price tariffs are present in different time blocks; the tariffs are higher during peak periods than during off-peak periods. A modification of this approach, known as Time-and-Level-of-Use (TLOU) [11] also considers the level of consumption at each time block. Other price-based programs are the real-time pricing, in which the tariffs are based on the real-time generation cost, and the critical peak pricing, which imposes very high tariffs during (very short) critical peak periods, e.g., during emergency conditions [8].
From the consumer side, the evolution in distributed energy resources, defined as demand- and supply-side resources that can be deployed throughout an electric distribution system to meet the energy and reliability needs [12], has changed the nature of the grid from passive to active [13]. These resources, such as roof solar photovoltaic (PV) panels and diesel generators with capacities ranging up to hundreds of MW, small-scale storage cells, electric vehicles, and smart meters, among other technologies, provide the consumers with the opportunity to generate and store energy. Different from utility-scale energy resources, such as coal-based thermal power plants, hydropower, and large-scale solar and nuclear plants, which are centralized and designed to supply energy for long distances, distributed energy resources are located close to the consumption points. This reduces the transmission and generation costs, particularly when the distributed resources are utilized to cover peak periods [14].
Distributed energy technologies allow for the formation of microgrids. A microgrid is defined as a small-scale power system with at least one distributed energy resource and one load, with clearly defined electrical boundaries, and a smart infrastructure that allows self-supply and islanding, providing security, reliability, autonomy, and resilience [14,15]. Microgrid technology lets electricity consumers improve their energy efficiency, including the use of clean energy. Some examples of microgrids include small distribution networks, manufacturing systems, data centers, smart buildings, and other systems that have onsite energy resources. For a comprehensive discussion about microgrids, including architecture, functions, challenges, and opportunities, see the research of Yoldaş et al. [16]. As discussed by Khodaei [15], microgrids offer improved reliability, resiliency, and efficiency, reduction of carbon emissions due to a higher power quality, the utilization of less costly RES, and supply of loads in remote areas. The management of a microgrid can be more sophisticated when compared to the distribution network (DN) due to their significantly smaller scale, allowing the potential to mitigate the effects of resource uncertainty and contingencies [17], particularly when load components are deferrable at some level [18]. By utilizing information and communication technologies, microgrids can coordinate the distributed energy resources and controllable loads efficiently, as well as switch over to islanded mode and operate autonomously for safety purposes [19,20]. Controllable distributed generators within a microgrid are much smaller than the utility-scale ones of the DN, allowing easier switching operations and a more flexible scheduling that can mitigate the impact of the intermittency of the RES [21,22]. Microgrids are connected to the DN through the point of common coupling, a special bus that allows electricity flow in both directions. A sophisticated energy management system (EMS) to monitor, control, and optimize the (onsite and grid) energy resources and find the system’s energy balance contributes to the above-mentioned benefits, making the microgrid technology attractive for the operation of diverse systems. For a comprehensive review of microgrids’ EMSs, see the research of Battula et al. [14].
The existence of microgrids that are capable of advanced energy-related decisions within the DN allows for the efficient implementation of DR policies. One of the main barriers that prevents achieving this implementation is the lack of effective collaboration and negotiation protocols between the DNO and the operators of the microgrids, considering that the former obtains profits for selling electricity to the latter, while the microgrids seek to minimize their operational costs. This paper proposes (i) a mathematical optimization framework to find the optimal operation of a DN with networked microgrids embedded into it and (ii) a collaboration scheme to coordinate the energy decisions between the DN and the microgrids. It is assumed that the DNO works closely with the operators of large-scale generators that power the network. Therefore, the decisions to be made include the unit commitment of the utility-scale generators, the optimal power flow problem for the DN, the incentives offered to the microgrids, and the bidirectional power flow between the DN and the microgrids to satisfy the network’s energy demand. The DNO uses a dynamic pricing scheme, and requests load reductions to the microgrids in exchange for economic incentives. The proposed model allows the microgrids to make their own decisions, considering the dynamic tariff policies and energy consumption reductions requested by the DNO. Additionally, the point of common coupling allows net metering; i.e., the microgrids can sell their power surplus to the DN. Each microgrid includes an EMS that takes into account the energy supply from the DN, the uncertainty related to the onsite renewable generators, the utilization of energy storage, and the DR policies to find the decisions that minimize its operational cost. The DN can utilize the microgrids’ energy resources and decisions as a reserve to improve the network’s flexibility. The proposed algorithm helps to carry out the collaboration between the DNO and the microgrids in terms of the interchange of energy, and the participation of the microgrids in the DR programs.
The main contributions of this paper are: (i) We introduce an optimization model that finds the optimal unit commitment and power flow of a DN with multiple embedded microgrids, without direct control nor complete information of the microgrids’ operations; the uncertainty related to the renewable sources is quantified and incorporated into the model via robust optimization. (ii) A collaboration and negotiation protocol between the DNO and the microgrids to allow the latter to reduce their energy usage, schedule the electricity usage under dynamic tariffs, and supply energy to the grid. (iii) A case study based on a modified instance of the IEEE 30-Bus System is analyzed to show the effectiveness of the proposed approach under different configurations, eliminating the necessity of load shedding and increasing the renewable energy share and power supplied by the microgrids to the distribution network.
The remainder of this paper is organized as follows: Previous work is reviewed in Section 2. In Section 3, the model for the optimal operation of the DNO with multiple embedded microgrids is presented, including the discussion of the uncertainty management and solution algorithm to solve the proposed model and the negotiation protocol. In Section 4, a case study and sensitivity analysis are presented, and the results are discussed in Section 5. Finally, Section 6 presents the conclusions.

2. Literature Review

DSM policies have been studied widely, and algorithms have been proposed to optimize the operation of diverse systems under these conditions. Habibian et al. [23] study the ability of an industrial consumer of electricity to flexibly reduce demand and offer interruptible load reserve and its impacts in the market. Anjos et al. [11] propose a framework to optimize the operation of an electricity retailer that sells energy using a TLOU pricing structure, finding the optimal price in this scheme considering a competitive environment and flexible loads from the consumers. Laitos et al. [10] present an implementation of DR emergency programs to induce the energy consumption of a power system’s consumers, considering renewable energy sources via particle swarm optimization. Mohammad et al. [24] propose an EMS for residential buildings with onsite PV, energy storage, and electric vehicles utilized as another storage system while minimizing the operating cost, peak to average ratio, and the user discomfort (including the electric vehicles’ availability). Zator [25] presents a case study to measure the impact of the utilization of an EMS for homes with PV and heat and electric energy storage system installations on self-consumption of a single-family house with price-based DSM. Kumar et al. [26] propose a method to integrate customer-oriented and utility-oriented DSM policies into microgrids’ EMS, considering non-dispatchable energy resources. Vergara-Fernandez et al. [27] present a case study for a seawater reverse osmosis desalination plant, whose operation is optimized using mixed-integer linear programming (MILP) that includes DSM policies. Baldi et al. [28] propose an optimal control framework for heat, ventilation, and air conditioning units in interconnected commercial buildings equipped with PV generators, achieving DR management and thermal comfort conditions for users. On the other hand, Korkas et al. [29] propose an optimal control algorithm for demand response and thermal comfort management for microgrids with onsite renewables (namely PV and wind turbines) and energy storage, considering the thermal requirements based on the occupancy level of the buildings under study. The framework performs a local optimization for each building as well as an aggregate level management carried out by a central controller. Gomez-Herrera and Anjos [30] present an optimization scheme for the operation of residential smart buildings, balancing cost and users’ comfort while participating in DR programs. Ruiz Duarte et al. [31] present an optimization framework for manufacturing facilities including onsite energy generation and storage technologies, operating under DSM policies. The approaches presented in the above-mentioned research present opportunities to the DN in which these systems are embedded to increase the network’s efficiency.
The optimization of networked microgrids has also been studied. Wang et al. [32] develop a mathematical model that captures heterogeneous management systems with flexible demands to find the optimal operation of a centralized DN with embedded networked microgrids with electricity consumer and producer characteristics. Fang et al. [33] present a cooperative energy dispatch approach to coordinate multiple networked autonomous microgrids with intermittent renewable generators. In a similar approach, Fang et al. [34] propose a collective energy dispatch approach to coordinate distributed generation, energy storage, and critical demand across multiple autonomous microgrids, considering small-scale distributed renewable generators. In their research, Wang et al. [35] develop a mathematical programming algorithm to optimally coordinate the operation of networked microgrids and the DNO, each of them with different objectives, considering uncertainties related to RES. Similarly, Wang et al. [36] propose a decentralized management system for networked microgrids, each with individual objectives and operating autonomously, considering both grid-connected and islanded modes. Wang et al. [37] present an algorithm for the operation of networked microgrids, which are connected through a physical common bus, to find the optimal schedule for distributed generation, energy storage, and controllable loads of each microgrid. Kou et al. [38] propose a distributed economic model to achieve coordination among multiple microgrids under a stochastic energy management approach for multiple microgrids and the DNO. Lv and Ai [39] propose a dynamic energy management strategy for an energy system consisting of multiple microgrids, with interactions between the DNO and the microgrids, and among the microgrids, with large-scale integration of solar and wind RES. Haddadian and Noroozian [40] develop an approach for optimizing a DN, first in an integrated approach, and after that decomposing the DN into a set of networked microgrids, considering RES and energy storage. Wang and Wang [41] propose a networked microgrid framework for outage detection and service restoration via network reconfiguration and power support from the microgrids. Ma et al. [42] propose a distributed algorithm for online energy management for networked microgrids with high RES integration utilizing regret minimization and alternating the direction method of multipliers. Ferro et al. [43] propose a bi-level optimization framework for an aggregator that manages different local consumers and/or producers (prosumers). The DN under study include the participation of microgrids and smart buildings as prosumers with renewable sources and storage systems, which can receive incentives for DR participation in a network topology. The aggregator, represented as the higher-level stage, determines the reference values that the prosumers try to reach at the lower-level stage. The main differences between the presented literature review and this paper are (i) the assumption that each microgrid optimizes its operation and energy usage individually, resulting in electricity demand/surplus; and (ii) the negotiation protocol assumes that the DNO does not have direct control nor complete information of the microgrids’ internal systems but requires a day-ahead estimate of the energy usage, availability, and decisions about participating in DR programs.

3. Materials and Methods

The system under study consists of a DN that requires finding the optimal unit commitment to satisfy its demand. It consists of a set of large-scale conventional generators, G (indexed by g), and a set of large-scale renewable generators, R (indexed by r), that provide electricity to a set of buses, I (indexed by i), of which some are microgrids, M I , (indexed by m), and a set of transmission lines, E (indexed by e), that distribute the electric energy. The DNO requires electric energy to be provided to satisfy the demand of set I , by either using the large-scale generators G and R , or leveraging on the onsite production of the microgrids, during a planning horizon T (indexed by t). The electricity is supplied under a dynamic tariff (DT) policy in which the prices per unit of electricity are adjusted according to K levels of consumption (indexed by k). Additionally, the DNO can request that the microgrids reduce their loads in exchange for economic incentives. Each microgrid m M can satisfy its demand via a combination of onsite energy generation and storage technologies, as well as energy purchased from the DN. The microgrid can decide to accept a reduction in their energy consumption levels when the DNO requests them. Furthermore, the microgrid has the capability of selling the surplus of onsite produced energy to the DN. Each microgrid has different characteristics depending on the system it represents, and its goal is to find its optimal operational and energy decisions. The microgrids under study include residential facilities, manufacturing facilities, data centers, and agricultural facilities. The system is shown in Figure 1, in which the large-scale generators are on the left-hand side of the figure, the non-microgrids buses are in the middle, and the microgrid buses are on the right-hand side.
The mathematical formulation is presented below. Table 1, Table 2 and Table 3 summarize the nomenclature utilized in the formulation. In particular, Table 1 lists the sets and indices for the time horizon, buses, microgrids, transmission lines, conventional and renewable generators, and dynamic electricity tariffs. Table 2 presents the costs associated with unit commitment, electricity generation and ramping, electricity purchased, and load shedding; the sales price of electricity; electricity generation and ramping capacities, as well as electricity trades limits; capacity and susceptance of transmission lines; and electricity demand of the buses. Table 3 lists the decision variables, including unit commitment, power generation, and ramping, energy flow in the transmission lines, power sold to and bought from microgrids, load reductions, and load shedding.

3.1. The Distribution Network

The DN needs to find the day-ahead schedule of the conventional generation units, the power used from the renewable generators, the optimal power flow, the energy transactions with the microgrids, the application of DR policies, and the energy balance that minimize the cost.

3.1.1. Unit Commitment

Since utility-scale conventional generators require larger periods to change their status, the problem of finding their optimal scheduling problem has been well studied. The unit commitment problem refers to the decisions of starting up and shutting down generators. A modified version of the unit commitment problem presented in [44] is shown in Equations (1)–(8).
t = 1 T g u 0 ( 1 x g , t ) = 0 , g G
t = 1 T g d 0 x g , t = 0 , g G
t = t t + T g u 1 x g , t T g u · y g , t u , g G , t { T g u 0 + 1 , , T T g u + 1 }
t = t t + T g d 1 ( 1 x g , t ) T g d · y g , t d , g G , t { T g d 0 + 1 , , T T g d + 1 }
t = t T ( x g , t y g , t u ) 0 , g G , t { T T g u + 2 , , T }
t = t T ( 1 x g , t ) y g , t d 0 , g G , t { T T g d + 2 , , T }
y g , t u y g , t d = x g , t x g , t 1 , g G , t T
x g , t , y g , t u , y g , t d { 0 , 1 } , g G , t T
Constraints (1) and (2) represent the initial online/offline requirements for generation unit g, which consists of T g u 0 and T g d 0 periods, respectively. The minimum number of periods that the generator unit g must remain online/offline after the minimum initial online/offline periods is represented by (3) and (4) ( T g u and T g d periods, respectively). These constraints are modified to account for the minimum number of periods that the generator unit g must remain online/offline at the end of the time horizon, as presented on (5) and (6). The relationship between starting up/shutting down generator unit g, represented by binary variables y g , t u and y g , t d , respectively, and its commitment status, modeled as the binary variable x g , t , is modeled in (7). The binary requirements for the decision of the online/offline status as well as the starting up/shutting down the generator unit g are modeled in (8).

3.1.2. Conventional Power Generation

The amount of power that the utility-scale generators produce, modeled as p g , t G + , is bounded by the lower and upper limits ( P ̲ g G and P ¯ g G , respectively) determined by their technical capacities. Additionally, a generator can contribute only if it is turned on during a certain period. The power generation is shown in Equations (9)–(14).
P ̲ g G x g , t p g , t G + P ¯ g G x g , t , g G , t T
p g , t G + = p g , t 1 G + + p g , t u p g , t d , g G , t T
0 p g , t u R g u · x g , t 1 + R ˜ g u · y g , t u , g G , t T
0 p g , t d R g d · x g , t + R ˜ g d · y g , t d , g G , t T
0 p g , t G p g , t G + , g G , t T
p g , t G = p g , t G + p g , t G , g G , t T
A generator g can produce power only if it is online, as modeled in (9). Generator unit g can change its power supply depending on its ramping rates, which determine how fast the electricity generation can be increased ( p g , t u ) or decreased ( p g , t d ) per unit time, or when the generator is turned on/off; this is modeled in (10). The generator’s ramping rate limits are presented in equations (11) and (12), where R g u and R g d represent the increase/decrease ramping limits when the generator is online, whereas R ˜ g u and R ˜ g d represent the increase/decrease ramping rates when the generator has started up or is shut down, respectively. In the case that the power generation of unit g is more than the required amount, it is possible to curtail some of that power to achieve power balance; the power curtailed, modeled as p g , t G , cannot be larger than the power generated, as modeled in (13). The amount of power supplied to the DN is modeled as p g , t G in (14).

3.1.3. Renewable Power Generation

For utility-scale renewable generators, the renewable power, modeled as p r , t R , is considered to depend only on the renewable power availability, P ¯ r , t R . The amount of renewable power that can be generated is presented in Equation (15).
0 p r , t R P ¯ r , t R , r R , t T
For this paper, the available renewable power is treated as uncertain. The integration of this uncertainty in the proposed model is discussed in Section 3.3.

3.1.4. Transmission Lines

A network topology is utilized to model the DN. Transmission lines can be considered as the edges of the network, while the buses are analogous to the nodes. The power flow that is generated on a bus i can satisfy the demand on a bus j if there is a path that connects buses i and j through transmission lines. The initial node of transmission line e is modeled as i e and its end node is modeled as j e . The dynamics of transmission lines are modeled as the DC approximation for the optimal power flow problem, as presented in Equations (16)–(18).
F e · z e , t f e , t F e · z e , t , e E , t T
f e , t = B e · ( θ i e , t θ j e , t ) · z e , t , e E , t T
z e , t { 0 , 1 } , e E , t T
The flow in transmission line e, modeled as f e , t , is bounded by its physical capacity, F e , and the line switching status, z e , t , as modeled in (16). Kirchhoff’s law for the DC optimal power flow problem, which considers that the power flow in a transmission line within the microgrid is proportional to the difference of the voltage phase angles at the terminal buses of the line ( θ i e , t and θ j e , t ) and the line’s susceptance ( B e ), taking into account the decision of connecting/disconnecting the transmission line as modeled in (17). The connected/disconnected binary status of the transmission lines is modeled in (18). Note that (17) involves a nonlinear term. Then, it can be reformulated as presented in (19).
B e · ( θ i e , t θ j e , t ) ( 1 z e , t ) · F e f e , t B e · ( θ i e , t θ j e , t ) + ( 1 z e , t ) · F e , e E , t T

3.1.5. Electricity Bought from Microgrids

The DNO has to supply electricity to every bus in the network. In the case of microgrids, it is assumed that they have electricity production capacity through distributed generation. It can be the case that the microgrid produces more electricity than the requirements for its operation. In this case, the DNO can leverage this surplus of electricity to satisfy other buses in the network by purchasing it from the microgrids. This type of transaction can be carried out at different rates. The electricity traded under a contract limit is called firm, modeled as p m , t F , and any amount of electricity traded beyond the contract limit is considered non-firm, modeled as p m , t N . The trade prices for both firm and non-firm purchases are different. The constraints for the connection and trades between the DNO and the microgrids are modeled in the Equations (20)–(22).
k K p i m , k , t ( ) ( p m , t N + p m , t F ) = e E : j e = i m f e , t , m M , t T
0 p m , t F P ¯ m , t F , m M , t T
0 p m , t N P ¯ m , t N , m M , t T
The total power flow in the transmission lines between the DNO and microgrid m is modeled in (20), where the left-hand side represents the net trade calculated by subtracting the electricity supplied from the microgrid to the DN from the energy sold to the microgrid node, modeled as p i m , k , t ( ) . The limits of the power that can be traded under firm and non-firm contracts, P ¯ m , t F and P ¯ m , t N , are modeled in (21) and (22), respectively.

3.1.6. Demand Response

In this paper, two DR policies are considered: DT and load-reduction incentives. Particularly, the DT policy is modeled for the power sold to the buses. Therefore, the demand at each bus i is divided based on the | K | consumption levels of the DT policy and represented by D i , k , t . Moreover, if bus i represents a microgrid, the DNO may request a load reduction of q m , k , t units of electricity at a particular period t. If the DNO is still unable to satisfy that demand, it is possible to incur load shedding of s i , k , t units of electricity. The amount of electricity sold to bus i after the above-mentioned considerations is modeled in Equations (23)–(25).
p i , k , t ( ) = D i , k , t s i , k , t q i m , k , t , i I , k K , t T
0 q m , k , t Q ¯ m , k , t , m M , k K , t T
0 s i , k , t D i , k , t , i I , k K , t T
The electricity sold to bus i equals its demand minus the load that is shed, and, in the case of microgrid nodes, minus the amount that is requested to reduce, as modeled in (23). The maximum amount that can be requested to be reduced from microgrids at each DT level is modeled as Q ¯ m , k , t in (24). The load that can be shed from a certain bus cannot be greater than the demand of this bus, as modeled in (25).
Note that the design of the set of equations to model the power sold to each bus allows other pricing policies, such as real-time pricing and critical peak pricing, by modifying the price parameters on each time frame t, which are discussed in Section 3.1.8.

3.1.7. Power Balance

The power supply and demand balance should be satisfied at each of the nodes of the DN. Equation (26) represents the power balance at bus i for each time frame t. The power supply (in the left-hand side of the equation) is the sum of the power generated by conventional and renewable generators placed at bus i and the incoming power flow. The power demand (on the right-hand side of the equation) is the power sold to bus i and the sum of the outgoing power flow.
g : i g = i p g , t G + r : i r = i p r , t R + e : j e = i f e = k K p i , k , t ( ) + e : i e = i f e , i I , t T

3.1.8. Cost Management

The cost function consists of the production cost for conventional and renewable generators, the ramping and commitment processes of the conventional generators and load shedding, the revenues obtained for selling electricity under DT policy, as well as the incentives paid for load reduction requests accepted by the microgrids. In general, the unit commitment and power generation costs are nonlinear and nonconvex. A linear approximation of the cost function consisting of a piecewise-linear reformulation based on the AC optimal power problem flow is presented in this paper. The unit commitment cost associated with starting and shutting down the generators, c UC , is shown in Equation (27).
c UC = t T g G C g T u · y g , t u + C g T d · y g , t d ,
where C g T u and C g T d represent the cost for starting up and shutting down the conventional generators, respectively.
The power generation and ramping cost for both conventional and renewable generators, c PG , are shown in Equation (28).
c PG = t T g G C g P · p g , t G + + C g R u · p g , t u + C g R d · p g , t d + r R C r R · p r , t R ,
where C g P , C g R u , C g R d and C r R indicate the generation cost, ramping up, and raping down of conventional generators, as well as the generation cost of renewable generators, respectively.
The cost of the electricity purchased from the microgrids and the incentives paid for accepting load reductions, c μ G , are shown in Equation (29).
c μ G = t T m M C m , t F · p m , t F + C m , t N · p m , t N + k K C m , k , t I n · q m , k , t ,
where C m , t F , C m , t N and C m , k , t I n model the cost of energy purchased from the microgrids under firm and non-firm contracts, as well as the load reduction incentives per electricity unit, respectively.
The cost of load shedding from different buses of the DN, c Sh , is reflected in Equation (30). Note that, depending on the criticality level of a given bus, the load shedding penalty C i , k , t S can be adjusted, e.g., this parameter can be relatively low for a residential bus, and very high for a hospital.
c Sh = t T k K i I C i , k , t S · s i , k , t
Finally, the revenue obtained from selling energy at different DT levels, c DT , is shown in (31). Note that the definition of the DT parameters C i , k , t D T allows the model to be flexible to utilize different tariffs schemes, such as TLOU (used in this paper), real-time pricing (by using different values on each time frame), and critical peak pricing (by setting very high values in particular time frames).
c DT = t T k K i I C i , k , t D T · p i , k , t ( )

3.1.9. Optimization Model

The optimization model for the operation of the DNO consists of minimizing the total cost, subject to the operational constraints. In order to facilitate the discussion of the solution methodology, a compact form of the problem is presented in (32). Let vector x represent the unit commitment (binary) decision variables, x g , t , y g , t u , y g , t d ; vector y denote the power generation and power flow (continuous) decision variables, p g , t G + / G / G , p g , t u / d , p r , t R , f e , t , θ i e , t , p m , t F / N , p i , k , t ( ) , q m , t , and s i , t ; and vector z represents the line switching (binary) decision variables, z e , t .
min x , y , z c UC + c PG + c μ G + c Sh c DT s . t . Constraints ( 1 ) ( 16 ) , ( 18 ) ( 26 )

3.2. The Microgtrids

In this paper, it is assumed that each microgrid is capable of making its own operational and energy decisions. Optimal control techniques for microgrids are discussed in more detail by Minchala-Avila et al. [45]. It is out of the scope of this paper to get into the details for each type of microgrid but is assumed, without loss of generality, that their management team utilizes mathematical optimization to find their optimal operation. A general cost minimization model for the optimal operation and energy decisions of microgrids that considers DSM policies can be modeled as shown in (33)–(36).
min o , p | δ D S M f ( o , p )
s . t . g ( o , p ) 0
o O
p P ,
where o and p represent the operational and energy decisions, respectively, δ D S M represents the DSM policies used by the DNO, and constraints (34)–(36) represent the operational and energy constraints. In this paper, δ D S M is composed of the incentives offered for load reduction requests, C m , k , t I n ; the load reduction requested by the DNO, q ¯ m , k , t ; the DT levels, P ¯ i m , k , t D T ; the DT prices, C ¯ i m , k , t D T ; the price of the energy bought at firm and non-firm contracts, C m , t F and C m , t N , respectively; and the maximum amount of the firm and non-firm contracts, P ¯ m , t F and P ¯ m , t N , respectively.
To achieve the optimal coordination of operational and energetic goals, an EMS is assumed to be installed within each microgrid to jointly optimize the use of energy and its operations. Figure 2 presents the scheme for a microgrid with a central EMS that utilizes information from the operational requirements, as well as the available onsite energy technologies and the DN, finding the optimal operations.
It is assumed that each microgrid can generate and store its own electricity, and utilize the electricity from the DN to optimize its processes. The EMS considers information, such as costs and incentives offered by the DN, including the DSM policies, to optimally decide how much energy to purchase or sell to the DN. It is assumed that the microgrid’s EMS will produce and communicate the following information to the DNO:
  • The amount of energy consumed from the grid at each DT level, D i m , k , t .
  • The amount of energy sold to the grid under firm and non-firm contracts, p ^ m , t F / N .
  • The amount of electricity consumption willing to reduce, q ^ m , k , t .
The information obtained from each microgrid can be used by the DNO to find its optimal operation.
The microgrids considered in this paper are classified into three categories, based on the definition of energy-consuming sectors defined by the U.S. EIA [46]:
(1)
Residential microgrids: The residential sector consists of living quarters for private households. Residential microgrids include the operation of appliances such as space heating, water heating, air conditioning, lighting, refrigeration, cooking, and running a variety of other appliances. Some of the residential appliances allow certain flexibility in the time frames they need to be operated. As an example, a washing machine can be scheduled to operate during specific periods that optimize energy usage.
(2)
Industrial microgrids: The industrial sector consists of facilities and equipment used for producing, processing, or assembling goods. This sector encompasses agriculture, forestry, fishing, hunting, mining, oil and gas extraction, and construction. Industrial microgrids present flexibility in their operation as long as the demand of goods is satisfied during the planning horizon. A change in their operation results in a change in their energy consumption patterns.
(3)
Commercial microgrids: The commercial sector consists of service-providing facilities and equipment of businesses, governments, and other private and public organizations. Common uses of energy associated with this sector include space heating, water heating, air conditioning, lighting, refrigeration, cooking, and running a wide variety of other equipment. Commercial microgrids present a certain flexibility in their operation, as long as a certain quality of service level is satisfied. Similarly to the industrial microgrids, a change in their operation results in a change in their energy consumption patterns.

3.3. Solution Methodology

In this section, the solution approach for the DN alone is presented in Section 3.3.1, and the framework for the collaboration between the DN and the microgrids is presented in Section 3.3.2.

3.3.1. Optimal Operation of the Distribution Network

In mathematical optimization, uncertainty can be addressed by finding the optimal operation of the system under the most adverse scenario with robust optimization. It is suitable for assessments in the long term, e.g., to compare the performance of the system after applying certain DSM policies. It is also suitable for short-term decisions since it can be assumed that, if the system can perform successfully in the worst-case scenario, then it can perform successfully in any other scenario as well.
The uncertainty of the presented model comes from the renewable power availability, P ¯ r , t R , given that its intermittency depends on weather conditions. Nevertheless, it is assumed that the expected value, E [ P ¯ r , t R ] , and the variance, Var [ P ¯ r , t R ] , of the renewable power output at each period can be obtained. Then, an uncertainty set, U , can be constructed utilizing the mean and the standard deviation, without a loss of generality. For the purpose of this paper, the uncertainty set is constructed by a lower bound and an upper bound, defined as follows:
U = P ¯ r , t R [ E [ P ¯ r , t R ] Y ̲ r , t , E [ P ¯ r , t R ] + Y ¯ r , t ] , r R , t T
where Y ̲ r , t and Y ¯ r , t denote the deviations that define the lower and upper bounds for a renewable source r at period t, and are a function of the variance. For this paper, Y ̲ r , t = Y ¯ r , t = Var [ P ¯ r , t R ] 1 / 2 . To determine the specific solar availability, auxiliary variables υ ̲ r , t and υ ¯ r , t are used to measure the deviation from the mean value, and a budget of uncertainty, U, limits the total deviation. The uncertainty set (37) can be formulated as follows:
P ¯ r , t R = E [ P ¯ r , t R ] υ ̲ r , t · Y ̲ r , t + υ ¯ r , t · Y ¯ r , t , r R , t T 0 υ ̲ r , t , υ ¯ r , t 1 , r R , t T t T r R ( υ ̲ r , t + υ ¯ r , t ) U
By design, the budget of uncertainty allows the DNO to control the level of conservativeness utilized in the planning process. If the DNO sets U = 0 , the renewable power output parameters are fixed to the mean value at each period; on the other hand, if U = T · | R | , it allows each renewable source to vary during all the time frames during the planning horizon.
The problem of finding the optimal operation of the DN with networked microgrids can be analyzed in two stages: First, the DNO finds the optimal unit commitment decisions for day-ahead planning. After this, the conventional and renewable generators produce the energy under the operational constraints and distribute it via the power grid. Note that the unit commitment occurs before the availability of the intermittent renewable sources is revealed in the first stage, and the power flow and line switching decisions occur in the second stage, once the renewable availability is unfolded. The model stated in (32) can be rewritten to follow the robust optimization approach as follows:
min x c UC + max r U min y , z | x ^ c PG + c μ G + c Sh c DT s . t . Constraints ( 1 ) ( 16 ) , ( 18 ) ( 26 )
where vector r denotes the uncertain parameters of renewable availability, P ¯ r , t R . The optimization model finds the worst-case renewable availability by utilizing the operator max r U . This can be understood as follows: While the DNO tries to minimize the operational cost by controlling the unit commitment ( min x ) and power flow ( min y , z ), another “agent” (i.e., the uncontrollable weather conditions) tries to maximize the cost by adjusting the renewable power availability ( max r U ). The resulting problem is a two-stage MILP problem, in which the first stage has only binary variables, while the second stage has both continuous and binary variables. This structure can be solved utilizing a nested column- and-constraint-generation approach, as demonstrated in [31].

3.3.2. Coordination between DNO and the Microgrids

The DNO and the microgrids try to optimize their operations and there may exist conflicting objectives, considering that the DN business consists of selling electricity to the latter, while the microgrids seek to minimize their operational costs, including the energy purchase. It is assumed that the DNO does not have direct control of the operation of the microgrids nor complete information of the microgrids’ internal systems. Nevertheless, the DNO leverages the DSM policies to influence the behavior of the microgrids, but the microgrids have the last decision on their operation. In this paper, the collaboration process is considered as follows: First, each of the microgrids presents an estimation of their electricity demand to the DNO, without considering any reduction request or power traded to the DN. Second, the DNO uses the consumption estimated by the microgrids to assess its operation and requests load reduction to each of the microgrids, with their respective incentives. Third, the microgrids consider the load reduction requests and the possibility of trading power to the DNO under firm and non-firm contracts, resulting in a new estimation for demand. Finally, the DNO uses the information of the accepted requests and power available from the microgrids to find its optimal operation. This approach can be summarized in the following algorithm:
  • Let the following elements on δ D S M be equal to zero: C m , k , t I n , q ¯ m , k , t , C m , t F , C m , t N , P ¯ m , t F and P ¯ m , t N , and solve each of the microgrids optimization models to obtain D i m , k , t .
  • Let Q ¯ m , k , t = D i m , k , t , and P ¯ m , t F = P ¯ m , t N = 0 and solve the DNO optimization model in (38) to obtain q m , k , t .
  • Let all elements of δ D S M take their original values, and q ¯ m , k , t = q m , k , t , and solve each of the microgrids optimization models to obtain D i m , k , t , p ^ m , t F / N , and q ^ m , k , t .
  • Let Q ¯ m , k , t = 0 , P ¯ m , t F = p ^ m , t F and P ¯ m , t N = p ^ m , t N and solve the DNO optimization model shown in (38) with the updated values obtained from the microgrids, adding the following term to its optimal cost t T k K m M C m , k , t I n · q ^ m , k , t , which represents the cost paid to the microgrids for accepting the load reduction requests.

4. Results

The algorithms and formulations presented in this paper were implemented using the C++ API of CPLEX 20.10 via ILOG Concert Technology. The hardware utilized consists of a DELL machine with Intel(R) Core(TM) i7-10700 CPU 2.90 GHz (eight-core) processor and 32 GB RAM, and it runs on Microsoft Windows 10. The computational time is reported as CPU seconds. For testing the framework discussed above, a benchmark case was carried out and compared with different configurations to demonstrate the flexibility of the proposed approach. The following experiments were carried out: (1) The benchmark case in Section 4.2 was studied under the configuration presented in Section 4.1; (2) In the case presented in Section 4.3, the settings of the conventional generators set, G , were modified to improve the results obtained in the benchmark case; (3) In the case presented in Section 4.4, the size of the renewable energy generators in both the DN and the microgrids was doubled. The second and third experiments are examples of long-term planning decisions.

4.1. Case Settings

A modified version of the IEEE 30-Bus System [47] was utilized as the DN. The operation of the DN was evaluated for a daily horizon divided into | T | = 24 periods, each corresponding to one hour. Within the DN, | M | = 5 microgrids are included in buses 5, 11, 18, 26, and 29. Each of the microgrids have onsite technologies of energy storage and PV farms, as well as an EMS that finds its optimal operation considering both the operational and energy constraints, as described in Section 3.2. The microgrid in bus 5 is a manufacturing facility with multi-process production scheduling, which optimization model is based on [31]; the microgrids in buses 11 and 26 are two interconnected green data centers with job migration capabilities, which optimization model is based on [48]; the microgrid in bus 18 is a vertical farm with controlled environment, whose optimization model is based on [49]; finally, the microgrid in bus 29 is a smart residential building, based on [32]. A total of | G | = 12 conventional generators are located within the DN. The parameters of the generators, adjusted based on the Annual Energy Outlook 2021 published by the EIA [50], are presented in Table 4.
There are | R | = 3 large-scale solar PV facilities of 10, 15, and 20 hectares located in buses 6, 10 and 19, respectively, with a corresponding levelized cost of energy of USD 29.04, 30.43, and 23.92 per MWh [50]. The expected value and variance of solar availability are estimated for the weather conditions of July 15 in Tucson, Arizona, based on the information from the National Renewable Energy Laboratory [51]. A budget of uncertainty of U = 9 is utilized. The expected solar value and the uncertainty for the solar availability for the case of 10 hectares are reported in Table 5, in MW. The availability for the other two solar PV farms are adjusted proportionally to the size of the facilities.
The DNO utilizes | K | = 3 TLOU levels. The rates in USD per MW and levels for TLOU are presented in Table 6 [50].
The demand for the buses that do not contain microgrids is showed in Table 7 and Table 8, where Table 7 reflects the maximum demand during the day for each bus and Table 8 reflects the proportion of that demand during each hour to describe the demand curve. The demand of the nodes with microgrids are not reported in Table 7.
The cost in USD per MW reduced by the microgrids is C m , 1 , t I n = 8 , C m , 2 , t I n = 4 and C m , 3 , t I n = 2 . The price in USD per MW for power purchased from each microgrid under firm contract is C m , t F = 30 and C m , t F = 25 , for each microgrid and during each period. The cost of load shedding in USD per MW shed is C m , 1 , t S = 2000 , C m , 2 , t S = 1000 and C m , 3 , t S = 500 .
The above-described DN with three large-scale solar PV facilities and five microgrids embedded is illustrated in Figure 3.
The information of the microgrids is as follows: (1) the manufacturing microgrid, considered as an industrial microgrid, has five consecutive processes and is planning to produce 30,000 units during the planning horizon, with an onsite solar PV facility of 7.5 hectares; (2) the networks of two green data centers, considered as commercial microgrids located in separate buses, need to process 90 batch-type jobs, and one has PV facilities of 2.5 hectares on each green data center; (3) the vertical farm facility, considered as an industrial microgrid, has a growing area of 1.5 hectares and a solar PV facility with the same extension; (4) the residential building, considered as a residential microgrid, has 20 different appliances with flexible load capabilities and a non-flexible electric demand for each period, and a solar PV facility of 0.5 hectares.

4.2. Benchmark Experiment

The algorithm is tested with the setup mentioned in Section 4.1, and is utilized for a comparison with a sensitivity analysis in the following sections. As shown in Table 9, the microgrids and the DNO collaborated negotiating load reduction requests and electricity trades from the microgrids to the DN. The information of the table can be interpreted as follows: the first and second columns provide general information of the microgrids; the third and fourth columns represent the optimal cost of the microgrids before (Cost 1) and after (Cost 2) the negotiation with the DNO; the fifth and sixth columns present the load reduction requested by the DNO and accepted by the microgrids; the seventh column represents the power that the DNO purchases from the microgrids; finally, the last column represents the green energy coefficient (GEC), which is the proportion of renewable energy utilized in the microgrids. The cost and GECs of microgrids 2 and 4 are the same because they both belong to the same network of green data centers, and both indicators are calculated globally for the network. As can be seen, there is no collaboration between the DNO and the microgrids in this case.
The optimal unit commitment and power generation of the conventional generators (labeled from g1 to g12) are shown in Figure 4. Generators 3, 4, 7 and 10 were not utilized during the planning horizon.
The balance between supply, demand and load shedding is shown in Figure 5. There is a signification amount of load shedding, mainly on buses 19 (4% of total shed) and 30 (96%). The total load shed is 172.15 MW.
The GEC for the DN is 6.79%, while the global GEC is 9.22%, and the total cost for the DNO is USD 25,810.40, indicating that the penalties incurred due to load shedding are greater than the incomes obtained from the electricity sold. The algorithm solved this instance in 55.44 CPU seconds, of which 32.81 CPU seconds corresponds to the two rounds of the DNO problem.

4.3. Change in Location of Generators

The results of the benchmark experiment in Section 4.2 show that there is load shedding on buses 19 and 30. Additionally, this reflects that some generators are not utilized during the planning horizon. Therefore, in this experiment, generator 7 was moved from bus 13 to bus 19, and generator 10 was moved from bus 23 to bus 30. The rest of the parameters remain the same. The results of the microgrid collaboration are shown in Table 10. It can be observed that the only collaboration consists of 0.50 MW of power sold from the vertical farm to the DNO.
The optimal unit commitment and power generation of the conventional generators are shown in Figure 6. Generators 3, 4, 9, and 11 were not utilized during the planning horizon.
The balance between supply, demand and load shedding is shown in Figure 7. There is no load shed in this experiment due to the relocation of generators 7 and 10.
The GEC for the DN is increased to 7.01%, while the global GEC is 9.37%, as a result of relocation and utilization of generators 7 and 10, and the total cost for the DNO is −USD 276,667.00, representing the profits obtained by the energy sold and no penalties associated to load shedding. The algorithm solved this instance in 56.02 CPU seconds, of which 30.86 CPU seconds corresponds to the two rounds of the DNO problem.

4.4. Double Size of Renewable Generators

The results of the experiments in Section 4.2 and Section 4.3 reflect low GEC values. Therefore, in this instance, the size of the renewable generators is doubled. In the DN, this is represented by installing another three large-scale solar PV facilities of 10, 15, and 20 hectares in buses 8, 9, and 14, respectively. The relocation of generators, as described in Section 4.3, is also applied in this experiment. In the microgrids, the size of the solar farms is doubled. The rest of the parameters remain the same as in Section 4.3. The results of the microgrids collaboration are shown in Table 11. It can be observed that the collaboration between the DNO and the microgrids is more significant than in the previous cases due to the expansion of solar capacity of the microgrids, allowing the data centers and the vertical farm sell the surplus energy to the DN.
The optimal unit commitment and power generation of the conventional generators is reflected in Figure 8. Generators 3, 4, 9, and 11 are not utilized during the planning horizon. In general, the amount of power generated is reduced particularly in the periods in which there exists solar availability.
The balance between supply, demand and load shedding is shown in Figure 9. There is no load shed in this experiment and there exists a considerably higher share of renewable energy utilized.
The GEC for the DN increases to 12.88%, while the global GEC is 17.91%, as a result of adding new renewable generators, and the total cost for the DNO is −USD 286,109.00, corresponding to the incomes obtained from the energy sold, and a cost reduction due to the use of cheaper renewable energy. The algorithm solved this instance in 156.70 CPU seconds, of which 138.48 CPU seconds correspond to the two rounds of the DNO problem.

5. Discussion

The framework discussed in this paper presents a collaborative approach between the DN and the microgrids operating within it. The results show that the communication between both entities allows a cost reduction for both systems and a better integration of RES. In the benchmark experiment, the DN requests 38 load reductions to the microgrids, none of which is accepted; additionally, there is not power sold from the microgrids to the DN. The GEC of the DN was 6.79%, and rose up to 9.33% when the RES of the microgrids were considered. The configuration presented in the first case was not capable of satisfying all the demand and some load shedding was carried out, resulting in a total cost of USD 25,810.40, mainly due to shedding penalties. The second case represents the long-term decision of relocating two conventional generators. The total demand is satisfied and no penalty costs are incurred. The profits obtained were USD 276,667.00, and the GEC of the DN increased to 7.01%, and up to 9.37% after considering the microgrids. The third experiment represents the long-term decision of increasing the renewable energy capacities and relocating the generators as was carried out in the second experiment. In this case, the profits increased to USD 286,109.00, a growth of 3.41%, while the GEC of the DN reached 12.88%, and up to 17.91% when the microgrids are considered (an increment of 5.03%).

6. Conclusions

In order to substantially increase the share of renewable energy in the generation mix, large-scale and small-scale solutions should be taken into consideration. The microgrids, small- to medium-scale power systems, provide an opportunity to achieve higher renewable penetration due to the optimized energy management frameworks and the rapid response provided by small-scale generators. The collaboration between microgrids and the distribution network in which they are installed is crucial for the optimal operation of the whole system. Nevertheless, this collaboration is complex due to the conflicting objectives between the distribution network operator and the operators of the microgrids, considering that the former obtains profits for selling electricity to the latter, while the microgrids seek to minimize their operational costs. This paper presents a collaborative approach based on mathematical optimization to coordinate the energy decisions between the distribution network and the multiple microgrids embedded into it. The model considers the energy consumption patterns of the microgrids based on their optimal operation, which is influenced by the demand-side management policies implemented by the distribution network, such as dynamic tariffs, load reduction requests, and energy trades with the main grid.
The proposed model is applied to a case study utilizing the IEEE 30-Bus system with five microgrids embedded into it, including one manufacturing facility, two interconnected data centers located in different buses, one vertical farm, and one smart residential building. Each of these microgrids has onsite renewable generators and energy storage technologies. Both the microgrids and the distribution network obtain their optimal operation by utilizing robust optimization to include the uncertainty related to renewable intermittency. A sensitivity analysis was carried out to demonstrate the flexibility of the algorithm to evaluate the diverse decision that can be made by utilizing it. The numerical experiments indicate that a correct network configuration eliminates the necessity of load shedding, allows a significant utilization of power supplied by the microgrids (22.3 MW), and increases the renewable energy share from 12.88 to 17.91% when microgrids are taken into consideration.
The presented algorithm considers the independence of the microgrids to make their own decisions to optimize their operation under the environment of demand-side management policies proposed by the distribution network operator. The collaboration between the distribution network operator and the microgrids allows finding the optimal operation of the distribution network without direct control nor complete information of the microgrids’ internal systems. Future work may consider other types of negotiation schemes and demand-side management policies, as well as the impact on greenhouse gas emissions.

Author Contributions

Conceptualization, J.L.R.D. and N.F.; methodology, J.L.R.D. and N.F.; software, J.L.R.D. and N.F.; validation, J.L.R.D. and N.F.; formal analysis, J.L.R.D. and N.F.; investigation, J.L.R.D. and N.F.; resources, J.L.R.D. and N.F.; data curation, J.L.R.D. and N.F.; writing—original draft preparation, J.L.R.D. and N.F.; writing—review and editing, J.L.R.D. and N.F.; visualization, J.L.R.D. and N.F.; supervision, J.L.R.D. and N.F.; project administration, J.L.R.D. and N.F.; funding acquisition, not applicable. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All the data utilized in the experiments for this paper is reported in Section 4.

Acknowledgments

J. L. Ruiz Duarte was supported by the Mexican National Council of Science and Technology (CONACYT) and the Mexican Department of Energy (SENER).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ACAlternating current
CO 2 Carbon dioxide
DCDirect current
DNDistribution network
DNODistribution network operator
DRDemand response
DSMDemand-side management
DTDynamic tariff
EIAEnergy Information Administration of the United States
EMSEnergy management system
GECGreen energy coefficient
IEEEInstitute of Electrical and Electronics Engineers
kWhKilowatt-hour
MILPMixed-Integer Linear Programming
MWMegawatt
MWhMegawatt-hour
PVPhotovoltaic
RESRenewable energy sources
TLOUTime-and-Level-of-Use

References

  1. World Bank Open Data. Available online: https://data.worldbank.org/ (accessed on 21 December 2021).
  2. Monthly Energy Review. Available online: https://www.eia.gov/totalenergy/data/monthly/index.php (accessed on 21 December 2021).
  3. Sustainable Development Goals. Available online: https://sdgs.un.org/goals (accessed on 2 August 2019).
  4. SB-350 Clean Energy and Pollution Reduction Act of 2015. Available online: https://www.energy.ca.gov/programs-and-topics/topics/renewable-energy (accessed on 26 July 2019).
  5. Brown, H.E.; Suryanarayanan, S.; Natarajan, S.A.; Rajopadhye, S. Improving reliability of islanded distribution systems with distributed renewable energy resources. IEEE Trans. Smart Grid 2012, 3, 2028–2038. [Google Scholar] [CrossRef]
  6. Mehrasa, M.; Pouresmaeil, E.; Sepehr, A.; Pournazarian, B.; Marzband, M. Control technique for the operation of grid-tied converters with high penetration of renewable energy resources. Electr. Power Syst. Res. 2019, 166, 18–28. [Google Scholar] [CrossRef]
  7. Behrangrad, M. A review of demand side management business models in the electricity market. Renew. Sustain. Energy Rev. 2015, 47, 270–283. [Google Scholar] [CrossRef]
  8. Fattahi Meyabadi, A.; Deihimi, M.H. A review of demand-side management: Reconsidering theoretical framework. Renew. Sustain. Energy Rev. 2017, 80, 367–379. [Google Scholar] [CrossRef]
  9. Kanakadhurga, D.; Prabaharan, N. Demand side management in microgrid: A critical review of key issues and recent trends. Renew. Sustain. Energy Rev. 2022, 156, 111915. [Google Scholar] [CrossRef]
  10. Laitsos, V.M.; Bargiotas, D.; Daskalopulu, A.; Arvanitidis, A.I.; Tsoukalas, L.H. An incentive-based implementation of demand side management in power systems. Energies 2021, 14, 7994. [Google Scholar] [CrossRef]
  11. Anjos, M.F.; Brotcorne, L.; Gómez-Herrera, J.A. Optimal setting of time-and-level-of-use prices for an electricity supplier. Energy 2021, 225, 120517. [Google Scholar] [CrossRef]
  12. Ackermann, T.; Andersson, G.; Söder, L. Distributed generation: A definition. Electr. Power Syst. Res. 2001, 57, 195–204. [Google Scholar] [CrossRef]
  13. Golari, M.; Fan, N.; Wang, J. Two-stage stochastic optimal islanding operations under severe multiple contingencies in power grids. Electr. Power Syst. Res. 2014, 114, 68–77. [Google Scholar] [CrossRef]
  14. Battula, A.R.; Vuddanti, S.; Salkuti, S.R. Review of energy management system approaches in microgrids. Energies 2021, 14, 5459. [Google Scholar] [CrossRef]
  15. Khodaei, A. Resiliency-oriented microgrid optimal scheduling. IEEE Trans. Smart Grid 2014, 5, 1584–1591. [Google Scholar] [CrossRef]
  16. Yoldaş, Y.; Önen, A.; Muyeen, S.M.; Vasilakos, A.V.; Alan, İ. Enhancing smart grid with microgrids: Challenges and opportunities. Renew. Sustain. Energy Rev. 2017, 72, 205–214. [Google Scholar] [CrossRef]
  17. Hytowitz, R.B.; Hedman, K.W. Managing solar uncertainty in microgrid systems with stochastic unit commitment. Electr. Power Syst. Res. 2015, 119, 111–118. [Google Scholar] [CrossRef]
  18. Al Zahr, S.; Doumith, E.A.; Forestier, P. Smart energy: A collaborative demand response solution for smart neighborhood. In Smart Cities: A Data Analytics Perspective; Khan, M.A., Algarni, F., Quasim, M.T., Eds.; Springer: Cham, Switzerland, 2021; pp. 43–62. [Google Scholar]
  19. Guo, Y.; Zhao, C. Islanding-aware robust energy management for microgrids. IEEE Trans. Smart Grid 2018, 9, 1301–1309. [Google Scholar] [CrossRef]
  20. Wang, Y.; Chen, C.; Wang, J.; Baldick, R. Research on resilience of power systems under natural disasters—A review. IEEE Trans. Smart Grid 2016, 31, 1604–1613. [Google Scholar] [CrossRef]
  21. Alharbi, W.; Raahemifar, K. Probabilistic coordination of microgrid energy resources operation considering uncertainties. Electr. Power Syst. Res. 2015, 128, 1–10. [Google Scholar] [CrossRef] [Green Version]
  22. Mirzaei, M.A.; Yazdankhah, A.S.; Mohammadi-Ivatloo, B.; Marzband, M.; Shafie-khah, M.; Catalão, J.P.S. Stochastic network-constrained co-optimization of energy and reserve products in renewable energy integrated power and gas networks with energy storage system. J. Clean. Prod. 2019, 223, 747–758. [Google Scholar] [CrossRef] [Green Version]
  23. Habibian, M.; Zakeri, G.; Downward, A.; Anjos, M.F.; Ferris, M. Co-optimization of demand response and interruptible load reserve offers for a price-making major consumer. Energy Syst. 2020, 11, 45–71. [Google Scholar] [CrossRef] [Green Version]
  24. Mohammad, A.; Zuhaib, M.; Ashraf, I.; Alsultan, M.; Ahmad, S.; Sarwar, A.; Abdollahian, M. Integration of electric vehicles and energy storage system in home energy management system with home to grid capability. Energies 2021, 14, 8557. [Google Scholar] [CrossRef]
  25. Zator, S. Power scheduling scheme for DSM in smart homes with photovoltaic and energy storage. Energies 2021, 14, 8571. [Google Scholar] [CrossRef]
  26. Seshu Kumar, R.; Phani Raghav, L.; Koteswara Raju, D.; Singh, A.R. Impact of multiple demand side management programs on the optimal operation of grid-connected microgrids. Appl. Energy 2021, 301, 117466. [Google Scholar] [CrossRef]
  27. Vergara-Fernandez, L.; Aguayo, M.M.; Moran, L.; Obreque, C. A MILP-based operational decision-making methodology for demand-side management applied to desalinated water supply systems supported by a solar photovoltaic plant: A case study in agricultural industry. J. Clean. Prod. 2022, 334, 130123. [Google Scholar] [CrossRef]
  28. Baldi, S.; Karagevrekis, A.; Michailidis, I.T.; Kosmatopoulos, E.B. Joint energy demand and thermal comfort optimization in photovoltaic-equipped interconnected microgrids. Energy Convers. Manag. 2015, 101, 352–363. [Google Scholar] [CrossRef] [Green Version]
  29. Korkas, C.D.; Baldi, S.; Michailidis, I.; Kosmatopoulos, E.B. Occupancy-based demand response and thermal comfort optimization in microgrids with renewable energy sources and energy storage. Appl. Energy 2016, 163, 93–104. [Google Scholar] [CrossRef]
  30. Gomez-Herrera, J.A.; Anjos, M.F. Optimal collaborative demand-response planner for smart residential buildings. Energy 2018, 161, 370–380. [Google Scholar] [CrossRef] [Green Version]
  31. Ruiz Duarte, J.L.; Fan, N.; Jin, T. Multi-process production planning with variable renewable integration and demand response. Eur. J. Oper. Res. 2020, 281, 186–200. [Google Scholar] [CrossRef]
  32. Wang, S.; Gangammanavar, H.; Ekşioğlu, S.D.; Mason, S.J. Stochastic optimization for energy management in power systems with multiple microgrids. IEEE Trans. Smart Grid 2019, 10, 1068–1079. [Google Scholar] [CrossRef]
  33. Fang, X.; Yang, Q.; Wang, J.; Yan, W. Coordinated dispatch in multiple cooperative autonomous islanded microgrids. Appl. Energy 2016, 162, 40–48. [Google Scholar] [CrossRef]
  34. Fang, X.; Ma, S.; Yang, Q.; Zhang, J. Cooperative energy dispatch for multiple autonomous microgrids with distributed renewable sources and storages. Energy 2016, 99, 48–57. [Google Scholar] [CrossRef]
  35. Wang, Z.; Chen, B.; Wang, J.; Begovic, M.M.; Chen, C. Coordinated energy management of networked microgrids in distribution systems. IEEE Trans. Smart Grid 2015, 6, 45–53. [Google Scholar] [CrossRef]
  36. Wang, Z.; Chen, B.; Wang, J.; Kim, J. Decentralized energy management system for networked microgrids in grid-connected and islanded modes. IEEE Trans. Smart Grid 2016, 7, 1097–1105. [Google Scholar] [CrossRef]
  37. Wang, Z.; Chen, B.; Wang, J.; Chen, C. Networked microgrids for self-healing power systems. IEEE Trans. Smart Grid 2016, 7, 310–319. [Google Scholar] [CrossRef]
  38. Kou, P.; Liang, D.; Gao, L. Distributed empc of multiple microgrids for coordinated stochastic energy management. Appl. Energy 2017, 185, 939–952. [Google Scholar] [CrossRef]
  39. Lv, T.; Ai, Q. Interactive energy management of networked microgrids-based active distribution system considering large-scale integration of renewable energy resources. Appl. Energy 2016, 163, 408–422. [Google Scholar] [CrossRef]
  40. Haddadian, H.; Noroozian, R. Multi-microgrids approach for design and operation of future distribution networks based on novel technical indices. Appl. Energy 2017, 185, 650–663. [Google Scholar] [CrossRef]
  41. Wang, Z.; Wang, J. Service restoration based on ami and networked mgs under extreme weather events. IET Gener. Transm. Distrib. 2017, 11, 401–408. [Google Scholar] [CrossRef]
  42. Ma, W.; Wang, J.; Gupta, V.; Chen, C. Distributed energy management for networked microgrids using online admm with regret. IEEE Trans. Smart Grid 2018, 9, 847–856. [Google Scholar] [CrossRef]
  43. Ferro, G.; Minciardi, R.; Parodi, L.; Robba, M.; Rossi, M. Optimal control of multiple microgrids and buildings by an aggregator. Energies 2020, 13, 1058. [Google Scholar] [CrossRef] [Green Version]
  44. Ruiz Duarte, J.L.; Fan, N. Operations of a microgrid with renewable energy integration and line switching. Energy Syst. 2019, 10, 247–272. [Google Scholar] [CrossRef] [Green Version]
  45. Minchala-Avila, L.I.; Garza-Castañón, L.E.; Vargas-Martínez, A.; Zhang, Y. A review of optimal control techniques applied to the energy management and control of microgrids. Procedia Comput. Sci. 2015, 52, 780–787. [Google Scholar] [CrossRef] [Green Version]
  46. Glossary. Available online: https://www.eia.gov/tools/glossary/ (accessed on 6 April 2020).
  47. Power Systems Test Case Archive. Available online: https://labs.ece.uw.edu/pstca/ (accessed on 5 April 2020).
  48. Ruiz Duarte, J.L.; Fan, N. Operations of data centers with onsite renewables considering greenhouse gas emissions. Submitted.
  49. Ruiz Duarte, J.L.; Fan, N.; Zhang, Y.; Kacira, M. Flexible operations of vertical farms with solar renewable integration. Submitted.
  50. Annual Energy Outlook 2021. Available online: https://www.eia.gov/outlooks/aeo/ (accessed on 17 February 2022).
  51. National Solar Radiation Database. Available online: https://nsrdb.nrel.gov/ (accessed on 14 September 2019).
Figure 1. DN with multiple microgrids.
Figure 1. DN with multiple microgrids.
Energies 15 02350 g001
Figure 2. The EMS of microgrids.
Figure 2. The EMS of microgrids.
Energies 15 02350 g002
Figure 3. IEEE 30-Bus system as the DN with large-scale solar PV farms and five embedded microgrids.
Figure 3. IEEE 30-Bus system as the DN with large-scale solar PV farms and five embedded microgrids.
Energies 15 02350 g003
Figure 4. Unit commitment and generation of conventional generators for the benchmark experiment.
Figure 4. Unit commitment and generation of conventional generators for the benchmark experiment.
Energies 15 02350 g004
Figure 5. Balance between electricity supply and demand for the benchmark experiment.
Figure 5. Balance between electricity supply and demand for the benchmark experiment.
Energies 15 02350 g005
Figure 6. Unit commitment and generation of conventional generators when generators are relocated.
Figure 6. Unit commitment and generation of conventional generators when generators are relocated.
Energies 15 02350 g006
Figure 7. Balance between electricity supply and demand when generators are relocated.
Figure 7. Balance between electricity supply and demand when generators are relocated.
Energies 15 02350 g007
Figure 8. Unit commitment and generation of conventional generators when solar capacity is doubled.
Figure 8. Unit commitment and generation of conventional generators when solar capacity is doubled.
Energies 15 02350 g008
Figure 9. Balance between electricity supply and demand when solar capacity is doubled.
Figure 9. Balance between electricity supply and demand when solar capacity is doubled.
Energies 15 02350 g009
Table 1. Sets and indices.
Table 1. Sets and indices.
SymbolDescription
T Set of time frames in the planning horizon, indexed by t { 1 , 2 , , T } , T = | T |
I Set of buses, indexed by i
M I Set of microgrid buses, indexed by m (connected to bus i m I )
E Set of transmission lines, indexed by e (denoted by e = ( i e , j e ) , with tail bus i e and head bus j e )
G Set of conventional energy generators, indexed by g (connected to bus i g I )
R Set of renewable energy generators, indexed by r (connected to bus i r I )
K Set of consumption levels for dynamic tariff policy, indexed by k
Table 2. Model parameters.
Table 2. Model parameters.
SymbolDescription
P ̲ g G / P ¯ g G Lower/upper limits of power output for generation unit g
T g d 0 / u 0 Minimum number of periods generation unit g must be initially offline/online
T g d / u Minimum number of periods generation unit g must remain offline/online once it is shut down/started up
R ˜ g d / u Maximum shutdown/startup ramping rate for generation unit g G when is shut down/started up
R g d / u Maximum ramp-down/up rate for generation unit g between adjacent periods
C g T d / T u Fixed shutdown/startup cost for generation unit g
C g R d / R u Ramp-down/up unitary cost for generation unit g
C g P Generation cost function for generation unit g
F e Capacity of transmission line e
B e Susceptance of transmission line e
P ¯ m , t F / N Limit under Firm (F)/Non-firm (N) contract for power purchased from microgrid m, at period t
C m , t F / N Cost of power purchased from microgrid m, under Firm (F)/Non-firm (N) contract, at period t
D i , k , t Energy demand at bus i, at consumption level k, at period t
C m , k , t I n Load reduction incentive per unit offered to microgrid m, at consumption level k, at period t
C r R Production cost of renewable type r
C i , k , t D T Energy price for bus i at consumption level k, at period t
C i , k , t S Cost of load shedding per unit on bus i, at consumption level k, at time t
P ¯ r , t R Availability of renewable power of RES type r at period t
P ¯ i , k , t D T Amount of conventional energy available at bus i, at consumption level k, at period t
Q ¯ m , k , t Limit for load reduction requests for microgrid m, at consumption level k, at period t
Table 3. Decision variables.
Table 3. Decision variables.
SymbolDescription
x g , t { 0 , 1 } Binary variable indicating if generation unit g is committed ( x g , t = 1 ) or not ( x g , t = 0 ) at period t
y g , t u / d { 0 , 1 } Binary variable indicating if generation unit g is started up/shut down ( y g , t u / d = 1 ) or not ( y g , t u / d = 0 ) at period t
p g , t G + / G Electricity generation/curtailment of conventional generator g at period t
p g , t G Power supplied from generator g to the DN at period t
p g , t u / d Ramp up/down of conventional generator g at period t
p r , t R Electricity generation of renewable generator r at period t
f e , t Electricity flow on line e, at period t
z e , t { 0 , 1 } Binary variable indicating if line e is switched on ( z e , t = 1 ) or off ( z e , t = 0 ), at period t
θ i e , t Voltage angle of bus i on line e, at period t
p m , t F / N Electricity purchased from microgrid m, under Firm (F)/Non-firm (N) contract, at period t
p i , k , t ( ) Electricity sold to bus i, at consumption level k at period t
q m , k , t Amount of load reduction requested to microgrid m, at consumption level k, at period t
s i , k , t Amount of load shed from bus i, at consumption level k, at period t
Table 4. Parameters of conventional generators.
Table 4. Parameters of conventional generators.
i g P ̲ g G P ¯ g G T g d 0 T g u 0 T g d T g u R ˜ g d R ˜ g u R g d R g u C g Td C g Tu C g Rd C g Ru C g P
11280002440404040012510.7810.78107.83
21280001140404040012510.6610.66106.62
3128000114040404002499.759.7597.50
4128000484040404002499.739.7397.32
71280002440404040003.713.7137.11
127.550001125252525003.323.3233.21
137.550001125252525012511.2011.20111.97
154.530001115151515004.124.1241.77
228.2555004827.527.527.527.5012510.7810.78107.83
234.53000241515151502499.759.7597.50
2764000112020202002499.739.7397.32
28640004820202020003.323.3233.21
Table 5. Parameters of the solar photovoltaic facility of 10 hectares in MW.
Table 5. Parameters of the solar photovoltaic facility of 10 hectares in MW.
t1…78910111213141516171819…24
E [ P ¯ 1 , t R ] 01.244.909.0612.8416.0717.3116.4714.6411.637.782.240
Var [ P ¯ 1 , t R ] 1 / 2 00.481.141.051.252.569.762.631.711.431.441.870
Table 6. Parameters of TLOU.
Table 6. Parameters of TLOU.
LevelResidentialIndustrialCommercial
MWPeakOff-PeakMWPeakOff-PeakMWPeakOff-Peak
k = 1 [0,5)171.86122.05[0,15)99.5070.66[0,10)150.53106.90
k = 2 [5,10)183.01133.20[15,50)105.9577.12[10,30)160.30116.67
k = 3 [10,)189.83140.02[50,)109.9081.06[30,)166.26122.64
Table 7. Maximum demand for each bus.
Table 7. Maximum demand for each bus.
Bus123467891012131415
Demand021.72.47.6022.83005.811.206.28.2
Bus161719202122232425272830-
Demand3.599.52.217.503.28.700010.6-
Table 8. Hourly demand for each bus.
Table 8. Hourly demand for each bus.
Hour123456789101112
Demand0.750.730.690.660.650.650.680.740.830.890.920.94
Hour131415161718192021222324
Demand0.910.900.900.860.850.880.921.000.970.950.900.85
Table 9. Collaboration of DN and microgrids in the benchmark experiment.
Table 9. Collaboration of DN and microgrids in the benchmark experiment.
MicrogridTypeCost 1 (USD)Cost 2 (USD)L. R. ReqL. R. AccTrades (MW)GEC (%)
1Manufacturing7,480,520.007,480,520.00220021.28
2Data center24,952.8024,952.80130033.63
3Vertical farm1363.271363.2700042.73
4Data center24,952.8024,952.8030033.63
5Residential8933.858933.8500021.84
Table 10. Collaboration of DN and microgrids when generators are relocated.
Table 10. Collaboration of DN and microgrids when generators are relocated.
MicrogridTypeCost 1 (USD)Cost 2 (USD)L. R. ReqL. R. AccTrades (MW)GEC (%)
1Manufacturing7,480,510.007,480,510.00220021.28
2Data center24,952.3024,952.3010033.59
3Vertical farm1366.281363.27100.5045.18
4Data center24,952.3024,952.3020033.59
5Residential8933.528933.5220021.84
Table 11. Collaboration of DN and microgrids when solar capacity is doubled.
Table 11. Collaboration of DN and microgrids when solar capacity is doubled.
MicrogridTypeCost 1 (USD)Cost 2 (USD)L. R. ReqL. R. AccTrades (MW)GEC (%)
1Manufacturing7,427,570.007,427,570.00220028.23
2Data center21,211.8018,968.5107.3049.00
3Vertical farm1116.811076.14106.8967.87
4Data center21,281.2021,211.80108.1149.00
5Residential8370.918370.9100027.12
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ruiz Duarte, J.L.; Fan, N. Operation of a Power Grid with Embedded Networked Microgrids and Onsite Renewable Technologies. Energies 2022, 15, 2350. https://doi.org/10.3390/en15072350

AMA Style

Ruiz Duarte JL, Fan N. Operation of a Power Grid with Embedded Networked Microgrids and Onsite Renewable Technologies. Energies. 2022; 15(7):2350. https://doi.org/10.3390/en15072350

Chicago/Turabian Style

Ruiz Duarte, José Luis, and Neng Fan. 2022. "Operation of a Power Grid with Embedded Networked Microgrids and Onsite Renewable Technologies" Energies 15, no. 7: 2350. https://doi.org/10.3390/en15072350

APA Style

Ruiz Duarte, J. L., & Fan, N. (2022). Operation of a Power Grid with Embedded Networked Microgrids and Onsite Renewable Technologies. Energies, 15(7), 2350. https://doi.org/10.3390/en15072350

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop