Optimal Load Shedding for Maximizing Satisfaction in an Islanded Microgrid

A microgrid (MG) is a discrete energy system that can operate either in parallel with or independently from a main power grid. It is designed to enhance reliability, carbon emission reduction, diversification of energy sources, and cost reduction. When a power fault occurs in a grid, an MG operates in an islanded manner from the grid and protects its power generations and loads from disturbance by means of intelligent load shedding. A load shedding is a control procedure that results in autonomous decrease of the power demands of loads in an MG. In this study, we propose a load shedding algorithm for the optimization problem to maximize the satisfaction of system components. The proposed algorithm preferentially assigns the power to the subdemand with a high preference to maximize the satisfaction of power consumers. In addition, the algorithm assigns the power to maximize the power sale and minimize the power surplus for satisfaction of power suppliers. To verify the performance of our algorithm, we implement a multi-agent system (MAS) on top of a conventional development framework and assess the algorithm’s adaptability, satisfaction metric, and running time.


Introduction
Recently, microgeneration has been considered for use in managing the growth of electric power demands and CO 2 emission.The Consortium for Electric Reliability Technology Solutions (CERTS) has initiated research on the impact of connecting small generation units to low-voltage electricity distribution systems in order to enhance the reliability of a utility grid system.Accordingly, a microgrid (MG) has been introduced [1].An MG is a small-scale autonomous power system that operates either in parallel to or islanded from utility main grids.An MG allows for fast installation of electricity supply without the need for expensive transmission infrastructure investments.It is composed of on-site distributed generators (DGs), energy storages (ESs), and loads [2].DGs are power generators such as renewable energy sources (e.g., photovoltaic (PV), wind, or solar) and conventional generators such as microturbines or diesel generators.ESs include lithium-ion batteries.They can be either power generators that produce power or power consumers that expend power.Loads such as residential or commercial buildings are power consumers.
An MG has two operational modes: grid-connected and islanded.A grid-connected operation indicates a situation in which an MG is connected to the main grid.An islanded operation indicates a situation in which an MG is disconnected from the main grid when a fault occurs in the grid.The amount of produced power of DGs depends on the demands of loads.The difference between the generated power and load demand is controlled by the well-functioning main grid in the Energies 2017, 10, 45 2 of 13 grid-connected MG.On an islanded MG, because the main energy sources are renewable energy sources that have an uncertain nature, the amount of produced power may not controlled.Thus, the control is commonly accomplished by a load shedding system, not a mechanical controller.
Load shedding is a process in which some of the power demands are shed to protect MGs from potential dangers [3].Two architectures for load shedding, centralized and decentralized, have been investigated [4].In a centralized architecture, system components are controlled by a central control system, which can lead to optimal utilization of MG resources.In a decentralized architecture, individual components cooperatively or competitively share the power supply or power demand based on their individual capabilities.A decentralized architecture can easily integrate new components.However, real-time optimal utility under various operating conditions is not possible.Although considerable research on load shedding systems such as islanding detection and autonomous operation/management have been extensively studied [5], research on load shedding algorithms remain an open issue for further research.A load shedding algorithm determines the amount of power demand that can be reduced.In general, a load shedding algorithm finds the most stable equilibrium operating point for the system with the fewest shed demands.In other words, a load shedding algorithm is designed to make full use of limited generation resources and maximize the beneficial performance [6].In this study, we propose a load shedding algorithm for the optimization problem in a centralized architecture of an islanded MG.We use the satisfaction of an MG system player as performance criteria.
In this paper, the major contributions are: -Identifying the satisfaction of system components as performance criteria; -Defining the optimization problem to maximize the satisfaction; -Developing a load shedding algorithm for the optimization problem; -Implementing a multi-agent system (MAS) for the load shedding algorithm; -Performing experiments to verify performance of the implemented system.
The remainder of the paper is organized as follows.In Section 2, we describe related studies and system models.In Section 3, we propose a load shedding algorithm for the optimization problem to maximize player satisfaction.In Section 4, experimental results are discussed to verify the performance of our algorithm based on an implemented system operating in a conventional development framework.Finally, a conclusion and future research directions are provided in Section 5.

