Energy Management of Virtual Power Plant Considering Distributed Generation Sizing and Pricing

: The energy management of virtual power plants faces some fundamental challenges that make it complicated compared to conventional power plants, such as uncertainty in production, consumption, energy price, and availability of network components. Continuous monitoring and scaling of network gain status, using smart grids provides valuable instantaneous information about network conditions such as production, consumption, power lines, and network availability. Therefore, by creating a bidirectional communication between the energy management system and the grid users such as producers or energy applicants, it will a ﬀ ord a suitable platform to develop more e ﬃ cient vector of the virtual power plant. The paper is treated with optimal sizing of DG units and the price of their electricity sales to achieve security issues and other technical considerations in the system. The ultimate goal in this study to determine the active demand power required to increase system loading capability and to withstand disturbances. The e ﬀ ect of di ﬀ erent types of DG units in simulations is considered and then the e ﬃ ciency of each equipment such as converters, wind turbines, electrolyzers, etc., is achieved to minimize the total operation cost and losses, improve voltage proﬁles, and address other security issues and reliability. The simulations are done in three cases and compared with HOMER software to validate the ability of proposed model.


Introduction
For the first time, the concept of virtual power plant was presented in 1994, with the aims of distributed energy resources (DERs) observation, providing suitable interface for local components, activation of distributed control strategies, and optimal use of available capacity. In the distributed generation (DG) units, the load response and energy storage systems (ESSs) which are used as the unit entity, are called the 'virtual power plant' (VPP). According to the advantages of distributed generation resources and ESSs, a VPP can be an appropriate alternative for conventional fossil fuels. The need for modifying, changing, and displacing of energy consumption provides an efficient atmosphere to implement the VPPs [1,2]. Therefore, a VPP can be defined as a cluster of dispersed generating units, flexible loads, and storage systems that are grouped in order to operate as a single entity. The generating units in the VPP can employ both fossil and renewable energy sources. loads, distributed discharges equipment, wind farms, and storage systems connected through a grid. Uncertainties in the price of electrical energy and the production of wind farms have been solved by robust optimization. Linear programming and direct optimal load flow have been used to implement the energy management algorithm.
In [24], the impact of forecasting accuracy and technical constraints on the energy management of VPP with and without the presence of smart grid has been investigated. The constraints of the storage system, stochastic generation of distributed resources, energy transfer constraints, and responsive loads are considered in this model. In order to model the uncertainties in electrical energy prices and generation of solar stations, point prediction method has been used [25]. In this study, the possibility of contingencies in power management of the VPPs has not been discussed.
In this paper, the contributions are listed as below: • In this manuscript, with the development of the VPP model, an optimal formulation is proposed to manage the VPPs energy scheduling.

•
The point prediction method is used to model non-deterministic variables with uncertainties and PSO methodology is proposed to solve the energy management concerns.

•
Taking into account that the possibility of contingency occurrence and equipping the microgrid with the necessary instruments to monitor the situation quickly, the impact of those phenomena in a smart microgrid is investigated on the management process.

•
Proposing a deterministic model for VPP energy management to monitor the contingencies of microgrid.

•
The sizing and siting of VPP components are investigated under different scenarios with reliability evaluation.

Concept of Energy Management Problem
The VPPs consist of a set of DGs and responsive loads in a microgrid equipped with smart grid technology. Responsive loads as energy applicants, such as an industrial estate, colleges of a university, or a residential area can change their energy consumption in response to change in energy prices. The energy supply unit is responsible for energy management and can provide the energy required by the main grid and solar stations as a source of uncertain DGs. Beside each of the solar stations, there is an energy storage system that enables the energy management system to store or produce energy at appropriate times. The set of responsive loads, DERs and storage systems in the form of VPP, have the ability to buy energy from the main grid during low-cost energy hours and sell them at expensive periods [26,27]. It is assumed that the VPP does not play a role in determining market price. Energy management system is related to wholesale electricity market and energy price announced by the main network, is considered as market price. The VPP owns solar stations and does not pay for solar energy.
Using smart grid technology, there is a bidirectional relation between time and reality of energy management and users in the microgrid. The energy management system will receive information about energy suppliers, including the energy price and the availability of obtainable power several minutes before the energy delivery date. Based on this information, the energy management system determines the optimal consumption rate, storage, and energy exchange for each time period and subsequent periods to send decisions a few minutes before the deadline for energy delivery to applicants and energy suppliers. These decisions are made in such that the social welfare of the VPP will include the difference in the cost of energy supply from the income of the applicants. Notice that information about the availability of network components including production resources and transmission lines is provided to the energy management system according to the scheduling period. If anything else happens, the energy management system will modify the decisions taken in the previous step to distribute power. The microgrid is also assumed to run at the beginning of the planning period. Figure 1 shows the proposed microgrid structure that interacts with the electricity market. As shown in this figure, the considered system includes wind turbine units, solar cells, electrolytic reservoirs, hydrogen fuel storage tanks, fuel cells, batteries, and Direct current (DC)/Alternative current (AC) converters and loads (in the proposed model, the load types are considered as non-interruptible). In an autonomous microgrid, distributed energy sources should be able to balance their production and consumption with sufficient capacity and proper operation, while islanded from the upstream network. In the proposed model of this paper, resources will only sell energy on the day-ahead market to the upstream distribution network due to the compatibility with the objective function target (cost minimization/profit maximization).  Figure 1 shows the proposed microgrid structure that interacts with the electricity market. As shown in this figure, the considered system includes wind turbine units, solar cells, electrolytic reservoirs, hydrogen fuel storage tanks, fuel cells, batteries, and Direct current (DC) / Alternative current (AC) converters and loads (in the proposed model, the load types are considered as non-interruptible). In an autonomous microgrid, distributed energy sources should be able to balance their production and consumption with sufficient capacity and proper operation, while islanded from the upstream network. In the proposed model of this paper, resources will only sell energy on the day-ahead market to the upstream distribution network due to the compatibility with the objective function target (cost minimization/profit maximization). Where in the above figure we have:

Mathematical Formulation
• P :Output power of each solar unit (kW) • :Output power of each wind turbine (kW) • :Total power generated by renewable energy units (kW) • :Transmission power from renewable sources to electrolyzer (kW) • :Transmission power from renewable sources to battery (kW) • :Transmission power from renewable sources and battery to DC/AC converter (kW) • :Transmission power from electrolyzer to hydrogen tank (kW) • :Transmission power from battery to DC/AC converter (kW) • :Transmission power from hydrogen tank to fuel cell (kW) • :Transmission power from fuel cell to DC/AC converter (kW) • :Transmission power from DC/AC converter to load and upstream network (kW) • P :Wholesale power sold to up-stream network (kW)

Solar Unit
The data received from the sun to the array surface is converted to its output power using the equation  Where in the above figure we have: • P PV : Output power of each solar unit (kW) • P WT : Output power of each wind turbine (kW) • P Gen : Total power generated by renewable energy units (kW) • P Gen-ele : Transmission power from renewable sources to electrolyzer (kW) • P Gen-bat : Transmission power from renewable sources to battery (kW) • P Gen-inv : Transmission power from renewable sources and battery to DC/AC converter (kW) • P ele-Tank : Transmission power from electrolyzer to hydrogen tank (kW) • P bat-inv : Transmission power from battery to DC/AC converter (kW) • P tank-FC : Transmission power from hydrogen tank to fuel cell (kW) • P FC-inv : Transmission power from fuel cell to DC/AC converter (kW) • P inv-load : Transmission power from DC/AC converter to load and upstream network (kW) • P sale : Wholesale power sold to up-stream network (kW)

Solar Unit
The data received from the sun to the array surface is converted to its output power using the equation In Equation (1), G is the radiation power perpendicular to the array surface (W/m 2 ) and P PV,rated represents the nominal power of each array, which is obtained for the cost. η PV,conv is also equivalent to the efficiency of the DC/DC converter installed between each DC array and the corresponding bus. By verifying the vertical and horizontal component of the solar power, at any moment, the power can be transmitted (vertically) onto the surface of the array installed with the angle θ PV could be calculated according to (2) where G H (t) and G V (t) are the horizontal and vertical radiation rates in terms of (W/m 2 ), respectively.

Wind Turbine Unit
The wind turbine class used in this paper is BWC Excel-R/48. Output power (P WT ) in terms of wind speed (v W ) can be approximated by relation (3) [17] where v cutin , v cutout , and v rated are the cut-off wind speed, the high cut-off rate and the nominal speed (m/s) of the turbine, respectively. The parameter P WT,max , shows the maximum turbine output power (kW) and P f url , also have an output power at the high cut-off speed. In this paper, m is considered to be (3).