Related Work
Recent studies have explored the load shedding problem in an islanded MG.In Reference [2], fine-grained load shedding policies for high-stressed grids were proposed.High-stressed grids are power grids with a very large and nearly continuous supply-demand gap.Two algorithms have been proposed.In a distributed algorithm, the power throttling level is stochastically computed and the demand is shed based on the level.In a centralized algorithm, feeder-level groups are formed among the loads.In a group, the demand for the highest consumption is shed to the specific level of power consumption.Other demands in the group are then shed alternatively.After the cycling, if a supply demand gap remains, the next group is selected and the process continues until the gap is resolved.
In References [7][8][9], researchers solved the load shedding problem by using a game theoretic approach.They considered MG as a collection of individual components that compete for power.The three methodologies frame a power system in the context of a multiplayer nonzero-sum game, the bankruptcy problem, or the bargaining solution, respectively.
In References [10][11][12], a heuristic method was adopted to solve the load shedding problem.In References [10,11], load shedding strategies were proposed in which a V-shaped transfer function for binary particle swarm optimization or a binary firefly algorithm was used.The objective function Energies 2017, 10, 45 3 of 13 is to meet the maximum demand based on the priority.In Reference [12], a genetic algorithm was applied to minimize the shed demands at all stages and maximize the lowest swing frequency.Penalty functions and chromosomes with varying lengths were utilized in the algorithm to determine the optimal shed demands.
In References [13,14], an optimization solution for load shedding was used.In Reference [13], specifically, a control strategy for stratified optimization was presented.This strategy stratifies demands into different levels based on their importance and the load frequency regulation effect coefficient.Based on the levels, the optimization solution is provided to minimize the cost of load shedding and voltage fluctuation.In Reference [14], it was observed that in addition to frequency, the rate of change of frequency (RoCoF) was an important index to determine the demand reduction.Based on RoCoF, a load shedding sequence was determined in real time by evaluating operating profits according to the willingness to pay (WTP) of the loads.

System Model
In a conventional under-frequency load shedding scheme, when the frequency of MG falls below permissible levels, some power demands in the system are shed in predetermined steps [15].In the first step, an anticipated overload (L) is given by: L = Total load demands − Total power supplied Total power supplied .
In the second step, the total amount of demands to be shed is calculated to maintain the frequency above permissible levels for the maximum L below that of [16]: where LD is the total number of demands to be shed, d is the load reduction factor, f is the minimum permissible frequency, and f n is the nominal system frequency.The load reduction factor indicates the percentage change in load demand per percentage change in frequency.In the final step, the number of demands to be shed for each load must be determined.In this study, we focus on the final step of a load shedding scheme.
We next describe a system model for load shedding in an islanded MG.The system is defined as a pair (L, P), where P is a vector of the power supplied by N DGs and L is a vector of power demanded by M loads (i.e., P = (P 1 , P 2 , . . . ,P N ) and L = (L 1 , L 2 , . . . ,L M )).There are some stochastic variables in DGs like power generation of renewable energy sources (e.g., wind turbines (WT) and PV cells) [17].The stochastic models which consider uncertain nature of some variables such as generation of WT and PV are proposed [18][19][20].In the models, the amount of wind speed or solar irradiance is forecasted by a time series method.For the forecasted amount, a probability distribution function is assigned to obtain the amount of stochastic variable.In our system model, P includes the power supplied by dispatchable energy sources, WT, and PV cells (i.e., P n ∈ {P dis , P WT , P PV }).
The demand of the load m is divided into two types: critical load demand and controllable load demand.Critical load includes hospitals, water stations, and other infrastructures to be associated with basic needs of human life.When MG is disconnected from the main grid, controllable load demand is shed gradually for the continuous power supply to critical load.The power supplied by DGs is firstly assigned to the critical load demands, L m .A controllable load demand is composed of ).Let l k m denote the k-th subdemand of the load m.For each subdemand, a preference is assigned in the range of 1 to Pre max .The power is assigned gradually to the subdemand with high preference.When ∑ ∑ M m=1 L m > ∑ N n=1 P n , our load shedding algorithm is activated.A summary of important symbols used in this study is given in Table 1.