Electrolyzer Unit
The function of electrolyzer, is decomposition of water into hydrogen and oxygen by a simple electrolysis process. In this way, the direct current of electricity passes through the path between two electrodes that immersed in water, causing water to decompose into oxygen and hydrogen. Oxygen is produced in the cathode and hydrogen at the anode side. In the proposed model, this paper is used to reduce the energy consumption of a compressor-less design, due to the usage of a polymer electrolyte membrane (PEM) fuel cell. In view of the fact that in this type of pressure fluid, the pressure required by hydrogen is 1.2 bar, the developed is very flexible and can be easily added to the compressor model. The electrochemical interactions in the electrolysis apparatus of water are To model the electrolyzer, its efficiency is used as an input parameter. The thermal value of hydrogen is 3.4 kV per cubic meter, which, taking into account the efficiency of 90% for an electrolyzer to produce a kilogram of hydrogen, consumes energy at 41.99 kWh/kg as (5).

90
× 100 = 41.97kwh/kg (5) The hydrogen weight produced by the excess of energy produced from the system to the electrolyzeris obtained at 41.97 kWh/kg.

Fuel Cell Unit
PEMs have a relatively fast dynamic response, about 1 to 3 s. The power output of these fuel cells can be calculated as a function of the input power of the hydrogen as well as its efficiency (η FC ), which can be assumed to be constant. Therefore the output power extracted from fuel cell stacks (P FC-inv ) could be represented with (7) in which P tank-FC is gross productive power of fuel cells.

Hydrogen Tank Unit
The energy stored in the tank for each step t, can be calculated as where, ∆t is the length of each time step, P el_tank , represents the transmission power from the electrolyzer to the hydrogen tank and P tank_FC , can be transmitted from the hydrogen tank to the fuel cell. The term η storage also signifies the efficiency of the storage system, which can indicate losses due to leakage or pumping. However, the hydrogen in the tank will always have a high and low range limitations, as

Energy Storage (Battery) Unit
The battery source is used to provide the load in the absence of renewable energy sources. The difference between the power produced and the load power required indicates whether the battery should be charged or discharged. The amount of charge of the battery bank is obtained in time horizon t using the following Equation (10) where, E bat (t) represents the amount of battery electric energy at time t. By the way, η bat and η dis_bat , are the charge and discharge efficiency of the battery bank, respectively.

Reliability
There are several references which have provided several indicators for calculating the reliability of systems, including indicators such as loss of load expectation (LOLE), loss of energy expectation (LOEE), or expected energy not supplied (EENS), loss of power supply probability (LPSP), equivalent loss factor (ELF), and so on. The above indicators are defined by the following relationships. In this paper, Markov chain method is used to calculate the reliability parameters. Other required parameters are also derived from [28][29][30].
In the above relation, E[LOL(t)] is the mathematical expectation of the time offset in time t, which can be defined by In this case, P s is the probability of being in the state of s and T s in the case of being in this situation. S is the total set of possible situations for the system. Therefore Here E[LOE(t)] is the mathematical expectation of the amount of lost energy in the time interval t that can be defined by Q s is the amount of lost load (kWh) if positioned in s. The probability of loss of source (LPSP) is obtained according to (15).
In the above relation, it is assumed that D(t) is equal to the load demand (kWh) in time t. Finally, the equivalent load offset coefficient can be defined as (16) Since the ELF contains more information, both the number of determinations N (and the values of Q(t)), this paper is used as the main criterion of reliability. The maximum allowed for ELF in developed countries is equal to 0.0001, however, for an independent system to the network, this limit is considered to be 0.01.

Objective Function
Life cycle cost analysis evaluates the costs of covering all expenditures incurred during the activity period. The net present cost (NPC) is used as the charge of the system life cycle. The NPC includes initial installation costs, replacement costs, repairs, and maintenance of the equipment, the cost of the power failure, the cost of connecting to the grid. In NPC calculations, the costs are considered positive and earnings are considered negative.
All costs and expenses are assessed at a fixed interest rate throughout the year. In this type of assessment, in order to influence the increase rate in calculations at the end of the process of analysis and review of the system, it should be applied to the NPC, by calculating the real interest rate caused by inflation, according to certain relationships. The real interest rate is the difference between the nominal interest rate and the inflation rate. The NPC value of the equipment can be calculated according to the equation In the above statement, N will equip the unit, or capacity (kW or kg), CC is initial investment cost ($/unit), RC stands for cost of each replacement annual maintenance cost ($/unit-year) equipping at R project lifetime (in this study is 20 years). The cost of the initial purchase of hydrogen is considered at the cost of the tank investment and represents real interest, which can be calculated in terms of the nominal interest (ir nominal ) and the annual inflation rate (ir) in accordance with Equation (18). PWA and K are respectively the annual and constant payments current value, which are defined as y and L are the number of replacements and useful life of the equipment, respectively. For each component of the system, the initial cost, in accordance with the cost in year zero, is the replacement cost according to the need to replace that component at the end of its lifetime and the cost of O & MC for each year of the project. The cost of replacing a component varies for several reasons with its initial cost. One of these reasons is that when replacing a component due to the expiration of its life-span, all the components of that component that have been expended during the initial installation do not need to be replaced.