Proposed Load-Shedding Algorithm
In this section, we devise an optimization scheme for load shedding to maximize the satisfaction of the system player.The system player is divided into power suppliers such as DGs and ESs and the power consumers such as loads and ESs.In this study, we use DGs as power suppliers and loads as power consumers.The system optimum is achieved when the satisfaction of both the DG and load are maximized.Given x m and y n , which are the satisfaction indices of the load m and DG n in the load shedding, U Load (x m ) and U DG (y n ) indicate utility functions of the load and the DG at a specific time interval.The time-dependent utility models general situations in which the load and DG may have various load demands and may supply power at different time intervals.We assume that U Load (x m ) and U DG (y n ) are continuously differentiable, strictly concave, and increase over the range x m ≥ 0 and y n ≥ 0 [21].In this model, our optimization problem is given by: Our optimization problem is divided into two parts: that of the load and that of the DG.First, we define the utility function from the load point of view.The utility function of the load reflects the satisfaction of power allocation and the dissatisfaction of power shortage.We analyze the relationship between the preference of a subdemand and the satisfaction.The satisfaction is achieved when the power is assigned to the subdemand with high preference greater than that when the power is assigned to the subdemand with a relatively low preference.Thus, we assign the power to subdemands in the order of high preference to low preference.To obtain the numerical value of the load satisfaction, we calculate k , which is given by: subject to: where K min indicates min 1≤m≤M K m .The satisfaction index of the subdemand l k m in the load m is given by: Energies 2017, 10, 45 where w l k m is the preference weight of the subdemand l k m .It is given by: where Pre l k m is the preference of the subdemand l k m .The satisfaction index of the load m, x m is given by: subject to: where ϕ m is a binary variable taken to be 0 or 1 for the linear programming problem.Thus, the satisfaction value of the load m is given by: subject to: where The variable ¬ϕ m indicates the negation of ϕ m .In other words, ¬ϕ m = (ϕ m → 0 = 1).The variable α L m denotes the sensitivity to the preference of L m .In Equation (10), the first term indicates the satisfaction of power allocation for l 1 m ∼ l k m (or l k +1 m ) and the second term indicates the dissatisfaction for l k +1 m (or l k +2 m ) ∼ l K m m resulting from the power shortage.Thus, the utility function of the load is defined as follows: In the second part of our optimization problem, we define the utility function from the DG point of view.The utility function of the DG reflects the satisfaction of the power sale and the dissatisfaction of the power surplus.The uncertainty of power generation is one of the important features when a load selects its power supplier.It is because the amount of produced power by DGs is not estimated exactly in an islanded MG with renewable energy sources.For dealing with the uncertainty, we assign the uncertainty value to each DG according to the deviations of forecasted amounts of produced power in long-term historical data [17].Using k in Equation ( 4), the subdemands of M loads, ∑ k k=1 ∑ M m=1 l k m , are assigned to DGs in the order of low to high uncertainty.The satisfaction index of the DG n, y n is given by: where: and subject to: Energies 2017, 10, 45 6 of 13 Thus, the satisfaction value of the DG n is given by where α P n is the sensitivity to the surplus power of P n .The utility function of DG is defined as follows: Therefore, our optimization problem is described as follows: maximize

Performance Evaluation
An MAS is presented as an effective solution for a centralized or decentralized MG system [22].Thus, we implemented an MAS for the load shedding management of an islanded MG using a Java Agent Development Environment (JADE) framework (Telecom Italia Lab, Torino, Italy) [23].The MAS operating on JADE emulates distributed load shedding operations.To optimize load shedding, the MAS was linked with the IBM CPLEX Optimizer (IBM, Armonk, NY, USA) [24], which is a mathematical programming solver for linear programming.We configured our system with MG central control (MGCC), using DGs as power generators and loads as power consumers.The components were classified according to operational capabilities and assigned appropriate control agents: Ag MGCC , Ag DG , and Ag load [25].One major functionality of MAS is the communication among the distributed agents.By using the Foundation for Intelligent Physical Agents (FIPA) guidelines for agent-based communication, FIPA-Contract Net Protocol (CNP) performatives are standardized [26].The performatives included CFP (Call For Proposal), PROPOSE, REFUSE, ACCEPT, REJECT, FAILURE, and INFORM.Figure 1 shows the message flows when using the FIPA-CNP performatives.The MG operation follows an interval-based planning paradigm.First, the Ag MGCC sends a CFP message to Ag DG s and Ag load s to produce a load shedding plan of the next interval (e.g., interval t).The Ag DG s and Ag load s respond with PROPOSE messages to send their information to the Ag MGCC .Such information includes the amount of power generated or the number of power demands.Each Ag load additionally sends a list of subdemands with their preferences.If the supplied power is less than the demanded power, the Ag MGCC initiates the load shedding scheme to determine the number of shed demands of the loads.Once the number of shed demands is determined, the Ag MGCC sends ACCEPT messages to notify the Ag DG s and Ag load s of the results of the load shedding.For the Ag DG s, the amount of power to be supplied to the loads is given.For the Ag load s, the number of shed demands and the amount of power to be allocated to the loads are given.The notified load shedding plan is implemented at time interval t.The load shedding planning and implementation procedures repeat at each interval.
To show the functionality and feasibility of our algorithm, we configured the system with one Ag MGCC , five Ag DG s, and seven Ag load s.The modified CERTS MG consists of five generators, seven loads, and 23 breakers [27].A graphical representation of normal operation of the modified CERTS is shown in Figure 2.
demands is determined, the sends ACCEPT messages to notify the s and s of the results of the load shedding.For the s, the amount of power to be supplied to the loads is given.For the s, the number of shed demands and the amount of power to be allocated to the loads are given.The notified load shedding plan is implemented at time interval .The load shedding planning and implementation procedures repeat at each interval.To show the functionality and feasibility of our algorithm, we configured the system with one , five s, and seven s.The modified CERTS MG consists of five generators, seven loads, and 23 breakers [27].A graphical representation of normal operation of the modified CERTS is shown in Figure 2. Figure 3 shows the system in a distributed environment based on three PCs.The left subfigure shows the test environment of our laboratory.One is on the first PC, five s are on the second PC, and seven s are on the third PC.The Internet is used for inter-agent communications.The right subfigure shows the operational process on the implemented JADE system.In our system, for the sake of experiments, we assume that the critical load demands ∑ = 0 and ∑ = ∑ .The total controllable load demands ∑ is set to 360 kWh and the total supply ∑ is set to 50%, 65%, or 80% of the total load demands.For α and α , we set both of them to 1.The uncertainty values of DGs are set to the same value.The details of load data is given in Table 2 and the subdemands ( ) for each load are sorted by preference ( ) in ascending order from the value '1'.We set to 10.We divided the time in a day to 24 intervals for MG operation, and the length of an interval was one hour.The amount of demand of an for each interval follows a Poisson distribution.Figure 3 shows the system in a distributed environment based on three PCs.The left subfigure shows the test environment of our laboratory.One Ag MGCC is on the first PC, five Ag DG s are on the second PC, and seven Ag load s are on the third PC.The Internet is used for inter-agent communications.The right subfigure shows the operational process on the implemented JADE system.In our system, for the sake of experiments, we assume that the critical load demands ∑ 7 m=1 L m = 0 and ∑ 5 n=1 P n = ∑ 5 n=1 P n .The total controllable load demands ∑ 7 m=1 L m is set to 360 kWh and the total supply ∑ 5 n=1 P n is set to 50%, 65%, or 80% of the total load demands.For α L m and α P n , we set both of them to 1.The uncertainty values of DGs are set to the same value.The details of load data is given in Table 2 and the subdemands (l k m ) for each load are sorted by preference (Pre l k m ) in ascending order from the value '1'.We set Pre max to 10.We divided the time in a day to 24 intervals for MG operation, and the length of an interval was one hour.The amount of demand of an Ag load for each interval follows a Poisson distribution.
Energies 2017, 10, 45 8 of 13 360 kWh and the total supply ∑ is set to 50%, 65%, or 80% of the total load demands.For α and α , we set both of them to 1.The uncertainty values of DGs are set to the same value.The details of load data is given in Table 2 and the subdemands ( ) for each load are sorted by preference ( ) in ascending order from the value '1'.We set to 10.We divided the time in a day to 24 intervals for MG operation, and the length of an interval was one hour.The amount of demand of an for each interval follows a Poisson distribution.Figure 4 shows the functionality of our load shedding algorithm in three test scenarios.In experiments, the total power supply was set to 180 kWh for 50% of the total load demands (50% test scenario), 234 kWh for 65% (65% test scenario), and 288 kWh for 80% (80% test scenario).The subfigures show the number of shed demands of 3 for 24 time intervals.In the figures, the height of a bar indicates the number of power demands, and the violet color of a bar indicates the amount of allocated power to the load.The blue color of a bar indicates the number of shed demands of the load.At each interval, as the power demands of the load changes, the number of shed demands accordingly changes.In addition, as the power supply increases, the number of shed demands decreases.In the experiments, the power was assigned about 42%, 52%, or 77% of the demands to 3 for 50%, 65% or 80% test scenarios.The plug-and-play functionality is one key feature of multi-agent based optimization systems.Thus, to investigate the plug-and-play functionality of our system, we configure that DG (50 kWh) among five DGs is removed and added at intervals 11 and 15. Figure 5 shows that our system is feasible for plug-and-play operation when ∑ is changed from 288 to 238 kWh.  Figure 4 shows the functionality of our load shedding algorithm in three test scenarios.In experiments, the total power supply was set to 180 kWh for 50% of the total load demands (50% test scenario), 234 kWh for 65% (65% test scenario), and 288 kWh for 80% (80% test scenario).The subfigures show the number of shed demands of 3 Ag load for 24 time intervals.In the figures, the height of a bar indicates the number of power demands, and the violet color of a bar indicates the amount of allocated power to the load.The blue color of a bar indicates the number of shed demands of the load.At each interval, as the power demands of the load changes, the number of shed demands accordingly changes.In addition, as the power supply increases, the number of shed demands decreases.In the Energies 2017, 10, 45 9 of 13 experiments, the power was assigned about 42%, 52%, or 77% of the demands to 3 Ag load for 50%, 65% or 80% test scenarios.
The plug-and-play functionality is one key feature of multi-agent based optimization systems.Thus, to investigate the plug-and-play functionality of our system, we configure that DG P 4 (50 kWh) among five DGs is removed and added at intervals 11 and 15. Figure 5 shows that our system is feasible for plug-and-play operation when ∑ 5 n=1 P n is changed from 288 to 238 kWh.To verify the performance of our algorithm, we compared it with the max-min fairness algorithm.In the max-min fairness, a resource is allocated to customers in order of increasing demand and customers with unsatisfied demands split the remaining resources.For the performance comparison, we designed the satisfaction metric (SM) of as follows: where * indicates the number of subdemands that the power is allocated.The SM of is given by:   To verify the performance of our algorithm, we compared it with the max-min fairness algorithm.In the max-min fairness, a resource is allocated to customers in order of increasing demand and customers with unsatisfied demands split the remaining resources.For the performance comparison, we designed the satisfaction metric (SM) of as follows: where * indicates the number of subdemands that the power is allocated.The SM of is given by: To verify the performance of our algorithm, we compared it with the max-min fairness algorithm.In the max-min fairness, a resource is allocated to customers in order of increasing demand and Energies 2017, 10, 45 10 of 13 customers with unsatisfied demands split the remaining resources.For the performance comparison, we designed the satisfaction metric (SM) of Ag load m as follows: where k * indicates the number of subdemands that the power is allocated.The SM of Ag DG n is given by: where P n indicates the amount of power to be allocated to the loads.Figure 6 shows the SM for different load conditions.Each value in the charts indicates ∑ M m=1 SM m or ∑ N n=1 SM n at each interval.Figure 6a shows the performance from the load points of view.On average, the SM of our algorithm was about 39% better than that of max-min algorithm.In detail, our algorithm showed about 55.1%, 39.4%, and 22.5% better performance for 50%, 65% and 80% test scenarios.Figure 6b shows the SM from the DG points of view.On average, the SM of our algorithm was about 30% better than that of the max-min algorithm.In detail, our algorithm showed about 44.8%, 29.3%, and 15.8% better performance for the three test scenarios.As a result, our algorithm improves the degree of satisfaction more than does the max-min fair algorithm from both the DG and load points of view.
Energies 2017, 10, 45 10 of 13 = , where indicates the amount of power to be allocated to the loads.Figure 6 shows the SM for different load conditions.Each value in the charts indicates ∑ or ∑ at each interval.Figure 6a shows the performance from the load points of view.On average, the SM of our algorithm was about 39% better than that of max-min algorithm.In detail, our algorithm showed about 55.1%, 39.4%, and 22.5% better performance for 50%, 65% and 80% test scenarios.Figure 6b shows the SM from the DG points of view.On average, the SM of our algorithm was about 30% better than that of the max-min algorithm.In detail, our algorithm showed about 44.8%, 29.3%, and 15.8% better performance for the three test scenarios.As a result, our algorithm improves the degree of satisfaction more than does the max-min fair algorithm from both the DG and load points of view.Figure 7 shows benefits for loads and DGs for different load conditions.On the load side, because the load demand is shed by force, we consider the benefits in terms of the satisfaction when the subdemands with high preferences are met.Figure 7a shows the benefit of load satisfaction when the SM of our algorithm is compared with the SM of the max-min algorithm.In detail, our algorithm shows about 153%, 132%, and 116% better benefits for 50%, 65%, and 80% test scenarios.Figure 7b shows the benefit of DGs' profit.To generate the power, DGs spend generation costs such as the maintenance cost, fuel cost, and emission penalty.Thus, we consider the benefit in terms of the profit when the generated power is maximally assigned to load demands.In detail, our algorithm shows about 159%, 141%, and 123% better benefits for 50%, 65%, and 80% test scenarios.Thus, as the power shortage becomes worse, the benefit of our algorithm increases in terms of load satisfaction and DGs' profit.
Figure 8 shows the running time of our load shedding algorithm.The dominant factor for the running time of our system is the time required to run the proposed optimization algorithm.For the efficient optimization calculation, we adopted the IBM CPLEX Optimizer in our MAS system.For our experiments, we varied the amount of total power supply.The running time of our system was up to 30 ms in the three test scenarios.In detail, the running times of our system were about 13.9 ms, 12.4 ms, and 11.3 ms for 50%, 65% and 80% test scenarios.Thus, we show that our optimization algorithm does not deteriorate performance in terms of running time.Figure 7 shows benefits for loads and DGs for different load conditions.On the load side, because the load demand is shed by force, we consider the benefits in terms of the satisfaction when the subdemands with high preferences are met.Figure 7a shows the benefit of load satisfaction when the SM of our algorithm is compared with the SM of the max-min algorithm.In detail, our algorithm shows about 153%, 132%, and 116% better benefits for 50%, 65%, and 80% test scenarios.Figure 7b shows the benefit of DGs' profit.To generate the power, DGs spend generation costs such as the maintenance cost, fuel cost, and emission penalty.Thus, we consider the benefit in terms of the profit when the generated power is maximally assigned to load demands.In detail, our algorithm shows about 159%, 141%, and 123% better benefits for 50%, 65%, and 80% test scenarios.Thus, as the power shortage becomes worse, the benefit of our algorithm increases in terms of load satisfaction and DGs' profit.
Figure 8 shows the running time of our load shedding algorithm.The dominant factor for the running time of our system is the time required to run the proposed optimization algorithm.For the efficient optimization calculation, we adopted the IBM CPLEX Optimizer in our MAS system.For our experiments, we varied the amount of total power supply.The running time of our system was up to Energies 2017, 10, 45 11 of 13 30 ms in the three test scenarios.In detail, the running times of our system were about 13.9 ms, 12.4 ms, and 11.3 ms for 50%, 65% and 80% test scenarios.Thus, we show that our optimization algorithm does not deteriorate performance in terms of running time.