Interruption Power Cost
The interruption cost of electricity equals the cost of damages caused by power outages to consumers. The cost of interrupting the supply of electrical energy required by load is estimated by different methods. For example, this can be calculated based on the willingness of the customer to pay for the expansion of the network or the losses incurred by the industries due to the interruption created in their production process. If the amount of lost hope is defined annually with relation (21), Then, the net present value of the load loss can be obtained according to (22) In fact, the Equation (22), is equal to the average loss due to the disconnection of every kW hour of charge ($/kWh). Regarding the type of load (non-removable, etc.) mentioned previously, the average loss due to the interruption of each kWh is different.

Power Sailing Revenue
Since the goal is to minimize the operation costs of objective function, we introduce negative revenue from the sale of electricity to the upstream network in the calculations. The net present value of electricity sales to the upstream network is In the above-mentioned equation, C sale is the revenues generated by selling electricity power in (kW/hour) to the upstream network, while it depends on the time of the power exchange and the price of energy at that period. Given the costs and income mentioned above, the objective function is defined as (24).
where i represents the desired equipment and x is a vector of optimization variables.

Constraints
At any given interval, the total production capacity of the hybrid production system should be able to meet the total demand for the load, taking into account the terms and conditions of reliability, which is calculated by the equation Accordingly, P Load (t), P sale (t), P Gen_inv (t), P FC_inv (t) represent the total load demand (disconnect-able and non-interruptible), high power upstream network, the power transmitted from DGs to battery and the transient power from fuel cell stacks to the DC/AC converter, respectively. Since 10% of the total demand per hour is considered as a removable load and the remaining 90% is assumed as an irremovable load, Therefore, when (Loss_power (t)/P Load (t) < 0.1), load interruptions are not included in the reliability calculations of the system. Consequently, when (Loss_power (t)/P Load (t) > 0.1), the ELF index in the reliability calculations of the system is computed. Thus, until the Equation (33) is correct, the load procurement can be continued. (26) and the capacity to be sold to the upstream network should not exceed a certain limit, which is determined by prior agreements The energy stored in the hydrogen tank and the battery should be within the following limits: The last two constraints imply that the energy stored in the tank and the battery at the end of the year should not be less than the energy stored at the beginning of the year. This ensures that the reliability calculations are performed for the worst possible situation.
The proposed formulation is optimized using particle swarm optimization (PSO) algorithm. The capacity of DGs is constant while the number of DGs is variable. In the other words our variables should be optimized are the number of DGs using PSO. By optimizing the variable, the objective function is concluded which is the total cost of VPP. The interested readers are referred to [31] for more studying about PSO.

Optimization Algorithm
Particle swarm optimization is an evolutionary and population-based algorithm that presented introduced by Eberhart and Kennedy [32]. Each particle in PSO is a candidate solution in search space of the problem. There are two main parts in these candidate solutions-current position (X i ) and current velocity (V i )-described as where n is the dimension of problem (solution) and k is the iteration index. New position of each particle is updated using current position and new velocity. New velocity also is produced by four factors including: current velocity, current position, best previous position of particle (Pbest), and best position in all particles of all iterations (GBest). Then, the new velocity is represented as where ω is particle inertia coefficient, c 1 and c 2 are accelerations coefficient, respectively. r 1 and r 2 are uniform random numbers between 0 and 1. New positions are obtained as where pbest i,j , is the particle best of position i with dimension j and gbest j is the dimension j of G best . The PSO algorithm is depicted in Figure 2.
, ,..., , ,..., Wheren is the dimension of problem (solution) and k is the iteration index. New position of each particle is updated using current position and new velocity. New velocity also is produced by four factors including: current velocity, current position, best previous position of particle (Pbest), and best position in all particles of all iterations (GBest). Then, the new velocity is represented as where ω is particle inertia coefficient, c1 and c2 are accelerations coefficient, respectively. r1 and r2 are uniform random numbers between 0 and 1. New positions are obtained as where i , j pbest , is the particle best of position i with dimension j and j gbest is the dimension j of Gbest.The PSO algorithm is depicted in Figure 2.