Conclusions
We developed an optimization algorithm for the load shedding problem in an islanded MG.We achieved the system optimum by maximizing the satisfaction of both the DG and load.To maximize satisfaction, we designed the utility function to reflect the satisfaction of power allocation and the dissatisfaction of power shortage.To verify the performance, we implemented an MAS for the load shedding algorithm on top of a JADE framework.Regarding the adaptability of our algorithm, we showed that, as the power demands of the load changes, the number of shed demands accordingly

Conclusions
We developed an optimization algorithm for the load shedding problem in an islanded MG.We achieved the system optimum by maximizing the satisfaction of both the DG and load.To maximize satisfaction, we designed the utility function to reflect the satisfaction of power allocation and the dissatisfaction of power shortage.To verify the performance, we implemented an MAS for the load shedding algorithm on top of a JADE framework.Regarding the adaptability of our algorithm, we

Conclusions
We developed an optimization algorithm for the load shedding problem in an islanded MG.We achieved the system optimum by maximizing the satisfaction of both the DG and load.To maximize satisfaction, we designed the utility function to reflect the satisfaction of power allocation and the dissatisfaction of power shortage.To verify the performance, we implemented an MAS for the load shedding algorithm on top of a JADE framework.Regarding the adaptability of our algorithm, we showed that, as the power demands of the load changes, the number of shed demands accordingly changes.To verify performance, we compared the proposed algorithm with the max-min fairness algorithm in terms of SM.To analyze the feasibility, we measured the running time of our algorithm and revealed that the optimization problem is solved quickly and efficiently.For a future study, we plan to identify an optimal strategy for the load to determine the preference of subdemands for maximizing the profit of the load in the long run.In addition, we plan to continue to do research on how to implement our system in a real world system.For example, the conventional distribution systems are changed into multi-microgrid systems that are usually faster, more controllable and reliable.We plan to investigate the energy management system for a multi-microgrid.