Simulation Results and Discussion
The capacity of each wind turbine is considered as 7.5 kW, 1 kW for each solar array, power of each electrolyzer and fuel cell are measured to 1 kW, then the size of each hydrogen tank is 1 Kg and the amount of energy stored per battery is 9.6 kWh. Tables 1 and 2 contain the specifications of the

Simulation Results and Discussion
The capacity of each wind turbine is considered as 7.5 kW, 1 kW for each solar array, power of each electrolyzer and fuel cell are measured to 1 kW, then the size of each hydrogen tank is 1 Kg and the amount of energy stored per battery is 9.6 kWh. Tables 1 and 2 contain the specifications of the system equipment, which are modeled. Also, the real interest rate (u) is 0.08 ir. Considering the useful life of the project is 20 years, the PWA (ir, R) value is 9.818. The maximum power available at upstream network is 500 kW. Also, the energy prices sold to the high-voltage network during different hours of the day are from the spring and autumn periods are those written in accordance with Table 3. The seasonal coefficient of summer and winter is defined as 1.3 and 0.8. The interruption cost for a removable load equivalent is 0.1 ($/kWh) and for non-removable loads equivalent is 0.5 ($/kWh). The price of electricity sales to domestic grid users (which is considered as the maximum profit for these prices) is at different times according to Table 4.
Transmission lines and transformer investment costs are 5% of equipment costs, as well as miscellaneous network costs. The cost cannot be evaluated precisely such as financial transaction fees, accounting costs, transportation, advertising and marketing, etc. Then there is 38% of the cost of installing scattered units and the cost of repairing and maintaining grid equipment, and 5% of the cost of repair and maintenance of dispersed production units. In this section, in order to validate the results of the developed software, initially, the results of the size and optimization of the three hybrid systems as shown in  They are compared with the output results extracted from the HOMER software. Hybrid systems consist of solar panels, wind turbines, and energy storage systems. In the hybrid system number 1, the hydrogen tank, in system 2, the battery and in system number 3, the battery and hydrogen tank are used as an energy storage device. systems. In the hybrid system number 1, the hydrogen tank, in system 2, the battery and in system number 3, the battery and hydrogen tank are used as an energy storage device.     The results obtained from HOMER software are compared with the results obtained from the software developed in MATLAB environment, shown in Table 5.  According to Table 5, the results obtained by HOMER software are very close to the results obtained from the developed software and confirm the accuracy of the developed coding. Also, the results of the review of the three hybrid systems show that the use of battery as a storage device instead of an electrolyzer assembly, hydrogen tank, and fuel cell will increase system costs. Therefore, this price rising is justified by the high speed of the battery in pursuit of charge. Table 6, also shows the optimal number of energy resources extended by considering the virtual power plant, the amount of revenue, and costs of the proposed grid.  The results obtained from HOMER software are compared with the results obtained from the software developed in MATLAB environment, shown in Table 5.  According to Table 5, the results obtained by HOMER software are very close to the results obtained from the developed software and confirm the accuracy of the developed coding. Also, the results of the review of the three hybrid systems show that the use of battery as a storage device instead of an electrolyzer assembly, hydrogen tank, and fuel cell will increase system costs. Therefore, this price rising is justified by the high speed of the battery in pursuit of charge. Table 6, also shows the optimal number of energy resources extended by considering the virtual power plant, the amount of revenue, and costs of the proposed grid. The results obtained from HOMER software are compared with the results obtained from the software developed in MATLAB environment, shown in Table 5. According to Table 5, the results obtained by HOMER software are very close to the results obtained from the developed software and confirm the accuracy of the developed coding. Also, the results of the review of the three hybrid systems show that the use of battery as a storage device instead of an electrolyzer assembly, hydrogen tank, and fuel cell will increase system costs. Therefore, this price rising is justified by the high speed of the battery in pursuit of charge. Table 6, also shows the optimal number of energy resources extended by considering the virtual power plant, the amount of revenue, and costs of the proposed grid. Now, considering the uncertainty in wind energy, for the 95% confidence level, we solve the problem again with the result of optimization presented in Table 6. then, Figure 6 represents the total cost of the system in terms of the number of repeat steps.  Now, considering the uncertainty in wind energy, for the 95% confidence level, we solve the problem again with the result of optimization presented in Table 6. then, Figure 6 represents the total cost of the system in terms of the number of repeat steps.  Table 7 shows the optimal number of sources of energy, the amount of income and expenses of the proposed grid.   Table 7 shows the optimal number of sources of energy, the amount of income and expenses of the proposed grid. The output results show that the inclusion of uncertainty in the output power of wind units increases the cost, and this growth in costs is justified by increasing the reliability level of wind turbine output, which increases the level of reliability in the output of wind turbines, correspondingly. The average reliability of the system will be reduced from 0.053 to 0.042 during the year, or increase the reliability of the system. Now we are examining and simulating system scenarios.