Figure 1 .
Figure 1.Message flow among agents for our load shedding system.

Figure 1 .
Figure 1.Message flow among agents for our load shedding system.

Figure 2 .
Figure 2. Modified Consortium for Electric Reliability Technology Solutions (CERTS) MG.

Figure 2 .
Figure 2. Modified Consortium for Electric Reliability Technology Solutions (CERTS) MG.

Figure 3 .
Figure 3. Implemented system in a distributed environment.(a) Test environment of our laboratory; and (b) operational process on the implemented system.

Figure 3 .
Figure 3. Implemented system in a distributed environment.(a) Test environment of our laboratory; and (b) operational process on the implemented system.

Figure 4 .
Figure 4. Number of shed demands of a specific in three test scenarios.(a) When the 50% test scenario is used; (b) when the 65% test scenario is used; and (c) when the 80% test scenario is used.

Figure 4 .
Figure 4. Number of shed demands of a specific Ag load in three test scenarios.(a) When the 50% test scenario is used; (b) when the 65% test scenario is used; and (c) when the 80% test scenario is used.

Figure 4 .
Figure 4. Number of shed demands of a specific in three test scenarios.(a) When the 50% test scenario is used; (b) when the 65% test scenario is used; and (c) when the 80% test scenario is used.

Figure 8 .
Figure 8.Running time of our load shedding system.

Figure 8 .
Figure 8.Running time of our load shedding system.

Figure 8 .
Figure 8.Running time of our load shedding system.

Table 1 .
Summary of important symbols.DG: distributed generator.

Table 2 .
Details of load data for modified Consortium for Electric Reliability Technology Solutions (CERTS) MG.

Table 2 .
Details of load data for modified Consortium for Electric Reliability Technology Solutions (CERTS) MG.