Cut-Off Wind Effect
In Table 8, the effect of cut-off wind turbine speed on the optimum size and cost of the proposed grid is shown. As it can be seen, while the number of wind turbines with low cut-off speed is decreased, the total cost increases by growing the cut-off speed. This indicates that more wind turbines working with low cut-off speed means less power will be generated. Therefore, other DGs have to operate at a higher power level and this raises the cost of operation. It is because the input of wind turbines is cheap wind energy, while some of those DGs use hydrogen fuel or heat as input power.

PV Investment
In Table 9, the effect of initial investment cost of solar panels on the optimum size and cost of the system is presented. According to Table 9, as the investment costs increase, the number of solar units and prices decrease. This is because the cheap and common energy of the sun will be less used and should be directed towards other energies like wind. Since the capacity of wind resources is also limited, the utilization of costly resources is required, which results in increasing the costs of exploitation.

DC/AC Converter
In Table 10, the effect of the DC/AC converter efficiency on the optimal size and cost of system execution is shown. According to Table 10, it is expected to decrease costs by increasing the efficiency of the converter. As the converter efficiency rises, the power losses are reduced, correspondingly. Then, with lower energy production, all demands will be supplied sufficiently.

Effect of Selling Electricity to the Upstream Network
In Table 11, the effect of the electricity exchanging price to the upstream network in non-peak hours is represented. This issue has been obtained in terms of the optimal size assessment of DGs considering grid connected operation mode. According to Table 11, the increase in electricity sales prices to the upstream network will increase the sales of electricity to the upstream network, which will reduce the total operation cost of system and it will raise the revenue.

Reliability Effects
In Table 12, the effect of increasing system reliability (ELF(t) or reliability index reduction) has been shown on the optimum size and cost of performing grid. According to Tables 13 and 14, increasing system reliability increases system costs, as expected. As a total deduction in this paper, solar panels and wind turbines are used as power generators and from hydrogen batteries and tanks as energy storage. The proposed system consists of two types of loads which are interruptible loads such as little residential areas and non-interruptible loads like the hospitals. Then, the expression of the mathematical model of each of the components of the proposed micro-network energy system was presented. After expressing the model of the proposed grid elements and components, the objective function was proposed to minimize (maximize) 20-year costs of the system along with its constraints.
The proposed algorithm and some the state of art methods are applied on system No. 1, and compared as shown in Tables 13 and 14. The parameter considered for evaluation consists of total operation cost, expected energy not supplied (EENS), loss of power supply probability (LPSP), and the numbers of DGs crate the VPP. As it is observed, the proposed method has the lowest operation cost and the best reliability indices among other approaches.

Conclusions
From the studies made in this paper, the following results can be deduced: • This major weakness greatly reduces the reliability of energy systems. The solution to this problem is the use of supportive production systems or energy storage systems, in which the battery and hydrogen tank are used as a storage system in order to provide optimal reliability.
• Compared to the use of storage systems, combining different energy sources that have complementary production characteristics (such as wind and sun) are considered as a convenient and inexpensive way to improve system reliability. The proposed grid has such a structure to overcome the above problems.

•
Comparison of the results of determining the optimal capacity of the three hybrid systems suggests that using the battery as a storage device in a hybrid system instead of a fuel cell would increase system costs which increase the cost of the system by increasing the system's ability to track the load.

•
The results show that simultaneous use of fuel and fuel cells as a storage medium in the hybrid system, in addition to reducing costs, also increases the system's ability to track the load. In addition, considering the disruptive load model as a virtual power plant reduces costs.

•
The results of determining the size of resources distributed energy microgrid suggested that in view of the uncertainty in the power output of wind turbines increased costs (loss of income) of the system during its life, is that the increased cost (reduced revenue) is justified by increasing the level of reliability in the power output of wind turbines.

•
Ability to cut off removable loads in case of necessity increases the reliability of the system and (possibly) reduces the cost of the grid.